Sensitivity of kinematics-based model predictions to optimization criteria in static lifting tasks

Biomechanics of the Human Knee Joint in Maximum Voluntary Isometric Flexion: Study of Changes in Applied Moment, Agonist-Antagonist Participations, Joint Center, and Flexion Angle

Estimation of the knee joint strength by maximum voluntary isometric contraction (MVIC) is a common practice to assess strength, coordination, safety to return to work or engage in sports after an injury, and to evaluate the efficacy of treatment modalities and rehabilitation strategies. In this study, we utilize a previously validated coupled finite element-musculoskeletal model of the lower extremity to explore the sensitivity of output measures (posterior cruciate ligament [PCL]/muscle/contact forces and passive moments) in knee MVIC flexion exercises at seated position. To do so, at three knee flexion angles (KFA), input measures (resistance moment and contribution moments of quadriceps and gastrocnemii) were varied at four levels each using the Taguchi design of experiment. Our findings reveal significant increases in PCL forces with KFA (p < 0.01), net MVIC moment (p < 0.01), and resistance moment of quadriceps (p < 0.01). In contrast, they drop at larger activity in gastrocnemii (p < 0.01). Tibiofemoral (TF) contact forces increase with the net MVIC moment (p < 0.01). The passive knee flexion moment, while highly dependent on the location at which computed, also increases with the net MVIC moment (p < 0.01). Changes in KFA, MVIC moment, and proportions thereof carried by quadriceps and/or gastrocnemii substantially affect biomechanics of the joint. Compared with level walking and stair ascent, slightly larger contact forces/stresses and much greater PCL forces are computed. This study improves our understanding of the knee joint behavior during MVIC in effective evaluation and rehabilitation interventions. Besides, it emphasizes the importance of positioning the joint center in model studies.

The ankle dorsiflexion kinetics demand to increase swing phase foot-ground clearance: implications for assistive device design and energy demands

BACKGROUND

The ankle is usually highly effective in modulating the swing foot's trajectory to ensure safe ground clearance but there are few reports of ankle kinetics and mechanical energy exchange during the gait cycle swing phase. Previous work has investigated ankle swing mechanics during normal walking but with developments in devices providing dorsiflexion assistance, it is now essential to understand the minimal kinetic requirements for increasing ankle dorsiflexion, particularly for devices employing energy harvesting or utilizing lighter and lower power energy sources or actuators.

METHODS

Using a real-time treadmill-walking biofeedback technique, swing phase ankle dorsiflexion was experimentally controlled to increase foot-ground clearance by 4 cm achieved via increased ankle dorsiflexion. Swing phase ankle moments and dorsiflexor muscle forces were estimated using AnyBody modeling system. It was hypothesized that increasing foot-ground clearance by 4 cm, employing only the ankle joint, would require significantly higher dorsiflexion moments and muscle forces than a normal walking control condition.

RESULTS

Results did not confirm significantly increased ankle moments with augmented dorsiflexion, with 0.02 N.m/kg at toe-off reducing to zero by the end of swing. Tibialis Anterior muscle force incremented significantly from 2 to 4 N/kg after toe-off, due to coactivation with the Soleus. To ensure an additional 4 cm mid swing foot-ground clearance, an estimated additional 0.003 Joules/kg is required to be released immediately after toe-off.

CONCLUSION

This study highlights the interplay between ankle moments, muscle forces, and energy demands during swing phase ankle dorsiflexion, offering insights for the design of ankle assistive technologies. External devices do not need to deliver significantly greater ankle moments to increase ankle dorsiflexion but, they should offer higher mechanical power to provide rapid bursts of energy to facilitate quick dorsiflexion transitions before reaching Minimum Foot Clearance event. Additionally, for ankle-related bio-inspired devices incorporating artificial muscles or humanoid robots that aim to replicate natural ankle biomechanics, the inclusion of supplementary Tibialis Anterior forces is crucial due to Tibialis Anterior and Soleus co-activation. These design strategies ensures that ankle assistive technologies are both effective and aligned with the biomechanical realities of human movement.

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.

Towards improving the accuracy of musculoskeletal simulation of dynamic three-dimensional spine rotations with optimizing model and algorithm

BACKGROUND

The accuracy of musculoskeletal simulations greatly relies on model structures and optimization algorithms. This study investigated the unclarified influence of accounting for several commonly-simplified different model components and optimization criteria on spinal musculoskeletal simulations.

METHODS

The study constructed a full-body musculoskeletal model with passive components of functional spinal units and spinal muscles subject-specifically refined. A muscle redundancy solver was built with 15 optimization criteria. Three-dimensional spine rotations and spinal muscle activities were measured using optical motion capture and electromyogram techniques when eight healthy volunteers performed standing, flexion/extension, lateral bending, and axial rotation. The effect of the model with four different conditions of the passive components and the sensitivity of the 15 optimization criteria on simulations were investigated.

RESULTS

Accounting for the refined passive components significantly improved the simulation accuracy. Different optimization criteria behaved distinctly for different motions. Generally minimizing the sum of squared muscle activations outperformed the others, with the highest averaged correlation coefficient (0.82) between the estimated erector spinae muscle activations and measured electromyography and with the estimated joint compression forces comparable to in vivo reference data.

CONCLUSION

This study highlights the importance of passive model components and proposes a suitable optimization framework for realistic spinal musculoskeletal simulations.

Biomechanical Properties of Paraspinal Muscles Influence Spinal Loading—A Musculoskeletal Simulation Study

Paraspinal muscles are vital to the functioning of the spine. Changes in muscle physiological cross-sectional area significantly affect spinal loading, but the importance of other muscle biomechanical properties remains unclear. This study explored the changes in spinal loading due to variation in five muscle biomechanical properties: passive stiffness, slack sarcomere length (SSL), in situ sarcomere length, specific tension, and pennation angle. An enhanced version of a musculoskeletal simulation model of the thoracolumbar spine with 210 muscle fascicles was used for this study and its predictions were validated for several tasks and multiple postures. Ranges of physiologically realistic values were selected for all five muscle parameters and their influence on L4-L5 intradiscal pressure (IDP) was investigated in standing and 36° flexion. We observed large changes in IDP due to changes in passive stiffness, SSL, in situ sarcomere length, and specific tension, often with interesting interplays between the parameters. For example, for upright standing, a change in stiffness value from one tenth to 10 times the baseline value increased the IDP only by 91% for the baseline model but by 945% when SSL was 0.4 μm shorter. Shorter SSL values and higher stiffnesses led to the largest increases in IDP. More changes were evident in flexion, as sarcomere lengths were longer in that posture and thus the passive curve is more influential. Our results highlight the importance of the muscle force-length curve and the parameters associated with it and motivate further experimental studies on in vivo measurement of those properties.

EMG optimization in OpenSim: A model for estimating lower back kinetics in gait

Participant-specific musculoskeletal models are needed to accurately estimate lower back internal kinetic demands and injury risk. In this study we developed the framework for incorporating an electromyography optimization (EMGopt) approach within OpenSim (https://simtk.org/projects/emg_opt_tool) and evaluated lower back demands estimated from the model during gait. Kinematic, external kinetic, and EMG data were recorded from six participants as they performed walking and carrying tasks on a treadmill. For evaluation, predicted lumbar vertebral joint forces were compared to those from a generic static optimization approach (SOpt) and to previous studies. Further, model-estimated muscle activations were compared to recorded EMG, and model sensitivity to day-to-day EMG variability was evaluated. Results showed the vertebral joint forces from the model were qualitatively similar in pattern and magnitude to literature reports. Compared to SOpt, the EMGopt approach predicted larger joint loads (p<.01) with muscle activations better matching individual participant EMG patterns. L5/S1 vertebral joint forces from EMGopt were sensitive to the expected variability of recorded EMG, but the magnitude of these differences (±4%) did not impact between-task comparisons. Despite limitations inherent to such models, the proposed musculoskeletal model and EMGopt approach appears well-suited for evaluating internal lower back demands during gait tasks.

Walking Biomechanics and Spine Loading in Patients With Symptomatic Lumbar Spinal Stenosis

Symptomatic lumbar spinal stenosis is a leading cause of pain and mobility limitation in older adults. It is clinically believed that patients with lumbar spinal stenosis adopt a flexed trunk posture or bend forward and alter their gait pattern to improve tolerance for walking. However, a biomechanical assessment of spine posture and motion during walking is broadly lacking in these patients. The purpose of this study was to evaluate lumbar spine and pelvic sagittal angles and lumbar spine compressive loads in standing and walking and to determine the effect of pain and neurogenic claudication symptoms in patients with symptomatic lumbar spinal stenosis. Seven participants with symptomatic lumbar spinal stenosis, aged 44–82, underwent a 3D opto-electronic motion analysis during standing and walking trials in asymptomatic and symptomatic states. Passive reflective marker clusters (four markers each) were attached to participants at T1, L1, and S2 levels of the spine, with additional reflective markers at other spinal levels, as well as the head, pelvis, and extremities. Whole-body motion data was collected during standing and walking trials in asymptomatic and symptomatic states. The results showed that the spine was slightly flexed during walking, but this was not affected by symptoms. Pelvic tilt was not different when symptoms were present, but suggests a possible effect of more forward tilt in both standing (p = 0.052) and walking (p = 0.075). Lumbar spine loading during symptomatic walking was increased by an average of 7% over asymptomatic walking (p = 0.001). Our results did not show increased spine flexion (adopting a trunk-flexed posture) and only indicate a trend for a small forward shift of the pelvis during both symptomatic walking and standing. This suggests that provocation of symptoms in these patients does not markedly affect their normal gait kinematics. The finding of increased spine loading with provocation of symptoms supports our hypothesis that spine loading plays a role in limiting walking function in patients with lumbar spinal stenosis, but additional work is needed to understand the biomechanical cause of this increase.

Sex-Dependent Estimation of Spinal Loads During Static Manual Material Handling Activities-Combined in vivo and in silico Analyses

Manual material handling (MMH) is considered as one of the main contributors to low back pain. While males traditionally perform MMH tasks, recently the number of females who undertake these physically-demanding activities is also increasing. To evaluate the risk of mechanical injuries, the majority of previous studies have estimated spinal forces using different modeling approaches that mostly focus on male individuals. Notable sex-dependent differences have, however, been reported in torso muscle strength and anatomy, segmental mass distribution, as well as lifting strategy during MMH. Therefore, this study aimed to use sex-specific models to estimate lumbar spinal and muscle forces during static MHH tasks in 10 healthy males and 10 females. Motion-capture, surface electromyographic from select trunk muscles, and ground reaction force data were simultaneously collected while subjects performed twelve symmetric and asymmetric static lifting (10 kg) tasks. AnyBody Modeling System was used to develop base-models (subject-specific segmental length, muscle architecture, and kinematics data) for both sexes. For females, female-specific models were also developed by taking into account for the female’s muscle physiological cross-sectional areas, segmental mass distributions, and body fat percentage. Males showed higher absolute L5-S1 compressive and shear loads as compared to both female base-models (25.3% compressive and 14% shear) and female-specific models (41% compressive and 23.6% shear). When the predicted spine loads were normalized to subjects’ body weight, however, female base-models showed larger loads (9% compressive and 16.2% shear on average), and female-specific models showed 2.4% smaller and 9.4% larger loads than males. Females showed larger forces in oblique abdominal muscles during both symmetric and asymmetric lifting tasks, while males had larger back extensor muscle forces during symmetric lifting tasks. A stronger correlation between measured and predicted muscle activities was found in females than males. Results indicate that female-specific characteristics affect the predicted spinal loads and must be considered in musculoskeletal models. Neglecting sex-specific parameters in these models could lead to the overestimation of spinal loads in females.

Knee flexion angle and muscle activations control the stability of an anterior cruciate ligament deficient joint in gait

Anterior cruciate ligament (ACL) is a primary structure and a commonly injured ligament of the knee joint. Some patients with ACL deficiency (ACLD) experience joint instability and require a reconstructive surgery to return to daily routines, some can adapt by limiting their activities while others, called copers, can return to high-level activities with no instability. We investigated the effects of alterations in the knee flexion angle (KFA) and muscle force activations on the stability and biomechanics of ACLD joints at 25, 50, and 75% periods of gait stance. ACLD joint stability is controlled by variations in both KFA and knee muscle forces. For the latter, a parameter called activity index is defined as the ratio of forces in ACL antagonists (quadriceps and gastrocnemii) to those in ACL agonists (hamstrings). Under a greater KFA (2-6° beyond the mean of reported values in healthy subjects), an ACLD joint regains its pre-injury stability levels. The ACLD joint stability also markedly improves at smaller quadriceps and larger hamstrings forces (activity indices of 2.0-3.6 at 25%) at the first half of stance and smaller gastrocnemii and larger hamstrings forces (activity indices of 0.1-1.1 at 50% and 0.1-1.2 at 75%) at the second half of stance. Activity index and KFA are both crucial when assessing the dynamic stability of an ACLD joint. These results are helpful in our understanding of the biomechanics and stability of ACLD joints towards improved prevention and treatment strategies.

A Dynamic Radiographic Imaging Study of Lumbar Intervertebral Disc Morphometry and Deformation In Vivo

Intervertebral discs are important structural components of the spine but also are significant sources of morbidity, especially for the “low back” lumbar region. Mechanical damage to, or degeneration of, the lumbar discs can diminish their structural integrity and elicit debilitating low back pain. Advancement of reparative or regenerative means to treat damaged or degenerated discs is hindered by a lack of basic understanding of the disc load-deformation characteristics in vivo. The current study presents an in vivo analysis of the morphometry and deformation of lumbar (L2-S1) intervertebral discs in 10 healthy participants while performing a common lifting act, using novel dynamic radiographic imaging of the lumbar vertebral body motion. Data analyses show uniquely different (p < 0.05) characteristics in morphometry, normal and shear strain patterns of the L5S1 discs, while the rest of lumbar discs exhibit great similarity. In particular shear strains in L2-L5 discs exhibited stronger linear correlations (R² ≥ 0.80) between strain changes and amount of lumbar flexion-extension motion compared to L5S1 (R² ≤ 0.5). The study therefore advances the state of knowledge on in vivo mechanical responses of the lumbar intervertebral discs during functional tasks.

Squat Lifting Imposes Higher Peak Joint and Muscle Loading Compared to Stoop Lifting

(1) Background: Yearly, more than 40% of the European employees suffer from work-related musculoskeletal disorders. Still, ergonomic guidelines defining optimal lifting techniques to decrease work-related musculoskeletal disorders (WMSDs) has not been unambiguously defined. Therefore, this study investigates if recommended squat lifting imposes lower musculoskeletal loading than stoop lifting while using a complex full body musculoskeletal OpenSim model. (2) Methods: Ten healthy participants lifted two different weights using both lifting techniques. 3D marker trajectories and ground reaction forces were used as input to calculate joint angles, moments and power using a full body musculoskeletal model with articulated lumbar spine. In addition, the muscle activity of nine different muscles was measured to investigate muscle effort when lifting. (3) Results: Peak moments and peak joint power in L5S1 were not different between the squat and the stoop, but higher peak moments and peak power in the hip, knee, elbow and shoulder were found during squat lifting. Moment impulses in L5S1 were higher during stoop lifting. This is reflected in higher peak electromyography (EMG) but lower muscle effort in prior described muscles during the squat. (4) Conclusions: Squat lifting imposes higher peak full body musculoskeletal loading but similar low back loading compared to stoop lifting, as reflected in peak moments, peak power, and peak EMG.

Predicting the influence of hip and lumbar flexibility on lifting motions using optimal control

Computational models of the human body coupled with optimization can be used to predict the influence of variables that cannot be experimentally manipulated. Here, we present a study that predicts the motion of the human body while lifting a box, as a function of flexibility of the hip and lumbar joints in the sagittal plane. We modeled the human body in the sagittal plane with joints actuated by pairs of agonist-antagonist muscle torque generators, and a passive hamstring muscle. The characteristics of a stiff, average and flexible person were represented by co-varying the lumbar range-of-motion, lumbar passive extensor-torque and the hamstring passive muscle-force. We used optimal control to solve for motions that simulated lifting a 10 kg box from a 0.3 m height. The solution minimized the total sum of the normalized squared active and passive muscle torques and the normalized passive hamstring muscle forces, over the duration of the motion. The predicted motion of the average lifter agreed well with experimental data in the literature. The change in model flexibility affected the predicted joint angles, with the stiffer models flexing more at the hip and knee, and less at the lumbar joint, to complete the lift. Stiffer models produced similar passive lumbar torque and higher hamstring muscle force components than the more flexible models. The variation between the motion characteristics of the models suggest that flexibility may play an important role in determining lifting technique.

A model-based approach for estimation of changes in lumbar segmental kinematics associated with alterations in trunk muscle forces

The kinematics information from imaging, if combined with optimization-based biomechanical models, may provide a unique platform for personalized assessment of trunk muscle forces (TMFs). Such a method, however, is feasible only if differences in lumbar spine kinematics due to differences in TMFs can be captured by the current imaging techniques. A finite element model of the spine within an optimization procedure was used to estimate segmental kinematics of lumbar spine associated with five different sets of TMFs. Each set of TMFs was associated with a hypothetical trunk neuromuscular strategy that optimized one aspect of lower back biomechanics. For each set of TMFs, the segmental kinematics of lumbar spine was estimated for a single static trunk flexed posture involving, respectively, 40° and 10° of thoracic and pelvic rotations. Minimum changes in the angular and translational deformations of a motion segment with alterations in TMFs ranged from 0° to 0.7° and 0 mm to 0.04 mm, respectively. Maximum changes in the angular and translational deformations of a motion segment with alterations in TMFs ranged from 2.4° to 7.6° and 0.11 mm to 0.39 mm, respectively. The differences in kinematics of lumbar segments between each combination of two sets of TMFs in 97% of cases for angular deformation and 55% of cases for translational deformation were within the reported accuracy of current imaging techniques. Therefore, it might be possible to use image-based kinematics of lumbar segments along with computational modeling for personalized assessment of TMFs.

Estimation of lumbar spinal loading and trunk muscle forces during asymmetric lifting tasks: application of whole-body musculoskeletal modelling in OpenSim

Large spinal compressive force combined with axial torsional shear force during asymmetric lifting tasks is highly associated with lower back injury (LBI). The aim of this study was to estimate lumbar spinal loading and muscle forces during symmetric lifting (SL) and asymmetric lifting (AL) tasks using a whole-body musculoskeletal modelling approach. Thirteen healthy males lifted loads of 7 and 12 kg under two lifting conditions (SL and AL). Kinematic data and ground reaction force data were collected and then processed by a whole-body musculoskeletal model. The results show AL produced a significantly higher peak lateral shear force as well as greater peak force of psoas major, quadratus lumborum, multifidus, iliocostalis lumborum pars lumborum, longissimus thoracis pars lumborum and external oblique than SL. The greater lateral shear forces combined with higher muscle force and asymmetrical muscle contractions may have the biomechanical mechanism responsible for the increased risk of LBI during AL. Practitioner Summary: Estimating lumbar spinal loading and muscle forces during free-dynamic asymmetric lifting tasks with a whole-body musculoskeletal modelling in OpenSim is the core value of this research. The results show that certain muscle groups are fundamentally responsible for asymmetric movement, thereby producing high lumbar spinal loading and muscle forces, which may increase risks of LBI during asymmetric lifting tasks.

MUSCLE FORCE MAGNITUDES IN THE HUMAN LEG FOR ISOMETRIC EXERCISES WITH VARIOUS RESULTANT FORCE DIRECTIONS AND JOINT ANGLES

Using a gradient-based numerical optimization routine, the force magnitudes required of 10 major sagittal plane leg muscles to create a constant magnitude isometric resultant force against a fixed surface at the toe directed anteriorly, posteriorly, superiorly, and inferiorly were quantitatively predicted for three sets of joint angles: a straight leg configuration, with the knee flexed, and with both the hip and knee flexed. Comparisons over the conditions studied for each individual system muscle found that the maximum variation occurred in knee and hip extensor forces (up to two orders of magnitude). Comparisons within the set of active muscles for each studied condition identified dominant muscles and muscle functions. All anteriorly-directed and posteriorly-directed resultants required a small number of muscles with common functions (mainly knee or hip extensors) and large force magnitudes (O(1000[Formula: see text]N)). In contrast, a large number of muscles, with wide-ranging synergistic and antagonistic functions, acting across multiple joints with relatively small magnitudes (O(100[Formula: see text]N)) were needed to create the superiorly-directed resultant with flexed hip and knee. With good correlation to experimentally measured trends in the interrelationships between leg joint angles and isometric forces, the systematic muscle force prediction and analysis presented in this work can be used to guide the design of targeted muscle strengthening exercises and study of muscle-specific injury.

Persons with unilateral transfemoral amputation experience larger spinal loads during level-ground walking compared to able-bodied individuals

BACKGROUND

Persons with lower limb amputation walk with increased and asymmetric trunk motion; a characteristic that is likely to impose distinct demands on trunk muscles to maintain equilibrium and stability of the spine. However, trunk muscle responses to such changes in net mechanical demands, and the resultant effects on spinal loads, have yet to be determined in this population.

METHODS

Building on a prior study, trunk and pelvic kinematics collected during level-ground walking from 40 males (20 with unilateral transfemoral amputation and 20 matched controls) were used as inputs to a kinematics-driven, nonlinear finite element model of the lower back to estimate forces in 10 global (attached to thorax) and 46 local (attached to lumbar vertebrae) trunk muscles, as well as compression, lateral, and antero-posterior shear forces at all spinal levels.

FINDINGS

Trunk muscle force and spinal load maxima corresponded with heel strike and toe off events, and among persons with amputation, were respectively 10-40% and 17-95% larger during intact vs. prosthetic stance, as well as 6-80% and 26-60% larger during intact stance relative to controls.

INTERPRETATION

During gait, larger spinal loads with transfemoral amputation appear to be the result of a complex pattern of trunk muscle recruitment, particularly involving co-activation of antagonistic muscles during intact limb stance; a period when these individuals are confident and likely to use the trunk to assist with forward progression. Given the repetitive nature of walking, repeated exposure to such elevated loading likely increases the risk for low back pain in this population.

The effect of parameters of equilibrium-based 3-D biomechanical models on extracted muscle synergies during isometric lumbar exertion

A hallmark of more advanced models is their higher details of trunk muscles represented by a larger number of muscles. The question is if in reality we control these muscles individually as independent agents or we control groups of them called "synergy". To address this, we employed a 3-D biomechanical model of the spine with 18 trunk muscles that satisfied equilibrium conditions at L4/5, with different cost functions. The solutions of several 2-D and 3-D tasks were arranged in a data matrix and the synergies were computed by using non-negative matrix factorization (NMF) algorithms. Variance accounted for (VAF) was used to evaluate the number of synergies that emerged by the analysis, which were used to reconstruct the original muscle activations. It was showed that four and six muscle synergies were adequate to reconstruct the input data of 2-D and 3-D torque space analysis. The synergies were different by choosing alternative cost functions as expected. The constraints affected the extracted muscle synergies, particularly muscles that participated in more than one functional tasks were influenced substantially. The compositions of extracted muscle synergies were in agreement with experimental studies on healthy participants. The following computational methods show that the synergies can reduce the complexity of load distributions and allow reduced dimensional space to be used in clinical settings.

Variations in trunk muscle activities and spinal loads following posterior lumbar surgery: A combined in vivo and modeling investigation

BACKGROUND

Iatrogenic injuries to paraspinal muscles during posterior lumbar surgery cause a reduction in their contractile cross-sectional area and thus presumably their postoperative activation. This study investigates the effect of such intraoperative injuries on postoperative patterns of muscle activations and spinal loads during various activities using a combined modeling and in vivo MR imaging approach.

METHODS

A three-dimensional, multi-joint, musculoskeletal model was used to estimate pre- and postoperative muscle forces and spinal loads under various activities in upright and flexed postures. According to our in vivo pre- and postoperative (~6 months) measurements in six patients using a 3-Tesla-MR scanner, physiological cross-sectional areas of multifidus and erector spinae were reduced in the postoperative model by 26 and 11%, respectively.

FINDINGS

Postoperative trunk extension strength was predicted to decrease by ~23% from 215 Nm in the intact model to 165 Nm in the postoperative model. Postoperative force in multifidus fascicles decreased by ~21-40% in flexion tasks and by ~14-35% in upright tasks. In contrast, the sum of the forces in all other intact and less injured extensor muscles slightly increased (by <6%) in the postoperative model. Postoperative L5-S1 compressive and shear loads varied slightly (by ~3%).

INTERPRETATION

Intraoperative injuries induced a shift in load-sharing from the most injured muscle (multifidus) toward other less injured and intact muscles during all simulated activities. Postoperative rehabilitation programs should therefore strengthen and facilitate (while avoiding muscle imbalance) not only the injured multifidus but also other intact and less injured trunk muscles that play a compensatory role after the operation.

An optimization-based method for prediction of lumbar spine segmental kinematics from the measurements of thorax and pelvic kinematics

Given measurement difficulties, earlier modeling studies have often used some constant ratios to predict lumbar segmental kinematics from measurements of total lumbar kinematics. Recent imaging studies suggested distribution of lumbar kinematics across its vertebrae changes with trunk rotation, lumbar posture, and presence of load. An optimization-based method is presented and validated in this study to predict segmental kinematics from measured total lumbar kinematics. Specifically, a kinematics-driven biomechanical model of the spine is used in a heuristic optimization procedure to obtain a set of segmental kinematics that, when prescribed to the model, were associated with the minimum value for the sum of squared predicted muscle stresses across all the lower back muscles. Furthermore, spinal loads estimated using the predicted kinematics by the present method were compared with those estimated using constant ratios. Predicted segmental kinematics were in good agreement with those obtained by imaging with an average error of ~10%. Compared with those obtained using constant ratios, predicted spinal loads using segmental kinematics obtained here were in general smaller. In conclusion, the proposed method offers an alternative tool for improving model-based estimates of spinal loads where image-based measurement of lumbar kinematics is not feasible.

Development and implementation of a coupled computational muscle force optimization bone shape adaptation modeling method

Improved methods to analyze and compare the muscle-based influences that drive bone strength adaptation can aid in the understanding of the wide array of experimental observations about the effectiveness of various mechanical countermeasures to losses in bone strength that result from age, disuse, and reduced gravity environments. The coupling of gradient-based and gradientless numerical optimization routines with finite element methods in this work results in a modeling technique that determines the individual magnitudes of the muscle forces acting in a multisegment musculoskeletal system and predicts the improvement in the stress state uniformity and, therefore, strength, of a targeted bone through simulated local cortical material accretion and resorption. With a performance-based stopping criteria, no experimentally based or system-based parameters, and designed to include the direct and indirect effects of muscles attached to the targeted bone as well as to its neighbors, shape and strength alterations resulting from a wide range of boundary conditions can be consistently quantified. As demonstrated in a representative parametric study, the developed technique effectively provides a clearer foundation for the study of the relationships between muscle forces and the induced changes in bone strength. Its use can lead to the better control of such adaptive phenomena.

Partitioning of knee joint internal forces in gait is dictated by the knee adduction angle and not by the knee adduction moment

Medial knee osteoarthritis is a debilitating disease. Surgical and conservative interventions are performed to manage its progression via reduction of load on the medial compartment or equivalently its surrogate measure, the external adduction moment. However, some studies have questioned a correlation between the medial load and adduction moment. Using a musculoskeletal model of the lower extremity driven by kinematics-kinetics of asymptomatic subjects at gait midstance, we aim here to quantify the relative effects of changes in the knee adduction angle versus changes in the adduction moment on the joint response and medial/lateral load partitioning. The reference adduction rotation of 1.6° is altered by ±1.5° to 3.1° and 0.1° or the knee reference adduction moment of 17Nm is varied by ±50% to 25.5Nm and 8.5Nm. Quadriceps, hamstrings and tibiofemoral contact forces substantially increased as adduction angle dropped and diminished as it increased. The medial/lateral ratio of contact forces slightly altered by changes in the adduction moment but a larger adduction rotation hugely increased this ratio from 8.8 to a 90 while in contrast a smaller adduction rotation yielded a more uniform distribution. If the aim in an intervention is to diminish the medial contact force and medial/lateral load ratio, a drop of 1.5° in adduction angle is much more effective (causing respectively 12% and 80% decreases) than a reduction of 50% in the adduction moment (causing respectively 4% and 13% decreases). Substantial role of changes in adduction angle is due to the associated alterations in joint nonlinear passive resistance. These findings explain the poor correlation between knee adduction moment and tibiofemoral compartment loading during gait suggesting that the internal load partitioning is dictated by the joint adduction angle.

Evaluation of knee joint muscle forces and tissue stresses-strains during gait in severe OA versus normal subjects

Osteoarthritis (OA) is the leading cause of pain and disability in the elderly with the knee being the most affected weight bearing joint. We used a musculoskeletal biomechanical model of the lower extremity including a detailed validated knee joint finite element model to compute lower extremity muscle forces and knee joint stresses-strains during the stance phase of gait. The model was driven by gait data on OA patients, and results were compared with those of the same model driven by data on normal controls. Additional analyses were performed with altered cartilage-menisci properties to evaluate the effects of deterioration during OA. In OA patients compared to normal subjects, muscle forces dropped at nearly all stance periods except mid-stance. Force in the anterior cruciate ligament remained overall the same. Total contact forces-stresses deceased by about 25%. Alterations in properties due to OA had negligible effects on muscle forces, but increased contact areas and cartilage strains and reduced contact pressures. Reductions in contact stresses and increases in tissue strains and transfer of load via menisci are partly due to the altered kinetics-kinematics of gait and partly due to deterioration in cartilage-menisci properties in OA patients.

Predictive equations for lumbar spine loads in load-dependent asymmetric one- and two-handed lifting activities

BACKGROUND

Asymmetric lifting activities are associated with low back pain.

METHODS

A finite element biomechanical model is used to estimate spinal loads during one- and two-handed asymmetric static lifting activities. Model input variables are thorax flexion angle, load magnitude as well as load sagittal and lateral positions while response variables are L4-L5 and L5-S1 disc compression and shear forces. A number of levels are considered for each input variable and all their possible combinations are introduced into the model. Robust yet user-friendly predictive equations that relate model responses to its inputs are established.

FINDINGS

Predictive equations with adequate goodness-of-fit (R(2) ranged from ~94% to 99%, P≤0.001) that relate spinal loads to task (input) variables are established. Contour plots are used to identify combinations of task variable levels that yield spine loads beyond the recommended limits. The effect of uncertainties in the measurements of asymmetry-related inputs on spinal loads is studied.

INTERPRETATION

A number of issues regarding the NIOSH asymmetry multiplier are discussed and it is concluded that this multiplier should depend on the trunk posture and be defined in terms of the load vertical and horizontal positions. Due to an imprecise adjustment of the handled load magnitude this multiplier inadequately controls the biomechanical loading of the spine. Ergonomists and bioengineers, faced with the dilemma of using either complex but more accurate models on one hand or less accurate but simple models on the other hand, have hereby easy-to-use predictive equations that quantify spinal loads under various occupational tasks.