Automated 2D and 3D finite element overclosure adjustment and mesh morphing using generalized regression neural networks.

Computer representations of three-dimensional (3D) geometries are crucial for simulating systems and processes in engineering and science. In medicine, and more specifically, biomechanics and orthopaedics, obtaining and using 3D geometries is critical to many workflows. However, while many tools exist to obtain 3D geometries of organic structures, little has been done to make them usable for their intended medical purposes. Furthermore, many of the proposed tools are proprietary, limiting their use. This work introduces two novel algorithms based on Generalized Regression Neural Networks (GRNN) and 4 processes to perform mesh morphing and overclosure adjustment. These algorithms were implemented, and test cases were used to validate them against existing algorithms to demonstrate improved performance. The resulting algorithms demonstrate improvements to existing techniques based on Radial Basis Function (RBF) networks by converting to GRNN-based implementations. Implementations in MATLAB of these algorithms and the source code are publicly available at the following locations: https://github.com/thor-andreassen/femors; https://simtk.org/projects/femors-rbf; https://www.mathworks.com/matlabcentral/fileexchange/120353-finite-element-morphing-overclosure-reduction-and-slicing.

Kinematic Performance of Medial Pivot Total Knee Arthroplasty.

BACKGROUND: Total knee arthroplasty (TKA) implants have continued to evolve to accommodate new understandings of knee mechanics. The medial-pivot implant is a newer design, which is intended to limit anterior-posterior translation in the medial compartment while allowing lateral compartment translation. However, evidence for a generalized medial-pivot characteristic across all activities is limited. The purpose of the study was to quantify and compare in vivo knee joint kinematics using high-speed stereo radiography during activities of daily living in patients who have undergone a TKA with a cruciate sacrificing medial-pivot implant to age-matched and sex-matched native controls. METHODS: Fifteen participants (7 patients, 4 women, mean age 70 years and 8 nonsymptomatic controls, 4 women, mean age 64 years) performed 6 functional tasks in high-speed stereo radiography: deep-knee lunge, chair rise, step down, gait, gait with 90° turn, and seated knee extension. Translational differences between groups (surgical versus control) were assessed for the medial and lateral condyle, while pivot location was normalized to subject-specific tibial plateau geometry. RESULTS: The surgical cohort displayed a more constrained medial condyle that provided greater stability of the medial compartment and did not result in the paradoxical anterior translation at mid-flexion angles during weight-bearing activities, but was associated with less condylar translation than native knees. Additionally, the transverse tibial pivot location occurs most commonly in the middle third of the tibial plateau and secondarily on the medial third. CONCLUSIONS: Some variability in pivot location occurs between activities and is more in nonsymptomatic, native knee controls.

Neural-driven activation of 3D muscle within a finite element framework: exploring applications in healthy and neurodegenerative simulations.

This paper presents a novel computational framework for neural-driven finite element muscle models, with an application to amyotrophic lateral sclerosis (ALS). The multiscale neuromusculoskeletal (NMS) model incorporates physiologically accurate motor neurons, 3D muscle geometry, and muscle fiber recruitment. It successfully predicts healthy muscle force and tendon elongation and demonstrates a progressive decline in muscle force due to ALS, dropping from 203 N (healthy) to 155 N (120 days after ALS onset). This approach represents a preliminary step towards developing integrated neural and musculoskeletal simulations to enhance our understanding of neurodegenerative and neurodevelopmental conditions through predictive NMS models.

Deciphering the "Art" in Modeling and Simulation of the Knee Joint: Assessing Model Calibration Workflows and Outcomes.

Model reproducibility is a point of emphasis for the National Institutes of Health (NIH) and in science, broadly. As the use of computational modeling in biomechanics and orthopedics grows, so does the need to assess the reproducibility of modeling workflows and simulation predictions. The long-term goal of the KneeHub project is to understand the influence of potentially subjective decisions, thus the modeler's "art", on the reproducibility and predictive uncertainty of computational knee joint models. In this paper, we report on the model calibration phase of this project, during which five teams calibrated computational knee joint models of the same specimens from the same specimen-specific joint mechanics dataset. We investigated model calibration approaches and decisions, and compared calibration workflows and model outcomes among the teams. The selection of the calibration targets used in the calibration workflow differed greatly between the teams and was influenced by modeling decisions related to the representation of structures, and considerations for computational cost and implementation of optimization. While calibration improved model performance, differences in the postcalibration ligament properties and predicted kinematics were quantified and discussed in the context of modeling decisions. Even for teams with demonstrated expertise, model calibration is difficult to foresee and plan in detail, and the results of this study underscore the importance of identification and standardization of best practices for data sharing and calibration.

Modeling of active skeletal muscles: a 3D continuum approach incorporating multiple muscle interactions.

Skeletal muscles have a highly organized hierarchical structure, whose main function is to generate forces for movement and stability. To understand the complex heterogeneous behaviors of muscles, computational modeling has advanced as a non-invasive approach to evaluate relevant mechanical quantities. Aiming to improve musculoskeletal predictions, this paper presents a framework for modeling 3D deformable muscles that includes continuum constitutive representation, parametric determination, model validation, fiber distribution estimation, and integration of multiple muscles into a system level for joint motion simulation. The passive and active muscle properties were modeled based on the strain energy approach with Hill-type hyperelastic constitutive laws. A parametric study was conducted to validate the model using experimental datasets of passive and active rabbit leg muscles. The active muscle model with calibrated material parameters was then implemented to simulate knee bending during a squat with multiple quadriceps muscles. A computational fluid dynamics (CFD) fiber simulation approach was utilized to estimate the fiber arrangements for each muscle, and a cohesive contact approach was applied to simulate the interactions among muscles. The single muscle simulation results showed that both passive and active muscle elongation responses matched the range of the testing data. The dynamic simulation of knee flexion and extension showed the predictive capability of the model for estimating the active quadriceps responses, which indicates that the presented modeling pipeline is effective and stable for simulating multiple muscle configurations. This work provided an effective framework of a 3D continuum muscle model for complex muscle behavior simulation, which will facilitate additional computational and experimental studies of skeletal muscle mechanics. This study will offer valuable insight into the future development of multiscale neuromuscular models and applications of these models to a wide variety of relevant areas such as biomechanics and clinical research.

Reproducibility in modeling and simulation of the knee: Academic, industry, and regulatory perspectives.

Stakeholders in the modeling and simulation (M&S) community organized a workshop at the 2019 Annual Meeting of the Orthopaedic Research Society (ORS) entitled "Reproducibility in Modeling and Simulation of the Knee: Academic, Industry, and Regulatory Perspectives." The goal was to discuss efforts among these stakeholders to address irreproducibility in M&S focusing on the knee joint. An academic representative from a leading orthopedic hospital in the United States described a multi-institutional, open effort funded by the National Institutes of Health to assess model reproducibility in computational knee biomechanics. A regulatory representative from the United States Food and Drug Administration indicated the necessity of standards for reproducibility to increase utility of M&S in the regulatory setting. An industry representative from a major orthopedic implant company emphasized improving reproducibility by addressing indeterminacy in personalized modeling through sensitivity analyses, thereby enhancing preclinical evaluation of joint replacement technology. Thought leaders in the M&S community stressed the importance of data sharing to minimize duplication of efforts. A survey comprised 103 attendees revealed strong support for the workshop and for increasing emphasis on computational modeling at future ORS meetings. Nearly all survey respondents (97%) considered reproducibility to be an important issue. Almost half of respondents (45%) tried and failed to reproduce the work of others. Two-thirds of respondents (67%) declared that individual laboratories are most responsible for ensuring reproducible research whereas 44% thought that journals are most responsible. Thought leaders and survey respondents emphasized that computational models must be reproducible and credible to advance knee M&S.

Knee pivot location in asymptomatic older adults.

Representative data of asymptomatic, native-knee kinematics is important when studying changes in knee function across the lifespan. High-speed stereo radiography (HSSR) provides a reliable measure of knee kinematics to <1 mm of translation and 1° of rotation, but studies often have limited statistical power to make comparisons between groups or measure the contribution of individual variability. The purpose of this study is to examine in vivo condylar kinematics to quantify the transverse center-of-rotation, or pivot, location across the flexion range and challenge the medial-pivot paradigm in asymptomatic knee kinematics. We quantified the pivot location during supine leg press, knee extension, standing lunge, and gait for 53 middle-aged and older adults (27 men; 26 women: 50.8 ± 7.0 yrs, 1.75 ± 0.1 m, 79.1 ± 15.4 kg). A central- to medial-pivot location was identified for all activities with increased knee flexion associated with posterior translation of the center-of-rotation. The association between knee angle and anterior-posterior center-of-rotation location was not as strong as the relation between medial-lateral and anterior-posterior location, excluding gait. The Pearson's correlation for gait was stronger between knee angle and anterior-posterior center-of-rotation location (P < 0.001) than medial-lateral and anterior-posterior location (P = 0.0122). Individual variability accounted for a measurable proportion in variance explained of center-of-rotation location. Unique to gait, the lateral translation of center-of-rotation location resulted in the anterior translation of center-of-rotation at <10° knee flexion. Furthermore, no association between vertical ground-reaction force and center-of-rotation was identified.

Three Dimensional Lower Extremity Musculoskeletal Geometry of the Visible Human Female and Male.

Models and simulations of human function impact medicine and medical technology. Particularly, musculoskeletal modeling provides an avenue for insight into the human body, which might not be otherwise possible. However, reaching the ultimate goal of functional multi-scale human models has been slowed by the lack of freely available datasets of anatomical models and geometries. Moreover, female-specific geometries have been neglected with a widespread emphasis on male geometry. To help realize this goal, we have developed and shared complete three-dimensional musculoskeletal geometries extracted from the National Libraries of Medicine Visible Human Female and Male cryosections. Muscle, bone, cartilage, ligament, and fat from the pelvis to the ankle were digitized and exported. These geometries provide a foundation for continued work in human musculoskeletal simulation with high-fidelity deformable tissues that enable a better understanding of normal function and the evaluation of pathologies and treatments. This work is novel as it includes both the male and female Visible Human specimens, outputs at multiple levels of post-processing for maximum data reuse, and is publicly available.

Assessment of reporting practices and reproducibility potential of a cohort of published studies in computational knee biomechanics.

Reproducible research serves as a pillar of the scientific method and is a foundation for scientific advancement. However, estimates for irreproducibility of preclinical science range from 75% to 90%. The importance of reproducible science has not been assessed in the context of mechanics-based modeling of human joints such as the knee, despite this being an area that has seen dramatic growth. Framed in the context of five experienced teams currently documenting knee modeling procedures, the aim of this study was to evaluate reporting and the perceived potential for reproducibility across studies the teams viewed as important contributions to the literature. A cohort of studies was selected by polling, which resulted in an assessment of nine studies as opposed to a broader analysis across the literature. Using a published checklist for reporting of modeling features, the cohort was evaluated for both "reporting" and their potential to be "reproduced," which was delineated into six major modeling categories and three subcategories. Logistic regression analysis revealed that for individual modeling categories, the proportion of "reported" occurrences ranged from 0.31, 95% confidence interval (CI) [0.23, 0.41] to 0.77, 95% CI: [0.68, 0.86]. The proportion of whether a category was perceived as "reproducible" ranged from 0.22, 95% CI: [0.15, 0.31] to 0.44, 95% CI: [0.35, 0.55]. The relatively low ratios highlight an opportunity to improve reporting and reproducibility of knee modeling studies. Ongoing efforts, including our findings, contribute to a dialogue that facilitates adoption of practices that provide both credibility and translation possibilities.

Effects of Weight-Bearing on Tibiofemoral, Patellofemoral, and Patellar Tendon Kinematics in Older Adults.

Quantification of natural knee kinematics is essential for the assessment of joint function in the diagnosis of pathologies. Combined measurements of tibiofemoral and patellofemoral joint kinematics are necessary because knee pathologies, such as progression of osteoarthritis and patellar instability, are a frequent concern in both articulations. Combined measurement of tibiofemoral and patellofemoral kinematics also enables calculation of important quantities, specifically patellar tendon angle, which partly determines the loading vector at the tibiofemoral joint and patellar tendon moment arm. The goals of this research were to measure the differences in tibiofemoral and patellofemoral kinematics, patellar tendon angle (PTA), and patellar tendon moment arm (PTMA) that occur during non-weight-bearing and weight-bearing activities in older adults. METHODS: High-speed stereo radiography was used to measure the kinematics of the tibiofemoral and patellofemoral joints in subjects as they performed seated, non-weight-bearing knee extension and two weight-bearing activities: lunge and chair rise. PTA and PTMA were extracted from the subject's patellofemoral and tibiofemoral kinematics. Kinematics and the root mean square difference (RMSD) between non-weight-bearing and weight-bearing activities were compared across subjects and activities. RESULTS: Internal rotation increased with weight-bearing (mean RMSD from knee extension was 4.2 ± 2.4° for lunge and 3.6 ± 1.8° for chair rise), and anterior translation was also greater (mean RMSD from knee extension was 2.2 ± 1.2 mm for lunge and 2.3 ± 1.4 mm for chair rise). Patellar tilt and medial-lateral translation changed from non-weight-bearing to weight-bearing. Changes of the patellar tendon from non-weight-bearing to weight-bearing were significant only for PTMA. CONCLUSIONS: While weight-bearing elicited changes in knee kinematics, in most degrees of freedoms, these differences were exceeded by intersubject differences. These results provide comparative kinematics for the evaluation of knee pathology and treatment in older adults.

Supine leg press as an alternative to standing lunge in high-speed stereo radiography.

The standing lunge is an activity commonly used to quantify in-vivo knee kinematics with fluoroscopy. The ability to perform the standing lunge varies between subjects and can necessitate movement accommodations to successfully complete the desired range of motion. We proposed a supine leg press as an alternative to the standing lunge that aimed to provide a similar evaluation of knee motion while increasing the measured range of motion. Tibiofemoral kinematics of 53 non-symptomatic adults (27 men, 26 women, 50.8 ± 7.0 yrs.) were calculated from the tracked high-speed stereo radiography (HSSR) images for supine leg press and standing lunge using CT-segmented bony geometries of the right lower limb. The supine leg press proved to be a useful alternative to the standing lunge while providing 46.2° greater range of motion in knee flexion. The difference in angle-matched kinematics across a 100° flexion range between the leg press and lunge was 0.70° in varus-valgus rotation, 1.5° in internal-external rotation, 1.0 mm in medial-lateral translation, 2.3 mm in anterior-posterior translation, and 0.46 mm in superior-inferior translation for men. The angle-matched difference for women across 100° was 0.58° in varus-valgus rotation, 2.4° internal-external rotation, 0.70 mm medial-lateral translation, 2.1 mm anterior-posterior translation, and 0.78 mm superior-inferior translation. The similar kinematics, while having a greater range of motion, and control of the applied load makes the supine leg press an alternative for quantifying in-vivo knee kinematics.

Integration of neural architecture within a finite element framework for improved neuromusculoskeletal modeling.

Neuromusculoskeletal (NMS) models can aid in studying the impacts of the nervous and musculoskeletal systems on one another. These computational models facilitate studies investigating mechanisms and treatment of musculoskeletal and neurodegenerative conditions. In this study, we present a predictive NMS model that uses an embedded neural architecture within a finite element (FE) framework to simulate muscle activation. A previously developed neuromuscular model of a motor neuron was embedded into a simple FE musculoskeletal model. Input stimulation profiles from literature were simulated in the FE NMS model to verify effective integration of the software platforms. Motor unit recruitment and rate coding capabilities of the model were evaluated. The integrated model reproduced previously published output muscle forces with an average error of 0.0435 N. The integrated model effectively demonstrated motor unit recruitment and rate coding in the physiological range based upon motor unit discharge rates and muscle force output. The combined capability of a predictive NMS model within a FE framework can aid in improving our understanding of how the nervous and musculoskeletal systems work together. While this study focused on a simple FE application, the framework presented here easily accommodates increased complexity in the neuromuscular model, the FE simulation, or both.

Apparatus for In Vivo Knee Laxity Assessment Using High-Speed Stereo Radiography.

Computational modeling is of growing importance in orthopedics and biomechanics as a tool to understand differences in pathology and predict outcomes from surgical interventions. However, the computational models of the knee have historically relied on in vitro data to create and calibrate model material properties due to the unavailability of accurate in vivo data. This work demonstrates the design and use of a custom device to quantify anterior-posterior (AP) and internal-external (IE) in vivo knee laxity, with an accuracy similar to existing in vitro methods. The device uses high-speed stereo radiography (HSSR) tracking techniques to accurately measure the resulting displacements of the femur, tibia, and patella bones during knee laxity assessment at multiple loads and knee flexion angles. The accuracy of the knee laxity apparatus was determined by comparing laxity data from two cadaveric specimens between the knee laxity apparatus and an existing in vitro robotic knee joint simulator. The accuracy of the knee laxity apparatus was within 1 mm (0.04 in.) for AP and 2.5 deg for IE. Additionally, two living subjects completed knee laxity testing to confirm the laboratory use of the novel apparatus. This work demonstrates the ability to use custom devices in HSSR to collect accurate data, in vivo, for calibration of computational models.

Automatic tracking of healthy joint kinematics from stereo-radiography sequences.

Kinematic tracking of healthy joints in radiography sequences is frequently performed by maximizing similarities between computed perspective projections of 3D computer models and corresponding objects' appearances in radiographic images. Significant human effort associated with manual tracking presents a major bottleneck in biomechanics research methods and limits the scale of target applications. The current work introduces a method for fully-automatic tracking of tibiofemoral and patellofemoral kinematics in stereo-radiography sequences for subjects performing dynamic activities. The proposed method involves the application of convolutional neural networks for annotating radiographs and a multi-stage optimization pipeline for estimating bone pose based on information provided by neural net predictions. Predicted kinematics are evaluated by comparing against manually-tracked trends across 20 distinct trials. Median absolute differences below 1.5 millimeters or degrees for 6 tibiofemoral and 3 patellofemoral degrees of freedom demonstrate the utility of our approach, which improves upon previous semi-automatic methods by enabling end-to-end automation. Implementation of a fully-automatic pipeline for kinematic tracking will benefit evaluation of human movement by enabling large-scale studies of healthy knee kinematics.

A Method to Compare Heterogeneous Types of Bone and Cartilage Meshes.

Accurately capturing the bone and cartilage morphology and generating a mesh remains a critical step in the workflow of computational knee joint modeling. Currently, there is no standardized method to compare meshes of different element types and nodal densities, making comparisons across research teams a significant challenge. The aim of this paper is to describe a method to quantify differences in knee joint bone and cartilages meshes, independent of bone and cartilage mesh topology. Bone mesh-to-mesh distances, subchondral bone boundaries, and cartilage thicknesses from meshes of any type of mesh are obtained using a series of steps involving registration, resampling, and radial basis function fitting after which the comparisons are performed. Subchondral bone boundaries and cartilage thicknesses are calculated and visualized in a common frame of reference for comparison. The established method is applied to models developed by five modeling teams. Our approach to obtain bone mesh-to-mesh distances decreased the divergence seen in selecting a reference mesh (i.e., comparing mesh A-to-B versus mesh B-to-A). In general, the bone morphology was similar across teams. The cartilage thicknesses for all models were calculated and the mean absolute cartilage thickness difference was presented, the articulating areas had the best agreement across teams. The teams showed disagreement on the subchondral bone boundaries. The method presented in this paper allows for objective comparisons of bone and cartilage geometry that is agnostic to mesh type and nodal density.

In vivo comparison of rotating platform and fixed bearing knee replacements during lunge and pivot activities.

BACKGROUND: The objective of this study was to discover whether notable differences in mobile and fixed-bearing kinematics occur during activity that promotes tibial rotation, and to compare these results with normal healthy kinematics. We hypothesized that rotating-platform knee replacements would exhibit greater rotation of the tibia relative to the fixed-bearing knee replacements. MATERIALS AND METHODS: The in vivo motion of the tibia relative to the femur was measured in subjects with posterior stabilized fixed-bearing (FB) and rotating-platform (RP) total knee arthroplasties using a high-speed stereo radiography system during a lunge and gait with a change in direction (pivot). RESULTS: The in vivo internal/external (IE) rotation and anterior/posterior translation of the tibia relative to the femur was similar between mobile and fixed-bearing total knee prostheses during two activities of daily living that included an activity that challenged tibial IE rotation. Measurements of IE rotation in participants with RP had higher variability and significantly greater range between maximum internal and external rotation compared with FB participants. The greater amount of variability of RP was not unlike the healthy knee. CONCLUSION: The pattern of IE rotation and AP translation for both RP and FB designs were similar to healthy kinematics but with less IE rotation. The RP implants more closely replicated the asymmetrical posterior condylar translation and range of IE rotation of the healthy knee during activity that challenged tibial IE rotation.

Deciphering the "Art" in Modeling and Simulation of the Knee Joint: Variations in Model Development.

The use of computational modeling to investigate knee joint biomechanics has increased exponentially over the last few decades. Developing computational models is a creative process where decisions have to be made, subject to the modelers' knowledge and previous experiences, resulting in the "art" of modeling. The long-term goal of the KneeHub project is to understand the influence of subjective decisions on the final outcomes and the reproducibility of computational knee joint models. In this paper, we report on the model development phase of this project, investigating model development decisions and deviations from initial modeling plans. Five teams developed computational knee joint models from the same dataset, and we compared each teams' initial uncalibrated models and their model development workflows. Variations in the software tools and modeling approaches were found, resulting in differences such as the representation of the anatomical knee joint structures in the model. The teams consistently defined the boundary conditions and used the same anatomical coordinate system convention. However, deviations in the anatomical landmarks used to define the coordinate systems were present, resulting in a large spread in the kinematic outputs of the uncalibrated models. The reported differences and similarities in model development and simulation presented here illustrate the importance of the "art" of modeling and how subjective decision-making can lead to variation in model outputs. All teams deviated from their initial modeling plans, indicating that model development is a flexible process and difficult to plan in advance, even for experienced teams.

Patellofemoral kinematics in healthy older adults during gait activities.

The patellofemoral (PF) joint is susceptible to many pathologies resulting from acute injury, chronic disease and complications following surgical treatment of the knee. The objectives of this study were to describe case series measurements of patellar motion in healthy older adults as they performed three gait activities, determine patellar tendon angle and moment arm, and show if these quantities were activity dependent. A stereo radiography system was utilized to obtain the 3D PF kinematics of seventeen healthy people over 55 years of age (8F/9M, 66 ± 7.9 years old, 75.7 ± 20.5 kg) as they performed level walking, a step down, and a pivot turn. For a similar portion of the gait cycle, patellar flexion (6.2° ± 5.8) and average range of motion (ROM) (11.0° ± 5.9°) for walking with a step down was greater compared to the other gait activities (gait ROM 6.9° ± 4.3°, pivot ROM 5.7° ± 3.3°), while the average range of motion for patella tilt was greater during walking with a pivot turn (8.6° ± 3.9°). However, each subject displayed distinct PF kinematic trends during all activities with a few notable exceptions. Importantly, the knee extensor mechanism characteristics of patellar tendon angle and moment arm showed considerable variation across subjects but were largely unaltered by changing activities. The variation between subjects and the different behavior of the patella during the step down and pivot emphasized the need for analysis of a range of activities to reveal individual response to pathology and treatment in patellar maltracking and osteoarthritis.

Lumbar Intervertebral Kinematics During an Unstable Sitting Task and Its Association With Standing-Induced Low Back Pain.

People developing transient low back pain during standing have altered control of their spine and hips during standing tasks, but the transfer of these responses to other tasks has not been assessed. This study used video fluoroscopy to assess lumbar spine intervertebral kinematics of people who do and do not develop standing-induced low back pain during a seated chair-tilting task. A total of 9 females and 8 males were categorized as pain developers (5 females and 3 males) or nonpain developers (4 females and 5 males) using a 2-hour standing exposure; pain developers reported transient low back pain and nonpain developers did not. Participants were imaged with sagittal plane fluoroscopy at 25 Hz while cyclically tilting their pelvises anteriorly and posteriorly on an unstable chair. Intervertebral angles, relative contributions, and anterior-posterior translations were measured for the L3/L4, L4/L5, and L5/S1 joints and compared between sexes, pain groups, joints, and tilting directions. Female pain developers experienced more extension in their L5/S1 joints in both tilting directions compared with female nonpain developers, a finding not present in males. The specificity in intervertebral kinematics to sex-pain group combinations suggests that these subgroups of pain developers and nonpain developers may implement different control strategies.

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.