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: Single-leg landing (SL) is an essential technique in sports such as basketball, soccer, and volleyball, which is often associated with a high risk of knee-related injury. The ankle motion pattern plays a crucial role in absorbing the load shocks during SL, but the effect on the knee joint is not yet clear. This work aims to explore the effects of different ankle plantarflexion angles during SL on the risk of knee-related injury. Methods: Thirty healthy male subjects were recruited to perform SL biomechanics tests, and one standard subject was selected to develop the finite element model of foot-ankle-knee integration. The joint impact force was used to evaluate the impact loads on the knee at various landing angles. The internal load forces (musculoskeletal modeling) and stress (finite element analysis) around the knee joint were simulated and calculated to evaluate the risk of knee-related injury during SL. To more realistically revert and simulate the anterior cruciate ligament (ACL) injury mechanics, we developed a knee musculoskeletal model that reverts the ACL ligament to a nonlinear short-term viscoelastic mechanical mechanism (strain rate-dependent) generated by the dense connective tissue as a function of strain. Results: As the ankle plantarflexion angle increased during landing, both the peak knee vertical impact force (p = 0.001) and ACL force (p = 0.001) decreased significantly. The maximum von Mises stress of ACL, meniscus, and femoral cartilage decreased as the ankle plantarflexion angle increased. The overall range of variation in ACL stress was small and was mainly distributed in the femoral and tibial attachment regions, as well as in the mid-lateral region. Conclusion: The current findings revealed that the use of larger ankle plantarflexion angles during landing may be an effective solution to reduce knee impact load and the risk of rupture of the medial femoral attachment area in the ACL. The findings of this study have the potential to offer novel perspectives in the optimized application of landing strategies, thus giving crucial theoretical backing for decreasing the risk of knee-related injury.

Previous research has primarily focused on evaluating the activity of individual muscles in dancers, often neglecting their synergistic interactions. Investigating the differences in lower limb muscle synergy during landing between dancers and healthy controls will contribute to a comprehensive understanding of their neuromuscular control patterns. This study enrolled 22 Latin dancers and 22 healthy participants, who performed a task involving landing from a 30 cm high platform. The data were collected using Vicon systems, force plates, and electromyography (EMG). The processed EMG data were subjected to non-negative matrix factorization (NNMF) for decomposition, followed by classification using K-means clustering algorithm and Pearson correlation coefficients. Three synergies were extracted for both Latin dancers and healthy participants. Synergy 1 showed increased contributions from the tibialis anterior (p < 0.001) and medial gastrocnemius (p = 0.024) in Latin dancers compared to healthy participants. Synergy 3 highlighted significantly greater contributions from the vastus lateralis in healthy participants compared to Latin dancers (p = 0.039). This study demonstrates that Latin dancers exhibit muscle synergies similar to those observed in healthy controls, revealing specific adjustments in the tibialis anterior and medial gastrocnemius muscles among dancers. This research illustrates how dancers optimize control strategies during landing tasks, offering a novel perspective for comprehensively understanding dancers’ neuromuscular control patterns.

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.