Effect of optimization criterion on spinal force estimates during asymmetric lifting

Optimization-based muscle force prediction models of the lumbar region are used in research and ergonomic practice, usually as a subroutine of a job analysis software package. Various optimization criteria have been put forward for use in rationally selecting a set of muscle forces to satisfy moment equilibrium, including the sum of cubed muscle contraction intensities and a double linear programming procedure for minimizing the spinal compression force and maximum muscle contraction intensity. A laboratory study was conducted to determine whether these two model formulations produce significantly different estimates of spinal forces for a dynamic asymmetric lift. Although statistically significant differences were found between the predictions of the two models, the difference in peak spinal compression force was only 1.1%.     

A new model for calculating muscle forces from electromyograms

A muscle model is described that uses electromyogram (EMG), muscle length and speed of contraction to predict muscle force. Physiological parameters are the Hill constants and the shape of the twitch response to a single stimulus. The model was incorporated in a jaw model of the rabbit and tested by predicting the bite force produced by the jaw muscles during mastication. The time course of the calculated force appeared to match the bite force, measured in vivo by a strain gauge, applied to the bone below the teeth. The variation in peak strain amplitude from cycle to cycle correlated with the variation predicted by the model. The peak amplitude of the integrated EMGs of individual jaw muscles showed an average correlation with peak strain of 0.41. Use of the sum of the available peak amplitudes, weighted according to their effect upon the bite force increased the correlation to 0.46; the model predicted bite forces showed a correlation of 0.57 with the strain. The increase in correlation was statistically significant. The muscle forces were calculated using a minimum number of easily obtainable constants.     

Investigation on masticatory muscular functionality following oral reconstruction – An inverse identification approach

The human masticatory system has received significant attention in the areas of biomechanics due to its sophisticated co-activation of a group of masticatory muscles which contribute to the fundamental oral functions. However, determination of each muscular force remains fairly challenging in vivo; the conventional data available may be inapplicable to patients who experience major oral interventions such as maxillofacial reconstruction, in which the resultant unsymmetrical anatomical structure invokes a more complex stomatognathic functioning system. Therefore, this study aimed to (1) establish an inverse identification procedure by incorporating the sequential Kriging optimization (SKO) algorithm, coupled with the patient-specific finite element analysis (FEA) in silico and occlusal force measurements at different time points over a course of rehabilitation in vivo; and (2) evaluate muscular functionality for a patient with mandibular reconstruction using a fibula free flap (FFF) procedure. The results from this study proved the hypothesis that the proposed method is of certain statistical advantage of utilizing occlusal force measurements, compared to the traditionally adopted optimality criteria approaches that are basically driven by minimizing the energy consumption of muscle systems engaged. Therefore, it is speculated that mastication may not be optimally controlled, in particular for maxillofacially reconstructed patients. For the abnormal muscular system in the patient with orofacial reconstruction, the study shows that in general, the magnitude of muscle forces fluctuates over the 28-month rehabilitation period regardless of the decreasing trend of the maximum muscular capacity. Such finding implies that the reduction of the masticatory muscle activities on the resection side might lead to non-physiological oral biomechanical responses, which can change the muscular activities for stabilizing the reconstructed mandible.     

The functional correlates of jaw‐muscle fiber architecture in tree‐gouging and nongouging callitrichid monkeys

Common (Callithrix jacchus) and pygmy (Cebuella pygmaea) marmosets and cotton‐top tamarins (Saguinus oedipus) share broadly similar diets of fruits, insects, and tree exudates. Marmosets, however, differ from tamarins in actively gouging trees with their anterior dentition to elicit tree exudates flow. Tree gouging in common marmosets involves the generation of relatively wide jaw gapes, but not necessarily relatively large bite forces. We compared fiber architecture of the masseter and temporalis muscles in C. jacchus (N = 18), C. pygmaea (N = 5), and S. oedipus (N = 13). We tested the hypothesis that tree‐gouging marmosets would exhibit relatively longer fibers and other architectural variables that facilitate muscle stretch. As an architectural trade‐off between maximizing muscle excursion/contraction velocity and muscle force, we also tested the hypothesis that marmosets would exhibit relatively less pinnate fibers, smaller physiologic cross‐sectional areas (PCSA), and lower priority indices (I) for force. As predicted, marmosets display relatively longer‐fibered muscles, a higher ratio of fiber length to muscle mass, and a relatively greater potential excursion of the distal tendon attachments, all of which favor muscle stretch. Marmosets further display relatively smaller PCSAs and other features that reflect a reduced capacity for force generation. The longer fibers and attendant higher contraction velocities likely facilitate the production of relatively wide jaw gapes and the capacity to generate more power from their jaw muscles during gouging. The observed functional trade‐off between muscle excursion/contraction velocity and muscle force suggests that primate jaw‐muscle architecture reflects evolutionary changes related to jaw movements as one of a number of functional demands imposed on the masticatory apparatus. Am J Phys Anthropol, 2009. © 2009 Wiley‐Liss, Inc.     

The functional significance of primate mandibular form.

A stress analysis of the primate mandible suggests that vertically deep jaws in the molar region are usually an adaptation to counter increased sagittal bending stress about the balancing-side mandibular corpus during unilateral mastication. This increased bending stress about the balancing side is caused by an increase in the amount of balancing-side muscle force. Furthermore, this increased muscle force will also cause an increase in dorso-ventral shear stress along the mandibular symphysis. Since increased symphyseal stress can be countered by symphyseal fusion and as increased bending stress can be countered by a deeper jaw, deep jaws and symphyseal fusion are often part of the same functional pattern. In some primates (e.g., Cercocebus albigena), deep jaws are an adaptation to counter bending in the sagittal plane during powerful incisor biting, rather than during unilateral mastication. The stress analysis of the primate mandible also suggests that jaws which are transversely thick in the molar region are an adaptation to counter increased torsion about the long axis of the working-side mandibular corpus during unilateral mastication. Increased torsion of the mandibular corpus can be caused by an increase in masticatory muscle force, an increase in the transverse component of the postcanine bite force and/or an increase in premolar use during mastication. Patterns of masticatory muscle force were estimated for galagos and macaques, demonstrating that the ratio of working-side muscle force to balancing-side muscle force is approximately 1.5:1 in macaques and 3.5:1 in galagos during unilateral isometric molar biting. These data support the hypothesis that mandibular symphyseal fusion is an adaptative response to maximize unilateral molar bite force by utilizing a greater percentage of balancing-side muscle force.     

The length dependence of muscle active force: considerations for parallel elastic properties.

The purpose of this study was to choose between two popular models of skeletal muscle: one with the parallel elastic component in parallel with both the contractile element and the series elastic component (model A), and the other in which it is in parallel with only the contractile element (model B). Passive and total forces were obtained at a variety of muscle lengths for the medial gastrocnemius muscle in anesthetized rats. Passive force was measured before the contraction (passive A) or was estimated for the fascicle length at which peak total force occurred (passive B). Fascicle length was measured with sonomicrometry. Active force was calculated by subtracting passive (A or B) force from peak total force at each fascicle or muscle length. Optimal length, that fascicle length at which active force is maximized, was 13.1 +/- 1.2 mm when passive A was subtracted and 14.0 +/- 1.1 mm with passive B (P < 0.01). Furthermore, the relationship between double-pulse contraction force and length was broader when calculated with passive B than with passive A. When the muscle was held at a long length, passive force decreased due to stress relaxation. This was accompanied by no change in fascicle length at the peak of the contraction and only a small corresponding decrease in peak total force. There is no explanation for the apparent increase in active force that would be obtained when subtracting passive A from the peak total force. Therefore, to calculate active force, it is appropriate to subtract passive force measured at the fascicle length corresponding to the length at which peak total force occurs, rather than passive force measured at the length at which the contraction begins.     

Mathematical Model for Skeletal Muscle to Simulate the Concentric and Eccentric Contraction

Skeletal muscles are responsible for the relative motion of the bones at the joints and provide the required strength. They exhibit highly nonlinear mechanical behaviour and are described by nonlinear hyperelastic constitutive relations. It is distinct from other biological soft tissue. Its hyperelastic or viscoelastic behaviour is modelled by using CE, SEE, and PEE. Contractile element simulates the behaviour of skeletal muscle when it is subjected to eccentric and concentric contraction. This research aims to estimate the stress induced in skeletal muscle in eccentric and concentric contraction with respect to the predefined strain. With the use of mathematical model for contraction of skeletal muscle for eccentric and concentric contraction, the stress induced in the skeletal muscle is estimated in this research. Mathematical model is developed for the muscle using EMG signals and Force-velocity relationship calculated. With the use of force-velocity of contraction of muscle, mathematical model is developed. This can be useful to understand the mechanical behaviour of skeletal muscles in eccentric and concentric contraction with clinical relevance. Authors are further working to develop the mathematical model with torsion force with proper activation function of muscle and experimentation for extraction of the anisotropic mechanical properties of skeletal muscle.     

Comparison of muscle weight and force ratios in New and Old World monkeys

Thin mandibles and small incisors found in New World monkeys as compared with Old World monkeys suggest that there may be differences in craniofacial loading patterns between these two groups, particularly in levels of mandibular corpus twisting (Hylander, 1975, 1979a; Eaglen, 1984; Bouvier, 1986a,b). This study examined the hypothesis that changes in the relative force contributions of the masticatory muscles were responsible for lowering torsion on the mandibular corpus in New World monkeys. Muscle weight and physiological cross-sections were compared using data from the literature (Schumacher, 1960: Turnbull, 1970; Cachel, 1979) as well as new data on adult male Cebus apella and Macaca mulatta. Both age and sex had an effect on muscle ratios. Mixed samples such as those used by Schumacher and Turnbull probably are not appropriate for drawing conclusions concerning species or group differences in muscle ratios. In addition, biomechanical conclusions based on muscle weight ratios alone to estimate muscle force may be misleading because fiber length inversely affects the amount of force a muscle can exert. A comparison of ratios based on physiological cross-section as an estimator of muscle force in New and Old World monkeys does not support the hypothesis that alterations in force contribution by individual masticatory muscles are responsible for minimizing mandibular corpus twisting in New World monkeys. Therefore, if twisting has been minimized in New World monkeys as suggested by their thin corpora, other changes in the craniofacial musculoskeletal complex, such as different muscle recruitment or pinnation patterns, may be responsible.     

Analysis of the tendinous structure in human masticatory muscles.

In the masticatory muscles, the development of bundles of the tendon was examined: they were composed of many collagen fibers and a few elastic fibers. In the masseter muscle, the property of the tendon differs in the distribution and size of collagen fibers and elastic fibers in comparison with those of other masticatory muscles. This difference is concerned with the kinetic force for the stress or the stretch of each tendon and muscle during jaw movement.     

Masticatory biomechanics and its relevance to early hominid phylogeny: an examination of palatal thickness using finite-element analysis

It has been proposed that morphological characters functionally related to mastication may be unreliable indicators of early hominid phylogeny. One hypothesis states that masticatory characters are highly prone to homoplasy. A second hypothesis states that such characters are likely to be morphologically integrated and thus violate the assumption of character independence implicit in all phylogenetic analyses. Evaluation of these hypotheses requires that masticatory features be accurately identified, but, to date, there have been relatively few attempts to test precisely which early hominid features are functionally related to chewing. This paper uses finite-element analysis to evaluate the functional relationships of a character--palatal thickness--that is one of several Paranthropus synapomorphies putatively related to mastication. A finite-element model of 145,680 elements was created from sixty-one 2-mm-thick CT scans of a Macaca fascicularis skull. The model was assigned the elastic properties of facial bone and loaded with muscle forces corresponding to the moment of centric occlusion during mastication. The model was constrained so as to produce a reaction force (corresponding to the bite force) at M(1). With a few exceptions, the strain patterns in the finite-element model compare well with those gathered from published and unpublished bone-strain experiments. The model was then modified to have a thick palate. The model was reloaded using an identical loading regime, and the strain patterns of the original and thick-palate models were compared. Although a thickened palate acts to reduce palatal strain, strains are elevated in other facial regions. This suggests that a thick palate would not have evolved in isolation as an adaptation to withstand masticatory stress. Rather, a thick palate may have evolved in concert with a suite of other facial features that share a stress-resistance function. This appears to be consistent with hypotheses positing that at least some facial features related to chewing evolved in an integrated fashion. More functional studies of other facial features are needed, as are formal studies of morphological integration.     

Stress and strain analysis of contractions during ramp distension in partially obstructed guinea pig jejunal segments

Previous studies have demonstrated morphological and biomechanical remodeling in the intestine proximal to an obstruction. The present study aimed to obtain stress and strain thresholds to initiate contraction and the maximal contraction stress and strain in partially obstructed guinea pig jejunal segments. Partial obstruction and sham operations were surgically created in mid-jejunum of male guinea pigs. The animals survived 2, 4, 7 and 14 days. Animals not being operated on served as normal controls. The segments were used for no-load state, zero-stress state and distension analyses. The segment was inflated to 10 cmH(2)O pressure in an organ bath containing 37°C Krebs solution and the outer diameter change was monitored. The stress and strain at the contraction threshold and at maximum contraction were computed from the diameter, pressure and the zero-stress state data. Young's modulus was determined at the contraction threshold. The muscle layer thickness in obstructed intestinal segments increased up to 300%. Compared with sham-obstructed and normal groups, the contraction stress threshold, the maximum contraction stress and the Young's modulus at the contraction threshold increased whereas the strain threshold and maximum contraction strain decreased after 7 days obstruction (P<0.05 and 0.01). In conclusion, in the partially obstructed intestinal segments, a larger distension force was needed to evoke contraction likely due to tissue remodeling. Higher contraction stresses were produced and the contraction deformation (strain) became smaller.     

Influence of occlusal stabilization splint on the asymmetric activity of masticatory muscles in patients with temporomandibular dysfunction

The aim of present study was to evaluate the symmetry of masticatory muscles' activity at various clenching levels in the intercuspal position in patients with functional disorders and in healthy subjects. The purpose was also to determine the effect of full-arch maxillary stabilization splint on the asymmetry of masticatory muscle activity in patients with temporomandibular dysfunction. In this study 6 TMD patients and 12 healthy subjects were investigated. Surface EMG recordings were obtained from left and right anterior temporal, left and right masseter and from the sub-mandibular group in the region of the anterior belly of the digastric muscle on the left and right side during clenching with the maximum 100% voluntary contraction (MVC) as well as during clenching at 50% and 25% of the maximum activity in the position of maximal intercuspation of teeth. In order to quantify asymmetrical masticatory muscle activity, the asymmetry index (AI) was calculated for each subject and for each muscle from the average anterior temporal, masseter and digastric potentials recorded during each test (100% MVC, 50% MVC and 25% MVC). In the group of patients EMG recordings were repeated during and after the splint therapy. The asymmetries of masticatory muscle activity was present in both groups, but in the group of TMD patients the asymmetry indices for anterior temporal muscle at 100% MVC (p = 0.049) and 50% MVC (p = 0.031) were significantly higher. Results have shown that the use of splint suppressed the asymmetry of all muscles, as during the splint therapy the asymmetry indices were lowered. After the therapy, the level of temporal muscle symmetry during submaximal clenching in the intercuspal position increased significantly (p = 0.046). This investigation points out that electromyography may be a valuable method of documenting that asymmetric activity of masticatory muscles improves after occlusal splint therapy in patients with TMD.     

Electromyography and mechanics of mastication in the albino rat.

The masticatory apparatus in the albino rat was studied by means of electromyography and subsequent estimation of muscular forces. The activity patterns of the trigeminal and suprahyoid musculature and the mandibular movements were recorded simultaneously during feeding. The relative forces of the individual muscles in the different stages of chewing cycles and biting were estimated on the basis of their physiological cross sections and their activity levels, as measured from integrated electromyograms. Workinglines and moment arms of these muscles were determined for different jaw positions. In the anteriorly directed masticatory grinding stroke the resultants of the muscle forces at each side are identical; they direct anteriorly, dorsally and slightly lingually and pass along the lateral side of the second molar. Almost the entire muscular resultant force is transmitted to the molars while the temporo-mandibular joint remains unloaded. A small transverse force, produced by the tense symphyseal cruciate ligaments balances the couple of muscle resultant and molar reaction force in the transverse plane. After each grinding stroke the mandible is repositioned for the next stroke by the overlapping actions of three muscle groups: the pterygoids and suprahyoids produce depression and forward shift, the suprahyoids and temporal backward shift and elevation of the mandible while the subsequent co-operation of the temporal and masseter causes final closure of the mouth and starting of the forward grinding movement. All muscles act in a bilaterally symmetrical fashion. The pterygoids contract more strongly, the masseter more weakly during biting than during chewing. The wide gape shifts the resultant of the muscle forces more vertically and moreposteriorly. The joint then becomes strongly loaded because the reaction forces are applied far anteriorly on the incisors. The charateristic angle between the almost horizontal biting force and the surface of the food pellet indicates that the lower incisors produce a chisel-like action. Tooth structure reflects chewing and biting forces. The transverse molar lamellae lie about parallel to the chewing forces whereas perpendicular loading of the occlusal surfaces is achieved by their inclination in the transverse plane. The incisors are loaded approximately parallel to their longitudinal axis, placement that avoids bending forces during biting. It is suggested that a predominantly protrusive musculature favors the effective force transmission to the lower incisors, required for gnawing. By grinding food across transversely oriented molar ridges the protrusive components of the muscles would be utilized best. From the relative weights of the masticatory muscles in their topographical relations with joints, molars and incisors it may be concluded that the masticatory apparatus is a construction adapted to optimal transmission of force from muscles to teeth.     

Mechanical work as predictor of force enhancement and force depression

The steady-state force following active muscle shortening or stretch differs from the maximum isometric force associated with the final length. This phenomenon proves that the isometric force production is not only dependent on current muscle length and length time derivative, but depends on the preceding contraction history. Isolated extensor digitorum longus and soleus muscles from mice (NMRI strain) were used to investigate the force produced by a muscle, and some parameters hypothetically influencing this history-dependent force modification. The muscles were pre-stimulated at a fixed length, then different stretch/shortening episodes were introduced, whereafter changes of the active force were recorded while the muscles were held isometrically to approach a steady-state force before de-stimulation. The mechanical work during active stretch and shortening was evaluated by integrating the product of force and ramp velocity over the length-varying period. The results show a negative linear correlation between the force modification and the mechanical work produced on or by the muscle, continuous between shortening and stretch. A corresponding modification of the passive force component following each stimulation was also observed. The conclusion is that the isometric force attained after stretch or shortening is well described by an asymptotic force which is determined by the mechanical work.     

Dynamics of myosin-driven skeletal muscle contraction: I. Steady-state force generation

Skeletal muscle contraction is a canonical example of motor-driven force generation. Despite the long history of research in this topic, a mechanistic explanation of the collective myosin force generation is lacking. We present a theoretical model of muscle contraction based on the conformational movements of individual myosins and experimentally measured chemical rate constants. Detailed mechanics of the myosin motor and the geometry of the sarcomere are taken into account. Two possible scenarios of force generation are examined. We find only one of the scenarios can give rise to a plausible contraction mechanism. We propose that the synchrony in muscle contraction is due to a force-dependent ADP release step. Computational results of a half sarcomere with 150 myosin heads can explain the experimentally measured force-velocity relationship and efficiency data. We predict that the number of working myosin motors increases as the load force is increased, thus showing synchrony among myosin motors during muscle contraction. We also find that titin molecules anchoring the thick filament are passive force generators in assisting muscle contraction.     

Assessment of intramuscular activation patterns using ultrasound M-mode strain

The intramuscular activation pattern can be connected to the motor unit recruitment strategy of force generation and fatigue resistance. Electromyography has earlier been used in several studies to quantify the spatial inhomogeneity of the muscle activation. We applied ultrasound M-mode strain to study the activation pattern through the tissue deformation. Correlation values of the strain at different force levels were used to quantify the spatial changes in the activation. The assessment was done including the biceps brachii muscle of 8 healthy subjects performing isometric elbow flexion contractions ranging from 0% to 80% of maximum voluntary contraction. The obtained results were repeatable and demonstrated consistent changes of the correlation values during force regulation, in agreement with previously presented EMG-results. Both intra-subject and inter-subject activation patterns of strain were considered along and transverse the fiber direction. The results suggest that ultrasound M-mode strain can be used as a complementary method to study intramuscular activation patterns with high spatial resolution.     

Passive and dynamic muscle architecture during transverse loading for gastrocnemius medialis in man

External forces from our environment impose transverse loads on our muscles. Studies in rats have shown that transverse loads result in a decrease in the longitudinal muscle force. Changes in muscle architecture during contraction may contribute to the observed force decrease. The aim of this study was to quantify changes in pennation angle, fascicle dimensions, and muscle thickness during contraction under external transverse load.Electrical stimuli were elicited to evoke maximal force twitches in the right calf muscles of humans. Trials were conducted with transverse loads of 2, 4.5, and 10 kg. An ultrasound probe was placed on the medial gastrocnemius in line with the transverse load to quantify muscle characteristics during muscle twitches.Maximum twitch force decreased with increased transverse muscle loading. The 2, 4.5, and 10 kg of transverse load showed a 9, 13, and 16% decrease in longitudinal force, respectively. Within the field of view of the ultrasound images, and thus directly beneath the external load, loading of the muscle resulted in a decrease in the muscle thickness and pennation angle, with higher loads causing greater decreases. During twitches the muscle transiently increased in thickness and pennation angle, as did fascicle thickness. Higher transverse loads showed a reduced increase in muscle thickness. Smaller increases in pennation angle and fascicle thickness strain also occurred with higher transverse loads.This study shows that increased transverse loading caused a decrease in ankle moment, muscle thickness, and pennation angle, as well as transverse deformation of the fascicles.     

Muscle extracellular matrix applies a transverse stress on fibers with axial strain

It is widely assumed that skeletal muscle contraction is isovolumic. This assumption has been verified at the single fiber and at the myofibril level. Model development and mechanical analyses often exploit this assumption when investigating skeletal muscle and evaluating muscle mechanical properties. This communication describes a method whereby individual muscle fibers and bundles of fibers, which include their constituent extracellular matrix (ECM), were tested to define the change in volume with axial strain. The results demonstrate that fibers are isovolumic, but bundles decrease in volume with strain. The loss of volume implicates a transverse force being applied to the fibers by the ECM. The nature and importance of this transverse force warrant further investigation.     

Normalizing upper trapezius EMG amplitude: Comparison of ramp and constant force procedures

The present study compared three procedures for normalization of upper trapezius surface electromyographic (EMG) amplitudes: (a) a ramp procedure (providing data in per cent of maximal voluntary contraction, MVC); (b) a constant force procedure based on two reference contractions (two-force procedure) (%MVC) and (c) a procedure expressing muscle activation in per cent of a reference voluntary electrical activity (%RVE). The study also evaluated the repeatability of the ramp and the RVE procedures and estimated the force exertion (%MVC) corresponding to the RVE. To illustrate the ergonomic effect of different normalization procedures, trapezius EMG during two work tasks was compared after normalization by the two-force and the RVE procedures. Fifteen subjects participated in the whole study. We found that force estimates obtained by the ramp procedure equation could be translated to force estimates obtained by the two-force procedure by the equation: %MVC_2force = − 0.6 + 0.9*%MVC_ramp, although with a considerable imprecision due to large inter-individual differences. In the ramp procedure, the intra-individual test-retest coefficient of variation (CV) depended on the force level; it was 45% at 5% MVC and 10% at 30% MVC. The CV of the RVE was 15%. The reference contraction used in the RVE procedure corresponded from 13–79% MVC (median 33%MVC). The load reducing effect of an ergonomic intervention was less obvious with the RVE procedure than with the two-force procedure due to a larger inter-individual variation. The advantages and disadvantages of the different procedures are discussed.