Subject-Specific Finite Element Models of the Tibia With Realistic Boundary Conditions Predict Bending Deformations Consistent With In Vivo Measurement

2019 ◽  
Vol 142 (2) ◽  
Author(s):  
Ifaz T. Haider ◽  
Michael Baggaley ◽  
W. Brent Edwards

Abstract Understanding the structural response of bone during locomotion may help understand the etiology of stress fracture. This can be done in a subject-specific manner using finite element (FE) modeling, but care is needed to ensure that modeling assumptions reflect the in vivo environment. Here, we explored the influence of loading and boundary conditions (BC), and compared predictions to previous in vivo measurements. Data were collected from a female participant who walked/ran on an instrumented treadmill while motion data were captured. Inverse dynamics of the leg (foot, shank, and thigh segments) was combined with a musculoskeletal (MSK) model to estimate muscle and joint contact forces. These forces were applied to an FE model of the tibia, generated from computed tomography (CT). Eight conditions varying loading/BCs were investigated. We found that modeling the fibula was necessary to predict realistic tibia bending. Applying joint moments from the MSK model to the FE model was also needed to predict torsional deformation. During walking, the most complex model predicted deformation of 0.5 deg posterior, 0.8 deg medial, and 1.4 deg internal rotation, comparable to in vivo measurements of 0.5–1 deg, 0.15–0.7 deg, and 0.75–2.2 deg, respectively. During running, predicted deformations of 0.3 deg posterior, 0.3 deg medial, and 0.5 deg internal rotation somewhat underestimated in vivo measures of 0.85–1.9 deg, 0.3–0.9 deg, 0.65–1.72 deg, respectively. Overall, these models may be sufficiently realistic to be used in future investigations of tibial stress fracture.

2019 ◽  
Vol 48 (4) ◽  
pp. 1181-1195 ◽  
Author(s):  
Yuyang Wei ◽  
Zhenmin Zou ◽  
Guowu Wei ◽  
Lei Ren ◽  
Zhihui Qian

AbstractThis paper aims to develop and validate a subject-specific framework for modelling the human hand. This was achieved by combining medical image-based finite element modelling, individualized muscle force and kinematic measurements. Firstly, a subject-specific human hand finite element (FE) model was developed. The geometries of the phalanges, carpal bones, wrist bones, ligaments, tendons, subcutaneous tissue and skin were all included. The material properties were derived from in-vivo and in-vitro experiment results available in the literature. The boundary and loading conditions were defined based on the kinematic data and muscle forces of a specific subject captured from the in-vivo grasping tests. The predicted contact pressure and contact area were in good agreement with the in-vivo test results of the same subject, with the relative errors for the contact pressures all being below 20%. Finally, sensitivity analysis was performed to investigate the effects of important modelling parameters on the predictions. The results showed that contact pressure and area were sensitive to the material properties and muscle forces. This FE human hand model can be used to make a detailed and quantitative evaluation into biomechanical and neurophysiological aspects of human hand contact during daily perception and manipulation. The findings can be applied to the design of the bionic hands or neuro-prosthetics in the future.


Author(s):  
Jason P. Halloran ◽  
Anthony J. Petrella ◽  
Paul J. Rullkoetter

The success of current total knee replacement (TKR) devices is contingent on the kinematics and contact mechanics during in vivo activity. Indicators of potential clinical performance of total joint replacement devices include contact stress and area due to articulations, and tibio-femoral and patello-femoral kinematics. An effective way of evaluating these parameters during the design phase or before clinical use is via computationally efficient computer models. Previous finite element (FE) knee models have generally been used to determine contact stresses and/or areas during static or quasi-static loading conditions. The majority of knee models intended to predict relative kinematics have not been able to determine contact mechanics simultaneously. Recently, however, explicit dynamic finite element methods have been used to develop dynamic models of TKR able to efficiently determine joint and contact mechanics during dynamic loading conditions [1,2]. The objective of this research was to develop and validate an explicit FE model of a TKR which includes tibio-femoral and patello-femoral articulations and surrounding soft tissues. The six degree-of-freedom kinematics, kinetics and polyethylene contact mechanics during dynamic loading conditions were then predicted during gait simulation.


2014 ◽  
Vol 136 (4) ◽  
Author(s):  
Robert E. Carey ◽  
Liying Zheng ◽  
Ameet K. Aiyangar ◽  
Christopher D. Harner ◽  
Xudong Zhang

In this paper, we present a new methodology for subject-specific finite element modeling of the tibiofemoral joint based on in vivo computed tomography (CT), magnetic resonance imaging (MRI), and dynamic stereo-radiography (DSX) data. We implemented and compared two techniques to incorporate in vivo skeletal kinematics as boundary conditions: one used MRI-measured tibiofemoral kinematics in a nonweight-bearing supine position and allowed five degrees of freedom (excluding flexion-extension) at the joint in response to an axially applied force; the other used DSX-measured tibiofemoral kinematics in a weight-bearing standing position and permitted only axial translation in response to the same force. Verification and comparison of the model predictions employed data from a meniscus transplantation study subject with a meniscectomized and an intact knee. The model-predicted cartilage-cartilage contact areas were examined against “benchmarks” from a novel in situ contact area analysis (ISCAA) in which the intersection volume between nondeformed femoral and tibial cartilage was characterized to determine the contact. The results showed that the DSX-based model predicted contact areas in close alignment with the benchmarks, and outperformed the MRI-based model: the contact centroid predicted by the former was on average 85% closer to the benchmark location. The DSX-based FE model predictions also indicated that the (lateral) meniscectomy increased the contact area in the lateral compartment and increased the maximum contact pressure and maximum compressive stress in both compartments. We discuss the importance of accurate, task-specific skeletal kinematics in subject-specific FE modeling, along with the effects of simplifying assumptions and limitations.


2014 ◽  
Vol 136 (2) ◽  
Author(s):  
Trent M. Guess ◽  
Antonis P. Stylianou ◽  
Mohammad Kia

Detailed knowledge of knee kinematics and dynamic loading is essential for improving the design and outcomes of surgical procedures, tissue engineering applications, prosthetics design, and rehabilitation. This study used publicly available data provided by the “Grand Challenge Competition to Predict in-vivo Knee Loads” for the 2013 American Society of Mechanical Engineers Summer Bioengineering Conference (Fregly et al., 2012, “Grand Challenge Competition to Predict in vivo Knee Loads,” J. Orthop. Res., 30, pp. 503–513) to develop a full body, musculoskeletal model with subject specific right leg geometries that can concurrently predict muscle forces, ligament forces, and knee and ground contact forces. The model includes representation of foot/floor interactions and predicted tibiofemoral joint loads were compared to measured tibial loads for two different cycles of treadmill gait. The model used anthropometric data (height and weight) to scale the joint center locations and mass properties of a generic model and then used subject bone geometries to more accurately position the hip and ankle. The musculoskeletal model included 44 muscles on the right leg, and subject specific geometries were used to create a 12 degrees-of-freedom anatomical right knee that included both patellofemoral and tibiofemoral articulations. Tibiofemoral motion was constrained by deformable contacts defined between the tibial insert and femoral component geometries and by ligaments. Patellofemoral motion was constrained by contact between the patellar button and femoral component geometries and the patellar tendon. Shoe geometries were added to the feet, and shoe motion was constrained by contact between three shoe segments per foot and the treadmill surface. Six-axis springs constrained motion between the feet and shoe segments. Experimental motion capture data provided input to an inverse kinematics stage, and the final forward dynamics simulations tracked joint angle errors for the left leg and upper body and tracked muscle length errors for the right leg. The one cycle RMS errors between the predicted and measured tibia contact were 178 N and 168 N for the medial and lateral sides for the first gait cycle and 209 N and 228 N for the medial and lateral sides for the faster second gait cycle. One cycle RMS errors between predicted and measured ground reaction forces were 12 N, 13 N, and 65 N in the anterior-posterior, medial-lateral, and vertical directions for the first gait cycle and 43 N, 15 N, and 96 N in the anterior-posterior, medial-lateral, and vertical directions for the second gait cycle.


Author(s):  
Mozammil Hussain ◽  
Raghu N. Natarajan ◽  
Gunnar B. J. Andersson ◽  
Howard S. An

Degenerative changes in the cervical spine due to aging are very common causes of neck pain in general population. Although many investigators have quantified the gross morphological changes in the disc with progressive degeneration, the biomechanical changes due to degenerative pathologies of the disc and its effect on the adjacent levels are not well understood. Despite many in vivo and in vitro techniques used to study such complex phenomena, the finite element (FE) method is still a powerful tool to investigate the internal mechanics and complex clinical situations under various physiological loadings particularly when large numbers of parameters are involved. The objective of the present study was to develop and validate a poroelastic FE model of a healthy C3-T1 segment of the cervical spine under physiologic moment loads. The model included the regional effect of change in the fixed charged density of proteoglycan concentration and change in the permeability and porosity due to change in the axial strain of disc tissues. The model was further modified to include various degrees of disc degeneration at the C5-C6 level. Outcomes of this study provided a better understanding on the progression of degeneration along the cervical spine by investigating the biomechanical response of the adjacent segments with an intermediate degenerated C5-C6 level.


2016 ◽  
Vol 138 (2) ◽  
Author(s):  
Yihwan Jung ◽  
Cong-Bo Phan ◽  
Seungbum Koo

Joint contact forces measured with instrumented knee implants have not only revealed general patterns of joint loading but also showed individual variations that could be due to differences in anatomy and joint kinematics. Musculoskeletal human models for dynamic simulation have been utilized to understand body kinetics including joint moments, muscle tension, and knee contact forces. The objectives of this study were to develop a knee contact model which can predict knee contact forces using an inverse dynamics-based optimization solver and to investigate the effect of joint constraints on knee contact force prediction. A knee contact model was developed to include 32 reaction force elements on the surface of a tibial insert of a total knee replacement (TKR), which was embedded in a full-body musculoskeletal model. Various external measurements including motion data and external force data during walking trials of a subject with an instrumented knee implant were provided from the Sixth Grand Challenge Competition to Predict in vivo Knee Loads. Knee contact forces in the medial and lateral portions of the instrumented knee implant were also provided for the same walking trials. A knee contact model with a hinge joint and normal alignment could predict knee contact forces with root mean square errors (RMSEs) of 165 N and 288 N for the medial and lateral portions of the knee, respectively, and coefficients of determination (R2) of 0.70 and −0.63. When the degrees-of-freedom (DOF) of the knee and locations of leg markers were adjusted to account for the valgus lower-limb alignment of the subject, RMSE values improved to 144 N and 179 N, and R2 values improved to 0.77 and 0.37, respectively. The proposed knee contact model with subject-specific joint model could predict in vivo knee contact forces with reasonable accuracy. This model may contribute to the development and improvement of knee arthroplasty.


2006 ◽  
Vol 39 ◽  
pp. S134 ◽  
Author(s):  
G. Bergmann ◽  
F. Graichen ◽  
A. Bender ◽  
M. Kääb ◽  
A. Rohlmann ◽  
...  

2018 ◽  
Vol 140 (7) ◽  
Author(s):  
Swithin S. Razu ◽  
Trent M. Guess

Computational models that predict in vivo joint loading and muscle forces can potentially enhance and augment our knowledge of both typical and pathological gaits. To adopt such models into clinical applications, studies validating modeling predictions are essential. This study created a full-body musculoskeletal model using data from the “Sixth Grand Challenge Competition to Predict in vivo Knee Loads.” This model incorporates subject-specific geometries of the right leg in order to concurrently predict knee contact forces, ligament forces, muscle forces, and ground contact forces. The objectives of this paper are twofold: (1) to describe an electromyography (EMG)-driven modeling methodology to predict knee contact forces and (2) to validate model predictions by evaluating the model predictions against known values for a patient with an instrumented total knee replacement (TKR) for three distinctly different gait styles (normal, smooth, and bouncy gaits). The model integrates a subject-specific knee model onto a previously validated generic full-body musculoskeletal model. The combined model included six degrees-of-freedom (6DOF) patellofemoral and tibiofemoral joints, ligament forces, and deformable contact forces with viscous damping. The foot/shoe/floor interactions were modeled by incorporating shoe geometries to the feet. Contact between shoe segments and the floor surface was used to constrain the shoe segments. A novel EMG-driven feedforward with feedback trim motor control strategy was used to concurrently estimate muscle forces and knee contact forces from standard motion capture data collected on the individual subject. The predicted medial, lateral, and total tibiofemoral forces represented the overall measured magnitude and temporal patterns with good root-mean-squared errors (RMSEs) and Pearson's correlation (p2). The model accuracy was high: medial, lateral, and total tibiofemoral contact force RMSEs = 0.15, 0.14, 0.21 body weight (BW), and (0.92 < p2 < 0.96) for normal gait; RMSEs = 0.18 BW, 0.21 BW, 0.29 BW, and (0.81 < p2 < 0.93) for smooth gait; and RMSEs = 0.21 BW, 0.22 BW, 0.33 BW, and (0.86 < p2 < 0.95) for bouncy gait, respectively. Overall, the model captured the general shape, magnitude, and temporal patterns of the contact force profiles accurately. Potential applications of this proposed model include predictive biomechanics simulations, design of TKR components, soft tissue balancing, and surgical simulation.


2019 ◽  
Author(s):  
◽  
Swithin Samuel Razu

"The goal of this dissertation is to develop a musculoskeletal model and corroborate model predictions to experimentally measured in vivo knee contact forces, in order to study the biomechanical consequences of two different total knee arthroplasty designs. The two main contributions of this dissertation are: (1) Corroboration to experimental data: The development of an EMG-driven, full-body, musculoskeletal model with subject-specific leg geometries including deformable contacts, ligaments, 6DOF knee joint, and a shoe-floor model that can concurrently predict muscle forces, ligament forces, and joint contact forces. Model predictions of tibiofemoral joint contact forces were evaluated against the subject-specific in vivo measurements from the instrumented TKR for three distinctly different styles of over ground gait. (2) Virtual surgery in TKA: The musculoskeletal modeling methodology was then used to develop a model for one healthy participant with a native knee and then virtually replacing the native knee with fixed-bearing and mobile-bearing total knee arthroplasty designs performing gait and step-up tasks. This approach minimized the biomechanical impact of variations in sex, geometry, implant size, design and positioning, ligament location and tension, and muscle forces found across patients. The differences in biomechanics were compared for the two designs. 1.2 Significance of this Research The world health organization ranks musculoskeletal disorders as the second largest contributor to disability worldwide. Conservative estimates put the national cost of direct care for musculoskeletal disease at $212.7 billion a year [1]. Many people who suffer from neuromuscular or musculoskeletal diseases may benefit from the insights gained from surgery simulations, since musculoskeletal reconstructions are commonly performed on these individuals. Improved surgical outcomes will benefit these individuals not only in the short-term, but also in the long-term, since their future rehabilitation needs may be reduced. For example, although total knee arthroplasty is a common surgical procedure for the treatment of osteoarthritis with over 700,000 procedures performed each year [2], many patients are unhappy with the ultimate results [3]. Ten to 30% of patients report [4] pain, dissatisfaction with function, and the need for further surgery such as revision after the initial surgery resulting in costs exceeding $11 billion [5]. Potentially, simulation studies that quantify the important biomechanical variables will reduce the need for revision surgeries in patients."--Introduction.


Sign in / Sign up

Export Citation Format

Share Document