Finite element analysis of plantar fascia under stretch-the relative contribution of windlass mechanism and Achilles tendon force

Stretching plays an important role in the treatment of plantar fasciitis. Information on the internal stresses/strains of the plantar fascia under stretch is useful in enhancing knowledge on the stretch mechanisms. Although direct measurement can monitor plantar fascia changes, it is invasive and gathers only localized information. The purpose of this paper was to construct a three-dimensional finite element model of the foot to calculate the stretch effects on plantar fascia and monitor its stress/strain distributions and concentrations. A three-dimensional foot model was developed and contained 26 bones with joint cartilages, 67 ligaments and a fan-like solid plantar fascia modeling. All tissues were idealized as linear elastic, homogeneous and isotropic whilst the plantar fascia was assigned as hyperelastic to represent its nonlinearity. The plantar fascia was monitored for its biomechanical responses under various stretch combinations: three toe dorsiflexion angles (windlass effect: 15 degrees , 30 degrees and 45 degrees ) and five Achilles tendon forces (100, 200, 300, 400 and 500N). Our results indicated that the plantar fascia strain increased as the dorsiflexion angles increased, and this phenomenon was enhanced by increasing Achilles tendon force. A stress concentration was found near the medial calcaneal tubercle, and the fascia stress was higher underneath the first foot ray and gradually decreased as it moved toward the fifth ray. The current model recreated the position of the foot when stretch is placed on the plantar fascia. The results provided a general insight into the mechanical and biomechanical aspects of the influences of windlass mechanism and Achilles tendon force on plantar fascia stress and strain distribution. These findings might have practical implications onto plantar fascia stretch approaches, and provide guidelines to its surgical release.     

Anatomic basis of plantar flap design: clinical applications

Closure of plantar defects with local rotation flaps was studied in 10 patients with 11 plantar defects. Ages ranged from 15 to 66 years, and the average defect was 3.0 X 3.6 cm. Two patients were diabetics. Etiology was variable and included trauma, tumors, and breakdown in patients with anesthetic plantar surfaces. Plantar flaps were designed superficial to the plantar fascia based on the proximal plantar subcutaneous plexus blood supply. Sensation was provided by including the medial calcaneal nerve territory within the flap and by performing a limited intraneural dissection of the medial and lateral plantar nerves. Flaps were medially based, although laterally based designs are also possible when sensation is absent. The follow-up period averaged 20.8 months. Patients with normal sensation preoperatively had full sensation postoperatively and were able to bear weight on the flap without limitation. There was minor breakdown in one patient with incomplete sensation. One patient developed a hematoma. Sensate plantar flaps can be designed superficial to the plantar fascia. These flaps are durable and allow normal weight-bearing on the reconstructed surface.     

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.     

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.     

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.     

Ultrasound elastographic assessment of plantar fascia in runners using rearfoot strike and forefoot strike

Forefoot strike is increasingly being adopted by runners because it can better attenuate impact than rearfoot strike. However, forefoot strike may overload the plantar fascia and alter the plantar fascia elasticity. This study aimed to use ultrasound elastography to investigate and compare shear wave elasticity of the plantar fascia between rearfoot strikers and forefoot strikers. A total of 35 participants (21 rearfoot strikers and 14 forefoot strikers), who were free of lower limb injuries and diseases, were recruited from a local running club. Individual foot strike patterns were identified through the measured plantar pressure during treadmill running. The B-Mode ultrasound images and shear wave elastographic images of the plantar fascia were collected from each runner. Two independent investigators reviewed the images and examined the plantar fascia qualitatively and quantitatively. The results demonstrated an overall good agreement between the investigators in the image review outcomes (ICC:0.96–0.98, κ: 0.89). There were no significant differences in the fascial thickness (p = 0.50) and hypoechogenicity on the gray-scale images (p = 0.54) between the two groups. Shear wave elastography showed that forefoot strikers exhibited reduced plantar fascia elasticity compared to rearfoot strikers (p = 0.01, Cohen’s d = 0.91). A less elastic fascial tissue was more easily strained under loading. Tissue overstrain is frequently related to the incidence of plantar fasciitis. While further study is needed for firm conclusions, runners using forefoot strike were encouraged to enhance their foot strength for better protection of the plantar fascia.     

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.     

Periostin differentially induces proliferation, contraction and apoptosis of primary Dupuytren's disease and adjacent palmar fascia cells

Dupuytren's disease, (DD), is a fibroproliferative condition of the palmar fascia in the hand, typically resulting in permanent contracture of one or more fingers. This fibromatosis is similar to scarring and other fibroses in displaying excess collagen secretion and contractile myofibroblast differentiation. In this report we expand on previous data demonstrating that POSTN mRNA, which encodes the extra-cellular matrix protein periostin, is up-regulated in Dupuytren's disease cord tissue relative to phenotypically normal palmar fascia. We demonstrate that the protein product of POSTN, periostin, is abundant in Dupuytren's disease cord tissue while little or no periostin immunoreactivity is evident in patient-matched control tissues. The relevance of periostin up-regulation in DD was assessed in primary cultures of cells derived from diseased and phenotypically unaffected palmar fascia from the same patients. These cells were grown in type-1 collagen-enriched culture conditions with or without periostin addition to more closely replicate the in vivo environment. Periostin was found to differentially regulate the apoptosis, proliferation, α smooth muscle actin expression and stressed Fibroblast Populated Collagen Lattice contraction of these cell types. We hypothesize that periostin, secreted by disease cord myofibroblasts into the extra-cellular matrix, promotes the transition of resident fibroblasts in the palmar fascia toward a myofibroblast phenotype, thereby promoting disease progression.     

Wide clinical spectrum in Zimmermann-Laband syndrome

Gingival fibromatosis can be present as an isolated form or be part of a genetic disease. The Zimmermann-Laband syndrome (ZLS) is a rare disorder inherited as an autosomal dominant fashion, clinically characterized by gingival fibromatosis, bulbous soft nose, thick floppy ears, nail dysplasia, joint hyperextensibility, hepatosplenomegaly, skeletal anomalies and occasional mental retardation. We studied a girl aged five years with clinical and radiological features of the ZLS, additionally she presented deafness not previously described in the ZLS, as only partial hearing loss was reported in some patients. The father presented some facial features suggestive of ZLS, nevertheless he did not have gingival fibromatosis or hypertrichosis. We suggest that this case supports that ZLS can be part a contiguous genes syndrome or be consequence ofa gene mutation with wide variable expression. The present report supports that ZLS has a wide clinical spectrum.     

Dermatan sulfate remodeling associated with advanced Dupuytren's contracture

Dermatan sulfate (DS) widespread as a component of extracellular matrix proteoglycans, is characterized by great bio-reactivity and remarkable structural heterogeneity due to distinct degrees of sulfation and glucuronosyl epimerization and different polymerization degrees. However, DS metabolism under various biological conditions is poorly known. Dupuytren's contracture is a benign fibromatosis leading to complex remodeling of the palmar fascia structure and properties. However, it remains unclear whether the disease affects the structure of DS, which is the major tissue glycosaminoglycan. Thus the aim of the study was to examine the structure of the total DS in Dupuytren's fascia. DS chains were extracted from 5 samples of normal fascia and 7 specimens of Dupuytren's tissue by papain digestion followed by fractionation with cetylpyridinium chloride. Then, DS structure analysis was performed comprising the evaluation of its molecular masses and sensitivity to hyaluronidase and chondroitinase B. Dupuytren's contracture is associated with significant remodeling of DS chain structure revealed by (1) a distinct profile of chain molecular masses characterized by the appearance of long size components as well as the increase in the content of small size chains; (2) a different glucuronosyl epimerization pattern connected with the enhanced content of glucuronate disaccharide blocks; (3) chain oversulfation. These structural alterations in total DS may modify the GAG interactions especially affecting collagen fibrillogenesis and growth factor availability. Thus, Dupuytren's contracture associated DS remodeling may promote the phenomena typical for advanced disease: apoptosis and reduction in cell number as well as the appearance of dense pseudotendinous collagen matrix.     

Matrix metalloproteinases and tissue inhibitors of metalloproteinases in sera and tissue of patients with Dupuytren's disease

Dupuytren's contracture is a fibroproliferative disorder characterized by progressive deposition of mature collagen fibers. In other fibrotic diseases affecting organs such as the liver, lung, heart, and skin, matrix metalloproteinases (MMPs) and their natural inhibitors, the tissue inhibitors of metalloproteinases (TIMPs), play an important role. In this study, serum concentrations of MMP-1, MMP-2, MMP-9, TIMP-1, and TIMP-2 were determined in 22 patients (five women and 17 men; average age, 67 +/- 11 years) with Dupuytren's disease using an enzyme-linked immunosorbent assay. Tissue samples were obtained for standard histological and immunohistochemical analyses. Sera and samples of palmar fascia from 20 patients (13 women and seven men; average age, 60 +/- 15 years) who had undergone hand surgery for carpal tunnel syndrome were used as the control group. Statistical analysis was performed using the Mann-Whitney test. Patients with Dupuytren's contracture presented with a TIMP-1 concentration of 437 +/- 160 ng/ml, a significantly higher TIMP-1 concentration than that seen in the control patients, who had a concentration of 321 +/- 70 ng/ml (p < 0.05). Patients with a proliferative active disease (n = 14) had a significantly higher TIMP-1 concentration (525 +/- 136 ng/ml) than patients (n = 8) with a contracture in the late involutional and residual phase (286 +/- 41 ng/ml; p < 0.05). There were no significant differences in the TIMP-2, MMP-1, MMP-2, and MMP-9 serum concentrations between patients with palmar fibromatosis and the control group. Patients with Dupuytren's disease had a significantly lower MMP-to-TIMP ratio (1.1 +/- 0.3; p < 0.05) than the control group (1.5 +/- 0.35). Patients with an active palmar fibromatosis presented a significantly (p < 0.05) reduced ratio (1 +/- 0.2) compared with those in later phases (1.4 +/- 0.3). TIMP-1 and TIMP-2 could be detected in tissue of patients with Dupuytren's contracture, with an accumulation in proliferative areas. MMPs could be detected locally in Dupuytren's tissue in a few patients, with less positive staining than for TIMPs. In the control group, there was just little or no staining for TIMPs and MMPs. The data indicate that the physiological balance between MMPs and their natural inhibitors is disturbed in patients with a proliferative active Dupuytren's disease. The decrease in the systemic MMP-to-TIMP ratio can cause increased synthesis and deposition of collagen, leading to palmar fibromatosis.     

Alterations in the extracellular matrix proteoglycan profile in Dupuytren's contracture affect the palmar fascia

Dupuytren's disease is a palmar fibromatosis associated with changes in fibroblast activity that also affect the metabolism of extracellular matrix components. In contrast to disease connected alterations in collagen and non-collagenous glycoproteins (mainly fibronectin), the metabolism of proteoglycans, being glycosaminoglycan modified glycoproteins, is poorly understood. Thus, the aim of the present study was the characterization of matrix proteoglycans (PGs) derived from normal fascia and Dupuytren's fascia. Extracted and purified PGs (particularly small PGs) were analysed for content, molecular mass, immunoreactivity and glycosaminoglycan chain structure. The matrix of normal fascia mainly contains decorin [small dermatan sulfate (DS) PG] with biglycan (another small DSPG) and large chondroitin sulfate(CS)/DSPG representing minor components. Dupuytren's disease is associated with the remodeling of matrix PG composition. The most prominent alteration is an accumulation of biglycan frequently bearing DS chains with higher molecular masses. Moreover, the amount of large CS/DSPG is increased. In contrast, decorin displays changes affecting mainly DS chain structure reflected in (i) an increase in some chain molecular masses, (ii) an enhanced content of iduronate disaccharide clusters, and (iii) oversulfation of disaccharide repeats. The PG alterations observed in Dupuytren's fascia may influence the matrix properties and contribute to disease progression.     

Changes in foot loading following plantar fasciotomy: a computer modeling study

Forward dynamic simulations of a toe-rise task were developed to explore the outcomes of plantar fasciotomy, a surgery commonly performed to relieve heel pain. The specific objectives of this study were to develop such a simulation, validate its predictions, and simulate rising on toes using a model from which the plantar fascia had been removed. Root-mean squared differences between the intact model and measurements of healthy subjects were found to be 0.009 body weights (BW) and 0.055 BW for the horizontal and vertical ground reaction forces and 7.1 mm, 11.3 mm, and 0.48 deg for the horizontal, vertical and rotational positions of the pelvis. Simulated plantar fasciotomy increased passive arch torques by 7.4%, increased metatarsal head contact forces by 18%, and resulted in greater toe flexor activity. These simulations may explain the mechanisms behind plantar fasciotomy complications when patients perform activities that require loading of the plantarflexors and the longitudinal arch.     

Evaluation of the rheumatoid arthritis susceptibility loci HLA-DRB1, PTPN22, OLIG3/TNFAIP3, STAT4 and TRAF1/C5 in an inception cohort

Introduction

The purpose of this study was to investigate the effectiveness of adalimumab in enthesitis and peripheral arthritis in patients with ankylosing spondylitis (AS).

Methods

Adults with active AS (Bath ankylosing spondylitis disease activity index [BASDAI] ≥ 4) received adalimumab 40 mg every other week with standard antirheumatic therapies in a 12-week, open-label study. Effectiveness in enthesitis was assessed using the Maastricht ankylosing spondylitis enthesitis score (MASES, 0-13) and by examining the plantar fascia in patients with enthesitis (≥ 1 inflamed enthesis) at baseline; effectiveness in peripheral arthritis was evaluated using tender and swollen joint counts (TJC, 0-46; SJC, 0-44) in patients with peripheral arthritis (≥ 1 swollen joint) at baseline. Overall effectiveness measures included Assessment of SpondyloArthritis International Society 20% response (ASAS20).

Results

Of 1,250 patients enrolled, 686 had enthesitis and 281 had peripheral arthritis. In 667 patients with MASES ≥ 1 at baseline, the median MASES was reduced from 5 at baseline to 1 at week 12. At week 12, inflammation of the plantar fascia ceased in 122 of 173 patients with inflammation at baseline. The median TJC in 281 patients with SJC ≥ 1 at baseline was reduced from 5 at baseline to 1 at week 12; the median SJC improved from 2 to 0. ASAS20 responses were achieved by 70.5% of 457 patients with no enthesitis and no arthritis; 71.0% of 512 patients with only enthesitis; 68.0% of 107 patients with only arthritis; and 66.7% of 174 patients with both.

Conclusions

Treatment with adalimumab improved enthesitis and peripheral arthritis in patients with active AS.

Trial registration

ClinicalTrials.gov NCT00478660.     

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.     

Real-time subject-specific monitoring of internal deformations and stresses in the soft tissues of the foot: a new approach in gait analysis

No technology is presently available to provide real-time information on internal deformations and stresses in plantar soft tissues of individuals during evaluation of the gait pattern. Because internal deformations and stresses in the plantar pad are critical factors in foot injuries such as diabetic foot ulceration, this severely limits evaluation of patients. To allow such real-time subject-specific analysis, we developed a hierarchal modeling system which integrates a two-dimensional gross structural model of the foot (high-order model) with local finite element (FE) models of the plantar tissue padding the calcaneus and medial metatarsal heads (low-order models). The high-order whole-foot model provides real-time analytical evaluations of the time-dependent plantar fascia tensile forces during the stance phase. These force evaluations are transferred, together with foot-shoe local reaction forces, also measured in real time (under the calcaneus, medial metatarsals and hallux), to the low-order FE models of the plantar pad, where they serve as boundary conditions for analyses of local deformations and stresses in the plantar pad. After careful verification of our custom-made FE solver and of our foot model system with respect to previous literature and against experimental results from a synthetic foot phantom, we conducted human studies in which plantar tissue loading was evaluated in real time during treadmill gait in healthy individuals (N = 4). We concluded that internal deformations and stresses in the plantar pad during gait cannot be predicted from merely measuring the foot-shoe force reactions. Internal loading of the plantar pad is constituted by a complex interaction between the anatomical structure and mechanical behavior of the foot skeleton and soft tissues, the body characteristics, the gait pattern and footwear. Real-time FE monitoring of internal deformations and stresses in the plantar pad is therefore required to identify elevated deformation/stress exposures toward utilizing it in gait laboratories to protect feet that are susceptible to injury.     

The first branch of the lateral plantar nerve and heel pain

The course and ramification pattern of the lateral plantar nerve was studied in serial sections from 4 fetal feet and in dissections from 34 adult feet with special reference to the so called first branch. This branch was found in all of the observed fetal and adult specimen. From its originating point the nerve runs immediately distally to the medial process of the calcaneal tuberosity in a lateral direction to the proximal part of the abductor digiti minimi muscle. During its course the FB gives two branches. One of them penetrates sometimes the insertion of the quadratus plantae muscle, whereas in adult feet it always sends fibres to the periosteum around the medial process of the calcaneal tuberosity and the long plantar ligament. The other innervates the flexor digitorum brevis muscle. The site of a possible entrapment is located between the abductor hallucis muscle and the medial head of the quadratus plantae muscle. There is strong indirect evidence that the nerve is of a mixed type consisting of sensory fibres for the calcaneal periosteum and the medial head of the quadratus plantae muscle. There is strong indirect evidence that the nerve is of a mixed type consisting of sensory fibres for the calcaneal periosteum and the long plantar ligament as well as motor fibres for the quadratus plantae, flexor digitorum brevis and abductor digiti minimi muscles, which may explain the characteristic pain complaints of the heel pain syndrome. The occurrence of a stiff fascia perforated by the nerve branch or a bursa around the insertion of the plantar aponeurosis as has been described by several authors and which was put forward as a possible aetiological factor could not be confirmed in our material.     

Towards patient-specific medializing calcaneal osteotomy for adult flatfoot: a finite element study

Clinically in medializing calcaneal osteotomy (MCO), foot and ankle surgeons are facing difficulties in choosing appropriate surgical parameters due to the individual differences in deformities among flatfoot patients. Traditional cadaveric studies have provided important information regarding the biomechanical effects of tendons, ligaments, and plantar fascia, but limitations have been reached when dealing with individual differences and tailoring patient-specific surgeries. Therefore, this study aimed at implementing the finite element (FE) method to investigate the effect of different MCO parameters to help foot and ankle surgeons performing patient-specific surgeries. This study constructed FE models of a flatfoot and a healthy foot based on computed tomography (CT) images. After validating the FE models with experimental measurements, differences in plantar stress were compared between two models and a criterion was established for evaluating the performance of surgical simulations. Four MCO parameters were then studied through FE simulations. Results suggested that the transverse angle, β, and translation distance, d, affected surgical performance. Therefore, special attentions may be recommended when choosing these two parameters clinically. However, the sagittal angle, α, and osteotomy position, p, were found to have less effect on the MCO performance.     

Three-dimensional finite element analysis of the foot during standing--a material sensitivity study

Information on the internal stresses/strains in the human foot and the pressure distribution at the plantar support interface under loading is useful in enhancing knowledge on the biomechanics of the ankle-foot complex. While techniques for plantar pressure measurements are well established, direct measurement of the internal stresses/strains is difficult. A three-dimensional (3D) finite element model of the human foot and ankle was developed using the actual geometry of the foot skeleton and soft tissues, which were obtained from 3D reconstruction of MR images. Except the phalanges that were fused, the interaction among the metatarsals, cuneiforms, cuboid, navicular, talus, calcaneus, tibia and fibula were defined as contact surfaces, which allow relative articulating movement. The plantar fascia and 72 major ligaments were simulated using tension-only truss elements by connecting the corresponding attachment points on the bone surfaces. The bony and ligamentous structures were embedded in a volume of soft tissues. The encapsulated soft tissue was defined as hyperelastic, while the bony and ligamentous structures were assumed to be linearly elastic. The effects of soft tissue stiffening on the stress distribution of the plantar surface and bony structures during balanced standing were investigated. Increases of soft tissue stiffness from 2 and up to 5 times the normal values were used to approximate the pathologically stiffened tissue behaviour with increasing stages of diabetic neuropathy. The results showed that a five-fold increase in soft tissue stiffness led to about 35% and 33% increase in the peak plantar pressure at the forefoot and rearfoot regions, respectively. This corresponded to about 47% decrease in the total contact area between the plantar foot and the horizontal support surface. Peak bone stress was found at the third metatarsal in all calculated cases with a minimal increase of about 7% with soft tissue stiffening.