Quantitative calculations of temporomandibular joint reaction forces--I. The importance of the magnitude of the jaw muscle forces

The effect of measurement errors on quantitative calculation of temporomandibular joint reaction force was investigated in a two-dimensional, two-muscle model. A computer program using the model incremented the magnitude of the bite force and muscle forces and the lengths of their moment arms, and calculated the joint reaction force at each increment. Computation of the joint reaction force is most sensitive to the relative lengths of the bite force and muscle forces moment arms. Absolute values for each muscle force are not required and errors in the magnitudes of the muscle forces have only a minor effect on calculation of the total joint reaction force.     

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.     

Predicted pattern of human muscle activity during clenching derived from a computer assisted model: symmetric vertical bite forces

A computer assisted three-dimensional model of the jaw, based on linear programming, is presented. The upper and lower attachments of the muscles of mastication have been measured on a single human skull and divided into thirteen independent units on each side--a total of 26 muscle elements. The direction (in three dimensions) and maximum forces that could be developed by each muscle element, the bite reaction and two joint reactions are included in the model. It is shown for symmetrical biting that a model which minimizes the sum of the muscle forces used to produce a given bite force activates muscles in a way which corresponds well with previous observations on human subjects. A model which minimizes the joint reactions behaves differently and is rejected. An analysis of the way the chosen model operates suggests that there are two types of jaw muscles, power muscles and control muscles. Power muscles (superficial masseter, medial pterygoid and some of temporalis) produce the bite force but tend to displace the condyle up or down the articular eminence. This displacement is prevented by control muscles (oblique temporalis and lateral pterygoid) which have very poor moment arms for generating usual bite forces, but are efficient for preventing condylar slide. The model incorporates the concept that muscles consist of elements which can contract independently. It predicts that those muscle elements with longer moment arms relative to the joint are the first to be activated and, as the bite force increases, a ripple of activity spreads into elements with shorter moment arms. In general, the model can be used to study the three-dimensional activity in any system of joints and muscles.     

Biomechanics of the masticatory apparatus of Pteropus giganteus (Megachiroptera)

Jaw mechanics in Pteropus were studied by means of a three-dimensional model. The model included several parameters of muscle architecture, combined with quantified movement and electromyographical data. Estimates of the nature of the applied forces that act upon the mandible during a chewing cycle, and subsequent estimates of reaction forces at the bite point and joints during the powerstroke, were thus obtained for different food consistencies. The resultant muscle force (relative to the palate) shifts from upward and slightly backward at large gapes to upward and markedly backward at the end of closing. The resultant simultaneously moves anteriorly. During the powerstroke it retains a constant position and orientation along the thickened anterior edge of the coronoid process. The early stages of opening are guided by the slope of the teeth and mandibular fossa; during the remaining part of opening the working line of the resultant crosses the skull behind the joint and thus acquires an opening moment. The bite force has downward and forward components, and a slight transverse component. For a given applied muscular force its magnitude is larger in more posteriorly positioned bite points. Both joints are loaded, the contralateral one more than the ipsilateral. Food consistency affects magnitude and orientation of the applied force, and hence, magnitude and orientation of the bite force and magnitude of the joint reaction forces. The magnitude of masseter activity relative to temporalis activity appears to be the key factor for the orientation of the bite force, and hence for the mechanical optimal position of the food. The adaptive value of the general topography of the masticatory muscles in Pteropus is discussed.     

Individual muscle force parameters and fiber operating ranges for elbow flexion-extension and forearm pronation-supination

We have quantified individual muscle force and moment contributions to net joint moments and estimated the operating ranges of the individual muscle fibers over the full range of motion for elbow flexion/extension and forearm pronation/supination. A three dimensional computer graphics model was developed in order to estimate individual muscle contributions in each degree of freedom over the full range of motion generated by 17 muscles crossing the elbow and forearm. Optimal fiber length, tendon slack length, and muscle specific tension values were adjusted within the literature range from cadaver studies such that the net isometric joint moments of the model approximated experimental joint moments within one standard deviation. Analysis of the model revealed that the muscles operate on varying portions of the ascending limb, plateau region, and descending limb of the force-length curve. This model can be used to further understand isometric force and moment contributions of individual muscles to net joint moments of the arm and forearm and can serve as a comprehensive reference for the forces and moments generated by 17 major muscles crossing the elbow and wrist.     

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.     

Estimates of forces exerted by the jaw muscles of some reptiles

Masses and muscle fibre lengths have been determined for the jaw muscles of Chrysemys (Chelonia). Varanus (Lacertilia) and Caiman (Crocodilia). Hence the force that each muscle can exert has been estimated. The position and direction of the line of action of each force have also been determined. The results are used to calculate the magnitude and direction of the bite forces that each reptile could exert, at particular positions in its mouth, and the magnitude and direction of the reactions at the jaw articulation.     

Velocity-dependent cost function for the prediction of force sharing among synergistic muscles in a one degree of freedom model

Prediction of accurate and meaningful force sharing among synergistic muscles is a major problem in biomechanics research. Given a resultant joint moment, a unique set of muscle forces can be obtained from this mathematically redundant system using nonlinear optimization. The classical cost functions for optimization involve a normalization of the muscle forces to the absolute force capacity of the target muscles, usually by the cross-sectional area or the maximal isometric force. In a one degree of freedom model this leads to a functional relationship between moment arms and the predicted muscle forces, such that for constant moment arms, or constant ratios of moment arms, agonistic muscle forces increase or decrease in unison. Experimental studies have shown however that the relationship between muscle forces is highly task-dependent often causing forces to increase in one muscle while decreasing in a functional agonist, likely because of the contractile conditions and contractile properties of the involved muscles. We therefore, suggest a modified cost function that accounts for the instantaneous contraction velocity of the muscles and its effect on the instantaneous maximal force. With this novel objective function, a task-dependent prediction of muscle force distribution is obtained that allows, even in a one degree of freedom system, the prediction of force sharing loops, and simultaneously increasing and decreasing forces for agonist pairs of muscles.     

The human mandible: Lever or link?

The mammalian mandible, and in particular the human mandible, is generally thought to function as a lever during biting. This notion, however, has not gone unchallenged. Various workers have suggested that the mandible does not function as a lever, and they base this proposition on essentially two assertions: (1) the resultant of the forces produced by the masticatory muscles always passes through the bite point; (2) the condylar neck and/or the temporomandibular joint is unsuited to withstand reaction forces during biting. A review of the electromyographic data and of the properties of the tissues of the temporomandibular joint do not support the non-lever hypothesis of mandibular function. In addition, an analysis of the strength of the condylar neck demonstrates that this structure is strong enough to withstand the expected reaction force during lever action. Ordinarily the human mandible and the forces that act upon it are analyzed solely in the lateral projection. Moments are then usually analyzed about the mandibular condyle; however, some workers have advocated taking moments about other points, e.g., the instantaneous center of rotation. Obviously it makes no difference as to what point is chosen since the moments about any point during equilibrium conditions are equal to zero. It is also useful to analyze the forces acting on the mandible in the frontal projection, particularly during unilateral biting. The electromyographic data suggest that during powerful unilateral molar biting the resultant adductor muscle force is passing between the bite point and the balancing (non-biting side) condyle. Therefore, in order for this system to be in equilibrium there must be a reaction force acting on the balancing condyle. This suggests that reaction forces are larger on the balancing side than on the working side, and possibly explains why individuals with a painful temporomandibular joint usually prefer to bite on the side of the diseased joint.     

Coactivation of jaw muscles: recruitment order and level as a function of bite force direction and magnitude

The aim of this study was to obtain insight into the coactivation behaviour of the jaw muscles under various a priori defined static loading conditions of the mandible. As the masticatory system is mechanically redundant, an infinite number of recruitment patterns is theoretically possible to produce a certain bite force. Using a three-component force transducer and a feedback method, subjects could be instructed to produce a bite force of specific direction and magnitude under simultaneous registration of the EMG activity of anterior and posterior temporal, masseter and digastric muscles on each side. Forces were measured at the second premolars. Vertical, anterior, posterior, lateral and medial force directions were examined; in each direction force levels between 50 N and maximal voluntary force were produced. The results show that for all muscles the bite force-EMG relationship obeys a straight-line fit for forces exceeding 50 N. The relationship varies with bite force direction, except in the case of the digastric muscles. Variation is small for the anterior temporal and large for the posterior temporal and masseter muscles. The relative activation of muscles for a particular force in a particular direction in unique, despite the redundancy.     

A parametric model of muscle moment arm as a function of joint angle: Application to the dorsiflexor muscle group in mice

A parametric model was developed to describe the relationship between muscle moment arm and joint angle. The model was applied to the dorsiflexor muscle group in mice, for which the moment arm was determined as a function of ankle angle. The moment arm was calculated from the torque measured about the ankle upon application of a known force along the line of action of the dorsiflexor muscle group. The dependence of the dorsiflexor moment arm on ankle angle was modeled as r=R sin(a+Δ), where r is the moment arm calculated from the measured torque and a is the joint angle. A least-squares curve fit yielded values for R, the maximum moment arm, and Δ, the angle at which the maximum moment arm occurs as offset from 90°. Parametric models were developed for two strains of mice, and no differences were found between the moment arms determined for each strain. Values for the maximum moment arm, R, for the two different strains were 0.99 and 1.14 mm, in agreement with the limited data available from the literature. While in some cases moment arm data may be better fitted by a polynomial, use of the parametric model provides a moment arm relationship with meaningful anatomical constants, allowing for the direct comparison of moment arm characteristics between different strains and species.     

Inter-joint coupling effects on muscle contributions to endpoint force and acceleration in a musculoskeletal model of the cat hindlimb

The biomechanical principles underlying the organization of muscle activation patterns during standing balance are poorly understood. The goal of this study was to understand the influence of biomechanical inter-joint coupling on endpoint forces and accelerations induced by the activation of individual muscles during postural tasks. We calculated induced endpoint forces and accelerations of 31 muscles in a 7 degree-of-freedom, three-dimensional model of the cat hindlimb. To test the effects of inter-joint coupling, we systematically immobilized the joints (excluded kinematic degrees of freedom) and evaluated how the endpoint force and acceleration directions changed for each muscle in 7 different conditions. We hypothesized that altered inter-joint coupling due to joint immobilization of remote joints would substantially change the induced directions of endpoint force and acceleration of individual muscles. Our results show that for most muscles crossing the knee or the hip, joint immobilization altered the endpoint force or acceleration direction by more than 90° in the dorsal and sagittal planes. Induced endpoint forces were typically consistent with behaviorally observed forces only when the ankle was immobilized. We then activated a proximal muscle simultaneous with an ankle torque of varying magnitude, which demonstrated that the resulting endpoint force or acceleration direction is modulated by the magnitude of the ankle torque. We argue that this simple manipulation can lend insight into the functional effects of co-activating muscles. We conclude that inter-joint coupling may be an essential biomechanical principle underlying the coordination of proximal and distal muscles to produce functional endpoint actions during motor tasks.     

Comparative analysis of methods for determining bite force in the spiny dogfish Squalus acanthias

Many studies have identified relationships between the forces generated by the cranial musculature during feeding and cranial design. Particularly important to understanding the diversity of cranial form amongst vertebrates is knowledge of the generated magnitudes of bite force because of its use as a measure of ecological performance. In order to determine an accurate morphological proxy for bite force in elasmobranchs, theoretical force generation by the quadratomandibularis muscle of the spiny dogfish Squalus acanthias was modeled using a variety of morphological techniques, and lever-ratio analyses were used to determine resultant bite forces. These measures were compared to in vivo bite force measurements obtained with a pressure transducer during tetanic stimulation experiments of the quadratomandibularis. Although no differences were found between the theoretical and in vivo bite forces measured, modeling analyses indicate that the quadratomandibularis muscle should be divided into its constituent divisions and digital images of the cross-sections of these divisions should be used to estimate cross-sectional area when calculating theoretical force production. From all analyses the maximum bite force measured was 19.57 N. This relatively low magnitude of bite force is discussed with respect to the ecomorphology of the feeding mechanism of S. acanthias to demonstrate the interdependence of morphology, ecology, and behavior in organismal design.     

Passive force characteristics of an architecturally complex muscle

In architecturally complex muscles with large attachment areas, it can be expected that during movement different muscle regions undergo different amounts of length excursions. As a consequence, the amount of passive force produced by the regions will differ. Therefore, we tested the hypothesis that during movement the vector of the passive force of such a muscle, which defines the magnitude, position and orientation of the resultant force of the various regions, has no fixed position, between the muscle's center of origin and insertion. As a model for an architecturally complex muscle we used the masseter muscle. It was expected that during jaw opening anterior muscle regions are more stretched than posterior regions, leading to an anterior shift of the passive force vector. A three-component force transducer was used to measure both the position and magnitude of passive force in the masseter muscle of 9 rabbits. Forces were recorded during repeated cycles of stepwise opening and closure of the jaw. The muscle exhibited a clear hysteresis: passive force measured during jaw opening was larger than that during jaw closing. With an increase of the jaw gape there was an approximately exponential increase of the magnitude of the passive muscle force, while simultaneously the passive force vector shifted anteriorly. Moment arm length of passive force increased by about 100%. This anterior shift contributed substantially to the increase of the passive muscle moment generated during jaw opening. It can be concluded that in architecturally complex muscles the increase of the passive resistance moment which is associated with muscle lengthening might not only be due to an increase of the magnitude of passive muscle force but also to an increase of the moment arm of this force.     

The effects of the antagonist muscle force on intersegmental loading during isokinetic efforts of the knee extensors

Hamstrings activation when acting as antagonists is considered very important for knee joint stability. However, the effect of hamstring antagonist activity on knee joint loading in vivo is not clear. Therefore, the purpose of this study was to examine the differences in antagonistic muscle force and their effect on agonist muscle and intersegmental forces during isokinetic eccentric and concentric efforts of the knee extensors. Ten males performed maximum isokinetic eccentric and concentric efforts of the knee extensors at 30 degrees s(-1). The muscular and tibiofemoral joint forces were then estimated using a two-dimensional model with and without including the antagonist muscle forces. The antagonist moment was predicted using an IEMG-moment model. The predicted antagonist force reached a maximum of 2.55 times body weight (BW) and 1.16 BW under concentric and eccentric conditions respectively. Paired t-tests indicated that these were significantly different (p<0.05). A one-way analysis of variance indicated that when antagonist forces are included in the calculations the patella tendon, compressive and posterior shear joint forces are significantly higher compared to those calculated without including the antagonist forces. The anterior shear force was not affected by antagonist activity. The antagonists produce considerable force throughout the range of motion and affect the joint forces exerted at the knee joint. Further, it appears that the antagonist effect depends on the type of muscle action examined as it is higher during concentric compared to eccentric efforts of the knee extensors.     

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.     

A Neuro-Muscular Elasto-Dynamic Model of the Human Arm Part 2: Musculotendon Dynamics and Related Stress Effects

In this paper we develop an elasto-dynamic model of the human arm that includes effects of neuro-muscular control uponelastic deformation in the limb.The elasto-dynamic model of the arm is based on hybrid parameter multiple body systemvariational projection principles presented in the companion paper.Though the technique is suitable for detailed bone and jointmodeling,we present simulations for simplified geometry of the bones,discretized as Rayleigh beams with elongation,whileallowing for large deflections.Motion of the upper extremity is simulated by incorporating muscle forces derived from aHill-type model of musculotendon dynamics.The effects of muscle force are modeled in two ways.In one approach,aneffective joint torque is calculated by multiplying the muscle force by a joint moment ann.A second approach models themuscle as acting along a straight line between the origin and insertion sites of the tendon.Simple arm motion is simulated byutilizing neural feedback and feedforward control.Simulations illustrate the combined effects of neural control strategies,models of muscle force inclusion,and elastic assumptions on joint trajectories and stress and strain development in the bone andtendon.     

Comparison between in vivo and theoretical bite performance: Using multi-body modelling to predict muscle and bite forces in a reptile skull

In biomechanical investigations, geometrically accurate computer models of anatomical structures can be created readily using computed-tomography scan images. However, representation of soft tissue structures is more challenging, relying on approximations to predict the muscle loading conditions that are essential in detailed functional analyses. Here, using a sophisticated multi-body computer model of a reptile skull (the rhynchocephalian Sphenodon), we assess the accuracy of muscle force predictions by comparing predicted bite forces against in vivo data. The model predicts a bite force almost three times lower than that measured experimentally. Peak muscle force estimates are highly sensitive to fibre length, muscle stress, and pennation where the angle is large, and variation in these parameters can generate substantial differences in predicted bite forces. A review of theoretical bite predictions amongst lizards reveals that bite forces are consistently underestimated, possibly because of high levels of muscle pennation in these animals. To generate realistic bites during theoretical analyses in Sphenodon, lizards, and related groups we suggest that standard muscle force calculations should be multiplied by a factor of up to three. We show that bite forces increase and joint forces decrease as the bite point shifts posteriorly within the jaw, with the most posterior bite location generating a bite force almost double that of the most anterior bite. Unilateral and bilateral bites produced similar total bite forces; however, the pressure exerted by the teeth is double during unilateral biting as the tooth contact area is reduced by half.