- 1Graduate School of Saitama Prefectural University, Graduate Course of Health and Social Services, Saitama, Japan
- 2Department of Biological Sciences, The University of Tokyo, Tokyo, Japan
- 3Department of Health Science, University of Human Arts and Sciences, Saitama, Japan
- 4Department of Health and Social Services, Saitama Prefectural University, Saitama, Japan
Kinetic multi-segment foot models have been proposed to evaluate the forces and moments generated in the foot during walking based on inverse dynamics calculations. However, these models did not consider the plantar aponeurosis (PA) despite its potential importance in generation of the ground reaction forces and storage and release of mechanical energy. This study aimed to develop a novel multi-segment foot model incorporating the PA to better elucidate foot kinetics. The foot model comprised three segments: the phalanx, forefoot, and hindfoot. The PA was modeled using five linear springs connecting the origins and the insertions via intermediate points. To demonstrate the efficacy of the foot model, an inverse dynamic analysis of human gait was performed and how the inclusion of the PA model altered the estimated joint moments was examined. Ten healthy men walked along a walkway with two force plates placed in series close together. The attempts in which the participant placed his fore- and hindfoot on the front and rear force plates, respectively, were selected for inverse dynamic analysis. The stiffness and the natural length of each PA spring remain largely uncertain. Therefore, a sensitivity analysis was conducted to evaluate how the estimated joint moments were altered by the changes in the two parameters within a range reported by previous studies. The present model incorporating the PA predicted that 13%–45% of plantarflexion in the metatarsophalangeal (MTP) joint and 8%–29% of plantarflexion in the midtarsal joints were generated by the PA at the time of push-off during walking. The midtarsal joint generated positive work, whereas the MTP joint generated negative work in the late stance phase. The positive and negative work done by the two joints decreased, indicating that the PA contributed towards transfer of the energy absorbed at the MTP joint to generate positive work at the midtarsal joint during walking. Although validation is limited due to the difficulty associated with direct measurement of the PA force in vivo, the proposed novel foot model may serve as a useful tool to clarify the function and mechanical effects of the PA and the foot during dynamic movements.
1 Introduction
The multi-segment foot model, such as the Oxford foot model (Carson et al., 2001; Wright et al., 2011; Li et al., 2022) and the Leardini foot model (Leardini et al., 2007; Deschamps et al., 2012; Watari et al., 2021), was developed to replace the conventional single-segment foot model for detailed in vivo evaluation of the foot segment kinematics during movements. Its clinical utility has been highly appreciated following assertions by several studies that foot kinematics are affected by age (Arnold et al., 2014; Deschamps et al., 2017), sex (Takabayashi et al., 2018; Sekiguchi et al., 2020), and deformities such as flat foot (Kim et al., 2020) or hallux valgus (Shin et al., 2019). However, detailed kinetics of the human foot during dynamic movements have not been similarly investigated, mainly due to the unavailability of a detailed foot model to estimate internal forces and moments within the foot segments. Kinetic assessments of the human foot during dynamic movements are particularly important since larger forces are applied to the foot and movements of the foot bones are greater in dynamic movements, possibly associated to the pathogenesis of foot injuries or disorders.
Recently, a few kinetic multi-segment foot models have been used to calculate the inter-foot segment moments based on an inverse dynamics analysis during gait (Bruening et al., 2012; Dixon et al., 2012; Saraswat et al., 2014; Kevin et al., 2017). However, these models lacked anatomical accuracy for two main reasons: 1) they did not calculate the detailed inertial properties of the divided foot, and 2) they did not consider forces generated by the plantar aponeurosis (PA). In previous studies, the mass and inertial tensor of each foot segment were determined arbitrarily (Dixon et al., 2012) or calculated by assuming a mathematical model such as a cylinder (Bruening et al., 2012; Saraswat et al., 2014; Kevin et al., 2017); however, such an assumption is unreasonable as the radius and height of each foot segment are not uniform. Although the influence of the inertia of the foot is expected to be relatively minor due to the small mass of the foot (Zelik and Honert, 2018), the acceleration in the most distal segment of the leg could be large, particularly in dynamic movements such as running and jumping. Therefore, accurate identification of the inertial parameters of each foot segment might be important for inverse dynamics analysis to calculate the joint moment of the foot.
The PA is an elastic band that supports the longitudinal foot arch, which consists of high density of collagen fibers. The occurrence of the windlass mechanism at push-off during walking and running has been proposed, wherein dorsiflexion of the metatarsophalangeal (MTP) joint winds the PA around the metatarsal heads, thereby increasing the rigidity of the foot (Hicks, 1954). Based on an inverse kinematic analysis, a study reported that the PA generates a tension force that is 1.5 times the bodyweight (Caravaggi et al., 2009), indicating that the forces generated by the PA are quite large, potentially having a major effect on the joint moments and forces computationally estimated based on an inverse dynamics analysis. Further, wire electromyographic analysis has shown that both the PA and the plantar intrinsic foot muscles contribute to increasing the foot rigidity in response to the magnitude of forces, such as stronger push-off during walking or running (Péter et al., 2015; Kelly et al., 2019; Farris et al., 2020). However, measuring the electromyography of the intrinsic muscles of the foot is difficult due to the invasive methods required. Dynamic finite element models of the human foot (e.g., Ito et al., 2022) can be used to estimate forces and moments generated by the foot muscles, but they are of limited utility due to large computational cost. Therefore, to evaluate foot function and elucidate the pathogenesis of foot injuries or disorders, it would be useful to establish a multi-segment kinetics model of the foot incorporating the PA that can be applied non-invasively during locomotion to estimate more accurate forces and moments in the foot.
In this study, we aimed to develop a novel multi-segment foot model to analyze foot kinematics and kinetics during dynamic movements by incorporating more accurate inertial parameters and the PA, and to demonstrate the efficacy of the developed foot model by applying it to gait analysis. Efforts have been previously made to incorporate the PA in a biomechanical foot model to estimate the mechanical contribution of the PA (Caravaggi et al., 2009; Chen et al., 2019; Farris et al., 2020; Welte et al., 2021). However, these studies estimated only the forces generated by the PA, but not the resultant joint moments within the foot segments generated during movements. For this purpose, a multi-segment foot model is necessary, but no studies have attempted to incorporate the PA in a multi-segment foot model for estimation of the resultant joint moments that can otherwise only be quantified using invasive techniques.
2 Materials and Methods
2.1 Model
To develop an anatomically accurate multi-segment foot model, computed tomography (CT) data of the foot of an adult male (age: 42 years, weight: 72 kg, height: 172 cm) were obtained (Ito et al., 2022). Three-dimensional surface models of the foot surface and skeleton were constructed (Figure 1A) using segmentation software (Analyze 9.0, Biomedical Imaging Resource, Mayo Clinic, Rochester, MN, United States). The human foot was represented as a chain comprising three segments (phalanx, forefoot, and hindfoot segments) based on the marker positions attached to anatomical landmarks that are widely used in kinematic studies (Leardini et al., 2007) (Table 1; Figure 1A). The joints between the tibia and hindfoot (ankle joint), hindfoot and forefoot (midtarsal joint), and forefoot and phalanx segments (MTP joint) were defined as the midpoints between the lateral and medial malleoli (ANKL and ANKM), the navicular tubercle (TN) and the fifth metatarsal base (VMB), and the heads of the first and fifth metatarsals (FMH and VMH), respectively. The coordinate systems of the three segments were defined as follows: the y and z axes of the hindfoot segment were defined as the normal vector of the plane defined by the superior (CA) and inferior (HE) points of the calcaneal tuberosity and the midtarsal joint, and the axis connecting the CA and HE, respectively; those of the forefoot segment were defined as the axis connecting the TN and VMB, and the normal vector of the plane defined by the TN, VMB, and MTP joints, respectively; those of the phalanx segment were defined as the axis connecting the FMH and VMH, and the normal vector of the plane defined by FMH, VMH, and the head of the first proximal phalanx (PM), respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/894731/fbioe-10-894731-HTML-r1/image_m/fbioe-10-894731-g001.jpg)
FIGURE 1. Multi-segment foot model incorporating the plantar aponeurosis (PA) for detailed kinematic and kinetic analyses of the foot. (A) Marker placements and the definitions of the coordinate systems of the phalanx, forefoot, and hindfoot segments. (B) The model of PA. The PA was modeled as five linear springs connecting the origin, via point, and insertion. (C) Segmentation of the phalanx, forefoot, and hindfoot segments. (D) Forces applied to the origin, via point, and insertions of the PA.
We represented the PA using five linear springs (PA1−5, from medial to lateral) connecting the hindfoot and phalanx segments via a point on the metatarsal plantar surface (Figure 1B). The positions of the origin, via point, and insertion of the PA were obtained from the CT data of the foot (Ito et al., 2022). The origin of the PA was defined as the plantar surface of the calcaneal tuberosity. The via points were defined as the plantar surfaces of the sesamoid of the first metatarsal and the second to fifth metatarsal heads. The insertion was defined as the plantar surface of the base of the first to fifth proximal phalanges. The via points were necessary to model the wrapping paths of PA slips around the metatarsal heads, analogous to the pully model in Caravaggi et al. (2009). To create a subject-specific model, the positions of the origin, via point, and insertion extracted from the CT data were described in the segment coordinate systems of the hindfoot, forefoot, and phalanx segments, the origins of which were HE, VMB, and VMH, respectively, and scaled by the length of the line connecting ANKL and VMH.
The length of the ith PA was calculated as the sum of the lengths from the origin to the via point and from the via point to insertion as:
The PA tension force (
where
where
Based on the free body diagrams in Figure 2, the Newton−Euler equations of the motion of the multi-segment foot model consisting of the phalanx, forefoot, and hindfoot segments can be written as:
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/894731/fbioe-10-894731-HTML-r1/image_m/fbioe-10-894731-g002.jpg)
FIGURE 2. Free body diagram of the phalanx (A), forefoot (B), and hindfoot segment (C). See Table 2 for notations.
The variables in the equations are presented in Table 2. The left superscripts P, F and H denote the coordinate systems in which the corresponding vector or matrix was represented.
2.2 Inertial Properties
To calculate the inertial properties (i.e., mass, inertial tensor around the center of mass (COM), and position of the center of mass in each segment), the foot surface was divided into segments by planes passing through ANKL and ANKM, TN and VMB, and FMH and VMH (Figure 1C). Therefore, the hindfoot, forefoot, and phalanx segments correspond to the calcaneus, talus, cuboid, and navicular; cuneiforms and metatarsals; and five phalanges, respectively. The inertial parameters of each segment were calculated using a computer-aided design software (Autodesk Inventor Professional, Autodesk, United States), assuming a homogeneous segment composition and a density of 1.1 g/cm3 (Winter, 1990). To create an individual-specific foot model, the relative mass of the foot segments and the relative inertial tensor around the segment COM were computed based on the inertial parameters from the CT data. Relative foot segment mass was defined as the mass of the foot segment as a percentage of the total foot mass. The relative inertial tensor around the segment COM was defined as the inertial tensor around the segment COM normalized by 5/3 power of each segment mass. Relative COM position was defined as the location of the COM, expressed as a percentage of the segment length from the proximal point of the segment, assuming that the COM is located on the line connecting the centers of the proximal and distal joints.
2.3 Participants
Ten males (age: 23.9 ± 3.0 years, height: 171.8 ± 5.1 cm, weight: 62.8 ±8.2 kg) without any deformity of the foot and lower extremity and with no history of orthopedic, neurological, and musculoskeletal disorders that are likely to affect gait were recruited for data collection. The number of participants was determined by referring to previous studies (Leardini et al., 2007; Portinaro et al., 2014). All participants provided written informed consent following a detailed explanation of the study’s purpose and the risks involved. The experimental procedures used in this study complied with the Declaration of Helsinki and were approved by the Ethics Committee on Human Experimentation at Saitama Prefectural University (No. 29508).
2.4 Experimental Procedure
Infrared-reflecting markers (diameters: 9.5 and 14 mm) were attached to 65 landmarks on the foot and the whole body, according to Table 1 and Figure 1A, and the Plug-in-Gait Full-body Ai model (Davis et al., 1991; Vicon Motion Systems Limited, 2016) (Supplementary Figure S1; Supplementary Table S1), respectively. The participants walked at a self-selected speed along a walkway with four force plates (two per side) placed in series close together. The three attempts in which the fore- and hindfoot contacted the front and rear force plates, respectively (Supplementary Figure S2), were selected for the inverse dynamic analysis.
Marker trajectories were collected using the Vicon Nexus 2.10.2, a three-dimensional motion analysis system (Vicon, Oxford, United Kingdom) with 20 infrared cameras at 100 Hz. Ground reaction force was collected from four force plates (Kistler Instrumente AG, Winterthur, Switzerland) at 1,000 Hz. All data were synchronized using the Vicon Workstation v4.5 software and saved for offline analysis.
2.5 Data Processing and Analysis
Marker trajectories and ground reaction forces were filtered using a zero-phased lag and fourth-order Butterworth filter with a cutoff frequency of 10 Hz. The 3D rotation angles of the ankle joint were described by the y-x-z Euler angle, and those of the midtarsal and MTP joints were described by z-x-y Euler angles. The rotational angles around the y-, x-, z-axes represent plantarflexion–dorsiflexion, inversion–eversion and adduction−abduction, respectively. The MTP joint in the present model is regarded as a hinge joint allowing only plantarflexion–dorsiflexion; hence the rotational angles around the other two axes are not presented.
For inverse dynamics, the total mass of each individual’s foot was estimated to be 0.0145 times the body weight of the individual (Winter, 1990). The mass, inertial tensor, and COM position of each foot segment were then calculated based on the relative mass, inertial tensor, and COM position obtained from the CT scan of the foot used to create the present multi-segment foot model (See Section 2.1). The ground reaction forces were applied to the foot segments based on the force plate data. The positional relationship between the center of pressure and the markers specified the segments on which the two force vectors were applied. The tension forces due to the PAs were applied to the insertions, via points, and origins of the PAs (Figure 1D). The joint moments were calculated by solving the equations for motions (Eqs 7–12) consecutively from the distal phalanx segment to the proximal hindfoot segment. The calculated joint moment vectors were transformed to the corresponding proximal segment coordinate systems to match them with the joint angles. The joint powers were calculated by multiplying the joint angular velocities and joint moments. The joint power curves were integrated into the positive and negative portions during one gait cycle, which represented energy generation and energy absorption, respectively. All data were analyzed using MATLAB 2018a (MathWorks, Natick, MA, United States).
In this study, the inverse dynamic analysis was performed with and without the PA to investigate how inclusion of the PA in the model alters the estimated joint moments and work. For this, the maximum and minimum values of the calculated joint moments during one gait cycle were compared between the two conditions. In addition, the PA contribution rate (%PA contribution) was calculated as the percentage of the joint moment generated by the PA with respect to the net joint moment, to quantify the contribution of the PA on joint moment generation at the time of push-off during walking. The changes in the positive and negative work generated or absorbed during walking were also quantified and compared.
2.6 Sensitivity Analysis
We estimated the relative mass of the segment and the inertia tensor of each segment from the CT data of a single male participant. Further, due to the uncertainty regarding the stiffness and the natural length of each PA spring, a sensitivity analysis was conducted to evaluate how the estimated joint moments were altered by the changes in these parameters within a reasonable range. For the mass and inertia of the segment, we made 0.5- and 1.5-fold changes in the mass and inertial parameters of the foot model to assess how changes in the mass and inertial parameters affected the joint moments, which was calculated using inverse dynamics. We made the 0.5- and 1.5-fold changes in the mass and inertial parameters because the standard deviations of the mass and inertial parameters were reportedly about a quarter of the respective mean values (Zatsiorsky, 2002) and approximately 95% of the population lies within two standard deviations. For the stiffness, we decreased and increased the natural PA length (default
2.7 Statistical Analysis
To test for possible statistical differences in the joint moment and power profiles between the models with and without PA in a continuous manner, we performed a one-dimensional Statistical Parametric Mapping (SPM) paired t-test (Pataky, 2012). To compare the statistical differences in the estimated peak joint moments, joint works, and %PA contributions, we performed a paired t-test if the normality was confirmed using the Kolmogorov–Smirnov test. If the normality was violated, we used a signed-rank test for statistical comparisons. All statistical tests were performed with a significance level of 5% using MATLAB 2018a (MathWorks, Natick, MA, United States).
3 Results
3.1 Inertial Properties of the Foot
The relative mass of the segment, relative COM position, and relative inertial tensor around the COM of the phalanx, forefoot, and hindfoot segments are shown in Table 3. The size-normalized position vectors of the insertions, via points, and origins of the PAs are presented in Table 4. The mean and standard deviations of the masses of the phalanx, forefoot, and hindfoot segments of the 10 participants were calculated to be 0.131 ± 0.017 kg, 0.386 ± 0.051 kg, and 0.393 ± 0.052 kg, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/894731/fbioe-10-894731-HTML-r1/image_m/fbioe-10-894731-t003.jpg)
TABLE 3. Relative segment mass, relative center of mass position, and relative inertial tensor around the center of mass of each phalanx, forefoot, and hindfoot segment.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/894731/fbioe-10-894731-HTML-r1/image_m/fbioe-10-894731-t004.jpg)
TABLE 4. Size-normalized position vectors of the insertions, via points, and origins of the plantar aponeuroses.
3.2 Foot Kinematics
The mean and standard deviations of the speed, cycle duration, and stance phase duration of the measured walking were 1.33 ± 0.17 m/s, 1.03 ± 0.05 s and 0.63 ± 0.05 s, respectively. The 3D joint angle profiles of the MTP, midtarsal, and ankle joints during walking are presented in Figure 3. The joint angles during quiet standing are also presented as a reference point. The joints angles were positive for eversion, dorsiflexion, and abduction. The ankle joint plantarflexed after heel-contact and dorsiflexed during single support phase so that the shank vaulted over the foot. The joint then plantarflexed while approaching toe-off and returned to its natural angle during the swing phase. Following heel contact, the midtarsal and MTP joints dorsiflexed and plantarflexed, respectively. During the single support phase, both joints remained largely unchanged; however, in the late stance phase, the midtarsal and MTP joints plantarflexed and dorsiflexed, respectively. Eversion−inversion and abduction−adduction of the two joints were relatively small, though slight eversion of the ankle and midtarsal joints and abduction of the ankle joint were observed during the stance phase.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/894731/fbioe-10-894731-HTML-r1/image_m/fbioe-10-894731-g003.jpg)
FIGURE 3. Mean joint angle profiles of the foot during walking. Mean (solid line) ± standard deviation (red band). Means and standard deviations of the joint angles during quiet standing were also plotted on the right side of each graph.
Figure 4 displays the changes in the PA length and force during walking. The length profiles were generally consistent with each other for the four lateral PAs. They were stretched after the heel-contact and remained stretched until 50% of the gait cycle. During push-off, they were sharply shortened and remained shortened during the swing phase. PA1 was more largely stretched at 50% of the gait cycle than the other four PAs, mainly due to the elongation of the distal portion of the PA connecting the via point and insertion. The maximum tension forces exerted by the five PAs were 0.23 ± 0.06 N/kg, 0.15 ± 0.04 N/kg, 0.16 ± 0.04 N/kg, 0.17 ± 0.04 N/kg, 0.16 ± 0.04 N/kg, respectively. The resultant maximum net PA force was 0.87 ± 0.20 N/kg.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/894731/fbioe-10-894731-HTML-r1/image_m/fbioe-10-894731-g004.jpg)
FIGURE 4. Mean plantar aponeurosis length and force profiles during walking. Mean (solid line) ± standard deviation (red band).
3.3 Foot Kinetics
The joint moment and joint power profiles calculated with and without incorporating the PA are presented in Figure 5. Plantarflexion moments were generated by all three joints during walking, particularly in the late stance phase. The magnitudes of the peak plantarflexion moments were 0.18 ± 0.04 Nm/kg, 0.98 ± 0.18 Nm/kg, and 1.58 ± 0.18 Nm/kg for the MTP, midtarsal, and ankle joints, respectively. Slight inversion and abduction moments were generated in phase with the plantarflexion moment by the midtarsal joint. The power generated by the midtarsal joint was calculated to be positive, as the joint plantarflexed while generating the plantarflexion moment. In contrast, the power generated by the MTP joint was negative because the joint dorsiflexed while generating the plantarflexion moment. Therefore, positive and negative work was performed by the midtarsal and MTP joints, respectively, during walking. The peak positive power generated by the midtarsal joint and peak negative power absorbed by the MTP joint were estimated as 1.81 ± 0.73 W/kg and −1.01 ± 0.48 W/kg, respectively, during walking. The power in the eversion−inversion and abduction−adduction directions was almost zero.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/894731/fbioe-10-894731-HTML-r1/image_m/fbioe-10-894731-g005.jpg)
FIGURE 5. Mean joint moment (A) and joint power (B) profiles during walking. Mean (solid red line = with PA, dotted blue line = without PA) ± standard deviation (red and blue bands, respectively). Color bar below each graph shows the results of the SPM analysis.
Due to the presence of the PA model, the plantarflexion moments and power of the MTP and midtarsal joints were significantly smaller in magnitude in the model with incorporating the PA than that without incorporating the PA (p < 0.01). The abduction and inversion moments of the midtarsal joint during the stance phase were also significantly smaller in the model with incorporating the PA than that without incorporating the PA (p < 0.01). The mean peak plantar flexion moments of the MTP were 0.07 Nm/kg smaller (p = 0.01), and the mean peak inversion, plantar flexion, and abduction moments of the midtarsal joint were 0.02 Nm/kg, 0.2 Nm/kg, and 0.04 Nm/kg smaller, respectively (p < 0.01, p < 0.01, and p = 0.01, respectively) when incorporating the PA model (Figure 6A). The mean negative work done by the MTP joint and the positive work done by the midtarsal joint were significantly reduced by 3.0 and 2.2 J, respectively (p < 0.01 and p < 0.01, respectively) (Figure 6B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/894731/fbioe-10-894731-HTML-r1/image_m/fbioe-10-894731-g006.jpg)
FIGURE 6. Comparison of the joint moment (A) and joint power (B) profiles between models with (blue) and without the plantar aponeurosis (red).
3.4 Sensitivity Analysis
Figure 7 compares the mean joint moment profiles of the three joints when the inertial parameters were altered. The change in the inertial parameters had virtually no effect on the calculated joint moments. Figure 8A compares the mean peak moments at toe-off of the midtarsal and MTP joints when the PA stiffness parameters were altered. The change in the PA stiffness altered the peak plantarflexion moments of the two joints but not the moments in the other directions. It was estimated that 13%–45% of the plantarflexion moment of the MTP joint and 8%–29% of plantarflexion in the midtarsal joints were generated by the PA at the time of push-off during walking (Figure 8B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/894731/fbioe-10-894731-HTML-r1/image_m/fbioe-10-894731-g008.jpg)
FIGURE 8. Comparisons of the mean peak moments (A) and %plantar aponeurosis (PA) contribution (B) at toe-off when the PA stiffness parameters were altered.
4 Discussion
In this study, we developed a novel multi-segment foot model incorporating the PA for a detailed inverse dynamic analysis of the foot segments. We then ran simulations with the model to address the biomechanical consequences of incorporating the PA on the estimated joint moments and power during human walking. We observed that ∼13%–45% of the plantarflexion moment in the MTP joint and 8%–29% of the plantarflexion moment in the midtarsal joints were generated by the PA at the time of push-off during walking (Figure 8). If the PA was not incorporated in the moment estimation, the joint moments generated by the foot muscles would be overestimated by these amounts. Therefore, we were able to perform a more precise estimation of the joint moments generated within the foot segments, which are not measurable non-invasively, during movements such as walking, running, and jumping. As such, our multi-segment foot model incorporating the PA contributes toward elucidating the basic biomechanics and motor control of the human foot, and to clarify the pathogenic mechanism and possible surgical or rehabilitative interventions for the treatment and prevention of foot pathologies.
We predicted that the PA was stretched in the early and mid-stance phase but shortened in the late stance phase (Figure 4). This is consistent with previous studies reporting the PA elongation profile during walking (Caravaggi et al., 2009, 2010; Welte et al., 2021). However, this study also demonstrated that the midtarsal joint generated positive mechanical work and that the MTP joint generated negative mechanical work in the second half of the stance phase during walking (Figure 5), and if the PA was incorporated in the inverse dynamic calculation, the positive and negative works done by the two joints were both reduced (Figure 6). This might be because the PA springs contributed toward transfer of the energy absorbed at the MTP joint to generate positive work at the midtarsal joint during walking. The human foot possesses a longitudinal arch with the PA spanning its plantar side, allowing stretch and recoil of the PA springs to store and release mechanical energy for generation of efficient locomotion (Ker et al., 1987; Kim and Voloshin, 1995; Stearne et al., 2016). Our study, which was based on inverse dynamics using the proposed multi-segment foot model, clarified the detailed energy recovery mechanism embedded in the human foot that possibly contributes to the reduced energy cost in human bipedal walking.
In the present study, the spring constant of the PA before the split was calculated as approximately 80 N/mm by substituting the mean natural length (171 mm) and mean body weight (63 kg) into Eq. 3, derived based on the published reports that the strain of the PA was approximately 0.07 and the force generated by the PA was approximately 1.5 times the body weight at the time of push-off during walking (see Materials and Methods). The stiffness value is less than a half of the PA stiffness previously obtained in vitro (204 N/mm; Kitaoka et al., 1994) and in vivo (170 N/mm; Gefen, 2003). However, if the PA is elongated by maximum 6 % and 12% as reported in Caravaggi et al. (2009) and Gefen (2003) during walking, the estimated tensile force generated by the PA will be 2093 and 4186 N, respectively, using 204 N/mm, and 1744 and 3488 N, respectively, using 170 N/mm, presumably too large for the estimated tensile force generated by the PA during walking. Therefore, we believe that the presently estimated stiffness of the PA is of reasonable accuracy, but this must be confirmed by further investigations.
Our study provided inertial parameters of the phalanx, forefoot, and hind foot segments based on the CT scan data of the foot, and we observed that the errors in the inertial parameters had virtually no impact on the joint moment calculation because the inertial forces and moments were much smaller than the others. Thus, the present dataset should serve as a useful reference for inertial parameters of the kinetic multi-segment foot model.
One limitation of the present study was that the kinematic measurements were possibly affected by skin marker artifacts (Nester et al., 2007; Shultz et al., 2011; Schallig et al., 2021). However, skin motion at the foot is generally regarded as relatively small compared to those of other parts of the body (Leardini et al., 2005). Therefore, this limitation should not have a major effect on the current results. Another limitation of our study was that we could not quantitatively validate the estimated PA lengths and forces against the corresponding measured data since in vivo direct measurement of the strain and tension generated by the PA during walking is technically impossible. Due to this limitation, we tried to estimate the PA parameters as precisely as possible based on available information and then conducted a sensitivity analysis to show how the uncertainty in the parameters of the PA contributes to the estimated PA forces and joint moments. For more precise estimation of the PA forces and joint moments, efforts should be made to better identify the parameters necessary to quantify the stiffness of the PA based on dissections of cadaver specimens (Guo et al., 2018; Sichting et al., 2020) or medical imaging of living persons using magnetic resonance imaging (Ehrmann et al., 2014; Shiotani et al., 2021), ultrasound imaging (Crofts et al., 2014; Boussouar et al., 2017), and shear wave elastography (Chino et al., 2019; Wang et al., 2019). A third limitation of the study is that the PA model was assumed to be identical for all participants, despite potentially large individual differences in the properties of the PA due to differences in age (Cheng et al., 2014), sex (Pascual Huerta and Alarcón García, 2007; Shiotani et al., 2019), and deformity (Angin et al., 2014, 2018). Such individual differences should also be incorporated in the estimation of the PA forces and joint moments in future studies. A fourth limitation of the study is that the joints of the present foot model were actuated by joint moments instead of muscle forces. Therefore, the model only allows the estimation of the joint moments representing the net effect of the muscles around each joint, but not the estimation of the force generated by each muscle during movements. For this, modeling of the paths of muscles in the foot is necessary, and this should also be investigated in future studies.
5 Conclusion
The present study proposed a novel multi-segment foot model incorporating the PA to analyze foot kinetics during dynamic movements and demonstrated the efficacy of the developed foot model by applying it to gait analysis. The present model incorporating the PA predicted that 13%–45% of plantarflexion in the MTP joint and 8%–29% of plantarflexion in the midtarsal joints were generated by the PA at the time of push-off during walking. The present model also demonstrated that the midtarsal joint generated positive work and that the MTP joint generated negative work in the late stance phase. The positive and negative work done by the two joints were both reduced, indicating that the PA contributed towards transfer of the energy absorbed at the MTP joint to generate positive work at the midtarsal joint during walking. The proposed novel foot model may serve as a useful tool to clarify the function and mechanical effects of the PA and the foot during dynamic movements.
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 author.
Ethics Statement
The studies involving human participants were reviewed and approved by Ethics Committee on Human Experimentation at Saitama Prefectural University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
YM, NO, TK, and NK contributed to conception and design of the study; YM and NO constructed the model and the computational framework; YM collected the data and performed the data analysis; NO and HH helped carrying out the data analysis; YM wrote the first draft of the manuscript. NO, HH, TK, and NK edited the manuscript. All authors contributed to revisions, and approved the final submitted version of the manuscript.
Funding
This study was supported by the Grant-in-Aid for JSPS Research Fellowship for Young Scientist (20J14897) to YM.
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 acknowledge the members of the Fundamental Science Lab, Saitama Prefectural University, for their assistance in the present study and experiments.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2022.894731/full#supplementary-material
References
Angin, S., Crofts, G., Mickle, K. J., and Nester, C. J. (2014). Ultrasound Evaluation of Foot Muscles and Plantar Fascia in Pes Planus. Gait Posture 40, 48–52. doi:10.1016/J.GAITPOST.2014.02.008
Angin, S., Mickle, K. J., and Nester, C. J. (2018). Contributions of Foot Muscles and Plantar Fascia Morphology to Foot Posture. Gait Posture 61, 238–242. doi:10.1016/J.GAITPOST.2018.01.022
Arnold, J. B., Mackintosh, S., Jones, S., and Thewlis, D. (2014). Differences in Foot Kinematics between Young and Older Adults during Walking. Gait Posture 39, 689–694. doi:10.1016/j.gaitpost.2013.09.021
Boussouar, A., Meziane, F., and Crofts, G. (2017). Plantar Fascia Segmentation and Thickness Estimation in Ultrasound Images. Comput. Med. Imaging Graph. 56, 60–73. doi:10.1016/J.COMPMEDIMAG.2017.02.001
Bruening, D. A., Cooney, K. M., and Buczek, F. L. (2012). Analysis of a Kinetic Multi-Segment Foot Model Part II: Kinetics and Clinical Implications. Gait Posture 35, 535–540. doi:10.1016/j.gaitpost.2011.11.012
Caravaggi, P., Pataky, T., Goulermas, J. Y., Savage, R., and Crompton, R. (2009). A Dynamic Model of the Windlass Mechanism of the Foot: Evidence for Early Stance Phase Preloading of the Plantar Aponeurosis. J. Exp. Biol. 212, 2491–2499. doi:10.1242/jeb.025767
Caravaggi, P., Pataky, T., Günther, M., Savage, R., and Crompton, R. (2010). Dynamics of Longitudinal Arch Support in Relation to Walking Speed: Contribution of the Plantar Aponeurosis. J. Anat. 217, 254–261. doi:10.1111/j.1469-7580.2010.01261.x
Carson, M. C., Harrington, M. E., Thompson, N., O’Connor, J. J., and Theologis, T. N. (2001). Kinematic Analysis of a Multi-Segment Foot Model for Research and Clinical Applications: A Repeatability Analysis. J. Biomechanics 34, 1299–1307. doi:10.1016/S0021-9290(01)00101-4
Chen, T. L.-W., Wong, D. W.-C., Wang, Y., Lin, J., and Zhang, M. (2019). Foot Arch Deformation and Plantar Fascia Loading during Running with Rearfoot Strike and Forefoot Strike: A Dynamic Finite Element Analysis. J. Biomechanics 83, 260–272. doi:10.1016/j.jbiomech.2018.12.007
Cheng, J.-W., Tsai, W.-C., and Yu, T.-Y. (2014). Gender-related Effect of Aging on the Sonographic Appearance of Plantar Fascia. J. Musculoskelet. Pain 22, 33–37. doi:10.3109/10582452.2014.883035
Chino, K., Lacourpaille, L., Sasahara, J., Suzuki, Y., and Hug, F. (2019). Effect of Toe Dorsiflexion on the Regional Distribution of Plantar Fascia Shear Wave Velocity. Clin. Biomech. 61, 11–15. doi:10.1016/J.CLINBIOMECH.2018.11.003
Crofts, G., Angin, S., Mickle, K. J., Hill, S., and Nester, C. J. (2014). Reliability of Ultrasound for Measurement of Selected Foot Structures. Gait Posture 39, 35–39. doi:10.1016/J.GAITPOST.2013.05.022
Davis, R. B., Õunpuu, S., Tyburski, D., and Gage, J. R. (1991). A Gait Analysis Data Collection and Reduction Technique. Hum. Mov. Sci. 10, 575–587. doi:10.1016/0167-9457(91)90046-Z
Kevin, D., Eerdekens, M., Desmet, D., Matricali, G. A., Wuite, S., and Staes, F. (2017). Estimation of Foot Joint Kinetics in Three and Four Segment Foot Models Using an Existing Proportionality Scheme: Application in Paediatric Barefoot Walking. J. Biomechanics 61, 168–175. doi:10.1016/j.jbiomech.2017.07.017
Deschamps, K., Staes, F., Bruyninckx, H., Busschots, E., Jaspers, E., Atre, A., et al. (2012). Repeatability in the Assessment of Multi-Segment Foot Kinematics. Gait Posture 35, 255–260. doi:10.1016/j.gaitpost.2011.09.016
Deschamps, K., Staes, F., Peerlinck, K., Van Geet, C., Hermans, C., Matricali, G. A., et al. (2017). 3D Multi-Segment Foot Kinematics in Children: A Developmental Study in Typically Developing Boys. Gait Posture 52, 40–44. doi:10.1016/j.gaitpost.2016.11.022
Dixon, P. C., Böhm, H., and Döderlein, L. (2012). Ankle and Midfoot Kinetics during Normal Gait: A Multi-Segment Approach. J. Biomechanics 45, 1011–1016. doi:10.1016/j.jbiomech.2012.01.001
Ehrmann, C., Maier, M., Mengiardi, B., Pfirrmann, C. W. A., and Sutter, R. (2014). Calcaneal Attachment of the Plantar Fascia: MR Findings in Asymptomatic Volunteers. Radiology 272, 807–814. doi:10.1148/radiol.14131410
Farris, D. J., Birch, J., and Kelly, L. (2020). Foot Stiffening during the Push-Off Phase of Human Walking Is Linked to Active Muscle Contraction, and Not the Windlass Mechanism. J. R. Soc. Interface. 17, 20200208. doi:10.1098/rsif.2020.0208
Gefen, A. (2003). The In Vivo Elastic Properties of the Plantar Fascia during the Contact Phase of Walking. Foot Ankle Int. 24, 238–244. doi:10.1177/107110070302400307
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. Biomechanics 76, 160–166. doi:10.1016/J.JBIOMECH.2018.05.032
Hicks, J. H. (1954). The Mechanics of the Foot. II. The Plantar Aponeurosis and the Arch. J. Anat. 88, 25–30. Available at: http://www.ncbi.nlm.nih.gov/pubmed/13129168.
Ito, K., Nakamura, T., Suzuki, R., Negishi, T., Oishi, M., Nagura, T., et al. (2022). Comparative Functional Morphology of Human and Chimpanzee Feet Based on Three-Dimensional Finite Element Analysis. Front. Bioeng. Biotechnol. 9, 1334. doi:10.3389/FBIOE.2021.760486
Kelly, L. A., Farris, D. J., Cresswell, A. G., and Lichtwark, G. A. (2019). Intrinsic Foot Muscles Contribute to Elastic Energy Storage and Return in the Human Foot. J. Appl. Physiology 126, 231–238. doi:10.1152/japplphysiol.00736.2018
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
Kim, E. J., Shin, H. S., Takatori, N., Yoo, H. J., Cho, Y. J., Yoo, W. J., et al. (2020). Inter‐segmental Foot Kinematics during Gait in Elderly Females According to the Severity of Hallux Valgus. J. Orthop. Res. 38, 2409–2418. doi:10.1002/jor.24657
Kim, W., and Voloshin, A. S. (1995). Role of Plantar Fascia in the Load Bearing Capacity of the Human Foot. J. Biomechanics 28, 1025–1033. doi:10.1016/0021-9290(94)00163-X
Kitaoka, H. B., Luo, Z. P., Growney, E. S., Berglund, L. J., and An, K.-N. (1994). Material Properties of the Plantar Aponeurosis. Foot Ankle Int. 15, 557–560. doi:10.1177/107110079401501007
Leardini, A., Benedetti, M. G., Berti, L., Bettinelli, D., Nativo, R., and Giannini, S. (2007). Rear-foot, Mid-foot and Fore-Foot Motion during the Stance Phase of Gait. Gait Posture 25, 453–462. doi:10.1016/j.gaitpost.2006.05.017
Leardini, A., Chiari, L., Croce, U. D., and Cappozzo, A. (2005). Human Movement Analysis Using Stereophotogrammetry. Gait Posture 21, 212–225. doi:10.1016/j.gaitpost.2004.05.002
Li, J., Song, Y., Xuan, R., Sun, D., Teo, E.-C., Bíró, I., et al. (2022). Effect of Long-Distance Running on Inter-segment Foot Kinematics and Ground Reaction Forces: A Preliminary Study. Front. Bioeng. Biotechnol. 10, 1–9. doi:10.3389/fbioe.2022.833774
Nester, C., Jones, R. K., Liu, A., Howard, D., Lundberg, A., Arndt, A., et al. (2007). Foot Kinematics during Walking Measured Using Bone and Surface Mounted Markers. J. Biomechanics 40, 3412–3423. doi:10.1016/j.jbiomech.2007.05.019
Pascual Huerta, J., and Alarcón García, J. M. (2007). Effect of Gender, Age and Anthropometric Variables on Plantar Fascia Thickness at Different Locations in Asymptomatic Subjects. Eur. J. Radiology 62, 449–453. doi:10.1016/J.EJRAD.2007.01.002
Pataky, T. C. (2012). One-dimensional Statistical Parametric Mapping in Python. Comput. Methods Biomechanics Biomed. Eng. 15, 295–301. doi:10.1080/10255842.2010.527837
Péter, A., Hegyi, A., Stenroth, L., Finni, T., and Cronin, N. J. (2015). EMG and Force Production of the Flexor Hallucis Longus Muscle in Isometric Plantarflexion and the Push-Off Phase of Walking. J. Biomechanics 48, 3413–3419. doi:10.1016/j.jbiomech.2015.05.033
Portinaro, N., Leardini, A., Panou, A., Monzani, V., and Caravaggi, P. (2014). Modifying the Rizzoli Foot Model to Improve the Diagnosis of Pes-Planus: Application to Kinematics of Feet in Teenagers. J. Foot Ankle Res. 7, 1–7. doi:10.1186/s13047-014-0057-2
Saraswat, P., MacWilliams, B. A., Davis, R. B., and D’Astous, J. L. (2014). Kinematics and Kinetics of Normal and Planovalgus Feet during Walking. Gait Posture 39, 339–345. doi:10.1016/j.gaitpost.2013.08.003
Schallig, W., Streekstra, G. J., Hulshof, C. M., Kleipool, R. P., Dobbe, J. G. G., Maas, M., et al. (2021). The Influence of Soft Tissue Artifacts on Multi-Segment Foot Kinematics. J. Biomechanics 120, 110359. doi:10.1016/j.jbiomech.2021.110359
Sekiguchi, Y., Kokubun, T., Hanawa, H., Shono, H., Tsuruta, A., and Kanemura, N. (2020). Foot Kinematics of Impact Absorption and Force Exertion during Depth-Jump Using a Multi-Segment Foot Model. J. Med. Biol. Eng. 40, 757–765. doi:10.1007/s40846-020-00560-5
Shin, H. S., Lee, J. H., Kim, E. J., Kyung, M. G., Yoo, H. J., and Lee, D. Y. (2019). Flatfoot Deformity Affected the Kinematics of the Foot and Ankle in Proportion to the Severity of Deformity. Gait Posture 72, 123–128. doi:10.1016/j.gaitpost.2019.06.002
Shiotani, H., Maruyama, N., Kurumisawa, K., Yamagishi, T., and Kawakami, Y. (2021). Human Plantar Fascial Dimensions and Shear Wave Velocity Change In Vivo as a Function of Ankle and Metatarsophalangeal Joint Positions. J. Appl. Physiology 130, 390–399. doi:10.1152/japplphysiol.00485.2020
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. Biomechanics 85, 198–203. doi:10.1016/j.jbiomech.2019.01.014
Shultz, R., Kedgley, A. E., and Jenkyn, T. R. (2011). Quantifying Skin Motion Artifact Error of the Hindfoot and Forefoot Marker Clusters with the Optical Tracking of a Multi-Segment Foot Model Using Single-Plane Fluoroscopy. Gait Posture 34, 44–48. doi:10.1016/j.gaitpost.2011.03.008
Sichting, F., Holowka, N. B., Ebrecht, F., and Lieberman, D. E. (2020). Evolutionary Anatomy of the Plantar Aponeurosis in Primates, Including Humans. J. Anat. 237, 85–104. doi:10.1111/JOA.13173
Stearne, S. M., McDonald, K. A., Alderson, J. A., North, I., Oxnard, C. E., and Rubenson, J. (2016). The Foot's Arch and the Energetics of Human Locomotion. Sci. Rep. 6, 1–10. doi:10.1038/srep19403
Takabayashi, T., Edama, M., Inai, T., and Kubo, M. (2018). Sex-related Differences in Coordination and Variability Among Foot Joints during Running. J. Foot Ankle Res. 11. doi:10.1186/s13047-018-0295-9
Vicon Motion Systems Limited (2016). Plug-in Gait Reference Guide. Available at: https://docs.vicon.com/display/Nexus212/Full+body+modeling+with+Plug-in+Gait (Accessed February 13, 2022).
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
Watari, R., Suda, E. Y., Santos, J. P. S., Matias, A. B., Taddei, U. T., and Sacco, I. C. N. (2021). Subgroups of Foot-Ankle Movement Patterns Can Influence the Responsiveness to a Foot-Core Exercise Program: A Hierarchical Cluster Analysis. Front. Bioeng. Biotechnol. 9. doi:10.3389/fbioe.2021.645710
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
Winter, D. A. (1990). Biomechanics and Motor Control of Human Movement. New York: John Wiley & Sons.
Wright, C. J., Arnold, B. L., Coffey, T. G., and Pidcoe, P. E. (2011). Repeatability of the Modified Oxford Foot Model during Gait in Healthy Adults. Gait Posture 33, 108–112. doi:10.1016/j.gaitpost.2010.10.084
Keywords: multi-segment foot model, inverse dynamics, plantar fascia, walking, healthy adults, motion analysis, foot kinematics, foot kinetics
Citation: Matsumoto Y, Ogihara N, Hanawa H, Kokubun T and Kanemura N (2022) Novel Multi-Segment Foot Model Incorporating Plantar Aponeurosis for Detailed Kinematic and Kinetic Analyses of the Foot With Application to Gait Studies. Front. Bioeng. Biotechnol. 10:894731. doi: 10.3389/fbioe.2022.894731
Received: 12 March 2022; Accepted: 01 June 2022;
Published: 24 June 2022.
Edited by:
Evie E. Vereecke, KU Leuven Kulak, BelgiumCopyright © 2022 Matsumoto, Ogihara, Hanawa, Kokubun and Kanemura. 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: Naohiko Kanemura, kanemura-naohiko@spu.ac.jp