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.     

Pre-impact lower extremity posture and brake pedal force predict foot and ankle forces during an automobile collision

BACKGROUND: The purpose of this study was to determine how a driver's foot and ankle forces during a frontal vehicle collision depend on initial lower extremity posture and brake pedal force. METHOD OF APPROACH: A 2D musculoskeletal model with seven segments and six right-side muscle groups was used. A simulation of a three-second braking task found 3647 sets of muscle activation levels that resulted in stable braking postures with realistic pedal force. These activation patterns were then used in impact simulations where vehicle deceleration was applied and driver movements and foot and ankle forces were simulated. Peak rearfoot ground reaction force (F(RF)), peak Achilles tendon force (FAT), peak calcaneal force (F(CF)) and peak ankle joint force (F(AJ)) were calculated. RESULTS: Peak forces during the impact simulation were 476 +/- 687 N (F(RF)), 2934 +/- 944 N (F(CF)) and 2449 +/- 918 N (F(AJ)). Many simulations resulted in force levels that could cause fractures. Multivariate quadratic regression determined that the pre-impact brake pedal force (PF), knee angle (KA) and heel distance (HD) explained 72% of the variance in peak FRF, 62% in peak F(CF) and 73% in peak F(AJ). CONCLUSIONS: Foot and ankle forces during a collision depend on initial posture and pedal force. Braking postures with increased knee flexion, while keeping the seat position fixed, are associated with higher foot and ankle forces during a collision.     

Roentgeno-functional study of the talocrural joint in the long-term follow-up of a fracture of the crural bone

On the basis of the results of clinicoroentgenological and tensographic investigations of 119 patients after traumas of the crural bones and ankle joint (2-36 yrs. ago) the authors showed the importance of roentgenofunctional investigation of the ankle joint, a patient being in a standing position with maximum flexion and extension of the foot. A specially designed footing was proposed. Of 77 patients after intra-articular fractures of the ankle bones various disorders in articular proportions, undetectable on routine roentgenography, were diagnosed in 29 by functional roentgenography. Articular changes on roentgenofunctional investigation were revealed in one patient only out of 42 patients with extra-articular fractures of the crural bones. Tensography showed disorders of foot biomechanics in all patients with subluxations in the ankle.     

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.     

Mechanically simulated muscle forces strongly stabilize intact and injured upper cervical spine specimens

Although muscles are assumed to be capable of stabilizing the spinal column in vivo, they have only rarely been simulated in vitro. Their effect might be of particular importance in unstable segments. The present study therefore tests the hypothesis that mechanically simulated muscle forces stabilize intact and injured cervical spine specimens. In the first step, six human occipito-cervical spine specimens were loaded intact in a spine tester with pure moments in lateral bending (+/- 1.5 N m), flexion-extension (+/- 1.5 N m) and axial rotation (+/- 0.5 N m). In the second step, identical flexibility tests were carried out during constant traction of three mechanically simulated muscle pairs: splenius capitits (5 N), semispinalis capitis (5 N) and longus colli (15 N). Both steps were repeated after unilateral and bilateral transection of the alar ligaments. The muscle forces strongly stabilized C0-C2 in all loading and injury states. This was most obvious in axial rotation, where a reduction of range of motion (ROM) and neutral zone to <50% (without muscles=100%) was observed. With increasing injury the normalized ROM (intact condition=100%) increased with and without muscles approximately to the same extend. With bilateral injury this increase was 125-132% in lateral bending, 112%-119% in flexion-extension and 103-116% in axial rotation. Mechanically simulated cervical spine muscles strongly stabilized intact and injured cervical spine specimens. Nevertheless, it could be shown that in vitro flexibility tests without muscle force simulation do not necessarily lead to an overestimation of spinal instability if the results are normalized to the intact state.     

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.     

Comparative study of plantar pressure during step exercise in different floor conditions

Mechanical load has been estimated during step exercise based on ground reaction force (GRF) obtained by force platforms. It is not yet accurately known whether these measures reflect foot contact forces once the latter depend on footwear and are potentially modified by the compliant properties of the step bench. The aim of the study was to compare maximal and mean plantar pressure (PP), and maximal GRF obtained by pressure insoles after performing seven movements both over two metal force platforms and over the step bench. Fifteen step-experienced females performed the movements at the cadences of 130 and 140 beats per minute. PP and GRF (estimated from PP) obtained for each floor condition were compared. Maximal PP ranged from 29.27 +/- 9.94 to 47.07 +/- 12.88 N/cm2 as for metal platforms, and from 28.20 +/- 9.32 to 43.00 +/- 13.80 N/cm2 as for the step bench. Mean PP ranged from 11.09 +/- 1.62 to 14.32 +/- 2.06 N/cm2 (platforms) and from 10.71 +/- 1.54 to 14.22 +/- 1.77 N/cm2 (step bench). GRF (normalized body weight) ranged from 1.43 +/- 0.14 to 2.41 +/- 0.24 BW (platforms) and from 1.38 +/- 0.14 to 2.36 +/- 0.19 BW (step bench). No significant statistical differences were obtained for most of the comparisons between the two conditions tested. The results suggest that metal force platform surfaces are suitable to assess mechanical load during this physical activity. The forces applied to the foot are similar to the softer step bench and the hard force platform surface. This may reflect the ability of the performers to adapt their movement patterns to normalize the impact forces in different floor conditions.     

Effect of upper and lower extremity control strategies on predicted injury risk during simulated forward falls: a study in healthy young adults

Fall-related wrist fractures are common at any age. We used a seven-link, sagittally symmetric, biomechanical model to test the hypothesis that systematically alterations in the configuration of the body during a forward fall from standing height can significantly influence the impact force on the wrists. Movement of each joint was accomplished by a pair of agonist and antagonist joint muscle torque actuators with assigned torque-angle, torque-velocity, and neuromuscular latency properties. Proportional-derivative joint controllers were used to achieve desired target body segment configurations in the pre- andor postground contact phases of the fall. Outcome measures included wrist impact forces and whole-body kinetic energy at impact in the best, and worst, case impact injury risk scenarios. The results showed that peak wrist impact force ranged from less than 1 kN to more than 2.5 kN, reflecting a fourfold difference in whole-body kinetic energy at impact (from less than 40 J to more than 160 J) over the range of precontact hip and knee joint angles used at impact. A reduction in the whole-body kinetic energy at impact was primarily associated with increasing negative work associated with hip flexion. Altering upper extremity configuration prior to impact significantly reduced the peak wrist impact force by up to 58% (from 919 N to 2212 N). Increased peak wrist impact forces associated greater shoulder flexion and less elbow flexion. Increasing postcontact arm retraction can reduce the peak wrist impact force by 28% (from 1491 N to 1078 N), but postcontact hip and knee rotations had a relatively small effect on the peak wrist impact force (8% reduction; from 1411 N to 1303 N). In summary, the choice of the joint control strategy during a forward fall can significantly affect the risk of wrist injury. The most effective strategy was to increase the negative work during hip flexion in order to dissipate kinetic energy thereby reducing the loss in potential energy prior to first impact. Extended hip or elbow configurations should be avoided in order to reduce forearm impact forces.     

Articulated external fixation of the ankle: minimizing motion resistance by accurate axis alignment

This study describes how an optimal single hinge axis position can be established for the application of articulated external fixation to the ankle joint. By deliberately introducing various amounts of relative mal-alignment between the optimal talocrural joint axis and the actual fixator hinge axis, it was possible to measure the corresponding amounts of additional resistance to joint motion. In a cadaveric study of six ankle specimens, we determined the instant axis of rotation of the talocrural joint from 3-D kinematic data. acquired by an electromagnetic motion tracking system. For each specimen, an optimal fixator hinge position was calculated from these motion data. Compared to the intact natural joint, aligning the fixator along the optimized axis position caused a moderate increase in energy (0.14 J) needed to rotate the ankle through a prescribed plantar/dorsiflexion range. However, malpositioning the hinge by 10 mm caused more than five times that amount of increase in motion resistance. While articulated external fixation with limited internal fixation can establish a favorable environment for the repair of severe injuries such as tibial pilon fractures, the large additional resistance to motion accompanying a malpositioned fixator axis suggests the development of untoward intra-articular forces that could act to disturb fragment alignment.     

The effect of static muscle forces on the fracture strength of the intact distal radius in vitro in response to simulated forward fall impacts

The distal radius fracture (DRF) is a particularly dominant injury of the wrist, commonly resulting from a forward fall on an outstretched hand. In an attempt to reduce the prevalence, costs, and potential long-term pain/deformities associated with this injury, in vivo and in vitro investigations have sought to classify the kinematics and kinetics of DRFs. In vivo forward fall work has identified a preparatory muscle contraction that occurs in the upper extremity prior to peak impact force. The present investigation constitutes the first attempt to systematically determine the effect of static muscle forces on the fracture threshold of the distal radius in vitro. Paired human cadaveric forearm specimens were divided into two groups, one that had no muscle forces applied (i.e., right arms) and the other that had muscle forces applied to ECU, ECRL, FCU and FCR (i.e., left arms), with magnitudes based on peak muscle forces and in vivo lower bound forward fall activation patterns. The specimens were secured in a custom-built pneumatic impact loading device and subjected to incremental impacts at pre-fracture (25 J) and fracture (150  J) levels. Similar fracture forces (6565 (866) N and 8665 (5133) N), impulses (47 (6) Ns and 57 (30) Ns), and energies (152 (38) J and 144 (45) J) were observed for both groups of specimens (p>0.05). Accordingly, it is suggested that, at the magnitudes presently simulated, muscle forces have little effect on the way the distal radius responds to forward fall initiated impact loading.     

An MRI-compatible foot-loading device for assessment of internal tissue deformation

It is well known that mechanical forces acting within the soft tissues of the foot can contribute to the formation of neuropathic ulcers in people with diabetes. Presently, only surface measurements of plantar pressure are used clinically to estimate risk status due to mechanical loading. It is currently not known how surface measurements relate to the three-dimensional (3-D) internal stress/strain state of the foot. This article describes the development of a foot-loading device that allows for the direct observation of the internal deformation of foot tissues under known forces. Ground reaction forces and plantar pressure distributions during normal walking were measured in ten healthy young adults. One instant in the gait cycle, when pressure under the metatarsal heads reached a peak, was extracted for simulation in an MR imager. T1-weighted 3-D gradient echo MRI sets were collected as the simulated walking ground reaction force was incrementally applied to the foot by the novel foot-loading device. The sub-metatarsal head soft-tissue thickness decreased rapidly at first and then reached a plateau. Peak plantar pressure measurements collected within the loading device (161+/-75kPa) were lower in magnitude and less focal than pressures measured during walking (492+/-91kPa). This finding implies that although the device successfully applied full peak walking ground reaction forces to the foot, they were not distributed in the same manner as during walking. Although not representative of gait, the data collected from this in vivo mechanical test are suitable for determination of foot tissue material properties or, when combined with finite element modeling, to examine the relationship between surface loading and internal stress.     

Scaling of plantar pressures in mammals

The interaction of the limbs with the substrate can teach a lot about an animal's gait mechanics. Unlike ground-reaction forces, plantar pressure distributions are rarely studied in animals, but they may provide more detailed information about the loading patterns and locomotor function of specific anatomical structures. With this study, we aim to describe pressures for a large and diverse sample of mammalian species, focusing on scaling effects. We collected dynamic plantar pressure distributions during voluntary walking in 28 mammal species. A dynamic classification of foot use was made, which distinguished between plantiportal, digitiportal and unguliportal animals. Analysis focused on scaling effects of peak pressures, peak forces and foot contact areas. Peak pressure for the complete mammal sample was found to scale to (mass)^1/2, higher than predicted assuming geometric similarity, and we found no difference between the different types of foot use. Only the scaling of peak force is dependent on the dynamic foot use. We conclude that plantar peak pressure rises faster with mass than expected, regardless of the type of foot use, and scales higher than in limb bones. These results might explain some anatomical and behavioural adaptations in graviportal animals.     

Natural gaits of the non-pathological flat foot and high-arched foot

There has been a controversy as to whether or not the non-pathological flat foot and high-arched foot have an effect on human walking activities. The 3D foot scanning system was employed to obtain static footprints from subjects adopting a half-weight-bearing stance. Based upon their footprints, the subjects were divided into two groups: the flat-footed and the high-arched. The plantar pressure measurement system was used to measure and record the subjects' successive natural gaits. Two indices were proposed: distribution of vertical ground reaction force (VGRF) of plantar and the rate of change of footprint areas. Using these two indices to compare the natural gaits of the two subject groups, we found that (1) in stance phase, there is a significant difference (p<0.01) in the distributions of VGRF of plantar; (2) in a stride cycle, there is also a significant difference (p<0.01) in the rate of change of footprint area. Our analysis suggests that when walking, the VGRF of the plantar brings greater muscle tension to the flat-footed while a smaller rate of change of footprint area brings greater stability to the high-arched.     

Multi-segment foot kinematics and ground reaction forces during gait of individuals with plantar fasciitis

Background

Clinically, plantar fasciitis (PF) is believed to be a result and/or prolonged by overpronation and excessive loading, but there is little biomechanical data to support this assertion. The purpose of this study was to determine the differences between healthy individuals and those with PF in (1) rearfoot motion, (2) medial forefoot motion, (3) first metatarsal phalangeal joint (FMPJ) motion, and (4) ground reaction forces (GRF).

Methods

We recruited healthy (n=22) and chronic PF individuals (n=22, symptomatic over three months) of similar age, height, weight, and foot shape (p>0.05). Retro-reflective skin markers were fixed according to a multi-segment foot and shank model. Ground reaction forces and three dimensional kinematics of the shank, rearfoot, medial forefoot, and hallux segment were captured as individuals walked at 1.35 ms⁻¹.

Results

Despite similarities in foot anthropometrics, when compared to healthy individuals, individuals with PF exhibited significantly (p<0.05) (1) greater total rearfoot eversion, (2) greater forefoot plantar flexion at initial contact, (3) greater total sagittal plane forefoot motion, (4) greater maximum FMPJ dorsiflexion, and (5) decreased vertical GRF during propulsion.

Conclusion

These data suggest that compared to healthy individuals, individuals with PF exhibit significant differences in foot kinematics and kinetics. Consistent with the theoretical injury mechanisms of PF, we found these individuals to have greater total rearfoot eversion and peak FMPJ dorsiflexion, which may put undue loads on the plantar fascia. Meanwhile, increased medial forefoot plantar flexion at initial contact and decreased propulsive GRF are suggestive of compensatory responses, perhaps to manage pain.     

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.     

Temporal characteristics of plantar shear distribution: relevance to diabetic patients

Diabetic foot ulcers are known to have a biomechanical etiology. Among the mechanical factors that cause foot lesions, shear stresses have been either neglected or underestimated. The purpose of this study was to determine various plantar pressure and shear variables in the diabetic and control groups and compare them. Fifteen diabetic patients with neuropathy and 20 non-diabetic subjects without foot symptoms were recruited. Subjects walked on a custom-built platform capable of measuring local normal and tangential forces simultaneously. Pressure-time integral quantities were increased by 54% (p=0.013) in the diabetic group. Peak AP and resultant shear magnitudes were found to be about 32% larger (p<0.05), even though diabetic subjects walked at a slower velocity. Lower AP and ML stress range (peak-to-peak) values were observed in the control subjects (p<0.05). Shear-time integral values were increased in the diabetic group by 61% and 132% for AP and resultant shear cases, respectively (p<0.05). Plantar shear is known to be a factor in callus formation and has previously been associated with higher ulcer incidence. During gait, shear stresses are induced with twice the frequency of pressure characteristically. Therefore, plantar shear should be investigated further from a broader perspective including the temporal specifications and fatigue failure characteristics of the affected plantar tissue.     

Dynamic plantar loading index: understanding the benefit of custom foot orthoses for painful pes cavus

The purpose of this study was to evaluate a new method showing how custom foot orthoses (CFO) improve dynamics of plantar loading. The method is based on the probability distribution of peak pressure time series and is quantified using the Regression Factor (RF). RF is a least square regression slope between the experimentally observed plantar pressure magnitude probability distribution and a modeled Gaussian shape. Plantar pressure data from a randomized controlled trial of 154 participants with painful Pes Cavus were retrospectively re-analyzed. The participants were randomized to an active treatment group given CFO or a control group given sham orthoses. The location of 2(nd) Peak pressure as a percentage of stance time (P(Loc2)) and its magnitude (P(M2)) was also calculated. In addition, plantar pressure data were collected on 23 healthy volunteers with normal foot alignment and no foot pain. Results demonstrated Pes Cavus had a significantly lower RF than healthy participants (0.30 v. 0.51; p<10(-7)). P(M2) was reduced in both active and control groups. However, RF and the P(Loc2) were only changed in the active group (p<0.005) without any significant change in the control group (p>0.5). This study suggests that painful Pes Cavus alters the shape of probability distribution of plantar loading during walking and CFO are an effective therapeutic solution that can significantly improve it. Further use of the RF index and 2(nd) peak pressure location as an outcome measure for treatment of foot and ankle deformities is suggested.     

An In Vivo Experimental Validation of a Computational Model of Human Foot

Reliable computational foot models offer an alternative means to enhance knowledge on the biomechanics of human foot.Model validation is one of the most critical aspects of the entire foot modeling and analysis process.This paper presents an invivo experiment combining motion capture system and plantar pressure measure platform to validate a three-dimensional finiteelement model of human foot.The Magnetic Resonance Imaging(MRI)slices for the foot modeling and the experimental datafor validation were both collected from the same volunteer subject.The validated components included the comparison of staticmodel predictions of plantar force,plantar pressure and foot surface deformation during six loading conditions,to equivalentmeasured data.During the whole experiment,foot surface deformation,plantar force and plantar pressure were recorded simultaneouslyduring six different loaded standing conditions.The predictions of the current FE model were in good agreementwith these experimental results.     

Functional Capacity Following Fractures of the Os Calcis

The long-term results, after applying the various methods of treatment for the different types of facture of the os calcis, were studied over a period of five years. Fractures with minimal displacement showed a satisfactory functional and economic end result. Those patients with the more severe types of fractures who had undergone arthrodesis early were free from pain early and were able to return to work earlier than those who underwent it as a delayed procedure.     

Analysis of the human and ape foot during bipedal standing with implications for the evolution of the foot

The ratio of the power arm (the distance from the heel to the talocrural joint) to the load arm (that from the talocrural joint to the distal head of the metatarsals), or RPL, differs markedly between the human and ape foot. The arches are relatively higher in the human foot in comparison with those in apes. This study evaluates the effect of these two differences on biomechanical effectiveness during bipedal standing, estimating the forces acting across the talocrural and tarsometatarsal joints, and attempts to identify which type of foot is optimal for bipedal standing. A simple model of the foot musculoskeletal system was built to represent the geometric and force relationships in the foot during bipedal standing, and measurements for a variety of human and ape feet applied. The results show that: (1) an RPL of around 40% (as is the case in the human foot) minimizes required muscle force at the talocrural joint; (2) the presence of an high arch in the human foot reduces forces in the plantar musculature and aponeurosis; and (3) the human foot has a lower total of force in joints and muscles than do the ape feet. These results indicate that the proportions of the human foot, and the height of the medial arch are indeed better optimized for bipedal standing than those of apes, further suggesting that their current state is to some extent the product of positive selection for enhanced bipedal standing during the evolution of the foot.