Computer simulation of stress distribution in the metatarsals at different inversion landing angles using the finite element method

Contribution of ankle motion pattern during landing to reduce the knee-related injury risk

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.

Analysis of stress response distribution in patients with lateral ankle ligament injuries: a study of neural control strategies utilizing predictive computing models

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.

Extended finite element analysis for the 2nd and 3rd metatarsals stress fracture during push-off

Metatarsal stress fractures (MSF), particularly the 2nd and 3rd MSF, are common injuries among athletes. Although there are several practices to reduce foot and ankle injuries, there is no injury prevention program specifically designed to minimize MSF. This is mainly due to the lack of information about the loadings/postures that cause MSF. Therefore, this study aimed to investigate dangerous loadings/postures potentially causing MSF during push-off (PO). The analysis was conducted with Finite Element Modelling (FEM), calibrated with the three-point bending test, and validated with peak plantar pressure (PPP) and fracture force measurement. Extended Finite Element Method was used for MSF simulation such that ten different foot and ankle configurations were designed, with five for each of the 2nd and 3rd MSF under pure vertical loadings. A more complex loading, ankle eversion/inversion during PO, was also examined for the MSF. The average error percentage for the calibration of the model with the three-point bending test was 3.05%. The average error percentages for the validation of the model with PPP and fracture force measurements were 18% and 30%, respectively. The outcomes of pure vertical loadings indicated the higher potential for the 2nd and 3rd MSF at 30% PO and 70% PO, respectively. The results of ankle eversion/inversion loadings represented that the most dangerous posture for MSF was 30° ankle eversion for the 3rd metatarsal at 70% PO. These results provide a guide, including what postures to avoid for the 2nd and 3rd MSF among people who are at high risk of MSF.

Effect of forefoot transverse arch stiffness on foot biomechanical response--based on finite element method

Background

The plantar vault, comprising the transverse and longitudinal arches of the human foot, is essential for impact absorption, elastic energy storage, and propulsion. Recent research underscores the importance of the transverse arch, contributing over 40% to midfoot stiffness. This study aimed to quantify biomechanical responses in the ankle-foot complex by varying the stiffness of the deep metatarsal transverse ligament (DTML).

Methods

Using CT image reconstruction, we constructed a complex three-dimensional finite element model of the foot and ankle joint complex, accounting for geometric complexity and nonlinear characteristics. The focus of our study was to evaluate the effect of different forefoot transverse arch stiffness, that is, different Young's modulus values of DTML (from 135 MPa to 405 MPa), on different biomechanical aspects of the foot and ankle complex. Notably, we analyzed their effects on plantar pressure distribution, metatarsal stress patterns, navicular subsidence, and plantar fascial strain.

Results

Increasing the stiffness of the DTML has significant effects on foot biomechanics. Specifically, higher DTML stiffness leads to elevate von Mises stress in the 1st, 2nd, and 3rd metatarsals, while concurrently reducing plantar pressure by 14.2% when the Young's modulus is doubled. This stiffening also impedes navicular bone subsidence and foot lengthening. Notably, a 100% increase in the Young's modulus of DTML results in a 54.1% decrease in scaphoid subsidence and a 2.5% decrease in foot lengthening, which collectively contribute to a 33.1% enhancement in foot longitudinal stiffness. Additionally, doubling the Young's modulus of DTML can reduce the strain stretch of the plantar fascia by 38.5%.

Conclusion

Preserving DTML integrity sustains the transverse arch, enhancing foot longitudinal stiffness and elastic responsiveness. These findings have implications for treating arch dysfunction and provide insights for shoe developers seeking to enhance propulsion.

Fifth Metatarsal Fractures

Fifth metatarsal features are the most common fractures in the foot. They have a long history that has resulted in many classification systems and little consensus on appropriate treatment. Although there is some agreement among experts, there are also many questions yet to be answered. There is a general consensus that dancer's fractures and zone 1 fractures can generally be treated nonoperatively. There is much more debate about zone 2 and 3 fractures and appropriate treatment guidelines. The authors review the current literature and give the recommendation for treatment based on their experience in a community-based private practice.

New Insights Optimize Landing Strategies to Reduce Lower Limb Injury Risk

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 rate-dependent) 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.

A combined musculoskeletal and finite element model of a foot to predict plantar pressure distribution

In this study, a combined subject-specific numerical and experimental investigation was conducted to explore the plantar pressure of an individual. The research utilized finite element (FE) and musculoskeletal modelling based on computed tomography (CT) images of an ankle-foot complex and three-dimensional gait measurements. Muscle forces were estimated using an individualized multi-body musculoskeletal model in five gait phases. The results of the FE model and gait measurements for the same subject revealed the highest stress concentration of 0.48 MPa in the forefoot, which aligns with previously-reported clinical observations. Additionally, the study found that the encapsulated soft tissue FE model with hyper-elastic properties exhibited higher stresses compared to the model with linear-elastic properties, with maximum ratios of 1.16 and 1.88 MPa in the contact pressure and von-Mises stress, respectively. Furthermore, the numerical simulation demonstrated that the use of an individualized insole caused a reduction of 8.3% in the maximum contact plantar pressure and 14.7% in the maximum von-Mises stress in the encapsulated soft tissue. Overall, the developed model in this investigation holds potential for facilitating further studies on foot pathologies and the improvement of rehabilitation techniques in clinical settings.

Can the Entire Function of the Foot Be Concentrated in the Forefoot Area during the Running Stance Phase? A Finite Element Study of Different Shoe Soles

The goal of this study was to use the finite element (FE) method to compare and study the differences between bionic shoes (BS) and normal shoes (NS) forefoot strike patterns when running. In addition, we separated the forefoot area when forefoot running as a way to create a small and independent area of instability. An adult male of Chinese descent was recruited for this investigation (age: 26 years old; body height: 185 cm; body mass: 82 kg) (forefoot strike patterns). We analyzed forefoot running under two different conditions through FE analysis, and used bone stress distribution feature classification and recognition for further analysis. The metatarsal stress values in forefoot strike patterns with BS were less than with NS. Additionally, the bone stress classification of features and the recognition accuracy rate of metatarsal (MT) 2, MT3 and MT5 were higher than other foot bones in the first 5%, 10%, 20% and 50% of nodes. BS forefoot running helped reduce the probability of occurrence of metatarsal stress fractures. In addition, the findings further revealed that BS may have important implications for the prevention of hallux valgus, which may be more effective in adolescent children. Finally, this study presents a post-processing method for FE results, which is of great significance for further understanding and exploration of FE results.

Effects of different contact angles during forefoot running on the stresses of the foot bones: a finite element simulation study

Introduction: The purpose of this study was to compare the changes in foot at different sole-ground contact angles during forefoot running. This study tried to help forefoot runners better control and improve their technical movements by comparing different sole-ground contact angles. Methods: A male participant of Chinese ethnicity was enlisted for the present study, with a recorded age of 25 years, a height of 183 cm, and a body weight of 80 kg. This study focused on forefoot strike patterns through FE analysis. Results: It can be seen that the peak von Mises stress of M1-5 (Metatarsal) of a (Contact angle: 9.54) is greater than that of b (Contact angle: 7.58) and c (Contact angle: 5.62) in the three cases. On the contrary, the peak von Mises stress of MC (Medial Cuneiform), IC (Intermediate Cuneiform), LC (Lateral Cuneiform), C (Cuboid), N (Navicular), T (Tarsal) in three different cases is opposite, and the peak von Mises stress of c is greater than that of a and b. The peak von Mises stress of b is between a and c. Conclusion: This study found that a reduced sole-ground contact angle may reduce metatarsal stress fractures. Further, a small sole-ground contact angle may not increase ankle joint injury risk during forefoot running. Hence, given the specialized nature of the running shoes designed for forefoot runners, it is plausible that this study may offer novel insights to guide their athletic pursuits.

Accurately and effectively predict the ACL force: Utilizing biomechanical landing pattern before and after-fatigue

BACKGROUND AND OBJECTIVE

As a fundamental exercise technique, landing can commonly be associated with anterior cruciate ligament (ACL) injury, especially during after-fatigue single-leg landing (SL). Presently, the inability to accurately detect ACL loading makes it difficult to recognize the risk degree of ACL injury, which reduces the effectiveness of injury prevention and sports monitoring. Increased risk of ACL injury during after-fatigue SL may be related to changes in ankle motion patterns. Therefore, this study aims to develop a highly accurate and easily implemented ACL force prediction model by combining deep learning and the explored relationship between ACL force and ankle motion pattern.

METHODS

First, 56 subjects' during before and after-fatigue SL data were collected to explore the relationship between the ankle initial contact angle (AIC), ankle range of motion (AROM) and peak ACL force (PAF). Then, the musculoskeletal model was developed to simulate and calculate the ACL force. Finally, the ACL force prediction model was constructed by combining the explored relationship and sparrow search algorithm (SSA) to optimize the extreme learning machine (ELM) and long short-term memory (LSTM).

RESULTS

There was almost a stronger linear relationship between the PAF and AIC (R = -0.70), AROM (R2 = -0.61). By substituting AIC and AROM as independent variables in the SSA-ELM prediction model, the model shows excellent prediction performance because of very strong correlation (R2 = 0.9992,  MSE = 0.0023,  RMSE = 0.0474). Based on the equal scaling by combining results of SSA-ELM and SSA-LSTM, the prediction model achieves excellent performance in ACL force prediction of the overall waveform (R2 = 0.9947,  MSE = 0.0076,  RMSE = 0.0873).

CONCLUSION

By increasing the AIC and AROM during SL, the lower limb joint energy dissipation can be increased and the PAF reduced, thus reducing the impact loads on the lower limb joints and reducing ACL injuries. The proposed ACL dynamic load force prediction model has low input variable demands (sagittal joint angles), excellent generalization capabilities and superior performance in terms of high accuracy. In the future, we plan to use it as an accurate ACL injury risk assessment tool to promote and apply it to a wider range of sports training and injury monitoring.

Effects of plantar fascia stiffness on the internal mechanics of idiopathic pes cavus by finite element analysis: implications for metatarsalgia

Metatarsalgia occurring in individuals with pes cavus is typically associated with abnormal loading patterns in the forefoot resulting from structural alterations. Simultaneously, the frequent overstress of the plantar fascia (PF) caused by the persistence of this foot deformity may further exacerbate the chronic pain induced by metatarsal overload. We aimed to investigate and quantify the effects of PF stiffness on the internal biomechanics of pes cavus using a computational modelling approach. A patient-specific finite element model of the foot-ankle complex using the actual three-dimensional geometry of idiopathic pes cavus bones and soft tissues was reconstructed. A sensitivity study was conducted to evaluate the effects of varying elastic modulus (0-700 MPa) of the PF on the metatarsal stress distribution, and force transmission through the metatarsophalangeal (MTP) and tarsometatarsal (TMT) joints in the pes cavus. The results indicated that variations in PF stiffness led to stress redistribution in the metatarsal region. Peak stress gradually reduced with decreasing stiffness until the PF was released, eventually resulting in a reduction of 22.39% compared to the reference value of 350 MPa. Furthermore, adjusting the PF stiffness to twice the reference value (700 MPa) increased the contact forces through the TMT and MTP joints by up to 23% and 116%, respectively. The reduction of PF stiffness alleviated focal metatarsal loading, and therefore, surgical fascia release can be considered to alleviate metatarsalgia in patients with pes cavus.

Predictive factors for the bone union disorder of intramedullary screw fixation in proximal fifth metatarsal bone fractures

INTRODUCTION

Although intramedullary screw fixation likely leads to successful union of Jones fractures compared to that of nonoperative treatments, bony union disorder after surgical treatment remains to be elucidated.

METHODS

Intramedullary screw fixation was performed for the surgical treatment of proximal fifth metatarsal stress fractures in this series. Between January 2008 and December 2019, the feet of 222 patients were investigated regarding the effective factors for postoperative bony union between the normal union group and the bony union disorder group according to the patients' physical status, radiological assessment, and screw size. The mean postoperative follow-up period was 11.1 months. Bone union disorder was defined as delayed union, nonunion, or a re-fracture recognized through a radiographic image.

RESULTS

The prevalence rate of union disorders occurred in 14% (31/222) of the patients. The risk of bone union disorder significantly increased when using a small-diameter screw (odds ratio 4.81, 95% confidence interval [CI] 1.62-14.2, p = 0.004) and non-bone graft procedures (odds ratio 3.13, 95% CI 1.22-8.02, p = 0.02). Screw length, preoperative Torg's classification, or patients' physical status did not affect postoperative bony union.

CONCLUSIONS

Approximately 14.0% of the patients in our study had postoperative bone union disorder. Small-diameter screws and non-bone graft procedures increased the risk of bone union disorder in the intramedullary screw fixation technique of fifth metatarsal bone stress fractures.

LEVEL OF EVIDENCE

Level 4, case series.

An investigation into the hammer toe effects on the lower extremity mechanics and plantar fascia tension: A case for a vicious cycle and progressive damage

Hammer toes are one of the common deformities of the forefoot that can lead to compensatory changes during walking in individuals with this condition. Predicting the adverse effects of tissue damage on the performance of other limbs is very important in the prevention of progressive damage. Finite element (FE) and musculoskeletal modeling can be helpful by allowing such effects to be studied in a way where the internal stresses in the tissue could be investigated. Hence, this study aims to investigate the effects of the hammer toe deformity on the lower extremity, especially on the plantar fascia functions. To compare the joint reactions of the hammer toe foot (HTF) and healthy foot (HF), two musculoskeletal models (MSM) of the feet of a healthy individual and that of a participant with hammer toe foot were developed based on gait analysis. A previously validated 3D finite element model which was constructed using Magnetic Resonance Imaging (MRI) of the diabetic participant with the hammer toe deformity was processed at five different events during the stance phase of gait. It was found that the hammer toe deformity makes dorsiflexion of the toes and the windlass mechanism less effective during walking. Specifically, the FE analysis results showed that plantar fascia (PF) in HTF compared to HF played a less dominant role in load bearing with both medial and lateral parts of PF loaded. Also, the results indicated that the stored elastic energy in PF was less in HTF than the HF, which can indicate a higher metabolic cost during walking. Internal stress distribution shows that the majority of ground reaction forces are transmitted through the lateral metatarsals in hammer toe foot, and the probability of fifth metatarsal fracture and also progressive deformity was subsequently increased. The MSM results showed that the joint reaction forces and moments in the hammer toe foot have deviated from normal, where the metatarsophalangeal joint reactions in the hammer toe were less than the values in the healthy foot. This can indicate a vicious cycle of foot deformity, leading to changes in body weight force transmission line, and deviation of joint reactions and plantar fascia function from normal. These in turn lead to increased internal stress concentration, which in turn lead to further foot deformities. This vicious cycle cause progressive damage and can lead to an increase in the risk of ulceration in the diabetic foot.

Development and Validation of a Subject-Specific Coupled Model for Foot and Sports Shoe Complex: A Pilot Computational Study

Nowadays, footwear serves an essential role in improving athletic performance and decreasing the risk of unexpected injuries in sports games. Finite element (FE) modeling is a powerful tool to reveal the biomechanical interactions between foot and footwear, and establishing a coupled foot-shoe model is the prerequisite. The purpose of this pilot study was to develop and validate a 3D FE coupled model of the foot and sports shoe complex during balanced standing. All major foot and shoe structures were constructed based on the participant’s medical CT images, and 3D gait analysis was conducted to define the loading and boundary conditions. Sensitivity analysis was applied to determine the optimum material property for shoe sole. Both the plantar and shoe sole areas were further divided into four regions for model validation, and the Bland–Altman method was used for consistency analysis between methods. The simulated peak plantar and sole pressure distribution showed good consistency with experimental pressure data, and the prediction errors were all less than 10% during balanced standing with only two exceptions (medial and lateral forefoot regions). Meanwhile, the Bland–Altman analysis demonstrated a good agreement between the two approaches. The sensitivity analysis suggested that shoe sole with Young’s modulus of 2.739 MPa presented the greatest consistency with the measured data in our scenario. The established model could be used for investing the complex biomechanical interactions between the foot and sports shoe and optimizing footwear design, after it has been fully validated in the subsequent works under different conditions.

A new method proposed to explore the feline's paw bones of contributing most to landing pattern recognition when landed under different constraints

Felines are generally acknowledged to have natural athletic ability, especially in jumping and landing. The adage “felines have nine lives” seems applicable when we consider its ability to land safely from heights. Traditional post-processing of finite element analysis (FEA) is usually based on stress distribution trend and maximum stress values, which is often related to the smoothness and morphological characteristics of the finite element model and cannot be used to comprehensively and deeply explore the mechanical mechanism of the bone. Machine learning methods that focus on feature pattern variable analysis have been gradually applied in the field of biomechanics. Therefore, this study investigated the cat forelimb biomechanical characteristics when landing from different heights using FEA and feature engineering techniques for post-processing of FEA. The results suggested that the stress distribution feature of the second, fourth metacarpal, the second, third proximal phalanx are the features that contribute most to landing pattern recognition when cats landed under different constraints. With increments in landing altitude, the variations in landing pattern differences may be a response of the cat's forelimb by adjusting the musculoskeletal structure to reduce the risk of injury with a more optimal landing strategy. The combination of feature engineering techniques can effectively identify the bone's features that contribute most to pattern recognition under different constraints, which is conducive to the grasp of the optimal feature that can reveal intrinsic properties in the field of biomechanics.

The contribution of the ligaments in progressive collapsing foot deformity: A comprehensive computational study

The contribution of each of the ligaments in preventing the arch loss, hindfoot valgus, and forefoot abduction seen in progressive collapsing foot deformity (PCFD) has not been well characterized. An improved understanding of the individual ligament contributions to the deformity would aid in selecting among available treatments, optimizing current surgical techniques, and developing new ones. In this study, we evaluated the contribution of each ligament to the maintenance of foot alignment using a finite element model of the foot reconstructed from computed tomography scan images. The collapsed foot was modeled by simulating the failure of all the ligaments involved in PCFD. The ligaments were removed one at a time to determine the impact of each ligament on foot alignment, and then restored one at a time to simulate isolated reconstruction. Our findings show that the failure of any one ligament did not immediately lead to deformity, but that combined failure of only a few (the plantar fascia, long plantar, short plantar, deltoid, and spring ligaments) could lead to significant deformity. The plantar fascia, deltoid, and spring ligaments were primarily responsible for the prevention of arch collapse, hindfoot valgus, and forefoot abduction, respectively. Moreover, to produce deformity, a considerable amount of attenuation in the spring, tibiocalcaneal, interosseous talocalcaneal, plantar naviculocuneiform, and first plantar tarsometatarsal ligaments, but only a small amount in the plantar fascia, long plantar, and short plantar ligaments was needed. The results of this study suggest that the ability of a ligament to prevent deformity may not correlate with its attenuation in a collapsed foot.

Orthosis and Foot Structure Affect the Fifth Metatarsal Principal Strains During Simulated Level Walking

BACKGROUND

Fractures of the proximal fifth metatarsal bone are common injuries in elite athletes and are associated with high rates of delayed union and nonunion. Structural features of the foot may increase fracture risk in some individuals, emphasizing the need for intervention strategies to prevent fracture. Although orthotic devices have shown promise in reducing fractures of the fifth metatarsal bone, the effect of orthosis on fifth metatarsal strains is not well understood.

PURPOSE

To quantify the effects of different foot orthotic constructs on principal tensile strains in the proximal fifth metatarsal bone during cadaveric simulations of level walking. An additional purpose was to investigate the relationships between structural features of the foot and corresponding strains on the fifth metatarsal bone during level walking.

STUDY DESIGN

Controlled laboratory study.

METHODS

A total of 10 midtibial cadaveric specimens were attached to a 6 degrees of freedom robotic gait simulator. Strain gauges were placed at the metaphyseal-diaphyseal junction (zone II) and the proximal diaphysis (zone III) during level walking simulations using 11 different foot orthotic configurations. Images of each specimen were used to measure structural features of the foot in an axially loaded position. The peak tensile strains were measured and reported relative to the sneaker-only condition for each orthotic condition and orthotic-specific association between structural features and principal strains of both zones.

RESULTS

In total, 2 of the 11 orthotic conditions significantly reduced strain relative to the sneaker-only condition in zone II. Further, 6 orthotic conditions significantly reduced strain relative to the sneaker-only condition in zone III. Increased zone II principal strain incurred during level walking in the sneaker-only condition showed a significant association with increases in the Meary's angle. Changes in zone III principal strain relative to the sneaker-only condition were significantly associated with increases in the Meary's angle and fourth-fifth intermetatarsal angle.

CONCLUSION

The use of orthotic devices reduced principal strain relative to the condition of a sneaker without any orthosis in zone II and zone III. The ability to reduce strain relative to the sneaker-only condition in zone III was indicated by increasing values of the Meary's angle and levels of the fourth-fifth intermetatarsal angle.

CLINICAL RELEVANCE

Clinicians can use characteristics of foot structure to determine the proper foot orthosis to potentially reduce stress fracture risk in high-risk individuals.

Workflow assessing the effect of Achilles tendon rupture on gait function and metatarsal stress: Combined musculoskeletal modeling and finite element analysis

Achilles tendon rupture (ATR) incidence has increased among badminton players in recent years. The foot internal stress was hard to obtain through experimental testing. The purpose of the current research is to develop a methodology that could improve the finite element model derived foot internal stress prediction for ATR clinical and rehabilitation applications. A subject-specific musculoskeletal model was combined with a 3D finite element model to predict the metatarsal stress. The 80% point during the push-off phase of walking was selected for the comparing between injured and uninjured sides. The surgical repaired Achilles tendon (AT) after 12 months was elongated by 5.5% than the uninjured tendon. At 80% point of stance phase, the ankle plantarflexion angle and AT force decreased by 39.6% and 21.9% on the injured side, respectively. The foot inversion degree increased by 22.9% and was accompanied by the redistribution of metatarsals von Mises stress. The stresses on the fourth and fifth metatarsals were increased by 59.5% and 85.9% on the injured side. The workflow is available to assess musculoskeletal disorders and obtain foot internal stress after ATR. The decreased ankle plantar flexor force may be affected by triceps surae muscle atrophy and weakened force transmission ability of elongated AT. The increased von Mises stress on fourth and fifth metatarsals accompanied by higher foot inversion may increase the ankle lateral sprain injury risk.

Asymptomatic foot and ankle structural injuries: a 3D imaging and finite element analysis of elite fencers

BACKGROUND

Fencing is a highly asymmetrical combat sport, that imposes high mechanical demands over repeated exposures on the musculoskeletal structures, a primary cause of injuries in fencers. However, there are limited epidemiological studies on the structural injuries of the foot and ankle in fencers. This study aimed to investigate foot and ankle structural injuries, and explore how metatarsophalangeal joint structural changes may affect the mechanisms of foot and ankle injuries in asymptomatic fencers.

METHODS

3D images of foot and ankle morphology using computed tomography were obtained from ten elite fencers. We then constructed finite element models of the first metatarsophalangeal joint in the foot of their trail legs. The validated models were used to simulate stress distribution changes from different ankle joint angles during lunging.

RESULTS

The findings showed that stress distribution changes at the medial and lateral sesamoid may have caused sesamoid fractures, and that habitual and concentrated stress on the metatarsal bones might have flattened the sesamoid groove. This process may damage the integrity of the first metatarsophalangeal joint, and consequently affect the efficiency of the windlass mechanism in fencers. During lunging, different ankle joint angles of the trail foot increased the total stress difference of the medial and lateral foot, and thus influenced the lunging quality and its stability.

CONCLUSIONS

Our findings revealed that the asymmetric nature of fencing might have caused asymptomatic foot and ankle structural injuries, and finite element analysis results indicated that this might increase the incidence of the serious injuries if unattended. Regular computed tomography examination should be introduced to monitor elite fencers' lower limb alterations, permitting unique angle adjustments in the trail foot without sacrificing technical or physiologic properties based on the exam results and reduce the lower limb injury risk.

Comparative Functional Morphology of Human and Chimpanzee Feet Based on Three-Dimensional Finite Element Analysis

To comparatively investigate the morphological adaptation of the human foot for achieving robust and efficient bipedal locomotion, we develop three-dimensional finite element models of the human and chimpanzee feet. Foot bones and the outer surface of the foot are extracted from computer tomography images and meshed with tetrahedral elements. The ligaments and plantar fascia are represented by tension-only spring elements. The contacts between the bones and between the foot and ground are solved using frictionless and Coulomb friction contact algorithms, respectively. Physiologically realistic loading conditions of the feet during quiet bipedal standing are simulated. Our results indicate that the center of pressure (COP) is located more anteriorly in the human foot than in the chimpanzee foot, indicating a larger stability margin in bipedal posture in humans. Furthermore, the vertical free moment generated by the coupling motion of the calcaneus and tibia during axial loading is larger in the human foot, which can facilitate the compensation of the net yaw moment of the body around the COP during bipedal locomotion. Furthermore, the human foot can store elastic energy more effectively during axial loading for the effective generation of propulsive force in the late stance phase. This computational framework for a comparative investigation of the causal relationship among the morphology, kinematics, and kinetics of the foot may provide a better understanding regarding the functional significance of the morphological features of the human foot.

Effect of Displacement Degree of Distal Chevron Osteotomy on Metatarsal Stress: A Finite Element Method

BACKGROUND

The stress of foot bone can effectively evaluate the functional damage caused by foot deformity and the results of operation. In this study, the finite element method was used to investigate the degree of displacement of distal chevron osteotomy on metatarsal stress and metatarsophalangeal joint load; Methods: Four finite element models of displacement were established by using the CT images of a patient with moderate hallux valgus (hallux valgus angle and intermetatarsal angle were 26.74° and 14.09°, respectively), and the validity of the model was verified. Each finite element model consisted of bones and various cartilage structures, ligaments, and plantar fascia, as well as encapsulated soft tissue. Except for soft tissue, the material properties of other parts were isotropic linear elastic material, and the encapsulated soft tissue was set as nonlinear hyperelastic material. The mesh was tetrahedral mesh. Link elements were used in ligament and plantar fascia. A ground reaction force with a half-body weight was applied at the bottom of the floor to simulate the ground reaction when standing. The upper surfaces of the encapsulated soft tissue, distal tibia, and distal fibula were fixed. The stress distribution of metatarsals and the stress of cartilage of the first metatarsophalangeal joint were compared and analyzed; Results: Compared with the hallux valgus without osteotomy, the stress of the first metatarsals and second metatarsals of 2-4 mm decreased, and the stress of the interarticular cartilage of the first metatarsophalangeal joint with 4 mm was reduced. In the case of 6 mm, the stress value between the first metatarsal and the first metatarsophalangeal joint increased, and 4 mm was the most suitable distance; Conclusions: Compared with the hallux valgus without osteotomy, the stress of the first metatarsals and second metatarsals of 2-4 mm decreased, and the stress of the interarticular cartilage of the first metatarsophalangeal joint with 4 mm was reduced. In the case of 6 mm, the stress value between the first metatarsal and the first metatarsophalangeal joint increased, and 4 mm was the most suitable distance. For the degree of displacement of the distal chevron osteotomy, the postoperative stability and the stress distribution of metatarsal bone should be considered. Factors such as hallux valgus angle, intermetatarsal angle, patient's age, body weight, and metatarsal width should be considered comprehensively. The factors affecting osteotomy need to be further explored. The degree of displacement of osteotomy can be evaluated by FE method before the operation, and the most suitable distance can be obtained.

The Influence of Heel Height on Strain Variation of Plantar Fascia During High Heel Shoes Walking-Combined Musculoskeletal Modeling and Finite Element Analysis

The therapeutic benefit of high heel shoes (HHS) for plantar fasciitis treatment is controversial. It has been suggested that plantar fascia strain can be decreased by heel elevation of shoes which helps in body weight redistribution throughout the length of the foot. Yet it is a fact that the repetitive tension caused by HHS wearing resulting in plantar fasciitis is a high-risk disease in HHS individuals who suffer heel and plantar pain. To explore the biomechanical function on plantar fascia under HHS conditions, in this study, musculoskeletal modeling (MsM) and finite element method (FEM) were used to investigate the effect of heel height on strain distribution of plantar fascia. Three-dimensional (3D) and one-dimensional (1D) finite element models of plantar fascia were generated to analyze the computed strain variation in 3-, 5-, and 7-cm heel heights. For validation, the computed foot contact pressure was compared with experimental measurement, and the strain value on 1D fascia was compared with previous studies. Results showed that the peak strain of plantar fascia was progressively increased on both 3D and 1D plantar fascia as heel elevated from 3 to 7 cm, and the maximum strain of plantar fascia occurs near the heel pain site at second peak stance. The 3D fascia model predicted a higher strain magnitude than that of 1D and provided a more reliable strain distribution on the plantar fascia. It is concluded that HHS with narrow heel support could pose a high risk on plantar fasciitis development, rather than reducing symptoms. Therefore, the heel elevation as a treatment recommendation for plantar fasciitis is questionable. Further studies of different heel support structures of shoes to quantify the effectiveness of heel elevation on the load-bearing mechanism of plantar fascia are recommended.

Avulsion injuries: an update on radiologic findings

Avulsion injuries result from the application of a tensile force to a musculoskeletal unit or ligament. Although injuries tend to occur more commonly in skeletally immature populations due to the weakness of their apophysis, adults may also be subject to avulsion fractures, particularly those with osteoporotic bones. The most common sites of avulsion injuries in adolescents and children are apophyses of the pelvis and knee. In adults, avulsion injuries commonly occur within the tendon due to underlying degeneration or tendinosis. However, any location can be involved in avulsion injuries. Radiography is the first imaging modality to diagnose avulsion injury, although advanced imaging modalities are occasionally required to identify subtle lesions or to fully delineate the extent of the injury. Ultrasonography has a high spatial resolution with a dynamic assessment potential and allows the comparison of a bone avulsion with the opposite side. Computed tomography is more sensitive for depicting a tiny osseous fragment located adjacent to the expected attachment site of a ligament, tendon, or capsule. Moreover, magnetic resonance imaging is the best imaging modality for the evaluation of soft tissue abnormalities, especially the affected muscles, tendons, and ligaments. Acute avulsion injuries usually manifest as avulsed bone fragments. In contrast, chronic injuries can easily mimic other disease processes, such as infections or neoplasms. Therefore, recognizing the vulnerable sites and characteristic imaging features of avulsion fractures would be helpful in ensuring accurate diagnosis and appropriate patient management. To this end, familiarity with musculoskeletal anatomy and mechanism of injury is necessary.

Single-Leg Landings Following a Volleyball Spike May Increase the Risk of Anterior Cruciate Ligament Injury More Than Landing on Both-Legs

Volleyball players often land on a single leg following a spike shot due to a shift in the center of gravity and loss of balance. Landing on a single leg following a spike may increase the probability of non-contact anterior cruciate ligament (ACL) injuries. The purpose of this study was to compare and analyze the kinematics and kinetics differences during the landing phase of volleyball players using a single leg (SL) and double-leg landing (DL) following a spike shot. The data for vertical ground reaction forces (VGRF) and sagittal plane were collected. SPM analysis revealed that SL depicted a smaller knee flexion angle (about 13.8°) and hip flexion angle (about 10.8°) during the whole landing phase, a greater knee and hip power during the 16.83–20.45% (p = 0.006) and 13.01–16.26% (p = 0.008) landing phase, a greater ankle plantarflexion angle and moment during the 0–41.07% (p < 0.001) and 2.76–79.45% (p < 0.001) landing phase, a greater VGRF during the 5.87–8.25% (p = 0.029), 19.75–24.14% (p = 0.003) landing phase when compared to DL. Most of these differences fall within the time range of ACL injury (30–50 milliseconds after landing). To reduce non-contact ACL injuries, a landing strategy of consciously increasing the hip and knee flexion, and plantarflexion of the ankle should be considered by volleyball players.

The biomechanical characteristics of a feline distal forelimb: A finite element analysis study

As a typical digitigrade mammal, the uniquely designed small distal limbs of the feline support two to three times of its body weight during daily movements. To understand how force transmission occurs in relation to the distal joint in a feline limb, which transfers bodyweight to the ground, it is necessary to examine the internal stress distribution of the distal joint limb in detail. Therefore, finite element models (FEM) of a healthy feline were established to predict the internal stress distribution of the distal limb. The FEM model included 23 bony components, various cartilaginous ligaments, as well as the encapsulated soft tissue of the paw. The FEM model was validated by comparison of paw pressure distribution, obtained from an experiment for balance standing. The results demonstrated a good agreement between the experimentally measured and numerically predicted pressure distribution in the feline paw. Additionally, higher stress levels were noted in the metacarpal segment, with smaller stresses observed in the phalanges portion including the proximal, middle, and distal segments. The raised metacarpal segment plays an important role in creating a stiff junction between the metacarpophalangeal (MCP) and wrist joint, stabilizing the distal limb. The paw pads help to optimize stress distribution in phalanx region. Findings from this study contribute to our understanding of feline distal forelimb biomechanical behavior. This information can be applied to bionic design of footwear since an optimal stiff junction and pressure distribution can be adapted to enhance injury relief and sports activities. Further developments may include progress, evaluation, and treatment of metatarsophalangeal joint injuries in human populations.

Anatomic Characteristics of Tissues Attached to the Fifth Metatarsal Bone

Background

Two types of stress, bending stress and traction stress, have been reported to be involved in the mechanism of Jones fracture. However, little is known about the risk factors for traction stress.

Purpose

To classify the attachment position of the peroneus brevis muscle (PB), peroneus tertius (PT), lateral band of the plantar aponeurosis (LB), and the long plantar ligament (LPL), focusing on the zone where a Jones fracture occurs (zone 2), and to compare the footprint area of each tissue type.

Study Design

Descriptive laboratory study.

Methods

This study examined 102 legs from 55 Japanese cadavers. Type classification was performed by focusing on the positional relationship between each tissue attachment and the zone where Jones fracture occurs (zone 2). The classifications were as follows: type I, attached proximal to the border between zones 1 and 2; type IIa, attached to the border between zones 1 and 2 with one attached part; and type IIb, attached across the border between zones 1 and 2 with two or more attached parts. The footprint areas of the PB, PT, LB, and LPL were compared between tissue types and within each attachment classification.

Results

The PB was recorded as type I in 41 feet (40.2%), type IIa in 56 feet (54.9%), and type IIb in 5 feet (4.9%); the PT was recorded as type IIa in 54 feet (60.0%) and type IIb in 36 feet (40.0%); and the LB was recorded as type I in 27 feet (26.5%) and type IIa in 75 feet (73.5%). The LPL did not attach to the fifth metatarsal bone. No significant difference was found in the footprint area between type I PB and type I LB.

Conclusion

The results indicate that type I, which attaches proximal to zone 2, occurs with PB and LB, and there was no significant difference in the footprint area between them. These findings suggest that type I is involved in traction stress. In the future, biomechanical research based on the results of this study will be necessary.

Clinical Relevance

The results of this study provide basic research for investigating the mechanism of Jones fracture and the cause of delayed healing.

Temporal Kinematic Differences between Forward and Backward Jump-Landing

Backward jump-landing during sports performance will result in dynamic postural instability with a greater risk of injury, and most research studies have focused on forward landing. Differences in kinematic temporal characteristics between single-leg and double-leg backward jump-landing are seldom researched and understood. The purpose of this study was to compare and analyze lower extremity kinematic differences throughout the landing phases of forward and backward jumping using single-leg and double-leg landings (FS and BS, FD and BD). Kinematic data were collected during the landing phases of FS and BS, FD and BD in 45 participants. Through statistical parametric mapping (SPM) analysis, we found that the BS showed smaller hip and knee flexion and greater vertical ground reactive force (VGRF) than the FS during 0–37.42% (p = 0.031), 16.07–32.11% (p = 0.045), and 23.03–17.32% (p = 0.041) landing phases. The BD showed smaller hip and knee flexion than the FD during 0–20.66% (p = 0.047) and 0–100% (p < 0.001) landing phases. Most differences appeared within a time frame during the landing phase at 30–50 ms in which non-contact anterior cruciate ligament (ACL) injuries are thought to occur and are consistent with the identification of risk in biomechanical analysis. A landing strategy that consciously increases the knee and hip flexion angles during backward landing should be considered for people as a measure to avoid injury during the performance of this type of physical activity.

Implication of obesity on motion, posture and internal stress of the foot: an experimental and finite element analysis

Obesity causes increased loading on the foot which can damage the soft tissue and bone ultimately leading to foot problems. Experimental and computational methods were used to analyse the chain of biomechanical changes in the lower limb due to obesity. The experimental study shows some changes in foot posture and gait where obese subjects were more likely to have pronated feet, smaller joint angles in the sagittal and frontal planes, smaller cadence, and smaller stride length. Anatomically correct finite element models generated on obese subjects showed increased and altered internal and plantar stress.  Altered foot posture was identified as a key indicator of increased internal stress indicating the importance of foot posture correction.

Cartilage Stiffness Effect on Foot Biomechanics of Chinese Bound Foot: A Finite Element Analysis

The purpose of this study is to investigate the effect of cartilage stiffness on inner foot biomechanics of Chinese bound foot while balanced standing using finite element method. A three-dimensional FE model of bound foot involving 28 bones, 72 ligaments, 5 plantar fascia, cartilages, and encapsulated soft tissue was constructed and validated. To conduct the sensitivity analysis of cartilage stiffness, the incremental Young’s modulus of 1, 5, 10, and 15 MPa were assigned to the cartilage. 25% of the body weight was applied to the Achilles tendon to adjust the anterior- posterior displacement of center of pressure agreeable with the measured result. As the Young’s modulus of cartilage increased, the peak von Mises stress in the fifth metatarsal increased obviously, while that in the calcaneus remains unchanged. The plantar fascia experienced reduced total tension with stiffer cartilage. The cartilage stiffening also caused a general increase of contact pressure at mid- and forefoot joints. Cartilage stiffening due to foot binding gave rise to risks of foot pain and longitude arch damage. Knowledge of this study contributes to the understanding of bound foot biomechanical behavior and demonstrating the mechanism of long-term injury and function damage in terms of weight-bearing due to foot binding.

Range limitation in hip internal rotation and fifth metatarsal stress fractures (Jones fracture) in professional football players

PURPOSE

To identify unknown risk factors associated with fifth metatarsal stress fracture (Jones fracture).

METHODS

A case-controlled study was conducted among male Japanese professional football (soccer) players with (N = 20) and without (N = 40) a history of Jones fracture. Injury history and physical examination data were reviewed, and the two groups were compared. Univariate and multivariate logistic regression controlling for age, leg dominance and body mass index were used to obtain odds ratios (ORs) and 95% confidence intervals (CIs) to describe the association between physical examination data and the presence or absence of Jones fractures.

RESULTS

From 2000 to 2014, among 162 professional football club players, 22 (13.6%; 21 Asians and one Caucasian) had a history of Jones fracture. Thirteen out of 22 (60%) had a Jones fracture in their non-dominant leg. The mean range of hip internal rotation (HIR) was restricted in players with a history of Jones fracture [25.9° ± 7.5°, mean ± standard deviation (SD)] compared to those without (40.4° ± 11.1°, P < 0.0001). Logistic regression analyses demonstrated that HIR limitation increased the risk of a Jones fracture (OR = 3.03, 95% CI 1.45-6.33, P = 0.003). Subgroup analysis using data prior to Jones fracture revealed a causal relationship, such that players with a restriction of HIR were at high risk of developing a Jones fracture [Crude OR (95% CI) = 6.66 (1.90-23.29), P = 0.003, Adjusted OR = 9.91 (2.28-43.10), P = 0.002]. In addition, right HIR range limitation increased the risks of developing a Jones fracture in the ipsilateral and the contralateral feet [OR = 3.11 (1.35-7.16) and 2.24 (1.22-4.12), respectively]. Similarly, left HIR range limitation increased the risks in the ipsilateral or the contralateral feet [OR (95% CI) = 4.88 (1.56-15.28) and 2.77 (1.08-7.08), respectively].

CONCLUSION

The restriction of HIR was associated with an increased risk of developing a Jones fracture. Since the HIR range is a modifiable factor, monitoring and improving the HIR range can lead to prevent reducing the occurrence of this fracture.

LEVEL OF EVIDENCE

III.

Fifth metatarsal stress fracture in elite male football players: an on-field analysis of plantar loading

Objective

Evaluate plantar loading during ‘on-field’ common football movements in players after fifth metatarsal (MT-5) stress fracture and compare with matched healthy players.

Methods

Fourteen elite male soccer players participated in the study conducted on a natural grass playing surface using firm ground football boots. Seven players who had suffered a primary stress fracture (MT-5 group) and seven matched healthy players (controls, CON) performed three common football movements while in-shoe plantar loading data were collected.

Results

Large between-group differences exist for maximal vertical force normalised to bodyweight (F_max) at the lateral toes (2-5) of the stance leg during a set-piece kick (MT-5: 0.2±0.06 bodyweight (BW), CON: 0.1±0.05 BW, effect size (ES) 1.4) and the curved run where the MT-5 group showed higher F_max with very large effect size at the lateral forefoot of the injured (closest to curve) limb when running a curve to receive a pass (MT-5 injured−CON=0.01 BW, ES 1.5). Small between-group differences were evident during straight-line running. However, between-limb analysis of MT-5 group showed significant unloading of the lateral forefoot region of the involved foot.

Conclusions

Elite male football players who have returned to play after MT-5 stress fracture display significantly higher maximum plantar force at the lateral forefoot and lateral toes (2-5) compared with healthy matched control players during two football movements (kick and curved run) with the magnitude of these differences being very large. These findings may have important implications for manipulating regional load during rehabilitation or should a player report lateral forefoot prodromal symptoms.

Subject-specific finite element modelling of the human foot complex during walking: sensitivity analysis of material properties, boundary and loading conditions

The objective of this study was to develop and validate a subject-specific framework for modelling the human foot. This was achieved by integrating medical image-based finite element modelling, individualised multi-body musculoskeletal modelling and 3D gait measurements. A 3D ankle-foot finite element model comprising all major foot structures was constructed based on MRI of one individual. A multi-body musculoskeletal model and 3D gait measurements for the same subject were used to define loading and boundary conditions. Sensitivity analyses were used to investigate the effects of key modelling parameters on model predictions. Prediction errors of average and peak plantar pressures were below 10% in all ten plantar regions at five key gait events with only one exception (lateral heel, in early stance, error of 14.44%). The sensitivity analyses results suggest that predictions of peak plantar pressures are moderately sensitive to material properties, ground reaction forces and muscle forces, and significantly sensitive to foot orientation. The maximum region-specific percentage change ratios (peak stress percentage change over parameter percentage change) were 1.935-2.258 for ground reaction forces, 1.528-2.727 for plantar flexor muscles and 4.84-11.37 for foot orientations. This strongly suggests that loading and boundary conditions need to be very carefully defined based on personalised measurement data.

Management of sports injuries of the foot and ankle: An update

Injuries to the foot in athletes are often subtle and can lead to a substantial loss of function if not diagnosed and treated appropriately. For these injuries in general, even after a diagnosis is made, treatment options are controversial and become even more so in high level athletes where limiting the time away from training and competition is a significant consideration. In this review, we cover some of the common and important sporting injuries affecting the foot including updates on their management and outcomes. Cite this article: Bone Joint J 2016;98-B:1299-1311.

GAIT ANALYSIS DRIVEN 2D FINITE ELEMENT MODEL OF THE NEUROPATHIC HINDFOOT

The diabetic foot is one of the most serious complications of diabetes mellitus and it can lead to foot ulcerations and amputations. Finite element analysis quantifies the loads developed in the different anatomical structures and describes how these affect foot tissue during foot–floor interaction. This approach for the diabetic subjects’ foot could provide valuable information in the process of plantar orthosis fabrication and fit. The purpose of this study was to develop two finite element models of the hindfoot, of healthy and diabetic neuropathic subjects. These models accounts for in vivo kinematics, kinetics, plantar pressure (PP) data and magnetic resonance images. These were acquired during gait analysis on 10 diabetic neuropathics and 10 healthy subjects. Validity of the models has been assessed through comparison between the peak PPs of simulated and experimental data: root mean square error (RMSE) in percentage of the experimental peak value was evaluated. Two different finite elements analysis were performed: subject-specific simulations in terms of both geometry and gait analysis, and by adopting the complete gait analysis dataset as boundary conditions. Model predicted plantar pressures were in good agreement with those experimentally measured. Best agreement was obtained in the subject-specific case (RMSE of 13%).

Constitutive formulation and numerical analysis of the biomechanical behaviour of forefoot plantar soft tissue

The aim of this work is to provide a numerical approach for the investigation of the mechanical behaviour of the forefoot soft tissues. The development of reliable numerical models of biological structures requires the definition of constitutive formulations that actually interpret the mechanical response of the constituent biological tissues and their structural arrangement. A specific visco-hyperelastic constitutive model is provided to account for the typical features of soft plantar tissue mechanics, as geometric and material non-linearity, almost-incompressible behaviour and time-dependent phenomena. Constitutive parameters are evaluated by the analysis of experimental data from compression and stress relaxation tests on tissue samples. A three-dimensional finite element model of the forefoot region is developed starting from the analysis of biomedical images, leading to the evaluation of overall structural response. The reliability of model and analyses is assessed by the comparison of experimental and numerical results pertaining to indentation tests. The numerical model developed allows to evaluate the mechanical response of plantar soft tissue in terms of stress and strain distribution.

MECHANICAL PROPERTIES OF THE HUMAN METATARSAL BONES

Despite the incidence of metatarsal fractures and the associated risk of significant disability, little is known about the biomechanical properties (strength and stiffness) of metatarsal bones. In most cases a single metatarsal bone (first, second and fifth) has been investigated. An extensive investigation of the biomechanical properties of the metatarsal bones is essential in the understanding and prevention of metatarsal injuries. Entire sets of metatarsal bones from four feet were tested. The first foot was used to fine-tune the testing set-ups. To measure the stiffness, each metatarsal bone was subjected to non-destructive four-point-bending in the sagittal and transverse planes, axial compression and torsion. Strain was measured at two locations. To measure the strength, each metatarsal bone was tested to failure in torsion. Significant differences (p < 0.0001) existed among the stiffness of the five metatarsal bones: (i) in torsion the first metatarsal bone was 2–3 times stiffer than the others; (ii) in four-point-bending and axial compression this difference was less pronounced than in torsion; (iii) differences were smaller among the other metatarsal bones; (iv) the second metatarsal bone was less stiff than the third and fourth in bending. The second, third and fourth metatarsal bones were stiffer in the sagittal than in the transverse plane (p < 0.0001). Conversely, there was no significant difference between the two planes of bending for the first and fifth bones. During destructive testing, all metatarsal bones exhibited a linear elastic behavior and brittle failure. The torsional strength at failure ranged between 1.9 Nm and 6.9 Nm. The first metatarsal bone was stronger than all the others. Stiffness in different loading conditions and failure were measured and compared for all metatarsal bones. These data corroborate previous biomechanical studies concerning the role and load sharing of the different metatarsal bones.

Effects of foot posture on fifth metatarsal fracture healing: a finite element study

The goal of this study was to evaluate the effects of maintaining different foot postures during healing of proximal fifth metatarsal fractures for each of 3 common fracture types. A 3-dimensional (3D) finite element model of a human foot was developed and 3 loading situations were evaluated, including the following: (1) normal weightbearing, (2) standing with the affected foot in dorsiflexion at the ankle, and (3) standing with the affected foot in eversion. Three different stages of the fracture-healing process were studied, including: stage 1, wherein the material interposed between the fractured edges was the initial connective tissue; stage 2, wherein connective tissue had been replaced by soft callus; and stage 3, wherein soft callus was replaced by mature bone. Thus, 30 3D finite element models were analyzed that took into account fracture type, foot posture, and healing stage. Different foot postures did not statistically significantly affect the peak-developed strains on the fracture site. When the fractured foot was everted or dorsiflexed, it developed a slightly higher strain within the fracture than when it was in the normal weightbearing position. In Jones fractures, eversion of the foot caused further torsional strain and we believe that this position should be avoided during foot immobilization during the treatment of fifth metatarsal base fractures. Tuberosity avulsion fractures and Jones fractures seem to be biomechanically stable fractures, as compared with shaft fractures. Our understanding of the literature and experience indicate that current clinical observations and standard therapeutic options are in accordance with the results that we observed in this investigation, with the exception of Jones fractures.