Medical, Pharma, Engineering, Science, Technology and Business

^{1}BGKT Consulting Ltd., BioEngineering, Huntington, New York, USA

^{2}Draper Laboratory, Research, Cambridge, Massachusetts, USA

- *Corresponding Author:
- Greene PR

BGKT Consulting Ltd

BioEngineering, Huntington

New York 11743, USA

**Tel:**+1 6319355666

**E-mail:**[email protected]

**Received date:** May 27, 2015; **Accepted date:** June 15, 2015; **Published date:** June 17, 2015

**Citation:** Greene PR, Coleman JD (2015) Elastic Spring Constants for Running Shoes: A Mathematical Model. J Comput Sci Syst Biol 08:215-218. doi: 10.4172/jcsb.1000191

**Copyright:** © 2015 Greene PR, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Computer Science & Systems Biology

**Background**: Running shoe compliance and track surface stiffness can reduce peak vertical foot forces. It is therefore of interest to measure directly the force-deflection curve for running shoes in the heel and forefoot areas. This study compares these measurements with similar work on track and field surfaces, and derives a mathematical stress-strain model useful over the entire force range.

**Methods**: Six different running shoes from 4 popular brands are measured to determine vertical spring stiffness. The heel and ball areas are tested with 3.8 and 5.1 cm (1.5 and 2.0 in.) diameter heels in the force range of 0 to 0.13 kN and 2.0 to 2.6 kN (0 to 30 lbf. and 450 to 600 lbf.). The results show a factor of 2 difference from one shoe to the next, holding test area, heel diameter and force range constant. Load increments are applied on a time scale of 0.1 seconds, comparable to typical foot contact times during running.

**Results**: The measured spring constants are essentially independent of plunger area, a useful simplification. For a given shoe, the ball area can be three times less compliant than the heel.

**Conclusion**: Heel spring constants at the high force levels fall in the range from 290 kN/m to 600 kN/m (20,000 to 42,000 lbf/ft.), and thus approximate optimal track stiffness. In terms of theory, an exponential function derives from the observation that the data fall along a straight line on semi-log co-ordinates. This mathematical model enables calculation of running shoe compressive response and spring constant at physiological force levels.

Running shoe compliance; Track stiffness; Compressive stress-strain; Orthopaedics; Force-deflection curve; Area independence

Compliance in the track or shoe is an important design parameter that can reduce foot force significantly, as reported in many technical reports [1-11]. The effects of compliance on vertical foot force, i.e. the vertical component of the ground reaction force, are reported by Aerts and de Clercq [1] . Running shoe and track compliance is reviewed by Frederick [12], as it influences locomotion at various speeds.

The rationale and motivation for this study is force reduction. The forces on the foot during running and jogging may contribute to soft tissue injuries of the calcaneous [3], Achilles tendon [4,13] anterior compartment (shin splints) [14], and knee ligaments and cartilage [11,15]. Excessive vertical force is of particular interest to amputees [16] in terms of prosthetic stability and comfort [17]. Shoe compliance measurements can vary according to the type of compression machine, indentor, and force range [2,12,18,19]. Using optimal compliance values for tracks and running shoes, peak vertical forces can be reduced without sacrificing speed, so one objective of this report is to compare vertical spring constants of typical running shoes with those of optimal running tracks [7,8,20]. Laboratory measurements of shoe dynamics in the heel and forefoot areas of the sole enable comparison with track stiffness measurements. This report investigates whether compliance results for a flat indentor are independent of plunger area, an important experimental and theoretical simplification. Lastly, a force-deflection equation is derived, facilitating calculation of shoe compliance over the entire physiological force range, applicable to the dynamics of walking, jogging, and running. In general, running shoes are the softest during walking, least compliant during sprinting periods. Serpell et al. [21] review physiological stiffness calculations of the leg and knee. Bo [22] discusses the difficult task of computing deflection of an elastic half space.

Nomenclature

A_{o}=plunger area [cm^{2} (sq. in.)]

ε=strain [dh/h_{0} (%)]

F=force=A_{o} × σ (ε) [lbs or N],

σ=stress=A [exp(α × ε)-1], [N/cm^{2} (lbf./in^{2})]

A=exponential material constant [N/cm^{2} (lbf./in^{2})]

α=exponential stiffening constant

ho=sole thickness [cm (inches)]

dh=vertical deflection [cm (mils.)]

A representative sample of 6 training shoes from 4 different manufacturers was measured to determine their spring constants. They were subjected to forces between 0 to 0.13 kN and 2.0 to 2.6 kN (0-30 and 450-600 lbf.) in the heel and ball (forefoot) regions of the shoe. The important tests were in the 2.0 to 2.6 kN range (450-600 lbf.), since this is the force level during most of the foot contact time during running. As measured, the force-deflection curve is non-linear; this is mathematically modeled here using basic stress-strain constitutive equations. One purpose of this report is to explore the feasibility of these experimental procedures together with basic equations, so several representative running shoes (N=6) are selected as typical samples, in order to demonstrate the utility of these techniques.

Sole thickness varies from the heel to ball region of the running shoe, so it was necessary to make separate measurements in each area. Two different size artificial heels were used, circular disks, with diameters D=3.8 cm and D=5.1 cm (1.5” and 2.0”), **Figure 1**, applied to the heel and forefoot of the shoe, to measure if an area independent index could be determined. The shoe uppers and soles were not cut in any way in order to perform these measurements, remaining intact along with the original liner inserts. The bottom of the circular aluminum heel, flat in profile, was teflon coated with a 1.6 mm (1/16”) thick sheet to reduce friction during the compression tests, and to not mar the shoe. As predicted by theory, an area independent index is experimentally possible, as shown in** Figure 2**, because the measured effective spring constants are independent of the heel diameter. An area independent index is basically a mechanical property independent of the indentor area, an important generalization.

The dynamic velocity component of the shoe sole, i.e. the dashpot component in the constituitive relation, was not measured directly. Typical foot contact times are the order of 0.10 to 0.15 sec during running, 0.3 to 0.8 sec during walking, jogging, accelerating, and turning [20,23-25]. The incremental loads in these experiments are applied on a time scale of 0.1 secs, so the elastic spring constants reported here are physiologically realistic, particular in the high force range.

The dynamic velocity component of the shoe sole, i.e. the dashpot component in the constituitive relation, was not measured directly. Typical foot contact times are the order of 0.10 to 0.15 sec during running, 0.3 to 0.8 sec during walking, jogging, accelerating, and turning [20,23-25]. The incremental loads in these experiments are applied on a time scale of 0.1 secs, so the elastic spring constants reported here are physiologically realistic, particular in the high force range. shoe on a time scale of 0.1 secs. An adjustable counter-weight is used, shown at the left in Figure 1, to null out and level the weight of the moment arm before the steel weights are applied. All shoes used in these experiments were men’s size 8-1/2. We recorded shoe stiffness constant k=dF/dx at the two different load ranges to characterize the non-linear material properties. An exponential constituitive relation, Fung [23] enables mathematical prediction of the stress-strain and load-deflection curves using the equation:

Eq. (1) F=A_{o} × σ (ε), where

A_{o}=plunger area [cm^{2} (sq. in.)]

σ=stress [N/cm^{2} (lbf./in^{2})]

ε=strain [dh/h_{0} (%)]

Figure 2 shows the results of all the force-deflection tests. For any given combination of shoe, test locale, plunger area, and force range, each test was repeated 4 times. The error bars indicate the standard deviation of these 4 repeated measurements. Overall there is a decrease in spring constant for increasing sole thickness, i.e. the shoe is more compliant as expected. Of the 6 running shoes tested, No. 5, Figure 2, has the smallest average spring constant in both the heel and ball area, and thus is the most compliant.

The heel value of k=332 kN/m (23,000 lbf./ft.) is close to the optimal value for running tracks [7,8]. The same weights, deflection gauge, and testing procedure used for the track and field experiments were also used for the shoe compressions. As shown in Figure 2, the ball area of the shoe is several times stiffer than the heel for all of the shoes tested. During running, a significant portion of the foot contact phase takes place with either both the heel and ball or just the ball of the foot in contact with the ground. Thus, the measured running shoes are several times stiffer than the optimal vertical compliance for track and field surfaces.

Following the theoretical results of Fung [26], the data are fit to an exponential stress-strain curve:

Eq. (2) σ (ε)=A [exp (α × ε)-1], where

σ=stress [N/cm^{2} (lbf./in^{2})]

ε=strain [dh / h_{0}]

A=material constant [N/cm^{2} (lbf./in^{2})]

α=exponential stiffening constant

Experimental results are shown in **Figures 3 and 4**. Figure 3 is a force-deflection graph for the ball region (squares) and the heel region (circles) of a running shoe. Typical load levels for walking, jogging, and running are shown, indicating that the effective stiffness does vary with the applied force level. Figure 4 is a semi-log stress-strain graph for the ball region of one of the six running shoes indicating that the sole material has an exponential stress-strain constituitive law as derived by Fung [26]. The +/-rms error bars are larger in the forefoot area (Figure 2) because the deflection in this area is approximately 1/3 that of the heel area.

In terms of practical clinical applications, excessive ground reaction forces on the foot during running and jogging may contribute to soft tissue injuries of the calcaneous, Achilles tendon, anterior compartment (shin splints), and knee ligaments and cartilage [11,14,16]. Excessive vertical force is of particular interest to amputees [16], in terms of shock transmission characteristics after foot strike [17], from the ground plain sequentially on upwards through the body. Thus, optimizing the running shoe spring constant to minimize foot force is a welldefined design objective. Other important design and experimental parameters include indentor area, hysteresis, and exponential stressstrain, as discussed below. Mathematical theory is important because it enables continuous calculation of the force-deflection and stress-strain curves over the entire force range, applicable to the forces and stresses typical of walking, jogging, and running, a useful generalization. These equations (Eq.1 and Eq. 2, Fung [26]) are an important economy, in terms of computing time and effort, because the exact elastic half-space calculations can be quite involved [3,22].

This study shows that compliance results for a flat indentor are independent of plunger area, an important experimental and theoretical simplification. An exponential force-deflection equation allows calculation of shoe compliance over the entire physiological force range, applicable to the dynamics of human locomotion at varying speeds. It remains to be seen with future research if we can confirm that compliant running shoes exhibit a 2% to 3% increase in top speed, as reported by McMahon and Greene [7,8] for tuned tracks.

**Area independence**

The data show that there is no essential difference in measured vertical stiffness constant between the small and large plungers, although the disk area changes by a factor of 77%. This effect can be predicted on theoretical grounds [26] as a result of the exponentially stiffening stress-strain curve of the polymer materials, a result which is also true for tendons stressed in tension[13]. Thus, a possible area independent index may be the exponential stiffening constant [26], typically 3<α<10 for running shoes, Eq. 2.

**Hysteresis**

Viscoelastic effects, i.e. hysteresis, with different trajectories for loading and unloading, are not measured for all the shoes reported here. Preliminary measurements, not presented here, indicated there is a 10% to 15% difference in effective stiffness, resulting from the viscous dashpot effects.

**Exponential stress-strain relation**

Following the theoretical results of Fung [26], Eq. (2), the data are fit to an exponential stress-strain curve σ (ε)=A [exp (α × ε)-1], as shown in Figures 3 and 4, with data from one of the shoes. All six shoes demonstrate similar exponential stress-strain curves. Exponential material constants range from 25<A<45 and 4<α<6 for the heel region. Figure 3 shows that the effective spring constant of the shoe varies considerably according to whether one is walking, jogging, or running.

**Shoe design for other sports**

When walking, jogging, accelerating [25], or running on turns [23,24], foot contact time is lengthened considerably from 0.15 secs. to 0.3-0.8 secs., although the forces remain at a high level. Thus, the dynamics of these longer impulse foot-strikes are quite different from the dynamics of sprinting, suggesting different design parameters for basketball, soccer, and football shoes [27].

**Study limitations**

Extrapolating these results to shoes of different sizes, i.e. scaling laws, is a challenging question, beyond the scope of this report, requiring further research. Sole stiffness is found to be proportional to the reciprocal of thickness h^{-1} by substituting Eq. 1 into Eq. 2 and differentiating: K=dF/dx ~ h^{-1}. This is a useful simplification. Since only 6 representative running shoes from various companies were used, this study is by no means exhaustive. For commercial reasons, we do not want to exclude any of the excellent manufacturers, nor do we wish to emphasize one brand over another, so the company names and models are deleted.

The vertical spring stiffness constants of six popular running shoes are reported. The measured average spring stiffness in the heel area for the 2.0 to 2.6 kN range (450-600 lbf.), loaded on a time scale of 0.1 seconds, is close to that for an optimal running track for some of the shoes reported here. Spring constant is found to be independent of plunger area, a useful theoretical and experimental simplification, which can be derived from the basic exponential stress-strain law. In the forefoot area of the shoe, the spring constant is larger than in the heel by a factor of 2 to 3, as shown in Figure 2.

From the theoretical point of view, it is shown that the stress-strain data falls along a straight line on semi-log co-ordinates, Figure 4, indicating an exponential constitutive relation [26]. Vertical stiffness constants are measured for the six shoes in both the heel and ball areas, both for the small and large plunger disks. The independence of plunger area is true in both regions. Thus, using basic equations, we can calculate the compression characteristics and vertical spring constant of a running shoe at any force level. This then is directly applicable to typical physiological forces and stresses experienced during locomotion [20]. Stiffness is found to be proportional to the reciprocal of the sole thickness h^{-1}, a useful simplification.

This report demonstrates that compliance results for a flat indentor are independent of plunger area, an important experimental and theoretical simplification. An exponential force-deflection equation enables calculation of shoe compliance over the entire physiological force range, applicable to the dynamics of walking, jogging, and running.

The authors would like to thank Tom McMahon, Peter Cavanagh, and Ned Frederick, for many helpful discussions. Special thanks to the “Norton Disk Doctor”, without which we might have lost this report entirely.

- Aerts P, De Clercq D (1993) Deformation characteristics of the heel region of the shod foot during a simulated heel strike: the effect of varying midsole hardness.J Sports Sci 11: 449-461.
- Cook SD, Kester MA, Brunet ME (1985) Shock absorption characteristics of running shoes.Am J Sports Med 13: 248-253.
- Erdemir A,Viveiros ML, Ul brecht JS, Cavanagh PR (2006) An inverse finite-element model of heel-pad indentation.J Biomech 39: 1279-1286.
- Herzog W (2012) Running injuries: is it a question of evolution, form, tissue properties, mileage, or shoes?Exerc Sport Sci Rev 40: 59-60.
- Ly QH,Alaoui A, Erlicher S, Baly L (2010) Towards a footwear design tool: influence of shoe midsole properties and ground stiffness on the impact force during running.J Biomech 43: 310-317.
- Logan S, Hunter I, J Ty Hopkins JT, Feland JB, Parcell AC (2010) Ground reaction force differences between running shoes, racing flats, and distance spikes in runners. J Sports Sci Med 9:147-153.
- McMahon TA, Greene PR (1979) The influence of track compliance on running.J Biomech 12: 893-904.
- McMahon TA, Greene PR (1978) Fast running tracks.Sci Am 239: 148-163.
- Morag E, Cavanagh PR (1999) Structural and functional predictors of regional peak pressures under the foot during walking.J Biomech 32: 359-370.
- Tessutti V,Ribeiro AP, Trombini-Souza F, Sacco IC (2012) Attenuation of foot pressure during running on four different surfaces: asphalt, concrete, rubber, and natural grass.J Sports Sci 30: 1545-1550.
- Zhang Y, Liu G, Xie SQ (2011) Biomechanical simulation of anterior cruciate ligament strain for sports injury prevention.ComputBiol Med 41: 159-163.
- Frederick EC (1984) Sports Shoes and Playing Surfaces, Human Kinetics Publ., Champaign, Ill., ISBN # 978-0931-2505-14, 198.
- McMahon TA (1984) Muscles Reflexes and Locomotion, Princeton Univ. Press, Princeton N.J, pp.359.
- Greene PR (1985) Constant strain increment for exponential tendons in the high-stress limit.J BiomechEng 107: 291.
- Greene PR, McMahon TA (1979) Reflex stiffness of man's anti-gravity muscles during kneebends while carrying extra weights.J Biomech 12: 881-891.
- Küçük H (2006) The effect of modeling cartilage on predicted ligament and contact forces at the knee.ComputBiol Med 36: 363-375.
- Beurskens R,Wilken JM, Dingwell JB (2014) Dynamic stability of superior vs. inferior body segments in individuals with transtibial amputation walking in destabilizing environments.J Biomech 47: 3072-3079.
- vanLeeuwen JL,Speth LA, Daanen HA (1990) Shock absorption of below-knee prostheses: a comparison between the SACH and the Multiflex foot.J Biomech 23: 441-446.
- Serpell BG, Ball NB, Scarvell JM, Smith PN (2012) A review of models of vertical, leg, and knee stiffness in adults for running, jumping or hopping tasks.J Sports Sci 30: 1347-1363.
- Bo J (1998) Using fredholm integral equation of the second kind to solve the vertical vibration of elastic plate on an elastic half space. Applied Mathematics and Mechanics 19: 157-162.
- Greene PR, McMahon TA (1979) Running in circles.Physiologist 22: S35-36.
- Greene PR,Monheit MA (1990) Optimal geometry for oval sprint tracks.J Biomech 23: 447-452.
- Fung YC (1967) Elasticity of soft tissues in simple elongation.Am J Physiol 213: 1532-1544.
- Price C, Cooper G, Jones R (2015)Themanipulation of midsole properties to alter impact characteristics in walking. Footwear Science 7: 9-16.
- Price C, Cooper G, Graham-Smith P, Jones R (2014) A mechanical protocol to replicate impact in walking footwear.Gait Posture 40: 26-31.
- Greene PR (1986) Predicting sprint dynamics from maximum-velocity measurements. Mathematical Biosciences80: 1-18.
- Zanetti EM, BignardiC, Franceschini G, Audenino AL (2013) Amateur football pitches: mechanical properties of the natural ground and of different artificial turf infills and their biomechanical implications. J Sports Sci 31:767-778.

Select your language of interest to view the total content in your interested language

- Advanced DNA Sequencing
- Algorithm
- Animal and Tissue Engineering
- Applications of Bioinformatics
- Artificial Intelligence Studies
- Artificial intelligence
- Artificial neural networks
- Big data
- Bioinformatics Algorithms
- Bioinformatics Databases
- Bioinformatics Modeling
- Bioinformatics Tools
- Biology Engineering
- Biostatistics: Current Trends
- Cancer Proteomics
- Chemistry of Biology
- Clinical Proteomics
- Cloud Computation
- Cluster analysis
- Comparative genomics
- Comparative proteomics
- Computational Chemistry
- Computational Sciences
- Computational drug design
- Computer Science
- Computer-aided design (CAD)
- Current Proteomics
- Data Mining Current Research
- Data algorithms
- Data mining applications in genomics
- Data mining applications in proteomics
- Data mining in drug discovery
- Data mining tools
- Data modelling and intellegence
- Data warehousing
- Ethics in Synthetic Biology
- Evolution of social network
- Evolutionary Optimisation
- Evolutionary algorithm
- Evolutionary algorithm in datamining
- Evolutionary computation
- Evolutionary science
- Experimental Physics
- Findings on Machine Learning
- Gene Synthesis
- Genome annotation
- Genomic data mining
- Genomic data warehousing
- Handover
- Human Proteome Project Applications
- Hybrid soft computing
- Industrial Biotechnology
- Knowledge modelling
- Machine Learninng
- Mapping of genomes
- Mass Spectrometry in Proteomics
- Mathematical Modeling
- Mathematics for Computer Science
- Meta genomics
- Microarray Proteomics
- Models of Science
- Molecular and Cellular Proteomics
- Multi Objective Programming
- Neural Network
- Ontology Engineering
- P4 medicine
- Physics Models
- Protein Sequence Analysis
- Proteogenomics
- Proteome Profiling
- Proteomic Analysis
- Proteomic Biomarkers
- Proteomics Clinical Applications
- Proteomics Research
- Proteomics Science
- Proteomics data warehousing
- Python for Bioinformatics
- Quantitative Proteomics
- Robotics Research
- Scientific Computing
- Simulation Computer Science
- Soft Computing
- Statistical data mining
- Studies on Computational Biology
- Swarm Robotics
- Swarm intelligence
- Synthetic Biology
- Synthetic Biology medicine
- Synthetic Biotechnology
- Synthetic Genomics
- Synthetic biology drugs
- Systems Biology
- Technologies in Computer Science
- Theoretical Chemistry
- Theoretical Computer Science
- Theoretical Issues in Ergonomics Science
- Theoretical Methods
- Theoretical and Applied Science

- International Conference on
**Metabolomics**and**Systems Biology**

June 11-12, 2018 London, UK - International Conference on
**Computational Biology**and**Bioinformatics**

Sep 05-06 2018 Tokyo, Japan - International Conference on Advancements in
**Bioinformatics**and**Drug Discovery**

November 26-27, 2018 Dublin, Ireland

- Total views:
**12624** - [From(publication date):

July-2015 - May 23, 2018] - Breakdown by view type
- HTML page views :
**8782** - PDF downloads :
**3842**

Peer Reviewed Journals

International Conferences 2018-19