Development of a finite element musculoskeletal model with the ability to predict contractions of three-dimensional muscles

Representation of realistic muscle geometries is needed for systematic biomechanical simulation of musculoskeletal systems. Most of the previous musculoskeletal models are based on multibody dynamics simulation with muscles simplified as one-dimensional (1D) line-segments without accounting for the large muscle attachment areas, spatial fibre alignment within muscles and contact and wrapping between muscles and surrounding tissues. In previous musculoskeletal models with three-dimensional (3D) muscles, contractions of muscles were among the inputs rather than calculated, which hampers the predictive capability of these models. To address these issues, a finite element musculoskeletal model with the ability to predict contractions of 3D muscles was developed. Muscles with realistic 3D geometry, spatial muscle fibre alignment and muscle-muscle and muscle-bone interactions were accounted for. Active contractile stresses of the 3D muscles were determined through an efficient optimization approach based on the measured kinematics of the lower extremity and ground force during gait. This model also provided stresses and strains of muscles and contact mechanics of the muscle-muscle and muscle-bone interactions. The total contact force of the knee predicted by the model corresponded well to the in vivo measurement. Contact and wrapping between muscles and surrounding tissues were evident, demonstrating the need to consider 3D contact models of muscles. This modelling framework serves as the methodological basis for developing musculoskeletal modelling systems in finite element method incorporating 3D deformable contact models of muscles, joints, ligaments and bones.     

Comparison of an EMG-based and a stress-based method to predict shoulder muscle forces

The estimation of muscle forces in musculoskeletal shoulder models is still controversial. Two different methods are widely used to solve the indeterminacy of the system: electromyography (EMG)-based methods and stress-based methods. The goal of this work was to evaluate the influence of these two methods on the prediction of muscle forces, glenohumeral load and joint stability after total shoulder arthroplasty. An EMG-based and a stress-based method were implemented into the same musculoskeletal shoulder model. The model replicated the glenohumeral joint after total shoulder arthroplasty. It contained the scapula, the humerus, the joint prosthesis, the rotator cuff muscles supraspinatus, subscapularis and infraspinatus and the middle, anterior and posterior deltoid muscles. A movement of abduction was simulated in the plane of the scapula. The EMG-based method replicated muscular activity of experimentally measured EMG. The stress-based method minimised a cost function based on muscle stresses. We compared muscle forces, joint reaction force, articular contact pressure and translation of the humeral head. The stress-based method predicted a lower force of the rotator cuff muscles. This was partly counter-balanced by a higher force of the middle part of the deltoid muscle. As a consequence, the stress-based method predicted a lower joint load (16% reduced) and a higher superior–inferior translation of the humeral head (increased by 1.2 mm). The EMG-based method has the advantage of replicating the observed cocontraction of stabilising muscles of the rotator cuff. This method is, however, limited to available EMG measurements. The stress-based method has thus an advantage of flexibility, but may overestimate glenohumeral subluxation.     

Maxillary expansion treatment using bone anchors: development and validation of a 3D finite element model

Objective: Develop a finite element (FE) model of a skull to perform biomechanical studies of maxillary expansion using bone anchors (BA).

Materials and methods: A skull model was developed and assigned material properties based on Hounsfield unit (HU) values of cone-beam computerized tomography (CBCT) images. A 3 mm diameter cylindrical BA was modelled and inserted in the palatal bone. A 4 mm transverse displacement was applied on the anchor. An evaluation on the effect on local stresses of BA implantation inclination angle was performed.

Results: Proper displacement results and strain–stress trends for the expansion process were present. Stress distribution patterns were similar as reported in the literature. No significant difference between BA inclination angles was found.

Conclusion: This work leads to a better understanding and prediction of craniofacial and maxillary bone remodelling during ME with BA treatments and is a first step towards the development of patient specific treatments.     

Influence of bicortical techniques in internal connection placed in premaxillary area by 3D finite element analysis

The aim of study was to evaluate the stress distribution in implant-supported prostheses and peri-implant bone using internal hexagon (IH) implants in the premaxillary area, varying surgical techniques (conventional, bicortical and bicortical in association with nasal floor elevation), and loading directions (0°, 30° and 60°) by three-dimensional (3D) finite element analysis. Three models were designed with Invesalius, Rhinoceros 3D and Solidworks software. Each model contained a bone block of the premaxillary area including an implant (IH, Ø4 × 10 mm) supporting a metal-ceramic crown. 178 N was applied in different inclinations (0°, 30°, 60°). The results were analyzed by von Mises, maximum principal stress, microstrain and displacement maps including ANOVA statistical test for some situations. Von Mises maps of implant, screws and abutment showed increase of stress concentration as increased loading inclination. Bicortical techniques showed reduction in implant apical area and in the head of fixation screws. Bicortical techniques showed slight increase stress in cortical bone in the maximum principal stress and microstrain maps under 60° loading. No differences in bone tissue regarding surgical techniques were observed. As conclusion, non-axial loads increased stress concentration in all maps. Bicortical techniques showed lower stress for implant and screw; however, there was slightly higher stress on cortical bone only under loads of higher inclinations (60°).     

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.     

Using 3D finite element models verified the importance of callus material and microstructure in biomechanics restoration during bone defect repair

Background: There is lack of further observations on the microstructure and material property of callus during bone defect healing and the relationships between callus properties and the mechanical strength. Methods: Femur bone defect model was created in rabbits and harvested CT data to reconstruct finite element models at 1 and 2 months. Three types of assumed finite element models were compared to study the callus properties, which assumed the material elastic property as heterogeneous (R-model), homogenous (H-model) or did not change from 1 to 2 months (U-model). Results: The apparent elastic moduli increased at 2 months (from 355.58 ± 132.67 to 1139.30 ± 967.43 MPa) in R-models. But there was no significant difference in apparent elastic moduli between R-models (355.58 ± 132.67 and 1139.30 ± 967.43 MPa) and H-models (344.79 ± 138.73 and 1001.52 ± 692.12 MPa) in 1 and 2 months. A significant difference of apparent elastic moduli was found between the R-model (1139.30 ± 967.43 MPa) and U-model group (207.15 ± 64.60 MPa) in 2 months. Conclusions: This study showed that the callus structure stability remodeled overtime to achieve a more effective structure, while the material quality of callus tissue is a very important factor for callus strength. At the meantime, this study showed an evidence that the material heterogeneity maybe not as important as it is in bone fracture model.     

Brief communication: comparing loading scenarios in lower first molar supporting bone structure using 3D finite element analysis

Finite element analysis (FEA) is a widespread technique to evaluate the stress/strain distributions in teeth or dental supporting tissues. However, in most studies occlusal forces are usually simplified using a single vector (i.e., point load) either parallel to the long tooth axis or oblique to this axis. In this pilot study we show how lower first molar occlusal information can be used to investigate the stress distribution with 3D FEA in the supporting bone structure. The LM₁ and the LP²‐LM¹ of a dried modern human skull were scanned by μCT in maximum intercuspation contact. A kinematic analysis of the surface contacts between LM₁ and LP²‐LM¹ during the power stroke was carried out in the occlusal fingerprint analyzer (OFA) software to visualize contact areas during maximum intercuspation contact. This information was used for setting the occlusal molar loading to evaluate the stress distribution in the supporting bone structure using FEA. The output was compared to that obtained when a point force parallel to the long axis of the tooth was loaded in the occlusal basin. For the point load case, our results indicate that the buccal and lingual cortical plates do not experience notable stresses. However, when the occlusal contact areas are considered, the disto‐lingual superior third of the mandible experiences high tensile stresses, while the medio‐lingual cortical bone is subjected to high compressive stresses. Developing a more realistic loading scenario leads to better models to understand the relationship between masticatory function and mandibular shape and structures. Am J Phys Anthropol, 2012. © 2011 Wiley Periodicals, Inc.     

Three-dimensional finite element analysis of the maxillary central incisor in two different situations of traumatic impact

Dental trauma is one of the most common events in dental practice. However, few studies have investigated the biomechanical characteristics of these injuries. The objective of this study was to analyse the stress distribution in the dentoalveolar structures of a maxillary central incisor subjected to two situations of impact loading. The following loading forces were applied using a 3D finite element model: a force of 2000 N acting at an angle of 90°on the buccal surface of the crown and a vertical 2000 N force acting in the cleidocranial direction on the incisal surface of the tooth. Harmful stresses were observed in both situations, causing damage to both the tooth and adjacent tissue. However, the damage found in soft tissues such as periodontal ligament and dental pulp was negligible. In conclusion, injuries resulting from the traumatic situations were more damaging to the integrity of the tooth and its associated hard-tissue structures.     

Mechanical significance of femoral head trabecular bone structure in Loris and Galago evaluated using micromechanical finite element models

Work on the interspecific and intraspecific variation of trabecular bone in the proximal femur of primates demonstrates important architectural variation between animals with different locomotor behaviors. This variation is thought to be related to the processes of bone adaptation whereby bone structure is optimized to the mechanical environment. Micromechanical finite element models were created for the proximal femur of the leaping Galago senegalensis and the climbing and quadrupedal Loris tardigradus by converting bone voxels from high-resolution X-ray computed tomography scans of the femoral head to eight-noded brick elements. The resulting models had approximately 1.8 million elements each. Loading conditions representing takeoff phase of a leap and more generalized load orientations were applied to the models, and the models were solved using the iterative "row-by-row" matrix-vector multiplication algorithm. The principal strain and Von Mises stress results for the leaping model were similar for both species at each load orientation. Similar hip joint reaction forces in the range of 4.9 x to 12 x body weight were calculated for both species under each loading condition, but the hip reaction values estimated for Loris were higher than predicted based on locomotor behavior. These results suggest that functional adaptation to hip joint loading may not fully explain the differences in femoral head trabecular bone structure in Galago and Loris. The finite element method represents a unique and useful tool for analyzing the functional adaptation of trabecular bone in a diversity of animals and for reconstructing locomotor behavior in extinct taxa.     

Validation performance comparison for finite element models of the human brain

The objective of this study was to compare the performance of six validated brain finite element (FE) models to localized brain motion validation data in five experimental configurations. Model performance was measured using the objective metric CORA (CORrelation and Analysis), where higher ratings represent better correlation. The KTH model achieved the highest average CORA rating, and the ABM received the highest average rating among models robustly validated against five cadaver impacts in three directions. This technique can be more frequently employed to build better models and, when associated limitations are well understood, to compare inter-model performance under similar conditions.     

Computer simulation of orthodontic tooth movement using CT image-based voxel finite element models with the level set method

Orthodontic tooth movement (OTM) is an adaptive biomechanical response of dentoalveolar components to orthodontic forces, in which remodeling of the alveolar bone occurs in response to changes in the surrounding mechanical environment. In this study, we developed a framework for OTM simulation by combining an image-based voxel finite element method, with a surface-tracking level set method using three-dimensional computer models. For a case study to demonstrate its capability of expressing clinical tooth movement, we observed displacement and rotation of the tooth under three types of force conditions. The simulation results demonstrate that the proposed simulation method has the potential to predict clinical OTM.     

Validation of a full body finite element model (THUMS) for running-type impacts to the lower extremity

The purpose of this study was to determine whether modifying an existing, highly biofidelic full body finite element model [total human model for safety (THUMS)] would produce valid amplitude and temporal shock wave characteristics as it travels proximally through the lower extremity. Modifying an existing model may be more feasible than developing a new model, in terms of cost, labour and expertise. The THUMS shoe was modified to more closely simulate the material properties of a heel pad. Relative errors in force and acceleration data from experimental human pendulum impacts and simulated THUMS impacts were 22% and 54%, respectively, across the time history studied. The THUMS peak acceleration was attenuated by 57.5% and took 19.7 ms to travel proximally along the lower extremity. Although refinements may be necessary to improve force and acceleration timing, the modified THUMS represented, to a certain extent, shock wave propagation and attenuation demonstrated by living humans under controlled impact conditions.     

3D finite element modeling of pelvic organ prolapse

Objectives: The purpose of this study is to develop a validated 3D finite element model of the pelvic floor system which can offer insights into the mechanics of anterior vaginal wall prolapse and have the ability to assess biomedical device treatment methods. The finite element results should accurately mimic the clinical findings of prolapse due to intra-abdominal pressure (IAP) and soft tissues impairment conditions. Methods: A 3D model of pelvic system was created in Creo Parametric 2.0 based on MRI Images, which included uterus, cervix, vagina, cardinal ligaments, uterosacral ligaments, and a simplified levator plate and rectum. The geometrical model was imported into ANSYS Workbench 14.5. Mechanical properties of soft tissues were based on experimental data of tensile test results from current literature. Studies were conducted for IAP loadings on the vaginal wall and uterus, increasing from lowest to extreme values. Results: Anterior vaginal wall collapse occurred at an IAP value corresponding to maximal valsalva and showed similar collapsed shape as clinical findings. Prolapse conditions exhibited high sensitivity to vaginal wall stiffness, whereas healthy tissues was found to support the vagina against prolapse. Ligament impairment was found to have only a secondary effect on prolapse.     

Computational modelling of the natural hip: a review of finite element and multibody simulations

Primary objective: The hip joint suffers from a high prevalence of degenerative conditions. Athough patient's well-being could be improved through early and more effective interventions, without a greater understanding of the mechanics of the hip, these developments cannot be attained. Thus, this review article summarises the current literature on this subject in order to provide a platform for future developments. To illustrate the influence computational simulations have had on the knowledge advancement in hip mechanics, we explored two methodological approaches: finite element (FE) analysis and multibody dynamics (MBD). Main outcomes and results: Notwithstanding the unique capabilities of FE and MBD, the former generally offers the micromechanics of the articulating surfaces whereas the latter the macromechanics of the skeleton, these two methodologies also provide the bulk of the literature regarding computational modelling of the musculoskeletal system. Although FE has provided significant knowledge on contact pressures and the effects of musculoskeletal geometries, in particular cartilage and bone shapes, MBD has afforded a wealth of understanding on the influence of gait patterns and muscle attachment locations on force magnitudes. Conclusions: These two computational techniques have, and will continue to, provide significant contributions towards the development of interventions. It is hoped that this article will help focus ongoing technological developments by highlighting areas of success, but also areas of under development.     

Volumetric locking free 3D finite element for modelling of anisotropic visco-hyperelastic behaviour of anterior cruciate ligament

Solids such as polymers, soft biological tissues display visco-hyperelastic, isochoric and finite deformation behaviour. The incompressibility constraint imposed severe restriction on the displacement field results in volumetric locking. Many techniques have been developed to address the issue such as reduced integration, mixed formulation, B-Bar and F-Bar methods, each of them with their own merits and demerits. In this work, we have developed a 3D finite element (hereby referred as J-Bar method) to counter volumetric locking in visco-hyperelastic solids. To validate the proposed J-Bar method, rheological characteristics of the human anterior cruciate ligament (ACL) were predicted and compared with the experimental results.     

Simulation of 3D tumor cell growth using nonlinear finite element method

We propose a novel parallel computing framework for a nonlinear finite element method (FEM)-based cell model and apply it to simulate avascular tumor growth. We derive computation formulas to simplify the simulation and design the basic algorithms. With the increment of the proliferation generations of tumor cells, the FEM elements may become larger and more distorted. Then, we describe a remesh and refinement processing of the distorted or over large finite elements and the parallel implementation based on Message Passing Interface to improve the accuracy and efficiency of the simulation. We demonstrate the feasibility and effectiveness of the FEM model and the parallelization methods in simulations of early tumor growth.     

Mechanical performance of endodontic restorations with prefabricated posts: sensitivity analysis of parameters with a 3D finite element model

Many studies have investigated the effect of different parameters of the endodontically restored tooth on its final strength, using in vitro tests and model simulations. However, the differences in the experimental set-up or modelling conditions and the limited number of parameters studied in each case prevent us from obtaining clear conclusions about the relative importance of each parameter. In this study, a validated 3D biomechanical model of the restored tooth was used for an exhaustive sensitivity analysis. The individual influence of 20 different parameters on the mechanical performance of an endodontic restoration with prefabricated posts was studied. The results bring up the remarkable importance of the loading angle on the final restoration strength. Flexural loads are more critical than compressive or tensile loads. Young's modulus of the post and its length and diameter are the most influential parameters for strength, whereas other parameters such as ferrule geometry or core and crown characteristics are less significant.