Gait Analysis & Ground Reaction Force (GRF) Calculator
Analyze human bipedal locomotion kinetics. Model vertical and anteroposterior ground reaction forces, loading peaks, valley unloading, stance/swing durations, vertical loading rates, and mechanical impulse.
Gait Analysis & Ground Reaction Force (GRF) Calculator: Free Gait Analysis & Ground Reaction Force (GRF) Calculator tool: calculate biomechanical gait parameters, vertical ground reaction force (GRF) peaks.
Subject & Locomotion Parameters
Kinetics & Biomechanical Metrics
Biomechanical Mechanics of Human Gait and Ground Reaction Forces
Terrestrial human locomotion is an exquisitely coordinated neuromusculoskeletal feat that balances forward propulsion against gravitational and energetic stability. Whenever a human foot makes contact with a force platform or the earth, the ground exerts an equal and opposite force vector upon the body, termed the Ground Reaction Force ($mathbf{F}_{ ext{GRF}}$) according to Newton's Third Law of Motion.
In three-dimensional motion capture and gait laboratories, the instantaneous ground reaction force is decomposed into three orthogonal Cartesian vector components:
where $F_z$ is the vertical component (resisting gravity and accelerating the body vertically), $F_y$ is the anteroposterior shear component (governing deceleration/braking and forward acceleration/propulsion), and $F_x$ is the mediolateral shear component (regulating lateral balance and turning).
The Inverted Pendulum Model of Walking
In level walking, the body's center of mass (CoM) vaults passively over a relatively straight stance leg, acting analogously to an inverted pendulum. The trajectory of the vertical force $F_z(t)$ exhibits a characteristic "M-shaped" bimodal curve:
- Initial Contact & Loading Response Peak ($F_{z1}$): Occurring at approximately 15% of the stance phase, this peak reaches $1.10$ to $1.25 imes ext{BW}$. As the heel strikes the ground, the downward momentum of the center of mass is arrested, requiring an upward acceleration ($F_z > mg$).
- Mid-Stance Unloading Valley ($F_{z,min}$): At approximately 50% of stance, the CoM passes over the upright support limb, reaching its maximum gravitational potential energy height and minimum forward kinetic velocity. The downward centrifugal acceleration reduces the vertical force below body weight ($0.65$ to $0.80 imes ext{BW}$).
- Terminal Stance Push-Off Peak ($F_{z2}$): At approximately 75% to 80% of stance, the plantarflexor muscle group (soleus, gastrocnemius) contracts forcefully to accelerate the body upward and forward into pre-swing, generating the second force peak ($1.15$ to $1.30 imes ext{BW}$).
- Toe-Off (Preswing, 100% Stance): The foot terminates ground contact as vertical force drops smoothly to zero.
The Spring-Mass (SLIP) Model of Running
When transitioning from walking to running (typically occurring at the Froude number transition speed $v approx sqrt{g cdot L_{ ext{leg}}} approx 2.0 ext{ m/s}$), the inverted pendulum mechanism breaks down. Running behaves as a Spring-Loaded Inverted Pendulum (SLIP), where kinetic energy and gravitational potential energy fluctuate in-phase. The stance phase compresses compliant muscle-tendon springs (primarily the Achilles tendon and quadriceps tendon), storing elastic strain energy during early stance and releasing it during push-off.
Key running kinetic characteristics include:
- Absence of Mid-Stance Valley: The vertical force curve becomes unimodal (a single large bell-shaped active peak reaching $2.2$ to $3.0 imes ext{BW}$).
- Impact Transient Spike (Rearfoot Strike): When striking on the calcaneus (heel), a sharp impact spike occurs within the first 25–40 milliseconds of contact ($1.5$ to $2.0 imes ext{BW}$), characterized by an exceptionally steep loading rate.
- Forefoot Attenuation: In barefoot or forefoot runners, eccentric contraction of the calf complex and controlled ankle dorsiflexion absorbs impact energy, smoothing the initial spike into a clean bell curve.
Vertical Loading Rate and Overuse Injury Biomechanics
The rate at which force rises upon impact is quantified by the Vertical Average Loading Rate (VALR) and Vertical Instantaneous Loading Rate (VILR):
Clinical research (e.g., Davis et al., Milner et al.) indicates that runners exhibiting VALR values exceeding $75 ext{ BW/s}$ possess significantly higher historical incidences of tibial stress fractures, Achilles tendinopathy, and plantar fasciitis. Interventions such as increasing cadence by 5% to 10% or transitioning to a midfoot strike pattern reduce loading rates by 20% to 40%.
Spatio-Temporal Relationships
- Cadence ($C$): Number of steps per minute. Step time is $T_{ ext{step}} = 60 / C$ (seconds).
- Stride Duration: In symmetrical gait, a full stride consists of two steps: $T_{ ext{stride}} = 2 imes T_{ ext{step}} = 120 / C$.
- Stance and Swing Durations: $$T_{ ext{stance}} = T_{ ext{stride}} imes left( rac{ ext{Duty Factor}}{100} ight), quad T_{ ext{swing}} = T_{ ext{stride}} - T_{ ext{stance}}$$
- Step and Stride Length: $$L_{ ext{step}} = rac{v}{C / 60} = rac{v cdot 60}{C} quad [ ext{meters}], quad L_{ ext{stride}} = 2 imes L_{ ext{step}}$$
- Mechanical Vertical Impulse ($J_z$): The time integral of vertical force over the stance duration: $$J_z = int_0^{T_{ ext{stance}}} F_z(t) , dt quad [ ext{N}cdot ext{s}]$$ In steady-state constant-velocity locomotion, the average vertical force over a full stride must equal body weight: $ar{F}_z imes T_{ ext{stance}} = m g imes T_{ ext{stride}}$.
Step-by-Step Worked Numerical Examples
Example 1: Biomechanical Analysis of Adult Walking
Given: A 70.0 kg adult walks at self-selected comfortable speed $v = 1.35 ext{ m/s}$ with a measured cadence of $C = 112 ext{ steps/min}$ and stance duty factor of $62%$.
Step 1: Compute Body Weight (BW):
$$ ext{BW} = m cdot g = 70.0 ext{ kg} imes 9.80665 ext{ m/s}^2 = 686.47 ext{ N}$$Step 2: Temporal Parameters:
$$T_{ ext{step}} = rac{60}{112} = 0.5357 ext{ s} = 535.7 ext{ ms}$$ $$T_{ ext{stride}} = 2 imes 0.5357 = 1.0714 ext{ s} = 1071.4 ext{ ms}$$ $$T_{ ext{stance}} = 1.0714 imes 0.62 = 0.6643 ext{ s} = 664.3 ext{ ms}$$ $$T_{ ext{swing}} = 1.0714 - 0.6643 = 0.4071 ext{ s} = 407.1 ext{ ms}$$Step 3: Step and Stride Length:
$$L_{ ext{step}} = rac{1.35 ext{ m/s}}{112 / 60 ext{ s}^{-1}} = rac{1.35}{1.8667} = 0.7232 ext{ m} approx 72.3 ext{ cm}$$ $$L_{ ext{stride}} = 2 imes 0.7232 = 1.446 ext{ m}$$Step 4: Peak Force Estimates:
$$F_{z1} approx 1.20 imes ext{BW} = 1.20 imes 686.47 = 823.8 ext{ N}$$ $$F_{z,min} approx 0.75 imes ext{BW} = 0.75 imes 686.47 = 514.9 ext{ N}$$ $$F_{z2} approx 1.25 imes ext{BW} = 1.25 imes 686.47 = 858.1 ext{ N}$$ Peak braking force: $F_{y, ext{brake}} approx -0.20 imes 686.47 = -137.3 ext{ N}$. Peak propulsive force: $F_{y, ext{prop}} approx +0.22 imes 686.47 = +151.0 ext{ N}$.Example 2: Running Biomechanics Comparison (Walking vs. Marathon Pace)
Now consider the same 70 kg athlete running at $v = 3.5 ext{ m/s}$ ($5:00 ext{ min/km}$ pace), cadence $C = 170 ext{ steps/min}$, duty factor $= 38%$.
1. Stride time: $T_{ ext{stride}} = 120 / 170 = 0.7059 ext{ s} = 705.9 ext{ ms}$.
2. Contact stance time: $T_{ ext{stance}} = 0.7059 imes 0.38 = 0.2682 ext{ s} = 268.2 ext{ ms}$ (less than half the walking contact time!).
3. Aerial flight time: Because stance is $38%$, the flight phase occupies $100% - 2 imes 38% = 24%$ of the stride ($169.4 ext{ ms}$).
4. Active peak force:
$$F_{z,max} approx 2.50 imes ext{BW} = 2.50 imes 686.47 = 1,716.2 ext{ N}$$
5. Vertical Loading Rate (VALR):
If the impact transient of $1.8 imes ext{BW}$ ($1,235.6 ext{ N}$) rises in $30 ext{ ms}$:
$$ ext{VALR} approx rac{1,235.6 ext{ N} imes 0.6}{0.018 ext{ s}} = rac{741.4}{0.018} = 41,189 ext{ N/s} approx 41.2 ext{ kN/s} approx 60.0 ext{ BW/s}$$
Locomotion Kinetic Comparison Across Modalities
| Locomotion Activity | Stance Duty Factor | Peak Vertical $F_z$ | Contact Time | Typical VALR | Mechanical Model |
|---|---|---|---|---|---|
| Slow Walking (0.8 m/s) | 65% – 68% | 1.05 – 1.12 BW | 750 – 900 ms | 5 – 8 BW/s | Rigid Inverted Pendulum |
| Comfortable Walk (1.35 m/s) | 60% – 63% | 1.15 – 1.25 BW | 600 – 700 ms | 10 – 15 BW/s | Compliant Inverted Pendulum |
| Jogging (2.8 m/s) | 40% – 45% | 2.00 – 2.40 BW | 280 – 330 ms | 45 – 70 BW/s | Spring-Mass (SLIP) |
| Distance Running (4.0 m/s) | 32% – 38% | 2.50 – 2.85 BW | 220 – 260 ms | 65 – 95 BW/s | Spring-Mass (SLIP) |
| Maximal Sprinting (9.0 m/s) | 18% – 22% | 3.50 – 4.50 BW | 90 – 120 ms | > 150 BW/s | High-Stiffness Ankle Piston |
Biomechanical FAQs
What is the Center of Pressure (CoP) and how is it calculated from a force plate?
The Center of Pressure (CoP) is the point of application of the net ground reaction force vector on the foot-ground interface. Multi-axis force plates contain four load cells positioned beneath each corner. By measuring vertical forces ($F_{z1..4}$) and moments about the axes ($M_x, M_y$), the CoP coordinates are computed as $x_{ ext{CoP}} = -M_y / F_z$ and $y_{ ext{CoP}} = M_x / F_z$. Tracking the CoP trajectory reveals excessive foot pronation, supination, or balance deficits.
Why does double support disappear as speed increases?
Double support (the phase where both feet simultaneously contact the ground) occupies approximately 20% to 24% of the walking cycle. As velocity increases, stance duration decreases quadratically while swing duration remains relatively invariant. When stance duration drops below 50% of the stride cycle, both feet cannot be on the ground at the same time, marking the fundamental physical threshold where double support vanishes and an aerial flight phase emerges.
How does shoe cushioning alter Ground Reaction Force kinetics?
Contrary to popular belief, thick cushioned running shoes do not substantially reduce the active peak vertical force ($F_{z,max}$), which is governed primarily by total body mass and center-of-mass trajectory. However, viscoelastic midsole foams (such as EVA or PEBA) deform upon contact, increasing the deceleration rise time of the foot. This attenuates the high-frequency impact transient and reduces the vertical loading rate (VALR), smoothing out shockwave transmission into the tibia.
Frequently Asked Questions
The Center of Pressure (CoP) is the point of application of the net ground reaction force vector on the foot-ground interface. Multi-axis force plates contain four load cells positioned beneath each corner. By measuring vertical forces ($F_{z1..4}$) and moments about the axes ($M_x, M_y$), the CoP coordinates are computed as $x_{ ext{CoP}} = -M_y / F_z$ and $y_{ ext{CoP}} = M_x / F_z$. Tracking the CoP trajectory reveals excessive foot pronation, supination, or balance deficits.
Double support (the phase where both feet simultaneously contact the ground) occupies approximately 20% to 24% of the walking cycle. As velocity increases, stance duration decreases quadratically while swing duration remains relatively invariant. When stance duration drops below 50% of the stride cycle, both feet cannot be on the ground at the same time, marking the fundamental physical threshold where double support vanishes and an aerial flight phase emerges.
Contrary to popular belief, thick cushioned running shoes do not substantially reduce the active peak vertical force ($F_{z,max}$), which is governed primarily by total body mass and center-of-mass trajectory. However, viscoelastic midsole foams (such as EVA or PEBA) deform upon contact, increasing the deceleration rise time of the foot. This attenuates the high-frequency impact transient and reduces the vertical loading rate (VALR), smoothing out shockwave transmission into the tibia.