Skip to main content
  • Research article
  • Open access
  • Published:

Effects of body size and load carriage on lower-extremity biomechanical responses in healthy women

Abstract

Background

Musculoskeletal injuries, such as stress fractures, are the single most important medical impediment to military readiness in the U.S. Army. While multiple studies have established race- and sex-based risks associated with a stress fracture, the role of certain physical characteristics, such as body size, on stress-fracture risk is less conclusive.

Methods

In this study, we investigated the effects of body size and load carriage on lower-extremity joint mechanics, tibial strain, and tibial stress-fracture risk in women. Using individualized musculoskeletal-finite-element-models of 21 women of short, medium, and tall statures (n = 7 in each group), we computed the joint mechanics and tibial strains while running on a treadmill at 3.0 m/s without and with a load of 11.3 or 22.7 kg. We also estimated the stress-fracture risk using a probabilistic model of bone damage, repair, and adaptation.

Results

Under all load conditions, the peak plantarflexion moment for tall women was higher than those in short women (p < 0.05). However, regardless of the load condition, we did not observe differences in the strains and the stress-fracture risk between the stature groups. When compared to the no-load condition, a 22.7-kg load increased the peak hip extension and flexion moments for all stature groups (p < 0.05). However, when compared to the no-load condition, the 22.7-kg load increased the strains and the stress-fracture risk in short and medium women (p < 0.05), but not in tall women.

Conclusion

These results show that women of different statures adjust their gait mechanisms differently when running with external load. This study can educate the development of new strategies to help reduce the risk of musculoskeletal injuries in women while running with external load.

Peer Review reports

Background

Overuse musculoskeletal injuries, such as stress fracture, account for 70% of all injuries among military recruits and active Soldiers in the U.S. Army [1]. The high incidence of stress fracture in recruits is attributed to strenuous exercises and other physical activities performed during basic combat training [2], specifically, walking and running while carrying a heavy load. Multiple studies have shown that stress fracture is multifactorial in nature and is influenced by modifiable factors, such as training intensity, training duration, and load carriage, as well as non-modifiable factors, such as race, sex, and body size [3,4,5]. For example, increased duration and frequency of training sessions as well as heavy load carriage can increase the risk of stress fracture [5]. Similarly, there is overwhelming evidence to support the fact that Whites compared to Blacks, women compared to men, and older individuals compared to younger ones are more susceptible to stress fracture [6, 7]. However, the role of certain physical characteristics, such as body size, on stress-fracture risk is less conclusive. For example, while Kelly et al. reported that shorter and lighter individuals have a higher incidence of stress fracture [8], Knapik et al. [9] and Välimäki et al. [10] reported that men who are taller are at an increased risk of stress fracture. In contrast, other studies have shown no association between incidence of stress fracture and body size [7, 11].

Body size influences an individual’s capacity for physical performance [12, 13] and the biomechanics of locomotion [14, 15]. For example, with increasing body height, cross-sectional area of muscle and bone increases [12], which will also lead to higher force development capacity (i.e., muscle strength) [13]. Similarly, individuals with greater body mass tend to present an extended knee posture during walking in order to mitigate body size-related increase in the knee moment [14, 15]. In addition, walking or running with a load causes alteration in gait kinematics and kinetics [16]. Specifically, walking and running with a load reduces stride length, increases cadence, increases hip and knee flexion angles and moments [17, 18], and increases tibial mechanical stress, which could potentially elevate the risk of tibial stress fracture [19]. However, to date, it is still unclear how load carriage influences the biomechanics of running and the associated risk of stress fracture as a function of body size. Such an investigation is critical, especially because there are considerable differences in the body size of enlisted U.S. female Soldiers [height: 148.0–178.1 cm, mass: 46.4–98.3 kg [20]], and regardless of their body size, each Soldier is required to complete the same physical activities while carrying the same load.

Recently, using individualized musculoskeletal-finite-element (M/FE) models, we determined the effects of load carriage on joint mechanics (i.e., kinematics and kinetics) and the spatiotemporal distribution of tibial stress in three women of different sizes, walking with or without load [21]. In the individualized M/FE models, we accounted for subject-specific features, such as body size, bone morphology, and bone mineral density distribution, in addition to motion characteristics. When compared to the predictions from a generic musculoskeletal model, joint reaction (i.e., contact) forces (JRFs) predicted by the individualized musculoskeletal models differed by as much as 22% in the knee and 26% in the ankle [21]. Moreover, tibial stresses predicted by the individualized models were 31% greater than those predicted by the generic finite element (FE) models. These results highlight the importance of individualized M/FE models in assessing the mechanical loads acting on different individuals performing the same activity.

In the current study, we used individualized M/FE models to investigate the effects of body size and load carriage on the joint kinematics and kinetics, tibial strain, and the risk of stress fracture in women when they ran at 3.0 m/s (i.e., a 9.0-min/mile pace) on a treadmill, without and with a load of 11.3 or 22.7 kg. The objective of this study is to determine the effects of women body size in joint kinematics and kinetics as well as strain in the tibia as a function of load carriage during level running.

Using a probabilistic stress-fracture prediction model [22,23,24], we also investigated the risk of stress fracture in these women resulting from running 4.8 km/day (i.e., 3.0 miles/day) for 100 consecutive days [22]. We hypothesized that the changes in the joint mechanics, specifically, those related to peak joint angles, forces, and moments, as well as the strain in the tibia and the likelihood of stress fracture in women running with a load are dependent on body size.

Methods

Imaging and motion-capture data

As reported previously [19], based on sample-size calculations, we enrolled 21 young, healthy women between the ages of 18 to 21 years [mean = 19.6 and standard deviation (SD) = 0.9 years], with a body mass index between 19 and 25 kg/m2 [mean = 22.2 and SD = 1.9 kg/m2]. We measured the height of each subject using a wall-mounted stadiometer and their body mass using a calibrated electronic scale (Table 1). Body mass index was calculated as mass (kg) divided by the square of height (m2). We assigned each subject to one of three stature groups determined a priori: short (n = 7; height, 149–160 cm; mass, 48–61 kg), medium (n = 7; height, 161–166 cm; mass, 54–64 kg), or tall (n = 7; height, 172–177 cm; mass, 59–72 kg). The subjects in the short, medium, and tall groups are less than the 30th percentile, between the 45th and 65th percentile, and greater than the 90th percentile of the U.S. female Soldier population, respectively [20]. All subjects were free of self-reported bone disorders or lower-limb injuries 3 months prior to data collection. For each subject, we obtained quantitative computed tomography (CT; GE Discovery Scanner, General Electric Medical System, Milwaukee, WI) scans of the left leg with an in-plane resolution of 0.49 × 0.49 mm2 and a slice thickness of 0.63 mm (Acquisition setting: 120 kVp and 200 mAs). Each CT scan also included a calibration phantom with known calcium hydroxyapatite concentrations (QRM, Moehrendorf, Germany).

Table 1 Anthropometric characteristics of 21 young, healthy subjects (n = 7 in each group) considered in the study. Subjects ranged from 18 to 21 years of age, with a body mass index (BMI) between 19 and 25 kg/m2. The data are presented as mean ± one standard deviation

We also acquired motion-capture data using an eight-camera motion analysis system (Vicon Nexus, Centennial, CO) that tracked 42 retroreflective markers secured on the subject’s body, including anatomical landmarks and segmental-tracking markers. Specifically, we placed the markers on the left and right anterior and posterior iliac crest, greater trochanters, femoral condyles, heel, and metatarsals. We collected the motion-capture data (at a frequency of 200 Hz) and ground reaction forces (GRF, frequency of 1000 Hz) while the subjects ran at 3.0 m/s on an instrumented treadmill (Bertec, Columbus, OH), carrying no load (0 kg, baseline model) or a load of 11.3 or 22.7 kg (i.e., 25 or 50 lb) using an adjustable weight vest (V-max, Rexburg, ID) [19]. These selected loads are representative of those carried by U.S. Soldiers during training and field operations [25]. We calculated stride frequency by using the time between successive ipsilateral foot contacts, and stride length by dividing the speed of the treadmill by the stride frequency while accounting for differences between ipsilateral foot contact positions. We received approval for the study from the Human Research Protection Office at the U.S. Army Medical Research and Development Command (Ft. Detrick, MD) and the Conjoint Health Research Ethics Board of the University of Calgary (Calgary, AB, Canada). All procedures were performed in accordance with applicable Department of Defense, U.S. Army, and U.S. Army Medical Research and Development Command human subject protection requirements. We obtained informed written consent from each participant prior to the study.

Individualized musculoskeletal analysis

Using individualized musculoskeletal models, we determined joint kinematics and kinetics for each subject while they ran without load, or with a load of 11.3 or 22.7 kg. We developed the individualized models by scaling a generic female musculoskeletal model available in the AnyBody Modeling System™ (AnyBody Technology, Aalborg, Denmark) [20, 21, 26] and then morphed the tibiae of the generic model with the subject-specific tibiae. To scale the generic model, we used the anthropometric measurements of the individual and an optimization scheme that minimized the errors between markers defined in the model and those tracked in the experiment [27]. The individualized musculoskeletal models consisted of seven rigid bodies in the lower extremity, including pelvis, thigh, shank, and foot, with 55 Hill-type muscles per leg. Moreover, in the model we represented the hip joints as spherical joints and the knee and ankle joints as revolute joints [21]. In the musculoskeletal analyses, we computed muscle activities and reaction forces over one gait cycle using an inverse dynamics approach by minimizing the sum of the cubed muscle activities. We defined muscle activity as the ratio of muscle force divided by the muscle strength. We performed the inverse dynamics analyses, using the most representative stride from among the multiple strides that were acquired during the motion-capture experiments [28], to compute joint kinematics and kinetics over one gait cycle for each of the 21 subjects under three different load conditions (i.e., 0, 11.3, and 22.7 kg).

Individualized finite element analysis

Similar to our previous work [21], we developed subject-specific three-dimensional (3-D) FE models of the left tibia. We created a 3-D FE model for each of the 21 subjects using 3-Matic (Materialise, Leuven, Belgium) by generating 10-noded quadratic tetrahedral (C3D10) elements (average edge length varied from 3.0 to 3.5 mm) from the surface geometry of the subject-specific tibia, which was developed from the CT images (Fig. 1). We considered the bone (Poisson’s ratio of 0.325) and intramedullary tissue (Poisson’s ratio of 0.167) regions as linear elastic, isotropic materials. We calculated inhomogeneous material properties of the bone by converting the intensity of the CT images (defined in Hounsfield units) into apparent bone densities, and then to elastic moduli based on the intensity of the calibration phantom.

Fig. 1
figure 1

Representative finite element mesh of the left tibia for short-, medium-, and tall-statured women

We specified the muscle forces as well as the bone forces and moments calculated from the individualized musculoskeletal analysis as loading conditions for the FE model. Specifically, we identified the muscle and ligament insertion points in the musculoskeletal model as FE constraint nodes, which we coupled to the outer surface of the tibial FE meshes. Through this procedure, on average, we specified 171 couplings between the muscle and the ligament nodes with the tibial FE meshes. We defined the bone-to-bone contact in a similar manner. Lastly, for each subject under each load condition, we performed a FE simulation (ABAQUS v6.17, Dassault Systèmes, Vélizy-Villacoublay, France) to determine the peak tibial strain for each element during one gait cycle. As the likelihood of stress fracture is dependent on the distribution of the peak strain in the bone [23, 24, 29], we also determined the 10th, 50th, and 90th percentile peak strains in the tibia during one gait cycle.

Probabilistic model of stress-fracture likelihood

We determined the likelihood of stress fracture in the tibia using a probabilistic model [23, 24, 29], which considers bone failure due to fatigue damage and accounts for bone repair and adaptation. We represented the fatigue life of a bone F using an inverse power-law relationship given as F = C × ε−6.6, where C denotes a constant obtained from in vitro fatigue experiments on bone samples [29, 30] and ε represents the peak von Mises equivalent strain determined from the FE analyses. We represented bone adaptation, which increases the tibial areal moment of inertia and decreases the tibial strain, by assuming a constant lamellar bone deposition rate of 4 μm/day [31]. We accounted for changes in the tibial strain due to bone adaptation by calculating the equivalent tibial strain Δ εe for each element in the FE model, using the following relationship

$$ \Delta {\upvarepsilon}_{\mathrm{e}}={\left(\frac{1}{{\mathrm{t}}_d}\;{\int}_0^{{\mathrm{t}}_d}{\left(\mathrm{R}.\Delta \upvarepsilon \right)}^{\mathrm{n}}\mathrm{dt}\right)}^{1/\mathrm{n}} $$
(1)

where td denotes the total time in days over which adaptation takes place, R denotes the adaptation strain ratio, and n represents a material constant. In this study, we assumed n = 6.6, based on cortical bone fatigue experiments reported by Carter et al. [30]. We represented the cumulative probability of bone repair Pr as [23, 24, 29].

$$ {\mathrm{P}}_{\mathrm{r}}=1-\exp \left[-{\left(\frac{\mathrm{t}}{26}\right)}^2\right] $$
(2)

where t represents the time in days. To account for variability in the fatigue life F within a tibia, we divided the bone into eight regions (N = 8) based on similar strain values obtained from the FE analyses, and then determined the cumulative probability of failure of each region Pn, with n = 1, 2, … ..., N, using the following relationship

$$ {\mathrm{P}}_{\mathrm{n}}=1-\exp \left[-\left(\frac{{\mathrm{V}}_{\mathrm{n}}}{{\mathrm{V}}_0}\right)\;{\left(\frac{\mathrm{t}}{{\mathrm{t}}_{\mathrm{f}}}\right)}^{1.2}\right] $$
(3)

where Vn denotes the volume of the nth elemental region, V0 denotes the volume of the bone sample used in fatigue experiments [22], and tf denotes the reference time to fatigue failure (defined in days). The reference time is a function of the material strength of the bone, the equivalent strain εe, and the number of loading cycles/day [23]. We represented the cumulative probability of fatigue failure of the entire tibia as

$$ {\mathrm{P}}_{\mathrm{f}}=1-{\prod}_{\mathrm{n}=1}^{\mathrm{N}}\left(1-{\mathrm{P}}_{\mathrm{n}}\right). $$

We combined the cumulative probabilities of repair Pr and fatigue failure Pf to determine the likelihood of stress fracture Ps in a tibia, which is given as

$$ {\mathrm{P}}_{\mathrm{s}}={\int}_0^{{\mathrm{t}}_d}{\mathrm{Q}}_{\mathrm{f}}\left(1-{\mathrm{P}}_{\mathrm{r}}\right)\kern0.28em \mathrm{dt} $$
(4)

where Qf denotes the time differential of Pf. We determined the risk of stress fracture for each subject by assuming that each subject ran 4.8 km/day (i.e., 3.0 miles/day) for 100 consecutive days [22]. We estimated the number of loading cycles/day by dividing the daily running distance by the subject’s stride length, which was measured during the motion-capture experiments.

Statistical analysis

We determined the relationships between dependent variables and predictor variables using a linear mixed-effects analysis. We considered the outputs from the individualized musculoskeletal analysis (i.e., peak joint angles, moments, and forces) and those from the individualized FE analysis (i.e., tibial strains), as well as the likelihood of stress fracture as dependent variables. We considered load and stature group, which were categorical factors, as predictor variables. We included a random intercept for each subject to account for the within-subject dependence. To determine the statistical significance of a predictor variable, we used a likelihood-ratio test. We then performed a post-hoc Tukey test on the linear mixed-effects model for pairwise comparisons and corrected the results for multiple comparisons using the Holm-Bonferroni correction. We separately performed between-stature-group comparisons of the kinematic and kinetic responses, tibial strain, and the likelihood of stress fracture for each load condition, using analysis of variance. We represent all data as mean ± one SD, unless otherwise noted. We tested for statistical significance using a criterion of p < 0.05, and performed all analyses in the statistical software package R [32, 33].

Results

Joint kinematics and kinetics

Figures 2, 3 and 4 show the mean and SD of the peak kinematic and kinetic parameters predicted by the individualized musculoskeletal model. Under the baseline condition, the average peak hip flexion angle ranged from 22.0 degrees [(SD = 2.8), tall women] to 26.5 degrees [(SD = 4.2), short women] (Fig. 2, top left). In tall women, when compared to the baseline condition, a load of 11.3 or 22.7 kg significantly increased their peak hip flexion angle to 24.7 degrees [(SD = 2.4), an increase of 12.3%] or 25.4 degrees [(SD = 1.8), 15.5%], respectively (p < 0.001). Similarly, in the medium women group, we observed a statistically significant increase in the peak hip flexion angle when they ran with load (8.2% for 11.3 kg and 13.4% for 22.7 kg; p < 0.01). In short women, when compared to the baseline condition, a 22.7-kg load significantly reduced peak knee flexion angle during stance [i.e., from 54.1 (SD = 6.2) to 50.9 (SD = 4.3) degrees, a 5.9% reduction, p < 0.05]. We did not observe statistically significant differences in peak hip extension and ankle plantarflexion between the three different loads in short, medium, and tall women (Fig. 2, right).

Fig. 2
figure 2

Comparison of peak joint angles at the hip, knee, and ankle under different loads during one gait cycle obtained from the individualized musculoskeletal analyses. Data are expressed as mean ± one standard deviation. *p < 0.05, **p < 0.01, ***p < 0.001

Fig. 3
figure 3

Comparison of peak joint moments at the hip, knee, and ankle under different loads during one gait cycle obtained from the individualized musculoskeletal analyses. Data are expressed as mean ± one standard deviation. *p < 0.05, **p < 0.01, ***p < 0.001

Fig. 4
figure 4

Comparison of peak ground reaction force (GRF) and peak joint reaction force (JRF) at the hip, knee, and ankle under different loads during one gait cycle obtained from the individualized musculoskeletal analyses. Data are expressed as mean ± one standard deviation. *p < 0.05, **p < 0.01, ***p < 0.001

When compared to the baseline condition, running with a load of 22.7 kg increased peak hip extension moment and flexion moment for the short, medium, and tall women (p < 0.05; Fig. 3). The peak knee extension moment increased from the baseline condition for tall women while running with a load of 11.3 or 22.7 kg (p < 0.001) but did not change for short and medium women. Running with a 22.7-kg load significantly increased the peak ankle plantarflexion moment in short women [i.e., from 126.3 (SD = 21.3) to 134.9 (SD = 19.3) N·m, p < 0.01] and medium women [i.e., from 151.3 (SD = 17.0) to 167.5 (SD = 18.0) N·m, p < 0.001], but not in tall women (Fig. 3).

As expected, the peak GRFs for the baseline condition were lower than the peak GRFs when running with a 11.3-kg or 22.7-kg load in short, medium, and tall women (p < 0.001; Fig. 4). Similarly, when running with a load, the peak JRF at the hip (p < 0.05) and the knee (p < 0.01) were higher than their corresponding values in the baseline condition. Moreover, running with a 22.7-kg load increased the average peak ankle JRFs by 6.7% for the short group and 8.7% for the medium group (p < 0.001).

We separately evaluated the effect of body size on joint kinematics and kinetics for each of the three load conditions. For the baseline condition, peak hip flexion moment and ankle plantarflexion moment for tall women were higher than the corresponding moments in short women by 25.9 and 36.8%, respectively (p < 0.05, Table 2). While running with an additional load of 11.3 or 22.7 kg, we observed differences in peak ankle plantarflexion moment (p < 0.01) and peak knee JRF (p < 0.05) between short and tall women. The peak hip and knee JRFs for tall women were higher than those in short women for the baseline conditions (p < 0.01).

Table 2 Peak joint angles, joint moments, and joint reaction forces at the hip, knee, and ankle when subjects ran without a load (0.0 kg), or with an additional load of 11.3 kg (25 lb) or 22.7 kg (50 lb), averaged within the group of short, medium, and tall women. The data are presented as mean ± one standard deviation

Tibial strains and stress-fracture risk

Figure 5 shows the predicted 10th, 50th, and 90th percentile peak tibial von Mises strain during the entire gait cycle in the tibial bone. The 10th percentile strain for the baseline condition was lower than the strain for the 11.3-kg and 22.7-kg load conditions for each of the three groups (p < 0.05). Compared to the baseline condition, an additional 22.7 kg significantly increased the 90th percentile von Mises strain in short women from 7183 (SD = 2069) to 7726 (SD = 2435) με (i.e., an increase of 7.6%) in short women and from 7531 (SD = 952) to 7968 (SD = 754) με (i.e., an increase of 5.8%) in medium women (p < 0.05; Fig. 5), but not in tall women. We did not observe any influence of body size on the tibial strain while running without or with a load of 11.3 and 22.7 kg.

Fig. 5
figure 5

Comparison of von Mises strains (10th, 50th, and 90th percentile) under different loads for the short, medium, and tall groups during one gait cycle. Data are expressed as mean ± one standard deviation. *p < 0.05, **p < 0.01, ***p < 0.001

When compared to the baseline condition, an additional load of 22.7 kg increased the likelihood of stress fracture in short and medium women by 9.7 and 7.4%, respectively (p < 0.001; Fig. 6). However, similar to the effect of body size on the strains, we did not observe statistically significant differences in the likelihood of stress fracture between the short, medium, and tall groups in any of the three load conditions.

Fig. 6
figure 6

Stress-fracture likelihood Ps when subjects ran 4.8 km/day for 100 consecutive days under different loads for the short, medium, and tall groups (n = 7 in each group). *p < 0.05, ***p < 0.001

Discussion

In this study, we developed individualized M/FE models for 21 women of short (n = 7), medium (n = 7), and tall (n = 7) statures and analyzed the effects of body size and load carriage on the kinematics and kinetics of lower extremity joints, tibial strain, and likelihood of stress fracture when these women ran without a load or with a load of 11.3 or 22.7 kg. When compared to the baseline (i.e., no-load) condition, we observed that the relative changes in the joint mechanics (i.e., kinematics and kinetics), tibial strain, and the risk of stress fracture when running with a load were dependent on body size (Figs. 2, 3, 4, 5 and 6). In partial agreement with our hypothesis, body size influenced the joint kinetics in women running with or without a load, but not the tibial strain and the likelihood of stress fracture.

Relative to the baseline condition, peak hip flexion angle increased with increasing load in medium and tall women (p < 0.01) but not in short women (Fig. 2). In addition, when compared to the baseline condition, peak hip flexion moment and extension moment of women running with a load of 22.7 kg increased for all stature groups (p < 0.05; Fig. 3). These results, which are in agreement with previous studies on women running with a load [21, 34], show engagement of hip muscles while performing a strenuous physical activity. Running with a load changed the kinematics and kinetics of the knee differently in the short, medium, and tall stature groups. For example, when compared to the baseline condition, running with a load of 22.7 kg reduced knee flexion angle during the stance phase in short women (p < 0.05), but not in tall women (Fig. 2). In contrast, with increasing load, the knee extension moment increased in tall women (p < 0.001) but not in short women (Fig. 3). When compared to the baseline condition, the reduction in the knee flexion angle in short women (Fig. 2) might reduce the ability of these women to dissipate forces while running with a 22.7-kg load. Increased knee stiffness, either due to a reduction in the knee flexion angle and/or increased knee flexion moment, is also considered as a risk factor for tibial stress fracture [35, 36].

Similar to the results of joint kinematics and kinetics, the change in the tibial strain (Fig. 5) and the likelihood of stress fracture with load carriage (Fig. 6), when compared to the baseline condition, was dependent on body size. For example, when compared to the baseline condition, running with 22.7-kg load increased the 50th and 90th percentile strain in the short and medium groups but not in the tall group. Similar to the change in the tibial strains, the likelihood of stress fracture in women resulting from running 4.8 km/day at 3.0 m/s for 100 consecutive days, increased with the load carried by women in the short and medium groups (p < 0.001). However, the differences in the strain and the likelihood of stress fracture for short, medium, and tall women under each of the three load conditions were not statistically significant, in spite of our model accounting for the important relationship between tibial size and stature [37]. It is interesting to note that, when compared to the other groups, the greater peak hip JRF for the baseline condition and peak knee JRF in the tall group for all load conditions (p < 0.05; Table 2) did not translate to higher strains or an increased risk of stress fracture for this group. We suspect that this observation could be due to stature-related differences in the geometry and material properties in the tibia, with the latter resulting from bone adaptation in response to mechanical forces generated during normal daily activity. These results also highlight the need for developing subject-specific M/FE models and combining the results of individuals with similar statures to perform group-level analyses.

The study has limitations. First, similar to our previous models [19, 21], we represented knee and ankle joints as revolute joints and did not include complex 3-D motions, including translation and rotation, of the joints. While inclusion of such complex motions can increase the fidelity of the model, revolute joints can adequately capture the kinematics and kinetics of the knee and ankle joints [38]. Therefore, we do not believe that incorporation of more complex 3-D motions in the musculoskeletal model would have changed the conclusions regarding the joint kinematics and kinetics as well as stress-fracture risk. Second, in the musculoskeletal model, we did not represent individualized muscle strength. Third, in the stress-fracture predictions, we did not consider any variation in gait mechanics due to muscle fatigue. Such an assumption may under-predict the likelihood of stress fracture, as fatigue may increase tibial strains and, thereby, the risk of stress fracture [39]. In addition, the probabilistic model assumes uniform bone remodeling and adaptation for all subjects. Although these assumptions are unlikely to affect the relative differences in the likelihood of stress fracture for the different load conditions for a subject, incorporation of muscle fatigue in the musculoskeletal model and an improved model of bone remodeling and adaptation [40] should be considered in the future to enhance the accuracy of stress-fracture predictions. Finally, we only quantified the effects of body size and load carriage for level running on a treadmill, without considering graded tests.

Conclusion

In summary, using individualized M/FE model, we evaluated the effect of body size and load carriage on lower-extremity joint kinematics and kinetics, tibial biomechanics, and the likelihood of tibial stress fracture in 21 women running with or without a load. The results show that women of different statures adjust their gait mechanisms differently when running with a load. We also observed that in women running with or without a load, body size influences the joint kinetics, but not the prediction of tibial strains and the likelihood of stress fracture. These findings will help in quantifying the effect of load carriage on bone health and ultimately reduce the risk of musculoskeletal injuries of the lower extremities, such as knee injury and tibial stress fracture, in women.

Availability of data and materials

The FE-generated data and related analyses will be made available through a written request to the corresponding author, including a summary of the planned research.

Abbreviations

3-D:

Three-dimensional

BMI:

Body mass index

CT:

Computed tomography

GRF:

Ground reaction force

JRF:

Joint reaction force

M/FE:

Musculoskeletal-finite-element

SD:

Standard deviation

References

  1. Molloy JM, Pendergrass TL, Lee IE, Chervak MC, Hauret KG, Rhon DI. Musculoskeletal injuries and United States Army readiness part I: overview of injuries and their strategic impact. Mil Med. 2020;185(9-10):e1461–71.

  2. Milgrom C, Giladi M, Stein M, Kashtan H, Margulies JY, Chisin R, et al. Stress fractures in military recruits. A prospective study showing an unusually high incidence. J Bone Joint Surg Br. 1985;67(5):732–5.

    Article  CAS  Google Scholar 

  3. Mattila VM, Niva M, Kiuru M, Pihlajamaki H. Risk factors for bone stress injuries: a follow-up study of 102,515 person-years. Med Sci Sports Exerc. 2007;39(7):1061–6.

    Article  Google Scholar 

  4. Moran DS, Israeli E, Evans RK, Yanovich R, Constantini N, Shabshin N, et al. Prediction model for stress fracture in young female recruits during basic training. Med Sci Sports Exerc. 2008;40(11 Suppl):S636–44.

    Article  Google Scholar 

  5. Foot Marching, Load Carriage, and Injury Risk: Army Public Health Center; https://apps.dtic.mil/dtic/tr/fulltext/u2/1010939.pdf; Accessed 01/24/2021.

  6. Bulathsinhala L, Hughes JM, McKinnon CJ, Kardouni JR, Guerriere KI, Popp KL, et al. Risk of stress fracture varies by race/ethnic origin in a cohort study of 1.3 million US Army soldiers. J Bone Miner Res. 2017;32(7):1546–53.

    Article  CAS  Google Scholar 

  7. Winfield AC, Moore J, Bracker M, Johnson CW. Risk factors associated with stress reactions in female marines. Mil Med. 1997;162(10):698–702.

    Article  CAS  Google Scholar 

  8. Kelly EW, Jonson SR, Cohen ME, Shaffer R. Stress fractures of the pelvis in female navy recruits: an analysis of possible mechanisms of injury. Mil Med. 2000;165(2):142–6.

    Article  CAS  Google Scholar 

  9. Knapik J, Montain SJ, McGraw S, Grier T, Ely M, Jones BH. Stress fracture risk factors in basic combat training. Int J Sports Med. 2012;33(11):940–6.

    Article  CAS  Google Scholar 

  10. Valimaki VV, Alfthan H, Lehmuskallio E, Loyttyniemi E, Sahi T, Suominen H, et al. Risk factors for clinical stress fractures in male military recruits: a prospective cohort study. Bone. 2005;37(2):267–73.

    Article  Google Scholar 

  11. Shaffer RA, Rauh MJ, Brodine SK, Trone DW, Macera CA. Predictors of stress fracture susceptibility in young female recruits. Am J Sports Med. 2006;34(1):108–15.

    Article  Google Scholar 

  12. Sumnik Z, Land C, Coburger S, Neu C, Manz F, Hrach K, et al. The muscle-bone unit in adulthood: influence of sex, height, age and gynecological history on the bone mineral content and muscle cross-sectional area. J Musculoskelet Neuronal Interact. 2006;6(2):195–200.

    CAS  PubMed  Google Scholar 

  13. Maughan RJ, Watson JS, Weir J. Strength and cross-sectional area of human skeletal muscle. J Physiol. 1983;338:37–49.

    Article  CAS  Google Scholar 

  14. Taboga P, Lazzer S, Fessehatsion R, Agosti F, Sartorio A, di Prampero PE. Energetics and mechanics of running men: the influence of body mass. Eur J Appl Physiol. 2012;112(12):4027–33.

    Article  Google Scholar 

  15. Hora M, Sládek V, Soumar L. stráníková K, Michálek T. Influence of body mass and lower limb length on knee flexion angle during walking in humans. Folia Zool. 2012;61(3–4):330–9.

  16. Polcyn AF, Bensel C, Harman E, Obusek J, Pandorf C, Frykman P. Effects of weight carried by soldiers: Combined analysis of four studies on maximal performance, physiology, and biomechanics. Natick: Army Research Institute of Environmental Medicine; 2002. Report No.: NATICK/TR-02/010

    Book  Google Scholar 

  17. Attwells RL, Birrell SA, Hooper RH, Mansfield NJ. Influence of carrying heavy loads on soldiers' posture, movements and gait. Ergonomics. 2006;49(14):1527–37.

    Article  Google Scholar 

  18. Fox BD, Judge LW, Dickin DC, Wang H. Biomechanics of military load carriage and resulting musculoskeletal injury: a review. J Orthopedics Orthopedic Surg. 2020;1(1):6–11.

    Google Scholar 

  19. Baggaley M, Esposito M, Xu C, Unnikrishnan G, Reifman J, Edwards WB. Effects of load carriage on biomechanical variables associated with tibial stress fractures in running. Gait Posture. 2020;77:190–4.

    Article  Google Scholar 

  20. Gordon CC, Blackwell CL, Bradtmiller B, Parham JL, Barrientos P, Paquette SP, et al. 2012 Anthropometric survey of U.S. Army personnel: methods and summary statistics. Technical report. U.S. Army NATICK soldier research, Center DaE; 2014. Report No.: NATICK/TR-15/007.

  21. Xu C, Silder A, Zhang J, Reifman J, Unnikrishnan G. A cross-sectional study of the effects of load carriage on running characteristics and tibial mechanical stress: implications for stress-fracture injuries in women. BMC Musculoskelet Disord. 2017;18(1):125–36.

    Article  Google Scholar 

  22. Edwards WB, Taylor D, Rudolphi TJ, Gillette JC, Derrick TR. Effects of running speed on a probabilistic stress fracture model. Clin Biomech. 2010;25(4):372–7.

    Article  Google Scholar 

  23. Taylor D, Kuiper JH. The prediction of stress fractures using a 'stressed volume' concept. J Orthop Res. 2001;19(5):919–26.

    Article  CAS  Google Scholar 

  24. Taylor D, Casolari E, Bignardi C. Predicting stress fractures using a probabilistic model of damage, repair and adaptation. J Orthop Res. 2004;22(3):487–94.

    Article  Google Scholar 

  25. Knapik J, Reynolds K. Load Carriage in Military Operations: A Review of Historical, Physiological, Biomechanical, and Medical Aspects. Military Quantitative Physiology. 2012. https://ke.army.mil/bordeninstitute/other_pub/loadcarriagePDF.pdf. Accessed 24 Jan 2021.

  26. Horsman K. The Twente lower extremity model. Consistent dynamic simulation of the human locomotor apparatus. Enschede: University of Twente; 2007.

    Google Scholar 

  27. Andersen MS, Damsgaard M, MacWilliams B, Rasmussen J. A computationally efficient optimisation-based method for parameter identification of kinematically determinate and over-determinate biomechanical systems. Comp Methods Biomechan Biomed Engineer. 2010;13(2):171–83.

    Article  CAS  Google Scholar 

  28. Sangeux M, Polak J. A simple method to choose the most representative stride and detect outliers. Gait Posture. 2015;41(2):726–30.

    Article  Google Scholar 

  29. Taylor D. Fatigue of bone and bones: an analysis based on stressed volume. J Orthop Res. 1998;16(2):163–9.

    Article  CAS  Google Scholar 

  30. Carter DR, Caler WE, Spengler DM, Frankel VH. Fatigue behavior of adult cortical bone: the influence of mean strain and strain range. Acta Orthop Scand. 1981;52(5):481–90.

    Article  CAS  Google Scholar 

  31. Turner CH, Forwood MR, Rho JY, Yoshikawa T. Mechanical loading thresholds for lamellar and woven bone formation. J Bone Miner Res. 1994;9(1):87–97.

    Article  CAS  Google Scholar 

  32. R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018.

    Google Scholar 

  33. Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67(1):1–48.

    Article  Google Scholar 

  34. Loverro KL, Hasselquist L, Lewis CL. Females and males use different hip and knee mechanics in response to symmetric military-relevant loads. J Biomech. 2019;95:109280.

    Article  Google Scholar 

  35. Milner CE, Ferber R, Pollard CD, Hamill J, Davis IS. Biomechanical factors associated with tibial stress fracture in female runners. Med Sci Sports Exerc. 2006;38(2):323–8.

    Article  Google Scholar 

  36. Milner CE, Hamill J, Davis I. Are knee mechanics during early stance related to tibial stress fracture in runners? Clin Biomech. 2007;22(6):697–703.

    Article  Google Scholar 

  37. Duyar I, Pelin C. Body height estimation based on tibia length in different stature groups. Am J Phys Anthropol. 2003;122(1):23–7.

    Article  Google Scholar 

  38. Konrath JM, Karatsidis A, Schepers HM, Bellusci G, de Zee M, Andersen MS. Estimation of the knee adduction moment and joint contact force during daily living activities using inertial motion capture. Sensors (Basel). 2019;19(7):1681.

    Article  Google Scholar 

  39. Milgrom C, Radeva-Petrova DR, Finestone A, Nyska M, Mendelson S, Benjuya N, et al. The effect of muscle fatigue on in vivo tibial strains. J Biomech. 2007;40(4):845–50.

    Article  Google Scholar 

  40. Beaupré GS, Orr TE, Carter DR. An approach for time-dependent bone modeling and remodeling-application: a preliminary remodeling simulation. J Orthop Res. 1990;8(5):662–70.

    Article  Google Scholar 

Download references

Acknowledgments

We acknowledge the high-performance computing resources made available by the U.S. Department of Defense High Performance Computing Modernization Program at the U.S. Army Research Laboratory.

Funding

This research was sponsored by the U.S. Army Network Science Initiative and a grant from the Defense Health Program managed by the Military Operational Medicine Program Area Directorate, U.S. Army Medical Research and Development Command, Ft. Detrick, Maryland.

Author information

Authors and Affiliations

Authors

Contributions

GU, CX, MB, WBE, and JR designed the study. MB and WBE performed the experiments and analyzed the experimental data. CX developed the computational model and performed simulations. CX and GU analyzed the simulation results and drafted the manuscript. JR contributed to data evaluation and manuscript preparation. JT and SK provided feedback on data interpretation. All authors edited the manuscript and reviewed the final version. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Jaques Reifman.

Ethics declarations

Ethics approval and consent to participate

We received approval for the study from the Human Research Protection Office (Protocol Number: A17952.2a; Approved on 27 December 2016) at the United States (U.S.) Army Medical Research and Development Command (Ft. Detrick, MD) and the Conjoint Health Research Ethics Board (Protocol Number: REB16–1475; Approved on 9 September 2016) of the University of Calgary (Calgary, AB, Canada). We obtained informed written consent from each participant prior to the study.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests. The opinions and assertions contained herein are the private views of the authors and are not to be construed as official or as reflecting the views of the U.S. Army, the U.S. Department of Defense (DoD), or The Henry M. Jackson Foundation (HJF) for the Advancement of Military Medicine, Inc. Any citations of commercial organizations and trade names in this report do not constitute an official U.S. Army, DoD, or HJF endorsement of approval of the products or services of these organizations. This paper has been approved for public release with unlimited distribution.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Unnikrishnan, G., Xu, C., Baggaley, M. et al. Effects of body size and load carriage on lower-extremity biomechanical responses in healthy women. BMC Musculoskelet Disord 22, 219 (2021). https://doi.org/10.1186/s12891-021-04076-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12891-021-04076-0

Keywords