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 strike pattern,step rate,and trunk posture combined gait modifications to reduce impact loading during running

Elevated impact loading can be detrimental to runners as it has been linked to the increased risk of tibial stress fracture and plantar fasciitis. The objective of this study was to investigate the combined effects of foot strike pattern, step rate, and anterior trunk lean gait modifications on impact loading in runners. Nineteen healthy runners performed 12 separate gait modification trials involving: three foot strike patterns (rearfoot, midfoot, and forefoot strike), two step rates (natural and 10% increased), and two anterior trunk lean postures (natural and 10-degree increased flexion). Overall, forefoot strike combined with increased step rate led to the lowest impact loading rates, and rearfoot strike combined with anterior trunk lean led to the highest impact loading rates. In addition, there were interaction effects between foot strike pattern and step rate on awkwardness and effort, such that it was both more natural and easier to transition to a combined gait modification involving forefoot strike and increased step rate than to an isolated gait modification involving either forefoot strike or increased step rate. These findings could help to inform gait modifications for runners to reduce impact loading and associated injury risks.     

Foot internal stress distribution during impact in barefoot running as function of the strike pattern

The aim of the present study is to examine the impact absorption mechanism of the foot for different strike patterns (rearfoot, midfoot and forefoot) using a continuum mechanics approach. A three-dimensional finite element model of the foot was employed to estimate the stress distribution in the foot at the moment of impact during barefoot running. The effects of stress attenuating factors such as the landing angle and the surface stiffness were also analyzed. We characterized rear and forefoot plantar sole behavior in an experimental test, which allowed for refined modeling of plantar pressures for the different strike patterns. Modeling results on the internal stress distributions allow predictions of the susceptibility to injury for particular anatomical structures in the foot.     

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.     

Acute changes in foot strike pattern and cadence affect running parameters associated with tibial stress fractures

Tibial stress fractures are a common and debilitating injury that occur in distance runners. Runners may be able to decrease tibial stress fracture risk by adopting a running pattern that reduces biomechanical parameters associated with a history of tibial stress fracture. The purpose of this study was to test the hypothesis that converting to a forefoot striking pattern or increasing cadence without focusing on changing foot strike type would reduce injury risk parameters in recreational runners. Running kinematics, ground reaction forces and tibial accelerations were recorded from seventeen healthy, habitual rearfoot striking runners while running in their natural running pattern and after two acute retraining conditions: (1) converting to forefoot striking without focusing on cadence and (2) increasing cadence without focusing on foot strike. We found that converting to forefoot striking decreased two risk factors for tibial stress fracture: average and peak loading rates. Increasing cadence decreased one risk factor: peak hip adduction angle. Our results demonstrate that acute adaptation to forefoot striking reduces different injury risk parameters than acute adaptation to increased cadence and suggest that both modifications may reduce the risk of tibial stress fractures.     

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.     

Joint contact loading in forefoot and rearfoot strike patterns during running

Research concerning forefoot strike pattern (FFS) versus rearfoot strike pattern (RFS) running has focused on the ground reaction force even though internal joint contact forces are a more direct measure of the loads responsible for injury. The main purpose of this study was to determine the internal loading of the joints for each strike pattern. A secondary purpose was to determine if converted FFS and RFS runners can adequately represent habitual runners with regards to the internal joint loading. Using inverse dynamics to calculate the net joint moments and reaction forces and optimization techniques to estimate muscle forces, we determined the axial compressive loading at the ankle, knee, and hip. Subjects consisted of 15 habitual FFS and 15 habitual RFS competitive runners. Each subject ran at a preferred running velocity with their habitual strike pattern and then converted to the opposite strike pattern. Plantar flexor muscle forces and net ankle joint moments were greater in the FFS running compared to the RFS running during the first half of the stance phase. The average contact forces during this period increased by 41.7% at the ankle and 14.4% at the knee joint during FFS running. Peak ankle joint contact force was 1.5 body weights greater during FFS running (p<0.05). There was no evidence to support a difference between habitual and converted running for joint contact forces. The increased loading at the ankle joint for FFS is an area of concern for individuals considering altering their foot strike pattern.     

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.     

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.     

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.     

Foot strike patterns of runners at the 15-km point during an elite-level half marathon

There are various recommendations by many coaches regarding foot landing techniques in distance running that are meant to improve running performance and prevent injuries. Several studies have investigated the kinematic and kinetic differences between rearfoot strike (RFS), midfoot strike (MFS), and forefoot strike (FFS) patterns at foot landing and their effects on running efficiency on a treadmill and over ground conditions. However, little is known about the actual condition of the foot strike pattern during an actual road race at the elite level of competition. The purpose of the present study was to document actual foot strike patterns during a half marathon in which elite international level runners, including Olympians, compete. Four hundred fifteen runners were filmed by 2 120-Hz video cameras in the height of 0.15 m placed at the 15.0-km point and obtained sagittal foot landing and taking off images for 283 runners. Rearfoot strike was observed in 74.9% of all analyzed runners, MFS in 23.7%, and FFS in 1.4%. The percentage of MFS was higher in the faster runners group, when all runners were ranked and divided into 50 runner groups at the 15.0-km point of the competition. In the top 50, which included up to the 69th place runner in actual order who passed the 15-km point at 45 minutes, 53 second (this speed represents 5.45 m x s(-1), or 15 minutes, 17 seconds per 5 km), RFS, MFS, and FFS were 62.0, 36.0, and 2.0%, respectively. Contact time (CT) clearly increased for the slower runners, or the placement order increased (r = 0.71, p < or = 0.05). The CT for RFS + FFS for every 50 runners group significantly increased with increase of the placement order. The CT for RFS was significantly longer than MFS + FFS (200.0 +/- 21.3 vs. 183.0 +/- 16 millisecond). Apparent inversion (INV) of the foot at the foot strike was observed in 42% of all runners. The percentage of INV for MFS was higher than for RFS and FFS (62.5, 32.0, and 50%, respectively). The CT with INV for MFS + FFS was significantly shorter than the CT with and without INV for RFS. Furthermore, the CT with INV was significantly shorter than push-off time without INV for RFS. The findings of this study indicate that foot strike patterns are related to running speed. The percentage of RFS increases with the decreasing of the running speed; conversely, the percentage of MFS increases as the running speed increases. A shorter contact time and a higher frequency of inversion at the foot contact might contribute to higher running economy.     

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.     

Midtarsal locking,the windlass mechanism,and running strike pattern: A kinematic and kinetic assessment

Changes in running strike pattern affect ankle and knee mechanics, but little is known about the influence of strike pattern on the joints distal to the ankle. The purpose of this study was to explore the effects of forefoot strike (FFS) and rearfoot strike (RFS) running patterns on foot kinematics and kinetics, from the perspectives of the midtarsal locking theory and the windlass mechanism. Per the midtarsal locking theory, we hypothesized that the ankle would be more inverted in early stance when using a FFS, resulting in decreased midtarsal joint excursions and increased dynamic stiffness. Associated with a more engaged windlass mechanism, we hypothesized that a FFS would elicit increased metatarsophalangeal joint excursions and negative work in late stance. Eighteen healthy female runners ran overground with both FFS and RFS patterns. Instrumented motion capture and a validated multi-segment foot model were used to analyze midtarsal and metatarsophalangeal joint kinematics and kinetics. During early stance in FFS the ankle was more inverted, with concurrently decreased midtarsal eversion (p < 0.001) and abduction excursions (p = 0.003) but increased dorsiflexion excursion (p = 0.005). Dynamic midtarsal stiffness did not differ (p = 0.761). During late stance in FFS, metatarsophalangeal extension was increased (p = 0.009), with concurrently increased negative work (p < 0.001). In addition, there was simultaneously increased midtarsal positive work (p < 0.001), suggesting enhanced power transfer in FFS. Clear evidence for the presence of midtarsal locking was not observed in either strike pattern during running. However, the windlass mechanism appeared to be engaged to a greater extent during FFS.     

Foot joint coupling variability differences between habitual rearfoot and forefoot runners prior to and following an exhaustive run

As joint coupling variability has been associated with running-related lower extremity injury, the purpose of this study was to identify how variability within the foot may be different between forefoot (FFS) and rearfoot strike (RFS) runners. Identifying typical variability in uninjured runners may contribute to understanding of ideal coordination associated with running foot strike patterns.Fifteen FFS and 15 RFS runners performed a maximal-effort 5 km treadmill run. A 7-segment foot model identified 6 functional articulations (rearfoot, medial and lateral midfoot and forefoot, and 1st metatarsophalangeal) for analysis. Beginning and end of the run motion capture data were analyzed. Vector coding was used to calculate 6 joint couples. Standard deviations of the coupling angles were used to identify variability within subphases of stance (loading, mid-stance, terminal, and pre-swing). Mixed between-within subjects ANOVAs compared differences between the foot strikes, pre and post run.Increased variability was identified within medial foot coupling for FFS and within lateral foot coupling for RFS during loading and mid-stance. The exhaustive run increased variability during mid-stance for both groups.Interpretation. Joint coupling variability profiles for FFS and RFS runners suggest different foot regions have varying coordination needs which should be considered when comparing the strike patterns.     

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.     

Development of a patient-specific anatomical foot model from structured light scan data

The use of anatomically accurate finite element (FE) models of the human foot in research studies has increased rapidly in recent years. Uses for FE foot models include advancing knowledge of orthotic design, shoe design, ankle–foot orthoses, pathomechanics, locomotion, plantar pressure, tissue mechanics, plantar fasciitis, joint stress and surgical interventions. Similar applications but for clinical use on a per-patient basis would also be on the rise if it were not for the high costs associated with developing patient-specific anatomical foot models. High costs arise primarily from the expense and challenges of acquiring anatomical data via magnetic resonance imaging (MRI) or computed tomography (CT) and reconstructing the three-dimensional models. The proposed solution morphs detailed anatomy from skin surface geometry and anatomical landmarks of a generic foot model (developed from CT or MRI) to surface geometry and anatomical landmarks acquired from an inexpensive structured light scan of a foot. The method yields a patient-specific anatomical foot model at a fraction of the cost of standard methods. Average error for bone surfaces was 2.53 mm for the six experiments completed. Highest accuracy occurred in the mid-foot and lowest in the forefoot due to the small, irregular bones of the toes. The method must be validated in the intended application to determine if the resulting errors are acceptable.     

A Coupling Analysis of the Biomechanical Functions of Human Foot Complex during Locomotion

This study represents a functional analysis of the human foot complex based on in-vivo gait measurements, finite element (FE) modeling and biological coupling theory, with the objective of achieving a comprehensive understanding of the impact attenuation and energy absorption functions of the human foot complex. A simplified heel pad FE model comprising reticular fiber structure and fat cells was constructed based on the foot pad Magnetic Resonance (MR) images. The model was then used to investigate the foot pad behaviors under impact during locomotion. Three-dimensional (3D) gait measurement and a 3D FE foot model comprising 29 bones, 85 ligaments and the plantar soft tissues were used to investigate the foot arch and plantar fascia deformations in mid-stance phase. The heel pad simulation results show that the pad model with fat cells (coupling model) has much stronger capacity in impact attenuation and energy storage than the model without fat cells (structure model). Furthermore, the FE simulation reproduced the deformations of the foot arch structure and the plantar fascia extension observed in the gait measurements, which reinforces the postulation that the foot arch structure also plays an important role in energy absorption during locomotion. Finally, the coupling mechanism of the human foot functions in impact attenuation and energy absorption was proposed.     

Dynamic Patterns of Forces and Loading Rate in Runners with Unilateral Plantar Fasciitis: A Cross-Sectional Study

Aim/Hypothesis

The etiology of plantar fasciitis (PF) has been related to several risk factors, but the magnitude of the plantar load is the most commonly described factor. Although PF is the third most-common injury in runners, only two studies have investigated this factor in runners, and their results are still inconclusive regarding the injury stage.

Objective

Analyze and compare the plantar loads and vertical loading rate during running of runners in the acute stage of PF to those in the chronic stage of the injury in relation to healthy runners.

Methods

Forty-five runners with unilateral PF (30 acute and 15 chronic) and 30 healthy control runners were evaluated while running at 12 km/h for 40 meters wearing standardized running shoes and Pedar-X insoles. The contact area and time, maximum force, and force-time integral over the rearfoot, midfoot, and forefoot were recorded and the loading rate (20–80% of the first vertical peak) was calculated. Groups were compared by ANOVAs (p<0.05).

Results

Maximum force and force-time integral over the rearfoot and the loading rate was higher in runners with PF (acute and chronic) compared with controls (p<0.01). Runners with PF in the acute stage showed lower loading rate and maximum force over the rearfoot compared to runners in the chronic stage (p<0.01).

Conclusion

Runners with PF showed different dynamic patterns of plantar loads during running over the rearfoot area depending on the injury stage (acute or chronic). In the acute stage of PF, runners presented lower loading rate and forces over the rearfoot, possibly due to dynamic mechanisms related to pain protection of the calcaneal area.     

The efficacy of a pneumatic compression device in the treatment of plantar fasciitis

Plantar fasciitis is a common foot disorder that affects more than two million Americans each year. Conservative management of plantar fasciitis is the first line of treatment by the medical specialist who frequently encounters this common foot complaint. In this study, 50 patients with plantar fasciitis of 4 weeks duration but less than 12 weeks were randomized into two groups. One group of 25 patients used the AirHeel (Aircast, Inc.) and the second group of 25 used the 1st Step prefabricated foot insert (Wrymark, Inc). Standardized weight-bearing radiographs were obtained in order to categorize the foot type: normal arch, pes planus, or pes cavus. Patients were initially evaluated and at 12 weeks they returned for a subsequent visit. The contact area of the foot with the AirHeel and 1st Step insert were similar, an increase of 27% and 26%, respectively, over the contact area barefoot. There was a noted difference in force reduction with the two devices. The AirHeel reduced the midstance force by 20.19%, as compared to the 1st Step insert which showed a 1.03% increase in midstance force. Patients with a higher initial pain score seemed to respond better initially to the AirHeel (p = 0.015) than the 1st Step insert (p = 0.035). This study demonstrates the benefits of offloading the foot at midstance with two devices to relieve the discomfort associated with plantar fasciitis. The Aircast AirHeel is a new modality for dynamic, functional treatment of proximal plantar fasciitis.     

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.