Musculoskeletal Modelling Approach Research Articles

Estimation of joint moments during ball throwing in simulated altered gravity conditions using a musculoskeletal modelling approach: preliminary results

Estimation of joint moments during ball throwing in simulated altered gravity conditions using a musculoskeletal modelling approach: preliminary results

Knee joint contact load associated by balance control for stance leg during taekwondo front kick: A musculoskeletal modeling approach

Knee joint contact load associated by balance control for stance leg during taekwondo front kick: A musculoskeletal modeling approach

Quantifying walking speeds in relation to ankle biomechanics on a real-time interactive gait platform: a musculoskeletal modeling approach in healthy adults.

Background: Given the inherent variability in walking speeds encountered in day-to-day activities, understanding the corresponding alterations in ankle biomechanics would provide valuable clinical insights. Therefore, the objective of this study was to examine the influence of different walking speeds on biomechanical parameters, utilizing gait analysis and musculoskeletal modelling. Methods: Twenty healthy volunteers without any lower limb medical history were included in this study. Treadmill-assisted gait-analysis with walking speeds of 0.8m/s and 1.1m/s was performed using the Gait Real-time Analysis Interactive Lab (GRAIL®). Collected kinematic data and ground reaction forces were processed via the AnyBody® modeling system to determine ankle kinetics and muscle forces of the lower leg. Data were statistically analyzed using statistical parametric mapping to reveal both spatiotemporal and magnitude significant differences. Results: Significant differences were found for both magnitude and spatiotemporal curves between 0.8m/s and 1.1m/s for the ankle flexion (p < 0.001), subtalar force (p < 0.001), ankle joint reaction force and muscles forces of the M. gastrocnemius, M. soleus and M. peroneus longus ( = 0.05). No significant spatiotemporal differences were found between 0.8m/s and 1.1m/s for the M. tibialis anterior and posterior. Discussion: A significant impact on ankle joint kinematics and kinetics was observed when comparing walking speeds of 0.8m/s and 1.1m/s. The findings of this study underscore the influence of walking speed on the biomechanics of the ankle. Such insights may provide a biomechanical rationale for several therapeutic and preventative strategies for ankle conditions.

Prediction of trunk muscle activation and spinal forces in adolescent idiopathic scoliosis during simulated trunk motion: A musculoskeletal modelling study

Due to lack of reference validation data, the common strategy in characterizing adolescent idiopathic scoliosis (AIS) by musculoskeletal modelling approach consists in adapting structure and parameters of validated body models of adult individuals with physiological alignments. Until now, only static postures have been replicated and investigated in AIS subjects. When aiming to simulate trunk motion, two critical factors need consideration: how distributing movement along the vertebral motion levels (lumbar spine rhythm), and if neglecting or accounting for the contribution of the stiffness of the motion segments (disc stiffness). The present study investigates the effect of three different lumbar spine rhythms and absence/presence of disc stiffness on trunk muscle imbalance in the lumbar region and on intervertebral lateral shear at different levels of the thoracolumbar/lumbar scoliotic curve, during simulated trunk motions in the three anatomical planes (flexion/extension, lateral bending, and axial rotation). A spine model with articulated ribcage previously developed in AnyBody software and adapted to replicate the spinal alignment in AIS subjects is employed. An existing dataset of 100 subjects with mild and moderate scoliosis is exploited. The results pointed out the significant impact of lumbar spine rhythm configuration and disc stiffness on changes in the evaluated outputs, as well as a relationship with scoliosis severity. Unfortunately, no optimal settings can be identified due to lack of reference validation data. According to that, extreme caution is recommended when aiming to adapt models of adult individuals with physiological alignments to adolescent subjects with scoliotic deformity.

How can saddle height changes the risk injuries of lumbar during the cycling: Kinematics and musculoskeletal modeling approach

How can saddle height changes the risk injuries of lumbar during the cycling: Kinematics and musculoskeletal modeling approach

Effect of obesity on spinal loads during load-reaching activities: A subject- and kinematics-specific musculoskeletal modeling approach

Obesity has been associated to increase the risk of low back disorders. Previous musculoskeletal models simulating the effect of body weight on intervertebral joint loads have assumed identical body postures for obese and normal-weight individuals during a given physical activity. Our recent kinematic-measurement studies, however, indicate that obese individuals adapt different body postures (segmental orientations) than normal-weight ones when performing load-reaching activities. The present study, therefore, used a subject- and kinematics-specific musculoskeletal modeling approach to compare spinal loads of nine normal-weight and nine obese individuals each performing twelve static two-handed load-reaching activities at different hand heights, anterior distances, and asymmetry angles (total of 12 tasks × 18 subjects = 216 model simulations). Each model incorporated personalized muscle architectures, body mass distributions, and full-body kinematics for each subject and task. Results indicated that even when accounting for subject-specific body kinematics obese individuals experienced significantly larger (by ∼38% in average) L5-S1 compression (2305 ± 468 N versus 1674 ± 337 N) and shear (508 ± 111 N versus 705 ± 150 N) loads during all reaching activities (p < 0.05 for all hand positions). This average difference of ∼38% was similar to the results obtained from previous modeling investigations that neglected kinematics differences between the two weight groups. Moreover, there was no significant interaction effect between body weight and hand position on the spinal loads; indicating that the effect of body weight on L5-S1 loads was not dependent on the position of hands. Postural differences alone appear, hence, ineffective in compensating the greater spinal loads that obese people experience during reaching activities.

Do Interlimb Knee Joint Loading Asymmetries Persist throughout Stance during Uphill Walking Following Total Knee Arthroplasty?

The purpose of this study was to determine differences in total (TCF), medial compartment (MCF), and lateral compartment (LCF) tibiofemoral joint compressive forces and related muscle forces between replaced and non-replaced limbs during level and uphill walking at an incline of 10°. A musculoskeletal modeling and simulation approach using static optimization was used to determine the muscle forces and TCF, MCF, and LCF for 25 patients with primary TKA. A statistical parametric mapping repeated-measures ANOVA was conducted on knee compressive forces and muscle forces using statistical parametric mapping. Greater TCF, MCF, and LCF values were observed throughout the loading response, mid-stance, and late stance during uphill walking. During level walking, knee extensor muscle forces were greater throughout the first 50% of the stance during level walking, yet greater during uphill walking during the last 50% of the stance. Conversely, knee flexor muscle forces were greater through the loading response and push-off phases of the stance. No between-limb differences were observed for compressive or muscle forces, suggesting that uphill walking may promote a more balanced loading of replaced and non-replaced limbs. Additionally, patients with TKA appear to rely on the hamstrings muscle group during the late stance for knee joint control, thus supporting uphill walking as an effective exercise modality to improve posterior chain muscle strength.

Joint loading topography during occupational tasks - A musculoskeletal modeling approach to substantiate ergonomic recommendations

BackgroundTo improve ergonomic recommendations and decrease work-related musculoskeletal disorders, thorough documentation of full-body, joint loading topography during occupational tasks is needed. Therefore, the purpose of this study was to document full-body internal joint loading topography in terms of estimated joint contact forces during occupational tasks. In addition, this internal loading topography was also compared to loading proxies (e.g., external joint moments) commonly used to assess injury risk during occupational tasks. Methods3D motion capture and ground reaction forces were measured while 20 participants performed ten occupational tasks. A musculoskeletal modeling workflow with a detailed spine was used to calculate internal joint loading in terms of contact forces, and their association with external joint moments was evaluated. FindingsLifting 10 kg from the ground imposed the highest full-body internal joint loading compared to all other lifting tasks, while lifting 10 kg from hip height to shoulder height imposed the lowest internal joint loading. Only during occupational tasks involving standing upright posture, loading proxies did correlate well with internal joint loading. InterpretationThe modeling workflow and the internal joint loading topography could inform ergonomic recommendations on optimized load distribution across different anatomical regions.

Effect of microgravity on mechanical loadings in lumbar spine at various postures: a numerical study

The aim of this study was to quantitatively analyze the mechanical change of spinal segments (disc, muscle, and ligament) at various postures under microgravity using a full-body musculoskeletal modeling approach. Specifically, in the lumbar spine, the vertebra were modeled as rigid bodies, the intervertebral discs were modeled as 6-degree-of-freedom joints with linear force-deformation relationships, the disc swelling pressure was deformation dependent, the ligaments were modeled as piecewise linear elastic materials, the muscle strength was dependent on its functional cross-sectional area. The neutral posture and the “fetal tuck” posture in microgravity (short as “Neutral 0G” and “Fetal Tuck 0G”, in our simulation, the G constant was set to 0 for simulating microgravity), and for comparison, the relaxed standing posture in 1G and 0G gravity (short as “Neutral 1G” and “Standing 0G”) were simulated. Compared to values at Neutral 1G, the mechanical response in the lower spine changed significantly at Neutral 0G. For example, the compressive forces on lumbar discs decreased 62–70%, the muscle forces decreased 55.7–92.9%, while disc water content increased 7.0–10.2%, disc height increased 2.1–3.0%, disc volume increased 6.4–9.3%, and ligament forces increased 59.5–271.3% at Neutral 0G. The fetal tuck 0G reversed these changes at Neutral 0G back toward values at Neutral 1G, with magnitudes much larger than those at Neutral 1G. Our results suggest that microgravity has significant influences on spinal biomechanics, alteration of which may increase the risks of disc herniation and degeneration, muscle atrophy, and/or ligament failure.

Effect of sprinting velocity on anterior cruciate ligament and knee load during sidestep cutting.

The purpose of the study was to investigate the effect of an increase in sprinting velocity on the anterior cruciate ligament (ACL) load, knee joint load, and activation of femoral muscles using the musculoskeletal modeling approach. Fourteen high school male athletes were recruited (age: 17.4 ± 0.7years, height: 1.75 ± 0.04m, weight: 73.3 ± 8.94kg), with the right foot dominant and physical activity level of about 3-4h per day. The kinematics, kinetics, and co-contraction index (CCI) of the extensors and flexors of the right leg's femoral muscles were calculated. The anterior cruciate ligament load was estimated using the musculoskeletal modeling method. In the results, it was observed that the anterior cruciate ligament load (p < 0.017) increased as sidestep cutting velocity increased, resulting in increased adduction (p < 0.017) and the internal rotation moment of the knee joint. This was significantly higher than when sprinting at a similar velocity. The co-contraction index result, which represents the balanced activation of the femoral extensor and flexor muscles, showed a tendency of decrement with increasing sprinting velocity during sidestep cutting (p < 0.017), whereas no significant differences were observed when running at different sprinting conditions. Therefore, we postulate that factors such as knee joint shear force, extended landing posture with increasing sprinting velocity, internal rotation moment, and femoral muscle activity imbalance influence the increase of anterior cruciate ligament load during a sidestep cutting maneuver.

Using musculoskeletal modelling to estimate knee joint loading pre and post high tibial osteotomy

BackgroundBoth medial knee osteoarthritis and associated varus alignment have been proposed to alter knee joint loading and consequently overloading the medial compartment. Individuals with knee osteoarthritis and varus deformity are candidates for coronal plane corrective surgery, high tibial osteotomy. This study evaluated knee loading and contact location for a control group, a pre-surgery cohort and the same cohort 12 months post-surgery using a musculoskeletal modelling approach. MethodsJoint kinematics during gait were measured in 30 knee osteoarthritis patients, before and after high tibial osteotomy, and 28 healthy adults. Using a musculoskeletal model that incorporated patient-specific mechanical tibial femoral angle, the resulting muscle, ligament, and contact forces were calculated and the medial - lateral condyle load distribution was analysed. FindingsSurgery changed medial compartment contact force throughout stance relative to pre-surgery. This reduction in medial compartment contact force pre- vs post-HTO is observed despite a significant increase in post-surgery walking speed compared to pre-HTO, where increased speed is typically associated with increased joint loading. InterpretationThis study has estimated the effects of high tibial osteotomy on knee loading using a generic model that incorporates a detailed knee model to better understand tibiofemoral contact loading. The findings support the aim of surgery to unload the medial knee compartment and lateralise joint contact forces.

Effects of auxetic shoe on lumbar spine kinematics and kinetics during gait and drop vertical jump by a combined in vivo and modeling investigation

The present study examined the effects of auxetic shoes on the biomechanics of the spine, as compared to barefoot and conventional shoe conditions, during gait and drop vertical jump (DVJ) activities using a combined in vivo and musculoskeletal modeling approach. Motion and force-plate data as well as electromyographic (EMG) activities of select trunk muscles of 11 individuals were collected during foregoing activities. In DVJ activity, two main phases of first landing (FL) and second landing (SL) were studied. In the FL phase of DVJ noticeable alternations were observed when auxetic shoes were used. That is, compared to the conventional footwear condition, smaller EMG activities in extensor muscles (by ~ 16–29%, p < 0.001), smaller anterior–posterior (AP) distance between the center of pressure of ground reaction force and heel (by ~ 19%, p = 0.002), generally larger maximal hip, knee, and ankle flexion angles (p < 0.005) and finally smaller maximal L5-S1 compression force and maximal external moment (by ~ 12 and 8%, respectively, p < 0.001) were obtained by wearing auxetic shoes. Our results, therefore, indicate that using auxetic shoes can reduce load on the lumbar spine during high-demanding activities such as vertical jump and thus may decrease the musculoskeletal risk of injuries during these activities.

Effect of lumbar muscle atrophy on the mechanical loading change on lumbar intervertebral discs

The objective of this study was to quantitatively analyze the effect of lumbar spinal muscle atrophy on the compressive (perpendicular to the upper surface of the disc) and shear (parallel to the upper surface of the disc in the anterior-posterior direction) forces change on lumbar intervertebral discs using a full body musculoskeletal modeling approach. Muscles atrophy was modeled with reduction of the functional cross-sectional area (FCSA) of the muscles. Compressive and shear forces under two levels of lumbar muscle atrophy (20% and 40%) at eight daily postures (lying on back, seating slouched, seating straight, standing, standing flexed (36°), standing lift a 20 kg weight close to chest, standing lift a 20 kg weight flexed (38°), and standing lift a 20 kg weight with arm stretched) were analyzed. There was small increase in compressive forces on lumbar discs with muscle atrophy at most postures except lying and sitting straight. The maximum increase of compressive forces on lumbar discs were 23 N (6%), 28 N (5%), 34 N (2%), 71 N (6%), 89 N (4%), and 190 N (10%) with 20% atrophy, and 66 N (19%), 77 N (12%), 98 N (6%), 169 N (14%), 256 N (12%), 501 N (24%) with 40% atrophy at seating slouched, standing, standing flexed, standing lift close, standing lift flexed, and standing stretched arm, respectively. The shear force did not change significantly on lumbar discs with muscle atrophy. This study is important for understanding the biomechanical mechanisms of how lumbar muscle atrophy may affect the lumbar IVD health.

Ground reaction force and joint moment estimation during gait using an Azure Kinect-driven musculoskeletal modeling approach

BackgroundGait analysis is burdened by time and equipment costs, interpretation, and accessibility of three-dimensional motion analysis systems. Evidence suggests growing adoption of gait testing in the shift toward evidence-based medicine. Further developments addressing these barriers will aid its efficacy in clinical practice. Previous research aiming to develop gait analysis systems for kinetics estimation using the Kinect V2 have provided promising results yet modified approaches using the latest hardware may further aid kinetics estimation accuracy Research questionCan a single Azure Kinect sensor combined with a musculoskeletal modeling approach provide kinetics estimations during gait similar to those obtained from marker-based systems with embedded force platforms? MethodsTen subjects were recruited to perform three walking trials at their normal speed. Trials were recorded using an eight-camera optoelectronic system with two embedded force plates and a single Azure Kinect sensor. Marker and depth data were both used to drive a musculoskeletal model using the AnyBody Modeling System. Predicted kinetics from the Azure Kinect-driven model, including ground reaction force (GRF) and joint moments, were compared to measured values using root meansquared error (RMSE), normalized RMSE, Pearson correlation, concordance correlation, and statistical parametric mapping ResultsHigh to very high correlations were observed for anteroposterior GRF (ρ = 0.889), vertical GRF (ρ = 0.940), and sagittal hip (ρ = 0.805) and ankle (ρ = 0.876) moments. RMSEs were 1.2 ± 2.2 (%BW), 3.2 ± 5.7 (%BW), 0.7 ± 0.1.3 (%BWH), and 0.6 ± 1.0 (%BWH) SignificanceThe proposed approach using the Azure Kinect provided higher accuracy compared to previous studies using the Kinect V2 potentially due to improved foot tracking by the Azure Kinect. Future studies should seek to optimize ground contact parameters and focus on regions of error between predicted and measured kinetics highlighted currently for further improvements in kinetic estimations.

Monitoring joint mechanics in anterior cruciate ligament reconstruction using depth sensor-driven musculoskeletal modeling and statistical parametric mapping

The incidence of anterior cruciate ligament injury and reconstruction (ACLR) may set the stage for the development of early onset osteoarthritis in these patients. Development of accessible quantitative motion capture methodologies for recurrent monitoring of knee joint loading during daily activities following ACLR is necessary. This study aimed to compare lower extremity kinetics between ACLR affected limbs, ACLR unaffected limbs, and dominant limbs of healthy control subjects during over-ground gait and stair ascent using a single depth sensor-driven musculoskeletal modeling approach. No meaningful differences were found between groups during over-ground gait in any kinetic variables. When subjected to a stair ascent task, both ACLR limbs showed greater hip extension and internal rotation moments compared to control subjects at approximately 72–79% stance. This was coincident with greater knee flexion moments in both ALCR limbs compared to control. The absence of differences during over-ground gait but presence of compensatory strategies during stair ascent, suggests task dependent recovery in this cohort who were tested at least 1-year following surgery. Importantly, this was determined using a portable low-cost motion capture method which may be attractive to professionals in sports medicine for recurrent monitoring following ACLR.

Single leg hop for distance symmetry masks lower limb biomechanics: time to discuss hop distance as decision criterion for return to sport after ACL reconstruction?

BackgroundWe evaluated the lower limb status of athletes after anterior cruciate ligament reconstruction (ACLR) during the propulsion and landing phases of a single leg hop for distance (SLHD) task after...

Force Production Patterns of Muscles Surrounding Knee During Running and Cutting Maneuvers: A Musculoskeletal Modeling Approach

OBJECTIVES The purpose of this study was to investigate the force production patterns of individual muscles surrounding the knee during running (RUN) and cutting (CUT) tasks.METHODS Thirteen women (24.2±3.5 yrs, 162.8±6.0 cm, 55.3±6.2 kg) performed a series of running and cutting tasks. Running and cutting motions were recorded using a motion capture system and ground reaction force (GRF) was recorded using a force plate. Three-dimensional knee angle, ground reaction force, and knee joint moment were calculated using Visual3D software. OpenSim musculoskeletal modeling software was used to calculate the force of individual muscles including the medial hamstring, biceps femoris long head, biceps femoris short head, rectus femoris, vastus medialis, vastus lateralis, gastrocnemius medialis, and gastrocnemius lateralis. All data were analyzed for loading response (or weight acceptance), mid-stance, and final push-off periods, respectively and were compared between two tasks.RESULTS At loading response: external rotation angle, medial and vertical GRFs, and valgus moment for the CUT task were greater than those of the RUN task. Compared to the RUN task, the CUT task showed: 1) an increase in lateral hamstring muscle force at weight acceptance, 2) a decrease in hamstring muscle force and an increase in medial vastus muscle force at mid-stance, and 3) an increase in lateral gastrocnemius muscle force at final push-off.CONCLUSIONS Selective force production patterns of muscles surrounding the knee seem to offset the external load caused by the cutting motion. We anticipate that our results will provide basic data for future training programs designed to prevent noncontact knee injuries.

From Stoop to Squat: A Comprehensive Analysis of Lumbar Loading Among Different Lifting Styles.

Lifting up objects from the floor has been identified as a risk factor for low back pain, whereby a flexed spine during lifting is often associated with producing higher loads in the lumbar spine. Even though recent biomechanical studies challenge these assumptions, conclusive evidence is still lacking. This study therefore aimed at comparing lumbar loads among different lifting styles using a comprehensive state-of-the-art motion capture-driven musculoskeletal modeling approach. Thirty healthy pain-free individuals were enrolled in this study and asked to repetitively lift a 15 kg-box by applying 1) a freestyle, 2) a squat and 3) a stoop lifting technique. Whole-body kinematics were recorded using a 16-camera optical motion capture system and used to drive a full-body musculoskeletal model including a detailed thoracolumbar spine. Continuous as well as peak compressive, anterior-posterior shear and total loads (resultant load vector of the compressive and shear load vectors) were calculated based on a static optimization approach and expressed as factor body weight (BW). In addition, lumbar lordosis angles and total lifting time were calculated. All parameters were compared among the lifting styles using a repeated measures design. For each lifting style, loads increased towards the caudal end of the lumbar spine. For all lumbar segments, stoop lifting showed significantly lower compressive and total loads (−0.3 to −1.0BW) when compared to freestyle and squat lifting. Stoop lifting produced higher shear loads (+0.1 to +0.8BW) in the segments T12/L1 to L4/L5, but lower loads in L5/S1 (−0.2 to −0.4BW). Peak compressive and total loads during squat lifting occurred approximately 30% earlier in the lifting cycle compared to stoop lifting. Stoop lifting showed larger lumbar lordosis range of motion (35.9 ± 10.1°) than freestyle (24.2 ± 7.3°) and squat (25.1 ± 8.2°) lifting. Lifting time differed significantly with freestyle being executed the fastest (4.6 ± 0.7 s), followed by squat (4.9 ± 0.7 s) and stoop (5.9 ± 1.1 s). Stoop lifting produced lower total and compressive lumbar loads than squat lifting. Shear loads were generally higher during stoop lifting, except for the L5/S1 segment, where anterior shear loads were higher during squat lifting. Lifting time was identified as another important factor, considering that slower speeds seem to result in lower loads.

Kinetic changes associated with extended knee landings following anterior cruciate ligament reconstruction in females

ObjectivesTo determine the relationship between knee flexion excursion symmetry and lower extremity kinematics, kinetics, and muscle, joint, and ligament forces in females 1–3 years after ACL reconstruction. DesignCross-sectional. SettingLaboratory. ParticipantsTwenty-one, college-aged females. Main outcome measuresLower extremity kinetics and kinematics, including estimated muscle, tibiofemoral, and ligament forces were assessed using 3D motion analysis and a musculoskeletal modeling approach. Participants demonstrating greater than 10% asymmetry in knee flexion excursion were classified as landing with an “extended knee”. Group and between-limb differences were compared. ResultsTen participants were classified as landing with an “extended knee” on the involved limb, while eleven exhibited a symmetric landing pattern. Participants landing with an “extended knee” demonstrated reduced knee extension moment and quadriceps force in the involved limb (p < 0.05). ConclusionsThese findings indicate that an “extended knee” landing pattern was associated with reduced knee extension moment and quadriceps muscle force in females 1–3 years after ACL reconstruction. This may represent an altered strategy that clinicians may choose to identify and address during rehabilitation.

A review of musculoskeletal modelling of human locomotion.

Locomotion through the environment is important because movement provides access to key resources, including food, shelter and mates. Central to many locomotion-focused questions is the need to understand internal forces, particularly muscle forces and joint reactions. Musculoskeletal modelling, which typically harnesses the power of inverse dynamics, unites experimental data that are collected on living subjects with virtual models of their morphology. The inputs required for producing good musculoskeletal models include body geometry, muscle parameters, motion variables and ground reaction forces. This methodological approach is critically informed by both biological anthropology, with its focus on variation in human form and function, and mechanical engineering, with a focus on the application of Newtonian mechanics to current problems. Here, we demonstrate the application of a musculoskeletal modelling approach to human walking using the data of a single male subject. Furthermore, we discuss the decisions required to build the model, including how to customize the musculoskeletal model, and suggest cautions that both biological anthropologists and engineers who are interested in this topic should consider.