A 3-D finite element study of the influence of crown-implant ratio on stress distribution

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Root canal filling: fracture strength of fiber-reinforced composite-restored roots and finite element analysis

The aims of this study were to evaluate the effect of root canal filling techniques on root fracture resistance and to analyze, by finite element analysis (FEA), the expansion of the endodontic sealer in two different root canal techniques. Thirty single-rooted human teeth were instrumented with rotary files to a standardized working length of 14 mm. The specimens were embedded in acrylic resin using plastic cylinders as molds, and allocated into 3 groups (n=10): G(lateral) - lateral condensation; G(single-cone) - single cone; G(tagger) - Tagger's hybrid technique. The root canals were prepared to a length of 11 mm with the #3 preparation bur of a tapered glass fiber-reinforced composite post system. All roots received glass fiber posts, which were adhesively cemented and a composite resin core was built. All groups were subjected to a fracture strength test (1 mm/min, 45°). Data were analyzed statistically by one-way ANOVA with a significance level of 5%. FEA was performed using two models: one simulated lateral condensation and Tagger's hybrid technique, and the other one simulated the single-cone technique. The second model was designed with an amount of gutta-percha two times smaller and a sealer layer two times thicker than the first model. The results were analyzed using von Mises stress criteria. One-way ANOVA indicated that the root canal filling technique affected the fracture strength (p=0.004). The G(lateral) and G(tagger) produced similar fracture strength values, while G(single-cone) showed the lowest values. The FEA showed that the single-cone model generated higher stress in the root canal walls. Sealer thickness seems to influence the fracture strength of restored endodontically treated teeth.

Complex modal analysis of a vertical rotor through finite element method

Natural frequencies were analyzed (axial, torsional and flexural) and frequency response of a vertical rotor with a hard disk at the edge through the classical and complex modal analysis. The mathematical modeling was based on the theory of Euler-Bernoulli beam. The equation that rules the movement was obtained through the Lagrangian formulation. The model considered the effects of bending, torsion and axial deformation of the shaft, besides the gravitational and gyroscopic effects. The finite element method was used to discretize the structure into hollow cylindrical elements with 12 degrees of freedom. Mass, stiffness and gyroscopic matrices were explained consistently. This type of tool, based on the use of complex coordinates to describe the dynamic behavior of rotating shaft, allows the decomposition of the system in two submodes, backward and forward. Thus, it is possible to clearly visualize that the orbit and direction of the precessional motion around the line of the rotating shaft is not deformed. A finite element program was developed using Matlab ®, and numerical simulations were performed to validate this model.

Three-Dimensional Finite Element Analysis of Stress Distribution on Different Bony Ridges With Different Lengths of Morse Taper Implants and Prosthesis Dimensions

This finite element analysis (FEA) compared stress distribution on different bony ridges rehabilitated with different lengths of morse taper implants, varying dimensions of metal-ceramic crowns to maintain the occlusal alignment. Three-dimensional FE models were designed representing a posterior left side segment of the mandible: group control, 3 implants of 11 mm length; group 1, implants of 13 mm, 11 mm and 5 mm length; group 2, 1 implant of 11 mm and 2 implants of 5 mm length; and group 3, 3 implants of 5 mm length. The abutments heights were 3.5 mm for 13- and 11-mm implants (regular), and 0.8 mm for 5-mm implants (short). Evaluation was performed on Ansys software, oblique loads of 365N for molars and 200N for premolars. There was 50% higher stress on cortical bone for the short implants than regular implants. There was 80% higher stress on trabecular bone for the short implants than regular implants. There was higher stress concentration on the bone region of the short implants neck. However, these implants were capable of dissipating the stress to the bones, given the applied loads, but achieving near the threshold between elastic and plastic deformation to the trabecular bone. Distal implants and/or with biggest occlusal table generated greatest stress regions on the surrounding bone. It was concluded that patients requiring short implants associated with increased proportions implant prostheses need careful evaluation and occlusal adjustment, as a possible overload in these short implants, and even in regular ones, can generate stress beyond the physiological threshold of the surrounding bone, compromising the whole system.

Assessment of nose protector for sport activities: finite element analysis

There has been a significant increase in the number of facial fractures stemming from sport activities in recent years, with the nasal bone one of the most affected structures. Researchers recommend the use of a nose protector, but there is no standardization regarding the material employed. Clinical experience has demonstrated that a combination of a flexible and rigid layer of ethylene vinyl acetate (EVA) offers both comfort and safety to practitioners of sports. The aim of the present study was the investigation into the stresses generated by the impact of a rigid body on the nasal bone on models with and without an EVA protector. For such, finite element analysis was employed. A craniofacial model was constructed from images obtained through computed tomography. The nose protector was modeled with two layers of EVA (1 mm of rigid EVA over 2 mm of flexible EVA), following the geometry of the soft tissue. Finite element analysis was performed using the LS Dyna program. The bone and rigid EVA were represented as elastic linear material, whereas the soft tissues and flexible EVA were represented as hyperelastic material. The impact from a rigid sphere on the frontal region of the face was simulated with a constant velocity of 20 m s-1 for 9.1 mu s. The model without the protector served as the control. The distribution of maximal stress of the facial bones was recorded. The maximal stress on the nasal bone surpassed the breaking limit of 0.130.34 MPa on the model without a protector, while remaining below this limit on the model with the protector. Thus, the nose protector made from both flexible and rigid EVA proved effective at protecting the nasal bones under high-impact conditions.

FINITE ELEMENT ANALYSIS OF MAXILLARY BONE STRESS CAUSED BY ARAMANY CLASS IV OBTURATOR PROSTHESES

Statement of problem. The retention of an Aramany Class IV removable partial dental prosthesis can be compromised by a lack of support. The biomechanics of this obturator prosthesis result in an unusual stress distribution on the residual maxillary bone. Purpose. This study evaluated the biomechanics of an Aramany Class IV obturator prosthesis with finite element analysis and a digital 3-dimensional (3-D) model developed from a computed tomography scan; bone stress was evaluated according to the load placed on the prosthesis. Material and methods. A 3-D model of an Aramany Class IV maxillary resection and prosthesis was constructed. This model was used to develop a finite element mesh. A 120 N load was applied to the occlusal and incisal platforms corresponding to the prosthetic teeth. Qualitative analysis was based on the scale of maximum principal stress; values obtained through quantitative analysis were expressed in MPa. Results. Under posterior load, tensile and compressive stresses were observed; the tensile stress was greater than the compressive stress, regardless of the bone region, and the greatest compressive stress was observed on the anterior palate near the midline. Under an anterior load, tensile stress was observed in all of the evaluated bone regions; the tensile stress was greater than the compressive stress, regardless of the bone region. Conclusions. The Aramany Class IV obturator prosthesis tended to rotate toward the surgical resection when subjected to posterior or anterior loads. The amount of tensile and compressive stress caused by the Aramany Class IV obturator prosthesis did not exceed the physiological limits of the maxillary bone tissue. (J Prosthet Dent 2012;107:336-342)

Finite element methods for the Stokes system based on a Zienkiewicz type N-simplex

Hermite interpolation is increasingly showing to be a powerful numerical solution tool, as applied to different kinds of second order boundary value problems. In this work we present two Hermite finite element methods to solve viscous incompressible flows problems, in both two- and three-dimension space. In the two-dimensional case we use the Zienkiewicz triangle to represent the velocity field, and in the three-dimensional case an extension of this element to tetrahedra, still called a Zienkiewicz element. Taking as a model the Stokes system, the pressure is approximated with continuous functions, either piecewise linear or piecewise quadratic, according to the version of the Zienkiewicz element in use, that is, with either incomplete or complete cubics. The methods employ both the standard Galerkin or the Petrov–Galerkin formulation first proposed in Hughes et al. (1986) [18], based on the addition of a balance of force term. A priori error analyses point to optimal convergence rates for the PG approach, and for the Galerkin formulation too, at least in some particular cases. From the point of view of both accuracy and the global number of degrees of freedom, the new methods are shown to have a favorable cost-benefit ratio, as compared to velocity Lagrange finite elements of the same order, especially if the Galerkin approach is employed.

Finite element models of coseismic deformation due to the 2009 L'Aquila (Italy) and 2008 Wenchuan(China) earthquakes

The topic of my Ph.D. thesis is the finite element modeling of coseismic deformation imaged by DInSAR and GPS data. I developed a method to calculate synthetic Green functions with finite element models (FEMs) and then use linear inversion methods to determine the slip distribution on the fault plane. The method is applied to the 2009 L’Aquila Earthquake (Italy) and to the 2008 Wenchuan earthquake (China). I focus on the influence of rheological features of the earth's crust by implementing seismic tomographic data and the influence of topography by implementing Digital Elevation Models (DEM) layers on the FEMs. Results for the L’Aquila earthquake highlight the non-negligible influence of the medium structure: homogeneous and heterogeneous models show discrepancies up to 20% in the fault slip distribution values. Furthermore, in the heterogeneous models a new area of slip appears above the hypocenter. Regarding the 2008 Wenchuan earthquake, the very steep topographic relief of Longmen Shan Range is implemented in my FE model. A large number of DEM layers corresponding to East China is used to achieve the complete coverage of the FE model. My objective was to explore the influence of the topography on the retrieved coseismic slip distribution. The inversion results reveals significant differences between the flat and topographic model. Thus, the flat models frequently adopted are inappropriate to represent the earth surface topographic features and especially in the case of the 2008 Wenchuan earthquake.

Mixed finite-element methods for elliptic convection-dominated problems arising in semiconductor physics

In this work we develop and analyze an adaptive numerical scheme for simulating a class of macroscopic semiconductor models. At first the numerical modelling of semiconductors is reviewed in order to classify the Energy-Transport models for semiconductors that are later simulated in 2D. In this class of models the flow of charged particles, that are negatively charged electrons and so-called holes, which are quasi-particles of positive charge, as well as their energy distributions are described by a coupled system of nonlinear partial differential equations. A considerable difficulty in simulating these convection-dominated equations is posed by the nonlinear coupling as well as due to the fact that the local phenomena such as "hot electron effects" are only partially assessable through the given data. The primary variables that are used in the simulations are the particle density and the particle energy density. The user of these simulations is mostly interested in the current flow through parts of the domain boundary - the contacts. The numerical method considered here utilizes mixed finite-elements as trial functions for the discrete solution. The continuous discretization of the normal fluxes is the most important property of this discretization from the users perspective. It will be proven that under certain assumptions on the triangulation the particle density remains positive in the iterative solution algorithm. Connected to this result an a priori error estimate for the discrete solution of linear convection-diffusion equations is derived. The local charge transport phenomena will be resolved by an adaptive algorithm, which is based on a posteriori error estimators. At that stage a comparison of different estimations is performed. Additionally a method to effectively estimate the error in local quantities derived from the solution, so-called "functional outputs", is developed by transferring the dual weighted residual method to mixed finite elements. For a model problem we present how this method can deliver promising results even when standard error estimator fail completely to reduce the error in an iterative mesh refinement process.

Finite-element discretization of 3D energy-transport equations for semiconductors

In this thesis a mathematical model was derived that describes the charge and energy transport in semiconductor devices like transistors. Moreover, numerical simulations of these physical processes are performed. In order to accomplish this, methods of theoretical physics, functional analysis, numerical mathematics and computer programming are applied. After an introduction to the status quo of semiconductor device simulation methods and a brief review of historical facts up to now, the attention is shifted to the construction of a model, which serves as the basis of the subsequent derivations in the thesis. Thereby the starting point is an important equation of the theory of dilute gases. From this equation the model equations are derived and specified by means of a series expansion method. This is done in a multi-stage derivation process, which is mainly taken from a scientific paper and which does not constitute the focus of this thesis. In the following phase we specify the mathematical setting and make precise the model assumptions. Thereby we make use of methods of functional analysis. Since the equations we deal with are coupled, we are concerned with a nonstandard problem. In contrary, the theory of scalar elliptic equations is established meanwhile. Subsequently, we are preoccupied with the numerical discretization of the equations. A special finite-element method is used for the discretization. This special approach has to be done in order to make the numerical results appropriate for practical application. By a series of transformations from the discrete model we derive a system of algebraic equations that are eligible for numerical evaluation. Using self-made computer programs we solve the equations to get approximate solutions. These programs are based on new and specialized iteration procedures that are developed and thoroughly tested within the frame of this research work. Due to their importance and their novel status, they are explained and demonstrated in detail. We compare these new iterations with a standard method that is complemented by a feature to fit in the current context. A further innovation is the computation of solutions in three-dimensional domains, which are still rare. Special attention is paid to applicability of the 3D simulation tools. The programs are designed to have justifiable working complexity. The simulation results of some models of contemporary semiconductor devices are shown and detailed comments on the results are given. Eventually, we make a prospect on future development and enhancements of the models and of the algorithms that we used.

Fluid-structure interaction by a coupled lattice Boltzmann-finite element approach

In this thesis, a strategy to model the behavior of fluids and their interaction with deformable bodies is proposed. The fluid domain is modeled by using the lattice Boltzmann method, thus analyzing the fluid dynamics by a mesoscopic point of view. It has been proved that the solution provided by this method is equivalent to solve the Navier-Stokes equations for an incompressible flow with a second-order accuracy. Slender elastic structures idealized through beam finite elements are used. Large displacements are accounted for by using the corotational formulation. Structural dynamics is computed by using the Time Discontinuous Galerkin method. Therefore, two different solution procedures are used, one for the fluid domain and the other for the structural part, respectively. These two solvers need to communicate and to transfer each other several information, i.e. stresses, velocities, displacements. In order to guarantee a continuous, effective, and mutual exchange of information, a coupling strategy, consisting of three different algorithms, has been developed and numerically tested. In particular, the effectiveness of the three algorithms is shown in terms of interface energy artificially produced by the approximate fulfilling of compatibility and equilibrium conditions at the fluid-structure interface. The proposed coupled approach is used in order to solve different fluid-structure interaction problems, i.e. cantilever beams immersed in a viscous fluid, the impact of the hull of the ship on the marine free-surface, blood flow in a deformable vessels, and even flapping wings simulating the take-off of a butterfly. The good results achieved in each application highlight the effectiveness of the proposed methodology and of the C++ developed software to successfully approach several two-dimensional fluid-structure interaction problems.

Finite element modeling of the cement-bone interface beneath the tibial tray of a total knee arthroplasty: Study of the stress shielding effect (modellazione agli elementi finiti dell'interfaccia cemento-osso sotto il piatto tibiale di protesi totale di ginocchio: Analisi dello stress shielding)

The goal of this thesis was the study of the cement-bone interface in the tibial component of a cemented total knee prosthesis. One of the things you can see in specimens after in vivo service is that resorption of bone occurs in the interdigitated region between bone and cement. A stress shielding effect was investigated as a cause to explain bone resorption. Stress shielding occurs when bone is loaded less than physiological and therefore it starts remodeling according to the new loading conditions. µCT images were used to obtain 3D models of the bone and cement structure and a Finite Element Analysis was used to simulate different kind of loads. Resorption was also simulated by performing erosion operations in the interdigitated bone region. Finally, 4 models were simulated: bone (trabecular), bone with cement, and two models of bone with cement after progressive erosions of the bone. Compression, tension and shear test were simulated for each model in displacement-control until 2% of strain. The results show how the principal strain and Von Mises stress decrease after adding the cement on the structure and after the erosion operations. These results show that a stress shielding effect does occur and rises after resorption starts.

Influence of Smoothing on Voxel-Based Mesh Accuracy in Micro-Finite Element

The interest in automatic volume meshing for finite element analysis (FEA) has grown more since the appearance of microfocus CT (μCT), due to its high resolution, which allows for the assessment of mechanical behaviour at a high precision. Nevertheless, the basic meshing approach of generating one hexahedron per voxel produces jagged edges. To prevent this effect, smoothing algorithms have been introduced to enhance the topology of the mesh. However, whether smoothing also improves the accuracy of voxel-based meshes in clinical applications is still under question. There is a trade-off between smoothing and quality of elements in the mesh. Distorted elements may be produced by excessive smoothing and reduce accuracy of the mesh. In the present work, influence of smoothing on the accuracy of voxel-based meshes in micro-FE was assessed. An accurate 3D model of a trabecular structure with known apparent mechanical properties was used as a reference model. Virtual CT scans of this reference model (with resolutions of 16, 32 and 64 μm) were then created and used to build voxel-based meshes of the microarchitecture. Effects of smoothing on the apparent mechanical properties of the voxel-based meshes as compared to the reference model were evaluated. Apparent Young’s moduli of the smooth voxel-based mesh were significantly closer to those of the reference model for the 16 and 32 μm resolutions. Improvements were not significant for the 64 μm, due to loss of trabecular connectivity in the model. This study shows that smoothing offers a real benefit to voxel-based meshes used in micro-FE. It might also broaden voxel-based meshing to other biomechanical domains where it was not used previously due to lack of accuracy. As an example, this work will be used in the framework of the European project ContraCancrum, which aims at providing a patient-specific simulation of tumour development in brain and lungs for oncologists. For this type of clinical application, such a fast, automatic, and accurate generation of the mesh is of great benefit.

Finite Element Modeling of Shell Shape in the Freshwater Turtle Pseudemys concinna Reveals a Trade-Off between Mechanical Strength and Hydrodynamic Efficiency

Aquatic species can experience different selective pressures on morphology in different flow regimes. Species inhabiting lotic regimes often adapt to these conditions by evolving low-drag (i.e., streamlined) morphologies that reduce the likelihood of dislodgment or displacement. However, hydrodynamic factors are not the only selective pressures influencing organismal morphology and shapes well suited to flow conditions may compromise performance in other roles. We investigated the possibility of morphological trade-offs in the turtle Pseudemys concinna. Individuals living in lotic environments have flatter, more streamlined shells than those living in lentic environments; however, this flatter shape may also make the shells less capable of resisting predator-induced loads. We tested the idea that ‘‘lotic’’ shell shapes are weaker than ‘‘lentic’’ shell shapes, concomitantly examining effects of sex. Geometric morphometric data were used to transform an existing finite element shell model into a series of models corresponding to the shapes of individual turtles. Models were assigned identical material properties and loaded under identical conditions, and the stresses produced by a series of eight loads were extracted to describe the strength of the shells. ‘‘Lotic’’ shell shapes produced significantly higher stresses than ‘‘lentic’’ shell shapes, indicating that the former is weaker than the latter. Females had significantly stronger shell shapes than males, although these differences were less consistent than differences between flow regimes. We conclude that, despite the potential for many-to-one mapping of shell shape onto strength, P. concinna experiences a trade-off in shell shape between hydrodynamic and mechanical performance. This trade-off may be evident in many other turtle species or any other aquatic species that also depend on a shell for defense. However, evolution of body size may provide an avenue of escape from this trade-off in some cases, as changes in size can drastically affect mechanical performance while having little effect on hydrodynamic performance.

Evaluation of the Design of Komendant's Cycloidal, Thin Shell Roofs At the Kimbell Art Museum Using Finite Element Analysis

The analysis of Komendant's design of the Kimbell Art Museum was carried out in order to determine the effectiveness of the ring beams, edge beams and prestressing in the shells of the roof system. Finite element analysis was not available to Komendant or other engineers of the time to aid them in the design and analysis. Thus, the use of this tool helped to form a new perspective on the Kimbell Art Museum and analyze the engineer's work. In order to carry out the finite element analysis of Kimbell Art Museum, ADINA finite element analysis software was utilized. Eight finite element models (FEM-1 through FEM-8) of increasing complexity were created. The results of the most realistic model, FEM-8, which included ring beams, edge beams and prestressing, were compared to Komendant's calculations. The maximum deflection at the crown of the mid-span surface of -0.1739 in. in FEM-8 was found to be larger than Komendant's deflection in the design documents before the loss in prestressing force (-0.152 in.) but smaller than his prediction after the loss in prestressing force (-0.3814 in.). Komendant predicted a larger longitudinal stress of -903 psi at the crown (vs. -797 psi in FEM-8) and 37 psi at the edge (vs. -347 psi in FEM-8). Considering the strength of concrete of 5000 psi, the difference in results is not significant. From the analysis it was determined that both FEM-5, which included prestressing and fixed rings, and FEM-8 can be successfully and effectively implemented in practice. Prestressing was used in both models and thus served as the main contribution to efficiency. FEM-5 showed that ring and edge beams can be avoided, however an architect might find them more aesthetically appropriate than rigid walls.