Multibody dynamics-based musculoskeletal modeling for gait analysis: a systematic review.

Beyond qualitative assessment, gait analysis involves the quantitative evaluation of various parameters such as joint kinematics, spatiotemporal metrics, external forces, and muscle activation patterns and forces. Utilizing multibody dynamics-based musculoskeletal (MSK) modeling provides a time and cost-effective non-invasive tool for the prediction of internal joint and muscle forces. Recent advancements in the development of biofidelic MSK models have facilitated their integration into clinical decision-making processes, including quantitative diagnostics, functional assessment of prosthesis and implants, and devising data-driven gait rehabilitation protocols. Through an extensive search and meta-analysis of over 116 studies, this PRISMA-based systematic review provides a comprehensive overview of different existing multibody MSK modeling platforms, including generic templates, methods for personalization to individual subjects, and the solutions used to address statically indeterminate problems. Additionally, it summarizes post-processing techniques and the practical applications of MSK modeling tools. In the field of biomechanics, MSK modeling provides an indispensable tool for simulating and understanding human movement dynamics. However, limitations which remain elusive include the absence of MSK modeling templates based on female anatomy underscores the need for further advancements in this area.

Positional Analysis of Assisting Muscles for Handling-Assisted Exoskeletons.

In order to better design handling-assisted exoskeletons, it is necessary to analyze the biomechanics of human hand movements. In this study, Anybody Modeling System (AMS) simulation was used to analyze the movement state of muscles during human handling. Combined with surface electromyography (sEMG) experiments, specific analysis and verification were carried out to obtain the position of muscles that the human body needs to assist during handling. In this study, the simulation and experiment were carried out for the manual handling process. A treatment group and an experimental group were set up. This study found that the vastus medialis muscle, vastus lateralis muscle, latissimus dorsi muscle, trapezius muscle, deltoid muscle and triceps brachii muscle require more energy in the process of handling, and it is reasonable and effective to combine sEMG signals with the simulation of the musculoskeletal model to analyze the muscle condition of human movement.

Evaluation of ground reaction forces and centers of pressure predicted by AnyBody Modeling System during load reaching/handling activities and effects of the prediction errors on model-estimated spinal loads.

Full-body and lower-extremity human musculoskeletal models require feet ground reaction forces (GRFs) and centers of pressure (CoPs) as inputs to predict muscle forces and joint loads. GRFs/CoPs are traditionally measured via floor-mounted forceplates that are usually restricted to research laboratories thus limiting their applicability in real occupational and clinical setups. Alternatively, GRFs/CoPs can be estimated via inverse dynamic approaches as also implemented in the Anybody Modeling System (AnyBody Technology, Aalborg, Denmark). The accuracy of Anybody in estimating GRFs/CoPs during load-handling/reaching activities and the effect of its prediction errors on model-estimated spinal loads remain to be investigated. Twelve normal- and over-weight individuals performed total of 480 static load-handling/reaching activities while measuring (by forceplates) and predicting (by AnyBody) their GRFs/CoPs. Moreover, the effects of GRF/CoP prediction errors on the estimated spinal loads were evaluated by inputting measured or predicted GRFs/CoPs into subject-specific musculoskeletal models. Regardless of the subject groups (normal-weight or overweight) and tasks (load-reaching or load-handling), results indicated great agreements between the measured and predicted GRFs (normalized root-mean-squared error, nRMSEs < 14% and R2 > 0.90) and between their model-estimated spinal loads (nRMSEs < 14% and R2 > 0.83). These agreements were good but relatively less satisfactory for CoPs (nRMSEs < 17% and 0.57 < R2 < 0.68). The only exception, requiring a more throughout investigation, was the situation when the ground-foot contact was significantly reduced during the activity. It appears that occupational/clinical investigations performed in real workstation/clinical setups with no access to forceplates may benefit from the AnyBody GRF/CoP prediction tools for a wide range of load-reaching/handling activities.

Exploring the impact of body mass change on fatigue and activity of the muscular system during daily routine.

BACKGROUND: Correct body weight is one of the factors that allows you to maintain a properly functioning body. Abnormal body weight can cause muscle tissue remodelling, affecting activity and muscle fatigue. Changes in the muscular system can cause occurrence of functional limitations. OBJECTIVE: To determine the effect of weight change on fatigue and activity of the muscular system during daily activities. METHODS: The evaluation of musculoskeletal functioning was based on the results of computer simulations conducted in the AnyBody Modeling System. The following activities were analysed: standing, sitting down and getting up from a chair, holding and lifting an object, and walking. The simulations of the activities were carried out using averaged kinematic data, and by changing body mass in the range of 50 kg to 100 kg by increments of 2 kg, to map different nutritional status from excessive thinness to extreme obesity. Identification of loads in the musculoskeletal system was based on solving an inverse dynamics problem and then the estimation of muscle force values using static optimization. The simulation results allowed to determine the value of muscle fatigue and the level of muscle activity. RESULTS: For activities (i.e., standing, walking, sitting down and getting up from a chair) it was observed that the value of muscle fatigue increases with increasing body mass. However, for activities that cause more load on the musculoskeletal system, i.e. lifting and holding an object, the highest value of muscle fatigue was observed in underweight individuals. CONCLUSION: The change in body weight alters the functioning of the muscular system and thus the ability to perform activities. It was shown that in case of underweight, overweight or obese people, abnormal body weight can be the reason for occurrence of difficulties in performing the activities of lifting and holding a 20 kg object, as well as walking.

Effect of strengthening and weakening of abdominal and dorsal muscles on lumbar spine loads in parents of disabled children.

Pain in the lower part of the back is one of the most common chronic illnesses globally. This work aimed to determine the impact of the reinforcement of particular groups of abdominal and dorsal muscles on the loads exerted on the lumbar section of the spine in 30 mothers of children with motor disabilities. An optical Ariel Performance Analysis System recorded and processed the kinematics data of everyday activities. Tests investigating the effects of the strengthening or weakening of abdominal and dorsal muscles on loading in the lumbar section of the spine utilized the AnyBody Modelling System. Input data for the simulations included mean values of body positions, while the effects of strengthening or weakening of muscles were simulated in the muscle forces model by introducing different values for muscle physiological cross-sectional area (PCSA). Simulations used decreasing or increasing PCSA values of abdominal muscles and the erector spinae. The analysis involved component and resultant force values on the lumbosacral joint (L5-S1) of the spine and intra-abdominal pressure values. The highest reduction of the resultant reaction value in L5-S1 was observed in the simulations that increased the PCSA of the transverse abdominal (TrA). Indeed, a double increase in the TrA cross-section caused a reduction of the resultant reaction in L5-S1 by 30% and the anterior-posterior and proximal-distal forces by approximately 20-30%. Increased PCSA of the erector spinae exerted higher loads on the spine. These results indicate that strengthening weakened abdominal muscles, particularly TrA, in parents of children with motor disabilities reduces lower spinal loads during daily activities.

Model-Based Biomechanical Exoskeleton Concept Optimization for a Representative Lifting Task in Logistics.

Occupational exoskeletons are a promising solution to prevent work-related musculoskeletal disorders (WMSDs). However, there are no established systems that support heavy lifting to shoulder height. Thus, this work presents a model-based analysis of heavy lifting activities and subsequent exoskeleton concept optimization. Six motion sequences were captured in the laboratory for three subjects and analyzed in multibody simulations with respect to muscle activities (MAs) and joint forces (JFs). The most strenuous sequence was selected and utilized in further simulations of a human model connected to 32 exoskeleton concept variants. Six simulated concepts were compared concerning occurring JFs and MAs as well as interaction loads in the exoskeleton arm interfaces. Symmetric uplifting of a 21 kg box from hip to shoulder height was identified as the most strenuous motion sequence with highly loaded arms, shoulders, and back. Six concept variants reduced mean JFs (spine: >70%, glenohumeral joint: >69%) and MAs (back: >63%, shoulder: >59% in five concepts). Parasitic loads in the arm bracing varied strongly among variants. An exoskeleton design was identified that effectively supports heavy lifting, combining high musculoskeletal relief and low parasitic loads. The applied workflow can help developers in the optimization of exoskeletons.

Ergonomic Assessment of a Lower-Limb Exoskeleton through Electromyography and Anybody Modeling System.

The aim of this study was to determine the muscle load reduction of the upper extremities and lower extremities associated with wearing an exoskeleton, based on analyses of muscle activity (electromyography: EMG) and the AnyBody Modeling System (AMS). Twenty healthy males in their twenties participated in this study, performing bolting tasks at two working heights (60 and 85 cm). The muscle activities of the upper trapezius (UT), middle deltoid (MD), triceps brachii (TB), biceps brachii (BB), erector spinae (ES), biceps femoris (BF), rectus femoris (RF), and tibialis anterior (TA) were measured by EMG and estimated by AMS, respectively. When working at the 60 cm height with the exoskeleton, the lower extremity muscle (BF, RF, TA) activities of EMG and AMS decreased. When working at the 85 cm height, the lower extremity muscle activity of EMG decreased except for TA, and those of AMS decreased except for RF. The muscle activities analyzed by the two methods showed similar patterns, in that wearing the exoskeleton reduced loads of the lower extremity muscles. Therefore, wearing an exoskeleton can be recommended to prevent an injury. As the results of the two methods show a similar tendency, the AMS can be used.

Hand Motion Capture from a 3D Leap Motion Controller for a Musculoskeletal Dynamic Simulation.

The AnyBody Modeling System™ (AMS) is a musculoskeletal software simulation solution using inverse dynamics analysis. It enables the determination of muscle and joint forces for a given bodily motion. The recording of the individual movement and the transfer into the AMS is a complex and protracted process. Researches indicated that the contactless, visual Leap Motion Controller (LMC) provides clinically meaningful motion data for hand tracking. Therefore, the aim of this study was to integrate the LMC hand motion data into the AMS in order to improve the process of recording a hand movement. A Python-based interface between the LMC and the AMS, termed ROSE Motion, was developed. This solution records and saves the data of the movement as Biovision Hierarchy (BVH) data and AnyScript vector files that are imported into the AMS simulation. Setting simulation parameters, initiating the calculation automatically, and fetching results is implemented by using the AnyPyTools library from AnyBody. The proposed tool offers a rapid and easy-to-use recording solution for elbow, hand, and finger movements. Features include animation, cutting/editing, exporting the motion, and remote controlling the AMS for the analysis and presentation of musculoskeletal simulation results. Comparing the motion tracking results with previous studies, covering problems when using the LMC limit the correctness of the motion data. However, fast experimental setup and intuitive and rapid motion data editing strengthen the use of marker less systems as the herein presented compared to marker based motion capturing.

Relevance of a novel external dynamic distraction device for treating back pain.

Low back pain is a common, expensive, and disabling condition in industrialized countries. There is still no consensus for its ideal management. Believing in the beneficial effect of traction, we developed a novel external dynamic distraction device. The purpose of this work was to demonstrate that external distraction allows limiting the pressure exerted in standing-up position on the lower intervertebral discs. Numerical and cadaveric studies were used as complementary approaches. Firstly, we implemented the device into a numerical model of a validated musculoskeletal software (Anybody Modeling System) and we calculated the lower disc pressure while traction forces were applied. Secondly, we performed an anatomical study using a non-formalin preserved cadaver placed in a sitting position. A pressure sensor was placed in the lower discs under fluoroscopic control through a Jamshidi needle. The intradiscal pressure was then measured continuously at rest while applying a traction force of 200 N. Both numerical and cadaveric studies demonstrated a decrease in intradiscal pressures after applying a traction force with the external device. Using the numerical model, we showed that tensile forces below 500 N in total were sufficient. The application of higher forces seems useless and potentially deleterious. External dynamic distraction device is able to significantly decrease the intradiscal pressure in a sitting or standing position. However, the therapeutic effects need to be proven using clinical studies.

Assessing the efficiency of exoskeletons in physical strain reduction by biomechanical simulation with AnyBody Modeling System.

Introduction: Recently, many industrial exoskeletons for supporting workers in heavy physical tasks have been developed. However, the efficiency of exoskeletons with regard to physical strain reduction has not been fully proved, yet. Several laboratory and field studies have been conducted, but still more data, that cannot be obtained solely by behavioral experiments, are needed to investigate effects on the human body. Methods: This paper presents an approach to extend laboratory and field research with biomechanical simulations using the AnyBody Modeling System. Based on a dataset recorded in a laboratory experiment with 12 participants using the exoskeleton Paexo Shoulder in an overhead task, the same situation was reproduced in a virtual environment and analyzed with biomechanical simulation. Results: Simulation results indicate that the exoskeleton substantially reduces muscle activity and joint reaction forces in relevant body areas. Deltoid muscle activity and glenohumeral joint forces in the shoulder were decreased between 54 and 87%. Simultanously, no increases of muscle activity and forces in other body areas were observed. Discussion: This study demonstrates how a simulation framework could be used to evaluate changes in internal body loads as a result of wearing exoskeletons. Biomechanical simulation results widely agree with experimental measurements in the previous laboratory experiment and supplement such by providing an insight into effects on the human musculoskeletal system. They confirm that Paexo Shoulder is an effective device to reduce physical strain in overhead tasks. The framework can be extended with further parameters, allowing investigations for product design and evaluation.

Torus obstacle method as a wrapping approach of the deltoid muscle group for humeral abduction in musculoskeletal simulation.

Musculoskeletal models of the shoulder complex are valuable research aids to investigate tears of the supraspinatus and the resulting mechanical impact during abduction of the humerus. One of the major contributors to this motion is the deltoid muscle group and for this, an accurate modeling of the lines of action is indispensable. The aim of this work was to utilize a torus obstacle wrapping approach for the deltoids of an existing shoulder model and assess the feasibility of the approach during humeral abduction. The shoulder model from the AnyBody™ modeling system was used as a platform. The size of the tori is based on a magnetic resonance imaging (MRI) approach and several kinematic couplings are implemented to determine the trajectories of the tori during abduction. To assess the model behavior, the moment arms of the virtual muscle elements and the resultant glenohumeral joint reaction force (GHJF) were compared with reference data from the literature during abduction of the humerus in the range 20°-120°. The root mean square error for the anterior, lateral and posterior part between the simulated muscle elements and reference data from the literature was 3.9, 1.7 and 5.8 mm, respectively. The largest deviation occurred on the outer elements of the muscle groups, with 12.6, 10.4 and 20.5 mm, respectively. During abduction, there is no overlapping of the muscle elements and these are in continuous contact with the torus obstacles, thus enabling a continuous force transmission. This results in a rising trend of the resultant GHJF. The torus obstacle approach as a wrapping method for the deltoid muscles provides a guided muscle pathing by simultaneously approximating the curvature of the deltoid muscle. The results from the comparison of the simulated moment arms and the resultant GHJF are in accordance with those in the literature in the range 20°-120° of abduction. Although this study shows the strength of the torus obstacle as a wrapping approach, the method of fitting the tori according to MRI data was not suitable. A cadaver study is recommended to better validate and mathematically describe the torus approach.

A numerical study to determine the effect of strengthening and weakening of the transversus abdominis muscle on lumbar spine loads.

This work aimed to determine the influence of strengthening or weakening of the transversus abdominis (TrA) muscle on loads in the lumbar spine using musculoskeletal modelling methods. The input kinematic data of two positions (a standing position and a position during a sitting-down task) were angles in the elbow joint (0°;4°), shoulder joint (0°;3°), hip joint (0°;75°) and knee joint (0°;69°) as well as the torso tilt angle (0°;32°). It was shown that a change in the TrA physiological cross section area (PCSA) has a crucial impact on lumbar spine loads (2xTrA PCSA causes a reduction in the force in joint L5-S1 by 11% for a standing position and by 25% for a sitting-down position) and value of intra-abdominal pressure.

Regarding loads after spinal fusion, every level should be seen separately: a musculoskeletal analysis.

INTRODUCTION: The number of spinal fusion surgeries is steadily increasing and biomechanical consequences are still in debate. The aim of this study is to provide biomechanical insights into the sagittal balance of the spine and to compare spinal load before and after spinal fusion. METHOD: The joint reaction forces of 52 patients were analyzed in proximo-distal and antero-posterior direction from the levels T12-L1 to L5-S1 using musculoskeletal simulations. RESULTS: In 104 simulations, pre-surgical forces were equal to post-surgical. The levels L4-L5 and T12-L1, however, showed increased spinal forces compression forces with higher sagittal displacement. Improved restauration of sagittal balance was accompanied by lower spinal load. AP shear stress, interestingly decreased with sagittal imbalance. CONCLUSION: Imbalanced spines have a risk of increased compression forces at Th12-L1. L4-L5 always has increased spinal loads. These slides can be retrieved under Electronic Supplementary Material.

Musculoskeletal modeling of human lower limb during normal walking, one-legged forward hopping and side jumping: Comparison of measured EMG and predicted muscle activity patterns.

This study focused on comparing muscle activities predicted by the Musculoskeletal Modeling System with EMG from ten healthy subjects who performed normal walking, one-legged forward hopping and side jumping. Eight EMG electrodes measured the activity of eight right leg muscles. Specific thresholds per muscle were applied on the EMG prior comparison. These thresholds were determined by equalizing the duration of EMG to AMS muscle activity. Three graph variables, number of onsets, offsets and hills were used to quantify the level of agreement by using Cohen׳s kappa analysis. The Pearson correlation coefficient was also calculated as a result comparison. Overall, visual inspection showed comparable activity patterns. However, when quantifying them some differences became apparent. The mean level of agreement of all tests was <0.20, meaning poor agreement. Pearson correlation showed better agreement compared to kappa analysis. In general, a more prescribed motion like FH and SJ showed a better agreement than NW. This explorative study shows that there are distinct differences between the model and EMG pattern. Those differences can be attributed to inevitable modeling limitation within the AMS framework like miscalculating the knee net moment, absence of co-contraction, simplified knee joint. Moreover, the delay between EMG and AMS has a clear effect on the comparison and this delay is obviously missing in the model. Despite those differences, this study can serve as a baseline measurement allowing progress in scientific work in order to reduce uncertainties with the aim to generate more reliable and robust musculoskeletal models in a valid manner.