Finite Element Stress Analysis of Edentulous Mandibles with Different Bone Types Supporting Multiple-Implant Superstructures

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

Influence of Platform and Abutment Angulation on Peri-Implant Bone. A Three-Dimensional Finite Element Stress Analysis

The aim of this study was to evaluate stress distribution on the pen-implant bone, simulating the influence of Nobel Select implants with straight or angulated abutments on regular and switching platform in the anterior maxilla, by means of 3-dimensional finite element analysis. Four mathematical models of a central incisor supported by external hexagon implant (13 mm x 5 mm) were created varying the platform (R, regular or S. switching) and the abutments (S, straight or A, angulated 15 degrees). The models were created by using Mimics 13 and Solid Works 2010 software programs. The numerical analysis was performed using ANSYS Workbench 10.0. Oblique forces (100 N) were applied to the palatine surface of the central incisor. The bone/implant interface was considered perfectly integrated. Maximum (sigma(max)) and minimum (sigma(min)) principal stress values were obtained. For the cortical bone the highest stress values (sigma(max)) were observed in the RA (regular platform and angulated abutment, 51 MPa), followed by SA (platform switching and angulated abutment, 44.8 MPa), RS (regular platform and straight abutment, 38.6 MPa) and SS (platform switching and straight abutment, 36.5 MPa). For the trabecular bone, the highest stress values (sigma(max)) were observed in the RA (6.55 MPa), followed by RS (5.88 MPa), SA (5.60 MPa), and SS (4.82 MPa). The regular platform generated higher stress in the cervical periimplant region on the cortical and trabecular bone than the platform switching, irrespective of the abutment used (straight or angulated).

Non-linear refined finite element elastic analysis of steel frames with generalized transverse member loads

In the finite element modelling of steel frames, external loads usually act along the members rather than at the nodes only. Conventionally, when a member is subjected to these transverse loads, they are converted to nodal forces which act at the ends of the elements into which the member is discretised by either lumping or consistent nodal load approaches. For a contemporary geometrically non-linear analysis in which the axial force in the member is large, accurate solutions are achieved by discretising the member into many elements, which can produce unfavourable consequences on the efficacy of the method for analysing large steel frames. Herein, a numerical technique to include the transverse loading in the non-linear stiffness formulation for a single element is proposed, and which is able to predict the structural responses of steel frames involving the effects of first-order member loads as well as the second-order coupling effect between the transverse load and the axial force in the member. This allows for a minimal discretisation of a frame for second-order analysis. For those conventional analyses which do include transverse member loading, prescribed stiffness matrices must be used for the plethora of specific loading patterns encountered. This paper shows, however, that the principle of superposition can be applied to the equilibrium condition, so that the form of the stiffness matrix remains unchanged with only the magnitude of the loading being needed to be changed in the stiffness formulation. This novelty allows for a very useful generalised stiffness formulation for a single higher-order element with arbitrary transverse loading patterns to be formulated. The results are verified using analytical stability function studies, as well as with numerical results reported by independent researchers on several simple structural frames.

Bearing capacity of circular footings over rock mass by using axisymmetric quasi lower bound finite element limit analysis

The ultimate bearing capacity of a circular footing, placed over rock mass, is evaluated by using the lower bound theorem of the limit analysis in conjunction with finite elements and nonlinear optimization. The generalized Hoek-Brown (HB) failure criterion, but by keeping a constant value of the exponent, alpha = 0.5, was used. The failure criterion was smoothened both in the meridian and pi planes. The nonlinear optimization was carried out by employing an interior point method based on the logarithmic barrier function. The results for the obtained bearing capacity were presented in a non-dimensional form for different values of GSI, m(i), sigma(ci)/(gamma b) and q/sigma(ci). Failure patterns were also examined for a few cases. For validating the results, computations were also performed for a strip footing as well. The results obtained from the analysis compare well with the data reported in literature. Since the equilibrium conditions are precisely satisfied only at the centroids of the elements, not everywhere in the domain, the obtained lower bound solution will be approximate not true. (C) 2015 Elsevier Ltd. All rights reserved.

Reply to Discussion on ``Bearing capacity of circular footings over rock mass by using axisymmetric quasi lower bound finite element limit analysis'' Comput. Geotech. 70 (2015) 138-149]

A discussion has been provided for the comments raised by the discusser (Clausen, 2015)1] on the article recently published by the authors (Chakraborty and Kumar, 2015). The effect of exponent alpha for values of GSI approximately smaller than 30 becomes more critical. On the other hand, for greater values of GSI, the results obtained by the authors earlier remain primarily independent of alpha and can be easily used. (C) 2015 Elsevier Ltd. All rights reserved.

Slope stability analysis for filled slopes using finite element limit analysis method

This paper uses the finite element upper and lower bound limit analysis to assess the stability of slopes mostly found in embankment cases where frictional materials are filled on purely cohesive undrained clay. For comparison purposes, the commonly used stability assessment method, limit equilibrium method (LEM) is also employed. The final results for both methods are then presented in the form of comprehensive chart solutions for the convenience of practicing engineers during preliminary slope designs. The failure mechanism will also be discussed in this paper. Ultimately, it should be noted that finite element limit analysis method holds the upper hand as its prior assumptions are not required. Thus, the obtained failure mechanism from the slope stability analysis will be more realistic. Hence, it will provide a better understanding for the slope failure surface. Therefore, engineers should design more carefully when the LEM is applied to the slopes with frictional materials filled on purely cohesive undrained clay. © 2014 American Society of Civil Engineers.

Slope stability charts for two-layered purely cohesive soils based on finite-element limit analysis methods

Stability charts for soil slopes, first produced in the first half of the twentieth century, continue to be used extensively as design tools, and draw the attention of many investigators. This paper uses finite-element upper and lower bound limit analysis to assess the short-term stability of slopes in which the slopematerial and subgrade foundation material have two distinctly different undrained strengths. The stability charts are proposed, and the exact theoretical solutions are bracketed to within 4.2% or better. In addition, results from the limit-equilibrium method (LEM) have been used for comparison. Differences of up to 20% were found between the numerical limit analysis and LEM solutions. It also shown that the LEM sometimes leads to errors, although it is widely used in practice for slope stability assessments.

Slope stability analysis for fill slopes using finite element limit analysis

This paper investigates the stability of fill slopes often found in embankment cases where frictional fill materials are placed on purely cohesive undrained clay with increasing strength. By using finite element upper and lower bound limit analysis for this investigation, the limit load can be truly bounded. It is known that two-dimensional analysis yields a more conservative result due to plain strain condition when compared to three-dimensional analysis. Therefore, this paper will focus on three-dimensional (3D) slope stability analysis and for comparison purposes two-dimensional analysis results will be employed. In fact, the final results are presented in the form of comprehensive chart solutions for the convenience of practicing engineers during preliminary slope design. The failure mechanism will also be discussed in order to further illustrate the situation during failure. It should be highlighted that the failure mechanisms are obtained through the numerical method itself and no prior assumptions are required, therefore, are more realistic and able to provide a better understanding for the slope failure surfaces.

Identifying vulnerable carotid plaques in vivo using high resolution magnetic resonance imaging-based finite element analysis

Object. Individuals with carotid atherosclerosis develop symptoms following rupture of vulnerable plaques. Biomechanical stresses within this plaque may increase vulnerability to rupture. In this report the authors describe the use of in vivo carotid plaque imaging and computational mechanics to document the magnitude and distribution of intrinsic plaque stresses. Methods. Ten (five symptomatic and five asymptomatic) individuals underwent plaque characterization magnetic resonance (MR) imaging. Plaque geometry and composition were determined by multisequence review. Intrinsic plaque stress profiles were generated from 3D meshes by using finite element computational analysis. Differences in principal (shear) stress between normal and diseased sections of the carotid artery and between symptomatic and asymptomatic plaques were noted. Results. There was a significant difference in peak principal stress between diseased and nondiseased segments of the artery (mean difference 537.65 kPa, p < 0.05). Symptomatic plaques had higher mean stresses than asymptomatic plaques (627.6 kPa compared with 370.2 kPa, p = 0.05), which were independent of luminal stenosis and plaque composition. Conclusions. Significant differences in plaque stress exist between plaques from symptomatic individuals and those from asymptomatic individuals. The MR imaging-based computational analysis may therefore be a useful aid to identification of vulnerable plaques in vivo.

Comparison of theoretical fixation stability of three devices employed in medial opening wedge high tibial osteotomy: a finite element analysis.

BACKGROUND Medial open wedge high tibial osteotomy is a well-established procedure for the treatment of unicompartmental osteoarthritis and symptomatic varus malalignment. We hypothesized that different fixation devices generate different fixation stability profiles for the various wedge sizes in a finite element (FE) analysis. METHODS Four types of fixation were compared: 1) first and 2) second generation Puddu plates, and 3) TomoFix plate with and 4) without bone graft. Cortical and cancellous bone was modelled and five different opening wedge sizes were studied for each model. Outcome measures included: 1) stresses in bone, 2) relative displacement of the proximal and distal tibial fragments, 3) stresses in the plates, 4) stresses on the upper and lower screw surfaces in the screw channels. RESULTS The highest load for all fixation types occurred in the plate axis. For the vast majority of the wedge sizes and fixation types the shear stress (von Mises stress) was dominating in the bone independent of fixation type. The relative displacements of the tibial fragments were low (in μm range). With an increasing wedge size this displacement tended to increase for both Puddu plates and the TomoFix plate with bone graft. For the TomoFix plate without bone graft a rather opposite trend was observed.For all fixation types the occurring stresses at the screw-bone contact areas pulled at the screws and exceeded the allowable threshold of 1.2 MPa for at least one screw surface. Of the six screw surfaces that were studied, the TomoFix plate with bone graft showed a stress excess of one out of twelve and without bone graft, five out of twelve. With the Puddu plates, an excess stress occurred in the majority of screw surfaces. CONCLUSIONS The different fixation devices generate different fixation stability profiles for different opening wedge sizes. Based on the computational simulations, none of the studied osteosynthesis fixation types warranted an intransigent full weight bearing per se. The highest fixation stability was observed for the TomoFix plates and the lowest for the first generation Puddu plate. These findings were revealed in theoretical models and need to be validated in controlled clinical settings.

Comparison of beam theory and finite-element analysis with in vivo bone strain data from the alligator cranium

The mechanical behavior of the vertebrate skull is often modeled using free-body analysis of simple geometric structures and, more recently, finite-element (FE) analysis. In this study, we compare experimentally collected in vivo bone strain orientations and magnitudes from the cranium of the American alligator with those extrapolated from a beam model and extracted from an FE model. The strain magnitudes predicted from beam and FE skull models bear little similarity to relative and absolute strain magnitudes recorded during in vivo biting experiments. However, quantitative differences between principal strain orientations extracted from the FE skull model and recorded during the in vivo experiments were smaller, and both generally matched expectations from the beam model. The differences in strain magnitude between the data sets may be attributable to the level of resolution of the models, the material properties used in the FE model, and the loading conditions (i.e., external forces and constraints). This study indicates that FE models and modeling of skulls as simple engineering structures may give a preliminary idea of how these structures are loaded, but whenever possible, modeling results should be verified with either in vitro or preferably in vivo testing, especially if precise knowledge of strain magnitudes is desired. (c) 2005 Wiley-Liss, Inc.

Tooth fracture risk analysis based on a new finite element dental structure models using micro-CT data

The finite element (FE) analysis is an effective method to study the strength and predict the fracture risk of endodontically-treated teeth. This paper presents a rapid method developed to generate a comprehensive tooth FE model using data retrieved from micro-computed tomography (μCT). With this method, the inhomogeneity of material properties of teeth was included into the model without dividing the tooth model into different regions. The material properties of the tooth were assumed to be related to the mineral density. The fracture risk at different tooth portions was assessed for root canal treatments. The micro-CT images of a tooth were processed by a Matlab software programme and the CT numbers were retrieved. The tooth contours were obtained with thresholding segmentation using Amira. The inner and outer surfaces of the tooth were imported into Solidworks and a three-dimensional (3D) tooth model was constructed. An assembly of the tooth model with the periodontal ligament (PDL) layer and surrounding bone was imported into ABAQUS. The material properties of the tooth were calculated from the retrieved CT numbers via ABAQUS user's subroutines. Three root canal geometries (original and two enlargements) were investigated. The proposed method in this study can generate detailed 3D finite element models of a tooth with different root canal enlargements and filling materials, and would be very useful for the assessment of the fracture risk at different tooth portions after root canal treatments.