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

Two versus three magnesium screws for osteosynthesis of mandibular condylar head fractures: a finite element analysis

OBJECTIVES

Previous finite element analyses (FEA) have shown promising results for using two titanium screws in treating mandibular condylar head fractures but limited mechanical stability of a two-screw osteosynthesis with magnesium screws. Given the potential benefits of magnesium screws in terms of biocompatibility and resorption, this study aimed to compare two- and three-screw osteosynthesis solutions for a right condylar head fracture (AO CMF type p) with magnesium screws with a FEA.

MATERIALS AND METHODS

A previously validated finite element model simulating a 350 N bite on the contralateral molars was used to analyze von Mises stress within the screws, fragment deformation, and fracture displacement. All screws were modeled with uniform geometric specifications mirroring the design of Medartis MODUS® Mandible Hexadrive cortical screws.

RESULTS

The three-screw configuration demonstrated lower values for all three parameters compared to the two-screw scenario. There was a 30% reduction in maximum von Mises stress for the top screw and a 46% reduction for the bottom screw.

CONCLUSIONS

Fracture treatment with three magnesium screws could be a valuable and sufficiently stable alternative to the established treatment with titanium screws. Further studies on screw geometry could help improve material stability under mechanical loading, enhancing the performance of magnesium screws in clinical applications.

CLINICAL RELEVANCE

The use of magnesium screws for osteosynthesis of mandibular condylar head fractures offers the benefit of reducing the need for second surgery for hardware removal. Clinical data is needed to determine whether the advantages of resorbable screw materials outweigh potential drawbacks in condylar head fracture treatment.

Prediction of jaw opening function after mandibular reconstruction using subject-specific musculoskeletal modelling

BACKGROUND

Mandibular reconstruction patients often suffer abnormalities in the mandibular kinematics. In silico simulations, such as musculoskeletal modelling, can be used to predict post-operative mandibular kinematics. It is important to validate the mandibular musculoskeletal model and analyse the factors influencing its accuracy.

OBJECTIVES

To investigate the jaw opening-closing movements after mandibular reconstruction, as predicted by the subject-specific musculoskeletal model, and the factors influencing its accuracy.

METHODS

Ten mandibular reconstruction patients were enrolled in this study. Cone-beam computed tomography images, mandibular movements, and surface electromyogram signals were recorded preoperatively. A subject-specific mandibular musculoskeletal model was established to predict surgical outcomes using patient-averaged muscle parameter changes as model inputs. Jaw bone geometry was replaced by surgical planning results, and the muscle insertion sites were registered based on the non-rigid iterative closest point method. The predicted jaw kinematic data were validated based on 6-month post-operative measurements. Correlations between the prediction accuracy and patient characteristics (age, pathology and surgical scope) were further analysed.

RESULTS

The root mean square error (RMSE) for lower incisor displacement was 31.4%, and the error for peak magnitude of jaw opening was 4.9 mm. Age, post-operative infection and radiotherapy influenced the prediction accuracy. The amount of masseter detachment showed little correlation with jaw opening.

CONCLUSION

The mandibular musculoskeletal model successfully predicted short-range jaw opening functions after mandibular reconstruction. It provides a novel surgical planning method to predict the risk of developing trismus.

Muscle and joint mechanics during maximum force biting following total temporomandibular joint replacement surgery

Total temporomandibular joint replacement (TMJR) surgery is the established treatment for severe temporomandibular joint disorders. While TMJR surgery is known to increase mouth-opening capacity, reduce pain and improve quality of life, little is known about post-surgical jaw function during activities of daily living such as biting and chewing. The aim of this study was to use subject-specific 3D bite force measurements to evaluate the magnitude and direction of joint loading in unilateral total TMJR patients and compare these data to those in healthy control subjects. An optoelectronic tracking system was used to measure jaw kinematics while biting a rubber sample for 5 unilateral total TMJR patients and 8 controls. Finite element simulations driven by the measured kinematics were employed to calculate the resultant bite force generated when compressing the rubber between teeth during biting tasks. Subject-specific musculoskeletal models were subsequently used to calculate muscle and TMJ loading. Unilateral total TMJR patients generated a bite force of 249.6 ± 24.4 N and 164.2 ± 62.3 N when biting on the contralateral and ipsilateral molars, respectively. In contrast, controls generated a bite force of 317.1 ± 206.6 N. Unilateral total TMJR patients biting on the contralateral molars had a significantly higher lateral TMJ force direction (median difference: 63.6°, p = 0.028) and a significantly lower ratio of working TMJ force to bite force (median difference: 0.17, p = 0.049) than controls. Results of this study may guide TMJ prosthesis design and evaluation of dental implants.

In-silico study of type 'B' condylar head fractures and evaluating the influence of two positional screw distance in two-screw osteosynthesis construct

Clinical fixation screws are common in clinical practices to fix mandibular condyle fractures. Evidence suggests significance of 'working length' that is, distance between proximal and distal fixation screws in proximity to the fracture in orthopaedic implant design. In pursuit of stable implant-bone construct, this study aims to investigate the biomechanical performance of each configuration considered in the study and provide an optimal working length between the screws for clinical reference. Finite element models of virtually designed broken condyle as type 'B' were simulated and analysed in ANSYS Workbench. Screws are implanted according to previous literature at five varied distances 'd' maintaining five different ratios with the fracture length 'D'. Based on a literature review, boundary conditions, muscle traction forces and non-linear contacts were assigned to obtain precise results. Each case is considered an individual configuration and von Mises distribution, microstrain in bone, screw-bone interface micromotion and fracture dislocation were evaluated for all these configurations. Stress-shielding phenomenon is observed for maximum von Mises stresses in bone. Microstrain concentration was significant in cancellous bone in the vicinity of the screw around the fracture line. Configurations were compared based on the stress-strain along with micromotion to support the required amount of osseointegration between implant and bone. Presented data from all five conditions supported the assumption that under physiological loading conditions, the D3 configuration provided stability for fracture healing. Further research on screw shapes, diameters and material properties, or investigating the direction of forces within the screws could provide further insight into this topic.

Designing customized temporomandibular fossa prosthesis based on envelope surface of condyle movement: validation via in silico musculoskeletal simulation

Objective: This study presents an innovative articular fossa prosthesis generated by the envelope surface of condyle movement, and compares its mandible movements, muscle activities, and joint reaction forces with two temporomandibular joint (TMJ) prostheses using multibody musculoskeletal simulation. Methods: A healthy 23-year-old female was recruited for this study. Cone-beam computed tomographic (CBCT) was performed to reconstruct the mandibular bone geometry. A customized TMJ fossa prosthesis was designed based on the subject-specific envelope surface of condyle movement (ESCM). Mandibular kinematics and jaw-closing muscle electromyography (EMG) were simultaneously recorded during maximum jaw opening-closing movements. To validate our prosthesis design, a mandibular musculoskeletal model was established using flexible multibody dynamics and the obtained kinematics and EMG data. The Biomet fossa prosthesis and the ellipsoidal fossa prosthesis designed by imitating the lower limb prostheses were used for comparison. Simulations were performed to analyze the effects of different fossa prostheses on jaw opening-closing motions, mandibular muscle activation, and contact forces. Results: The maximum opening displacement for the envelope-based fossa prosthesis was greater than those for Biomet and ellipsoidal prostheses (36 mm, 35 mm, and 33 mm, respectively). The mandibular musculoskeletal model with ellipsoidal prosthesis led to dislocation near maximal jaw opening. Compared to Biomet, the envelope-based fossa reduced the digastric and lateral pterygoid activation at maximal jaw opening. It also reduced the maximal resistance to condylar sliding on the intact side by 63.2 N. Conclusion: A customized TMJ fossa prosthesis was successfully developed using the ESCM concept. Our study of musculoskeletal multibody modeling has highlighted its advantages and potential. The artificial fossa design successfully achieved a wider condylar range of motion. It also reduced the activation of jaw opening muscles on the affected side and resistance on the intact side. This study showed that an ESCM-based approach may be useful for optimizing TMJ fossa prostheses design.

Reverse dynamics analysis of contact force and muscle activities during the golf swing after total hip arthroplasty

There are no reports on hip kinetics including contact forces and muscle activities during the golf swing after total hip arthroplasty (THA). The aim of this study was to identify the characteristics of three-dimensional dynamics during the golf swing. Ten unilateral primary THA patients participated in motion capture test of their driver golf swing. The driver swing produced approximately 20-30° of rotation in both lead and trail replaced hips. The mean hip contact forces (HCFs) of lead and trail replaced hips were 5.1 and 6.6 × body weight, respectively. Left and right THAs showed similar HCFs of lead and trail hips. More than 60% of the Percent maximum voluntary isometric contraction was found in bilateral iliopsoas muscles in all unilateral THA. Three factors [female sex, lower modified Harris Hip Score, and higher HCF of surgical side] were associated with the golf-related replacement hip pain. Golf is an admissible sport after THA because driver swings do not contribute excessive rotation or contact forces to hip prostheses. HCF could be reduced through swing adjustments, which may allow patients with golf-related replacement hip pain to develop a comfortable golf game free from pain.

Two-screw osteosynthesis of the mandibular condylar head with different screw materials: a finite element analysis

This study compared the biomechanical behavior of titanium, magnesium, and polylactic acid screws for two-screw osteosynthesis of mandibular condylar head fractures using finite element analysis. Von Mises stress distribution, fracture displacement, and fragment deformation were evaluated. Titanium screws performed the best in terms of carrying the highest load, resulting in the least fracture displacement and fragment deformation. Magnesium screws showed intermediate results, while PLA screws were found to be unsuitable with stress values exceeding their tensile strength. These findings suggest that magnesium alloys could be considered a suitable alternative to titanium screws in mandibular condylar head osteosynthesis.

Two-versus three-screw osteosynthesis of the mandibular condylar head: A finite element analysis

Titanium screws are commonly used for osteosynthesis of mandibular condylar head fractures. Evidence suggests that the insertion of three screws may result in better fracture stability. Two screws only, on the other hand, could reduce adverse effects, mainly bone resorption. This study aimed to investigate the biomechanical differences in mandibular condylar head osteosynthesis with two versus three titanium screws using finite element analysis. A finite element model of the mandible with a right type P condylar head fracture fixed with two or three titanium screws was analyzed in ANSYS Mechanical. The geometry of the model assembly was constructed in ANSYS Spaceclaim. Biomechanical load boundary conditions were obtained from a validated musculoskeletal model in AnyBody Modeling System™. The preprocessing of the finite element model and mapping of the obtained boundary conditions was done in docq VIT. Fracture displacement, fragment deformation, von Mises stress distribution, and reaction forces within the screws were evaluated in ANSYS for three different loading scenarios. Finite element analysis showed similar results when comparing osteosynthesis with two versus three titanium screws for all three loading scenarios. Contralateral molar loading resulted in the highest stress on both the fracture and the screws with the maximum von Mises stress being found at the condylar neck. Stress concentration within the screws was found in the fracture gap and was higher in the lateral fragment. In all scenarios, maximum von Mises stress values were smaller when forces were distributed among three screws. However, stability was also adequate when two screws were used. Mandibular condylar head osteosynthesis with two titanium screws appears to provide sufficient fracture stability. Further clinical studies are needed to clarify the implications of these results.

Impact of Cleft Palate on Tongue-Based Upper Airway Obstruction in Pierre Robin Sequence: Implications for Mandibular Distraction Osteogenesis and Timing of Cleft Palate Repair

INTRODUCTION

The objectives of this study are to analyze the impact of cleft palate (CP) on upper airway obstruction using polysomnography in patients with Pierre Robin Sequence (PRS) undergoing mandibular distraction osteogenesis (MDO) and subsequent CP repair.

METHODS

A single-surgeon, retrospective chart review was performed of all patients with nonsyndromic PRS treated with MDO. Severity of upper airway obstruction was evaluated pre- and post-distraction via polysomnography. Details of MDO and CP repair were collected and any complications recorded.

RESULTS

Twenty-one nonsyndromic PRS patients with CP and 6 patients without CP met inclusion criteria. There was no significant difference in predistraction apnea-hypopnea index between the 2 groups (53.4 ± 42.1 versus 34.4 ± 18.9; P = 0.3). Patients with CP had significantly higher predistraction SpO2 saturation (94.5 ± 1.6% versus 91.0% ± 4.8%; P = 0.01), SpO2 nadir (74.5% ± 9.1% versus 63.6% ± 11.6%; P = 0.03) and lower percentage time spent below 90% SpO2 (6.0% ± 7.2% versus 23.6% ± 29.9%; P = 0.04). The rate of oronasal fistulas formation was 38%. The time between MDO and CP repair was shorter for patients with complications (250.8 ± 3.3 versus 370.8 ± 191.9 days; P = 0.08). Five patients experienced relapse of respiratory difficulties after CP repair.

CONCLUSIONS

The presence of CP in nonsyndromic PRS patients decreases the severity of obstructive sleep apnea by oxygen parameters on PSG. Palatal fistulas and relapse of respiratory distress are common complications of CP repair following MDO. Delaying CP repair may help to decrease complication rates.

Electromyography Based Validation of a Musculoskeletal Hand Model

Regarding the prevention of injuries and rehabilitation of the human hand, musculoskeletal simulations using an inverse dynamics approach allow for insights of the muscle recruitment and thus acting forces on the hand. Currently, several hand models from various research groups are in use, which are mainly validated by the comparison of numerical and anatomical moment arms. In contrast to this validation and model-building technique by cadaver studies, the aim of the present study is to further validate a recently published hand model [1] by analyzing numerically calculated muscle activities in comparison to experimentally measured electromyographical signals of the muscles. Therefore, the electromyographical signals of 10 hand muscles of five test subjects performing seven different hand movements were measured. The kinematics of these tasks were used as input for the hand model, and the numerical muscle activities were computed. To analyze the relationship between simulated and measured activities, the time difference of the muscle on- and off-set points were calculated, which resulted in a mean on- and off-set time difference of 0.58 s between the experimental data and the model. The largest differences were detected for movements that mainly addressed the wrist. One major issue comparing simulated and measured muscle activities of the hand is cross-talk. Nevertheless, the results show that the hand model fits the experiment quite accurately despite some limitations and is a further step towards patient-specific modelling of the upper extremity.

3D Printing Experimental Validation of the Finite Element Analysis of the Maxillofacial Model

Contacts used in finite element (FE) models were considered as the best simulation for interactions in the temporomandibular joint (TMJ). However, the precision of simulations should be validated through experiments. Three-dimensional (3D) printing models with the high geometric and loading similarities of the individuals were used in the validation. This study aimed to validate the FE models of the TMJ using 3D printing models. Five asymptomatic subjects were recruited in this study. 3D models of mandible, disc, and maxilla were reconstructed according to cone-beam CT (CBCT) image data. PLA was chosen for 3D printing models from bottom to top. Five pressure forces corresponding to the central occlusion were applied to the 3D printing models. Ten strain rosettes were distributed on the mandible to record the horizontal and vertical strains. Contact was used in the FE models with the same geometries, material properties, loadings, and boundary conditions as 3D printing models to simulate the interaction of the disc-condyle, disc-temporal bone, and upper-lower dentition. The differences of the simulated and experimental results for each sample were less than 5% (maximum 4.92%) under all five loadings. In conclusion, it was accurate to use contact to simulate the interactions in TMJs and upper-lower dentition.

Time-dependent in silico modelling of orthognathic surgery to support the design of biodegradable bone plates

Orthognathic surgery is performed to realign the jaws of a patient through several osteotomies. The state-of-the-art bone plates used to maintain the bone fragments in place are made of titanium. The presence of these non-degradable plates can have unwanted side effects on the long term (e.g. higher infection risk) if they are not removed. Using a biodegradable material such as magnesium may be a possible solution to this problem. However, biodegradation leads to a decrease of mechanical strength, therefore a time-dependent computational approach can help to evaluate the performance of such plates. In the present work, a computational framework has been developed to include biodegradation and bone healing algorithms in a finite element model. It includes bone plates and the mandible, which are submitted to masticatory loads during the early healing period (two months following the surgery). Two different bone plate designs with different stiffnesses have been tested. The stiff design exhibited good mechanical stability, with maximum Von Mises stress being less than 40% of the yield strength throughout the simulation. The flexible design shows high stresses when the bone healing has not started in the fracture gaps, indicating possible failure of the plate. However, this design leads to a higher bone healing quality after two months, as more cartilage is formed due to higher strains exerted in fracture gaps. We therefore conclude that in silico modelling can support tuning of the design parameters to ensure mechanical stability and while promoting bone healing.

A Biomechanical Analysis of Muscle Force Changes After Bilateral Sagittal Split Osteotomy

A basic procedure affecting maxillofacial geometry is the bilateral sagittal split osteotomy. During the surgery, the bony segments are placed in a new position that provides the correct occlusion. Changes in the geometry of the mandible will affect the surrounding structures and will have a significant impact on the functioning of the masticatory system. As a result of the displacement of the bone segment, the biomechanical conditions change, i.e., the load and the position of the muscles. The primary aim of this study was to determine the changes in the values of the muscular forces caused by mandible geometry alteration. The study considered the translation and rotation of the distal segment, as well as rotations of the proximal segments in three axes. Calculations were performed for the unilateral, static loading of a model based on rigid body mechanics. Muscles were modeled as spring elements, and a novel approach was used to determine muscle stiffness. In addition, an attempt was made, based on the results obtained for single displacements separately, to determine the changes in muscle forces for geometries with complex displacements. Based on the analysis of the results, it was shown that changes in the geometry of the mandibular bone associated with the bilateral sagittal split osteotomy will have a significant effect on the values of the masticatory muscle forces. Displacement of the distal segment has the greatest effect from -21.69 to 26.11%, while the proximal segment rotations affected muscle force values to a less extent, rarely exceeding 1%. For Yaw and Pitch rotations, the opposite effect of changes within one muscle is noticed. Changes in muscle forces for complex geometry changes can be determined with a high degree of accuracy by the appropriate summation of results obtained for simple cases.

A new musculoskeletal AnyBody™ detailed hand model

Musculoskeletal research questions regarding the prevention or rehabilitation of the hand can be addressed using inverse dynamics simulations when experiments are not possible. To date, no complete human hand model implemented in a holistic human body model has been fully developed. The aim of this work was to develop, implement, and validate a fully detailed hand model using the AnyBody Modelling System (AMS) (AnyBody, Aalborg, Denmark). To achieve this, a consistent multiple cadaver dataset, including all extrinsic and intrinsic muscles, served as a basis. Various obstacle methods were implemented to obtain with the correct alignment of the muscle paths together with the full range of motion of the fingers. These included tori, cylinders, and spherical ellipsoids. The origin points of the lumbrical muscles within the tendon of the flexor digitorum profundus added a unique feature to the model. Furthermore, the possibility of an entire patient-specific scaling based on the hand length and width were implemented in the model. For model validation, experimental datasets from the literature were used, which included the comparison of numerically calculated moment arms of the wrist, thumb, and index finger muscles. In general, the results displayed good comparability of the model and experimental data. However, the extrinsic muscles showed higher accordance than the intrinsic ones. Nevertheless, the results showed, that the proposed developed inverse dynamics hand model offers opportunities in a broad field of applications, where the muscles and joint forces of the forearm play a crucial role.

Effect of suprahyoid muscles on mouth opening with and without lateral pterygoid muscle: 3D inverse dynamic model analysis

Objective: To evaluate the contribution of suprahyoid muscles to mouth opening in different hyoid bone positions.Methods: The jaw-opening and closing movements were imposed on the 3D inverse dynamic jaw model with and without the lateral pterygoid muscle (LPM). The activation of the muscles was evaluated for eight different positions of the hyoid bone.Results: The muscles with 100% activation provided maximum mouth opening (MMO). When the hyoid bone was replaced, the muscles could achieve MMO at the 135º, 180º, and 225º points with the LPM. Mouth opening was estimated to be 36.8 mm in the absence of the LPM. A jaw opening greater than 36.8 mm was seen when the hyoid bone was repositioned at the 90º, 180º, and 135º points.Discussion: The contribution of suprahyoid muscles to mouth opening varies in different hyoid bone positions, with the inferior and anterior positions having a positive impact.

Muscle Contributions to Upper-Extremity Movement and Work From a Musculoskeletal Model of the Human Shoulder

Musculoskeletal models enable movement scientists to examine muscle function by computing the mechanical work done by muscles during motor tasks. To estimate muscle work accurately requires a model that is physiologically plausible. Previous models of the human shoulder have coupled scapula movement to humeral movement. While coupled movement produces a stereotypical scapulohumeral rhythm, it cannot model shrugging or independent movement of the scapula and humerus. The artificial coupling of humeral elevation to scapular rotation permits muscles that cross the glenohumeral joint, such as the rotator-cuff muscles and deltoids, to do implausible work to elevate and rotate the scapula. In reality, the motion of the scapula is controlled by thoracoscapular muscles, yet the roles of these muscles in shoulder function remains unclear. To elucidate the roles of the thoracoscapular muscles, we developed a shoulder model with an accurate scapulothoracic joint and includes scapular muscles to drive its motion. We used the model to compute the work done by the thoracoscapular muscles during shrugging and arm elevation. We found that the bulk of the work done in upper-extremity tasks is performed by the largest muscles of the shoulder: trapezius, deltoids, pectoralis major, and serratus-anterior. Trapezius and serratus anterior prove to be important synergists in performing upward-rotation of the scapula. We show that the large thoracoscapular muscles do more work than glenohumeral muscles during arm-elevation tasks. The model, experimental data and simulation results are freely available on SimTK.org to enable anyone to explore our results and to perform further studies in OpenSim 4.0.

A Dynamic Jaw Model With a Finite-Element Temporomandibular Joint

The masticatory region is an important human motion system that is essential for basic human tasks like mastication, speech or swallowing. An association between temporomandibular disorders (TMDs) and high temporomandibular joint (TMJ) stress has been suggested, but in vivo joint force measurements are not feasible to directly test this assumption. Consequently, biomechanical computer simulation remains as one of a few means to investigate this complex system. To thoroughly examine orofacial biomechanics, we developed a novel, dynamic computer model of the masticatory system. The model combines a muscle driven rigid body model of the jaw region with a detailed finite element model (FEM) disk and elastic foundation (EF) articular cartilage. The model is validated using high-resolution MRI data for protrusion and opening that were collected from the same volunteer. Joint stresses for a clenching task as well as protrusive and opening movements are computed. Simulations resulted in mandibular positions as well as disk positions and shapes that agree well with the MRI data. The model computes reasonable disk stress patterns for dynamic tasks. Moreover, to the best of our knowledge this model presents the first ever contact model using a combination of EF layers and a FEM body, which results in a clear decrease in computation time. In conclusion, the presented model is a valuable tool for the investigation of the human TMJ and can potentially help in the future to increase the understanding of the masticatory system and the relationship between TMD and joint stress and to highlight potential therapeutic approaches for the restoration of orofacial function.

Muscle-tendon unit parameter estimation of a Hill-type musculoskeletal model based on experimentally obtained subject-specific torque profiles

The aim of this study was to generate a subject-specific musculoskeletal muscle model, based on isometric and isovelocity measurements of the whole lower extremity. A two-step optimisation procedure is presented for optimising the muscle-tendon parameters for isometric and isovelocity joint torque profiles. A significant decrease in the differences between the experimental and model-predicted joint torques was obtained for both isometric and isovelocity cases. However, only small to moderate differences between an isometric based optimisation and an isovelocity-based optimisation were observed. Depending on the specific purpose of the model, it may be beneficial to use an isometric-only rather than isovelocity-based optimisation.

In vivo prediction of temporomandibular joint disc thickness and position changes for different jaw positions

Temporomandibular joint disorders (TMD) are common dysfunctions of the masticatory region and are often linked to dislocation or changes of the temporomandibular joint (TMJ) disc. Magnetic resonance imaging (MRI) is the gold standard for TMJ imaging but standard clinical sequences do not deliver a sufficient resolution and contrast for the creation of detailed meshes of the TMJ disc. Additionally, bony structures cannot be captured appropriately using standard MRI sequences due to their low signal intensity. The objective of this study was to enable researchers to create high resolution representations of all structures of the TMJ and consequently investigate morphological as well as positional changes of the masticatory system. To create meshes of the bony structures, a single computed tomography (CT) scan was acquired. In addition, a high-resolution MRI sequence was produced, which is used to collect the thickness and position change of the disc for various static postures using bite blocks. Changes in thickness of the TMJ disc as well as disc translation were measured. The newly developed workflow successfully allows researchers to create high resolution models of all structures of the TMJ for various static positions, enabling the investigation of TMJ disc translation and deformation. Discs were thinnest in the lateral part and moved mainly anteriorly and slightly medially. The procedure offers the most comprehensive picture of disc positioning and thickness changes reported to date. The presented data can be used for the development of a biomechanical computer model of TMJ anatomy and to investigate dynamic and static loads on the components of the system, which could be useful for the prediction of TMD onset.

Agreement of measured and calculated muscle activity during highly dynamic movements modelled with a spherical knee joint

The inclusion of muscle forces into the analysis of joint contact forces has improved their accuracy. But it has not been validated if such force and activity calculations are valid during highly dynamic multidirectional movements. The purpose of this study was to validate calculated muscle activation of a lower extremity model with a spherical knee joint for running, sprinting and 90°-cutting. Kinematics, kinetics and lower limb muscle activation of ten participants were investigated in a 3D motion capture setup including EMG. A lower extremity rigid body model was used to calculate the activation of these muscles with an inverse dynamics approach and a cubic cost function. Correlation coefficients were calculated to compare measured and calculated activation. The results showed good correlation of the modelled and calculated data with a few exceptions. The highest average correlations were found during walking (r = 0.81) and the lowest during cutting (r = 0.57). Tibialis anterior had the lowest average correlation (r = 0.33) over all movements while gastrocnemius medius had the highest correlation (r = 0.9). The implementation of a spherical knee joint increased the agreement between measured and modelled activation compared to studies using a hinge joint knee. Although some stabilizing muscles showed low correlations during dynamic movements, the investigated model calculates muscle activity sufficiently.

Ex-vivo and in vitro validation of an innovative mandibular condyle implant concept

PURPOSE

The purpose of this study is to pre-validate an inovative implant concept, and to compare the behavior of the mandibular condyle against a commercial Biomet implant in an ex vivo model and present results of the first cadaveric studies.

MATERIALS AND METHODS

Three experimental cadaveric condyles were tested under three conditions: one intact, another with the Biomet model, and one with the innovative concept. The condyle was tested with a reaction of 300 N in all situations and the principal strains were measured. Before the geometry of the cadaveric condyle was reconstructed from a microCT scan, and a finite element model was created. Finally, a procedure was carried out with the new implant by two expert surgeons on a two cadaveric head model.

RESULTS

In vitro the mandible condyle presents a linear behavior until maximum load. The strain measured with Biomet implant indicates a strain shielding effect in the proximal region, inducing bone loss in the long term. The lingual side of the Biomet implanted condyle presents an increase of +44% in strain.

CONCLUSION

The new concept was evaluated and showed a similar behavior to the intact model, and better behavior than the Biomet. The innovative concept proves that it is possible to avoid screws for a TMJ fixation and improve the TMJ alloplastic behavior.

Construction and evaluation of a model for wheelchair propulsion in an individual with tetraplegia

Upper limb overuse injuries are common in manual wheelchair users with spinal cord injury. Patient-specific in silico models enhance experimental biomechanical analyses by estimating in vivo shoulder muscle and joint contact forces. Current models exclude deep shoulder muscles that have important roles in wheelchair propulsion. Freely accessible patient-specific models have not been generated for persons with tetraplegia, who have a greater risk for shoulder pain and injury. The objectives of this work were to (i) construct a freely accessible, in silico, musculoskeletal model capable of generating patient-specific dynamic simulations of wheelchair propulsion and (ii) establish proof-of-concept with data obtained from an individual with tetraplegia. Constructed with OpenSim, the model features muscles excluded in existing models. Shoulder muscle forces and activations were estimated via inverse dynamics. Mean absolute error of estimated muscle activations and fine-wire electromyography (EMG) recordings was computed. Mean muscle activation for five consecutive stroke cycles demonstrated good correlation (0.15-0.17) with fine-wire EMG. These findings, comparable to other studies, suggest that the model is capable of estimating shoulder muscle forces during wheelchair propulsion. The additional muscles may provide a greater understanding of shoulder muscle contribution to wheelchair propulsion. The model may ultimately serve as a powerful clinical tool. Graphical abstract ᅟ.

Three-dimensional temporomandibular joint muscle attachment morphometry and its impacts on musculoskeletal modeling

In musculoskeletal models of the human temporomandibular joint (TMJ), muscles are typically represented by force vectors that connect approximate muscle origin and insertion centroids (centroid-to-centroid force vectors). This simplification assumes equivalent moment arms and muscle lengths for all fibers within a muscle even with complex geometry and may result in inaccurate estimations of muscle force and joint loading. The objectives of this study were to quantify the three-dimensional (3D) human TMJ muscle attachment morphometry and examine its impact on TMJ mechanics. 3D muscle attachment surfaces of temporalis, masseter, lateral pterygoid, and medial pterygoid muscles of human cadaveric heads were generated by co-registering measured attachment boundaries with underlying skull models created from cone-beam computerized tomography (CBCT) images. A bounding box technique was used to quantify 3D muscle attachment size, shape, location, and orientation. Musculoskeletal models of the mandible were then developed and validated to assess the impact of 3D muscle attachment morphometry on joint loading during jaw maximal open-close. The 3D morphometry revealed that muscle lengths and moment arms of temporalis and masseter muscles varied substantially among muscle fibers. The values calculated from the centroid-to-centroid model were significantly different from those calculated using the 'Distributed model', which considered crucial 3D muscle attachment morphometry. Consequently, joint loading was underestimated by more than 50% in the centroid-to-centroid model. Therefore, it is necessary to consider 3D muscle attachment morphometry, especially for muscles with broad attachments, in TMJ musculoskeletal models to precisely quantify the joint mechanical environment critical for understanding TMJ function and mechanobiology.

Estimation of muscle activation during different walking speeds with two mathematical approaches compared to surface EMG

BACKGROUND

Muscle force estimation could improve clinical gait analysis by enhancing insight into causes of impairments and informing targeted treatments. However, it is not currently standard practice to use muscle force models to augment clinical gait analysis, partly, because robust validations of estimated muscle activations, underpinning force modelling processes, against recorded electromyography (EMG) are lacking.

RESEARCH QUESTION

Therefore, in order to facilitate future clinical use, this study sought to validate estimated lower limb muscle activation using two mathematical models (static optimisation SO, computed muscle control CMC) against recorded muscle activations of ten healthy participants.

METHODS

Participants walked at five speeds. Visual agreement in activation onset and offset as well as linear correlation (r) and mean absolute error (MAE) between models and EMG were evaluated.

RESULTS

MAE between measured and recorded activations were variable across speeds (SO vs EMG 15-68%, CMC vs EMG 13-69%). Slower speeds resulted in smaller deviations (mean MAE < 30%) than faster speeds. Correlation was high (r > 0.5) for only 11/40 (CMC) and 6/40 (SO) conditions (muscles X speeds) compared to EMG.

SIGNIFICANCE

Modelling approaches do not yet show sufficient consistency of agreement between estimated and recorded muscle activation to support recommending immediate clinical adoption of muscle force modelling. This may be because assumptions underlying muscle activation estimations (e.g. muscles' anatomy and maximum voluntary contraction) are not yet sufficiently individualizable. Future research needs to find timely and cost efficient ways to scale musculoskeletal models for better individualisation to facilitate future clinical implementation.

An Investigation of Two Finite Element Modeling Solutions for Biomechanical Simulation Using a Case Study of a Mandibular Bone

The method used in biomechanical modeling for finite element method (FEM) analysis needs to deliver accurate results. There are currently two solutions used in FEM modeling for biomedical model of human bone from computerized tomography (CT) images: one is based on a triangular mesh and the other is based on the parametric surface model and is more popular in practice. The outline and modeling procedures for the two solutions are compared and analyzed. Using a mandibular bone as an example, several key modeling steps are then discussed in detail, and the FEM calculation was conducted. Numerical calculation results based on the models derived from the two methods, including stress, strain, and displacement, are compared and evaluated in relation to accuracy and validity. Moreover, a comprehensive comparison of the two solutions is listed. The parametric surface based method is more helpful when using powerful design tools in computer-aided design (CAD) software, but the triangular mesh based method is more robust and efficient.

Risk of subacromial shoulder disorder in airport baggage handlers: combining duration and intensity of musculoskeletal shoulder loads

Musculoskeletal shoulder load among baggage handlers measured by combining duration and intensity based on biomechanical and epidemiological information may be a stronger predictor of subacromial shoulder disorders than baggage handler seniority. In 2012, a cohort of baggage handlers employed at Copenhagen Airport in 1990-2012, and a cohort of unskilled otherwise employed men answered a survey. Self-reported information on work tasks during employment in the airport in combination with work task specific biomechanically modelled forces in the shoulder joint was used to estimate shoulder load. Exposure measures were accumulated shoulder abduction moment, accumulated shoulder compression force, accumulated supraspinatus force and baggage handler seniority. The outcome was subacromial shoulder disorder registered in the Danish National Patient Register. When analyses were adjusted by all confounders except age, exposure variables showed close to significant associations with subacromial shoulder disorder. Results could not confirm our hypothesis that combined information on work task duration and shoulder load intensity was stronger associated with subacromial shoulder disorder than seniority. Practitioner Summary: In this study we sought to identify if the exposure to work-related musculoskeletal shoulder loading including duration and intensity among baggage handlers was associated with subacromial shoulder disorder. We found that there was an association but this was not stronger than that between baggage handler seniority and subacromial shoulder disorder.

Occlusal loading during biting from an experimental and simulation point of view

OBJECTIVES

Occlusal loading during clenching and biting is achieved by the action of the masticatory system, and forms the basis for the evaluation of the functional performance of prosthodontic and maxillofacial components. This review provides an overview of (i) current bite force measurement techniques and their limitations and (ii) the use of computational modelling to predict bite force. A brief simulation study highlighting the challenges of current computational dental models is also presented.

METHODS

Appropriate studies were used to highlight the development and current bite force measurement methodologies and state-of-the-art simulation for computing bite forces using biomechanical models.

RESULTS

While a number of strategies have been developed to measure occlusal forces in three-dimensions, the use of strain-gauges, piezo-electric sensors and pressure sheets remain the most widespread. In addition to experimental-based measurement techniques, bite force may be also estimated using computational models of the masticatory system. Simulations of different bite force models clearly show that the use of three-dimensional force measurements enriches the evaluation of masticatory functional performance.

SIGNIFICANCE

Hence, combining computational modelling with three-dimensional force measurement techniques can significantly improve the evaluation of masticatory system and the functional performance of prosthodontic components.

Gait training assist system of a lower limb prosthetic visualizing muscle activation pattern using a color-depth sensor

Some unilateral lower-limb amputees load the intact limb more than the prosthetic limb. This can cause chronic pains, fatigue, lumbago, and joint diseases, including knee osteoarthritis. To avoid and counteract these symptoms it is necessary to improve their asymmetric gait. Increasing the function of the hip abductor muscle is important to maintaining symmetrical weight distribution. Therefore, the purpose of this study is to develop a training assist system, which estimates and visualizes an abductor muscle by using a color-depth sensor. To estimate the muscle activation, first, the floor reaction force is calculated using a simple dynamic model. Then, the hip torque is calculated using joint angles. The floor reaction force and, the muscle length are calculated based on a human musculoskeletal model. Muscle activity is estimated by these parameters. Evaluation experiments of this proposed method were performed on healthy persons and unilateral trans femoral amputees, and the effectiveness of this proposed algorithm has been confirmed.

Introduction to Force-Dependent Kinematics: Theory and Application to Mandible Modeling

Knowledge of the muscle, ligament, and joint forces is important when planning orthopedic surgeries. Since these quantities cannot be measured in vivo under normal circumstances, the best alternative is to estimate them using musculoskeletal models. These models typically assume idealized joints, which are sufficient for general investigations but insufficient if the joint in focus is far from an idealized joint. The purpose of this study was to provide the mathematical details of a novel musculoskeletal modeling approach, called force-dependent kinematics (FDK), capable of simultaneously computing muscle, ligament, and joint forces as well as internal joint displacements governed by contact surfaces and ligament structures. The method was implemented into the anybody modeling system and used to develop a subject-specific mandible model, which was compared to a point-on-plane (POP) model and validated against joint kinematics measured with a custom-built brace during unloaded emulated chewing, open and close, and protrusion movements. Generally, both joint models estimated the joint kinematics well with the POP model performing slightly better (root-mean-square-deviation (RMSD) of less than 0.75 mm for the POP model and 1.7 mm for the FDK model). However, substantial differences were observed when comparing the estimated joint forces (RMSD up to 24.7 N), demonstrating the dependency on the joint model. Although the presented mandible model still contains room for improvements, this study shows the capabilities of the FDK methodology for creating joint models that take the geometry and joint elasticity into account.

The influence of the human TMJ eminence inclination on predicted masticatory muscle forces

Aim of this paper was to investigate the change in masticatory muscle forces and temporomandibular joint (TMJ) reaction forces simulated by inverse dynamics when thesteepness of the anterior fossa slope was varied. We used the model by de Zee et al. (2007) created in AnyBody™. The model was equipped with 24musculotendon actuators. Mandibular movement was governed by thetrajectory of theincisal point. The TMJ was modelled as a planar constraint canted 5°medially and thecaudal inclination relative to the occlusal plane was varied from 10° to 70°. Our models showed that for the two simulated movements (empty chewing and unilateral clenching) the joint reaction forces were smallest for the eminence inclination of 30° and 40° and highest for 70°. The muscle forces were relatively insensitive to change of the eminence inclination for the angles between 20° and 50°. This did not hold for the pterygoid muscle, for which the muscle forces increased continually with increasing fossa inclination. For empty chewing the muscle force reached smaller values than for clenching. During clenching, the muscle forces changed by up to 200N.

Metatarsal Loading During Gait—A Musculoskeletal Analysis

Detailed knowledge of the loading conditions within the human body is essential for the development and optimization of treatments for disorders and injuries of the musculoskeletal system. While loads in the major joints of the lower limb have been the subject of extensive study, relatively little is known about the forces applied to the individual bones of the foot. The objective of this study was to use a detailed musculoskeletal model to compute the loads applied to the metatarsal bones during gait across several healthy subjects. Motion-captured gait trials and computed tomography (CT) foot scans from four healthy subjects were used as the inputs to inverse dynamic simulations that allowed the computation of loads at the metatarsal joints. Low loads in the metatarsophalangeal (MTP) joint were predicted before terminal stance, however, increased to an average peak of 1.9 times body weight (BW) before toe-off in the first metatarsal. At the first tarsometatarsal (TMT) joint, loads of up to 1.0 times BW were seen during the early part of stance, reflecting tension in the ligaments and muscles. These loads subsequently increased to an average peak of 3.0 times BW. Loads in the first ray were higher compared to rays 2–5. The joints were primarily loaded in the longitudinal direction of the bone.

Masticatory biomechanics in the rabbit: a multi-body dynamics analysis

Multi-body dynamics is a powerful engineering tool which is becoming increasingly popular for the simulation and analysis of skull biomechanics. This paper presents the first application of multi-body dynamics to analyse the biomechanics of the rabbit skull. A model has been constructed through the combination of manual dissection and three-dimensional imaging techniques (magnetic resonance imaging and micro-computed tomography). Individual muscles are represented with multiple layers, thus more accurately modelling muscle fibres with complex lines of action. Model validity was sought through comparing experimentally measured maximum incisor bite forces with those predicted by the model. Simulations of molar biting highlighted the ability of the masticatory system to alter recruitment of two muscle groups, in order to generate shearing or crushing movements. Molar shearing is capable of processing a food bolus in all three orthogonal directions, whereas molar crushing and incisor biting are predominately directed vertically. Simulations also show that the masticatory system is adapted to process foods through several cycles with low muscle activations, presumably in order to prevent rapidly fatiguing fast fibres during repeated chewing cycles. Our study demonstrates the usefulness of a validated multi-body dynamics model for investigating feeding biomechanics in the rabbit, and shows the potential for complementing and eventually reducing in vivo experiments.

Finite element analysis of the human mastication cycle

The aim of this paper is to propose a biomechanical model that could serve as a tool to overcome some difficulties encountered in experimental studies of the mandible. One of these difficulties is the inaccessibility of the temporomandibular joint (TMJ) and the lateral pterygoid muscle. The focus of this model is to study the stresses in the joint and the influence of the lateral pterygoid muscle on the mandible movement. A finite element model of the mandible, including the TMJ, was built to simulate the process of unilateral mastication. Different activation patterns of the left and right pterygoid muscles were tried. The maximum stresses in the articular disc and in the whole mandible during a complete mastication cycle were reached during the instant of centric occlusion. The simulations show a great influence of the coordination of the right and left lateral pterygoid muscles on the movement of the jaw during mastication. An asynchronous activation of the lateral pterygoid muscles is needed to achieve a normal movement of the jaw during mastication.

Modeling the biomechanics of swine mastication--an inverse dynamics approach

A novel reconstructive alternative for patients with severe facial structural deformity is Le Fort-based, face-jaw-teeth transplantation (FJTT). To date, however, only ten surgeries have included underlying skeletal and jaw-teeth components, all yielding sub-optimal results and a need for a subsequent revision surgery, due to size mismatch and lack of precise planning. Numerous studies have proven swine to be appropriate candidates for translational studies including pre-operative planning of transplantation. An important aspect of planning FJTT is determining the optimal muscle attachment sites on the recipient's jaw, which requires a clear understanding of mastication and bite mechanics in relation to the new donated upper and/or lower jaw. A segmented CT scan coupled with data taken from literature defined a biomechanical model of mandible and jaw muscles of a swine. The model was driven using tracked motion and external force data of one cycle of chewing published earlier, and predicted the muscle activation patterns as well as temporomandibular joint (TMJ) reaction forces and condylar motions. Two methods, polynomial and min/max optimization, were used for solving the muscle recruitment problem. Similar performances were observed between the two methods. On average, there was a mean absolute error (MAE) of <0.08 between the predicted and measured activation levels of all muscles, and an MAE of <7 N for TMJ reaction forces. Simulated activations qualitatively followed the same patterns as the reference data and there was very good agreement for simulated TMJ forces. The polynomial optimization produced a smoother output, suggesting that it is more suitable for studying such motions. Average MAE for condylar motion was 1.2mm, which reduced to 0.37 mm when the input incisor motion was scaled to reflect the possible size mismatch between the current and original swine models. Results support the hypothesis that the model can be used for planning of facial transplantation.

Strengths and limitations of a musculoskeletal model for an analysis of simulated meat cutting tasks

This study assessed the capacity of a musculoskeletal model to predict the relative muscle activation changes as a function of the workbench height and the movement direction during a simulated meat cutting task. Seven subjects performed a cutting task alternating two cutting directions for 20 s at four different workbench heights. Kinematics, electromyography (EMG), and cutting force data were collected and used to drive a musculoskeletal model of the shoulder girdle. The model predicted the muscle forces exerted during the task. Both the recorded and computed activation of the muscles was then compared by means of cross-correlation and by comparison of muscle activation trends with respect to the workstation parameters, i.e. cutting direction and workbench height. The results indicated that cutting movements involving arm flexion are preferable to movement requiring internal arm rotation and abduction. The optimal bench height for meat cutting tasks should be between 20 and 30 cm below the worker's elbow height. The present study underlines a beneficial use of musculoskeletal models for adjusting workstation parameters.

Muscle load in reaching movements performed by a wheelchair user: a case study

PURPOSE

The aim of this study was to analyse the load on the shoulder muscles during reaching movements that are specific to wheelchair users in relation to the risk of impingement.

METHOD

Three activities of daily living were performed: putting a book on a shelf in front and at the side and putting a pack of water bottles on a table. The AnyBody shoulder model was used to calculate the activity and forces of the shoulder muscles.

RESULTS

Handling the pack of bottles caused the highest forces in the deltoideus, trapezius, serratus anterior and rotator cuff muscles. For handling the book, the highest forces were found in the deltoideus (scapular part) and the serratus anterior, especially during the put phase.

CONCLUSIONS

Handling heavy objects such as a pack of bottles or a wheelchair produces high forces on the rotator cuff muscles and can lead to early fatigue. Therefore, these activities seem to be associated with a high risk of developing impingement syndrome. Implications for Rehabilitation In a single patient, this study demonstrates that the load on the rotator cuff is high during reaching movements. Handling a pack of water bottles, which resembles wheelchair handling, represents an activity associated with a high risk of developing impingement syndrome. Shoulder muscles must be trained in a balanced way to provide stabilization at the shoulder joint and prevent fatigue.

The morphology of the mouse masticatory musculature

The mouse has been the dominant model organism in studies on the development, genetics and evolution of the mammalian skull and associated soft-tissue for decades. There is the potential to take advantage of this well studied model and the range of mutant, knockin and knockout organisms with diverse craniofacial phenotypes to investigate the functional significance of variation and the role of mechanical forces on the development of the integrated craniofacial skeleton and musculature by using computational mechanical modelling methods (e.g. finite element and multibody dynamic modelling). Currently, there are no detailed published data of the mouse masticatory musculature available. Here, using a combination of micro-dissection and non-invasive segmentation of iodine-enhanced micro-computed tomography, we document the anatomy, architecture and proportions of the mouse masticatory muscles. We report on the superficial masseter (muscle, tendon and pars reflecta), deep masseter, zygomaticomandibularis (anterior, posterior, infraorbital and tendinous parts), temporalis (lateral and medial parts), external and internal pterygoid muscles. Additionally, we report a lateral expansion of the attachment of the temporalis onto the zygomatic arch, which may play a role in stabilising this bone during downwards loading. The data presented in this paper now provide a detailed reference for phenotypic comparison in mouse models and allow the mouse to be used as a model organism in biomechanical and functional modelling and simulation studies of the craniofacial skeleton and particularly the masticatory system.

The importance of accurate muscle modelling for biomechanical analyses: a case study with a lizard skull

Computer-based simulation techniques such as multi-body dynamics analysis are becoming increasingly popular in the field of skull mechanics. Multi-body models can be used for studying the relationships between skull architecture, muscle morphology and feeding performance. However, to be confident in the modelling results, models need to be validated against experimental data, and the effects of uncertainties or inaccuracies in the chosen model attributes need to be assessed with sensitivity analyses. Here, we compare the bite forces predicted by a multi-body model of a lizard (Tupinambis merianae) with in vivo measurements, using anatomical data collected from the same specimen. This subject-specific model predicts bite forces that are very close to the in vivo measurements and also shows a consistent increase in bite force as the bite position is moved posteriorly on the jaw. However, the model is very sensitive to changes in muscle attributes such as fibre length, intrinsic muscle strength and force orientation, with bite force predictions varying considerably when these three variables are altered. We conclude that accurate muscle measurements are crucial to building realistic multi-body models and that subject-specific data should be used whenever possible.

Measuring functional outcome after total hip replacement with subject-specific hip joint loading

Total hip replacement is an often-performed orthopedic surgical procedure; the amount of procedures undertaken will increase since our life expectancy is growing. In order to optimize function, hip biomechanics should be restored to as near normal as possible. The goal of this pilot study was to determine whether or not it is feasible to compute the vectorial hip reaction force pathways on the head of the prosthesis and the force angles relative to the cup of the prosthesis that occur during gait in total hip replacement patients, serving as an objective measurement of the functional outcome following hip replacement. A three-dimensional gait analysis, measuring ground reaction forces and kinematics, was performed. The data retrieved from the gait analysis was used as the input for the musculoskeletal model to compute vectorial joint reaction forces for data processing. To evaluate the position and orientation of the joint reaction forces, the force path, as well as the force angles for the operated and non-operated joint, has been calculated during the stance phase of the specific leg. The force path for subject 2 on the non-operated side is only located in the posterior-lateral quarter, as is the force path for subject 1. In contrast to this subject, the force path for subject 2 at the operated hip joint can be found only within the anterior quarter of the head of the implant, where it is nearly equally distributed in the medio-lateral half of the prosthesis head. The force-inclination angles on the cup of subject 1, with respect to the plane of the socket face, indicates that the force vector is mainly positioned in the same quadrant when compared with subject 2 (in a cup-fixed coordinate system). The force-anteversion angle behaves similarly to the force-inclination angle, even when the effects are not as pronounced. The proposed methods in this article are aiming to define two functional outcomes of total hip replacement that are related to wear and rim loading. It is accepted that wear is not only a function of time, but a function of use. Owing to the methods listed in this article, we are able to determine a) the applied force and b) the sliding distance (force pathway) in a subject-specific manner. The computed hip-reaction force angles and the distance to the rim cup are a measurement for cup or rim loading, and occurs in the so-called safe-zones. This method may well give us insight into the biomechanical situation during gait, after receiving total hip replacement, that we need to fully understand the mechanisms acting on a hip joint and to prove a possible increase of functional outcome after receiving total hip replacement.

On validation of multibody musculoskeletal models

We review the opportunities to validate multibody musculoskeletal models in view of the current transition of musculoskeletal modelling from a research topic to a practical simulation tool in product design, healthcare and other important applications. This transition creates a new need for justification that the models are adequate representations of the systems they simulate. The need for a consistent terminology and established standards is identified and knowledge from fields with a more progressed state-of-the-art in verification and validation is introduced. A number of practical steps for improvement of the validation of multibody musculoskeletal models are pointed out and directions for future research in the field are proposed. It is hoped that a more structured approach to model validation can help to improve the credibility of musculoskeletal models.

Biomechanical characteristics of regenerated cortical bone in the canine mandible

To test the mechanical properties of regenerate cortical bone created using mandibular bone transport (MBT) distraction, five adult male American foxhound dogs underwent unilateral distraction of the mandible with a novel MBT device placed to linearly repair a 30-35 mm bone defect. The animals were sacrificed 12 weeks after the beginning of the consolidation period. Fourteen cylindrical specimens were taken from the inner (lingual) and outer (buccal) plates of the reconstructed mandible and 21 control specimens were removed from the contralateral aspect of the mandible. The mechanical properties of the 35 cylindrical cortical bone specimens were assessed by using a non-destructive pulse ultrasound technique. Results showed that all of the cortical mechanical properties exhibit higher numerical values on the control side than the MBT regenerate side. In addition, both densities and the elastic moduli in the direction of maximum stiffness of the regenerate cortical bone specimens are higher on the lingual side than the buccal side. Interestingly, there is no statistical difference between elastic modulus (E(1) and E(2) ) in orthogonal directions throughout the 35 cortical specimens. The data suggest that not only is the regenerate canine cortical bone heterogeneous, but the elastic mechanical properties tend to approximate transverse isotropy at a tissue level, as opposed to control cortical bone, which is orthotropic. In addition, the elastic mechanical properties are higher not only on the control side but also in the lingual anatomical position, suggesting a stress shielding effect from the presence of the reconstruction plate.

Craniofacial biomechanics: an overview of recent multibody modelling studies

Multibody modelling is underutilised in craniofacial analyses, particularly when compared to other computational methods such as finite element analysis. However, there are many potential applications within this area, where bony movements, muscle forces, joint kinematics and bite forces can all be studied. This paper provides an overview of recent, three-dimensional, multibody modelling studies related to the analysis of skulls. The goal of this paper is not to offer a critical review of past studies, but instead intends to inform the reader of what has been achieved with multibody modelling.

Development of a mandibular motion simulator for total joint replacement

PURPOSE

The purpose of this study was to develop a motion simulator capable of recreating and recording the full range of mandibular motions in a cadaveric preparation for an intact temporomandibular joint (TMJ) and after total joint replacement.

MATERIAL AND METHODS

A human cadaver head was used. Two sets of tracking balls were attached to the forehead and mandible, respectively. Computed tomographic (CT) scan was performed and 3-dimensional CT models of the skull were generated. The cadaver head was then dissected to attach the muscle activation cables and mounted onto the TMJ simulator. Realistic jaw motions were generated through the application of the following muscle forces: lateral pterygoid muscle, suprahyoid depressors (geniohyoid, mylohyoid, and digastric muscles), and elevator muscles. To simulate muscle contraction, cables were inserted into the mandible at the center area of each muscle's attachment. To provide a minimum mouth closing force at the initial position, the elevator muscles were combined at the anterior mandible. During mandibular movement, each motion was recorded using a high-resolution laser scanner. The right TMJ of the same head was reconstructed with a total TMJ prosthesis. The same forces were applied and the jaw motions were recorded again. CT scan was performed and 3-dimensional CT models of the skull with TMJ prosthesis were generated.

RESULTS

Mandibular motions, before and after TMJ replacement, with and without lateral pterygoid muscle reattachment, were re-created in a cadaveric preparation. The laser-scanned data during the mandibular motion were used to drive 3-dimensional CT models. A movie for each mandibular motion was subsequently created for motion path analysis. Compared with mandibular motion before TMJ replacement, mandibular lateral and protrusive motions after TMJ replacement, with and without lateral pterygoid muscle reattachment, were greatly limited. The jaw motion recorded before total joint replacement was applied to the mandibular and prostheses models after total TMJ replacement. The condylar component was observed sinking into the fossa during jaw motion.

CONCLUSION

A motion simulator capable of re-creating and recording full range of mandibular motions in a cadaveric preparation has been developed. It can be used to simulate mandibular motions for the intact TMJ and total joint prosthesis, and to re-create and record their full range of mandibular motions. In addition, the full range of the recorded motion can be re-created as motion images in a computer. These images can be used for motion path analysis and to study the causation of limited range of motion after total joint replacement and strategies for improvement.

Influence of abutment material on the fracture strength and failure modes of abutment-fixture assemblies when loaded in a bio-faithful simulation

OBJECTIVES

The aim of the present study was to evaluate differences in the ultimate fracture resistance of titanium and zirconia abutments.

MATERIAL AND METHODS

Twenty titanium fixtures were embedded in 20 resin mandible section simulators to mimic osseointegrated implants in the premolar area. The embedded implants were then randomly divided into two groups. Afterwards, specimens in group A (n=10) were connected to titanium abutments (TiDesign™ 3.5/4.0, 5.5, 1.5 mm), while specimens in group B (n=10) were connected to zirconia abutments (ZirDesign ™ 3.5/4.0, 5.5, 1.5 mm). Both groups were loaded to failure in a dynamometric testing machine. Fractured samples were then analyzed by scanning electron microscopy (SEM).

RESULTS

Group A showed a significantly higher fracture strength than that observed in group B. Group A failures were observed at the screw that connects the abutment with the implant while the abutment connection hexagons were plastically bent by the applied load. Group B failures were a result of abutment fractures. SEM analysis showed that in group A the screw failure was driven by crack nucleation, coalescence and propagation, while in group B, the SEM analysis of failed surfaces showed the conchoidal fracture profile characteristic of brittle materials.

CONCLUSIONS

The strength of both tested systems is adequate to resist physiologic chewing forces in the premolar area. Conversely, the titanium and zirconia failure modes evaluated here occurred at unphysiological loads. In addition, because the abutments were tested without crowns, the presented data have limited direct transfer to the clinical situation.

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.

[Validating a loading protocol for weakened mandibular models]

INTRODUCTION

The purpose of this study was to validate an experimental model of weakened synthetic mandible loaded in a device simulating masticatory forces by studying the reproducibility of the fracture threshold.

MATERIALS AND METHODS

Both sides of five synthetic mandibles were weakened by progressive synthetic bone resection. The synthetic mandibles were submitted, after each resection, to static chewing exercises in the simulator. The threshold fracture was recorded.

RESULTS

The fracture threshold mean resection height was 19.33 mm with a percentage of error of 6.72%.

DISCUSSION

The low percentage of error allowed us to validate this experimental model. A comparison with other experimental or finite element models confirmed our protocol and results. We plan to apply this protocol to fresh weakened mandibles.