Predicting muscle patterns for hemimandibulectomy models

The effect of simulated radiation induced fibrosis on tongue protrusion

Radiation therapy (RT) is an important adjuvant and primary treatment modality for head and neck cancers. A severe side effect of RT is fibrosis or scarring of muscle tissues of the oral cavity including the tongue. Previous studies have demonstrated that increased radiation doses to the oral cavity structures have led to decrements in function, hypothesized to result from changes in muscle tissue properties that affect the tongue's function. To understand the complex relationship between tongue muscle fibrosis and tongue function, the current study used a virtual biomechanical model of the tongue. Fibrosis parameters including density (high, low), area (large, small) and location (946 node centres) were systematically varied in the model to test its impact on a target tongue tip motion (protrusion). The impact of fibrosis lesion parameters on three directional components of the tip (anterior-inferior, lateral-medial, and superior-inferior) were analyzed using multi linear regression models. Increases in density and area of fibrosis significantly predicted tongue protrusion movements compared to baseline. In the anterior-posterior direction, reductions in the tongue protrusion were observed. In the inferior-superior direction, the tongue height remained above baseline for the majority of cases. In the lateral-medial direction, ipsilateral deviations were observed. The location of fibrosis modulated these three main effects by either amplifying the observed effect or minimizing it. The findings support the hypothesis that changes in muscle tissue properties because of fibrosis impact tongue function. Increases in density and area of fibrosis impact key muscles in the target motion. The range of modulating effects of the lesion location (i.e., either amplifying or minimizing certain impact patterns) highlights the intricacy of tongue anatomy/soft tissue biomechanics and may suggest that lesions in any location will compromise the tongue's movement.

To what extent can mastication functionality be restored following mandibular reconstruction surgery? A computer modeling approach

STATEMENT OF PROBLEM

Advanced cases of head and neck cancer involving the mandible often require surgical removal of diseased sections and subsequent replacement with donor bone. During the procedure, the surgeon must make decisions regarding which bones or tissues to resect. This requires balancing tradeoffs related to issues such as surgical access and post-operative function; however, the latter is often difficult to predict, especially given that long-term functionality also depends on the impact of post-operative rehabilitation programs.

PURPOSE

To assist in surgical decision-making, we present an approach for estimating the effects of reconstruction on key aspects of post-operative mandible function.

MATERIAL AND METHODS

We develop dynamic biomechanical models of the reconstructed mandible considering different defect types and validate them using literature data. We use these models to estimate the degree of functionality that might be achieved following post-operative rehabilitation.

RESULTS

We find significant potential for restoring mandibular functionality, even in cases involving large defects. This entails an average trajectory error below 2 mm, bite force comparable to a healthy individual, improved condyle mobility, and a muscle activation change capped at a maximum of 20%.

CONCLUSION

These results suggest significant potential for adaptability in the masticatory system and improved post-operative rehabilitation, leading to greater restoration of jaw function.

Management of atypical occlusal discrepancy after condylar resection: A clinical report

This clinical report describes the prosthetic management of occlusion for a patient who had received condylar resection as part of cancer treatment. Previous reports have identified that patients with unrepaired segmental resection of the mandible experienced a frontal plane rotation of the mandible toward the nonsurgical side. In contrast, because of preservation of temporomandibular muscles and their attachments, the mandible rotated toward the surgical side, and occlusal contacts were limited to a pair of molars on that side. Manual manipulation and instructions for muscular stretching and massages were provided to reduce muscular tension. A mandibular guidance prosthesis was fabricated and gradually adjusted to guide the mandible progressively toward a normal position. These treatments helped improve general comfort, mastication, occlusion, and the gradual rotation of the mandible toward a normal position.

Computational models and their applications in biomechanical analysis of mandibular reconstruction surgery

Advanced head and neck cancers involving the mandible often require surgical removal of the diseased parts and replacement with donor bone or prosthesis to recreate the form and function of the premorbid mandible. The degree to which this reconstruction successfully replicates key geometric features of the original bone critically affects the cosmetic and functional outcomes of speaking, chewing, and breathing. With advancements in computational power, biomechanical modeling has emerged as a prevalent tool for predicting the functional outcomes of the masticatory system and evaluating the effectiveness of reconstruction procedures in patients undergoing mandibular reconstruction surgery. These models offer cost-effective and patient-specific treatment tailored to the needs of individuals. To underscore the significance of biomechanical modeling, we conducted a review of 66 studies that utilized computational models in the biomechanical analysis of mandibular reconstruction surgery. The majority of these studies employed finite element method (FEM) in their approach; therefore, a detailed investigation of FEM has also been provided. Additionally, we categorized these studies based on the main components analyzed, including bone flaps, plates/screws, and prostheses, as well as their design and material composition.

Muscle-driven forward dynamic active hybrid model of the lumbosacral spine: combined FEM and multibody simulation

Most spine models belong to either the musculoskeletal multibody (MB) or finite element (FE) method. Recently, coupling of MB and FE models has increasingly been used to combine advantages of both methods. Active hybrid FE-MB models, still rarely used in spine research, avoid the interface and convergence problems associated with model coupling. They provide the inherent ability to account for the full interplay of passive and active mechanisms for spinal stability. In this paper, we developed and validated a novel muscle-driven forward dynamic active hybrid FE-MB model of the lumbosacral spine (LSS) in ArtiSynth to simultaneously calculate muscle activation patterns, vertebral movements, and internal mechanical loads. The model consisted of the rigid vertebrae L1-S1 interconnected with hyperelastic fiber-reinforced FE intervertebral discs, ligaments, facet joints, and force actuators representing the muscles. Morphological muscle data were implemented via a semi-automated registration procedure. Four auxiliary bodies were utilized to describe non-linear muscle paths by wrapping and attaching the anterior abdominal muscles. This included an abdominal plate whose kinematics was optimized using motion capture data from upper body movements. Intra-abdominal pressure was calculated from the forces of the abdominal muscles compressing the abdominal cavity. For the muscle-driven approach, forward dynamics assisted data tracking was used to predict muscle activation patterns that generate spinal postures and balance the spine without prescribing accurate spinal kinematics. During calibration, the maximum specific muscle tension and spinal rhythms resulting from the model dynamics were evaluated. To validate the model, load cases were simulated from -10° extension to +30° flexion with weights up to 20 kg in both hands. The biomechanical model responses were compared with in vivo literature data of intradiscal pressures, intra-abdominal pressures, and muscle activities. The results demonstrated high agreement with this data and highlight the advantages of active hybrid modeling for the LSS. Overall, this new self-contained tool provides a robust and efficient estimation of LSS biomechanical responses under in vivo similar loads, for example, to improve pain treatment by spinal stabilization therapies.

Calibration and validation of a novel hybrid model of the lumbosacral spine in ArtiSynth-The passive structures

In computational biomechanics, two separate types of models have been used predominantly to enhance the understanding of the mechanisms of action of the lumbosacral spine (LSS): Finite element (FE) and musculoskeletal multibody (MB) models. To combine advantages of both models, hybrid FE-MB models are an increasingly used alternative. The aim of this paper is to develop, calibrate, and validate a novel passive hybrid FE-MB open-access simulation model of a ligamentous LSS using ArtiSynth. Based on anatomical data from the Male Visible Human Project, the LSS model is constructed from the L1-S1 rigid vertebrae interconnected with hyperelastic fiber-reinforced FE intervertebral discs, ligaments, and facet joints. A mesh convergence study, sensitivity analyses, and systematic calibration were conducted with the hybrid functional spinal unit (FSU) L4/5. The predicted mechanical responses of the FSU L4/5, the lumbar spine (L1-L5), and the LSS were validated against literature data from in vivo and in vitro measurements and in silico models. Spinal mechanical responses considered when loaded with pure moments and combined loading modes were total and intervertebral range of motions, instantaneous axes and centers of rotation, facet joint contact forces, intradiscal pressures, disc bulges, and stiffnesses. Undesirable correlations with the FE mesh were minimized, the number of crisscrossed collagen fiber rings was reduced to five, and the individual influences of specific anatomical structures were adjusted to in vitro range of motions. Including intervertebral motion couplings for axial rotation and nonlinear stiffening under increasing axial compression, the predicted kinematic and structural mechanics responses were consistent with the comparative data. The results demonstrate that the hybrid simulation model is robust and efficient in reproducing valid mechanical responses to provide a starting point for upcoming optimizations and extensions, such as with active skeletal muscles.

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.

A comprehensive and volumetric musculoskeletal model for the dynamic simulation of the shoulder function

We present a volumetric and extensive finite element model of the shoulder usable in the context of inverse control, in which the scapula is left unconstrained on the ribcage. Such a model allows for exploring various shoulder movements, which are essential for making patient-specific decisions. The proposed model consists of 23 volumetric muscles parts modelled using the finite element method. The glenohumeral, acromioclavicular and sternoclavicular joints are modelled with soft ball-socket constraints. The musculoskeletal model can be controlled by a tracking-based algorithm, finding the excitations values in the muscles needed to follow some target points. The moment arms obtained during abduction and rotation are compared with the literature, which includes results from cadaveric data and a fine FE model of the rotator cuff and the deltoid. We simulated the paralysis of serratus anterior, a main reason of scapular winging, and compared it with its physiological counterpart. A deficiency in the range of motion as well as a reduction in upward rotation were observed, which both corroborate clinical observations. This is one of the most comprehensive model of the shoulder, which can be used to study complex pathologies of the shoulder and their impact on functional outcome such as range-of-motion.

Fast Forward-Dynamics Tracking Simulation: Application to Upper Limb and Shoulder Modeling

OBJECTIVE

Musculoskeletal simulation can be used to estimate muscle forces in clinical movement studies. However, such simulations typically only target movement measurements and are not applicable to force exertion tasks which are commonly used in rehabilitation therapy. Simulations can also produce nonphysiological joint forces or be too slow for real-time clinical applications, such as rehabilitation with real-time feedback. The objective of this study is to propose and evaluate a new formulation of forward-dynamics assisted tracking simulation that incorporates measured reaction forces as targets or constraints without any additional computational cost.

METHODS

We illustrate our method with idealized proof-of-concept models and evaluate it with two upper limb cases: Tracking of hand reaction forces during an isometric force-generation task and constraining glenohumeral joint reaction forces for stability during arm elevation.

RESULTS

We show that the addition of reaction force optimization terms within our simulations generates plausible muscle force predictions for these tasks, which are strongly related to reaction forces in addition to movement. Execution times for all models tested were not different when run with or without the reaction force optimization term, ensuring that the simulations are fast enough for real-time clinical applications.

CONCLUSION

Our novel reaction force optimization term leads to more realistic shoulder reaction forces, without any additional computational costs.

SIGNIFICANCE

Our formulation is not only valuable for shoulder simulations, but could be used in various clinical situations (e.g., for different joints and rehabilitation therapy tasks) where the direction and/or magnitude of reaction forces are of interest.

Three-dimensional longitudinal changes in craniofacial growth in untreated hemifacial microsomia patients with cone-beam computed tomography

INTRODUCTION

The purpose of this study was to evaluate the concept that the affected and contralateral sides do not grow at the same rate in patients with hemifacial microsomia. Changes in the cranial base, maxilla, mandible, and occlusal plane were evaluated on 3-dimensional images from cone-beam computed tomography data in untreated patients.

METHODS

Six patients were classified as having mandibular Pruzansky/Kaban type I, IIA, or IIB hemifacial microsomia. Cone-beam computed tomography (MercuRay; Hitachi, Tokyo, Japan) scans were taken before orthodontic treatment during both growth and postpuberty periods.

RESULTS

The cranial base as defined by the position of the mastoid process was in a different position between the affected and contralateral control sides. The nasomaxillary length or height was shorter on the affected side for all 6 patients with hemifacial microsomia regardless of its severity, and it grew less than on the contralateral control side in 5 of the 6 patients. The occlusal plane angle became more inclined in 4 of the 6 patients. The mandibular ramus was shorter on the affected side in all patients and grew less on the affected side in 5 of the 6 patients. The mandibular body grew slower, the same, or faster than on the control side.

CONCLUSIONS

The cranial base, position of the condyle, lengths of the condyle and ramus, and positions of the gonial angle and condyle can vary between the affected and contralateral control sides of patients with hemifacial microsomia, with the ramus and nasomaxillary length usually growing slower than they grow on the control side. These results suggest that many factors affect the growth rate of the craniofacial region and, specifically, the mandible in patients with hemifacial microsomia.

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.

Physics-based deformable tongue visualization

In this paper, a physics-based framework is presented to visualize the human tongue deformation. The tongue is modeled with the Finite Element Method (FEM) and driven by the motion capture data gathered during speech production. Several novel deformation visualization techniques are presented for in-depth data analysis and exploration. To reveal the hidden semantic information of the tongue deformation, we present a novel physics-based volume segmentation algorithm. This is accomplished by decomposing the tongue model into segments based on its deformation pattern with the computation of deformation subspaces and fitting the target deformation locally at each segment. In addition, the strain energy is utilized to provide an intuitive low-dimensional visualization for the high-dimensional sequential motion. Energy-interpolation-based morphing is also equipped to effectively highlight the subtle differences of the 3D deformed shapes without any visual occlusion. Our experimental results and analysis demonstrate the effectiveness of this framework. The proposed methods, though originally designed for the exploration of the tongue deformation, are also valid for general deformation analysis of other shapes.

Automatic prediction of tongue muscle activations using a finite element model

Computational modeling has improved our understanding of how muscle forces are coordinated to generate movement in musculoskeletal systems. Muscular-hydrostat systems, such as the human tongue, involve very different biomechanics than musculoskeletal systems, and modeling efforts to date have been limited by the high computational complexity of representing continuum-mechanics. In this study, we developed a computationally efficient tracking-based algorithm for prediction of muscle activations during dynamic 3D finite element simulations. The formulation uses a local quadratic-programming problem at each simulation time-step to find a set of muscle activations that generated target deformations and movements in finite element muscular-hydrostat models. We applied the technique to a 3D finite element tongue model for protrusive and bending movements. Predicted muscle activations were consistent with experimental recordings of tongue strain and electromyography. Upward tongue bending was achieved by recruitment of the superior longitudinal sheath muscle, which is consistent with muscular-hydrostat theory. Lateral tongue bending, however, required recruitment of contralateral transverse and vertical muscles in addition to the ipsilateral margins of the superior longitudinal muscle, which is a new proposition for tongue muscle coordination. Our simulation framework provides a new computational tool for systematic analysis of muscle forces in continuum-mechanics models that is complementary to experimental data and shows promise for eliciting a deeper understanding of human tongue function.

Positional changes of the masseter and medial pterygoid muscles after surgical mandibular advancement procedures: an MRI study

This study evaluated whether surgical mandibular advancement procedures induced a change in the direction and the moment arms of the masseter (MAS) and medial pterygoid (MPM) muscles. Sixteen adults participated in this study. The sample was divided in two groups: Group I (n=8) with a mandibular plane angle (mpa) <39° and Group II (n=8) with an mpa >39°. Group I patients were treated with a bilateral sagittal split osteotomy (BSSO). Those in Group II were treated with a BSSO combined with a Le Fort I osteotomy. Pre- and postoperative direction and moment arms of MAS and MPM were compared in these groups. Postsurgically, MAS and MPM in Group II showed a significantly more vertical direction in the sagittal plane. Changes of direction in the frontal plane and changes of moment arms were insignificant in both groups. This study demonstrated that bimaxillary surgery in patients with an mpa >39° leads to a significant change of direction of MAS and MPM in the sagittal plane.

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.