The effect of heel-pad thickness and loading protocol on measured heel-pad stiffness and a standardized protocol for inter-subject comparability

Current poisson's ratio values of finite element models are too low to consider soft tissues nearly-incompressible: illustration on the human heel region

Tissues' nearly incompressibility was well reported in the literature but little effort has been made to compare volume variations computed by simulations with in vivo measurements. In this study, volume changes of the fat pad during controlled indentations of the human heel region were estimated from segmented medical images using digital volume correlation. The experiment was reproduced using finite element modelling with several values of Poisson's ratio for the fat pad, from 0.4500 to 0.4999. A single value of Poisson's ratio could not fit all the indentation cases. Estimated volume changes were between 0.9% - 11.7%.

Assessing reliability and validity of different stiffness measurement tools on a multi-layered phantom tissue model

Changes in the mechanical properties (i.e., stiffness) of soft tissues have been linked to musculoskeletal disorders, pain conditions, and cancer biology, leading to a rising demand for diagnostic methods. Despite the general availability of different stiffness measurement tools, it is unclear as to which are best suited for different tissue types and the related measurement depths. The study aimed to compare different stiffness measurement tools' (SMT) reliability on a multi-layered phantom tissue model (MPTM). A polyurethane MPTM simulated the four layers of the thoracolumbar region: cutis (CUT), subcutaneous connective tissue (SCT), fascia profunda (FPR), and erector spinae (ERS), with varying stiffness parameters. Evaluated stiffness measurement tools included Shore Durometer, Semi-Electronic Tissue Compliance Meter (STCM), IndentoPRO, MyotonPRO, and ultrasound imaging. Measurements were made by two independent, blinded examiners. Shore Durometer, STCM, IndentoPRO, and MyotonPRO reliably detected stiffness changes in three of the four MPTM layers, but not in the thin (1 mm thick) layer simulating FPR. With ultrasound imaging, only stiffness changes in layers thicker than 3 mm could be measured reliably. Significant correlations ranging from 0.70 to 0.98 (all p < 0.01) were found. The interrater reliability ranged from good to excellent (ICC(2,2) = 0.75-0.98). The results are encouraging for researchers and clinical practitioners as the investigated stiffness measurement tools are easy-to-use and comparatively affordable.

Analysis and Assessment through Mechanical Static Compression Tests of Damping Capacity in a Series of Orthosis Plantar Materials Used as Supports

High plantar pressure is the cause of multiple types of foot injuries and one of the main reasons for patient visits in podiatry and traumatology. Therefore, there is a need to acquire specific tools to address such injuries. The aim of this study was to determine the absorption capacity of selected materials applied as plantar supports and their response to pressure. The study had a cross-sectional design. A total of 21 materials were chosen and different material families were assessed, including ethylene-vinyl acetate, polyurethane foams, and polyethylene foams. Static compression tests were performed to analyze each material. The system is ideally suited for lower-force applications, small components, biomedical applications, and lower-strength materials. Damping was determined using mathematical calculations performed on the study data. It was found that materials with a low Shore A, or soft materials, exhibited worse absorption capacity than harder materials. Ethyl-vinyl acetates had good absorption capacity, polyurethane foams had a poor absorption capacity, and soft materials provided better adaption to impact. The results suggested that damping is not determined by the hardness of the material, and materials within the same family exhibit different damping capabilities.

A novel simplified biomechanical assessment of the heel pad during foot plantarflexion

The heel pad (HP) which is located below the calcaneus comprises a composition of morphometrical and morphological arrangements of soft tissues that are influenced by factors such as gender, age and obesity. It is well known that HP pain and Achilles tendonitis consist of discomfort, pain and swelling symptoms that usually develop from excessive physical activities such as walking, jumping and running. The purpose of this study was to develop biomechanical techniques to evaluate the function and characteristics of the HP. Ten healthy participants (five males and five females) participated in this laboratory-based study, each performing a two-footed heel raise to mimic the toe-off phase during human locomotion. Twenty-six (3 mm) retroreflective markers were attached to the left and right heels (thirteen markers on each heel). Kinematic data was captured using three-dimensional motion analysis cameras synchronised with force plates. Descriptive and multivariate statistical tests were used in this study. In addition, a biomechanical technique that utilises only six markers from 26 markers to assess HP deformation and function has been developed and used in this study. Overall HP displacement was significantly higher in males on the most lateral part of the right heel (p < 0.05). No significant differences were evident when comparing the non-dominant and dominant heels during the baseline, unloading and loading phases (p > 0.05). Findings from this study suggested that biomechanical outputs expressed as derivatives from tracked HP marker movements can morphologically and morphometrically characterise HP soft tissue deformation changes. The outcome of this study highlights the importance of 3D motion analysis being used as a potential prospective intervention to quantify the function / characteristics of the heel pad soft tissues.

Adaptation to Coriolis force perturbations of postural sway requires an asymmetric two-leg model

In the companion paper (Bakshi A, DiZio P, Lackner JR. J Neurophysiol. In press, 2019), we reported how voluntary forward-backward sway in a rotating room generated medial-lateral Coriolis forces that initially deviated intended body sway paths. Pure fore-aft sway was gradually restored over per-rotation trials, and a negative aftereffect occurred during postrotation sway. Force plate recordings showed that subjects learned to compensate for the Coriolis forces by executing a bimodal torque, the distribution of which was asymmetric across the two legs and of opposite sign for forward vs. backward sway. To explain these results, we have developed an asymmetric, nonparallel-leg, inverted pendulum model to characterize upright balance control in two dimensions. Fore-aft and medial-lateral sway amplitudes can be biomechanically coupled or independent. Biomechanical coupling occurs when Coriolis forces orthogonal to the direction of movement perturb sway about the ankles. The model includes a mechanism for alternating engagement/disengagement of each leg and for asymmetric drive to the ankles to achieve adaptation to Coriolis force-induced two-dimensional sway. The model predicts the adaptive control underlying the adaptation of voluntary postural sway to Coriolis forces. A stability analysis of the model generates parameter values that match those measured experimentally, and the parameterized model simulations reproduce the experimentally observed sway trajectories. NEW & NOTEWORTHY This paper presents a novel nonparallel leg model of postural control that correctly predicts the perturbations of voluntary sway that occur in a rotating environment and the adaptive changes that occur to restore faithful movement trajectories. This engaged leg model (ELM) predicts the asymmetries in force distribution and their patterns between the two legs to restore accurate movement trajectories. ELM has clinical relevance for pathologies that generate postural asymmetries and for altered gravitoinertial force conditions.

Mechanical Modeling of Healthy and Diseased Calcaneal Fat Pad Surrogates

The calcaneal fat pad is a major load bearing component of the human foot due to daily gait activities such as standing, walking, and running. Heel and arch pain pathologies such as plantar fasciitis, which over one third of the world population suffers from, is a consequent effect of calcaneal fat pad damage. Also, fat pad stiffening and ulceration has been observed due to diabetes mellitus. To date, the biomechanics of fat pad damage is poorly understood due to the unavailability of live human models (because of ethical and biosafety issues) or biofidelic surrogates for testing. This also precludes the study of the effectiveness of preventive custom orthotics for foot pain pathologies caused due to fat pad damage. The current work addresses this key gap in the literature with the development of novel biofidelic surrogates， which simulate the in vivo and in vitro compressive mechanical properties of a healthy calcaneal fat pad. Also, surrogates were developed to simulate the in vivo mechanical behavior of the fat pad due to plantar fasciitis and diabetes. A four-part elastomeric material system was used to fabricate the surrogates, and their mechanical properties were characterized using dynamic and cyclic load testing. Different strain (or displacement) rates were tested to understand surrogate behavior due to high impact loads. These surrogates can be integrated with a prosthetic foot model and mechanically tested to characterize the shock absorption in different simulated gait activities, and due to varying fat pad material property in foot pain pathologies (i.e., plantar fasciitis, diabetes, and injury). Additionally, such a foot surrogate model, fitted with a custom orthotic and footwear, can be used for the experimental testing of shock absorption characteristics of preventive orthoses.

In-vivo viscous properties of the heel pad by stress-relaxation experiment based on a spherical indentation

Identifying the viscous properties of the plantar soft tissue is crucial not only for understanding the dynamic interaction of the foot with the ground during locomotion, but also for development of improved footwear products and therapeutic footwear interventions. In the present study, the viscous and hyperelastic material properties of the plantar soft tissue were experimentally identified using a spherical indentation test and an analytical contact model of the spherical indentation test. Force-relaxation curves of the heel pads were obtained from the indentation experiment. The curves were fit to the contact model incorporating a five-element Maxwell model to identify the viscous material parameters. The finite element method with the experimentally identified viscoelastic parameters could successfully reproduce the measured force-relaxation curves, indicating the material parameters were correctly estimated using the proposed method. Although there are some methodological limitations, the proposed framework to identify the viscous material properties may facilitate the development of subject-specific finite element modeling of the foot and other biological materials.

Material properties of the heel fat pad across strain rates

The complex structural and material behaviour of the human heel fat pad determines the transmission of plantar loading to the lower limb across a wide range of loading scenarios; from locomotion to injurious incidents. The aim of this study was to quantify the hyper-viscoelastic material properties of the human heel fat pad across strains and strain rates. An inverse finite element (FE) optimisation algorithm was developed and used, in conjunction with quasi-static and dynamic tests performed to five cadaveric heel specimens, to derive specimen-specific and mean hyper-viscoelastic material models able to predict accurately the response of the tissue at compressive loading of strain rates up to 150 s⁻¹. The mean behaviour was expressed by the quasi-linear viscoelastic (QLV) material formulation, combining the Yeoh material model (C10=0.1MPa, C30=7MPa, K=2GPa) and Prony׳s terms (A1=0.06, A2=0.77, A3=0.02 for τ1=1ms, τ2=10ms, τ3=10s). These new data help to understand better the functional anatomy and pathophysiology of the foot and ankle, develop biomimetic materials for tissue reconstruction, design of shoe, insole, and foot and ankle orthoses, and improve the predictive ability of computational models of the foot and ankle used to simulate daily activities or predict injuries at high rate injurious incidents such as road traffic accidents and underbody blast.

Dynamic material characterization of the human heel pad based on in vivo experimental tests and numerical analysis

A numerical-experimental, proof-of-concept approach is described to characterize the mechanical material behavior of the human heel pad under impact conditions similar to a heel strike while running. A 3D finite-element model of the right foot of a healthy female subject was generated using magnetic resonance imaging. Based on quasi-static experimental testing of the subject's heel pad, force-displacement data was obtained. Using this experimental data as well as a numerical optimization algorithm, an inverse finite-element analysis and the 3D model, heel pad hyperelastic (long-term) material parameters were determined. Applying the same methodology, based on the dynamic experimental data from the impact test and obtained long-term parameters, linear viscoelastic parameters were established with a Prony series. Model validation was performed employing quasi-static and dynamic force-displacement data. Coefficients of determination when comparing model to experimental data during quasi-static and dynamic (initial velocity: 1480mm/s) procedure were R(2) = 0.999 and R(2) = 0.990, respectively. Knowledge of these heel pad material parameters enables realistic numerical analysis to evaluate internal stress and strain in the heel pad during different quasi-static or dynamic load conditions.

The relationship between the mechanical properties of heel-pad and common clinical measures associated with foot ulcers in patients with diabetes

AIM

The present study aims at investigating the correlation between the mechanical properties of the heel-pad of people with type-2 diabetes and the clinical parameters used to monitor their health and ulceration risk.

METHODS

A new device for the in-vivo testing of plantar soft tissues was built and pilot-tested. This device consists of an ultrasound probe connected in series with a dynamometer. Loading is applied manually using a ball-screw actuator. A total of 35 volunteers with type-2 diabetes were recruited and the thickness, stiffness of their heel-pads as well as the energy absorbed during loading were assessed. The participants with diabetes also underwent blood tests and measurements of Ankle Brachial Index and Vibration Perception Threshold.

RESULTS

Pearson correlation analysis revealed strong correlations between triglycerides and heel-pad stiffness (r=0.675, N=27, p<0.001) and between triglycerides and energy (r=-0.598, N=27, p=0.002). A correlation of medium strength was found between Fasting Blood Sugar (FBS) and stiffness (r=0.408, N=29, p=0.043).

CONCLUSIONS

People with type-2 diabetes and high levels of triglycerides and FBS are more likely to have stiffer heel-pads. Increased stiffness could limit the tissues' ability to evenly distribute loads making them more vulnerable to trauma and ulceration.

A stereologic study of the plantar fat pad in young and aged rats

Plantar fat pad (PFP) is a tissue structure that absorbs the initial impact of walking and running and ultimately bears body weight at standing. This study was designed to quantify the histomorphological changes of the PFP in aged rats. The most medial PFP was dissected from the hind feet of young rats (4 months old, n = 6) and aged rats (24 months old, n = 6). Histological structure and cellular senescence of PFP were analyzed stereologically and histomorphometrically. Immunohistochemistry of matrix metalloproteinase 9 (MMP9) was also performed on PFP tissue sections. Compared with young rats, the thickness of epidermis, dermis and septa of the PFP were significantly reduced in the aged rats. The total volume of adipose tissue in the PFP of aged rats was only about 65% of that in the young rats. The microvascular density and the number of fat pad units (FPU), a cluster of adipocytes enclosed by elastin septa, in the PFP were unchanged in the aged rats. In the aged rats, the number of adipocytes per FPU was reduced but the number of simple adipocyte clusters, without surrounding septa, was increased. The shift of the types of adipocyte clusters in the aged PFP was accompanied by degradation of elastin fibers and increased expression of MMP9. In conclusion, the PFP, particularly the elastic septa, degenerates significantly in aged rats and this may contribute to the pathology of PFP-related diseases.

Internal strain estimation for quantification of human heel pad elastic modulus: A phantom study

Shock absorption is the most important function of the human heel pad. However, changes in heel pad elasticity, as seen in e.g. long-distance runners, diabetes patients, and victims of Falanga torture are affecting this function, often in a painful manner. Assessment of heel pad elasticity is usually based on one or a few strain measurements obtained by an external load-deformation system. The aim of this study was to develop a technique for quantitative measurements of heel pad elastic modulus based on several internal strain measures from within the heel pad by use of ultrasound images. Nine heel phantoms were manufactured featuring a combination of three heel pad stiffnesses and three heel pad thicknesses to model the normal human variation. Each phantom was tested in an indentation system comprising a 7MHz linear array ultrasound transducer, working as the indentor, and a connected load cell. Load-compression data and ultrasound B-mode images were simultaneously acquired in 19 compression steps of 0.1mm each. The internal tissue displacement was for each step calculated by a phase-based cross-correlation technique and internal strain maps were derived from these displacement maps. Elastic moduli were found from the resulting stress-strain curves. The elastic moduli made it possible to distinguish eight of nine phantoms from each other according to the manufactured stiffness and showed very little dependence of the thickness. Mean elastic moduli for the three soft, the three medium, and the three hard phantoms were 89kPa, 153kPa, and 168kPa, respectively. The combination of ultrasound images and force measurements provided an effective way of assessing the elastic properties of the heel pad due to the internal strain estimation.

Investigations on the viscoelastic behaviour of a human healthy heel pad: In vivo compression tests and numerical analysis

The aim of this study was to investigate the viscoelastic behaviour of the human heel pad by comparing the stress–relaxation curves obtained from a compression device used on an in vivo heel pad with those obtained from a three-dimensional computer-based subject-specific heel pad model subjected to external compression. The three-dimensional model was based on the anatomy revealed by magnetic resonance imaging of a 31-year-old healthy female. The calcaneal fat pad tissue was described with a viscohyperelastic model, while a fibre-reinforced hyperelastic model was formulated for the skin. All numerical analyses were performed to interpret the mechanical response of heel tissues, with loading conditions and displacement rate in agreement with experimental tests. The heel tissues showed a non-linear, viscoelastic behaviour described by characteristic hysteretic curves, stress–relaxation and viscous recovery phenomena. The reliability of the investigations was validated by the interpretation of the mechanical response of heel tissues under the application of three pistons with diameter of 15, 20 and 40 mm, at the same displacement rate of about 1.7 mm/s. The maximum and minimum relative errors were found to be less than 0.95 and 0.064, respectively.

HYPER-ELASTIC PARAMETER ESTIMATION OF HUMAN HEEL-PAD: A FINITE ELEMENT AND EVOLUTIONARY BASED ALGORITHM

The heel-pad as a biological shock absorber has an important role in the initial contact phase of gait cycle dissipating the impact forces resulted in locomotion. An axisymmetric finite element model of human heel-pad has been generated and the heel-pad experimental data deduced from a published force-deflection graph of the same specimen (Iain R. Spears, Janice E. Miller-Young), Iterative identification task has been used to extract nonlinear material properties describing hyper-elastic behavior of heel-pad. The genetic algorithm was incorporated into estimation process using an interface program. Two parameters of hyper-elastic materials potential energy function represented by Mooney–Rivlin were determined by using the genetic algorithm technique to minimize the displacement error between the experimental data and the corresponding finite element results after a considerable number of iterations. The result can be used for design and construction of synthetic heel-pad and therapeutic foot wear as well as insoles, especially for diabetic patients.

Biomechanical behaviour of heel pad tissue experimental testing, constitutive formulation, and numerical modelling

This paper deals with the constitutive formulation of heel pad tissue and presents a procedure for identifying constitutive parameters using experimental data, with the aim of developing a computational approach for investigating the actual biomechanical response. The preliminary definition of constitutive parameters was developed using a visco-hyperelastic formulation, considering experimental data from in vitro compression tests on specimens of fat pad tissue and data from in vivo tests to identify the actual trend of tissue stiffness. The discrepancy between model results and experimental data was evaluated on the basis of a specific cost function, adopting a stochastic/deterministic procedure. The parameter evaluation was upgraded by considering experimental tests performed on the fat pad tissues of a cadaveric foot using in situ indentation tests at 0.01 and 350 mm/s strain rates. The constitutive formulation was implemented in a numerical model. The comparison of data from in situ tests and numerical results led to an optimal domain of parameters based on an admissible discrepancy criterion. Numerical results evaluated for different sets of parameters inside the domain are reported and compared with experimental data for a reliability evaluation of the proposed procedure.

Diabetic effects on microchambers and macrochambers tissue properties in human heel pads

BACKGROUND

The study attempted to highlight the differences of mechanical properties in microchambers and macrochambers between patients with type 2 diabetes mellitus and age-matched healthy volunteers.

METHODS

A total of 29 heels in 18 diabetic patients and 28 heels in 16 age-matched healthy participants were examined by a loading device consisting of a 10-MHz compact linear-array ultrasound transducer, a Plexiglas cylinder, and a load cell. Subjects in both groups were on average about 55 years old with a body mass index of approximately 25 kg/m(2). A stepping motor was used to progressively load the transducer on the tested heels at a velocity of 6mm/s from zero to the maximum stress of 78 kPa. Unloaded thickness, strain, and elastic modulus in microchambers, macrochambers and heel pads were measured.

FINDINGS

Microchambers strain in diabetic patients was significantly greater than that in healthy subjects (0.291 (SD 0.14) vs. 0.104 (SD 0.057); P<0.001). Macrochambers strain in diabetic patients was significantly less than that in healthy subjects (0.355 (SD 0.098) vs. 0.450 (SD 0.092); P=0.001). Microchambers stiffness in diabetic patients was significantly less than that in healthy persons (393 (SD 371)kPa vs. 1140 (SD 931)kPa; P<0.001). Macrochambers stiffness in diabetic patients was significantly greater than that in healthy persons (239 (SD 77)kPa vs. 181 (SD 42)kPa; P=0.001).

INTERPRETATION

Heel pad tissue properties are altered heterogeneously in people with diabetes. Increased macrochambers but decreased microchambers stiffness may cause diminished cushioning capacities in diabetic heels.

The Diabetic Foot Load Monitor: A Portable Device for Real-Time Subject-Specific Measurements of Deep Plantar Tissue Stresses During Gait

Elevated stresses in deep plantar tissue of diabetic neuropathic patients were associated with an increased risk for foot ulceration, but only interfacial foot pressures are currently measured to evaluate susceptibility to ulcers. The goals of this study were to develop a real-time patient-specific plantar tissue stress monitor based on the Hertz contact theory. The biomechanical model for stress calculations considers the heel and metatarsal head pads, where most ulcers occur. For calculating stress concentrations around the bone-pad interface, plantar tissue is idealized as elastic and incompressible semi-infinite bulk (with properties measured by indentation), which is penetrated by a rigid sphere with the bone’s radius of curvature (from X-ray). Hertz’s theory is used to solve the bone-pad mechanical interactions, after introducing correction coefficients to consider large deformations. Foot-shoe forces are measured to solve and display the principal compressive, tensile, and von Mises plantar tissue stresses in real time. Our system can be miniaturized in a handheld computer, allowing plantar stress monitoring in the patient’s natural environment. Small groups of healthy subjects (N=6) and diabetic patients (N=3) participated in an evaluation study in which the differences between free walking and treadmill walking were examined. We also compared gait on a flat surface to gait on an ascending/descending slope of 3.5deg and when ascending/descending stairs. Peak internal compression stress was about threefold greater than the interface pressure at the calcaneus region. Subjects who were inexperience in treadmill walking displayed high gait-cycle variability in the internal stresses as well as poor foot loading. There was no statistical difference between gait on a flat surface and gait when ascending/descending a slope. Internal stresses under the calcaneus during gait on a flat surface, however, were significantly higher than when ascending/descending stairs. We conclude that the present stress monitor is a promising tool for real-time patient-specific evaluation of deep tissue stresses, providing valuable information in the effort to protect diabetic patients from foot ulceration. Clinical studies are now underway to identify which stress parameters can distinguish between diabetic and normal subjects; these parameters may be used for establishing injury threshold criteria.

Effects of varying material properties on the load deformation characteristics of heel cushions

Various insole materials were used in attenuation of heel-strike impact. This study presented a compression test to investigate the deformation characteristics of common heel cushions. There were two materials (thermoplastic elastomer "TPE" and silicone) with three hardness and six thickness being analyzed. They underwent consecutive loading-unloading cycles with a load control mode. The displacement of material thickness was recorded during cyclic compression being applied and released from 0 to 1050 N. The energy input, return and dissipation were evaluated based on the load deformation curves when new and after repeated compression. The TPE recovered more deformed energy and thickness than the silicone after the first loading cycle. The silicone would preserve more strain energy with increasing its hardness for the elastic recovery in the unloading process. The deformed energy was decreased as the original thickness did not completely recover under cyclic tests. The reduction in hysteresis area was gradually converged within 20 cycles. The silicone attenuated more impact energy in the initial cycles, but its energy dissipation was reduced after repeated loading. To increase hardness or thickness should be considered to improve resilience or accommodate persistent compression without flattening. The careful selection of cushion materials is imperative to meet individual functional demands.

Microchambers and macrochambers in heel pads: are they functionally different?

The heel pad consists of a superficial microchamber layer and a deep macrochamber layer. This study highlights the different biomechanical behaviors between the microchamber and macrochamber layers using ultrasonography. The heel pad in each left foot of six healthy volunteers aged ∼25 yr old was measured with a device consisting of a 10-MHz linear-array ultrasound transducer and a load cell. The testing heels were loaded on the ultrasound transducer with a loading velocity of ∼0.5 cm/s and were withdrawn when the specified maximum stress (158 kPa) was reached. Unloaded tissue thickness, end-loaded thickness, deformation proportion, average deformation, and rebound rates and elastic modulus of the microchamber and macrochamber layers were assessed. The unloaded thickness of the microchamber layer was ∼30% of the macrochamber layer. The microchamber layer also had significantly less unloaded thickness, end-loaded thickness, mean deformation rate, mean rebound rate, and deformation proportion than the macrochamber layer. A significant difference between the unloaded and end-loaded thickness in the macrochamber layer was observed. The average soft tissue deformation rate was significantly different from the rebound rate in the microchamber layer. A similar trend was detected in the macrochamber layer. The elastic modulus of the microchamber layer was 450 kPa (SD 240), which was nearly 10 times of that in the macrochamber layer. In conclusion, ultrasound can identify the heterogeneous tissue properties of the heel pad. The macrochamber layer responds to loading with large deformation, and the microchamber layer has a high degree of tissue stiffness.

The potential influence of the heel counter on internal stress during static standing: A combined finite element and positional MRI investigation

Confinement of the heel due to the counter of the shoe is believed to influence heel pad biomechanics. Using a two-dimensional finite element model of the heel pad and shoe during a simulation of static standing, the aim of this study was to quantify the potential effect of confinement on internal heel pad stress. Non-weightbearing MRI and weightbearing MRI with plantar pressure and ground reaction force data were recorded for a single subject. The non-weightbearing MRI was used to create two FE models of the heel pad, using either homogeneous or composite material properties. The composite model included a distinction in material properties between fat pad and skin. Vertical and medial-lateral forces, as measured on the subject's heel, were applied to the models and vertical compressive strains for both models were comparable with those observed by weightbearing MRI. However, only for the composite model was the predicted plantar pressure distribution comparable with measured data. The composite model was therefore used in further analyses. In this composite model, the internal stresses were located mainly in the skin and were predominantly tensile in nature, whereas the stress state in the fat pad approached hydrostatic conditions. A representation of a running shoe, including an insole, midsole and heel counter was then added to the composite heel pad to form the shod model. In order to investigate the counter effect, the load was applied to the shod model with and without the heel counter. The effect of the counter on peak stress was to elevate compression (0-50%), reduce tension (22-34%) and reduce shear (22-28%) in the skin. In addition, the counter reduced both compressive (20-40%) and shear (58-80%) stress in the fat pad and tension in the fat pad remained negligible. Taken together the results indicate that a well-fitted counter works in sympathy with the internal structure of the heel pad and could be an effective reducer of heel pad stress. However, further research needs to be undertaken to assess the long-term effects on the soft-tissues, practicalities of achieving good fit and behavior under dynamic events.