Subject-specific geometrical detail rather than cost function formulation affects hip loading calculation

A three-dimensional musculoskeletal model of the pelvis and lower limb of Australopithecus afarensis

OBJECTIVES

Musculoskeletal modeling is a powerful approach for studying the biomechanics and energetics of locomotion. Australopithecus (A.) afarensis is among the best represented fossil hominins and provides critical information about the evolution of musculoskeletal design and locomotion in the hominin lineage. Here, we develop and evaluate a three-dimensional (3-D) musculoskeletal model of the pelvis and lower limb of A. afarensis for predicting muscle-tendon moment arms and moment-generating capacities across lower limb joint positions encompassing a range of locomotor behaviors.

MATERIALS AND METHODS

A 3-D musculoskeletal model of an adult A. afarensis pelvis and lower limb was developed based primarily on the A.L. 288-1 partial skeleton. The model includes geometric representations of bones, joints and 35 muscle-tendon units represented using 43 Hill-type muscle models. Two muscle parameter datasets were created from human and chimpanzee sources. 3-D muscle-tendon moment arms and isometric joint moments were predicted over a wide range of joint positions.

RESULTS

Predicted muscle-tendon moment arms generally agreed with skeletal metrics, and corresponded with human and chimpanzee models. Human and chimpanzee-based muscle parameterizations were similar, with some differences in maximum isometric force-producing capabilities. The model is amenable to size scaling from A.L. 288-1 to the larger KSD-VP-1/1, which subsumes a wide range of size variation in A. afarensis.

DISCUSSION

This model represents an important tool for studying the integrated function of the neuromusculoskeletal systems in A. afarensis. It is similar to current human and chimpanzee models in musculoskeletal detail, and will permit direct, comparative 3-D simulation studies.

EMG-Driven Musculoskeletal Model Calibration With Wrapping Surface Personalization

Muscle forces and joint moments estimated by electromyography (EMG)-driven musculoskeletal models are sensitive to the wrapping surface geometry defining muscle-tendon lengths and moment arms. Despite this sensitivity, wrapping surface properties are typically not personalized to subject movement data. This study developed a novel method for personalizing OpenSim cylindrical wrapping surfaces during EMG-driven model calibration. To avoid the high computational cost of repeated OpenSim muscle analyses, the method uses two-level polynomial surrogate models. Outer-level models specify time-varying muscle-tendon lengths and moment arms as functions of joint angles, while inner-level models specify time-invariant outer-level polynomial coefficients as functions of wrapping surface parameters. To evaluate the method, we used walking data collected from two individuals post-stroke and performed four variations of EMG-driven lower extremity model calibration: 1) no calibration of scaled generic wrapping surfaces (NGA), 2) calibration of outer-level polynomial coefficients for all muscles (SGA), 3) calibration of outer-level polynomial coefficients only for muscles with wrapping surfaces (LSGA), and 4) calibration of cylindrical wrapping surface parameters for muscles with wrapping surfaces (PGA). On average compared to NGA, SGA reduced lower extremity joint moment matching errors by 31%, LSGA by 24%, and PGA by 12%, with the largest reductions occurring at the hip. Furthermore, PGA reduced peak hip joint contact force by 47% bodyweight, which was the most consistent with published in vivo measurements. The proposed method for EMG-driven model calibration with wrapping surface personalization produces physically realistic OpenSim models that reduce joint moment matching errors while improving prediction of hip joint contact force.

Comparative validation of two patient-specific modelling pipelines for predicting knee joint forces during level walking

Over the past few years, the use of computer models and simulations tailored to the patient's physiology to assist clinical decision-making has increased enormously.While several pipelines to develop personalized models exist, their adoption on a large scale is still limited due to the required niche computational skillset and the lengthy operations required. Novel toolboxes, such as STAPLE, promise to streamline and expedite the development of image-based skeletal lower limb models. STAPLE-generated models can be rapidly generated, with minimal user input, and present similar joint kinematics and kinetics compared to models developed employing the established INSIGNEO pipeline. Yet, it is unclear how much the observed discrepancies scale up and affect joint contact force predictions. In this study, we compared image-based musculoskeletal models developed (i) with the INSIGNEO pipeline and (ii) with a semi-automated pipeline that combines STAPLE and nmsBuilder, and assessed their accuracy against experimental implant data.Our results showed that both pipelines predicted similar total knee joint contact forces between one another in terms of profiles and average values, characterized by a moderately high level of agreement with the experimental data. Nonetheless, the Student t-test revealed statistically significant differences between both pipelines. Of note, the STAPLE-based pipeline required considerably less time than the INSIGNEO pipeline to generate a musculoskeletal model (i.e., 60 vs 160 min). This is likely to open up opportunities for the use of personalized musculoskeletal models in clinical practice, where time is of the essence.

On the estimation of hip joint loads through musculoskeletal modeling

Noninvasive estimation of joint loads is still an open challenge in biomechanics. Although musculoskeletal modeling represents a solid resource, multiple improvements are still necessary to obtain accurate predictions of joint loads and to translate such potential into practical utility. The present study, focused on the hip joint, is aimed at reviewing the state-of-the-art literature on the estimation of hip joint reaction forces through musculoskeletal modeling. Our literature inspection, based on well-defined selection criteria, returned seventeen works, which were compared in terms of methods and results. Deviations between predicted and in vivo measured hip joint loads, taken from the OrthoLoad database, were assessed through quantitative deviation indices. Despite the numerous modeling and computational improvements made over the last two decades, predicted hip joint loads still deviate from their experimental counterparts and typically overestimate them. Several critical aspects have emerged that affect muscle force estimation, hence joint loads. Among them, the physical fidelity of the musculoskeletal model, with its parameters and geometry, plays a crucial role. Also, predicted joint loads are markedly affected by the selected muscle recruitment strategy, which reflects the underlying motor control policy. Practical guidelines for researchers interested in noninvasive estimation of hip joint loads are also provided.

Biomechanical compensations during a stand-to-sit maneuver using transfemoral osseointegrated prostheses: A case series

BACKGROUND

Patients with transfemoral amputation and socket prostheses are at a heightened risk of developing musculoskeletal overuse injuries, commonly due to altered joint biomechanics. Osseointegrated prostheses, which involve direct anchorage of the prosthesis to the residual limb through a bone anchored prosthesis, are a novel alternative to sockets yet their biomechanical effect is largely unknown.

METHODS

Four patients scheduled to undergo unilateral transfemoral prosthesis osseointegration completed two data collections (baseline with socket prosthesis and 12-months after prosthesis osseointegration) in which whole-body kinematics and ground reaction forces were collected during stand-to-sit tasks. Trunk, pelvis, and hip kinematics, and the surrounding muscle forces, were calculated using subject-specific musculoskeletal models developed in OpenSim. Peak joint angles and muscle forces were compared between timepoints using Cohen's d effect sizes.

FINDINGS

Compared to baseline with socket prostheses, patients with osseointegrated prostheses demonstrated reduced lateral trunk bending (d = 1.46), pelvic obliquity (d = 1.09), and rotation (d = 1.77) toward the amputated limb during the stand to sit task. This was accompanied by increased amputated limb hip flexor, abductor, and rotator muscle forces (d> > 0.8).

INTERPRETATION

Improved lumbopelvic movement patterns and stabilizing muscle forces when using an osseointegrated prosthesis indicate that this novel prosthesis type likely reduces the risk of the development and/or progression of overuse injuries, such as low back pain and osteoarthritis. We attribute the increased muscle hip muscle forces to the increased load transmission between the osseointegrated prosthesis and residual limb, which allows a greater eccentric ability of the amputated limb to control lowering during the stand-to-sit task.

Thin patient-specific clavicle fracture fixation plates can mechanically outperform commercial plates: An in silico approach

Current fixation plates used to operatively stabilize clavicular fractures are suboptimal, leading to reoperation rates up to 53%. Plate irritation, which can be caused by a poor geometric fit and plate thickness, has been found to be an important factor for reoperation. Moreover, muscle attachment sites (MAS) have to be detached occasionally. To improve current surgical clavicle fracture treatment with plate osteosynthesis, a change in plate design is required. The goal of this study was to design a patient-specific clavicle fracture fixation plate that includes geometrical optimization and stiffness determination. Its biomechanical performance has been compared with a commercial plate by examining the geometric fit, anatomical outline, stresses and interfragmentary motion using a finite element analysis with physiological loading and boundary conditions. Evaluation showed a better geometrical fit of the patient-specific plate as well as an improved fracture reduction. Displacements between fracture fragments were lower in case of the patient-specific plate, both when a fracture gap and no fracture gap were present. Results indicate a superior mechanical performance in terms of all investigated outcomes of the patient-specific plate compared to the commercial plate, while better aligning with the patient-specific geometry and without the need for MAS release. Due to the patient-specific geometry and reduced thickness, these fixation plates are expected to decrease the operation time, as intraoperative contouring will become irrelevant, and to decrease reoperation rates as implant irritation will be minimized.

Electromyography-Assisted Neuromusculoskeletal Models Can Estimate Physiological Muscle Activations and Joint Moments Across the Neck Before Impacts

Knowledge of neck muscle activation strategies before sporting impacts is crucial for investigating mechanisms of severe spinal injuries. However, measurement of muscle activations during impacts is experimentally challenging and computational estimations are not often guided by experimental measurements. We investigated neck muscle activations before impacts with the use of electromyography (EMG)-assisted neuromusculoskeletal models. Kinematics and EMG recordings from four major neck muscles of a rugby player were experimentally measured during rugby activities. A subject-specific musculoskeletal model was created with muscle parameters informed from MRI measurements. The model was used in the calibrated EMG-informed neuromusculoskeletal modeling toolbox and three neural solutions were compared: (i) static optimization (SO), (ii) EMG-assisted (EMGa), and (iii) MRI-informed EMG-assisted (EMGaMRI). EMGaMRI and EMGa significantly (p < 0.01) outperformed SO when tracking cervical spine net joint moments from inverse dynamics in flexion/extension (RMSE = 0.95, 1.14, and 2.32 N·m) but not in lateral bending (RMSE = 1.07, 2.07, and 0.84 N·m). EMG-assisted solutions generated physiological muscle activation patterns and maintained experimental cocontractions significantly (p < 0.01) outperforming SO, which was characterized by saturation and nonphysiological "on-off" patterns. This study showed for the first time that physiological neck muscle activations and cervical spine net joint moments can be estimated without assumed a priori objective criteria before impacts. Future studies could use this technique to provide detailed initial loading conditions for theoretical simulations of neck injury during impacts.

Femoral version deformities alter joint reaction forces in dysplastic hips during gait

Developmental dysplasia of the hip (DDH) causes hip instability and early-onset osteoarthritis. The focus on pathomechanics in DDH has centered on the shallow acetabulum, however there is growing awareness of the role of femoral deformities in joint damage. The objective of this study was to determine the influence of femoral version (FV) on the muscle and joint reaction forces (JRFs) of dysplastic hips during gait. Magnetic resonance images, in-vivo gait data, and musculoskeletal models were used to calculate JRFs and simulate changes due to varying FV deformities. Rotation about the long axis of the femur was added in the musculoskeletal models to simulate FV values from -5° (relative retroversion) to + 35° (increased anteversion). In our simulations, FV deformities caused the largest changes to the anteroposterior and resultant JRFs. From a normal FV of 15°, a 15° increase in femoral anteversion caused JRFs to be less posterior in early stance (Δ = 0.43 ± 0.22 xbodyweight) and more anterior in late stance (Δ = 0.60 ± 14 xbodyweight). Relative retroversion caused anteroposterior changes that were similar to anteversion in early stance but opposite in late stance. Resultant JRFs experienced the largest changes during late stance where anteversion raised the peak by 0.48 ± 0.15 xbodyweight and relative retroversion lowered the peak by 0.32 ± 0.30 xbodyweight. Increasing anteversion increased hip flexor and abductor muscle forces, which caused the changes in JRFs. Identifying how FV deformities influence hip joint loading can elucidate their role in the mechanisms of hip degeneration in patients with DDH.

A Musculoskeletal Model Customized for Sagittal and Frontal Knee Kinematics with Improved Knee Joint Stability

Current lower limb musculoskeletal (MSK) models focus on sagittal plane kinematics. However, abnormal gait is typically associated with sagittal plane motions crossing into other planes, limiting the use of current MSK models. The purpose of this study was twofold, first, to extend the capability of a full-body MSK model from the literature to include frontal knee plane kinematics during healthy gait, and second, to propose and implement a realistic muscle discretization technique. Two MSK model constructs were derived - the first construct (Knee2_SM) allowed two degrees of freedom (sagittal and coronal) at the knee and the second construct (Knee2_MM) implemented multi-line elements for all the lower limb muscles in conjunction with two knee degrees of freedom. Motion analysis data of normal gait cycle from 10 healthy adults were used to compare joint kinematics, muscle moment arms, muscle forces, and muscle activations, between new constructs and the original model. Knee varus-valgus trajectories were estimated with the mean peak values ranging from 9.49° valgus to 1.57° varus. Knee2_MM predicted significant difference (p < 0.05) in moment arms and forces in those muscles responsible for medial-lateral stability of the knee. The simulated muscle activations generated by Knee2_MM model matched more closely to the experimental EMG when qualitatively compared. This study enhances the capability of sagittal plane full-body MSK model to incorporate knee varus-valgus motion while keeping the joint stability intact and improving muscle prediction.

A Systematic Review of the Associations Between Inverse Dynamics and Musculoskeletal Modeling to Investigate Joint Loading in a Clinical Environment

The assessment of knee or hip joint loading by external joint moments is mainly used to draw conclusions on clinical decision making. However, the correlation between internal and external loads has not been systematically analyzed. This systematic review aims, therefore, to clarify the relationship between external and internal joint loading measures during gait. A systematic database search was performed to identify appropriate studies for inclusion. In total, 4,554 articles were identified, while 17 articles were finally included in data extraction. External joint loading parameters were calculated using the inverse dynamics approach and internal joint loading parameters by musculoskeletal modeling or instrumented prosthesis. It was found that the medial and total knee joint contact forces as well as hip joint contact forces in the first half of stance can be well predicted using external joint moments in the frontal plane, which is further improved by including the sagittal joint moment. Worse correlations were found for the peak in the second half of stance as well as for internal lateral knee joint contact forces. The estimation of external joint moments is useful for a general statement about the peak in the first half of stance or for the maximal loading. Nevertheless, when investigating diseases as valgus malalignment, the estimation of lateral knee joint contact forces is necessary for clinical decision making because external joint moments could not predict the lateral knee joint loading sufficient enough. Dependent on the clinical question, either estimating the external joint moments by inverse dynamics or internal joint contact forces by musculoskeletal modeling should be used.

MRI-based anatomical characterisation of lower-limb muscles in older women

The ability of muscles to produce force depends, among others, on their anatomical features and it is altered by ageing-associated weakening. However, a clear characterisation of these features, highly relevant for older individuals, is still lacking. This study hence aimed at characterising muscle volume, length, and physiological cross-sectional area (PCSA) and their variability, between body sides and between individuals, in a group of post-menopausal women. Lower-limb magnetic resonance images were acquired from eleven participants (69 (7) y. o., 66.9 (7.7) kg, 159 (3) cm). Twenty-three muscles were manually segmented from the images and muscle volume, length and PCSA were calculated from this dataset. Personalised maximal isometric force was then calculated using the latter information. The percentage difference between the muscles of the two lower limbs was up to 89% and 22% for volume and length, respectively, and up to 84% for PCSA, with no recognisable pattern associated with limb dominance. Between-subject coefficients of variation reached 36% and 13% for muscle volume and length, respectively. Generally, muscle parameters were similar to previous literature, but volumes were smaller than those from in-vivo young adults and slightly higher than ex-vivo ones. Maximal isometric force was found to be on average smaller than those obtained from estimates based on linear scaling of ex-vivo-based literature values. In conclusion, this study quantified for the first time anatomical asymmetry of lower-limb muscles in older women, suggesting that symmetry should not be assumed in this population. Furthermore, we showed that a scaling approach, widely used in musculoskeletal modelling, leads to an overestimation of the maximal isometric force for most muscles. This heavily questions the validity of this approach for older populations. As a solution, the unique dataset of muscle segmentation made available with this paper could support the development of alternative population-based scaling approaches, together with that of automatic tools for muscle segmentation.

Botulinum toxin injections minimally affect modelled muscle forces during gait in children with cerebral palsy

BACKGROUND

Children with cerebral palsy (CP) present altered gait patterns and electromyography (EMG) activity compared to typically developing children. To temporarily reduce muscular activity and to correct the abnormal muscle force balance, Botulinum Toxin type A (BTX-A) injections are used.

RESEARCH QUESTION

What is the effect of BTX-A injections on dynamic muscle forces during gait, when calculated using an EMG-constrained approach?.

METHODS

Retrospective data of ten typically developing (TD) and fourteen children with spastic diplegic CP were used for musculoskeletal modeling and dynamic simulations of gait, before and after BTX-A treatment. Individual muscle forces were calculated using an EMG-constrained optimization, in which EMG of eight muscles was used as muscle excitation signal to constrain the muscle activation patterns. Paired t-tests were used to compare average modelled muscle forces in different phases of the gait cycle pre- and post-BTX-A, summarized in the muscle profile score. Two-sample t-tests were used to determine significant differences between TD and pre- and post-BTX-A modelled muscle forces.

RESULTS

For most muscles, the force was decreased in CP compared to TD children in all phases of the gait cycle, both before and after BTX-A treatment. Differences in muscle forces before and after BTX-A treatment were limited, with only few significant differences between pre- and post-BTX-A. Compared to a standard static optimization approach, imposing the EMG activity increased modelled muscle forces for most muscles.

SIGNIFICANCE

Our findings indicate that BTX-A treatment has a limited effect on the muscle balance in CP children. Besides that, the use of EMG-constrained optimization is recommended when studying muscle balance in children with CP.

Dysplastic hip anatomy alters muscle moment arm lengths, lines of action, and contributions to joint reaction forces during gait

Developmental dysplasia of the hip (DDH) is characterized by abnormal bony anatomy, which causes detrimental hip joint loading and leads to secondary osteoarthritis. Hip joint loading depends, in part, on muscle-induced joint reaction forces (JRFs), and therefore, is influenced by hip muscle moment arm lengths (MALs) and lines of action (LoAs). The current study used subject-specific musculoskeletal models and in-vivo motion analysis to quantify the effects of DDH bony anatomy on dynamic muscle MALs, LoAs, and their contributions to JRF peaks during early (~17%) and late-stance (~52%) of gait. Compared to healthy hips (N = 15, 16-39 y/o), the abductor muscles in patients with untreated DDH (N = 15, 16-39 y/o) had smaller abduction MALs (e.g. anterior gluteus medius, 35.3 vs. 41.6 mm in early stance, 45.4 vs. 52.6 mm late stance, p ≤ 0.01) and more medially-directed LoAs. Abduction-adduction and rotation MALs also differed for major hip flexors such as rectus femoris and iliacus. The altered MALs in DDH corresponded to higher hip abductor forces, medial JRFs (1.26 vs. 0.87 × BW early stance, p = 0.03), and resultant JRFs (5.71 vs. 4.97 × BW late stance, p = 0.05). DDH anatomy not only affected hip muscle force generation in the primary plane of function, but also their out-of-plane mechanics, which collectively elevated JRFs. Overall, hip muscle MALs and their contributions to JRFs were significantly altered by DDH bony anatomy. Therefore, to better understand the mechanisms of joint degeneration and improve the efficacy of treatments for DDH, the dynamic anatomy-force relationships and multi-planar functions of the whole hip musculature must be collectively considered.

Predicting longitudinal changes in joint contact forces in a juvenile population: scaled generic versus subject-specific musculoskeletal models

Subject-specific musculoskeletal model use in clinical settings is limited due to development-associated time and effort burdens together with potential medical imaging unavailability. As an alternative, this study investigated consistency in estimating longitudinal changes in joint contact forces (JCF) between scaled generic and subject-specific models. For 11 children, joint kinematics and JCF were calculated using subject-specific and scaled generic models. JCF changes estimated by both models were strongly correlated for the hip and knee although JCF estimates varied between models. Findings suggest that within specified limits of accuracy, scaled generic models are sensitive enough to detect JCF changes consistent with subject-specific models.

Influence of musculotendon geometry variability in muscle forces and hip bone-on-bone forces during walking

Inverse dynamics problems are usually solved in the analysis of human gait to obtain reaction forces and moments at the joints. However, these actions are not the actual forces and moments supported by the joint structure, because they do not consider the forces of the muscles acting across the joint. Therefore, to analyse bone-on bone forces it is necessary to estimate those muscle forces. Usually, this problem is addressed by means of optimization algorithms. One of the parameters required to solve this problem is the musculotendon geometry. These data are usually taken from cadavers or MRI data from several subjects, different from the analysed subject. Then, the model is scaled to the subject morphology. This procedure constitutes a source of error. The goals of this work were two. First, to perform a sensitivity analysis of the influence of muscle insertion locations on the muscle forces acting on the hip joint and on the hip joint bone-on-bone forces. Second, to compare the hip joint bone-on-bone forces during gait cycle obtained through muscle insertion locations taken from a musculoskeletal model template and a scaling procedure to those obtained from a subject-specific model using an MRI of the subject. The problem was solved using OpenSim. Results showed that anatomical variability should be analysed from two perspectives. One the one hand, throughout the gait cycle, in a global way. On the other hand, at a characteristic instant of the gait cycle. Variations of ±1 cm in the position of the attachment points of certain muscles caused variations of up to 14.21% in averaged deviation of the muscle forces and 58.96% in the peak force in the modified muscle and variations up to 57.23% in the averaged deviation of the muscle force and up to 117.23% in the peak force in the rest of muscles. Then, the influence of that variability on muscle activity patterns and hip bone-on-bone forces could be described more precisely. A biomechanical analysis of a subject-specific musculoskeletal model was carried out. Using MRI data, variations up to 5 cm in the location of the insertion points were introduced. These modifications showed significant differences between the baseline model and the customized model: within the range [-12%, 10%] for muscle forces and around 35% of body weight for hip bone-on-bone forces.

Subject-specific geometry affects acetabular contact pressure during gait more than subject-specific loading patterns

Finite element modeling (FEM) can predict hip cartilage contact mechanics. This study investigated how subject-specific boundary conditions and joint geometry affect acetabular cartilage contact mechanics using a multi-scale workflow. For two healthy subjects, musculoskeletal models calculated subject-specific hip kinematics and loading, which were used as boundary conditions for FEM. Cartilage contact mechanics were predicted using either generic or subject-specific FEM and boundary conditions. A subject-specific mesh resulted in a more lateral contact. Effects of subject-specific boundary conditions varied between both subjects. Results highlight the complex interplay between loading and kinematics and their effect on cartilage contact mechanics.

Musculoskeletal models with generic and subject-specific geometry estimate different joint biomechanics in dysplastic hips

Optimizing the geometric complexity of musculoskeletal models is important for reliable yet feasible estimation of joint biomechanics. This study investigated the effects of subject-specific model geometry on hip joint reaction forces (JRFs) and muscle forces in patients with developmental dysplasia of the hip (DDH) and healthy controls. For nine DDH and nine control subjects, three models were created with increasingly subject-specific pelvis geometry, hip joint center locations and muscle attachments. Hip JRFs and muscle forces during a gait cycle were compared among the models. For DDH subjects, resultant JRFs from highly specific models including subject-specific pelvis geometry, joint locations and muscle attachments were not significantly different compared to models using generic geometry in early stance, but were significantly higher in late stance (p = 0.03). Estimates from moderately specific models using CT-informed scaling of generic pelvis geometry were not significantly different from low specificity models using generic geometry scaled with skin markers. For controls, resultant JRFs in early stance from highly specific models were significantly lower than moderate and low specificity models (p ≤ 0.02) with no significant differences in late stance. Inter-model JRF differences were larger for DDH subjects than controls. Inter-model differences for JRF components and muscle forces were similar to resultant JRFs. Incorporating subject-specific pelvis geometry significantly affects JRF and muscle force estimates in both DDH and control groups, which may be especially important for reliable estimation of pathomechanics in dysplastic hips.

Refining muscle geometry and wrapping in the TLEM 2 model for improved hip contact force prediction

Musculoskeletal models represent a powerful tool to gain knowledge on the internal forces acting at the joint level in a non-invasive way. However, these models can present some errors associated with the level of detail in their geometrical representation. For this reason, a thorough validation is necessary to prove the reliability of their predictions. This study documents the development of a generic musculoskeletal model and proposes a working logic and simulation techniques for identifying specific model features in need of refinement; as well as providing a quantitative validation for the prediction of hip contact forces (HCF). The model, implemented in the AnyBody Modeling System and based on the cadaveric dataset TLEM 2.0, was scaled to match the anthropometry of a patient fitted with an instrumented hip implant and to reproduce gait kinematics based on motion capture data. The relative contribution of individual muscle elements to the HCF and joint moments was analyzed to identify critical geometries, which were then compared to muscle magnetic resonance imaging (MRI) scans and, in case of inconsistencies, were modified to better match the volumetric scans. The predicted HCF showed good agreement with the overall trend and timing of the measured HCF from the instrumented prosthesis. The average root mean square error (RMSE), calculated for the total HCF was found to be 0.298*BW. Refining the geometries of the muscles thus identified reduced RMSE on HCF magnitudes by 17% (from 0.359*BW to 0.298*BW) over the whole gait cycle. The detailed study of individual muscle contributions to the HCF succeeded in identifying muscles with incorrect anatomy, which would have been difficult to intuitively identify otherwise. Despite a certain residual over-prediction of the final hip contact forces in the stance phase, a satisfactory level of geometrical accuracy of muscle paths has been achieved with the refinement of this model.

Advancing quantitative techniques to improve understanding of the skeletal structure-function relationship

Although all functional movement arises from the interplay between the neurological, skeletal, and muscular systems, it is the skeletal system that forms the basic framework for functional movement. Central to understanding human neuromuscular development, along with the genesis of musculoskeletal pathologies, is quantifying how the human skeletal system adapts and mal-adapts to its mechanical environment. Advancing this understanding is hampered by an inability to directly and non-invasively measure in vivo strains, stresses, and forces on bone. Thus, we traditionally have turned to animal models to garner such information. These models enable direct in vivo measures that are not available for human subjects, providing information in regards to both skeletal adaptation and the interplay between the skeletal and muscular systems. Recently, there has been an explosion of new imaging and modeling techniques providing non-invasive, in vivo measures and estimates of skeletal form and function that have long been missing. Combining multiple modalities and techniques has proven to be one of our most valuable resources in enhancing our understanding of the form-function relationship of the human skeletal, muscular, and neurological systems. Thus, to continue advancing our knowledge of the structural-functional relationship, validation of current tools is needed, while development is required to limit the deficiencies in these tools and develop new ones.

Longitudinal joint loading in patients before and up to one year after unilateral total hip arthroplasty

Abnormal kinematics and kinetics have been reported in hip osteoarthritis (OA) patients before and after total hip arthroplasty (THA). These changes can affect the loading of the ipsilateral hip, as well as the contralateral hip and knee joint. As it is not clear how hip and knee loading evolves in THA patients during the first year after surgery, the goal of this study is to define how joint loading changes in patients before and at three evaluation times after THA surgery. Musculoskeletal modelling in combination with gait analysis data was used to calculate hip and knee contact forces in 14 patients before and 3-, 6- and 12-months after unilateral THA, as well as in 18 healthy controls. Results showed that bilateral hip and knee loading were decreased compared to controls, both before and after THA surgery. Loading symmetry was altered compared to controls at 3-months post-surgery for the hip and at all evaluation times, except for 6-months post-surgery, for the knee, with ipsilateral joint loading decreased compared to the contralateral side. To conclude, 12-months after THA joint loading was not normalized, with both hip and knee loading in patients decreased compared to controls. Therefore, no overloading of the ipsi- or contralateral hip and knee joint was found before and up to one year after unilateral THA.

Bioinspired Technologies to Connect Musculoskeletal Mechanobiology to the Person for Training and Rehabilitation

Musculoskeletal tissues respond to optimal mechanical signals (e.g., strains) through anabolic adaptations, while mechanical signals above and below optimal levels cause tissue catabolism. If an individual's physical behavior could be altered to generate optimal mechanical signaling to musculoskeletal tissues, then targeted strengthening and/or repair would be possible. We propose new bioinspired technologies to provide real-time biofeedback of relevant mechanical signals to guide training and rehabilitation. In this review we provide a description of how wearable devices may be used in conjunction with computational rigid-body and continuum models of musculoskeletal tissues to produce real-time estimates of localized tissue stresses and strains. It is proposed that these bioinspired technologies will facilitate a new approach to physical training that promotes tissue strengthening and/or repair through optimal tissue loading.