Spinal Loading Research Articles

Associations of back muscle endurance with occupational back muscle activity and spinal loading among subsistence farmers and office workers in Rwanda.

Over the course of the physical activity transition, machines have largely replaced skeletal muscle as the source of work for locomotion and other forms of occupational physical activity in industrial environments. To better characterize this transition and its effect on back muscles and the spine, we tested to what extent typical occupational activities of rural subsistence farmers demand higher magnitudes and increased variability of back muscle activity and spinal loading compared to occupational activities of urban office workers in Rwanda, and whether these differences were associated with back muscle endurance, the dominant risk factor for back pain. Using electromyography, inertial measurement units, and OpenSim musculoskeletal modeling, we measured back muscle activity and spinal loading continuously while participants performed occupational activities for one hour. We measured back muscle endurance using electromyography median frequency analysis. During occupational work, subsistence farmers activate their back muscles and load their spines at 390% higher magnitudes and with 193% greater variability than office workers. Partial correlations accounting for body mass show magnitude and variability response variables are positively associated with back muscle endurance (R = 0.39-0.90 [P < 0.001-0.210] and R = 0.54-0.72 [P = 0.007-0.071], respectively). Body mass is negatively correlated with back muscle endurance (R = -0.60, P = 0.031), suggesting higher back muscle endurance may be also partly attributable to having lower body mass. Because higher back muscle endurance is a major factor that prevents back pain, these results reinforce evidence that under-activating back muscles and under-loading spines at work increases vulnerability to back pain and may be an evolutionary mismatch. As sedentary occupations become more common, there is a need to study the extent to which occupational and leisure time physical activities that increase back muscle endurance helps prevent back pain.

Effects of a soft back exoskeleton on lower lumbar spine loads during manual materials handling: a musculoskeletal modelling study

The aim of this study was to append a passive soft back exoskeleton to a validated musculoskeletal model and assess its effectiveness in reducing lumbar loads. Fifteen participants lifted a box, with and without wearing a CORFOR® exoskeleton. A full body OpenSim model was used to estimate lumbar joint moments and reaction forces, as well as low back muscles forces. Wearing the exoskeleton reduced the peak flexion moment, muscles forces, as well as peak compressive and shear forces. This musculoskeletal modelling study shows that wearing the exoskeleton may reduce lumbar spine loads and may contribute to prevent low back disorders.

Assessment of Musculoskeletal Loads among Office Workers Due to Predicted BMI Increase

The increasing prevalence of adult obesity is a major global health problem. Causes include sedentary lifestyles and physical inactivity, as well as the place and style of work. This work aimed to assess the impact of an increase in BMI, resulting from a sedentary lifestyle, on the functioning of the human musculoskeletal system during daily activities. Standing and sitting down were analyzed using numerical simulations in the AnyBody Modeling System. The results showed that by 2050, an increase in body weight will lead to an 11% increase in lumbar spine loads and up to 14% increase in knee and hip joint loads compared to 2015. Increased body weight will also increase muscle fatigue, requiring the body to exert more force to perform the same activities. The BMI change prediction model shows that the average man will be overweight or obese in the coming years, which will likely increase musculoskeletal loads.

Correlation between sagittal balance and thoracolumbar elastic energy parameters in 42 spines subject to spondylolisthesis or spinal stenosis and 21 normal spines

The curvature of the lumbar spine plays a critical role in maintaining spinal function, stability, weight distribution, and load transfer. We have developed a mathematical model of the lumbar spine curve by introducing a novel mechanism: minimization of the elastic bending energy of the spine with respect to two biomechanical parameters: dimensionless lumbosacral spinal curvature cLS and dimensionless curvature increment along the spine CI. While most of the biomechanical studies focus on a particular segment of the spine, the distinction of the presented model is that it describes the shape of the thoracolumbar spine by considering it as a whole (non-locally) and thus includes interactions between the different spinal levels in a holistic approach. From radiographs, we have assessed standard geometrical parameters: lumbar lordosis LL, pelvic incidence PI, pelvic tilt PT, sacral slope ψ0 and sagittal balance parameter SB = sagittal vertical axis (SVA)/sacrum-bicoxofemoral distance (SFD) of 42 patients with lumbar spinal stenosis (SS) or degenerative spondylolisthesis (SL) and 21 radiologically normal subjects. SB statistically significantly correlated with model parameters cL5 (r = −0.34, p = 0.009) and −CI (r = 0.33, p = 0.012) but not with standard geometrical parameters. A statistically significant difference with sufficient statistical power between the patients and the normal groups was obtained for cLS, CI, and SB but not for standard geometrical parameters. The model provides a possibility to predict changes in the thoracolumbar spine shape in surgery planning and in assessment of different spine pathologies.

The effect of a soft active back support exosuit on trunk motion and thoracolumbar spine loading during squat and stoop lifts

Back support exosuits aim to reduce tissue demands and thereby risk of injury and pain. However, biomechanical analyses of soft active exosuit designs have been limited. The objective of this study was to evaluate the effect of a soft active back support exosuit on trunk motion and thoracolumbar spine loading in participants performing stoop and squat lifts of 6 and 10 kg crates, using participant-specific musculoskeletal models. The exosuit did not change overall trunk motion but affected lumbo-pelvic motion slightly, and reduced peak compressive and shear vertebral loads at some levels, although shear increased slightly at others. This study indicates that soft active exosuits have limited kinematic effects during lifting, and can reduce spinal loading depending on the vertebral level. These results support the hypothesis that a soft exosuit can assist without limiting trunk movement or negatively impacting skeletal loading and have implications for future design and ergonomic intervention efforts.

The effects of spinal stabilization exercises on functional exercise capacity in individuals with transtibial amputation: A randomized trial.

Optimization of mobility activities with spinal stabilization exercises and reducing spinal loads in individuals with transtibial amputation are recommended by researchers; however, the effect of spinal stabilization exercises on functional exercise capacity has not been investigated. This study aims to investigate the effect of spinal stabilisation exercises on functional exercise capacity in individuals with transtibial amputation. Randomized controlled trial with concealed allocation, assessor blinding. Eighteen individuals with transtibial amputation were included. Conventional physiotherapy program was applied to individuals in control group, and spinal stabilization exercise training was applied to individuals in intervention group together with conventional physiotherapy program for 8 weeks. The primary outcome was maximal oxygen uptake during 6-Minute Stepper Test. The secondary outcomes were maximal energy expenditure, step count and fatigue during 6-Minute Stepper Test, perceived mobility level with prosthesis, and stabilization strength of deep spinal muscles. In the within-group comparisons, the change in the functional exercise capacity, step count, and stabilization strength of deep spinal muscles were found to be statistically significant in both groups. In the between-group comparison, the changes in the functional exercise capacity, perceived mobility level, and stabilization strength of deep spinal muscles were found to be significant in favor of the group in which spinal stabilization exercises were applied. Along with improvement in proximal control of movement, improvements were observed in some parameters that are indicative of physical capacity. Considering the described developments, exercise programs focused on spinal stabilization in individuals with amputation are becoming a serious alternative in clinical practice.

Facet deflection and strain are dependent on axial compression and distraction in C5-C7 spinal segments under constrained flexion.

Facet fractures are frequently associated with clinically observed cervical facet dislocations (CFDs); however, to date there has only been one experimental study, using functional spinal units (FSUs), which has systematically produced CFD with concomitant facet fracture. The role of axial compression and distraction on the mechanical response of the cervical facets under intervertebral motions associated with CFD in FSUs has previously been shown. The same has not been demonstrated in multi-segment lower cervical spine specimens under flexion loading (postulated to be the local injury vector associated with CFD). This study investigated the mechanical response of the bilateral inferior C6 facets of thirteen C5-C7 specimens (67±13 yr, 6 male) during non-destructive constrained flexion, superimposed with each of five axial conditions: (1) 50 N compression (simulating weight of the head); (2-4) 300, 500, and 1000 N compression (simulating the spectrum of intervertebral compression resulting from neck muscle bracing prior to head-first impact and/or externally applied compressive forces); and, (5) 2 mm of C6/C7 distraction (simulating the intervertebral distraction present during inertial loading of the cervical spine by the weight of the head). Linear mixed-effects models (α = 0.05) assessed the effect of axial condition. Increasing amounts of intervertebral compression superimposed on flexion rotations, resulted in increased facet surface strains (range of estimated mean difference relative to Neutral: maximum principal = 77 to 110 με, minimum principal = 126 to 293 με, maximum shear = 203 to 375 με) and angular deflection of the bilateral inferior C6 facets relative to the C6 vertebral body (range of estimated mean difference relative to Neutral = 0.59° to 1.47°). These findings suggest increased facet engagement and higher load transfer through the facet joint, and potentially a higher likelihood of facet fracture under the compressed axial conditions.

The effects of setup parameters on the measured kinetic output of cervical disc prostheses

Mechanical testing machines are used to evaluate kinematics, kinetics, wear, and efficacy of spinal implants. The simulation of "physiological" spinal loading conditions necessitates the simultaneous use of multiple actuators. The challenge in achieving a desired loading profile lies in achieving close synchronization of these actuators. Errors in load application can be attributed to both the control system and the intrinsic sample response. Moreover, the presence of friction in the setup can have an impact on the measured outcome. The optimization of setup parameters can substantially improve the ability to simulate spinal loading conditions and obtain reliable data on implant performance. In this study, a reproducible kinematic test protocol was developed to evaluate the sensitivity of the kinetic response (i.e., measured loads, moments, and stiffnesses) of a cervical disc prosthesis to several testing parameters. In this context, five ceramic ball and socket sample implants were mounted in a 6 DOF material testing machine and tested with a constant axial compressive force of 100 N in two motion modes: 1) flexion-extension (±7.5°) and 2) lateral bending (±6°). Parameters including rotation rate, slider friction, friction between the samples' articulating surfaces, and moment arm were considered to determine their effects on measured kinetic parameters. The sensitivity analysis indicated that all setup parameters except friction between the samples' articulating surfaces had a substantial effect on the results. The findings were then compared to predictions from a free body diagram to determine the optimal setup parameters. Consequently, the setup with the lowest rotation rate and employing passive sliders yielded results that were consistent with the free body diagram. This study demonstrated the significance of a comprehensive setup evaluation for reliable and reproducible testing of spinal implants, also for comparison between labs.

Cadence matters: Influence of cadence on spinal load during running

BackgroundRunning exposes the body to physiological and mechanical stresses that generate musculoskeletal injuries, such as low back pain due to large spinal loading. Increasing running cadence may reduce impact forces and spinal shrinkage. Research questionThis study aimed to determine the relationship between spinal loading and running cadence. MethodsThis cross-sectional study included 15 runners from the local community (36 ± 11 years; 23 ± 2 kg.m−2, and 8 ± 9 years of running experience) who ran for 30 min (R30) and 60 min (R60) at a constant speed (10 km.h−1). The spinal loading was assessed via fine stature variation measurements before the run (baseline) at R30 and R60. Cadence was monitored via a wristwatch. The cadence ranged from 150 to 180 steps.min−1. A t-test was used to compare stature loss between R30 and R60 (relative to baseline), and a stepwise linear regression equation was used to identify the relationship between cadence and stature variation in each instant. ResultsThere was a stature loss throughout the race (R30 = 5.27 ± 1.92 mm and R60 =7.51 ± 2.51 mm). A linear regression analysis revealed a negative relationship between stature loss and cadence, indicating that running at a faster cadence produces smaller spinal loading than running at slower cadences after R60 (R2 = 0.38; p<0.05). SignificanceIncreasing running cadence might cause less spinal loading than running with a slower cadence, which may reduce the risk of injury and back disorders in runners.

PUBLIC AND PATIENT ENGAGEMENT IN RESEARCH TO DEVELOP TECHNOLOGIES FOR PREDICTING SPINE LOADS IN VIVO

BackgroundOur current research aims to develop technologies to predict spinal loads in vivo using a combination of imaging and modelling methods. To ensure the project's success and inform future applications of the technology, we sought to understand the opinions and perspectives of patients and the public.MethodsA 90-minute public and patient involvement event was developed in collaboration with Exeter Science Centre and held on World Spine Day 2023. The event involved a brief introduction to the project goals followed by an interactive questionnaire to gauge the participants’ background knowledge and interest. The participants then discussed five topics: communication, future directions of the research, concerns about the research protocol, concerns about data, and interest in the project team and research process. A final questionnaire was used to determine their thoughts about the event.ResultsTwelve adults attended the event, many motivated by their experience or interest in back pain. A thematic analysis was used to review participant comments on the research project, identifying the need to relate the research to everyday life, present risks in various ways, and be transparent about funding and data sharing. In terms of future applications, participants felt the technology should be used to understand normal spine behaviour, prevent problems, and improve treatment. Participants agreed that they had got something positive out of engaging in the event.ConclusionEngagement with public and patient stakeholders is an essential activity that can generate vital information to inform and add value to technology development projects.Conflicts of interestNo conflicts of interestSources of fundingEPSRC grants EP/V036602/1 (Meakin, Holsgrove &amp; Javadi) and EP/V032275/1 (Holt &amp; Williams).

Cyclic Mechanism Affects Lumbar Spine Creep Response.

This study aims to explore how cyclic loading influences creep response in the lumbar spine under combined flexion-compression loading. Ten porcine functional spinal units (FSUs) were mechanically tested in cyclic or static combined flexion-compression loading. Creep response between loading regimes was compared using strain-time histories and linear regression. High-resolution computed tomography (µCT) visualized damage to FSUs. Statistical methods, ANCOVA and ANOVA, assessed differences in behavior between loading regimes. Cyclic and static loading regimes exhibited distinct creep response patterns and biphasic response. ANCOVA and ANOVA analyses revealed significant differences in slopes of creep behavior in both linear phases. Cyclic tests consistently showed endplate fractures in µCT imaging. The study reveals statistically significant differences in creep response between cyclic and static loading regimes in porcine lumbar spinal units under combined flexion-compression loading. The observed biphasic behavior suggests distinct phases of tissue response, indicating potential shifts in load transfer mechanisms. Endplate fractures in cyclic tests suggest increased injury risk compared to static loading. These findings underscore the importance of considering loading conditions in computational models and designing preventive measures for occupations involving repetitive spinal loading.

Modelling the female torso and breast during physical activity: Implications on spinal loading

ABSTRACT Methods of modelling the female torso during physical activity often neglect the position and movement of the breast. This novel investigation compares three female torso modelling approaches that differ in complexity (integrated breast, fixed breast, dynamic breast) to determine the effect on spinal joint moments during running and jumping. The commonly used integrated breast model distributed breast mass within the torso, the fixed breast model attached the mass of the breasts to fixed positions on the anterior of the torso, and a new dynamic breast model enabled relative motion between the breasts and anterior torso. Key findings demonstrated minimal differences in lumbar spine moments (<0.05 Nm/kg; 4%) between integrated breast and fixed breast models but greater differences, up to 0.86 Nm/kg (68%) during running and 0.89 Nm/kg (82%) during jumping, when breast motion was included. Thoracic spine moments revealed similar patterns with minimal differences (<0.05 Nm/kg; 11%) between integrated breast and fixed breast models and greater differences, up to 0.48 Nm/kg (92%) during running and 0.63 Nm/kg (66%) during jumping, with the dynamic breast model. Future female musculoskeletal models should consider including breast mass and motion to avoid mis-representing spinal loading in females during running and jumping.

Surgical treatment of L5 spondylolysis in an athlete using custom-made implant

Background: Spondylolysis is one of the most common causes of lower back pain in children and adolescents who are professionally involved in sports. It is noted that spondylolysis is observed more often when practicing a number of sports that are associated with repeated axial load and/or hyperextension of the lumbar spine with rotation. In most cases, the treatment of spondylolysis, including cases of its occurrence in professional athletes, is conservative. Surgical treatment is indicated only if conservative treatment is ineffective or if symptoms progress. One of the most common methods of surgical treatment of spondylolysis is to restore the integrity of the arch using various metal structures. The use of additive methods for the manufacture of individual implants currently allows the manufacture of personalized implants with a number of advantages. The article describes the first experience of using an individual implant for surgical treatment of spondylolysis and provides a brief review of the literature. Clinical case description: A clinical case is presented involving the treatment of a 16-year-old female patient who is a professional gymnast. The report includes a description of the patient’s medical history, clinical manifestations, and specialized diagnostic methods. The preoperative planning, design of a custom implant, the surgical procedure, and long-term treatment outcomes are detailed. A brief literature review highlights the results of conservative treatment, the main indications and methods of surgical therapy for spondylolysis, and justifies the use of a custom-made implant for its surgical treatment. CONCLUSION: For the surgical treatment of L5 spondylolysis and restoration of vertebral arch integrity without limiting motion on L5-S1 level, the use of a custom-made implant is possible. The use of customized implants may improve outcomes in cases where spondylolysis is combined with abnormalities and individual characteristics of the vertebral bone structures, including the patient’s sports activity.

Effects of laptop screen height on neck and shoulder muscle fatigue and spine loading for office workers.

Due to the unfavourable neck-shoulder muscle loads caused by poor posture, the people who use the laptop for a long time may face the risk of neck and shoulder injuries. The purpose of this study investigates the impact of the screen height on the muscle activation of head flexion, neck and shoulder, and the cervical spine torque to provide the favorite screen height for laptop user. Twelve healthy young participants completed a15-minute task of the reading at the four different screen heights. sEMG signals of the splenius capitis (SC) and upper trapezius (UT) were measured and calculated the root mean square (RMS) and mean power frequency (MPF) to determine muscle fatigue. The different height of laptop users was simulated and the forces on the spine of users at different screen heights were analyzed by Jack. Adjusting the height of the laptop screen can effectively reduce head flexion and muscle activity of SC and UT, and has a positive effect on reducing fatigue of SC, but has no significant effect on UT. Adjusting the height of the laptop screen can delay the occurrence of SC muscle fatigue to a certain extent. The joint analysis of sEMG spectrum and amplitude reports that the screen heights of D15 and D45 have the highest and the lowest frequency of fatigue, respectively. At the same time, the moment of spineT1/T2 and spineL4/L5 decrease with the increase of screen height.

Spinal loads during dynamic full flexion and return to standing posture in different age and sex groups: A musculoskeletal model study

During forward flexion, spine motion varies due to age and sex differences. Previous studies showed that lumbar/pelvis range of flexion (RoF) and lumbo-pelvic ratio (L/P) are age/sex dependent. How variation of these parameters affects lumbar loading in a normal population requires further assessment. We aimed to estimate lumbar loads during dynamic flexion-return cycle and the differences in peak loads (compression) and corresponding trunk inclinations due to variation in lumbar/pelvis RoF and L/P.Based on in vivo L/P (0.11–3.44), temporal phases of flexion (early, middle, and later), the lumbar (45-55°) and hip (60-79°) RoF; full flexion-return cycles of six seconds were reconstructed for three age groups (20–35, 36–50 and 50+ yrs.) in both sexes. Six inverse dynamic analyses were performed with a 50th percentile model, and differences in peak loads and corresponding trunk inclinations were calculated.Peak loads at L4-L5 were 179 N higher in younger males versus females, but 228 N and 210 N lower in middle-aged and older males, respectively, compared to females. Females exhibited higher trunk inclinations (6°-20°) than males across all age groups. Age related differences in L4-L5 peak loads and corresponding trunk inclinations were found up to 415 N and 19° in males and 152 N and 13° in females. With aging, peak loads were reduced in males but were found non-monotonic in females, whereas trunk inclinations at peak loads were reduced in both sexes from young to middle/old age groups.In conclusion, lumbar loading and corresponding trunk inclinations varied notably due to age/sex differences. Such data may help distinguishing normal or pathological condition of the lumbar spine.

Musculoskeletal spine modeling in large patient cohorts: how morphological individualization affects lumbar load estimation.

Introduction: Achieving an adequate level of detail is a crucial part of any modeling process. Thus, oversimplification of complex systems can lead to overestimation, underestimation, and general bias of effects, while elaborate models run the risk of losing validity due to the uncontrolled interaction of multiple influencing factors and error propagation. Methods: We used a validated pipeline for the automated generation of multi-body models of the trunk to create 279 models based on CT data from 93 patients to investigate how different degrees of individualization affect the observed effects of different morphological characteristics on lumbar loads. Specifically, individual parameters related to spinal morphology (thoracic kyphosis (TK), lumbar lordosis (LL), and torso height (TH)), as well as torso weight (TW) and distribution, were fully or partly considered in the respective models according to their degree of individualization, and the effect strengths of these parameters on spinal loading were compared between semi- and highly individualized models. T-distributed stochastic neighbor embedding (T-SNE) analysis was performed for overarching pattern recognition and multiple regression analyses to evaluate changes in occurring effects and significance. Results: We were able to identify significant effects (p < 0.05) of various morphological parameters on lumbar loads in models with different degrees of individualization. Torso weight and lumbar lordosis showed the strongest effects on compression (β ≈ 0.9) and anterior-posterior shear forces (β ≈ 0.7), respectively. We could further show that the effect strength of individual parameters tended to decrease if more individual characteristics were included in the models. Discussion: The induced variability due to model individualization could only partly be explained by simple morphological parameters. Our study shows that model simplification can lead to an emphasis on individual effects, which needs to be critically assessed with regard to in vivo complexity. At the same time, we demonstrated that individualized models representing a population-based cohort are still able to identify relevant influences on spinal loading while considering a variety of influencing factors and their interactions.

Spinal joint moment prediction following simulated breast surgery using a female whole-body musculoskeletal model

This study aimed to use a musculoskeletal model to predict changes in spinal moments following simulated breast surgery. A female full body musculoskeletal model with a fully articulated thoracolumbar spine and independent moveable breast segments was customised for this study. Key findings suggest that the simulated removal of breast tissue (750 g to 1501 g) can reduce the magnitude of lumbar spine extensor moments by >0.05 Nm/kg during walking and jogging. A customised female whole-body musculoskeletal model is capable of providing a first approximation of changes in spinal loading following simulated breast surgery.

Advanced subject-specific neck musculoskeletal modeling unveils sex differences in muscle moment arm and cervical spine loading

Neck pain and injuries are growing healthcare burdens with women having a higher incidence rate and poorer treatment outcomes than males. A better understanding of sex differences in neck biomechanics, foundational for more targeted, effective prevention or treatment strategies, calls for more advanced subject-specific musculoskeletal modeling. Current neck musculoskeletal models are based on generic anatomy, lack subject specificity beyond anthropometric scaling, and are unable to accurately reproduce neck strengths exhibited in vivo without arbitrary muscle force scaling factors or residual torque actuators. In this work, subject-specific neck musculoskeletal models of 23 individuals (11 male, 12 female) were constructed by integrating multi-modality imaging and biomechanical measurements. Each model simulated maximal voluntary neck static exertions in three postures: neck flexion in a neutral posture, flexion in a 40° extended posture, and extension in a 40° flexed posture. Quantitative model validation showed close agreement between model-predicted muscle activation and EMG measurement. The models unveiled that (1) males have greater moment arms in one flexor muscle group and five extensor muscle groups, (2) females exhibited higher cervical spinal compression per unit exertion force in the flexed posture, and (3) the variability of compression force was much greater in females in all three exertions but most notably in the extension with a flexed “dropped head” position. These insights illuminated a plausible pathway from sex differences in neck biomechanics to sex disparities in the risk and prevalence of neck pain.

Muscle Strength Identification Based on Isokinetic Testing and Spine Musculoskeletal Modeling.

Subject-specific spinal musculoskeletal modeling can help understand the spinal loading mechanism during human locomotion. However, existing literature lacks methods to identify the maximum isometric strength of individual spinal muscles. In this study, a muscle strength identification method combining isokinetic testing and musculoskeletal simulations was proposed, and the influence of muscle synergy and intra-abdominal pressure (IAP) on identified spinal muscle strength was further discussed. A multibody dynamic model of the spinal musculoskeletal system was established and controlled by a feedback controller. Muscle strength parameters were adjusted based on the measured isokinetic moments, and muscle synergy vectors and the IAP piston model were further introduced. The results of five healthy subjects showed that the proposed method successfully identified the subject-specific spinal flexor/extensor strength. Considering the synergistic activations of antagonist muscles improved the correlation between the simulated and measured spinal moments, and the introduction of IAP slightly increased the identified spinal extensor strength. The established method is beneficial for understanding spinal loading distributions for athletes and patients with sarcopenia.

Analysis of lumbar spine loading during walking in patients with chronic low back pain and healthy controls: An OpenSim-Based study.

Low back pain (LBP) is one of the most prevalent and disabling disease worldwide. However, the specific biomechanical changes due to LBP are still controversial. The purpose of this study was to estimate the lumbar and lower limb kinematics, lumbar moments and loads, muscle forces and activation during walking in healthy adults and LBP. A total of 18 healthy controls and 19 patients with chronic LBP were tested for walking at a comfortable speed. The kinematic and dynamic data of the subjects were collected by 3D motion capture system and force plates respectively, and then the motion simulation was performed by OpenSim. The OpenSim musculoskeletal model was used to calculate lumbar, hip, knee and ankle joint angle variations, lumbar moments and loads, muscle forces and activation of eight major lumbar muscles. In our results, significant lower lumbar axial rotation angle, lumbar flexion/extension and axial rotation moments, as well as the muscle forces of the four muscles and muscle activation of two muscles were found in patients with LBP than those of the healthy controls (p < 0.05). This study may help providing theoretical support for the evaluation and rehabilitation treatment intervention of patients with LBP.