Next Article in Journal
Numerical and Experimental Investigations of Humeral Greater Tuberosity Fractures with Plate Fixation under Different Shoulder Rehabilitation Activities
Next Article in Special Issue
Comparison of Spine–Pelvis Kinematics Variability during Sit-to-Stand and Stand-to-Sit in People with & without Chronic Low Back Pain: A Vector Coding and Statistical Parametric Mapping Approach
Previous Article in Journal
How Useful Are Strain Rates for Estimating the Long-Term Spatial Distribution of Earthquakes?
Previous Article in Special Issue
The Relationship between Personality and Postural Control in Young Adults—A Pilot Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Analytical Model to Predict Foot Sole Temperature: Implications to Insole Design for Physical Activity in Sport and Exercise

1
Department of Mechanics, Marvdasht Branch, Islamic Azad University, Marvdasht 73711-13119, Iran
2
Centre for Biomechanics and Rehabilitation Technologies, Staffordshire University, Leek Road, Stoke-on-Trent ST4 2DF, UK
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(13), 6806; https://doi.org/10.3390/app12136806
Submission received: 27 May 2022 / Revised: 20 June 2022 / Accepted: 28 June 2022 / Published: 5 July 2022
(This article belongs to the Special Issue Biomechanics in Sport Performance and Injury Preventing)

Abstract

:
Foot sole temperature, besides its importance in thermal comfort, can be considered an important factor in identifying tissue injuries due to heavy activities or diseases. Hyperthermia, which is a raise in the foot temperature, increases the risk of diabetic ulcers considerably. In this study, a model is proposed to predict the foot sole temperature with acceptable accuracy. This model for the first time considers both the thermal and mechanical properties of the shoe sole, the intensity of the activity, the ambient condition, and sweating, which are involved in the thermal interaction between the sole of the foot and footwear. Furthermore, the proposed model provides the opportunity to estimate the contributions of different parameters in foot thermal regulation by describing the interaction of activity, duration, and intensity as well as sweating in influencing the foot sole temperature. In doing so it takes into account the relative importance of heat capacitance and the thermal conductivity. The results of this study revealed that sweating is not as effective in cooling the ball area of the foot while it is the principal contributor to thermal regulation in the arch area. The model also showed the importance of trapped air in keeping the foot warm, especially in cold conditions. Based on the simulation results, in selecting the shoe sole, and in addition to the conductivity, the thermal capacity of the sole of the shoe needs to be considered. The developed analytical model allowed the investigation of the contribution of all the involved parameters in foot thermal regulation and has shown that a different foot temperature can be achieved when the amount of material versus air is changed in the insole design. This can have practical implications in the insole design for a variety of conditions such as hypo and hyper-thermia in physical activities in sports and exercise settings.

1. Introduction

In the design of footwear, usually, aesthetics and comfort are the main important factors considered, and thermal characteristics are not commonly taken into account. Hence, limited studies can be found that have a focus on the thermal characteristics of the footwear and the foot thermal regulation in a shoe [1,2,3,4,5]. The majority of the studies about footwear heat transfer deal with keeping feet warm in cold conditions [1,4,6,7,8,9]. Generally, feet are kept warm using thermal insulation or an electric heater pad to inhibit or compensate for the heat loss from the foot sole to the environment. However, the other aspect of foot thermal regulation involves cooling. The human body’s mechanical efficiency is not more than 20~25% [10] and the rest of the energy dissipates to the environment in the form of heat. Consequently, the foot as a body thermal radiator [11] is responsible for dissipating heat. It was observed that in vigorous physical activities such as running, the temperature of the midsole can reach 50 °C or higher during a summer daytime [12].
Cooling feet is important for two reasons: first, thermal comfort is related mainly to skin temperature [13], and for feet, this comfort depends more on skin temperature rather than moisture [14]. The second and even more important reason, is that excess foot sole temperature is associated with diabetic foot ulcers [15,16,17].
Few experimental methods to measure the foot temperature and footwear thermal insulation assessment exist [2,3,5,18,19]. To assess the footwear thermal insulation, a plastic bag filled with hot water is usually placed inside the shoe to occupy the shoe space perfectly, and then, the temperatures of the shoe’s different parts are measured during the time [5,20]. Some other experiments involve measuring the temperature in situ, by thermocouples during walking or running [2,3,18,21,22,23]. A comprehensive study in this field was conducted as a series of experiments by Shimazaki et al. [2,3,23,24].
Despite the relatively developed models to assess the mechanical behavior of footwear [25,26,27], there is currently no accepted model to justify the experimental results and predict the foot sole temperature. Recently, Nemati et al. [28] proposed an analytical method to predict the foot sole temperature during running and jogging at different speeds. In their model, they accounted for both the thermal and mechanical properties of the shoe sole. The model was then validated against experimental data and found to have a very good agreement between the measured foot sole temperatures and the predicted values. That model enables the designer to compare, for example, the softness of EVA12 and the superiority of the thermal properties of EVA08. The model has an analytical solution that can be used to accurately predict the arch area temperature for the initial duration of activity as up to 20 min during slow walking at 3 km/h and for up to 15 min for running at 9 km/h. However, after these times, the model over-predicts the foot sole temperature. Nemati et al. claimed that this may be due to the fact that the effect of sweating was not considered in their study [28]. Since the foot is capsulated inside the shoe, it is not possible to exchange heat to ambient environment by radiation or convection, efficiently. Therefore, sweating is the only effective temperature-controlling mechanism which relies on latent heat. While sweating is an effective factor in cooling feet, further care is required in modeling to assess the effect of sweating on the foot temperature.
The main aim of this study was to develop a model which can predict the temperature at the ball area and arch area in the first instance and to assess the accuracy of the model in predicting the foot sole temperature.

2. Model Development

The presented mathematical model is based on the following experimental scenario in which candidates wear the shoes without socks and stand still for a long time (around 10 min). At this time, the foot plantar surface and shoe sole reach thermal equilibrium, i.e., the shoe sole and foot plantar skin have the same temperature at their interface. Then, jogging commences for around 30 min. More detail of the experiment can be found in Section 3.
Based on the above, a transient conductive heat transfer equation can be considered for the shoe sole as follows [29]:
k s d 2 T d x 2 + q G e n = ρ s c p s d T d τ
where, ρ s and c p s are shoe sole density and specific heat properties, respectively. q G e n is the average heat generation rate per unit volume of the shoe sole due to viscous work during the stance phase of gait (Section 2.2). T is the temperature at any point along the shoe sole thickness at any time, τ .
Since the person has started the exercise from a standstill condition (please refer to Equation (12) in Section 2.3), therefore, there is a linear temperature distribution as the initial condition of Equation (1), i.e.,:
T ( x ,   τ = 0 ) = ( T T f 0 ) x l + T f 0 , 0 x l
The boundary conditions at the inner surface of the shoe sole and the outer surface of the shoe sole are, respectively:
[ q m e t ( τ ) q s w = k s d T d x ] x = 0
[ h ( T T ) = k s d T d x ] x = l
where, as discussed earlier, q m e t and h are the metabolism heat flux as a function of time and convective heat transfer coefficient, respectively, that will be discussed later. q s w is the evaporative cooling rate due to sweating.

2.1. Convective Heat Transfer from Shoe Sole to Air

During walking/jogging, the outsole of the shoe is mainly cooled by air, forced convective heat transfer ( q c o n v in Figure 1) and the air average velocity over the shoe is equal to the person walking/jogging velocity (i.e. walking or jogging under zero wind condition). Considering, the outsole as a flat surface, the average convective heat transfer coefficient (h) can be calculated as [30]:
h = 0.037 k a L   R e 4 / 5 P r 1 / 3
where L is the shoe outsole length. In the above equation:
R e = ρ a · V · L μ
P r = μ a · c p a k a

2.2. Heat Generation in Shoe Sole due to Periodic Loading–Unloading Condition

During the jogging, the shoe sole is periodically under the loading-unloading condition. The shoe sole is often constructed from viscoelastic materials (e.g., Polyurethane Elastomers, EVA foam, etc.) because of their shock-absorbent behavior [24,31]. For jogging activity, under the harmonic loading–unloading condition, viscous heat generates inside the shoe sole ( q G e n ) .
For this harmonic loading as shown in Figure 2, the volumetric generated heat during one period of the harmonic strain excitation is [32]:
U d = π ε 0 2 E
ε 0 = σ 0 π ( E ) 2 + ( E ) 2
E and E are loss modulus and storage modulus, respectively. In addition, ε 0 and σ 0 and are strain and stress amplitudes, respectively.
Walking is also a harmonic motion. An example of the time-dependent gait motion pattern is presented in Figure 3, in which, the contact force (normalized by the person’s weight) and the contact area for walking speed V = 3 km/h can be considered. It is known that the maximum contact force value depends on the gait speed and it is higher at higher velocities. However, the average contact area is independent of gait speed [2]. Furthermore, the maximum contact force can be beyond the body’s weight at foot landing during running [33].
The stress amplitude, σ 0 , is equal to the maximum contact force divided by the average contact area. The average contact area is defined as the average foot sole area that is in contact with the shoe sole during walking.
Considering the above explanation, the volumetric generated heat can be calculated by dividing the volumetric generated heat (Equation (9)) by the period time ( τ S t a n c e + τ S w i n g ) :
q G e n = U d 2 ( τ S t a n c e + τ S w i n g ) . A c o n t A
For a harmonic loading, (Figure 2) a full period of a harmonic strain excitation includes both compression and expansion, while during walking, as shown in Figure 3, only the compressive loading is included (one-half of the harmonic cycle). Consequently, factor 2 appears in the denominator of Equation (10). Moreover, based on Figure 3, since the contact area, A c o n t , is different from the outsole area (the heat transfer area), Equation (10) is multiplied by A c o n t A .

2.3. Metabolism Rate Assessment during Walking/Jogging

The described experiment [3] was comprised of two steps. After wearing shoes and 10 min of rest, participants walked for 30 min. So, at the end of the first 10 min, it can be expected that the foot and shoe have reached a thermal equilibrium (Figure 4). In the absence of any mechanical work due to walking, and for this steady-state condition, the governing equation (Equation (1)) is reduced to Equation (11) with the following boundary condition.
k s d 2 T d T 2 = 0 T ( x = 0 ) = T f 0 T ( x = l ) = T G
where k s is the shoe sole thermal conductivity and T is temperature. l is the shoe sole thickness. T f 0 and T G are the average foot sole temperature and ground temperature, respectively. Equation (11) is an ordinary differential equation with a solution in the form of:
T ( x ) = ( T T f 0 ) x l + T f 0 , 0 x l
This would be the initial condition at the begging of exercise.
Because of the thermal equilibrium and steady-state condition, the dissipated heat flux from the foot sole is equal to (Figure 4):
q m e t 0 = k s T f 0 T G l
Starting from the rest condition to walking, the metabolism rate increases sharply by factor F:
q m e t ( τ ) = F . q m e t 0
The metabolism intensification factor (F) is a function of the gait speed and is adopted from Kipp’s experimental correlation [34].
F ( V ) = 1.3 0.009 V 2 + 0.002 V + 0.082 0.009 V r e f 2 + 0.002 V r e f + 0.082
V r e f is 0.833 m/s (3 km/h).

2.4. Latent Heat Loss due to Sweating during Walking/Jogging

The evaporation of sweat is an effective cooling mechanism for thermal regulation. Neither the sweating rate, nor the evaporation rate is uniform over the skin surface. It will be shown later that each part of foot sole has a different sweating evaporation. Cooling due to sweating is more considerable in the arc zone and is negligible in the ball of the foot zone. In the absence of more precise equations, the following equation is used to simulate the sweating–cooling effect:
q s w = h f g h e × 4.7 × 10 5 × ( P s a t P ) ( T f T f 0 )
Equation (16) was adopted from [35]. h f g = 2.43 × 10 6 J/kg is the water latent heat [36]. According to Equation (16), the sweating rate increases by increasing the temperature elevation ( T f T f 0 ) . In the above equation, P is the vapour pressure in vicinity of the ball zone (the area that sweating is considerable). T f is the foot sole skin temperature and P s a t is the saturate vapour pressure at T f . It was assumed that the vapour pressure at the skin is at the saturated pressure [37,38]. So, when the vapor pressure reaches the saturated pressure, evaporation stops. P s a t can be estimated by the Antoine equation [36]:
P s a t = 100   exp ( 18.956 4030.18 235 + T f )
However, the vapor pressure can be defined based on relative humidity, RH:
P = R H · P s a t
and Equation (16) can be represented as:
q s w = h f g h e × 4.7 × 10 5 × P s a t ( 1 R H ) ( T f T f 0 )
h e is the evaporative heat transfer coefficient, which is proportional to the square root of velocity [39]:
h e = 0.000192 × L R × V 0.5
where LR is the Lewis ratio and approximately equals 16.5 × 10 3 K/Pa.

3. Heat Transfer Modeling of a Shod Foot

The results of the current simulation were compared against the results from the literature [3]. The experiment was comprised of two stages [3]. In that experiment, first, participants wore the shoe with no socks, and after 10 min of rest, they did a jogging exercise for 30 min [3]. Foot skin temperatures were measured at various points every minute by thermocouples [3]. Participants were seven healthy adult males (1.72 ± 0.07 m height, 61.8 ± 4.7 kg weight, and 23.8 ± 4.3 years old) [3]. The footwear was a running shoe with a woven upper (Figure 5) [3]. The foot sole inside length is L = 25 cm (1 cm less than the shoe outsole length due to shoe thickness) [3]. The experiments were performed in a climate chamber with a constant air temperature of 25 °C, 58.5% RH, no solar radiation, and wind for three different speeds, i.e., 3, 6, and 9 Km/h [3]. More details can be followed in [3].
The shoe sole was made-up of EVA08, and its properties are presented in Table 1:
Thermo-physical properties of air are also presented in Table 2.
For this experiment, the average contact area is 0.006 m2 for all gait speeds and the recorded maximum contact force is reported in Table 3 for each velocity. The average bare feet temperature at the beginning of the jogging was T f 0 = 34 °C
The equations were solved numerically using the Crank–Nicolson scheme. Since the set of equations is nonlinear, they had to be solved iteratively. To validate the numerical method, the problem described in [28] was solved numerically and compared with the proposed analytical solution in [28]. The maximum considered difference was less than 0.05%. The foot sole temperature rises ( T f T f 0 ) were calculated for three speeds with and without considering the sweating evaporative heat loss. The results are compared with experimental measurements based on [3], as presented in Figure 6.
A good agreement can be observed between the experimental and numerical results (Figure 6). According to Figure 6, the maximum increase in foot sole temperature after 30 min is 6 °C for a gait speed of 3 km/h. This value is 8 °C and 11.5 °C, for 6 km/h and 9 km/h, respectively.
Figure 6 shows that for a gait speed of 3 km/h, the effect of sweating on cooling is not strong and all curves, i.e., the numerical results with and without sweating and experimental measurements for ball and arch, are approximately overlapping. However, at higher speeds, especially for 9 km/h, the effect of sweating is very considerable. For a gait speed of 9 km/h, after the first 15 min, sweating as a thermoregulation mechanism comes into play and cools the arch area. The predicted temperature-time trend with sweating much closely replicate the temperature measured under the arch area while temperatures predicted without sweating more closely match the temperature-time trend at the ball of the foot.
For a gait speed of 6 km/h, sweating-induced cooling is triggered at a later time around 20 min, while at this speed the sweating at the ball zone seems to be negligible (Figure 6b).

4. Illustrative Example

The results of this study show that sweating plays a principal role in the thermal regulation of the arch area. As an example of the application of the developed method, the effect of a porous turf-like insole design on foot temperature can be investigated. A schematic of the fibrous insole is shown in Figure 7. This insole is comprised of longitudinal fibers arranged in an aligned layout on a base. The thickness of the base is 1 mm and the length of the fiber is 10 mm. The fiber diameter is 1.2 mm (center to center distance) with the variable pitch.
The insole is to be located inside the specified shoe for medical purposes. It was assumed that the insole material is thermoplastic polyurethane with the properties mentioned in Table 4.
To consider the thermal properties of the insole, it was considered as a porous media, and properties were assumed as a volume-weighted average of air at 35 °C (Table 2) and thermoplastic polyurethane. To find the effective thermal conductivity, two different constant temperatures were assumed for the top and bottom of the insole and the 3D steady-state conduction heat transfer equation (Laplacian of temperature) was solved for two conjugated domains (air and thermoplastic polyurethane). Then, this was equated to a 1D heat transfer with the same temperature difference, the same geometry, and the same heat transfer rate to find the effective thermal conductivity. Results are presented in Table 5.
Simulations were performed at speeds of 3 and 6 km/h. To find the extreme values of temperature, sweating was omitted. Results are presented in Figure 8.
Based on Figure 8, the presence of the insole causes a sharp temperature increase at the beginning of the activity, which is mainly due to trapped air between the fibers. This shows the importance of the presence of still air in keeping the foot warm. It was interesting to note that the difference between the “without insole” condition and the insole with a pitch of 1.5 mm decreases in the course of jogging. This shows the importance of insole heat capacity ( ρ e f f c p e f f ), which prevents the temperature rise and keeps the foot sole temperature cool. So, in the design of an insole, both thermal conductivity and heat capacity must be considered.

5. Conclusions

An increase in the foot temperature raises the risk of plantar soft tissue injuries. So, the design of a comfortable shoe sole with proper thermal specifications is very important. However, the insole design is challenging as the effective parameters and the phenomena that regulated the foot temperature are complex. In this study, a comprehensive model was developed to predict the temperature of the barefoot sole in a shoe during walking and jogging. The results were validated against the experimental measurements. The effects of both the thermal and the mechanical properties of shoe soles were considered. Furthermore, the cooling effect due to sweat evaporation came into account as was presented at three different gait speeds. It was revealed that sweating is an effective parameter in foot arch zone thermal regulation while it is negligible at the ball of the foot area.
The presented method enables designers and researchers to investigate the interaction of different parameters and their impact before doing a large series of experiments. The study revealed the importance of a trapped layer of air in warming the foot sole and showed the importance of shoe sole thermal capacity besides its conductivity.
As an illustrative example, the effect of wearing a porous turf-like insole concept was investigated using the developed analytical model, where the foot sole temperature during the different levels of activities was predicted. This has shown that a different foot temperature can be achieved when the amount of material versus the air changes in the insole design. This developed model can potentially decrease the cost of experiments considerably by allowing an approximation of the insole design before producing the prototype.
As with any modelling approach, there are certain limitations that need to be taken into account. Since participants in the experiment were barefooted, the effect of sock properties such as material, thickness, and or humidity absorption can be considered in future studies.

Author Contributions

Conceptualization, H.N. and R.N.; methodology, H.N.; software, H.N.; writing—original draft preparation, H.N.; writing—review and editing, H.N and R.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

A Area
c p Specific heat
E Storage modulus
E Loss modulus
F Metabolism intensification factor
h Average convective heat transfer coefficient
h e Evaporative heat transfer coefficient
h f g Water latent heat
k Thermal conductivity
L Shoe length
LRLewis ratio
P Pressure
P r Prandtl number
q Heat flux
q Volumetric heat generation rate
R e Reynolds number
R H Relative humidity
T Temperature
l Thickness
U d Dissipated energy per unit of volume in one period
V Gait speed
x Cartesian coordinate axis
Greek letters
ε 0 Strain amplitude
μ Viscosity
ρ Density
σ 0 Stress amplitude
τ Time
Subscripts
contContact
convConvection
fFoot sole in the jogging condition
f0Foot sole in the static steady condition
GGround
genGeneration
IInsole
metMetabolism
sShoe sole
satSaturated
swSweating
Ambient

References

  1. Havenith, G.; Ouzzahra, Y.; Gerrett, N.; Faulkner, S. Assessing the lower temperature limit for comfort in footwear. In Proceedings of the 15th International Conference on Environmental Ergonomics, Queenstown, New Zealand, 11–15 February 2013. [Google Scholar]
  2. Shimazaki, Y.; Murata, M. Effect of gait on formation of thermal environment inside footwear. Appl. Ergon. 2015, 49, 55–62. [Google Scholar] [CrossRef] [PubMed]
  3. Shimazaki, Y.; Matsutani, T.; Satsumoto, Y. Evaluation of thermal formation and air ventilation inside footwear during gait: The role of gait and fitting. Appl. Ergon. 2016, 55, 234–240. [Google Scholar] [CrossRef] [PubMed]
  4. Kopitar, D.; Akalovic, J.; Skenderi, Z. Thermal resistance of leather and membranes for summer desert military footwear under different climate conditions. In Functional Textiles and Clothing; Springer: Singapore, 2019; pp. 3–9. [Google Scholar]
  5. Puszkarz, A.K.; Usupov, A. The Study of Footwear Thermal Insulation Using Thermography and the Finite Volume Method. Int. J. Thermophys. 2019, 40, 45. [Google Scholar] [CrossRef] [Green Version]
  6. Rintamäki, H.; Rissanen, S.; Mäkinen, T.; Peitso, A. Finger temperatures during military field training at 0 to −29 °C. J. Therm. Biol. 2004, 29, 857–860. [Google Scholar] [CrossRef]
  7. Covill, D.; Guan, Z.; Bailey, M.; Pope, D. Finite element analysis of the heat transfer in footwear. In The Engineering of Sport 7; Springer: Paris, France, 2008; pp. 247–254. [Google Scholar]
  8. Covill, D.; Guan, Z.; Bailey, M.; Raval, H. Development of thermal models of footwear using finite element analysis. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 2011, 225, 268–281. [Google Scholar] [CrossRef]
  9. Liu, G.; Su, Y.; Zhu, W.; Tian, M.; Li, J.; Tian, Y. Numerical simulation of heat transfer in electrically heated footwear in a severely cold environment. Build. Environ. 2022, 207, 108429. [Google Scholar] [CrossRef]
  10. James, S. Encyclopaedia of Occupational Health and Safety, 4th ed.; International Labour Office: Geneva, Switzerland, 1999. [Google Scholar]
  11. Taylor, N.A.; Machado-Moreira, C.A.; van den Heuvel, A.M.; Caldwell, J.N. Hands and feet: Physiological insulators, radiators and evaporators. Eur. J. Appl. Physiol. 2014, 114, 2037–2060. [Google Scholar] [CrossRef]
  12. Kinoshita, H.; Bates, B.T. The effect of environmental temperature on the properties of running shoes. J. Appl. Biomech. 1996, 12, 258–268. [Google Scholar] [CrossRef]
  13. Bulcao, C.F.; Frank, S.M.; Raja, S.N.; Tran, K.M.; Goldstein, D.S. Relative contribution of core and skin temperatures to thermal comfort in humans. J. Therm. Biol. 2000, 25, 147–150. [Google Scholar] [CrossRef]
  14. Arezes, P.; Neves, M.M.; Teixeira, S.; Leão, C.P.; Cunha, J. Testing thermal comfort of trekking boots: An objective and subjective evaluation. Appl. Ergon. 2013, 44, 557–565. [Google Scholar] [CrossRef]
  15. Ersen, A.; Adams, L.S.; Myers, R.T.; Hirschman, G.B.; Lavery, L.A.; Yavuz, M. Temperature Regulating Shoes for Prevention of Diabetic Foot Ulcers. Diabetes 2018, 67, 41-LB. [Google Scholar] [CrossRef]
  16. Houghton, V.J.; Bower, V.M.; Chant, D.C. Is an increase in skin temperature predictive of neuropathic foot ulceration in people with diabetes? A systematic review and meta-analysis. J. Foot Ankle Res. 2013, 6, 31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Yavuz, M.; Ersen, A.; Monga, A.; Lavery, L.A.; Garrett, A.G.; Salem, Y.; Hirschman, G.B.; Myers, R. Temperature-and Pressure-Regulating Insoles for Prevention of Diabetic Foot Ulcers. J. Foot Ankle Surg. 2020, 59, 685–688. [Google Scholar] [CrossRef] [PubMed]
  18. Taiar, R.; Rebay, M.; Vannozzi, G.; Sanna, G.; Cappozzo, A. Evolution of the in-Shoe Temperature during Walking and Running; Andreas, H., Ed.; ACTA Press: Anaheim, CA, USA, 2013; pp. 130–133. [Google Scholar]
  19. Carbonell, L.; Priego Quesada, J.I.; Retorta, P.; Benimeli, M.; Cibrián Ortiz De Anda, R.M.; Salvador Palmer, R.; González Peña, R.J.; Galindo, C.; Pino Almero, L.; Blasco, M.C. Thermographic quantitative variables for diabetic foot assessment: Preliminary results. Comput. Methods Biomech. Biomed. Eng. Imaging Vis. 2019, 7, 660–666. [Google Scholar] [CrossRef]
  20. Dessing, O.; Jansen, A.J.; Leihitu, C.; Overhage, D. Experimental study of heat dissipation in indoor sports shoes. Procedia Eng. 2014, 72, 575–580. [Google Scholar] [CrossRef] [Green Version]
  21. Rebay, M.; Arfaoui, A.; Leopol, A.; Perin, J.; Taiar, R. Heat transfer in athletic shoes during the running. In Proceedings of the 5th IASME/WSEAS International Conference on Heat Transfer, Thermal Engineering and Environment, Athens, Greece, 25–27 August 2007. [Google Scholar]
  22. Irzmańska, E. The microclimate in protective fire fighter footwear: Foot temperature and air temperature and relative humidity. Autex. Res. J. 2016, 16, 75–79. [Google Scholar] [CrossRef] [Green Version]
  23. Shimazaki, Y.; Aisaka, K. Novel Thermal Analysis Model of the Foot-Shoe Sole Interface during Gait Motion. Multidiscip. Digit. Publ. Inst. Proc. 2018, 2, 278. [Google Scholar]
  24. Shimazaki, Y.; Nozu, S.; Inoue, T. Shock-absorption properties of functionally graded EVA laminates for footwear design. Polym. Test. 2016, 54, 98–103. [Google Scholar] [CrossRef]
  25. Naemi, R.; Chockalingam, N. Development of a method for quantifying the midsole reaction model parameters. Comput. Methods Biomech. Biomed. Eng. 2013, 16, 1273–1277. [Google Scholar] [CrossRef]
  26. Naemi, R.; Chockalingam, N. Mathematical models to assess foot-ground interaction: An overview. Med. Sci. Sports Exerc. 2013, 45, 1524–1533. [Google Scholar] [CrossRef]
  27. Naemi, R.; Chatzistergos, P.E.; Chockalingam, N. A mathematical method for quantifying in vivo mechanical behaviour of heel pad under dynamic load. Med. Biol. Eng. Comput. 2016, 54, 341–350. [Google Scholar] [CrossRef] [PubMed]
  28. Nemati, H.; Moghimi, M.; Naemi, R. A mathematical model to investigate heat transfer in footwear during walking and jogging. J. Therm. Biol. 2021, 97, 102778. [Google Scholar] [CrossRef] [PubMed]
  29. Hahn, D.W.; Özisik, M.N. Heat Conduction; John Wiley & Sons: New York, NY, USA, 2012. [Google Scholar]
  30. Bergman, T.L.; Incropera, F.P.; Lavine, A.S.; DeWitt, D.P. Introduction to Heat Transfer; John Wiley & Sons: New York, NY, USA, 2011. [Google Scholar]
  31. Even-Tzur, N.; Weisz, E.; Hirsch-Falk, Y.; Gefen, A. Role of EVA viscoelastic properties in the protective performance of a sport shoe: Computational studies. Bio-Med. Mater. Eng. 2006, 16, 289–299. [Google Scholar]
  32. Roylance, D. Engineering Viscoelasticity; Department of Materials Science and Engineering–Massachusetts Institute of Technology: Cambridge, MA, USA, 2001; pp. 1–37. Available online: http://web.mit.edu/course/3/3.11/www/modules/visco.pdf (accessed on 26 May 2022).
  33. Chuckpaiwong, B.; Nunley, J.A.; Mall, N.A.; Queen, R.M. The effect of foot type on in-shoe plantar pressure during walking and running. Gait Posture 2008, 28, 405–411. [Google Scholar] [CrossRef] [PubMed]
  34. Kipp, S. Why Does Metabolic Rate Increase Curvilinearly with Running Velocity. Master’s Thesis, University of Colorado, Boulder, CO, USA, 2017. [Google Scholar]
  35. Gagge, A.P. An effective temperature scale based on a simple model of human physiological regulatory response. Ashrae Trans. 1971, 77, 247–262. [Google Scholar]
  36. Parsons, K. Human Thermal Environments: The Effects of Hot, Moderate, and Cold Environments on Human Health, Comfort and Performance; CRC Press: New York, NY, USA, 2007. [Google Scholar]
  37. Arens, E.A.; Zhang, H. The skin’s role in human thermoregulation and comfort. In Thermal and Moisture Transport in Fibrous Materials; Pan, N., Gibson, P., Eds.; Woodhead Publishing Ltd.: Cambridge, UK, 2006; pp. 560–602. [Google Scholar]
  38. Foda, E.; Sirén, K. A new approach using the Pierce two-node model for different body parts. Int. J. Biometeorol. 2011, 55, 519–532. [Google Scholar] [CrossRef] [PubMed]
  39. Gonzalez, R.R.; Cheuvront, S.N.; Montain, S.J.; Goodman, D.A.; Blanchard, L.A.; Berglund, L.G.; Sawka, M.N. Expanded prediction equations of human sweat loss and water needs. J. Appl. Physiol. 2009, 107, 379–388. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Different forms of heat transfer to/from shoe sole.
Figure 1. Different forms of heat transfer to/from shoe sole.
Applsci 12 06806 g001
Figure 2. Harmonic strain excitation and stress response.
Figure 2. Harmonic strain excitation and stress response.
Applsci 12 06806 g002
Figure 3. Shoe contact area and contact force for V = 3.0 km/h modified from [2].
Figure 3. Shoe contact area and contact force for V = 3.0 km/h modified from [2].
Applsci 12 06806 g003
Figure 4. Thermal conditions for 10 min rest in the experiment.
Figure 4. Thermal conditions for 10 min rest in the experiment.
Applsci 12 06806 g004
Figure 5. The running shoe used for the experiment (adapted from [3]).
Figure 5. The running shoe used for the experiment (adapted from [3]).
Applsci 12 06806 g005
Figure 6. Comparison of foot sole temperature (numerical simulation and experimental measurements). (a) Gait speed 3 km/h; (b) Gait speed 6 km/h; and (c) Gait speed 9 km/h.
Figure 6. Comparison of foot sole temperature (numerical simulation and experimental measurements). (a) Gait speed 3 km/h; (b) Gait speed 6 km/h; and (c) Gait speed 9 km/h.
Applsci 12 06806 g006
Figure 7. Schematic of the porous turf-like insole.
Figure 7. Schematic of the porous turf-like insole.
Applsci 12 06806 g007
Figure 8. The foot sole temperature during shod gait at 3 and 6 km/h, using two porous turf-like insole designs with pitch of 2 mm and 1.5 mm and comparing with no insole condition.
Figure 8. The foot sole temperature during shod gait at 3 and 6 km/h, using two porous turf-like insole designs with pitch of 2 mm and 1.5 mm and comparing with no insole condition.
Applsci 12 06806 g008
Table 1. EVA08 thermal–mechanical properties based on what was reported in [23].
Table 1. EVA08 thermal–mechanical properties based on what was reported in [23].
c p s
( J · k g 1 · K 1 )
k s
( W · m 1 · K 1 )
ρ s
( k g · m 3 )
E
(Pa)
E
(Pa)
Hardness
25850.11240 1.55 × 10 6 12.6 × 10 6 A45
Table 2. Air thermo-physical properties at 25 °C based on what was reported in [30].
Table 2. Air thermo-physical properties at 25 °C based on what was reported in [30].
Temperature
(°C)
c p a
( J · k g 1 · K 1 )
k a
( W · m 1 · K 1 )
ρ a
( k g · m 3 )
μ a
( Pa · s )
251006.960.02611.171 1.84 × 10 5
35 *10070.02691.135---
* is used in Section 4.
Table 3. Maximum contact force in percent of body weight (%B.W) for different gait speeds as reported in [2].
Table 3. Maximum contact force in percent of body weight (%B.W) for different gait speeds as reported in [2].
Speed
(km/h)
369
Max. Contact force
%B.W.
100110170
Table 4. Thermoplastic polyurethane properties.
Table 4. Thermoplastic polyurethane properties.
c p I
( J · k g 1 · K 1 )
k I
( W · m 1 · K 1 )
ρ I
( k g · m 3 )
17000.2581130
Table 5. Effective properties of the insole.
Table 5. Effective properties of the insole.
Pitch
(mm)
c p e f f
( J · k g 1 · K 1 )
k e f f
( W · m 1 · K 1 )
ρ e f f
( k g · m 3 )
1.51386.70.145619.6
21248.10.096393.9
2.51184.00.073289.5
31149.20.060232.7
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nemati, H.; Naemi, R. An Analytical Model to Predict Foot Sole Temperature: Implications to Insole Design for Physical Activity in Sport and Exercise. Appl. Sci. 2022, 12, 6806. https://doi.org/10.3390/app12136806

AMA Style

Nemati H, Naemi R. An Analytical Model to Predict Foot Sole Temperature: Implications to Insole Design for Physical Activity in Sport and Exercise. Applied Sciences. 2022; 12(13):6806. https://doi.org/10.3390/app12136806

Chicago/Turabian Style

Nemati, Hossain, and Roozbeh Naemi. 2022. "An Analytical Model to Predict Foot Sole Temperature: Implications to Insole Design for Physical Activity in Sport and Exercise" Applied Sciences 12, no. 13: 6806. https://doi.org/10.3390/app12136806

APA Style

Nemati, H., & Naemi, R. (2022). An Analytical Model to Predict Foot Sole Temperature: Implications to Insole Design for Physical Activity in Sport and Exercise. Applied Sciences, 12(13), 6806. https://doi.org/10.3390/app12136806

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop