Showing papers on "Extended finite element method published in 2003"

Finite Element Methods for Maxwell's Equations

Abstract: We survey finite element methods for approximating the time harmonic Maxwell equations. We concentrate on comparing error estimates for problems with spatially varying coefficients. For the conforming edge finite element methods, such estimates allow, at least, piecewise smooth coefficients. But for Discontinuous Galerkin (DG) methods, the state of the art of error analysis is less advanced (we consider three DG families of methods: Interior Penalty type, Hybridizable DG, and Trefftz type methods). Nevertheless, DG methods offer significant potential advantages compared to conforming methods.

New crack-tip elements for XFEM and applications to cohesive cracks

Abstract: An extended finite element method scheme for a static cohesive crack is developed with a new formulation for elements containing crack tips. This method can treat arbitrary cracks independent of the mesh and crack growth without remeshing. All cracked elements are enriched by the sign function so that no blending of the local partition of unity is required. This method is able to treat the entire crack with only one type of enrichment function, including the elements containing the crack tip. This scheme is applied to linear 3-node triangular elements and quadratic 6-node triangular elements. To ensure smooth crack closing of the cohesive crack, the stress projection normal to the crack tip is imposed to be equal to the material strength. The equilibrium equation and the traction condition are solved by the Newton–Raphson method to obtain the nodal displacements and the external load simultaneously. The results obtained by the new extended finite element method are compared to reference solutions and show excellent agreement. Copyright © 2003 John Wiley & Sons, Ltd.

Dynamic crack propagation based on loss of hyperbolicity and a new discontinuous enrichment

Abstract: SUMMARY A methodology is developed for switching from a continuum to a discrete discontinuity where the governing partial dierential equation loses hyperbolicity. The approach is limited to rate-independent materials, so that the transition occurs on a set of measure zero. The discrete discontinuity is treated by the extendednite element method (XFEM) whereby arbitrary discontinuities can be incorporated in the model without remeshing. Loss of hyperbolicity is tracked by a hyperbolicity indicator that enables both the crack speed and crack direction to be determined for a given material model. A new method was developed for the case when the discontinuity ends within an element; it facilitates the modelling of crack tips that occur within an element in a dynamic setting. The method is applied to several dynamic crack growth problems including the branching of cracks. Copyright ? 2003 John Wiley & Sons, Ltd.

A mixed multiscale finite element method for elliptic problems with oscillating coefficients

Abstract: The recently introduced multiscale finite element method for solving elliptic equations with oscillating coefficients is designed to capture the large-scale structure of the solutions without resolving all the fine-scale structures. Motivated by the numerical simulation of flow transport in highly heterogeneous porous media, we propose a mixed multiscale finite element method with an over-sampling technique for solving second order elliptic equations with rapidly oscillating coefficients. The multiscale finite element bases are constructed by locally solving Neumann boundary value problems. We provide a detailed convergence analysis of the method under the assumption that the oscillating coefficients are locally periodic. While such a simplifying assumption is not required by our method, it allows us to use homogenization theory to obtain the asymptotic structure of the solutions. Numerical experiments are carried out for flow transport in a porous medium with a random log-normal relative permeability to demonstrate the efficiency and accuracy of the proposed method.

Spectral Element Method in Structural Dynamics

Abstract: Preface. Part One Introduction to the Spectral Element Method and Spectral Analysis of Signals. 1 Introduction. 1.1 Theoretical Background. 1.2 Historical Background. 2 Spectral Analysis of Signals. 2.1 Fourier Series. 2.2 Discrete Fourier Transform and the FFT. 2.3 Aliasing. 2.4 Leakage. 2.5 Picket-Fence Effect. 2.6 Zero Padding. 2.7 Gibbs Phenomenon. 2.8 General Procedure of DFT Processing. 2.9 DFTs of Typical Functions. Part Two Theory of Spectral Element Method. 3 Methods of Spectral Element Formulation. 3.1 Force-Displacement Relation Method. 3.2 Variational Method. 3.3 State-Vector Equation Method. 3.4 Reduction from the Finite Models. 4 Spectral Element Analysis Method. 4.1 Formulation of Spectral Element Equation. 4.2 Assembly and the Imposition of Boundary Conditions. 4.3 Eigenvalue Problem and Eigensolutions. 4.4 Dynamic Responses with Null Initial Conditions. 4.5 Dynamic Responses with Arbitrary Initial Conditions. 4.6 Dynamic Responses of Nonlinear Systems. Part Three Applications of Spectral Element Method. 5 Dynamics of Beams and Plates. 5.1 Beams. 5.2 Levy-Type Plates. 6 Flow-Induced Vibrations of Pipelines. 6.1 Theory of Pipe Dynamics. 6.2 Pipelines Conveying Internal Steady Fluid. 6.3 Pipelines Conveying Internal Unsteady Fluid. Appendix 6.A: Finite Element Matrices: Steady Fluid. Appendix 6.B: Finite Element Matrices: Unsteady Fluid. 7 Dynamics of Axially Moving Structures. 7.1 Axially Moving String. 7.2 Axially Moving Bernoulli-Euler Beam. 7.3 Axially Moving Timoshenko Beam. 7.4 Axially Moving Thin Plates. Appendix 7.A: Finite Element Matrices for Axially Moving String. Appendix 7.B: Finite Element Matrices for Axially Moving Bernoulli-Euler Beam. Appendix 7.C: Finite Element Matrices for Axially Moving Timoshenko Beam. Appendix 7.D: Finite Element Matrices for Axially Moving Plate. 8 Dynamics of Rotor Systems. 8.1 Governing Equations. 8.2 Spectral Element Modeling. 8.3 Finite Element Model. 8.4 Numerical Examples. Appendix 8.A: Finite Element Matrices for the Transverse Bending Vibration. 9 Dynamics of Multi-Layered Structures. 9.1 Elastic-Elastic Two-Layer Beams. 9.2 Elastic-Viscoelastic-elastic-Three-Layer (PCLD) Beams. Appendix 9.A: Finite Element Matrices for the Elastic-Elastic Two-Layer Beam. Appendix 9.B: Finite Element Matrices for the Elastic-VEM-Elastic Three-Layer Beam. 10 Dynamics of Smart Structures. 10.1 Elastic-Piezoelectric Two-Layer Beams. 10.2 Elastic-Viscoelastic-Piezoelctric Three-Layer (ACLD) Beams. 11 Dynamics of Composite Laminated Structures. 11.1 Theory of Composite Mechanics. 11.2 Equations of Motion for Composite Laminated Beams. 11.3 Dynamics of Axial-Bending-Shear Coupled Composite Beams. 11.4 Dynamics of Bending-Torsion-Shear Coupled Composite Beams. Appendix 11.A: Finite Element Matrices for Axial-Bending-Shear Coupled Composite Beams. Appendix 11.B: Finite Element Matrices for Bending-Torsion-Shear Coupled Composite Beams. 12 Dynamics of Periodic Lattice Structures. 12.1 Continuum Modeling Method. 12.2 Spectral Transfer Matrix Method. 13 Biomechanics: Blood Flow Analysis. 13.1 Governing Equations. 13.2 Spectral Element Modeling: I. Finite Element. 13.3 Spectral Element Modeling: II. Semi-Infinite Element. 13.4 Assembly of Spectral Elements. 13.5 Finite Element Model. 13.6 Numerical Examples. Appendix 13.A: Finite Element Model for the 1-D Blood Flow. 14 Identification of Structural Boundaries and Joints. 14.1 Identification of Non-Ideal Boundary Conditions. 14.2 Identification of Joints. 15 Identification of Structural Damage. 15.1 Spectral Element Modeling of a Damaged Structure. 15.2 Theory of Damage Identification. 15.3 Domain-Reduction Method. 16 Other Applications. 16.1 SEM-FEM Hybrid Method. 16.2 Identification of Impact Forces. 16.3 Other Applications. References. Index.

New Cartesian grid methods for interface problems using the finite element formulation

Abstract: New finite element methods based on Cartesian triangulations are presented for two dimensional elliptic interface problems involving discontinuities in the coefficients. The triangulations in these methods do not need to fit the interfaces. The basis functions in these methods are constructed to satisfy the interface jump conditions either exactly or approximately. Both non-conforming and conforming finite element spaces are considered. Corresponding interpolation functions are proved to be second order accurate in the maximum norm. The conforming finite element method has been shown to be convergent. With Cartesian triangulations, these new methods can be used as finite difference methods. Numerical examples are provided to support the methods and the theoretical analysis.

Electromagnetic Modeling by Finite Element Methods

Abstract: PREFACE MATHEMATICAL PRELIMINARIES Introduction The Vector Notation Vector Derivation The Gradient The Divergence The Rotational Second-Order Operators Application of Operators to More than One Function Expressions in Cylindrical and Spherical Coordinates MAXWELL EQUATIONS, ELECTROSTATICS, MAGNETOSTATICS, AND MAGNETODYNAMIC FIELDS Introduction The EM Quantities Local Form of the Equations The Anisotropy The Approximation of Maxwell's Equations The Integral Form of Maxwell's Equations Electrostatic Fields Magnetostatic Fields Magnetodynamic Fields BRIEF PRESENTATION OF THE FINITE ELEMENT METHOD Introduction The Galerkin Method - Basic Concepts A First-Order Finite Element Program Generalization of the Finite Element Method Numerical Integration Some 2D Finite Elements Coupling Different Finite Elements Calculation of Some Terms in the Field Equation A Simplified 2D Second-Order Finite Element Program THE FINITE ELEMENT METHOD APPLIED TO 2D ELECTROMAGNETIC CASES Introduction Some Static Cases Application to 2D Eddy Current Problems Axi-Symmetric Application Advantages and Limitation of 2D Formulations Non-Linear Applications Geometric Repetition of Domains Thermal Problems Voltage-Fed Electromagnetic Devices Static Examples Dynamic Examples COUPLING OF FIELD AND ELECTRICAL CIRCUIT EQUATIONS Introduction Electromagnetic Equations Equations for Different Conductor Configurations Connections Between Electromagnetic Devices and External Feeding Circuits Examples MOVEMENT MODELING FOR ELECTRICAL MACHINES Introduction The Macro-Element The Moving Band The Skew Effect in Electrical Machines Using 2D Simulation Examples INTERACTION BETWEEN ELECTROMAGNETIC AND MECHANICAL FORCES Introduction Methods Based on Direct Formulations Methods Based on the Force Density Electrical Machine Vibrations Originated by Magnetic Forces Example of Coupling Between the Field and Circuit Equations, Including Mechanical Transients IRON LOSSES Introduction Eddy Current Losses Hysteresis Anomalous or Excess Losses Total Iron Losses The Jiles-Atherton Model The Inverse Jiles-Atherton Model Including Iron Losses in Finite Element Calculations BIBLIOGRAPHY INDEX

On the construction of blending elements for local partition of unity enriched finite elements

Abstract: For computational efficiency, partition of unity enrichments are preferably localized to the sub-domains where they are needed. It is shown that an appropriate construction of the elements in the blending area, the regionwhere the enriched elements blend to unenriched elements, is often crucial for good performance of local partition of unity enrichments. An enhanced strain formulation is developed which leads to good performance; the optimal rate of convergence is achieved. For polynomial enrichments, it is shown that a proper choice of the finite element shape functions and partition of unity shape functions also improves the accuracy and convergence. The methods are illustrated by several examples. The examples deal primarily with the signed distance function enrichment for treating discontinuous derivatives inside an element, but other enrichments are also considered. Results show that both methods provide optimal rates of convergence.

An Extended Finite Element Method for Two-Phase Fluids

Abstract: An extended finite element method with arbitrary interior discontinuous gradients is applied to two-phase immiscible flow problems. The discontinuity in the derivative of the velocity field is introduced by an enrichment with an extended basis whose gradient is discontinuous across the interface. Therefore, the finite element approximation can capture the discontinuities at the interface without requiring the mesh to conform to the interface, eliminating the need for remeshing. The equations for incompressible flow are solved by a fractional step method where the advection terms are stabilized by a characteristic Galerkin method. The phase interfaces are tracked by level set functions which are discretized by the same finite element mesh and are updated via a stabilized conservation law. The method is demonstrated in several examples

Finite volume methods for the simulation of skeletal muscle

Abstract: Since it relies on a geometrical rather than a variational framework, many find the finite volume method (FVM) more intuitive than the finite element method (FEM). We show that the FVM allows one to interpret the stress inside a tetrahedron as a simple "multidimensional force" pushing on each face. Moreover, this interpretation leads to a heuristic method for calculating the force on each node, which is as simple to implement and comprehend as masses and springs. In the finite volume spirit, we also present a geometric rather than interpolating function definition of strain. We use the FVM and a quasi-incompressible, transversely isotropic, hyperelastic constitutive model to simulate contracting muscle tissue. B-spline solids are used to model fiber directions, and the muscle activation levels are derived from key frame animations.

An extended finite element method with higher-order elements for curved cracks

Abstract: A finite element method for linear elastic fracture mechanics using enriched quadratic interpolations is presented. The quadratic finite elements are enriched with the asymptotic near tip displacement solutions and the Heaviside function so that the finite element approximation is capable of resolving the singular stress field at the crack tip as well as the jump in the displacement field across the crack face without any significant mesh refinement. The geometry of the crack is represented by a level set function which is interpolated on the same quadratic finite element discretization. Due to the higher-order approximation for the crack description we are able to represent a crack with curvature. The method is verified on several examples and comparisons are made to similar formulations using linear interpolants.

Extended finite element method for quasi‐brittle fracture

Abstract: A methodology for the simulation of quasi-static cohesive crack propagation in quasi-brittle materials is presented. In the framework of the recently proposed extended finite element method, the partition of unity property of nodal shape functions has been exploited to introduce a higher-order displacement discontinuity in a standard finite element model. In this way, a cubic displacement discontinuity, able to reproduce the typical cusp-like shape of the process zone at the tip of a cohesive crack, is allowed to propagate without any need to modify the background finite element mesh. The effectiveness of the proposed method has been assessed by simulating mode-I and mixed-mode experimental tests. Copyright © 2003 John Wiley & Sons, Ltd.

Hierarchic finite element bases on unstructured tetrahedral meshes

Abstract: The problem of constructing hierarchic bases for finite element discretization of the spaces H1, H(curl), H(div) and L2 on tetrahedral elements is addressed A simple and efficient approach to ensuring conformity of the approximations across element interfaces is described Hierarchic bases of arbitrary polynomial order are presented It is shown how these may be used to construct finite element approximations of arbitrary, non-uniform, local order approximation on unstructured meshes of curvilinear tetrahedral elements

Vector level sets for description of propagating cracks in finite elements

Abstract: A new level set method is developed for describing surfaces that are frozen behind a moving front, such as cracks. In this formulation, the level set is described in two dimensions by a three-tuple: the sign of the level set function and the components of the closest point projection to the surface. The update of the level set is constructed by geometric formulas, which are easily implemented. Results are given for growth of lines in two dimensions that show the method is very accurate. The method combines very naturally with the extended finite element method (XFEM) where the discontinuous enrichment for cracks is best described in terms of level set functions. Examples of crack growth simulations obtained by combining this level set method with the extended finite element method are given. Copyright © 2003 John Wiley & Sons, Ltd.

Stress intensity factor analysis of interface cracks using X‐FEM

Abstract: The extended finite element method (X-FEM) proposed by Belytschko et al. (International Journal for Numerical Methods in Engineering 1999; 45: 602; 1999; 46: 131; 2001; 50: 993) uses interpolation functions based on the concept of partition of unity, and considers the asymptotic solution and the discontinuity of displacement fields near a crack independently of the finite element mesh. This paper describes the application of X-FEM to stress analyses of structures containing interface cracks between dissimilar materials. In X-FEM, an interface crack can be modelled by locally changing an interpolation function in the element near a crack. The energy release rate should be separated into individual stress intensity factors, K1 and K2, because the stress field around the interface crack has mixed modes coupled with mode-I and mode-II. For this purpose, various evaluation methods used in conjunction with numerical methods such as FEM and BEM are reviewed. These methods are examined in numerical examples of elastostatic analyses of structures containing interface cracks using X-FEM. The numerical results show that X-FEM is an effective method for performing stress analyses and evaluating stress intensity factors in problems related to bi-material fractures. Copyright © 2003 John Wiley & Sons, Ltd.

Brittle fracture in polycrystalline microstructures with the extended finite element method

Abstract: A two-dimensional numerical model of microstructural effects in brittle fracture is presented, with an aim towards the understanding of toughening mechanisms in polycrystalline materials such as ceramics. Quasi-static crack propagation is modelled using the extended finite element method (X-FEM) and microstructures are simulated within the framework of the Potts model for grain growth. In the X-FEM, a discontinuous function and the two-dimensional asymptotic crack-tip displacement fields are added to the finite element approximation to account for the crack using the notion of partition of unity. This enables the domain to be modelled by finite elements with no explicit meshing of the crack surfaces. Hence, crack propagation can be simulated without any user-intervention or the need to remesh as the crack advances. The microstructural calculations are carried out on a regular lattice using a kinetic Monte Carlo algorithm for grain growth. We present a novel constrained Delaunay triangulation algorithm with grain boundary smoothing to create a finite element mesh of the microstructure. The fracture properties of the microstructure are characterized by assuming that the critical fracture energy of the grain boundary (Gcgb) is different from that of the grain interior (Gci). Numerical crack propagation simulations for varying toughness ratios Gcgb/Gci are presented, to study the transition from the intergranular to the transgranular mode of crack growth. This study has demonstrated the capability of modelling crack propagation through a material microstructure within a finite element framework, which opens-up exciting possibilities for the fracture analysis of functionally graded material systems. Copyright © 2003 John Wiley & Sons, Ltd.

Mesh independent modelling of cracks by using higher order shape functions

Abstract: A mesh independent crack modelling approach based on displacement approximation with higher order shape functions is proposed. The Heaviside step function based local enrichment method, known as eXtended Finite Element Method, is modified by replacing the step function with a higher order shape functions approximation. Polynomial B-spline approximation functions are used in the present paper. An advantage of the proposed method is that its implementation only involves integration of the products of original shape functions and their derivatives and does not require modification of the integration domains. A volume integral based expression is proposed to calculate the effective surface area of the crack modelled by using an approximate step function. It is shown to give the actual crack surface area in the limit of the approximate step function approaching the Heaviside function. The convergence and accuracy of the method is illustrated in examples of transverse and oblique (with respect to loading direction) crack problems in rectangular plates. Uniaxial tension of a unidirectional composite with an open hole is considered. Hoop stress relaxation due to longitudinal splitting is successfully modelled by the method proposed and compared to direct modelling by using ANSYS software. Published in 2002 by John Wiley & Sons, Ltd.