Coupled motions under compressive load in intact and ACL-deficient knees: a cadaveric study

Knowledge of the coupled motions, which develop under compressive loading of the knee, is useful to determine which degrees of freedom should be included in the study of tibiofemoral contact and also to understand the role of the anterior cruciate ligament (ACL) in coupled motions. The objectives of this study were to measure the coupled motions of the intact knee and ACL-deficient knee under compression and to compare the coupled motions of the ACL-deficient knee with those of the intact knee. Ten intact cadaveric knees were tested by applying a 1600 N compressive load and measuring coupled internal-external and varus-valgus rotations and anterior-posterior and medial-lateral translations at 0 deg, 15 deg, and 30 deg of flexion. Compressive loads were applied along the functional axis of axial rotation, which coincides approximately with the mechanical axis of the tibia. The ACL was excised and the knees were tested again. In the intact knee, the peak coupled motions were 3.8 deg internal rotation at 0 deg flexion changing to -4.9 deg external rotation at 30 deg of flexion, 1.4 deg of varus rotation at 0 deg flexion changing to -1.9 deg valgus rotation at 30 deg of flexion, 1.4 mm of medial translation at 0 deg flexion increasing to 2.3 mm at 30 deg of flexion, and 5.3 mm of anterior translation at 0 deg flexion increasing to 10.2 mm at 30 deg of flexion. All changes in the peak coupled motions from 0 deg to 30 deg flexion were statistically significant (p<0.05). In ACL-deficient knees, there was a strong trend (marginally not significant, p=0.07) toward greater anterior translation (12.7 mm) than that in intact knees (8.0 mm), whereas coupled motions in the other degrees of freedom were comparable. Because the coupled motions in all four degrees of freedom in the intact knee and ACL-deficient knee are sufficiently large to substantially affect the tibiofemoral contact area, all degrees of freedom should be included when either developing mathematical models or designing mechanical testing equipment for study of tibiofemoral contact. The increase in coupled anterior translation in ACL-deficient knees indicates the important role played by the ACL in constraining anterior translation during compressive loading.     

Dynamic finite element knee simulation for evaluation of knee replacement mechanics

In vitro pre-clinical testing of total knee replacement (TKR) devices is a necessary step in the evaluation of new implant designs. Whole joint knee simulators, like the Kansas knee simulator (KKS), provide a controlled and repeatable loading environment for comparative evaluation of component designs or surgical alignment under dynamic conditions. Experimental testing, however, is time and cost prohibitive for design-phase evaluation of tens or hundreds of design variations. Experimentally-verified computational models provide an efficient platform for analysis of multiple components, sizes, and alignment conditions. The purpose of the current study was to develop and verify a computational model of a dynamic, whole joint knee simulator. Experimental internal-external and valgus-varus laxity tests, followed by dynamic deep knee bend and gait simulations in the KKS were performed on three cadaveric specimens. Specimen-specific finite element (FE) models of posterior-stabilized TKR were created from magnetic resonance images and CAD geometry. The laxity data was used to optimize mechanical properties of tibiofemoral soft-tissue structures on a specimen-specific basis. Each specimen was subsequently analyzed in a computational model of the experimental KKS, simulating both dynamic activities. The computational model represented all joints and actuators in the experimental setup, including a proportional-integral-derivative (PID) controller to drive quadriceps actuation. The computational model was verified against six degree-of-freedom patellofemoral (PF) and tibiofemoral (TF) kinematics and actuator loading during both deep knee bend and gait activities, with good agreement in trends and magnitudes between model predictions and experimental kinematics; differences were less than 1.8 mm and 2.2° for PF and TF translations and rotations. The whole joint FE simulator described in this study can be applied to investigate a wide range of clinical and research questions.     

A novel non-invasive method for measuring knee joint laxity in four dof: In vitro proof-of-concept and validation

Knee joint laxity or instability is a common problem that may have detrimental consequences for patients. Unfortunately, assessment of knee joint laxity is limited by current methodologies resulting in suboptimal diagnostics and treatment. This paper presents a novel method for accurately measuring non-invasive knee joint laxity in four degrees-of-freedom (DOF). An arthrometer, combining a parallel manipulator and a six-axis force/moment sensor, was developed to be used in combination with a low-dose biplanar x-ray system and 3D image data to reconstruct tibiofemoral position and orientation of laxity measurements.As proof-of-concept, four cadaveric knees were tested in the device. Each cadaveric knee was mounted in the device at approximately 30° of flexion and twelve monoplanar anteroposterior, mediolateral and internal/external load cases were applied. Additionally, four biplanar load cases were applied, consisting of different combinations of anteroposterior and internal/external loads. The arthrometer was limited to four DOF to address the specific measurements. For validation purposes, the pose reconstructions of tibia and femur were compared with pose reconstructions of bone pin marker frames mounted on each bone.The measurements from the arthrometer in terms of translation and rotations displayed comparable values to what have previously been presented in the literature. Furthermore, the measurements revealed coupled motions in multiple planes, demonstrating the importance of multi DOF laxity measurements. The validation displayed an average mean difference for translations of 0.08 mm and an average limit of agreement between −1.64 mm and 1.80 mm. The average mean difference for rotations was 0.10° and the limit of agreement was between −0.85° and 1.05°.The presented method eliminates several limitations present in current methods and may prove a valuable tool for assessing knee joint laxity.     

A new device to measure the structural properties of the femur-anterior cruciate ligament-tibia complex

Previous studies of biomechanical properties of femur-anterior cruciate ligament-tibia complex (FATC) utilized a wide variety of testing methodologies, particularly with respect to ligament orientation relative to loading direction. A new device was designed and built to test the anterior-posterior displacement of the intact porcine knee at 30 and 90 deg of flexion, as well as the tensile properties of the FATC at any loading direction and flexion angle. Tensile tests were performed with the knees at 30 and 90 deg of flexion with the loading direction along either the axis of the tibia (tibial axis) or the axis of the anterior cruciate ligament (ligament axis). The results showed that the stiffness, ultimate load and energy absorbed were all significantly increased when the FATC was tested along the ligament axis. This study demonstrates the importance of alignment in the evaluation of the biomechanical characteristics of the femur-ACL-tibia complex.     

Ligament balancing in TKA: evaluation of a force-sensing device and the influence of patellar eversion and ligament release

Ligament balancing in total knee arthroplasty may have an important influence on joint stability and prosthesis lifetime. In order to provide quantitative information and assistance during ligament balancing, a device that intraoperatively measures knee joint forces and moments was developed. Its performance and surgical advantages were evaluated on six cadaver specimens mounted on a knee joint loading apparatus allowing unconstrained knee motion as well as compression and varus-valgus loading. Four different experiments were performed on each specimen. (1) Knee joints were axially loaded. Comparison between applied and measured compressive forces demonstrated the accuracy and reliability of in situ measurements (1.8N). (2) Assessment of knee stability based on condyle contact forces or varus-valgus moments were compared to the current surgical method (difference of varus-valgus loads causing condyle lift-off). The force-based approach was equivalent to the surgical method while the moment-based, which is considered optimal, showed a tendency of lateral imbalance. (3) To estimate the importance of keeping the patella in its anatomical position during imbalance assessment, the effect of patellar eversion on the mediolateral distribution of tibiofemoral contact forces was measured. One fourth of the contact force induced by the patellar load was shifted to the lateral compartment. (4) The effect of minor and major medial collateral ligament releases was biomechanically quantified. On average, the medial contact force was reduced by 20% and 46%, respectively. Large variation among specimens reflected the difficulty of ligament release and the need for intraoperative force monitoring. This series of experiments thus demonstrated the device's potential to improve ligament balancing and survivorship of total knee arthroplasty.     

Experimental validation of a tibiofemoral model for analyzing joint force distribution

A computational model of the tibiofemoral joint utilizing the discrete element analysis method has been developed and validated with human cadaveric knees. The computational method can predict load distributions to within a root mean square error (RMSE) of 3.6%. The model incorporates subject-specific joint geometry and the health of the subjects’ articular cartilage to determine the cartilage stiffness. It also includes the collateral and cruciate ligaments and utilizes stiffness values derived from literature for these elements. Comparisons of the total load, peak load, and peak load location for axial, varus, and valgus loading conditions confirmed that there was less than 4% RMSE between the analytical and experimental results. The model presented in this paper can generate results with minimal computational time and it can be used as a non-invasive method for characterizing and monitoring subject-specific knee loading patterns.     

In vivo tibiofemoral cartilage deformation during the stance phase of gait

The knowledge of articular cartilage contact biomechanics in the knee joint is important for understanding the joint function and cartilage pathology. However, the in vivo tibiofemoral articular cartilage contact biomechanics during gait remains unknown. The objective of this study was to determine the in vivo tibiofemoral cartilage contact biomechanics during the stance phase of treadmill gait. Eight healthy knees were magnetic resonance (MR) scanned and imaged with a dual fluoroscopic system during gait on a treadmill. The tibia, femur and associated cartilage were constructed from the MR images and combined with the dual fluoroscopic images to determine in vivo cartilage contact deformation during the stance phase of gait. Throughout the stance phase of gait, the magnitude of peak compartmental contact deformation ranged between 7% and 23% of the resting cartilage thickness and occurred at regions with thicker cartilage. Its excursions in the anteroposterior direction were greater in the medial tibiofemoral compartment as compared to those in the lateral compartment. The contact areas throughout the stance phase were greater in the medial compartment than in the lateral compartment. The information on in vivo tibiofemoral cartilage contact biomechanics during gait could be used to provide physiological boundaries for in vitro testing of cartilage. Also, the data on location and magnitude of deformation among non-diseased knees during gait could identify where loading and later injury might occur in diseased knees.     

On the estimate of the two dominant axes of the knee using an instrumented spatial linkage

This article presents the validation of a technique to assess the appropriateness of a 2 degree-of-freedom model for the human knee, and, in which case, the dominant axes of flexion/extension and internal/external longitudinal rotation are estimated. The technique relies on the use of an instrumented spatial linkage for the accurate detection of passive knee kinematics, and it is based on the assumption that points on the longitudinal rotation axis describe nearly circular and planar trajectories, whereas the flexion/extension axis is perpendicular to those trajectories through their centers of rotation. By manually enforcing a tibia rotation while bending the knee in flexion, a standard optimization algorithm is used to estimate the approximate axis of longitudinal rotation, and the axis of flexion is estimated consequently. The proposed technique is validated through simulated data and experimentally applied on a 2 degree-of-freedom mechanical joint. A procedure is proposed to verify the fixed axes assumption for the knee model. The suggested methodology could be possibly valuable in understanding knee kinematics, and in particular for the design and implant of customized hinged external fixators, which have shown to be effective in knee dislocation treatment and rehabilitation.     

Simulating dynamic activities using a five-axis knee simulator

This work describes the design and capabilities of the Purdue Knee Simulator: Mark II and a sagittal-plane model of the machine. This five-axis simulator was designed and constructed to simulate dynamic loading activities on either cadaveric knee specimens or total knee prostheses mounted on fixtures. The purpose of the machine was to provide a consistent, realistic loading of the knee joint, allowing the kinematics and specific loading of the structures of the knee to be determined based on condition, articular geometry, and simulated activity. The sagittal-plane model of the knee simulator was developed both to predict the loading at the knee from arbitrary inputs and to generate the necessary inputs required to duplicate specified joint loading. Measured tibio-femoral compressive force and quadriceps tension were shown to be in good agreement with the predicted loads from the model. A controlled moment about the ankle-flexion axis was also shown to change the loading on the quadriceps.     

A new method to measure post-traumatic joint contractures in the rabbit knee

A new device and method to measure rabbit knee joint angles are described. The method was used to measure rabbit knee joint angles in normal specimens and in knee joints with obvious contractures. The custom-designed and manufactured gripping device has two clamps. The femoral clamp sits on a pinion gear that is driven by a rack attached to a materials testing system. A 100 N load cell in series with the rack gives force feedback. The tibial clamp is attached to a rotatory potentiometer. The system allows the knee joint multiple degrees-of-freedom (DOF). There are two independent DOF (compression-distraction and internal-external rotation) and two coupled motions (medial-lateral translation coupled with varus-valgus rotation; anterior-posterior translation coupled with flexion-extension rotation). Knee joint extension-flexion motion is measured, which is a combination of the materials testing system displacement (converted to degrees of motion) and the potentiometer values (calibrated to degrees). Internal frictional forces were determined to be at maximum 2% of measured loading. Two separate experiments were performed to evaluate rabbit knees. First, normal right and left pairs of knees from four New Zealand White (NZW) rabbits were subjected to cyclic loading. An extension torque of 0.2 Nm was applied to each knee. The average change in knee joint extension from the first to the fifth cycle was 1.9 deg +/- 1.5 deg (mean +/- sd) with a total of 49 tests of these eight knees. The maximum extension of the four left knees (tested 23 times) was 14.6 deg +/- 7.1 deg, and of the four right knees (tested 26 times) was 12.0 deg +/- 10.9 deg. There was no significant difference in the maximum extension between normal left and right knees. In the second experiment, nine skeletally mature NZW rabbits had stable fractures of the femoral condyles of the right knee that were immobilized for five, six or 10 weeks. The left knee served as an unoperated control. Loss of knee joint extension (flexion contracture) was demonstrated for the experimental knees using the new methodology where the maximum extension was 35 deg +/- 9 deg, compared to the unoperated knee maximum extension of 11 deg +/- 7 deg, 10 or 12 weeks after the immobilization was discontinued. The custom gripping device coupled to a materials testing machine will serve as a measurement test for future studies characterizing a rabbit knee model of post-traumatic joint contractures.     

The effect of hamstring muscle compensation for anterior laxity in the ACL-deficient knee during gait

The hamstring muscles have been recognized as an important element in compensating for the loss of stability in the ACL-deficient knee, but it is still not clear whether the hamstring muscle force can completely compensate for the loss of ACL, and the consequences of increased hamstring muscle force. A two-dimensional anatomical knee model in the sagittal plane was developed to examine the effect of various levels of hamstring muscle activation on restraining anterior tibial translation in the ACL-deficient knee during level walking. The model included the tibiofemoral and patellofemoral joints, four major ligaments, the medial capsule, and five muscle units surrounding the knee. Simulations were conducted to determine anterior tibial translation and internal joint loading at a single selected position when the knee was under a peak external flexion moment during early stance phase of gait. Incremental hamstring muscle forces were applied to the modeled normal and the ACL-deficient knees. Results of simulations showed that the ACL injury increased the anterior tibial translation by 11.8mm, while 56% of the maximal hamstring muscle force could reduce the anterior translation of the tibia to a normal level during the stance phase of gait. The consequences of increased hamstring muscle force included increased quadriceps muscle force and joint contact force.     

The effect of the frontal plane tibiofemoral angle and varus knee moment on the contact stress and strain at the knee cartilage

Subject-specific models were developed and finite element analysis was performed to observe the effect of the frontal plane tibiofemoral angle on the normal stress, Tresca shear stress and normal strain at the surface of the knee cartilage. Finite element models were created for three subjects with different tibiofemoral angle and physiological loading conditions were defined from motion analysis and muscle force mathematical models to simulate static single-leg stance. The results showed that the greatest magnitude of the normal stress, Tresca shear stress and normal strain at the medial compartment was for the varus aligned individual. Considering the lateral knee compartment, the individual with valgus alignment had the largest stress and strain at the cartilage. The present investigation is the first known attempt to analyze the effects of tibiofemoral alignment during single-leg support on the contact variables of the cartilage at the knee joint. The method could be potentially used to help identify individuals most susceptible to osteoarthritis and to prescribe preventive measures.     

Dynamic joint and muscle forces during knee isokinetic exercise

Isokinetic exercise has been commonly used in knee rehabilitation, conditioning and research in the past two decades. Although many investigators have used various experimental and theoretical approaches to study the muscle and joint force involved in isokinetic knee extension and flexion exercises, only a few of these studies have actually distinguished between the tibiofemoral joint forces and muscle forces. Therefore, the objective of this study was to specify, via an eletromyography(EMG)-driven muscle force model of the knee, the magnitude of the tibiofemoral joint and muscle forces acting during isokinetic knee extension and flexion exercises. Fifteen subjects ranging from 21 to 36 years of age volunteered to participate in this study. A Kin Com exercise machine (Chattecx Corporation, Chattanooga, TN, U.S.A.) was used as the loading device. An EMG-driven muscle force model was used to predict muscle forces, and a biomechanical model was used to analyze two knee joint constraint forces; compression and shear force. The methods used in this study were shown to be valid and reliable (r > 0.84 andp < 0.05). The effects on the tibiofemoral joint force during knee isokinetic exercises were compared with several functional activities that were investigated by earlier researchers. The muscle forces generated during knee isokinetic exercise were also obtained. Based on the findings obtained in this study, several therapeutic justifications for knee rehabilitation are proposed.     

'Equivalent geometry' of the knee and the prediction of tensions along the cruciates: an experimental study

In this paper we introduce the concept of the functional (or equivalent) geometry of the knee, which is an attempt to reduce the natural knee with its complex geometry, frictional resistance and deformable cartilage into a two-dimensional joint comprising rigid femur and tibia in frictionless contact. An apparatus and method are described to measure the slope of the tangent to the surfaces of the 'equivalent' bones at their 'point' of contact. An antero posterior force of +/-300-500 N and axial compressive load of twice body weight were applied on cadaveric knee joints. The corresponding displacement of the tibia in the saggital plane was measured firstly with both cruciates intact and then when each was severed in turn. From the data obtained both the slope of the tangent mentioned above and the tensions developed along the cruciates under the influence of the forces applied were calculated. The results showed that the functional geometry of the knee in the saggital plane can be represented by a convex femur and a concave tibia. The tensions along the cruciates calculated on the basis of the experimental measurements were nearly always lower than the antero posterior force applied, and although this corroborated the trend demonstrated in a previous theoretical analysis, they were lower still. The reason for this may be the deformation of the cartilage under load, thus modifying the geometry of contact resulting in a more concave tibia of the 'equivalent' knee joint, than that of the rigid model used in the theoretical analysis.     

Quantified in vitro tibiofemoral contact during bodyweight back squats

Squats are a common lower extremity task used in strength and conditioning, balance training, and rehabilitation. It is important to understand how slight alterations in lower extremity kinematics during a squat affect the internal joint loading of the knee. This study directly quantified tibiofemoral contact throughout the in vitro simulation of a bodyweight back squat performed two ways: a heel squat (knees in line with toes) and a toe squat (knees anterior to the toes) at peak knee flexion. Three cadaveric right lower extremities were instrumented and positioned into the University of Texas Joint Load Simulator. Kinematics, kinetics, and predicted muscle forces from a 20-year-old athletic male performing the two back squats were used as inputs for the in vitro simulations. The quantified tibiofemoral contact area, peak pressure, net force, and center of pressure location were significantly different between squat types (p > 0.05). Net contact area on the tibial plateau at peak knee flexion was significantly larger in the heel versus toe squat (599 ± 80 mm² vs. 469 ± 125 mm²; p < 0.05). Peak lateral pressure was significantly higher in the heel versus toe squat (2.73 ± 0.54 MPa vs. 0.87 ± 0.56 MPa; p < 0.05). Results suggest the heel squat generates an even load distribution, which is less likely to affect joint degeneration. Future in vitro simulations should quantify the effects lower extremity kinematics, kinetics, and individual muscle forces have on tibiofemoral contact parameters during common athletic tasks.     

The coupled motion of the femur and patella during in vivo weightbearing knee flexion

The movement of the knee joint consists of a coupled motion between the tibiofemoral and patellofemoral articulations. This study measured the six degrees-of-freedom kinematics of the tibia, femur, and patella using dual-orthogonal fluoroscopy and magnetic resonance imaging. Ten normal knees from ten living subjects were investigated during weightbearing flexion from full extension to maximum flexion. The femoral and the patellar motions were measured relative to the tibia. The femur externally rotated by 12.9 deg and the patella tilted laterally by 16.3 deg during the full range of knee flexion. Knee flexion was strongly correlated with patellar flexion (R(2)=0.91), posterior femoral translation was strongly correlated to the posterior patellar translation (R(2)=0.87), and internal-external rotation of the femur was correlated to patellar tilt (R(2)=0.73) and medial-lateral patellar translation (R(2)=0.63). These data quantitatively indicate a kinematic coupling between the tibia, femur, and patella, and provide base line information on normal knee joint kinematics throughout the full range of weightbearing flexion. The data also suggest that the kinematic coupling of tibia, femur, and patella should be considered when investigating patellar pathologies and when developing surgical techniques to treat knee joint diseases.     

A method of calculating physiologically relevant joint reaction forces during forward dynamic simulations of movement from an existing knee model

In the commonly used SIMM software, which includes a complete musculoskeletal model of the lower limbs, the reaction forces at the knee are computed. These reaction forces represent the bone-on-bone contact forces and the soft tissue forces (e.g. ligaments) other than muscles acting at the joint. In the knee model integrated into this software, a patellotibial joint rather than a patellofemoral joint is defined, and a force acting along the direction of the patellar ligament is not included. Although this knee model results in valid kinematics and muscle moment arms, the reaction forces at the knee calculated do not represent physiologic knee joint reaction forces. Hence our objectives were to develop a method of calculating physiologic knee joint reaction forces using the knee model incorporated into the SIMM software and to demonstrate the differences in the forces returned by SIMM and the physiologic forces in an example. Our method converts the anatomically fictional patellotibial joint into a patellofemoral joint and computes the force in an inextensible patellar ligament. In our example, the rectus femoris was fully excited isometrically, with the knee and hip flexed to 90 degrees . The resulting SIMM tibiofemoral joint reaction force was primarily shear, because the quadriceps force was applied to the tibia via the fictional patellotibial joint. In contrast the physiologic tibiofemoral joint reaction force was primarily compression, because the quadriceps force was applied through the patellar ligament. This result illustrates that the physiologic knee joint reaction forces are profoundly different than the forces returned by SIMM. However physiologic knee joint reaction forces can be computed with postprocessing of SIMM results.     

Subject-specific knee joint geometry improves predictions of medial tibiofemoral contact forces

Estimating tibiofemoral joint contact forces is important for understanding the initiation and progression of knee osteoarthritis. However, tibiofemoral contact force predictions are influenced by many factors including muscle forces and anatomical representations of the knee joint. This study aimed to investigate the influence of subject-specific geometry and knee joint kinematics on the prediction of tibiofemoral contact forces using a calibrated EMG-driven neuromusculoskeletal model of the knee. One participant fitted with an instrumented total knee replacement walked at a self-selected speed while medial and lateral tibiofemoral contact forces, ground reaction forces, whole-body kinematics, and lower-limb muscle activity were simultaneously measured. The combination of generic and subject-specific knee joint geometry and kinematics resulted in four different OpenSim models used to estimate muscle–tendon lengths and moment arms. The subject-specific geometric model was created from CT scans and the subject-specific knee joint kinematics representing the translation of the tibia relative to the femur was obtained from fluoroscopy. The EMG-driven model was calibrated using one walking trial, but with three different cost functions that tracked the knee flexion/extension moments with and without constraint over the estimated joint contact forces. The calibrated models then predicted the medial and lateral tibiofemoral contact forces for five other different walking trials. The use of subject-specific models with minimization of the peak tibiofemoral contact forces improved the accuracy of medial contact forces by 47% and lateral contact forces by 7%, respectively compared with the use of generic musculoskeletal model.     

Repeated application of incremental landing impact loads to intact knee joints induces anterior cruciate ligament failure and tibiofemoral cartilage deformation and damage: A preliminary cadaveric investigation

Anterior cruciate ligament (ACL) injury is a major problem worldwide and prevails during high-impact activities. It is not well-understood how the extent and distribution of cartilage damage will arise from repetitive landing impact loads that can lead to ACL failure. This study seeks to investigate the sole effect of repetitive incremental landing impact loads on the induction of ACL failure, and extent and distribution of tibiofemoral cartilage damage in cadaveric knees. Five cadaveric knees were mounted onto a material testing system at 70° flexion to simulate landing posture. A motion-capture system was used to track rotational and translational motions of the tibia and femur, respectively. Each specimen was compressed at a single 10 Hz haversine to simulate landing impact. The compression trial was successively repeated with increasing actuator displacement till a significant compressive force drop was observed. All specimens underwent ACL failure, which was confirmed via magnetic resonance scans and dissection. Volume analysis, thickness measurement and histological techniques were employed to assess cartilage lesion status. For each specimen, the highest peak compressive force (1.9–7.8 kN) was at the final trial in which ACL failure occurred; corresponding posterior femoral displacement (7.6–18.0 mm) and internal tibial rotation (0.6°–4.7°) were observed. Significant compressive force drop (79.8–90.9%) was noted upon ACL failure. Considerable cartilage deformation and damage were found in exterior, posterior and interior femoral regions with substantial volume reduction in lateral compartments. Repeated application of incremental landing impact loads can induce both ACL failure and cartilage damage, which may accelerate the risk of developing osteoarthritis.