Model and in vivo studies on human trunk load partitioning and stability in isometric forward flexions

To resolve the trunk redundancy to determine muscle forces, spinal loads, and stability margin in isometric forward flexion tasks, combined in vivo-numerical model studies was undertaken. It was hypothesized that the passive resistance of both the ligamentous spine and the trunk musculature plays a crucial role in equilibrium and stability of the system. Fifteen healthy males performed free isometric trunk flexions of approximately 40 degrees and approximately 65 degrees +/- loads in hands while kinematics by skin markers and EMG activity of trunk muscles by surface electrodes were measured. A novel kinematics-based approach along with a nonlinear finite element model were iteratively used to calculate muscle forces and internal loads under prescribed measured postures and loads considered in vivo. Stability margin was investigated using nonlinear, linear buckling, and perturbation analyses under various postures, loads and alterations in ligamentous stiffness. Flexion postures significantly increased activity in extensor muscles when compared with standing postures while no significant change was detected in between flexed postures. Compression at the L5-S1 substantially increased from 570 and 771 N in upright posture, respectively, for +/-180 N, to 1912 and 3308 N at approximately 40 degrees flexion, and furthermore to 2332 and 3850 N at approximately 65 degrees flexion. Passive ligamentous/muscle components resisted up to 77% of the net moment. In flexion postures, the spinal stability substantially improved due both to greater passive stiffness and extensor muscle activities so that, under 180 N, no muscle stiffness was required to maintain stability. The co-activity of abdominal muscles and the muscle stiffness were of lesser concern to maintain stability in forward flexion tasks as compared with upright tasks. An injury to the passive system, on one hand, required a substantial compensatory increase in active muscle forces which further increased passive loads and, hence, the risk of injury and fatigue. On the other hand, it deteriorated the system stability which in turn could require greater additional muscle activation. This chain of events would place the entire trunk active-passive system at higher risks of injury, fatigue and instability.     

Trunk response analysis under sudden forward perturbations using a kinematics-driven model

Accurate quantification of the trunk transient response to sudden loading is crucial in prevention, evaluation, rehabilitation and training programs. An iterative dynamic kinematics-driven approach was used to evaluate the temporal variation of trunk muscle forces, internal loads and stability under sudden application of an anterior horizontal load. The input kinematics is hypothesized to embed basic dynamic characteristics of the system that can be decoded by our kinematics-driven approach. The model employs temporal variation of applied load, trunk forward displacement and surface EMG of select muscles measured on two healthy and one chronic low-back pain subjects to a sudden load. A finite element model accounting for measured kinematics, nonlinear passive properties of spine, detailed trunk musculature with wrapping of global extensor muscles, gravity load and trunk biodynamic characteristics is used to estimate the response under measured sudden load. Results demonstrate a delay of ～200 ms in extensor muscle activation in response to sudden loading. Net moment and spinal loads substantially increase as muscles are recruited to control the trunk under sudden load. As a result and due also to the trunk flexion, system stability significantly improves. The reliability of the kinematics-driven approach in estimating the trunk response while decoding measured kinematics is demonstrated. Estimated large spinal loads highlight the risk of injury that likely further increases under larger perturbations, muscle fatigue and longer delays in activation.     

Seated whole body vibrations with high-magnitude accelerations--relative roles of inertia and muscle forces

Reliable computation of spinal loads and trunk stability under whole body vibrations with high acceleration contents requires accurate estimation of trunk muscle activities that are often overlooked in existing biodynamic models. A finite element model of the spine that accounts for nonlinear load- and direction-dependent properties of lumbar segments, complex geometry and musculature of the spine, and dynamic characteristics of the trunk was used in our iterative kinematics-driven approach to predict trunk biodynamics in measured vehicle's seat vibrations with shock contents of about 4g (g: gravity acceleration of 9.8m/s(2)) at frequencies of about 4 and 20Hz. Muscle forces, spinal loads and trunk stability were evaluated for two lumbar postures (erect and flexed) with and without coactivity in abdominal muscles. Estimated peak spinal loads were substantially larger under 4Hz excitation frequency as compared to 20Hz with the contribution of muscle forces exceeding that of inertial forces. Flattening of the lumbar lordosis from an erect to a flexed posture and antagonistic coactivity in abdominal muscles, both noticeably increased forces on the spine while substantially improving trunk stability. Our predictions clearly demonstrated the significant role of muscles in trunk biodynamics and associated risk of back injuries. High-magnitude accelerations in seat vibration, especially at near-resonant frequency, expose the vertebral column to large forces and high risk of injury by significantly increasing muscle activities in response to equilibrium and stability demands.     

Elevation and orientation of external loads influence trunk neuromuscular response and spinal forces despite identical moments at the L5–S1 level

A wide range of loading conditions involving external forces with varying magnitudes, orientations and locations are encountered in daily activities. Here we computed the effect on trunk biomechanics of changes in force location (two levels) and orientation (5 values) in 4 subjects in upright standing while maintaining identical external moment of 15 Nm, 30 N m or 45 Nm at the L5–S1. Driven by measured kinematics and gravity/external loads, the finite element models yielded substantially different trunk neuromuscular response with moderate alterations (up to 24% under 45 Nm moment) in spinal loads as the load orientation varied. Under identical moments, compression and shear forces at the L5–S1 as well as forces in extensor thoracic muscles progressively decreased as orientation of external forces varied from downward gravity (90°) all the way to upward (−25°) orientation. In contrast, forces in local lumbar muscles followed reverse trends. Under larger horizontal forces at a lower elevation, lumbar muscles were much more active whereas extensor thoracic muscle forces were greater under smaller forces at a higher elevation. Despite such differences in activity pattern, the spinal forces remained nearly identical (<6% under 45 Nm moment). The published recorded surface EMG data of extensor muscles trend-wise agreed with computed local muscle forces as horizontal load elevation varied but were overall different from results in both local and global muscles when load orientation altered. Predictions demonstrate the marked effect of external force orientation and elevation on the trunk neuromuscular response and spinal forces and questions attempts to estimate spinal loads based only on consideration of moments at a spinal level.     

Computation of trunk equilibrium and stability in free flexion-extension movements at different velocities

Velocity of movement has been suggested as a risk factor for low-back disorders. The effect of changes in velocity during unconstrained flexion-extension movements on muscle activations, spinal loads, base reaction forces and system stability was computed. In vivo measurements of kinematics and ground reaction forces were initially carried out on young asymptomatic subjects. The collected kinematics of three subjects representing maximum, mean and minimum lumbar rotations were subsequently used in the kinematics-driven model to compute results during the entire movements at three different velocities. Estimated spinal loads and muscle forces were significantly larger in fastest pace as compared to slower ones indicating the effect of inertial forces. Spinal stability was improved in larger trunk flexion angles and fastest movement. Partial or full flexion relaxation of global extensor muscles occurred only in slower movements. Some local lumbar muscles, especially in subjects with larger lumbar flexion and at slower paces, also demonstrated flexion relaxation. Results confirmed the crucial role of movement velocity on spinal biomechanics. Predictions also demonstrated the important role on response of the magnitude of peak lumbar rotation and its temporal variation.     

Trunk strength,muscle activity and spinal loads in maximum isometric flexion and extension exertions: A combined in vivo-computational study

Determination of the trunk maximum voluntary exertion moment capacity and associated internal spinal forces could serve in proper selection of workers for specific occupational task requirements, injury prevention and treatment outcome evaluations. Maximum isometric trunk exertion moments in flexion and extension along with surface EMG of select trunk muscles are measured in 12 asymptomatic subjects. Subsequently and under individualized measured harness-subject forces, kinematics and upper trunk gravity, an iterative kinematics-driven finite element model is used to compute muscle forces and spinal loads in 4 of these subjects. Different co-activity and intra-abdominal pressure levels are simulated. Results indicate significantly larger maximal resistant moments and spinal compression/shear forces in extension exertions than flexion exertions. The agonist trunk muscles reach their maximum force generation (saturation) to greater extent in extension exertions compared to flexion exertions. Local lumbar extensor muscles are highly active in extension exertions and generate most of the internal spinal forces. The maximum exertion attempts produce large spinal compression and shear loads that increase with the antagonist co-activity level but decrease with the intra-abdominal pressure. Intra-abdominal pressure decreases agonist muscle forces in extension exertions but generally increase them in flexion exertions.     

Estimation of trunk muscle forces using the finite element method and in vivo loads measured by telemeterized internal spinal fixation devices.

Direct quantitative measurement of muscle forces is not possible. Forces in the trunk muscles were estimated for standing and flexion of the upper body using three-dimensional, nonlinear finite element models of the lumbar spine with and without an internal spinal fixation device. Muscle forces assumed were two pairs dorsally and one ventrally, each representing several muscles. Muscle forces in the model with internal fixators were varied in discrete steps until the implant loads calculated closely corresponded to those measured in a patient with an instrumented implant. The calculated angles between adjacent lumbar vertebrae were compared with corresponding values measured on X-ray films of a patient as well as with literature values and served as a second criterion for predicting muscle forces. For the model without an implant, the muscle forces of the first model were slightly varied until the lumbar spine shape and the intradiscal pressure were physiological. The abdomen was shown to have a considerable supporting function for flexion.     

Analysis of large compression loads on lumbar spine in flexion and in torsion using a novel wrapping element

Axial compression on the spine could reach large values especially in lifting tasks which also involve large rotations. Experimental and numerical investigations on the spinal multi motion segments in presence of physiological compression loads cannot adequately be carried out due to the structural instability and artefact loads. To circumvent these problems, a novel wrapping cable element is used in a nonlinear finite element model of the lumbosacral spine (L1-S1) to investigate the role of moderate to large compression loads on the lumbar stiffness in flexion and axial moments/rotations. The compression loads up to 2,700 N was applied with no instability or artefact loads. The lumbar stiffness substantially increased under compression force, flexion moment, and axial torque when applied alone. The presence of compression preloads significantly stiffened the load-displacement response under flexion and axial moments/rotations. This stiffening effect was much more pronounced under larger preloads and smaller moments/rotations. Compression preloads also increased intradiscal pressure, facet contact forces, and maximum disc fibre strain at different levels. Forces in posterior ligaments were, however, diminished with compression preload. The significant increase in spinal stiffness, hence, should be considered in biomechanical studies for accurate investigation of the load partitioning, system stability, and fixation systems/disc prostheses.     

A frontal plane model of the lumbar spine subjected to a follower load: implications for the role of muscles

Compression on the lumbar spine is 1000 N for standing and walking and is higher during lifting. Ex vivo experiments show it buckles under a vertical load of 80-100 N. Conversely, the whole lumbar spine can support physiologic compressive loads without large displacements when the load is applied along a follower path that approximates the tangent to the curve of the lumbar spine. This study utilized a two-dimensional beam-column model of the lumbar spine in the frontal plane under gravitational and active muscle loads to address the following question: Can trunk muscle activation cause the path of the internal force resultant to approximate the tangent to the spinal curve and allow the lumbar spine to support compressive loads of physiologic magnitudes? The study identified muscle activation patterns that maintained the lumbar spine model under compressive follower load, resulting in the minimization of internal shear forces and bending moments simultaneously at all lumbar levels. The internal force resultant was compressive, and the lumbar spine model, loaded in compression along the follower load path, supported compressive loads of physiologic magnitudes with minimal change in curvature in the frontal plane. Trunk muscles may coactivate to generate a follower load path and allow the ligamentous lumbar spine to support physiologic compressive loads.     

Role of trunk muscles in generating follower load in the lumbar spine of neutral standing posture

Recently, experimental results have demonstrated that the load carrying capacity of the human spine substantially increases under the follower load condition. Thus, it is essential to prove that a follower load can be generated in vivo by activating the appropriate muscles in order to demonstrate the possibility that the stability of the spinal column could be maintained through a follower load mechanism. The aim of this study was to analyze the coordination of the trunk muscles in order to understand the role of the muscles in generating the follower load. A three-dimensional finite element model of the lumbar spine was developed from T12 to S1 and 117 pairs of trunk muscles (58 pairs of superficial muscles and 59 pairs of deep muscles) were considered. The follower load concept was mathematically represented as an optimization problem. The muscle forces required to generate the follower load were predicted by solving the optimization problem. The corresponding displacements and rotations at all nodes were estimated along with the follower forces, shear forces, and joint moments acting on those nodes. In addition, the muscle forces and the corresponding responses were investigated when the activations of the deep muscles or the superficial muscles were restricted to 75% of the maximum activation, respectively. Significantly larger numbers of deep muscles were involved in the generation of the follower load than the number of superficial muscles, regardless of the restriction on muscle activation. The shear force and the resultant joint moment are more influenced by the change in muscle activation in the superficial muscles. A larger number of deep trunk muscles were activated in order to maintain the spinal posture in the lumbar spine. In addition, the deep muscles have a larger capability to reduce the shear force and the resultant joint moment with respect to the perturbation of the external load or muscle fatigue compared to the superficial muscles.     

A comparative study of two trunk biomechanical models under symmetric and asymmetric loadings

Despite recent advances in modeling of the human spine, simplifying assumptions are still required to tackle complexities. Such assumptions need to be scrutinized to assess their likely impacts on predictions. A comprehensive comparison of muscle forces and spinal loads estimated by a single-joint (L5–S1) optimisation-assisted EMG-driven (EMGAO) and a multi-joint Kinematics-driven (KD) model of the spine under symmetric (symmetric trunk flexion from neutral upright to maximum forward flexion) and asymmetric (holding a load at various heights in the right hand) activities is carried out. Regardless of the task simulated, the KD model predicted greater activities in extensor muscles as compared to the EMGAO model. Such differences in the symmetric tasks was due mainly to the distinct approaches to resolve the redundancy while in the asymmetric tasks they were due also to the different methods used to estimate joint moments. Shear and compression forces were generally higher in the KD model. Differences in predictions between these modeling approaches varied depending on the task simulated and the joint considered in the single-joint EMGAO model. The EMGAO model should incorporate a multi-joint strategy to satisfy equilibrium at different levels while the KD model should benefit from recorded EMG activities of the antagonistic muscles to supplement input measured kinematics.     

Determination of trunk muscle forces for flexion and extension by using a validated finite element model of the lumbar spine and measured in vivo data

Muscle forces stabilize the spine and have a great influence on spinal loads. But little is known about their magnitude. In a former in vitro experiment, a good agreement with intradiscal pressure and fixator loads measured in vivo could be achieved for standing and extension of the lumbar spine. However, for flexion the agreement between in vitro and in vivo measurements was insufficient. In order to improve the determination of trunk muscle forces, a three-dimensional nonlinear finite element model of the lumbar spine with an internal fixation device was created and the same loads were applied as in a previous in vitro experiment. An extensive adaptation process of the model was performed for flexion and extension angles up to 20 degrees and -15 degrees, respectively. With this validated computer model intra-abdominal pressure, preload in the fixators, and a combination of hip- and lumbar flexion angle were varied until a good agreement between analytical and in vivo results was reached for both, intradiscal pressure and bending moments in the fixators. Finally, the fixators were removed and the muscle forces for the intact lumbar spine calculated. A good agreement with the in vivo results could only be achieved at a combination of hip- and lumbar flexion. For the intact spine, forces of 170, 100 and 600 N are predicted in the m. erector spinae for standing, 5 degrees extension and 30 degrees flexion, respectively. The force in the m. rectus abdominus for these body positions is less than 25 N. For more than 10 degrees extension the m. erector spinae is unloaded. The finite element method together with in vivo data allows the estimation of trunk muscle forces for different upper body positions in the sagittal plane. In our patients, flexion of the upper body was most likely a combination of hip- and lumbar spine bending.     

Coupled objective function to study the role of abdominal muscle forces in lifting using the kinematics-driven model

To circumvent the existing shortcoming of optimisation algorithms in trunk biomechanical models, both agonist and antagonist trunk muscle stresses to different powers are introduced in a novel objective function to evaluate the role of abdominal muscles in trunk stability and spine compression. This coupled objective function is introduced in our kinematics-driven finite element model to estimate muscle forces and to identify the role of abdominal muscles in upright standing while lifting symmetrically a weight at different heights. Predictive equations for the compression and buckling forces are developed. Results are also compared with the conventional objective function that neglects abdominal muscle forces. An overall optimal solution involving smaller spinal compression but higher trunk stability is automatically attained when choosing muscle stress powers at and around 3. Results highlight the internal oblique muscle as the most efficient abdominal muscle during the tasks investigated. The estimated relative forces in abdominal muscles are found to be in good agreement with the respective ratios of recorded electromyography activities.     

The effect of control strategies for an active back-support exoskeleton on spine loading and kinematics during lifting

With mechanical loading as the main risk factor for LBP, exoskeletons (EXO) are designed to reduce the load on the back by taking over part of the moment normally generated by back muscles. The present study investigated the effect of an active exoskeleton, controlled using three different control modes (INCLINATION, EMG & HYBRID), on spinal compression forces during lifting with various techniques.Ten healthy male subjects lifted a 15 kg box, with three lifting techniques (free, squat & stoop), each of which was performed four times, once without EXO and once each with the three different control modes. Using inverse dynamics, we calculated L5/S1 joint moments. Subsequently, we estimated spine forces using an EMG-assisted trunk model.Peak compression forces substantially decreased by 17.8% when wearing the EXO compared to NO EXO. However, this reduction was partly, by about one third, attributable to a reduction of 25% in peak lifting speed when wearing the EXO. While subtle differences in back load patterns were seen between the three control modes, no differences in peak compression forces were found. In part, this may be related to limitations in the torque generating capacity of the EXO. Therefore, with the current limitations of the motors it was impossible to determine which of the control modes was best. Despite these limitations, the EXO still reduced both peak and cumulative compression forces by about 18%.     

Effect of load position on muscle forces, internal loads and stability of the human spine in upright postures

A novel kinematics-based approach coupled with a non-linear finite element model was used to investigate the effect of changes in the load position and posture on muscle activity, internal loads and stability margin of the human spine in upright standing postures. In addition to 397 N gravity, external loads of 195 and 380 N were considered at different lever arms and heights. Muscle forces, internal loads and stability margin substantially increased as loads displaced anteriorly away from the body. Under same load magnitude and location, adopting a kyphotic posture as compared with a lordotic one increased muscle forces, internal loads and stability margin. An increase in the height of a load held at a fixed lever arm substantially diminished system stability thus requiring additional muscle activations to maintain the same margin of stability. Results suggest the importance of the load position and lumbar posture in spinal biomechanics during various manual material handling operations.     

Effect of instruction,surface stability,and load intensity on trunk muscle activity

The aim of this study was to assess the effect of verbal instruction, surface stability, and load intensity on trunk muscle activity levels during the free weight squat exercise. Twelve trained males performed a free weight squat under four conditions: (1) standing on stable ground lifting 50% of their 1-repetition maximum (RM), (2) standing on a BOSU balance trainer lifting 50% of their 1-RM, (3) standing on stable ground lifting 75% of their 1-RM, and (4) receiving verbal instructions to activate the trunk muscles followed by lifting 50% of their 1-RM. Surface EMG activity from muscles rectus abdominis (RA), external oblique (EO), transversus abdominis/internal oblique (TA/IO), and erector spinae (ES) were recorded for each condition and normalized for comparisons. Muscles RA, EO, and TA/IO displayed greater peak activity (39–167%) during squats with instructions compared to the other squat conditions (P = 0.04–0.007). Peak EMG activity of muscle ES was greater for the 75% 1-RM condition than squats with instructions or lifting 50% of 1-RM (P = 0.04–0.02). The results indicate that if the goal is to enhance EMG activity of the abdominal muscles during a multi-joint squat exercise then verbal instructions may be more effective than increasing load intensity or lifting on an unstable surface. However, in light of other research, conscious co-activation of the trunk muscles during the squat exercise may lead to spinal instability and hazardous compression forces in the lumbar spine.     

Effects of lumbo-pelvic rhythm on trunk muscle forces and disc loads during forward flexion: A combined musculoskeletal and finite element simulation study

Previous in-vivo studies suggest that the ratio of total lumbar rotation over pelvic rotation (lumbo-pelvic rhythm) during trunk sagittal movement is essential to evaluate spinal loads and discriminate between low back pain and asymptomatic population. Similarly, there is also evidence that the lumbo-pelvic rhythm is key for evaluation of realistic muscle and joint reaction forces and moments predicted by various computational musculoskeletal models. This study investigated the effects of three lumbo-pelvic rhythms defined based on in-vivo measurements on the spinal response during moderate forward flexion (60°) using a combined approach of musculoskeletal modeling of the upper body and finite element model of the lumbosacral spine. The muscle forces and joint loads predicted by the musculoskeletal model, together with the gravitational forces, were applied to the finite element model to compute the disc force and moment, intradiscal pressure, annular fibers strain, and load-sharing. The results revealed that a rhythm with high pelvic rotation and low lumbar flexion involves more global muscles and increases the role of the disc in resisting spinal loads, while its counterpart, with low pelvic rotation, recruits more local muscles and engages the ligaments to lower the disc loads. On the other hand, a normal rhythm that has balanced pelvic and lumbar rotations yields almost equal disc and ligament load-sharing and results in more balanced synergy between global and local muscles. The lumbo-pelvic rhythm has less effect on the intradiscal pressure and annular fibers strain. This work demonstrated that the spinal response during forward flexion is highly dependent on the lumbo-pelvic rhythm. It is therefore, essential to adapt this parameter instead of using the default values in musculoskeletal models for accurate prediction of muscle forces and joint reaction forces and moments. The findings provided by this work are expected to improve knowledge of spinal response during forward flexion, and are clinically relevant towards low back pain treatment and disc injury prevention.     

Effect of intervertebral translational flexibilities on estimations of trunk muscle forces,kinematics, loads,and stability

Due to the complexity of the human spinal motion segments, the intervertebral joints are often simulated in the musculoskeletal trunk models as pivots thus allowing no translational degrees of freedom (DOFs). This work aims to investigate, for the first time, the effect of such widely used assumption on trunk muscle forces, spinal loads, kinematics, and stability during a number of static activities. To address this, the shear deformable beam elements used in our nonlinear finite element (OFE) musculoskeletal model of the trunk were either substantially stiffened in translational directions (SFE model) or replaced by hinge joints interconnected through rotational springs (HFE model). Results indicated that ignoring intervertebral translational DOFs had in general low to moderate impact on model predictions. Compared with the OFE model, the SFE and HFE models predicted generally larger L4–L5 and L5–S1 compression and shear loads, especially for tasks with greater trunk angles; differences reached ～15% for the L4–L5 compression, ～36% for the L4–L5 shear and ～18% for the L5–S1 shear loads. Such differences increased, as location of the hinge joints in the HFE model moved from the mid-disc height to either the lower or upper endplates. Stability analyses of these models for some select activities revealed small changes in predicted margin of stability. Model studies dealing exclusively with the estimation of spinal loads and/or stability may, hence with small loss of accuracy, neglect intervertebral translational DOFs at smaller trunk flexion angles for the sake of computational simplicity.     

Effect of load position on muscle forces,internal loads and stability of the human spine in upright postures

A novel kinematics-based approach coupled with a non-linear finite element model was used to investigate the effect of changes in the load position and posture on muscle activity, internal loads and stability margin of the human spine in upright standing postures. In addition to 397 N gravity, external loads of 195 and 380 N were considered at different lever arms and heights. Muscle forces, internal loads and stability margin substantially increased as loads displaced anteriorly away from the body. Under same load magnitude and location, adopting a kyphotic posture as compared with a lordotic one increased muscle forces, internal loads and stability margin. An increase in the height of a load held at a fixed lever arm substantially diminished system stability thus requiring additional muscle activations to maintain the same margin of stability. Results suggest the importance of the load position and lumbar posture in spinal biomechanics during various manual material handling operations.