Abstract

People with unilateral transtibial amputations (TTA) have greater risks of bilateral hip osteoarthritis, related to asymmetric biomechanics compared to people without TTA. Running is beneficial for physical health and is gaining popularity. However, people with TTA may not have access to running-specific prostheses (RSPs), which are designed for running, and may instead run using their daily-use prosthesis (DUP). Differences in joint loading may result from prosthesis choice; thus, it is important to characterize changes in peak and impulsive hip joint contact loading during running. Six people with and without TTA ran at 3.5 m/s while ground reaction forces, kinematics, and electromyography were collected. People with TTA ran using their own RSP and DUP. Musculoskeletal models incorporating prosthesis type of each individual were used to quantify individual muscle forces and hip joint contact forces (HJCFs) during running. People using RSPs had smaller bilateral peak hip joint contact forces compared to when wearing DUPs during stance and swing, and a smaller impulse over the entire gait cycle. Greater amputated leg peak hip joint contact forces for people wearing DUPs compared to RSPs occurred with greater forces from the ipsilateral gluteus maximus during stance. People with TTA also had greater bilateral peak hip joint contact forces during swing compared to people without TTA, which occurred with greater peak gluteus medius forces. Running with more compliant RSPs may be beneficial for long-term joint health by reducing peak and impulsive hip loading compared to DUPs.

Introduction

People with a transtibial amputation (TTA) are increasingly interested and participating in physical activities like running [1], which has numerous emotional, social, and physical health benefits [13]. However, limited access to appropriate prostheses is a barrier for many people with TTA who wish to run. Running-specific prostheses (RSPs) are designed for running and sprinting tasks with large capacity for energy storage and return. However, RSPs can be difficult to obtain as they are expensive and many insurance providers in the United States do not cover the cost of the devices. Therefore, if a person with TTA desires to run, he or she may do so using their daily-use prosthesis (DUP), which are stiffer devices that are not designed for highly dynamic activities like RSPs.

People with TTA have greater risks and rates of bilateral hip osteoarthritis (OA) compared to the age-matched population without TTA (14% versus 1%) [47], but there is conflicting evidence as to whether the intact or amputated side hip of people with TTA is at greater risk [5,7,8]. OA is a degenerative and painful disease characterized by the breakdown of articular cartilage [9]. The causality of OA in people with TTA is not fully understood, but asymmetric, large, and repeated joint contact forces are risk factors for the development of the disease in people with and without TTA [4,1013]. For people with TTA, asymmetric movement biomechanics and muscle strength further increase OA risk factors [4,14]. Greater device energy storage and return in RSPs reduce hip joint kinetic asymmetry compared to DUPs [15]. Specifically, using compliant RSPs decreases amputated leg positive hip work (i.e., the time integral of positive hip joint power) compared to DUPs during running. Smaller positive hip work from the amputated leg reduces hip joint work asymmetry commonly observed in people with TTA. RSPs are generally more compliant compared to DUPs (RSP: 18–42 kN/m, DUP: 49–72 kN/m) [16,17], although these can vary greatly depending on manufacturer, design, and prescription. Repeated, large joint contact forces are also contributing factors to high rates of hip OA, and are likely affected by the presence of an amputation and device type, but have not yet been quantified in people with TTA during running [4,11,12,18]. Joint contact impulse, which is the time integral of joint contact force, is descriptive of cumulative load. Quantifying impulse over the gait cycle can reveal additional differences in loading than peak force alone [19].

During walking, the gluteus maximus, gluteus medius, vasti, and hamstrings are substantial contributors to hip joint contact forces (HJCFs) [20], in addition to the rectus femoris and iliopsoas providing smaller contributions. People with TTA have reduced amputated leg quadriceps strength and altered activity in the amputated leg gluteus medius and quadriceps during locomotion compared to people without TTA [2123]. Thus, HJCFs are likely also different between people with and without TTA. Furthermore, characterizing the effect of device type on HJCFs is important for understanding potential effects of prosthesis choice on OA risk factors. Investigating HJCFs becomes especially important for people with TTA during activities with repeated and large magnitude joint loading, such as running. Joint contact forces are nearly impossible to measure directly, and joint intersegmental force computations do not include contributions from muscle forces [24,25]. However, muscle forces substantially contribute to joint contact forces [26]. Alternatively, musculoskeletal modeling and simulation can be used to quantify internal muscle force and HJCF results.

The purpose of this study was to determine how the presence of a unilateral transtibial amputation and the type of prosthesis (DUP versus RSP) affect peak HJCFs and impulses during running. We hypothesized that people with TTA would have greater peak HJCFs during the stance and swing phases and hip joint contact impulses over the running gait cycle compared to people without TTA due to documented hip joint kinetic compensations in this population [15,27]. We also hypothesized that using an RSP would reduce peak HJCFs and impulses in the amputated leg compared with using a DUP, due to prior evidence that use of an RSP reduces the magnitude of hip compensations (i.e., hip joint work) in the amputated leg. Similarly, given hip compensatory strategies for people with TTA, we expected that the amputated leg would have greater peak HJCF compared to the intact leg. To assist in interpretation of HJCF and impulse results, we also examined muscle forces that largely contribute to HJCFs, including the gluteus maximus, gluteus medius, vasti, hamstrings, rectus femoris, and iliopsoas.

Methods

Six people with unilateral TTA (5 M/1F, 1.78 m ± 0.057 m, 73.4 kg ± 8.14 kg, 30 ± 7 years) and six people without TTA (5 M/1F, 1.76 m ± 0.062 m, 76.7 kg ± 10.5 kg, 28 ± 5 years) (Table 1) provided written informed consent to participate in the experimental protocol approved by the Institutional Review Board. Participants were free from injury at the time of testing.

Table 1

Participant prosthesis type and etiology

ParticipantDaily-use prosthesisRunning-specific prosthesisTime since amputationEtiology
P1Össur Proflex XCFreedom Innovations Catapult1 year 5 monthsTraumatic
P2Ottobock Advantage DP2Össur FlexRun16 years 5 monthsTraumatic
P3Freedom Innovations RenegadeFreedom Innovations Catapult9 yearsTraumatic
P4Fillauer All-ProFreedom Innovations Catapult8 years 6 monthsTraumatic
P5Freedom Innovations SilhouetteEndolite Running Prosthesis2 years 5 monthsTraumatic
P6Freedom Innovations RenegadeÖssur FlexRun8 yearsTraumatic
ParticipantDaily-use prosthesisRunning-specific prosthesisTime since amputationEtiology
P1Össur Proflex XCFreedom Innovations Catapult1 year 5 monthsTraumatic
P2Ottobock Advantage DP2Össur FlexRun16 years 5 monthsTraumatic
P3Freedom Innovations RenegadeFreedom Innovations Catapult9 yearsTraumatic
P4Fillauer All-ProFreedom Innovations Catapult8 years 6 monthsTraumatic
P5Freedom Innovations SilhouetteEndolite Running Prosthesis2 years 5 monthsTraumatic
P6Freedom Innovations RenegadeÖssur FlexRun8 yearsTraumatic

Experimental Protocol.

Participants ran on an instrumented treadmill (Bertec, Inc., Columbus, OH) at 3.5 m/s (7:40 min/mile) while GRFs were collected at 2000 Hz. Participants were instrumented with a set of 56 (people without TTA and people with TTA wearing DUPs) or 63 (people with TTA wearing RSPs) active markers. Markers on the DUP were placed similarly to the intact leg. The RSP had seven individual markers on the prosthesis to characterize deformation of the carbon-fiber device (four physical and three digitized) [15]. Kinematics were collected with an optoelectric motion capture system at 100 Hz (3DInvestigator, Northern Digital, Inc., ON, Canada).

Surface electromyography (EMG) signals were collected at 2000 Hz (Bortec, Inc., Canada) on a set of four muscles including the bilateral vastus lateralis, biceps femoris long head, rectus femoris, and medial gastrocnemius (bilateral for people without TTA and intact leg of people with TTA). EMG electrodes were placed parallel to the direction of the muscle fiber (Ag/AgCl, 22 mm interelectrode distance (Blue Sensor N-00-S, Ambu, Inc., Denmark)).

Model Development.

Musculoskeletal models were developed for each participant by scaling generic models in opensim v3.3 to match each participant's size and mass [2830]. Each model had varying degrees-of-freedom (DOF) and number of Hill-type musculotendon actuators depending on the presence of TTA and the type of prosthesis (Fig. 1). Hill-type musculotendon actuators incorporated activation and deactivation dynamics as a simplified first-order dynamic model and contraction dynamics with force–length–velocity properties [31]. For each model, lumbar motion was modeled with three rotational DOFs relative to the pelvis, the upper arm had three rotational DOFs at the shoulder, and elbow and forearm motion were each modeled with 1 DOF revolute joints. Each DOF of the torso and arms were controlled with idealized torque actuators. Hip motion was modeled with three rotational DOFs, and the knee was modeled with one rotational DOF [28]. RSPs and DUPs were modeled with 4 and 1 sagittal plane DOF, respectively, and actuated with idealized torque actuators that reflected the net joint moments computed from inverse dynamics (Fig. 1). The ankle joint of the intact leg and legs of people without TTA was modeled as 1 DOF. Center of mass locations for the residual leg in models of people with TTA were adjusted using experimentally collected residual leg dimensions and defined as a right frustum [32]. Prosthesis mass was measured and applied to the model for each participant. Inertial properties of the RSP were estimated using previous data [33] according to the shape of RSP. DUP inertial properties were estimated using participant-specific DUP mass in combination with the geometry of a biological foot [34,35]. The pylon was modeled as a hollow tube.

Fig. 1
Musculoskeletal modeling and simulation workflow in opensim 3.3 for people with a transtibial amputation(TTA) using RSPs and DUPs and for people without TTA. A total of 108 running simulations were developed.
Fig. 1
Musculoskeletal modeling and simulation workflow in opensim 3.3 for people with a transtibial amputation(TTA) using RSPs and DUPs and for people without TTA. A total of 108 running simulations were developed.
Close modal

Running Simulations.

Experimental GRFs and kinematics were low-pass filtered using a fourth order bidirectional Butterworth filter at 6 Hz. An inverse kinematics solution constrained to model DOFs was computed using a least squares optimization approach [36] (C-Motion, Germantown, MD). A residual reduction algorithm was used to ensure dynamic consistency between the inverse kinematics solution and the ground reaction forces [37]. A computed muscle control algorithm (CMC) was then used to determine individual muscle forces by minimizing the sum of squared muscle controls, which reproduced the running kinematics [38] (Fig. 1). Three gait cycles on each leg of each participant were simulated for a total of 108 running simulations.

Individual muscle forces determined from CMC were used to calculate the three-dimensional vector magnitude of the total hip joint contact force (e.g., Ref. [39]) over the stance and swing phases of running. HJCF and muscle force results were low-pass filtered using a bidirectional fourth order Butterworth with a 6 Hz cutoff frequency and normalized by body weight for each participant. Peak stance and swing phase values of HJCF of the intact and amputated legs of people with TTA using RSPs and DUPs and people without TTA were compared in addition to total gait cycle impulse, which was calculated as the time integral of HJCF from heel strike to the following heel strike. Peak muscle forces of the vasti, hamstrings, gluteus maximus, gluteus medius, rectus femoris, and iliopsoas were compared during the stance and swing phases. Vastus medialis, vastus intermedius, and vastus lateralis were combined to represent total vasti muscle force. Biceps femoris long head, semitendinosus, and semimembranosus were combined to represent total hamstrings muscle group force.

Electromyography signals were demeaned, bandpass filtered (20 Hz–500 Hz, bidirectional fourth order Butterworth), full-wave rectified, and low-pass filtered (6 Hz cutoff, bidirectional fourth-order Butterworth) to create a linear envelope [40]. EMG was normalized to the maximum activation of the muscle during the trial for each participant. Muscle activations from the running simulations were similarly, low-pass filtered (6 Hz cutoff, bidirectional fourth order Butterworth) and normalized to peak activation (Fig. 1). EMG and simulation muscle activations were visually compared to ensure similar onset and offset times [41]. Residual forces and moments from the residual reduction algorithm and kinematic tracking errors were also evaluated as a metric of simulation quality.

Statistical Analysis.

Differences in stance and swing phase peak HJCFs and muscle forces, as well as total gait cycle hip joint contact impulse, were each assessed in rstatisticalcomputingsoftware v 1.1.153 using three separate, two-factor analyses of variance (ANOVAs) [42] with factors of (1) side and (2) prosthesis condition. Separate ANOVAs were used to include within- and between-subjects comparisons, which has been previously used in similar research designs [43]. The side factor had two levels (amputated/left versus intact/right) and the prosthesis condition factor had two levels (RSP versus DUP, RSP versus No TTA, DUP versus No TTA) (Table 2). Tukey's correction for multiple comparisons was used for post hoc pairwise comparisons when significant main or interaction effects were found (ɑ = 0.05).

Table 2

Explanation of statistical design

ANOVADependent variableWithin-subjects factorsBetween-subjects factors
DUP versus no TTAStance/swing phase peak HJCF/muscle force, total gait cycle impulseSide: two levels (R/Intact vs. L/Amputated)Prosthesis condition: two levels (Person with TTA using DUP/Person without TTA)
RSP versus no TTAStance/swing phase peak HJCF/muscle force, total gait cycle impulseSide: two levels (R/intact versus L/amputated)Prosthesis condition: two levels (person with TTA using RSP/person without TTA)
DUP versus RSPStance/swing phase peak HJCF/muscle force, total gait cycle impulseProsthesis condition: two levels (RSP/DUP)None
Side: two levels (intact versus amputated)
ANOVADependent variableWithin-subjects factorsBetween-subjects factors
DUP versus no TTAStance/swing phase peak HJCF/muscle force, total gait cycle impulseSide: two levels (R/Intact vs. L/Amputated)Prosthesis condition: two levels (Person with TTA using DUP/Person without TTA)
RSP versus no TTAStance/swing phase peak HJCF/muscle force, total gait cycle impulseSide: two levels (R/intact versus L/amputated)Prosthesis condition: two levels (person with TTA using RSP/person without TTA)
DUP versus RSPStance/swing phase peak HJCF/muscle force, total gait cycle impulseProsthesis condition: two levels (RSP/DUP)None
Side: two levels (intact versus amputated)

RSP/DUP = running-specific or daily-use prosthesis. No TTA = people without an amputation.

Results

Simulation Quality.

Simulations reproduced the inverse kinematics solution, with the largest root-mean-squared (RMS) rotational errors in arm pronation/supination (i.e., intact side arm pronation = 7.04 ± 13.32 deg), and the largest translational error in vertical translation of the pelvis for people wearing RSPs (1.52 ± 3.02 cm) (Appendix Table 6). Residual forces and moments at the pelvis were less than 1% overall bodyweight/bodyweight-height (Appendix Table 7). Experimental EMG and simulated muscle activations were similar, although simulated gastrocnemius and vastus lateralis peak activations were delayed compared to EMG signals (Fig. 2).

Fig. 2
Experimental EMG signals and simulated muscle activations for people using RSPs, DUPs, and for people without a transtibial amputation (TTA). EMG is displayed as group averages±standard deviation. Vertical dashed lines represent toe-off for the amputated leg and vertical solid lines represent toe-off for intact legs or the legs of people without an amputation.
Fig. 2
Experimental EMG signals and simulated muscle activations for people using RSPs, DUPs, and for people without a transtibial amputation (TTA). EMG is displayed as group averages±standard deviation. Vertical dashed lines represent toe-off for the amputated leg and vertical solid lines represent toe-off for intact legs or the legs of people without an amputation.
Close modal
Table 6

Mean (standard deviation) RMS tracking errors for each degree-of-freedom in the models of people without a transtibial amputation (no TTA) and for people with TTA wearing DUPs and RSPs

DUPRSPNo TTA
AmputatedIntactAmputatedIntact
Hip flexion0.12 (0.67)0.11 (0.59)0.06 (0.28)−0.01 (0.06)1.55 (1.24)
Hip adduction0.06 (0.18)−0.05 (0.18)0.07 (0.17)0.05 (0.24)−0.57 (0.53)
Hip rotation−0.01 (0.04)0.05 (0.07)−0.01 (0.04)−0.04 (0.21)−0.78 (0.44)
Knee angle−0.28 (2.50)0.03 (1.78)−0.48 (1.54)−0.55 (1.02)1.07 (0.53)
Ankle/prosthesis angle0.15 (0.30)0.22 (0.34)RSP4: −0.09 (0.37)0.24 (0.35)0.61 (0.89)
RSP3: 0.18 (0.77)
RSP2: −0.18 (0.77)
RSP1: 0.08 (0.30)
Arm flexion0.23 (2.97)−1.60 (7.21)1.43 (1.85)0.70 (1.19)−0.32 (0.23)
Arm adduction−2.44 (4.60)−0.57 (3.70)−1.53 (4.96)1.23 (2.50)0.01 (0.13)
Arm rotation−1.26 (3.34)−1.06 (2.03)−0.58 (1.32)0.35 (1.23)−0.09 (0.07)
Elbow flexion2.63 (9.23)1.30 (2.21)−0.43 (1.90)0.67 (2.62)0.22 (0.13)
Pronation/supination−1.75 (7.24)0.65 (2.79)3.24 (9.65)7.04 (13.32)−0.09 (1.00)
Pelvis vertical−0.00 (5.18)1.52 (3.02)−0.19 (0.14)
Pelvis anterior0.77 (1.73)1.18 (1.07)−0.06 (0.05)
Pelvis mediolateral−0.06 (3.83)−0.35 (1.93)0.00 (0.03)
Pelvis tilt0.35 (1.55)0.06 (0.16)−1.10 (1.02)
Pelvis list−0.43 (1.41)−0.38 (0.87)0.02 (0.22)
Pelvis rotation−0.35 (1.07)−0.22 (0.43)0.02 (0.25)
Lumbar extension0.24 (2.71)−0.01 (0.22)1.31 (1.08)
Lumbar bending−0.28 (1.05)−0.43 (1.05)0.09 (0.22)
Lumbar rotation0.02(0.21)−0.01 (0.09)−0.04 (0.27)
DUPRSPNo TTA
AmputatedIntactAmputatedIntact
Hip flexion0.12 (0.67)0.11 (0.59)0.06 (0.28)−0.01 (0.06)1.55 (1.24)
Hip adduction0.06 (0.18)−0.05 (0.18)0.07 (0.17)0.05 (0.24)−0.57 (0.53)
Hip rotation−0.01 (0.04)0.05 (0.07)−0.01 (0.04)−0.04 (0.21)−0.78 (0.44)
Knee angle−0.28 (2.50)0.03 (1.78)−0.48 (1.54)−0.55 (1.02)1.07 (0.53)
Ankle/prosthesis angle0.15 (0.30)0.22 (0.34)RSP4: −0.09 (0.37)0.24 (0.35)0.61 (0.89)
RSP3: 0.18 (0.77)
RSP2: −0.18 (0.77)
RSP1: 0.08 (0.30)
Arm flexion0.23 (2.97)−1.60 (7.21)1.43 (1.85)0.70 (1.19)−0.32 (0.23)
Arm adduction−2.44 (4.60)−0.57 (3.70)−1.53 (4.96)1.23 (2.50)0.01 (0.13)
Arm rotation−1.26 (3.34)−1.06 (2.03)−0.58 (1.32)0.35 (1.23)−0.09 (0.07)
Elbow flexion2.63 (9.23)1.30 (2.21)−0.43 (1.90)0.67 (2.62)0.22 (0.13)
Pronation/supination−1.75 (7.24)0.65 (2.79)3.24 (9.65)7.04 (13.32)−0.09 (1.00)
Pelvis vertical−0.00 (5.18)1.52 (3.02)−0.19 (0.14)
Pelvis anterior0.77 (1.73)1.18 (1.07)−0.06 (0.05)
Pelvis mediolateral−0.06 (3.83)−0.35 (1.93)0.00 (0.03)
Pelvis tilt0.35 (1.55)0.06 (0.16)−1.10 (1.02)
Pelvis list−0.43 (1.41)−0.38 (0.87)0.02 (0.22)
Pelvis rotation−0.35 (1.07)−0.22 (0.43)0.02 (0.25)
Lumbar extension0.24 (2.71)−0.01 (0.22)1.31 (1.08)
Lumbar bending−0.28 (1.05)−0.43 (1.05)0.09 (0.22)
Lumbar rotation0.02(0.21)−0.01 (0.09)−0.04 (0.27)

Tracking errors are in degrees except for pelvis translations, which are in centimeters. The RSP degrees-of-freedom represent the most proximal point (RSP4) to the most distal point (RSP1) along the prosthesis.

Table 7

Mean (standard deviation) RMS residual forces and moments applied to the pelvis

RSPDUPNo TTA
Fx (%BW)0.64 (0.53)0.21 (0.83)0.60 (0.62)
Fy (%BW)−0.50 (1.24)−0.49 (4.55)−0.08 (1.25)
Fz (%BW)−0.09 (0.33)−0.14 (0.79)−0.41 (0.72)
Mx (%BW*height)−0.51 (1.06)−0.31 (1.03)−0.33 (1.88)
My (%BW*height)−0.04 (0.11)−0.04 (0.29)−0.10 (0.17)
Mz (%BW*height)0.44 (0.84)0.51 (2.18)−0.91 (1.69)
RSPDUPNo TTA
Fx (%BW)0.64 (0.53)0.21 (0.83)0.60 (0.62)
Fy (%BW)−0.50 (1.24)−0.49 (4.55)−0.08 (1.25)
Fz (%BW)−0.09 (0.33)−0.14 (0.79)−0.41 (0.72)
Mx (%BW*height)−0.51 (1.06)−0.31 (1.03)−0.33 (1.88)
My (%BW*height)−0.04 (0.11)−0.04 (0.29)−0.10 (0.17)
Mz (%BW*height)0.44 (0.84)0.51 (2.18)−0.91 (1.69)

Forces (Fx, Fy, and Fz) are normalized to percentage of body weight (BW), and moments (Mx, My, Mz) are normalized to percent of body weight-height (in meters).

Peak Hip Joint Contact Forces.

During the stance phase, there was a significant main effect of side, and a significant interaction of side × prosthesis condition when comparing people wearing RSPs and people without TTA, in addition to significant main effects of side and prosthesis condition when comparing people wearing RSPs to people wearing DUPs (p ≤ 0.04). The intact leg of people wearing RSPs had a greater peak HJCF compared to the amputated leg (p = 0.02) (Fig. 3, Table 3). Similarly, the intact leg of people wearing DUPs had a greater stance phase peak HJCF compared to the amputated leg (p = 0.05) (Fig. 3, Table 3). The amputated leg of people wearing DUPs had greater stance phase peak HJCF compared to the amputated leg of people wearing RSPs (p < 0.01) (Table 3). Overall stance phase peak HJCF was greater for people wearing DUPs compared to people wearing RSPs (p < 0.01) and intact legs had greater peak HJCFs compared to amputated legs, regardless of prosthesis type (p = 0.03).

Fig. 3
Hip joint contact force for people without TTA and people with TTA wearing RSPs and DUPs while running at 3.5 m/s. Contact force was normalized by body weight and analyzed during the stance and swing phases. Vertical dashed lines represent toe-off for the amputated leg and vertical solid lines represent toe-off for intact legs or the legs of people without an amputation.
Fig. 3
Hip joint contact force for people without TTA and people with TTA wearing RSPs and DUPs while running at 3.5 m/s. Contact force was normalized by body weight and analyzed during the stance and swing phases. Vertical dashed lines represent toe-off for the amputated leg and vertical solid lines represent toe-off for intact legs or the legs of people without an amputation.
Close modal
Table 3

Average (standard deviation) stance and swing phase peak hip joint contact force (HJCF) and total running gait cycle hip joint contact impulse for people with and without a transtibial amputation (TTA) during running

Running-SpecificDaily-Use
IntactAmputatedIntactAmputatedNo TTA
Stance phase peak HJCF (bodyweight multiple)9.14(0.81)*8.37(0.99)10.14 (1.14)*9.68(1.61)9.12(1.13)
Swing phase peak HJCF (bodyweight multiple)7.12(1.06)7.91(1.44)8.2(1.23)8.66(2.04)6.27(0.97)
Total running gait cycle impulse (bodyweight × seconds)3.83(0.58)*3.45(0.41)4.10(0.64)*3.91(0.82)3.56(0.41)
Running-SpecificDaily-Use
IntactAmputatedIntactAmputatedNo TTA
Stance phase peak HJCF (bodyweight multiple)9.14(0.81)*8.37(0.99)10.14 (1.14)*9.68(1.61)9.12(1.13)
Swing phase peak HJCF (bodyweight multiple)7.12(1.06)7.91(1.44)8.2(1.23)8.66(2.04)6.27(0.97)
Total running gait cycle impulse (bodyweight × seconds)3.83(0.58)*3.45(0.41)4.10(0.64)*3.91(0.82)3.56(0.41)

Significant differences compared to the amputated leg of the same prosthesis type are indicated by ‘*’. Significant differences compared to the amputated leg of people wearing DUPs are indicated by “†”. Significant differences compared to people without TTA are indicated by “○”. Significant differences compared to the intact leg of people wearing DUPs are indicated by “▲”.

During the swing phase, there were significant main effects of group for all ANOVAs (p ≤ 0.01) (Table 3). People with TTA wearing DUPs and RSPs had greater peak HJCF compared to people without TTA (p ≤ 0.01) (Fig. 3, Table 3) and people wearing DUPs had greater peak HJCFs compared to people wearing RSPs (p < 0.01). In addition, the intact leg of people wearing DUPs had greater swing phase peak HJCF compared to wearing RSPs (p < 0.01), and the amputated leg of people with TTA had greater peak HJCFs compared to the leg of people without TTA regardless of prosthesis type (p ≤ 0.03) (Fig. 3, Table 3). Finally, the intact leg of people wearing DUPs had greater swing phase peak HJCFs compared to the leg of people without TTA (p = 0.02) (Fig. 3, Table 3).

Hip Joint Contact Impulse.

There was a significant main effect of side and an interaction effect of side × prosthesis condition when comparing people without TTA to people wearing RSPs (p < 0.04). There were also significant main effects of side and prosthesis condition when comparing people with TTA wearing DUPs and RSPs (p < 0.03). The intact leg of people wearing RSPs and DUPs had a greater hip joint contact impulse compared to the amputated leg (p = 0.02, p < 0.01, respectively) (Table 3) over the total gait cycle. In addition, people wearing DUPs had greater overall hip joint contact impulse compared to when they were wearing RSPs (group effect, p < 0.01) and the intact leg had greater hip joint contact impulse compared to the amputated leg, regardless of prosthesis type (p = 0.03).

Gluteus Maximus.

There was a significant main effect of side for people without TTA compared to people using a DUP (p = 0.04) and a significant main effect of prosthesis condition when comparing people using a DUP to people without TTA and to people using an RSP (p ≤ 0.02) over stance. People using DUPs had greater stance phase leg gluteus maximus force compared to people without TTA and compared to people using RSPs (main effect, p ≤ 0.02) (Table 4). The amputated leg of people wearing DUPs had greater peak gluteus maximus force compared to people without TTA (p = 0.03) (Fig. 4, Table 4).

Fig. 4
Simulated muscle forces over the running gait cycle. Average muscle force for people with a transtibial amputation (TTA) using RSPs and DUPs. Average±standard deviation of muscle force for people without TTA is shown. Vertical dashed lines represent toe-off for the amputated leg and vertical solid lines represent toe-off for intact legs or the leg of a person without an amputation.
Fig. 4
Simulated muscle forces over the running gait cycle. Average muscle force for people with a transtibial amputation (TTA) using RSPs and DUPs. Average±standard deviation of muscle force for people without TTA is shown. Vertical dashed lines represent toe-off for the amputated leg and vertical solid lines represent toe-off for intact legs or the leg of a person without an amputation.
Close modal
Table 4

Average peak (standard deviation) muscle forces normalized by body weight during the stance phase of running

Stance phase peak muscle force (bodyweight multiple)
RSPDUP
IntactAmputatedIntactAmputatedNo TTA
Gluteus maximus1.74 (0.56)1.91 (0.45)2.05 (0.47)2.53 (0.78)*1.57 (0.35)
Gluteus medius2.88 (0.28)2.48 (0.42)2.99 (0.70)2.80 (0.75)2.82 (0.44)
Vasti7.71 (1.65)*†3.72 (1.56)*7.21 (1.18)*†4.31 (1.68)*4.74 (1.08)
Hamstrings2.46 (0.11)*†2.03 (0.27)*2.19 (0.27)2.17 (0.41)1.86 (0.37)
Rectus Femoris1.01 (0.39)0.80 (0.17)1.09 (0.39)0.91 (0.20)1.06 (0.61)
Iliopsoas1.67 (0.26)1.77 (0.36)1.73 (0.45)†1.91 (0.38)1.69 (0.24)
Stance phase peak muscle force (bodyweight multiple)
RSPDUP
IntactAmputatedIntactAmputatedNo TTA
Gluteus maximus1.74 (0.56)1.91 (0.45)2.05 (0.47)2.53 (0.78)*1.57 (0.35)
Gluteus medius2.88 (0.28)2.48 (0.42)2.99 (0.70)2.80 (0.75)2.82 (0.44)
Vasti7.71 (1.65)*†3.72 (1.56)*7.21 (1.18)*†4.31 (1.68)*4.74 (1.08)
Hamstrings2.46 (0.11)*†2.03 (0.27)*2.19 (0.27)2.17 (0.41)1.86 (0.37)
Rectus Femoris1.01 (0.39)0.80 (0.17)1.09 (0.39)0.91 (0.20)1.06 (0.61)
Iliopsoas1.67 (0.26)1.77 (0.36)1.73 (0.45)†1.91 (0.38)1.69 (0.24)

Vasti are the sum of the vastus lateralis, vastus intermedius, and vastus medialis. Hamstrings are the sum of the biceps femoris long head, semitendinosus, semimembranosus. “*” indicates difference compared to people without TTA. “†” indicates difference compared to the amputated leg of the same prosthesis condition.

For swing phase, there was a significant prosthesis condition effect when comparing people wearing RSPs to people wearing DUPs (p = 0.02). People wearing DUPs had greater peak gluteus maximus force compared to people wearing RSPs (Table 5).

Table 5

Average peak (standard deviation) muscle forces normalized by body weight during the swing phase of running

Swing phase peak muscle force (bodyweight multiple)
RSPDUP
IntactAmputatedIntactAmputatedNo TTA
Gluteus maximus0.83 (0.44)1.00 (0.60)1.34 (0.65)1.34 (0.44)0.84 (0.46)
Gluteus medius1.94 (0.65)*1.89 (0.53)*2.40 (0.57)*2.15 (0.40)*1.03 (0.35)
Vasti2.00 (0.88)†1.08 (0.49)1.99 (0.97)1.61 (0.39)1.12 (0.38)
Hamstrings2.30 (0.27)†1.77 (0.41)2.06 (0.25)2.16 (0.18)1.88 (0.34)
Rectus Femoris1.72 (0.33)†1.11 (0.29)1.71 (0.45)1.58 (0.23)1.52 (0.34)
Iliopsoas1.59 (0.27)1.83 (0.23)1.74 (0.34)1.79 (0.22)1.66 (0.20)
Swing phase peak muscle force (bodyweight multiple)
RSPDUP
IntactAmputatedIntactAmputatedNo TTA
Gluteus maximus0.83 (0.44)1.00 (0.60)1.34 (0.65)1.34 (0.44)0.84 (0.46)
Gluteus medius1.94 (0.65)*1.89 (0.53)*2.40 (0.57)*2.15 (0.40)*1.03 (0.35)
Vasti2.00 (0.88)†1.08 (0.49)1.99 (0.97)1.61 (0.39)1.12 (0.38)
Hamstrings2.30 (0.27)†1.77 (0.41)2.06 (0.25)2.16 (0.18)1.88 (0.34)
Rectus Femoris1.72 (0.33)†1.11 (0.29)1.71 (0.45)1.58 (0.23)1.52 (0.34)
Iliopsoas1.59 (0.27)1.83 (0.23)1.74 (0.34)1.79 (0.22)1.66 (0.20)

Vasti are the sum of the vastus lateralis, vastus intermedius, and vastus medialis. Hamstrings are the sum of the biceps femoris long head, semitendinosus, semimembranosus. “*” indicates difference compared to people without TTA. “†” indicates difference compared to the amputated leg of the same prosthesis condition.

Hamstrings.

During stance phase, there were significant main effects of side and prosthesis condition when comparing people wearing RSPs to people without TTA (p = 0.02, p = 0.05) and a significant main effect of side when comparing people wearing RSPs and DUPs (p = 0.03). There was significant side × prosthesis condition interaction when comparing people wearing RSPs to people without TTA (p < 0.01) and to people wearing DUPs (p = 0.04). People wearing RSPs had greater stance phase peak hamstrings force compared to people without TTA (p = 0.048), and the intact leg had greater peak hamstrings force compared to the amputated leg (p < 0.01) and compared to people without TTA (p = 0.04). Finally, the intact leg had greater peak hamstrings force compared to the amputated leg regardless of prosthesis type (p = 0.03) (Fig. 4, Table 4).

During swing phase, there was a significant main effect of side and interaction of side × prosthesis condition when comparing people wearing RSPs to people without TTA and people wearing DUPs (p ≤ 0.04). The intact leg of people wearing RSPs had greater swing phase peak hamstrings force compared to the amputated leg (p < 0.01) (Table 5).

Iliopsoas.

In stance, there was a significant effect of side when comparing people without TTA to people wearing DUPs (p = 0.02) (Table 3). The amputated leg of people wearing DUPs had greater stance phase peak iliopsoas muscle force compared to the intact leg (p = 0.04) (Fig. 4, Table 4).

There were no significant differences in iliopsoas peak force during swing phase.

Gluteus Medius and Rectus Femoris.

There were no significant differences in peak gluteus medius or rectus femoris force for prosthesis condition or side during the stance phase.

During swing phase, there was a significant effect of prosthesis condition when comparing gluteus medius force of people without TTA and people with TTA. People with TTA had greater peak swing phase gluteus medius force when wearing DUPs (group effect, p < 0.01) and RSPs (p < 0.01). The amputated and intact leg of people with TTA had greater peak force compared to people without TTA (p ≤ 0.05).

For the rectus femoris during swing phase, there was a significant main effect of side when comparing people wearing RSPs to people without TTA and to people wearing DUPs (p ≤ 0.03) and a significant interaction of side × prosthesis condition when comparing people without TTA to people wearing RSPs. The intact leg had greater peak rectus femoris force during swing phase compared to the amputated leg (p < 0.01) (Table 5).

Vasti.

There were significant main effects of side (p < 0.01) when comparing all groups and a significant side × prosthesis condition interaction for people without TTA compared to people with TTA (p < 0.01) during stance. The intact leg of people wearing DUPs and RSPs had greater stance phase peak vasti force compared to the amputated leg (p < 0.01) and compared to people without TTA (p ≤ 0.04). The intact leg had greater overall force compared to the amputated leg regardless of prosthesis type (p < 0.01) (Fig. 4, Table 4).

During swing phase, there was a significant effect of prosthesis condition when comparing people without TTA to people wearing DUPs (p = 0.04) and a significant effect of side when comparing people wearing RSPs and DUPs (p = 0.03). There was a significant side × prosthesis condition effect when comparing people without TTA to people wearing RSPs (p = 0.05). People wearing DUPs had greater peak vasti force during swing phase compared to people without TTA and the intact leg of people with TTA had greater peak vasti force compared to the amputated leg (p = 0.03) (Table 5).

Discussion

There were no differences in stance phase peak HJCF and total gait cycle impulse between people with and without TTA, contrary to our expectations. Using an RSP significantly reduced bilateral HJCF and impulse and resulted in smaller amputated leg peak HJCF compared to the use of a DUP during stance phase, which supported our hypothesis. However, people with TTA had greater intact leg stance phase peak HJCF and total gait cycle impulse compared to the amputated leg regardless of prosthesis type, which did not support our hypothesis. The intact leg of people wearing RSPs had a smaller peak HJCF during swing phase compared to the use of a DUP. Smaller amputated side vasti muscle force for people with TTA regardless of prosthesis type and greater amputated side gluteus maximus force when using a DUP both likely contributed to the resulting asymmetric HJCFs during stance phase. During swing phase, larger peak gluteus medius force for people with TTA likely contributed to greater peak HJCF for people with TTA compared to people without TTA. Muscular weakness, overloading of a joint, and asymmetric loading are all risk factors for the onset of OA [4,9]. Our results suggest that these risk factors are greater in people with TTA compared to people without TTA, especially given asymmetric running biomechanics [15].

Hip joint contact impulse is important to consider as cumulative loading is associated with the development of OA [44]. Greater impulsive joint loading indicates higher cumulative loading, as a high impulse over many gait cycles results in greater cumulative load. The intact leg of people with TTA had a greater total gait cycle impulse compared to the amputated leg (Table 3), which may be detrimental for long-term joint health [44]. Reduced overall hip joint contact impulse with the use of an RSP may be beneficial for long-term joint health during running. Changes in peak HJCF are driven by differences in muscle forces, which we also observed across conditions and legs.

During stance phase, the use of a DUP resulted in greater peak gluteus maximus force compared to people without TTA or people using RSPs (p ≤ 0.03). The gluteus maximus is important for providing body support and propulsion during the stance phase of running [45], and also contributes to positive hip joint work. Smaller propulsion from DUPs compared to RSPs may result in greater required gluteus maximus force as a result of limited energy storage and return capabilities of DUPs. Greater gluteus maximus force may also contribute to greater amputated side positive hip work when using DUPs compared to RSPs, which has been previously observed in these participants [15]. Smaller gluteus maximus peak activation in the amputated leg of people wearing RSPs compared to DUPs likely contributes to smaller peak amputated HJCF (Fig. 3), especially as gluteus maximus is a large contributor to HJCF during walking [20].

Gluteus medius is important for providing body support during the stance phase of running for people without TTA [45], is important for medial/lateral balance control, and is the largest contributor to HJCF during walking [20] for people without TTA. People with TTA characteristically favor weight-bearing on the intact leg, which results in associated poor amputated side hip abductor strength [46]. While no significant differences in gluteus medius peak force were found during stance phase, there was greater overall bilateral peak gluteus medius force during swing phase shortly after toe-off in people with TTA compared to people without TTA, regardless of prosthesis type. The role of the gluteus medius during the swing phase of running is not as well-understood as it generally has greater levels of activation during stance. However, different components of the gluteus medius control both hip internal and external rotation, in addition to hip abduction. Altered swing limb trajectories in people with TTA may be indicative of greater force requirements from this muscle.

The iliopsoas contributes to leg swing initiation by transferring power from the trunk to the amputated leg during walking [47] and may produce greater force in people using DUPs due to the lack of biarticular gastrocnemius function, and associated smaller propulsion compared to RSPs and biological legs (Fig. 4, Table 4). Observed greater peak force from the amputated side iliopsoas for people wearing DUPs (Fig. 4) may contribute to greater amputated side HJCF near the end of stance phase (Fig. 3).

The vasti are also important to consider in interpreting HJCF, even though they do not span the hip joint. The vasti muscles, which are uni-articular knee extensors, substantially contribute to HJCF through dynamic coupling by accelerating the knee and hip joints into extension [20,48]. Smaller vasti force output during stance phase likely contributes to the observed smaller stance phase peak HJCF and impulse in the amputated leg compared to the intact leg for people with TTA, particularly as the hip is extending during this period (Figs. 3 and 4). Supporting this notion, people with TTA have atrophied quadriceps resulting in smaller force output [49,50]. While we did not account for asymmetric muscle strength in the musculoskeletal model, the optimization predicted reduced force requirements from the amputated leg vasti muscles. Thus, reduced thigh muscle strength is consistent with our observed reduction in simulated amputated leg vasti force (Fig. 4, Table 4), and the large difference between amputated and intact legs, regardless of prosthesis type. This asymmetry in vasti force likely contributes to peak HJCF and impulse asymmetry (Figs. 3 and 4).

People with TTA included in our study used their own DUPs and RSPs, which was intended to capture the habituated running mechanics of people with TTA. However, variation due to prosthesis design is important to consider and is a potential limitation to this study. For example, participant 4 used a DUP (Table 1) that was much more responsive than a traditional DUP design as it allowed for greater deflection, thus storing and returning more energy. The use of this device decreased amputated leg peak HJCF relative to the intact leg (Fig. 5). In contrast, participant 2 used a DUP with a geometry that had a smaller deflection range (Ottobock Advantage DP2) and resulted in amputated leg peak HJCF that was much greater than the intact leg. Further, individual experience may play a significant role in running mechanics as well as device choice. For example, in addition to a more responsive DUP, Participant 4 was a consistent runner, with frequent access to prosthetic adjustments. Individual differences in running biomechanics and prosthesis choice likely resulted in variable results in magnitude and timing of hip joint loading and muscle forces across participants. The effect of specific device characteristics and running experience on running biomechanics should be further investigated in future work.

Fig. 5
Hip joint contact force (HJCF) during running at 3.5 m/s for participants with a transtibial amputation (TTA) wearing RSPs/DUPs. Average hip joint contact force (±one standard deviation) for people without TTA is shaded in gray.
Fig. 5
Hip joint contact force (HJCF) during running at 3.5 m/s for participants with a transtibial amputation (TTA) wearing RSPs/DUPs. Average hip joint contact force (±one standard deviation) for people without TTA is shaded in gray.
Close modal

Our overall loading magnitude results (peak HJCF: 8-11 times body weight (BW)) were in the range of multiple previous experimental and modeling studies, which have examined HJCF. Specifically, in vivo HJCF measurements from instrumented hip implants during running have ranged between 5 and 6 times bodyweight [51,52]. However, results from instrumented implants were collected at a lower running velocity (2.2 m/s) and may be influenced by surgical approaches that alter gait biomechanics [53]. Musculoskeletal modeling has also been used to examine HJCF during running at 3.33 m/s using a static optimization approach with peak loads estimated at approximately 9 times bodyweight during the stance phase and approximately 5 times bodyweight during swing phase [54], similar to our results.

While model muscle activations and EMG were similar, we observed timing shifts in the vastus lateralis and gastrocnemius (Fig. 2), and this result should be acknowledged as a limitation. Changes in timing between model muscle activations and EMG may be due to limitations both in experimental muscle activity measurement, modeling of muscle activation and deactivation dynamics, and the movement optimization algorithm. EMG signals are sensitive to sensor placement and movement artifact. We followed recommended SENIAM placement guidelines [55]; however, cross-talk between muscles may have occurred. Further, EMG signals reflect activity at the muscle belly, rather than the muscle activation occurring along the length of the muscle. Timing shifts may also be a result of differences between the simulations and the participants in electromechanical delay. Electromechanical delay in vivo has been observed at levels exceeding 100 ms [56]. In addition, the muscle recruitment problem that is solved within CMC may result in different muscle activity results across agonist muscles. For example, semitendinosus, semimembranosus, and the biceps femoris long head are biarticular hip extensors, and thus may be activated to varying levels in a group, which may be different than the singular EMG signal that is experimentally collected. Even with some shifts in timing, overall shape of the simulated muscle activation and EMG signals was comparable. Further, the simulations tracked the experimental running kinematics well with small residual forces and kinematic tracking errors, giving us confidence in simulation results. In addition, peak HJCFs from our simulations were similar to previously collected experimental and simulated values.

Musculoskeletal modeling and simulation are useful to estimate in vivo loading and its implications for injury risk. Incorporating muscle forces as a part of estimating hip joint contact forces is important given their substantial contributions to joint loading. In a previous study of knee joint contact loading, muscular forces added approximately two more bodyweight multiples of force during walking compared to simulations only accounting for intersegmental loading [57]. These results suggest that joint contact load differences at the hip for people with TTA are largely affected by changes in muscular contraction, rather than changes in prosthesis mass. Muscles adapt to different device and skeletal characteristics, and it is challenging to separate the effects of device alone from the effects of muscles and their associated joint kinetics that have adapted to using a specific device. However, specifying relative contributions to HJCF from muscular load, prostheses, and skeletal characteristics may be an important area of future work to inform device design and movement training.

Conclusion

Runners with TTA using RSPs had smaller overall peak HJCFs during the stance phase and smaller intact leg peak HJCFs during the swing phase compared to when using DUPs. People with TTA also had greater peak HJCFs during swing compared to people without TTA. The presence of an amputation resulted in greater stance phase intact leg peak HJCF and impulse over the gait cycle compared to the amputated leg and greater overall peak HJCFs during swing compared to people without TTA, which may be detrimental to long-term hip joint health. However, the use of an RSP decreased overall hip joint force impulse compared to the use of a DUP during the gait cycle, which may be especially important for repeated, high impact activities such as running. The use of an RSP also resulted in smaller peak HJCFs just after toe-off in the amputated leg compared to the use of a DUP, which may be a result of the lighter RSP. HJCFs are dependent on muscle force contributions from leg muscles. The use of a DUP resulted in greater amputated leg gluteus maximus peak force compared to people without TTA and compared to people using RSPs during stance, which may reflect the lack of energy storage and return of the DUP compared to the RSP. Swing phase peak gluteus medius forces from both legs were greater for people with TTA compared to people without TTA. Vasti atrophy in the amputated leg likely also plays an important role in the asymmetric HJCF observed for people with TTA regardless of the prosthesis type they use.

Acknowledgment

This work was supported by the Office of the Assistant Secretary of Defense for Health Affairs through the Orthotics and Prosthetics Outcomes Research Program under Award No. W81XWH-15-1-0518. Opinions, interpretations, conclusions, and recommendations are those of the authors and are not necessarily endorsed by the Department of Defense. Additional support was provided by the 2018 American Society of Biomechanics Graduate Student Grant-In-Aid.

Appendix

References

1.
Power
,
L.
,
2016
,
Amputee Sport and Physical Activity Survey
,
England
,
Lingfield, UK
.
2.
Singh
,
R.
,
Hunter
,
J.
, and
Philip
,
A.
,
2007
, “
The Rapid Resolution of Depression and Anxiety Symptoms After Lower Limb Amputation
,”
Clin. Rehab.
,
21
(
8
), pp.
754
759
.10.1177/0269215507077361
3.
Kavanagh
,
T.
,
1983
, “
Exercise and the Heart
,”
Ann. Acad. Med. Singapore
,
12
(
3
), pp.
331
337.
https://pubmed.ncbi.nlm.nih.gov/6378048/
4.
Morgenroth
,
D. C.
,
Segal
,
A. D.
,
Zelik
,
K. E.
,
Czerniecki
,
J. M.
,
Klute
,
G. K.
,
Adamczyk
,
P. G.
,
Orendurff
,
M. S.
,
Hahn
,
M. E.
,
Collins
,
S. H.
, and
Kuo
,
A. D.
,
2011
, “
The Effect of Prosthetic Foot Push-Off on Mechanical Loading Associated With Knee Osteoarthritis in Lower Extremity Amputees
,”
Gait Posture
,
34
(
4
), pp.
502
507
.10.1016/j.gaitpost.2011.07.001
5.
Struyf
,
P. A.
,
van Heugten
,
C. M.
,
Hitters
,
M. W.
, and
Smeets
,
R. J.
,
2009
, “
The Prevalence of Osteoarthritis of the Intact Hip and Knee Among Traumatic Leg Amputees
,”
Arch. Phys. Med. Rehab.
,
90
(
3
), pp.
440
446
.10.1016/j.apmr.2008.08.220
6.
Norvell
,
D. C.
,
Czerniecki
,
J. M.
,
Reiber
,
G. E.
,
Maynard
,
C.
,
Pecoraro
,
J. A.
, and
Weiss
,
N. S.
,
2005
, “
The Prevalence of Knee Pain and Symptomatic Knee Osteoarthritis Among Veteran Traumatic Amputees and Nonamputees
,”
Arch. Phys. Med. Rehab.
,
86
(
3
), pp.
487
493
.10.1016/j.apmr.2004.04.034
7.
Kulkarni
,
J.
,
Adams
,
J.
,
Thomas
,
E.
, and
Silman
,
A.
,
1998
, “
Association Between Amputation, Arthritis and Osteopenia in British Male War Veterans With Major Lower Limb Amputations
,”
Clin. Rehab.
,
12
(
4
), pp.
348
353
.10.1191/026921598672393611
8.
Burke
,
M. J.
,
Roman
,
V.
, and
Wright
,
V.
,
1978
, “
Bone and Joint Changes in Lower Limb Amputees
,”
Ann. Rheumatic Diseases
,
37
(
3
), pp.
252
254
.10.1136/ard.37.3.252
9.
Felson
,
D. T.
,
2000
, “
Osteoarthritis: New Insights Part 1: The Disease and Its Risk Factors
,”
Ann. Intern. Med.
,
133
(
8
), pp.
635
639
.10.7326/0003-4819-133-8-200010170-00016
10.
Radin
,
E. L.
,
Burr
,
D. B.
,
Caterson
,
B.
,
Fyhrie
,
D.
,
Brown
,
T. D.
, and
Boyd
,
R. D.
,
1991
, “
Mechanical Determinants of Osteoarthritis
,”
Semin. Arthritis Rheum.
,
21
(
3
), pp.
12
21
.10.1016/0049-0172(91)90036-Y
11.
Ganz
,
R.
,
Leunig
,
M.
,
Leunig-Ganz
,
K.
, and
Harris
,
W. H.
,
2008
, “
The Etiology of Osteoarthritis of the Hip
,”
Clin. Orthop. Relat. Res.
,
466
(
2
), pp.
264
272
.10.1007/s11999-007-0060-z
12.
Hurwitz
,
D. E.
,
Foucher
,
K. C.
,
Sumner
,
D. R.
,
Andriacchi
,
T. P.
,
Rosenberg
,
A. G.
, and
Galante
,
J. O.
,
1998
, “
Hip Motion and Moments During Gait Relate Directly to Proximal Femoral Bone Mineral Density in Patients With Hip Osteoarthritis
,”
J. Biomech.
,
31
(
10
), pp.
919
925
.10.1016/S0021-9290(98)00097-9
13.
Spector
,
T. D.
,
Harris
,
P. A.
,
Hart
,
D. J.
,
Cicuttini
,
F. M.
,
Nandra
,
D.
,
Etherington
,
J.
,
Wolman
,
R. L.
, and
Doyle
,
D. V.
,
1996
, “
Risk of Osteoarthritis Associated With Long-Term Weight-Bearing Sports: A Radiologic Survey of the Hips and Knees in Female Ex-Athletes and Population Controls
,”
Arthritis Rheumatism
,
39
(
6
), pp.
988
995
.10.1002/art.1780390616
14.
Lloyd
,
C. H.
,
Stanhope
,
S. J.
,
Davis
,
I. S.
, and
Royer
,
T. D.
,
2010
, “
Strength Asymmetry and Osteoarthritis Risk Factors in Unilateral Trans-Tibial, Amputee Gait
,”
Gait Posture
,
32
(
3
), pp.
296
300
.10.1016/j.gaitpost.2010.05.003
15.
Sepp
,
L. A.
,
Baum
,
B. S.
,
Nelson-Wong
,
E.
, and
Silverman
,
A. K.
,
2020
, “
Joint Work and Ground Reaction Forces During Running With Daily-Use and Running-Specific Prostheses
,”
J. Biomech.
,
101
, p.
109629
.10.1016/j.jbiomech.2020.109629
16.
Beck
,
O. N.
,
Taboga
,
P.
, and
Grabowski
,
A. M.
,
2016
, “
Characterizing the Mechanical Properties of Running-Specific Prostheses
,”
PLoS One
,
11
(
12
), p.
e0168298
.10.1371/journal.pone.0168298
17.
Van Jaarsveld
,
H. W. L.
,
Grootenboer
,
H. J.
,
De Vries
,
J.
, and
Koopman
,
H. F.
,
1990
, “
Stiffness and Hysteresis Properties of Some Prosthetic Feet
,”
Prosthet. Orthot. Int.
,
14
(
3
), pp.
117
124
.10.3109/03093649009080337
18.
Wesseling
,
J.
,
Bastick
,
A. N.
,
Ten Wolde
,
S.
,
Kloppenburg
,
M.
,
Lafeber
,
F. P.
,
Bierma-Zeinstra
,
S. M.
, and
Bijlsma
,
J. W.
,
2015
, “
Identifying Trajectories of Pain Severity in Early Symptomatic Knee Osteoarthritis: A 5-Year Followup of the Cohort Hip and Cohort Knee (CHECK) Study
,”
J. Rheumatol.
,
42
(
8
), pp.
1470
1477
.10.3899/jrheum.141036
19.
Prince
,
F.
,
Allard
,
P.
,
Therrien
,
R.
, and
McFadyen
,
B. J.
,
1992
, “
Running Gait Impulse Asymmetries in Below-Knee Amputees
,”
Prosthet. Orthot. Int.
,
16
(
1
), pp.
19
24
.10.3109/03093649209164303
20.
Correa
,
T. A.
,
Crossley
,
K. M.
,
Kim
,
H. J.
, and
Pandy
,
M. G.
,
2010
, “
Contributions of Individual Muscles to Hip Joint Contact Force in Normal Walking
,”
J. Biomech.
,
43
(
8
), pp.
1618
1622
.10.1016/j.jbiomech.2010.02.008
21.
Isakov
,
E.
,
Keren
,
O.
, and
Benjuya
,
N.
,
2000
, “
Trans-Tibial Amputee Gait: Time-Distance Parameters and EMG Activity
,”
Prosthet. Orthot. Int.
,
24
(
3
), pp.
216
220
.10.1080/03093640008726550
22.
Sepp
,
L. A.
,
Nelson-Wong
,
E.
,
Baum
,
B. S.
, and
Silverman
,
A. K.
,
2020
, “
Running-Specific Prostheses Reduce Lower-Limb Muscle Activity Compared to Daily-Use Prostheses in People With Unilateral Transtibial Amputations
,”
J. Electromyogr. Kinesiol.
,
55
, p.
102462
.10.1016/j.jelekin.2020.102462
23.
Fey
,
N. P.
,
Klute
,
G. K.
, and
Neptune
,
R. R.
,
2011
, “
The Influence of Energy Storage and Return Foot Stiffness on Walking Mechanics and Muscle Activity in Below-Knee Amputees
,”
Clinical Biomech.
,
26
(
10
), pp.
1025
1032
.10.1016/j.clinbiomech.2011.06.007
24.
Hurley
,
G.
,
McKenney
,
R.
,
Robinson
,
M.
,
Zadravec
,
M.
, and
Pierrynowski
,
M.
,
1990
, “
The Role of the Contralateral Limb in Below-Knee Amputee Gait
,”
Prosthet. Orthot. Int.
,
14
(
1
), pp.
33
42
.10.3109/03093649009080314
25.
Zajac
,
F. E.
,
Neptune
,
R. R.
, and
Kautz
,
A.
,
2002
, “
Biomechanics and Muscle Coordination of Human Walking Part I: Introduction to Concepts, Power Transfer, Dynamics and Simulations
,”
Gait Posture
,
16
, pp.
215
232
.10.1016/S0966-6362(02)00068-1
26.
Lu
,
T. W.
,
O'Connor
,
J. J.
,
Taylor
,
S. J. G.
, and
Walker
,
P. S.
,
1997
, “
Validation of a Lower Limb Model With In Vivo Femoral Forces Telemetered From Two Subjects
,”
J. Biomech.
,
31
(
1
), pp.
63
69
.10.1016/S0021-9290(97)00102-4
27.
Baum
,
B. S.
,
Hobara
,
H.
,
Koh
,
K.
,
Kwon
,
H. J.
,
Miller
,
R. H.
, and
Shim
,
J. K.
,
2019
, “
Amputee Locomotion: Joint Moment Adaptations to Running Speed Using Running-Specific Prostheses After Unilateral Transtibial Amputation
,”
Am. J. Phys. Med. Rehab.
, 98(3), pp.
182
190
.10.1097/PHM.0000000000000905
28.
Rajagopal
,
A.
,
Dembia
,
C. L.
,
DeMers
,
M. S.
,
Delp
,
D. D.
,
Hicks
,
J. L.
, and
Delp
,
S. L.
,
2016
, “
Full-Body Musculoskeletal Model for Muscle-Driven Simulation of Human Gait
,”
IEEE Trans. Biomed. Eng.
,
63
(
10
), pp.
2068
2079
.10.1109/TBME.2016.2586891
29.
Lai
,
A. K. M.
,
Arnold
,
A. S.
, and
Wakeling
,
J. M.
,
2017
, “
Why Are Antagonist Muscles Co-Activated in My Simulation? A Musculoskeletal Model for Analysing Human Locomotor Tasks
,”
Ann. Biomed. Eng.
,
45
(
12
), pp.
2762
2713
.10.1007/s10439-017-1920-7
30.
Delp
,
S. L.
,
Loan
,
J. P.
,
Hoy
,
M. G.
,
Zajac
,
F. E.
,
Topp
,
E. L.
, and
Rosen
,
J. M.
,
1990
, “
An Interactive Graphics-Based Model of the Lower Extremity to Study Orthopaedic Surgical Procedures
,”
IEEE Trans. Biomed. Eng.
,
37
(
8
), pp.
757
767
.10.1109/10.102791
31.
Millard
,
M.
,
Uchida
,
T.
,
Seth
,
A.
, and
Delp
,
S. L.
,
2013
, “
Flexing Computational Muscle: Modeling and Simulation of Musculotendon Dynamics
,”
ASME J. Biomech. Eng.
,
135
(
2
), p.
021005
.10.1115/1.4023390
32.
Hanavan
,
E. P.
,
1964
, A Mathematical Model of the Human Body, Air Force Aerospace Medical Research Lab, Wright-Patterson AFB, OH, No.
AFIT-GA-PHYS-64-3
.https://pubmed.ncbi.nlm.nih.gov/14243640/
33.
Baum
,
B. S.
,
Schultz
,
M. P.
,
Tian
,
A.
,
Shefter
,
B.
,
Wolf
,
E. J.
,
Kwon
,
H. J.
, and
Shim
,
J. K.
,
2013
, “
Amputee Locomotion: Determining the Inertial Properties of Running-Specific Prostheses
,”
Arch. Phys. Med. Rehab.
,
94
(
9
), pp.
1776
1783
.10.1016/j.apmr.2013.03.010
34.
Miller
,
D. I.
,
1987
, “
Resultant Lower Extremity Joint Moments in Below-Knee Amputees During Running Stance
,”
J. Biomech.
,
20
(
5
), pp.
529
541
.10.1016/0021-9290(87)90253-3
35.
Sanderson
,
D.
, and
Martin
,
P.
,
1996
, “
Joint Kinetics in Unilateral Below-Knee Amputee Patients During Running
,”
Arch. Phys. Med. Rehab.
,
77
(
12
), pp.
1279
1285
.10.1016/S0003-9993(96)90193-8
36.
Lu
,
T. W.
, and
O'Connor
,
J. J.
,
1999
, “
One Position Estimation From Skin Marker co - Ordinates Using Global Optimisation With Joint Constraints
,”
J. Biomech.
,
32
(
2
), pp.
129
134
.10.1016/S0021-9290(98)00158-4
37.
Delp
,
S. L.
,
Anderson
,
F. C.
,
Arnold
,
A. S.
,
Loan
,
P.
,
Habib
,
A.
,
John
,
C. T.
,
Guendelman
,
E.
, and
Thelen
,
D. G.
,
2007
, “
OpenSim: Open Source to Create and Analyze Dynamic Simulations of Movement
,”
IEEE Trans. Bio-Med. Eng.
,
54
(
11
), pp.
1940
1950
.10.1109/TBME.2007.901024
38.
Thelen
,
D. G.
, and
Anderson
,
F. C.
,
2006
, “
Using Computed Muscle Control to Generate Forward Dynamic Simulations of Human Walking From experimental data
,”
J. Biomech.
,
39
(
6
), pp.
1107
1115
.10.1016/j.jbiomech.2005.02.010
39.
Steele
,
K. M.
,
DeMers
,
M. S.
,
Schwartz
,
M. H.
, and
Delp
,
S. L.
,
2012
, “
Compressive Tibiofemoral Force During Crouch Gait
,”
Gait Posture
, 35(4), pp.
556
560
.10.1016/j.gaitpost.2011.11.023
40.
De Luca
,
C. J.
,
Donald Gilmore
,
L.
,
Kuznetsov
,
M.
, and
Roy
,
S. H.
,
2010
, “
Filtering the Surface EMG Signal: Movement Artifact and Baseline Noise Contamination
,”
J. Biomech.
,
43
(
8
), pp.
1573
1579
.10.1016/j.jbiomech.2010.01.027
41.
Hicks
,
J. L.
,
Uchida
,
T. K.
,
Seth
,
A.
,
Rajagopal
,
A.
, and
Delp
,
S. L.
,
2015
, “
Is my Model Good Enough? Best Practices for Verification and Validation of Musculoskeletal Models and Simulations of Movement
,”
ASME J. Biomech. Eng.
,
137
(
2
), p.
020905
.10.1115/1.4029304
42.
Pinheiro
,
J.
,
Bates
,
D.
,
DebRoy
,
S.
, and
Sarkar
,
D.
, and
2013
, “
R Development Core Team, NLME: Linear and Nonlinear Mixed Effects Models-R Package Version 3.1-108
,” R Foundation for Statistical Computing, Vienna, Austria.
43.
Pickle
,
N. T.
,
Wilken
,
J. M.
,
Aldridge
,
J. M.
,
Neptune
,
R. R.
, and
Silverman
,
A. K.
,
2014
, “
Whole-Body Angular Momentum During Stair Walking Using Passive and Powered Lower-Limb Prostheses
,”
J. Biomech.
,
47
(
13
), pp.
3380
3389
.10.1016/j.jbiomech.2014.08.001
44.
Miller
,
R. H.
,
2017
, “
Joint Loading in Runners Does Not Initiate Knee Osteoarthritis
,”
Exerc. Sport Sci. Rev
,
45
(
2
), pp.
87
95
.10.1249/JES.0000000000000105
45.
Hamner
,
S. R.
,
Seth
,
A.
, and
Delp
,
S. L.
,
2010
, “
Muscle Contributions to Propulsion and Support During Running
,”
J. Biomech.
,
43
(
14
), pp.
2709
2716
.10.1016/j.jbiomech.2010.06.025
46.
Nadollek
,
H.
,
Brauer
,
S.
, and
Isles
,
R.
,
2002
, “
Outcomes After Trans-Tibial Amputation: The Relationship Between Quiet Stance Ability, Strength of Hip Abductor Muscles and Gait
,”
Physiother. Res. Int. J. Res. Clin. Phys. Ther.
,
7
(
4
), pp.
203
214
.10.1002/pri.260
47.
Silverman
,
A. K.
, and
Neptune
,
R. R.
,
2012
, “
Muscle and Prosthesis Contributions to Amputee Walking Mechanics: A Modeling Study
,”
J. Biomech.
,
45
(
13
), pp.
2271
2278
.10.1016/j.jbiomech.2012.06.008
48.
Zajac
,
F. E.
, and
Gordon
,
M. E.
,
1989
, “
Determining Muscle's Force and Action in Multi-Articular Movement
,”
Exer. Sport Sci. Rev.
,
17
(
1
), pp.
187
230
.https://pubmed.ncbi.nlm.nih.gov/2676547/
49.
Renström
,
P.
,
Grimby
,
G.
, and
Larsson
,
E.
,
1983
, “
Thigh Muscle Strength in Below-Knee Amputees
,”
Scand. J. Rehab. Med.
,
9
, pp.
163
173
. https://pubmed.ncbi.nlm.nih.gov/6585938/
50.
Isakov
,
E.
,
Burger
,
H.
,
Gregorič
,
M.
, and
Marinček
,
C.
,
1996
, “
Isokinetic and Isometric Strength of the Thigh Muscles in Below-Knee Amputees
,”
Clin. Biomech.
,
11
(
4
), pp.
233
235
.10.1016/0268-0033(95)00078-X
51.
Bergmann
,
G.
,
Graichen
,
F.
, and
Rohlmann
,
A.
,
1993
, “
Hip Joint Loading During Walking and running, measured in two patients
,”
J. Biomech.
,
26
(
8
), pp.
969
990
.10.1016/0021-9290(93)90058-M
52.
Bergmann
,
G.
,
Bender
,
A.
,
Dymke
,
J.
,
Duda
,
G.
, and
Damm
,
P.
,
2016
, “
Standardized Loads Acting in Hip Implants
,”
PLoS One
,
11
(
5
), p.
e0155612
.10.1371/journal.pone.0155612
53.
Ewen
,
A. M.
,
Stewart
,
S.
,
Gibson
,
A. S. C.
,
Kashyap
,
S. N.
, and
Caplan
,
N.
,
2012
, “
Post-Operative Gait Analysis in Total Hip Replacement Patients—A Review of Current Literature and Meta-Analysis
,”
Gait Posture
,
36
(
1
), pp.
1
6
.10.1016/j.gaitpost.2011.12.024
54.
Giarmatzis
,
G.
,
Jonkers
,
I.
,
Wesseling
,
M.
,
Van Rossom
,
S.
, and
Verschueren
,
S.
,
2015
, “
Loading of Hip Measured by Hip Contact Forces at Different Speeds of Walking and Running
,”
J. Bone Miner. Res.
,
30
(
8
), pp.
1431
1440
.10.1002/jbmr.2483
55.
Hermens
,
H. J.
,
Freriks
,
B.
,
Merletti
,
R.
,
Stegeman
,
D.
,
Blok
,
J.
,
Rau
,
G.
,
Disselhorst-Klug
,
C.
, and
Hägg
,
G.
,
1999
, “
European Recommendations for Surface Electromyography
,”
Roessingh Res. Dev.
,
8
(
2
), pp.
13
54
. http://www.seniam.org/pdf/contents8.PDF
56.
Blackburn
,
J. T.
,
Bell
,
D. R.
,
Norcross
,
M. F.
,
Hudson
,
J. D.
, and
Engstrom
,
L. A.
,
2009
, “
Comparison of Hamstring Neuromechanical Properties Between Healthy Males and Females and the Influence of Musculotendinous Stiffness
,”
J. Electromyogr. Kinesiol.
,
19
(
5
), pp.
e362
e369
.10.1016/j.jelekin.2008.08.005
57.
Lerner
,
Z. F.
,
Haight
,
D. J.
,
DeMers
,
M. S.
,
Board
,
W. J.
, and
Browning
,
R. C.
,
2014
, “
The Effects of Walking Speed on Tibiofemoral Loading Estimated Via Musculoskeletal Modeling
,”
J. Appl. Biomech.
,
30
(
2
), pp.
197
205
.10.1123/jab.2012-0206