- 1Key Laboratory of Bionic Engineering, Jilin University, Changchun, China
- 2Orthopaedic Medical Center, The Second Hospital of Jilin University, Changchun, China
- 3School of Mechanical, Aerospace and Civil Engineering, University of Manchester, Manchester, United Kingdom
Plantar fascia plays an important role in human foot biomechanics; however, the morphology and mechanical properties of plantar fascia in patients with flexible flatfoot are unknown. In this study, 15 flexible flatfeet were studied, each plantar fascia was divided into 12 positions, and the morphologies and mechanical properties in the 12 positions were measured in vivo with B-mode ultrasound and shear wave elastography (SWE). Peak pressures under the first to fifth metatarsal heads (MH) were measured with FreeStep. Statistical analysis included 95% confidence interval, intragroup correlation coefficient (ICC1,1), one-way analysis of variance (one-way ANOVA), and least significant difference. The results showed that thickness and Young’s modulus of plantar fascia were the largest at the proximal fascia (PF) and decreased gradually from the proximal end to the distal end. Among the five distal branches (DB) of the fascia, the thickness and Young’s modulus of the second and third DB were larger. The peak pressures were also higher under the second and third MH. This study found a gradient distribution in that the thickness and Young’s modulus gradient decreased from the proximal end to the distal end of plantar fascia in the longitudinal arch of flexible flatfeet. In the transverse arch, the thickness and Young’s modulus under the second and third DB were larger than those under the other three DB in flexible flatfoot, and the peak pressures under the second and third MH were also larger than those under the other three MH in patients with flexible flatfoot. These findings deepen our understanding of the changes of biomechanical properties and may be meaningful for the study of pathological mechanisms and therapy for flexible flatfoot.
Introduction
Plantar fascia is a ligament that attaches the calcaneus to metatarsals (Orchard, 2012). It plays an important role in passive force transmission (Stecco et al., 2013). Its main task is to stabilize the arch of the foot and reduce the influence of ground reaction force on metatarsal heads (MH) and the longitudinal foot arch (Hicks, 1954; Ker et al., 1987; McKeon et al., 2015). There is a close relationship between plantar fascia and foot function, and studies have shown that when plantar fascia changes, it will produce clinical problems, for example, heel pain (Wearing et al., 2006). Thus, research on plantar fascia has a broad interest.
During the past decades, numerous studies on plantar fascia have been conducted. Guo et al. (2018) found a certain relationship between the mechanical tension of plantar fascia and fiber morphology. Chen et al. (2019a,b) found that people who used forefoot strike were more likely to suffer from plantar fasciitis. Tas and Cetin (2019a) focused on the relationship between plantar pressure distribution and the morphology and mechanical properties of plantar fascia. Welte et al. (2021) revealed the effect of plantar fascia extensibility on the windlass mechanism of plantar fascia. Wang et al. (2019) illustrated the morphology and mechanical properties of plantar fascia in normal feet. These studies strengthen the understanding of the mechanical properties of plantar fascia in normal feet. However, to the author’s knowledge, the morphology and mechanical properties of the whole plantar fascia of flexible flatfeet have not been reported to date.
Flatfoot is a common foot posture abnormality, with the highest incidence of 78% (Sung, 2016), and is characterized by a low medial longitudinal arch (Pehlivan et al., 2009). Flatfoot can be divided into rigid flatfoot and flexible flatfoot. Rigid flatfoot means that the medial longitudinal arch is always missing in both load-bearing and nonload-bearing positions. Flexible flatfoot means that the medial longitudinal arch is missing only in the load-bearing position, while in the nonload-bearing position, it is the same as that of a normal foot (Carr et al., 2016). The abnormal structural changes of the flexible flatfoot under load will gradually lead to changes in the morphology and mechanical properties of the plantar fascia, which may lead to plantar fasciitis and other diseases. The changes in the morphology and mechanical properties of plantar fascia will in turn affect the foot kinematics of patients with flatfoot, resulting in clinical symptoms such as patellar tendinopathy and medial tibial stress syndrome (Kohls-Gatzoulis et al., 2004; Van der Worp et al., 2011; Hamstra-Wright et al., 2015). Studies have shown that the potential cause of plantar fasciitis is the abnormal morphology and mechanical properties of plantar fascia (Wearing et al., 2006; Wu et al., 2011).
Therefore, the objective of this study was to investigate the morphology and mechanical properties of plantar fascia of patients with flexible flatfoot by B-mode ultrasound and shear wave elastography (SWE) in vivo. A comprehensive analysis was conducted combined with plantar pressure measurement. The results of the study may provide a meaningful reference and basis for analysis of the pathological mechanism and rehabilitation in patients with flexible flatfoot as well as more accurate definitions for foot finite element models.
Methods
Ethics Statement
This study was based on the principles outlined in the Helsinki Declaration, which was approved by the Ethics Committee of the Second Hospital of Jilin University (No. 2020085). All volunteers who participated in the study signed written informed consent agreements.
Selection of Research Subjects
The subjects of this experiment were patients with flexible flatfeet. They had the typical characteristics in that the medial longitudinal arch was missing only in the load-bearing position, while in the nonload-bearing position, it was the same as that of a normal foot (Carr et al., 2016). An intelligent scanner was employed to confirm the diagnosis and severity of flatfeet, of which the diagnostic principle was the arch index proposed by Cavanagh et al. (1987). The arch index was widely accepted and adopted (Wearing et al., 2004; Wong et al., 2012; Nirenberg et al., 2020; Wang et al., 2020). The inclusion criteria were as follows: 1) healthy male, 20–30 years old; 2) the diagnosis being flexible flatfeet; and 3) no history of other foot diseases. The exclusion criteria were as follows: 1) rigid flatfeet; 2) a history of foot trauma or surgery; 3) presence of systemic diseases that may affect plantar fascia, such as rheumatoid arthritis, diabetes, and gout; and 4) the presence of diseases that affect local plantar fascia, such as calcaneal spur or nodular fasciitis and plantar fibromatosis. Finally, 10 volunteers with 15 flexible flatfeet were included, and the basic characteristics of the volunteers were age, 26.2 ± 1.6 years; weight, 65.2 ± 2.2 kg; height, 175.2 ± 2.7 cm.
Test Device and Procedure
The subjects were asked to avoid intense sports 1 week before the test. B-mode and SWE mode of an Aixplorer ultrasonic scanner (Aixplorer ultrasonic imager, Aix-en-Provence, France) were used to measure the thickness and Young’s modulus of plantar fascia, respectively. The linear transducer frequency was 10–2 MHz for this study. The sampling depth was adjusted according to the positions of plantar fascia. It was set at 1.5–2.5 cm to include the whole plantar fascia, and the mechanical index was 1.0 in this study. During measurement, each subject lay prone on the examination bed, with the lower limbs straight and the feet hanging naturally (Haen et al., 2017) on the edge of the examination bed (Figure 1). The upper body and legs were relaxed.
In order to observe the entire changes in the plantar fascia, it was divided into four main regions: proximal fascia (PF), middle fascia (MF), five branches of fascia (BF1-BF5), and five distal branches (DB1-DB5), 12 positions in total (Figure 2). The PF was measured at a point 1 cm away from the insertion to the calcaneus. The location of the five DBs was defined as the farthest end where the plantar fascia has not been fused with joint capsule. The ultrasonic transducer was parallel to the plantar fascia, and the thickness of plantar fascia was measured in the middle of every position. Subsequently, the elastic measurement via SWE was performed. The width of the square-shaped elastography window (region of interest, ROI) was as large as possible, and the height was set to include the complete plantar fascia. Q-BoxTMTrace was used to measure Young’s modulus (maximum, minimum, and average, in kPa) of plantar fascia with a length of 1 cm at each position, and Young’s modulus scale was adjusted to 0–600 kPa (Wang et al., 2019). Additionally, the mean Young’s modulus value was used for the data analysis in the study. At each position, Young’s modulus and thickness of plantar fascia were measured three times.
Plantar Pressure Measurement
FreeStep (Sensor Medica, Italy) was employed to detect the plantar pressure of the subjects during level walking. Subjects were requested to walk normally, without rushing, acceleration, or deceleration. Data were collected barefoot at a self-selected speed (Rao et al., 2011; Hillstrom et al., 2013) along a 2 m walkway, and the walking velocity was 1.33 ± 0.97 m/s. The peak pressures under the first to fifth MH were measured during the push-off stage.
Statistical Analysis
IBM Statistical Package for the Social Sciences (SPSS) statistical software version 26.0 (SPSS Inc., Chicago, IL, United States) was used to analyze all the data. The 95% confidence interval (95% CI) and intragroup correlation coefficient (ICC1,1) were used to measure and evaluate the reliability of plantar fascia thickness and Young’s modulus. Generally, the values of ICC1,1 in the ranges of 0–0.40, 0.41–0.6, 0.61–0.79, and 0.8–1.0, respectively, indicate poor, medium, good, and excellent reliability. At the same time, the one-way analysis of variance (one-way ANOVA) was used to compare the differences between different positions of plantar fascia. If the result of one-way ANOVA was p < 0.05, least significant difference was used to compare the differences between every two positions of plantar fascia. For least significant difference, P values we used had been corrected by the number of pairwise comparisons. Statistical difference was defined as p < 0.05. In order to better understand the spatial distribution in thickness and Young’s modulus of plantar fascia, an exponential function (first-order exponential decay) was used to fit and analyze the variation trend of plantar fascia from the calcaneal to the five DB.
Results
Intragroup Correlation Results of Thickness and Young’s Modulus
The intragroup correlation results of the thickness and Young’s modulus of the 15 flatfeet are listed in Table 1. The ICC1,1 ranged from 0.976 to 0.995 and the corresponding 95% CI was 0.938, 0.998 for thickness of plantar fascia. The ICC1,1 ranged from 0.985 to 0.999 and the 95% CI was 0.961, 1.000 for Young’s modulus of plantar fascia.
Distribution Pattern of Thickness and Young’s Modulus of Plantar Fascia
The results of thickness and Young’s modulus of plantar fascia of 15 flexible flatfeet are shown in Figure 3. The results showed that both the thickness and Young’s modulus of plantar fascia decreased gradually from the proximal end to the distal end. Among the five DB, the thickness and Young’s modulus of the second and third branches were larger than the other three.
The one-way ANOVA was used to compare the differences between different positions of plantar fascia. If the result of one-way ANOVA was p < 0.05, least significant difference was used to compare the differences between every two positions of plantar fascia. The one-way ANOVA results showed that the differences in thickness and Young’s modulus between different positions were statistically significant (p < 0.05). Least significant difference results showed that, in terms of plantar fascia thickness, PF > MF > all the five BFs > all the five DBs. Among the five DBs, DB2 and DB3 > DB1 and DB4 > DB5. The differences were statistically significant (p < 0.05). There was no statistical difference between DB2 and DB3, and there was also no statistical difference between DB1 and DB4 (Table 2). For Young’s modulus, PF > MF > all the five BF > the corresponding position of DB. Among the five DBs, DB2 > DB4 and DB5; DB3 > DB1 and DB4 and DB5; DB1 > DB5. All the differences were statistically significant (p < 0.05). There was no statistical difference between DB2 and DB3, and no statistical difference was found between DB2 and DB1. There was also no statistical difference between DB4 and DB5 (Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727940/fbioe-09-727940-HTML/image_m/fbioe-09-727940-t002.jpg)
TABLE 2. P value of least significant difference results between different positions in thickness of plantar fascia.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727940/fbioe-09-727940-HTML/image_m/fbioe-09-727940-t003.jpg)
TABLE 3. P value of least significant difference results between different positions in Young’s modulus of plantar fascia.
Peak Pressure Distribution Under Five MHs
The peak pressure under five MHs of 15 flexible flatfeet is shown in Figure 4. The pressures under the second and third MH were higher than those under the other three MH, and the differences were statistically significant (p < 0.05) (Table 4). This distribution pattern is similar to the thickness and Young’s modulus in the five DBs.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727940/fbioe-09-727940-HTML/image_m/fbioe-09-727940-g004.jpg)
FIGURE 4. The peak pressure under the five metatarsal heads (MH1–MH5) (A), the thickness (B), and Young’s modulus (C) in the five distal branches of the plantar fascia in 15 flatfeet.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727940/fbioe-09-727940-HTML/image_m/fbioe-09-727940-t004.jpg)
TABLE 4. P value of least significant difference results between different positions in peak pressure.
Spatial Distribution of Plantar Fascia Thickness and Young’s Modulus
In order to better understand the spatial distribution in thickness and Young’s modulus of plantar fascia, an exponential function (first-order exponential decay) was used to fit and analyze the variation trend of plantar fascia from the calcaneal to the five DB.
The spatial distribution of plantar fascia thickness and Young’s modulus of foot #1 is shown in Figure 5. The results showed that the thickness and Young’s modulus of plantar fascia were the largest at the calcaneus tubercle, and the thickness and Young’s modulus of five fascial bundles gradually decreased as plantar fascia extended from the calcaneus to the five toes. The spatial distribution of thickness and Young’s modulus in the other 14 flexible flatfeet also showed a similar tendency. The thickness and Young’s modulus of plantar fascia of 15 flatfeet at PF and five DBs are shown in Figure 6.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727940/fbioe-09-727940-HTML/image_m/fbioe-09-727940-g005.jpg)
FIGURE 5. The curve of the thickness and Young’s modulus of the plantar fascia from the calcaneus to the five distal branches in foot #1: (A) thickness of the first branch, (B) Young’s modulus of the first branch, (C) thickness of the second branch, (D) Young’s modulus of the second branch, (E) thickness of the third branch, (F) Young’s modulus of the third branch, (G) thickness of the forth branch, (H) Young’s modulus of the forth branch, (I) thickness of the fifth branch, and (J) Young’s modulus of the fifth branch.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727940/fbioe-09-727940-HTML/image_m/fbioe-09-727940-g006.jpg)
FIGURE 6. Distribution of plantar fascia thickness (A) and Young’s modulus (B) at proximal fascia (PF) and the five distal branches in 15 flatfeet.
Discussion
This study investigated the morphology and mechanical properties of plantar fascia of patients with flexible flatfoot by B-mode ultrasound and ultrasonic elastography in vivo. A comprehensive analysis was conducted combined with plantar pressure measurements.
In order to evaluate the accuracy of the data, the repeatability of the thickness and Young’s modulus data was analyzed in all 15 flexible flatfeet. The results showed that all the values of ICC1,1 were more than 0.9, which indicated that the data of the study had good reliability. At the same time, a previous study reported that B-mode ultrasound was a reliable and reproducible method for detecting the thickness of plantar fascia and SWE mode was a reliable and reproducible method for detecting the elasticity of plantar fascia (Wang et al., 2019).
Plantar fasciitis is one of the most common foot musculoskeletal diseases in primary diagnosis and treatment institutions (Thing et al., 2012; Young, 2012), and it is more likely to occur in patients with flatfoot than with normal foot (Riddle et al., 2003). It is characterized by heel pain after rest because it mainly affects the plantar fascia inserted into the calcaneus (Huang et al., 2000). The pain can also extend along the length of the plantar fascia (Thomas et al., 2016; Babatunde et al., 2019). In this study, the maximum Young’s modulus of proximal plantar fascia was 387.1kPa, while that of a normal foot was about 300kPa (Wang et al., 2019). Studies showed that there was a positive correlation between Young’s modulus and tendon force (Yeh et al., 2013; Yeh et al., 2016). Thus, the increased Young's modulus of the proximal plantar fascia indicates that the plantar fascia bears greater stress, leading more easily to the degeneration of plantar fascia (Huffer et al., 2017). The increase of Young’s modulus in plantar fascia near calcaneus attachment in patients with flatfoot may provide a theoretical explanation for the high incidence of plantar fasciitis in patients with flatfoot.
The results of the study showed that the plantar fascia of the flexible flatfoot was spatially dependent from proximal to distal, and the thickness and Young’s modulus of the five branches decreased gradually from proximal to distal. The differences between different parts were statistically significant. This feature of gradient changes is consistent with the results in normal plantar fascia (Wang et al., 2019). In the finite element model, the plantar fascia is often regarded as a linear elastic material, and the whole plantar fascia has the same Young’s modulus (Phan et al., 2021). Thus, the spatial distribution feature (different Young’s modulus in different regions) obtained in this study is helpful to define more accurate material properties for flatfeet finite element models to achieve more meaningful simulation results.
However, among the DB, Wang et al. (2019) showed that the thickness and Young’s modulus between the five branches of the normal plantar fascia were the greatest under the first MH, while this study showed that the thickness and Young’s modulus under the second and third MH were greater in patients with flexible flatfoot. At the same time, this study showed that the peak pressures under the second and third MH were greater than that under the fourth and fifth MH, which was consistent with the results of Buldt et al. (2018) and Hillstrom et al. (2013). It is speculated that this result may be due to the difference in the degree of collapse of the medial and lateral longitudinal arches in patients with flexible flatfoot. These results indicate that, in patients with flexible flatfoot, the degree of collapse of the medial longitudinal arch is more than that of the lateral arch, resulting in higher force and higher pressure on the medial side in the push-off phase. The stronger pressure stimulates plantar fascia, leading to its degeneration (Wearing et al., 2006). Shiotani et al. (2019) also noted that plantar fascia is mechanically stretched, so the morphology and mechanical properties of plantar fascia may be adapted to stress accumulation.
The center of pressure (COP) is defined as the centroid of the pressure distribution at a series of moments in time as the ground reaction is applied over the plantar surface of the foot (Cho and Choi, 2005). It was found that the peak pressures under the second and third metatarsals were higher than those under the other metatarsals. Thus, the COP would move laterally from the first MH. These results were the same as those of Han et al. (2011). They found that, in the normal foot, the trajectory of the COP moved from the lateral heel, moved medially in forefoot, and then ended at the big toe. In flatfeet, the COP moved straight from the heel to the toe without medial shifting in the forefoot. There was a tendency for the COP in flatfoot to shift laterally in the forefoot than the COP in normal foot. These results also confirm our inference; that is, the medial longitudinal arch collapses more than the lateral arch in flatfoot, which leads to the higher force and higher pressure under the second and third metatarsals and the COP moving outward.
Morphologic and mechanical properties of the plantar fascia may be important factors affecting the plantar pressure distribution because the primary task of the plantar fascia is to stabilize the foot arches (McKeon et al., 2015). Studies (Tas and Cetin, 2019a) also show that there is a significant positive correlation between plantar pressure distribution and the thickness of plantar fascia. Higher plantar pressure may lead to plantar fascia hypertrophy. Foot orthoses could modify tissue loading by altering kinematics, kinetics, muscle activity, and sensory feedback (Mills et al., 2010), and they have been demonstrated to have a good therapeutic effect in plantar fasciitis (Buchbinder, 2004). The changes in morphology and mechanical properties of plantar fascia and peak pressure of the forefoot in patients with flexible flatfoot found in our study may provide the basis for the development of new foot orthoses for flexible flatfoot.
There were limitations in this study. Firstly, the sample size was limited to 15 cases, and there was no grading according to mild, moderate, and severe flexible flatfeet. However, the results showed that although the sample size is small and there may be some differences in disease degree among participants, the spatial distribution characteristics of thickness and Young’s modulus of plantar fascia in all 15 flexible flat feet were similar, which indicated that the spatial distribution characteristics are less affected by the disease severity, and the research results may have broad representative significance. Secondly, this study did not include the control group, but our group has previously conducted and published one study on the morphology and mechanical properties of plantar fascia in normal feet (Wang et al., 2019). In addition, the age, height, and body weight of the volunteers who participated in this study are similar to those in our previous published work. Therefore, we cited and employed the published data (normal foot data) as the healthy control group in this study (Wang et al., 2019). Thirdly, though SWE has been used to evaluate the material properties of plantar fascia (Shiotani et al., 2019; Tas and Cetin, 2019b), studies have shown that the shear wave velocity of layered tissue is affected by its thickness and surrounding tissue properties (Helfenstein-Didier et al., 2016; Martin et al., 2018, 2019; Sadeghi and Cortes, 2020); especially when the thickness of the relevant tissue is equal to or less than the wavelength, SWE is no longer applicable (Li et al., 2018). The thickness of plantar fascia measured in this study is millimeter, which is far greater than the wavelength. In addition, the results of Helfenstein-Didier et al. (2016) in measuring the human Achilles tendon show that there is a high correlation between the shear modulus measured by SWE and the new guided wave technology-phase velocity mode, even considering the influence of thickness. Therefore, although the results of the differences between different positions in plantar fascia as well as between patients with flexible flat feet and healthy volunteers in this study may not be affected, it is necessary to explore the influence of thickness on the properties of plantar fascia materials by using guided wave technology in the future.
Conclusion
This study found a gradient distribution in that the thickness and Young’s modulus gradient decreased from the proximal end to the distal end of plantar fascia in the longitudinal arch of flexible flatfeet. In the transverse arch, the thickness and Young’s modulus under the second and third DB were larger than those under the other three DB in flexible flatfoot, and the peak pressures under the second and third MH were also larger than those under the other three MH in patients with flexible flatfoot. These findings deepen our understanding of the changes of biomechanical properties and may be meaningful for the study of pathological mechanisms and therapy for flexible flatfoot.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding authors.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of the Second Hospital of Jilin University (No. 2020085). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
ZQ and ZJ were responsible for the experiments and manuscript preparation. JW, FC, and LR participated in discussions and revisions. JL and LR worked as supervisors for all procedures.
Funding
This research was supported by the project of National Natural Science Foundation of China (No. 52175270), the Key Project of the National Natural Science Foundation of China (No. 91848204), National Natural Science Foundation of China (No. 52005209), the Interdisciplinary Research Funding Program for Doctoral Students of Jilin University (No. 101832020DJX049), and China Postdoctoral Science Foundation (No. 2021T140260 and No. 2021M691206).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors would like to thank all the volunteers who contributed to this study.
Abbreviations
SWE, shear wave elastography; PF, proximal fascia; MF, middle fascia; BF, branches of fascia; DB, distal branches; ROI, region of interest; 95% CI, 95% confidence interval; ICC1,1, intragroup correlation coefficient; one-way ANOVA, one-way analysis of variance; MH, metatarsal head; COP, center of pressure.
References
Babatunde, O. O., Legha, A., Littlewood, C., Chesterton, L. S., Thomas, M. J., Menz, H. B., et al. (2019). Comparative Effectiveness of Treatment Options for Plantar Heel Pain: a Systematic Review with Network Meta-Analysis. Br. J. Sports Med. 53, 182–194. doi:10.1136/bjsports-2017-098998
Buldt, A. K., Forghany, S., Landorf, K. B., Levinger, P., Murley, G. S., and Menz, H. B. (2018). Foot Posture Is Associated with Plantar Pressure during Gait: A Comparison of normal, Planus and Cavus Feet. Gait & Posture 62, 235–240. doi:10.1016/j.gaitpost.2018.03.005
Carr, J. B., Yang, S., and Lather, L. A. (2016). Pediatric Pes Planus: A State-Of-The-Art Review. Pediatrics 137, e20151230. doi:10.1542/peds.2015-1230
Cavanagh, P. R., Rodgers, M. M., and liboshi, A. (1987). Pressure Distribution under Symptom-free Feet during Barefoot Standing. Foot & Ankle 7, 262–278. doi:10.1177/107110078700700502
Chen, T. L.-W., Agresta, C. E., Lipps, D. B., Provenzano, S. G., Hafer, J. F., Wong, D. W.-C., et al. (2019a). Ultrasound Elastographic Assessment of Plantar Fascia in Runners Using Rearfoot Strike and Forefoot Strike. J. Biomech. 89, 65–71. doi:10.1016/j.jbiomech.2019.04.013
Chen, T. L.-W., Wong, D. W.-C., Wang, Y., Lin, J., and Zhang, M. (2019b). Foot Arch Deformation and Plantar Fascia Loading during Running with Rearfoot Strike and Forefoot Strike: A Dynamic Finite Element Analysis. J. Biomech. 83, 260–272. doi:10.1016/j.jbiomech.2018.12.007
Cho, W. H., and Choi, H. (2005). Center of Pressure (COP) during the Postural Balance Control of High-Heeled Woman. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2005, 2761–2764. doi:10.1109/IEMBS.2005.1617044
Guo, J., Liu, X., Ding, X., Wang, L., and Fan, Y. (2018). Biomechanical and Mechanical Behavior of the Plantar Fascia in Macro and Micro Structures. J. Biomech. 76, 160–166. doi:10.1016/j.jbiomech.2018.05.032
Haen, T. X., Roux, A., Soubeyrand, M., and Laporte, S. (2017). Shear Waves Elastography for Assessment of Human Achilles Tendon's Biomechanical Properties: an Experimental Study. J. Mech. Behav. Biomed. Mater. 69, 178–184. doi:10.1016/j.jmbbm.2017.01.007
Hamstra-Wright, K. L., Bliven, K. C. H., and Bay, C. (2015). Risk Factors for Medial Tibial Stress Syndrome in Physically Active Individuals Such as Runners and Military Personnel: a Systematic Review and Meta-Analysis. Br. J. Sports Med. 49, 362–369. doi:10.1136/bjsports-2014-093462
Han, J. T., Koo, H. M., Jung, J. M., Kim, Y. J., and Lee, J. H. (2011). Differences in Plantar Foot Pressure and COP between Flat and Normal Feet during Walking. J. Phys. Ther. Sci. 23, 683–685. doi:10.1589/jpts.23.683
Helfenstein-Didier, C., Andrade, R. J., Brum, J., Hug, F., Tanter, M., Nordez, A., et al. (2016). In Vivoquantification of the Shear Modulus of the Human Achilles Tendon during Passive Loading Using Shear Wave Dispersion Analysis. Phys. Med. Biol. 61, 2485–2496. doi:10.1088/0031-9155/61/6/2485
Hicks, J. H. (1954). The Mechanics of the Foot. II. The Plantar Aponeurosis and the Arch. J. Anat. 88, 25–30.
Hillstrom, H. J., Song, J., Kraszewski, A. P., Hafer, J. F., Mootanah, R., Dufour, A. B., et al. (2013). Foot Type Biomechanics Part 1: Structure and Function of the Asymptomatic Foot. Gait & Posture 37, 445–451. doi:10.1016/j.gaitpost.2012.09.007
Huang, H. H., Qureshi, A. A., and Biundo, J. J. (2000). Sports and Other Soft Tissue Injuries, Tendinitis, Bursitis, and Occupation-Related Syndromes. Curr. Opin. Rheumatol. 12, 150–154. doi:10.1097/00002281-200003000-00009
Huffer, D., Hing, W., Newton, R., and Clair, M. (2017). Strength Training for Plantar Fasciitis and the Intrinsic Foot Musculature: A Systematic Review. Phys. Ther. Sport 24, 44–52. doi:10.1016/j.ptsp.2016.08.008
Ker, R. F., Bennett, M. B., Bibby, S. R., Kester, R. C., and Alexander, R. M. (1987). The spring in the Arch of the Human Foot. Nature 325, 147–149. doi:10.1038/325147a0
Kohls-Gatzoulis, J., Angel, J. C., Singh, D., Haddad, F., Livingstone, J., and Berry, G. (2004). Tibialis Posterior Dysfunction: a Common and Treatable Cause of Adult Acquired Flatfoot. BMJ 329, 1328–1333. doi:10.1136/bmj.329.7478.1328
Martin, J. A., Brandon, S. C. E., Keuler, E. M., Hermus, J. R., Ehlers, A. C., Segalman, D. J., et al. (2018). Gauging Force by Tapping Tendons. Nat. Commun. 9, 1592. doi:10.1038/s41467-018-03797-6
Martin, J. A., Schmitz, D. G., Ehlers, A. C., Allen, M. S., and Thelen, D. G. (2019). Calibration of the Shear Wave Speed-Stress Relationship in Ex Vivo Tendons. J. Biomech. 90, 9–15. doi:10.1016/j.jbiomech.2019.04.015
Mckeon, P. O., Hertel, J., Bramble, D., and Davis, I. (2015). The Foot Core System: a New Paradigm for Understanding Intrinsic Foot Muscle Function. Br. J. Sports Med. 49, 290. doi:10.1136/bjsports-2013-092690
Mills, K., Blanch, P., Chapman, A. R., Mcpoil, T. G., and Vicenzino, B. (2010). Foot Orthoses and Gait: a Systematic Review and Meta-Analysis of Literature Pertaining to Potential Mechanisms. Br. J. Sports Med. 44, 1035–1046. doi:10.1136/bjsm.2009.066977
Nirenberg, M., Ansert, E., Campbell, J., and Curran, M. (2020). Forensic Implications of Foot Arch index Comparison between Dynamic Bare Footprints and Shoe Insole Foot Impressions. Sci. Justice 60, 375–380. doi:10.1016/j.scijus.2020.03.001
Pehlivan, O., Cilli, F., Mahirogullari, M., Karabudak, O., and Koksal, O. (2009). Radiographic Correlation of Symptomatic and Asymptomatic Flexible Flatfoot in Young Male Adults. Int. Orthopaedics (Sicot) 33, 447–450. doi:10.1007/s00264-007-0508-5
Phan, P. K., Vo, A. T. N., Bakhtiarydavijani, A., Burch, R., Smith, B., Ball, J. E., et al. (2021). In Silico Finite Element Analysis of the Foot Ankle Complex Biomechanics: A Literature Review. J. Biomech. Eng. 143. doi:10.1115/1.4050667
Rao, S., Song, J., Kraszewski, A., Backus, S., Ellis, S. J., Md, J. T. D., et al. (2011). The Effect of Foot Structure on 1st Metatarsophalangeal Joint Flexibility and Hallucal Loading. Gait & Posture 34, 131–137. doi:10.1016/j.gaitpost.2011.02.028
Riddle, D. L., Pulisic, M., Pidcoe, P., and Johnson, R. E. (2003). Risk Factors for Plantar Fasciitis. The J. Bone Jt. Surgery-American Volume 85, 872–877. doi:10.2106/00004623-200305000-00015
Sadeghi, S., and Cortes, D. H. (2020). Measurement of the Shear Modulus in Thin-Layered Tissues Using Numerical Simulations and Shear Wave Elastography. J. Mech. Behav. Biomed. Mater. 102, 103502. doi:10.1016/j.jmbbm.2019.103502
Shiotani, H., Yamashita, R., Mizokuchi, T., Naito, M., and Kawakami, Y. (2019). Site- and Sex-Differences in Morphological and Mechanical Properties of the Plantar Fascia: A Supersonic Shear Imaging Study. J. Biomech. 85, 198–203. doi:10.1016/j.jbiomech.2019.01.014
Stecco, C., Corradin, M., Macchi, V., Morra, A., Porzionato, A., Biz, C., et al. (2013). Plantar Fascia Anatomy and its Relationship with Achilles Tendon and Paratenon. J. Anat. 223, 665–676. doi:10.1111/joa.12111
Sung, P. S. (2016). The Ground Reaction Force Thresholds for Detecting Postural Stability in Participants with and without Flat Foot. J. Biomech. 49, 60–65. doi:10.1016/j.jbiomech.2015.11.004
Tas, S., and Cetin, A. (2019a). An Investigation of the Relationship between Plantar Pressure Distribution and the Morphologic and Mechanic Properties of the Intrinsic Foot Muscles and Plantar Fascia. Gait Posture 72, 217–221. doi:10.1016/j.jbiomech.2015.11.004
Taş, S., and Çetin, A. (2019b). Mechanical Properties and Morphologic Features of Intrinsic Foot Muscles and Plantar Fascia in Individuals with Hallux Valgus. Acta Orthopaedica et Traumatologica Turcica 53, 282–286. doi:10.1016/j.aott.2019.03.009
Thing, J., Maruthappu, M., and Rogers, J. (2012). Diagnosis and Management of Plantar Fasciitis in Primary Care. Br. J. Gen. Pract. 62, 443–444. doi:10.3399/bjgp12X653769
Thomas, M. J., Menz, H. B., and Mallen, C. D. (2016). Plantar Heel Pain. BMJ 353, i2175. doi:10.1136/bmj.i2175
Van Der Worp, H., Van Ark, M., Roerink, S., Pepping, G.-J., Van Den Akker-Scheek, I., and Zwerver, J. (2011). Risk Factors for Patellar Tendinopathy: a Systematic Review of the Literature. Br. J. Sports Med. 45, 446–452. doi:10.1136/bjsm.2011.084079
Wang, K., Liu, J., Wu, J., Qian, Z., Ren, L., and Ren, L. (2019). Noninvasive In Vivo Study of the Morphology and Mechanical Properties of Plantar Fascia Based on Ultrasound. Ieee Access 7, 53641–53649. doi:10.1109/Access.2019.2909409
Wang, Y.-T., Chen, J.-C., and Lin, Y.-S. (2020). Effects of Artificial Texture Insoles and Foot Arches on Improving Arch Collapse in Flat Feet. Sensors 20, 3667. doi:10.3390/s20133667
Wearing, S. C., Hills, A. P., Byrne, N. M., Hennig, E. M., and Mcdonald, M. (2004). The Arch index: a Measure of Flat or Fat Feet? Foot Ankle Int. 25, 575–581. doi:10.1177/107110070402500811
Wearing, S. C., Smeathers, J. E., Urry, S. R., Hennig, E. M., and Hills, A. P. (2006). The Pathomechanics of Plantar Fasciitis. Sports Med. 36, 585–611. doi:10.2165/00007256-200636070-00004
Welte, L., Kelly, L. A., Kessler, S. E., Lieberman, D. E., D'andrea, S. E., Lichtwark, G. A., et al. (2021). The Extensibility of the Plantar Fascia Influences the Windlass Mechanism during Human Running. Proc. R. Soc. B. 288, 20202095. doi:10.1098/rspb.2020.2095
Wong, C. K., Weil, R., and De Boer, E. (2012). Standardizing Foot-type Classification Using Arch index Values. Physiother. Can. 64, 280–283. doi:10.3138/ptc.2011-40
Wu, C.-H., Chang, K.-V., Mio, S., Chen, W.-S., and Wang, T.-G. (2011). Sonoelastography of the Plantar Fascia. Radiology 259, 502–507. doi:10.1148/radiol.11101665
Yeh, C.-L., Kuo, P.-L., Gennisson, J.-L., Brum, J., Tanter, M., and Li, P.-C. (2016). Shear Wave Measurements for Evaluation of Tendon Diseases. IEEE Trans. Ultrason. Ferroelect., Freq. Contr. 63, 1906–1921. doi:10.1109/TUFFC.2016.2591963
Yeh, C.-L., Kuo, P.-L., and Li, P.-C. (2013). Correlation between the Shear Wave Speed in Tendon and its Elasticity Properties. 2013 Ieee Int. Ultrason. Symp. (Ius), 9–12. doi:10.1109/Ultsym.2013.0003
Keywords: flexible flatfoot, plantar fascia, shear wave elastography, morphology properties, mechanical properties
Citation: Qian Z, Jiang Z, Wu J, Chang F, Liu J, Ren L and Ren L (2021) Morphology and Mechanical Properties of Plantar Fascia in Flexible Flatfoot: A Noninvasive In Vivo Study. Front. Bioeng. Biotechnol. 9:727940. doi: 10.3389/fbioe.2021.727940
Received: 20 June 2021; Accepted: 16 August 2021;
Published: 15 September 2021.
Edited by:
Fang Pu, Beihang University, ChinaCopyright © 2021 Qian, Jiang, Wu, Chang, Liu, Ren and Ren. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lei Ren, lei.ren@manchester.ac.uk; Jing Liu, jingliu@jlu.edu.cn