Improved single- and multi-contact life-time testing of dental restorative materials using key characteristics of the human masticatory system and a force/position-controlled robotic dental wear simulator

This paper presents a new in vitro wear simulator based on spatial parallel kinematics and a biologically inspired implicit force/position hybrid controller to replicate chewing movements and dental wear formations on dental components, such as crowns, bridges or a full set of teeth. The human mandible, guided by passive structures such as posterior teeth and the two temporomandibular joints, moves with up to 6 degrees of freedom (DOF) in Cartesian space. The currently available wear simulators lack the ability to perform these chewing movements. In many cases, their lack of sufficient DOF enables them only to replicate the sliding motion of a single occlusal contact point by neglecting rotational movements and the motion along one Cartesian axis. The motion and forces of more than one occlusal contact points cannot accurately be replicated by these instruments. Furthermore, the majority of wear simulators are unable to control simultaneously the main wear-affecting parameters, considering abrasive mechanical wear, which are the occlusal sliding motion and bite forces in the constraint contact phase of the human chewing cycle. It has been shown that such discrepancies between the true in vivo and the simulated in vitro condition influence the outcome and the quality of wear studies. This can be improved by implementing biological features of the human masticatory system such as tooth compliance realized through the passive action of the periodontal ligament and active bite force control realized though the central nervous system using feedback from periodontal preceptors. The simulator described in this paper can be used for single- and multi-occlusal contact testing due to its kinematics and ability to exactly replicate human translational and rotational mandibular movements with up to 6 DOF without neglecting movements along or around the three Cartesian axes. Recorded human mandibular motion and occlusal force data are the reference inputs of the simulator. Experimental studies of wear using this simulator demonstrate that integrating the biological feature of combined force/position hybrid control in dental material testing improves the linearity and reduces the variability of results. In addition, it has been shown that present biaxially operated dental wear simulators are likely to provide misleading results in comparative in vitro/in vivo one-contact studies due to neglecting the occlusal sliding motion in one plane which could introduce an error of up to 49% since occlusal sliding motion D and volumetric wear loss V(loss) are proportional.     

Application and validation of a three-dimensional mathematical model of the human masticatory system in vivo.

A previously described three-dimensional mathematical model of the human masticatory system, predicting maximum possible bite forces in all directions and the recruitment patterns of the masticatory muscles necessary to generate these forces, was validated in in vivo experiments. The morphological input parameters to the model for individual subjects were collected using MRI scanning of the jaw system. Experimental measurements included recording of maximum voluntary bite force (magnitude and direction) and surface EMG from the temporalis and masseter muscles. For bite forces with an angle of 0, 10 and 20 degrees relative to the normal to the occlusal plane the predicted maximum possible bite forces were between 0.9 and 1.2 times the measured ones and the average ratio of measured to predicted maximum bite force was close to unity. The average measured and predicted muscle recruitment patterns showed no striking differences. Nevertheless, some systematic differences, dependent on the bite force direction, were found between the predicted and the measured maximum possible bite forces. In a second series of simulations the influence of the direction of the joint reaction forces on these errors was studied. The results suggest that they were caused primarily by an improper determination of the joint force directions.     

Bite forces and their resultants during forceful intercuspal clenching in humans

We aimed to develop a method of gathering complete information on the system of bite forces acting on the dental arches during clenching with the teeth in maximum intercuspation. Further, we attempted to reduce this system into an equivalent wrench—a force–couple system comprising a single force and a single couple acting along a unique line of action. We investigated the normative distribution of the bite forces and the location and orientation of their resultant wrench in 30 young adults (18–23 yr) with natural dentitions. The number of detected occlusal contacts varied from 12 to 46 (mean: 26.1; SD: 8.4), and was significantly greater for the molars than the premolar and anterior teeth, as were the bite-force magnitudes at individual occlusal contacts (1.2–218.4 N); those resulted in the antero-posteriorly slanted bite-force distribution. The magnitude of the bite-force resultants varied from 246.9 to 2091.9 N, and the points at which the resultant wrench axes intersected the mandibular occlusal plane were located 21.3–37.6 mm posterior to the incisal point and less than 8.9 mm from the midline bilaterally. The bite-force resultant was slightly inclined anteriorly from the perpendicular direction to the mandibular occlusal plane. Our method of using pressure-sensitive films to obtain information on all parameters needed to mechanically define a force (such as magnitude, direction, and point of application) is novel. To our knowledge, this is the first study investigating the system of bite forces during forceful intercuspal clenching in six degrees-of-freedom.     

Factors influencing the anterior component of occlusal force

We hypothesized that the anterior component of occlusal force (ACF) generated by mandibular molars was a function of molar inclination, height of the transverse condylar axis above the occlusal plane, steepness of the occlusal plane, gape, molar root dimensions, interproximal tooth contact force when not biting, and bite force. Our research aim was to identify those biomechanical factors which determine ACF. Mandibular second molars were axially loaded with a 90 N force (10 mm second molar gape) in 15 subjects, and the resulting ACF was measured at the mandibular first molar-second premolar contact using a recording technique based on interproximal frictional forces. Morphologic measurements were obtained from lateral cephalometric radiographs of each subject and included: Frankfort mandibular plane angle, occlusal plane angle, angles formed by the longitudinal axis of the second molar and the occlusal and mandibular planes, perpendicular distance from the top of the condyle to the occlusal plane, and second molar root width and root length. For ten subjects, ACF resulting from axial loads of 50, 100, 150, and 200 N was measured. For ten subjects, ACF resulting from an axial load of 50 N and second molar gapes of 10 mm, 14 mm, 18 mm, and 22 mm were measured. ACF increased with increasing gape and increased proportionally to increasing bite force. Correlation and stepwise regression analyses revealed that ACF varies with interproximal tooth contact force when not biting (contact ‘tightness’) and molar root width (model R² = 0.71, p < 0.01). The hypothesis that ACF is a function of bite force, gape, molar root width, and interproximal contact tightness has been supported, and the hypothesis that ACF is a function of molar inclination, occlusal plane steepness, condylar axis height, and root length was rejected.     

A novel computational method to determine subject-specific bite force and occlusal loading during mastication

The evaluation of three-dimensional occlusal loading during biting and chewing may assist in development of new dental materials, in designing effective and long-lasting restorations such as crowns and bridges, and for evaluating functional performance of prosthodontic components such as dental and/or maxillofacial implants. At present, little is known about the dynamic force and pressure distributions at the occlusal surface during mastication, as these quantities cannot be measured directly. The aim of this study was to evaluate subject-specific occlusal loading forces during mastication using accurate jaw motion measurements. Motion data was obtained from experiments in which an individual performed maximal effort dynamic chewing cycles on a rubber sample with known mechanical properties. A finite element model simulation of one recorded chewing cycle was then performed to evaluate the deformation of the rubber. This was achieved by imposing the measured jaw motions on a three-dimensional geometric surface model of the subject’s dental impressions. Based on the rubber’s deformation and its material behaviour, the simulation was used to compute the resulting stresses within the rubber as well as the contact pressures and forces on the occlusal surfaces. An advantage of this novel modelling approach is that dynamic occlusal pressure maps and biting forces may be predicted with high accuracy and resolution at each time step throughout the chewing cycle. Depending on the motion capture technique and the speed of simulation, the methodology may be automated in such a way that it can be performed chair-side. The present study demonstrates a novel modelling methodology for evaluating dynamic occlusal loading during biting or chewing.     

A three-dimensional mathematical model of the human masticatory system predicting maximum possible bite forces

A three-dimensional mathematical model of the human masticatory system, containing 16 muscle forces and two joint reaction forces, is described. The model allows simulation of static bite forces and concomitant joint reaction forces for various bite point locations and mandibular positions. The system parameters for the model were obtained from a cadaver head. Maximum possible bite forces were computed using optimization techniques; the optimization criterion we used was the minimizing of the relative activity of the most active muscle. The model predicts that at each specific bite point, bite forces can be generated in a wide range of directions, and that the magnitude of the maximum bite force depends on its direction. The relationship between bite force direction and its maximum magnitude depends on bite point location and mandibular position. In general, the direction of the largest possible bite force does not coincide with the direction perpendicular to the occlusal plane.     

Three-dimensional finite element analysis of glass-ceramic dental crowns

Because of the improved esthetic potential of glass-ceramic crowns as dental restorations, they are sometimes preferred over metal-ceramic crowns for restoration of anterior teeth. Because of their relatively high strength, these ceramic crowns are also frequently used for restoration of posterior teeth. However, due to the larger magnitude of biting forces on posterior teeth, intraoral fracture of all-ceramic crowns tends to occur more frequently in posterior crowns (Moffa, 1988). The objective of this study was to determine the relative influence of load orientation and the occlusal thickness of posterior ceramic crowns on the stress distribution which develops under these loading and design conditions. Three-dimensional finite element models for a molar crown were developed to determine the stress distribution under simulated applied loads. Glass-ceramic crowns with occlusal thicknesses of 0.5, 1.5, and 3.0 mm were considered. The largest principal tensile stresses induced in ceramic due to a distributed load of 600 N applied in a cuspal region were approximately 12 and 182 MPa for vertical and horizontal loading orientations, respectively. Stresses which developed in the facial and lingual marginal regions were primarily compressive under vertical loads. However, tensile stresses developed when the load was applied horizontally. Differences in stress distribution within crowns with the three occlusal thicknesses occurred only near the site of loading. Because of the relatively large failure rates of ceramic crowns in the posterior regions, these restorations should be strengthened by improvement in design, composition, and thermal processing conditions. Before any significant progress is made in these areas, these restorations should be used for the anterior teeth. The results of this study suggest that orientation of the applied load has a more important effect on development of large tensile stresses than the occlusal thickness of ceramic.     

Efficiency as a predictor of human jaw design in the sagittal plane

The effects of changing the direction of the bite force and of the mandibular joint reaction have been studied with a mathematical model assisted by a computer using the technique of linear programming. We conclude the following: In the sagittal plane the long axes of lower molars are each tilted in the direction that most efficiently converts muscle force into work at the bite point rather than in the direction that would maximize static bite force. These genetically determined angles are referred to as the most 'work efficient' angles. Collectively they lead to the appearance of the curve of Spee associated with the postcanines. Given the most work efficient angle of the first molar, the model indicates for bite forces generated in this direction the joint reaction is least when tilted forward from the vertical at between 20 degrees and 30 degrees. The joint reaction is normal to the articular surface of the condyle which is itself tilted forward 20-30 degrees from the occlusal plane. We conclude the condyle and articular eminence are remodelled to the angle that minimizes the joint reaction. The direction of the bite force may be controlled via neuronal circuitry connecting mechanoreceptors of the periodontal ligament with motor nerves supplying the jaw-closing muscles. The height of the occlusal plane in the molar region has little effect on jaw efficiency.     

Brief communication: comparing loading scenarios in lower first molar supporting bone structure using 3D finite element analysis

Finite element analysis (FEA) is a widespread technique to evaluate the stress/strain distributions in teeth or dental supporting tissues. However, in most studies occlusal forces are usually simplified using a single vector (i.e., point load) either parallel to the long tooth axis or oblique to this axis. In this pilot study we show how lower first molar occlusal information can be used to investigate the stress distribution with 3D FEA in the supporting bone structure. The LM₁ and the LP²‐LM¹ of a dried modern human skull were scanned by μCT in maximum intercuspation contact. A kinematic analysis of the surface contacts between LM₁ and LP²‐LM¹ during the power stroke was carried out in the occlusal fingerprint analyzer (OFA) software to visualize contact areas during maximum intercuspation contact. This information was used for setting the occlusal molar loading to evaluate the stress distribution in the supporting bone structure using FEA. The output was compared to that obtained when a point force parallel to the long axis of the tooth was loaded in the occlusal basin. For the point load case, our results indicate that the buccal and lingual cortical plates do not experience notable stresses. However, when the occlusal contact areas are considered, the disto‐lingual superior third of the mandible experiences high tensile stresses, while the medio‐lingual cortical bone is subjected to high compressive stresses. Developing a more realistic loading scenario leads to better models to understand the relationship between masticatory function and mandibular shape and structures. Am J Phys Anthropol, 2012. © 2011 Wiley Periodicals, Inc.     

A three-dimensional finite element study on the biomechanical behavior of an FGBM dental implant in surrounding bone

Dental implants made of functionally graded biomaterials (FGBM) have been receiving increasing attention due to their unique advantage of being able to simultaneously satisfy biocompatibility, strength, corrosion resistance, etc., which a single composition with a uniform structure cannot satisfy. This paper investigates the biomechanical behavior of a threaded FGBM dental implant/surrounding bone system under static and harmonic occlusal forces by using a three-dimensional finite element method. The implant is a mixture of a bioceramic and a biometal with a smooth gradient in both the material composition and properties in the longitudinal direction. The interaction of the implant and the supporting bone tissues is considered. Three contact conditions at the implant-bone interface are used to model different osseointegration stages. A comprehensive parametric study is conducted to highlight the influence of the material properties, the volume fraction index, the occlusal force orientation, and the osseointegration quality on the maximum von-Mises stress, deformation distribution, natural frequencies, and harmonic response.     

The role of cutaneous feedback for anticipatory grip force adjustments during object movements and externally imposed variation of the direction of gravity

Grip force adjustments to changes of object loading induced by external changes of the direction of gravity during discrete arm movements with a grasped object were analyzed during normal and anesthetized finger sensibility. Two subjects were seated upright in a rotatable chair and rotated backwards into a horizontal position during discrete movements with a hand-held instrumented object. The movement direction varied from vertical to horizontal inducing corresponding changes in the direction of gravity, but the orientation of the movement in relation to the body remained unaffected. During discrete vertical movements a maximum of load force occurs early in upward and late in downward movements; during horizontal movements two load force peaks result from both acceleratory and deceleratory phases of the movement. During performance with normal finger sensibility grip force was modulated in parallel with fluctuations of load force during vertical and horizontal movements. The grip force profile adopted to the varying load force profile during the transition from the vertical to the horizontal position. The maximum grip force occurred at the same time of maximum load force irrespective of the movement plane. During both subjects' first experience of digital anesthesia the object slipped from the grasp during rotation to the horizontal plane. During the following trials with anesthetized fingers subjects substantially increased their grip forces, resulting in elevated force ratios between maximum grip and load force. However, grip force was still modulated with the movement-induced load fluctuations and maximum grip force coincided with maximum load force during vertical and horizontal movements. This implies that the elevated force ratio between maximum grip and load force does not alter the feedforward system of grip force control. Cutaneous afferent information from the grasping digits seems to be important for the economic scaling of the grip force magnitude according to the actual loading conditions and for reactive grip force adjustments in response to load perturbations. However, it plays a subordinate role for the precise anticipatory temporal coupling between grip and load forces during voluntary object manipulation.     

Retention of telescopic denture in elderly patients with maximum partially edentulous arch

The telescopic denture is one possible solution for prosthetic rehabilitation for a partially dentate arch with few remaining teeth, which is the reason why it is usually considered in elderly patients. The system of double crowns consists of the male component of attachment or patrix, which is cemented to the abutment tooth, and the female component of attachment or matrix, which is the removable part of the restoration. Retention force represents resistance to those forces which tend to separate the partial denture and supporting tissue in the occlusal direction. Its optimal value is 5–9 N per each patrix and it should strictly be between these limits as a stronger or weaker force will endanger the functionality of the denture or damage the supporting tissue of the abutment teeth. Objective: The goal of this study was to establish the exact value of the individual and total retention force of the completed telescopic system. Method: The Bredent dynamometer was used for measuring these values. The study involved measuring the individual force of separation of 50 patrix–matrix components, and total force of 20 partial dentures of various designs retained by double crowns made in the Dental Laboratory of the Clinic for Prosthetic Dentistry, Faculty of Stomatology, Belgrade University. The research was carried out in the period from January 2006 to October 2006. Results and conclusion: The results showed different values of separation force – from very low values to those exceeding optimal limits. A reduction in the force of separation is necessary in the case of excessive tightness. If there was looseness between patrix and matrix, the cylindrical components of double crowns need to be lined by a suitable composite material.     

The role of cutaneous feedback for anticipatory grip force adjustments during object movements and externally imposed variation of the direction of gravity

Grip force adjustments to changes of object loading induced by external changes of the direction of gravity during discrete arm movements with a grasped object were analyzed during normal and anesthetized finger sensibility. Two subjects were seated upright in a rotatable chair and rotated backwards into a horizontal position during discrete movements with a hand-held instrumented object. The movement direction varied from vertical to horizontal inducing corresponding changes in the direction of gravity, but the orientation of the movement in relation to the body remained unaffected. During discrete vertical movements a maximum of load force occurs early in upward and late in downward movements; during horizontal movements two load force peaks result from both acceleratory and deceleratory phases of the movement. During performance with normal finger sensibility grip force was modulated in parallel with fluctuations of load force during vertical and horizontal movements. The grip force profile adopted to the varying load force profile during the transition from the vertical to the horizontal position. The maximum grip force occurred at the same time of maximum load force irrespective of the movement plane. During both subjects' first experience of digital anesthesia the object slipped from the grasp during rotation to the horizontal plane. During the following trials with anesthetized fingers subjects substantially increased their grip forces, resulting in elevated force ratios between maximum grip and load force. However, grip force was still modulated with the movement-induced load fluctuations and maximum grip force coincided with maximum load force during vertical and horizontal movements. This implies that the elevated force ratio between maximum grip and load force does not alter the feedforward system of grip force control. Cutaneous afferent information from the grasping digits seems to be important for the economic scaling of the grip force magnitude according to the actual loading conditions and for reactive grip force adjustments in response to load perturbations. However, it plays a subordinate role for the precise anticipatory temporal coupling between grip and load forces during voluntary object manipulation.     

Friction and adhesion in the tarsal and metatarsal scopulae of spiders

Friction and adhesion forces of the ventral surface of tarsi and metatarsi were measured in the bird spider Aphonopelma seemanni (Theraphosidae) and the hunting spider Cupiennius salei (Ctenidae). Adhesion measurements revealed no detectable attractive forces when the ventral surfaces of the leg segments were loaded and unloaded against the flat smooth glass surface. Strong friction anisotropy was observed: friction was considerably higher during sliding in the distal direction. Such anisotropy is explained by an anisotropic arrangement of microtrichia on setae: only the setal surface facing in the distal direction of the leg is covered by the microtrichia with spatula-like tips. When the leg is pushed, the spatula-shaped tips of microtrichia contact the substrate, whereas, when the leg is pulled over a surface, setae bend in the opposite direction and contact the substrate with their spatulae-lacking sides. In an additional series of experiments, it was shown that desiccation has an effect on the friction force. Presumably, drying of the legs results in reduction of the flexibility of the setae, microtrichia, spatulae, and underlying cuticle; this diminishes the ability to establish proper contact with the substrate and thus reduces the contact forces.     

Determination of muscle loading at the hip joint for use in pre-clinical testing

The stability of joint endoprostheses depends on the loading conditions to which the implant-bone complex is exposed. Due to a lack of appropriate muscle force data, less complex loading conditions tend to be considered in vitro. The goal of this study was to develop a load profile that better simulates the in vivo loading conditions of a "typical" total hip replacement patient and considers the interdependence of muscle and joint forces. The development of the load profile was based on a computer model of the lower extremities that has been validated against in vivo data. This model was simplified by grouping functionally similar hip muscles. Muscle and joint contact forces were computed for an average data set of up to four patients throughout walking and stair climbing. The calculated hip contact forces were compared to the average of the in vivo measured forces. The final derived load profile included the forces of up to four muscles at the instances of maximum in vivo hip joint loading during both walking and stair climbing. The hip contact forces differed by less than 10% from the peak in vivo value for a "typical" patient. The derived load profile presented here is the first that is based on validated musculoskeletal analyses and seems achievable in an in vitro test set-up. It should therefore form the basis for further standardisation of pre-clinical testing by providing a more realistic approximation of physiological loading conditions.     

Contact fracture of full-ceramic crowns subjected to occlusal loads

The mechanisms contributing to failure of full dental ceramic crowns under occlusal loads were studied using a unique optical approach. Model specimens comprising triple-layered crowns (veneer, core and substrate) were developed with both flat and curved occlusal surfaces and then subjected to simulated quasi-static occlusal loading using a spherical indenter. Deformation within the specimens during loading was analyzed by means of digital image correlation (DIC). Finite element models were also developed and used to examine the mechanics of contact. Results of the experiments with flat dental crowns indicated three typical modes of failure, i.e. cone cracks, plastic yielding and radial cracks. Fracture of the specimens with curved dental crowns was complicated by contributions from competing and multiple modes of failure. Both experimental and numerical results conclude that the dominant fracture mode in the full-ceramic crowns was radial cracking in the core beneath the contact area. However, displacement fields obtained using DIC showed that debonding developed near the shoulder of the crown, particularly during off-axis loading, and initiated under substantially lower occlusal loads than those required for crack initiation.     

Pattern of anterior cruciate ligament force in normal walking

The goal of this study was to calculate and explain the pattern of anterior cruciate ligament (ACL) loading during normal level walking. Knee-ligament forces were obtained by a two-step procedure. First, a three-dimensional (3D) model of the whole body was used together with dynamic optimization theory to calculate body-segmental motions, ground reaction forces, and leg-muscle forces for one cycle of gait. Joint angles, ground reaction forces, and muscle forces obtained from the gait simulation were then input into a musculoskeletal model of the lower limb that incorporated a 3D model of the knee. The relative positions of the femur, tibia, and patella and the forces induced in the knee ligaments were found by solving a static equilibrium problem at each instant during the simulated gait cycle. The model simulation predicted that the ACL bears load throughout stance. Peak force in the ACL (303 N) occurred at the beginning of single-leg stance (i.e., contralateral toe off). The pattern of ACL force was explained by the shear forces acting at the knee. The balance of muscle forces, ground reaction forces, and joint contact forces applied to the leg determined the magnitude and direction of the total shear force acting at the knee. The ACL was loaded whenever the total shear force pointed anteriorly. In early stance, the anterior shear force from the patellar tendon dominated the total shear force applied to the leg, and so maximum force was transmitted to the ACL at this time. ACL force was small in late stance because the anterior shear forces supplied by the patellar tendon, gastrocnemius, and tibiofemoral contact were nearly balanced by the posterior component of the ground reaction.     

实验性牙周炎合并咬合创伤模型的研究进展

不正常的咬合接触关系或者过大的咬合力,造成咀嚼系统各部位的病理性损害或者适应性变化称为咬合创伤,而牙周炎是一种慢性炎症性疾病,是由牙菌斑生物膜引起的牙周组织的感染性疾病,导致牙周支持组织的破坏-牙周袋形成和炎症、进行性的附着丧失和牙槽骨吸收。近年来有学者对咬合创伤作为一类独立疾病进行研究,也有学者将咬合创伤作为牙周炎的一个重要的局部促进因素进行研究。牙周炎和咬合创伤的关系至今不是很明确,关于咬合创伤和牙周炎的相互关系,虽然早在20世纪初已开始研究,却直到70年代才有严格对照的动物实验研究。大鼠牙周炎合并咬合创伤的模型在牙周炎和咬合创伤关系的相关研究中起着十分重要的作用。本文从牙周炎和咬合创伤复合模型的建立方面,综述了各种建模方法的研究进展,以期为牙周炎合并咬合创伤的深入研究提供参考。     

Characteristics of human fingertips in the shearing direction

The shearing strain of the human fingertip plays an important role in the determination of the optimal grasping force and in the perception of texture. Most research concerned with the mechanical impedance of the human fingertips has treated the orthogonal direction to the tip surface, and little attention has been paid to the tangential direction. This paper describes impedance characteristics of the human fingertips in the tangential directions to the tip surface. In the experiment, step and ramp shearing forces were individually applied to the tips of the thumb, middle finger, and little finger. Dynamics of the fingertips were represented by the Kelvin model. Experimental results show that each fingertip had different properties with respect to the shearing strain versus the applied force, and that the thumb had the strongest shearing stiffness among these three digits. Moreover, the shearing stiffness depended on the direction of the applied force, and the stiffness in the pointing direction was stronger than that in the perpendicular direction. As the contact force in the orthogonal direction to the fingertip surface was increased, the shearing stiffness and viscosity increased without regard to the load speed of the shearing force. Furthermore, it is shown that the average strain rate of the fingertip in the tangential direction to the fingertip surface became slower and converged to a constant value with higher contact forces.     

The influence of patellofemoral joint contact geometry on the modeling of three dimensional patellofemoral joint forces

The purpose of this study was to determine the influence of patellofemoral joint contact geometry on the modeling of three-dimensional patellofemoral joint forces. To achieve this goal, patellofemoral joint reaction forces (PFJRFs) that were measured from an in-vitro cadaveric set-up were compared to PFJRFs estimated from a computer model that did not consider patellofemoral joint contact geometry. Ten cadaver knees were used in this study. Each was mounted on a custom jig that was fixed to an Instron frame. Quadriceps muscle loads were accomplished using a pulley system and weights. The force in the patellar ligament was obtained using a buckle transducer. To quantify the magnitude and direction of the PFJRF, a six-axis load cell was incorporated into the femoral fixation system so that a rigid body assumption could be made. PFJRF data were obtained at 0 degrees , 20 degrees , 40 degrees and 60 degrees of knee flexion. Following in vitro testing, SIMM modeling software was used to develop computational models based on the three-dimensional coordinates (Microscribe digitizer) of individual muscle and patellar ligament force vectors obtained from the cadaver knees. The overall magnitude of the PFJRF estimated from the computer generated models closely matched the direct measurements from the in vitro set-up (Pearson's correlation coefficient, R(2)=0.91, p<0.001). Although the computational model accurately estimated the posteriorly directed forces acting on the joint, some discrepancies were noted in the forces acting in the superior and lateral directions. These differences however, were relatively small when expressed as a total of the overall PFJRF magnitude.