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.     

A new transducer for facet force measurement in the lumbar spine: benchmark and in vitro test results.

A new transducer capable of direct measurement of time-dependent loads in human lumbar facet joints was developed and tested. The transducer was comprised of a force-sensitive resistor (FSR) in series with a pressure-sensitive film. A wide range of experiments revealed the performance attributes and limitations of the FSR. The output signal of the FSR is actually sensitive to both force and area of contact independently. Therefore, a pressure-sensitive film was used to quantify the contact area. At least two transformation equations were calculated for each FSR corresponding to known contact areas. Each equation was a linearization of the log of the FSR output vs the log of the applied ramp loads. Coefficients of determination (CD) were calculated for small (21 mm2) and large (32 mm2) contact areas, and were found to exceed 0.900 for all data. The average of nine cycles was nearly linear for some FSRs (CD of 0.999). FSR output signal and contact area were recorded in cadaveric lumbar facets under ramp load. The appropriate transformation equation was determined by a linear interpolation between benchmark equations based on the contact area measured in vitro. Facet force measurements compared well with those of other researchers. The transducer was found to be quite easy to use.     

Material and functional properties of articular cartilage and patellofemoral contact mechanics in an experimental model of osteoarthritis

The purposes of this study were to determine the in situ functional and material properties of articular cartilage in an experimental model of joint injury, and to quantify the corresponding in situ joint contact mechanics. Experiments were performed in the anterior cruciate ligament (ACL) transected knee of the cat and the corresponding, intact contralateral knee, 16 weeks following intervention. Cartilage thickness, stiffness, effective Young’s modulus, and permeability were measured and derived from six locations of the knee. The total contact area and peak pressures in the patellofemoral joint were obtained in situ using Fuji Pressensor film, and comparisons between experimental and contralateral joint were made for corresponding loading conditions. Total joint contact area and peak pressure were increased and decreased significantly (=0.01), respectively, in the experimental compared to the contralateral joint. Articular cartilage thickness and stiffness were increased and decreased significantly (=0.01), respectively, in the experimental compared to the contralateral joint in the four femoral and patellar test locations. Articular cartilage material properties (effective Young’s modulus and permeability) were the same in the ACL-transected and intact joints. These results demonstrate for the first time the effect of changes in articular cartilage properties on the load transmission across a joint. They further demonstrate a substantial change in the joint contact mechanics within 16 weeks of ACL transection. The results were corroborated by theoretical analysis of the contact mechanics in the intact and ACL-transected knee using biphasic contact analysis and direct input of cartilage properties and joint surface geometry from the experimental animals. We conclude that the joint contact mechanics in the ACL-transected cat change within 16 weeks of experimental intervention.     

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.     

Load sharing between solid and fluid phases in articular cartilage: I--Experimental determination of in situ mechanical conditions in a porcine knee.

The in situ mechanical conditions of cartilage in the articulated knee were quantified during joint loading. Six porcine knees were subjected to a 445 N compressive load while cartilage deformations and contact pressures were measured. From roentgenograms, cartilage thickness before and during loading allowed the calculation of tissue deformation on the lateral femoral condyle at different times during the loading process. Contact pressures on the articular surface were measured with miniature fiber-optic pressure transducers. Results showed that the medial side of the lateral femoral condyle had higher contact pressures, as well as deformations. To begin to correlate the pressures and resulting deformations, the intrinsic material properties of the cartilage on the lateral condyle were obtained from indentation tests. Data from four normal control specimens indicated that the aggregate modulus of the medial side was significantly higher than in other areas of the condyle. These experimental measures of the in situ mechanical conditions of articular cartilage can be combined with theoretical modeling to obtain valuable information about the relative contributions of the solid and fluid phases to supporting the applied load on the cartilage surface (see Part II).     

An improved method for measuring tibiofemoral contact areas in total knee arthroplasty: a comparison of K-scan sensor and Fuji film.

A computerised, real time, thin-film pressure transducer method is used to measure tibiofemoral contact area in total knee arthroplasty (TKA) devices that is easier and more reliable and reproducible as compared to the Fuji pressure-sensitive film technique. Many authors have suggested that contact areas and pressures within TKA devices can be a predictor of wear and failure of the ultra-high molecular weight polyethylene (UHMWPE) tibial insert. In this study, two contact area measurement techniques (Fuji pressure-sensitive Film and K-scan sensor system) were compared using a custom TKA testing jig designed for freedom of movement so that in any loading configuration the component found and seated in its own "home" position. The K-scan system was used to measure contact areas of one TKA design at several angles from 0 to 110 degrees flexion with loads equating to 4, 4.5, and 5 times body weight. For comparison, four ranges of Fuji film were used to measure areas at the same flexion angles but at 5 times body weight only. Contact areas measured with the Fuji films were 11-36% (p < 0.05) lower than those measured by the K-scan sensor.     

Contact area and pressure distribution in the feline patellofemoral joint under physiologically meaningful loading conditions.

The purpose of this study was to determine contact area and mean and peak pressures in the healthy feline patellofemoral joint over the complete range of possible applied force. Furthermore, we wanted to improve upon the repeatability of previous measurements while maximizing the physiological relevance of the results obtained. The patellae and femora were secured in a loading frame approximating an in situ loading configuration. Low- and medium-grade Fuji film was used to assess patellofemoral contact area and pressure distribution, respectively. Constant force was applied to the patellofemoral joints for 2s (short duration trials) or 5min (long duration trials). For the short duration trials, contact area was shown to increase logarithmically with the force applied. In contrast, mean and peak pressures increased linearly with force. Furthermore, the rate of increase of peak pressure with force was approximately three times greater than that of mean pressure. For the long duration trials, contact area increased up to 33% compared to the short duration trials. This effect could no longer be detected with our approach after an unloading period of 5-10s. Increasing contact area is one mechanism that the feline patellofemoral joint may use to regulate the pressures experienced by the cartilage as the force applied to the joint increases. The attenuation of external forces inside a joint is achieved by the specific geometry of the articulating surfaces and the viscoelastic properties of the articular cartilage. It likely represents a natural protection of joints to high external load magnitudes.     

Comparison of nonlinear mechanical properties of bovine articular cartilage and meniscus

Nonlinear, linear and failure properties of articular cartilage and meniscus in opposing contact surfaces are poorly known in tension. Relationships between the tensile properties of articular cartilage and meniscus in contact with each other within knee joints are also not known. In the present study, rectangular samples were prepared from the superficial lateral femoral condyle cartilage and lateral meniscus of bovine knee joints. Tensile tests were carried out with a loading rate of 5 mm/min until the tissue rupture. Nonlinear properties of the toe region, linear properties in larger strains, and failure properties of both tissues were analysed. The strain-dependent tensile modulus of the toe region, Young's modulus of the linear region, ultimate tensile stress and toughness were on average 98.2, 8.3, 4.0 and 1.9 times greater (p<0.05) for meniscus than for articular cartilage. In contrast, the toe region strain, yield strain and failure strain were on average 9.4, 3.1 and 2.3 times greater (p<0.05) for cartilage than for meniscus. There was a significant negative correlation between the strain-dependent tensile moduli of meniscus and articular cartilage samples within the same joints (r=−0.690, p=0.014). In conclusion, the meniscus possesses higher nonlinear and linear elastic stiffness and energy absorption capability before rupture than contacting articular cartilage, while cartilage has longer nonlinear region and can withstand greater strains before failure. These findings point out different load carrying demands that both articular cartilage and meniscus have to fulfil during normal physiological loading activities of knee joints.     

An articular cartilage contact model based on real surface geometry

Abnormal, excessive stresses acting on articular joint surfaces are speculated to be one of the causes for joint degeneration. However, articular surface stresses have not been studied systematically, since it is technically difficult to measure in vivo contact areas and pressures in dynamic situations. Therefore, we implemented a numerical model of articular surface contact using accurate surface geometries. The model was developed for the cat patellofemoral joint. We demonstrated that small misalignments of the patella relative to the femur change the joint contact mechanics substantially for a given external load. These results suggest that misalignment might be studied as one of the factors causing articular cartilage disorder and joint degeneration.     

Adsorption, lubrication, and wear of lubricin on model surfaces: polymer brush-like behavior of a glycoprotein

Using a surface force apparatus, we have measured the normal and friction forces between layers of the human glycoprotein lubricin, the major boundary lubricant in articular joints, adsorbed from buffered saline solution on various hydrophilic and hydrophobic surfaces: i), negatively charged mica, ii), positively charged poly-lysine and aminothiol, and iii), hydrophobic alkanethiol monolayers. On all these surfaces lubricin forms dense adsorbed layers of thickness 60–100 nm. The normal force between two surfaces is always repulsive and resembles the steric entropic force measured between layers of end-grafted polymer brushes. This is the microscopic mechanism behind the antiadhesive properties showed by lubricin in clinical tests. For pressures up to ~6 atm, lubricin lubricates hydrophilic surfaces, in particular negatively charged mica (friction coefficient μ = 0.02–0.04), much better than hydrophobic surfaces (μ > 0.3). At higher pressures, the friction coefficient is higher (μ > 0.2) for all surfaces considered and the lubricin layers rearrange under shear. However, the glycoprotein still protects the underlying substrate from damage up to much higher pressures. These results support recent suggestions that boundary lubrication and wear protection in articular joints are due to the presence of a biological polyelectrolyte on the cartilage surfaces.     

Glenohumeral joint cartilage contact in the healthy adult during scapular plane elevation depression with external humeral rotation

The shoulder (glenohumeral) joint has the greatest range of motion of all human joints; as a result, it is particularly vulnerable to dislocation and injury. The ability to non-invasively quantify in-vivo articular cartilage contact patterns of joints has been and remains a difficult biomechanics problem. As a result, little is known about normal in-vivo glenohumeral joint contact patterns or the consequences that surgery has on altering them. In addition, the effect of quantifying glenohumeral joint contact patterns by means of proximity mapping, both with and without cartilage data, is unknown. Therefore, the objectives of this study are to (1) describe a technique for quantifying in-vivo glenohumeral joint contact patterns during dynamic shoulder motion, (2) quantify normal glenohumeral joint contact patterns in the young healthy adult during scapular plane elevation depression with external humeral rotation, and (3) compare glenohumeral joint contact patterns determined both with and without articular cartilage data. Our results show that the inclusion of articular cartilage data when quantifying in-vivo glenohumeral joint contact patterns has significant effects on the anterior–posterior contact centroid location, the superior–inferior contact centroid range of travel, and the total contact path length. As a result, our technique offers an advantage over glenohumeral joint contact pattern measurement techniques that neglect articular cartilage data. Likewise, this technique may be more sensitive than traditional 6-Degree-of-Freedom (6-DOF) joint kinematics for the assessment of overall glenohumeral joint health. Lastly, for the shoulder motion tested, we found that glenohumeral joint contact was located on the anterior–inferior glenoid surface.     

Two-dimensional strain fields on the cross-section of the bovine humeral head under contact loading

The objective of this study was to provide a detailed experimental assessment of the two-dimensional cartilage strain distribution on the cross-section of immature and mature bovine humeral heads subjected to contact loading at a relatively rapid physiological loading rate. Six immature and six mature humeral head specimens were loaded against glass and strains were measured at the end of a 5s loading ramp on the textured articular cross-section using digital image correlation analysis. The primary findings indicate that elevated tensile and compressive strains occur near the articular surface, around the center of the contact region. Few qualitative or quantitative differences were observed between mature and immature joints. Under an average contact stress of approximately 1.7 MPa, the peak compressive strains averaged -0.131+/-0.048, which was significantly less than the relative change in cartilage thickness, -0.104+/-0.032 (p<0.05). The peak tensile strains were significantly smaller in magnitude, at 0.0325+/-0.013. These experimental findings differ from a previous finite element analysis of articular contact, which predicted peak strains at the cartilage-bone interface even when accounting for the porous-hydrated nature of the tissue, its depth-dependent inhomogeneity, and the disparity between its tensile and compressive properties. These experimental results yield new insights into the local mechanical environment of the tissue and cells, and suggest that further refinements are needed in the modeling of contacting articular layers.     

An augmented lagrangian method for sliding contact of soft tissue

Despite the importance of sliding contact in diarthrodial joints, only a limited number of studies have addressed this type of problem, with the result that the mechanical behavior of articular cartilage in daily life remains poorly understood. In this paper, a finite element formulation is developed for the sliding contact of biphasic soft tissues. The augmented Lagrangian method is used to enforce the continuity of contact traction and fluid pressure across the contact interface. The resulting method is implemented in the commercial software COMSOL Multiphysics. The accuracy of the new implementation is verified using an example problem of sliding contact between a rigid, impermeable indenter and a cartilage layer for which analytical solutions have been obtained. The new implementation's capability to handle a complex loading regime is verified by modeling plowing tests of the temporomandibular joint (TMJ) disc.     

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).     

Biomechanical properties of human articular cartilage under compressive loads

The function of articular cartilage is to support and distribute loads and to provide lubrication in the diarthrodial joints. Cartilage function is described by proper mechanical and rheological properties, strain and depth-dependent, which are not completely assessed. Unconfined and confined compression are commonly used to evaluate the Young's modulus (E) and the aggregate modulus (H(A)), respectively. The Poisson's ratio (nu) can be calculated indirectly from the equilibrium compression data, or using the biphasic indentation technique; it has recently been optically evaluated by using video microscopy during unconfined compression. The transient response of articular cartilage during confined compression depends on its permeability k; a constant value of k can be easily identified by a simple analytical model of confined compression tests, whereas more complex models or direct measurements (permeation tests) are needed to study the permeability dependence on deformation. A poroelastic finite element model of articular cartilage was developed for this purpose. The elastic parameters (E,nu) of the model were evaluated performing unconfined compression creep tests on human articular cartilage disks, whereas k was identified from the confined test response. Our combined experimental and computational method can be used to identify the parameters that define the permeability dependence on deformation, as a function of depth from articular surface.     

Load-sharing in the lumbosacral spine in neutral standing & flexed postures – A combined finite element and inverse static study

Understanding load-sharing in the spine during in-vivo conditions is critical for better spinal implant design and testing. Previous studies of load-sharing that considered actual spinal geometry applied compressive follower load, with or without moment, to simulate muscle forces. Other studies used musculoskeletal models, which include muscle forces, but model the discs by simple beams or spherical joints and ignore the articular facet joints.This study investigated load-sharing in neutral standing and flexed postures using a detailed Finite Element (FE) model of the ligamentous lumbosacral spine, where muscle forces, gravity loads and intra-abdominal pressure, as predicted by a musculoskeletal model of the upper body, are input into the FE model. Flexion was simulated by applying vertebral rotations following spine rhythm measured in a previous in-vivo study, to the musculoskeletal model. The FE model predicted intradiscal pressure (IDP), strains in the annular fibers, contact forces in the facet joints, and forces in the ligaments. The disc forces and moments were determined using equilibrium equations, which considered the applied loads, including muscle forces and IDP, as well as forces in the ligaments and facet joints predicted by the FE model. Load-sharing was calculated as the portion of the total spinal load carried along the spine by each individual spinal structure. The results revealed that spinal loads which increased substantially from the upright to the flexed posture were mainly supported by the discs in the upright posture, whereas the ligaments’ contribution in resisting shear, compression, and moment was more significant in the flexed posture.     

Measurement of finger joint angles and maximum finger forces during cylinder grip activity

Finger joint angles and finger forces during maximal cylindrical grasping were measured using multi-camera photogrammetry and pressure-sensitive sheets, respectively. The experimental data were collected from four healthy subjects gripping cylinders of five different sizes. For joint angles, an image analysis system was used to digitize slides showing markers. During the calibration of the camera system, both the nonlinear least square and the direct linear transform methods were applied and compared, the former providing the fewer errors; it was used to determine joint angles. Data were collected from the pressure-sensitive grip films by using the same image analysis system as used in the collection of the joint angle data. The method of using pressure-sensitive sheets provided an estimation of the weighted centre of the phalangeal forces. Results indicate that finger flexion angles at the metacarpophalangeal and proximal interphalangeal joints gradually increase as cylinder diameter decreases, but that at the distal interphalangeal joint the angle remains constant throughout all cylinder sizes. It was also found that most of the radio-ulnar deviation and the axial rotation angles at the finger joints deviate from zero, but the deviations are small. For the force measurement, it was found that total finger force increases as cylinder size decreases, and the phalangeal force centres are not located at the mid-points of the phalanges. The data obtained in this experiment would be useful for muscle force predictions and for the design of handles.     

Validation of radiocarpal joint contact models based on images from a clinical MRI scanner

This study was undertaken to assess magnetic resonance imaging (MRI)-based radiocarpal surface contact models of functional loading in a clinical MRI scanner for future in vivo studies, by comparison with experimental measures from three cadaver forearm specimens. Experimental data were acquired using a Tekscan sensor during simulated light grasp. Magnetic resonance (MR) images were used to obtain model geometry and kinematics (image registration). Peak contact pressures (PPs) and average contact pressures (APs), contact forces and contact areas were determined in the radiolunate and radioscaphoid joints. Contact area was also measured directly from MR images acquired with load and compared with model data. Based on the validation criteria (within 25% of experimental data), out of the six articulations (three specimens with two articulations each), two met the criterion for AP (0%, 14%); one for peak pressure (20%); one for contact force (5%); four for contact area with respect to experiment (8%, 13%, 19% and 23%), and three contact areas met the criterion with respect to direct measurements (14%, 21% and 21%). Absolute differences between model and experimental PPs were reasonably low (within 2.5 MPa). Overall, the results indicate that MRI-based models generated from 3T clinical MR scanner appear sufficient to obtain clinically relevant data.     

Finger joint impedance during tapping on a computer keyswitch

We studied the dynamic behavior of finger joints during the contact period of tapping on a computer keyswitch, to characterize and parameterize joint function with a lumped-parameter impedance model. We tested the hypothesis that the metacarpophalangeal (MCP) and interphalangeal (IP) joints act similarly in terms of kinematics, torque, and energy production when tapping. Fifteen human subjects tapped with the index finger of the right hand on a computer keyswitch mounted on a two-axis force sensor, which measured forces in the vertical and sagittal planes. Miniature fiber-optic goniometers mounted across the dorsal side of each joint measured joint kinematics. Joint torques were calculated from endpoint forces and joint kinematics using an inverse dynamic algorithm. For each joint, a linear spring and damper model was fitted to joint torque, position, and velocity during the contact period of each tap (22 per subject on average). The spring-damper model could account for over 90% of the variance in torque when loading and unloading portions of the contact were separated, with model parameters comparable to those previously measured during isometric loading of the finger. The finger joints functioned differently, as illustrated by energy production during the contact period. During the loading phase of contact the MCP joint flexed and produced energy, whereas the proximal and distal IP joints extended and absorbed energy. These results suggest that the MCP joint does work on the interphalangeal joints as well as on the keyswitch.