1. Introduction
Over the years, researchers have been fascinated by the gait and postural control dynamics and its complicated responses to potential perturbations. In fields such as physics, mathematics, and biology, perturbation refers to a small change or disturbance applied to a system or process, which can lead to alterations in its behaviour, properties, or outcomes [
1,
2]. These disturbances can arise from external forces [
3], such as environmental conditions or physical interactions with objects in the surroundings [
4]. Alternatively, they can originate from internal factors, including physiological changes or neural signals within the organism [
5]. These disturbances might manifest intentionally, such as during experimental manipulations [
6], or they could occur accidentally due to unforeseen events during daily life. It is worth mentioning that external perturbations are not limited to mechanical disturbances. Mechanical disturbances address any physical force or movement that disrupts the balance or proper operation of the system [
7]. These may include vibration, impact, pressure, temperature changes, or external forces like wind or water flow. They can also manifest as sensory disturbances, including hearing and vision [
8].
Roeles, et al. [
8] demonstrated that mechanical perturbations like contralateral sway and treadmill belt deceleration perturbation significantly altered the gait patterns in both young and older adults. In contrast, sensory perturbations (visual, i.e., room darkening, and auditory, i.e., unexpected loud noises) had no significant effect on the gait pattern [
8]. Similar conclusions were reached in the studies conducted by Thies, et al. [
9] and Rogers, et al. [
10]. These findings suggest that an analysis of responses to mechanical perturbations is more justified. To date, many perturbation methods have been proposed in the literature. These include placing obstacles underfoot [
11], using slippery or unstable surfaces [
12], moving plates [
13], and unexpected treadmill accelerations and decelerations [
4,
11,
14]. The field of gait perturbation research, particularly involving both one- and dual-belt treadmill setups, is gaining increasing attention and is presently the most commonly used equipment in such studies [
15]. Semaan, et al. [
16] demonstrated that treadmill gait is a reliable approximation of overground walking. Notwithstanding the divergence between treadmill gait perturbations and those observed in the ground gait condition [
17], Rieger, et al. [
18] demonstrated encouraging outcomes regarding motor skills transfer from an artificial context to real life. Furthermore, treadmills offer the advantage of high repeatability in perturbations and walking conditions [
19]. Modern dual-belt treadmills, such as GRAIL or CAREN, facilitate the generation of many perturbation types across various gait phases and planes. This comprehensive capability enables a detailed analysis of individual gait characteristics in response to perturbations.
Research focused on perturbations is strictly related to the adaptive functions of the human neuromuscular system that allow the maintenance of balance and a straight position [
4,
6,
20,
21,
22]. The elderly, particularly vulnerable to falls, are a primary focus due to their diminished ability to adapt to external disturbances, with nearly a quarter of this population at risk [
23,
24,
25]. Understanding the perturbation response is increasingly important as the global population ages. Numerous studies have examined kinematics, kinetics, spatiotemporal parameters, and muscle activity during gait perturbations [
15]. Predicting responses in young, healthy individuals provides a baseline for future research and potential training or rehabilitation protocols. In addition, young and healthy individuals can safely handle more severe perturbations under laboratory conditions than older adults [
26]. Determining the effects of high-intensity perturbations, such as those induced by the Motek GRAIL system, on kinematic, kinetic, and spatiotemporal parameters is crucial for understanding the body’s response limits and gait pattern changes.
Few studies have investigated perturbations during different gait cycle phases, such as Initial Contact (IC) [
27,
28], Mid Stance (MS) [
29,
30], or Pre-Swing (PS) [
31]. However, Błażkiewicz and Hadamus [
4] investigated how external perturbations induced by treadmill belt acceleration and deceleration during these phases affect gait regularity in young adults. Their findings revealed that, during treadmill deceleration in the Pre-Swing phase, the center of mass (CoM) displacement displayed a consistent pattern in the anterior–posterior (AP) and vertical directions. Conversely, the least regular CoM trajectories were observed during treadmill deceleration in Mid Stance in the AP direction. However, no paper has been found that directly compares the effects of perturbations during all three phases (Initial Contact, Mid Stance, and Pre-Swing) on altered lower limb joint angles and torques compared to undisturbed gait. To address this issue, the purpose of this study was to examine the effect of the acceleration of one belt of a treadmill during three different phases of the gait cycle (Initial Contact, Mid Stance, and Pre-Swing) on kinematic and kinetic parameters and relate these changes to unperturbed gait.
2. Materials and Methods
2.1. Participants
The study group consisted of 21 healthy young females with a mean age of 21.38 ± 1.32 years, body weight of 61.38 ± 6.48 kg, and body height of 165.9 ± 4.53 cm. Eligibility criteria included no muscular or neurological disorders, no lower limb injuries within the last six months, engaging in physical activity at least twice weekly, prior experience with treadmill walking, and right leg dominance [
32]. Leg dominance was determined using the kicking test, where participants indicated their preferred lower limb for kicking a ball. In this study, all individuals had the right lower limb dominant. Exclusion criteria encompassed inexperience with treadmill walking, balance issues, or medications affecting the nervous system. All participants provided written informed consent to participate in the study. The study received approval from the University Review Committee (no. SKE01-15/2023) and followed the ethical guidelines and principles of the Declaration of Helsinki.
2.2. Measurement Protocol
The kinematic (joint angles) and kinetic (joint torques and ground reaction forces) parameters of the perturbed and unperturbed gait were measured in the Interactive Gait Real-Time Analysis Laboratory (GRAIL, Motek Medical B.V., Amsterdam, The Netherlands). The GRAIL consisted of a 2200 mm-long and 2 × 500 mm-wide dual-split-belt treadmill (1000 Hz), a motion capture system (Vicon Metrics Ltd., Oxford, UK) with ten Bonita cameras operating at 100 Hz, three video cameras, synchronised virtual reality environments (circular screen with 5 m diameter, 180°, and 3 m height) with three projectors, and a safety harness suspended from the ceiling (
Figure 1). The accompanying D-Flow software (Motek Medical B.V., Amsterdam, The Netherlands) was used for model adjustment, initiating perturbations, and data collection.
When the participants arrived, anthropometric measurements were collected, and 25 reflective markers were placed on their bodies according to Human Body Model 2 (HBM2) guidelines (
Figure 2). Throughout the study, participants wore their own comfortable sports shoes, which they typically use for walking on a treadmill or exercising at the gym.
Subsequently, participants stepped on the treadmill and aligned in a T-pose to achieve an accurate body model fit in D-Flow 3.26 software. The test protocol included three trials. In each trial, the walking speed was constant at 1.2 m/s. Each trial lasted 80 s and contained five perturbations. The perturbations were introduced at the 30th, 40th, 50th, 60th, and 70th seconds of each trial. The magnitude of the perturbations was five on a scale of 1 to 5. This setting resulted in a change in treadmill belt speed by about 0.5~0.6 m per second. Thus, during the perturbation, the velocity of the left lane of the treadmill was 1.7~1.8 m/s (
Figure 3). This setup was in line with that proposed by Sloot, et al. [
33]. The duration of the perturbation was ~0.82 s.
Perturbations were introduced only on the left lane of the treadmill, affecting the left limb, which was the non-dominant limb (as described in
Section 2.1). Using a high perturbation magnitude facilitated a clearer estimation of peak loads, as the body’s response to significant unexpected forces was more pronounced, making it easier to observe and measure resulting joint loads.
The first trial included perturbations generated only during the Initial Contact (IC) phases. The participants then had a break of about 5 min, during which they did not step off the treadmill. This time was used to save data and set up D-Flow software for the next attempt. Another trial consisted of perturbations during the Mid Stance (MS) phase. Again, this trial was followed by a break of about 5 min. Subsequently, the participants followed the last trial, where perturbations were generated in the Pre-Swing (PS) phase (
Figure 4). Therefore, five perturbations were recorded separately for IC, MS, and PS in each of the three trials.
Perturbations were configured in the D-flow software, allowing for the selection of perturbation intensity, event type, treadmill lane, and generation time interval. Since gait phases are defined as intervals rather than specific points in time, generating a perturbation at, for example, the 30th second during the MS phase (which, according to Perry, et al. [
34], falls at 10–30% of the gait cycle) could cause a perturbation to occur at 20% of the gait cycle for one person and at 25% for another.
Perturbation presence was most visible for the vertical component of the ground reaction force (vGRF) (
Figure 5A). Thus, the definition of the gait cycle, both with and without perturbations, was based on this variable. For gait without and with perturbations, the gait cycle was defined from the first contact of the left foot with the ground to its re-contact with the ground. Additionally, a visual verification was performed using motion data from OpenSim to ensure that the gait cycle was accurately identified, especially during perturbations. For each participant, the undisturbed gait cycle from the second trial, which occurred around the 20th second of gait, was used for analysis, as participants were more accustomed to walking on a treadmill (
Figure 5B). For perturbed gait, only the first perturbation gait cycle in each trial was considered, typically occurring around the 30th second of gait (
Figure 5C).
Therefore, for analysis, there were 21 gait cycles without perturbations for the left lower limb, 21 gait cycles with perturbations in the IC phase, 21 gait cycles with perturbations in the MS phase, and 21 gait cycles with perturbations in the PS phase.
2.3. Data Processing
Data obtained from D-Flow 3.26 software were saved as *.mox files, and then imported into Matlab 2021a (MathWorks, Natick, MA, USA) using a toolbox developed by Feldhege, et al. [
35]. Subsequently, they were transformed into OpenSim 4.2 input files: *.trc, which contained positions of markers placed on a subject at different times during a motion capture trial, and *.mot, which included ground reaction force data [
36]. The gait2354_model, a three-dimensional, 23-degree-of-freedom computer model of the human musculoskeletal system with 54 musculotendon actuators representing 27 muscles in the lower extremities and torso, was adjusted to participants’ anthropometry using the Scale Tool. Joint angles were determined using the Inverse Kinematics Tool (IK) based on marker positions during motion. Joint torques were obtained using the Inverse Dynamics Tool (ID) and ground reaction forces (GRFs).
The first perturbed gait cycle for the left leg for each of the three trials, for each person, was determined in the OpenSim software, considering the visualisation of inverse kinematics and GRF. As previously described, it was around the 30th second of gait. The same process was applied to gait without perturbations, but only for the second trial. The gait cycle, which occurred around the 20th second of the trial, was included for each individual. The gait cycle length for each perturbed trial and the gait without perturbations was estimated using the Euclidean distance equation based on the position of the left heel marker at the heel strike, as follows:
where
x1, y1, z1 denote the position of the heel marker at the beginning of the gait cycle, and
x2, y2, z2 denote the position of the heel marker at the end of the gait cycle, which was taken into consideration. For subsequent statistical analysis, the stride length, and extreme values (minimal and maximal) of kinematic and kinetic parameters were extracted for gait cycles both with and without perturbations.
2.4. Statistical Analysis
Statistical analysis was performed using PQStat 2021 software v.1.8.2.238 (PQStat Software, Poznań, Poland). Normal distribution of the data was assessed using the Shapiro–Wilk test, indicating normal distributions for extreme values of joint angles and stride length in both unperturbed and perturbed gait conditions (IC, MS, and PS). Extreme values of joint torques and stride length during MS perturbation did not follow a normal distribution.
ANOVA with Tukey’s HSD post hoc test was used to identify statistically significant differences in the extreme values of joint angles achieved in gait cycles with perturbations in IC, MS, and PS and gait cycles without perturbation (N). Friedman’s ANOVA test with Bonferroni’s post hoc test was employed to detect statistically significant differences in the extreme values of joint torques and stride length achieved in the gait cycle with perturbations in IC, MS, and PS, and the gait cycle without perturbation (N).
ANOVA and post hoc tests, with a non-parametric approach using statistical parametric mapping (SPM) [
37] was employed on continuous curves, specifically focusing on angles and muscle torques at the lower limb joints. The SPM1D package “(
https://www.spm1d.org) (accessed on 14 February 2024)” facilitated these analyses. For each SPM ANOVA, a statistical parametric map SPM{F} was created by calculating the univariate F-Statistic at each point of the gait curve [
37]. If SPM{F} crossed the threshold line corresponding to a confidence level equal to 0.95, post hoc SPM{t} maps were calculated for comparing each pair of independent variables. When the SPM{t} map crossed the critical threshold, a significant difference was found between the examined pair of motion tasks. This manuscript only reported comparisons of curves that contained perturbations in the IC, MS, and PS phases, relative to curves recorded during unperturbed gait. All analyses were performed in Matlab R2021a. A similar approach was used in the paper [
38].
4. Discussion
An understanding of the influence of external perturbations on gait dynamics can assist clinicians and researchers in the design of personalised rehabilitation programmes aimed at improving balance, stability, and overall gait function. These programmes can reduce the risk of falls and enhance mobility and quality of life for individuals with neuromuscular impairments or age-related declines in gait performance. This study aimed to provide support for this field by identifying the most vulnerable phase of the gait cycle to perturbation and the highest overloads that occur during the perturbation. This study employed a perturbation protocol using the Grail system (Motek Medical BV, Amsterdam, The Netherlands), which allowed for the precise timing and intensity control of perturbations. The findings indicate that all perturbations significantly influenced the examined parameters, probably due to the application of high perturbation force [
4]. A consistent significant alteration observed across all perturbations was the reduction in full knee extension, as evidenced by the knee joint angle curves (
Figure 6B) and extreme values for the knee extension angle (
Table 2). All this suggests a critical role of the knee joint in adapting to perturbations, a trend similarly noted by Shokouhi, et al. [
39], who examined adjustments in lower limb joint power and work in response to trip-induced perturbations. In contrast to Shokouhi et al.’s [
39] trip-like perturbation, which involved rapidly decelerating the belt under the dominant foot, this study utilised accelerations. Additionally, Błażkiewicz and Hadamus [
4] assessed gait regularity during external perturbations caused by treadmill belt acceleration and deceleration across different phases of the gait cycle. They found that perturbations caused by treadmill belt deceleration contribute to higher irregularities, with the most irregular behaviour observed during the Pre-Swing phase (vertical direction) among those caused by acceleration. In this paper, an analysis of the kinetic and kinematic values of the perturbed gait indicated that perturbations induced during the Initial Contact (IC) phase were the most intense. Perturbations in the Pre-Swing (PS) phase also had a significant impact on joint angular values, but they did not result in substantial changes in torque values, thus making them less likely to cause overloads. Perturbations in the Mid Stance (MS) phase were determined to have the least significance. These findings will be discussed sequentially, starting with the most impactful perturbations.
4.1. Perturbations in the Initial Contact Phase
The Initial Contact (IC) phase, marking the beginning of the gait cycle, plays a fundamental role in locomotion [
34]. It requires good co-ordination, balance, and effective shock absorption capabilities. The perturbation induced in the IC phase affects the entire gait cycle, which is probably the reason why it accounts for the most significant changes observed in the examined parameters.
The most noteworthy discovery revolved around the alteration in knee joint torques, as the areas showing significant changes covered 55.86% of the gait cycle (
Figure 6E). Specifically, the peak knee flexion torque reached 1.15 ± 0.29 Nm/kg, marking a staggering 248.48% increase compared to unperturbed gait. It was the most substantial kinetic shift across all perturbations and joints. Additionally, there was a noteworthy rise in peak knee extension torque −0.66 ± 0.21 Nm/kg in the presence of perturbation, compared to the −0.31 ± 0.03 Nm/kg recorded during unperturbed gait. Moreover, the perturbation during the Initial Contact (IC) phase induced substantial alterations in knee kinematics, affecting up to 72.34% of the gait cycle. However, there was only a minor increase of 12.26% in extreme flexion values.
Perturbations significantly altered the ankle joint angle curve values in 77.38% of the gait cycle area compared to a typical gait pattern. These alterations were most notable in the extremes of plantar flexion, exhibiting a substantial increase from −23.91 ± 9.42 deg to −4.96 ± 2.25 deg. Despite this pronounced shift, the torque exerted by the plantar flexors experienced a notable decrease of 27.72% compared to the torque observed during unperturbed gait.
Perturbations also significantly affected the kinematic parameters associated with the hip joint. These alterations impacted 61% of the gait cycle, particularly visible during the late stance and all swing phases. Notably, perturbation during the Initial Contact phase induced the most substantial changes in the extreme angular values of the hip joint, with hip flexion increasing by 29.71% and hip extension by 36.62%. In terms of torque, there was a noteworthy 78.69% increase in hip flexion torque.
It is worth noting that reviews by Taylor, et al. [
15], as well as Ferreira, et al. [
11], have underscored that heel strike perturbations are prevalent. This study corroborates this finding, demonstrating that perturbations during the IC phase significantly impacted all examined parameters except torques related to dorsiflexion and hip extension.
Notably, knee flexion torque exhibited a nearly threefold increase, emphasising the pivotal role of the knee joint and quadriceps muscle in responding to such perturbations. This observation is particularly relevant given that even short periods of immobilisation can lead to rapid muscle weakness [
40], potentially heightening the risk of falls triggered by perturbations.
4.2. Perturbations in the Pre-Swing Phase
The Pre-Swing phase occurs at 50–60% of the gait cycle, immediately preceding the swing phase. Perturbation induced in this phase had the most significant impact on increasing step length, artificially accelerating the swing phase, and forcing participants to actively decelerate the limb. However, concerns arise as the gait cycle for perturbations in the PS phase was disrupted even before the perturbation occurred.
This disruption does not seem attributable to mistimed perturbation induction, as Golyski, et al. [
29] reported an average delay of 56 ms in perturbation induction by the GRAIL and D-Flow systems, nor does it appear to stem from changes in percentage intervals due to gait cycle shortening or lengthening. A more plausible explanation seems to be that the perturbation in the PS phase was induced last, prompting participants to adopt a different gait pattern in anticipation, as observed in studies by Swart, et al. [
41]. The hypothesis that participants adopted a different gait strategy is also supported by the kinematic curves of the knee joint, showing a lack of full extension compared to unperturbed gait. This relationship is particularly visible in the area of [7–35.4]% of the gait cycle (
Figure 6E). Significant changes also occur in the hip joint [16–45]% and in the ankle joint [9–17]%, and [43–48]% in areas before the disturbance. The perturbation in the PS phase significantly affected all extreme kinematic values, except for the hip extension angle. The greatest changes compared to unperturbed gait were observed for plantar flexion −20.93 ± 8.22 vs. −4.96 ± 2.25 deg and knee extension 3.4 ± 3.81 vs. −3.28 ± 1.95 deg.
As for the analysis of kinetic parameters, changes in the gait pattern appeared in the [5–36]% of the gait cycle for the ankle torque and within [13–21.18]%, and [36.31–52.14]% of the gait cycle for the knee joint torque curves. However, changes in the kinetic parameters of the hip joint are directly related to the moment of perturbation, appearing, respectively, in the areas of [55–60]% and [72.41–95]% of the gait cycle. The most significant alteration in extreme values occurred at the hip joint, with the hip extension torque increasing by 73.02%. This represented the most substantial change in this parameter among all examined perturbations. These results are consistent with Vlutters, et al. [
42], who also identified the hip joint as a primary control for this type of disturbance.
Regarding the perturbation during the Pre-Swing (PS) phase, the most notable changes were noted for the hip extension torque. However, the existence of significant alterations preceding the induced perturbation is intriguing and challenging to explain. These pre-perturbation changes are possibly attributable to anticipatory muscle activation in preparation for the impending perturbation. But this issue requires further research.
4.3. Perturbations in the Mid Stance Phase
The perturbations induced during the Mid Stance phase had minimal impact on the parameters studied compared to the previously mentioned perturbations. Notably, only this perturbation did not alter stride length compared to unperturbed gait. However, this finding contrasts with the Madehkhaksar, et al. [
43] report, which observed a significantly shorter step length under similar conditions.
The most common strategy adopted in response to this perturbation involved rising onto the toes (
Figure 3). This is evident in the notable deviations observed in ankle plantar flexion occurring within the intervals [23.47–51.17; 57.81–59.22]% of the gait cycle. Moreover, extreme torque values for plantar flexion increased by 118.18% compared to those recorded during undisturbed gait, marking the most significant change in this parameter of all the perturbations. Furthermore, there was a significant increase in extreme torque values for knee flexion—93.94%—compared to unaffected gait. Aligned with the findings presented in this study are the results of Taborri, et al. [
38]. Taborri, et al. [
38] demonstrated a noticeable decrease in ankle dorsiflexion and plantar flexion when locomotion was perturbed. However, it is worth noting that the authors analysed ranges of motion rather than angles.
Before perturbation, as for perturbations in the Pre-Swing phase, some areas showed significant differences from the unperturbed gait cycle. These included changes in [0–9]% and [5–60]% of the gait cycle for ankle and knee kinematics, respectively. For kinetic parameters, the change regions were [2–4.22; 6.63–20]% and [6–22.1]% of the gait cycle for the knee and hip joint, respectively (
Figure 6). It seems that the explanation for this phenomenon may be similar to that for perturbation in the Pre-Swing phase, and it relates to the premature activation of muscles in subconscious anticipation of a perturbation, as this type of perturbation was studied as the second one in the experiment. Santuz, et al. [
44], by applying the muscle synergy theory, demonstrated that humans are able to modify their motor control strategies, especially in terms of timing, when walking in unsteady conditions.
4.4. Study Limitations
This study has some limitations. First, due to the technical limitations of the treadmill, the perturbations did not occur at the same time of the gait cycle, only in the same interval, which could lead to different muscle activations and, thus, different behaviours. However, as presented, the results were pretty consistent within the perturbations. Another limitation concerns the analysis of the perturbed limb only. In subsequent studies, it would also be worthwhile to consider the recovery limb. Given that there was a delay in the triggering of the perturbation due to limitations in the configuration of the treadmill device, the frame in which the perturbation was presented could not be detected, making point analysis impossible and forcing an interval analysis (SPM). It is not known whether the extreme values were after the perturbation or at the time of its occurrence. Thus, some cautions in interpreting the results should be taken into account, considering that, depending on the moment of the perturbation within the gait cycle, the response is different.
4.5. Future Research
In future research, it would be beneficial to investigate the response to perturbations of both young and elderly individuals through comprehensive assessments of kinematic and kinetic parameters. However, challenges may arise due to differences in treadmill walking experience and the need to adjust the walking speed accordingly. Another avenue for future research could be to compare the impact of perturbations on the Motek GRAIL system to that of more affordable tools, such as a moving platform or other similar devices that are less expensive to maintain. Furthermore, the literature on induced perturbations uses the terms slip, trip, or stumble misleadingly, as authors use the same terms to describe different types of perturbations. Therefore, it would be beneficial for future research to provide a clear and detailed explanation for all induced perturbations and propose a common terminology.
5. Conclusions
The research findings suggest that each type of examined one-belt-acceleration perturbation presents a challenge to balance capabilities. However, the perturbation in the IC phase was associated with the highest loads, particularly in terms of the extreme values of the knee flexion torque, which were three times higher compared to normal gait.
Perturbations in the PS phase primarily influenced the hip joint, with the most significant change observed in the extension torque. The perturbations that occurred in the MS phase mainly affected the ankle joint, with a notable influence on the plantar flexion torque.