Finite element lumbar spine facet contact parameter predictions are affected by the cartilage thickness distribution and initial joint gap size

Current finite element modeling techniques utilize geometrically inaccurate cartilage distribution representations in the lumbar spine. We hypothesize that this shortcoming severely limits the predictive fidelity of these simulations. Specifically, it is unclear how these anatomically inaccurate cartilage representations alter range of motion and facet contact predictions. In the current study, cadaveric vertebrae were serially sectioned, and images were taken of each slice in order to identify the osteochondral interface and the articulating surface. A series of custom-written algorithms were utilized in order to quantify each facet joint's three-dimensional cartilage distribution using a previously developed methodology. These vertebrae-dependent thickness cartilage distributions were implemented on an L1 through L5 lumbar spine finite element model. Moments were applied in three principal planes of motion, and range of motion and facet contact predictions from the variable thickness and constant thickness distribution models were determined. Initial facet gap thickness dimensions were also parameterized. The data indicate that the mean and maximum cartilage thickness increased inferiorly from L1 to L5, with an overall mean thickness value of 0.57 mm. Cartilage distribution and initial facet joint gap thickness had little influence on the lumbar range of motion in any direction, whereas the mean contact pressure, total contact force, and total contact area predictions were altered considerably. The data indicate that range of motion predictions alone are insufficient to establish model validation intended to predict mechanical contact parameters. These data also emphasize the need for the careful consideration of the initial facet joint gap thickness with respect to the spinal condition being studied.     

Specimen-specific predictions of contact stress under physiological loading in the human hip: validation and sensitivity studies

Hip osteoarthritis may be initiated and advanced by abnormal cartilage contact mechanics, and finite element (FE) modeling provides an approach with the potential to allow the study of this process. Previous FE models of the human hip have been limited by single specimen validation and the use of quasi-linear or linear elastic constitutive models of articular cartilage. The effects of the latter assumptions on model predictions are unknown, partially because data for the instantaneous behavior of healthy human hip cartilage are unavailable. The aims of this study were to develop and validate a series of specimen-specific FE models, to characterize the regional instantaneous response of healthy human hip cartilage in compression, and to assess the effects of material nonlinearity, inhomogeneity and specimen-specific material coefficients on FE predictions of cartilage contact stress and contact area. Five cadaveric specimens underwent experimental loading, cartilage material characterization and specimen-specific FE modeling. Cartilage in the FE models was represented by average neo-Hookean, average Veronda Westmann and specimen- and region-specific Veronda Westmann hyperelastic constitutive models. Experimental measurements and FE predictions compared well for all three cartilage representations, which was reflected in average RMS errors in contact stress of less than 25 %. The instantaneous material behavior of healthy human hip cartilage varied spatially, with stiffer acetabular cartilage than femoral cartilage and stiffer cartilage in lateral regions than in medial regions. The Veronda Westmann constitutive model with average material coefficients accurately predicted peak contact stress, average contact stress, contact area and contact patterns. The use of subject- and region-specific material coefficients did not increase the accuracy of FE model predictions. The neo-Hookean constitutive model underpredicted peak contact stress in areas of high stress. The results of this study support the use of average cartilage material coefficients in predictions of cartilage contact stress and contact area in the normal hip. The regional characterization of cartilage material behavior provides the necessary inputs for future computational studies, to investigate other mechanical parameters that may be correlated with OA and cartilage damage in the human hip. In the future, the results of this study can be applied to subject-specific models to better understand how abnormal hip contact stress and contact area contribute to OA.     

Finite element simulations of a focal knee resurfacing implant applied to localized cartilage defects in a sheep model

Articular resurfacing metal implants have recently been tested in animal models to treat full thickness localized articular cartilage defects, showing promising results. However, the mechanical behavior of cartilage surrounding the metal implant has not been studied yet as it is technically challenging to measure in vivo contact areas, pressures, stresses and deformations from the metal implant. Therefore, we implemented a detailed numerical finite element model by approximating one of the condyles of the sheep tibiofemoral joint and created a defect of specific size to accommodate the implant. Using this model, the mechanical behavior of the surrounding of metal implant was studied. The model showed that the metal implant plays a significant role in the force transmission. Two types of profiles were investigated for metal implant. An implant with a double-curved profile, i.e., a profile fully congruent with the articular surfaces in the knee, gives lower contact pressures and stresses at the rim of the defect than the implant with unicurved spherical profile. The implant should be placed at a certain distance into the cartilage to avoid damage to opposing biological surface. Too deep positions, however, lead to high shear stresses in the cartilage edges around the implant. Mechanical sealing was achieved with a wedge shape of the implant, also useful for biochemical sealing of cartilage edges at the defect.     

Effects of idealized joint geometry on finite element predictions of cartilage contact stresses in the hip

Computational models may have the ability to quantify the relationship between hip morphology, cartilage mechanics and osteoarthritis. Most models have assumed the hip joint to be a perfect ball and socket joint and have neglected deformation at the bone-cartilage interface. The objective of this study was to analyze finite element (FE) models of hip cartilage mechanics with varying degrees of simplified geometry and a model with a rigid bone material assumption to elucidate the effects on predictions of cartilage stress. A previously validated subject-specific FE model of a cadaveric hip joint was used as the basis for the models. Geometry for the bone-cartilage interface was either: (1) subject-specific (i.e. irregular), (2) spherical, or (3) a rotational conchoid. Cartilage was assigned either a varying (irregular) or constant thickness (smoothed). Loading conditions simulated walking, stair-climbing and descending stairs. FE predictions of contact stress for the simplified models were compared with predictions from the subject-specific model. Both spheres and conchoids provided a good approximation of native hip joint geometry (average fitting error ～0.5 mm). However, models with spherical/conchoid bone geometry and smoothed articulating cartilage surfaces grossly underestimated peak and average contact pressures (50% and 25% lower, respectively) and overestimated contact area when compared to the subject-specific FE model. Models incorporating subject-specific bone geometry with smoothed articulating cartilage also underestimated pressures and predicted evenly distributed patterns of contact. The model with rigid bones predicted much higher pressures than the subject-specific model with deformable bones. The results demonstrate that simplifications to the geometry of the bone-cartilage interface, cartilage surface and bone material properties can have a dramatic effect on the predicted magnitude and distribution of cartilage contact pressures in the hip joint.     

Load sharing between solid and fluid phases in articular cartilage: II--Comparison of experimental results and u-p finite element predictions.

Experimental measurements in conjunction with theoretical predictions were used to determine the extent of load supported by the fluid phase of cartilage at the articular surface. The u-p finite element model was used to simulate the loading of six separate porcine knee joints and to predict surface deformations of the cartilage layer on the lateral femoral condyle. Representative geometry for the condyle, contact pressures, and intrinsic material properties of the cartilage layer were supplied from experimental measures (see Part I). The u-p finite element predictions for surface deformations of the cartilage layer were obtained for several load partitioning states between the solid and fluid phases of cartilage at the articular surface. These were then compared to actual surface deformations obtained experimentally. It appeared from the comparison that approximately 75 percent of the applied load was borne by the fluid phase at the articular surface under this loading regime. This was qualitatively in agreement with the hypothesis that an applied load to articular joints is partitioned at the surface to the two phases according to the surface area ratios of the solid and fluid phases. It appeared that the solid phase was shielded from the total applied stress on the articular surface by the fluid and could be a reason for the excellent durability of the tissue under the demanding conditions in a diarthrodial joint.     

Development of a statistical shape model of the patellofemoral joint for investigating relationships between shape and function

Patellofemoral (PF)-related pathologies, including joint laxity, patellar maltracking, cartilage degradation and anterior knee pain, affect nearly 25% of the population. Researchers have investigated the influence of articular geometry on kinematics and contact mechanics in order to gain insight into the etiology of these conditions. The purpose of the current study was to create a three-dimensional statistical shape model of the PF joint and to characterize relationships between PF shape and function (kinematics and contact mechanics). A statistical shape model of the patellar and femoral articular surfaces and their relative alignment was developed from magnetic resonance images. Using 15 shape parameters, the model characterized 97% of the variation in the training set. The first three shape modes primarily described variation in size, patella alta-baja and depth of the sulcus groove. A previously verified finite element model was used to predict kinematics and contact mechanics for each subject. Combining the shape and joint mechanics data, a statistical shape-function model was developed that established quantitative relations of how changes in the shape of the PF joint influence mechanics. The predictive capability of the shape-function model was evaluated by comparing statistical model and finite element predictions, resulting in kinematic root mean square errors of less than 3° and 2.5 mm. The key results of the study are dually in the implementation of a novel approach linking statistical shape and finite element models and the relationships elucidated between PF articular geometry and mechanics.     

Anisotropic hydraulic permeability in compressed articular cartilage

The extent to which articular cartilage hydraulic permeability is anisotropic is largely unknown, despite its importance for understanding mechanisms of joint lubrication, load bearing, transport phenomena, and mechanotransduction. We developed and applied new techniques for the direct measurement of hydraulic permeability within statically compressed adult bovine cartilage explant disks, dissected such that disk axes were perpendicular to the articular surface. Applied pressure gradients were kept small to minimize flow-induced matrix compaction, and fluid outflows were measured by observation of a meniscus in a glass capillary under a microscope. Explant disk geometry under radially unconfined axial compression was measured by direct microscopic observation. Pressure, flow, and geometry data were input to a finite element model where hydraulic permeabilities in the disk axial and radial directions were determined. At less than 10% static compression, near free-swelling conditions, hydraulic permeability was nearly isotropic, with values corresponding to those of previous studies. With increasing static compression, hydraulic permeability decreased, but the radially directed permeability decreased more dramatically than the axially directed permeability such that strong anisotropy (a 10-fold difference between axial and radial directions) in the hydraulic permeability tensor was evident for static compression of 20-40%. Results correspond well with predictions of a previous microstructurally-based model for effects of tissue mechanical deformations on glycosaminoglycan architecture and cartilage hydraulic permeability. Findings inform understanding of structure-function relationships in cartilage matrix, and suggest several biomechanical roles for compression-induced anisotropic hydraulic permeability in articular cartilage.     

Variability of a three-dimensional finite element model constructed using magnetic resonance images of a knee for joint contact stress analysis

Magnetic resonance (MR) imaging has been widely used to evaluate the thickness and volume of articular cartilage both in vivo and in vitro. While morphological information on the cartilage can be obtained using MR images, image processing for extracting geometric boundaries of the cartilage may introduce variations in the thickness of the cartilage. To evaluate the variability of using MR images to construct finite element (FE) knee cartilage models, five investigators independently digitized the same set of MR images of a human knee. The topology of cartilage thickness was determined using a minimal distance algorithm. Less than 8 percent variation in cartilage thickness was observed from the digitized data. The effect of changes in cartilage thickness on contact stress analysis was then investigated using five FE models of the knee. One FE model (average FE model) was constructed using the mean values of the digitized contours of the cartilage, and the other four were constructed by varying the thickness of the average FE model by +/- 5 percent and +/- 10 percent, respectively. The results demonstrated that under axial tibial compressive loading (up to 1,400 N), variations of cartilage thickness caused by digitization of MR images may result in a difference of approximately 10 percent in peak contact stresses (surface pressure, von Mises stress, and hydrostatic pressure) in the cartilage. A reduction of cartilage thickness caused increases of contact stresses, while an increase of cartilage thickness reduced contact stresses. Furthermore, the effect of variation of material properties of the cartilage on contact stress analysis was investigated. The peak contact stress increased almost linearly with the Young's modulus of the cartilage. The peak von Mises stress was dramatically reduced when the Poisson,s ratio was increased from 0.05 to 0.49 under an axial compressive load of 1,400 N, while peak hydrostatic pressure was dramatically increased. Peak surface pressure was also increased with the Poisson's ratio, but with a lower magnitude compared to von Mises stress and hydrostatic pressure. In conclusion, the imaging process may cause 10 percent variations in peak contact stress, and the predicted stress distribution is sensitive to the accuracy of the material properties of the cartilage model, especially to the variation of Poisson's ratio.     

Modelling cartilage mechanobiology

The growth, maintenance and ossification of cartilage are fundamental to skeletal development and are regulated throughout life by the mechanical cues that are imposed by physical activities. Finite element computer analyses have been used to study the role of local tissue mechanics on endochondral ossification patterns, skeletal morphology and articular cartilage thickness distributions. Using single-phase continuum material representations of cartilage, the results have indicated that local intermittent hydrostatic pressure promotes cartilage maintenance. Cyclic tensile strains (or shear), however, promote cartilage growth and ossification. Because single-phase material models cannot capture fluid exudation in articular cartilage, poroelastic (or biphasic) solid/fluid models are often implemented to study joint mechanics. In the middle and deep layers of articular cartilage where poroelastic analyses predict little fluid exudation, the cartilage phenotype is maintained by cyclic fluid pressure (consistent with the single-phase theory). In superficial articular layers the chondrocytes are exposed to tangential tensile strain in addition to the high fluid pressure. Furthermore, there is fluid exudation and matrix consolidation, leading to cell 'flattening'. As a result, the superficial layer assumes an altered, more fibrous phenotype. These computer model predictions of cartilage mechanobiology are consistent with results of in vitro cell and tissue and molecular biology experiments.     

Verification of predicted specimen-specific natural and implanted patellofemoral kinematics during simulated deep knee bend

Verified computational models represent an efficient method for studying the relationship between articular geometry, soft-tissue constraint, and patellofemoral (PF) mechanics. The current study was performed to evaluate an explicit finite element (FE) modeling approach for predicting PF kinematics in the natural and implanted knee. Experimental three-dimensional kinematic data were collected on four healthy cadaver specimens in their natural state and after total knee replacement in the Kansas knee simulator during a simulated deep knee bend activity. Specimen-specific FE models were created from medical images and CAD implant geometry, and included soft-tissue structures representing medial–lateral PF ligaments and the quadriceps tendon. Measured quadriceps loads and prescribed tibiofemoral kinematics were used to predict dynamic kinematics of an isolated PF joint between 10° and 110° femoral flexion. Model sensitivity analyses were performed to determine the effect of rigid or deformable patellar representations and perturbed PF ligament mechanical properties (pre-tension and stiffness) on model predictions and computational efficiency.Predicted PF kinematics from the deformable analyses showed average root mean square (RMS) differences for the natural and implanted states of less than 3.1° and 1.7 mm for all rotations and translations. Kinematic predictions with rigid bodies increased average RMS values slightly to 3.7° and 1.9 mm with a five-fold decrease in computational time. Two-fold increases and decreases in PF ligament initial strain and linear stiffness were found to most adversely affect kinematic predictions for flexion, internal–external tilt and inferior–superior translation in both natural and implanted states. The verified models could be used to further investigate the effects of component alignment or soft-tissue variability on natural and implant PF mechanics.     

Effects of inserting a pressensor film into articular joints on the actual contact mechanics.

Fuji film has been widely used in studies aimed at obtaining the contact mechanics of articular joints. Once sealed for practical use in biological joints, Fuji Pressensor film has a total effective thickness of 0.30 mm, which is comparable to the cartilage thickness in the joints of many small animals. The average effective elastic modulus of Fuji film is approximately 100 MPa in compression, which is larger by a factor of 100-300 compared to that of normal articular cartilage. Therefore, inserting a Pressensor film into an articular joint will change the contact mechanics of the joint. The measurement precision of the Pressensor film has been determined systematically; however, the changes in contact mechanics associated with inserting the film into joints have not been investigated. This study was aimed at quantifying the changes in the contact mechanics associated with inserting sealed Fuji Pressensor film into joints. Spherical and cylindrical articular joint contact mechanics with and without Pressensor film and for varying degrees of surface congruency were analyzed and compared by using finite element models. The Pressensor film was taken as linearly elastic and the cartilage was assumed to be biphasic, composed of a linear elastic solid phase and an inviscid fluid phase. The present analyses showed that measurements of the joint contact pressures with Fuji Pressensor film will change the maximum true contact pressures by 10-26 percent depending on the loading, geometry of the joints, and the mechanical properties of cartilage. Considering this effect plus the measurement precision of the film (approximately 10 percent), the measured joint contact pressures in a joint may contain errors as large as 14-28 percent.     

Changes in articular cartilage mechanics with meniscectomy: A novel image-based modeling approach and comparison to patterns of OA

Meniscectomy is a significant risk factor for osteoarthritis, involving altered cell synthesis, central fibrillation, and peripheral osteophyte formation. Though changes in articular cartilage contact pressure are known, changes in tissue-level mechanical parameters within articular cartilage are not well understood. Recent imaging research has revealed the effects of meniscectomy on the time-dependent deformation of physiologically loaded articular cartilage. To determine tissue-level cartilage mechanics that underlie observed deformation, a novel finite element modeling approach using imaging data and a contacting indenter boundary condition was developed. The indenter method reproduces observed articular surface deformation and avoids assumptions about tangential stretching. Comparison of results from an indenter model with a traditional femur-tibia model verified the method, giving errors in displacement, solid and fluid stress, and strain below 1% (RMS) and 7% (max.) of the absolute maximum of the parameters of interest. Indenter finite element models using real joint image data showed increased fluid pressure, fluid exudation, loss of fluid load support, and increased tensile strains centrally on the tibial condyle after meniscectomy-patterns corresponding to clinical observations of cartilage matrix damage and fibrillation. Peripherally there was decreased consolidation, which corresponds to reduced contact and fluid pressure in this analysis. Clinically, these areas have exhibited advance of the subchondral growth front, biological destruction of the cartilage matrix, cartilage thinning, and eventual replacement of the cartilage via endochondral ossification. Characterizing the changes in cartilage mechanics with meniscectomy and correspondence with observed tissue-level effects may help elucidate the etiology of joint-level degradation seen in osteoarthritis.     

C4–C5 segment finite element model development,validation, and load-sharing investigation

Detailed cervical spine models are necessary to better understand cervical spine response to loading, improve our understanding of injury mechanisms, and specifically for predicting occupant response and injury in auto crash scenarios. The focus of this study was to develop a C4–C5 finite element model with accurate representations of each tissue within the segment. This model incorporates more than double the number of elements of existing models, required for accurate prediction of response. The most advanced material data available were then incorporated using appropriate nonlinear constitutive models to provide accurate predictions of response at physiological levels of loading. This tissue-scale segment model was validated against a wide variety of experimental data including different modes of loading (axial rotation, flexion, extension, lateral bending, and translation), and different load levels. In general, the predicted response of the model was within the single standard deviation response corridors for both low and high load levels. Importantly, this model demonstrates that appropriate refinement of the finite element mesh, representation at the tissue level, and sufficiently detailed material properties and constitutive models provide excellent response predictions without calibration of the model to experimental data. Load sharing between the disc, ligaments, and facet joints was investigated for various modes of loading, and the dominant load-bearing structure was found to correlate with typical anatomical injury sites for these modes of loading. The C4–C5 model forms the basis for the development of a full cervical spine model. Future studies will focus on tissue-level injury prediction and dynamic response.     

Influence of cervical disc degeneration after posterior surgical techniques in combined flexion-extension--a nonlinear analytical study

Laminectomy and facetectomy are surgical techniques used for decompression of the cervical spinal stenosis. Recent in vitro and finite element studies have shown significant cervical spinal instability after performing these surgical techniques. However, the influence of degenerated cervical disk on the biomechanical responses of the cervical spine after these surgical techniques remains unknown. Therefore, a three-dimensional nonlinear finite element model of the human cervical spine (C2-C7) was created. Two types of disk degeneration grades were simulated. For each grade of disk degeneration, the intact as well as the two surgically altered models simulating C5 laminectomy with or without C5-C6 total facetectomies were exercised under flexion and extension. Intersegmental rotational motions, internal disk annulus, cancellous and cortical bone stresses were obtained and compared to the normal intact model. Results showed that the cervical rotational motion decreases with progressive disk degeneration. Decreases in the rotational motion due to disk degeneration were accompanied by higher cancellous and cortical bone stress. The surgically altered model showed significant increases in the rotational motions after laminectomies and facetectomies when compared to the intact model. However, the percentage increases in the rotational motions after various surgical techniques were reduced with progressive disk degeneration.     

Predictive modelling of cervical disc implant wear

This study presents a chain of simulations aimed at estimating the wear in a cervical disc implant and providing insight into the in vivo biomechanical performance of the implant. The simulation chain can start with determining a representative maximum range of motion (ROM) of a person's head. The ROM is used as motion input to a kinematic simulation of the cervical spine containing a disc implant. The cervical spine geometry is obtained from computed tomography (CT) scans and converted to STL format using reverse engineering software. The time histories of the loads imposed by the adjacent vertebrae on the implant, as well as the vertebral relative rotations can be extracted from the kinematic simulation. Alternatively, force and motion profiles prescribed by wear test protocols (e.g. ISO 18192-1 and ASTM F2423-05) can be used. The force and motion profiles are applied as boundary conditions to a non-linear finite element model (FEM) of the implant to determine the time-varying contact stress and slip velocity distributions at the interface between the two halves of the implant. The stresses and slip velocities are used in a linear wear model to estimate the wear rate distribution at the FEM's nodal points where contact occurs. Reverse engineering software is used to triangulate the contact surface so that the total wear volume can be calculated. The simulation chain's predicted wear rate shows good agreement with in vitro results in the literature. The simulation chain is thereby demonstrated to be suitable for comparative pre-experimental studies of spinal implant designs.     

Influence of aspect ratios on the creep behaviour of articular cartilage in indentation

Indentation tests are commonly used to determine the mechanical behaviour of articular cartilage with varying properties, thickness, and geometry. This investigation evaluated the effect of changing geometric parameters on the properties determined from creep indentation tests. Finite element analyses simulated the indentation behaviour of two models, an excised cylindrical specimen of cartilage with either normal and repair qualities and an osteochondral defect represented as a cylindrical region of repair cartilage integrated with a surrounding layer of normal tissue. For each model, the ratios of indenter radius to cartilage height (a/h=0.5,1.5) and cartilage radius to indenter radius (r/a=2,5) were varied. The vertical displacement of the cartilage under the indenter obtained through finite element analysis was fitted to a numerical algorithm to determine the aggregate modulus, permeability, and Poisson's ratio. Indentation behaviours of cartilage specimens for either model with a/h=1.5 were not affected by r/a for values of 2 and 5. Aggregate modulus was not greatly affected by the geometric changes studied. Permeability was affected by changes in the ratio of specimen to indenter radii for a/h=0.5. These findings suggest that experimental configurations of excised cylindrical specimens, also representing osteochondral defects with no or unknown degree of integration, where the cartilage layer has a/h=0.5 should not have r/a values on the order of 2 for confidence in the mechanical properties determined. Indentation of osteochondral defects where the repair cartilage is fully integrated to the surrounding cartilage can be performed with confidence for all cases tested.     

A finite element model of the human knee joint for the study of tibio-femoral contact

As a step towards developing a finite element model of the knee that can be used to study how the variables associated with a meniscal replacement affect tibio-femoral contact, the goals of this study were 1) to develop a geometrically accurate three-dimensional solid model of the knee joint with special attention given to the menisci and articular cartilage, 2) to determine to what extent bony deformations affect contact behavior, and 3) to determine whether constraining rotations other than flexion/extension affects the contact behavior of the joint during compressive loading. The model included both the cortical and trabecular bone of the femur and tibia, articular cartilage of the femoral condyles and tibial plateau, both the medial and lateral menisci with their horn attachments, the transverse ligament, the anterior cruciate ligament, and the medial collateral ligament. The solid models for the menisci and articular cartilage were created from surface scans provided by a noncontacting, laser-based, three-dimensional coordinate digitizing system with an root mean squared error (RMSE) of less than 8 microns. Solid models of both the tibia and femur were created from CT images, except for the most proximal surface of the tibia and most distal surface of the femur which were created with the three-dimensional coordinate digitizing system. The constitutive relation of the menisci treated the tissue as transversely isotropic and linearly elastic. Under the application of an 800 N compressive load at 0 degrees of flexion, six contact variables in each compartment (ie., medial and lateral) were computed including maximum pressure, mean pressure, contact area, total contact force, and coordinates of the center of pressure. Convergence of the finite element solution was studied using three mesh sizes ranging from an average element size of 5 mm by 5 mm to 1 mm by 1 mm. The solution was considered converged for an average element size of 2 mm by 2 mm. Using this mesh size, finite element solutions for rigid versus deformable bones indicated that none of the contact variables changed by more than 2% when the femur and tibia were treated as rigid. However, differences in contact variables as large as 19% occurred when rotations other than flexion/extension were constrained. The largest difference was in the maximum pressure. Among the principal conclusions of the study are that accurate finite element solutions of tibio-femoral contact behavior can be obtained by treating the bones as rigid. However, unrealistic constraints on rotations other than flexion/extension can result in relatively large errors in contact variables.     

Biomechanical properties of knee articular cartilage

Structure and properties of knee articular cartilage are adapted to stresses exposed on it during physiological activities. In this study, we describe site- and depth-dependence of the biomechanical properties of bovine knee articular cartilage. We also investigate the effects of tissue structure and composition on the biomechanical parameters as well as characterize experimentally and numerically the compression-tension nonlinearity of the cartilage matrix. In vitro mechano-optical measurements of articular cartilage in unconfined compression geometry are conducted to obtain material parameters, such as thickness, Young's and aggregate modulus or Poisson's ratio of the tissue. The experimental results revealed significant site- and depth-dependent variations in recorded parameters. After enzymatic modification of matrix collagen or proteoglycans our results show that collagen primarily controls the dynamic tissue response while proteoglycans affect more the static properties. Experimental measurements in compression and tension suggest a nonlinear compression-tension behavior of articular cartilage in the direction perpendicular to articular surface. Fibril reinforced poroelastic finite element model was used to capture the experimentally found compression-tension nonlinearity of articular cartilage.     

Quantitation of articular surface topography and cartilage thickness in knee joints using stereophotogrammetry

An analytical stereophotogrammetry (SPG) technique has been developed based upon some of the pioneering work of Selvik [Ph.D. thesis, University of Lund, Sweden (1974)] and Huiskes and coworkers [J. Biomechanics 18, 559-570 (1985)], and represents a fundamental step in the construction of biomechanical models of diarthrodial joints. Using this technique, the precise three-dimensional topography of the cartilage surfaces of various diarthrodial joints has been obtained. The system presented in this paper delivers an accuracy of 90 microns in the least favorable conditions with 95% coverage using the same calibration method as Huiskes et al. (1985). In addition, a method has been developed, using SPG, to quantitatively map the cartilage thickness over the entire articular surface of a joint with a precision of 134 microns (95% coverage). In the present study, our SPG system has been used to quantify the topography, including surface area, of the articular surfaces of the patella, distal femur, tibial plateau, and menisci of the human knee. Furthermore, examples of cartilage thickness maps and corresponding thickness data including coefficient of variation, minimum, maximum, and mean cartilage thickness are also provided for the cartilage surfaces of the knee. These maps illustrate significant variations over the joint surfaces which are important in the determination of the stresses and strains within the cartilage during diarthrodial joint function. In addition, these cartilage surface topographies and thickness data are essential for the development of anatomically accurate finite element models of diarthrodial joints.(ABSTRACT TRUNCATED AT 250 WORDS)     

Dynamic measurement of internal solid displacement in articular cartilage using ultrasound backscatter

Mechanics of articular cartilage can be represented using poroelastic theories where fluid and solid displacements are viscously coupled to create a time-dependent spatially heterogeneous behavior. In recent models of this tissue, finite element methods have been used to predict tissue deformation as a function of time for adult articular cartilage bearing a characteristic depth-dependent structure and composition. However, current experimental methods are limited in providing verification of these predictions. The current study presents an apparatus for imaging the radial displacement profile of cartilage in unconfined compression using an ultrasound technique called elastography. We acquired ultrasound A-scans across the lateral diameter of full-thickness cartilage disks containing a thin layer of underlying bone, during axial compression. Elastography was then applied to correlate temporally sequential A-scans to estimate the solid radial displacement profile in articular cartilage while it undergoes compression and stress-relaxation. Both time-dependent and depth-dependent solid radial displacement profiles were obtained with a precision better than 0.2 micro The results generally agree with predictions of poroelastic models, demonstrating lateral expansion with an effective Poisson's ratio just after completion of the compression phase of the mechanical tests reaching values from 0.18 to 0.4 (depending on compression speed), followed by contraction to lower values. A more restricted movement was observed at both the articular surface and near to the subchondral bone than at regions midway between these two locations.