Three-dimensional dynamical capabilities of the human masticatory muscles

Clinical finite element analysis of mandibular displacement model treated with Twin-block appliance

INTRODUCTION

The mechanical distribution of the mandible is an important factor that affects functional orthosis during Twin-block (TB) appliance correction. Changes in the mandible before and after TB appliance correction are also key factors in maintaining the therapeutic effect. Finite element analysis, a powerful numerical, analytical tool, is widely used to predict the stress and strain distribution of the craniofacial bone that orthodontics generates.

METHODS

The sample was a 14-year-old male patient with Class II malocclusion during growth. A cone-beam computed tomography scan was undertaken at pretreatment and posttreatment. In the Finite element analysis of the pretreatment model, the remote displacement model of the mandible was established with the sella point as the center. A mandibular model under TB appliance loading was established. Its mandibular displacement and von Mises stress were compared before and after loading. Three-dimensional registration was conducted on the pretreatment and posttreatment models to measure the sagittal displacement of the centrosome.

RESULTS

The force on the mandible occurred mainly in the condyle neck and medial mandible after the TB appliance moved the mandible. After displacement, the posterior upper margin of the condyle was farther away from the articular fossa. Three-dimensional registration results showed that new bone had formed behind and above the condyle after TB appliance treatment.

CONCLUSION

The TB appliance provides additional advantages in treating skeletal Class II malocclusions by helping to reduce the burden on the temporomandibular joint and promoting the adaptive reconstruction of the mandible.

Physiological oculo-auricular-facial-mandibular synkinesis elicited in humans by gaze deviations

Integrated motor behaviors involving ocular motion-associated movements of the head, neck, pinna, and parts of the face are commonly seen in animals orienting to a visual target. A number of coordinated movements have also been observed in humans making rapid gaze shifts to horizontal extremes, which may be vestiges of these. Since such integrated mechanisms point to a non-pathological co-activation of several anatomically separate cranial circuits in humans, it is important to see how the different pairs of integrative motor behaviors with a common trigger (i.e., ocular motion) manifest in relation to one another. Here, we systematically examined the pattern of eye movement-induced recruitment of multiple cranial muscles in humans. Simultaneous video-oculography and bilateral surface electromyograms of transverse auricular, temporalis, frontalis, and masseter muscles were recorded in 15 healthy subjects (8 females; 29.3±5.2 years) while they made head-fixed, horizontal saccadic, pursuit and optokinetic eye movements. Potential chin laterotrusion linked to contractions of masticator muscles was captured with a yaw-fixed accelerometer. Our findings objectively show an orchestrated aural-facial-masticatory muscle response to a range of horizontal eye movements (prevalence of 21-93%). These responses were most prominent during eccentric saccades. We further reveal distinctions between the various observed activation patterns in terms of their profile (transient or sustained), laterality (with respect to direction of gaze) and timing (with respect to saccade onset). Possible underlying neural substrates, their atavistic behavioral significance, and potential clinical applications for monitoring sensory attention and designing attention-directed hearing aids in the future are discussed.

Biomechanical analysis of occlusal modes on the periodontal ligament while orthodontic force applied

OBJECTIVE

The study objective was to investigate four common occlusal modes by using the finite element (FE) method and to conduct a biomechanical analysis of the periodontal ligament (PDL) and surrounding bone when orthodontic force is applied.

MATERIALS AND METHODS

A complete mandibular FE model including teeth and the PDL was established on the basis of cone-beam computed tomography images of an artificial mandible. In the FE model, the left and right mandibular first premolars were not modeled because both canines required distal movement. In addition, four occlusal modes were simulated: incisal clench (INC), intercuspal position (ICP), right unilateral molar clench (RMOL), and right group function (RGF). The effects of these four occlusal modes on the von Mises stress and strain of the canine PDLs and bone were analyzed.

RESULTS

Occlusal mode strongly influenced the distribution and value of von Mises strain in the canine PDLs. The maximum von Mises strain values on the canine PDLs were 0.396, 1.811, 0.398, and 1.121 for INC, ICP, RMOL, and RGF, respectively. The four occlusal modes had smaller effects on strain distribution in the cortical bone, cancellous bone, and miniscrews.

CONCLUSION

Occlusal mode strongly influenced von Mises strain on the canine PDLs when orthodontic force was applied.

CLINICAL RELEVANCE

When an FE model is used to analyze the biomechanical behavior of orthodontic treatments, the effect of muscle forces caused by occlusion must be considered.

Patient-specific finite element models of the human mandible: Lack of consensus on current set-ups

The use of finite element analysis (FEA) has increased rapidly over the last decennia and has become a popular tool to design implants, osteosynthesis plates and prostheses. With increasing computer capacity and the availability of software applications, it has become easier to employ the FEA. However, there seems to be no consensus on the input variables that should be applied to representative FEA models of the human mandible. This review aims to find a consensus on how to define the representative input factors for a FEA model of the human mandible. A literature search carried out in the PubMed and Embase database resulted in 137 matches. Seven papers were included in this current study. Within the search results, only a few FEA models had been validated. The material properties and FEA approaches varied considerably, and the available validations are not strong enough for a general consensus. Further validations are required, preferably using the same measuring workflow to obtain insight into the broad array of mandibular variations. A lot of work is still required to establish validated FEA settings and to prevent assumptions when it comes to FEA applications.

Mandible Integrity and Material Properties of the Periodontal Ligament during Orthodontic Tooth Movement: A Finite-Element Study

We used the finite-element method (FEM) to investigate the effects of jawbone model integrity and the material properties of the periodontal ligament (PDL) on orthodontic tooth movement. Medical imaging software and computer-aided design software were used to create finite-element models of a partial and complete mandibles based on dental cone beam computed tomography images of the human skull. Additionally, we exerted an orthodontic force on the canine crown in the direction of an orthodontic miniscrew under a lower molar root to compare the von Mises strain on the canine PDL in three models: a partial mandible model under orthodontic force (Model 1), a complete mandible model under orthodontic force (Model 2), and a complete mandible model under orthodontic force with clench occlusion in the intercuspal position (ICP; Model 3). Additionally, in the complete mandible model under orthodontic force with ICP occlusion, we analyzed the effects of a PDL with a low (Model 4), moderate (Model 5), and high (Model 6) linear elastic modulus and a PDL a bilinear elastic modulus (Model 7). The simulation results for mandible integrity indicated that the maximum von Mises strains on the canine PDL for Models 1, 2, and 3 were 0.461, 0.394, and 1.811, respectively. Moreover, for the models with different PDL material properties, the maximum von Mises strains on the canine PDLs for Models 4, 5, 6, and 7 were 6.047, 2.594, 0.887, and 1.811, respectively. When the FEM was used to evaluate tooth movement caused by orthodontic force, the transformation of a complete mandible model into a partial mandible model or alteration of the elastic modulus of the PDL influenced the biomechanical responses of the PDL. Additionally, the incorporation of daily ICP occlusion resulted in a larger effect.

Ultrasonographic characterization of lingual structures pertinent to oral, periodontal, and implant surgery

OBJECTIVES

Increased applications of ridge augmentation in the lingual posterior mandible call for an urgent need to study its anatomy. Therefore, our first aim was to validate ultrasound in measuring the mandibular lingual structures in human cadavers. Secondarily, to test its feasibility in imaging the lingual nerve in live humans.

MATERIALS AND METHODS

Nine fresh un-embalmed fully/partially edentulous cadaver heads were utilized for aim 1. Three areas in the lingual mandible were imaged (mandibular premolar, molar, and retromolar). Immediately after, biopsies were harvested from each site. The thickness of the mucosa, mylohyoid muscle, and lingual nerve diameter was measured via ultrasound and statistically compared to histology. Similarly, the lingual nerve in live humans was also imaged.

RESULTS

None of the differences between the ultrasound and histology measurements reached statistical significance (p > .05). The mean mucosal thickness via ultrasound and histology was 1.45 ± 0.49 and 1.39 ± 0.50 mm, 5 mm lingual to the mylohyoid muscle attachment. At 10 mm beyond the attachment, the ultrasound and histologic values were 1.54 ± 0.48 and 1.37 ± 0.49, respectively. The mean muscle thickness measured via ultrasound and histology was 2.31 ± 0.56 and 2.25 ± 0.47 mm, at the 5 mm distance. At the 10 mm distance, the measurements were 2.46 ± 0.56 and 2.36 ± 0.5 mm, respectively. The mean ultrasonic lingual nerve diameter was 2.38 ± 0.44 mm, versus 2.43 ± 0.42 mm, with histology. The lingual nerve diameter on 19 live humans averaged to 2.01 ± 0.35 mm (1.4-3.1 mm).

CONCLUSIONS

Within its limitations, ultrasound accurately measured mandibular lingual soft tissue structures on cadavers, and the lingual nerve on live humans.

Modelling of the forces acting on the human stomatognathic system during dynamic symmetric incisal biting of foodstuffs

A major stage in the preparation of a computational model of the human stomatognathic system is the determination of the values of the forces for the adopted loading configuration. In physiological conditions, food is a factor having a significant effect on the values of the loads acting on the stomatognathic system. Considering that the act of mastication is a complex process, this research undertook to determine the forces (bite forces, muscular forces and temporomandibular joint reaction forces) acting on the stomatognathic system during the dynamic symmetric incisal biting of selected foodstuffs. The investigations were divided into two stages: (1) experimental tests and (2) numerical simulations. In the first stage, classic force-displacement characteristic curves (F_i-Δh) were determined for the food while in the second stage, the curves were used as a dynamic stomatognathic system model load function. One of the most important results of this research is that the food characteristic in the form of a force-displacement function has been shown to have a significant effect not only on the values of the muscular forces and the temporomandibular joint reaction forces, but also on their curves during the dynamic loading of the stomatognathic system. The analysis of the results indicates that F_i-Δh has an effect on not only the (active and passive) forces, but also on other parameters, such as stress, deformation, displacement, and probably the rigidity of the muscles.

Challenges in computational modelling of food breakdown and flavour release

A dynamic, three dimensional (3D) computational model that predicts the breakdown of food and the release of tastants and aromas could enhance the understanding of how food is perceived during consumption. This model could also shorten the development process of new foods because many virtual foods could be assessed, and discarded if unsuitable, before any physical prototyping is required. The construction and testing of a complete 3D model of mastication presents many challenges including an accurate representation of: the anatomical movements of the oral cavity (including the teeth, tongue, cheeks and palates), the breakdown behaviour of the food, the interactions between comminuted food and saliva as the bolus is formed, the release and transport of taste and aromas and how these physical and chemical processes are perceived by a person. These challenges are discussed in reference to previous experimental and simulation work and using results of new applications of a coupled biomechanical-smoothed particle hydrodynamics (B-SPH) model. The B-SPH model is demonstrated to simulate several complicated aspects of mastication including: (1) the sensitivity of particle size to changes in the movements of the jaw and tongue; (2) large strain behaviour of food due to softening by heating; (3) interactions between solid and liquid food components; (3) the release of tastants into the saliva; and (4) the transport of tastants to the taste buds. These applications show the possibilities of a model to viably simulate mastication, but highlight the many modelling and experimental challenges that remain.

Influence of mastication and edentulism on mandibular bone density

The aim of this study was to demonstrate that external loading due to daily activities, including mastication, speech and involuntary open-close cycles of the jaw contributes to the internal architecture of the mandible. A bone remodelling algorithm that regulates the bone density as a function of stress and loading cycles is incorporated into finite element analysis. A three-dimensional computational model is constructed on the basis of computerised tomography (CT) images of a human mandible. Masticatory muscle activation involved during clenching is modelled by static analysis using linear optimisation. Other loading conditions are approximated by imposing mandibular flexure. The simulations predict that mandibular bone density distribution results in a tubular structure similar to what is observed in the CT images. Such bone architecture is known to provide the bone optimum strength to resist bending and torsion during mastication while reducing the bone mass. The remodelling algorithm is used to simulate the influence of edentulism on mandibular bone loss. It is shown that depending on the location and number of missing teeth, up to one-third of the mandibular bone mass can be lost due to lack of adequate mechanical stimulation.

Masticatory loadings and cranial deformation in Macaca fascicularis: a finite element analysis sensitivity study

Biomechanical analyses are commonly conducted to investigate how craniofacial form relates to function, particularly in relation to dietary adaptations. However, in the absence of corresponding muscle activation patterns, incomplete muscle data recorded experimentally for different individuals during different feeding tasks are frequently substituted. This study uses finite element analysis (FEA) to examine the sensitivity of the mechanical response of a Macaca fascicularis cranium to varying muscle activation patterns predicted via multibody dynamic analysis. Relative to the effects of varying bite location, the consequences of simulated variations in muscle activation patterns and of the inclusion/exclusion of whole muscle groups were investigated. The resulting cranial deformations were compared using two approaches; strain maps and geometric morphometric analyses. The results indicate that, with bite force magnitude controlled, the variations among the mechanical responses of the cranium to bite location far outweigh those observed as a consequence of varying muscle activations. However, zygomatic deformation was an exception, with the activation levels of superficial masseter being most influential in this regard. The anterior portion of temporalis deforms the cranial vault, but the remaining muscles have less profound effects. This study for the first time systematically quantifies the sensitivity of an FEA model of a primate skull to widely varying masticatory muscle activations and finds that, with the exception of the zygomatic arch, reasonable variants of muscle loading for a second molar bite have considerably less effect on cranial deformation and the resulting strain map than does varying molar bite point. The implication is that FEA models of biting crania will generally produce acceptable estimates of deformation under load as long as muscle activations and forces are reasonably approximated. In any one FEA study, the biological significance of the error in applied muscle forces is best judged against the magnitude of the effect that is being investigated.

Developing a musculoskeletal model of the primate skull: predicting muscle activations, bite force, and joint reaction forces using multibody dynamics analysis and advanced optimisation methods

An accurate, dynamic, functional model of the skull that can be used to predict muscle forces, bite forces, and joint reaction forces would have many uses across a broad range of disciplines. One major issue however with musculoskeletal analyses is that of muscle activation pattern indeterminacy. A very large number of possible muscle force combinations will satisfy a particular functional task. This makes predicting physiological muscle recruitment patterns difficult. Here we describe in detail the process of development of a complex multibody computer model of a primate skull (Macaca fascicularis), that aims to predict muscle recruitment patterns during biting. Using optimisation criteria based on minimisation of muscle stress we predict working to balancing side muscle force ratios, peak bite forces, and joint reaction forces during unilateral biting. Validation of such models is problematic; however we have shown comparable working to balancing muscle activity and TMJ reaction ratios during biting to those observed in vivo and that peak predicted bite forces compare well to published experimental data. To our knowledge the complexity of the musculoskeletal model is greater than any previously reported for a primate. This complexity, when compared to more simple representations provides more nuanced insights into the functioning of masticatory muscles. Thus, we have shown muscle activity to vary throughout individual muscle groups, which enables them to function optimally during specific masticatory tasks. This model will be utilised in future studies into the functioning of the masticatory apparatus.

Combining geometric morphometrics and functional simulation: an emerging toolkit for virtual functional analyses

The development of virtual methods for anatomical reconstruction and functional simulation of skeletal structures offers great promise in evolutionary and ontogenetic investigations of form-function relationships. Key developments reviewed here include geometric morphometric methods for the analysis and visualization of variations in form (size and shape), finite element methods for the prediction of mechanical performance of skeletal structures under load and multibody dynamics methods for the simulation and prediction of musculoskeletal function. These techniques are all used in studies of form and function in biology, but only recently have they been combined in novel ways to facilitate biomechanical modelling that takes account of variations in form, can statistically compare performance, and relate performance to form and its covariates. Here we provide several examples that illustrate how these approaches can be combined and we highlight areas that require further investigation and development before we can claim a mature theory and toolkit for a statistical biomechanical framework that unites these methods.

Current computational modelling trends in craniomandibular biomechanics and their clinical implications

Computational models of interactions in the craniomandibular apparatus are used with increasing frequency to study biomechanics in normal and abnormal masticatory systems. Methods and assumptions in these models can be difficult to assess by those unfamiliar with current practices in this field; health professionals are often faced with evaluating the appropriateness, validity and significance of models which are perhaps more familiar to the engineering community. This selective review offers a foundation for assessing the strength and implications of a craniomandibular modelling study. It explores different models used in general science and engineering and focuses on current best practices in biomechanics. The problem of validation is considered at some length, because this is not always fully realisable in living subjects. Rigid-body, finite element and combined approaches are discussed, with examples of their application to basic and clinically relevant problems. Some advanced software platforms currently available for modelling craniomandibular systems are mentioned. Recent studies of the face, masticatory muscles, tongue, craniomandibular skeleton, temporomandibular joint, dentition and dental implants are reviewed, and the significance of non-linear and non-isotropic material properties is emphasised. The unique challenges in clinical application are discussed, and the review concludes by posing some questions which one might reasonably expect to find answered in plausible modelling studies of the masticatory apparatus.

Masticatory motion is controlled in humans by a limited set of muscle synergies

The masticatory motion, whereby food introduced into the mouth is processed into a bolus suitable for swallowing, can be divided into successive masticatory cycles, each comprising downward and subsequent upward movements of the mandible. The present study deals with the problem of the existence of muscle synergies in mastication, that is whether some of the muscles involved in mastication receive common motor drives, rather than controlled individually. Evidence for muscle synergy during mastication is scarce, partly due to the difficulties in simultaneous recording of the electromyographic (EMG) activities from all the muscles involved. Thus, we analyzed the variability of the mandibular motion during mastication rather than to examine the EMG patterns, based on the hypothesis that a motion elicited by a limited set of muscle synergies can be approximated as a superposition of the same number of independent motions. Mandibular motion paths were recorded from 8 healthy males (25-31 years), who chewed gum or gummy candy. A morphometric technique, which describes the shape of a closed curve by using normalized elliptic Fourier descriptors and reduces the variance of the shape by using principal component analysis, was applied to analyze the variability of the mandibular motion paths. We found three independent variations of the motion paths, whose linear combinations accounted for an average of 93% (range, 88-96%) of the total variance. The extracted variations were similar among the subjects. These findings provide indirect evidence for the existence of a limited set of muscle synergies for mastication in humans.

A dynamic model of jaw and hyoid biomechanics during chewing

Our understanding of human jaw biomechanics has been enhanced by computational modelling, but comparatively few studies have addressed the dynamics of chewing. Consequently, ambiguities remain regarding predicted jaw-gapes and forces on the mandibular condyles. Here, we used a new platform to simulate unilateral chewing. The model, based on a previous study, included curvilinear articular guidance, a mobile hyoid apparatus, and a compressible food bolus. Muscles were represented by Hill-type actuators with drive profiles tuned to produce target jaw and hyoid movements. The cycle duration was 732 ms. At maximum gape, the lower incisor-point was 20.1mm down, 5.8mm posterior, and 2.3mm lateral to its initial, tooth-contact position. Its maximum laterodeviation to the working-side during closing was 6.1mm, at which time the bolus was struck. The hyoid's movement, completed by the end of jaw-opening, was 3.4mm upward and 1.6mm forward. The mandibular condyles moved asymmetrically. Their compressive loads were low during opening, slightly higher on the working-side at bolus-collapse, and highest bilaterally when the teeth contacted. The model's movements and the directions of its condylar forces were consistent with experimental observations, resolving seeming discordances in previous simulations. Its inclusion of hyoid dynamics is a step towards modelling mastication.

A call to action for bioengineers and dental professionals: directives for the future of TMJ bioengineering

The world's first TMJ Bioengineering Conference was held May 25-27, 2006, in Broomfield, Colorado. Presentations were given by 34 invited speakers representing industry, academics, government agencies such as NIH, and private practice, which included surgeons, engineers, biomedical scientists, and patient advocacy leaders. Other attendees included documentary film makers and FDA officials. The impetus for the conference was that the field of TMJ research has been lacking continuity, with no open forum available for surgeons, scientists, and bioengineers to exchange scientific and clinical ideas and identify common goals, strengths, and capabilities. The goal was thus to plant the seeds for establishing a forum for multidisciplinary and interdisciplinary interactions. The collective wisdom and interactions brought about by a melting pot of these diverse individuals has been pooled and is disseminated in this article, which offers specific directives to bioengineers, basic scientists, and medical and dental professionals including oral and maxillofacial surgeons, pain specialists, orthodontists, prosthodontists, endocrinologists, rheumatologists, immunologists, radiologists, neurologists, and orthopaedic surgeons. A primary goal of this article was to attract researchers across a breadth of research areas to lend their expertise to a significant clinical problem with a dire need for new talent. For example, researchers with expertise in finite element modeling will find an extensive list of clinically significant problems. Specific suggestions for TMJ research were presented by the leading organizations for TMJ surgeons and TMJ patients, and further research needs were identified in a series of group discussions. The specific needs identified at the conference and presented here will be essential for those who endeavor to engage in TMJ research, especially in the areas of tissue engineering and biomechanics. Collectively, it is our hope that many of the questions and directives presented here find their way into the proposals of multidisciplinary teams across the world with new and promising approaches to diagnose, prevent and treat TMJ disorders.

Validation of a musculo-skeletal model of the mandible and its application to mandibular distraction osteogenesis

Mandibular distraction osteogenesis will lead to a change in muscle coordination and load transfer to the temporomandibular joints (TMJ). The objective of this work is to present and validate a rigid-body musculo-skeletal model of the mandible based on inverse dynamics for calculation of the muscle activations, muscle forces and TMJ reaction forces for different types of clenching tasks and dynamic tasks. This approach is validated on a symmetric mandible model and an application will be presented where the TMJ reaction forces during unilateral clenching are estimated for a virtual distraction patient with a shortened left ramus. The mandible model consists of 2 rigid segments and has 4 degrees-of-freedom. The model was equipped with 24 hill-type musculotendon actuators. During the validation experiment one subject was asked to do several tasks while measuring EMG activity, bite force and kinematics. The bite force and kinematics were used as input for the simulations of the same tasks after which the estimated muscle activities were compared with the measured muscle activities. This resulted in an average correlation coefficient of 0.580 and an average of the Mean Absolute Error of 0.109. The virtual distraction model showed a large difference in the TMJ reaction forces between left and right compared with the symmetric model for the same loading case. The present work is a step in the direction of building patient-specific mandible models, which can assess the mechanical effects on the TMJ before mandibular distraction osteogenesis surgery.

Combined finite-element and rigid-body analysis of human jaw joint dynamics

The jaw joint plays a crucial role in human mastication. It acts as a guidance for jaw movements and as a fulcrum for force generation. The joint is subjected to loading which causes tensions and deformations in its cartilaginous structures. These are assumed to be a major determinant for development, maintenance and also degeneration of the joint. To analyze the distribution of tensions and deformations in the cartilaginous structures of the jaw joint during jaw movement, a dynamical model of the human masticatory system has been constructed. Its movements are controlled by muscle activation. The articular cartilage layers and articular disc were included as finite-element (FE) models. As this combination of rigid-body and FE modeling had not been applied to musculoskeletal systems yet, its benefits and limitations were assessed by simulating both unloaded and loaded jaw movements. It was demonstrated that joint loads increase with muscle activation, irrespective of the external loads. With increasing joint load, the size of the stressed area of the articular surfaces was enlarged, whereas the peak stresses were much less affected. The results suggest that the articular disc enables distribution of local contact stresses over a much wider area of the very incongruent articular surfaces by transforming compressive principal stress into shear stress.

Using the finite element method to model the biomechanics of the asymmetric mandible before, during and after skeletal correction by distraction osteogenesis

An approach was developed to evaluate the load transfer mechanism in the temporomandibular joint (TMJ) area before, during and after mandibular ramus elongation by distraction osteogenesis (DO). In a concerted approach using computer tomography, magnetic resonance imaging (MRI), and finite element analysis, three-dimensional numerical models based on a young male patient, with a dento-facial deformity were generated. The magnitude and direction of the muscle forces acting on the mandible were assessed using both values derived from the muscles volume and cross-section as retrieved from the MRI-scan data-sets and taken from the literature. The resistance of the soft tissue envelope towards elongation during the DO-phase was also included. The finite element analyses showed that before skeletal correction by DO the load transfer was asymmetrical with high peak stresses in the affected joint. Following ramus elongation a more symmetrical loading in TMJs was predicted. The reaction forces in the TMJs during DO were low.

Three-Dimensional Finite Element Analysis of the Mandible and Temporomandibular Joint on Simulated Occlusal Forces before and after Vertical Ramus Elongation by Distraction Osteogenesis

Distraction osteogenesis has recently become a mainstay for treatment of mandibular hypoplasia. Thorough knowledge about changes in the temporomandibular joint (TMJ) and the surrounding parts of the mandible and the skull after mandibular distraction is still lacking. The purpose of the current study was to investigate the stress distribution in the mandible and the TMJ before and after skeletal correction by intraoral unilateral vertical mandibular ramus distraction, using a finite element (FE) model. The FE models were based on computed tomography scans and magnetic resonance imaging scans of a patient with unilateral hypoplasia of the right mandibular ramus caused by juvenile idiopathic arthritis. The character of stress distribution in the mandible and TMJ before and after skeletal correction by 15 mm of vertical distraction of the mandibular ramus was analyzed quantitatively and compared during centric occlusion. Before the distraction osteogenesis treatment, the condyles, articular discs, and glenoid fossa regions are loaded with a different stress pattern. The affected right condyle, disc, and fossa are loaded diffusely and externally in comparison with the anterior and with centralized loading on the normal left side. After unilateral mandibular distraction osteogenesis, the load became more centric and symmetrical. The results suggest that correction of the mandibular deformity by distraction osteogenesis tends to normalize the stress patterns in the TMJ.

A method for material parameter determination for the human mandible based on simulation and experiment

In cranio-maxillofacial surgery planning and implant design, it is important to know the elastic response of the mandible to load forces as they occur, e.g., in biting. The goal of the present study is to provide a method for a quantitative determination of material parameters for the human jaw bone, whose values can, e.g., be used to devise a prototype plastic model for the mandible. Non-destructive load experiments are performed on a cadaveric mandible using a specially designed test bed. The identical physiological situation is simulated in a computer program. The underlying mathematical model is based on a two component, linear elastic material law. The numerical realization of the model, difficult due to the complex geometry and morphology of the mandible, is via the finite element (FE) method. Combining the validated simulation with the results of the tests, an inverse problem for the determination of Young's modulus and the Poisson ratio of both cortical and cancellous bone can then be solved.

Functional significance of the coupling between head and jaw movements

When humans open or close the jaw they also move the head. Unintentionally, it rotates backwards when the jaw opens and returns upon jaw closure. We hypothesized that this mutual movement coupling is related to the muscles in the floor of the mouth. A biomechanical model was applied to comprehend the functional significance of this movement coupling. As the jaw opened the jaw opening muscles shortened and became less forceful. Meanwhile they had to stretch the jaw closing muscles. The simulations showed that a simultaneous head extension facilitated jaw opening. A possible functional significance for the coupling between head and jaw movements is that it can extend jaw gape. Head extension can contribute to a wider jaw gape by on the one hand a reduced shortening of the jaw opening muscles and on the other hand by a reorientation of these muscles so that they obtain a more favorable position for jaw opening.

Neuromuscular objectives of the human masticatory apparatus during static biting

OBJECTIVE

The central nervous system controls the muscles of mastication and may dictate muscle outputs according to a biologically important objective. This study tested the hypotheses that (a) the effective sagittal TMJ eminence morphology, and (b) the outputs of the masticatory muscles during static biting, are consistent with minimisation of joint loads or minimisation of muscle effort.

DESIGN

Numerical modelling predicted effective eminence morphology (from sagittal plane directions of TMJ force for centred loading over a range from molar to incisor biting) and TMJ and muscle forces during static unilateral biting in seven subjects. In vivo effective eminence morphology was measured from jaw tracking recorded from each subject. Muscle activities during biting tasks on first molar and incisor teeth were measured by electromyography using surface or indwelling electrodes.

RESULTS

Subject-specific predicted effective eminence morphology correlated with in vivo data (0.85< or =R2< or =0.99). Mixed and random coefficient analysis of covariance indicated good agreement between predicted and measured muscle outputs for all muscles of mastication investigated. Individual linear regression analysis showed that modelled muscle outputs accurately predicted EMG data, with average errors of 8% for molar and 15% for incisor biting.

CONCLUSIONS

Effective sagittal eminence morphology was consistent with minimisation of joint loads for all subjects. Masticatory muscle outputs during unilateral biting were consistent with minimisation of joint loads or minimisation of muscle effort, or both, depending on the subject. These results are believed to be the first to test model predictions of muscle output during biting for all muscles of mastication.

Muscle and temporomandibular joint forces associated with chincup loading predicted by numerical modeling

Development of the components of the temporomandibular joint (TMJ) is thought to reflect joint loading. The aims of this project were to test 3 hypotheses: whether effective eminence morphology, masticatory muscle forces, and predicted TMJ forces during chincup loading of the mandible were consistent with the objectives of minimization of joint loads (MJL) or muscle effort (MME), or both. Regression relationships of MJL model-predicted versus measured eminence shapes in 9 subjects indicated a high degree of correlation (mean slope = 0.99, compared with perfect-match slope = 1.00). Model predictions of muscle output during chincup loading of the mandible were tested by comparison with data gathered in 6 subjects. Midsagittal plane chin loads were applied over a range of 60 degrees while bilateral masticatory muscle surface electromyography was quantified. The regression relationships of predicted versus measured masseter and anterior digastric muscle outputs indicated that model predictions were highly correlated (mean slope (masseter muscle) = 1.02; mean slope (digastric muscle) = 0.96). TMJ forces predicted by modeling showed intersubject differences of up to 34% for similar chincup loading conditions. Intrasubject variation in TMJ forces was as high as 57%, depending on chin load angle. The results demonstrated that TMJ eminence shape and masticatory muscle forces were consistent with objectives of both MJL and MME. Variation in TMJ forces depended on the subject and the direction of chincup loading.

A biomechanical model for analysis of muscle force, power output and lower jaw motion in fishes

Fish skulls are complex kinetic systems with movable components that are powered by muscles. Cranial muscles for jaw closing pull the mandible around a point of rotation at the jaw joint using a third-order lever mechanism. The present study develops a lever model for the jaw of fishes that uses muscle design and the Hill equation for nonlinear length-tension properties of muscle to calculate dynamic power output. The model uses morphometric data on skeletal dimensions and muscle proportions in order to predict behavior and force transmission mediated by lever action. The computer model calculates a range of dynamic parameters of jaw function including muscle force, torque, effective mechanical advantage, jaw velocity, bite duration, bite force, work and power. A complete list of required morphometrics is presented and a software program (MandibLever 2.0) is available for implementing lever analysis. Results show that simulations yield kinematics and timing profiles similar to actual fish feeding events. Simulation of muscle properties shows that mandibles reach their peak velocity near the start of jaw closing, peak force at the end of jaw closing, and peak power output at about 25% of the closing cycle time. Adductor jaw muscles with different mechanical designs must have different contractile properties and/or different muscle activity patterns to coordinate jaw closing. The effective mechanical advantage calculated by the model is considerably lower than the mechanical advantage estimated from morphological lever ratios, suggesting that previous studies of morphological lever ratios have overestimated force and underestimated velocity transmission to the mandible. A biomechanical model of jaw closing can be used to interpret the mechanics of a wide range of jaw mechanisms and will enable studies of the functional results of developmental and evolutionary changes in skull morphology and physiology.

Human masticatory muscle forces during static biting

Muscle forces determine joint loads, but the objectives governing the mix of muscle forces involved are unknown. This study tested the hypothesis that masticatory muscle forces exerted during static biting are consistent with objectives of minimization of joint loads (MJL) or muscle effort (MME). To do this, we compared numerical model predictions with data measured from six subjects. Biting tasks which produced moments on molar and incisor teeth were modeled based on MJL or MME. The slope of predicted vs. electromyographic (EMG) data for an individual was compared with a perfect match slope of 1.00. Predictions based on MME matched best with EMG activity for molar biting (slopes, 0.89-1.16). Predictions from either or both models matched EMG results for incisor biting (best-match slopes, 0.95-1.07). Muscle forces during isometric biting appear to be consistent with objectives of MJL or MME, depending on the individual, biting location, and moment.

Dynamics of the human masticatory system

In this review, the movement characteristics of the human masticatory system are discussed from a biomechanical perspective. The discussion is based upon the three fundamental laws of mechanics applied to the various anatomical structures that are part of the masticatory system. An analysis of the forces and torques applied to the mandible by muscles, joints, articular capsules, and teeth is used to assess the determinants of jaw movement. The principle of relating the interplay of forces to the center of gravity of the lower jaw, in contrast to a hinge axis near its joints, is introduced. It is evident that the muscles are the dominant determinants of jaw movement. The contributions of the individual muscles to jaw movements can be derived from the orientation of their lines of action with respect to the center of gravity of the lower jaw. They cause the jaw to accelerate with six degrees of freedom. The ratio between linear and angular accelerations is subtly dependent on the mass and moments of inertia of the jaw, and the structures that are more or less rigidly attached to it. The effects of articular forces must be taken into account, especially if the joints are loaded asymmetrically. The muscles not only move the jaw but also maintain articular stability during midline movements. Passive structures, such as the ligaments, become dominant only when the jaw reaches its movement boundaries. These ligaments are assumed to prevent joint dislocation during non-midline movements.

Validated numerical modeling of the effects of combined orthodontic and orthognathic surgical treatment on TMJ loads and muscle forces

Investigations of the changes in the mechanics of the craniomandibular system as a result of treatment have been limited by the lack of validated models of this system. The aims of this project were to (1) validate numerical model predictions of temporomandibular joint (TMJ) eminence morphology and muscle forces produced during molar biting and (2) use the validated models to calculate the changes in TMJ and muscle forces as a consequence of treatment involving orthognathic surgery. Ten volunteers participated; their combined orthodontic and orthognathic surgical treatments were completed. Three-dimensional anatomical data from each subject were used in computer models to predict the sagittal TMJ eminence morphology and joint and muscle forces for each subject, consistent with the neuromuscular objectives of minimizing joint loads and muscle effort. The actual shape of the eminence in each subject was measured with jaw tracking. Surface electromyographic recordings were a measure of the muscle forces involved in static molar biting. Model predictions were compared with measured data from the subjects for eminence shape (R(2) = 0.96) and for muscle activity ratios (R(2) = 0.98). The strength of these relationships validated the models for use in calculating changes in joint loads and muscle forces after treatment. The results suggested that the mechanics of the masticatory system are affected by the combined treatments. The TMJ loads increased in 8 subjects. The average increases in condylar and muscle forces were 4% relative to the applied bite force, but in 1 case the increases were up to 20%. Therefore, although average increases in the forces were small, some persons may experience biologically significant increases in joint and muscle forces as a result of treatment.

The three-dimensional active envelope of jaw border movement and its determinants

The sagittal and frontal active envelope of border movement is applied regularly as a clinical tool in functional examinations of the human masticatory system. In contrast, the three-dimensional movement area has hardly been examined. Furthermore, the determinants of this area are not established unambiguously. In the present study, the three-dimensional envelope of incisor movement was predicted with a three-dimensional mathematical model of the human masticatory system, which included the morphology of the system and the fine architecture of its muscles. With this model, the influence of the temporomandibular ligaments and the passive muscle tensions on the envelope were estimated. The predicted three-dimensional active envelope of border movements was limited in horizontal directions, predominantly by the temporomandibular ligaments. The passive tensions of the masticatory muscles influenced, although marginally, its vertical extension. It appeared unlikely that, in a normal situation, active muscle tensions (casu quo muscle reflexes) contribute to the shape of the envelope.

A method to predict muscle control in the kinematically and mechanically indeterminate human masticatory system

A method is proposed to generate muscle activation patterns for goal-directed movements of the human masticatory system. This system is special because apart from a larger amount of muscles than degrees of freedom its joints do not restrict its movements a priori. Therefore, each muscle is able to influence all six degrees of freedom which makes the system kinematically and mechanically indeterminate. Furthermore, its working space is principally determined by the dynamical properties of its muscles and not by passive constraints. The presented method determines the contribution of each degree of freedom to a movement of a reference point on the mandible. It avails of straightforward mathematical techniques like Linear Programming. It does not require a separate trajectory planning step. It was applied in a six degrees of freedom dynamical mathematical model of the human masticatory system. This model which was based upon rigid-body dynamics incorporating skull morphology and muscle architecture including dynamical properties. Movements were exclusively defined by a goal position of the mandibular reference point. The method proved to be robust in generating muscle activation patterns for both feasible and infeasible movement tasks. Generally, they were accomplished faster than habitually observed. If the task was infeasible the movement stopped at the outer boundary of the working space at the side of the unreachable goal. The method, therefore, enables to explore the working space of the mandible and the factors that are relevant for its boundaries.

Three-dimensional finite element analysis of the human temporomandibular joint disc

A three-dimensional finite element model of the articular disc of the human temporomandibular joint has been developed. The geometry of the articular cartilage and articular disc surfaces in the joint was measured using a magnetic tracking device. First, polynomial functions were fitted through the coordinates of these scattered measurements. Next, the polynomial description was transformed into a triangulated description to allow application of an automatic mesher. Finally, a finite element mesh of the articular disc was created by filling the geometry with tetrahedral elements. The articulating surfaces of the mandible and skull were modeled by quadrilateral patches. The finite element mesh and the patches were combined to create a three-dimensional model in which unrestricted sliding of the disc between the articulating surfaces was allowed. Simulation of statical joint loading at the closed jaw position predicted that the stress and strain distributions were located primarily in the intermediate zone of the articular disc with the highest values in the lateral part. Furthermore, it was predicted that considerable deformations occurred for relatively small joint loads and that relatively large variations in the direction of joint loading had little influence on the distribution of the deformations.