In-vitro results of rapid maxillary expansion on adults compared with finite element simulations

This study was mainly performed to investigate the effects of high maxillary expansion forces on the skull with fresh and thiel-fixed human skulls. The maxillary suture was not weakened except in one experiment. This study compares the strain measured on the zygomatic process of the skull with the results of a finite element model generated for this purpose. An increasing transversal force was applied on the alveolar process (teeth) until rupture. Strain on the zygomatic process, maxilla displacement and the expanding forces were registered.The results of this study show linear material behaviour of the skull before rupture. The highest stress during the experiments and FE simulation was observed on the alveolar process.Conclusions of this study are the necessity of the existence of appropriate models and that female specimens seem to rupture at a lower force than male ones. Both male and female specimens show a similar linear behaviour in the force/strain curve within each gender group. The probability of maxillary suture opening in adults during ultra-rapid maxillary expansion with tooth anchorage is very low. Complications and unwanted rupture could occur.     

Physical validation of a patient-specific contact finite element model of the ankle

A validation study was conducted to determine the extent to which computational ankle contact finite element (FE) results agreed with experimentally measured tibio-talar contact stress. Two cadaver ankles were loaded in separate test sessions, during which ankle contact stresses were measured with a high-resolution (Tekscan) pressure sensor. Corresponding contact FE analyses were subsequently performed for comparison. The agreement was good between FE-computed and experimentally measured mean (3.2% discrepancy for one ankle, 19.3% for the other) and maximum (1.5% and 6.2%) contact stress, as well as for contact area (1.7% and 14.9%). There was also excellent agreement between histograms of fractional areas of cartilage experiencing specific ranges of contact stress. Finally, point-by-point comparisons between the computed and measured contact stress distributions over the articular surface showed substantial agreement, with correlation coefficients of 90% for one ankle and 86% for the other. In the past, general qualitative, but little direct quantitative agreement has been demonstrated with articular joint contact FE models. The methods used for this validation enable formal comparison of computational and experimental results, and open the way for objective statistical measures of regional correlation between FE-computed contact stress distributions from comparison articular joint surfaces (e.g., those from an intact versus those with residual intra-articular fracture incongruity).     

Investigation of effective intrusion and extrusion force for maxillary canine using finite element analysis

Abstract

Orthodontic tooth movement is mainly regulated by the biomechanical responses of loaded periodontal ligament (PDL). We investigated the effective intervals of orthodontic force in pure maxillary canine intrusion and extrusion referring to PDL hydrostatic stress and logarithmic strain. Finite element analysis (FEA) models, including a maxillary canine, PDL and alveolar bone, were constructed based on computed tomography (CT) images of a patient. The material properties of alveolar bone were non-uniformly defined using HU values of CT images; PDL was assumed to be a hyperelastic–viscoelastic material. The compressive stress and tensile stress ranging from 0.47 to 12.8?kPa and 18.8 to 51.2?kPa, respectively, were identified as effective for tooth movement; a strain 0.24% was identified as the lower limit of effective strain. The stress/strain distributions within PDL were acquired in canine intrusion and extrusion using FEA; root apex was the main force-bearing area in intrusion–extrusion movements and was more prone to resorption. Owing to the distinction of PDL biomechanical responses to compression and tension, the effective interval of orthodontic force was substantially lower in canine intrusion (80–90?g) than in canine extrusion (230–260?g). A larger magnitude of force remained applicable in canine extrusion. This study revised and complemented orthodontic biomechanical behaviours of tooth movement with intrusive–extrusive force and could further help optimize orthodontic treatment.     

Variability of a three-dimensional finite element model constructed using magnetic resonance images of a knee for joint contact stress analysis

Magnetic resonance (MR) imaging has been widely used to evaluate the thickness and volume of articular cartilage both in vivo and in vitro. While morphological information on the cartilage can be obtained using MR images, image processing for extracting geometric boundaries of the cartilage may introduce variations in the thickness of the cartilage. To evaluate the variability of using MR images to construct finite element (FE) knee cartilage models, five investigators independently digitized the same set of MR images of a human knee. The topology of cartilage thickness was determined using a minimal distance algorithm. Less than 8 percent variation in cartilage thickness was observed from the digitized data. The effect of changes in cartilage thickness on contact stress analysis was then investigated using five FE models of the knee. One FE model (average FE model) was constructed using the mean values of the digitized contours of the cartilage, and the other four were constructed by varying the thickness of the average FE model by +/- 5 percent and +/- 10 percent, respectively. The results demonstrated that under axial tibial compressive loading (up to 1,400 N), variations of cartilage thickness caused by digitization of MR images may result in a difference of approximately 10 percent in peak contact stresses (surface pressure, von Mises stress, and hydrostatic pressure) in the cartilage. A reduction of cartilage thickness caused increases of contact stresses, while an increase of cartilage thickness reduced contact stresses. Furthermore, the effect of variation of material properties of the cartilage on contact stress analysis was investigated. The peak contact stress increased almost linearly with the Young's modulus of the cartilage. The peak von Mises stress was dramatically reduced when the Poisson,s ratio was increased from 0.05 to 0.49 under an axial compressive load of 1,400 N, while peak hydrostatic pressure was dramatically increased. Peak surface pressure was also increased with the Poisson's ratio, but with a lower magnitude compared to von Mises stress and hydrostatic pressure. In conclusion, the imaging process may cause 10 percent variations in peak contact stress, and the predicted stress distribution is sensitive to the accuracy of the material properties of the cartilage model, especially to the variation of Poisson's ratio.     

An anatomically based patient-specific finite element model of patella articulation: towards a diagnostic tool

A 3D anatomically based patient-specific finite element (FE) model of patello-femoral (PF) articulation is presented to analyse the main features of patella biomechanics, namely, patella tracking (kinematics), quadriceps extensor forces, surface contact and internal patella stresses. The generic geometries are a subset from the model database of the International Union of Physiological Sciences (IUPS) (http://www.physiome.org.nz) Physiome Project with soft tissue derived from the widely used visible human dataset, and the bones digitised from an anatomically accurate physical model with muscle attachment information. The models are customised to patient magnetic resonance images using a variant of free-form deformation, called 'host-mesh' fitting. The continuum was solved using the governing equation of finite elasticity, with the multibody problem coupled through contact mechanics. Additional constraints such as tissue incompressibility are also imposed. Passive material properties are taken from the literature and implemented for deformable tissue with a non-linear micro-structurally based constitutive law. Bone and cartilage are implemented using a 'St-Venant Kirchoff' model suitable for rigid body rotations. The surface fibre directions have been estimated from anatomy images of cadaver muscle dissections and active muscle contraction was based on a steady-state calcium-tension relation. The 3D continuum model of muscle, tendon and bone is compared with experimental results from the literature, and surgical simulations performed to illustrate its clinical assessment capabilities (a Maquet procedure for reducing patella stresses and a vastus lateralis release for a bipartite patella). Finally, the model limitations, issues and future improvements are discussed.     

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.     

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.     

Stent expansion in curved vessel and their interactions: a finite element analysis

Coronary restenosis after angioplasty has been reduced by stenting procedure, but in-stent restenosis (ISR) has not been eliminated yet, especially in tortuous vessels. In this paper, we proposed a finite element method (FEM) to study the expansion of a stent in a curved vessel (the CV model) and their interactions. A model of the same stent in a straight vessel (the SV model) was also studied and mechanical parameters of both models were researched and compared, including final lumen area, tissue prolapse between stent struts and stress distribution. Results show that in the CV model, the vessel was straightened by stenting and a hinge effect can be observed at extremes of the stent. The maximum tissue prolapse of the CV model was more severe (0.079 mm) than the SV model (0.048 mm); and the minimum lumen area of the CV was decreased (6.10 mm(2)), compared to that of the SV model (6.28 mm(2)). Tissue stresses of the highest level were concentrated in the inner curvature of the CV model. The simulations offered some explanations for the clinical results of ISR in curved vessels and gave design suggestions of the stent and balloon for tortuous vessels. This FEM provides a tool to study mechanisms of stents in curved vessels and can improve new stent designs especially for tortuous vessels.     

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.     

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)     

A finite element model of myocardial infarction using a composite material approach

Computational models are effective tools to study cardiac mechanics under normal and pathological conditions. They can be used to gain insight into the physiology of the heart under these conditions while they are adaptable to computer assisted patient-specific clinical diagnosis and therapeutic procedures. Realistic cardiac mechanics models incorporate tissue active/passive response in conjunction with hyperelasticity and anisotropy. Conventional formulation of such models leads to mathematically-complex problems usually solved by custom-developed non-linear finite element (FE) codes. With a few exceptions, such codes are not available to the research community. This article describes a computational cardiac mechanics model developed such that it can be implemented using off-the-shelf FE solvers while tissue pathologies can be introduced in the model in a straight-forward manner. The model takes into account myocardial hyperelasticity, anisotropy, and active contraction forces. It follows a composite tissue modeling approach where the cardiac tissue is decomposed into two major parts: background and myofibers. The latter is modelled as rebars under initial stresses mimicking the contraction forces. The model was applied in silico to study the mechanics of infarcted left ventricle (LV) of a canine. End-systolic strain components, ejection fraction, and stress distribution attained using this LV model were compared quantitatively and qualitatively to corresponding data obtained from measurements as well as to other corresponding LV mechanics models. This comparison showed very good agreement.     

Stability analysis and finite element simulation of bone remodeling model

Bone remodeling is widely viewed as a dynamic process--maintaining bone structure through a balance between the opposed activities of osteoblast and osteoclast cells--in which the stability problem is often pointed out. By an analytical approach, we present a bone remodeling model applied to n unit-elements in order to analyze the stationary states and the condition of their stability. In addition, this theory has been simulated in a computer model using the Finite Element Method (FEM) to show a relationship between the bone remodeling process and the stability analysis.     

Development of a finite element/multi-body model of a newborn infant for restraint analysis and design

A finite element/multi-body model of a newborn infant has been developed by researchers at the University of Windsor. The geometry of this model is derived from a Nita newborn hospital training mannequin. It consists of 17 parts: eight upper and lower limb segments, the torso, head, and a seven-segment neck with seven translational and eight rotational joints. Anthropometry is consistent with hospital growth charts, measurements requested from health professionals and data from the open literature. The biomechanical properties of the model (i.e. joint stiffnesses) are implementations of data identified in the open literature. The model has been validated with respect to studies of the biomechanics of shaken baby syndrome, infant falls and the Q0 anthropomorphic testing device. A significant conclusion of this study is that the kinetics of the Q0 neck is not biofidelic. This model is currently used in an analysis of airway patency for infants in modern automotive child restraints.