Three-dimensional finite element stress analysis of the dentate human mandible.

The biomechanical events which accompany functional loading of the human mandible are not fully understood. The techniques normally used to record them are highly invasive. Computer modelling offers a promising alternative approach in this regard, with the additional ability to predict regional stresses and strains in inaccessible locations. In this study, we built two three-dimensional finite element (FE) models of a human mandible reconstructed from tomographs of a dry dentate jaw. The first model was used for a complete mechanical characterization of physical events. It also provided comparative data for the second model, which had an increased vertical corpus depth. In both cases, boundary conditions included rigid restraints at the first right molar and endosteal cortical surfaces of the articular eminences of temporal bones. Groups of parallel multiple vectors simulated individual masticatory muscle loads. The models were solved for displacements, stresses, strains, and forces. The simulated muscle loads in the first model deformed the mandible helically upward and toward its right (working) side. The highest principal stresses occurred at the bite point, anterior aspects of the coronoid processes, symphyseal region, and right and left sides of the mandibular corpus. In general, the observed principal stresses and strains were highest on the periosteal cortical surface and alveolar bone. At the symphyseal region, maximum principal stresses and strains were highest on the lower lingual mandibular aspect, whereas minimum principal stresses and strains were highest on its upper labial side. Subcondylar principal strains and condylar forces were higher on the left (balancing or nonbiting) side than on the right mandibular side, with condylar forces more concentrated on the anteromedial aspect of the working-side condyle and on the central and lateral aspects of the left. When compared with in vivo strain data from macaques during comparable biting events, the predictive strain values from the first model were qualitatively similar. In the second model, the reduced tensile stress on the working-side, and decreased shear stress bilaterally, confirmed that lower stresses occurred on the lower mandibular border with increased jaw depth. Our results suggested that although the mandible behaved in a beam-like manner, its corpus acted more like a combination of open and closed cross sections due to the presence of tooth sockets, at least for the task modelled.(ABSTRACT TRUNCATED AT 400 WORDS)     

The effect of boundary condition on the biomechanics of a human pelvic joint under an axial compressive load: a three-dimensional finite element model

The finite element (FE) model of the pelvic joint is helpful for clinical diagnosis and treatment of pelvic injuries. However, the effect of an FE model boundary condition on the biomechanical behavior of a pelvic joint has not been well studied. The objective of this study was to study the effect of boundary condition on the pelvic biomechanics predictions. A 3D FE model of a pelvis using subject-specific estimates of intact bone structures, main ligaments and bone material anisotropy by computed tomography (CT) gray value was developed and validated by bone surface strains obtained from rosette strain gauges in an in vitro pelvic experiment. Then three FE pelvic models were constructed to analyze the effect of boundary condition, corresponding to an intact pelvic joint, a pelvic joint without sacroiliac ligaments and a pelvic joint without proximal femurs, respectively. Vertical load was applied to the same pelvis with a fixed prosthetic femoral stem and the same load was simulated in the FE model. A strong correlation coefficient (R(2)=0.9657) was calculated, which indicated a strong correlation between the FE analysis and experimental results. The effect of boundary condition changes on the biomechanical response depended on the anatomical location and structure of the pelvic joint. It was found that acetabulum fixed in all directions with the femur removed can increase the stress distribution on the acetabular inner plate (approximately double the original values) and decrease that on the superior of pubis (from 7 MPa to 0.6 MPa). Taking sacrum and ilium as a whole, instead of sacroiliac and iliolumber ligaments, can influence the stress distribution on ilium and pubis bone vastly. These findings suggest pelvic biomechanics is very dependent on the boundary condition in the FE model.     

Bone ingrowth around porous-coated acetabular implant: a three-dimensional finite element study using mechanoregulatory algorithm

Fixation of uncemented implant is influenced by peri-prosthetic bone ingrowth, which is dependent on the mechanical environment of the implant–bone structure. The objective of the study is to gain an insight into the tissue differentiation around an acetabular component. A mapping framework has been developed to simulate appropriate mechanical environment in the three-dimensional microscale model, implement the mechanoregulatory tissue differentiation algorithm and subsequently assess spatial distribution of bone ingrowth around an acetabular component, quantitatively. The FE model of implanted pelvis subjected to eight static load cases during a normal walking cycle was first solved. Thereafter, a mapping algorithm has been employed to include the variations in implant–bone relative displacement and host bone material properties from the macroscale FE model of implanted pelvis to the microscale FE model of the beaded implant–bone interface. The evolutionary tissue differentiation was observed in each of the 13 microscale models corresponding to 13 acetabular regions. The total implant–bone relative displacements, averaged over each region of the acetabulum, were found to vary between 10 and 60 \(\upmu \hbox {m}\). Both the linear elastic and biphasic poroelastic models predicted similar mechanoregulatory peri-prosthetic tissue differentiation. Considerable variations in bone ingrowth (13–88 %), interdigitation depth (0.2–0.82 mm) and average tissue Young’s modulus (970–3430 MPa) were predicted around the acetabular cup. A progressive increase in the average Young’s modulus, interdigitation depth and decrease in average radial strains of newly formed tissue layer were also observed. This scheme can be extended to investigate tissue differentiation for different surface texture designs on the implants.     

Considerations for reporting finite element analysis studies in biomechanics

Simulation-based medicine and the development of complex computer models of biological structures is becoming ubiquitous for advancing biomedical engineering and clinical research. Finite element analysis (FEA) has been widely used in the last few decades to understand and predict biomechanical phenomena. Modeling and simulation approaches in biomechanics are highly interdisciplinary, involving novice and skilled developers in all areas of biomedical engineering and biology. While recent advances in model development and simulation platforms offer a wide range of tools to investigators, the decision making process during modeling and simulation has become more opaque. Hence, reliability of such models used for medical decision making and for driving multiscale analysis comes into question. Establishing guidelines for model development and dissemination is a daunting task, particularly with the complex and convoluted models used in FEA. Nonetheless, if better reporting can be established, researchers will have a better understanding of a model's value and the potential for reusability through sharing will be bolstered. Thus, the goal of this document is to identify resources and considerate reporting parameters for FEA studies in biomechanics. These entail various levels of reporting parameters for model identification, model structure, simulation structure, verification, validation, and availability. While we recognize that it may not be possible to provide and detail all of the reporting considerations presented, it is possible to establish a level of confidence with selective use of these parameters. More detailed reporting, however, can establish an explicit outline of the decision-making process in simulation-based analysis for enhanced reproducibility, reusability, and sharing.     

Using the finite element method to model the biomechanics of the asymmetric mandible before, during and after skeletal correction by distraction osteogenesis

An approach was developed to evaluate the load transfer mechanism in the temporomandibular joint (TMJ) area before, during and after mandibular ramus elongation by distraction osteogenesis (DO). In a concerted approach using computer tomography, magnetic resonance imaging (MRI), and finite element analysis, three-dimensional numerical models based on a young male patient, with a dento-facial deformity were generated. The magnitude and direction of the muscle forces acting on the mandible were assessed using both values derived from the muscles volume and cross-section as retrieved from the MRI-scan data-sets and taken from the literature. The resistance of the soft tissue envelope towards elongation during the DO-phase was also included. The finite element analyses showed that before skeletal correction by DO the load transfer was asymmetrical with high peak stresses in the affected joint. Following ramus elongation a more symmetrical loading in TMJs was predicted. The reaction forces in the TMJs during DO were low.     

A fluid-structure interaction problem in biomechanics: prestressed vibrations of the eye by the finite element method.

The object of this work has been to develop a mechanical and numerical model of the eye submitted to vibrations, and in particular, to calculate the influence of intraocular pressure (IOP) on the eye resonance frequencies. Our mechanical model of the eye relies upon the theory of the mechanics of continuous media. The numerical model results from a model analysis of the vibrations of the eye using a finite element method (FEM) for discretization. The eye can be schematically represented as a prestressed shell, filled by an inviscid barotropic compressible fluid, which leads us to formulate and solve a problem of vibrations of a coupled fluid-structure system. The corneoscleral shell has been modeled as a thin and thick shell, taking into account material nonlinearities in the thick case. Numerical results obtained for the attached eye demonstrate a fair sensitivity of the resonance frequencies to the variations of the IOP; thus, founding the interest of the surveillance of the resonance frequency of the eye.     

A three-dimensional finite element analysis of the upper tibia

A three-dimensional finite element model of the proximal tibia has been developed to provide a base line for further modeling of prosthetic resurfaced tibiae. The geometry for the model was developed by digitizing coronal and transverse sections made with the milling machine, from one fresh tibia of average size. The load is equally distributed between the medial and lateral compartments over contact areas that were reported in the literature. An indentation test has been used to measure the stiffness and the ultimate strength of cancellous bone in four cadaver tibiae. These values provided the statistical basis for characterising the inhomogeneous distribution of the cancellous bone properties in the proximal tibia. All materials in the model were assumed to be linearly elastic and isotropic. Mechanical properties for the cortical bone and cartilage have been taken from the literature. Results have been compared with strain gage tests and with a two-dimensional axisymmetric finite element model both from the literature. Qualitative comparison between trabecular alignment, and the direction of the principal compressive stresses in the cancellous bone, showed a good relationship. Maximum stresses in the cancellous bone and cortical bone, under a load which occurs near stance phase during normal gait, show safety factors of approximately eight and twelve, respectively. The load sharing between the cancellous bone and the cortical bone has been plotted for the first 40 mm distally from the tibial eminence.     

Modeling leaflet correction techniques in aortic valve repair: A finite element study

In aortic valve sparing surgery, cusp prolapse is a common cause of residual aortic insufficiency. To correct cusp pathology, native leaflets of the valve frequently require adjustment which can be performed using a variety of described correction techniques, such as central or commissural plication, or resuspension of the leaflet free margin. The practical question then arises of determining which surgical technique provides the best valve performance with the most physiologic coaptation. To answer this question, we created a new finite element model with the ability to simulate physiologic function in normal valves, and aortic insufficiency due to leaflet prolapse in asymmetric, diseased or sub-optimally repaired valves. The existing leaflet correction techniques were simulated in a controlled situation, and the performance of the repaired valve was quantified in terms of maximum leaflets stress, valve orifice area, valve opening and closing characteristics as well as total coaptation area in diastole. On the one hand, the existing leaflet correction techniques were shown not to adversely affect the dynamic properties of the repaired valves. On the other hand, leaflet resuspension appeared as the best technique compared to central or commissural leaflet plication. It was the only method able to achieve symmetric competence and fix an individual leaflet prolapse while simultaneously restoring normal values for mechanical stress, valve orifice area and coaptation area.     

RETRACTED: Modelling subcortical bone in finite element analyses: A validation and sensitivity study in the macaque mandible

This article has been retracted: please see Elsevier Policy on Article Withdrawal (http://www.elsevier.com/locate/withdrawalpolicy).This article has been retracted at the request of the authors who unanimously wish to retract this paper because of errors in the validation protocol and data, and the finite element analysis data and results.     

Mechanics of the right whale mandible: Full scale testing and finite element analysis

In an effort to better understand the mechanics of ship-whale collision and to reduce the associated mortality of the critically endangered North Atlantic right whale, a comprehensive biomechanical study has been conducted by the Woods Hole Oceanographic Institution and the University of New Hampshire. The goal of the study is to develop a numerical modeling tool to predict the forces and stresses during impact and thereby the resulting mortality risk to whales from ship strikes.Based on post-mortem examinations, jaw fracture was chosen as a fatal endpoint for the whales hit by a vessel. In this paper we investigate the overall mechanical behavior of a right whale mandible under transverse loading and develop a finite element analysis model of the bone. The equivalent elastic modulus of the cortical component of right whale mandible is found by comparing full-scale bending tests with the results of numerical modeling. The finite element model of the mandible can be used in conjunction with a vessel-whale collision event model to predict bone fracture for various ship strike scenarios.     

Mixed-mode loading of the cement-bone interface: a finite element study

While including the cement-bone interface of complete cemented hip reconstructions is crucial to correctly capture their response, its modelling is often overly simplified. In this study, the mechanical mixed-mode response of the cement-bone interface is investigated, taking into account the effects of the well-defined microstructure that characterises the interface. Computed tomography-based plain strain finite element analyses models of the cement-bone interface are built and loaded in multiple directions. Periodic boundaries are considered and the failure of the cement and bone fractions by cracking of the bulk components are included. The results compare favourably with experimental observations. Surprisingly, the analyses reveal that under shear loading no failure occurs and considerable normal compression is generated to prevent interface dilation. Reaction forces, crack patterns and stress fields provide more insight into the mixed-mode failure process. Moreover, the cement-bone interface analyses provide details which can serve as a basis for the development of a cohesive law.     

A three-dimensional inverse finite element analysis of the heel pad

Quantification of plantar tissue behavior of the heel pad is essential in developing computational models for predictive analysis of preventive treatment options such as footwear for patients with diabetes. Simulation based studies in the past have generally adopted heel pad properties from the literature, in return using heel-specific geometry with material properties of a different heel. In exceptional cases, patient-specific material characterization was performed with simplified two-dimensional models, without further evaluation of a heel-specific response under different loading conditions. The aim of this study was to conduct an inverse finite element analysis of the heel in order to calculate heel-specific material properties in situ. Multidimensional experimental data available from a previous cadaver study by Erdemir et al. ("An Elaborate Data Set Characterizing the Mechanical Response of the Foot," ASME J. Biomech. Eng., 131(9), pp. 094502) was used for model development, optimization, and evaluation of material properties. A specimen-specific three-dimensional finite element representation was developed. Heel pad material properties were determined using inverse finite element analysis by fitting the model behavior to the experimental data. Compression dominant loading, applied using a spherical indenter, was used for optimization of the material properties. The optimized material properties were evaluated through simulations representative of a combined loading scenario (compression and anterior-posterior shear) with a spherical indenter and also of a compression dominant loading applied using an elevated platform. Optimized heel pad material coefficients were 0.001084 MPa (μ), 9.780 (α) (with an effective Poisson's ratio (ν) of 0.475), for a first-order nearly incompressible Ogden material model. The model predicted structural response of the heel pad was in good agreement for both the optimization (<1.05% maximum tool force, 0.9% maximum tool displacement) and validation cases (6.5% maximum tool force, 15% maximum tool displacement). The inverse analysis successfully predicted the material properties for the given specimen-specific heel pad using the experimental data for the specimen. The modeling framework and results can be used for accurate predictions of the three-dimensional interaction of the heel pad with its surroundings.     

How plate positioning impacts the biomechanics of the open wedge tibial osteotomy; a finite element analysis

A numerical model of the medial open wedge tibial osteotomy based on the finite element method was developed. Two plate positions were tested numerically. In a configuration, (a), the plate was fixed in a medial position and (b) in an anteromedial position. The simulation took into account soft tissues preload, muscular tonus and maximal gait load.The maximal stresses observed in the four structural elements (bone, plate, wedge, screws) of an osteotomy with plate in medial position were substantially higher (1.13-2.8 times more) than those observed in osteotomy with an anteromedial plate configuration. An important increase (1.71 times more) of the relative micromotions between the wedge and the bone was also observed. In order to avoid formation of fibrous tissue at the bone wedge interface, the osteotomy should be loaded under 18.8% (approximately 50 kg) of the normal gait load until the osteotomy interfaces union is achieved.     

Development of a three-dimensional finite element model of a human tibia using experimental modal analysis.

The modal analysis of a human tibia consisted of characterizing its dynamic behavior by determining natural frequency, damping ratio and mode shapes. Two methods were used to perform the modal analysis: (1) a finite element method (structural model); (2) an experimental modal analysis (modal model). The experimental modal model was used to optimize the structural model. After optimization, differences in results between the two models were found to be due only to mechanical properties and mass distribution. The influences of boundary conditions and geometric properties (such as inertia and length) were eliminated by the finite element model itself. The percent relative error between the two methods was approximately 3%, corresponding to the standard deviation of the measured frequencies. For the frequency range considered, the mode shapes were bending modes in two different vibration planes (latero-medial and sagittal), with a slight torsion effect due to the twisted geometry of the tibia.     

A three-dimensional finite element model for arterial clamping

Clamp induced injuries of the arterial wall may determine the outcome of surgical procedures. Thus, it is important to investigate the underlying mechanical effects. We present a three-dimensional finite element model, which allows the study of the mechanical response of an artery-treated as a two-layer tube-during arterial clamping. The important residual stresses, which are associated with the load-free configuration of the artery, are also considered. In particular, the finite element analysis of the deformation process of a clamped artery and the associated stress distribution is presented. Within the clamping area a zone of axial tensile peak-stresses was identified, which (may) cause intimal and medial injury. This is an additional injury mechanism, which clearly differs from the commonly assumed wall damage occurring due to compression between the jaws of the clamp. The proposed numerical model provides essential insights into the mechanics of the clamping procedure and the associated injury mechanisms. It allows detailed parameter studies on a virtual clamped artery, which can not be performed with other methodologies. This approach has the potential to identify the most appropriate clamps for certain types of arteries and to guide optimal clamp design.     

Validating a voxel-based finite element model of a human mandible using digital speckle pattern interferometry

Finite element analysis is a powerful tool for predicting the mechanical behaviour of complex biological structures like bones, but to be confident in the results of an analysis, the model should be validated against experimental data. In such validation experiments, the strains in the loaded bones are usually measured with strain gauges glued to the bone surface, but the use of strain gauges on bone can be difficult and provides only very limited data regarding surface strain distributions. This study applies the full-field strain measurement technique of digital speckle pattern interferometry to measure strains in a loaded human mandible and compares the results with the predictions of voxel-based finite element models of the same specimen. It is found that this novel strain measurement technique yields consistent, reliable measurements. Further, strains predicted by the finite element analysis correspond well with the experimental data. These results not only confirm the usefulness of this technique for future validation studies in the field of bone mechanics, but also show that the modelling approach used in this study is able to predict the experimental results very accurately.     

Development of an accurate three-dimensional finite element knee model

This paper presents the development of a detailed articulating three-dimensional finite-element model of the human knee, derived from MRI scan images. The model utilises precise material models and many contact interfaces in order to produce a realistic kinematic response. The behaviour of the model was examined within two fields of biomechanical simulations: general life and car-crash. These simulations were performed with the non-linear explicit dynamic code PAM-SAFE trade mark. The knee model produced results that compared favourably with existing literature. Such a model (together with other joint models that could be constructed using the same techniques) would be a valuable tool for examining new designs of prosthesis and mechanisms of injury.     

A three-dimensional non-linear finite element study of the effect of cement-prosthesis debonding in cemented femoral total hip components.

A three-dimensional non-linear finite element analysis of a cemented femoral component in which the component was partially debonded from the cement mantle was used to assess the effects of debonding on stresses in the cement. Three cases of partial cement-metal debonding were modelled with debonding of the proximal portion of the implant down to a horizontal plane which was 35, 62.5, or 82.5 mm below the prosthesis collar. Each situation was studied under loads simulating both gait and stairclimbing. Also, complete debonding between the implant and the surrounding cement mantle was modeled for loads simulating gait. Under stair climbing loads with partial cement-mental debonding, hoop stresses of 13-18 MPa were observed in the cement at the cement-metal interface at the proximal postero-medial corner of the implant. Similarly, in stair climbing, the maximum principal stresses in the cement were also adjacent to the proximal postero-medial region of the implant. These stresses were compressive and increased from 15 MPa with fully bonded interfaces to 48 MPa with debonding down to 82.5 mm below the prosthesis collar. Under gait loads, complete debonding caused high compressive stresses up to 34.9 MPa in the cement distal to the prosthesis tip. Thus, cement failure subsequent to prosthesis debonding is likely in the proximal region in a partially debonded implant due to stair climbing loads and is likely below the prosthesis tip in a fully debonded implant due to gait loading.