The three-dimensional determination of internal loads in the lower extremity

Generation of Customized Bone Implants from CT Scans Using FEA and AM

Additive manufacturing (AM) allows the creation of customized designs for various medical devices, such as implants, casts, and splints. Amongst other AM technologies, fused filament fabrication (FFF) facilitates the production of intricate geometries that are often unattainable through conventional methods like subtractive manufacturing. This study aimed to develop a methodology for substituting a pathological talus bone with a personalized one created using additive manufacturing. The process involved generating a numerical parametric solid model of the specific anatomical region using computed tomography (CT) scans of the corresponding healthy organ from the patient. The healthy talus served as a mirrored template to replace the defective one. Structural simulation of the model through finite element analysis (FEA) helped compare and select different materials to identify the most suitable one for the replacement bone. The implant was then produced using FFF technology. The developed procedure yielded commendable results. The models maintained high geometric accuracy, while significantly reducing the computational time. PEEK emerged as the optimal material for bone replacement among the considered options and several specimens of talus were successfully printed.

Lower extremity joint contact force symmetry during walking and running, 2-7 years post-ACL reconstruction

Premature osteoarthritis after anterior cruciate ligament reconstruction (ACLR) is common among athletes. Reduced knee contact forces after ACLR likely contribute to the multifactorial etiology of the disease. Whether this reduction is accompanied by compensatory increases in joint contact forces (JCF) at adjacent or contralateral joints is unclear. It is also unclear if compensatory effects depend on the task demands. Thus, we compared hip, knee, and ankle JCF symmetry between individuals with reconstruction and a matched control group during walking and running. Thirty participants (19 females), 2-7 years post-unilateral ACLR (mean = 47.8 months), and 30 controls matched on sex, mass, and activity level were recruited. Limb symmetry indices of peak contact forces and force impulses were calculated for each joint during walking and running, and analyzed using two-factor (group, activity) analysis of variances. Lower ACLR group peak knee JCF (p = 0.009) and knee JCF impulse (p = 0.034) during walking and running were observed. An interaction of group and activity was observed for peak hip JCF, with ACLR participants demonstrating greater involved limb peak hip JCF during running (p = 0.012). Ankle JCF and ground reaction force symmetry indices were not different between groups or across tasks. Decreased knee and increased ipsilateral peak hip JCF during running suggests that proximal adaptations exist at 2-7 years after ACLR, particularly during activities with increased task demand. Clinical significance: Knee and hip JCF asymmetry at 2-7 years after ACLR may underscore a need for clinical strategies and follow-up assessments to identify and target such outcomes.

Quantifying demands on the hamstrings during high-speed running: A systematic review and meta-analysis

INTRODUCTION

Hamstring strain injury (HSI) remains a performance, economic, and player availability burden in sport. High-speed running (HSR) is cited as a common mechanism for HSI. While evidence exists regarding the high physical demands on the hamstring muscles in HSR, meta-analytical synthesis of related activation and kinetic variables is lacking.

METHODS

A systematic search of Medline, Embase, Scopus, CINAHL, SportDiscus, and Cochrane library databases was conducted in accordance with the PRISMA 2020 guidelines. Studies reporting hamstring activation (electromyographic [EMG]) or hamstring muscle/related joint kinetics were included where healthy adult participants ran at or beyond 60% of maximum speed (activation studies) or 4 m per second (m/s) (kinetic studies).

RESULTS

A total of 96 studies met the inclusion criteria. Run intensities were categorized as "slow," "moderate," or "fast" in both activation and kinetic based studies with appropriate relative, and raw measures, respectively. Meta-analysis revealed pooled mean lateral hamstring muscle activation levels of 108.1% (95% CI: 84.4%-131.7%) of maximal voluntary isometric contraction (MVIC) during "fast" running. Meta-analysis found swing phase peak knee flexion internal moment and power at 2.2 Newton meters/kilogram (Nm/kg) (95% CI: 1.9-2.5) and 40.3 Watts/kilogram (W/kg) (95% CI: 31.4-49.2), respectively. Hip extension peak moment and power was estimated as 4.8 Nm/kg (95% CI: 3.9-5.7) and 33.1 W/kg (95% CI: 17.4-48.9), respectively.

CONCLUSIONS

As run intensity/speed increases, so do the activation and kinetic demands on the hamstrings. The presented data will enable clinicians to incorporate more objective measures into the design of injury prevention and return-to-play decision-making strategies.

Effect of gait types and external weight carrying strategies on the femoral neck strains during stair descent

Gait and weight carrying method may change the femoral neck load during stair descent. Applying specific gait and weight carrying methods may reduce the femoral neck load during stair descent, which may reduce hip pain, hip pain related falls and fall related fractures for the older population. The purpose of this study was to test the effect of different gait types (step-over-step v.s. step-by-step) and external weight carrying strategies (ipsilateral v.s. contralateral side) on the femoral neck load, discover which method could reduce the femoral neck load effectively. Seventeen healthy adults from 50 to 70 yrs old were recruited. The kinematic and kinetic analysis, musculoskeletal modelling method were used to estimate the joint and muscle loads for the lower extremities. Finite element analysis was used with the femur model to calculate the femoral neck strains during stair descent with different gait types and weight carrying strategies. The compressive strains were reduced for step-by-step gait method than step-over-step (p<0.015, 12.3-17.4% decrease of strains), the tensile strains were significantly increased for the trailing leg of step-by-step than the leading leg (p<0.001, 24.7% increase of strains). Contralateral weight carrying increased compressive and tensile strains than ipsilateral (p<0.001, 9.9-24.5% increase of strains) in most conditions. Applying step-by-step method and avoiding contralateral side weight carrying could be effective to reduce femoral neck strains. These outcomes could be helpful for the older population to reduce the risks of hip pain, femoral neck pain or pain related falls and fractures.

Peak tibial acceleration should not be used as indicator of tibial bone loading during running

Peak tibial acceleration (PTA) is a widely used indicator of tibial bone loading. Indirect bone loading measures are of interest to reduce the risk of stress fractures during running. However, tibial compressive forces are caused by both internal muscle forces and external ground reaction forces. PTA might reflect forces from outside the body, but likely not the compressive force from muscles on the tibial bone. Hence, the strength of the relationship between PTA and maximum tibial compression forces in rearfoot-striking runners was investigated. Twelve runners ran on an instrumented treadmill while tibial acceleration was captured with accelerometers. Force plate and inertial measurement unit data were spatially aligned with a novel method based on the centre of pressure crossing a virtual toe marker. The correlation coefficient between maximum tibial compression forces and PTA was 0.04 ± 0.14 with a range of -0.15 to +0.28. This study showed a very weak and non-significant correlation between PTA and maximum tibial compression forces while running on a level treadmill at a single speed. Hence, PTA as an indicator for tibial bone loading should be reconsidered, as PTA does not provide a complete picture of both internal and external compressive forces on the tibial bone.  .

Independently ambulatory children with spina bifida experience near-typical knee and ankle joint moments and forces during walking

BACKGROUND

Spina bifida, a neurological defect, can result in lower-limb muscle weakness. Altered ambulation and reduced musculoskeletal loading can yield decreased bone strength in individuals with spina bifida, yet individuals who remain ambulatory can exhibit normal bone outcomes.

RESEARCH QUESTION

During walking, how do lower-limb joint kinematics and moments and tibial forces in independently ambulatory children with spina bifida differ from those of children with typical development?

METHODS

We retrospectively analyzed data from 16 independently ambulatory children with spina bifida and 16 children with typical development and confirmed that tibial bone strength was similar between the two groups. Plantar flexor muscle strength was measured by manual muscle testing, and 14 of the children with spina bifida wore activity monitors for an average of 5 days. We estimated tibial forces at the knee and ankle using motion capture data and musculoskeletal simulations. We used Statistical Parametric Mapping t-tests to compare lower-limb joint kinematic and kinetic waveforms between the groups with spina bifida and typical development. Within the group with spina bifida, we examined relationships between plantar flexor muscle strength and peak tibial forces by calculating Spearman correlations.

RESULTS

Activity monitors from the children with spina bifida reported typical daily steps (9656 [SD 3095]). Despite slower walking speeds (p = 0.004) and altered lower-body kinematics (p < 0.001), children with spina bifida had knee and ankle joint moments and forces similar to those of children with typical development, with no detectable differences during stance. Plantar flexor muscle weakness was associated with increased compressive knee force (p = 0.002) and shear ankle force (p = 0.009).

SIGNIFICANCE

High-functioning, independently ambulatory children with spina bifida exhibited near-typical tibial bone strength and near-typical step counts and tibial load magnitudes. Our results suggest that the tibial forces in this group are of sufficient magnitudes to support the development of normal tibial bone strength.

Contributions of individual muscle forces to hip, knee, and ankle contact forces during the stance phase of running: a model-based study

Knowledge of muscle forces' contributions to the joint contact forces can assist in the evaluation of muscle function, joint injury prevention, treatment of gait disorders, and arthroplasty planning. This study's objective was to evaluate the contributions of human lower limb muscles to the hip, knee, and ankle joint contact forces during the stance phase of running. A total of 25 muscles (or groups) were investigated based on the OpenSim framework along the anterior-posterior, superoinferior, and mediolateral components of each joint coordinate system. It was revealed that, during the running stance phase, the gluteus medius, gluteus maximus, and iliopsoas mainly contributed to the hip contact force. The soleus, vastus group, and rectus femoris primarily contributed to the knee contact force, while the peroneus, soleus, gluteus medius, and gastrocnemius mainly contributed to the ankle joint force; some muscles simultaneously offloaded the joints during the stance phase. The distributive pattern of the individual muscle functions contributing to the joint load may substantially differ during the running and walking stance phases. This study's findings may further provide suggestive information for the design of lower limb joint prosthesis, the study of the biomechanics of pathologic walking and running, and the progression of joint osteoarthritis.

Differences in Muscle Demand and Joint Contact Forces Between Running and Skipping

Skipping has been proposed as a viable cross-training exercise to running due to its lower knee contact forces and higher whole-body energy expenditure. However, how individual muscle forces, energy expenditure, and joint loading are affected by differences in running and skipping mechanics remains unclear. The purpose of this study was to compare individual muscle forces, energy expenditure, and lower extremity joint contact forces between running and skipping using musculoskeletal modeling and simulations of young adults (n = 5) performing running and skipping at 2.5 m·s−1 on an instrumented treadmill. In agreement with previous work, running had greater knee and patella contact forces than skipping which was accompanied by greater knee extensor energetic demand. Conversely, skipping had greater ankle contact forces and required greater energetic demand from the uniarticular ankle plantarflexors. There were no differences in hip contact forces between gaits. These findings further support skipping as a viable alternative to running if the primary goal is to reduce joint loading at the commonly injured patellofemoral joint. However, for those with ankle injuries, skipping may not be a viable alternative due to the increased ankle loads. These findings may help clinicians prescribe activities most appropriate for a patient’s individual training or rehabilitation goals.

Peak and Per-Step Tibial Bone Stress During Walking and Running in Female and Male Recreational Runners

BACKGROUND

Athletes, especially female athletes, experience high rates of tibial bone stress injuries (BSIs). Knowledge of tibial loads during walking and running is needed to understand injury mechanisms and design safe running progression programs.

PURPOSE

To examine tibial loads as a function of gait speed in male and female runners.

STUDY DESIGN

Controlled laboratory study.

METHODS

Kinematic and kinetic data were collected on 40 recreational runners (20 female, 20 male) during 4 instrumented gait speed conditions on a treadmill (walk, preferred run, slow run, fast run). Musculoskeletal modeling, using participant-specific magnetic resonance imaging and motion data, was used to estimate tibial stress. Peak tibial stress and stress-time impulse were analyzed using 2-factor multivariate analyses of variance (speed*sex) and post hoc comparisons (α = .05). Bone geometry and tibial forces and moments were examined.

RESULTS

Peak compression was influenced by speed (P < .001); increasing speed generally increased tibial compression in both sexes. Women displayed greater increases in peak tension (P = .001) and shear (P < .001) than men when transitioning from walking to running. Further, women displayed greater peak tibial stress overall (P < .001). Compressive and tensile stress-time impulse varied by speed (P < .001) and sex (P = .006); impulse was lower during running than walking and greater in women. A shear stress-time impulse interaction (P < .001) indicated that women displayed greater impulse relative to men when changing from a walk to a run. Compared with men, women displayed smaller tibiae (P < .001) and disproportionately lower tibial forces (P≤ .001-.035).

CONCLUSION

Peak tibial stress increased with gait speed, with a 2-fold increase in running relative to walking. Women displayed greater tibial stress than men and greater increases in stress when shifting from walking to running. Sex differences appear to be the result of smaller bone geometry in women and tibial forces that were not proportionately lower, given the womens' smaller stature and lower mass relative to men.

CLINICAL RELEVANCE

These results may inform interventions to regulate running-related training loads and highlight a need to increase bone strength in women. Lower relative bone strength in women may contribute to a sex bias in tibial BSIs, and female runners may benefit from a slower progression when initiating a running program.

Finite element analysis of femoral neck strains during stair ascent and descent

For older population, a better understanding of the hip joint loading environment is needed for the prevention of hip pain, and the reduction of the stress fractures and fall risks. Using the motion analysis and inverse dynamics methods, combined with musculoskeletal modelling, static optimization, and finite element (FE) femur model, the difference of femoral neck strains between stair ascent vs. descent, young vs. older populations was compared. A two-way repeated-measures MANOVA was applied to test the effect of age and stair direction on the femoral neck strains. The strains at the femoral neck cross-section were greater for stair descent than ascent for both age groups (mostly P = 0.001 to 0006) but there was no difference between age groups. In this study, femoral neck strains represented bone loading environment in more direct ways than joint reaction forces/moments or joint contact forces, the risk of hip pain, falls and stress fractures might be greater during stair descent than ascent. Possible preventative methods to reduce these risks should be developed in the future studies.

Determination of the correlation between muscle forces obtained from OpenSim and muscle activities obtained from electromyography in the elderly

Measurement of muscle forces related to aging can help to better identify the gait impairment mechanisms in the elderly. To this end, musculoskeletal modeling has been developed to estimate muscle forces. This study aimed to check the validity of OpenSim modeling (i.e., computed muscle control) approach in elderly subjects. Kinematic and kinetic data and Electromyography (EMG) signals for four different muscles were collected in nine healthy elderly males during walking. Dynamic simulation was done within OpenSim. Correlation analysis was performed to quantitatively compare the maximum estimated muscle forces with maximum measured muscle activities during the first double limb support, single limb support, and the second double limb support phases. The area-time plots of OpenSim and EMG data during gait cycle were obtained for qualitative assessment. In quantitative assessment, a low to moderate correlation was observed for the peak of muscle force and muscle activation of four muscles during sub phases of gait. The muscle forces pattern from OpenSim was found to be relatively similar to the muscle activity pattern from EMG especially for Gastrocnemius Medialis. A low to moderate consistency between OpenSim and EMG in the elderly can be explained by using a single mathematical estimation approach.

Surface acceleration transmission during drop landings in humans

The purpose of this study was to quantify the magnitude and frequency content of surface-measured accelerations at each major human body segment from foot to head during impact landings. Twelve males performed two single leg drop landings from each of 0.15 m, 0.30 m, and 0.45 m. Triaxial accelerometers (2000 Hz) were positioned over the: first metatarsophalangeal joint; distal anteromedial tibia; superior to the medial femoral condyle; L5 vertebra; and C6 vertebra. Analysis of acceleration signal power spectral densities revealed two distinct components, 2-14 Hz and 14-58 Hz, which were assumed to correspond to time domain signal joint rotations and elastic wave tissue deformation, respectively. Between each accelerometer position from the metatarsophalangeal joint to the L5 vertebra, signals exhibited decreased peak acceleration, increased time to peak acceleration, and decreased power spectral density integral of both the 2-14 Hz and 14-58 Hz components, with no further attenuation beyond the L5 vertebra. This resulted in peak accelerations close to vital organs of less than 10% of those at the foot. Following landings from greater heights, peak accelerations measured distally were greater, as was attenuation prior to the L5 position. Active and passive mechanisms within the lower limb therefore contribute to progressive attenuation of accelerations, preventing excessive accelerations from reaching the torso and head, even when distal accelerations are large.

Measuring femoral neck loads in healthy young and older adults during stair ascent and descent

Understanding the hip loading environment for daily activities is useful for hip fracture prevention, rehabilitation, and the design of osteogenic exercises. Seventeen older adults (50–70 yrs) and twenty young adults (18–30 yrs) were recruited. A rigid body model combined with a musculoskeletal model was used to estimate lower extremity loading. An elliptical cross-section model of the femoral neck was used to estimate femoral neck stress during stair ascent and descent. Two peaks were identified in the stress curves, corresponding to the peaks in the vertical ground reaction force. During stair ascent, significantly higher tension on the superior femoral neck was found for the young group at peak 1 (young: 13.5±6.1 MPa, older: 4.2±6.5 MPa, p<0.001). Also during stair ascent, significantly higher compression on the posterior femoral neck was found for the older group at peak 2 (young: -11.4±4.9 MPa, old: -18.1±8.6 MPa, p = 0.006). No significant difference was found for stair descent. Components of stress (muscle vs. reaction forces; axial forces vs. bending moments) were also examined for each trial of stair ascent and descent. The stresses and their components provided loading magnitude and locations of higher stress on the femoral neck during stair ascent and descent. Understanding femoral neck stresses may be used to help prevent hip fractures, reduce pain, improve rehabilitation, and design osteogenic exercises.

Biomechanical evaluation of total ankle arthroplasty. Part I: Joint loads during simulated level walking

In total ankle arthroplasty, the interaction at the joint between implant and bone is driven by a complex loading environment. Unfortunately, little is known about the loads at the ankle during daily activities since earlier attempts use two- or three-dimensional models to explore simplified joint mechanics. Our goal was to develop a framework to calculate multi-axial loads at the joint during simulated level walking following total ankle arthroplasty. To accomplish this, we combined robotic simulations of level walking at one-quarter bodyweight in three cadaveric foot and ankle specimens with musculoskeletal modeling to calculate the multi-axial forces and moments at the ankle during the stance phase. The peak compressive forces calculated were between 720 and 873 N occurring around 77%-80% of stance. The peak moment, which was the internal moment for all specimens, was between 6.1 and 11.6 N m and occurred between 72% and 88% of the stance phase. The peak moment did not necessarily occur with the peak force. The ankle joint loads calculated in this study correspond well to previous attempts in the literature; however, our robotic simulator and framework provide an opportunity to resolve the resultant three-dimensional forces and moments as others have not in previous studies. The framework may be useful to calculate ankle joint loads in cadaveric specimens as the first step in evaluating bone-implant interactions in total ankle replacement using specimen specific inputs. This approach also provides a unique opportunity to evaluate changes in joint loads and kinematics following surgical interventions of the foot and ankle.

Allowed Activities After Primary Total Knee Arthroplasty and Total Hip Arthroplasty

It is essential for total knee arthroplasty patients to return to their previous level of activity to maintain a healthy lifestyle. This article reviews the current recommendations regarding return to physical activity after total knee arthroplasty and trying to find the balance between levels of activity and prosthetic joint preservation. In general, most total joint replacement patients are able to return to their previous level of activity and to a lesser extent to sports. This article discuss patients' actual levels of activity including their return to work and sport and the factors that influence meeting their expectations for surgery.

Subject-specific muscle properties from diffusion tensor imaging significantly improve the accuracy of musculoskeletal models

Musculoskeletal modelling is an important platform on which to study the biomechanics of morphological structures in vertebrates and is widely used in clinical, zoological and palaeontological fields. The popularity of this approach stems from the potential to non-invasively quantify biologically important but difficult-to-measure functional parameters. However, while it is known that model predictions are highly sensitive to input values, it is standard practice to build models by combining musculoskeletal data from different sources resulting in 'generic' models for a given species. At present, there are little quantitative data on how merging disparate anatomical data in models impacts the accuracy of these functional predictions. This issue is addressed herein by quantifying the accuracy of both subject-specific human limb models containing individualised muscle force-generating properties and models built using generic properties from both elderly and young individuals, relative to experimental muscle torques obtained from an isokinetic dynamometer. The results show that subject-specific models predict isokinetic muscle torques to a greater degree of accuracy than generic models at the ankle (root-mean-squared error - 7.9% vs. 49.3% in elderly anatomy-based models), knee (13.2% vs. 57.3%) and hip (21.9% vs. 32.8%). These results have important implications for the choice of musculoskeletal properties in future modelling studies, and the relatively high level of accuracy achieved in the subject-specific models suggests that such models can potentially address questions about inter-subject variations of muscle functions. However, despite relatively high levels of overall accuracy, models built using averaged generic muscle architecture data from young, healthy individuals may lack the resolution and accuracy required to study such differences between individuals, at least in certain circumstances. The results do not wholly discourage the continued use of averaged generic data in musculoskeletal modelling studies but do emphasise the need for to maximise the accuracy of input values if studying intra-species form-function relationships in the musculoskeletal system.

Multibody Kinematics Optimization for the Estimation of Upper and Lower Limb Human Joint Kinematics: A Systematized Methodological Review

Multibody kinematics optimization (MKO) aims to reduce soft tissue artefact (STA) and is a key step in musculoskeletal modeling. The objective of this review was to identify the numerical methods, their validation and performance for the estimation of the human joint kinematics using MKO. Seventy-four papers were extracted from a systematized search in five databases and cross-referencing. Model-derived kinematics were obtained using either constrained optimization or Kalman filtering to minimize the difference between measured (i.e., by skin markers, electromagnetic or inertial sensors) and model-derived positions and/or orientations. While hinge, universal, and spherical joints prevail, advanced models (e.g., parallel and four-bar mechanisms, elastic joint) have been introduced, mainly for the knee and shoulder joints. Models and methods were evaluated using: (i) simulated data based, however, on oversimplified STA and joint models; (ii) reconstruction residual errors, ranging from 4 mm to 40 mm; (iii) sensitivity analyses which highlighted the effect (up to 36 deg and 12 mm) of model geometrical parameters, joint models, and computational methods; (iv) comparison with other approaches (i.e., single body kinematics optimization and nonoptimized kinematics); (v) repeatability studies that showed low intra- and inter-observer variability; and (vi) validation against ground-truth bone kinematics (with errors between 1 deg and 22 deg for tibiofemoral rotations and between 3 deg and 10 deg for glenohumeral rotations). Moreover, MKO was applied to various movements (e.g., walking, running, arm elevation). Additional validations, especially for the upper limb, should be undertaken and we recommend a more systematic approach for the evaluation of MKO. In addition, further model development, scaling, and personalization methods are required to better estimate the secondary degrees-of-freedom (DoF).

Model-Based Estimation of Individual Muscle Force Based on Measurements of Muscle Activity in Forearm Muscles During Isometric Tasks

OBJECTIVE

Several forward dynamics estimation approaches have been proposed to estimate individual muscle force. However, characterization of the estimation error that arises when measurements are available only from a subset of the muscles involved in the movement under analysis, as is the case of the forearm muscles, has been limited. Our objectives were: first, to quantify the accuracy of forward-dynamics muscle force estimators for forearm muscles; and second, to develop a muscle force estimator that is accurate even when measurements are available only from a subset of muscles acting on a given joint or segment.

METHODS

We developed a neuromusculoskeletal (NMSK) estimator that integrates forward dynamics estimation with a neural model of muscle cocontraction to estimate individual muscle force during isometric contractions, suitable to operate when measurements are not available for all muscles. We developed a computational framework to assess the effect of physiological variability in muscle cocontraction, cross-talk, and measurement error on the estimator accuracy using a sensitivity analysis. We thus compared the performance of our estimator with that of a standard estimator that neglects the contribution of unmeasured muscles.

RESULTS

The NMSK estimator reduces the estimation error by 25% in average noise conditions. Moreover, the NMSK estimator is robust against physiological variability in muscle cocontraction and outperforms the standard estimator even when the validity of the neural model is compromised.

CONCLUSION AND SIGNIFICANCE

In isometric tasks, the NMSK estimator reduces muscle force estimation error compared to a standard estimator, and may enable future applications involving estimation of forearm muscle force during coordinated movements.

A 2-Year Prospective Cohort Study of Overuse Running Injuries: The Runners and Injury Longitudinal Study (TRAILS)

BACKGROUND

The National Center for Injury Prevention and Control, noting flaws in previous running injury research, called for more rigorous prospective designs and comprehensive analyses to define the origin of running injuries.

PURPOSE

To determine the risk factors that differentiate recreational runners who remain uninjured from those diagnosed with an overuse running injury during a 2-year observational period.

STUDY DESIGN

Cohort study; Level of evidence, 2.

METHODS

Inclusion criteria were running a minimum of 5 miles per week and being injury free for at least the past 6 months. Data were collected at baseline on training, medical and injury histories, demographics, anthropometrics, strength, gait biomechanics, and psychosocial variables. Injuries occurring over the 2-year observation period were diagnosed by an orthopaedic surgeon on the basis of predetermined definitions.

RESULTS

Of the 300 runners who entered the study, 199 (66%) sustained at least 1 injury, including 73% of women and 62% of men. Of the injured runners, 111 (56%) sustained injuries more than once. In bivariate analyses, significant ( P ≤ .05) factors at baseline that predicted injury were as follows: Short Form Health Survey-12 mental component score (lower mental health-related quality of life), Positive and Negative Affect Scale negative affect score (more negative emotions), sex (higher percentage of women were injured), and knee stiffness (greater stiffness was associated with injury); subsequently, knee stiffness was the lone significant predictor of injury (odds ratio = 1.18) in a multivariable analysis. Flexibility, quadriceps angle, arch height, rearfoot motion, strength, footwear, and previous injury were not significant risk factors for injury.

CONCLUSION

The results of this study indicate the following: (1) among recreational runners, women sustain injuries at a higher rate than men; (2) greater knee stiffness, more common in runners with higher body weights (≥80 kg), significantly increases the odds of sustaining an overuse running injury; and (3) contrary to several long-held beliefs, flexibility, arch height, quadriceps angle, rearfoot motion, lower extremity strength, weekly mileage, footwear, and previous injury are not significant etiologic factors across all overuse running injuries.

Treadmill Running Ameliorates Destruction of Articular Cartilage and Subchondral Bone, Not Only Synovitis, in a Rheumatoid Arthritis Rat Model

We analyzed the influence of treadmill running on rheumatoid arthritis (RA) joints using a collagen-induced arthritis (CIA) rat model. Eight-week-old male Dark Agouti rats were randomly divided into four groups: The control group, treadmill group (30 min/day for 4 weeks from 10-weeks-old), CIA group (induced CIA at 8-weeks-old), and CIA + treadmill group. Destruction of the ankle joint was evaluated by histological analyses. Morphological changes of subchondral bone were analyzed by μ-CT. CIA treatment-induced synovial membrane invasion, articular cartilage destruction, and bone erosion. Treadmill running improved these changes. The synovial membrane in CIA rats produced a large amount of tumor necrosis factor-α and Connexin 43; production was significantly suppressed by treadmill running. On μ-CT of the talus, bone volume fraction (BV/TV) was significantly decreased in the CIA group. Marrow star volume (MSV), an index of bone loss, was significantly increased. These changes were significantly improved by treadmill running. Bone destruction in the talus was significantly increased with CIA and was suppressed by treadmill running. On tartrate-resistant acid phosphate and alkaline phosphatase (TRAP/ALP) staining, the number of osteoclasts around the pannus was decreased by treadmill running. These findings indicate that treadmill running in CIA rats inhibited synovial hyperplasia and joint destruction.

Femoral Neck Stress in Older Adults During Stair Ascent and Descent

A detailed understanding of the hip loading environment is needed to help prevent hip fractures, minimize hip pain, rehabilitate hip injuries, and design osteogenic exercises for the hip. The purpose of this study was to compare femoral neck stress during stair ascent and descent and to identify the contribution of muscles and reaction forces to the stress environment in mature adult subjects (n = 17; age: 50-65 y). Motion analysis and inverse dynamics were combined with musculoskeletal modeling and optimization, then used as input to an elliptical femoral neck cross-sectional model to estimate femoral neck stress. Peak stress values at the 2 peaks of the bimodal stress curves (stress vs time plot) were compared between stair ascent and descent. Stair ascent had greater compressive stress than descent during the first peak at the anterior (ascent: -18.0 [7.9] MPa, descent: -12.9 [5.4] MPa, P < .001) and posterior (ascent: -34.4 [10.9] MPa, descent: -27.8 [10.1] MPa, P < .001) aspects of the femoral neck cross section. Stair descent had greater tensile stress during both peaks at the superior aspect (ascent: 1.3 [7.0] MPa, descent: 24.8 [9.7] MPa, peak 1: P < .001; ascent: 15.7 [6.1] MPa, descent: 18.0 [8.4] MPa, peak 2: P = .03) and greater compressive stress during the second peak at the inferior aspect (ascent: -43.8 [9.7] MPa, descent: -51.1 [14.3] MPa, P = .004). Understanding this information can provide a more comprehensive view of bone loading at the femoral neck for older population.

A systematic review of approaches to modelling lower limb muscle forces during gait: Applicability to clinical gait analyses

Computational methods to estimate muscle forces during walking are becoming more common in biomechanical research but not yet in clinical gait analysis. This systematic review aims to identify the current state-of-the-art, examine the differences between approaches, and consider applicability of the current approaches in clinical gait analysis. A systematic database search identified studies including estimated muscle force profiles of the lower limb during healthy walking. These were rated for quality and the muscle force profiles digitised for comparison. From 13.449 identified studies, 22 were finally included which used four modelling approaches: static optimisation, enhanced static optimisation, forward dynamics and EMG-driven. These used a range of different musculoskeletal models, muscle-tendon characteristics and cost functions. There is visually broad agreement between and within approaches about when muscles are active throughout the gait cycle. There remain, considerable differences (CV 7%-151%, range of timing of peak forces in gait cycle 1%-31%) in patterns and magnitudes of force between and within modelling approaches. The main source of this variability is not clear. Different musculoskeletal models, experimental protocols, and modelling approaches will clearly have an effect as will the variability of joint kinetics between healthy individuals. Limited validation of modelling approaches, particularly at the level of individual participants, makes it difficult to conclude if any of the approaches give consistently better estimates than others. While muscle force modelling has clear potential to enhance clinical gait analyses future research is needed to improve validation, accuracy and feasibility of implementation in clinical practice.

Vacuum level effects on knee contact force for unilateral transtibial amputees with elevated vacuum suspension

The elevated vacuum suspension system (EVSS) has demonstrated unique health benefits for amputees, but the effect of vacuum pressure values on knee contact force (KCF) is still unclear. The objective of this study was to investigate the effect of vacuum levels on KCF for unilateral transtibial amputees (UTA) using the EVSS. Three-dimensional gait was modeled for 9 UTA with five vacuum levels (0-20inHg [67.73kPa], 5inHg [16.93kPa] increments) and 9 non-amputees based on kinematic and ground reaction force data. The results showed that the vacuum level effects were significant for peak axial KCF, which had a relatively large value at 0 and 20inHg (67.73kPa). The intact limb exhibited a comparable peak axial KCF to the non-amputees at 15inHg (50.79kPa). At moderate vacuum levels (5inHg [16.93kPa] to 15inHg [50.79kPa]), co-contraction of quadriceps and hamstrings at peak axial KCF was similar for the intact limb, but was smaller for the residual limb comparing with the non-amputees. The intact limb showed a similar magnitude of quadriceps and hamstrings force at 15inHg (50.79kPa) to the non-amputees, but the muscle coordination patterns varied between the residual and intact limbs. These findings indicate that a proper vacuum level may partially compensate for the lack of ankle plantarflexor and reduce the knee loading. Of the tested vacuum levels, 15inHg (50.79kPa) appears most favorable, although additional analyses with more amputees are suggested to confirm these results prior to establishing clinical guidelines.

Lower extremity joint loads in habitual rearfoot and mid/forefoot strike runners with normal and shortened stride lengths

Our purpose was to compare joint loads between habitual rearfoot (hRF) and habitual mid/forefoot strikers (hFF), rearfoot (RFS) and mid/forefoot strike (FFS) patterns, and shorter stride lengths (SLs). Thirty-eight hRF and hFF ran at their normal SL, 5% and 10% shorter, as well as with the opposite foot strike. Three-dimensional ankle, knee, patellofemoral (PF) and hip contact forces were calculated. Nearly all contact forces decreased with a shorter SL (1.2-14.9% relative to preferred SL). In general, hRF had higher PF (hRF-RFS: 10.8 ± 1.4, hFF-FFS: 9.9 ± 2.0 BWs) and hip loads (axial hRF-RFS: -9.9 ± 0.9, hFF-FFS: -9.6 ± 1.0 BWs) than hFF. Many loads were similar between foot strike styles for the two groups, including axial and lateral hip, PF, posterior knee and shear ankle contact forces. Lateral knee and posterior hip contact forces were greater for RFS, and axial ankle and knee contact forces were greater for FFS. The tibia may be under greater loading with a FFS because of these greater axial forces. Summarising, a particular foot strike style does not universally decrease joint contact forces. However, shortening one's SL 10% decreased nearly all lower extremity contact forces, so it may hold potential to decrease overuse injuries associated with excessive joint loads.

Model-Based Estimation of Ankle Joint Stiffness

We address the estimation of biomechanical parameters with wearable measurement technologies. In particular, we focus on the estimation of sagittal plane ankle joint stiffness in dorsiflexion/plantar flexion. For this estimation, a novel nonlinear biomechanical model of the lower leg was formulated that is driven by electromyographic signals. The model incorporates a two-dimensional kinematic description in the sagittal plane for the calculation of muscle lever arms and torques. To reduce estimation errors due to model uncertainties, a filtering algorithm is necessary that employs segmental orientation sensor measurements. Because of the model's inherent nonlinearities and nonsmooth dynamics, a square-root cubature Kalman filter was developed. The performance of the novel estimation approach was evaluated in silico and in an experimental procedure. The experimental study was conducted with body-worn sensors and a test-bench that was specifically designed to obtain reference angle and torque measurements for a single joint. Results show that the filter is able to reconstruct joint angle positions, velocities and torque, as well as, joint stiffness during experimental test bench movements.

Immediate effects of modified landing pattern on a probabilistic tibial stress fracture model in runners

BACKGROUND

Tibial stress fracture is a common injury in runners. This condition has been associated with increased impact loading. Since vertical loading rates are related to the landing pattern, many heelstrike runners attempt to modify their footfalls for a lower risk of tibial stress fracture. Such effect of modified landing pattern remains unknown. This study examined the immediate effects of landing pattern modification on the probability of tibial stress fracture.

METHODS

Fourteen experienced heelstrike runners ran on an instrumented treadmill and they were given augmented feedback for landing pattern switch. We measured their running kinematics and kinetics during different landing patterns. Ankle joint contact force and peak tibial strains were estimated using computational models. We used an established mathematical model to determine the effect of landing pattern on stress fracture probability.

FINDINGS

Heelstrike runners experienced greater impact loading immediately after landing pattern switch (P<0.004). There was an increase in the longitudinal ankle joint contact force when they landed with forefoot (P=0.003). However, there was no significant difference in both peak tibial strains and the risk of tibial stress fracture in runners with different landing patterns (P>0.986).

INTERPRETATION

Immediate transitioning of the landing pattern in heelstrike runners may not offer timely protection against tibial stress fracture, despite a reduction of impact loading. Long-term effects of landing pattern switch remains unknown.

Bone stress in runners with tibial stress fracture

BACKGROUND

Combinations of smaller bone geometry and greater applied loads may contribute to tibial stress fracture. We examined tibial bone stress, accounting for geometry and applied loads, in runners with stress fracture.

METHODS

23 runners with a history of tibial stress fracture & 23 matched controls ran over a force platform while 3-D kinematic and kinetic data were collected. An elliptical model of the distal 1/3 tibia cross section was used to estimate stress at 4 locations (anterior, posterior, medial and lateral). Inner and outer radii for the model were obtained from 2 planar x-ray images. Bone stress differences were assessed using two-factor ANOVA (α=0.05). Key contributors to observed stress differences between groups were examined using stepwise regression.

FINDINGS

Runners with tibial stress fracture experienced greater anterior tension and posterior compression at the distal tibia. Location, but not group, differences in shear stress were observed. Stepwise regression revealed that anterior-posterior outer diameter of the tibia and the sagittal plane bending moment explained >80% of the variance in anterior and posterior bone stress.

INTERPRETATION

Runners with tibial stress fracture displayed greater stress anteriorly and posteriorly at the distal tibia. Elevated tibial stress was associated with smaller bone geometry and greater bending moments about the medial-lateral axis of the tibia. Future research needs to identify key running mechanics associated with the sagittal plane bending moment at the distal tibia as well as to identify ways to improve bone geometry in runners in order to better guide preventative and rehabilitative efforts.

Effect of lower-limb joint models on subject-specific musculoskeletal models and simulations of daily motor activities

Understanding the validity of using musculoskeletal models is critical, making important to assess how model parameters affect predictions. In particular, assumptions on joint models can affect predictions from simulations of movement, and the identification of image-based joints is unavoidably affected by uncertainty that can decrease the benefits of increasing model complexity. We evaluated the effect of different lower-limb joint models on muscle and joint contact forces during four motor tasks, and assessed the sensitivity to the uncertainties in the identification of anatomical four-bar-linkage joints. Three MRI-based musculoskeletal models having different knee and ankle joint models were created and used for the purpose. Model predictions were compared against a baseline model including simpler and widely-adopted joints. In addition, a probabilistic analysis was performed by perturbing four-bar-linkage joint parameters according to their uncertainty. The differences between models depended on the motor task analyzed, and there could be marked differences at peak loading (up to 2.40 BW at the knee and 1.54 BW at the ankle), although they were rather small over the motor task cycles (up to 0.59 BW at the knee and 0.31 BW at the ankle). The model including more degrees of freedom showed more discrepancies in predicted muscle activations compared to measured muscle activity. Further, including image-based four-bar-linkages was robust to simulate walking, chair rise and stair ascent, but not stair descent (peak standard deviation of 2.66 BW), suggesting that joint model complexity should be set according to the imaging dataset available and the intended application, performing sensitivity analyses.

The development of a segment-based musculoskeletal model of the lower limb: introducing FreeBody

Traditional approaches to the biomechanical analysis of movement are joint-based; that is the mechanics of the body are described in terms of the forces and moments acting at the joints, and that muscular forces are considered to create moments about the joints. We have recently shown that segment-based approaches, where the mechanics of the body are described by considering the effect of the muscle, ligament and joint contact forces on the segments themselves, can also prove insightful. We have also previously described a simultaneous, optimization-based, musculoskeletal model of the lower limb. However, this prior model incorporates both joint- and segment-based assumptions. The purpose of this study was therefore to develop an entirely segment-based model of the lower limb and to compare its performance to our previous work. The segment-based model was used to estimate the muscle forces found during vertical jumping, which were in turn compared with the muscular activations that have been found in vertical jumping, by using a Geers' metric to quantify the magnitude and phase errors. The segment-based model was shown to have a similar ability to estimate muscle forces as a model based upon our previous work. In the future, we will evaluate the ability of the segment-based model to be used to provide results with clinical relevance, and compare its performance to joint-based approaches. The segment-based model described in this article is publicly available as a GUI-based Matlab® application and in the original source code (at www.msksoftware.org.uk).

Unilateral total hip replacement patients with symptomatic leg length inequality have abnormal hip biomechanics during walking

BACKGROUND

Symptomatic leg length inequality accounts for 8.7% of total hip replacement related claims made against the UK National Health Service Litigation authority. It has not been established whether symptomatic leg length inequality patients following total hip replacement have abnormal hip kinetics during gait.

METHODS

Hip kinetics in 15 unilateral total hip replacement patients with symptomatic leg length inequality during gait was determined through multibody dynamics and compared to 15 native hip healthy controls and 15 'successful' asymptomatic unilateral total hip replacement patients.

FINDING

More significant differences from normal were found in symptomatic leg length inequality patients than in asymptomatic total hip replacement patients. The leg length inequality patients had altered functions defined by lower gait velocity, reduced stride length, reduced ground reaction force, decreased hip range of motion, reduced hip moment and less dynamic hip force with a 24% lower heel-strike peak, 66% higher mid-stance trough and 37% lower toe-off peak. Greater asymmetry in hip contact force was also observed in leg length inequality patients.

INTERPRETATION

These gait adaptions may affect the function of the implant and other healthy joints in symptomatic leg length inequality patients. This study provides important information for the musculoskeletal function and rehabilitation of symptomatic leg length inequality patients.

Force-controlled dynamic wear testing of total ankle replacements

Currently, our knowledge of wear performance in total ankle replacements is limited. The aim of this study is to develop a scenario for force-controlled testing and wear testing of total ankle replacements. A force-controlled wear test was developed: based on cadaver measurements, the passive stabilization (ligaments and soft tissue) of the ankle joint was characterized and a restraint model for ankle stabilization was developed. Kinematics and kinetics acting at the replaced ankle joint were defined based on literature data and gait analysis. Afterwards, force-controlled wear testing was carried out on a mobile, three-component, total ankle replacement design. Wear was assessed gravimetrically and wear particles were analyzed. Wear testing resulted in a mean wear rate of 18.2±1.4mm(3)/10(6) cycles. Wear particles showed a mean size of 0.23μm with an aspect ratio of 1.61±0.96 and a roundness of 0.62±0.14. Wear testing of total ankle replacement shows that a relevant wear mass is generated with wear particles in a biologically relevant size range. The developed wear test provides a basis for future wear testing of total ankle replacements.

Muscle synergies may improve optimization prediction of knee contact forces during walking

The ability to predict patient-specific joint contact and muscle forces accurately could improve the treatment of walking-related disorders. Muscle synergy analysis, which decomposes a large number of muscle electromyographic (EMG) signals into a small number of synergy control signals, could reduce the dimensionality and thus redundancy of the muscle and contact force prediction process. This study investigated whether use of subject-specific synergy controls can improve optimization prediction of knee contact forces during walking. To generate the predictions, we performed mixed dynamic muscle force optimizations (i.e., inverse skeletal dynamics with forward muscle activation and contraction dynamics) using data collected from a subject implanted with a force-measuring knee replacement. Twelve optimization problems (three cases with four subcases each) that minimized the sum of squares of muscle excitations were formulated to investigate how synergy controls affect knee contact force predictions. The three cases were: (1) Calibrate+Match where muscle model parameter values were calibrated and experimental knee contact forces were simultaneously matched, (2) Precalibrate+Predict where experimental knee contact forces were predicted using precalibrated muscle model parameters values from the first case, and (3) Calibrate+Predict where muscle model parameter values were calibrated and experimental knee contact forces were simultaneously predicted, all while matching inverse dynamic loads at the hip, knee, and ankle. The four subcases used either 44 independent controls or five synergy controls with and without EMG shape tracking. For the Calibrate+Match case, all four subcases closely reproduced the measured medial and lateral knee contact forces (R2 ≥ 0.94, root-mean-square (RMS) error &lt; 66 N), indicating sufficient model fidelity for contact force prediction. For the Precalibrate+Predict and Calibrate+Predict cases, synergy controls yielded better contact force predictions (0.61 &lt; R2 &lt; 0.90, 83 N &lt; RMS error &lt; 161 N) than did independent controls (-0.15 &lt; R2 &lt; 0.79, 124 N &lt; RMS error &lt; 343 N) for corresponding subcases. For independent controls, contact force predictions improved when precalibrated model parameter values or EMG shape tracking was used. For synergy controls, contact force predictions were relatively insensitive to how model parameter values were calibrated, while EMG shape tracking made lateral (but not medial) contact force predictions worse. For the subject and optimization cost function analyzed in this study, use of subject-specific synergy controls improved the accuracy of knee contact force predictions, especially for lateral contact force when EMG shape tracking was omitted, and reduced prediction sensitivity to uncertainties in muscle model parameter values.

Why don't most runners get knee osteoarthritis? A case for per-unit-distance loads

Peak knee joint contact forces ("loads") in running are much higher than they are in walking, where the peak load has been associated with the initiation and progression of knee osteoarthritis. However, runners do not have an especially high risk of osteoarthritis compared with nonrunners. This paradox suggests that running somehow blunts the effect of very high peak joint contact forces, perhaps to provide a load per unit distance (PUD) traveled that is relatively low.

PURPOSE

This study aimed to compare peak and PUD knee joint loads between human walking and running.

METHODS

Fourteen healthy adults walked and ran at self-selected speeds. Ground reaction force and motion capture data were measured and combined with inverse dynamics and musculoskeletal modeling to estimate the peak knee joint loads, PUD knee joint loads, and the impulse of the knee joint contact force for each gait with a matched-pair (within-subject) design.

RESULTS

The peak load was three times higher in running (8.02 vs 2.72 body weight, P < 0.001), but the PUD load did not differ between running and walking (0.80 vs 0.75 body weight per meter, P = 0.098). The impulse of the joint contact force was greater for running than for walking (1.30 vs 1.04 body weight per second, P < 0.001). The peak load increased with increasing running speed, whereas the PUD load decreased with increasing speed.

CONCLUSIONS

Compared with walking, the relatively short duration of ground contact and relatively long length of strides in running seem to blunt the effect of high peak joint loads, such that the PUD loads are no higher than that in walking. Waveform features other than or in addition to the peak value should be considered when studying joint loading and injuries.

Femoral shaft strains during daily activities: Implications for atypical femoral fractures

BACKGROUND

Atypical femoral fractures are low-energy fractures initiating in the lateral femoral shaft. We hypothesized that atypical femoral fracture onset is associated with daily femoral strain patterns. We examined femoral shaft strains during daily activities.

METHODS

We analyzed earlier calculations of femoral strain during walking, sitting and rising from a chair, stair ascent, stair descent, stepping up, and squatting based on anatomically consistent musculoskeletal and finite-element models from a single donor and motion recordings from a body-matched volunteer. Femoral strains in the femoral shaft were extracted for the different activities and compared. The dependency between femoral strains in the lateral shaft and kinetic parameters was studied using multi-parametric linear regression analysis.

FINDINGS

Tensile strain in the lateral femoral shaft varied from 327 με (squatting) to 2004 με (walking). Walking and stair descent imposed tensile loading on the lateral shaft, whereas the other activities mainly imposed tensile loads on the anterior shaft. The multi-parametric linear regression showed a moderately strong correlation between tensile strains in the lateral shaft and the motion kinetic (joint moments and ground reaction force) in the proximal (R(2)=0.60) and the distal shaft (R(2)=0.46).

INTERPRETATION

Bone regions subjected to tensile strains are associated with atypical femoral fractures. Walking is the daily activity that induces the highest tensile strain in the lateral femoral shaft. The kinetics of motion explains 46%-50% of the tensile strain variation in the lateral shaft, whereas the unexplained part is likely to be attributed to the way joint moments are decomposed into muscle forces.

Hip contact forces in asymptomatic total hip replacement patients differ from normal healthy individuals: Implications for preclinical testing

BACKGROUND

Preclinical durability testing of hip replacement implants is standardised by ISO-14242-1 (2002) which is based on historical inverse dynamics analysis using data obtained from a small sample of normal healthy individuals. It has not been established whether loading cycles derived from normal healthy individuals are representative of loading cycles occurring in patients following total hip replacement.

METHODS

Hip joint kinematics and hip contact forces derived from multibody modelling of forces during normal walking were obtained for 15 asymptomatic total hip replacement patients and compared to 38 normal healthy individuals and to the ISO standard for pre-clinical testing.

FINDINGS

Hip kinematics in the total hip replacement patients were comparable to the ISO data and the hip contact force in the normal healthy group was also comparable to the ISO cycles. Hip contact forces derived from the asymptomatic total hip replacement patients were comparable for the first part of the stance period but exhibited 30% lower peak loads at toe-off.

INTERPRETATION

Although the ISO standard provides a representative kinematic cycle, the findings call into question whether the hip joint contact forces in the ISO standard are representative of those occurring in the joint following total hip replacement.

Changes in tibiofemoral forces due to variations in muscle activity during walking

Muscles induce large forces in the tibiofemoral joint during walking and thereby influence the health of tissues like articular cartilage and menisci. It is possible to walk with a wide variety of muscle coordination patterns, but the effect of varied muscle coordination on tibiofemoral contact forces remains unclear. The goal of this study was to determine the effect of varied muscle coordination on tibiofemoral contact forces. We developed a musculoskeletal model of a subject walking with an instrumented knee implant. Using an optimization framework, we calculated the tibiofemoral forces resulting from muscle coordination that reproduced the subject's walking dynamics. We performed a large set of optimizations in which we systematically varied the coordination of muscles to determine the influence on tibiofemoral force. Model-predicted tibiofemoral forces arising with minimum muscle activation matched in vivo forces measured during early stance, but were greater than in vivo forces during late stance. Peak tibiofemoral forces during late stance could be reduced by increasing the activation of the gluteus medius, uniarticular hip flexors, and soleus, and by decreasing the activation of the gastrocnemius and rectus femoris. These results suggest that retraining of muscle coordination could substantially reduce tibiofemoral forces during late stance.

Assessment of external and internal loads in the triple jump via inverse dynamics simulation

The triple jump is a demanding athletics event that, after an approach run, consists of three consecutive phases: the hop, the bound, and the jump. During the involved three take-off actions a jumper is exposed to increased risk of injury due to the high impact forces from the ground and powerful muscle/tendon efforts, which are further reflected in the internal loads of the lower limb joints. While external ground reactions can possibly be measured using force platforms, in vivo measurements of the internal loads are practically not feasible. The purpose of the paper is to present the development of an effective formulation for the inverse dynamics simulation of the triple jump, based on the jumper dynamical model and non-invasive kinematic recordings of the movement. The developed simulation model serves for the analysis of all the triple jump phases, irrespective of whether the jumper is in flight or in contact with the ground with one of his feet, and is focused on effective assessment of the external reactions on the supporting leg as well as the muscle forces and joint reaction forces in the leg. Some numerical results of inverse dynamics simulation of the triple jump are reported.