Estimation of muscle response using three-dimensional musculoskeletal models before impact situation: a simulation study

When car crash experiments are performed using cadavers or dummies, the active muscles' reaction on crash situations cannot be observed. The aim of this study is to estimate muscles' response of the major muscle groups using three-dimensional musculoskeletal model by dynamic simulations of low-speed sled-impact. The three-dimensional musculoskeletal models of eight subjects were developed, including 241 degrees of freedom and 86 muscles. The muscle parameters considering limb lengths and the force-generating properties of the muscles were redefined by optimization to fit for each subject. Kinematic data and external forces measured by motion tracking system and dynamometer were then input as boundary conditions. Through a least-squares optimization algorithm, active muscles' responses were calculated during inverse dynamic analysis tracking the motion of each subject. Electromyography for major muscles at elbow, knee, and ankle joints was measured to validate each model. For low-speed sled-impact crash, experiment and simulation with optimized and unoptimized muscle parameters were performed at 9.4 m/h and 10 m/h and muscle activities were compared among them. The muscle activities with optimized parameters were closer to experimental measurements than the results without optimization. In addition, the extensor muscle activities at knee, ankle, and elbow joint were found considerably at impact time, unlike previous studies using cadaver or dummies. This study demonstrated the need to optimize the muscle parameters to predict impact situation correctly in computational studies using musculoskeletal models. And to improve accuracy of analysis for car crash injury using humanlike dummies, muscle reflex function, major extensor muscles' response at elbow, knee, and ankle joints, should be considered.     

The sensitivity of endpoint forces produced by the extrinsic muscles of the thumb to posture

This study utilizes a biomechanical model of the thumb to estimate the force produced at the thumb-tip by each of the four extrinsic muscles. We used the principle of virtual work to relate joint torques produced by a given muscle force to the resulting endpoint force and compared the results to two separate cadaveric studies. When we calculated thumb-tip forces using the muscle forces and thumb postures described in the experimental studies, we observed large errors. When relatively small deviations from experimentally reported thumb joint angles were allowed, errors in force direction decreased substantially. For example, when thumb posture was constrained to fall within ±15° of reported joint angles, simulated force directions fell within experimental variability in the proximal–palmar plane for all four muscles. Increasing the solution space from ±1° to an unbounded space produced a sigmoidal decrease in error in force direction. Changes in thumb posture remained consistent with a lateral pinch posture, and were generally consistent with each muscle’s function. Altering thumb posture alters both the components of the Jacobian and muscle moment arms in a nonlinear fashion, yielding a nonlinear change in thumb-tip force relative to muscle force. These results explain experimental data that suggest endpoint force is a nonlinear function of muscle force for the thumb, support the continued use of methods that implement linear transformations between muscle force and thumb-tip force for a specific posture, and suggest the feasibility of accurate prediction of lateral pinch force in situations where joint angles can be measured accurately.     

Force, work and power output of lower limb muscles during human maximal-effort countermovement jumping.

The purpose of this study was to simulate human maximal-effort countermovement jumping with a three-dimensional neuromusculoskeletal model. The specific aim was to investigate muscle force, work and power output of major lower limb muscles during the motion. A neuromusculoskeletal model that has nine rigid body segments, 20 degrees of freedom, 32 Hill-type lower limb muscles was developed. The neural activation input signal was represented by a series of step functions with step duration of 0.05 s. The excitation-contraction dynamics of the contractile element, the tissues around the joints to limit the joint range of motion, as well as the foot-ground interaction were implemented. A simulation was started from a standing posture. Optimal pattern of the activation input signal was searched through numerical optimization with a goal of maximizing the height reached by the mass center of body after jumping up. As a result, feasible kinematics, ground reaction force profile and muscle excitation profile were generated. It was found that monoarticular muscles had major contributions of mechanical work and power output, whereas biarticular muscles had minor contributions. Hip adductors, abductors and external rotator muscles were vigorously activated, although their mechanical work and power output was minor because of their limited length change during the motion. Joint flexor muscles such as m. iliopsoas, m. biceps femoris short head and m. tibialis anterior were activated in the beginning of the motion with an effect of facilitating the generation of a countermovement.     

Posture-dependent trunk extensor EMG activity during maximum isometrics exertions in normal male and female subjects.

Posture-dependent trunk function data are important for appropriate normalization of submaximal trunk exertions, and is also necessary to define a more precise and specific use for strength testing in the prevention and diagnosis of spinal disorders. The aim of the current study was to quantify maximal effort trunk muscle extensor activity and trunk isometric extension torque over a functional range of sagittal standing postures. Twenty healthy, young adult male and female subjects performed isometric extension tasks over a sagittal posture range of -20 degrees extension to +50 degrees flexion, in 10 degrees increments. Erector spinae muscle activity was recorded bilaterally at the level of L3 using surface EMG electrodes. Isometric trunk extension torque was measured using a trunk dynamometer. EMG and trunk torque differed significantly between genders, but there were no differences between male and female subjects when the data were normalized with respect to the upright posture. For the combined male and female population, upright posture normalized L3 EMG activity (EMGn) and trunk extension torque (Tn) increased 1.7-fold and 3.5-fold, respectively, over the 70 degrees range of sagittal postures examined. The ratio (Tn/EMGn) increased two-fold (0.83 to 1.67) from -20 degrees extension to +50 degrees flexion, indicating that the neuromuscular efficiency increases with flexion. Trunk extension torque normalized with respect to the upright posture was linearly and positively correlated (r = 0.59, P < 0.001) to similarly normalized L3 EMG activity. This relatively weak correlation suggests that trunk muscle synergism and/or intrinsic muscle length-tension relationships are also modulated by posture. This study provides data that can be used to estimate trunk extensor muscle function over a broad range of sagittal postures. Our findings indicate that appropriate postural normalization of trunk extensor EMG activity is necessary for studies where submaximal trunk exertions are performed over a range of upright postures.     

An examination of the flexion-relaxation phenomenon in the cervical spine in lumbo-pelvic sitting

The flexion-relaxation phenomenon (FRP) is well documented at end-range lumbar spine flexion in both standing and sitting however, the FRP has been insufficiently investigated in cervico-thoracic musculature. The aim of this study was to determine whether the FRP occurs during forward flexion of the neck, in lumbo-pelvic sitting, amongst a pain-free population. Surface electromyography (EMG) was used to measure muscle activation in 20 (10 men, 10 women) asymptomatic subjects in selected cervico-thoracic muscles during four, 5-s phases (upright posture, forward flexion, full flexion and return to upright) while subjects were positioned in lumbo-pelvic sitting. Spinal kinematics were simultaneously measured using an electromagnetic motion tracking device. No FRP was observed in upper trapezius or thoracic erector spinae (T4). When using visual methods to determine the presence/absence of the FRP, five subjects were believed to show evidence of the FRP in the cervical erector spinae. However, when using various non-visual criteria to determine the existence of the FRP, substantial variations (0–13 subjects) were evident. We recommend that criteria based upon relatively large differences in muscle activation should be considered when defining the FRP. These findings are of significance for future investigations examining specific cervical pain disorders.     

Measurement of mandibular movements in patients with temporomandibular disorders and in asymptomatic subjects

The aim of the study was to investigate the range of mandibular movements and to analyze the difference in range of mouth opening, right and left lateral movements, and protrusive movement between patients with clinical diagnoses of temporomandibular disorders and asymptomatic subjects (control group) in a young male population. A total of 240 subjects, aged 19-28, were included in the study. The TMD sample comprised 180 patients (60 patients with muscle disorders; 60 patients with disc displacement with reduction; and 60 patients with muscle disorders and disc displacement with reduction) and was compared with 60 healthy control subjects. All participants were evaluated by the attending dentists at baseline by means of a physical examination of the masticatory system and a history questionnaire which included the Research Diagnostic Criteria for Temporomandibular Disorders (RDC/TMD) Axis I measures. Analysis of variance (ANOVA) with the post hoc Bonferroni criteria showed significant difference in ranges of mandibular movements between and within the groups of asymptomatic subjects and TMD patients for active mouth opening (p = 0.001), right lateral movement (p = 0.002), left lateral movement (p = 0.006), and protrusive movement (p = 0.05). It has been found that there are statistically significant differences in the range of mandibular movements that separate asymptomatic subjects and patients with muscle disorders and disc displacements with reduction in this young male population. However, we cannot conclude that measurements of active mandibular movements can discriminate one group (TMD patients) from the other (asymptomatic subjects), because the mean ranges of these active movements between the groups were measured in clinically "normal" values.     

Upper extremity biomechanics in computer tasks differ by gender

ObjectivesThis laboratory study examined gender differences in upper extremity postures, applied forces, and muscle activity when a computer workstation was adjusted to individual anthropometry according to current guidelines.MethodsFifteen men and 15 women completed five standardized computer tasks: touch-typing, completing a form, editing text, sorting and resizing graphical objects and navigating intranet pages. Subjects worked at a height-adjustable workstation with the keyboard on top of the work surface and the mouse to the right. Subjects repeated the text editing task with the mouse in two other locations: a “high” mouse position, which simulated using a keyboard drawer with the mouse on the primary work surface, and “center” mouse position with the mouse between the keyboard and the body, centered with the body’s center line. Surface electromyography measured muscle activity; electrogoniometric and magnetic motion analysis system measured wrist, forearm and upper arm postures; load-cells measured typing forces; and a force-sensing mouse measured applied forces.ResultsRelative forces applied to the keyboard, normalized muscle activity of two forearm muscles, range of motion for the wrist and shoulder joints and external rotation of the shoulder were higher for women (p < 0.05). When subjects were dichotomized instead by anthropometry (either large/small shoulder width or arm length), the differences in forces, muscle activity of the shoulder and wrist posture and shoulder posture became more pronounced with smaller subjects having higher values. Postural differences between the genders increased in the high mouse position and decreased in the center mouse location.ConclusionsWhen a workstation is adjusted per current guidelines differences in upper extremity force, muscle activity and postural factors still exist between genders. However, these were often stronger when subjects were grouped by anthropometry suggesting that perhaps the computer input devices themselves should be scaled to be more in proportion with the anthropometry and strength of the user.     

The effects of accentuated eccentric loading on strength, muscle hypertrophy, and neural adaptations in trained individuals

The purpose of this study was to compare the strength and neuromuscular adaptations for dynamic constant external resistance (DCER) training and dynamic accentuated external resistance (DAER) training (resistance training employing an accentuated load during eccentric actions). Male subjects active in resistance training were assigned to either a DCER training group (n = 10) or a DAER training group (n = 8) for 9 weeks. Subjects in the DCER group performed 4 sets of 10 repetitions with a load of 75% concentric 1 repetition maximum (RM). Subjects in the DAER group performed 3 sets of 10 repetitions with a concentric load of 75% of 1RM and an eccentric load of approximately 120% of concentric 1RM. Three measures reflecting adaptation of elbow flexors and extensors were recorded pretraining and posttraining: concentric 1RM, muscle cross-sectional area (CSA), and specific tension. Strength was assessed at midtraining periods. No significant changes in muscle CSA were observed in either group. Both training groups experienced significant increases in concentric 1RM and specific tension of both the elbow flexors and extensors, but compared with DCER training, DAER training produced significantly greater increases in concentric 1RM of the elbow extensors. These results suggest that, for some exercises, DAER training may be more effective than DCER training in developing strength within a 9-week training phase. However, for trained subjects, neither protocol is effective in eliciting muscle hypertrophy.     

Static optimization underestimates antagonist muscle activity at the glenohumeral joint: A musculoskeletal modeling study

Static optimization is commonly employed in musculoskeletal modeling to estimate muscle and joint loading; however, the ability of this approach to predict antagonist muscle activity at the shoulder is poorly understood. Antagonist muscles, which contribute negatively to a net joint moment, are known to be important for maintaining glenohumeral joint stability. This study aimed to compare muscle and joint force predictions from a subject-specific neuromusculoskeletal model of the shoulder driven entirely by measured muscle electromyography (EMG) data with those from a musculoskeletal model employing static optimization. Four healthy adults performed six sub-maximal upper-limb contractions including shoulder abduction, adduction, flexion, extension, internal rotation and external rotation. EMG data were simultaneously measured from 16 shoulder muscles using surface and intramuscular electrodes, and joint motion evaluated using video motion analysis. Muscle and joint forces were calculated using both a calibrated EMG-driven neuromusculoskeletal modeling framework, and musculoskeletal model simulations that employed static optimization. The EMG-driven model predicted antagonistic muscle function for pectoralis major, latissimus dorsi and teres major during abduction and flexion; supraspinatus during adduction; middle deltoid during extension; and subscapularis, pectoralis major and latissimus dorsi during external rotation. In contrast, static optimization neural solutions showed little or no recruitment of these muscles, and preferentially activated agonistic prime movers with large moment arms. As a consequence, glenohumeral joint force calculations varied substantially between models. The findings suggest that static optimization may under-estimate the activity of muscle antagonists, and therefore, their contribution to glenohumeral joint stability.     

Vibration-Induced Frontal Postural Reactions in Humans Standing on Unstable Supports of Various Types

The authors studied the influence of the direction of support stability the vibration-induced frontal postural reactions of healthy humans during unilateral vibration of three lateral muscles, namely, the long peroneal muscle, tensor muscle of fascia lata, and abdominal external oblique muscle. The subjects stood on a movable blotter- shaped support. Its base was cylindrical or spheric; its height was 24 cm, and its base radius was 40 cm. The platform turn angle and sagittal and frontal horizontal shift of the upper part of the subjects' bodies were recorded. Reacting to vibrations, the subjects tilted their bodies contralaterally to the vibrated muscle, irrespective of the support type. These reactions were strongest when the subjects stood on a stable support and were weakest when the support was unstable. The reaction values were also influenced by the distance between the vibrated muscle and general center of gravity of the body. As a rule, all other things being equal, vibration-induced reaction strengthened when the stimulated muscle was closer to the center of gravity. The authors concluded that vibration-induced frontal reactions of a standing human depend on the support properties and reflect the features of the participation of the lateral muscles in the standing posture regulation. The information received from the lateral muscular receptors by the control system may be used in different manners, depending on the internal body–support interaction model formed by the CNS.     

Interaction between resistance training and flexibility training in healthy young adults

To test the hypothesis that increases in muscle strength and flexibility are developed by specific training programs, 43 healthy young adults were tested before and after 4 different interventions conducted twice a week for 12 weeks: (a) resistance training only (n = 13); (b) flexibility training only (n = 11); (c) resistance and flexibility training (n = 9); and (d) no intervention (n = 10). There was no change in either strength or flexibility in the control group (p > 0.05). Resistance training improved muscle strength either alone (+14%; effect size = 0.53; p < 0.001) or in combination with flexibility training (+16%; effect size = 0.66; p = 0.032), but did not change flexibility (p = 0.610). Flexibility increased with specific training alone (+33%; p < 0.001) or in combination with resistance training (+18%; p < 0.001). In conclusion, in young, healthy subjects, resistance training alone did not increase flexibility, but resistance training did not interfere with the increase in joint range of motion during flexibility training. These results support the concept that specific training should be employed in order to increase either muscle strength or flexibility.     

Prediction of antagonistic muscle forces using inverse dynamic optimization during flexion/extension of the knee.

This paper examined the feasibility of using different optimization criteria in inverse dynamic optimization to predict antagonistic muscle forces and joint reaction forces during isokinetic flexion/extension and isometric extension exercises of the knee. Both quadriceps and hamstrings muscle groups were included in this study. The knee joint motion included flexion/extension, varus/valgus, and internal/external rotations. Four linear, nonlinear, and physiological optimization criteria were utilized in the optimization procedure. All optimization criteria adopted in this paper were shown to be able to predict antagonistic muscle contraction during flexion and extension of the knee. The predicted muscle forces were compared in temporal patterns with EMG activities (averaged data measured from five subjects). Joint reaction forces were predicted to be similar using all optimization criteria. In comparison with previous studies, these results suggested that the kinematic information involved in the inverse dynamic optimization plays an important role in prediction of the recruitment of antagonistic muscles rather than the selection of a particular optimization criterion. Therefore, it might be concluded that a properly formulated inverse dynamic optimization procedure should describe the knee joint rotation in three orthogonal planes.     

Predicting optimal electrical stimulation for repetitive human muscle activation.

Functional electrical stimulation is the use of electrical currents to activate paralyzed muscles to produce functional movements. Muscle force output must meet or exceed the external load to maintain a posture or produce movements. A mathematical force-fatigue modeling system that predicts muscle force responses during repetitive electrical stimulation has been developed in our laboratory to help identify stimulation patterns that optimize force output for individual subjects. This study tests how well this model predicts the number of contractions that can be maintained above a required force level (successful contractions) during repetitive activation of a muscle. Healthy human quadriceps muscles were tested isometrically on 12 subjects. Data were first collected and used to parameterize the model. Next, the model was used to predict the number of successful contractions that were produced by trains with frequencies ranging from 5 to 100 Hz while the pulse durations and amplitudes were held constant. Finally, three clinically relevant stimulation frequencies were selected and tested to verify the model's predictions. Under these conditions, the model accurately predicted the number of successful contractions for clinically relevant stimulation frequencies. Furthermore, the model appears to have the potential to identify the stimulation frequency that maximizes muscle force output and minimizes fatigue for each subject.     

Age-related changes in posture response under a continuous and unexpected perturbation

Aging is a critical factor to influence the functional performance during daily life. Without an appropriate posture control response when experiencing an unexpected external perturbation, fall may occur. A novel six-degree-of freedom platform with motion control protocol was designed to provide a real-life simulation of unexpected disturbance in order to discriminate the age-related changes of the balance control and the recovery ability. Twenty older adults and 20 healthy young adults participated in the study. The subjects stood barefoot on the novel movable platform, data of the center of mass (COM) excursion, joint rotation angle and electromyography (EMG) were recorded and compared. The results showed that the older adults had similar patterns of joint movement and COM excursion as the young adults during the balance reactive-recovery. However, larger proximal joint rotation in elderly group induced larger COM sway envelop and therefore loss of the compensatory strategy of posture recovery. The old adults also presented a lower muscle power. In order to keep an adequate joint stability preventing from falling, the EMG activity was increased, but the asymmetric pattern might be the key reason of unstable postural response. This novel design of moveable platform and test protocol comprised the computerized dynamic posturography (CDP) demonstrate its value to assess the possible sensory, motor, and central adaptive impairments to balance control and could be the training tool for posture inability person.     

Inter-individual variation in vertebral kinematics affects predictions of neck musculoskeletal models

Experimental studies have found significant variation in cervical intervertebral kinematics (IVK) among healthy subjects, but the effect of this variation on biomechanical properties, such as neck strength, has not been explored. The goal of this study was to quantify variation in model predictions of extension strength, flexion strength and gravitational demand (the ratio of gravitational load from the weight of the head to neck muscle extension strength), due to inter-subject variation in IVK. IVK were measured from sagittal radiographs of 24 subjects (14F, 10M) in five postures: maximal extension, mid-extension, neutral, mid-flexion, and maximal flexion. IVK were defined by the position (anterior-posterior and superior-inferior) of each cervical vertebra with respect to T1 and its angle with respect to horizontal, and fit with a cubic polynomial over the range of motion. The IVK of each subject were scaled and incorporated into musculoskeletal models to create models that were identical in muscle force- and moment-generating properties but had subject-specific kinematics. The effect of inter-subject variation in IVK was quantified using the coefficient of variation (COV), the ratio of the standard deviation to the mean. COV of extension strength ranged from 8% to 15% over the range of motion, but COV of flexion strength was 20–80%. Moreover, the COV of gravitational demand was 80–90%, because the gravitational demand is affected by head position as well as neck strength. These results indicate that including inter-individual variation in models is important for evaluating neck musculoskeletal biomechanical properties.     

Crouched postures reduce the capacity of muscles to extend the hip and knee during the single-limb stance phase of gait

Many children with cerebral palsy walk in a crouch gait that progressively worsens over time, decreasing walking efficiency and leading to joint degeneration. This study examined the effect of crouched postures on the capacity of muscles to extend the hip and knee joints and the joint flexions induced by gravity during the single-limb stance phase of gait. We first characterized representative mild, moderate, and severe crouch gait kinematics based on a large group of subjects with cerebral palsy (N=316). We then used a three-dimensional model of the musculoskeletal system and its associated equations of motion to determine the effect of these crouched gait postures on (1) the capacity of individual muscles to extend the hip and knee joints, which we defined as the angular accelerations of the joints, towards extension, that resulted from applying a 1N muscle force to the model, and (2) the angular acceleration of the joints induced by gravity. Our analysis showed that the capacities of almost all the major hip and knee extensors were markedly reduced in a crouched gait posture, with the exception of the hamstrings muscle group, whose extension capacity was maintained in a crouched posture. Crouch gait also increased the flexion accelerations induced by gravity at the hip and knee throughout single-limb stance. These findings help explain the increased energy requirements and progressive nature of crouch gait in patients with cerebral palsy.     

Support surface related changes in feedforward and feedback control of standing posture

The aim of the study was to investigate the effect of different support surfaces on feedforward and feedback components of postural control. Nine healthy subjects were exposed to external perturbations applied to their shoulders while standing on a rigid platform, foam, and wobble board with eyes open or closed.Electrical activity of nine trunk and leg muscles and displacements of the center of pressure were recorded and analyzed during the time frames typical of feedforward and feedback postural adjustments. Feedforward control of posture was characterized by earlier activation of anterior muscles when the subjects stood on foam compared to a wobble board or a firm surface. In addition, the magnitude of feedforward muscle activity was the largest when the foam was used. During the feedback control, anterior muscles were activated prior to posterior muscles irrespective of the nature of surface. Moreover, the largest muscle activity was seen when the supporting surface was foam. Maximum CoP displacement occurred when subjects were standing on a rigid surface.Altering support surface affects both feedforward and feedback components of postural control. This information should be taken into consideration in planning rehabilitation interventions geared towards improvement of balance.     

An EMG-assisted model of loads on the lumbar spine during asymmetric trunk extensions

An EMG-assisted, low-back, lifting model is presented which simulates spinal loading as a function of dynamic, asymmetric, lifting exertions. The purpose of this study has been to develop a model which overcomes the limitations of previous models including static or isokinetic mechanics, inaccurate predictions of muscle coactivity, static interpretation of myoelectric activity, and physiologically unrealistic or variable muscle force per unit area. The present model predicts individual muscle forces from processed EMG data, normalized as a function of trunk angle and asymmetry, and modified to account for muscle length and velocity artifacts. The normalized EMGs are combined with muscle cross-sectional area and intrinsic strength capacity as determined on a per subject basis, to represent tensile force amplitudes. Dynamic internal and external force vectors are employed to predict trunk moments, spinal compression, lateral and anterior shear forces. Data from 20 subjects performing a total of 2160 exertions showed good agreement between predicted and measured values under all trunk angle, asymmetry, velocity, and acceleration conditions. The design represents a significant step toward accurate, fully dynamic modeling of the low-back in multiple dimensions. The benefits of such a model are the insights provided into the effects of motion induced, muscle co-activity on spinal loading in multiple dimensions.     

Head and shoulder posture affect scapular mechanics and muscle activity in overhead tasks

Forward head and rounded shoulder posture (FHRSP) is theorized to contribute to alterations in scapular kinematics and muscle activity leading to the development of shoulder pain. However, reported differences in scapular kinematics and muscle activity in those with forward head and rounded shoulder posture are confounded by the presence of shoulder pain. Therefore, the purpose of this study was to compare scapular kinematics and muscle activity in individuals free from shoulder pain, with and without FHRSP. Eighty volunteers were classified as having FHRSP or ideal posture. Scapular kinematics were collected concurrently with muscle activity from the upper and lower trapezius as well as the serratus anterior muscles during a loaded flexion and overhead reaching task using an electromagnetic tracking system and surface electromyography. Separate mixed model analyses of variance were used to compare three-dimensional scapular kinematics and muscle activity during the ascending phases of both tasks. Individuals with FHRSP displayed significantly greater scapular internal rotation with less serratus anterior activity, during both tasks as well as greater scapular upward rotation, anterior tilting during the flexion task when compared with the ideal posture group. These results provide support for the clinical hypothesis that FHRSP impacts shoulder mechanics independent of shoulder pain.     

Sarcomere length organization as a design for cooperative function amongst all lumbar spine muscles

The functional design of spine muscles in part dictates their role in moving, loading, and stabilizing the lumbar spine. There have been numerous studies that have examined the isolated properties of these individual muscles. Understanding how these muscles interact and work together, necessary for the prediction of muscle function, spine loading, and stability, is lacking. The objective of this study was to measure sarcomere lengths of lumbar muscles in a neutral cadaveric position and predict the sarcomere operating ranges of these muscles throughout full ranges of spine movements. Sarcomere lengths of seven lumbar muscles in each of seven cadaveric donors were measured using laser diffraction. Using published anatomical coordinate data, superior muscle attachment sites were rotated about each intervertebral joint and the total change in muscle length was used to predict sarcomere length operating ranges. The extensor muscles had short sarcomere lengths in a neutral spine posture and there were no statistically significant differences between extensor muscles. The quadratus lumborum was the only muscle with sarcomere lengths that were optimal for force production in a neutral spine position, and the psoas muscles had the longest lengths in this position. During modeled flexion the extensor, quadratus lumborum, and intertransversarii muscles lengthened so that all muscles operated in the approximate same location on the descending limb of the force-length relationship. The intrinsic properties of lumbar muscles are designed to complement each other. The extensor muscles are all designed to produce maximum force in a mid-flexed posture, and all muscles are designed to operate at similar locations of the force-length relationship at full spine flexion.