Next Article in Journal
Udder Health Monitoring for Prevention of Bovine Mastitis and Improvement of Milk Quality
Next Article in Special Issue
The Effect of Arch Stiffness on the Foot–Ankle Temporal Kinematics during Gait Termination: A Statistical Nonparametric Mapping Study
Previous Article in Journal
Efficient and Economic Heparin Recovery from Porcine Intestinal Mucosa Using Quaternary Ammonium-Functionalized Silica Gel
Previous Article in Special Issue
The Implications of Sports Biomechanics Studies on the Research and Development of Running Shoes: A Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Running Velocity and Longitudinal Bending Stiffness Influence the Asymmetry of Kinematic Variables of the Lower Limb Joints

1
Faculty of Sports Science, Ningbo University, Ningbo 315211, China
2
Faculty of Engineering, University of Szeged, 6720 Szeged, Hungary
3
Savaria Institute of Technology, Faculty of Informatics, Eötvös Loránd University, 9700 Szombathely, Hungary
*
Authors to whom correspondence should be addressed.
Bioengineering 2022, 9(11), 607; https://doi.org/10.3390/bioengineering9110607
Submission received: 11 September 2022 / Revised: 10 October 2022 / Accepted: 19 October 2022 / Published: 23 October 2022

Abstract

:
Running-related limb asymmetries suggest specific sports injuries and recovery circumstances. It is debatable if running speed affected asymmetry, and more research is required to determine how longitudinal bending stiffness (LBS) affected asymmetry. The purpose of this study was to investigate the influence of running velocity and LBS on kinematic characteristics of the hip, knee, ankle, metatarsophalangeal joint (MTP) and the corresponding asymmetry. Kinematic (200 Hz) running stance phase data were collected bilaterally for 16 healthy male recreational runners (age: 23.13 ± 1.17, height: 175.2 ± 1.6 cm, body mass: 75.7 ± 3.6 kg, BMI: 24.7 ± 1.3 kg/m2) running on a force plate at three different velocities (10, 12 and 14 km/h) and three increasing-LBS shoes in a randomized order. The symmetry angle (SA) was calculated to quantify gait asymmetry magnitude at each running velocity and LBS. Changes in running velocity and LBS led to differences in kinematic variables between the hip, knee, ankle and MTP (p < 0.05). Significant changes in SA caused by running velocity were found in the knee flexion angle (p = 0.001) and flexion angle peak velocity (p < 0.001), ankle plantarflexion angle (p = 0.001) and plantarflexion angle peak velocity (p = 0.043) and MTP dorsiflexion angle (p = 0.001) and dorsiflexion angle peak velocity (p = 0.019). A significant change in the SA caused by LBS was found in the MTP dorsiflexion peak angle velocity (p = 0.014). There were interaction effects between running velocity and LBS on the MTP plantarflexion angle (p = 0.033) and plantarflexion angle peak velocity (p = 0.038). These findings indicate the existence of bilateral lower limb asymmetry. Meanwhile, it was proved that running velocity and LBS can influence the asymmetry of lower limb joints. Additionally, there was an interaction between running velocity and LBS on the asymmetry of the lower limb. These findings can provide some information for sports injuries, such as metatarsal stress fractures and anterior cruciate ligament injuries. They can also provide some useful information for running velocities and running shoes.

1. Introduction

Longitudinal bending stiffness (LBS) is a footwear property that is known as one key consideration for footwear performance development [1,2]. A suitable LBS is also thought to be a crucial element in comfort and athletic performance [2,3]. The fast development of footwear technology has increased awareness of LBS. The majority of studies have also found a link between LBS and sports injuries [1,4,5]. Increasing the forefoot bending stiffness of footwear may help with injury prevention in particular by preventing excessive forefoot extension during sports activities to reduce the incidence of injuries such as turf toe [4]. Studies suggest that the stiffness of the forefoot bending may have an effect on the prevalence of metatarsal stress fractures. Stiff shoes cause the center of pressure under the foot to shift more anteriorly and alter the peak pressures acting on various foot regions [1]. Performance in long-distance running is influenced by running economy (RE) [6]. In a study by Roy and Stefanyshyn et al. [7], approximately 1% metabolic energy savings were observed when participants ran in a stiff midsole. In addition, studies had shown that stiff carbon-fiber plates may reduce the energetic cost of running by 4% [8]. These suggested that enhancing the midsole’s LBS would enhance RE.
The bulk of research overlooked the asymmetry between the bilateral lower limbs and instead used the unilateral dominant leg to represent the overall performance of the bilateral lower limbs [9] in order to simplify data collection and processing. The human body’s symmetry, however, is not perfect. There are asymmetries between the dominant and non-dominant legs, according to studies [10]. Lower limb asymmetry was not only caused by genetics and hormones, but also by biomechanical variables [11,12]. The differences in the load-related dynamics of the bilateral limbs may lead to asymmetry [12]. The existence of lower limb asymmetry was demonstrated in an experiment that included jumping–landing exercises. The knee moment of the dominant leg was greater than that of the non-dominant limb [13]. Additionally, it has been shown that knee flexion moment asymmetry predicts re-injury in athletes who had an anterior cruciate ligament reconstruction [14]. In earlier studies, the symmetry angle (SA), which reflects the symmetry of kinematic and kinetic variables of lower extremity joints, was frequently used [15]. In comparison to pre-fatigue, the SA of the knee flexion angle, hip flexion angle and hip extension angle was significantly higher in post-fatigue. In other words, after running, the joint asymmetry of the lower extremities becomes worse [16].
The biomechanical asymmetry of the lower limb can provide some information related to sports injuries. Idiopathic scoliosis may be caused by pelvic and hip angle asymmetry [17]. Kotwicki et al. [18], in a comparison of hip range of motion (ROM) asymmetry between scoliosis and normal adolescents, found that the scoliosis group showed greater asymmetry. For evaluating athletes’ return to the field following anterior cruciate ligament surgery and recuperation, the asymmetry of knee ROM provides important reference data [19]. Although asymmetry has often been often considered a manifestation of pathology, for some lower limb joint movements, the asymmetry in the range of variation of bilateral lower limbs in healthy people remains to be investigated [20].
Lower limb asymmetry has been shown to increase with walking velocity, suggesting that running may result in even more asymmetry. Biomechanical asymmetry is not detrimental during walking tasks; however, greater biomechanical demands are imposed on the musculoskeletal system during running [9]. Mo et al. discovered that among recreational runners, SA altered nonlinearly and displayed an approximately U-shaped trend across velocities [21]. Running velocity, in contrast, had no impact on the asymmetry of the kinematic characteristics of the joints in the lower limbs, according to Jiang et al. [11]. The negative work performed by the metatarsophalangeal joint (MTP) decreased with the increase in LBS and significantly changed the mechanical properties of the ankle and knee [22]. The MTP dorsiflexion angle and dorsiflexion angle velocity reduced with an increase in LBS to lessen the chance of forefoot injury [1,23]. Research on the kinematic variations of lower limb joints regarding LBS is very extensive. Research on the impact of LBS on lower limb asymmetry is lacking, however. Additionally, the impact of running speed on the kinematics of the lower limbs has always been debatable.
As a result, the purpose of this study was to explore the difference and asymmetry in lower limb kinematic variables when wearing increasing-LBS shoes at different running velocities. We analyzed the influence of running velocity and LBS on angles, angle peak velocities and the SA of the hip, knee, ankle and MTP.

2. Materials and Methods

2.1. Participants

The sample size was calculated using G*Power 3.1 (Franz Faul, Germany) for univariate analysis of variance for detecting a medium Cohen’s effect size (d = 0.4), α error probability = 0.05 and power (1 − β) = 0.95. Based on these parameters, it was estimated that a minimum of 14 participants would be required for this study [24]. A total of 16 healthy males (age: 23.13 ± 1.17, height: 175.2 ± 1.6 cm, body mass: 75.7 ± 3.6 kg, BMI: 24.7 ± 1.3 kg/m2) who were recreational runners (no formal running competition, training at least 3 times a week) were recruited [11,25]. The recruitment criteria for recreational runners in this experiment were running for at least 6 months and running a minimum distance of 10 km per week and having the right-side limb as the dominant limb [26]. The dominant limb was defined as the preferred leg when kicking a ball. All participants were free from health problems and/or neuromuscular disorders and/or known gait impairments, and had had no lower limb injuries in the previous six months. All participants were rearfoot strikers and were recruited from Ningbo University for this study. Before the experiment, all participants gave written consent. The study was approved by the Ethics Committee of the Research Institute at Ningbo University.

2.2. The Experimental Process

The general process of the experiment is shown in Figure 1b. Before the formal test, participants had 10 min to warm up and familiarize themselves with experimental settings. In the first step of the formal test, each participant was asked to stand on a force platform to collect static coordinates by standing parallel to the Y-axis of the force platform with arms crossed over shoulders and eyes looking forward until the full static coordinates were captured. All participants were allowed three trials to familiarize themselves with the test maneuvers before the formal test. During the test, participants were asked to wear Shoe 1 (S1), Shoe 2 (S2) and Shoe 3 (S3) at 10 ± 5% km/h (V1), 12 ± 5% km/h (V2) and 14 ± 5% km/h (V3), respectively, over a 10 m track [5,11,15,27]. The shoe information is shown in Table 1. Additionally, participants were asked to complete a full gait cycle on a 2 m force platform (Kistler, Winterthur, Switzerland) located in the middle of the track. Five trials were achieved to gather eligible data on the dominant leg, in which the running speed of the participant had less than 5% variance and was within 5% of the predefined running speed. The full gait cycle was defined as the time from the right heel strike to the left forefoot coming off the ground in this test. The LBS values of the shoes were measured by a rotational axis material-testing machine (Instron ElectroPuls E1000, Norwood, MA, USA). The force platform recorded the ground reaction force at 1000 Hz to distinguish a complete gait cycle. An eight-camera motion capture system (Vicon Metrics Ltd., Oxford, United Kingdom) was used to record running kinematic data during the stance phase at a frequency of 200 Hz. A threshold of 20 N on the vertical ground reaction force was applied to identify the initial foot contact and toe-off. To manage running velocity, Brower timing lights (Brower Timing System, Draper, UT, USA) were used. Before the experiment, participants were asked to apply 38 reflective markers (diameter: 14 mm) on their bodies. The specific positions of the markers are shown in Figure 1a.

2.3. Data Analysis

This study focused on the sagittal planes of the hip, knee, ankle and MTP. Dominant variation in the sagittal plane is reported to occur during running [28]. Marker trajectories were filtered by zero-latency fourth-order Butterworth low-pass filters at 12 Hz. The C3D file data were converted to formats recognized in OpenSim 4.3 (.mot and .trc) by Matlab R2016a (The MathWorks, Natick, MA, USA), and then imported into OpenSim for data processing [29]. A musculoskeletal model in OpenSim (gait 2392) was used. The model was scaled using the participant’s marker point location and weight in a static calibration. The static weight of each marker was manually adjusted according to the root mean square (RMS) error value (less than 0.02) between the experimental and virtual markers in the model until it was adjusted to the appropriate position before applying the scaled model to the data calculation. The joint angles were calculated using the inverse kinematics (IK) calculation tool in OpenSim, and the results were optimized using least squares to minimize the error between the experimental and virtual markers.
SA was used to evaluate the biomechanical symmetry of the participants’ dominant limb and non-dominant limb. SA can be calculated as follows:
SA = ( 45 ° arctan ( X l e f t / X r i g h t ) ) 90 ° × 100 %
If
( 45 ° arctan ( X l e f t / X r i g h t ) )
Then
SA = ( 45 ° arctan ( X l e f t / X r i g h t ) 180 ° ) 90 ° × 100 %
Xleft represents the kinematic variables of the left lower limb, and Xright represents the kinematic variables of the right lower limb. A score of 0% suggests perfect symmetry and 100% suggests perfect asymmetry between the right and left leg [30].

2.4. Statistical Analysis

SPSS 26.0 (SPSS, Chicago, IL, USA) software was used for statistical analysis.
Descriptive statistics were provided as means and standard deviations (SDs). Tests for normality and homogeneity of variances (Shapiro–Wilk and Levene’s, respectively) were conducted on all SA data before the analysis. Univariate ANOVA was performed to determine the effects of running velocity and LBS on bilateral joint angle and angle peak velocity asymmetry during running. In the presence of interaction, simple effect comparisons were performed. Pairwise comparisons with Bonferroni were used post hoc to further analyze significant effects of running velocity, LBS and interaction. Effect sizes were calculated using partial eta2 (η2), with the relative magnitude of any differences expressed as a standard criterion: small effect size (0.01 < η2 < 0.05), medium effect size (0.06 < η2 < 0.14) and large effect size (η2 > 0.14) [31,32,33]. Paired t-tests assessed differences in joint angles and angle peak velocities between legs at different running velocities and in increasing-LBS shoes. This was primarily used to investigate if running at different velocities while wearing increasing-LBS shoes caused variations in the kinematic characteristics of bilateral lower limbs. The level of statistical significance was set at p < 0.05.
Joint angles and joint angle velocities during the running stance phase were compared using one-dimension statistical parameter mapping (SPM1d). For SPM1d, kinematic data for each step were time-normalized to the stance phase (101 data points per stance phase). Biomechanical differences between the right and left lower limb of each running velocity and LBS were obtained by statistically examining the entire time series using SPM1d with post hoc paired t-tests. All SPM1d analyses were conducted in MATLAB R2016a (The MathWorks, MA, USA) using the open-source software package spm1d 8 [34] (www.spm1d.org, accessed on 2 May 2022).

3. Results

The effects of running velocity and LBS on kinematic variables of lower limb joints and SA were shown in Table 2 and Table 3, respectively. Figure 2 and Figure 3 show the angle and angle velocity changes in the hip, knee, ankle and MTP.

3.1. The Influence of Running Velocity and LBS on Lower Limb Joints’ Kinematic Variables

3.1.1. Hip

There were significant differences in the hip extension angle when the participants ran at V1 wearing S1, S2 and S3 (p < 0.05). The extension angle and extension angle peak velocity were significantly different when they ran at V2 wearing S3 (p < 0.05).

3.1.2. Knee

The flexion angle was significantly different when the participants ran at V3 wearing S1, S2 and S3 (p < 0.05). The extension angle was significantly different when they ran at V1 in S1, S2 and S3 (p < 0.05). It was significantly different when they ran at V2 in S1 and S3 (p < 0.05). There was a significant difference when they ran at V3 in S3 (p < 0.05). The flexion angle peak velocity was significantly different when they ran at V1 wearing S3 (p < 0.05). The extension angle peak velocity was significantly different when they ran at V2 and V3 wearing S1 and at V1 and V2 wearing S3 (p < 0.05).

3.1.3. Ankle

The ankle dorsiflexion angle was significantly different when the participants ran at V1, V2 and V3 wearing S1, S2 and S3 (p < 0.05). The plantarflexion angle was significantly different when they ran at V1 wearing S1 and S2 and at V2 wearing S2 (p < 0.05). The dorsiflexion angle peak velocity was significantly different when they ran at V1 wearing S3 and at V2 wearing S2 (p < 0.05). The plantarflexion angle peak velocity was significantly different when they ran at V1 wearing S2 and at V2 wearing S2 and S3 (p < 0.05).

3.1.4. MTP

The dorsiflexion angle was significantly different when the participants ran at V1 wearing S3, at V2 wearing S1 and S3 and at V3 wearing S1, S2 and S3 (p < 0.05). The plantarflexion angle was significantly different when they ran at V2 and V3 wearing S2 and S3, respectively (p < 0.05). The dorsiflexion angle peak velocity was significantly different when they ran at V3 wearing S2 and S3 (p < 0.05). The plantarflexion angle peak velocity was significantly different when they ran at V1 wearing S1 and ran at V3 wearing S3 (p < 0.05).

3.1.5. SPM Results

The SPM results are displayed as:
(1) Hip angles were significantly different in the phases of 0–27%, 46–53% and 58–100% when the participants ran at V1 and wore S1; the phases of 11–28% and 37–55% when they ran at V1 and wore S3; the phase of 0–85% when they ran at V2 and wore S1; the phase of 0–75% when they ran at V2 and wore S2; the phase of 18–70% when they ran at V2 and wore S3; and the phase of 42–70% when they ran at V3 and wore S2. Hip angle velocities were significantly different in the phase of 35–37% when they ran at V1 and wore S1; the phases of 21–22% and 76–77% when they ran at V1 and wore S2; the phases of 35–38%, 62–65%, 66–76%, 78–81%, 83–85% and 88–90% when they ran at V1 and wore S3; the phases of 38–42%, 65–68% and 71–72% when they ran at V2 and wore S1; the phases of 90–94% and 97–100% when they ran at V2 and wore S2; the phases of 65–77% and 80–82% when they ran at V2 and wore S3; the phases of 43–47% and 58–63% when they ran at V3 and wore S1; the phase of 44–52% when they ran at V3 and wore S2; and the phase of 45–55% when they ran at V3 and wore S3.
(2) Knee angles were significantly different in the phase of 98–100% when the participants ran at V1 and wore S3; the phase of 11–65% when they ran at V2 and wore S1; the phase of 0–3% when they ran at V2 and wore S2; and the phases of 0–8% and 88–100% when they ran at V3 and wore S1. Knee angle velocities were significantly different in the phases of 21–22% and 34–35% when they ran at V1 and wore S2; the phases of 44–45%, 58–68% and 71–73% when they ran at V1 and wore S3; the phases of 40–44% and 64–76% when they ran at V2 and wore S1; the phase of 41–49% when they ran at V3 and wore S2; and the phase of 26–31% when they ran at V3 and wore S3.
(3) Ankle angles were significantly different in the phases of 0–77%, 45–54% and 57–100% when the participants ran at V1 and wore S1; the phase of 0–11% when they ran at V1 and wore S3; the phase of 13–44% when they ran at V2 and wore S1; the phase of 94–100% when they ran at V2 and wore S2; the phases of 1–13%, 15–17% and 22–26% when they ran at V2 and wore S3; and the phase of 5–9% when they ran at V3 and wore S1. Ankle angle velocities were significantly different in the phases of 22–25%, 27–31% and 84–85% when they ran at V1 and wore S1; the phases of 74–75% and 77–82% when they ran at V1 and wore S2; the phases of 13–14%, 52–53% and 57–65% when ran at V2 and wore S1; the phase of 88–91% when they ran at V2 and wore S2; the phase of 42–48% when they ran at V2 and wore S3; and the phase of 37–40% when they ran at V3 and wore S3.
(4) MTP angles were significantly different in the phase of 0–100% when the participants ran at V1 and wore S2; the phase of 80–94% when they ran at V1 and wore S3; the phase of 75–100% when they ran at V2 and wore S1; the phase of 0–100% when they ran at V2 and S2; the phases of 4–12% and 58–100% when they ran at V2 and wore S3; the phase of 85–95% when they ran at V3 and wore S1; the phase of 6–100% when they ran at V3 and wore S2; and the phase of 72–93% when they ran at V3 and wore S3. MTP angle velocities were significantly different in the phase of 32–41% when they ran at V1 and wore S3; the phase of 35–48% when they ran at V2 and wore S3; the phase of 0–11% when they ran at V3 and wore S1; the phase of 10–13% when they ran at V3 and wore S2; and the phase of 0–10% when they ran at V3 and wore S3.

3.2. The Influence of Running Velocity and LBS on SA

3.2.1. Running Velocity

The knee flexion angle and flexion angle peak velocity, ankle plantarflexion angle and plantarflexion angle peak velocity, MTP dorsiflexion angle and dorsiflexion angle peak velocity were significantly different due to running velocity (p < 0.05). There was no interaction effect with LBS. This suggests that running velocity affected the SA of lower limb joints, regardless of the type of shoes participants wore. The statistical difference was specific as (1) the knee flexion angle was significantly different when the participants ran at V3 compared to V1 and V2 (p < 0.05). The knee flexion angle peak velocity was significantly different at V1, V2 and V3 (p < 0.05). (2) The ankle plantarflexion angle in V3 was significantly different from that in V1 and V2 (p < 0.05). The plantarflexion peak velocity in V1 was also significantly different from that in V3 (p < 0.05). (3) The MTP dorsiflexion angle was significantly different at V1, V2 and V3 (p < 0.05). The plantarflexion angle peak velocity was significantly different between V1 and V3 (p < 0.05).

3.2.2. LBS

LBS had a significant effect on the MTP dorsiflexion peak velocity, and there was no interaction effect with running velocity (p < 0.05). This indicates that running velocity would not affect the influence of LBS on SA. MTP plantarflexion angle peak velocity was significantly different when worn in S1 and S3.

3.2.3. The Interaction between Running Velocity and LBS

Interaction effects were observed in MTP. There was an interaction effect between running velocity and LBS in both the MTP plantarflexion angle and plantarflexion angle peak velocity. When ran at V2, the SA corresponding to the plantarflexion angle was significantly different between S1 and S2 shoes (p < 0.05). The SA of the plantarflexion angle peak velocity was significantly different between the participants ran at V1 and V2 when wearing S3 (p < 0.05). In addition, there was also a significant difference in the plantarflexion angle peak velocity between S2 and S3 when the participants ran at V1 (p < 0.05).

4. Discussion

This study aimed to investigate how bilateral lower limb asymmetry evolves when wearing increasing-LBS shoes at different running velocities. Asymmetry was present as evidenced by differences in angle and angle peak velocity between the hip, knee, ankle, MTP of both lower limbs and the SPM results. It also further confirmed the previous research [9,11] showing that differences exist in bilateral lower limbs. On this basis, it was vital to investigate the changes in the SA. In the exploration of the SA, we found that running velocity influenced the knee flexion angle and flexion angle peak velocity, ankle plantarflexion angle and plantarflexion angle peak velocity, MTP dorsiflexion angle and dorsiflexion angle peak velocity. Furthermore, LBS had a significant influence on the MTP dorsiflexion angle peak velocity. The interaction between running velocity and LBS existed in the MTP plantarflexion angle and plantarflexion angle peak velocity. It was found that the difference in kinematic variables of lower limb joints was not completely consistent with the change in SA. This was also reflected in the study by Jiang et al. and Gao et al. [11,12].
Studies on the influence of the running velocity and LBS on kinematic variables suggested that there were differences in the hip. However, the SA of the hip did not change significantly with running velocity and LBS. In addition, the SA of the hip was generally smaller than that of other joints, suggesting a higher symmetry. Hannah et al. pointed out that a high level of symmetry was observed for the sagittal plane of the hip [35].
The difference existed both In knee angle and angle peak velocity. In addition, the influence of running velocity on the SA of the knee mainly focused on the knee flexion angle and knee flexion angle peak velocity. In their investigation, Rodolfo et al. [36] found that there were differences in the strength of the knee flexors in the dominant and non-dominant legs. According to Ayako Higashihara et al. [37], the characteristics of the electromyographic activity of the biceps femoris and semitendinosus muscles changed as running velocity increased. Additionally, the biceps femoris activated much earlier during the stance phase than the semitendinosus. The difference between flexor and extensor muscles may be the cause of the alteration in knee kinematics, which results in the alteration in the SA with velocity. The majority of injuries sustained during running occur around or close to the knee, especially in overuse conditions, including patellar tendinitis, meniscal tears, tibial stress fractures and patellofemoral pain syndrome [38]. The knee kinematic difference is frequently used by medical professionals, physiotherapists or strength and conditioning specialists to quantify the functional deficit caused by a knee injury and/or surgery, like an anterior cruciate ligament injury, to track the effectiveness of sports rehabilitation programs and establish baseline data for readiness to return to play after an injury [39].
The influence of running velocity on the SA of the ankle was primarily focused on the plantarflexion angle and plantarflexion angle peak velocity. As the largest articular joint of the distal lower limb, differences in side-to-side force production may influence joint kinematics, potentially increasing the injury risk [40]. Additionally, the changes in asymmetries might be seen as a compensatory mechanism to keep the gait stable when running [11]. As a consequence, changes in the difference and SA of the ankle plantarflexion angle and plantarflexion angle peak velocity may serve as a mechanism of self-protection during running. When the participants ran at V1, the SA corresponding to the ankle plantarflexion angle and plantarflexion angle peak velocity were generally smaller than when the participants ran at V2 and V3. This may imply that slower velocities are preferable for recreational runners.
Running velocity influenced the SA of the MTP dorsiflexion angle and dorsiflexion angle peak velocity significantly. Additionally, the SA generally increased with the increase in running velocity. Limb dominance is believed to be task-dependent. It has been proposed that the functions of the lower limbs (such as those seen in muscle activity and joint kinematics) lead to a local asymmetry [41]. Asymmetry in the plantar pressure of the MTP as well as functional differences between the dominant and non-dominant legs may be to blame for changes in the SA of the MTP dorsiflexion angle and dorsiflexion angle peak velocity. In the study of plantar pressure symmetry by Gao et al., it was found that the metatarsal pressure and the corresponding SA increased with increasing velocity [12]. Furthermore, greater dorsiflexion was one of the factors contributing to elevated metatarsal pressure [42]. Changes in the MTP dorsiflexion angle and dorsiflexion angle peak velocity may be related to changes in the SA due to metatarsal pressure. MTP serves as the foundation of support once the heel has lifted while running [43]. Moving at a higher speed requires the MTP to create significantly greater propulsion [44]. Therefore, the MTP is significantly affected by running velocity. The asymmetry of the MTP may be the underlying cause of metatarsal stress fractures [45]. Additionally, we discovered that, typically speaking, the SA rose as running velocity rose. This might imply that casual runners are more suited to running at a slower pace. Additionally, the MTP dorsiflexion angle peak velocity was mostly where LBS had an effect on the SA. Furthermore, the SA increased with the increase in LBS. Willwacher et al. [46] demonstrated that increasing the longitudinal bending stiffness in running shoes led to a significant anterior shift of all lower limb lever arms, and this effect was greater at the more distal joints. This could be the cause of LBS’s effects in this investigation, which were limited to the MTP’s SA. Gait asymmetry in healthy individuals has been referred to as functional asymmetry, which describes limb asymmetry based on the distinction of fundamental functions, namely stabilization or propulsion [20]. Differences in muscle function and strength could be the source of the asymmetry of the bilateral MTP. Interestingly, there were interaction effects between LBS and running velocity on the MTP plantarflexion angle and plantarflexion angle peak velocity. The result for the MTP plantarflexion angle showed that when the participants ran at V2, the SA corresponding to S2 was significantly smaller than that corresponding to S1. In addition, the result for MTP plantarflexion angle peak velocity showed that the SA corresponding to S2 was significantly smaller than that corresponding to S3 when the participants ran at V1. When they ran in S3, the SA corresponding to V2 was significantly smaller than that corresponding to V1. In their study, Cornelis et al. [47] found that recreational runners performed better when they ran at velocities between 10 and 13 km/h, which is broadly comparable to the improved symmetry in our results for running at low and moderate velocities. Increased LBS can protect the MTP by limiting the plantarflexion of the MTP. However, larger LBS required greater muscle strength to overcome the stiff carbon-fiber plate [8]. When the threshold is exceeded, running performance may deteriorate [48]. This suggests that there may be an optimal LBS. In addition, the results of this study show that the influence of LBS on the SA was nearly U-shaped when run at low or moderate velocity. This further indicates that there was an optimal LBS for the influence of LBS on the SA of the MTP. This could mean that recreational runners would benefit more from using shoes with modest LBS when running at low-to-moderate velocities since they would have better symmetry and a better chance of avoiding sports injuries.
Although our study is very novel, there are some limitations. Only kinematic data of lower extremity joints were discussed in this study. The asymmetry of the dynamic variables of the lower limb joints was not discussed. The following research can explore the asymmetry of lower limb dynamics variables and expand the sample to include novice runners and competitive runners.

5. Conclusions

These findings indicate the existence of bilateral lower limb asymmetry. Meanwhile, it has been established that running velocity and LBS affect the asymmetry of the lower limb. Additionally, there was an interaction between running velocity and LBS on the asymmetry of lower limb. These findings can provide some information for sports injuries, such as metatarsal stress fractures and anterior cruciate ligament injuries. They can also provide some useful suggestions for running velocities and running shoes.

Author Contributions

Q.L., N.A. and H.C. conceived the presented idea, developed the framework and wrote the manuscript. Y.S., G.F., J.L. and Y.G. provided critical feedback and contributed to the final version. All authors were involved in the final direction of the paper and contributed to the final version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Major Program of the National Natural Science Foundation of China (19ZDA352), Zhejiang Provincial Key Research and Development Program of China (2021C03130), Zhejiang Provincial Natural Science Foundation of China for Distinguished Young Scholars (LR22A020002), Philosophy and Social Sciences Project of Zhejiang Province, China (22QNYC10ZD, 22NDQN223YB), Educational science planning project of Zhejiang Province (2021SCG083), the Fundamental Research Funds for the Provincial Universities of Zhejiang (SJWY2022014), Ningbo Natural Science Foundation (20221JCGY010532; 20221JCGY010607), and K. C. Wong Magna Fund in Ningbo University.

Institutional Review Board Statement

This research was approved by the Ethics Committee of Ningbo University (approval number: RAGH20220420).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data that support the findings of this study are available on reasonable request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stefanyshyn, D.J.; Wannop, J.W. The influence of forefoot bending stiffness of footwear on athletic injury and performance. Footwear Sci. 2016, 8, 51–63. [Google Scholar] [CrossRef]
  2. Park, S.K.; Lam, W.K.; Yoon, S.; Lee, K.K.; Ryu, J. Effects of forefoot bending stiffness of badminton shoes on agility, comfort perception and lower leg kinematics during typical badminton movements. Sports Biomech. 2017, 16, 374–386. [Google Scholar] [CrossRef] [PubMed]
  3. Jiang, C. The Effect of Basketball Shoe Collar on Ankle Stability: A Systematic Review and Meta-Analysis. Phys. Act. Health 2020, 4, 11–18. [Google Scholar] [CrossRef] [Green Version]
  4. Wannop, J.W.; Stefanyshyn, D.J. Special issue: Bending stiffness: Performance and injury effects. Footwear Sci. 2016, 8, 49–50. [Google Scholar] [CrossRef] [Green Version]
  5. Chen, H.; Shao, E.; Sun, D.; Xuan, R.; Baker, J.S.; Gu, Y. Effects of footwear with different longitudinal bending stiffness on biomechanical characteristics and muscular mechanics of lower limbs in adolescent runners. Front. Physiol. 2022, 8, 1499. [Google Scholar] [CrossRef] [PubMed]
  6. Cheung, R.T.; Ngai, S.P. Effects of footwear on running economy in distance runners: A meta-analytical review. J. Sci. Med. Sport 2016, 19, 260–266. [Google Scholar] [CrossRef]
  7. Roy, J.P.; Stefanyshyn, D.J. Shoe midsole longitudinal bending stiffness and running economy, joint energy, and EMG. Med. Sci. Sports Exerc. 2006, 38, 562–569. [Google Scholar] [CrossRef] [Green Version]
  8. Day, E.; Hahn, M. Optimal footwear longitudinal bending stiffness to improve running economy is speed dependent. Footwear Sci. 2020, 12, 3–13. [Google Scholar] [CrossRef]
  9. Radzak, K.N.; Putnam, A.M.; Tamura, K.; Hetzler, R.K.; Stickley, C.D. Asymmetry between lower limbs during rested and fatigued state running gait in healthy individuals. Gait Posture 2017, 51, 268–274. [Google Scholar] [CrossRef]
  10. Yen, S.C.; Olsavsky, L.C.; Cloonan, C.M.; Llanos, A.R.; Dwyer, K.J.; Nabian, M.; Farjadian, A.B. An examination of lower limb asymmetry in ankle isometric force control. Hum. Mov. Sci. 2018, 57, 40–49. [Google Scholar] [CrossRef]
  11. Xinyan, J.; Hairong, C.; Dong, S.; Baker, J.S.; Yaodong, G. Running speed does not influence the asymmetry of kinematic variables of the lower limb joints in novice runners. Acta Bioeng. Biomech. 2021, 23, 69–81. [Google Scholar] [CrossRef]
  12. Gao, Z.; Mei, Q.; Xiang, L.; Baker, J.S.; Fernandez, J.; Gu, Y. Effects of limb dominance on the symmetrical distribution of plantar loading during walking and running. Proc. Inst. Mech. Eng. Part P J. Sports Eng. Technol. 2020, 236, 17–23. [Google Scholar] [CrossRef]
  13. Pappas, E.; Carpes, F.P. Lower extremity kinematic asymmetry in male and female athletes performing jump-landing tasks. J. Sci. Med. Sport 2012, 15, 87–92. [Google Scholar] [CrossRef]
  14. Paterno, M.V.; Schmitt, L.C.; Ford, K.R.; Rauh, M.J.; Myer, G.D.; Huang, B.; Hewett, T.E. Biomechanical measures during landing and postural stability predict second anterior cruciate ligament injury after anterior cruciate ligament reconstruction and return to sport. Am. J. Sports Med. 2010, 38, 1968–1978. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Girard, O.; Morin, J.-B.; Ryu, J.; Read, P.; Townsend, N. Running velocity does not influence lower limb mechanical asymmetry. Front. Sports Act. Living 2019, 1, 36. [Google Scholar] [CrossRef] [Green Version]
  16. Gao, Z.; Mei, Q.; Fekete, G.; Baker, J.S.; Gu, Y. The Effect of Prolonged Running on the Symmetry of Biomechanical Variables of the Lower Limb Joints. Symmetry 2020, 12, 720. [Google Scholar] [CrossRef]
  17. Karatel, M.; Yagci, G.; Yakut, Y. Investigation of multidirectional hip range of motion and hip motion asymmetry in individuals with idiopathic scoliosis with different curve patterns. J. Bodyw. Mov. Ther. 2021, 27, 77–83. [Google Scholar] [CrossRef]
  18. Kotwicki, T.; Walczak, A.; Szulc, A. Trunk rotation and hip joint range of rotation in adolescent girls with idiopathic scoliosis: Does the “dinner plate” turn asymmetrically? Scoliosis 2008, 3, 1. [Google Scholar] [CrossRef] [Green Version]
  19. Renner, K.E.; Franck, C.T.; Miller, T.K.; Queen, R.M. Limb asymmetry during recovery from anterior cruciate ligament reconstruction. J. Orthop. Res.® 2018, 36, 1887–1893. [Google Scholar] [CrossRef] [Green Version]
  20. Sadeghi, H.; Allard, P.; Prince, F.; Labelle, H. Symmetry and limb dominance in able-bodied gait: A review. Gait Posture 2000, 12, 34–45. [Google Scholar] [CrossRef]
  21. Mo, S.; Lau, F.O.Y.; Lok, A.K.Y.; Chan, Z.Y.S.; Zhang, J.H.; Shum, G.; Cheung, R.T.H. Bilateral asymmetry of running gait in competitive, recreational and novice runners at different speeds. Hum. Mov. Sci. 2020, 71, 102600. [Google Scholar] [CrossRef] [PubMed]
  22. Healey, L.A.; Hoogkamer, W. Longitudinal bending stiffness does not affect running economy in Nike Vaporfly Shoes. J. Sport Health Sci. 2022, 11, 285–292. [Google Scholar] [CrossRef] [PubMed]
  23. Ortega, J.A.; Healey, L.A.; Swinnen, W.; Hoogkamer, W. Energetics and Biomechanics of Running Footwear with Increased Longitudinal Bending Stiffness: A Narrative Review. Sports Med. 2021, 51, 873–894. [Google Scholar] [CrossRef] [PubMed]
  24. Tartibian, B.; Malandish, A.; Afsargharehbagh, R.; Eslami, R.; Sheikhlou, Z. Assessment of hepatic and lipid profiles following 12 weeks of aerobic exercise in overweight postmenopausal women. Int. J. Basic Sci. Med. 2018, 3, 159–167. [Google Scholar] [CrossRef]
  25. Delp, S.L.; Loan, J.P.; Hoy, M.G.; Zajac, F.E.; Topp, E.L.; Rosen, J.M. An interactive graphics-based model of the lower extremity to study orthopaedic surgical procedures. IEEE Trans. Biomed. Eng. 1990, 37, 757–767. [Google Scholar] [CrossRef]
  26. Saragiotto, B.T.; Yamato, T.P.; Lopes, A.D. What do recreational runners think about risk factors for running injuries? A descriptive study of their beliefs and opinions. J. Orthop. Sports Phys. Ther. 2014, 44, 733–738. [Google Scholar] [CrossRef] [Green Version]
  27. Monte, A.; Baltzopoulos, V.; Maganaris, C.N.; Zamparo, P. Gastrocnemius Medialis and Vastus Lateralis in vivo muscle-tendon behavior during running at increasing speeds. Scand. J. Med. Sci. Sports 2020, 30, 1163–1176. [Google Scholar] [CrossRef]
  28. Sinclair, J.; Selfe, J. Sex differences in knee loading in recreational runners. J. Biomech. 2015, 48, 2171–2175. [Google Scholar] [CrossRef] [Green Version]
  29. Zhou, H.; Xu, D.; Quan, W.; Liang, M.; Ugbolue, U.C.; Baker, J.S.; Gu, Y. A Pilot Study of Muscle Force between Normal Shoes and Bionic Shoes during Men Walking and Running Stance Phase Using Opensim. Actuators 2021, 10, 274. [Google Scholar] [CrossRef]
  30. Zifchock, R.A.; Davis, I.; Higginson, J.; Royer, T. The symmetry angle: A novel, robust method of quantifying asymmetry. Gait Posture 2008, 27, 622–627. [Google Scholar] [CrossRef]
  31. Lu, Z.; Li, X.; Xuan, R.; Song, Y.; Bíró, I.; Liang, M.; Gu, Y. Effect of Heel Lift Insoles on Lower Extremity Muscle Activation and Joint Work during Barbell Squats. Bioengineering 2022, 9, 301. [Google Scholar] [CrossRef] [PubMed]
  32. Flores, N.; Delattre, N.; Berton, E.; Rao, G. Does an increase in energy return and/or longitudinal bending stiffness shoe features reduce the energetic cost of running? Eur. J. Appl. Physiol. 2019, 119, 429–439. [Google Scholar] [CrossRef] [PubMed]
  33. Day, E.M.; Hahn, M.E. Does running speed affect the response of joint level mechanics in non-rearfoot strike runners to footwear of varying longitudinal bending stiffness? Gait Posture 2021, 84, 187–191. [Google Scholar] [CrossRef]
  34. Pataky, T.C. One-dimensional statistical parametric mapping in Python. Comput. Methods Biomech. Biomed. Eng. 2012, 15, 295–301. [Google Scholar] [CrossRef] [PubMed]
  35. Hannah, R.; Morrison, J.; Chapman, A. Kinematic symmetry of the lower limbs. Arch. Phys. Med. Rehabil. 1984, 65, 155–158. [Google Scholar] [PubMed]
  36. Dellagrana, R.A.; Diefenthaeler, F.; Carpes, F.P.; Hernandez, S.G.; de Campos, W. Evidence for isokinetic knee torque asymmetries in male long distance-trained runners. Int. J. Sports Phys. Ther. 2015, 10, 514. [Google Scholar]
  37. Higashihara, A.; Ono, T.; Kubota, J.; Okuwaki, T.; Fukubayashi, T. Functional differences in the activity of the hamstring muscles with increasing running speed. J. Sports Sci. 2010, 28, 1085–1092. [Google Scholar] [CrossRef]
  38. Thakre, P.S.; Zakiuddin, K.S.; Khan, I.A.; Faizan, M.S. Finite element analysis of tibia bone. Int. J. Biomed. Eng. Technol. 2021, 35, 318–339. [Google Scholar] [CrossRef]
  39. Dempster, J.; Dutheil, F.; Ugbolue, U.C. The Prevalence of Lower Extremity Injuries in Running and Associated Risk Factors: A Systematic Review. Phys. Act. Health 2021, 5, 133–145. [Google Scholar] [CrossRef]
  40. Tiwari, A.; Singh, O.P.; Bhatia, D. A review on motor neuron disabilities and treatments. Int. J. Biomed. Eng. Technol. 2021, 37, 154–175. [Google Scholar] [CrossRef]
  41. Wafai, L.; Zayegh, A.; Woulfe, J.; Aziz, S.M.; Begg, R. Identification of foot pathologies based on plantar pressure asymmetry. Sensors 2015, 15, 20392–20408. [Google Scholar] [CrossRef] [PubMed]
  42. Nagel, A.; Fernholz, F.; Kibele, C.; Rosenbaum, D. Long distance running increases plantar pressures beneath the metatarsal heads: A barefoot walking investigation of 200 marathon runners. Gait Posture 2008, 27, 152–155. [Google Scholar] [CrossRef] [PubMed]
  43. Pg Damit, D.N.F.; Senanayake, S.M.N.A.; Malik, O.A.; Tuah, N.J. Integrated neuromuscular fatigue analysis system for soldiers’ load carriage trial using DWT. Int. J. Biomed. Eng. Technol. 2021, 35, 1–18. [Google Scholar] [CrossRef]
  44. Merry, C.; Baker, J.S.; Dutheil, F.; Ugbolue, U.C. Do Kinematic Study Assessments Improve Accuracy & Precision in Golf Putting? A Comparison between Elite and Amateur Golfers: A Systematic Review and Meta-Analysis. Phys. Act. Health 2022, 6, 108–123. [Google Scholar]
  45. Xu, D.; Quan, W.; Zhou, H.; Sun, D.; Baker, J.S.; Gu, Y. Explaining the differences of gait patterns between high and low-mileage runners with machine learning. Sci. Rep. 2022, 12, 2981. [Google Scholar] [CrossRef]
  46. Willwacher, S.; König, M.; Braunstein, B.; Goldmann, J.-P.; Brüggemann, G.-P. The gearing function of running shoe longitudinal bending stiffness. Gait Posture 2014, 40, 386–390. [Google Scholar] [CrossRef]
  47. De Ruiter, C.J.; Van Daal, S.; Van Dieën, J.H. Individual optimal step frequency during outdoor running. Eur. J. Sport Sci. 2020, 20, 182–190. [Google Scholar] [CrossRef]
  48. Xiang, L.; Mei, Q.; Wang, A.; Shim, V.; Fernandez, J.; Gu, Y. Evaluating function in the hallux valgus foot following a 12-week minimalist footwear intervention: A pilot computational analysis. J. Biomech. 2022, 132, 110941. [Google Scholar] [CrossRef]
Figure 1. (a) The front, side and back positions of markers. Blue dots: markers. (b) Illustration of experiment design for collecting the kinematics data during the running stance phase.
Figure 1. (a) The front, side and back positions of markers. Blue dots: markers. (b) Illustration of experiment design for collecting the kinematics data during the running stance phase.
Bioengineering 09 00607 g001
Figure 2. Descriptive results of angles between left lower limb and right lower limb SPM during running at V1, V2 and V3. L = left lower limb, R = right lower limb.
Figure 2. Descriptive results of angles between left lower limb and right lower limb SPM during running at V1, V2 and V3. L = left lower limb, R = right lower limb.
Bioengineering 09 00607 g002
Figure 3. Descriptive results of angle velocities between left lower limb and right lower limb SPM during running at V1, V2 and V3. L = left lower limb, R = right lower limb.
Figure 3. Descriptive results of angle velocities between left lower limb and right lower limb SPM during running at V1, V2 and V3. L = left lower limb, R = right lower limb.
Bioengineering 09 00607 g003
Table 1. The information of shoes.
Table 1. The information of shoes.
S1S2S3
LBS value (Nm/rad)2.75.08.6
AppearanceBioengineering 09 00607 i001Bioengineering 09 00607 i002Bioengineering 09 00607 i003
Table 2. The influence of running velocity and LBS on lower limb joints’ kinematic variables.
Table 2. The influence of running velocity and LBS on lower limb joints’ kinematic variables.
Joint KinematicsS1p-ValueS2p-ValueS3p-Value
LRLRLR
HipExtension angle (◦)V149.07 ± 10.1144.01 ± 8.920.03348.27 ± 9.1141.90 ± 7.830.01047.56 ± 9.8340.28 ± 8.50<0.001
V248.71 ± 14.8446.05 ± 10.950.09453.79 ± 2.4650.15 ± 8.820.21353.84 ± 10.3950.80 ± 9.410.037
V366.87 ± 7.2961.32 ± 10.410.14267.94 ± 7.5866.06 ± 9.410.45163.99 ± 8.3467.77 ± 7.390.287
Extension angle velocity (rad/s)V16.47 ± 0.446.16 ± 0.970.4335.97 ± 0.236.18 ± 0.460.1306.05 ± 0.335.94 ± 0.490.595
V26.99 ± 0.666.99 ± 0.560.9777.12 ± 0.606.99 ± 0.480.6997.61 ± 0.486.90 ± 0.650.004
V39.75 ± 0.889.37 ± 0.970.5199.49 ± 0.609.42 ± 0.500.8189.13 ± 0.939.45 ± 0.520.434
KneeFlexion angle (◦)V129.10 ± 4.1927.36 ± 4.340.13526.98 ± 2.4527.75 ± 2.670.21528.77 ± 3.9728.06 ± 4.850.689
V226.20 ± 3.0326.37 ± 4.630.86425.22 ± 3.5225.95 ± 1.410.57326.49 ± 6.0225.49 ± 2.440.653
V318.13 ± 2.5715.58 ± 4.220.01219.01 ± 2.0715.01 ± 2.350.00321.38 ± 3.0515.99 ± 2.860.007
Extension angle (◦)V125.06 ± 8.8119.03 ± 9.290.00522.16 ± 7.3719.23 ± 5.520.03423.04 ± 7.6215.91 ± 8.35<0.001
V222.56 ± 11.3617.80 ± 10.780.00220.16 ± 4.5518.28 ± 9.880.48820.70 ± 12.1417.66 ± 11.390.014
V322.22 ± 10.7816.42 ± 10.890.00723.97 ± 10.7322.70 ± 10.660.63221.21 ± 9.7625.48 ± 6.600.192
Flexion angle velocity (rad/s)V14.99 ± 1.064.12 ± 1.470.0584.28 ± 1.144.32 ± 1.280.9384.51 ± 1.113.37 ± 1.280.003
V29.39 ± 1.678.30 ± 1.830.0638.93 ± 1.668.31 ± 0.580.2769.58 ± 2.728.02 ± 0.560.122
V39.30 ± 2.798.38 ± 2.540.4829.72 ± 2.418.13 ± 1.900.23610.73 ± 2.188.12 ± 2.180.065
Extension angle velocity (rad/s)V14.99 ± 1.064.12 ± 1.470.0584.28 ± 1.144.32 ± 1.280.9384.51 ± 1.113.37 ± 1.280.003
V24.67 ± 1.913.91 ± 1.450.0084.33 ± 0.883.86 ± 1.100.2054.53 ± 2.063.62 ± 1.470.002
V35.38 ± 2.094.28 ± 1.870.0275.98 ± 2.205.13 ± 1.560.2515.32 ± 2.015.56 ± 0.990.700
AnkleDorsiflexion angle (◦)V124.30 ± 2.1018.25 ± 1.57<0.00122.34 ± 1.8918.53 ± 0.980.00123.87 ± 1.7118.78 ± 0.92<0.001
V219.61 ± 2.2317.28 ± 1.660.00422.09 ± 3.4117.66 ± 1.520.01723.99 ± 2.2119.74 ± 1.160.001
V321.43 ± 3.5317.34 ± 1.720.01123.04 ± 3.1019.34 ± 1.830.03022.54 ± 2.3519.14 ± 1.620.004
Plantarflexion angle (◦)V142.45 ± 4.6239.90 ± 4.010.04843.59 ± 4.4238.43 ± 4.790.00541.11 ± 5.9138.12 ± 4.970.096
V240.01 ± 8.0239.89 ± 5.680.96746.45 ± 1.6040.94 ± 6.340.04744.88 ± 5.3544.01 ± 4.920.491
V343.53 ± 8.8034.87 ± 6.390.07347.37 ± 7.9545.66 ± 9.130.71446.00 ± 4.7743.67 ± 7.910.552
Dorsiflexion angle velocity (rad/s)V16.12 ± 0.715.71 ± 0.790.3237.37 ± 2.137.24 ± 2.050.7648.23 ± 2.577.10 ± 1.870.038
V26.14 ± 0.726.07 ± 1.430.8978.50 ± 2.716.69 ± 1.340.0488.67 ± 2.948.28 ± 2.790.616
V38.32 ± 2.997.31 ± 1.590.3439.75 ± 3.688.80 ± 2.030.29310.19 ± 3.659.26 ± 3.330.530
Plantarflexion angle velocity (rad/s)V18.52 ± 0.527.80 ± 0.660.1278.69 ± 0.657.38 ± 0.510.0038.14 ± 0.807.45 ± 0.840.138
V210.07 ± 1.539.15 ± 0.700.13411.26 ± 1.108.77 ± 0.620.00111.41 ± 1.139.65 ± 0.46<0.001
V313.20 ± 2.6411.49 ± 2.160.17214.73 ± 2.3713.55 ± 1.820.32115.29 ± 3.2412.73 ± 1.970.124
MTPDorsiflexion angle (◦)V115.03 ± 3.8215.90 ± 3.890.61414.91 ± 2.0517.20 ± 4.220.05911.56 ± 1.6713.40 ± 1.330.046
V213.63 ± 4.1520.11 ± 1.88<0.00118.10 ± 3.3019.45 ± 2.180.42013.10 ± 0.8618.60 ± 2.55<0.001
V316.95 ± 3.6825.28 ± 2.29<0.00115.85 ± 6.1024.57 ± 4.200.02013.67 ± 4.2222.17 ± 1.89<0.001
Plantarflexion angle (◦)V114.05 ± 4.0915.65 ± 4.420.53413.29 ± 3.3712.60 ± 3.470.6218.86 ± 3.678.88 ± 2.340.988
V214.84 ± 5.3713.08 ± 6.000.60812.82 ± 2.5315.30 ± 2.760.0119.89 ± 1.8412.72 ± 2.700.050
V317.71 ± 5.4019.89 ± 7.040.45413.53 ± 7.4020.37 ± 3.020.0139.77 ± 1.6414.83 ± 5.490.020
Dorsiflexion angle velocity (rad/s)V13.53 ± 0.723.53 ± 0.490.9852.81 ± 0.733.78 ± 2.460.2052.83 ± 0.683.35 ± 0.630.214
V23.53 ± 0.723.53 ± 0.490.9852.81 ± 0.733.78 ± 2.460.2052.83 ± 0.683.35 ± 0.630.214
V35.18 ± 2.086.05 ± 1.030.2265.33 ± 1.907.99 ± 1.550.0043.82 ± 1.166.96 ± 2.530.001
Plantarflexion angle velocity (rad/s)V16.66 ± 1.808.44 ± 1.580.0398.00 ± 2.767.13 ± 1.680.3604.84 ± 2.614.23 ± 1.320.618
V27.57 ± 1.018.14 ± 3.590.6629.50 ± 2.428.83 ± 1.990.5876.59 ± 0.997.65 ± 2.050.155
V312.84 ± 4.6215.06 ± 3.570.23611.31 ± 4.7915.73 ± 3.140.0648.29 ± 1.6411.16 ± 2.820.022
L = left lower limb, R = right lower limb. Data in bold in the table indicate statistical significance, p < 0.05. Value: mean ± SD.
Table 3. The influence of running velocity and LBS on SA.
Table 3. The influence of running velocity and LBS on SA.
Joint KinematicsSymmetry Angle [%]VSV·S
S1S2S3pFη2pFη2pFη2
HipExtensionV14.92 ± 2.755.22 ± 3.855.27 ± 2.400.0772.6500.0610.9200.0840.0020.2151.4800.068
V22.49 ± 1.404.90 ± 3.352.58 ± 1.91
V34.93 ± 4.163.03 ± 1.944.42 ± 3.18
Extension angle velocityV15.35 ± 3.441.95 ± 1.062.50 ± 2.020.7020.3550.0090.1192.1890.0510.0782.1810.097
V22.04 ± 1.773.09 ± 2.553.32 ± 2.32
V34.36 ± 3.912.32 ± 2.433.43 ± 2.54
KneeFlexionV13.27 ± 2.42 c1.91 ± 0.99 c4.95 ± 2.74 c0.0017.7530.1610.0672.7980.0650.8520.3370.016
V23.25 ± 1.98 c4.02 ± 2.61 c6.48 ± 3.54 c
V36.31 ± 6.74 ab8.14 ± 5.19 ab9.26 ± 7.91 ab
ExtensionV112.01 ± 9.346.69 ± 3.3113.80 ± 9.100.5720.5620.0140.8920.1140.0030.4490.9340.044
V29.54 ± 8.9312.58 ± 10.837.48 ± 6.75
V312.93 ± 1.1811.88 ± 11.7513.05 ± 12.54
Flexion angle velocityV125.18 ± 2.67 bc28.23 ± 2.11 bc26.40 ± 2.92 bc<0.00172.8450.6430.9780.0220.0010.6520.6170.030
V25.74 ± 5.20 ac4.21 ± 3.73 ac6.66 ± 3.03 ac
V312.73 ± 10.76 ab12.04 ± 7.63 ab10.86 ± 10.55 ab
Extension angle velocityV18.85 ± 8.777.26 ± 6.0711.59 ± 8.010.0782.6310.0610.8880.1190.0030.6940.5580.027
V27.00 ± 3.577.36 ± 7.577.47 ± 5.92
V310.43 ± 5.9513.26 ± 9.8710.91 ± 9.14
AnkleDorsiflexionV19.07 ± 3.996.02 ± 3.027.49 ± 3.070.3701.0080.0240.7310.3140.0080.1421.7720.080
V24.39 ± 2.447.71 ± 5.636.05 ± 3.64
V37.73 ± 4.106.39 ± 5.155.07 ± 3.83
PlantarflexionV12.58 ± 1.87 c4.35 ± 3.16 c4.24 ± 2.44 c0.0017.7550.1610.1302.0950.0490.1841.5940.073
V25.24 ± 6.58 c5.14 ± 4.71 c1.95 ± 2.08 c
V311.58 ± 5.95 ab8.20 ± 5.96 ab5.47 ± 7.22 ab
Dorsiflexion angle velocityV15.45 ± 3.483.54 ± 2.925.61 ± 3.930.1531.9200.0450.3251.1410.0270.4590.9150.043
V26.79 ± 4.327.37 ± 7.115.85 ± 3.90
V310.70 ± 5.976.80 ± 4.495.47 ± 8.78
Plantarflexion angle velocityV13.62 ± 2.11 c5.32 ± 3.40 c4.92 ± 2.77 c0.0433.2790.0750.6700.4030.0100.6250.6550.031
V25.34 ± 3.717.73 ± 4.305.17 ± 2.19
V39.00 ± 6.28 a7.31 ± 4.56 a7.06 ± 8.08 a
MTPDorsiflexionV17.32 ± 7.34 bc4.25 ± 4.20 bc6.10 ± 4.67 bc<0.00116.7480.2930.6260.4720.0120.1231.8720.085
V213.09 ± 7.61 ac6.94 ± 4.32 ac10.65 ± 3.98 ac
V313.00 ± 5.83 ab17.67 ± 7.83 ab15.84 ± 7.42 ab
PlantarflexionV111.41 ± 8.897.50 ± 5.0814.78 ± 9.210.3031.2110.0290.2171.5560.0370.0332.7570.120
V219.47 ± 14.11 f7.10 ± 3.35 f9.57 ± 6.91
V312.65 ± 7.7416.57 ± 11.0315.09 ± 7.67
Dorsiflexion angle velocityV17.48 ± 6.07 ce9.36 ± 8.60 c9.85 ± 7.54 cd0.0194.1640.0930.0144.5160.1000.4600.9130.043
V27.00 ± 5.03 e6.35 ± 6.5615.50 ± 5.77 d
V311.11 ± 8.16 ae14.04 ± 8.92 a17.40 ± 7.99 ad
Plantarflexion angle velocityV19.92 ± 5.768.30 ± 7.78 g18.77 ± 12.51 gh0.6900.3730.0090.7870.2400.0060.0382.6620.116
V214.69 ± 12.9510.44 ± 6.368.05 ± 5.42 h
V311.35 ± 9.9015.58 ± 8.0512.40 ± 6.88
V = running velocity, S = shoes, V·S = the interaction between running velocity and shoes. a = significantly different with V1, b = significantly different with V2, c = significantly different with V3, d = significantly different with S1, e = significantly different with S3, f = have a significant difference between V2/S1 and V2/S2, g = a significant difference between V1/S2 and V1/S3, h = a significant difference between S3/V1 and S3/V2. Data in bold in the table indicate statistical significance, p < 0.05. Value: mean ± SD.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, Q.; Chen, H.; Song, Y.; Alla, N.; Fekete, G.; Li, J.; Gu, Y. Running Velocity and Longitudinal Bending Stiffness Influence the Asymmetry of Kinematic Variables of the Lower Limb Joints. Bioengineering 2022, 9, 607. https://doi.org/10.3390/bioengineering9110607

AMA Style

Liu Q, Chen H, Song Y, Alla N, Fekete G, Li J, Gu Y. Running Velocity and Longitudinal Bending Stiffness Influence the Asymmetry of Kinematic Variables of the Lower Limb Joints. Bioengineering. 2022; 9(11):607. https://doi.org/10.3390/bioengineering9110607

Chicago/Turabian Style

Liu, Qian, Hairong Chen, Yang Song, Nykytiuk Alla, Gusztáv Fekete, Jianpeng Li, and Yaodong Gu. 2022. "Running Velocity and Longitudinal Bending Stiffness Influence the Asymmetry of Kinematic Variables of the Lower Limb Joints" Bioengineering 9, no. 11: 607. https://doi.org/10.3390/bioengineering9110607

APA Style

Liu, Q., Chen, H., Song, Y., Alla, N., Fekete, G., Li, J., & Gu, Y. (2022). Running Velocity and Longitudinal Bending Stiffness Influence the Asymmetry of Kinematic Variables of the Lower Limb Joints. Bioengineering, 9(11), 607. https://doi.org/10.3390/bioengineering9110607

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