The extensibility of the plantar fascia influences the windlass mechanism during human running

The arch of the human foot is unique among hominins as it is compliant at ground contact but sufficiently stiff to enable push-off. These behaviours are partly facilitated by the ligamentous plantar fascia whose role is central to two mechanisms. The ideal windlass mechanism assumes that the plantar fascia has a nearly constant length to directly couple toe dorsiflexion with a change in arch shape. However, the plantar fascia also stretches and then shortens throughout gait as the arch-spring stores and releases elastic energy. We aimed to understand how the extensible plantar fascia could behave as an ideal windlass when it has been shown to strain throughout gait, potentially compromising the one-to-one coupling between toe arc length and arch length. We measured foot bone motion and plantar fascia elongation using high-speed X-ray during running. We discovered that toe plantarflexion delays plantar fascia stretching at foot strike, which probably modifies the distribution of the load through other arch tissues. Through a pure windlass effect in propulsion, a quasi-isometric plantar fascia''s shortening is delayed to later in stance. The plantar fascia then shortens concurrently to the windlass mechanism, likely enhancing arch recoil at push-off.     

Deformation and stress distribution of the human foot after plantar ligaments release: A cadaveric study and finite element analysis

The majority of foot deformities are related to arch collapse or instability, especially the longitudinal arch. Although the relationship between the plantar fascia and arch height has been previously investigated, the stress distribution remains unclear. The aim of this study was to explore the role of the plantar ligaments in foot arch biomechanics. We constructed a geometrical detailed three-dimensional (3-D) finite element (FE) model of the human foot and ankle from computer tomography images. The model comprised the majority of joints in the foot as well as bone segments, major ligaments, and plantar soft tissue. Release of the plantar fascia and other ligaments was simulated to evaluate the corresponding biomechanical effects on load distribution of the bony and ligamentous structures. These intrinsic ligaments of the foot arch were sectioned to simulate different pathologic situations of injury to the plantar ligaments, and to explore bone segment displacement and stress distribution. The validity of the 3-D FE model was verified by comparing results with experimentally measured data via the displacement and von Mise stress of each bone segment. Plantar fascia release decreased arch height, but did not cause total collapse of the foot arch. The longitudinal foot arch was lost when all the four major plantar ligaments were sectioned simultaneously. Plantar fascia release was compromised by increased strain applied to the plantar ligaments and intensified stress in the midfoot and metatarsal bones. Load redistribution among the centralized metatarsal bones and focal stress relief at the calcaneal insertion were predicted. The 3-D FE model indicated that plantar fascia release may provide relief of focal stress and associated heel pain. However, these operative procedures may pose a risk to arch stability and clinically may produce dorsolateral midfoot pain. The initial strategy for treating plantar fasciitis should be non-operative.     

Parametric study of orthopedic insole of valgus foot on partial foot amputation

Orthopedic insole was important for partial foot amputation (PFA) to achieve foot balance and avoid foot deformity. The inapposite insole orthosis was thought to be one of the risk factors of reamputation for foot valgus patient, but biomechanical effects of internal tissues on valgus foot had not been clearly addressed. In this study, plantar pressure on heel and metatarsal regions of PFA was measured using F-Scan. The three-dimensional finite element (FE) model of partial foot evaluated different medial wedge angles (MWAs) (0.0°–10.0°) of orthopedic insole on valgus foot. The effect of orthopedic insole on the internal bone stress, the medial ligament tension of ankle, plantar fascia tension, and plantar pressure was investigated. Plantar pressure on medial heel region was about 2.5 times higher than that of lateral region based on the F-Scan measurements. FE-predicted results showed that the tension of medial ankle ligaments was the lowest, and the plantar pressure was redistributed around the heel, the first metatarsal, and the lateral longitudinal arch regions when MWA of orthopedic insole ranged from 7.5° to 8.0°. The plantar fascias maintained about 3.5% of the total load bearing on foot. However, the internal stresses from foot bones increased. The simulation in this study would provide the suggestion of guiding optimal design of orthopedic insole and therapeutic planning to pedorthist.     

Stress analysis of the standing foot following surgical plantar fascia release

Plantar fascia release is a surgical alternative for patients who suffer chronic heel pain due to plantar fasciitis and are unaffected by conservative treatment. A computational (finite element) model for analysis of the structural behavior of the human foot during standing was utilized to investigate the biomechanical effects of releasing the plantar fascia. The model integrates a system of five planar structures in the directions of the foot rays. It was built according to accurate geometric data of MRI, and includes linear and non-linear elements that represent bony, cartilaginous, ligamentous and fatty tissues. The model was successfully validated by comparing its resultant ground reactions with foot-ground pressure measurements and its predicted displacements with those observed in radiological tests. Simulation of plantar fascia release (partial or total) was accomplished by gradually removing parts of the fascia in the model. The results showed that total fascia release causes extensive arch deformation during standing, which is greater than normal deformation by more than 2.5mm. Tension stresses carried by the long plantar ligaments increased significantly, and may exceed the normal average stress by more than 200%. Since the contribution of the plantar fascia to the foot's load-bearing ability is of major importance, its release must be very carefully considered, and the present model may be used to help surgeons decide upon the desired degree of release.     

Muscle-driven finite element simulation of human foot movements

This paper describes a finite element scheme for realistic muscle-driven simulation of human foot movements. The scheme is used to simulate human ankle plantar flexion. A three-dimensional anatomically detailed finite element model of human foot and lower leg is developed and the idea of generating natural foot movement based entirely on the contraction of the plantar flexor muscles is used. The bones, ligaments, articular cartilage, muscles, tendons, as well as the rest soft tissues of human foot and lower leg are included in the model. A realistic three-dimensional continuum constitutive model that describes the biomechanical behaviour of muscles and tendons is used. Both the active and passive properties of muscle tissue are accounted for. The materials for bones and ligaments are considered as homogeneous, isotropic and linearly elastic, whereas the articular cartilage and the rest soft tissues (mainly fat) are defined as hyperelastic materials. The model is used to estimate muscle tissue deformations as well as stresses and strains that develop in the lower leg muscles during plantar flexion of the ankle. Stresses and strains that develop in Achilles tendon during such a movement are also investigated.     

Muscle-driven finite element simulation of human foot movements

This paper describes a finite element scheme for realistic muscle-driven simulation of human foot movements. The scheme is used to simulate human ankle plantar flexion. A three-dimensional anatomically detailed finite element model of human foot and lower leg is developed and the idea of generating natural foot movement based entirely on the contraction of the plantar flexor muscles is used. The bones, ligaments, articular cartilage, muscles, tendons, as well as the rest soft tissues of human foot and lower leg are included in the model. A realistic three-dimensional continuum constitutive model that describes the biomechanical behaviour of muscles and tendons is used. Both the active and passive properties of muscle tissue are accounted for. The materials for bones and ligaments are considered as homogeneous, isotropic and linearly elastic, whereas the articular cartilage and the rest soft tissues (mainly fat) are defined as hyperelastic materials. The model is used to estimate muscle tissue deformations as well as stresses and strains that develop in the lower leg muscles during plantar flexion of the ankle. Stresses and strains that develop in Achilles tendon during such a movement are also investigated.     

Analysis of the main passive soft tissues associated with adult acquired flatfoot deformity development: A computational modeling approach

Adult acquired flatfoot deformity (AAFD) is a pathology with a wide range of treatment options. Physicians decide the best treatment based on their experience, so the process is entirely subjective. A better understanding of soft tissue stress and its contribution in supporting the plantar arch could help to guide the clinical decision. Traditional experimental trials cannot consistently evaluate the contribution of each tissue. Therefore, in this research a 3-Dimensional FE foot model was reconstructed from a normal patient in order to measure the stress of the passive stabilizers of the arch, and its variation in different scenarios related with intermediate stages of AAFD development. All bones, the plantar fascia (PF), cartilages, plantar ligaments and the spring ligament (SL) were included, respecting their anatomical distribution and biomechanical characteristics. An AAFD evaluation scenario was simulated. The relative contribution of each tissue was obtained comparing each result with a normal case. The results show that PF is the main tissue that prevents the arch elongation, while SL mainly reduces the foot pronation. Long and short plantar ligaments play a secondary role in this process. The stress increment on both PF and SL when one of two fails suggests that these tissues complement each other. These findings support the theory that regards the tibialis posterior tendon as a secondary actor in the arch maintenance, compared with the PF and the SL, because this tendon is overstretched by the hindfoot pronation around the talonavicular joint. This approach could help to improve the understanding of AAFD.     

Foot arch deformation and plantar fascia loading during running with rearfoot strike and forefoot strike: A dynamic finite element analysis

Forefoot strike becomes popular among runners because it facilitates better impact attenuation. However, forefoot strike may overload the plantar fascia and impose risk of plantar fasciitis. This study aimed to examine and compare the foot arch deformation and plantar fascia tension between different foot strike techniques in running using a computational modelling approach. A three-dimensional finite element foot model was reconstructed from the MRI of a healthy runner. The foot model included twenty bones, bulk soft tissue, ligaments, tendons, and plantar fascia. The time-series data of segmental kinematics, foot muscle force, and ankle joint reaction force were derived from a musculoskeletal model of the same participant based on the motion capture analysis and input as the boundary conditions for the finite element analysis. Rearfoot strike and forefoot strike running were simulated using a dynamic explicit solver. The results showed that, compared to rearfoot strike, forefoot strike reduced the foot arch height by 9.12% and increased the medial longitudinal arch angle by 2.06%. Forefoot strike also increased the plantar connective tissues stress by 18.28–200.11% and increased the plantar fascia tensile force by 18.71–109.10%. Although it is currently difficult to estimate the threshold value of stress or force that results in injury, forefoot strike runners appeared to be more vulnerable to plantar fasciitis.     

Effect of muscle contraction levels on the force-length relationship of the human Achilles tendon during lengthening of the triceps surae muscle-tendon unit

Findings from animal experiments are sometimes contradictory to the idea that the tendon structure is a simple elastic spring in series with muscle fibers, and suggest influence of muscle contraction on the tendon mechanical properties. The purpose of the present study was to investigate the influence of muscle contraction levels on the force-length relationship of the human Achilles tendon during lengthening of the triceps surae muscle-tendon unit. For seven subjects, ankle dorsiflexion was performed without (passive condition) and with contraction of plantar flexor muscles (eccentric conditions, at 3 contraction levels) on an isokinetic dynamometer. Deformation of the Achilles tendon during each trial was measured using ultrasonography. The Achilles tendon force corresponding to the tendon elongation of 10mm in the passive condition was significantly smaller than those in the eccentric conditions (p<0.05 or p<0.01). Within the eccentric conditions, the Achilles tendon force corresponding to the tendon elongation of 10mm was significantly greater in the maximal contraction level than those in submaximal eccentric conditions (p<0.05 or p<0.01). In addition, the tendon stiffness was greater in higher contraction levels (p<0.05 or p<0.01). Present results suggest that the human tendon structure is not a simple elastic spring in series with muscle fibers.     

Finite element modelling of radial shock wave therapy for chronic plantar fasciitis

Therapeutic use of high-amplitude pressure waves, or shock wave therapy (SWT), is emerging as a popular method for treating musculoskeletal disorders. However, the mechanism(s) through which this technique promotes healing are unclear. Finite element models of a shock wave source and the foot were constructed to gain a better understanding of the mechanical stimuli that SWT produces in the context of plantar fasciitis treatment. The model of the shock wave source was based on the geometry of an actual radial shock wave device, in which pressure waves are generated through the collision of two metallic objects: a projectile and an applicator. The foot model was based on the geometry reconstructed from magnetic resonance images of a volunteer and it comprised bones, cartilage, soft tissue, plantar fascia, and Achilles tendon. Dynamic simulations were conducted of a single and of two successive shock wave pulses administered to the foot. The collision between the projectile and the applicator resulted in a stress wave in the applicator. This wave was transmitted into the soft tissue in the form of compression–rarefaction pressure waves with an amplitude of the order of several MPa. The negative pressure at the plantar fascia reached values of over 1.5 MPa, which could be sufficient to generate cavitation in the tissue. The results also show that multiple shock wave pulses may have a cumulative effect in terms of strain energy accumulation in the foot.     

Force output during and following active stretches of rat plantar flexor muscles: effect of velocity of ankle rotation

During the development of force deficits by repeated stretches, velocity-sensitive changes in the extra force produced during and after subsequent stretching has not been studied. In the present study, repeated dorsiflexion of the foot of rats with maximally contracting plantar flexor muscles was performed at two angular velocities [0.87 (slow muscle stretch) and 10.47rads(-1) (fast muscle stretch)] to examine the active force of the muscles during and following dorsiflexion. Dorsiflexion was performed 30 times with a rest period of 3min between the stretches to minimize muscle fatigue. The ability of rat plantar flexor muscles to produce additional force during the stretch was not velocity sensitive. In contrast, repeated dorsiflexion with fast muscle stretches, but not with slow muscle stretches, resulted in an increase in the force decay with time following the stretches (i.e. increased stress relaxation), as indicated by a change in the time constant of force decay during stress relaxation. Apparently, the stress-relaxation of rat plantar flexor muscles is sensitive to angular velocity of ankle movements; repeated fast, but not slow dorsiflexion, alters the stress relaxation process of active skeletal muscles exposed to stretches which create a force deficit. The change in time constant of force decay during stress relaxation in response to a series of repeated stretches might provide information on the sarcomere length distribution in skeletal muscles.     

Effect of stretching training on the viscoelastic properties of human tendon structures in vivo.

The purpose of this study was to examine whether stretching training altered the viscoelastic properties of human tendon structures in vivo. Eight men performed the stretching training for 3 wk. Before and after the stretching training, the elongation of the tendon and aponeurosis of medial gastrocnemius muscle was directly measured by ultrasonography while the subjects performed ramp isometric plantar flexion up to the voluntary maximum, followed by a ramp relaxation. The relationship between the estimated muscle force (Fm) and tendon elongation (L) during the ascending phase was fitted to a linear regression, the slope of which was defined as stiffness of tendon structures. The percentage of the area within the Fm-L loop to the area beneath the curve during ascending phase was calculated as an index representing hysteresis. To assess the flexibility, the passive torque of the plantar flexor muscles was measured during the passive stretch from 0 degrees (anatomic position) to 25 degrees of dorsiflexion with a constant velocity of 5 degrees/s. The slope of the linear portion of the passive torque-angle curve during stretching was defined as flexibility index. Flexibility index decreased significantly after stretching training (-13.4 +/- 4.6%). On the other hand, the stretching training produced no significant change in stiffness but significantly decreased hysteresis from 19.9 +/- 11.7 to 12.5 +/- 9.5%. The present results suggested that stretching training affected the viscosity of tendon structures but not the elasticity.     

Mechanical stress redistribution in the calcaneus after autologous bone harvesting

The calcaneus is a desirable site for harvesting autologous bone for use in foot surgery. However, fracture of the calcaneus is a serious complication associated with bone harvesting from this site. Currently it is unknown how much bone may be safely harvested from the calcaneus without inducing a fracture. The purpose of this study was to investigate the effect of progressive bone removal from the calcaneus onto the mechanical stress redistribution of the foot, and therefore on the increase in fracture risk. Different loads were applied on the talus to evaluate the calcaneus stress distribution at different situations. Because of the potential increase in mechanical stress in the calcaneus, secondary to contraction of the Achilles tendon, we also evaluated the mechanical behavior properties of the foot with increasing traction force in the Achilles tendon. A three-dimensional (3D) finite element (FE) model developed from CT images obtained from a healthy individual was used to compute displacement, tension and compression stresses in six situations, including intact foot, and five depth of the bone block removed, with a maximum depth of 7.5 mm. The results from these simulations indicated that when the maximum load was applied at the Achilles tendon, the tension stress increased from 42.16 MPa in the intact foot to 86.28 MPa with maximum bone harvesting. Furthermore, as the volume of bone extracted from the calcaneus increases, there is a redistribution of stresses that differs significantly from the intact foot. In fact, although the maximum stress was not significantly affected by increasing the volume of bone harvested-except when increasing the Achilles tendon force-, stresses did increase in areas of the calcaneus is vulnerable to injury, leading to an increase in fracture risk.     

The contraction-induced increase in Achilles tendon moment arm: A three-dimensional study

The present study aimed to re-examine the influence of the isometric plantarflexors contraction on the Achilles tendon moment arm (ATMA) and the factors influencing the ATMA in three-dimensions. A series of coronal magnetic resonance images of the right ankle were recorded at foot positions of 10° of dorsiflexion, neutral position, and 10° of plantarflexion for the rest condition and the plantarflexors contraction condition at 30% maximal voluntary effort. The shortest distance between the talocrural joint axis and the line of action of the Achilles tendon force projected to the orthogonal plane of the talocrural joint axis was determined as the ATMA. The ATMA determined in the contraction condition was significantly greater by 8 mm than that determined in the rest condition. The talocrural joint axis was displaced anteriorly by 3 mm and distally by 2 mm due to the muscle contraction. As the same time, the line of action of the Achilles tendon force was displaced posteriorly by 5 mm and medially by 2 mm. These linear displacements of the talocrural joint axis and the line of action of the Achilles tendon force accounted for the difference in the ATMAs between the two conditions by 35.9 and 62.4%, respectively. These angular displacements accounted for the total of 0.4% increase in the ATMA. These results confirm the previous findings reported in two-dimensional studies and found that the linear displacement of the line of action of the Achilles tendon force is the primary source of the contraction-induced increase in the ATMA.     

Finite element modeling of the first ray of the foot: a tool for the design of interventions

Disorders of the first ray of the foot (defined as the hard and soft tissues of the first metatarsal, the sesamoids, and the phalanges of the great toe) are common, and therapeutic interventions to address these problems range from alterations in footwear to orthopedic surgery. Experimental verification of these procedures is often lacking, and thus, a computational modeling approach could provide a means to explore different interventional strategies. A three-dimensional finite element model of the first ray was developed for this purpose. A hexahedral mesh was constructed from magnetic resonance images of the right foot of a male subject. The soft tissue was assumed to be incompressible and hyperelastic, and the bones were modeled as rigid. Contact with friction between the foot and the floor or footwear was defined, and forces were applied to the base of the first metatarsal. Vertical force was extracted from experimental data, and a posterior force of 0.18 times the vertical force was assumed to represent loading at peak forefoot force in the late-stance phase of walking. The orientation of the model and joint configuration at that instant were obtained by minimizing the difference between model predicted and experimentally measured barefoot plantar pressures. The model were then oriented in a series of postures representative of push-off, and forces and joint moments were decreased to zero simultaneously. The pressure distribution underneath the first ray was obtained for each posture to illustrate changes under three case studies representing hallux limitus, surgical arthrodesis of the first ray, and a footwear intervention. Hallux limitus simulations showed that restriction of metatarsophalangeal joint dorsiflexion was directly related to increase and early occurrence of hallux pressures with severe immobility increasing the hallux pressures by as much as 223%. Modeling arthrodesis illustrated elevated hallux pressures when compared to barefoot and was dependent on fixation angles. One degree change in dorsiflexion and valgus fixation angles introduced approximate changes in peak hallux pressure by 95 and 22 kPa, respectively. Footwear simulations using flat insoles showed that using the given set of materials, reductions of at least 18% and 43% under metatarsal head and hallux, respectively, were possible.     

Dynamic creep and pre-conditioning of the Achilles tendon in-vivo

Warm-up exercises are often advocated prior to strenuous exercise, but the warm-up duration and effect on muscle–tendon behavior are not well defined. The gastrocnemius–Achilles tendon complexes of 18 subjects were studied to quantify the dynamic creep response of the Achilles tendon in-vivo and the warm-up dose required for the Achilles tendon to achieve steady-state behavior. A custom testing chamber was used to determine each subject's maximum voluntary contraction (MVC) during an isometric ankle plantar flexion effort. The subject's right knee and ankle were immobilized for one hour. Subjects then performed over seven minutes of cyclic isometric ankle plantar flexion efforts equal to 25–35% of their MVC at a frequency of 0.75 Hz. Ankle plantar flexion effort and images from dual ultrasound probes located over the gastrocnemius muscle–Achilles tendon and the calcaneus–Achilles tendon junction were acquired for eight seconds at the start of each sequential minute of the activity. Ultrasound images were analyzed to quantify the average relative Achilles tendon strain at 25% MVC force (ε_25%MVC) for each minute. The ε_25%MVC increased from 0.3% at the start of activity to 3.3% after seven minutes, giving a total dynamic creep of ~3.0%. The ε_25%MVC increased by more than 0.56% per minute for the first five minutes and increased by less than 0.13% per minute thereafter. Therefore, following a period of inactivity, a low intensity warm-up lasting at least six minutes or producing 270 loading cycles is required for an Achilles tendon to reach a relatively steady-state behavior.     

Quantifying rearfoot-forefoot coordination in human walking

A method is proposed to facilitate the quantification and interpretation of inter-joint/-segment coordination. This technique is illustrated using rearfoot-forefoot kinematic data. We expand existing vector coding techniques and introduce a set of operational terms through which the coordinative patterns between the rearfoot segment and the forefoot segment are summarized: in-phase, anti-phase, rearfoot phase and forefoot phase. The literature on foot mechanics has characterized the stable foot at pushoff by a decreasing medial longitudinal arch angle in the sagittal plane, which is accompanied by forefoot pronation and concurrent rearfoot supination-in other words, anti-phase motion. Nine skin markers were placed on the rearfoot and forefoot segments according to a multi-segment foot model. Three healthy subjects performed standing calibration and walking trials (1.35ms(-1)), while a three-dimensional motion capture system acquired their kinematics. Rearfoot-forefoot joint angles were derived and the arch angle was inferred from the sagittal plane. Coupling angles of rearfoot and forefoot segments were derived and categorized into one of the four coordination patterns. Arch kinematics were consistent with the literature; in stance, the arch angle reached peak dorsiflexion, and then decreased rapidly. However, anti-phase coordination was not the predominant pattern during mid- or late stance. These preliminary data suggest that the coordinative interactions between the rearfoot and the forefoot are more complicated than previously described. The technique offers a new perspective on coordination and may provide insight into deformations of underlying tissues, such as the plantar fascia.     

Biomechanical analysis of foot with different foot arch heights: a finite element analysis

Clinically, different foot arch heights are associated with different tissue injuries to the foot. To investigate the possible factors contributing to the difference in foot arch heights, previous studies have mostly measured foot pressure in either low-arched or high-arched feet. However, little information exists on stress variation inside the foot with different arch heights. Therefore, this study aimed to implement the finite element (FE) method to analyse the influence of different foot arches. This study established a 3D foot FE model using software ANSYS 11.0. After validating the FE model, this study created low-arched, high-arched and normal-arched foot FE models. The FE analysis found that both the stress and strain on the plantar fascia and metatarsal were higher in the high-arched foot, whereas the stress and strain on the calcaneous, navicular and cuboid were higher in low-arched foot. Additionally, forefoot pressure was increased with an increase in arch height.     

Elastic properties of human Achilles tendon are correlated to muscle strength.

The purpose of this study was to investigate whether the mechanical properties of the Achilles tendon were correlated to muscle strength in the triceps surae in humans. Twenty-four men and twelve women exerted maximal voluntary isometric plantar flexion (MVIP) torque. The elongation (DeltaX) and strain of the Achilles tendon (epsilon), the proximal part of which is the composite of the gastrocnemius tendon and the soleus aponeurosis, at MVIP were determined from the displacement of the distal myotendinous junction of the medial gastrocnemius using ultrasonography. The Achilles tendon force at MVIP (F) was calculated from the MVIP torque and the Achilles tendon moment arm. There were no significant differences in either the F-DeltaX or F-epsilon relationships between men and women. DeltaX and epsilon were 9.8 +/- 2.6 mm and 5.3 +/- 1.6%, respectively, and were positively correlated to F (r = 0.39, P < 0.05; r = 0.39, P < 0.05), which meant that subjects with greater muscle strength could store more elastic energy in the tendon. The regression y-intercepts for the F-DeltaX (P < 0.01) and F-epsilon (P < 0.05) relationship were significantly positive. These results might indicate that the Achilles tendon was stiffer in subjects with greater muscle strength, which may play a role in reducing the probability of tendon strain injuries. It was suggested that the Achilles tendon of subjects with greater muscle strength did not impair the potential for storing elastic energy in tendons and may be able to deliver the greater force supplied from a stronger muscle more efficiently. Furthermore, the difference in the Achilles tendon mechanical properties between men and women seemed to be correlated to the difference in muscle strength rather than gender.