Specimen-Specific Natural, Pathological, and Implanted Knee Mechanics Using Finite Element Modeling

There is an increasing incidence of knee pain and injury among the population, and increasing demand for higher knee function in total knee replacement designs. As a result, clinicians and implant manufacturers are interested in improving patient outcomes, and evaluation of knee mechanics is essential for better diagnosis and repair of knee pathologies. Common knee pathologies include osteoarthritis (degradation of the articulating surfaces), patellofemoral pain, and cruciate ligament injury and/or rupture. The complex behavior of knee motion presents unique challenges in the diagnosis of knee pathology and restoration of healthy knee function. Quantifying knee mechanics is essential for developing successful rehabilitation therapies and surgical treatments. Researchers have used in-vitro and in-vivo experiments to quantify joint kinematics and loading, but experiments can be costly and time-intensive, and contact and ligament mechanics can be difficult to measure directly. Computational modeling can complement experimental studies by providing cost-effective solutions for quantifying joint and soft tissue forces. Musculoskeletal models have been used to measure whole-body motion, and predict joint and muscle forces, but these models can lack detail and accuracy at the joint-level. Finite element modeling provides accurate solutions of the internal stress/strain behavior of bone and soft tissue using subject-specific geometry and complex contact and material representations. While previous FE modeling has been used to simulate injury and repair, models are commonly based on literature description or average knee behavior. The research presented in this dissertation focused on developing subject-specific representations of the TF and PF joints including calibration and validation to experimental data for healthy, pathological, and implanted knee conditions. A combination of in-vitro experiment and modeling was used to compare healthy and cruciate-deficient joint mechanics, and develop subject-specific computational representations. Insight from in-vitro testing supported in-vivo simulations of healthy and implanted subjects, in which PF mechanics were compared between two common patellar component designs and the impact of cruciate ligament variability on joint kinematics and loads was assessed. The suite of computational models developed in this dissertation can be used to investigate knee pathologies to better inform clinicians on the mechanisms surrounding injury, support the diagnosis of at-risk patients, explore rehabilitation and surgical techniques for repair, and support decision-making for new innovative implant designs

Development and Validation of a Tibiofemoral Joint Finite Element Model and Subsequent Gait Analysis of Intact ACL and ACL Deficient Individuals

Osteoarthritis (OA) is a degenerative condition of articular cartilage that affects more than 25 million people in the US. Joint injuries, like anterior cruciate ligament (ACL) tears, can lead to OA due to a change in articular cartilage loading. Gait analysis combined with knee joint finite element modeling (FEM) has been used to predict the articular cartilage loading. To predict the change of articular cartilage loading during gait due to various ACL injuries, a tibiofemoral FEM was developed from magnetic resonance images (MRIs) of a 33 year male, with no prior history of knee injuries. The FEM was validated for maximum contact pressure and anterior tibial translation using cadaver knee studies. The FEM was used to model gait of knees with an intact ACL, anteromedial (AM) bundle injury, posterolateral (PL) bundle injury, complete ACL injury, AM deficiency, PL deficiency, complete ACL rupture, as well as a bone-patellar tendon-bone (BPTB) graft. Generally, the predicted maximum contact pressure and contact area increased for all the ACL injuries when compared to intact ACLs. While an increase in maximum contact pressure and contact area is an indication of an increased risk of the development of OA, the percent of increase was typically small suggesting that walking is a safe activity for individuals with ACL injuries

Subject Specific Computational Models of the Knee to Predict Anterior Cruciate Ligament Injury

Knee joint is a complex joint involving multiple interactions between cartilage, bone, muscles, ligaments, tendons and neural control. Anterior Cruciate Ligament (ACL) is one ligament in the knee joint that frequently gets injured during various sports or recreational activities. ACL injuries are common in college level and professional athletes especially in females and the injury rate is growing in epidemic proportions despite significant increase in the research focusing on neuromuscular and proprioceptive training programs. Most ACL injuries lead to surgical reconstruction followed by a lengthy rehabilitation program impacting the health and performance of the athlete. Furthermore, the athlete is still at the risk of early onset of osteoarthritis. Regardless of the gender disparity in the ACL injury rates, a clear understanding of the underlying injury mechanisms is required in order to reduce the incidence of these injuries. Computational modeling is a resourceful and cost effective tool to investigate the biomechanics of the knee. The aim of this study was twofold. The first aim was to develop subject specific computational models of the knee joint and the second aim to gain an improved understanding of the ACL injury mechanisms using the subject specific models. We used a quasi-static, multi-body modeling approach and developed MRI based tibio-femoral computational knee joint models. Experimental joint laxity and combined loading data was obtained using five cadaveric knee specimens and a state-of-the-art robotic system. Ligament zero strain lengths and insertion points were optimized using joint laxity data. Combined loading and ACL strain data were used for model validations. ACL injury simulations were performed using factorial design approach comprising of multiple factors and levels to replicate a large and rich set of loading states. This thesis is an extensive work covering all the details of the ACL injury project explained above and highlighting the importance of 1) computational modeling in inj

Subject Specific Computational Models of the Knee to Predict Anterior Cruciate Ligament Injury

Knee joint is a complex joint involving multiple interactions between cartilage, bone, muscles, ligaments, tendons and neural control. Anterior Cruciate Ligament (ACL) is one ligament in the knee joint that frequently gets injured during various sports or recreational activities. ACL injuries are common in college level and professional athletes especially in females and the injury rate is growing in epidemic proportions despite significant increase in the research focusing on neuromuscular and proprioceptive training programs. Most ACL injuries lead to surgical reconstruction followed by a lengthy rehabilitation program impacting the health and performance of the athlete. Furthermore, the athlete is still at the risk of early onset of osteoarthritis. Regardless of the gender disparity in the ACL injury rates, a clear understanding of the underlying injury mechanisms is required in order to reduce the incidence of these injuries. Computational modeling is a resourceful and cost effective tool to investigate the biomechanics of the knee. The aim of this study was twofold. The first aim was to develop subject specific computational models of the knee joint and the second aim to gain an improved understanding of the ACL injury mechanisms using the subject specific models. We used a quasi-static, multi-body modeling approach and developed MRI based tibio-femoral computational knee joint models. Experimental joint laxity and combined loading data was obtained using five cadaveric knee specimens and a state-of-the-art robotic system. Ligament zero strain lengths and insertion points were optimized using joint laxity data. Combined loading and ACL strain data were used for model validations. ACL injury simulations were performed using factorial design approach comprising of multiple factors and levels to replicate a large and rich set of loading states. This thesis is an extensive work covering all the details of the ACL injury project explained above and highlighting the importance of 1) computational modeling in inj

Articular human joint modelling

Copyright @ Cambridge University Press 2009.The work reported in this paper encapsulates the theories and algorithms developed to drive the core analysis modules of the software which has been developed to model a musculoskeletal structure of anatomic joints. Due to local bone surface and contact geometry based joint kinematics, newly developed algorithms make the proposed modeller different from currently available modellers. There are many modellers that are capable of modelling gross human body motion. Nevertheless, none of the available modellers offer complete elements of joint modelling. It appears that joint modelling is an extension of their core analysis capability, which, in every case, appears to be musculoskeletal motion dynamics. It is felt that an analysis framework that is focused on human joints would have significant benefit and potential to be used in many orthopaedic applications. The local mobility of joints has a significant influence in human motion analysis, in understanding of joint loading, tissue behaviour and contact forces. However, in order to develop a bone surface based joint modeller, there are a number of major problems, from tissue idealizations to surface geometry discretization and non-linear motion analysis. This paper presents the following: (a) The physical deformation of biological tissues as linear or non-linear viscoelastic deformation, based on spring-dashpot elements. (b) The linear dynamic multibody modelling, where the linear formulation is established for small motions and is particularly useful for calculating the equilibrium position of the joint. This model can also be used for finding small motion behaviour or loading under static conditions. It also has the potential of quantifying the joint laxity. (c) The non-linear dynamic multibody modelling, where a non-matrix and algorithmic formulation is presented. The approach allows handling complex material and geometrical nonlinearity easily. (d) Shortest path algorithms for calculating soft tissue line of action geometries. The developed algorithms are based on calculating minimum ‘surface mass’ and ‘surface covariance’. An improved version of the ‘surface covariance’ algorithm is described as ‘residual covariance’. The resulting path is used to establish the direction of forces and moments acting on joints. This information is needed for linear or non-linear treatment of the joint motion. (e) The final contribution of the paper is the treatment of the collision. In the virtual world, the difficulty in analysing bodies in motion arises due to body interpenetrations. The collision algorithm proposed in the paper involves finding the shortest projected ray from one body to the other. The projection of the body is determined by the resultant forces acting on it due to soft tissue connections under tension. This enables the calculation of collision condition of non-convex objects accurately. After the initial collision detection, the analysis involves attaching special springs (stiffness only normal to the surfaces) at the ‘potentially colliding points’ and motion of bodies is recalculated. The collision algorithm incorporates the rotation as well as translation. The algorithm continues until the joint equilibrium is achieved. Finally, the results obtained based on the software are compared with experimental results obtained using cadaveric joints

Two-dimensional dynamics model of the lower limb to include viscoelastic knee ligaments

A dynamic, 2D, anatomical knee joint model has been developed to simulate knee reactions to external input forces. A deformable contact area approach is used to find contact forces and moments, and a method of applying nonlinear viscoelastic ligament strain rate response was also developed and implemented on the model to account for the effects of viscoelasticity on the ligament fibers. The ligaments were then tested for various deficiencies to identify their effects on the natural frequency of the knee. Internal knee forces from ligaments, muscles, and contacting surfaces are modeled and then numerically found for different exercises. Static and dynamic equations for knee motion are developed. These equations are then transformed into differential algebraic equation (DAE) systems for modeling various exercises. The DAE systems and the model simulations are performed using Matlab solver ODE15S, and predicted data from the model is compared to data published in literature for validation

Development of a Human Tibiofemoral Joint Finite Element Model to Investigate the Effects of Obesity and Knee Malalignment on Joint Contact Pressure

Obesity is a known risk factor for osteoarthritis (OA). Excess body weight generates greater joint contact forces at the knee; however, obese individuals alter their gait to decrease joint contact forces. Knee malalignment has been identified as a strong mediating factor between obesity and knee OA progression. Excess body weight acting on a varus malaligned knee would have an additive effect on cartilage stress and could cause stress levels to exceed the threshold limit for damage and loss of cartilage matrix. A finite element (FE) model of the human tibiofemoral joint was developed and validated in order to investigate changes in cartilage pressure due to obesity and knee varus malalignment. The results of this analysis show that obese loading conditions caused greater contact pressure in both the lateral and medial tibiofemoral compartments at most phases of stance. Increased contact pressure applied cyclically during daily activities could make obese individuals more susceptible to OA. Varus malalignment increased medial contact pressure as expected, but lateral contact pressure also increased during midstance for both normal weight and obese load conditions. These results suggest that varus malaligned individuals could be susceptible to OA development in both tibiofemoral compartments due to the overall increase in joint contact pressure. As a qualitative tool, the FE model functioned well in highlighting changes in joint contact pressure due to the addition of obesity or varus malalignment. Further work can be done to increase confidence in the quantitative outputs of the model by using more sophisticated material models for soft tissue structures and incorporating the patellofemoral joint into the FE model

Development of a planar multi-body model of the human knee joint

The aim of this work is to develop a dynamic model for the biological human knee joint. The model is formulated in the framework of multibody systems methodologies, as a system of two bodies, the femur and the tibia. For the purpose of describing the formulation, the relative motion of the tibia with respect to the femur is considered. Due to their higher stiffness compared to that of the articular cartilages, the femur and tibia are considered as rigid bodies. The femur and tibia cartilages are considered to be deformable structures with specific material characteristics. The rotation and gliding motions of the tibia relative to the femur can not be modeled with any conventional kinematic joint, but rather in terms of the action of the knee ligaments and potential contact between the bones. Based on medical imaging techniques, the femur and tibia profiles in the sagittal plane are extracted and used to define the interface geometric conditions for contact. When a contact is detected, a continuous non-linear contact force law is applied which calculates the contact forces developed at the interface as a function of the relative indentation between the two bodies. The four basic cruciate and collateral ligaments present in the knee are also taken into account in the proposed knee joint model, which are modeled as non-linear elastic springs. The forces produced in the ligaments, together with the contact forces, are introduced into the system’s equations of motion as external forces. In addition, an external force is applied on the center of mass of the tibia, in order to actuate the system mimicking a normal gait motion. Finally, numerical results obtained from computational simulations are used to address the assumptions and procedures adopted in this study.Fundação para a Ciência e a Tecnologia (FCT

Development and validation of a computational model of the knee joint for the evaluation of surgical treatments for osteoarthritis

A three-dimensional (3D) knee joint computational model was developed and validated to predict knee joint contact forces and pressures for different degrees of malalignment. A 3D computational knee model was created from high-resolution radiological images to emulate passive sagittal rotation (full-extension to 658-flexion) and weight acceptance. A cadaveric knee mounted on a six-degree-of-freedom robot was subjected to matching boundary and loading conditions. A ligamenttuning process minimised kinematic differences between the robotically loaded cadaver specimen and the finite element (FE) model. The model was validated by measured intra-articular force and pressure measurements. Percent full scale error between FE-predicted and in vitro-measured values in the medial and lateral compartments were 6.67% and 5.94%, respectively, for normalised peak pressure values, and 7.56% and 4.48%, respectively, for normalised force values. The knee model can accurately predict normalised intra-articular pressure and forces for different loading conditions and could be further developed for subject-specific surgical planning