Stress distribution of metatarsals during forefoot strike versus rearfoot strike: A finite element study

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.

Case study: The influence of Achilles tendon rupture on knee joint stress during counter-movement jump - Combining musculoskeletal modeling and finite element analysis

BACKGROUND

Presently, the current research concerning Achilles tendon rupture repair (ATR) is predominantly centered on the ankle joint, with a paucity of evidence regarding its impact on the knee joint. ATR has the potential to significantly impede athletic performance and increase tibiofemoral contact forces in athletes. The purpose of this study was to prognosticate the distribution of stress within the knee joint during a countermovement jump through the use of a simulation method that amalgamated a musculoskeletal model of a patient who underwent Achilles tendon rupture repair with a finite element model of the knee joint.

METHODS

A male elite badminton player who had suffered an acute Achilles tendon rupture in his right leg one year prior was selected as our study subject. In order to analyze his biomechanical data, we employed both the OpenSim musculoskeletal model and finite element model to compute various parameters such as joint angles, joint moments, joint contact forces, and the distribution of knee joint stress.

RESULTS

During the jumping phase, a significantly lower knee extension angle (p < 0.001), ankle dorsiflexion angle (p = 0.002), peak vertical ground reaction force (p < 0.001), and peak tibiofemoral contact force (p = 0.009) were observed on the injured side than on the uninjured side. During the landing phase, the ankle range of motion (ROM) was significantly lower on the injured side than on the uninjured side (p = 0.009), and higher peak vertical ground reaction forces were observed (p = 0.012). Additionally, it is logical that an injured person will put higher load on the uninjured limb, but the finite element analysis indicated that the stresses on the injured side of medial meniscus and medial cartilage were significantly greater than the uninjured side.

CONCLUSIONS

An Achilles tendon rupture can limit ankle range of motion and lead to greater joint stress on the affected area during countermovement jumps, especially during the landing phase. This increased joint stress may also transfer more stress to the soft tissues of the medial knee, thereby increasing the risk of knee injury. It is worth noting that this study only involves the average knee flexion angle and load after ATR in one athlete. Caution should be exercised when applying the conclusions, and in the future, more participants should be recruited to establish personalized knee finite element models to validate the results.

Asymmetries of foot strike patterns during running in high-level female and male soccer players

BACKROUND

Foot strike pattern (FSP) is defined by the way the foot makes initial ground contact and is influenced by intrinsic and extrinsic factors. This study investigated the effect of running speed on asymmetries of FSP.

METHODS

Seventeen female and nineteen male soccer players performed an incremental running test on an instrumented treadmill starting at 2.0 m/s until complete exhaustion. Force plate data were used to categorize foot strikes into rearfoot (RFS) and non-rearfoot strikes. Additionally, peak vertical ground reaction force (peakGRF) and stride time were calculated. The symmetry index (SI) was used to quantify lateral asymmetries between legs.

RESULTS

The SI indicated asymmetries of the rate of RFS (%RFS) of approximately 30% at slow running speed which decreased to 4.4% during faster running speed (p = 0.001). There were minor asymmetries in peakGRF and stride time at each running stage. Running speed influenced %RFS (p < 0.001), peakGRF (p < 0.001) and stride time (p < 0.001). Significant interaction effects between running speed and sex were shown for %RFS (p = 0.033), peakGRF (p < 0.001) and stride time (p = 0.041).

CONCLUSION

FSP of soccer players are asymmetric at slower running speed, but symmetry increases with increasing speed. Future studies should consider that FSP are non-stationary and influenced by running speed but also differ between legs.

Effects of Midsole Hardness on the Mechanical Response Characteristics of the Plantar Fascia during Running

High long-term stress on the plantar fascia (PF) is the main cause of plantar fasciitis. Changes in the midsole hardness (MH) of running shoes are an important factor leading to the alteration of the PF. This study aims to establish a finite-element (FE) model of the foot-shoe, and investigates the effects of midsole hardness on PF stress and strain. The FE foot-shoe model was built in ANSYS using computed-tomography imaging data. Static structural analysis was used to simulate the moment of running push and stretch. Plantar stress and strain under different MH levels were quantitatively analyzed. A complete and valid 3D FE model was established. With an increase in MH from 10 to 50 Shore A, the overall stress and strain of the PF were decreased by approximately 1.62%, and the metatarsophalangeal (MTP) joint flexion angle was decreased by approximately 26.2%. The height of the arch descent decreased by approximately 24.7%, but the peak pressure of the outsole increased by approximately 26.6%. The established model in this study was effective. For running shoes, increasing the MH reduces the stress and strain of PF, but also imposes a higher load on the foot.

The influence of running shoe with different carbon-fiber plate designs on internal foot mechanics: A pilot computational analysis

A carbon-fiber plate (CFP) embedded into running shoes is a commonly applied method to improve running economy, but little is known in regard the effects of CFP design features on internal foot mechanics. This study aimed to explore how systematic changes in CFP geometrical variations (i.e., thickness and location) can alter plantar pressure and strain under the forefoot as well as metatarsal stress state through computational simulations. A foot-shoe finite element (FE) model was built and different CFP features including three thicknesses (1 mm, 2 mm, and 3 mm) and three placements (high-loaded (just below the insole), mid-loaded (in between the midsole), and low-loaded (just above the outsole)) were further modulated within the shoe sole. Simulations were conducted at the impact peak instant during forefoot strike running. Compared with the no-CFP shoe, peak plantar pressure and compressive strain under the forefoot consistently decreased when the CFP thickness increased, and the low-loaded conditions were found more effective (peak pressure decreased up to 31.91% and compressive strain decreased up to 18.61%). In terms of metatarsal stress, CFP designs resulted in varied effects and were dependent on their locations. Specifically, high-loaded CFP led to relatively higher peak metatarsal stress without the reduction trend as thickness increased (peak stress increased up to 12.91%), while low-loaded conditions showed a gradual reduction in peak stress, decreasing by 0.74%. Therefore, a low-loaded thicker CFP should be considered to achieve the pressure-relief effects of running shoes without the expense of increased metatarsal stress.

Biomechanical Analysis of Foot–Ankle Complex during Jogging with Rearfoot Strike versus Forefoot Strike

Background and Aim. In order to reduce foot and ankle injuries induced by jogging, two-foot strike patterns, rearfoot strike (RFS), and forefoot strike (FFS), were adopted and compared. First, RFS jogging and FFS jogging were experimentally studied, so as to acquire kinematic and kinetic data, including foot strike angle, knee flexion angle, and ground reaction force (GRF). Then, a 3D finite element model of foot–ankle complex was reconstructed from the scanned 2D-stacked images. Biomechanical characteristics, including plantar pressure, stress of metatarsals, midfoot bone, calcaneus and cartilage, and tensile force of plantar fascia and ligaments, were obtained. The results showed that RFS jogging and FFS jogging had a similar change trend and a close peak value of GRF. Since possessing more momentum in the push stage and less momentum in the brake stage, FFS jogging could be in favor of a higher jogging speed. However, FFS jogging produced larger metatarsal stress in the 5th metatarsal and much larger tensile force of plantar fascia, which might cause metatarsal fracture and heel pain. While RFS jogging produced larger plantar pressure in the hindfoot area, larger calcaneus stress, and much larger tarsal navicular stress, which might cause heel tissue injury, calcaneus damage, and stress fracture of naviculocuneiform joint. In addition, talocrural and talocalcaneal joint cartilage could bear jogging loads, as the peak contact pressure were both small in RFS jogging and FFS jogging. Therefore, jogging with rearfoot or FFS pattern should be chosen according to the health condition of foot–ankle parts.

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.

Finite element modelling for footwear design and evaluation: A systematic scoping review

Finite element modelling has become an efficient tool for an in-depth understanding of the foot, footwear biomechanics and footwear optimization. The aim of this paper was to provide an updated overview in relation to the footwear finite element (FE) analysis published since 2000. The paper will attempt to outline the main challenges and research gaps that need confronting in the further development of realistic and accurate models for clinical and industrial applications. English databases of the Web of Science and PubMed were used to search (‘finite element’ OR ‘FEA’ OR ‘computational model’) AND (‘shoe’ OR ‘footwear’) until 16 December 2021. Articles that conducted FE analyses on the whole foot and footwear structures were included in this review. Twelve articles met the eligibility criteria, and were grouped into three categories for further analysis, (1) finite element modelling of the foot and high-heeled shoes; (2) finite element modelling of the foot and boot; (3) finite element modelling of the foot and sports shoe. Even though most of the existing foot-shoe FE analyses were performed under certain simplifications and assumptions, they have provided essential contributions in identifying the mechanical response of the foot in casual or athletic footwear. Further to this, the results have provided information in relation to optimizing footwear design to enhance functional performance. Nevertheless, further simulations still present several challenges, including reliable data information for geometry reconstruction, the balance between accurate details and computational cost, accurate representations of material properties, realistic boundary and loading conditions, and thorough model validation. In addition, some research gaps in terms of the coverage of footwear design, the consideration of insole/orthosis and socks, and the internal and external validity of the FE design should be fully covered.

3D printed auxetic heel pads for patients with diabetic mellitus

More than 422 million people worldwide suffered from diabetes mellitus (DM) in 2021. Diabetic foot is one the most critical complications resultant of DM. Foot ulceration and infection are frequently arisen, which are associated with changes in the mechanical properties of the plantar soft tissues, peripheral arterial disease, and sensory neuropathy. Diabetic insoles are currently the mainstay in reducing the risk of foot ulcers by reducing the magnitude of the pressure on the plantar Here, we propose a novel pressure relieving heel pad based on a circular auxetic re-entrant honeycomb structure by using three-dimensional (3D) printing technology to minimize the pressure on the heel, thus reducing the occurrence of foot ulcers. Finite element models (FEMs) are developed to evaluate the structural changes of the developed circular auxetic structure upon exertion of compressive forces. Moreover, the effects of the internal angle of the re-entrant structure on the peak contact force and the mean pressure acting on the heel as well as the contact area between the heel and the pads are investigated through a finite element analysis (FEA). Based on the result from the validated FEMs, the proposed heel pad with an auxetic structure demonstrates a distinct reduction in the peak contact force (∼10%) and the mean pressure (∼14%) in comparison to a conventional diabetic insole (PU foam). The characterized result of the designed circular auxetic structure not only provides new insights into diabetic foot protection, but also the design and development of various impact resistance products.

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.

A Review of Finite Element Models of Ligaments in the Foot and Considerations for Practical Application

PURPOSE

Finite element (FE) modeling has been used as a research tool for investigating underlying ligaments biomechanics and orthopedic applications. However, FE models of the ligament in the foot have been developed with various configurations, mainly due to their complex 3D geometry, material properties, and boundary conditions. Therefore, the purpose of this review was to summarize the current state of finite element modeling approaches that have been used in the ?eld of ligament biomechanics, to discuss their applicability to foot ligament modeling in a practical setting, and also to acknowledge current limitations and challenges.

METHODS

A comprehensive literature search was performed. Each article was analyzed in terms of the methods used for: (a) ligament geometry, (b) material property, (c) boundary and loading condition related to its application, and (d) model verification and validation.

RESULTS

Of the reviewed studies, 80% of the studies used simplified representations of ligament geometry, the non-linear mechanical behavior of ligaments was taken into account in only 19.2% of the studies, 33% of included studies did not include any kind of validation of the FE model.

CONCLUSION

Further refinement in the functional modeling of ligaments, the micro-structure level characteristics, nonlinearity, and time-dependent response, may be warranted to ensure the predictive ability of the models.

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.

Do non-rearfoot runners experience greater second metatarsal stresses than rearfoot runners?

Stress fracture of the second metatarsal is a common and problematic injury for runners. The choice of foot strike pattern is known to affect external kinetics and kinematics but its effect on internal loading of the metatarsals is not well understood. Models of various complexities can be used to investigate the effects of running characteristics on metatarsal stresses. This study aimed to compare second metatarsal stress between habitual rearfoot and non-rearfoot strikers during barefoot running, using a novel participant-specific finite element model, including accurate metatarsal and soft tissue geometry. Synchronised force and kinematic data were collected during barefoot overground running from 20 participants (12 rearfoot strikers). Stresses were calculated using a previously evaluated and published 3D finite element model. Non-rearfoot strikers demonstrated greater external loading and joint contact forces than rearfoot runners, but there were no differences in stresses between groups. Additionally, the study allowed for a qualitative assessment of bone geometries and stresses. No correlation was found between bone volume and stresses, however, there was found to be a large variation in metatarsal shapes, possibly accounting for the lack of difference in stresses. This emphasises the importance of bone geometry when estimating bone stress and supports the suggestion that external forces should not be assumed to be representative of internal loading.

Associations between changes in loading pattern, deformity, and internal stresses at the foot with hammer toe during walking; a finite element approach

Over the past decade, Finite Element (FE) modelling has been used as a method to understand the internal stresses within the diabetic foot. Foot deformities such as hammer toe have been associated with increased risk of foot ulcers in diabetic patients. Hence the aim of this study is to investigate the influence of hammer toe deformity on internal stresses during walking. A 3D finite element model of the human foot was constructed based on capturing Magnetic Resonance Imaging (MRI) of a diabetic neuropathic volunteer exhibiting hammer toe. 3D gait measurements and a multi-body musculoskeletal model for the same participant were used to define muscle forces. FE simulations were run at five different instances during the stance phase of gait. Peak plantar pressure and pressure distribution results calculated from the model showed a good agreement with the experimental measurement having less than 11% errors. Maximum von Mises internal stresses in the forefoot hard tissue were observed at the 3^rd and 5^th metatarsals and 4^th proximal phalanx. Moreover, presence of hammer toe deformity was found to shift the location of maximum internal stresses on the soft tissue to the forefoot by changing the location of centre of pressure with internal stress 1.64 times greater than plantar pressure. Hammer toe deformity also showed to reduce the involvement of the first phalanx in internal/external load-bearing during walking. The findings of this study support the association between changes in loading pattern, deformity, and internal stresses in the soft tissue that lead to foot ulceration.

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.

A Narrative Review of Metatarsal Bone Stress Injury in Athletic Populations: Etiology, Biomechanics, and Management

Metatarsal bone stress injuries (BSIs) are common in athletic populations. BSIs are overuse injuries that result from an accumulation of microdamage that exceeds bone remodeling. Risk for metatarsal BSI is multifactorial and includes factors related to anatomy, biology, and biomechanics. In this article, anatomic factors including foot type, metatarsal length, bone density, bone geometry, and intrinsic muscle strength, which each influence how the foot responds to load, are discussed. Biologic factors such as low energy availability and impaired bone metabolism influence the quality of the bone. Finally, the influence of biomechanical loads to bone such as peak forces, load rates, and loading cycles are reviewed. General management of metatarsal BSI is discussed, including acute care, rehabilitation, treatment of refractory metatarsal BSI, and evaluation of healing/return to sport. Finally, we identify future research priorities and emerging treatments for metatarsal BSI.

Plantar load characteristics among runners with different strike patterns during preferred speed

Objectives

This study aimed to compare the plantar loads between habitual rearfoot strike (RFS) and non-RFS (NRFS) during running under the participant’s preferred speed.

Methods

A total of 66 (36 RFS, 30 NRFS) healthy amateur male runners were included in our study. In-shoe pressure sensors were utilised to the test plantar loads when participants were running using their preferred foot strike pattern and running speed (RFS: 3.2 ± 0.3 m/s; NRFS: 3.4 ± 0.4 m/s).

Results

Results indicated that running speed has a significant effect on the total contact area [F (1, 64) = 7.061, P = 0.01, η² = 0.101], which also affects midfoot and forefoot regions. No significant difference was found on the total maximum force, force-time-integral, peak pressure (PP) and pressure-time-integral (PTI), but the total contact area of RFS was higher than that of NRFS runners [F (1, 64) = 77.406, P < 0.001, η² = 0.551]. Plantar loads were mainly focused on the heel and midfoot for RFS runners in all variables, and NRFS runners experienced increased PP and PTI in medial forefoot regions.

Conclusion

Habitual runners tend to adjust their contact area according to the running speed through midfoot and forefoot regions. RFS runners remain susceptible to high impact force on the heel and midfoot, and NRFS runners experience high impact force in the first metatarsal regions. Therefore, runners should note this situation to avoid running-related injuries.

Comparison of plantar loads among runners with different strike patterns

This study aimed to explore the plantar loading variables between habitual rearfoot strike (RFS) and non-rearfoot strike (NRFS) during running. 78 healthy males participated in this study (41 RFS, 37 NRFS). In-shoe pressure sensors were used to measure plantar loading while the participants were running on a 15 m indoor runway with their preferred foot strike pattern (FSP) at 12.0 ± 5% km/h. Results indicate that force and pressure parameters were much higher in the rearfoot and midfoot regions during RFS running and relatively greater in forefoot region during NRFS running. However, compared with NRFS running, the contact area, maximum force and force-time-integrals during RFS running on total foot were 21.44% (P < 0.001, ES = 2.29), 13.99% (P = 0.006, ES = 0.64) and 21.27% (P < 0.001, ES = 0.85) higher, respectively. Total foot peak pressure and pressure-time-integral between two FSPs were similar. Higher loads in the rearfoot region may transmit to the knee joint and result in patellofemoral joint injuries. NRFS runners' higher loads in forefoot seem to be ralated to metatarsal stress fractures and compensatory damage to the Achilles tendon. Therefore, runners should choose proper FSPs according to their unique physical conditions.

The effect of running on foot muscles and bones: A systematic review

Despite the widespread evidence of running as a health-preserving exercise, little is known concerning its effect on the foot musculature and bones. While running may influence anatomical foot adaptation, it remains unclear to what extent these adaptations occur. The aim of this paper is to provide a systematic review of the studies that investigated the effects of running and the adaptations that occur in foot muscles and bones. The search was performed following the PRISMA guidelines. Relevant keywords were used for the search through PubMed/MEDLINE, Scopus and SPORTDiscus. The methodological quality of intervention studies was assessed using the Downs and Black checklist. For cross-sectional studies, the Newcastle-Ottawa scale was used. Sixteen studies were found meeting the inclusion criteria. In general, the included studies were deemed to be of moderate methodological quality. Although results of relevant literature are limited and somewhat contradictory, the outcome suggests that running may increase foot muscle volume, muscle cross-sectional area and bone density, but this seems to depend on training volume and experience. Future studies conducted in this area should aim for a standard way of reporting foot muscle/bone characteristics. Also, herein, suggestions for future research are provided.