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

A Systematic Review of Finite Element Analysis in Running Footwear Biomechanics: Insights for Running-Related Musculoskeletal Injuries

This study presented a systematic review of recent advancements in the application of finite element (FE) methods to running and running shoe biomechanics. It focused on outlining the general approach to build foot-running shoe FE models, exploring their current applications and challenges, and providing directions for future research. The review also aimed to highlight the gap between theoretical mechanical responses in simulations and real-world manifestations of running-related musculoskeletal injuries (RRMI). A comprehensive search of electronic databases, including Web of Science, PubMed, and Scopus, identified 12 eligible articles for inclusion in this review. Current studies have examined the effects of various running shoe design features and conditions on the mechanical response of internal foot tissues using foot-running shoe FE models. These models have gradually evolved from simplified local representations to more realistic and comprehensive models, with the incorporation of experimental data enhancing simulation accuracy. However, to further improve simulation outcomes, key advancements are proposed to reduce development time and enhance model robustness. These include high-fidelity 3D model development, personalized shape transformation, AI-driven automated reconstruction, comprehensive dynamic running simulations, and improved validation methods. More importantly, future research needs to bridge the gap between FE simulations and RRMI risk by addressing the complexities of bone fracture criteria and conducting localized assessments of bone properties. Overall, this review provided valuable insights for biomedical engineers, medical professionals, and researchers, facilitating more accurate investigations of foot-running shoe FE models. Ultimately, these advancements aim to improve footwear design and training programs to reduce the risk of RRMI.

Testing the effects of footwear on biomechanics of human body: A review

Studies show the high potential of shoes to impact human movement and reduce the risk of injuries during normal and high-demanding activities. This review will delve into the existing literature on mechanical and biomechanical tests of shoes and their effects on the human body. Mechanical tests mainly include compression, bending, torsional flexibility, and impact tests. Biomechanical tests, on the other hand, study the kinematics and kinetics of the human body while performing different tasks. The primary goal of this review is to highlight the importance of isolating parameters in shoe design and testing to achieve optimal results in providing comfort, support, and injury prevention. Key conclusions include the influence of lattice structures on shoe stiffness and stress distribution, the effectiveness of composite loofah sponge for vibration damping, the benefits of Poron insoles for impact attenuation, the potential injury risk reduction with auxetic shoes, and the need for future research on mechanical tests, parameter investigation, and optimization of shoe sole structures.

The influence of simulated worn shoe and foot inversion on heel internal biomechanics during running impact: A subject-specific finite element analysis

This study explored how systematic changes in running shoe degradation and foot inversion alter the distribution and peak value of heel pressure and calcaneus stress, as well as the total stress-concentration exposure (TSCE) within the calcaneal bone. A foot-shoe finite element model was employed and three shoe wear conditions (new shoe (CON), moderate worn shoe (MWSC), excessive worn shoe (EWSC)) coupled with three foot inversion angles (0°, 10°, 20°) were further modulated. Simulations were conducted at the impact peak instant during running. Compared to CON0, heel pressure during neutral landings shifted medially and increased with progressive shoe wear, peaking under EWSC0. This shift expanded the high-pressure area by 1.333 cm2 and raised peak pressure by 24.42 %. Foot inversion landings exhibited an opposite trend: increased shoe wear promoted balanced pressure distribution, centralizing the load and eliminating high-pressure areas under EWSC10, where peak pressure was 11.36 % lower than CON10. Calcaneus stress during neutral landings, initially concentrated on the medial calcaneal surface and inferior tuberosity, intensified with wear, expanding high-stress area by 5.276 cm2 and increasing peak stress by 22.79 % under EWSC0. For foot inversion, the high-stress region shifted to the inferior tuberosity, with wear reducing peak stress by 10.41 % and eliminating high-stress area in EWSC10 compared to CON10. TSCE analysis revealed that EWSC10 had the lowest stress exposure (0 %kPa) across all conditions. Worn-out shoes would exacerbate heel internal biomechanics, while these effects may be mitigated by foot inversion, likely due to the formation of a relatively flat and larger contact area between the lateral sole and the ground.

Curved carbon-plated shoe may further reduce forefoot loads compared to flat plate during running

Using a curved carbon-fiber plate (CFP) in running shoes may offer notable performance benefit over flat plates, yet there is a lack of research exploring the influence of CFP geometry on internal foot loading during running. The objective of this study was to investigate the effects of CFP mechanical characteristics on forefoot biomechanics in terms of plantar pressure, bone stress distribution, and contact force transmission during a simulated impact peak moment in forefoot strike running. We employed a finite element model of the foot-shoe system, wherein various CFP configurations, including three stiffnesses (stiff, stiffer, and stiffest) and two shapes (flat plate (FCFP) and curved plate (CCFP)), were integrated into the shoe sole. Comparing the shoes with no CFP (NCFP) to those with CFP, we consistently observed a reduction in peak forefoot plantar pressure with increasing CFP stiffness. This decrease in pressure was even more notable in a CCFP demonstrating a further reduction in peak pressure ranging from 5.51 to 12.62%, compared to FCFP models. Both FCFP and CCFP designs had a negligible impact on reducing the maximum stress experienced by the 2nd and 3rd metatarsals. However, they greatly influenced the stress distribution in other metatarsal bones. These CFP designs seem to optimize the load transfer pathway, enabling a more uniform force transmission by mainly reducing contact force on the medial columns (the first three rays, measuring 0.333 times body weight for FCFP and 0.335 for CCFP in stiffest condition, compared to 0.373 in NCFP). We concluded that employing a curved CFP in running shoes could be more beneficial from an injury prevention perspective by inducing less peak pressure under the metatarsal heads while not worsening their stress state compared to flat plates.

Cushioning mechanism of the metatarsals during landing for the skateboarding ollie maneuver

Skateboarding is an Olympic event with frequent jumping and landing, where the cushioning effect by the foot structure (from the arch, metatarsals, etc.) and damping performance by sports equipment (shoes, insoles, etc.) can greatly affect an athlete's sports performance and lower the risk of limb injury. Skateboarding is characterized by the formation of a "man-shoe-skateboard system," which makes its foot cushioning mechanism different from those of other sports maneuvers, such as basketball vertical jump and gymnastics broad jump. Therefore, it is necessary to clarify the cushioning mechanism of the foot structure upon landing on a skateboard. To achieve this, a multibody finite element model of the right foot, shoe, and skateboard was created using Mimics, Geomagic, and ANSYS. Kinetic data from the ollie maneuver were used to determine the plantar pressure and Achilles tendon force at three characteristics (T1, T2, and T3). The stress and strain on the foot and metatarsals (MT1-5) were then simulated. The simulation results had an error of 6.98% compared to actual measurements. During landing, the force exerted on the internal soft tissues tends to increase. The stress and strain variations were highest on MT2, MT3, and MT4. Moreover, the torsion angle of MT1 was greater than those of the other metatarsals. Additionally, the displacements of MT2, MT3, and MT4 were higher than those of the other parts. This research shows that skateboarders need to absorb the ground reaction force through the movements of the MTs for ollie landing. The soft tissues, bones, and ligaments in the front foot may have high risks of injury. The developed model serves as a valuable tool for analyzing the foot mechanisms in skateboarding; furthermore, it is crucial to enhance cushioning for the front foot during the design of skateboard shoes to reduce potential injuries.

Modeling and Analysis of Foot Function in Human Gait Using a Two-Degrees-of-Freedom Inverted Pendulum Model with an Arced Foot

Gait models are important for the design and control of lower limb exoskeletons. The inverted pendulum model has advantages in simplicity and computational efficiency, but it also has the limitations of oversimplification and lack of realism. This paper proposes a two-degrees-of-freedom (DOF) inverted pendulum walking model by considering the knee joints for describing the characteristics of human gait. A new parameter, roll factor, is defined to express foot function in the model, and the relationships between the roll factor and gait parameters are investigated. Experiments were conducted to verify the model by testing seven healthy adults at different walking speeds. The results demonstrate that the roll factor has a strong relationship with other gait kinematics parameters, so it can be used as a simple parameter for expressing gait kinematics. In addition, the roll factor can be used to identify walking styles with high accuracy, including small broken step walking at 99.57%, inefficient walking at 98.14%, and normal walking at 99.43%.

The Influence of Materials on Footwear Behaviour: A Finite Element Simulation Study

The objective of this study was to analyse the influence of materials and their position within the upper assembly on the behaviour of casual footwear using finite element simulation tools. The study was carried out on three models of casual footwear, which are identical in terms of design lines, varying only in the materials of the upper assembly, namely calfskin leather (M1), knitted fabric (M2), and combination of knitted fabric and calfskin leather (M3). The footwear models were designed according to the design constraints specific to casual footwear. The foot was reconstructed based on the shoe last obtained based on anthropometric data. Material definition, 3D models editing, setting up analysis conditions, and constraints were performed using the Ansys 17.2 software. Gait biomechanics were taken into account to define the loading model, force distribution, force values, and constraints. The study evaluates footwear behaviour in terms of directional deformation (Z axis), equivalent von Mises stress, and equivalent elastic strain distribution. This paper explores a methodology that has the potential to enhance the footwear design and manufacturing process, providing designers with information about the deformations and stress distribution on upper parts of the footwear product.

Effects of heel apex position, apex angle and rocker radius on plantar pressure in the heel region

Introduction

Rocker shoes and insoles are used to prevent diabetic foot ulcers in persons with diabetes mellitus and loss of protective sensation, by reducing the plantar pressure in regions with high pressure values (>200 kPa) (e.g., hallux, metatarsal heads and heel). However, forefoot rocker shoes that reduce pressure in the forefoot inadvertently increase pressure in the heel. No studies focused on mitigating the negative effects on heel pressure by optimizing the heel rocker midsole, yet. Therefore, we analyze the effect of different heel rocker parameters on the heel plantar pressure.

Methods

In-shoe pressure was measured, while 10 healthy participants walked with control shoe and 10 different heel rocker settings. Peak pressure was determined in 7 heel masks, for all shoes. Generalized estimating equations was performed to test the effect of the different shoes on the peak pressure in the different heel masks.

Results

In the proximal heel, a rocker shoe with distal apex position, small rocker radius and large apex angle (100°), shows the largest significant decrease in peak pressure compared to rocker shoes with more proximally located apex positions. In the midheel and distal heel, the same rocker shoes or any other rocker shoes, analyzed in this study, do not reduce the PP more than 2 % compared to the control shoe. For the midheel and distal heel region with high pressure values (>200 kPa), rocker shoes alone are not the correct option to reduce the pressure to below 200 kPa.

Conclusion

When using rocker shoes to reduce the pressure in the forefoot, a heel rocker midsole with a distal apex position, small rocker radius and apex angle of 100°, mitigates the negative effects on proximal heel pressure. For the midheel and distal heel, other footwear options as an addition or instead of rocker shoes are needed to reduce the pressure.

Effects of Barefoot and Shod Conditions on the Kinematics and Kinetics of the Lower Extremities in Alternating Jump Rope Skipping-A One-Dimensional Statistical Parameter Mapping Study

PURPOSE

To explore the difference in the biomechanics of the lower extremity during alternating jump rope skipping (AJRS) under barefoot and shod conditions.

METHODS

Fourteen experienced AJRS participants were randomly assigned to wear jump rope shoes or be barefoot (BF) during the AJRS at a self-selected speed. The Qualisys motion capture system and Kistler force platform were used to synchronously collect the ground reaction forces and trajectory data of the hip, knee, ankle, and metatarsophalangeal (MTP) joints. One-dimensional statistical parameter mapping was used to analyze the kinematics and kinetics of the lower extremity under both conditions using paired t-tests.

RESULTS

Wearing shoes resulted in a significant decrease in the ROM (p < 0.001) and peak angular velocity (p < 0.001) of the MTP joint during the landing phase. In addition, the MTP joint power (p < 0.001) was significantly larger under shod condition at 92-100% of the landing phase. Moreover, wearing shoes reduced the peak loading rate (p = 0.002).

CONCLUSION

The findings suggest that wearing shoes during AJRS could provide better propulsion during push-off by increasing the MTP plantarflexion joint power. In addition, our results emphasize the significance of the ankle and MTP joint by controlling the ankle and MTP joint angle.

Comparative Efficacy of Vibration foam Rolling and Cold Water Immersion in Amateur Basketball Players after a Simulated Load of Basketball Game

To compare the efficacy of different recovery strategies (sitting; cold water immersion, CWI; vibration foam rolling, VFR) on the lower extremities of amateur basketball players after the simulated load of a basketball game, we assessed the power, agility, and dynamic balance before and after interventions. Ten amateur basketball players alternately underwent 12 min of sitting, 12 min of CWI at 5 °C, and 12 min of VFR. The power, agility, and dynamic balance were measured immediately post-warm-up, immediately post-game, immediately post-intervention, 1 h after interventions, and 24 h after interventions. To simulate the load of a basketball game, specific movements were designed and implemented. Jump height was measured using a Kistler force plate. Reaction time and dynamic balance score were assessed using the Pavigym agility response system and the Y balance test, respectively. The data were analyzed with a two-way repeated measures analysis of variance (ANOVA). The results showed that the vertical jump height significantly decreased after the CWI intervention compared to the CON and VFR groups (p < 0.001). At 1 h after the intervention, the vertical jump height in the CON group showed delayed recovery compared to the CWI and VFR groups (p = 0.007; p < 0.001). At 24 h after the intervention, the vertical jump height in the CWI group further increased and was significantly different from the CON and VFR groups (p < 0.001; p = 0.005). Additionally, reaction times significantly increased immediately after the CWI intervention (p = 0.004) but showed further recovery at 24 h compared to the CON group (p < 0.001). The dynamic balance score significantly rebounded after the CWI intervention compared to the CON group (p = 0.021), with further improvement at 24 h (p < 0.001). CWI initially showed negative effects, but over time, its recovery effect was superior and more long-lasting. VFR had the best immediate effect on lower limb recovery after the game.

The effect of foot posture on static balance, ankle and knee proprioception in 18-to-25-year-old female student: a cross-sectional study

BACKGROUND & PURPOSE

Afferent input from the sole affects postural stability. Cutaneous reflexes from the foot are important to posture and gait. Lower-limb afferents alone provide enough information to maintain upright stance and are critical in perceiving postural sway. Altered feedback from propreoceptive receptors alters gait and patterns of muscle activation. The position and posture of the foot and ankle may also play an important role in proprioceptive input.Therefore, the current research aims to compare static balance and ankle and knee proprioception in people with and without flexible flatfeet.

METHODOLOGY

91 female students between the ages of 18 and 25 voluntarily participated in this study, of which 24 were in the flexible flatfoot group and 67 were in the regular foot group after evaluating the longitudinal arch of the foot. The position sense of ankle and knee joints were measured using the active reconstruction test of the ankle and knee angle; Static balance was measured using the Sharpened Romberg test. Data were non-normally distributed. Accordingly, non-parametric tests were applied. The Kruskal-Wallis test was applied to compare differences between groups in variables.

RESULT

Kruskal-Wallis test showed a significant difference between two groups of flat feet and normal feet in the variables of static balance and position sense of ankle plantarflexion, ankle dorsiflexion, and knee flexion (p ≤ 0.05). A significant correlation was found between static balance and sense of ankle and knee position in the group with normal feet. The analysis of the regression line also showed that ankle and knee position sense could predict the static balance score in the regular foot group (ankle dorsiflexion position sense 17% (R² = 0.17), ankle plantarflexion position sense 17% (R² = 0.17) and knee flexion position sense 46% (R² = 0.46) explain of changes in static balance).

DISCUSSION & CONCLUSION

Flexible flatfoot soles can cause loss of balance and sense of joint position; therefore, according to this preliminary study, clinicians must be aware and should take into account this possible deficit in the management of these patients.

Applications of Finite Element Modeling in Biomechanical Analysis of Foot Arch Deformation: A Scoping Review

Excessive foot arch deformation is associated with plantar tissue overload and ligamentous injury pathologies. Finite element (FE) analysis, as an effective tool for modeling and simulation, has been utilized clinically for providing insights into arch biomechanics. This systematic scoping review aimed to summarize the current state of computational modeling techniques utilized in arch biomechanics from 2000 onwards and outline the main challenges confronting the further development of accurate models in clinical conditions. English-language searches of the electronic databases were conducted in the Web of Science, PubMed, and Scopus until July 2022. Articles that investigated arch deformation mechanisms by FE modeling were included. The methodological quality was assessed utilizing the Methodological Quality Assessment of Subject-Specific Finite Element Analysis Used in Computational Orthopedics (MQSSFE). Seventeen articles were identified in this systematic scoping review, mostly focusing on constructing models for specific pathological conditions, such as progressive collapsing foot deformity, valgus foot, and posterior tibial tendon dysfunction. However, given the complexity of the arch problem, geometrical simplifications regarding the balance between accurate detail and computational cost and assumptions made in defining modeling parameters (material properties and loading and boundary conditions) may bring challenges to the accuracy and generalizability of models applied to clinical settings. Overall, advances in computational modeling techniques have contributed to reliable foot deformation simulation and analysis in modern personalized medicine.

Assessment of the Efficiency of Measuring Foot and Ankle Edema with a 3D Portable Scanner

Background: To prospectively evaluate the reliability of a portable optical scanner compared to the water displacement technique for volumetric measurements of the foot and ankle and to compare the acquisition time associated with these two methods. Methods: Foot volume was measured in 29 healthy volunteers (58 feet, 24 females and 5 males) by a 3D scanner (UPOD-S 3D Laser Full-Foot Scanner®) and by water displacement volumetry. Measurements were performed on both feet, up to a height of 10 cm above the ground. The acquisition time for each method was evaluated. The Kolmogorov-Smirnov test, Lin's Concordance Correlation Coefficient, and a Student's t-test were performed. Results: Mean foot volume was 869.7 +/- 165.1 cm³ (3D scanner) versus 867.9 +/- 155.4 cm³ (water-displacement volumetry) (p < 10^-5). The concordance of measurements was 0.93, indicative of a high correlation between the two techniques. Volumes were 47.8 cm³ lower when using the 3D scanner versus water volumetry. After statistically correcting this underestimation, the concordance was improved (0.98, residual bias = -0.03 +/- 35.1 cm³). The mean examination time was 4.2 +/- 1.7 min (3D optical scanner) versus 11.1 +/- 2.9 min (water volumeter) (p < 10^-4). Conclusions: Ankle/foot volumetric measurements performed using this portable 3D scanner are reliable and fast and can be used in clinical practice and research.

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.

Silicone Elastomeric-Based Materials of Soft Pneumatic Actuator for Lower-Limb Rehabilitation: Finite Element Modelling and Prototype Experimental Validation

This study describes the basic design, material selection, fabrication, and evaluation of soft pneumatic actuators (SPA) for lower-limb rehabilitation compression therapy. SPAs can be a promising technology in proactive pressure delivery, with a wide range of dosages for treating venous-related diseases. However, the most effective design and material selection of SPAs for dynamic pressure delivery have not been fully explored. Therefore, a SPA chamber with two elastomeric layers was developed for this study, with single-side inflation. The 3D deformation profiles of the SPA chamber using three different elastomeric rubbers were analyzed using the finite element method (FEM). The best SPA-compliant behavior was displayed by food-grade silicone A10 Shore with a maximum deformation value of 25.34 mm. Next, the SPA chamber was fabricated using A10 Shore silicone and experimentally validated. During the simulation in FEM, the air pressure was applied on the inner wall of the chamber (i.e., the affected area). This is to ensure the applied pressure was evenly distributed in the inner wall while the outer wall of the chamber remained undeformed for all compression levels. During the inflation process, pressure will be applied to the SPA chamber, causing exerted pressure on the skin which is then measured for comparison. The simulation and experimental results show an excellent agreement of pressure transmission on the skin for the pressure range of 0–120 mmHg, as depicted in the Bland–Altman plots. The findings exhibited promising results in the development of the SPA chamber using low-cost and biocompatible food-grade silicone.