The biomechanical function of the patellar tendon during in-vivo weight-bearing flexion

Few studies have investigated the function of the patellar tendon in-vivo. This study quantified the three-dimensional (3D) kinematics of the patellar tendon during weight-bearing flexion. Eleven subjects were imaged using magnetic resonance (MR). Sagittal plane images were outlined to create a 3D model of the patella, tibia, and femur and included the attachment sites of the patellar tendon. Each attachment site was divided into central, medial, and lateral thirds. Next, the subjects were imaged using fluoroscopy from two orthogonal directions while performing a single-leg lunge. The models and fluoroscopic images were used to reproduce the motion of the patella, tibia, and femur. The apparent elongation, sagittal plane angle, and coronal plane angle of each third of the patellar tendon were measured from the relative motion of the attachment sites. All three portions of the patellar tendon deformed similarly with flexion. The length of the patellar tendon significantly from full extension to 30 degrees . From 30 degrees -110 degrees , no significant change in the length of the patellar tendon was observed. The patellar tendon was oriented anteriorly at flexion angles less than 60 degrees and posteriorly thereafter. From full extension to 60 degrees , the medial orientation of the patellar tendon decreased significantly with flexion. These data may have important implications for anterior cruciate ligament reconstruction using patellar tendon autografts and for the design of rehabilitation regimens for patients of patellar tendon repair.     

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.     

A planar model of the knee joint to characterize the knee extensor mechanism

A simple planar static model of the knee joint was developed to calculate effective moment arms for the quadriceps muscle. A pathway for the instantaneous center of rotation was chosen that gives realistic orientations of the femur relative to the tibia. Using the model, nonlinear force and moment equilibrium equations were solved at one degree increments for knee flexion angles from 0 (full extension) to 90 degrees, yielding patellar orientation, patellofemoral contact force and patellar ligament force and direction with respect to both the tibial insertion point and the tibiofemoral contact point. The computer-derived results from this two-dimensional model agree with results from more complex models developed previously from experimentally obtained data. Due to our model's simplicity, however, the operation of the patellar mechanism as a lever as well as a spacer is clearly illustrated. Specifically, the thickness of the patella was found to increase the effective moment arm significantly only at flexions below 35 degrees even though the actual moment arm exhibited an increase throughout the flexion range. Lengthening either the patella or the patellar ligament altered the force transmitted from the quadriceps to the patellar ligament, significantly increasing the effective moment arm at flexions greater than 25 degrees. We conclude that the levering action of the patella is an essential mechanism of knee joint operation at moderate to high flexion angles.     

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

Influence of patella alta on knee extensor mechanics

The purpose of this study was to compare the knee extensor mechanics in persons with and without patella alta. Thirteen subjects with patella alta and 14 subjects with normal patellar position participated in the study. Sagittal and axial MR images of the knee were acquired at 0°, 20°, 40°, and 60° of knee flexion. Measurements of actual moment arm, patellar ligament/quadriceps tendon force ratio, quadriceps effective moment arm, and joint reaction force/quadriceps force ratio were obtained. There were no differences between groups in terms of actual moment arm. However, subjects with patella alta had significantly larger patellar ligament/quadriceps tendon force ratios (1.04±0.02 vs. 0.92±0.02) and quadriceps effective moment arms (4.40±0.09 vs. 4.00±0.09 cm) when averaged across the range of knee flexion angles tested. There was no difference in the joint reaction force/quadriceps force ratio between groups. The observed differences in knee extensor mechanics suggest that individuals with patella alta have a more efficient knee extensor mechanism and would be expected to generate similar joint reaction forces per unit quadriceps force compared to subjects with normal patellar position. Therefore, persons with patella alta may experience less patellofemoral joint reaction force to overcome the same knee flexion moment in the range of 0°–60° of knee flexion.     

In vitro forces in the normal and cruciate-deficient knee during simulated squatting motion

Three orthogonal components of the tibiofemoral and patellofemoral forces were measured simultaneously for knees with intact cruciate ligaments (nine knees), following anterior cruciate ligament resection (six knees), and subsequent posterior cruciate ligament resection (six knees). The knees were loaded using an experimental protocol that modeled static double-leg squat. The mean compressive tibial force increased with flexion angle. The mean anteroposterior tibial shear force acted posteriorly on the tibia below 50 deg flexion and anteriorly above 55 deg. Mediolateral shear forces were low compared to the other force components and tended to be directed medially on both the patella and tibia. The mean value of the ratio of the resultant tibial force divided by the quadriceps force decreased with increasing flexion angle and was between 0.6 and 0.7 above 70 deg flexion. The mean value of the ratio of the resultant tibiofemoral contact force divided by the resultant patellofemoral contact force decreased with increasing flexion and was between 0.8 and 1.0 above 55 deg flexion. Cruciate ligament resection resulted in no significant changes in the patellar contact forces. Following resection of the anterior cruciate ligament, the tibial anteroposterior shear force was directed anteriorly over all flexion angles tested. Subsequent resection of the posterior cruciate ligament resulted in an approximately 10 percent increase in the quadriceps tendon and tibial compressive force.     

In vivo tracking of the human patella.

The purpose of this study was to describe the dynamic, in vivo, three-dimensional tracking pattern of the patella for one normal male subject. Intracortical pins were inserted into the patella, tibia, and femur. The subject performed seated and squatting knee flexion/extension, and maximum voluntary quadriceps contractions. In addition, the vastus medialis oblique was subjected to maximal electrical stimulation. Motions of the markers attached to the intracortical pins were analyzed using an automated video system. Patellar and tibial motions were determined relative to a femoral reference system. While the tibia flexed 50 degrees from full extension (seated condition), the patella flexed 30.3 degrees, tilted laterally 10.3 degrees, and shifted laterally 8.6 mm. In general, these results show qualitative agreement with the data collected from cadaveric specimens [van Kampen and Huiskes, J. orthop. Res. 8, 372-382 (1990)]. The differences present may reflect different passive constraints to patellar motions, and different relative loading of the individual quadriceps components, in our study compared to the cadaveric study. Only small differences were found between patellar motions in the seated and squatting conditions. Differences in patellar displacements produced by (1) maximal electrical stimulation of the vastus medialis oblique, and (2) maximum voluntary quadriceps contraction, at 30 degrees knee flexion and full extension, may reflect the dominant influence of passive constraints, and the vastus lateralis, on normal patellar motions. Further in vivo study of patellar tracking seems warranted to evaluate surgical and conservative interventions for patellofemoral disorders.     

Moment arm of the patellar tendon in the human knee

The moment arm of the knee-extensor mechanism is described by the moment arm of the patellar tendon calculated with respect to the screw axis of the tibia relative to the femur. The moment arm may be found once the line of action of the patellar tendon and the position and orientation of the screw axis are known. In this study, the orientation of the patellar tendon and the position and orientation of the finite screw axis of the tibia relative to the femur were calculated from measurements of the three-dimensional positions of the bones obtained from fresh cadaver specimens. Peak values of the patellar tendon moment arm ranged from 4-6 cm for the six knees tested; the moment arm was maximum near 45 degrees of knee flexion. The moment arm of the patellar tendon was nearly equal to the shortest (perpendicular) distance between the line of action of the patellar tendon and the axis of rotation of the knee at all flexion angles, except near full extension. Near full extension, the angle between the patellar tendon and the screw axis was significantly less than 90 degrees, and the magnitude of the moment arm was then less than the perpendicular distance between these two lines. The patellar tendon moment arm remained roughly constant across individuals when normalized by femoral condyle width, suggesting that anatomical differences play a large role in determining the moment arm of the extensor mechanism.     

In vivo assessment of the interaction of patellar tendon tibial shaft angle and anterior cruciate ligament elongation during flexion

A potential cause of non-contact anterior cruciate ligament (ACL) injury is landing on an extended knee. In line with this hypothesis, studies have shown that the ACL is elongated with decreasing knee flexion angle. Furthermore, at low flexion angles the patellar tendon is oriented to increase the anterior shear component of force acting on the tibia. This indicates that knee extension represents a position in which the ACL is taut, and thus may have an increased propensity for injury, particularly in the presence of excessive force acting via the patellar tendon. However, there is very little in vivo data to describe how patellar tendon orientation and ACL elongation interact during flexion. Therefore, this study measured the patellar tendon tibial shaft angle (indicative of the relative magnitude of the shear component of force acting via the patellar tendon) and ACL length in vivo as subjects performed a quasi-static lunge at varying knee flexion angles. Spearman rho rank correlations within each individual revealed that flexion angles were inversely correlated to both ACL length (rho = −0.94 ± 0.07, mean ± standard deviation, p < 0.05) and patellar tendon tibial shaft angle (rho = −0.99 ± 0.01, p < 0.05). These findings indicate that when the knee is extended, the ACL is both elongated and the patellar tendon tibial shaft angle is increased, resulting in a relative increase in anterior shear force on the tibia acting via the patellar tendon. Therefore, these data support the hypothesis that landing with the knee in extension is a high risk scenario for ACL injury.     

Experimental determination of forces transmitted through the patello-femoral joint

Tensions in the quadriceps tendon and infrapatellar ligament were measured as a function of flexion angle in eight cadaver knees using a load cell of a materials tester to determine the quadriceps force and a spring balance to quantify the patellar tendon force. The ratio between the tensions in the quadriceps tendon and the patellar tendon (FQ/FP) ranged from 1.55 at 70 degrees of flexion to 0.86 at 10 degrees of flexion. The patello-femoral joint reaction (PFJR) force for extension against resistance was maximal at 60 degrees. No change in the quadriceps force required to extend the knee occurred with changes of the Q-angle of +/- 5 degrees. This study demonstrates that FQ does not equal FP as several authors have reported (Bandi, 1972; Barry, 1979; Ficat and Hungerford, 1977; Hungerford and Barry, 1979; Reilly and Martens, 1972; Smidt, 1973). Furthermore, the difference in FQ and FP influences both the magnitude and direction of PFJR. Studies that assess the influence of surgical procedures which alter the patello-femoral joint or the extensor mechanism must take these differences into account.     

The effect of design variables of condylar total knees on the joint forces in step climbing based on a computer model.

The ability to climb a steep step or rise from a low chair after total knee replacement may be enhanced if the required force in the quadriceps muscle is reduced. This can potentially be achieved if the total knee produces a large lever arm measured from the femoral-tibial contact point to the patellar ligament. A reduced quadriceps force would also reduce the patello-femoral force and the femoral-tibial contact force. The contact point location is likely to be a function of the geometry of the femoral and tibial components in the sagittal plane, including the relative distal and posterior radii of the femoral profile, the location of the bottom-of-the-dish of the tibial surface, the radius of the tibial surface, and the presence or absence of the posterior cruciate ligament. A three-dimensional model of the knee was developed including the quadriceps and various ligaments. In the study, the motion was confined to flexion extension and displacement in the sagittal plane. The quadriceps was assumed to be the only muscle acting. A standard software package (Pro/Mechanica) was used for the analysis. For a femoral component with a smaller distal radius, there was 12% reduction in the quadriceps muscle force and up to 11% reduction in the patello-femoral force from about 100 up to 60 degrees flexion. However, apart from that, there were less than 10% differences in all the forces as a function of all of the design variables studied. This was attributed to the relatively small changes in the lever arm of the patella tendon, since the tendon moves in an anterior-posterior direction along with the femur. An additional factor explaining the results was the change in the anterior-posterior contact point as controlled by the forces in the patella tendon and in the soft tissues. The results imply that for a standard condylar replacement knee, the muscle and contact forces are not greatly affected by the geometrical design variables.     

Influence of patellar position on the knee extensor mechanism in normal and crouched walking

Patella alta is common in cerebral palsy, especially in patients with crouch gait. Correction of patella alta has been advocated in the treatment of crouch, however the appropriate degree of correction and the implications for knee extensor function remain unclear. Therefore, the goal of this study was to assess the impact of patellar position on quadriceps and patellar tendon forces during normal and crouch gait. To this end, a lower extremity musculoskeletal model with a novel 12 degree of freedom knee joint was used to simulate normal gait in a healthy child, as well as mild (23 deg min knee flexion in stance), moderate (41 deg), and severe (67 deg) crouch gait in three children with cerebral palsy. The simulations revealed that quadriceps and patellar tendon forces increase dramatically with crouch, and are modulated by patellar position. For example with a normal patellar tendon position, peak patellar tendon forces were 0.7 times body weight in normal walking, but reached 2.2, 3.2 and 5.4 times body weight in mild, moderate and severe crouch. Moderate patella alta acted to reduce quadriceps and patellar tendon loads in crouch gait, due to an enhancement of the patellar tendon moment arms with alta in a flexed knee. In contrast, patella baja reduced the patellar tendon moment arm in a flexed knee and thus induced an increase in the patellar tendon loads needed to walk in crouch. Functionally, these results suggest that patella baja could also compromise knee extensor function for other flexed knee activities such as chair rise and stair climbing. The findings are important to consider when using surgical approaches for correcting patella alta in children who exhibit crouch gait patterns.     

A detailed and validated three dimensional dynamic model of the patellofemoral joint

A detailed 3D anatomical model of the patellofemoral joint was developed to study the tracking, force, contact and stability characteristics of the joint. The quadriceps was considered to include six components represented by 15 force vectors. The patellar tendon was modeled using four bundles of viscoelastic tensile elements. Each of the lateral and medial retinaculum was modeled by a three-bundle nonlinear spring. The femur and patella were considered as rigid bodies with their articular cartilage layers represented by an isotropic viscoelastic material. The geometrical and tracking data needed for model simulation, as well as validation of its results, were obtained from an in vivo experiment, involving MR imaging of a normal knee while performing isometric leg press against a constant 140 N force. The model was formulated within the framework of a rigid body spring model and solved using forth-order Runge-Kutta, for knee flexion angles between zero and 50 degrees. Results indicated a good agreement between the model predictions for patellar tracking and the experimental results with RMS deviations of about 2 mm for translations (less than 0.7 mm for patellar mediolateral shift), and 4 degrees for rotations (less than 3 degrees for patellar tilt). The contact pattern predicted by the model was also consistent with the results of the experiment and the literature. The joint contact force increased linearly with progressive knee flexion from 80 N to 210 N. The medial retinaculum experienced a peak force of 18 N at full extension that decreased with knee flexion and disappeared entirely at 20 degrees flexion. Analysis of the patellar time response to the quadriceps contraction suggested that the muscle activation most affected the patellar shift and tilt. These results are consistent with the recent observations in the literature concerning the significance of retinaculum and quadriceps in the patellar stability.     

A mathematical model of the patellofemoral joint

A mathematical model of the patellofemoral joint taking into account movements and forces in the sagittal plane is described. The system parameters of the model are the locations of the attachments of the quadriceps muscle and the patellar ligament, the length of the patellar ligament, the dimensions of the patella and the geometry of the articulating surfaces. They were obtained from ten autopsy knees. The model enables calculation of the relative position of the patella, patellar ligament and quadriceps tendon, the location of the patellofemoral contact point and the magnitude of the patellofemoral compression force and the force in the patellar ligament as a function of the location of the tibial tuberosity at different flexion-extension angles of the knee. The model is validated by comparing model data with experimentally determined data.     

In vivo patellar tracking induced by individual quadriceps components in individuals with patellofemoral pain

Patellofemoral pain is a common knee disorder with a multi-factorial etiology related to abnormal patellar tracking. Our hypothesis was that the pattern of three-dimensional rotation and translation of the patella induced by selective activation of individual quadriceps components would differ between subjects with patellofemoral pain and healthy subjects. Nine female subjects with patellofemoral pain and seven healthy female subjects underwent electrical stimulation to selectively activate individual quadriceps components (vastus medialis obliquus, VMO; vastus medialis lateralis, VML; vastus lateralis, VL) with the knee at 0° and 20° flexion, while three-dimensional patellar tracking was recorded. Normalized direction of rotation and direction of translation characterized the relative amplitudes of each component of patellar movement. VMO activation in patellofemoral pain caused greater medial patellar rotation (distal patellar pole rotates medially in frontal plane) at both knee positions (p<0.01), and both VMO and VML activation caused increased anterior patellar translation (p<0.001) in patellofemoral pain compared to healthy subjects at 20° knee flexion. VL activation caused more lateral patellar translation (p<0.001) in patellofemoral pain compared to healthy subjects. In healthy subjects the 3-D mechanical action of the VMO is actively modulated with knee flexion angle while such modulation was not observed in PFP subjects. This could be due to anatomical differences in the VMO insertion on the patella and medial quadriceps weakness. Quantitative evaluation of the influence of individual quadriceps components on patellar tracking will aid understanding of the knee extensor mechanism and provide insight into the etiology of patellofemoral pain.     

3-D anatomically based dynamic modeling of the human knee to include tibio-femoral and patello-femoral joints

An anatomical dynamic model consisting of three body segments, femur, tibia and patella, has been developed in order to determine the three-dimensional dynamic response of the human knee. Deformable contact was allowed at all articular surfaces, which were mathematically represented using Coons' bicubic surface patches. Nonlinear elastic springs were used to model all ligamentous structures. Two joint coordinate systems were employed to describe the six-degrees-of-freedom tibio-femoral (TF) and patello-femoral (PF) joint motions using twelve kinematic parameters. Two versions of the model were developed to account for wrapping and nonwrapping of the quadriceps tendon around the femur. Model equations consist of twelve nonlinear second-order ordinary differential equations coupled with nonlinear algebraic constraint equations resulting in a Differential-Algebraic Equations (DAE) system that was solved using the Differential/Algebraic System Solver (DASSL) developed at Lawrence Livermore National Laboratory. Model calculations were performed to simulate the knee extension exercise by applying non-linear forcing functions to the quadriceps tendon. Under the conditions tested, both "screw home mechanism" and patellar flexion lagging were predicted. Throughout the entire range of motion, the medial component of the TF contact force was found to be larger than the lateral one while the lateral component of the PF contact force was found to be larger than the medial one. The anterior and posterior fibers of both anterior and posterior cruciate ligaments, ACL and PCL, respectively, had opposite force patterns: the posterior fibers were most taut at full extension while the anterior fibers were most taut near 90 degrees of flexion. The ACL was found to carry a larger total force than the PCL at full extension, while the PCL carried a larger total force than the ACL in the range of 75 degrees to 90 degrees of flexion.     

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 1: Theoretical Construct

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.