A mass-flowing muscle model with shape restrictive soft tissues: correlation with sonoelastography

Prediction of jaw opening function after mandibular reconstruction using subject-specific musculoskeletal modelling

BACKGROUND

Mandibular reconstruction patients often suffer abnormalities in the mandibular kinematics. In silico simulations, such as musculoskeletal modelling, can be used to predict post-operative mandibular kinematics. It is important to validate the mandibular musculoskeletal model and analyse the factors influencing its accuracy.

OBJECTIVES

To investigate the jaw opening-closing movements after mandibular reconstruction, as predicted by the subject-specific musculoskeletal model, and the factors influencing its accuracy.

METHODS

Ten mandibular reconstruction patients were enrolled in this study. Cone-beam computed tomography images, mandibular movements, and surface electromyogram signals were recorded preoperatively. A subject-specific mandibular musculoskeletal model was established to predict surgical outcomes using patient-averaged muscle parameter changes as model inputs. Jaw bone geometry was replaced by surgical planning results, and the muscle insertion sites were registered based on the non-rigid iterative closest point method. The predicted jaw kinematic data were validated based on 6-month post-operative measurements. Correlations between the prediction accuracy and patient characteristics (age, pathology and surgical scope) were further analysed.

RESULTS

The root mean square error (RMSE) for lower incisor displacement was 31.4%, and the error for peak magnitude of jaw opening was 4.9 mm. Age, post-operative infection and radiotherapy influenced the prediction accuracy. The amount of masseter detachment showed little correlation with jaw opening.

CONCLUSION

The mandibular musculoskeletal model successfully predicted short-range jaw opening functions after mandibular reconstruction. It provides a novel surgical planning method to predict the risk of developing trismus.

Determining a musculoskeletal system's pre-stretched state using continuum-mechanical forward modelling and joint range optimization

The subject-specific range of motion (RoM) of a musculoskeletal joint system is balanced by pre-tension levels of individual muscles, which affects their contraction capability. Such an inherent pre-tension or pre-stretch of muscles is not measureable with in vivo experiments. Using a 3D continuum mechanical forward simulation approach for motion analysis of the musculoskeletal system of the forearm with 3 flexor and 2 extensor muscles, we developed an optimization process to determine the muscle fibre pre-stretches for an initial arm position, which is given human dataset. We used RoM values of a healthy person to balance the motion in extension and flexion. The performed sensitivity study shows that the fibre pre-stretches of the m. brachialis, m. biceps brachii and m. triceps brachii with 91 % dominate the objective flexion ratio, while m. brachiradialis and m. anconeus amount 7.8 % and 1.2 % . Within the multi-dimensional space of the surrogate model, 3D sub-spaces of primary variables, namely the dominant muscles and the global objective, flexion ratio, exhibit a path of optimal solutions. Within this optimal path, the muscle fibre pre-stretch of two flexors demonstrate a negative correlation, while, in contrast, the primary extensor, m. triceps brachii correlates positively to each of the flexors. Comparing the global optimum with four other designs along the optimal path, we saw large deviations, e.g., up to 15∘ in motion and up to 40% in muscle force. This underlines the importance of accurate determination of fibre pre-stretch in muscles, especially, their role in pathological muscular disorders and surgical applications such as free muscle or tendon transfer.

Muscle Strength Identification Based on Isokinetic Testing and Spine Musculoskeletal Modeling

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

Construction and Simulation of Biomechanical Model of Human Hip Joint Muscle-Tendon Assisted by Elastic External Tendon by Hill Muscle Model

Based on the Hill muscle model (HMM), a biomechanical model of human hip muscle tendon assisted by elastic external tendon (EET) was preliminarily established to investigate and analyze the biomechanical transition between the hip joint (HJ) and related muscle tendons. Using the HMM, the optimal muscle fiber length and muscle force scaling variables were introduced by means of constrained optimization problems and were optimized. The optimized HMM was constructed with human parameters of 170 cm and 70 kg. The biomechanical model simulation test of the hip muscle tendon was performed in the automatic dynamic analysis of mechanical systems (ADAMS) software to analyze and optimize the changes in the root mean square error (RMSE), biological moment, muscle moment distribution coefficient (MDC), muscle moment, muscle force, muscle power, and mechanical work of the activation curves of the hip major muscle, iliopsoas muscle, rectus femoris muscle, and hamstring muscle under analyzing the optimized HMM and under different EET auxiliary stiffnesses from the joint moment level, joint level, and muscle level, respectively. It was found that the trends of the output joint moment of the optimized HMM and the biological moment of the human HJ were basically the same, r 2 = 0.883 and RMSE = 0.18 Nm/kg, and the average metabolizable energy consumption of the HJ was (243.77 ± 1.59) J. In the range of 35%∼65% of gait cycle (GC), the auxiliary moment showed a significant downward trend with the increase of EET stiffness, when the EET stiffness of the human body was less than 200 Nm/rad, the biological moment of the human HJ gradually decreased with the increase of EET stiffness, and the MDC of the iliopsoas and hamstring muscles gradually decreased; when the EET stiffness was greater than 200 Nm/rad, the increase of the total moment of the extensor muscles significantly increased, the MDC of the gluteus maximus and rectus muscles gradually increased, and the gluteus maximus and hamstring muscle moments and muscle forces gradually increased; the results show that the optimized muscle model based on Hill can reflect the law of human movement and complete the simulation test of HJ movements, which provides a new idea for the analysis of energy migration in the musculoskeletal system of the lower limb.

Embodiment of intra-abdominal pressure in a flexible multibody model of the trunk and the spinal unloading effects during static lifting tasks

The role of intra-abdominal pressure (IAP) in spinal load reduction has remained controversial, partly because previous musculoskeletal models did not introduce the pressure generating mechanism. In this study, an integrated computational methodology is proposed to combine the IAP change with core muscle activations. An ideal gas relationship was introduced to calculate pressure distribution within the abdominal cavity. Additionally, based on flexible multibody dynamics, a muscle membrane element was developed by incorporating the muscular fiber deformation, inter-fiber stiffness, and volume constancy. This element was then utilized in discretizing the diaphragm and transversus abdominis, forming an IAP-muscle coupling system of the abdominal cavity. Based on this methodology, a forward dynamic simulation of spinal flexion was presented to examine the unloading effect of abdominal breathing. The results confirm that core muscle contraction during the abdominal breathing cycle can substantially reduce the forces of spinal compression together with trunk extensor muscles, and this effect is more pronounced when the IAP increase is produced by contraction of the transversus abdominis. This unloading effect still holds even with the co-activation of other abdominal muscles, providing a potential choice when designing trunk movements during weight-lifting tasks.