Simulated increases in body weight led to increased displacement, von Mises stress, and contact pressure in finite element models of the extended and flexed knee. Contact shifted to locations of typical medial osteoarthritis lesions in the extended knee models. Obesity is commonly associated with increased risk of osteoarthritis (OA). The effects of increases in body weight and other loads on the stresses and strains within a joint can be calculated using finite element (FE) models. The specific effects for different individuals can be calculated using subject-specific FE models which take individual geometry and forces into account. Model results can then be used to propose mechanisms by which damage within the joint may initiate.Summary Statement
Introduction
Osteoarthritis (OA) is a degenerative, chronic disease of the articular cartilage that affects more than 150 million people [1]. In the knee, OA can begin as either isolated medial OA or isolated lateral OA. Previous research [2,3] shows medial OA and lateral OA have characteristic cartilage lesion locations and progression patterns as well as flexion angles associated with lesion development, indicating strong involvement of mechanical factors in disease initiation. Therefore, it is important to investigate these mechanical factors. Previous studies combined data sets (geometry, motion, load) from separate sources. The aim of the current work was to use a consistent multi-modal approach. A finite element (FE) model of a healthy knee in full extension was created using magnetic resonance imaging (MRI) and motion analysis data from the same subject (female, 24 yrs). MRI data was obtained using a 3T MRI scanner (Philips Medical Systems/Achieva). Surface geometries of the tibia, femur, and associated cartilage were then semi-automatically segmented and processed (Mimics 12.5; Geomagic Studio 11; SolidWorks 2009). Motion data was collected at 100 Hz (Vicon 612) during level walking and subsequently applied to a lower limb model (AnyBody Version 3.0) to calculate muscle forces. Both sets of data were then combined to create a subject-specific FE model (ANSYS 11.0) which was solved to determine relative contact areas, pressures, and deformations in the medial and lateral tibiofemoral compartments.Background
Method
Iterative finite element (FE) models are used to simulate bone remodelling that takes place due to the surgical insertion of an implant or to simulate fracture healing. In such simulations element material properties are calculated after each iteration of solving the model. New material properties are calculated based on the results derived by the model during the last iteration. Once the FE model has gone through a number of such iterations it is often necessary to assess the remodelling that has taken place. The method widely used to do this is to analyse element Young's modulus plots taken at particular sections through the model. Although this method gives relevant information which is often helpful when comparing different implants, the information is rather abstract and is difficult to compare with patient data which is commonly in the form of radiographs. The authors suggest a simple technique that can be used to generate synthetic radiograph images from FE models. These images allow relatively easy comparisons of FE derived information with patient radiographs. Another clear advantage of this technique is that clinicians (who are familiar with reading radiographs) are able to understand and interpret them readily. To demonstrate the technique a three dimensional (3D) model of the proximal tibia implanted with an Oxford Unicompartmental Knee replacement was created based on CT data obtained from a cadaveric tibia. The model's initial element material properties were calculated from the same CT data set using a relationship between radiographic density and Young's modulus. The model was subject to simplified loading conditions and solved over 365 iterations representing one year of in vivo remodelling. After each iteration the element material properties were recalculated based on previously published remodelling rules. Next, synthetic anteroposterior radiographs were generated by back calculating radiographic densities from material properties of the model after 365 iterations. A 3D rectangular grid of sampling points which encapsulated the model was defined. For each of the elements in the FE model radiographic densities were back calculated based on the same relationships used to calculate material properties from radiographic densities. The radiographic density of each element was assigned to all the sampling grid points within the element. The 3D array of radiographic densities was summed in the anteroposterior direction thereby creating a 2D array of radiographic densities. This 2D array was plotted giving an image analogous to anteroposterior patient radiographs. Similar to a patient radiograph denser material appeared lighter while less dense material appeared darker. The resulting synthetic radiographs were compared to patient radiographs and found to have similar patterns of dark and light regions. The synthetic radiographs were relatively easy to produce based on the FE model results, represented FE results in a manner easily comparable to patient radiographs, and represented FE results in a clinician friendly manner.
The intact femur geometry was derived from a CT dataset of a cadaveric femur and CT numbers were converted into a realistic distribution of material properties. The FE intact mesh was based on an experimentally validated mesh of a human femur. The femur was segmented into 22 neck sections. The loading condition was modelled to represent an instant at 10% of gait where all muscle forces were included. The femoral neck regions were compared between the models to evaluate the effect of notch sizes on stress distribution. Maximum tensile stresses were compared to the ultimate tensile stress (UTS) of cortical and cancellous bone.
Three-dimensional motion of the lower limbs was measured using gait analysis. Transverse plane kinematics, including hip rotation and foot progression angles were recorded.
In finite element (FE) analysis of long bones it is now common practice to calculate the material properties based on CT data. Although a unique material property is calculated for each element, assigning each element an individual material property results in excessively large models. To avoid this, it is usual to group the elements based on their material properties and to assign each group a single material property (Zannoni 1998). No study has analysed the effect the number of material properties used in a long bone FE model has on the accuracy of the results. The aim of this study was to evaluate the variation in the calculated mechanical environment as a function of the number of material properties used in an FE model. An FE mesh of a cadaveric human tibia containing 47,696 ten-node tetrahedron elements and 75,583 nodes was created using CT scans. Material properties were calculated for each element of the mesh based on previous work (Rho 1995, 1996). Eleven FE models were created by varying the number of groups (1, 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024) the elements were divided into. A single material property was assigned to each group. All models were subject to an axial point load of 300N applied on the medial condyle of the tibial plateau while the distal end was fixed. The variation in maximum and minimum principal strains and deflections, at 17 well distributed surface nodes and at 65 randomly distributed nodes within the bone were plotted against the number of element groups. The total strain energy was also plotted against the number of groups. The errors for strain, deflection, and total strain energy were calculated for each model assuming that the model using 1024 element groups was accurate. The parameter to converge with the least number of element groups was the total strain energy. At 512 element groups the error was less than 0.001% (0.7% for the two material model). The next to converge were the displacements. Using 512 materials the maximum error in displacement at the surface nodes was 0.001% (4.7% for the 2 material model), while for the internal nodes the maximum error was 0.53% (36.7% for the 2 material model). The least convergence occurred for principal strains. The maximum errors when 512 materials were used were 1.06% (57.7% for the 2 material model) and 3.02% (104.5% for the 2 material model) for the surface and the internal nodes respectively. This study demonstrates the relationship between the accuracy of calculated mechanical environment and the number of material properties assigned to the model. While this study will allow the analyst to make an informed decision on the number of material properties for modelling the human tibia it also helps examine the validity of previous studies which, usually due to limited resources, used fewer material properties.
Finite element (FE) analysis is widely used to calculate stresses and strains within human bone in order to improve implant designs. Although validated FE models of the human femur have been created (Lengsfeld et al., 1998), no equivalent yet exists for the tibia. The aim of this study was to create such an FE model, both with and without the tibial component of a knee replacement, and to validate it against experimental Results: A set of reference axes was marked on a cleaned, fresh frozen cadaveric human tibia. Seventeen triaxial stacked strain rosettes were attached along the bone, which was then subjected to nine axial loading conditions, two four-point bending loading conditions, and a torsional loading condition using a materials testing machine (MTS 858). Deflections and strain readings were recorded. Axial loading was repeated after implantation of a knee replacement (medial tibial component, Biomet Oxford Unicompartmental Phase 3). The intact tibia was CT scanned (GE HiSpeed CT/i) and the images used to create a 3D FE mesh. The CT data was also used to map 600 transversely isotropic material properties (Rho, 1996) to individual elements. All experiments were simulated on the FE model. Measured principal strains and displacements were compared to their corresponding FE values using regression analysis. Experimental results were repeatable (mean coefficients of variation for intact and implanted tibia, 5.3% and 3.9%). They correlated well with those of the FE analysis (R squared = 0.98, 0.97, 0.97, and 0.99 for axial (intact), axial (implanted), bending, torsional loading). For each of the load cases the intersects of the regression lines were small in comparison to the maximum measured strains (<
1.5%). While the model was more rigid than the bone under torsional loading (slope =0.92), the opposite was true for axial (slope = 1.14 (intact) 1.24 (implanted)) and bending (slope = 1.06) loads. This is probably due to a discrepancy in the material properties of the model.
More than 100,000 anterior cruciate ligament reconstructions are performed annually in the USA. The hamstrings and the patellar tendons are the most frequently used graft tissues. Up to ten percent of these grafts are deemed to have failed, generating considerable discussion in the literature regarding the ideal graft choice. A three-dimensional computational model, taking into account both material and geometrical non-linearities, would be useful in predicting the performance of different graft tissues and fixations. Unfortunately, the mechanical characteristics and parameters needed for such a model are complex and largely unknown. The aim of this study is to develop a method for measuring the geometrical properties needed as input for a three-dimensional tendon model. A laser-based, non-contact technique is used to generate a series of cross-sectional profiles along the length of the tendon. Unlike previously proposed methods, it is able to detect concavities and can be constructed using equipment commonly found in an engineering laboratory. A laser line generator (Stocker-Yale Lasiris SNF, Quebec, Canada) projects a horizontal line onto the sample. Images of the line are acquired with a digital video camera (Basler A631fc, Germany) as the tendon is rotated. These images are reassembled into 2-D slices using MatLab software. Multiple cross-sections can be combined to create three dimensional geometries. The new method was validated on objects of known shape (circular and hexagonal cylinders). The cross-sectional area measurement was found to be accurate to within 2.5%. The method was repeatable to within 1.7%. Six bovine flexor tendons have been analysed; concavities were evident in four of these. This method could be adapted to determine the surface geometries of other long and slender objects.
This study assesses the functional in vivo kinematics of Advanced Medial Pivot (AMP) TKR and compares it to kinematics of the normal knee.
Finite element (FE) analysis is widely used to calculate stresses and strains within human bone in order to improve implant designs. Although validated FE models of the human femur have been created (Lengsfeld et al., 1998), no equivalent yet exists for the tibia. The aim of this study was to create such an FE model, both with and without the tibial component of a knee replacement, and to validate it against experimental results. A set of reference axes was marked on a cleaned, fresh frozen cadaveric human tibia. Seventeen triaxial stacked strain rosettes were attached along the bone, which was then subjected to nine axial loading conditions, two four-point bending loading conditions, and a torsional loading condition using a materials testing machine (MTS 858). Deflections and strain readings were recorded. Axial loading was repeated after implantation of a knee replacement (medial tibial component, Biomet Oxford Unicompartmental Phase 3). The intact tibia was CT scanned (GE HiSpeed CT/i) and the images used to create a 3D FE mesh. The CT data was also used to map 600 transversely isotropic material properties (Rho, 1996) to individual elements. All experiments were simulated on the FE model. Measured principal strains and displacements were compared to their corresponding FE values using regression analysis. Experimental results were repeatable (mean coeffi-cients of variation for intact and implanted tibia, 5.3% and 3.9%). They correlated well with those of the FE analysis (R squared = 0.98, 0.97, 0.97, and 0.99 for axial (intact), axial (implanted), bending, torsional loading). For each of the load cases the intersects of the regression lines were small in comparison to the maximum measured strains (<
1.5%). While the model was more rigid than the bone under torsional loading (slope =0.92), the opposite was true for axial (slope = 1.14 (intact) 1.24 (implanted)) and bending (slope = 1.06) loads. This is probably due to a discrepancy in the material properties of the model.
Patellofemoral pain is a significant problem for patients with Total Knee Replacements (TKRs). It is hypothesized that pain is related to high patellofemoral forces (PFF). The aim of this study is to validate a model to estimate PFF after TKR, using a combination of non-invasive measurement and theoretical modeling. Experiments were performed on four cadaver knee specimens to compare the PFF and the quadriceps force (QF) estimated by a model, with those measured using force transducers. Each knee was tested in its initial state and after implantation of three Scorpio designs: Cruciate Retaining (CR), Posterior Stabilised (PS), and the Posterior Stabilised Mobile Bearing (PS+). Each knee was extended/flexed under a simulated quadriceps load with 3 kg hung from the distal tibia. Relative movement of the bones was measured using a Vicon 612 motion analysis system. A 6DOF force transducer was used to measure PFFs and a uni-axial transducer was used to measure QFs. A fluoroscope simultaneously captured images of the leg extension activity. Parameters measured from the images were used as inputs to the model. The measured and estimated PFF and QF matched closely between 20o and 80o of knee flexion for the TKRs. At higher flexion angles, the model overestimated the PFF by a maximum of 23N (7.6% max) for the PFF and by 31N for the QF (10.3% max). The estimated and measured Patellar Flexion Angles (PFA) were within 3.5o throughout the flexion range. The model accurately predicts sagittal plane patellar kinematics and kinetics, using only fluoroscopy and externally measured forces as inputs. However, the model has a limitation in assuming that the extending moment is only due to the quadriceps.
Award for the best student biomaterials paper (US$ 2,000); a proper certificate