A finite element approach for predicting the ultimate rotation capacity of RC beams

This paper presents a numerical approach to model the complex failure mechanisms that define the ultimate rotational capacity of reinforced concrete beams. The behavior in tension and compression is described by a constitutive damage model derived from a combination of two specific damage models [1]. The nonlinear behavior of the compressed region is treated by the compressive damage model based on the Drucker-Prager criterion written in terms of the effective stresses. The tensile damage model employs a failure criterion based on the strain energy associated with the positive part the effective stress tensor. This model is used to describe the behavior of very thin bands of strain localization, which are embedded in finite elements to represent multiple cracks that occur in the tensioned region [2]. The softening law establishes dissipation energy compatible with the fracture energy of the concrete. The reinforcing steel bars are modeled by truss elements with elastic-perfect plastic behavior. It is shown that the resulting approach is able to predict the different stages of the collapse mechanism of beams with distinct sizes and reinforcement ratios. The tensile damage model and the finite element embedded crack approach are able to describe the stiffness reduction due to concrete cracking in the tensile zone. The truss elements are able to reproduce the effects of steel yielding and, finally, the compressive damage model is able to describe the non-linear behavior of the compressive zone until the complete collapse of the beam due to crushing of concrete. The proposed approach is able to predict well the plastic rotation capacity of tested beams [3], including size-scale effects.

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.

Finite element analysis of concrete cracking at early age

The study of the early age concrete properties is becoming more important, as the thermal effects and the shrinkage, even in the first hours, could generate cracks, increasing the permeability of the structure and being able to induce problems of durability and functionality in the same ones. The detailed study of the stresses development during the construction process can be decisive to keep low the cracking levels. In this work a computational model, based on the finite element method, was implemented to simulate the early age concrete behavior and, specially, the evaluation of the cracking risk. The finite element analysis encloses the computational modeling of the following phenomena: chemical, thermal, moisture diffusion and mechanical which occur at the first days after the concrete cast. The developed software results were compared with experimental values found in the literature, demonstrating an excellent approach for all the implemented analysis.

Finite element and dimensional analysis algorithm for the prediction of mechanical properties of bulk materials and thin films

In this work, the applicability of a new algorithm for the estimation of mechanical properties from instrumented indentation data was studied for thin films. The applicability was analyzed with the aid of both three-dimensional finite element simulations and experimental indentation tests. The numerical approach allowed studying the effect of the substrate on the estimation of mechanical properties of the film, which was conducted based on the ratio h(max)/l between maximum indentation depth and film thickness. For the experimental analysis, indentation tests were conducted on AISI H13 tool steel specimens, plasma nitrated and coated with TiN thin films. Results have indicated that, for the conditions analyzed in this work, the elastic deformation of the substrate limited the extraction of mechanical properties of the film/substrate system. This limitation occurred even at low h(max)/l ratios and especially for the estimation of the values of yield strength and strain hardening exponent. At indentation depths lower than 4% of the film thickness, the proposed algorithm estimated the mechanical properties of the film with accuracy. Particularly for hardness, precise values were estimated at h(max)/l lower than 0.1, i.e. 10% of film thickness. (C) 2010 Published by Elsevier B.V.

Finite element modelling of three-dimensional convection problems in fluid-saturated porous media heated from below

We use the finite element method to model three-dimensional convective pore-fluid flow in fluid-saturated porous media when they are heated from below. In particular, we employ the particle-tracking technique to mimic the trajectories of particles in three-dimensional fluid flow problems. The related numerical results demonstrated that: (1) The progressive asymptotic approach procedure, which was previously developed for the finite element modelling of two-dimensional convective pore-fluid flow problems, is equally applicable to the finite element modelling of three-dimensional convective pore-fluid flow in fluid-saturated porous media heated from below. (2) The perturbation of gravity at different planes has a significant effect on the pattern of three-dimensional convective pore-fluid flow and therefore, may influence the pattern of orebody formation and mineralization in three-dimensional hydrothermal systems. Copyright (C) 2001 John Wiley & Sons, Ltd.

Finite element modelling of rock alteration and metamorphic process in hydrothermal systems

We use the finite element method to simulate the rock alteration and metamorphic process in hydrothermal systems. In particular, we consider the fluid-rock interaction problems in pore-fluid saturated porous rocks. Since the fluid rock interaction takes place at the contact interface between the pore-fluid and solid minerals, it is governed by the chemical reaction which usually takes place very slowly at this contact interface, from the geochemical point of view. Due to the relative slowness of the rate of the chemical reaction to the velocity of the pore-fluid flow in the hydrothermal system to be considered, there exists a retardation zone, in which the conventional static theory in geochemistry does not hold true. Since this issue is often overlooked by some purely numerical modellers, it is emphasized in this paper. The related results from a typical rock alteration and metamorphic problem in a hydrothermal system have shown not only the detailed rock alteration and metamorphic process, but also the size of the retardation zone in the hydrothermal system. Copyright (C) 2001 John Wiley & Sons, Ltd.

Estimating Local Part Thickness in Midplane Meshes for Finite Element Analysis

Within the development of motor vehicles, crash safety (e.g. occupant protection, pedestrian protection, low speed damageability), is one of the most important attributes. In order to be able to fulfill the increased requirements in the framework of shorter cycle times and rising pressure to reduce costs, car manufacturers keep intensifying the use of virtual development tools such as those in the domain of Computer Aided Engineering (CAE). For crash simulations, the explicit finite element method (FEM) is applied. The accuracy of the simulation process is highly dependent on the accuracy of the simulation model, including the midplane mesh. One of the roughest approximations typically made is the actual part thickness which, in reality, can vary locally. However, almost always a constant thickness value is defined throughout the entire part due to complexity reasons. On the other hand, for precise fracture analysis within FEM, the correct thickness consideration is one key enabler. Thus, availability of per element thickness information, which does not exist explicitly in the FEM model, can significantly contribute to an improved crash simulation quality, especially regarding fracture prediction. Even though the thickness is not explicitly available from the FEM model, it can be inferred from the original CAD geometric model through geometric calculations. This paper proposes and compares two thickness estimation algorithms based on ray tracing and nearest neighbour 3D range searches. A systematic quantitative analysis of the accuracy of both algorithms is presented, as well as a thorough identification of particular geometric arrangements under which their accuracy can be compared. These results enable the identification of each technique’s weaknesses and hint towards a new, integrated, approach to the problem that linearly combines the estimates produced by each algorithm.

Strength prediction of single- and double-lap joints by standard and extended finite element modelling

The structural integrity of multi-component structures is usually determined by the strength and durability of their unions. Adhesive bonding is often chosen over welding, riveting and bolting, due to the reduction of stress concentrations, reduced weight penalty and easy manufacturing, amongst other issues. In the past decades, the Finite Element Method (FEM) has been used for the simulation and strength prediction of bonded structures, by strength of materials or fracture mechanics-based criteria. Cohesive-zone models (CZMs) have already proved to be an effective tool in modelling damage growth, surpassing a few limitations of the aforementioned techniques. Despite this fact, they still suffer from the restriction of damage growth only at predefined growth paths. The eXtended Finite Element Method (XFEM) is a recent improvement of the FEM, developed to allow the growth of discontinuities within bulk solids along an arbitrary path, by enriching degrees of freedom with special displacement functions, thus overcoming the main restriction of CZMs. These two techniques were tested to simulate adhesively bonded single- and double-lap joints. The comparative evaluation of the two methods showed their capabilities and/or limitations for this specific purpose.

The use of the boundary element method in the analysis of single lap joints

The most common techniques for stress analysis/strength prediction of adhesive joints involve analytical or numerical methods such as the Finite Element Method (FEM). However, the Boundary Element Method (BEM) is an alternative numerical technique that has been successfully applied for the solution of a wide variety of engineering problems. This work evaluates the applicability of the boundary elem ent code BEASY as a design tool to analyze adhesive joints. The linearity of peak shear and peel stresses with the applied displacement is studied and compared between BEASY and the analytical model of Frostig et al., considering a bonded single-lap joint under tensile loading. The BEM results are also compared with FEM in terms of stress distributions. To evaluate the mesh convergence of BEASY, the influence of the mesh refinement on peak shear and peel stress distributions is assessed. Joint stress predictions are carried out numerically in BEASY and ABAQUS®, and analytically by the models of Volkersen, Goland, and Reissner and Frostig et al. The failure loads for each model are compared with experimental results. The preparation, processing, and mesh creation times are compared for all models. BEASY results presented a good agreement with the conventional methods.

Reappraisal of the effective elastic thickness for the sub-Andes using 3-D finite element flexural modelling, gravity and geological constraints

P>Estimates of effective elastic thickness (T(e)) for the western portion of the South American Plate using, independently, forward flexural modelling and coherence analysis, suggest different thermomechanical properties for the same continental lithosphere. We present a review of these T(e) estimates and carry out a critical reappraisal using a common methodology of 3-D finite element method to solve a differential equation for the bending of a thin elastic plate. The finite element flexural model incorporates lateral variations of T(e) and the Andes topography as the load. Three T(e) maps for the entire Andes were analysed: Stewart & Watts (1997), Tassara et al. (2007) and Perez-Gussinye et al. (2007). The predicted flexural deformation obtained for each T(e) map was compared with the depth to the base of the foreland basin sequence. Likewise, the gravity effect of flexurally induced crust-mantle deformation was compared with the observed Bouguer gravity. T(e) estimates using forward flexural modelling by Stewart & Watts (1997) better predict the geological and gravity data for most of the Andean system, particularly in the Central Andes, where T(e) ranges from greater than 70 km in the sub-Andes to less than 15 km under the Andes Cordillera. The misfit between the calculated and observed foreland basin subsidence and the gravity anomaly for the Maranon basin in Peru and the Bermejo basin in Argentina, regardless of the assumed T(e) map, may be due to a dynamic topography component associated with the shallow subduction of the Nazca Plate beneath the Andes at these latitudes.

Comparison of different designs of implant-retained overdentures and fixed full-arch implant-supported prosthesis on stress distribution in edentulous mandible - A computed tomography-based three-dimensional finite element analysis

A finite element analysis was used to compare the effect of different designs of implant-retained overdentures and fixed full-arch implant-supported prosthesis on stress distribution in edentulous mandible. Four models of an human mandible were constructed. In the OR (O'ring) group, the mandible was restored with an overdenture retained by four unsplinted implants with O'ring attachment; in the BC (bar-clip) -C and BC groups, the mandibles were restored with overdentures retained by four splinted implants with bar-clip anchor associated or not with two distally placed cantilevers, respectively; in the FD (fixed denture) group, the mandible was restored with a fixed full-arch four-implant-supported prosthesis. Models were supported by the masticatory muscles and temporomandibular joints. A 100-N oblique load was applied on the left first molar. Von Mises (σvM), maximum (σmax) and minimum (σmin) principal stresses (in MPa) analyses were obtained. BC-C group exhibited the highest stress values (σvM=398.8, σmax=580.5 and σmin=-455.2) while FD group showed the lowest one (σvM=128.9, σmax=185.9 and σmin=-172.1). Within overdenture groups, the use of unsplinted implants reduced the stress level in the implant/prosthetic components (59.4% for σvM, 66.2% for σmax and 57.7% for σmin versus BC-C group) and supporting tissues (maximum stress reduction of 72% and 79.5% for σmax, and 15.7% and 85.7% for σmin on the cortical and trabecular bones, respectively). Cortical bone exhibited greater stress concentration than the trabecular bone for all groups. The use of fixed implant dentures and removable dentures retained by unsplinted implants to rehabilitate edentulous mandible reduced the stresses in the periimplant bone tissue, mucosa and implant/prosthetic components. © 2013 Elsevier Ltd.

Free tools and strategies for the generation of 3D finite element meshes: Modeling of the cardiac structures

The Finite Element Method is a well-known technique, being extensively applied in different areas. Studies using the Finite Element Method (FEM) are targeted to improve cardiac ablation procedures. For such simulations, the finite element meshes should consider the size and histological features of the target structures. However, it is possible to verify that some methods or tools used to generate meshes of human body structures are still limited, due to nondetailed models, nontrivial preprocessing, or mainly limitation in the use condition. In this paper, alternatives are demonstrated to solid modeling and automatic generation of highly refined tetrahedral meshes, with quality compatible with other studies focused on mesh generation. The innovations presented here are strategies to integrate Open Source Software (OSS). The chosen techniques and strategies are presented and discussed, considering cardiac structures as a first application context. © 2013 E. Pavarino et al.

Bone stress and strain after use of a miniplate for molar protraction and uprighting: a 3-dimensional finite element analysis

The aim of this study was to use the finite element method to evaluate the distribution of stresses and strains on the local bone tissue adjacent to the miniplate used for anchorage of orthodontic forces. Methods: A 3-dimensional model composed of a hemimandible and teeth was constructed using dental computed tomographic images, in which we assembled a miniplate with fixation screws. The uprighting and mesial movements of the mandibular second molar that was anchored with the miniplate were simulated. The miniplate was loaded with horizontal forces of 2, 5, and 15 N. A moment of 11.77 N.mm was also applied. The stress and strain distributions were analyzed, and their correlations with the bone remodeling criteria and miniplate stability were assessed. Results: When orthodontic loads were applied, peak bone strain remained within the range of bone homeostasis (100-1500 mu m strain) with a balance between bone formation and resorption. The maximum deformation was found to be 1035 mu m strain with a force of 5 N. At a force of 15 N, bone resorption was observed in the region of the screws. Conclusions: We observed more stress concentration around the screws than in the cancellous bone. The levels of stress and strain increased when the force was increased but remained within physiologic levels. The anchorage system of miniplate and screws could withstand the orthodontic forces, which did not affect the stability of the miniplate.

Finite element simulation of a local scale air quality model over complex terrain

[EN]In this paper we propose a finite element method approach for modelling the air quality in a local scale over complex terrain. The area of interest is up to tens of kilometres and it includes pollutant sources. The proposed methodology involves the generation of an adaptive tetrahedral mesh, the computation of an ambient wind field, the inclusion of the plume rise effect in the wind field, and the simulation of transport and reaction of pollutants. We apply our methodology to simulate a fictitious pollution episode in La Palma island (Canary Island, Spain)...