Finite element simulations of ductile rupture in a constrained metal foil

A numerical study of the ductile rupture in a metal foil constrained between two stiff ceramic blocks is performed. The finite element analysis is carried out under the conditions of mode I, plane strain, small-scale yielding. The rate-independent version of the Gurson model that accounts for the ductile failure mechanisms of microvoid nucleation, growth and coalescence is employed to represent the behavior of the metal foil. Different distributions of void nucleating sites in the metal foil are considered for triggering the initiation of discrete voids. The results clearly show that far-field triaxiality-induced cavitation is the dominant failure mode when the spacing of the void nucleating sites is large. On the contrary, void coalescence near the notch tip is found to be the operative failure mechanism when closely spaced void nucleating sites are considered.

A finite element assessment of flexural strength of prestressed concrete beams with fiber reinforcement

This paper presents an assessment of the flexural behavior of 15 fully/partially prestressed high strength concrete beams containing steel fibers investigated using three-dimensional nonlinear finite elemental analysis. The experimental results consisted of eight fully and seven partially prestressed beams, which were designed to be flexure dominant in the absence of fibers. The main parameters varied in the tests were: the levels of prestressing force (i.e, in partially prestressed beams 50% of the prestress was reduced with the introduction of two high strength deformed bars instead), fiber volume fractions (0%, 0.5%, 1.0% and 1.5%), fiber location (full depth and partial depth over full length and half the depth over the shear span only). A three-dimensional nonlinear finite element analysis was conducted using ANSYS 5.5 [Theory Reference Manual. In: Kohnke P, editor. Elements Reference Manual. 8th ed. September 1998] general purpose finite element software to study the flexural behavior of both fully and partially prestressed fiber reinforced concrete beams. Influence of fibers on the concrete failure surface and stress-strain response of high strength concrete and the nonlinear stress-strain curves of prestressing wire and deformed bar were considered in the present analysis. In the finite element model. tension stiffening and bond slip between concrete and reinforcement (fibers., prestressing wire, and conventional reinforcing steel bar) have also been considered explicitly. The fraction of the entire volume of the fiber present along the longitudinal axis of the prestressed beams alone has been modeled explicitly as it is expected that these fibers would contribute to the mobilization of forces required to sustain the applied loads across the crack interfaces through their bridging action. A comparison of results from both tests and analysis on all 15 specimens confirm that, inclusion of fibers over a partial depth in the tensile side of the prestressed flexural structural members was economical and led to considerable cost saving without sacrificing on the desired performance. However. beams having fibers over half the depth in only the shear span, did not show any increase in the ultimate load or deformational characteristics when compared to plain concrete beams. (C) 2002 Published by Elsevier Science Ltd.

Reconstructing Solid Model from 2D Scanned Images of Biological Organs for Finite Element Simulation

This work presents a methodology to reconstruct 3D biological organs from image sequences or other scan data using readily available free softwares with the final goal of using the organs (3D solids) for finite element analysis. The methodology deals with issues such as segmentation, conversion to polygonal surface meshes, and finally conversion of these meshes to 3D solids. The user is able to control the detail or the level of complexity of the solid constructed. The methodology is illustrated using 3D reconstruction of a porcine liver as an example. Finally, the reconstructed liver is imported into the commercial software ANSYS, and together with a cyst inside the liver, a nonlinear analysis performed. The results confirm that the methodology can be used for obtaining 3D geometry of biological organs. The results also demonstrate that the geometry obtained by following this methodology can be used for the nonlinear finite element analysis of organs. The methodology (or the procedure) would be of use in surgery planning and surgery simulation since both of these extensively use finite elements for numerical simulations and it is better if these simulations are carried out on patient specific organ geometries. Instead of following the present methodology, it would cost a lot to buy a commercial software which can reconstruct 3D biological organs from scanned image sequences.

Uniformly convergent finite element methods for singularly perturbed parabolic partial differential equations

This thesis is concerned with uniformly convergent finite element methods for numerically solving singularly perturbed parabolic partial differential equations in one space variable. First, we use Petrov-Galerkin finite element methods to generate three schemes for such problems, each of these schemes uses exponentially fitted elements in space. Two of them are lumped and the other is non-lumped. On meshes which are either arbitrary or slightly restricted, we derive global energy norm and L2 norm error bounds, uniformly in the diffusion parameter. Under some reasonable global assumptions together with realistic local assumptions on the solution and its derivatives, we prove that these exponentially fitted schemes are locally uniformly convergent, with order one, in a discrete L∞norm both outside and inside the boundary layer. We next analyse a streamline diffusion scheme on a Shishkin mesh for a model singularly perturbed parabolic partial differential equation. The method with piecewise linear space-time elements is shown, under reasonable assumptions on the solution, to be convergent, independently of the diffusion parameter, with a pointwise accuracy of almost order 5/4 outside layers and almost order 3/4 inside the boundary layer. Numerical results for the above schemes are presented. Finally, we examine a cell vertex finite volume method which is applied to a model time-dependent convection-diffusion problem. Local errors away from all layers are obtained in the l2 seminorm by using techniques from finite element analysis.

Evaluation of distortions in laser welded shipbuilding parts using local-global finite element approach

In recognition of the differences of scale between the welding pool and the heat affected zone along the welding line on one hand, and the overall size of the components being welded on the other, a local-global finite element approach was developed for the evaluation of distortions in laser welded shipbuilding parts. The approach involves the tandem use of a 'local' and a 'global' step. The local step involves a three-dimensional finite element model for the simulation of the laser welding process using the Sysweld finite element code, which takes into account thermal, metallurgical, and mechanical aspects. The simulation of the laser welding process was performed using a non-linear heat transfer analysis, based on a keyhole formation model, and a coupled transient thermomechanical analysis, which takes into account metallurgical transformations using the temperature dependent material properties and the continuous cooling transformation diagram. The size and shape of the keyhole used in the local finite element analysis was evaluated using a keyhole formation model and the Physica finite volume code. The global step involves the transfer of residual plastic strains and the stiffness of the weld obtained from the local model to the global analysis, which then provides the predicted distortions for the whole part. This newly developed methodology was applied to the evaluation of global distortions due to laser welding of stiffeners on a shipbuilding part. The approach has been proved reliable in comparison with experiments and of practical industrial use in terms of computing time and storage.

Finite element modelling of crack detection tests

Four non-destructive tests for determining the length of fatigue cracks within the solder joints of a 2512 surface mount resistor are investigated. The sensitivity of the tests is obtained using finite element analysis with some experimental validation. Three of the tests are mechanically based and one is thermally based. The mechanical tests all operate by applying different loads to the PCB and monitoring the strain response at the top of the resistor. The thermal test operates by applying a heat source underneath the PCB, and monitoring the temperature response at the top of the resistor. From the modelling work done, two of these tests have shown to be sensitive to cracks. Some experimental results are presented but further work is required to fully validate the simulation results.

Mixed-dimensional coupling in finite element models

In many finite element analysis models it would be desirable to combine reduced- or lower-dimensional element types with higher-dimensional element types in a single model. In order to achieve compatibility of displacements and stress equilibrium at the junction or interface between the differing element types, it is important in such cases to integrate into the analysis some scheme for coupling the element types. A novel and effective scheme for establishing compatibility and equilibrium at the dimensional interface is described and its merits and capabilities are demonstrated. Copyright (C) 2000 John Wiley & Sons, Ltd.

Finite element investigation of the structural behaviour of deck slabs in composite bridges

This research studies the structural behaviour of bridge deck slabs under static patch loads in steel–concrete composite bridges and investigates compressive membrane action (CMA) in concrete bridge decks slabs, which governs the structural behaviour. A non-linear 3D finite element analysis models was developed using ABAQUS 6.5 software packages. Experimental data from one-span composite bridge structures are used to validate and calibrate the proposed FEM models. A series of parametric studies is conducted. The analysis results are discussed and conclusions on the behaviour of the bridge decks are presented.

Predicting revision risk for aseptic loosening of femoral components in total flip arthroplasty in individual patients - A finite element study

Advances in surgical procedure, prosthesis design, and biomaterials performance have considerably increased the longevity of total joint replacements. Preoperative planning is another step in joint replacement that may have the potential to improve clinical outcome for the individual patient, but has remained relatively consistent for a longtime. One means of advancing this aspect of joint replacement surgery may be to include predictive computer simulation into the planning process. In this article, the potential of patient-specific finite element analysis in preoperative assessment is investigated. Seventeen patient-specific finite element models of cemented Charnley reconstructions were created, of which six were early (

Evaluation of cement stresses in finite element analyses of cemented orthopaedic implants

Stress analysis of the cement fixation of orthopaedic implants to bone is frequently? carried out using finite element analysis. However, the stress distribution in the cement laver is usually intricate, and it is difficult to report it in a way that facilitates comparison of implants for pre-clinical testing. To study this problem, and make recommendations for stress reporting, a finite element analysis of a hip prosthesis implanted into a synthetic composite femur is developed. Three cases are analyzed: a fully bonded implant, a debonded implant, and a debonded implant where the cement is removed distal to the stein tip. In addition to peak stresses, and contour and vector plots, a stressed volume and probability-of-failure analysis is reported. It is predicted that the peak stress is highest for the debonded stem, and that removal of the distal cement more than halves this peak stress. This would suggest that omission of the distal cement is good for polished prostheses (as practiced for the Exeter design). However; if the percentage of cement stressed above a certain threshold (say 3 MPa) is considered, then the removal of distal cement is shown to be disadvantageous because a higher volume of cement is stressed to above the threshold. Vector plots clearly demonstrate the different load transfer for bonded and debonded prostheses: A bonded stein generates maximum tensile stresses in the longitudinal direction, whereas a debonded stem generates most tensile stresses in the hoop direction, except near the tip where tensile longitudinal stresses occur due to subsidence of the stein. Removal of the cement distal to the tip allows greater subsidence but alleviates these large stresses at the tip, albeit at the expense of increased hoop stresses throughout the mantle. It is concluded that a thorough analysis of cemented implants should not report peak stress, which can be misleading, but rather stressed volume, and that vector plots should be reported if a precise analysis of the load transfer mechanism is required.

Finite element modelling of composite structures under crushing load

This paper details the theory and implementation of a composite damage model, addressing damage within a ply (intralaminar) and delamination (interlaminar), for the simulation of crushing of laminated composite structures. It includes a more accurate determination of the characteristic length to achieve mesh objectivity in capturing intralaminar damage consisting of matrix cracking and fibre failure, a load-history dependent material response, an isotropic hardening nonlinear matrix response, as well as a more physically-based interactive matrix-dominated damage mechanism. The developed damage model requires a set of material parameters obtained from a combination of standard and non-standard material characterisation tests. The fidelity of the model mitigates the need to manipulate, or "calibrate", the input data to achieve good agreement with experimental results. The intralaminar damage model was implemented as a VUMAT subroutine, and used in conjunction with an existing interlaminar damage model, in Abaqus/Explicit. This approach was validated through the simulation of the crushing of a cross-ply composite tube with a tulip-shaped trigger, loaded in uniaxial compression. Despite the complexity of the chosen geometry, excellent correlation was achieved with experimental 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.

The design and development of an integrated system for object-oriented finite element computing

Finite Element Method (FEM) is widely used in Science and Engineering since 1960’s. The vast majority of FEM software is procedure-oriented. However, this conventional style of designing FEM software encounters problems in maintenance, reuse, and expansion of the software. Recently the object-oriented finite element method attracts the attention of lots of researchers, and now there is a growing interest in this method. In this paper, the object-oriented finite element (OOFE) is briefly introduced. Then the design and development of an integrated OOFE system is described. A comparison of the integrated OOFE system and a procedure-oriented system shows that our OOFE system has many advantages.

Finite element simulation of deformation processes of micro-components

The aim of this paper is to improve the understanding of deformation of micro medical needle and thread during assembly and then to develop an economical and flexible deformation method. Therefore, the swaging process is computationally simulated with the finite element method in this paper. A commercially available explicit nonlinear finite element analysis code, LS-Dyna, is used to model the 3-D deformation and contact problem. As the firmness of the assembly on the needle depends on the contact force and friction, the contact and the slide between the needle and thread are taken into account in the simulation. The general surface-to-surface contact algorithm (STS) is used to simulate the contact. The paper provides an insight into the deformation of the micro products.