Some aspects of rat femorotibial joint microanatomy as demonstrated by high-resolution magnetic resonance imaging

High-resolution magnetic resonance images (MRI) of the right femorotibial joint of normal Han:Wistar rats were acquired using a 4.7 Tesla magnet and a single-turn solenoid radio frequency coil (built in-house). Some anatomical findings of the rat femorotibial joint, which have not been reported previously using MRI, are described. The separation of patellar ligament and crural fascia was feasible on MRI. This separation would not be seen on images of lower resolution and its presence on high-resolution images could be mistaken for artefact due to the magic angle effect. Band-like fibrous structures exist in the infra-patellar fat pad, which might be mistaken as ligaments within the femorotibial joint. On sagittal MRI a vessel was seen inserted on the central part of the caudal surface of the patellar ligament. Subcutaneous fascia/cutaneous muscles (panniculus carnosus) could also be demonstrated with MRI in the femorotibial joint area.     

Ex Vivo Pathomechanics of the Canine Pond-Nuki Model

Background

Transection of the canine cranial cruciate ligament (CCL) is a well-established osteoarthritis (OA) model. The effect of CCL loss on contact pressure and joint alignment has not been quantified for stifle loading in standing. The purposes of the study were to measure femorotibial contact areas and stresses and joint alignment following transection of the CCL in an ex vivo model. We hypothesized that transection of the CCL would lead to abnormal kinematics, as well as alterations in contact mechanics of the femorotibial joint.

Methodology/Principal Findings

Eight canine hindlimbs were tested in a servo-hydraulic materials testing machine using a custom made femoral jig. Contact area and pressure measurements, and femorotibial rotations and translations were measured in the normal and the CCL–deficient stifle in both standing and deep flexion angles.We found that at standing angle, transection of the CCL caused cranial translation and internal rotation of the tibia with a concurrent caudal shift of the contact area, an increase in peak pressure and a decrease in contact area. These changes were not noted in deep flexion. At standing, loss of CCL caused a redistribution of the joint pressure, with the caudal region of the compartment being overloaded and the rest of the joint being underloaded.

Conclusion

In the Pond-Nuki model alterations in joint alignment are correlated with shifting of the contact points to infrequently loaded areas of the tibial plateau. The results of this study suggest that this cadaveric Pond-Nuki model simulates the biomechanical changes previously reported in the in-vivo Pond-Nuki model.     

In vivo kinematic evaluation and design considerations related to high flexion in total knee arthroplasty

In designing a posterior-stabilized total knee arthroplasty (TKA) it is preferable that when the cam engages the tibial spine the contact point of the cam move down the tibial spine. This provides greater stability in flexion by creating a greater jump distance and reduces the stress on the tibial spine. In order to eliminate edge loading of the femoral component on the posterior tibial articular surface, the posterior femoral condyles need to be extended. This provides an ideal femoral contact with the tibial articular surface during high flexion angles. To reduce extensor mechanism impingement in deep flexion, the anterior margin of the tibial articular component should be recessed. This provides clearance for the patella and patella tendon. An in vivo kinematic analysis that determined three dimensional motions of the femorotibial joint was performed during a deep knee bend using fluoroscopy for 20 subjects having a TKA designed for deep flexion. The average weight-bearing range-of-motion was 125 degrees . On average, TKA subjects experienced 4.9 degrees of normal axial rotation and all subjects experienced at least -4.4 mm of posterior femoral rollback. It is assumed that femorotibial kinematics can play a major role in patellofemoral kinematics. In this study, subjects implanted with a high-flexion TKA design experienced kinematic patterns that were similar to the normal knee. It can be hypothesized that forces acting on the patella were not substantially increased for TKA subjects compared with the normal subjects.     

Development and validation of a weight-bearing finite element model for total knee replacement

Total knee arthroplasty (TKA) is a successful procedure for osteoarthritis. However, some patients (19%) do have pain after surgery. A finite element model was developed based on boundary conditions of a knee rig. A 3D-model of an anatomical full leg was generated from magnetic resonance image data and a total knee prosthesis was implanted without patella resurfacing. In the finite element model, a restarting procedure was programmed in order to hold the ground reaction force constant with an adapted quadriceps muscle force during a squat from 20° to 105° of flexion. Knee rig experimental data were used to validate the numerical model in the patellofemoral and femorotibial joint. Furthermore, sensitivity analyses of Young’s modulus of the patella cartilage, posterior cruciate ligament (PCL) stiffness, and patella tendon origin were performed. Pearson’s correlations for retropatellar contact area, pressure, patella flexion, and femorotibial ap-movement were near to 1. Lowest root mean square error for retropatellar pressure, patella flexion, and femorotibial ap-movement were found for the baseline model setup with Young’s modulus of 5 MPa for patella cartilage, a downscaled PCL stiffness of 25% compared to the literature given value and an anatomical origin of the patella tendon. The results of the conducted finite element model are comparable with the experimental results. Therefore, the finite element model developed in this study can be used for further clinical investigations and will help to better understand the clinical aspects after TKA with an unresurfaced patella.     

Validation of the Tekscan system for statistic and dynamic pressure measurements of the human femorotibial joint]

In vitro dynamic pressure measurements in the healthy and pathologically altered knee joint help to improve our understanding of the loading pattern on femorotibial surfaces. The aim of the study was to evaluate a piezoresistive pressure measuring system. A human cadaveric knee was mounted in a material-testing machine (Bionix 858) using a specially designed knee-holding device. Axial loading of the knee, flexed at 20o, at 500 N, 1000N and 1500 N was then carried out. For the static investigations, the piezoresistive measuring system (Tekscan), was compared with the FUJI measuring system. In addition, dynamic measurements were also performed with the Tekscan System. With the exception of the lateral compartment at a load of 1500 N, no differences in maximum pressures were observed between the two systems. Nor were there any differences with regard to contact surfaces, either in the medial or lateral compartment (p > 0.05). However, the reproducibility of the data was significantly higher with the Tekscan System (p < 0.01). Dynamic pressure measurements obtained with the knee flexed 20 to 90o showed that the lateral contact area shifted from anterior to posterior, while the medial contact area remained virtually unchanged. The Tekscan System proved to be more reliable than the FUJI System, and permits simultaneous measurements in both compartments. The Tekscan System is suitable for dynamic measurement of the femorotibial joint, and permits measurements to be made under more physiological conditions.     

Contact pressure in the facet joint during sagittal bending of the cadaveric cervical spine

The facet joint contributes to the normal biomechanical function of the spine by transmitting loads and limiting motions via articular contact. However, little is known about the contact pressure response for this joint. Such information can provide a quantitative measure of the facet joint's local environment. The objective of this study was to measure facet pressure during physiologic bending in the cervical spine, using a joint capsule-sparing technique. Flexion and extension bending moments were applied to six human cadaveric cervical spines. Global motions (C2-T1) were defined using infra-red cameras to track markers on each vertebra. Contact pressure in the C5-C6 facet was also measured using a tip-mounted pressure transducer inserted into the joint space through a hole in the postero-inferior region of the C5 lateral mass. Facet contact pressure increased by 67.6 ± 26.9 kPa under a 2.4 Nm extension moment and decreased by 10.3 ± 9.7 kPa under a 2.7 Nm flexion moment. The mean rotation of the overall cervical specimen motion segments was 9.6 ± 0.8° and was 1.6 ± 0.7° for the C5-C6 joint, respectively, for extension. The change in pressure during extension was linearly related to both the change in moment (51.4 ± 42.6 kPa/Nm) and the change in C5-C6 angle (18.0 ± 108.9 kPa/deg). Contact pressure in the inferior region of the cervical facet joint increases during extension as the articular surfaces come in contact, and decreases in flexion as the joint opens, similar to reports in the lumbar spine despite the difference in facet orientation in those spinal regions. Joint contact pressure is linearly related to both sagittal moment and spinal rotation. Cartilage degeneration and the presence of meniscoids may account for the variation in the pressure profiles measured during physiologic sagittal bending. This study shows that cervical facet contact pressure can be directly measured with minimal disruption to the joint and is the first to provide local pressure values for the cervical joint in a cadaveric model.     

Changes in knee function associated with treadmill ambulation

A comparison of level walking, on a walkway and on a treadmill, was performed using ten normal subjects. Motion about the knee was measured using a triaxial electrogoniometer, and foot-floor contact patterns were recorded by means of four foot switches attached to the sole of each shoe. On the walkway, the data were collected with the subject moving at a comfortable walking speed. The treadmill was then set at the average velocity obtained on the walkway. Knee joint rotation in the coronal and transverse planes did not change significantly between the walkway and the treadmill. In the sagittal plane, significant differences were found for total motion (p less than 0.01), swing phase motion (p less than 0.01), knee position at heel strike (p less than 0.05), and maximum swing phase extension (p less than 0.01). A comparison of the foot-floor contact patterns between walkway and treadmill ambulation revealed reduced heel contact time, with an increase in toe contact while on the treadmill. It was concluded that sagittal plane knee kinematics during level treadmill walking differ significantly from level overground walking.     

The effect of resistance level and stability demands on recruitment patterns and internal loading of spine in dynamic flexion and extension using a simple trunk model

The effects of external resistance on the recruitment of trunk muscles in sagittal movements and the coactivation mechanism to maintain spinal stability were investigated using a simple computational model of iso-resistive spine sagittal movements. Neural excitation of muscles was attained based on inverse dynamics approach along with a stability-based optimisation. The trunk flexion and extension movements between 60° flexion and the upright posture against various resistance levels were simulated. Incorporation of the stability constraint in the optimisation algorithm required higher antagonistic activities for all resistance levels mostly close to the upright position. Extension movements showed higher coactivation with higher resistance, whereas flexion movements demonstrated lower coactivation indicating a greater stability demand in backward extension movements against higher resistance at the neighbourhood of the upright posture. Optimal extension profiles based on minimum jerk, work and power had distinct kinematics profiles which led to recruitment patterns with different timing and amplitude of activation.     

The effect of resistance level and stability demands on recruitment patterns and internal loading of spine in dynamic flexion and extension using a simple trunk model

The effects of external resistance on the recruitment of trunk muscles in sagittal movements and the coactivation mechanism to maintain spinal stability were investigated using a simple computational model of iso-resistive spine sagittal movements. Neural excitation of muscles was attained based on inverse dynamics approach along with a stability-based optimisation. The trunk flexion and extension movements between 60° flexion and the upright posture against various resistance levels were simulated. Incorporation of the stability constraint in the optimisation algorithm required higher antagonistic activities for all resistance levels mostly close to the upright position. Extension movements showed higher coactivation with higher resistance, whereas flexion movements demonstrated lower coactivation indicating a greater stability demand in backward extension movements against higher resistance at the neighbourhood of the upright posture. Optimal extension profiles based on minimum jerk, work and power had distinct kinematics profiles which led to recruitment patterns with different timing and amplitude of activation.     

Estimates of muscle function in human gait depend on how foot-ground contact is modelled

Computational analyses of leg-muscle function in human locomotion commonly assume that contact between the foot and the ground occurs at discrete points on the sole of the foot. Kinematic constraints acting at these contact points restrict the motion of the foot and, therefore, alter model calculations of muscle function. The aim of this study was to evaluate how predictions of muscle function obtained from musculoskeletal models are influenced by the model used to simulate ground contact. Both single- and multiple-point contact models were evaluated. Muscle function during walking and running was determined by quantifying the contributions of individual muscles to the vertical, fore-aft and mediolateral components of the ground reaction force (GRF). The results showed that two factors – the number of foot-ground contact points assumed in the model and the type of kinematic constraint enforced at each point – affect the model predictions of muscle coordination. Whereas single- and multiple-point contact models produced similar predictions of muscle function in the sagittal plane, inconsistent results were obtained in the mediolateral direction. Kinematic constraints applied in the sagittal plane altered the model predictions of muscle contributions to the vertical and fore-aft GRFs, while constraints applied in the frontal plane altered the calculations of muscle contributions to the mediolateral GRF. The results illustrate the sensitivity of calculations of muscle coordination to the model used to simulate foot-ground contact.     

Estimates of muscle function in human gait depend on how foot-ground contact is modelled

Computational analyses of leg-muscle function in human locomotion commonly assume that contact between the foot and the ground occurs at discrete points on the sole of the foot. Kinematic constraints acting at these contact points restrict the motion of the foot and, therefore, alter model calculations of muscle function. The aim of this study was to evaluate how predictions of muscle function obtained from musculoskeletal models are influenced by the model used to simulate ground contact. Both single- and multiple-point contact models were evaluated. Muscle function during walking and running was determined by quantifying the contributions of individual muscles to the vertical, fore-aft and mediolateral components of the ground reaction force (GRF). The results showed that two factors--the number of foot-ground contact points assumed in the model and the type of kinematic constraint enforced at each point--affect the model predictions of muscle coordination. Whereas single- and multiple-point contact models produced similar predictions of muscle function in the sagittal plane, inconsistent results were obtained in the mediolateral direction. Kinematic constraints applied in the sagittal plane altered the model predictions of muscle contributions to the vertical and fore-aft GRFs, while constraints applied in the frontal plane altered the calculations of muscle contributions to the mediolateral GRF. The results illustrate the sensitivity of calculations of muscle coordination to the model used to simulate foot-ground contact.     

In vivo patellofemoral forces in high flexion total knee arthroplasty

This study compares the in vivo patellofemoral contact forces generated in high flexion fixed bearing posterior cruciate retaining Nexgen CR-Flex (PCR) and high flexion posterior stabilized Nexgen LPS-Flex (LPS) TKAs with that of normal knees from full knee extension to maximum weight bearing flexion. Ten patients with the PCR total knee arthroplasty (TKA), ten with the LPS TKA and seven patients having normal knees were fluoroscoped while performing a deep knee bend activity. In vivo femorotibial kinematics, obtained from 3D-to-2D registration technique, and patellar kinematics obtained by direct measurements from the fluoroscopic images were entered into a 3D inverse dynamics mathematical model to determine the in vivo contact forces at the knee. The variation in the patellofemoral and quadriceps forces with flexion were found to be similar across the three groups-increasing from full extension to 90 degrees of flexion, reaching a maximum between 90 degrees and 120 degrees of flexion and then decreasing until maximum flexion. At maximum knee flexion, these forces were found to be significantly lower in the normal knees than in the TKAs. The patellar ligament to quadriceps force ratio decreased with the increase in knee flexion while the patellofemoral to quadriceps force ratio increased. A strong correlation was found to exist between the patellofemoral forces, the femorotibial contact forces and the forces in the extensor mechanism. The PCR TKA in this study exhibited greater resemblance to the normal patients with respect to the patellofemoral forces than the LPS TKA though significant differences in the two implant types were not observed.     

Effects of variation on extensor elements and operative procedures in patellofemoral disorders

The influence of both insertion and strength/elasticity of each extensor in patellofemoral disorders was fully investigated through a two-dimensional mathematical model analysis in a horizontal plane, in combination with experimental design theory for analyzing mutually correlated influences. In the model, patellofemoral joint profiles projected on a horizontal plane have been expressed as spline functions. Each muscle of the quadriceps has been represented as a string pulled by the respective force; fascias and tendons have been represented by springs. Nonlinear equations have been constructed to represent the forces involved, and then solved by numerical iteration. An analysis of variance was performed on the data derived from a series of simulations, obtaining the following results. The strength of most extensors has been shown to have an influence on the increase in lateral contact force but not patellar translation. The tibial tubercle position has significant influence on both patellar translation and lateral contact force. The quadriceps' insertion on the femur has no influence on patellar translation. The insertion of each extensor on the patella has been shown to have a strong effect on patellar translation but not on contact force.     

Physiologically based boundary conditions in finite element modelling

Finite element analysis has been used extensively in the study of bone loading and implant performance, such as in the femur. The boundary conditions applied vary widely, generally producing excessive femoral deformation, and although it has been shown that the muscle forces influence femoral deflections and loading, little consideration has been given to the displacement constraints. It is hypothesised that careful application of physiologically based constraints can produce physiological deformation, and therefore straining, of the femur. Joint contact forces and a complete set of muscle forces were calculated based on the geometry of the Standardised Femur using previously validated musculoskeletal models. Five boundary condition cases were applied to a finite element model of the Standardised Femur: (A) diaphyseally constrained with hip contact and abductor forces; (B) case A plus vasti forces; (C) case A with complete set of muscle forces; (D) distally constrained with all muscle forces; (E) physiological constraints with all muscle forces. It was seen that only the physiological boundary conditions, case E, produced physiological deflections (< 2.0mm) of the femoral head in both the coronal and sagittal planes, which resulted in minimal reaction forces at the constrained nodes. Strains in the mid-diaphysis varied by up to 600 micro-strain under walking loads and 1000 micro-strain under stair climbing loads. The mode of loading, as indicated by the strain profiles on the cortex also varied substantially under these boundary conditions, which has important consequences for studies that examine localised bone loading such as fracture or bone remodelling simulations.     

Contact Profiles in Eight European Countries and Implications for Modelling the Spread of Airborne Infectious Diseases

Background

For understanding the spread of infectious diseases it is crucial to have knowledge of the patterns of contacts in a population during which the infection can be transmitted. Besides contact rates and mixing between age groups, the way individuals distribute their contacts across different locations may play an important role in determining how infections spread through a population.

Methods and Findings

Representative surveys were performed in eight countries to assess the number of social contacts (talking to another person at close distance either with or without physical contact), using a diary approach in which participants recorded individual contacts. The overall sample size was 7290 respondents. We analyzed the reported numbers of contacts per respondent in six different settings (household, work, school, leisure, transportation and others) to define different contact profiles. The identification of the profiles and classification of respondents according to these profiles was conducted using a two-step cluster analysis algorithm as implemented in SPSS.We identified seven distinct contact profiles: respondents having (1) mixed: contacts predominantly at school, during transportation and leisure time, (2) contacts during leisure time, (3) contacts mainly in the household (large family), (4) contacts at work, (5) contacts solely at school, (6) contacts in other places and finally (7) respondents having a low number of contacts in any setting. Similar contact profiles can be found in all eight European countries which participated in the study. The distributions of respondents across the profiles were similar in all countries. The profiles are dominated by work, school and household contacts. But also contacts during leisure activities play an important role in the daily lives of a large fraction of individuals. A surprisingly large number of individuals has only few contacts in all locations. There was a distinct age-dependence in the distribution of the population across contact profiles.

Conclusions

In contrast with earlier studies that focussed on the contribution of different age groups to the spread of an infectious disease, our results open up the opportunity to analyze how an infection spreads between locations and how locations as work or school are interconnected via household contacts. Mathematical models that take these local contact patterns into account can be used to assess the effect of intervention measures like school closure and cancelling of leisure activities on the spread of influenza.     

The effects of tibial component inclination on bone stress after unicompartmental knee arthroplasty

Unlike the case with total knee arthroplasty, the femorotibial angle (FTA) after unicompartmental knee arthroplasty (UKA) does not directly depend on the inclination of the tibial component when the height of the joint line is maintained. This study analyzed the effects of the inclination of the tibial component in the coronal plane on the contact pressure of the implant-bone surface and the stresses on the proximal tibia. A two-dimensional, coronal plane model of the proximal tibia was subjected to finite-element analysis. Sixteen patterns of finite-element models of equal FTA were developed in which the inclination of tibial components ranged from 5 degrees valgus to 10 degrees varus in increments of 1 degrees. Stress concentration at the proximal medial diaphyseal cortex gradually increased as the inclination changed from valgus to varus. Maximum contact pressure on the metal-bone interface similarly changed and shifted from the lateral edge to the medial edge of the implant as the inclination changed to varus. It was found that even without changing FTA, the inclination of the tibial component might affect stress concentration and contact pressure in the proximal tibia after UKA. The results suggested that slight valgus inclination of the tibial component might be preferable to varus and even to 0 degrees (square) inclination so far as the stress distribution is concerned.     

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.     

Are clinical parameters sufficient to model gait patterns in patients with cerebral palsy using a multilinear approach?

The aim of this study was to evaluate whether clinical parameters are sufficient using, a multilinear regression model, to reproduce the sagittal plane joint angles (hip, knee, and ankle) in cerebral palsy gait. A total of 154 patients were included. The two legs were considered (308 observations). Thirty-six clinical parameters were used as regressors (range of motion, muscle strength, and spasticity of the lower). From the clinical gait analysis, the joint angles of the sagittal plane were selected. Results showed that clinical parameter does not provide sufficient information to recover joint angles and/or that the multilinear regression model is not an appropriate solution.     

Finger joint impedance during tapping on a computer keyswitch

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

A model of temporomandibular joint function in anthropoid primates based on condylar movements during mastication

The hypothesis that the shape of the bony temporomandibular joint (TMJ) is functionally related to sagittal sliding of the condyle during mastication is tested, and a model of the relation of sagittal sliding to mandibular size, TMJ shape, and diet is developed. Sagittal sliding is defined as fore-aft motion of the condyle during mandibular translation and/or angular rotation. Ascending ramus height is used as a structural correlate of the distance between the condyle and the mandibular axis of rotation (CR). Cineradiographic data on sagittal sliding and gape during mastication in Ateles spp., Macaca fascicularis, Papio anubis, and Pan troglodytes in conjunction with comparative data on mandibular size and TMJ shape are used to evaluate the hypothesis. The results show that 1) linear and angular gape are highly positively correlated with sagittal sliding, 2) pure mandibular translation is rare during mastication, 3) the CR is rarely if ever located at the condyle during mastication, 4) angular gape should be standardized in interindividual comparisons of sagittal sliding, and 5) the height of the ascending ramus (and by inference the CR-to-condyle distance) is highly positively correlated with absolute sagittal sliding. Sagittal sliding relative to the length of the articular eminence was the variable used to explore the relation between TMJ shape and sliding. This variable standardized absolute sagittal sliding relative to joint size. The relative depth and orientation of the articular eminence were not correlated with relative sagittal sliding. The anteroposterior curvature of the condyle was highly negatively correlated with relative sagittal sliding. Flat condyles are associated with large amounts of relative sagittal sliding. A flat condyle increases joint contact area, which reduces joint stress. A flat condyle also increases joint congruence, and this may facilitate the combined sliding and rolling motion of the condyle when the sliding motion is relatively large. The shape of the entoglenoid process was also positively correlated with relative sagittal sliding. A relatively large entoglenoid process may help to guide sagittal sliding and prevent excessive mediolateral sliding of the condyle. The functional model makes a number of predictions about the correlations between food consistency and food object size, mandibular size, TMJ shape, and sagittal sliding of the condyle during mastication and incision. Am J Phys Anthropol 109:67–88, 1999. © 1999 Wiley-Liss, Inc.