Abstract:
A method is disclosed for determining and predicting temporomandibular joint (TMJ) motion and to use the predictive model to prepare improved TMJ replacement joints, improved crown, cap or bridge contouring, improved teeth alignments and placements, and improved mouth guard and dentures.

Description:
RELATED APPLICATIONS 
   This application is a 35 U.S.C.§371 nationalization of PCT/US05/04465 filed 14 Feb. 2005, which claims priority to United States Provisional Patent Application Ser. No. 60/544,725 filed 13 Feb. 2004. 

   BACKGROUND OF THE INVENTION 
   1. Field of the Invention 
   The present invention relates to a system and a method for determining a guide plane of a temporomandibular joint (TMJ) in an animal including a human and using the guide plane to design dental prostheses, to aid in dental reconstruction surgery, to aid in orthodontic configurations and final teeth placement, to aid in cap design and contouring and to aid in the design of mouth guards to protect and properly align the jaw. 
   More particularly, the present invention relates to a system and a method for determining a guide plane of an TMJ in an animal including a human and using the guide plane to design dental prostheses, to aid in dental reconstruction surgery, to aid in orthodontic configurations and final teeth placement, to aid in cap design and contouring and to the design of mouth guards to protect and properly align the jaw, where the system includes obtaining digital images at a fixed focal length of a lower and upper jaw, scanning the digital image into 3-D images, transporting the 3-D images to a 3-D software program, using a derived mathematical function and determining contact planes of functional teeth surfaces (wear facets) and a glenoid fossa, Zola&#39;s tubercle. The invention also relates to methods for using the determined teeth data to design prosthesis, to design caps, to design mouth guards, to improve orthodontic configuration to produce improved teeth alignment, to augment teeth placement to relieve stress on the TM joint. 
   2. Description of the Related Art 
   Analysis of human skulls is a delicate and difficult process. Much of the significant work that has been done in this area of study, morphometrics, has been conducted using analog techniques, i.e., direct measurements, using instrument such as calipers, dividers, rulers, grids and/or custom designed “guides”. 
   To better understand the functional relationship of dental occlusions and the bony morphology of the skull and mandible, direct measurement and/or analog models were insufficient at supplying necessary mathematical formula to develop a predictive model evidencing the relationship. 
   After determining that classical analysis using analog models and protocols derived thereform did not yield the desired mathematical described relationship, the inventors sought a more detailed measurements analytical approach. First, the inventors studied radiographic techniques and attempted correlation the radiographic data to photographic evidence. However, the distortional aspect of the radiograph data rendered such techniques useless. The inventors also attempted using 3-D, corrected CT Scan data, but this technique proved difficult to manage 
   Thus, there is a need in the art for an improved method for determining the motion of the temporomandibular joint allowing for improved dental, medical and rehabilitative prostheses, mouth guards, teeth alignment, engaging surface contouring and other similar methods. 
   SUMMARY OF THE INVENTION 
   The present invention provides a method for analyzing jaw motion including the steps of determining planes associated with wear surfaces of at least two teeth on a lower jaw, at least one of the teeth located on an ipsilateral posterior side of the jaw and at least another one of the teeth located on a contralateral anterior side of the jaw, calculating normals vectors to the two planes directed from the planes towards the upper jaw, determining position vectors between a gliding temporomandibular joint (TMJ) pivot point and a center point of each of the planes and a rotating temporomandibular joint (TMJ) pivot point. 
   The present invention also provides a method implemented on a computer for predicting motion of a temporomandibular joint including the steps of generating a 3-D image of an upper jaw, a lower jaw, upper and lower teeth and the temporomandibular joint (TMJ). After the image is created, a right or left TMJ pivot point is located and planar wear surfaces are determined for at least two teeth on the lower and upper jaw, at least one of the teeth located on an ipsilateral posterior side of the jaw and at least another one of the teeth located on a contralateral anterior side of the jaw. Once the pivot point and the planes are located, normal vectors to the planes directed towards the upper jaw are determined and pointing vectors are determine between the pivot point and an opposite side TMJ pivot and between the pivot point and the two normal vectors. With the vectors determined, creating a TMJ trajectory given by equation (1): 
                   {           P     TMJ   ⁢           ⁢   x     ′               P     TMJ   ⁢           ⁢   y     ′               P     TMJ   ⁢           ⁢   z     ′           }     =       [           ⁢             cos   ⁢           ⁢   Φ     +       E     R   ⁢           ⁢   x     2     ⁡     (     1   -     cos   ⁢           ⁢   Φ       )                   E     R   ⁢           ⁢   x       ⁢       E     R   ⁢           ⁢   y       ⁡     (     1   -     cos   ⁢           ⁢   Φ       )         +       E     R   ⁢           ⁢   z       ⁢   sin   ⁢           ⁢   Φ                 E     R   ⁢           ⁢   x       ⁢       E     R   ⁢           ⁢   z       ⁡     (     1   -     cos   ⁢           ⁢   Φ       )         -       E     R   ⁢           ⁢   y       ⁢   sin   ⁢           ⁢   Φ                     E     R   ⁢           ⁢   x       ⁢       E     R   ⁢           ⁢   y       ⁡     (     1   -     cos   ⁢           ⁢   Φ       )         -       E     R   ⁢           ⁢   z       ⁢   sin   ⁢           ⁢   Φ               cos   ⁢           ⁢   Φ     +       E     R   ⁢           ⁢   y     2     ⁡     (     1   -     cos   ⁢           ⁢   Φ       )                   E     R   ⁢           ⁢   y       ⁢       E     R   ⁢           ⁢   z       ⁡     (     1   -     cos   ⁢           ⁢   Φ       )         +       E     R   ⁢           ⁢   x       ⁢   sin   ⁢           ⁢   Φ                     E     R   ⁢           ⁢   x       ⁢       E     R   ⁢           ⁢   z       ⁡     (     1   -     cos   ⁢           ⁢   Φ       )         +       E     R   ⁢           ⁢   y       ⁢   sin   ⁢           ⁢   Φ                 E     R   ⁢           ⁢   y       ⁢       E     R   ⁢           ⁢   z       ⁡     (     1   -     cos   ⁢           ⁢   Φ       )         -       E     R   ⁢           ⁢   x       ⁢   sin   ⁢           ⁢   Φ               cos   ⁢           ⁢   Φ     +       E     R   ⁢           ⁢   z     2     ⁡     (     1   -     cos   ⁢           ⁢   Φ       )               ⁢           ]     ⁢     {           P     TMJ   ⁢           ⁢   x                 P     TMJ   ⁢           ⁢   y                 P     TMJ   ⁢           ⁢   z             }               (   1   )               
where P′ TMJx , P′ TMJy , and P′ TMJz  are components of an opposite TMJ trajectory vector P′ TMJ , Φ is any rotation about the TMJ pivot point, P TMJx , P TMJy , and P TMJz  are components of a TMJ position vector P TMJ , and a 3×3 rotation matrix, where E rx , E ry  and E Rz  are components of a vector E R  given by the formula:
 
                     N   -&gt;     P1     =         P   -&gt;     1     ×       N   -&gt;     1               (   2   )                   N   -&gt;     P2     =         P   -&gt;     2     ×       N   -&gt;     2               (   3   )                   E   -&gt;     R     =           N   -&gt;     P1     ×       N   -&gt;     P2                  N   -&gt;     P1     ×       N   -&gt;     P2                      (   4   )               
where N 1  and N 2  are normal vectors to the two planes and P 1  and P 2  are the position vectors from the TMJ pivot point to the normals and N P1  and N P2  are vector cross products as defined in equations 2 and 3. The TMJ trajectory permits a medical professional to design improved replacement TMJs from the wear surfaces of two lower or upper jaw teeth. The TMJ trajectory permits improved contouring of crowns, caps and bridges. The TMJ trajectory permits improved teeth placement augmentations to improve bite, TMJ functionality and to maximize TMJ motion and to reduce or correct teeth induced stress on the TMJ or on TMJ motion.
 
   The present invention also provides improved TMJ replacement joints designed using the formula of equation (1). 
   The present invention also provides improved crowns, caps or bridges designed using the formula of equation (1). 
   The present invention also provides improved mouth guards designed using the formula of equation (1). 
   The present invention also provides improved teeth alignments using the formula of equation (1). 
   The present invention also provides improved denture designs using the formula of equation (1). 

   
     DESCRIPTION OF THE DRAWINGS 
     The invention can be better understood with reference to the following detailed description together with the appended illustrative drawings in which like elements are numbered the same: 
       FIG. 1  depicts a photograph of a Sony Mavica MVC-CD1OOO 000 digital camera set a fixed focal length from a skull; 
       FIG. 2  depicts the image being digitally scanned using a Minolta Vivid 3-D laser scanner; 
       FIG. 3  depicts importation of the digitized image into a Dual Processor Dell workstation and assembled into 3-D image processing software; 
       FIG. 4  depicts the scanned surfaces selected for analysis from the skull image; 
       FIG. 5  depicts the digital composite of the skull and the vectors required by the invention to determine and TMJ trajectory equation; 
       FIG. 6  depicts specific areas of interest of a medial wall of a glenoid fossa on a moving or gliding side, a posterior tooth contact on an ipsilateral side and an anterior tooth on an opposite or contralateral side; 
       FIG. 7  depicts spherically oriented locations of selected planes of contact; 
       FIG. 8  depicts areas selected for analysis by the image processing software; 
       FIG. 9  depicts a plot of a mathematical model showing excellent correlation; 
       FIG. 10  depicts a diagram of the TMJ pivot and the two wear surface planes of interest; 
       FIG. 11  depicts a diagram of the position, normal and TMJ vectors; 
       FIG. 12  depicts a diagram illustrating the jaw rigid body motion; 
       FIG. 13  depicts that the rotation vector must be in the planes defined by the position and normal vectors to the wear surfaces; 
       FIGS. 14A&amp;B  depict opposite TMJ trajectory constrained to the surface of a sphere centered at the opposite TMJ pivot point with a radius of |P TMJ |; 
       FIGS. 15A&amp;B  depict test case study for the A4 skull, TMJ and tooth locations, and TMJ trajectories; 
       FIGS. 16A&amp;B  depict an actual case study for the A4 skull, TMJ and tooth locations, and TMJ trajectories; 
       FIGS. 17A&amp;B  depict an actual case study for the A10 skull, TMJ and tooth locations, and TMJ trajectories; 
       FIGS. 18A&amp;B  depict an actual case study for the A12 skull, TMJ and tooth locations, and TMJ trajectories; 
       FIGS. 19A&amp;B  depict an actual case study for the A14 skull, TMJ and tooth locations, and TMJ trajectories; 
       FIGS. 20A&amp;B  depict an actual case study for the A skull, TMJ and tooth locations, and TMJ trajectories; 
       FIGS. 21A&amp;B  depict an actual case study for the Anat skull, TMJ and tooth locations, and TMJ trajectories; 
       FIGS. 22A&amp;B  depict an actual case study for the C13 skull, TMJ and tooth locations, and TMJ trajectories; 
       FIG. 23  depicts that the TMJ trajectory normal vector is constrained to the plane of the circular trajectory path; 
       FIG. 24  depicts that multiple possibilities of surface orientations exist between the TMJ and the mandibular cavity; and 
       FIGS. 25A&amp;B  depict mandibular cavity normal vector correlation for skulls A4, A10, A12, A14, A, Anat and C13 in tabulated form. 
   

   DETAILED DESCRIPTION OF THE INVENTION 
   The inventors have found that a mathematical model can be constructed from wear patterns on at least two teeth on opposite side of the jaw to determine and understand temporomandibular joint (TMJ) motion. The determination and understanding of TMJ motion allows the inventors to construct improved dental prostheses such as artificial TMJs, improved caps and/or replacement teeth, or the like and to design improved teeth alignment for orthodontic configurations and improved tooth augmentation to improve or change TMJ dynamics. 
   The present invention broadly relates to a predictive mathematical formula based on morphometric analysis, analog and/or digital, of a TMJ and teeth of an animal including a human, where the mathematical formula is capable of describing or predicting complex TMJ motion patterns as evidenced by tooth wear and/or jaw wear. Data derived from the formula can then be used to construct remedial appliances or prostheses for the head and neck, and for dental, medical and/or physical therapeutic applications. 
   The present invention also broadly relates to a method including the step of determining a set of vectors derived from planes of wear of at least two teeth on opposite sides of a jaw and a TMJ pivot point, supplying the vectors to the predictive mathematical formula, deriving a TMJ trajectory model and designing remedial appliances or prostheses for the head and neck, and for dental, medical and/or physical therapeutic applications from the TMJ trajectory model. 
   Laser Scan and Data Analysis 
   The inventors investigated the use of a laser scanner technique using a Minolta Laser Scanner made available for initial studies at the research laboratories of Dr. John Kappelman at the University of Texas at Austin, Department of Physical Anthropology. A skull sample was selected from the Department&#39;s collection and a trial laser scan was made to test the applicability of the apparatus and to verify that this technique could provide data sufficient to determine the desired relationships and ultimately reduce the relationships to a universal mathematical formula that can be used to predict TMJ motion. From data derived from this formula, improved dental, medical and physical therapeutic treatments of the jaw and teeth can be designed as well as improved prostheses, cap, teeth alignment, teeth resurfacing or other corrective dental procedures. 
   The sample was scanned using a Minolta Vivid 3-D scanner. The scan was transferred to a digital processing unit, a Dell dual processor work station and was imported into a digital analysis software program called Rapidform available from INUS Technology, Inc. for data analysis and direct measurement determinations. 
   The trial scan was used to revise the basic thesis and protocols of this invention. We performed a random series of scans of relevant dry specimens in a laboratory facility. The selected specimens came from a private collection of 7 skulls. The seven skulls had been acquired over a period of time and were not collected specifically for this study. Thus, the selection process was random. In addition to these seven skulls, a random selection of skulls from the University of Texas Dental Branch at Houston was also selected. 
   The selected skulls had to have a fall compliment of teeth and were selected in three categories: (a) youth (deciduous or mixed teeth); (b) adolescent (all permanent teeth with little wear); and (c) adult (permanent teeth in varying stages of wear). Ethnicity was not considered (or known) in the selection process. A group of 24 skulls was selected that fit the basic criteria, grouped as to age and then a non-professional staff member was asked to pick any 6 skulls from the 24 skulls, while the skulls were all closed in their skull cases, i.e., the skulls were not visible. 
   With skulls in hand, each skull was photographed at a standardized, fixed, focal length and then scanned from the inferior surface using the Minolta Vivid 3-D scanner. The scanning process surveyed the base of the skull, the glenoid fossa, the occipital tubercle, the surfaces of the teeth and the hard palate. No mandibles were scanned. 
   Multiple scans were made of each skull to obtain the surfaces required to complete the study. These scans were assembled into single composite 3-D images which could be manipulated as needed. The composite 3-D images were compared to the physical specimens for accuracy of detail and measurement. 
   These images were then used, with the appropriate software, to make linear, angular, and complex spherical studies of multiple surfaces in the 3-D composite images. Of particular interest in the early analyses was to find a mathematical solution for correlated angular function of the wear surfaces of teeth (facets) and the guide (glide) angle of the medial aspect of the glenoid fossa, Zola&#39;s Tubercle. See Zola A and Rothschild E. A.,  J. Pros Dent ., “Posterior Condyle Positions in Unimpeded Jaw Movements,” 1962, September-October, Vol. II, No. 5, pp 873-881. 
   The use of non-linear mathematic analysis for solving for angular relationships with 6 degrees of freedom has demonstrated a positive relationship to these areas shaped by dental occlusal function. This work is original; no comparable work is seen in the research literature. 
   Protocol 
   Skulls used in the study were photographed at a fixed focal length using a Sony Mavica MVC-CD1000 digital camera as shown in  FIG. 1  and then digitally scanned using a Minolta Vivid 3-D laser scanner as shown in  FIG. 2 . The scans were imported into a Dual Processor Dell workstation and assembled into manipulatable 3-D images using INUS Technology, Inc., Rapidform software as shown in  FIG. 3 . 
   Specific areas of interest on scanned surfaces were selected for analysis as shown in FIG.  4 —skull and as shown in FIG.  5 —digital composite. The areas identified had been previously marked in dynamic function using an inked red silk ribbon. The specific areas of interest were a medial wall of a glenoid fossa on a moving or gliding side, a posterior tooth contact on an ipsilateral side and an anterior tooth on an opposite or contralateral side as shown in  FIG. 6 . 
   Areas of coordinated functional contacts were mapped by the imaging program which, then, provided detailed mapping coordinates of these specific areas as shown in FIG.  7 —data chart. The data were analyzed using linear and non-linear geometric techniques for complex mathematical analyses. The first goal was to make a mathematical model that would transpose the scanned data in a computer model to plot the selected points of function and then to demonstrate the accuracy of these plots as shown in FIGS.  8  and  9 —skull and plot photo&#39;s. This series of analyses showed the greater majority of the digital data correlated with the scanned data and the physical model of the skull. Where discrepancies were found, we assumed that the selected scan was inconsistent with the anatomical model and reassessment of the scanned images corrected the discrepancy. The study continued with more than 80% of the data being accurate. 
   In conclusion, this series of studies substantially proved that the digital scans provide sufficient data to locate contact planes of functioning surfaces on teeth (wear facets) and the glenoid fossa (Zola&#39;s tubercle). See Zola A and Rothschild E. A.,  J Pros Dent ., “Posterior Condyle Positions in Unimpeded Jaw Movements,” 1962, September-October, Vol. II, No. 5, pp 873-881. 
   Derivation of Formula 
   Melvin Moss posited the “Functional Matrix” in several papers in the 1970&#39;s. In Dental Clinics of North America, (Moss Melvin L, DCNA,  A Functional Cranial Analysis of CentriC Relation,  1975,0 I,. Vol., L9, No. 3 pp 431-442), Melvin clearly described the functional relationship between muscles and the joints they operate. Melvin described the adaptation of the boney matrix to the dynamics and strength of muscle function. 
   This concept was then applied to the dynamics of motion and bony architecture of the TMJ. Zola described a functional guiding plane on the medial wall of the glenoid fossa-now called Zola&#39;s tubercle. Zola A and Rothschild E. A., J Pros Dent.,  Posterior Condyle Positions in Unimpeded Jaw Movements, r 962, September-October, Vol. II, No. 5, PP 873-881. 
   The inventors have found that the motion of the mandible in the chewing function is controlled by the strength of the major muscles used in food mastication and in the guide planes of contacting teeth and the angle of guidance in the TM joint, more specifically, Zola&#39;s tubercle can be used to construct a predictive TMJ joint operational model. The contact of the teeth in forceful use (mastication and/or bruxism or thegosis) moderate and coordinate the power of muscle contraction and these combined forces influence the pressure within the gliding TM joint which effects modification of the boney surfaces, through adaptive change, to be in harmony with the guiding teeth. 
   The techniques of measurement, data collection and confirmation were described. Using data from physical examination of human dry specimen skulls (digital laser scan images) described above, a mathematical formula to solve the complex relationships of contact and function was determined using non-linear vector mathematical analysis. 
   To confirm that a single mathematical model could determine the motion of an TMJ, two separate problems had to be solved and the results related. 
   The derivation of the mathematical equation or mathematical model that permits an understanding and determination of the functioning of the TMJ under normal and forceful use can then be used to construct improved TMJ prostheses, improved crown, cap and bridge designs, improve teeth placement augmentation (orthodontic) and improve mouth guards to correct, prevent or reduce damage to the TMJ. The model can also be used to determine proper reconstruction and/or augmentation surgical procedures. 
   The basic realization leading to the mathematical model was that jaw motion about the TMJ pivot is constrained by wear surfaces on two specific teeth as shown in  FIG. 10 . The geometry is defined by position vectors (P 1  and P 2 ) and normal vectors (N 1  and N 2 ) and a TMJ vector (P TMJ ) as shown in  FIG. 11 . The jaw rigid body motion is defined by a rotation about a directional vector at the TMJ pivot as shown in  FIG. 12 . The rotation vector must be in the planes defined by the position vectors and the normal vectors to the wear surfaces as shown in  FIG. 13 . 
   The equations that model the motion are given below: 
                     N   -&gt;     P1     =         P   -&gt;     1     ×       N   -&gt;     1               (   2   )                   N   -&gt;     P2     =         P   -&gt;     2     ×       N   -&gt;     2               (   3   )                   E   -&gt;     R     =           N   -&gt;     P1     ×       N   -&gt;     P2                  N   -&gt;     P1     ×       N   -&gt;     P2                      (   4   )                 {           P     TMJ   ⁢           ⁢   x     ′               P     TMJ   ⁢           ⁢   y     ′               P     TMJ   ⁢           ⁢   z     ′           }     =       [           ⁢             cos   ⁢           ⁢   Φ     +       E     R   ⁢           ⁢   x     2     ⁡     (     1   -     cos   ⁢           ⁢   Φ       )                   E     R   ⁢           ⁢   x       ⁢       E     R   ⁢           ⁢   y       ⁡     (     1   -     cos   ⁢           ⁢   Φ       )         +       E     R   ⁢           ⁢   z       ⁢   sin   ⁢           ⁢   Φ                 E     R   ⁢           ⁢   x       ⁢       E     R   ⁢           ⁢   z       ⁡     (     1   -     cos   ⁢           ⁢   Φ       )         -       E     R   ⁢           ⁢   y       ⁢   sin   ⁢           ⁢   Φ                     E     R   ⁢           ⁢   x       ⁢       E     R   ⁢           ⁢   y       ⁡     (     1   -     cos   ⁢           ⁢   Φ       )         -       E     R   ⁢           ⁢   z       ⁢   sin   ⁢           ⁢   Φ               cos   ⁢           ⁢   Φ     +       E     R   ⁢           ⁢   y     2     ⁡     (     1   -     cos   ⁢           ⁢   Φ       )                   E     R   ⁢           ⁢   y       ⁢       E     R   ⁢           ⁢   z       ⁡     (     1   -     cos   ⁢           ⁢   Φ       )         +       E     R   ⁢           ⁢   x       ⁢   sin   ⁢           ⁢   Φ                     E     R   ⁢           ⁢   x       ⁢       E     R   ⁢           ⁢   z       ⁡     (     1   -     cos   ⁢           ⁢   Φ       )         +       E     R   ⁢           ⁢   y       ⁢   sin   ⁢           ⁢   Φ                 E     R   ⁢           ⁢   y       ⁢       E     R   ⁢           ⁢   z       ⁡     (     1   -     cos   ⁢           ⁢   Φ       )         -       E     R   ⁢           ⁢   x       ⁢   sin   ⁢           ⁢   Φ               cos   ⁢           ⁢   Φ     +       E     R   ⁢           ⁢   z     2     ⁡     (     1   -     cos   ⁢           ⁢   Φ       )               ⁢           ]     ⁢     {           P     TMJ   ⁢           ⁢   x                 P     TMJ   ⁢           ⁢   y                 P     TMJ   ⁢           ⁢   z             }               (   1   )               
where P′ TMJx , P′ TMJy , and P′ TMJz  are components of an opposite TMJ trajectory vector P′ TMJ , Φ is any rotation about the TMJ pivot point, P TMJx , P TMJy , and P TMJz  are components of the TMJ position vector P TMJ , and a 3×3 rotation matrix set forth in brackets. In the rotation matrix, E rx , E ry , and E Rz  are components of a vector E R , N 1  and N 2  are the normals to the planes and P 1  and P 2  are the position vectors from the TMJ pivot to the normals and N P1  and N P2  are vector cross products as described in. The opposite TMJ trajectory is constrained to the surface of a sphere centered at the opposite TMJ pivot point with a radius of |P TMJ | as shown in  FIGS. 14A and 14B . Input data provided via Excel spread sheet: location coordinates for TMJ; location coordinates for constraint teeth, and normal vectors for wear surfaces on constraint teeth. TMJ trajectories plotted on y-z side view along the y reference axis. Trajectories are reported as displacement from the jaw center configuration and left TMJ and right TMJ plotted separately.
 
   A mathematical model was constructed to solve for anyone of three missing guide planes given the other two. In the initial study, the givens were two tooth contacts on opposite sides of the arch and the solution was for the TMJ guide plane. The accuracy of the mathematical model was calculated at 0.99611%, an error factor of 0.088102% (less than 1%), and the solution for the guide angle of Zola&#39;s tubercle was consistent with the analog model. 
   Initial conclusion that can be drawn from this study: (a) There are cross arch tooth contacts which guide the motion of the mandible; (b) The tooth contacts modify muscle function—speed, trajectory and power; (c) TMJ boney surfaces adapt to this muscle function to guide the jaw in harmony with the guidance of the teeth in forceful use of the masticatory system. 
   Utilizing the A4 skull data to build simplified trajectory test case: TMJ and tooth location unchanged; normal wear surfaces are horizontal; and TMJ trajectory normal vectors are simple normalization of relative position vectors X and Y components from opposite TMJ and horizontal motion normal vector Z component should be zero. TMJ trajectories should result in a pure horizontal motion. TMJ trajectory normal vectors should coincide with relative position vector X and Y components. Referring now to  FIGS. 15A&amp;B , test data for the A4 skull, TMJ and tooth locations, and TMJ trajectories of skull sample A4. The TMJ trajectory test case results are as follows: calculated TMJ trajectories are purely horizontal; calculated trajectory normal vectors coincide exactly with relative position vector constraints—edata is the normal vector input data—normal is the calculated value: 
           edat   =              0.9961111           0.088102           0                          normal   =              0.996111           0.008102           0                    
The mathematical trajectory model yields the expected results for the simplified test case.
 
   Referring now to  FIGS. 16A&amp;B , actual case data for the A4 skull, TMJ and tooth locations, and TMJ trajectories of skull sample A4. The TMJ trajectory test case results are as follows: calculated TMJ trajectories are purely horizontal; calculated trajectory normal vectors coincide exactly with relative position vector constraints. 
   Referring now to  FIGS. 17A&amp;B , actual case data for the A4 skull, TMJ and tooth locations, and TMJ trajectories of skull sample A10. The TMJ trajectory test case results are as follows: calculated TMJ trajectories are purely horizontal; calculated trajectory normal vectors coincide exactly with relative position vector constraints. 
   Referring now to  FIGS. 18A&amp;B , actual case data for the A4 skull, TMJ and tooth locations, and TMJ trajectories of skull sample A12. The TMJ trajectory test case results are as follows: calculated TMJ trajectories are purely horizontal; calculated trajectory normal vectors coincide exactly with relative position vector constraints. 
   Referring now to  FIGS. 19A&amp;B , actual case data for the A4 skull, TMJ and tooth locations, and TMJ trajectories of skull sample A14. The TMJ trajectory test case results are as follows: calculated TMJ trajectories are purely horizontal; calculated trajectory normal vectors coincide exactly with relative position vector constraints. 
   Referring now to  FIGS. 20A&amp;B , actual case data for the A4 skull, TMJ and tooth locations, and TMJ trajectories of skull sample A. The TMJ trajectory test case results are as follows: calculated TMJ trajectories are purely horizontal; calculated trajectory normal vectors coincide exactly with relative position vector constraints. 
   Referring now to  FIGS. 21A&amp;B , actual case data for the A4 skull, TMJ and tooth locations, and TMJ trajectories of skull sample Anat. The TMJ trajectory test case results are as follows: calculated TMJ trajectories are purely horizontal; calculated trajectory normal vectors coincide exactly with relative position vector constraints. 
   Referring now to  FIGS. 22A&amp;B , actual case data for the A4 skull, TMJ and tooth locations, and TMJ trajectories of skull sample C13. The TMJ trajectory test case results are as follows: calculated TMJ trajectories are purely horizontal; calculated trajectory normal vectors coincide exactly with relative position vector constraints. 
   TMJ Physiological Supposition 
   The skull structure conforms to TMJ motion, adjusting over time. Because of this conformity, the inventors expected a correlation to exist between the predicted TMJ trajectories within the mandibular cavity. The data showed that localized skull surfaces conform to the predicted TMJ trajectory and that normal vector to surfaces were theorized to be correlated to normal vector calculated from the TMJ trajectory. 
   The TMJ trajectory normal vector is constrained to the plane of the circular trajectory path as shown in  FIG. 23 . 
   View Parallel TMJ Trajectory 
   Multiple possibilities of surface orientations exist between the TMJ and the mandibular cavity as shown in  FIG. 24 . Trajectory normal correlation to measured normal must adjust the measurement to account for that vector component out of plane from the trajectory as shown in Figure. 
   Trajectory Normal Vector 
   The trajectory normal vector is defined via vector mathematics. Calculate the instantaneous trajectory tangent vector and trajectory normal vector is normal to both tangent and the rotation orientation vectors: 
             T   TMJ     =         ⅆ     P   TMJ         ⅆ   Φ                ⅆ     P   TMJ         ⅆ   Φ                            N   TMJ     =         T   TMJ     ×     E   R                T   TMJ     ×     E   R                    
Measured Mandibular Cavity Normal Vector
 
   The measured normal vector is adjusted to remove vector components normal to the calculated trajectory plane. Calculate the normal vector components normal to the calculated TMJ trajectory plane. Subtract from the measured normal and renormalized to unit length: 
             N   projected     =       (       E   R     ·     N   measured       )     *     E   R                     N   adjusted     =       N   measured     -     N   projected                     N   adjusted     =       N   adjusted            N   adjusted                  
Mandidular Cavity Normal Vector Correlation
 
   Referring now to  FIGS. 25A&amp;B , mandibular cavity normal vector correlation for skulls A4, A10, A12, A14, A, Anat and C13 are shown in tabulated form. The data indicates that trajectories are highly variable from individual to individual and from jaw to jaw on any particular individual. The data does not support a correlation with mandibular cavity wear surface measured normal vectors. However, the determined variation is proposed to be within selected data collection sites which are to be re-evaluated. The mathematical model is proposed to be accurate and valid. 
   All references cited herein are incorporated by reference. While this invention has been described fully and completely, it should be understood that, within the scope of the appended claims, the invention may be practiced otherwise than as specifically described. Although the invention has been disclosed with reference to its preferred embodiments, from reading this description those of skill in the art may appreciate changes and modification that may be made which do not depart from the scope and spirit of the invention as described above and claimed hereafter.