A Three-Dimensional Musculoskeletal Model of the Human Knee Joint. Part 1: Theoretical Construction

A three-dimensional model of the knee is developed to study the interactions between the muscles, ligaments, and bones during activity. The geometry of the distal femur, proximal tibia, and patella is based on cadaver data reported for an average-size knee. The shapes of the femoral condyles are represented by high-order polynomials: the tibial plateaux and patellar facets are approximated as flat surfaces. The contacting surfaces of the femur and tibia are modeled as deformable, while those of the femur and patella are assumed to be rigid. Interpenetration of the femur and tibia is taken into account by modeling cartilage as a thin, linear, elastic layer, mounted on rigid bone. Twelve elastic elements describe the geometry and mechanical properties of the cruciate ligaments, the collateral ligaments, and the posterior capsule. The model is actuated by thirteen musculotendinous units, each unit modeled as a three-element muscle in series with tendon. The path of each muscle is approximated as a straight line, except where it contacts and wraps around bone and other muscles; changes in muscle paths are taken into account using data obtained from MRI. In the first part of this paper, the model is used to simulate passive knee flexion. Quantitative comparisons of the model results with experimental data reported in the literature indicate that the relative movements of the bones and the geometry of the ligaments and muscles in the model are similar to those evident in the real knee. In Part II, the model is used to describe knee-ligament function during anterior-posterior draw, axial rotation, and isometric knee-extension and knee-flexion exercises.     

A Three-Dimensional Musculoskeletal Model of the Human Knee Joint. Part 2: Analysis of Ligament Function

A three-dimensional model of the knee is used to study ligament function during anterior-posterior (a-p) draw, axial rotation, and isometric contractions of the extensor and flexor muscles. The geometry of the model bones is based on cadaver data. The contacting surfaces of the femur and tibia are modeled as deformable; those of the femur and patella are assumed to be rigid. Twelve elastic elements are used to describe the geometry and mechanical properties of the cruciate ligaments, the collateral ligaments, and the posterior capsule. The model is actuated by thirteen musculotendinous units, each unit represented as a three-element muscle in series with tendon. The calculations show that the forces applied during a-p draw are substantially different from those applied by the muscles during activity. Principles of knee-ligament function based on the results of in vitro experiments may therefore be overstated. Knee-ligament forces during straight a-p draw are determined solely by the changing geometry of the ligaments relative to the bones: ACL force decreases with increasing flexion during anterior draw because the angle between the ACL and the tibial plateau decreases as knee flexion increases; PCL force increases with increasing flexion during posterior draw because the angle between the PCL and the tibial plateau increases. The pattern of ligament loading during activity is governed by the geometry of the muscles spanning the knee: the resultant force in the ACL during isometric knee extension is determined mainly by the changing orientation of the patellar tendon relative to the tibia in the sagittal plane; the resultant force in the PCL during isometric knee flexion is dominated by the angle at which the hamstrings meet the tibia in the sagittal plane.     

A Three-Dimensional Musculoskeletal Model of the Human Knee Joint. Part 2: Analysis of Ligament Function

A three-dimensional model of the knee is used to study ligament function during anterior-posterior (a-p) draw, axial rotation, and isometric contractions of the extensor and flexor muscles. The geometry of the model bones is based on cadaver data. The contacting surfaces of the femur and tibia are modeled as deformable; those of the femur and patella are assumed to be rigid. Twelve elastic elements are used to describe the geometry and mechanical properties of the cruciate ligaments, the collateral ligaments, and the posterior capsule. The model is actuated by thirteen musculotendinous units, each unit represented as a three-element muscle in series with tendon. The calculations show that the forces applied during a-p draw are substantially different from those applied by the muscles during activity. Principles of knee-ligament function based on the results of in vitro experiments may therefore be overstated. Knee-ligament forces during straight a-p draw are determined solely by the changing geometry of the ligaments relative to the bones: ACL force decreases with increasing flexion during anterior draw because the angle between the ACL and the tibial plateau decreases as knee flexion increases; PCL force increases with increasing flexion during posterior draw because the angle between the PCL and the tibial plateau increases. The pattern of ligament loading during activity is governed by the geometry of the muscles spanning the knee: the resultant force in the ACL during isometric knee extension is determined mainly by the changing orientation of the patellar tendon relative to the tibia in the sagittal plane; the resultant force in the PCL during isometric knee flexion is dominated by the angle at which the hamstrings meet the tibia in the sagittal plane.     

A dynamic model of the knee and lower limb for simulating rising movements

A two-dimensional dynamical model of the human body was developed and used to simulate muscle and knee-ligament loading during a fast rising movement. The hip, ankle, and toes were each modeled as a simple hinge joint. Relative movements of the femur, tibia, and patella in the sagittal plane were described using a more detailed representation of the knee. The geometry of the model bones was adapted from cadaver data. Eleven elastic elements described the geometric and mechanical properties of the knee ligaments and joint capsule. The patella was assumed to be massless. Smooth hypersurfaces were constructed and used to calculate the position and orientation of the patella during a forward integration of the model. Each hypersurface was formed by applying the principle of static equilibrium to approximate patellofemoral mechanics during the simulation. The model was actuated by 22 musculotendinous units, each unit represented as a three-element muscle in series with tendon. A first-order process was assumed to model muscle excitation-contraction dynamics. Dynamic optimization theory was used to calculate the pattern of muscle excitations that produces a coordinated rising movement from an initial squatting position in minimum time. The calculations support the contention that squatting is a relatively safe exercise for rehabilitation following ACL reconstruction. ACL forces remain less than 20 N for the duration of the task.     

Two-dimensional dynamic modelling of human knee joint

A mathematical dynamic model of the two-dimensional representation of the knee joint is presented. The profiles of the joint surfaces are determined from X-ray films and they are represented by polynomials. The joint ligaments are modelled as nonlinear elastic springs of realistic stiffness properties. Nonlinear equations of motion coupled with nonlinear constraint conditions are solved numerically. Time derivatives are approximated by Newmark difference formulae and the resulting nonlinear algebraic equations are solved employing the Newton-Raphson iteration scheme. Several dynamic loads are applied to the center of mass of the tibia and the ensuing motion is investigated. Numerical results on ligament forces, contact point locations between femur and tibia, and the orientation of tibia relative to femur are presented. The results are shown to be consistent with the anatomy of the knee joint.     

Bone adaptation to altered loading after spinal cord injury: a study of bone and muscle strength

Bone loss from the paralysed limbs after spinal cord injury (SCI) is well documented. Under physiological conditions, bones are adapted to forces which mainly emerge from muscle pull. After spinal cord injury (SCI), muscles can no longer contract voluntarily and are merely activated during spasms. Based on the Ashworth scale, previous research has suggested that these spasms may mitigate bone losses. We therefore wished to assess muscle forces after SCI with a more direct measure and compare it to measures of bone strength. We hypothesized that the bones in SCI patients would be in relation to the loss of muscle forces. Six male patients with SCI 6.4 (SD 4.3) years earlier and 6 age-matched, able-bodied control subjects were investigated. Bone scans from the right knee were obtained by pQCT. The knee extensor muscles were electrically stimulated via the femoral nerve, isometric knee extension torque was measured and patellar tendon force was estimated. Tendon force upon electrical stimulation in the SCI group was 75% lower than in the control subjects (p<0.01). Volumetric bone mineral density of the patella and of the proximal tibia epiphysis were 50% lower in the SCI group than in the control subjects (p<0.01). Cortical area was lower by 43% in the SCI patients at the proximal tibia metaphysis, and by 33% at the distal femur metaphysis. No group differences were found in volumetric cortical density. Close curvilinear relationships were found between stress and volumetric density for the tibia epiphysis (r(2)=0.90) and for the patella (r(2)=0.91). A weaker correlation with the tendon force was found for the cortical area of the proximal tibia metaphysis (r(2)=0.63), and none for the distal femur metaphysis. These data suggest that, under steady state conditions after SCI, epiphyseal bones are well adapted to the muscular forces. For the metaphysis of the long bones, such an adaptation appears to be less evident. The reason for this remains unclear.     

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

A procedure for estimating the relevant forces in the human knee using a four-bar mechanism

Knee injuries, especially those that affect the cruciate and lateral ligaments, are one of the most serious and frequent pathologies that affect the lower human extremity. Hence, the aim of this study is to develop a dynamic model for the lower extremity capable of estimating forces, forces in the cruciate and collateral ligaments and those normal to the articular cartilage, generated in the knee. The proposed model considers a four-bar mechanism in the knee, a spherical joint in the pelvis and a revolute one in the ankle. The four-bar mechanism is obtained by a synthesis process. The dynamic model includes the inertial properties of the femur, tibia, patella and the foot, the ground reaction force and the most important muscles in the knee. Muscle forces are estimated using an optimisation technique. Results from the application of the model on a real human task are presented.     

Computational wear simulation of patellofemoral articular cartilage during in vitro testing

Though changes in normal joint motions and loads (e.g., following anterior cruciate ligament injury) contribute to the development of knee osteoarthritis, the precise mechanism by which these changes induce osteoarthritis remains unknown. As a first step toward identifying this mechanism, this study evaluates computational wear simulations of a patellofemoral joint specimen wear tested on a knee simulator machine. A multibody dynamic model of the specimen mounted in the simulator machine was constructed in commercial computer-aided engineering software. A custom elastic foundation contact model was used to calculate contact pressures and wear on the femoral and patellar articular surfaces using geometry created from laser scan and MR data. Two different wear simulation approaches were investigated--one that wore the surface geometries gradually over a sequence of 10 one-cycle dynamic simulations (termed the "progressive" approach), and one that wore the surface geometries abruptly using results from a single one-cycle dynamic simulation (termed the "non-progressive" approach). The progressive approach with laser scan geometry reproduced the experimentally measured wear depths and areas for both the femur and patella. The less costly non-progressive approach predicted deeper wear depths, especially on the patella, but had little influence on predicted wear areas. Use of MR data for creating the articular and subchondral bone geometry altered wear depth and area predictions by at most 13%. These results suggest that MR-derived geometry may be sufficient for simulating articular cartilage wear in vivo and that a progressive simulation approach may be needed for the patella and tibia since both remain in continuous contact with the femur.     

A constraint-based approach to modelling the mobility of the human knee joint

A model of knee mobility able to predict the range and pattern of movement in the unloaded joint was proposed by Wilson et al. (J. Biomech. 31 (1998) 1127-1136). The articular surfaces in the lateral and medial compartments and isometric fascicles in three of the knee ligaments were represented as five constraints on motion between the femur and tibia in a single degree-of-freedom parallel spatial mechanism. The path of movement of the bones during passive flexion was found by solving the forward kinematics of the mechanism using an iterative method. The present paper shows that such a mechanism-based solution approach can lead to an underestimation of the flexion range. This is due to the mechanism reaching a 'stationary configuration' and 'locking'. A new, constraint-based approach to the solution of the model joint displacement is proposed. It avoids the representation of ligaments and articular surfaces by kinematically equivalent chains of one degree-of-freedom pairs which are prone to singularities. It relies instead on a numerical solution of five non-linear constraint equations to find the relative positions of the bones at a series of flexion angles. The method is successful both in its ability to predict motion through a physiological range and in its efficiency with a solution rate forty times faster than the original algorithm. The new approach may be extended to include more complex joint surface geometry, allowing a study of the effects of articular surface shape and ligament arrangement on joint kinematics.     

Control exerted by ligaments

The function of the ligaments as local controllers, independent of the central nervous system, in maintaining the integrity of the joint is demonstrated by modelling the human knee in the sagittal plane, and studying its anterior-posterior motion. In addition to the ligaments, the model includes the characteristic geometry of the joint surface and some muscle groups. The connecting reaction forces at the point of contact between the tibia and the femur are considered to be constraint forces due to three different surface motions--gliding, rolling and combined gliding and rolling. It is demonstrated that the ligamentous structure maintains these holonomic and nonholonomic constraints that describe the joint motion, and that stability of the knee joint is provided mainly by ligaments. Muscular structures further stabilize and contribute to joint movement. Computer simulation of rolling movement of the knee is presented to illustrate the importance of the ligaments for joint integrity and stability.     

The orientation of the distal part of the quadriceps femoris muscle as a function of the knee flexion-extension angle

Lateral view radiographs of ten autopsy knees were used to determine the orientation of the patellar ligament, patella and quadriceps tendon relative to tibia and femur at different flexion-extension angles (0-120 degrees) of the knee. The results show a linear relationship between the angle of flexion and the movement of the patellar ligament relative to the tibia and of the movement of the patella relative to tibia and femur. There is a non-linear relationship between angle of flexion and the movement of the quadriceps tendon relative to the patellar ligament, patella and femur. The angular changes between patella and patellar ligament are negligible. The complicated movements of the distal part of the quadriceps femoris muscle may significantly influence biomechanical parameters such as the forces acting at the patella and tibial tuberosity.     

Influence of patellofemoral articular geometry and material on mechanics of the unresurfaced patella

Patellar resurfacing during knee replacement is still under debate, with several studies reporting higher incidence of anterior knee pain in unresurfaced patellae. Congruency between patella and femur impacts the mechanics of the patellar cartilage and strain in the underlying bone, with higher stresses and strains potentially contributing to cartilage wear and anterior knee pain. The material properties of the articulating surfaces will also affect load transfer between femur and patella. The purpose of this study was to evaluate the mechanics of the unresurfaced patella and compare with natural and resurfaced conditions in a series of finite element models of the patellofemoral joint. In the unresurfaced analyses, three commercially available implants were compared, in addition to an 'ideal' femoral component which replicated the geometry, but not the material properties, of the natural femur. Hence, the contribution of femoral component material properties could be assessed independently from geometry changes. The ideal component tracked the kinematics and patellar bone strain of the natural knee, but had consistently inferior contact mechanics. In later flexion, compressive patellar bone strain in unresurfaced conditions was substantially higher than in resurfaced conditions. Understanding how femoral component geometry and material properties in unresurfaced knee replacement alters cartilage contact mechanics and bone strain may aid in explaining why the incidence of anterior knee pain is higher in the unresurfaced population, and ultimately contribute to identifying criteria to pre-operatively predict which patients are suited to an unresurfaced procedure and reducing the incidence of anterior knee pain in the unresurfaced patient population.     

Articular-surface-based automatic anatomical coordinate systems for the knee bones

Increasing use of patient-specific surgical procedures in orthopaedics means that patient-specific anatomical coordinate systems (ACSs) need to be determined. For knee bones, automatic algorithms constructing ACSs exist and are assumed to be more reliable than manual methods, although both approaches are based on non-unique numerical reconstructions of true bone geometries. Furthermore, determining the best algorithms is difficult, as algorithms are evaluated on different datasets. Thus, in this study, we developed 3 algorithms, each with 3 variants, and compared them with 5 from the literature on a dataset comprising 24 lower-limb CT-scans. To evaluate algorithms’ sensitivity to the operator-dependent reconstruction procedure, the tibia, patella and femur of each CT-scan were each reconstructed once by three different operators.Our algorithms use principal inertia axis (PIA), cross-sectional area, surface normal orientations and curvature data to identify the bone region underneath articular surfaces (ASs). Then geometric primitives are fitted to ASs, and the ACSs are constructed from the geometric primitive points and/or axes. For each bone type, the algorithm displaying the least inter-operator variability is identified. The best femur algorithm fits a cylinder to posterior condyle ASs and a sphere to the femoral head, average axis deviations: 0.12°, position differences: 0.20 mm. The best patella algorithm identifies the AS PIAs, average axis deviations: 0.91°, position differences: 0.19 mm. The best tibia algorithm finds the ankle AS center and the 1st PIA of a layer around a plane fitted to condyle ASs, average axis deviations: 0.38°, position differences: 0.27 mm.     

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.     

Effects of joint unloading and reloading on human cartilage morphology and function, muscle cross-sectional areas, and bone density - a quantitative case report

Recent studies have shown that thinning of human cartilage occurs with unloading, but no data are available on the effect of remobilization (after immobilization) on knee joint cartilage status in humans. We examined a 36-year-old patient after 6 weeks of unilateral immobilization. Knee joint cartilage morphology (patella and tibia), patellar cartilage deformation, and thigh muscle cross-sectional areas were assessed with quantitative MR imaging and bone density with peripheral quantitative computed tomography (pQCT) during 24 months of remobilization. The immobilized limb displayed lower muscle cross-sectional areas (MCSA) of the knee extensors (-36%), lower bone density of the femur and tibia (-12/-6%), lower patellar cartilage thickness (-14%), but no side differences of tibial cartilage thickness. During remobilization, side differences decreased to -4% for knee extensor MCSAs, to -6%/-3% for femoral and tibial BMD, and to -8% for patellar cartilage thickness. No change was observed in tibial cartilage. Patellar deformation decreased from 9% to 4% after 15 months. In conclusion, we observed substantial changes of thigh MCSAs, but little (patella) to no (tibia) change in cartilage thickness during remobilization. These preliminary results indicate that human cartilage macro-morphology may be less adaptive to variations of the mechanical loading than muscle and bone.     

Ligaments and articular contact guide passive knee flexion

The aim of this study was to test the hypothesis that the coupled features of passive knee flexion are guided by articular contact and by the isometric fascicles of the ACL, PCL and MCL. A three-dimensional mathematical model of the knee was developed, in which the articular surfaces in the lateral and medial compartments and the isometric fascicles in the ACL, PCL and MCL were represented as five constraints in a one degree-of-freedom parallel spatial mechanism. Mechanism analysis techniques were used to predict the path of motion of the tibia relative to the femur. Using a set of anatomical parameters obtained from a cadaver specimen, the model predicts coupled internal rotation and ab/adduction with flexion. These predictions correspond well to measurements of the cadaver specimen’s motion. The model also predicts posterior translation of contact on the tibia with flexion. Although this is a well-known feature of passive knee flexion, the model predicts more translation than has been reported from experiments in the literature. Modelling of uncertainty in the anatomical parameters demonstrated that the discrepancy between theoretical predictions and experimental measurement can be attributed to parameter sensitivity of the model. This study shows that the ligaments and articular surfaces work together to guide passive knee motion. A principal implication of the work is that both articular surface geometry and ligament geometry must be preserved or replicated by surgical reconstruction and replacement procedures to ensure normal knee kinematics and by extension, mechanics.     

A three-dimensional anatomical model of the human patello-femoral joint,for the determination of patello-femoral motions and contact characteristics

The object of this study is to develop a three-dimensional mathematical model of the patello-femoral joint, which is modelled as two rigid bodies representing a moving patella and a fixed femur. Two-point contact was assumed between the femur and patella at the medial and lateral sides and in the analysis, the femoral and patellar articular surfaces were mathematically represented using Coons' bicubic surface patches. Model equations include six equilibrium equations and eleven constraints: six contact conditions, four geometric compatibility conditions, and the condition of a rigid patellar ligament; the model required the solution of a system of 17 nonlinear equations in 17 unknowns, its response describing the six-degress-of-freedom patellar motions and the forces acting on the patella. Patellar motions are described by six motion parameters representing the translations and rotations of the patella with respect to the femur. The forces acting on the patella include the medial and lateral component of patello-femoral contact and the patellar ligament force, all of which were represented as ratios to the quadriceps tendon force. The model response also includes the locations of the medial and lateral contact points on the femur and the patella. A graphical display of its response was produced in order to visualize better the motion of the components of the extensor mechanism.Model calculations show good agreement with experimental results available from the literature. The patella was found to move distally and posteriorly on the femoral condyles as the knee was flexed from full extension. Results indicate that the relative orientation of the patellar ligament with respect to the patella remains unchanged during this motion. The model also predicts a patellar flexion which always lagged knee flexion.Our calculations show that as the angle of knee flexion increased, the lateral contact point moved distally on the femur without moving significantly either medially or laterally. The medial contact point also moved distally on the femur but moved medially from full extension to about 40° of knee flexion, then laterally as the knee flexion angle increased. The lateral contact point on the patella did not change significantly in the medial and lateral direction as the knee was flexed; however, this point moved proximally toward the basis of the patella with knee flexion. The medial contact point also moved proximally on the patella with knee flexion, and in a similar manner the medial contact point on the patella moved distally with flexion from full extension to about 40° of flexion. However, as the angle of flexion increased, the medial contact point did not move significantly in the medial-lateral direction.Model calculations also show that during the simulated knee extension exercise, the ratio of the force in the patellar ligament to the force in the quadriceps tendon remains almost unchanged for the first 30° of knee flexion, then decreases as the angle of knee flexion increases. Furthermore, model results show that the lateral component of the patello-femoral contact force is always greater than the medial component, both components increasing with knee flexion.     

Deformation of the distal femur: a contribution towards the pathogenesis of osteochondrosis dissecans in the knee joint.

Osteochondrosis dissecans (OD) is a process of subchondral bone necrosis occurring predominantly in young individuals at specific sites. The aetiology of this disease remains controversial with mechanical processes due to trauma and/or ischaemic factors being proposed. This study aims at explaining the aetiology of OD in the knee joint as a result of the particular deformation of the condyles. A finite element analysis of the distal third of the femur was performed. A three-dimensional model was developed based on computed tomography scans of a normal femur, consisting of cortical bone, cancellous bone and articular cartilage. This model was subjected to physiological loads at 0, 30, 60 and 90 degrees of knee flexion. A complex deformation was found within each condyle as well as between the two condyles. Both medial and lateral condyles are deformed in the medio-lateral direction and at the same time compressed between the patella and the tibia in the antero-posterior direction. This effect is highest at 60 degrees of knee flexion. In both planes, the medial condyle is distorted more than the lateral one. Strain concentration in the subchondral bone facing the patella varies with flexion, especially for angles exceeding 60 degrees. The deformation of the femur in the predominant locus of OD in the medial condyle exceeds that of the lateral condyle considerably. The analysis shows that repeated vigorous exercise including extreme knee flexion may produce rapidly changing strains which in turn could ultimately be responsible for local subchondral bone collapse.