A finite element method for nonlinear elliptic problems

We present a Galerkin method with piecewise polynomial continuous elements for fully nonlinear elliptic equations. A key tool is the discretization proposed in Lakkis and Pryer, 2011, allowing us to work directly on the strong form of a linear PDE. An added benefit to making use of this discretization method is that a recovered (finite element) Hessian is a byproduct of the solution process. We build on the linear method and ultimately construct two different methodologies for the solution of second order fully nonlinear PDEs. Benchmark numerical results illustrate the convergence properties of the scheme for some test problems as well as the Monge–Amp`ere equation and the Pucci equation.

A finite element method for second order nonvariational elliptic problems

We propose a numerical method to approximate the solution of second order elliptic problems in nonvariational form. The method is of Galerkin type using conforming finite elements and applied directly to the nonvariational (nondivergence) form of a second order linear elliptic problem. The key tools are an appropriate concept of “finite element Hessian” and a Schur complement approach to solving the resulting linear algebra problem. The method is illustrated with computational experiments on three linear and one quasi-linear PDE, all in nonvariational form.

A finite element-spherical harmonics model for radiative transfer in inhomogeneous clouds. Part II. some applications

EVENT has been used to examine the effects of 3D cloud structure, distribution, and inhomogeneity on the scattering of visible solar radiation and the resulting 3D radiation field. Large eddy simulation and aircraft measurements are used to create realistic cloud fields which are continuous or broken with smooth or uneven tops. The values, patterns and variance in the resulting downwelling and upwelling radiation from incident visible solar radiation at different angles are then examined and compared to measurements. The results from EVENT confirm that 3D cloud structure is important in determining the visible radiation field, and that these results are strongly influenced by the solar zenith angle. The results match those from other models using visible solar radiation, and are supported by aircraft measurements of visible radiation, providing confidence in the new model.

Noether-type discrete conserved quantities arising from a finite element approximation of a variational problem

In this work, we prove a weak Noether-type Theorem for a class of variational problems that admit broken extremals. We use this result to prove discrete Noether-type conservation laws for a conforming finite element discretisation of a model elliptic problem. In addition, we study how well the finite element scheme satisfies the continuous conservation laws arising from the application of Noether’s first theorem (1918). We summarise extensive numerical tests, illustrating the conservation of the discrete Noether law using the p-Laplacian as an example and derive a geometric-based adaptive algorithm where an appropriate Noether quantity is the goal functional.

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.

Finite-Element Analysis of Stress on Dental Implant Prosthesis

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

NUMERICAL SIMULATION of CONVECTION-DIFFUSION PROBLEMS BY THE CONTROL-VOLUME-BASED FINITE-ELEMENT METHOD

This work presents a numerical study of the tri-dimensional convection-diffusion equation by the control-volume-based on finite-element method using quadratic hexahedral elements. Considering that the equation governing this problem in its main variable may represent several properties, including temperature, turbulent kinetic energy, viscous dissipation rate of the turbulent kinetic energy, specific dissipation rate of the turbulent kinetic energy, or even the concentration of a contaminant in a given medium, among others, the wide applicability of this problem is thus evidenced. Three cases of temperature distributions will be studied specifically in this work, in addition to one case of pollutant dispersion upon analysis of the concentration of a contaminant in a fixed flow point. Some comparisons will be carried out against works found in the open literature, while others will be done according to each phenomenon characteristics.

Finite Element Modeling for Development and Optimization of a Bone Plate for Mandibular Fracture in Dogs

This study aimed to develop a plate to treat fractures of the mandibular body in dogs and to validate the project using finite elements and biomechanical essays. Mandible prototypes were produced with 10 oblique ventrorostral fractures (favorable) and 10 oblique ventrocaudal fractures (unfavorable). Three groups were established for each fracture type. Osteosynthesis with a pure titanium plate of double-arch geometry and blocked monocortical screws offree angulanon were used. The mechanical resistance of the prototype with unfavorable fracture was lower than that of the fcworable fracture. In both fractures, the deflection increased and the relative stiffness decreased proportionally to the diminishing screw number The finite element analysis validated this plate study, since the maximum tension concentration observed on the plate was lower than the resistance limit tension admitted by the titanium. In conclusion, the double-arch geometry plate fixed with blocked monocortical screws has sufficient resistance to stabilize oblique,fractures, without compromising mandibular dental or neurovascular structures. J Vet Dent 24 (7); 212 - 221, 2010

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

Passivity Versus Unilateral Angular Misfit: Evaluation of Stress Distribution on Implant-Supported Single Crowns: Three-Dimensional Finite Element Analysis

The aim of this study was to evaluate the effect of unilateral angular misfit of 100 Km on stress distribution of implant-supported single crowns with ceramic veneering and gold framework by three-dimensional finite element analysis. Two three-dimensional models representing a maxillary section of premolar region were constructed: group 1 (control)-crown completely adapted to the implant and group 2-crown with unilateral angular misfit of 100 Km. A vertical force of 100 N was applied on 2 centric points of the crown. The von Mises stress was used as an analysis criterion. The stress values and distribution in the main maps (204.4 MPa for group 1 and 205.0 MPa for group 2) and in the other structures (aesthetic veneering, framework, retention screw, implant, and bone tissue) were similar for both groups. The highest stress values were observed between the first and second threads of the retention screw. Considering the bone tissue, the highest stress values were exhibited in the peri-implant cortical bone. The unilateral angular misfit of 100 Km did not influence the stress distribution on the implant-supported prosthesis under static loading.

Comparison of Single-Standing or Connected Implants on Stress Distribution in Bone of Mandibular Overdentures: A Two-Dimensional Finite Element Analysis

This study aimed to compare the influence of single-standing or connected implants on stress distribution in bone of mandibular overdentures by means of two-dimensional finite element analysis. Two finite element models were designed using software (ANSYS) for 2 situations: bar-clip (BC) group-model of an edentulous mandible supporting an overdenture over 2 connected implants with BC system, and o'ring (OR) group-model of an edentulous mandible supporting an overdenture over 2 single-standing implants with OR abutments. Axial loads (100 N) were applied on either central (L1) or lateral (L2) regions of the models. Stress distribution was concentrated mostly in the cortical bone surrounding the implants. When comparing the groups, BC (L1, 52.0 MPa and L2, 74.2 MPa) showed lower first principal stress values on supporting tissue than OR (L1, 78.4 MPa and L2, 76.7 MPa). Connected implants with BC attachment were more favorable on stress distribution over peri-implant-supporting tissue for both loading conditions.

Stress Analysis in Simulation Models With or Without Implant Threads Representation

Purpose: This study aimed to evaluate the influence of implants with or without threads representation on the outcome of a two-dimensional finite element (FE) analysis. Materials and Methods: Two-dimensional FE models that reproduced a frontal section of edentulous mandibular posterior bone were constructed using a standard crown/implant/screw system representation. To evaluate the effect of implant threads, two models were created: a model in which the implant threads were accurately simulated (precise model) and a model in which implants with a smooth surface (press-fit implant) were used (simplified model). An evaluation was performed on ANSYS software, in which a load of 133 N was applied at a 30-degree angulation and 2 mm off-axis from the long axis of the implant on the models, The Von Mises stresses were measured. Results: The precise model (1.45 MPa) showed higher maximum stress values than the simplified model (1.2 MPa). Whereas in the cortical bone, the stress values differed by about 36% (292.95 MPa for the precise model and 401.14 MPa for the simplified model), in trabecular bone (19.35 MPa and 20.35 MPa, respectively), the stress distribution and stress values were similar. Stress concentrations occurred around the implant neck and the implant apex. Conclusions: Considering implant and cortical bone analysis, remarkable differences in stress values were found between the models. Although the models showed different absolute stress values, the stress distribution was similar. INT J ORAL MAXILLOFAC IMPLANTS 2009;24:1040-1044

Effect of Passive Fit Absence in the Prosthesis/Implant/Retaining Screw System: A Two-Dimensional Finite Element Analysis

The misfit between prostheses and implants is a clinical reality, but the level that can be accepted without causing mechanical or biologic problem is not well defined. This study investigates the effect of different levels of unilateral angular misfit prostheses in the prosthesis/implant/retaining screw system and in the surrounding bone using finite element analysis. Four models of a two-dimensional finite element were constructed: group I (control), prosthesis that fit the implant; groups 2 to 4, prostheses with unilateral angular misfit of 50, 100, and 200 mu m, respectively. A load of 133 N was applied with a 30-degree angulation and off-axis at 2 mm from the long axis of the implant at the opposite direction of misfit on the models. Taking into account the increase of the angular misfit, the stress maps showed a gradual increase of prosthesis stress and uniform stress in the implant and trabecular bone. Concerning the displacement, an inclination of the system due to loading and misfit was observed. The decrease of the unilateral contact between prosthesis and implant leads to the displacement of the entire system, and distribution and magnitude alterations of the stress also occurred.

Finite Element Analysis to Compare Complete Denture and Implant-Retained Overdentures With Different Attachment Systems

This finite element analysis compared stress distribution on complete dentures and implant-retained overdentures with different attachment systems. Four models of edentulous mandible were constructed: group A (control), complete denture; group B, overdenture retained by 2 splinted implants with bar-clip system; group C, overdenture retained by 2 unsplinted implants with o'ring system; and group D, overdenture retained by 2 splinted implants with bar-clip and 2 distally placed o'ring system. Evaluation was performed on Ansys software, with 100-N vertical load applied on central incisive teeth. The lowest maximum general stress value (in megapascal) was observed in group A (64.305) followed by groups C (119.006), D (258.650), and B (349.873). The same trend occurred it) supporting tissues with the highest stress value for cortical bone. Unsplinted implants associated with the o'ring attachment system showed the lowest maximum stress values among all overdenture groups. Furthermore, o'ring system also improved stress distribution when associated with bar-clip system.