Two-dimensional isostatic meshes in the finite element method

In a Finite Element (FE) analysis of elastic solids several items are usually considered, namely, type and shape of the elements, number of nodes per element, node positions, FE mesh, total number of degrees of freedom (dot) among others. In this paper a method to improve a given FE mesh used for a particular analysis is described. For the improvement criterion different objective functions have been chosen (Total potential energy and Average quadratic error) and the number of nodes and dof's of the new mesh remain constant and equal to the initial FE mesh. In order to find the mesh producing the minimum of the selected objective function the steepest descent gradient technique has been applied as optimization algorithm. However this efficient technique has the drawback that demands a large computation power. Extensive application of this methodology to different 2-D elasticity problems leads to the conclusion that isometric isostatic meshes (ii-meshes) produce better results than the standard reasonably initial regular meshes used in practice. This conclusion seems to be independent on the objective function used for comparison. These ii-meshes are obtained by placing FE nodes along the isostatic lines, i.e. curves tangent at each point to the principal direction lines of the elastic problem to be solved and they should be regularly spaced in order to build regular elements. That means ii-meshes are usually obtained by iteration, i.e. with the initial FE mesh the elastic analysis is carried out. By using the obtained results of this analysis the net of isostatic lines can be drawn and in a first trial an ii-mesh can be built. This first ii-mesh can be improved, if it necessary, by analyzing again the problem and generate after the FE analysis the new and improved ii-mesh. Typically, after two first tentative ii-meshes it is sufficient to produce good FE results from the elastic analysis. Several example of this procedure are presented.

Analysis of the moisture effect on the mechanical behaviour of an adhesively bonded joint under proportional multi-axial loads

The objective of the study is to identify the 3D behaviour of an adhesive in an assembly, and to take into account the effect of ageing in a marine environment. To that end, three different tests were employed. Gravimetric analyses were used to determine the water diffusion kinetics in the adhesive. Bulk tensile tests were performed to highlight the effects of humid ageing on the adhesive behaviour. Modified Arcan tests were performed for several ageing times to obtain the experimental database which was necessary to identify constitutive models. A Mahnken-Schlimmer type model was determined for the unaged state according to a procedure developed in a previous study. This identification used inverse techniques. It was based on the unaged modified Arcan results and on a coupling between an optimisation routine and finite-element analysis. Then, a global inverse identification procedure was developed. Its aim was to relate the unaged parameters to the moisture concentration and overcome the difficulties usually associated with ageing of bonded assemblies in a humid environment: a non-uniformity of the stress state and a gradient of mechanical properties in the adhesive. This procedure was similar to the one used in the first part but needed modified Arcan results for several ageing times. It also required an initial assumption for the evolution of the Mahnken-Schlimmer parameters with the moisture concentration.

Modelling adhesive joints with cohesive zone models: effect of the cohesive law shape of the adhesive layer

Adhesively-bonded joints are extensively used in several fields of engineering. Cohesive Zone Models (CZM) have been used for the strength prediction of adhesive joints, as an add-in to Finite Element (FE) analyses that allows simulation of damage growth, by consideration of energetic principles. A useful feature of CZM is that different shapes can be developed for the cohesive laws, depending on the nature of the material or interface to be simulated, allowing an accurate strength prediction. This work studies the influence of the CZM shape (triangular, exponential or trapezoidal) used to model a thin adhesive layer in single-lap adhesive joints, for an estimation of its influence on the strength prediction under different material conditions. By performing this study, guidelines are provided on the possibility to use a CZM shape that may not be the most suited for a particular adhesive, but that may be more straightforward to use/implement and have less convergence problems (e.g. triangular shaped CZM), thus attaining the solution faster. The overall results showed that joints bonded with ductile adhesives are highly influenced by the CZM shape, and that the trapezoidal shape fits best the experimental data. Moreover, the smaller is the overlap length (LO), the greater is the influence of the CZM shape. On the other hand, the influence of the CZM shape can be neglected when using brittle adhesives, without compromising too much the accuracy of the strength predictions.

Optimization study of hybrid spot-welded/bonded single-lap joints

Joining of components with structural adhesives is currently one of the most widespread techniques for advanced structures (e.g., aerospace or aeronautical). Adhesive bonding does not involve drilling operations and it distributes the load over a larger area than mechanical joints. However, peak stresses tend to develop near the overlap edges because of differential straining of the adherends and load asymmetry. As a result, premature failures can be expected, especially for brittle adhesives. Moreover, bonded joints are very sensitive to the surface treatment of the material, service temperature, humidity and ageing. To surpass these limitations, the combination of adhesive bonding with spot-welding is a choice to be considered, adding a few advantages like superior static strength and stiffness, higher peeling and fatigue strength and easier fabrication, as fixtures during the adhesive curing are not needed. The experimental and numerical study presented here evaluates hybrid spot-welded/bonded single-lap joints in comparison with the purely spot-welded and bonded equivalents. A parametric study on the overlap length (LO) allowed achieving different strength advantages, up to 58% compared to spot-welded joints and 24% over bonded joints. The Finite Element Method (FEM) and Cohesive Zone Models (CZM) for damage growth were also tested in Abaqus® to evaluate this technique for strength prediction, showing accurate estimations for all kinds of joints.

Numerical evaluation of three-dimensional scarf repairs in carbon-epoxy structures

The widespread employment of carbon-epoxy laminates in high responsibility and severely loaded applications introduces an issue regarding their handling after damage. Repair of these structures should be evaluated, instead of their disposal, for cost saving and ecological purposes. Under this perspective, the availability of efficient repair methods is essential to restore the strength of the structure. The development and validation of accurate predictive tools for the repairs behaviour are also extremely important, allowing the reduction of costs and time associated to extensive test programmes. Comparing with strap repairs, scarf repairs have the advantages of a higher efficiency and the absence of aerodynamic disturbance. This work reports on a numerical study of the tensile behaviour of three-dimensional scarf repairs in carbon-epoxy structures, using a ductile adhesive (Araldite® 2015). The finite elements analysis was performed in ABAQUS® and Cohesive Zone Modelling was used for the simulation of damage onset and growth in the adhesive layer. Trapezoidal cohesive laws in each pure mode were used to account for the ductility of the specific adhesive mentioned. A parametric study was performed on the repair width and scarf angle. The use of over-laminating plies covering the repaired region at the outer or both repair surfaces was also tested as an attempt to increase the repairs efficiency. The obtained results allowed the proposal of design principles for repairing composite structures.

Biomechanics Studies in Dentistry: Bioengineering Applied in Oral Implantology

The application of engineering knowledge in dentistry has helped the understanding of biomechanics aspects related to osseointegrated implants. Several techniques have been used to evaluate the biomechanical load oil implants comprising the use of photoelastic stress analysis, finite element stress analysis, and strain-gauge analysis. Therefore, the purpose of this Study was to describe engineering methods used in dentistry to evaluate the biomechanical behavior of osseointegrated implants. Photoelasticity provides good qualitative information oil the overall location and concentration of stresses but produces limited quantitative information. The method serves as ail important tool for determining the critical stress points in a material and is often used for determining stress concentration factors in irregular geometries. The application of strain-gauge method oil dental implants is based oil the use of electrical resistance strain gauges and its associated equipment and provides both in vitro and vivo measurements strains under static and dynamic loads. However, strain-gauge method provides only the data regarding strain at the gauge. Finite element analysis can Simulate stress using a computer-created model to calculate stress, strain, and displacement. Such analysis has the advantage of allowing several conditions to be changed easily and allows measurement of stress distribution around implants at optional points that are difficult to examine clinically All the 3 methodologies call be useful to evaluate biomechanical implant behavior close to the clinical condition but the researcher should have enough knowledge in model fabrication (experimental delineation) and results analysis.

Influence of convergence angle of tooth preparation on the fracture resistance of Y-TZP-based all-ceramic restorations

Objective: To investigate the influence of the convergence angle of tooth preparation on the fracture load of Y-TZP-based ceramic (YZ-Vita YZ) substructure (SB) veneered with a feldspathic porcelain (VM9-Vita VM9). Methods: Finite element stress analysis (FEA) was performed to examine the stress distribution of the system. Eighty YZ SB were fabricated using a CAD-CAM system and divided into four groups (n = 20), according to the total occlusal convergence (TOC) angle: G6-6° TOC; G12-12° TOC; G20-20° TOC; and G20MOD-20° TOC with modified SB. All SB were veneered with VM9, cemented in a fiber reinforced epoxy resin die, and loaded to failure. Half of the specimens from each group (n = 10) were cyclic fatigued (106 cycles) before testing. Failure analysis was performed to determine the fracture origin. Data were statistically analyzed using Anova and Tukey's tests (α = 0.05). Results: The greatest mean load to fracture value was found for the G20MOD, which was predicted by the FEA. Cyclic fatigue did not significantly affect the load of fracture. Catastrophic failure originating from the internal occlusal surface of the SB was the predominant failure mode, except for G20MOD. Significance: The YZ-VM9 restorations resisted greater compression load than the usual physiological occlusal load, regardless of the TOC angle of preparations. Yet, the G20MOD design produced the best performance among the experimental conditions evaluated. © 2013 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

Influência da angulação do preparo coronário na resistência à fratura de restaurações cerâmicas à base de Y-TZP

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Effect Of Platform Connection And Abutment Material On Stress Distribution In Single Anterior Implant-supported Restorations: A Nonlinear 3-dimensional Finite Element Analysis.

Although various abutment connections and materials have recently been introduced, insufficient data exist regarding the effect of stress distribution on their mechanical performance. The purpose of this study was to investigate the effect of different abutment materials and platform connections on stress distribution in single anterior implant-supported restorations with the finite element method. Nine experimental groups were modeled from the combination of 3 platform connections (external hexagon, internal hexagon, and Morse tapered) and 3 abutment materials (titanium, zirconia, and hybrid) as follows: external hexagon-titanium, external hexagon-zirconia, external hexagon-hybrid, internal hexagon-titanium, internal hexagon-zirconia, internal hexagon-hybrid, Morse tapered-titanium, Morse tapered-zirconia, and Morse tapered-hybrid. Finite element models consisted of a 4×13-mm implant, anatomic abutment, and lithium disilicate central incisor crown cemented over the abutment. The 49 N occlusal loading was applied in 6 steps to simulate the incisal guidance. Equivalent von Mises stress (σvM) was used for both the qualitative and quantitative evaluation of the implant and abutment in all the groups and the maximum (σmax) and minimum (σmin) principal stresses for the numerical comparison of the zirconia parts. The highest abutment σvM occurred in the Morse-tapered groups and the lowest in the external hexagon-hybrid, internal hexagon-titanium, and internal hexagon-hybrid groups. The σmax and σmin values were lower in the hybrid groups than in the zirconia groups. The stress distribution concentrated in the abutment-implant interface in all the groups, regardless of the platform connection or abutment material. The platform connection influenced the stress on abutments more than the abutment material. The stress values for implants were similar among different platform connections, but greater stress concentrations were observed in internal connections.

Three-Dimensional Stress Distribution in the Human Periodontal Ligament in Masticatory, Parafunctional, and Trauma Loads: Finite Element Analysis

Background: The presence of the periodontal ligament (PDL) makes it possible to absorb and distribute loads produced during masticatory function and other tooth contacts into the alveolar process via the alveolar bone proper. However, several factors affect the integrity of periodontal structures causing the destruction of the connective matrix and cells, the loss of fibrous attachment, and the resorption of alveolar bone. Methods: The purpose of this study was to evaluate the stress distribution by finite element analysis in a PDL in three-dimensional models of the upper central incisor under three different load conditions: 100 N occlusal loading at 45 degrees (model 1: masticatory load); 500 N at the incisal edge at 45 degrees (model 2: parafunctional habit); and 800 N at the buccal surface at 90 degrees (model 3: trauma case). The models were built from computed tomography scans. Results: The stress distribution was quite different among the models. The most significant values (harmful) of tensile and compressive stresses were observed in models 2 and 3, with similarly distinct patterns of stress distributions along the PDL. Tensile stresses were observed along the internal and external aspects of the PDL, mostly at the cervical and middle thirds. Conclusions: The stress generation in these models may affect the integrity of periodontal structures. A better understanding of the biomechanical behavior of the PDL under physiologic and traumatic loading conditions might enhance the understanding of the biologic reaction of the PDL in health and disease. J Periodontol 2009;80:1859-1867.

Finite-Element Analysis of Stress on Dental Implant Prosthesis

Background: Understanding how clinical variables affect stress distribution facilitates optimal prosthesis design and fabrication and may lead to a decrease in mechanical failures as well as improve implant longevity. Purpose: In this study, the many clinical variations present in implant-supported prosthesis were analyzed by 3-D finite element method. Materials and Method: A geometrical model representing the anterior segment of a human mandible treated with 5 implants supporting a framework was created to perform the tests. The variables introduced in the computer model were cantilever length, elastic modulus of cancellous bone, abutment length, implant length, and framework alloy (AgPd or CoCr). The computer was programmed with physical properties of the materials as derived from the literature, and a 100N vertical load was used to simulate the occlusal force. Images with the fringes of stress were obtained and the maximum stress at each site was plotted in graphs for comparison. Results: Stresses clustered at the elements closest to the loading point. Stress increase was found to be proportional to the increase in cantilever length and inversely proportional to the increase in the elastic modulus of cancellous bone. Increasing the abutment length resulted in a decrease of stress on implants and framework. Stress decrease could not be demonstrated with implants longer than 13 mm. A stiffer framework may allow better stress distribution. Conclusion: The relative physical properties of the many materials involved in an implant-supported prosthesis system affect the way stresses are distributed.

Finite-Element Analysis of Stress on Dental Implant Prosthesis

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

Influence of Cusp Inclination on Stress Distribution in Implant-Supported Prostheses. A Three-Dimensional Finite Element Analysis

Purpose: The aim of this study was to assess the influence of cusp inclination on stress distribution in implant-supported prostheses by 3D finite element method.Materials and Methods: Three-dimensional models were created to simulate a mandibular bone section with an implant (3.75 mm diameter x 10 mm length) and crown by means of a 3D scanner and 3D CAD software. A screw-retained single crown was simulated using three cusp inclinations (10 degrees, 20 degrees, 30 degrees). The 3D models (model 10d, model 20d, and model 30d) were transferred to the finite element program NeiNastran 9.0 to generate a mesh and perform the stress analysis. An oblique load of 200 N was applied on the internal vestibular face of the metal ceramic crown.Results: The results were visualized by means of von Mises stress maps. Maximum stress concentration was located at the point of application. The implant showed higher stress values in model 30d (160.68 MPa). Cortical bone showed higher stress values in model 10d (28.23 MPa).Conclusion: Stresses on the implant and implant/abutment interface increased with increasing cusp inclination, and stresses on the cortical bone decreased with increasing cusp inclination.

Stress Analysis in Platform-Switching Implants: A 3-Dimensional Finite Element Study

The aim of this study was to evaluate the influence of the platform-switching technique on stress distribution in implant, abutment, and pen-implant tissues, through a 3-dimensional finite element study. Three 3-dimensional mandibular models were fabricated using the Solid Works 2006 and InVesalius software. Each model was composed of a bone block with one implant 10 mm long and of different diameters (3.75 and 5.00 mm). The UCLA abutments also ranged in diameter from 5.00 mm to 4.1 mm. After obtaining the geometries, the models were transferred to the software FEMAP 10.0 for pre- and postprocessing of finite elements to generate the mesh, loading, and boundary conditions. A total load of 200 N was applied in axial (0 degrees), oblique (45 degrees), and lateral (90) directions. The models were solved by the software NeiNastran 9.0 and transferred to the software FEMAP 10.0 to obtain the results that were visualized through von Mises and maximum principal stress maps. Model A (implants with 3.75 mm/abutment with 4.1 mm) exhibited the highest area of stress concentration with all loadings (axial, oblique, and lateral) for the implant and the abutment. All models presented the stress areas at the abutment level and at the implant/abutment interface. Models B (implant with 5.0 mm/abutment with 5.0 mm) and C (implant with 5.0 mm/abutment with 4.1 mm) presented minor areas of stress concentration and similar distribution pattern. For the cortical bone, low stress concentration was observed in the pen-implant region for models B and C in comparison to model A. The trabecular bone exhibited low stress that was well distributed in models B and C. Model A presented the highest stress concentration. Model B exhibited better stress distribution. There was no significant difference between the large-diameter implants (models B and C).