Development and validation of a weight-bearing finite element model for total knee replacement

Total knee arthroplasty (TKA) is a successful procedure for osteoarthritis. However, some patients (19%) do have pain after surgery. A finite element model was developed based on boundary conditions of a knee rig. A 3D-model of an anatomical full leg was generated from magnetic resonance image data and a total knee prosthesis was implanted without patella resurfacing. In the finite element model, a restarting procedure was programmed in order to hold the ground reaction force constant with an adapted quadriceps muscle force during a squat from 20° to 105° of flexion. Knee rig experimental data were used to validate the numerical model in the patellofemoral and femorotibial joint. Furthermore, sensitivity analyses of Young’s modulus of the patella cartilage, posterior cruciate ligament (PCL) stiffness, and patella tendon origin were performed. Pearson’s correlations for retropatellar contact area, pressure, patella flexion, and femorotibial ap-movement were near to 1. Lowest root mean square error for retropatellar pressure, patella flexion, and femorotibial ap-movement were found for the baseline model setup with Young’s modulus of 5 MPa for patella cartilage, a downscaled PCL stiffness of 25% compared to the literature given value and an anatomical origin of the patella tendon. The results of the conducted finite element model are comparable with the experimental results. Therefore, the finite element model developed in this study can be used for further clinical investigations and will help to better understand the clinical aspects after TKA with an unresurfaced patella.     

Development and validation of a subject-specific finite element model of a human clavicle

This study aimed to develop and validate a finite element (FE) model of a human clavicle which can predict the structural response and bone fractures under both axial compression and anterior–posterior three-point bending loads. Quasi-static non-injurious axial compression and three-point bending tests were first conducted on a male clavicle followed by a dynamic three-point bending test to fracture. Then, two types of FE models of the clavicle were developed using bone material properties which were set to vary with the computed tomography image density of the bone. A volumetric solid FE model comprised solely of hexahedral elements was first developed. A solid-shell FE model was then created which modelled the trabecular bone as hexahedral elements and the cortical bone as quadrilateral shell elements. Finally, simulations were carried out using these models to evaluate the influence of variations in cortical thickness, mesh density, bone material properties and modelling approach on the biomechanical responses of the clavicle, compared with experimental data. The FE results indicate that the inclusion of density-based bone material properties can provide a more accurate reproduction of the force–displacement response and bone fracture timing than a model with uniform bone material properties. Inclusion of a variable cortical thickness distribution also slightly improves the ability of the model to predict the experimental response. The methods developed in this study will be useful for creating subject-specific FE models to better understand the biomechanics and injury mechanism of the clavicle.     

Development and validation of a finite element model to simulate the opening of a medial opening wedge high tibial osteotomy

Medial opening wedge high tibial osteotomy (MOWHTO) is a surgical procedure intended to alter the coronal and sagittal plane alignment of the lower limb to primarily relieve the symptoms of osteoarthritis in the medial compartment of the knee. The purpose of this work was to develop and validate a finite element model to simulate the opening of a high tibial osteotomy and determine whether a pilot hole at the cortical hinge reduces the risk of lateral cortical fracture. Fifteen models were reconstructed from CT images of eight cadaveric specimens. The validated models indicated that the addition of the pilot hole increased the stresses and likelihood of a type-I and type-II fractures during the opening of a medial open wedge high tibial osteotomy compared to the no-hole condition.     

Physical validation of a patient-specific contact finite element model of the ankle

A validation study was conducted to determine the extent to which computational ankle contact finite element (FE) results agreed with experimentally measured tibio-talar contact stress. Two cadaver ankles were loaded in separate test sessions, during which ankle contact stresses were measured with a high-resolution (Tekscan) pressure sensor. Corresponding contact FE analyses were subsequently performed for comparison. The agreement was good between FE-computed and experimentally measured mean (3.2% discrepancy for one ankle, 19.3% for the other) and maximum (1.5% and 6.2%) contact stress, as well as for contact area (1.7% and 14.9%). There was also excellent agreement between histograms of fractional areas of cartilage experiencing specific ranges of contact stress. Finally, point-by-point comparisons between the computed and measured contact stress distributions over the articular surface showed substantial agreement, with correlation coefficients of 90% for one ankle and 86% for the other. In the past, general qualitative, but little direct quantitative agreement has been demonstrated with articular joint contact FE models. The methods used for this validation enable formal comparison of computational and experimental results, and open the way for objective statistical measures of regional correlation between FE-computed contact stress distributions from comparison articular joint surfaces (e.g., those from an intact versus those with residual intra-articular fracture incongruity).     

Lumbar spine finite element model for healthy subjects: development and validation

Finite element (FE) method is a proven powerful and efficient tool to study the biomechanics of the human lumbar spine. However, due to the large inter-subject variability of geometries and material properties in human lumbar spines, concerns existed on the accuracy and predictive power of one single deterministic FE model with one set of spinal geometry and material properties. It was confirmed that the combined predictions (median or mean value) of several distinct FE models can be used as an improved prediction of behavior of human lumbar spine under identical loading and boundary conditions. In light of this improved prediction, five FE models (L1-L5 spinal levels) of the human lumbar spine were developed based on five healthy living subjects with identical modeling method. The five models were extensively validated through experimental and computational results in the literature. Mesh convergence and material sensitivity analysis were also conducted. We have shown that the results from the five FE models developed in this paper were consistent with the experimental data and simulation results from the existing literature. The validated modeling method introduced in this study can be used in modeling dysfunctional lumber spines such as disc degeneration and scoliosis in future work.     

Accommodation of the human lens capsule using a finite element model based on nonlinear regionally anisotropic biomembranes

Accommodation of the eyes, the mechanism that allows humans to focus their vision on near objects, naturally diminishes with age via presbyopia. People who have undergone cataract surgery, using current surgical methods and artificial lens implants, are also left without the ability to accommodate. The process of accommodation is generally well known; however the specific mechanical details have not been adequately explained due to difficulties and consequences of performing in vivo studies. Most studies have modeled the mechanics of accommodation under assumptions of a linearly elastic, isotropic, homogenous lens and lens capsule. Recent experimental and numerical studies showed that the lens capsule exhibits nonlinear elasticity and regional anisotropy. In this paper we present a numerical model of human accommodation using a membrane theory based finite element approach, incorporating recent findings on capsular properties. This study seeks to provide a novel perspective of the mechanics of accommodation. Such findings may prove significant in seeking biomedical solutions to restoring loss of visual power.     

Implementation of reflex loops in a biomechanical finite element model

In the field of biomechanics, the offer of models which are more and more realistic requires to integrate a physiological response, in particular, the controlled muscle bracing and the reflexes. The following work aims to suggest a unique methodology which couples together a sensory and motor loop with a finite element model. Our method is applied to the study of the oscillation of the elbow in the case of a biceps brachial stretch reflex. The results obtained are promising in the purpose of the development of reactive human body models.     

Finite element models of total shoulder replacement: Application of boundary conditions

We present a one-dimensional (1D) fluid dynamic model that can predict blood flow and blood pressure during exercise using data collected at rest. To facilitate accurate prediction of blood flow, we developed an impedance boundary condition using morphologically derived structured trees. Our model was validated by computing blood flow through a model of large arteries extending from the thoracic aorta to the profunda arteries. The computed flow was compared against measured flow in the infrarenal (IR) aorta at rest and during exercise. Phase contrast-magnetic resonance imaging (PC-MRI) data was collected from 11 healthy volunteers at rest and during steady exercise. For each subject, an allometrically-scaled geometry of the large vessels was created. This geometry extends from the thoracic aorta to the femoral arteries and includes the celiac, superior mesenteric, renal, inferior mesenteric, internal iliac and profunda arteries. During rest, flow was simulated using measured supraceliac (SC) flow at the inlet and a uniform set of impedance boundary conditions at the 11 outlets. To simulate exercise, boundary conditions were modified. Inflow data collected during steady exercise was specified at the inlet and the outlet boundaries were adjusted as follows. The geometry of the structured trees used to compute impedance was scaled to simulate the effective change in the cross-sectional area of resistance vessels and capillaries due to exercise. The resulting computed flow through the IR aorta was compared to measured flow. This method produces good results with a mean difference between paired data to be 1.1 ± 7 cm³ s^? 1 at rest and 4.0 ± 15 cm³ s^? 1 at exercise. While future work will improve on these results, this method provides groundwork with which to predict the flow distributions in a network due to physiologic regulation.     

The effects of total ankle arthroplasty on postural stability and loading symmetry in quiet stance

Ankle osteoarthritis is a debilitating condition affecting about 1% of the population with approximately 50,000 new instances annually. One treatment is total ankle arthroplasty (TAA), however, its effects on balance are not well understood. This study analyzed balance over a two-year period following TAA. 408 subjects (177 left, 231 right ankles) diagnosed with end-stage ankle osteoarthritis performed quiet standing trials while center of pressure (COP) data were collected. Data were compared across three time points (pre-op, 1-year, and 2-years post-op) and between surgical and non-surgical limbs using a linear mixed model with significance set at P = 0.05. COP excursions in the feet-together condition were not significantly different between limbs after 2 years in anteroposterior or mediolateral directions (P = 0.06, 0.08) after being significantly different between limbs in the anteroposterior (P = 0.014) and mediolateral direction (P < 0.001) pre-op. The vertical ground reaction force significantly decreased across time in the non-surgical limb, while reciprocally increasing in the surgical limb (P < 0.001). After 2 years, no significant difference in vertical ground reaction force between limbs existed (P = 0.20). Limb asymmetry indices decreased at each time point in both conditions (all P < 0.001) and were not significantly different from zero after 2 years in the feet-together condition (P = 0.290). In conclusion, surgical limb balance improved compared to pre-op, resulting in increased symmetry between limbs after 2 years. Vertical ground reaction forces on both limbs converge and limb asymmetry indices approach zero two years post-op. Differences in the COP excursion-loading symmetry relationship between limbs could be useful for identifying instability in other pathologies.     

A strain-based finite element model for calcification progression in aortic valves

Calcific aortic valve disease (CAVD) is a serious disease affecting the aging population. A complex interaction between biochemicals, cells, and mechanical cues affects CAVD initiation and progression. In this study, motivated by the progression of calcification in regions of high strain, we developed a finite element method (FEM) based spatial calcification progression model. Several cardiac cycles of transient structural FEM simulations were simulated. After each simulation cycle, calcium deposition was placed in regions of high circumferential strain. Our results show the radial expansion of calcification as spokes starting from the attachment region, agreeing very well with the reported clinical data.     

Development and validation of an atlas-based finite element brain model

Traumatic brain injury is a leading cause of disability and injury-related death. To enhance our ability to prevent such injuries, brain response can be studied using validated finite element (FE) models. In the current study, a high-resolution, anatomically accurate FE model was developed from the International Consortium for Brain Mapping brain atlas. Due to wide variation in published brain material parameters, optimal brain properties were identified using a technique called Latin hypercube sampling, which optimized material properties against three experimental cadaver tests to achieve ideal biomechanics. Additionally, falx pretension and thickness were varied in a lateral impact variation. The atlas-based brain model (ABM) was subjected to the boundary conditions from three high-rate experimental cadaver tests with different material parameter combinations. Local displacements, determined experimentally using neutral density targets, were compared to displacements predicted by the ABM at the same locations. Error between the observed and predicted displacements was quantified using CORrelation and Analysis (CORA), an objective signal rating method that evaluates the correlation of two curves. An average CORA score was computed for each variation and maximized to identify the optimal combination of parameters. The strongest relationships between CORA and material parameters were observed for the shear parameters. Using properties obtained through the described multiobjective optimization, the ABM was validated in three impact configurations and shows good agreement with experimental data. The final model developed in this study consists of optimized brain material properties and was validated in three cadaver impacts against local brain displacement data.     

Development of structural and material clavicle response corridors under axial compression and three point bending loading for clavicle finite element model validation

Clavicle injuries were frequently observed in automotive side and frontal crashes. Finite element (FE) models have been developed to understand the injury mechanism, although no clavicle loading response corridors yet exist in the literature to ensure the model response biofidelity. Moreover, the typically developed structural level (e.g., force–deflection) response corridors were shown to be insufficient for verifying the injury prediction capacity of FE model, which usually is based on strain related injury criteria. Therefore, the purpose of this study is to develop both the structural (force vs deflection) and material level (strain vs force) clavicle response corridors for validating FE models for injury risk modeling. 20 Clavicles were loaded to failure under loading conditions representative of side and frontal crashes respectively, half of which in axial compression, and the other half in three point bending. Both structural and material response corridors were developed for each loading condition. FE model that can accurately predict structural response and strain level provides a more useful tool in injury risk modeling and prediction. The corridor development method in this study could also be extended to develop corridors for other components of the human body.     

Construction of a risk model through the fusion of experimental data and finite element modeling: Application to car crash-induced TBI

This article introduces a new approach for the construction of a risk model for the prediction of Traumatic Brain Injury (TBI) as a result of a car crash. The probability of TBI is assessed through the fusion of an experiment-based logistic regression risk model and a finite element (FE) simulation-based risk model. The proposed approach uses a multilevel framework which includes FE simulations of vehicle crashes with dummy and FE simulations of the human brain. The loading conditions derived from the crash simulations are transferred to the brain model thus allowing the calculation of injury metrics such as the Cumulative Strain Damage Measure (CSDM). The framework is used to propagate uncertainties and obtain probabilities of TBI based on the CSDM injury metric. The risk model from FE simulations is constructed from a support vector machine classifier, adaptive sampling, and Monte-Carlo simulations. An approach to compute the total probability of TBI, which combines the FE-based risk assessment as well as the risk prediction from the experiment-based logistic regression model is proposed. In contrast to previous published work, the proposed methodology includes the uncertainty of explicit parameters such as impact conditions (e.g., velocity, impact angle), and material properties of the brain model. This risk model can provide, for instance, the probability of TBI for a given assumed crash impact velocity.     

Finite element analysis of retroacetabular osteolytic defects following total hip replacement

Periprosthetic osteolysis in the retroacetabular region with cancellous bone loss is a recognized phenomenon in the long-term follow-up of total hip replacement. The effects on load transfer in the presence of defects are less well known. A finite element model incorporating a retroacetabular defect behind a cementless component was validated against a 4th generation sawbone pelvis. Computational predictions of surface strain and von Mises stresses were closely correlated to experimental findings. The presence of a cancellous defect increased von Mises stress in the cortical bone of the medial wall of the pelvis. At a load of 600 N this was under the predicted failure stress for cortical bone. Increases in the cup size relative to the acetabulum caused increased stress in the cortical bone of the lateral wall of the pelvis, adjacent to the acetabulum. We are confident that our modeling approach can be applied to patient specific defects to predict pelvis stress with large loads and a range of activities.     

Parametric convergence sensitivity and validation of a finite element model of the human lumbar spine

The primary objective of this study was to generate a finite element model of the human lumbar spine (L1–L5), verify mesh convergence for each tissue constituent and perform an extensive validation using both kinematic/kinetic and stress/strain data. Mesh refinement was accomplished via convergence of strain energy density (SED) predictions for each spinal tissue. The converged model was validated based on range of motion, intradiscal pressure, facet force transmission, anterolateral cortical bone strain and anterior longitudinal ligament deformation predictions. Changes in mesh resolution had the biggest impact on SED predictions under axial rotation loading. Nonlinearity of the moment-rotation curves was accurately simulated and the model predictions on the aforementioned parameters were in good agreement with experimental data. The validated and converged model will be utilised to study the effects of degeneration on the lumbar spine biomechanics, as well as to investigate the mechanical underpinning of the contemporary treatment strategies.     

Generalized finite element solution to one-dimensional flux problems

A finite element numerical solution to the general one-dimensional flow equation is derived in a form that provides a convenient and general means to simulate a wide variety of one-dimensional flow techniques of interest to biological scientists, e.g., ultracentrifugation, electrophoresis, chromatography, etc. Diverse physical models defined in terms of column geometry, solute interactions, and the dependence of transport parameters on column position, time, or concentrations of one or more solutes, can be accommodated. A particularly useful aspect of the formulation is that a wide variety of boundary conditions can be simply applied to the end result, without rederivation of the solution for each new case. The numerical solution is expressed as matrix equations that are sufficiently general so that incorporation of particular models can be effected by substitution of appropriate quantities into the final result.