- 1Faculty of Sports Science, Ningbo University, Ningbo, China
- 2Faculty of Engineering, University of Pannonia, Veszprem, Hungary
- 3School of Mechanical and Aerospace Engineering, Nanyang Technological University, Singapore, Singapore
- 4Department of Biomedical Engineering, Faculty of Engineering, The Hong Kong Polytechnic University, Kowloon, Hong Kong SAR, China
Background: Ankle sprains are prevalent in sports, often causing complex injuries to the lateral ligaments. Among these, anterior talofibular ligament (ATFL) injuries constitute 85%, and calcaneofibular ligament (CFL) injuries comprise 35%. Despite conservative treatment, some ankle sprain patients develop chronic lateral ankle instability (CLAI). Thus, this study aimed to investigate stress response and neural control alterations during landing in lateral ankle ligament injury patients.
Method: This study recruited twenty individuals from a Healthy group and twenty CLAI patients performed a landing task using relevant instruments to collect biomechanical data. The study constructed a finite element (FE) foot model to examine stress responses in the presence of laxity of the lateral ankle ligaments. The lateral ankle ligament was modeled as a hyperelastic composite structure with a refined representation of collagen bundles and ligament laxity was simulated by adjusting material parameters. Finally, the validity of the finite element model is verified by a high-speed dual fluoroscopic imaging system (DFIS).
Result: CLAI patients exhibited earlier Vastus medialis (p < 0.001) and tibialis anterior (p < 0.001) muscle activation during landing. The FE analysis revealed that with laxity in the ATFL, the peak von Mises stress in the fifth metatarsal was 20.74 MPa, while with laxity in the CFL, it was 17.52 MPa. However, when both ligaments were relaxed simultaneously, the peak von Mises stress surged to 21.93 MPa. When the ATFL exhibits laxity, the CFL is subjected to a higher stress of 3.84 MPa. Conversely, when the CFL displays laxity, the ATFL experiences a peak von Mises stress of 2.34 MPa.
Conclusion: This study found that changes in the laxity of the ATFL and the CFL are linked to shifts in metatarsal stress levels, potentially affecting ankle joint stability. These alterations may contribute to the progression towards CLAI in individuals with posterolateral ankle ligament injuries. Additionally, significant muscle activation pattern changes were observed in CLAI patients, suggesting altered neural control strategies post-ankle ligament injury.
1 Introduction
Lateral Ankle Ligament (LAL) injuries are among the most prevalent injuries affecting the ankle joint, with occurrences noted across various populations and levels of athletic activity (Swanson et al., 2022). According to research statistics, LAL injuries represent approximately 85% or more of all ankle injuries (Ferran and Maffulli, 2006). Among these, a prevalent type of LAL injury is the anterior talofibular ligament (ATFL) injury, often occurring concomitantly with injuries to the calcaneofibular ligament (CFL) or both the calcaneofibular ligament and posterior talofibular ligament (PTFL) (DiGiovanni et al., 2004; Maffulli et al., 2008). Injuries to the ATFL represent 85% of these injuries, while injuries to the CFL account for 35%. Extending from the anterior margin of the lateral fibular ankle to the neck of the talus, ATFL is also recognized as the structurally weakest among the lateral ankle ligaments (Khawaji and Soames, 2015). Despite a 60%–80% healing rate with conservative treatment for ankle sprains, cases where injuries fail to heal properly may lead to chronic lateral ankle instability (CLAI) (Hintermann et al., 2002). In this case, altered kinematics of the tibial talonavicular joint and increased cartilage contact stresses lead to degenerative joint disease (Bischof et al., 2010).
Injury to the LAL not only induces ligament structural alterations but may also exert diverse impacts on the neuromuscular system. Research indicates that following LAL injury, patients may experience diminished neuromuscular recruitment, weakened muscle strength, and impaired motor coordination (Gribble et al., 2013). These consequences can manifest as challenges in daily activities, including uncoordinated gait patterns during walking and running, or even difficulties in maintaining postural stability. Damage to the neuromuscular system not only impacts the patient’s quality of life but also directly influences their motor performance and athletic capabilities (Xu et al., 2023a). Diminished neuromuscular recruitment is a common consequence following LAL injury (Doherty et al., 2014). Ligament injuries can result in reduced joint proprioception, affecting the nervous system’s perception and responsiveness, consequently diminishing the sensitivity and responsiveness of muscles to nerve signals (Kunugi et al., 2017). This can decrease muscle activation, disrupting normal motor control and movement coordination. Moreover, LAL injuries may contribute to a decline in muscle strength. Because ligament injuries induce joint instability, the adjacent muscles may undergo functional deterioration, resulting in diminished muscle strength (Lin et al., 2019; Han et al., 2022). This reduction in muscle strength can exacerbate joint instability, perpetuating a cycle of vulnerability to injury (Xu et al., 2023b).
Studies have shown that the LAL plays a crucial role in maintaining ankle stability, and injury to it may lead to impaired joint stability (Gribble et al., 2016). Patients may experience ankle laxity and lateral instability following injury, increasing the risk of re-sprain or recurrence (Xu et al., 2024). Normally, the LAL maintains stability by supporting the ankle and limiting joint motion. However, once the LAL is damaged, the joint’s stability will be compromised. Secondly, the loss of ankle stability not only affects the patient’s daily activities but may also increase the risk of joint damage (Hintermann et al., 2004). Due to the lack of joint stability, patients are more susceptible to external impact or torque during sports or activities, increasing the likelihood of re-injury (Bae et al., 2015). Ankle instability may also damage the periarticular soft tissues, exacerbating joint dysfunction and pain symptoms. Chronic ankle instability can ultimately result in cartilage degeneration and the eventual development of ankle osteoarthritis (Valderrabano et al., 2006). Long-term follow-up studies have reported that osteoarthritis occurs in 13%–78% of patients with ankle instability over 10 years (Valderrabano et al., 2009). Additionally, recent studies involving patients have identified LAL injuries as the primary cause of ankle osteoarthritis following ligamentous lesions (Löfvenberg et al., 1994).
The finite element (FE) method plays an important role in the field of biomechanics, enabling in-depth analysis of the behavior of complex joints and tissues under clinically relevant loading conditions. It not only accurately simulates real-life scenarios but also provides insights beyond traditional biomechanical studies. When analyzing ankle stability, finite element modeling emerges as a powerful tool for studying the mechanical behavior of the lateral ligament of the ankle. The material behavior of both the ankle and the lateral collateral ligament typically exhibits nonlinearity and nonuniformity, particularly under conditions involving significant deformations and high strain rates (Telfer et al., 2014). Finite element models can integrate these material properties, thereby enabling a more precise simulation of the stress-strain response of the ligaments, as well as the deformation behavior under varying loading conditions (Peng et al., 2023). By meticulously modeling the geometry and material properties of the ankle joint, finite element simulations offer researchers valuable insights into the mechanisms governing ankle stability. Additionally, they unveil the internal stress distribution and deformation patterns of the ligaments under stress conditions (Mabrouk et al., 2022).
There have been previous finite element studies on the ankle joint. Mangwani et al. utilized finite elements to analyze the impacts of ankle ligament injuries, including lateral, syndesmotic, and medial injuries (Talbott et al., 2023). They also investigated the effects of stepwise repairs for each injury on joint displacements and contact stresses, providing valuable insights into optimal repair strategies and prognosis (Halloran et al., 2023). Furthermore, finite element modeling can optimize techniques such as determining the number and placement of bone anchors in lateral ligament reconstruction, evaluating the impact of fibrous bands or hamstring reinforcement on lateral ligament reconstruction, and dynamically stabilizing screws in the joint (Shin et al., 2012). It can also quantify the effects of bone alignment, such as heel pronation and correction. Additionally, finite element analysis can serve as a treatment and reattachment strategy for insertional Achilles tendinopathy (Tits and Ruffoni, 2021).
Previous research has not explored the changes in joint stability and stress response distribution in patients with ankle ligament injuries during landing, understanding these changes is crucial for developing effective rehabilitation strategies and preventive measures for ankle ligament injuries. Therefore, this study aimed to explore the stress response and neural control changes during landing in individuals with lateral ankle ligament injuries and to evaluate the impact of varying degrees of ligament laxity on metatarsal stress. Our hypothesis posited that metatarsal stress would escalate with increased ligament laxity.
2 Method
The main framework of this study consists of the following components: 1) Collection of subject biomechanical data. 2) Processing data using Matlab. 3) Construction of finite element models. 4) Validation of finite element model results using a high-speed dual fluoroscopic imaging system (DFIS). First, kinematic and kinetic data of the participants were collected using Kistler force plates and Vicon, while sEMG signals were gathered using EMG sensors. Next, the data were preprocessed using MATLAB. Then, a finite element model was constructed, and the data were imported. Finally, the results of the finite element model were validated using DFIS.
3 Subjects
Twenty healthy subjects and twenty subjects with CLAI were recruited for this study (Table 1). All CLAI subjects were treated conservatively for 6 months before the experiment; however, they continued to exhibit symptoms of pain, instability and decreased proprioception. The diagnosis of CLAI was confirmed through a clinical examination by a foot and ankle orthopedic surgeon, and the clinical manifestations of lateral ankle injury were confirmed by magnetic resonance imaging (MRI). The criteria for subject screening were: 1) no history of foot and ankle surgeries; 2) no additional ankle pathology was identified by MRI. Reasons for exclusion included peroneal tendinopathy, and intra-articular small bone or cartilage degeneration (Caputo et al., 2009).
Prior to engaging in the study, all individuals were thoroughly informed on the study’s purpose, requirements, and procedures. This information was provided after each participant had signed the informed consent form. The study protocol received approval from the Ningbo University Scientific Research Ethics Committee (Approval Number: RAGH20240405).
3.1 Biomechanics parameters collection and processing
This study utilized a Kistler force platform and an eight-camera Vicon motion capture system (Oxford Metrics Ltd., Oxford, United Kingdom) to synchronize the kinetic and kinematic data collection. Kinematic data collection was performed at frequencies of 100 Hz, while kinetic data collection was conducted at frequencies of 200 Hz (Ward et al., 2009). Markers were positioned according to the gait 2,392 (Figure 1A), and EMG sensors were placed following the guidelines provided by SENIAM (Figure 1C) (Hermens et al., 2000). This study minimized the impedance of the skin-electrode interface by shaving the hair near the skin and cleaning the area with alcohol. Muscle activation was measured using six electromyography (EMG) sensors from Delsys (Boston, MA, United States) (Figure 2). All subjects were instructed to perform landings from two consecutive steps to collect the corresponding kinetic data. Data was successfully collected 20 times for each subject, recording only the data from stable landings (Figure 1B). The three-dimensional marker trajectories and ground reaction force data were collected using Vicon Nexus 2.14.0 and exported as C3D format files. These data were then processed in MATLAB (MathWorks, MA, United States), which involved coordinate system conversion, low-pass filtering, data extraction, and format conversion for both the kinematic and ground reaction force data. This study utilized the OpenSim software (Stanford, CA, United States) for finite element analysis to calculate biomechanical parameters as the boundary condition (Cheung et al., 2009).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438194/fphys-15-1438194-HTML/image_m/fphys-15-1438194-g001.jpg)
Figure 1. Illustrative representation depicting the placement of reflective marker points on the musculoskeletal model. (A) Visual representation of the placement of reflective markers on the constructed musculoskeletal model. (B) Depiction of the positioning of EMG electrodes on the lower limbs of human subjects. (C) Illustration outlining the procedure of the unanticipated landing biomechanics test. (D) Finite Element Boundary Condition Setting.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438194/fphys-15-1438194-HTML/image_m/fphys-15-1438194-g002.jpg)
Figure 2. Validation was conducted by comparing the patterns of measured muscle activation and simulated activation.
3.2 The process of obtaining and reconstructing geometric data
One of the subjects with chronic lateral ankle instability (CLAI) was chosen to undergo magnetic resonance imaging (MRI) and computed tomography (CT) with a 2 mm interval for this study. The resulting 2D images were then segmented using Mimics 21.0 (Materialise, Leuven, Belgium). Subsequently, the generated foot model was imported into Geomagic Studio 2021 (Geomagic, Inc., Research Triangle Park, NC, United States) to optimize the model. The imported components were then assembled into solid form using SolidWorks 2017 (SolidWorks Corporation, Waltham, MA, United States). Finally, the model’s contacts were meshed and modeled using Workbench 2021 (ANSYS, Inc. located in Canonsburg, PA, United States), and static finite element analysis was conducted (Figure 1D).
3.3 Composite representation of the LAL
In the finite element model, the lateral collateral ligament (LAL) of the ankle was initially simulated too homogeneously, and this simplistic representation could not accurately capture the laxity of the lateral ankle ligament. Setting the lateral ankle ligament to be isotropic may lead to inaccurate results (Siegler et al., 1988). Therefore, we utilized a refined hyperelastic composite to model the physiological structure of the LAL, incorporating collagen fibers for enhanced accuracy (Figure 3A). This composite LAL model comprises a proteoglycan matrix reinforced by collagen fibers, with a collagen fiber volume fraction set at 60% (Kumai et al., 2002). This value was determined based on the upper limit of the total cross-sectional area occupied by collagenous proto-fibers in the LAL.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438194/fphys-15-1438194-HTML/image_m/fphys-15-1438194-g003.jpg)
Figure 3. Finite element simulation workflows. (A) Composite representation of the anterior talofibular ligament; (B) Finite element model validation.
The strain energy formula for hyperelastic material is derived from the strain energy density function:
In this study, a Neo-Hookean model was employed. The volumetric component, denoted as
For the volumetric component
For the deviatoric component
Limited ability of collagen fibers in the ligament to withstand loads, so it is defined as (Weiss et al., 1996):
Therefore, the strain energy equation is determined as:
3.4 Boundary, loading condition and model validation
To simulate the real situation, in this study, the angle of the ankle joint was determined based on the position of the foot model. This was achieved by fixing the floor and adjusting the angle between the tibial axis and the longitudinal axis of the foot in the sagittal plane within the finite element model. The adjustment ensured alignment between the global coordinate system and the original coordinate system of OpenSim (Delp et al., 2007). The inertial forces experienced during landing were simulated by applying ankle moments and reaction forces to the talar sliding joint and tibiotalar joint surfaces, respectively. Additionally, joint forces of the MPJ were applied to the upper surface of the first metatarsal and proximal phalanx bone. All materials were assumed to be isotropic and linearly elastic, except for the encapsulated soft tissues, skin, and ligaments, for which properties were derived from previous studies (Table 2). In this study, fluoroscopic image data of the ankle joints of subjects in the grounded condition were acquired using a high-speed dual fluoroscopic imaging system (DFIS) (Ti-WISH-II, Ti-Motion Ltd., Shanghai City, CN). These computed ankle joint displacements were then compared with the results obtained from finite-element analysis to validate its accuracy (Roach et al., 2017). Specifically, we validated the accuracy of the finite element model of the foot by measuring the displacement of the navicular bone, which is commonly used in clinical practice to assess specific deformations of the foot.
First, high-resolution images of the ankle joint were obtained using DFIS imaging technology to ensure clarity and accuracy. These images were then imported into Rhinoceros software for detailed model adjustments. Next, the displacement of the navicular bone was calculated using the coordinate system calculator plugin in Rhinoceros (Li et al., 2008). This plugin accurately measures the spatial position changes of various points in the model to compute the actual displacement data of the navicular bone. Finally, these calculated results were compared with the finite element model for verification. The finite element model provides theoretical displacement data, and by comparing it with the actual calculated displacement data, the accuracy of the model can be validated (Figure 3B) (Nielsen et al., 2009; Koo and Li, 2016).
DFIS system generates high voltage through a high-voltage generator, and the X-ray tube is used to generate X-rays. The X-ray tube contains a cathode and an anode, when high voltage is applied to the X-ray tube, the cathode releases free electrons. These changes produce an image on the detector, forming an X-ray fluoroscopic image. The DFIS system is equipped with two X-ray sources and two detectors, which are positioned on each side of the object under observation. This setup allows for the simultaneous acquisition of X-ray fluoroscopic images from both directions, enabling reliable medical image data for clinical diagnosis.
4 Result
This study investigated the kinematic and kinetic changes in healthy individuals and patients with CLAI during landing. Additionally, we explored the stress distribution using a finite element model. The specific findings are detailed as follows:
4.1 Joint angle and moment
According to Figure 4, during the 0%–14% landing phase, the ankle joint angle in the sagittal plane shows significant changes in CLAI patients compared to the healthy group (p < 0.001). In the coronal plane, during the 19%–38% landing phase, the ankle joint angle changes in the healthy group are significantly greater than those in CLAI patients (p = 0.033). In the horizontal plane, during the initial 0%–3% landing phase, the ankle joint angle changes in CLAI patients are significantly greater than in the healthy group (p = 0.016). Additionally, during the 70%–100% landing phase, the ankle joint angle in CLAI patients shows significant changes compared to the healthy group (p < 0.001).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438194/fphys-15-1438194-HTML/image_m/fphys-15-1438194-g004.jpg)
Figure 4. Mean and standard deviation of the ankle and knee joints. The red line illustrates the SPM analysis findings between the CLAI and healthy groups. From top to bottom are sagittal plane, coronal plane and horizontal plane.
For the knee joint, in the sagittal plane, CLAI patients are greater than that of the healthy group. Significant changes in the knee joint angle of CLAI patients are observed during the 42%–100% landing phase (p < 0.001). In the coronal plane, during the 77%–100% landing phase, the knee joint angle changes in the healthy group are significantly greater than those in CLAI patients (p < 0.05). In the horizontal plane, during the 23%–31% landing phase, the knee joint angle changes in CLAI patients are significantly greater than in the healthy group (p < 0.001), while during the 52%–71% phase, the angle changes are significantly smaller than in the healthy group (p < 0.05) (Table 3).
The results indicate that in the sagittal plane, the peak ankle joint moment in the CLAI group is significantly smaller than in the healthy group (p < 0.01). In contrast, in the coronal plane, the peak ankle joint moment in the CLAI group is significantly higher than in the healthy group (p < 0.05). However, in the horizon plane, there was no significant difference between the two groups (p = 0.05). For the knee joint moment in the sagittal plane, the CLAI group showed significantly higher compared to the healthy group (p < 0.001). Similarly, in the coronal plane, the knee joint moment was significantly greater in the CLAI group than in the healthy group (p < 0.01). Differences were also observed in the horizon plane, with the CLAI group exhibiting higher knee joint moments than the healthy group (p < 0.05) (Table 3).
4.2 Muscle activation
Patients with chronic lateral ankle instability (CLAI) exhibit different patterns of muscle activation compared to the healthy group during landing. Based on Figure 5, the muscle activation of the vastus lateralis muscle during landing was higher in CLAI patients compared to the healthy group (p < 0.05). Additionally, the results revealed that during the initial 0%–12% phase of landing, the muscle activation of the vastus medialis muscle in CLAI patients preceded that of the healthy group (p < 0.001), with higher peak muscle activation observed as well. However, injury to the lateral ankle ligament during landing resulted in significant differences. Conversely, during the 8%–35% landing phase (p = 0.012), the healthy group exhibited significantly higher muscle activation in the soleus compared to the CLAI group. Regarding the medial gastrocnemius and lateral gastrocnemius muscles, the results demonstrated a similar disparity in muscle activation between the healthy group and the CLAI group. Muscle activation in the medial gastrocnemius exhibited an earlier peak activation in the healthy group (p < 0.001), along with a higher peak muscle activation. Similarly, activation of the lateral gastrocnemius muscle occurred earlier in the healthy group compared to the CLAI group (p < 0.001). Furthermore, significant differences were observed in the tibialis anterior muscle between the CLAI and healthy groups during the 22%–40% and 48%–61% landing phases (p < 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438194/fphys-15-1438194-HTML/image_m/fphys-15-1438194-g005.jpg)
Figure 5. Mean and standard deviation of the simulated muscle activations. The red line illustrates the SPM analysis findings between the CLAI and healthy groups.
4.3 Muscle force
Based on SPM1d analysis, according to Figure 6, significant differences were observed in the muscle force of the tibialis anterior muscle between the healthy group and CLAI patients during the 0%–36% landing phase (p < 0.05). Conversely, during the 81%–98% phase of landing, the CLAI group exhibited higher muscle force compared to the healthy group. For the soleus, the CLAI group was higher than the healthy group most of the time, notably demonstrating a significant difference during the 30%–40% phase (p < 0.05). However, it was lower than the healthy group during the 0%–19% and 93%–100% landing phases (p < 0.001, p < 0.05, respectively). The muscle force of the peroneus longus muscle in the CLAI group was comparable to that of the healthy group, except for the 32%–44% phase of the landing, where it was lower than the healthy group (p < 0.001). Meanwhile, the medial gastrocnemius muscle force was significantly higher in the healthy group than in the CLAI group during the 0%–16% landing phase (p < 0.001), while in the 42%–78% phase (p < 0.001) it was higher in the CLAI group than in the healthy group. Lateral gastrocnemius muscle force was essentially similar between the two groups, although differences were observed during the 0%–5% and 21%–32% landing phases (p = 0.034, p < 0.05, respectively).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438194/fphys-15-1438194-HTML/image_m/fphys-15-1438194-g006.jpg)
Figure 6. Mean and standard deviation of the simulated muscle forces. The red line illustrates the SPM analysis findings between the CLAI and healthy groups.
4.4 Stress distribution results
During the simulation of a CLAI patient’s landing, the thermograms generated from finite element analysis illustrate changes in stress distribution on the ligaments as different ligaments exhibit laxity (Figure 7). When the ATFL exhibits laxity, the peak von Mises stress in the CFL rises to 3.8359 MPa, while the peak von Mises stress within the ATFL itself is 1.7911 MPa. A similar situation was observed when the CFL exhibited laxity. The stress on the ATFL increased significantly to 2.3381 MPa compared to its laxity. The CFL did not experience any additional stress and only 2.1875 MPa. However, some differences emerge when the ATFL and CFL exhibit laxity. The stresses experienced by the ATFL and CFL do not increase, with the peak von Mises stress on the ATFL being 1.3351 MPa, and on the CFL being 1.9835 MPa. Figure 8 illustrates the effect of ligament laxity on the calcaneus and talus bones. Under conditions where the CFL exhibits laxity, the peak von Mises stress on the calcaneus is 7.6799 MPa. However, when the ATFL exhibits laxity, the stress on the calcaneus increases to 8.5861 MPa, and the overall stress distribution is higher than when the CFL exhibits laxity. When both the ATFL and CFL exhibit laxity, the peak von Mises stress on the calcaneus significantly increases to 9.1576 MPa. The stress distribution of the talus bone shows a similar trend. When the ATFL exhibits laxity, the peak von Mises stress on the talus bone is 4.8085 MPa, while under CFL laxity, the stress on the talus bone is slightly lower at 3.6159 MPa. When both ligaments exhibit laxity, the peak von Mises stress on the talus bone increases to 5.2019 MPa.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438194/fphys-15-1438194-HTML/image_m/fphys-15-1438194-g007.jpg)
Figure 7. Distribution of ligament von Mises stress response during ligament laxity. The results from left to right are anterior talofibular ligament laxity, calcaneofibular ligament laxity, and both ligaments laxity. (A) anterior talofibular ligament; (B) calcaneofibular ligament.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438194/fphys-15-1438194-HTML/image_m/fphys-15-1438194-g008.jpg)
Figure 8. Distribution of von Mises stress response during ligament laxity. The results from left to right are anterior talofibular ligament laxity, calcaneofibular ligament laxity, and both ligaments laxity. (A) calcaneus bone; (B) talus bone.
Figure 9 illustrates the changes in stress on different metatarsal bones when ligaments exhibit laxity. The stress variation is relatively small in the first metatarsal bone, with a peak von Mises stress of 11.9532 MPa observed when the ATFL exhibits laxity. Under CFL laxity conditions, the peak von Mises stress measures only 10.1714 MPa, slightly lower than observed with ATFL laxity. However, with simultaneous laxity in both ligaments, the peak von Mises stress reaches 13.898 MPa. However, the second metatarsal bone exhibits a different pattern. With the ATFL laxity, the peak von Mises stress measures 15.366 MPa, which increases to 17.936 MPa under CFL laxity. When laxity occurs in both ligaments simultaneously, the peak von Mises stress rises to 20.201 MPa. The third metatarsal bone demonstrates a behavior that resembles that of the second metatarsal bone, showing a peak von Mises stress of 10.183 MPa under ATFL laxity and 12.537 MPa under CFL laxity. With laxity present in both ligaments, the peak von Mises stress increases to 16.673 MPa. The variation in peak von Mises stress is less noticeable in the fourth metatarsal bone, with 13.548 MPa when the ATFL exhibits laxity, 10.506 MPa when the CFL exhibits laxity, and 15.096 MPa when both ligaments exhibit laxity. The fifth metatarsal bone exhibits the highest von Mises stress among the metatarsals. Under ATFL laxity, the peak von Mises stress measures 20.740 MPa, while under CFL laxity, it measures 17.521 MPa. Simultaneous laxity in both ligaments results in a peak von Mises stress of 21.931 MPa.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438194/fphys-15-1438194-HTML/image_m/fphys-15-1438194-g009.jpg)
Figure 9. Distribution of metatarsal von Mises stress response during ligament laxity. The results from left to right are anterior talofibular ligament laxity, calcaneofibular ligament laxity and both ligaments laxity. (A) first metatarsal; (B) second metatarsal; (C) third metatarsal; (D) fourth metatarsal; (E) fifth metatarsal.
5 Discussion
The purpose of this study was to examine alterations in stress response and neural control during landing among patients with lateral collateral ligament injuries of the ankle. We conducted a comparative analysis of metatarsal stress changes during landing in subjects exhibiting varying degrees of lateral collateral ligament laxity in the ankle, aiming to enhance ankle stability in patients. Additionally, we aimed to contribute to the development of more scientifically informed treatment protocols for the rehabilitation of such patients through further research. We hypothesized that metatarsal stress would vary with the degree of laxity of the lateral ligament of the ankle joint, thus affecting the stability of the ankle joint. The results of this study are consistent with our initial hypothesis.
The study found that compared to the healthy group, CLAI patients exhibited greater ankle joint angle changes in the sagittal plane. This difference may be due to their inability to effectively regulate ankle joint angles upon landing, leading to abnormal joint movement patterns. Additionally, there were differences in knee joint angles between the healthy group and CLAI patients, indicating distinct knee joint movement control strategies between the two groups. These findings are consistent with previous research, suggesting that CLAI patients display altered landing patterns compared to the healthy group (Brown et al., 2008). Gehring et al. (2014) research highlights significant deviations in knee joint angles during landing and cutting movements in the CLAI group compared to the healthy group, emphasizing the importance of considering unanticipated tasks when assessing actual LAL. The study revealed that the CLAI group employed different movement strategies during landing tasks than the healthy group. Further research is needed to verify the effectiveness of the observed changes in neural control strategies in reducing the risk of LAL.
This study observed that laxity of the lateral ankle ligaments increases stress on the metatarsals, leading to ankle instability. Research indicates that damage or laxity of the lateral ligaments alters the biomechanical properties of the ankle joint. Typically, the lateral ankle ligaments maintain stability by supporting and restricting joint movement (Bonnel et al., 2010; Jiang et al., 2024). Damage to these ligaments may result in increased pressure on the metatarsals, impairing their ability to effectively support ankle movement. The findings also corroborate previous research, showing that further laxity of the ATFL and CFL significantly alters metatarsal stress during simulated CLAI landings. Under CFL laxity, peak von Mises stress on the fifth metatarsal increases, and under ATFL laxity, the peak von Mises stress on the fifth metatarsal is even greater, reaching its peak when both ligaments are lax. Studies by Doherty indicate that lateral ankle ligament injuries restrict inversion and eversion movements of the ankle joint, leading to additional load on the metatarsals and increasing the risk of ankle injury (Doherty et al., 2016). Dobbe et al. observed increased metatarsal stress in patients with chronic lateral ankle instability, particularly during physical activity (Dobbe et al., 2020). Additionally, the study results demonstrate that ligament stress gradually increases upon landing with lateral ligament laxity, suggesting changes in the surrounding tissues and structures, consistent with previous findings. During this process, the mechanical properties of the ligaments may change, making them more sensitive to external stress and potentially increasing the risk of secondary injury when damaged ligaments are subjected to external stress (Barelds et al., 2018).
Furthermore, the study found that ligament laxity can induce muscle compensation. Results suggest that laxity of the ATFL may transfer stress to the CFL, and conversely, laxity of the CFL may increase stress on the ATFL. When two ligaments are lax, stress may transfer to other muscles to maintain the stability of the foot structure. This muscle compensation may involve muscles and ligaments surrounding the foot. These findings are significant for understanding muscle strength reconstruction during ligament injury rehabilitation and preventing further injury. Banerjee and Agarwal. (1998) research emphasizes the indispensable role of the talus and calcaneus in weight-bearing, providing anchorage for other crucial ligament structures, thereby tightly connecting the critical joints of the distal lower limb. Laxity of the lateral ankle ligaments may further damage the talus and calcaneus. Moreover, the study also found that with ligament laxity, stress between the calcaneus and talus changes, potentially leading to talus displacement and affecting joint stability. Therefore, maintaining the stability of the lateral ankle ligaments is crucial for protecting the talus and calcaneus and preserving the function and stability of the foot joints (Ferran et al., 2009; Kang and Jiang, 2024).
Another noteworthy discovery in the present study is the distinct muscle activation pattern observed in patients with chronic lateral ankle instability (CLAI) during landing compared to the healthy group, which corroborates earlier research. DeJong et al. (2020) demonstrated that individuals with CLAI showcased compensatory muscle activation in the proximal joint, while Rios et al. (2015) observed heightened muscle activation in CLAI patients during single-leg stance compared to healthy controls, particularly in the activation of proximal joint muscles during the standing phase. The current study noted increased muscle activation in the medial and lateral femoral muscles among CLAI patients compared to the healthy group. This indicates that individuals with CLAI tend to rely on heightened proximal muscle activation as a compensatory mechanism to improve motor control and mitigate neuromuscular deficits in the ankle joint. Previous research has similarly highlighted decreased neuromuscular recruitment as a common manifestation following LAL injury (Feger et al., 2014). The ligament injury may lead to decreased joint positional sensation, thereby reducing muscle sensitivity and responsiveness to nerve signals. The study findings indicated delayed activation of the medial gastrocnemius and lateral gastrocnemius muscles in the CLAI group compared to the healthy group. Furthermore, the soleus muscle exhibited not only earlier activation but also higher peak muscle activation. H. Kim’s analysis revealed that muscles like the peroneus longus, medial gastrocnemius, and lateral gastrocnemius were activated earlier during walking in the healthy group. Conversely, CLAI patients displayed delayed activation of the lateral gastrocnemius and soleus muscles during landing, potentially contributing to ankle instability (Kim et al., 2019). Moreover, CLAI patients exhibited greater tibialis anterior muscle strength and activation compared to the healthy group, indicating an augmented reliance on ankle dorsiflexion for stabilization during landing (Jie et al., 2024). The above study suggests that CLAI patients exhibit altered motor control patterns during landing compared to the healthy group. In the absence of adequate treatment and rehabilitation, CLAI patients may employ alternative strategies to compensate for ankle instability, potentially heightening the risk of secondary injury (Son et al., 2019).
Given the observed compensatory muscle activation patterns in patients with CLAI, targeted rehabilitation programs can be developed to address these neuromuscular deficits. For instance, strengthening exercises focused on the medial and lateral femoral muscles may help improve proximal muscle support and overall joint stability (Xu et al., 2022; Zhou and Ugbolue, 2024). Additionally, proprioceptive training can be incorporated to enhance joint positional sense and muscle responsiveness, thereby reducing the risk of further injury. By customizing rehabilitation protocols to address specific muscle activation delays and deficiencies identified in this study, doctors can optimize recovery and functional outcomes for individuals with CLAI (Watabe et al., 2021). Such strategies not only aim to restore normal muscle activation patterns but also to prevent compensatory mechanisms that could lead to secondary injuries.
This study also has several limitations. Firstly, only one patient with stabilized chronic lateral ankle instability was included in the development of the finite element model. Given inherent individual differences, the conclusions drawn from the study may vary. Additionally, while the ligaments were modeled as hyperelastic in this study, the overall rigidity of the model was only tailored to a single individual and did not fully encompass collective variations. Moreover, the boundary conditions of the models were uniform and did not simulate the actual process of landing injury. Variations in setup conditions such as material properties, mesh size, and mesh behavior could significantly influence the results. It’s important to acknowledge that this scenario may not entirely reflect real-world conditions. Meanwhile, only displacement of navicular height was used in the model validation. Subsequent experiments should include validation through plantar pressure distribution, joint contact stress, and contact area.
6 Conclusion
In this study, we observed the impact of lateral ankle ligament laxity on the biomechanical characteristics and stability of the foot. The results indicate that injury or laxity of the lateral ligaments may lead to increased stress on the metatarsals, thereby compromising the stability of the ankle joint. Specifically, further laxity of the ATFL and CFL significantly altered the distribution of stress on the metatarsals, resulting in a higher risk of ankle injury. Additionally, we observed that ligament laxity may trigger muscle compensation, further affecting the stability of the foot structure. These findings are crucial for understanding the rebuilding of muscle strength during the rehabilitation process following ligament injuries and for preventing recurrent injuries. Therefore, maintaining the stability of the lateral ankle ligaments is essential for preserving the function and stability of the foot joints.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving humans were approved by Ethical committee of Ningbo University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
ZZ: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Writing–original draft. HZ: Data curation, Investigation, Project administration, Writing–review and editing. TJ: Investigation, Methodology, Writing–original draft. DX: Data curation, Formal Analysis, Investigation, Writing–review and editing. E-CT: Data curation, Formal Analysis, Investigation, Writing–review and editing. MW: Data curation, Formal Analysis, Writing–review and editing. YG: Conceptualization, Formal Analysis, Funding acquisition, Project administration, Resources, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was funded by Zhejiang Provincial Natural Science Foundation of China for Distinguished Young Scholars (LR22A020002), Zhejiang Provincial Key Research and Development Program of China (2023C03197), Zhejiang Provincial Natural Science Foundation (LTGY23H040003), Ningbo key R&D Program (2022Z196), Zhejiang Province Exploring Public Welfare Projects (LTGY23H040003), Ningbo Natural Science Foundation (20221JCGY010532, 20221JCGY010607), Public Welfare Science & Technology Project of Ningbo, China (2021S134), and Zhejiang Rehabilitation Medical Association Scientific Research Special Fund (ZKKY2023001).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Bae J. Y., Park K. S., Seon J. K., Jeon I. (2015). Analysis of the effects of normal walking on ankle joint contact characteristics after acute inversion ankle sprain. Ann. Biomed. Eng. 43 (12), 3015–3024. doi:10.1007/s10439-015-1360-1
Banerjee K., Agarwal B. (1998). Estimation of age from epiphyseal union at the wrist and ankle joints in the capital city of India. Forensic Sci. Int. 98 (1-2), 31–39. doi:10.1016/s0379-0738(98)00134-0
Barelds I., van den Broek A. G., Huisstede B. M. (2018). Ankle bracing is effective for primary and secondary prevention of acute ankle injuries in athletes: a systematic review and meta-analyses. Sports Med. 48, 2775–2784. doi:10.1007/s40279-018-0993-2
Bischof J. E., Spritzer C. E., Caputo A. M., Easley M. E., DeOrio J. K., Nunley II J. A., et al. (2010). In vivo cartilage contact strains in patients with lateral ankle instability. J. Biomechanics 43 (13), 2561–2566. doi:10.1016/j.jbiomech.2010.05.013
Bonnel F., Toullec E., Mabit C., Tourné Y.Sofcot (2010). Chronic ankle instability: biomechanics and pathomechanics of ligaments injury and associated lesions. Orthop. Traumatology Surg. Res. 96 (4), 424–432. doi:10.1016/j.otsr.2010.04.003
Brown C., Padua D., Marshall S. W., Guskiewicz K. (2008). Individuals with mechanical ankle instability exhibit different motion patterns than those with functional ankle instability and ankle sprain copers. Clin. Biomech. 23 (6), 822–831. doi:10.1016/j.clinbiomech.2008.02.013
Caputo A. M., Lee J. Y., Spritzer C. E., Easley M. E., DeOrio J. K., Nunley J. A., et al. (2009). In vivo kinematics of the tibiotalar joint after lateral ankle instability. Am. J. Sports Med. 37 (11), 2241–2248. doi:10.1177/0363546509337578
Chen W.-M., Park J., Park S.-B., Shim V. P.-W., Lee T. (2012). Role of gastrocnemius–soleus muscle in forefoot force transmission at heel rise—a 3D finite element analysis. Am. J. Sports Med. 45 (10), 1783–1789. doi:10.1016/j.jbiomech.2012.04.024
Cheung J. T.-M., Yu J., Wong D. W.-C., Zhang M. (2009). Current methods in computer-aided engineering for footwear design. Footwear Sci. 1 (1), 31–46. doi:10.1080/19424280903002323
Cheung J. T.-M., Zhang M., Leung A. K.-L., Fan Y.-B. (2005). Three-dimensional finite element analysis of the foot during standing—a material sensitivity study. J. Biomechanics 38 (5), 1045–1054. doi:10.1016/j.jbiomech.2004.05.035
DeJong A. F., Koldenhoven R. M., Hart J. M., Hertel J. (2020). Gluteus medius dysfunction in females with chronic ankle instability is consistent at different walking speeds. Clin. Biomech. 73, 140–148. doi:10.1016/j.clinbiomech.2020.01.013
Delp S. L., Anderson F. C., Arnold A. S., Loan P., Habib A., John C. T., et al. (2007). OpenSim: open-source software to create and analyze dynamic simulations of movement. IEEE Trans. Biomed. Eng. 54 (11), 1940–1950. doi:10.1109/tbme.2007.901024
DiGiovanni B. F., Partal G., Baumhauer J. F. (2004). Acute ankle injury and chronic lateral instability in the athlete. Clin. Sports Med. 23 (1), 1–19. doi:10.1016/s0278-5919(03)00095-4
Dobbe A., Beaupre L. A., Almansoori K. A., Fung T.-S., Scharfenberger A. V. (2020). Functional outcomes of isolated infrasyndesmotic fibula fractures. Foot Ankle Orthop. 5 (1), 2473011419892227. doi:10.1177/2473011419892227
Doherty C., Bleakley C., Hertel J., Caulfield B., Ryan J., Delahunt E. (2016). Dynamic balance deficits in individuals with chronic ankle instability compared to ankle sprain copers 1 year after a first-time lateral ankle sprain injury. Knee Surg. Sports Traumatol. Arthrosc. 24 (4), 1086–1095. doi:10.1007/s00167-015-3744-z
Doherty C., Delahunt E., Caulfield B., Hertel J., Ryan J., Bleakley C. (2014). The incidence and prevalence of ankle sprain injury: a systematic review and meta-analysis of prospective epidemiological studies. Sports Med. 44 (1), 123–140. doi:10.1007/s40279-013-0102-5
Feger M. A., Donovan L., Hart J. M., Hertel J. (2014). Lower extremity muscle activation during functional exercises in patients with and without chronic ankle instability. PM&R 6 (7), 602–611. doi:10.1016/j.pmrj.2013.12.013
Ferran N. A., Maffulli N. (2006). Epidemiology of sprains of the lateral ankle ligament complex. Foot Ankle Clin. 11 (3), 659–662. doi:10.1016/j.fcl.2006.07.002
Ferran N. A., Oliva F., Maffulli N. (2009). Ankle instability. Sports Med. Arthrosc. Rev. 17 (2), 139–145. doi:10.1097/JSA.0b013e3181a3d790
Gehring D., Faschian K., Lauber B., Lohrer H., Nauck T., Gollhofer A. (2014). Mechanical instability destabilises the ankle joint directly in the ankle-sprain mechanism. Br. J. Sports Med. 48 (5), 377–382. doi:10.1136/bjsports-2013-092626
Gribble P. A., Bleakley C. M., Caulfield B. M., Docherty C. L., Fourchet F., Fong D. T.-P., et al. (2016). Evidence review for the 2016 International Ankle Consortium consensus statement on the prevalence, impact and long-term consequences of lateral ankle sprains. Br. J. Sports Med. 50 (24), 1496–1505. doi:10.1136/bjsports-2016-096189
Gribble P. A., Delahunt E., Bleakley C., Caulfield B., Docherty C., Fourchet F., et al. (2013). Selection criteria for patients with chronic ankle instability in controlled research: a position statement of the International Ankle Consortium. J. Orthop. Sports Phys. Ther. 43 (8), 585–591. doi:10.2519/jospt.2013.0303
Gu Y., Ren X., Li J., Lake M., Zhang Q., Zeng Y. (2010). Computer simulation of stress distribution in the metatarsals at different inversion landing angles using the finite element method. Int. Orthop. 34 (5), 669–676. doi:10.1007/s00264-009-0856-4
Halloran J. P., Abdollahi Nohouji N., Hafez M. A., Besier T. F., Chokhandre S. K., Elmasry S., et al. (2023). Assessment of reporting practices and reproducibility potential of a cohort of published studies in computational knee biomechanics. J. Orthop. Research® 41 (2), 325–334. doi:10.1002/jor.25358
Han S., Son S. J., Kim H., Lee H., Seeley M., Hopkins T. (2022). Prelanding movement strategies among chronic ankle instability, coper, and control subjects. Sports Biomech. 21 (4), 391–407. doi:10.1080/14763141.2021.1927163
Hermens H. J., Freriks B., Disselhorst-Klug C., Rau G. (2000). Development of recommendations for SEMG sensors and sensor placement procedures. J. Electromyogr. Kinesiol. 10 (5), 361–374. doi:10.1016/s1050-6411(00)00027-4
Hintermann B., Boss A., Schäfer D. (2002). Arthroscopic findings in patients with chronic ankle instability. Am. J. Sports Med. 30 (3), 402–409. doi:10.1177/03635465020300031601
Hintermann B., Valderrabano V., Boss A., Trouillier H. H., Dick W. (2004). Medial ankle instability: an exploratory, prospective study of fifty-two cases. Am. J. Sports Med. 32 (1), 183–190. doi:10.1177/0095399703258789
Jiang X., Sárosi J., Bíró I. (2024). Characteristics of lower limb running-related injuries in trail runners: a systematic review. Phys. Activity Health 8 (1), 137–147. doi:10.5334/paah.375
Jie T., Xu D., Zhang Z., Chon T.-E., Baker J. S., Zhou H., et al. (2024). Structural and organizational strategies of locomotor modules during landing in patients with chronic ankle instability. Bioengineering 11 (5), 518. doi:10.3390/bioengineering11050518
Kang Z. H., Jiang X. Y. (2024). The effect of running experience on muscle forces and knee joint reaction forces during running. Int. J. Biomed. Eng. Technol. 45 (3), 183–197. doi:10.1504/IJBET.2024.138969
Khawaji B., Soames R. (2015). The anterior talofibular ligament: a detailed morphological study. Foot 25 (3), 141–147. doi:10.1016/j.foot.2015.05.004
Kim H., Son S. J., Seeley M. K., Hopkins J. T. (2019). Altered movement strategies during jump landing/cutting in patients with chronic ankle instability. Scand. J. Med. Sci. Sports 29 (8), 1130–1140. doi:10.1111/sms.13445
Kitaoka H. B., Luo Z. P., Growney E. S., Berglund L. J., An K.-N. (1994). Material properties of the plantar aponeurosis. Foot Ankle Int. 15 (10), 557–560. doi:10.1177/107110079401501007
Koo T. K., Li M. Y. (2016). A guideline of selecting and reporting intraclass correlation coefficients for reliability research. J. Chiropr. Med. 15 (2), 155–163. doi:10.1016/j.jcm.2016.02.012
Kumai T., Takakura Y., Rufai A., Milz S., Benjamin M. (2002). The functional anatomy of the human anterior talofibular ligament in relation to ankle sprains. J. Anat. 200 (5), 457–465. doi:10.1046/j.1469-7580.2002.00050.x
Kunugi S., Masunari A., Yoshida N., Miyakawa S. (2017). Postural stability and lower leg muscle activity during a diagonal single-leg landing differs in male collegiate soccer players with and without functional ankle instability. J. Phys. Fit. Sports Med. 6 (4), 257–265. doi:10.7600/jpfsm.6.257
Li G., Van de Velde S. K., Bingham J. T. (2008). Validation of a non-invasive fluoroscopic imaging technique for the measurement of dynamic knee joint motion. J. Biomechanics 41 (7), 1616–1622. doi:10.1016/j.jbiomech.2008.01.034
Lin J.-Z., Lin Y.-A., Lee H.-J. (2019). Are landing biomechanics altered in elite athletes with chronic ankle instability. J. Sports Sci. Med. 18 (4), 653–662.
Löfvenberg R., Kärrholm J., Lund B. (1994). The outcome of nonoperated patients with chronic lateral instability of the ankle: a 20-year follow-up study. Foot Ankle Int. 15 (4), 165–169. doi:10.1177/107110079401500401
Mabrouk M. S., Marzouk S. Y., Afify H. M. (2022). A biomechanical analysis of prosthesis disc in lumbar spinal segment using three-dimensional finite element modelling. Int. J. Biomed. Eng. Technol. 39 (1), 1–21. doi:10.1504/ijbet.2022.10047884
Maffulli N., Ferran N. A. (2008). Management of acute and chronic ankle instability. JAAOS-Journal Am. Acad. Orthop. Surg. 16 (10), 608–615. doi:10.5435/00124635-200810000-00006
Nielsen R. G., Rathleff M. S., Simonsen O. H., Langberg H. (2009). Determination of normal values for navicular drop during walking: a new model correcting for foot length and gender. J. Foot Ankle Res. 2 (1), 12. doi:10.1186/1757-1146-2-12
Pailler-Mattei C., Bec S., Zahouani H. (2008). In vivo measurements of the elastic mechanical properties of human skin by indentation tests. Med. Eng. Phys. 30 (5), 599–606. doi:10.1016/j.medengphy.2007.06.011
Peng L., Yu L., Jia J., Gan Y., Ding A., Xiong P., et al. (2023). The effect of thickness and elastic modulus of the anterior talofibular ligament on anterior ankle joint stiffness: a subject-specific finite element study. Front. Bioeng. Biotechnol. 11 (2023), 1175347. doi:10.3389/fbioe.2023.1175347
Rios J. L., Gorges A. L., dos Santos M. J. (2015). Individuals with chronic ankle instability compensate for their ankle deficits using proximal musculature to maintain reduced postural sway while kicking a ball. Hum. Mov. Sci. 43, 33–44. doi:10.1016/j.humov.2015.07.001
Roach K. E., Foreman K. B., Barg A., Saltzman C. L., Anderson A. E. (2017). Application of high-speed dual fluoroscopy to study in vivo tibiotalar and subtalar kinematics in patients with chronic ankle instability and asymptomatic control subjects during dynamic activities. Foot Ankle Int. 38 (11), 1236–1248. doi:10.1177/1071100717723128
Shin J., Yue N., Untaroiu C. D. (2012). A finite element model of the foot and ankle for automotive impact applications. Ann. Biomed. Eng. 40 (12), 2519–2531. doi:10.1007/s10439-012-0607-3
Siegler S., Block J., Schneck C. D. (1988). The mechanical characteristics of the collateral ligaments of the human ankle joint. Foot Ankle 8 (5), 234–242. doi:10.1177/107110078800800502
Son S. J., Kim H., Seeley M. K., Hopkins J. T. (2019). Altered walking neuromechanics in patients with chronic ankle instability. J. Athl. Train. 54, 684–697. doi:10.4085/1062-6050-478-17
Swanson B., Hewitt M. A., Buckley S. E., Collins C., Hunt K. J. (2022). Epidemiology of ankle ligament injuries in US high school athletes. Foot Ankle Orthop. 7 (4), 2473011421S0096. doi:10.1177/2473011421S00961
Talbott H., Jha S., Gulati A., Brockett C., Mangwani J., Pegg E. C. (2023). Clinically useful finite element models of the natural ankle–A review. Clin. Biomech. 106, 106006. doi:10.1016/j.clinbiomech.2023.106006
Telfer S., Erdemir A., Woodburn J., Cavanagh P. R. (2014). What has finite element analysis taught us about diabetic foot disease and its management? A systematic review. PLoS ONE 9 (10), e109994. doi:10.1371/journal.pone.0109994
Tits A., Ruffoni D. (2021). Joining soft tissues to bone: insights from modeling and simulations. Bone Rep. 14, 100742. doi:10.1016/j.bonr.2020.100742
Valderrabano V., Hintermann B., Horisberger M., Fung T. S. (2006). Ligamentous posttraumatic ankle osteoarthritis. Am. J. Sports Med. 34 (4), 612–620. doi:10.1177/0363546505281813
Valderrabano V., Horisberger M., Russell I., Dougall H., Hintermann B. (2009). Etiology of ankle osteoarthritis. Clin. Orthop. Relat. Research® 467 (7), 1800–1806. doi:10.1007/s11999-008-0543-6
Ward S. R., Eng C. M., Smallwood L. H., Lieber R. L. (2009). Are current measurements of lower extremity muscle architecture accurate? Clin. Orthop. Relat. Res. 467 (4), 1074–1082. doi:10.1007/s11999-008-0594-8
Watabe T., Takabayashi T., Tokunaga Y., Yoshida T., Kubo M. (2021). Copers adopt an altered movement pattern compared to individuals with chronic ankle instability and control groups in unexpected single-leg landing and cutting task. J. Electromyogr. Kinesiol. 57, 102529. doi:10.1016/j.jelekin.2021.102529
Weiss J. A., Gardiner J. C., Bonifasi-Lista C. (2002). Ligament material behavior is nonlinear, viscoelastic and rate-independent under shear loading. J. Biomechanics 35 (7), 943–950. doi:10.1016/s0021-9290(02)00041-6
Weiss J. A., Maker B. N., Govindjee S. (1996). Finite element implementation of incompressible, transversely isotropic hyperelasticity. Comput. Methods Appl. Mech. Eng. 135 (1-2), 107–128. doi:10.1016/0045-7825(96)01035-3
Xiang L., Mei Q., Wang A., Shim V., Fernandez J., Gu Y. (2022). Evaluating function in the hallux valgus foot following a 12-week minimalist footwear intervention: a pilot computational analysis. J. Biomechanics 132, 110941. doi:10.1016/j.jbiomech.2022.110941
Xu D., Zhou H., Quan W., Gusztav F., Baker J. S., Gu Y. (2023a). Adaptive neuro-fuzzy inference system model driven by the non-negative matrix factorization-extracted muscle synergy patterns to estimate lower limb joint movements. Comput. Methods Programs Biomed. 242, 107848. doi:10.1016/j.cmpb.2023.107848
Xu D., Zhou H., Quan W., Gusztav F., Wang M., Baker J. S., et al. (2023b). Accurately and effectively predict the ACL force: utilizing biomechanical landing pattern before and after-fatigue. Comput. Methods Programs Biomed. 241, 107761. doi:10.1016/j.cmpb.2023.107761
Xu D., Zhou H., Quan W., Ma X., Chon T.-E., Fernandez J., et al. (2024). New insights optimize landing strategies to reduce lower limb injury risk. Cyborg Bionic Syst. 5, 0126. doi:10.34133/cbsystems.0126
Xu Y., Song B., Ming A., Zhang C., Ni G. (2022). Chronic ankle instability modifies proximal lower extremity biomechanics during sports maneuvers that may increase the risk of ACL injury: a systematic review. Front. Physiology 13, 1036267. doi:10.3389/fphys.2022.1036267
Keywords: lateral ankle ligament injury, finite element analysis, foot, ligament mechanics, composite material
Citation: Zhou Z, Zhou H, Jie T, Xu D, Teo E-C, Wang M and Gu Y (2024) Analysis of stress response distribution in patients with lateral ankle ligament injuries: a study of neural control strategies utilizing predictive computing models. Front. Physiol. 15:1438194. doi: 10.3389/fphys.2024.1438194
Received: 25 May 2024; Accepted: 10 July 2024;
Published: 24 July 2024.
Edited by:
Qipeng Song, Shandong Sport University, ChinaReviewed by:
Qiuxia Zhang, Soochow University, ChinaDusan Radivoje Mitic, University of Belgrade, Serbia
Copyright © 2024 Zhou, Zhou, Jie, Xu, Teo, Wang and Gu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Huiyu Zhou, zhouhuiyu@aliyun.com; Yaodong Gu, guyaodong@nbu.edu.cn