Application of domain decomposition and transmission line modelling techniques to 2D, time-domain, finite element problems

In this paper a recently published finite element method, which combines domain decomposition with a novel technique for solving nonlinear magnetostatic finite element problems is described. It is then shown how the method can be extended to, and optimised for, the solution of time-domain problems. © 1999 IEEE.

Finite element analysis of stress and strain distributions in InAs/GaAs quantum dots

In this paper, we perform systematic calculations of the stress and strain distributions in InAs/GaAs truncated pyramidal quantum dots (QDs) with different wetting layer (WL) thickness, using the finite element method (FEM). The stresses and strains are concentrated at the boundaries of the WL and QDs, are reduced gradually from the boundaries to the interior, and tend to a uniform state for the positions away from the boundaries. The maximal strain energy density occurs at the vicinity of the interface between the WL and the substrate. The stresses, strains and released strain energy are reduced gradually with increasing WL thickness. The above results show that a critical WL thickness may exist, and the stress and strain distributions can make the growth of QDs a growth of strained three-dimensional island when the WL thickness is above the critical value, and FEM can be applied to investigate such nanosystems, QDs, and the relevant results are supported by the experiments.

Finite-Element Time-Step Coupled Generator, Load, AVR and Brushless Exciter Modelling

This paper presented results from a details and comprehensive simulation using finite element method of the practical operation of an electrical machine. The results it displayed have been used in practice to design more efficient equipment.

A finite element solution scheme for an elastic–viscoplastic soil model

A Newton–Raphson solution scheme with a stress point algorithm is presented for the implementation of an elastic–viscoplastic soilmodel in a finite element program. Viscoplastic strain rates are calculated using the stress and volumetric states of the soil. Sub-incrementsof time are defined for each iterative calculation of elastic–viscoplastic stress changes so that their sum adds up to the time incrementfor the load step. This carefully defined ‘iterative time’ ensures that the correct amount of viscoplastic straining is accumulated overthe applied load step. The algorithms and assumptions required to implement the solution scheme are provided. Verification of the solutionscheme is achieved by using it to analyze typical boundary value problems.

Finite Element Application to Analysis, Testing, and Manufacturing of Satellite Structures

Designing satellite structures poses an ongoing challenge as the interaction between analysis, experimental testing, and manufacturing phases is underdeveloped. Finite Element Analysis for Satellite Structures: Applications to Their Design, Manufacture and Testing explains the theoretical and practical knowledge needed to perform design of satellite structures. By layering detailed practical discussions with fully developed examples, Finite Element Analysis for Satellite Structures: Applications to Their Design, Manufacture and Testing provides the missing link between theory and implementation.
Computational examples cover all the major aspects of advanced analysis; including modal analysis, harmonic analysis, mechanical and thermal fatigue analysis using finite element method. Test cases are included to support explanations an a range of different manufacturing simulation techniques are described from riveting to shot peening to material cutting. Mechanical design of a satellites structures are covered in three steps: analysis step under design loads, experimental testing to verify design, and manufacturing.
Stress engineers, lecturers, researchers and students will find Finite Element Analysis for Satellite Structures: Applications to Their Design, Manufacture and Testing a key guide on with practical instruction on applying manufacturing simulations to improve their design and reduce project cost, how to prepare static and dynamic test specifications, and how to use finite element method to investigate in more details any component that may fail during testing.

Reliability Assessment of 3-D Space Frame Structures Applying Stochastic Finite Element Analysis

Throughout design development of satellite structure, stress engineer is usually challenged with randomness in applied loads and material properties. To overcome such problem, a risk-based design is applied which estimates satellite structure probability of failure under static and thermal loads. Determining probability of failure can help to update initially applied factors of safety that were used during structure preliminary design phase. These factors of safety are related to the satellite mission objective. Sensitivity-based analysis is to be implemented in the context of finite element analysis (probabilistic finite element method or stochastic finite element method (SFEM)) to determine the probability of failure for satellite structure or one of its components.

Impact of finite element idealisation on the prediction of welded fuselage stiffened panel buckling

Lap joints are widely used in the manufacture of stiffened panels and influence local panel sub-component stability, defining buckling unit dimensions and boundary conditions. Using the Finite Element method it is possible to model joints in great detail and predict panel buckling behaviour with accuracy. However, when modelling large panel structures such detailed analysis becomes computationally expensive. Moreover, the impact of local behaviour on global panel performance may reduce as the scale of the modelled structure increases. Thus this study presents coupled computational and experimental analysis, aimed at developing relationships between modelling fidelity and the size of the modelled structure, when the global static load to cause initial buckling is the required analysis output. Small, medium and large specimens representing welded lap-joined fuselage panel structure are examined. Two element types, shell and solid-shell, are employed to model each specimen, highlighting the impact of idealisation on the prediction of welded stiffened panel initial skin buckling.

Finite element modeling of debonding failures in FRP-strengthened RC beams: A dynamic approach

This paper is concerned with the finite element simulation of debonding failures in FRP-strengthened concrete beams. A key challenge for such simulations is that common solution techniques such as the Newton-Raphson method and the arc-length method often fail to converge. This paper examines the effectiveness of using a dynamic analysis approach in such FE simulations, in which debonding failure is treated as a dynamic problem and solved using an appropriate time integration method. Numerical results are presented to show that an appropriate dynamic approach effectively overcomes the convergence problem and provides accurate predictions of test results.

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.

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)

Dynamic stationary response of reinforced plates by the boundary element method

A direct version of the boundary element method (BEM) is developed to model the stationary dynamic response of reinforced plate structures, such as reinforced panels in buildings, automobiles, and airplanes. The dynamic stationary fundamental solutions of thin plates and plane stress state are used to transform the governing partial differential equations into boundary integral equations (BIEs). Two sets of uncoupled BIEs are formulated, respectively, for the in-plane state ( membrane) and for the out-of-plane state ( bending). These uncoupled systems are joined to formamacro-element, in which membrane and bending effects are present. The association of these macro-elements is able to simulate thin-walled structures, including reinforced plate structures. In the present formulation, the BIE is discretized by continuous and/or discontinuous linear elements. Four displacement integral equations are written for every boundary node. Modal data, that is, natural frequencies and the corresponding mode shapes of reinforced plates, are obtained from information contained in the frequency response functions (FRFs). A specific example is presented to illustrate the versatility of the proposed methodology. Different configurations of the reinforcements are used to simulate simply supported and clamped boundary conditions for the plate structures. The procedure is validated by comparison with results determined by the finite element method (FEM).

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.

Finite element analysis and fracture resistance testing of a new intraradicular post

Objectives: The objective of the present study was to evaluate a prefabricated intraradicular threaded pure titanium post, designed and developed at the Sao Jose dos Campos School of Dentistry - UNESP, Brazil. This new post was designed to minimize stresses observed with prefabricated post systems and to improve cost-benefits. Materials and and methods: Fracture resistance testing of the post/core/root complex, fracture analysis by microscopy and stress analysis by the finite element method were used for post evaluation. The following four prefabricated metal post systems were analyzed: group 1, experimental post; group 2, modification of the experimental post; group 3, Flexi Post, and group 4, Para Post. For the analysis of fracture resistance, 40 bovine teeth were randomly assigned to the four groups (n=10) and used for the fabrication of test specimens simulating the situation in the mouth. The test specimens were subjected to compressive strength testing until fracture in an EMIC universal testing machine. After fracture of the test specimens, their roots were sectioned and analyzed by microscopy. For the finite element method, specimens of the fracture resistance test were simulated by computer modeling to determine the stress distribution pattern in the post systems studied. Results: The fracture test presented the following averages and standard deviation: G1 (45.63 +/- 8.77), G2 (49.98 +/- 7.08), G3 (43.84 +/- 5.52), G4 (47.61 +/- 7.23). Stress was homogenously distributed along the body of the intraradicular post in group 1, whereas high stress concentrations in certain regions were observed in the other groups. These stress concentrations in the body of the post induced the same stress concentration in root dentin. Conclusions: The experimental post (original and modified versions) presented similar fracture resistance and better results in the stress analysis when compared with the commercial post systems tested (08/2008PA/CEP).

Three-dimensional finite element analysis of MHD duct flow by the penalty function formulation

A numerical scheme based on the Finite Element Method (FEM) is presented to calculate the full solution of a three-dimensional steady magnetohydrodynamic (MHD) flow with moderately high Hartmann numbers and interaction parameters. An incompressible, viscous and electrically conducting liquid-metal is considered. Assuming a low magnetic Reynolds number, the solution method solves the coupled Navier-Stokes and Maxwell's equations through the use of a penalty function method. Results are presented for Hartmann numbers in the range 10(2)-10(3).