Three-dimensional finite element analysis of glass-ceramic dental crowns

Because of the improved esthetic potential of glass-ceramic crowns as dental restorations, they are sometimes preferred over metal-ceramic crowns for restoration of anterior teeth. Because of their relatively high strength, these ceramic crowns are also frequently used for restoration of posterior teeth. However, due to the larger magnitude of biting forces on posterior teeth, intraoral fracture of all-ceramic crowns tends to occur more frequently in posterior crowns (Moffa, 1988). The objective of this study was to determine the relative influence of load orientation and the occlusal thickness of posterior ceramic crowns on the stress distribution which develops under these loading and design conditions. Three-dimensional finite element models for a molar crown were developed to determine the stress distribution under simulated applied loads. Glass-ceramic crowns with occlusal thicknesses of 0.5, 1.5, and 3.0 mm were considered. The largest principal tensile stresses induced in ceramic due to a distributed load of 600 N applied in a cuspal region were approximately 12 and 182 MPa for vertical and horizontal loading orientations, respectively. Stresses which developed in the facial and lingual marginal regions were primarily compressive under vertical loads. However, tensile stresses developed when the load was applied horizontally. Differences in stress distribution within crowns with the three occlusal thicknesses occurred only near the site of loading. Because of the relatively large failure rates of ceramic crowns in the posterior regions, these restorations should be strengthened by improvement in design, composition, and thermal processing conditions. Before any significant progress is made in these areas, these restorations should be used for the anterior teeth. The results of this study suggest that orientation of the applied load has a more important effect on development of large tensile stresses than the occlusal thickness of ceramic.     

Three-dimensional finite element analysis on cochlear implantation electrode insertion

Biomechanics and Modeling in Mechanobiology - Studying the insertion process of cochlear implant (CI) electrode array (EA) is important to ensure successful, sufficient, and safe implantation. A...     

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)     

Three-dimensional finite element analysis of the composite and compomer onlays in primary molars

Resin onlay restoration is an esthetic alternative technique used for restoring extensively damaged primary molars. Understanding the behavior of materials under repeated functional stress and how the stress is transmitted to the remaining tooth structure is important. The aim of this study was to compare stresses in primary molars restored with indirect composite and compomer onlay. 3D frame models of the right mandibular and maxillary primary molars and the alveolar bone were created using computerized tomography images of a six-year-old girl. The enamel and dentine layers above the cement layer were unified to generate onlay restoration, and composite and compomer were used as restorative materials. The vertical occlusal load (100?N) was applied to the teeth in the occlusal contact areas. The von Mises stress distributions and normal stress distributions of the y-axis (parallel to the long axis of tooth) were evaluated. The occlusal stress is transmitted to the cervical part of healthy teeth by spreading it through the enamel layer. The composite and compomer restorative materials exhibited similar stress distribution patterns. An indirect technique creates a structure similar to the original morphological form, and it allows restorations to distribute high occlusal stresses and to minimize possible breakages.     

Three-dimensional finite element analysis of the foot during standing--a material sensitivity study

Information on the internal stresses/strains in the human foot and the pressure distribution at the plantar support interface under loading is useful in enhancing knowledge on the biomechanics of the ankle-foot complex. While techniques for plantar pressure measurements are well established, direct measurement of the internal stresses/strains is difficult. A three-dimensional (3D) finite element model of the human foot and ankle was developed using the actual geometry of the foot skeleton and soft tissues, which were obtained from 3D reconstruction of MR images. Except the phalanges that were fused, the interaction among the metatarsals, cuneiforms, cuboid, navicular, talus, calcaneus, tibia and fibula were defined as contact surfaces, which allow relative articulating movement. The plantar fascia and 72 major ligaments were simulated using tension-only truss elements by connecting the corresponding attachment points on the bone surfaces. The bony and ligamentous structures were embedded in a volume of soft tissues. The encapsulated soft tissue was defined as hyperelastic, while the bony and ligamentous structures were assumed to be linearly elastic. The effects of soft tissue stiffening on the stress distribution of the plantar surface and bony structures during balanced standing were investigated. Increases of soft tissue stiffness from 2 and up to 5 times the normal values were used to approximate the pathologically stiffened tissue behaviour with increasing stages of diabetic neuropathy. The results showed that a five-fold increase in soft tissue stiffness led to about 35% and 33% increase in the peak plantar pressure at the forefoot and rearfoot regions, respectively. This corresponded to about 47% decrease in the total contact area between the plantar foot and the horizontal support surface. Peak bone stress was found at the third metatarsal in all calculated cases with a minimal increase of about 7% with soft tissue stiffening.     

Biomechanical interactions of different mini-plate fixations and maxilla advancements in the Le Fort I Osteotomy: a finite element analysis

This study investigates the biomechanical interaction of different mini-plate fixation types (shapes/sizes and patterns) with segmental advancement levels on the Le Fort I osteotomy using the non-linear finite element (FE) approach. Nine models were generated under a standard 1-piece LeFort I osteotomy for advancement with 3, 6 and 9 mm distances and four mini-plates with three fixation patterns including LL, LI, and II patterns placed on the maxillae models by integrating computed tomography images and computer-aided design system. The axial and oblique occlusal forces were 250 N applied to each premolar/molar and 125 N applied at 30° inclination to the tooth long axis and from palatal to buccal, respectively. The relative micro-movement values between the two maxillary bone segments and maximum mini-plate stress increased obviously with maxilla advancement increment and the increasing trend can be fitted by exponential curve. The corresponding values in II mini-plate fixation presented apparently high values in all simulated cases. The mini-plate stress concentration locations were found at the bending regions to increase high fracture risk. The mini-plate yield strength can be mapped to a critical (limited) advancement for three types of fixations for safe consideration. This study concluded that L-shaped mini-plates with lateral fixation are recommended to provide better stability. The risk for mini-plate fracture and bone relapse increases when maxillary advancement is larger than a critical value of 5 mm in the Le Fort I osteotomy.     

Three-dimensional finite element modelling of muscle forces during mastication

This paper presents a three-dimensional finite element model of human mastication. Specifically, an anatomically realistic model of the masseter muscles and associated bones is used to investigate the dynamics of chewing. A motion capture system is used to track the jaw motion of a subject chewing standard foods. The three-dimensional nonlinear deformation of the masseter muscles are calculated via the finite element method, using the jaw motion data as boundary conditions. Motion-driven muscle activation patterns and a transversely isotropic material law, defined in a muscle-fibre coordinate system, are used in the calculations. Time-force relationships are presented and analysed with respect to different tasks during mastication, e.g. opening, closing, and biting, and are also compared to a more traditional one-dimensional model. The results strongly suggest that, due to the complex arrangement of muscle force directions, modelling skeletal muscles as conventional one-dimensional lines of action might introduce a significant source of error.     

Three-dimensional model of the human craniofacial skeleton: method and preliminary results using finite element analysis

The purpose of this study was to develop a three-dimensional finite element model of the craniofacial skeleton using a dry human skull. The model consisted of 2918 nodes and 1776 solid elements, and was used to investigate the biomechanical effect of a distally directed orthopaedic force on the craniofacial complex. The force was applied at the level of the maxillary first molar. The results indicated that in response to the force system applied: the nasomaxillary complex displaces in a backward and downward direction and rotates in clockwise sense; the nasomaxillary complex, including the zygomatic bone, experiences high stress levels in comparison with those at the remaining bones; the stress distribution in the maxillary basal bone area is relatively uniform; and the stress distribution across the opposing surface of the bony margins of the sutures is non-uniform.     

Three-dimensional finite element analysis of glenoid replacement prostheses: a comparison of keeled and pegged anchorage systems

Glenoid component loosening is the dominant cause of failure in total shoulder arthroplasty. It is presumed that loosening in the glenoid is caused by high stresses in the cement layer. Several anchorage systems have been designed with the aim of reducing the loosening rate, the two major categories being "keeled" fixation and "pegged" fixation. However, no three-dimensional finite element analysis has been performed to quantify the stresses in the cement or to compare the different glenoid prosthesis anchorage systems. The objective of this study was to determine the stresses in the cement layer and surrounding bone for glenoid replacement components. A three-dimensional model of the scapula was generated using CT data for geometry and material property definition. Keeled and pegged designs were inserted into the glenoid, surrounded by a 1-mm layer of bone cement. A 90 deg arm abduction load with a full muscle and joint load was applied, following van der Helm (1994). Deformations of the prosthesis, stresses in the cement, and stresses in the bone were calculated. Stresses were also calculated for a simulated case of rheumatoid arthritis (RA) in which bone properties were modified to reflect that condition. A maximum principal stress-based failure model was used to predict what quantity of the cement is at risk of failure at the levels of stress computed. The prediction is that 94 percent (pegged prosthesis) and 68 percent (keeled prosthesis) of the cement has a greater than 95 percent probability of survival in normal bone. In RA bone, however, the situation is reversed where 86 percent (pegged prosthesis) and 99 percent (keeled prosthesis) of the cement has a greater than 95 percent probability of survival. Bone stresses are shown to be not much affected by the prosthesis design, except at the tip of the central peg or keel. It is concluded that a "pegged" anchorage system is superior for normal bone, whereas a "keeled" anchorage system is superior for RA bone.     

Subject-specific finite element model of the pelvis: development, validation and sensitivity studies

A better understanding of the three-dimensional mechanics of the pelvis, at the patient-specific level, may lead to improved treatment modalities. Although finite element (FE) models of the pelvis have been developed, validation by direct comparison with subject-specific strains has not been performed, and previous models used simplifying assumptions regarding geometry and material properties. The objectives of this study were to develop and validate a realistic FE model of the pelvis using subject-specific estimates of bone geometry, location-dependent cortical thickness and trabecular bone elastic modulus, and to assess the sensitivity of FE strain predictions to assumptions regarding cortical bone thickness as well as bone and cartilage material properties. A FE model of a cadaveric pelvis was created using subject-specific computed tomography image data. Acetabular loading was applied to the same pelvis using a prosthetic femoral stem in a fashion that could be easily duplicated in the computational model. Cortical bone strains were monitored with rosette strain gauges in ten locations on the left hemipelvis. FE strain predictions were compared directly with experimental results for validation. Overall, baseline FE predictions were strongly correlated with experimental results (r2=0.824), with a best-fit line that was not statistically different than the line y=x (experimental strains = FE predicted strains). Changes to cortical bone thickness and elastic modulus had the largest effect on cortical bone strains. The FE model was less sensitive to changes in all other parameters. The methods developed and validated in this study will be useful for creating and analyzing patient-specific FE models to better understand the biomechanics of the pelvis.     

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.     

Effects of several temporomandibular disorders on the stress distributions of temporomandibular joint: a finite element analysis

The aim of this study was to evaluate stress distributions in the temporomandibular joints (TMJs) with temporomandibular disorders (TMDs) for comparison with healthy TMJs. A model of mandible and normal TMJs was developed according to CT images. The interfaces between the discs and the articular cartilages were treated as contact elements. Nonlinear cable elements were used to simulate disc attachments. Based on this model, seven models of various TMDs were established. The maximum stresses of the discs with anterior, posterior, medial and lateral disc displacement (ADD, PDD, MDD and LDD) were 12.09, 9.33, 10.71 and 6.07 times magnitude of the identically normal disc, respectively. The maximum stresses of the posterior articular eminences in ADD, PDD, MDD, LDD, relaxation of posterior attachments and disc perforation models were 21, 59, 46, 21, 13 and 15 times greater than the normal model, respectively. TMDs could cause increased stresses in the discs and posterior articular eminences.     

Statistical methods in finite element analysis

Finite element analysis (FEA) is a commonly used tool within many areas of engineering and can provide useful information in structural analysis of mechanical systems. However, most analyses within the field of biomechanics usually take no account either of the wide variation in material properties and geometry that may occur in natural tissues or manufacturing imperfections in synthetic materials. This paper discusses two different methods of incorporating uncertainty in FE models. The first, Taguchi's robust parameter design, uses orthogonal matrices to determine how to vary the parameters in a series of FE models, and provides information on the sensitivity of a model to input parameters. The second, probabilistic analysis, enables the distribution of a response variable to be determined from the distributions of the input variables. The methods are demonstrated using a simple example of an FE model of a beam that is assigned material properties and geometry over a range similar to an orthopaedic fixation plate. In addition to showing how each method may be used on its own, we also show how computational effort may be minimised by first identifying the most important input variables before determining the effects of imprecision.     

In vivo characterization of mechanical tissue properties of internal organs using endoscopic microscopy and inverse finite element analysis

Knowledge about mechanical tissue properties is required for functional modelling and simulating of tissue and organ responses to external mechanical stress. To get the right properties especially for functional modelling of organs, tissue properties have to be determined in vivo. There are only few described methods for characterization of internal organ's tissue mechanics that can be applied in vivo. We introduce and evaluate a method to determine mechanical tissue properties, especially those of lung tissue, endoscopically. Inverse finite element analysis (utilizing a Neo-Hookean model for hyperelastic materials) and image processing algorithms are used to determine the shear modulus of a soft tissue. The resulting values for shear moduli were normally distributed. The shear modulus of the artificial tissue sample was determined with a relative error of 0.47% compared to the value obtained by uniaxial tensile test.     

Three-dimensional simulation of calcium waves and contraction in cardiomyocytes using the finite element method

To investigate the characteristics and underlying mechanisms of Ca²⁺ wave propagation, we developed a three-dimensional (3-D) simulator of cardiac myocytes, in which the sarcolemma, myofibril, and Z-line structure with Ca²⁺ release sites were modeled as separate structures using the finite element method. Similarly to previous studies, we assumed that Ca²⁺ diffusion from one release site to another and Ca²⁺-induced Ca²⁺ release were the basic mechanisms, but use of the finite element method enabled us to simulate not only the wave propagation in 3-D space but also the active shortening of the myocytes. Therefore, in addition to the dependence of the Ca²⁺ wave propagation velocity on the sarcoplasmic reticulum Ca²⁺ content and affinity of troponin C for Ca²⁺, we were able to evaluate the influence of active shortening on the propagation velocity. Furthermore, if the initial Ca²⁺ release took place in the proximity of the nucleus, spiral Ca²⁺ waves evolved and spread in a complex manner, suggesting that this phenomenon has the potential for arrhythmogenicity. The present 3-D simulator, with its ability to study the interaction between Ca²⁺ waves and contraction, will serve as a useful tool for studying the mechanism of this complex phenomenon. cardiac muscle cell; excitation-contraction coupling; mechanoelectrical feedback; spiral wave; arrhythmia     

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°).     

Implication of obesity on motion,posture and internal stress of the foot: an experimental and finite element analysis

Obesity causes increased loading on the foot which can damage the soft tissue and bone ultimately leading to foot problems. Experimental and computational methods were used to analyse the chain of biomechanical changes in the lower limb due to obesity. The experimental study shows some changes in foot posture and gait where obese subjects were more likely to have pronated feet, smaller joint angles in the sagittal and frontal planes, smaller cadence, and smaller stride length. Anatomically correct finite element models generated on obese subjects showed increased and altered internal and plantar stress. Altered foot posture was identified as a key indicator of increased internal stress indicating the importance of foot posture correction.     

Potential of finite element analysis for cephalometric investigation

In view of concerns relating to the validity of traditional cephalometric appraisals, we undertook this study to use a potentially powerful method, finite element analysis, to compare cephalometric changes between two samples (subjects with and without orthodontic treatment). The derived data show varying sample contrasts depending on the particular finite element array included in the analysis. Thus, although finite element analysis facilitates rigorous morphometric analysis, further investigation is required before it can be applied universally in cephalometric studies.