Novel force-displacement control passive finite element models of the spine to simulate intact and pathological conditions; comparisons with traditional passive and detailed musculoskeletal models

Evaluation of ground reaction forces and centers of pressure predicted by AnyBody Modeling System during load reaching/handling activities and effects of the prediction errors on model-estimated spinal loads

Full-body and lower-extremity human musculoskeletal models require feet ground reaction forces (GRFs) and centers of pressure (CoPs) as inputs to predict muscle forces and joint loads. GRFs/CoPs are traditionally measured via floor-mounted forceplates that are usually restricted to research laboratories thus limiting their applicability in real occupational and clinical setups. Alternatively, GRFs/CoPs can be estimated via inverse dynamic approaches as also implemented in the Anybody Modeling System (AnyBody Technology, Aalborg, Denmark). The accuracy of Anybody in estimating GRFs/CoPs during load-handling/reaching activities and the effect of its prediction errors on model-estimated spinal loads remain to be investigated. Twelve normal- and over-weight individuals performed total of 480 static load-handling/reaching activities while measuring (by forceplates) and predicting (by AnyBody) their GRFs/CoPs. Moreover, the effects of GRF/CoP prediction errors on the estimated spinal loads were evaluated by inputting measured or predicted GRFs/CoPs into subject-specific musculoskeletal models. Regardless of the subject groups (normal-weight or overweight) and tasks (load-reaching or load-handling), results indicated great agreements between the measured and predicted GRFs (normalized root-mean-squared error, nRMSEs < 14% and R2 > 0.90) and between their model-estimated spinal loads (nRMSEs < 14% and R2 > 0.83). These agreements were good but relatively less satisfactory for CoPs (nRMSEs < 17% and 0.57 < R2 < 0.68). The only exception, requiring a more throughout investigation, was the situation when the ground-foot contact was significantly reduced during the activity. It appears that occupational/clinical investigations performed in real workstation/clinical setups with no access to forceplates may benefit from the AnyBody GRF/CoP prediction tools for a wide range of load-reaching/handling activities.

Muscle-driven forward dynamic active hybrid model of the lumbosacral spine: combined FEM and multibody simulation

Most spine models belong to either the musculoskeletal multibody (MB) or finite element (FE) method. Recently, coupling of MB and FE models has increasingly been used to combine advantages of both methods. Active hybrid FE-MB models, still rarely used in spine research, avoid the interface and convergence problems associated with model coupling. They provide the inherent ability to account for the full interplay of passive and active mechanisms for spinal stability. In this paper, we developed and validated a novel muscle-driven forward dynamic active hybrid FE-MB model of the lumbosacral spine (LSS) in ArtiSynth to simultaneously calculate muscle activation patterns, vertebral movements, and internal mechanical loads. The model consisted of the rigid vertebrae L1-S1 interconnected with hyperelastic fiber-reinforced FE intervertebral discs, ligaments, facet joints, and force actuators representing the muscles. Morphological muscle data were implemented via a semi-automated registration procedure. Four auxiliary bodies were utilized to describe non-linear muscle paths by wrapping and attaching the anterior abdominal muscles. This included an abdominal plate whose kinematics was optimized using motion capture data from upper body movements. Intra-abdominal pressure was calculated from the forces of the abdominal muscles compressing the abdominal cavity. For the muscle-driven approach, forward dynamics assisted data tracking was used to predict muscle activation patterns that generate spinal postures and balance the spine without prescribing accurate spinal kinematics. During calibration, the maximum specific muscle tension and spinal rhythms resulting from the model dynamics were evaluated. To validate the model, load cases were simulated from -10° extension to +30° flexion with weights up to 20 kg in both hands. The biomechanical model responses were compared with in vivo literature data of intradiscal pressures, intra-abdominal pressures, and muscle activities. The results demonstrated high agreement with this data and highlight the advantages of active hybrid modeling for the LSS. Overall, this new self-contained tool provides a robust and efficient estimation of LSS biomechanical responses under in vivo similar loads, for example, to improve pain treatment by spinal stabilization therapies.

Biomechanical evaluation and optimal design of a pedicle screw with double bent rods internal fixation system based on PE-PLIF fusion

Problems, such as broken screws, broken rods, and cage subsidence after clinical spinal fusion surgery affect the success rate of fusion surgery and the fixation effect of fusion segments, and these problems still affect the treatment and postoperative recovery of patients. In this study, we used the biomechanical finite element analysis method to analyze and study the fixation effect of three kinds of spinal internal fixation systems on L4-L5 lumbar spine segments in percutaneous endoscopic posterior lumbar interbody fusion (PE-PLIF). The three different fixation systems compared in this study include bilateral pedicle screw fixation (M1); bilateral pedicle screw with cross-link fixation (M2); bilateral pedicle screws with double bent rods fixation (M3). The internal fixation systems with different structures were analyzed with the help of Hypermesh, and Abaqus. It was found that the internal fixation system with double bent rods reduced screw stresses by 23.8 and 22.2% in right and left axial rotation than the traditional bilateral pedicle screw system, while titanium rod stresses were reduced by 9.6, 3.7, 9.6, and 2.9% in flexion, left and right lateral bending, and right axial rotation, respectively, and L5 upper endplate stresses were reduced by 35.5, 18.9, 38.4, 10.2, and 48.3% in flexion, left and right lateral bending, and left and right axial rotation, respectively. The spinal range of motion (ROM) of the M3 internal fixation system was less than that of the M1 and M2 internal fixation systems in left lateral bending, left lateral rotation, and right axial rotation, and the intact vertebral ROM was reduced by 93.7, 94.9, and 90.9%, respectively. The double bent rod structure of the spinal internal fixation system has better biomechanical properties, which can effectively reduce the risk of screw breakage, loosening, cage subsidence, and endplate collapse after fusion surgery.