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.

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.

Background: Understanding the internal load characteristics of the knee joint is essential for investigating unilateral knee injuries associated with running. This study examined the differences in the location and magnitude of von Mises stress in the internal structures of bilateral knee joints during the stance phase of gait following 10 km running at submaximal speeds. Methods: A healthy male recreational runner participated in this study. We employed a synergistic approach, integrating subject-specific knee joint angles, reaction forces, and moments derived from musculoskeletal modeling to inform and drive the finite element (FE) modeling of running. This methodology ensured a detailed and accurate representation of knee joint mechanics. The peak stresses of the bilateral knee menisci, tibial cartilage, and five main ligaments were quantified using a FE model during the stance phase. Results: The meniscus, tibial cartilage, anterior (ACL), posterior cruciate ligament (PCL), medial (MCL), lateral collateral ligament (LCL) and experienced larger loads in the non-dominant limb across most phases of stance. Additionally, fatigue elevated the peak loading on the non-dominant limb’s ACL, PCL, and LCL during the gait stance phase but diminished the load on these ligaments in the dominant knee joint. For Patellar ligament (PL), the non-dominant side showed maximal stress at initial contact, while the dominant side dominated during the remaining stance phases. Conclusions: This proof-of-concept substantially enhances our understanding of the impact of running-induced fatigue on bilateral knee loading. It provides a detailed analysis of factors leading to unilateral knee overload during extended running. These insights are essential in formulating targeted strategies to reduce injury risks.