The pathogenesis of musculoskeletal disorders is closely associated with the cumulative damage and fatigue failure behavior of fibrous connective tissues under long-term repetitive loading. However, significant technological challenges remain in real-time dynamic monitoring of ligament fatigue life, particularly the lack of efficient computational mechanics modeling frameworks and precise assessment tools adaptable to real-world movement scenarios. The multimodal integrated framework for ligament fatigue life assessment was proposed in this study. First, the high-accuracy subject-specific musculoskeletal models were developed based on individualized medical imaging data. A coupled hyperelastic-viscoelastic constitutive model was incorporated to accurately characterize the nonlinear mechanical behavior of ligamentous tissues and their fatigue damage evolution under cyclic loading. Furthermore, by integrating continuum damage mechanics theory, a time-dependent cumulative damage evolution equation was established to systematically quantify the coupling relationship between fatigue failure probability and dynamic mechanical loading. In the data-driven prediction module, an innovative deep-learning model that integrates kinematic-dynamic coupling was developed. By integrating wearable inertial measurement units, the model enables real-time inversion of ligament loading force-fatigue failure states and prediction of fatigue life. This approach effectively overcomes the limitations of traditional mechanical modeling in long-term, multi-scenario dynamic monitoring, achieving high-precision and minimally invasive fatigue life evaluation of ligaments. The proposed computational framework breaks the static-loading constraints of conventional fatigue testing, achieving the dynamic biomechanical analysis and fatigue life prediction under real movement conditions. This work not only provides novel theoretical insights into the mechanisms and modeling of ligament fatigue damage, but also provides a generalizable tool for biomechanical injury prevention, rehabilitation planning, and soft tissue fatigue analysis in the musculoskeletal system.

Background: During city running or marathon races, shifts in level ground and up-and-down slopes are regularly encountered, resulting in changes in lower limb biomechanics. The longitudinal bending stiffness of the running shoe affects the running performance. Purpose: This research aimed to investigate the biomechanical changes in the lower limbs when transitioning from level ground to an uphill slope under different longitudinal bending stiffness (LBS) levels in running shoes. Methods: Fifteen male amateur runners were recruited and tested while wearing three different LBS running shoes. The participants were asked to pass the force platform with their right foot at a speed of 3.3 m/s ± 0.2. Kinematics data and GRFs were collected synchronously. Each participant completed and recorded ten successful experiments per pair of shoes. Results: The range of motion in the sagittal of the knee joint was reduced with the increase in the longitudinal bending stiffness. Positive work was increased in the sagittal plane of the ankle joint and reduced in the keen joint. The negative work of the knee joint increased in the sagittal plane. The positive work of the metatarsophalangeal joint in the sagittal plane increased. Conclusion: Transitioning from running on a level surface to running uphill, while wearing running shoes with high LBS, could lead to improved efficiency in lower limb function. However, the higher LBS of running shoes increases the energy absorption of the knee joint, potentially increasing the risk of knee injuries. Thus, amateurs should choose running shoes with optimal stiffness when running.

As one of many fundamental sports techniques, the landing maneuver is also frequently used in clinical injury screening and diagnosis. However, the landing patterns are different under different constraints, which will cause great difficulties for clinical experts in clinical diagnosis. Machine learning (ML) have been very successful in solving a variety of clinical diagnosis tasks, but they all have the disadvantage of being black boxes and rarely provide and explain useful information about the reasons for making a particular decision. The current work validates the feasibility of applying an explainable ML (XML) model constructed by Layer-wise Relevance Propagation (LRP) for landing pattern recognition in clinical biomechanics. This study collected 560 groups landing data. By incorporating these landing data into the XML model as input signals, the prediction results were interpreted based on the relevance score (RS) derived from LRP. The interpretation obtained from XML was evaluated comprehensively from the statistical perspective based on Statistical Parametric Mapping (SPM) and Effect Size. The RS has excellent statistical characteristics in the interpretation of landing patterns between classes, and also conforms to the clinical characteristics of landing pattern recognition. The current work highlights the applicability of XML methods that can not only satisfy the traditional decision problem between classes, but also largely solve the lack of transparency in landing pattern recognition. We provide a feasible framework for realizing interpretability of ML decision results in landing analysis, providing a methodological reference and solid foundation for future clinical diagnosis and biomechanical analysis.

Background: Finding the best subset of gait features among biomechanical variables is considered very important because of its ability to identify relevant sports and clinical gait pattern differences to be explored under specific study conditions. This study proposes a new method of metaheuristic optimization-based selection of optimal gait features, and then investigates how much contribution the selected gait features can achieve in gait pattern recognition. Methods: Firstly, 800 group gait datasets performed feature extraction to initially eliminate redundant variables. Then, the metaheuristic optimization algorithm model was performed to select the optimal gait feature, and four classification algorithm models were used to recognize the selected gait feature. Meanwhile, the accuracy results were compared with two widely used feature selection methods and previous studies to verify the validity of the new method. Finally, the final selected features were used to reconstruct the data waveform to interpret the biomechanical meaning of the gait feature. Results: The new method finalized 10 optimal gait features (6 ankle-related and 4-related knee features) based on the extracted 36 gait features (85 % variable explanation) by feature extraction. The accuracy in gait pattern recognition among the optimal gait features selected by the new method (99.81 % ± 0.53 %) was significantly higher than that of the feature-based sorting of effect size (94.69 % ± 2.68 %), the sequential forward selection (95.59 % ± 2.38 %), and the results of previous study. The interval between reconstructed waveform-high and reconstructed waveform-low curves based on the selected feature was larger during the whole stance phase. Significance: The selected gait feature based on the proposed new method (metaheuristic optimization-based selection) has a great contribution to gait pattern recognition. Sports and clinical gait pattern recognition can benefit from population-based metaheuristic optimization techniques. The metaheuristic optimization algorithms are expected to provide a practical and elegant solution for sports and clinical biomechanical feature selection with better economy and accuracy.

Single-leg landing (SL) is often associated with a high injury risk, especially anterior cruciate ligament (ACL) injuries and lateral ankle sprain. This work investigates the relationship between ankle motion patterns (ankle initial contact angle [AICA] and ankle range of motion [AROM]) and the lower limb injury risk during SL, and proposes an optimized landing strategy that can reduce the injury risk. To more realistically revert and simulate the ACL injury mechanics, we developed a knee musculoskeletal model that reverts the ACL ligament to a nonlinear short-term viscoelastic mechanical mechanism (strain ratedependent) generated by the dense connective tissue as a function of strain. Sixty healthy male subjects were recruited to collect biomechanics data during SL. The correlation analysis was conducted to explore the relationship between AICA, AROM, and peak vertical ground reaction force (PVGRF), joint total energy dissipation (TED), peak ankle knee hip sagittal moment, peak ankle inversion angle (PAIA), and peak ACL force (PAF). AICA exhibits a negative correlation with PVGRF (r = -0.591) and PAF (r = -0.554), and a positive correlation with TED (r = 0.490) and PAIA (r = 0.502). AROM exhibits a positive correlation with TED (r = 0.687) and PAIA (r = 0.600). The results suggested that the appropriate increases in AICA (30° to 40°) and AROM (50° to 70°) may reduce the lower limb injury risk. This study has the potential to offer novel perspectives on the optimized application of landing strategies, thus giving the crucial theoretical basis for decreasing injury risk.