scispace - formally typeset
Search or ask a question

Showing papers on "Finite element method published in 1997"


Journal ArticleDOI
TL;DR: This paper studies a multiscale finite element method for solving a class of elliptic problems arising from composite materials and flows in porous media, which contain many spatial scales and proposes an oversampling technique to remove the resonance effect.

1,825 citations


Book ChapterDOI
01 Jan 1997
TL;DR: This chapter introduces the finite element method (FEM) as a tool for solution of classical electromagnetic problems and discusses the main points in the application to electromagnetic design, including formulation and implementation.
Abstract: This chapter introduces the finite element method (FEM) as a tool for solution of classical electromagnetic problems. Although we discuss the main points in the application of the finite element method to electromagnetic design, including formulation and implementation, those who seek deeper understanding of the finite element method should consult some of the works listed in the bibliography section.

1,820 citations


Journal ArticleDOI
TL;DR: This paper extends a discontinuous finite element discretization originally considered for hyperbolic systems such as the Euler equations to the case of the Navier?Stokes equations by treating the viscous terms with a mixed formulation, and finds the method is ideally suited to compute high-order accurate solution of theNavier?

1,750 citations


Journal ArticleDOI
TL;DR: In this article, the smoothed particle hydrodynamics (SPH) method is extended to model incompressible flows of low Reynolds number, and the results show that the SPH results exhibit small pressure fluctuations near curved boundaries.

1,696 citations


Book
01 Jan 1997
TL;DR: In this paper, the authors present a one-dimensional analysis of fiber-reinforced composite materials and their properties, including the properties of the components of a Lamina and their relationship with other components.
Abstract: Introduction and Mathematical Preliminaries Fiber-Reinforced Composite Materials. Vectors and Tensors. Matrices. Transformation of Vector and Tensor Components. Integral Relations. Equations of Anisotropic Elasticity Classification of Equations. Kinematics. Kinetics. Constitutive Equations. Equations of Thermoelasticity and Electroelasticity. Summary. Virtual Work Principles and Variational Methods Virtual Work. The Variational Operator and Functionals. Extrema of Functionals. Virtual Work Principles. Variational Methods. Summary. Introduction to Composite Materials Basic Concepts and Terminology. Constitutive Equations of a Lamina. Transformation of Stresses and Strains. Plane Stress Constitutive Relations. Classical and First-Order Theories of Laminated Composite Plates Introduction. An Overview of ESL Laminate Theories. The Classical Laminated Plate Theory. The First-Order Laminated Plate Theory. Stiffness Characteristics for Selected Laminates. One-Dimensional Analysis of Laminated Plates Introduction. Analysis of Laminated Beams Using CLPT. Analysis of Laminated Beams Using FSDT. Cylindrical Bending Using CLPT. Cylindrical Bending Using FSDT. Closing Remarks. Analysis of Specially Orthotropic Plates Using CLPT Introduction. Bending of Simply Supported Plates. Bending of Plates with Two Opposite Edges Simply Supported. Bending of Rectangular Plates with Various Boundary Conditions. Buckling of Simply Supported Plates Under Compressive Loads. Buckling of Rectangular Plates Under Inplane Shear Load. Vibration of Simply Supported Plates. Buckling and Vibration of Plates with Two Parallel Edges Simply Supported. Closure. Analytical Solutions of Rectangular Laminates Using CLPT Governing Equations in Terms of Displacements. Admissible Boundary Conditions for the Navier Solutions. Navier Solutions of Antisymmetric Cross-Ply Laminates. The Navier Solutions of Antisymmetric Angle-Ply Laminates. The LTvy Solutions. Analysis of Midplane Symmetric Laminates. Transient Analysis. Summary. Analytical Solutions of Rectangular Laminates Using FSDT Introduction. Simply Supported Antisymmetric Cross-Ply Laminates. Simply Supported Antisymmetric Angle-Ply Laminates. Antisymmetric Cross-Ply Laminates with Two Opposite Edges Simply Supported. Antisymmetric Angle-Ply Laminates with Two Opposite Edges Simply Supported. Transient Solutions. Summary. Finite Element Analysis of Composite Laminates Introduction. Laminated Beams and Plate Strips by CLPT. Timoshenko Beam/Plate Theory. Numerical Results for Beams and Plate Strips. Finite Element Models of Laminated Plates (CLPT). Finite Element Models of Laminated Plates (FSDT). Summary. Refined Theories of Laminated Composite Plates Introduction. A Third-Order Plate Theory. Higher-Order Laminate Stiffness Characteristics. The Navier Solutions. LTvy Solutions of Cross-Ply Laminates. Displacement Finite Element Model. Layerwise Theories and Variable Kinematic Models In troduction. Development of the Theory. Finite Element Model. Variable Kinematic Formulations. Nonlinear Analysis of Composite Laminates Introduction. Nonlinear Stiffness Coefficients. Solution Methods for Nonlinear Algebraic Equations. Computational Aspects and Numerical Examples. Closure. Index Most chapters include Exercises and References for Additional Reading

1,344 citations


Journal ArticleDOI
TL;DR: In this article, a three-phase topology optimization method was proposed to find the distribution of material phases that optimizes an objective function (e.g. thermoelastic properties) subject to certain constraints, such as elastic symmetry or volume fractions of the constituent phases, within a periodic base cell.
Abstract: Composites with extremal or unusual thermal expansion coefficients are designed using a three-phase topology optimization method. The composites are made of two different material phases and a void phase. The topology optimization method consists in finding the distribution of material phases that optimizes an objective function (e.g. thermoelastic properties) subject to certain constraints, such as elastic symmetry or volume fractions of the constituent phases, within a periodic base cell. The effective properties of the material structures are found using the numerical homogenization method based on a finite-element discretization of the base cell. The optimization problem is solved using sequential linear programming. To benchmark the design method we first consider two-phase designs. Our optimal two-phase microstructures are in fine agreement with rigorous bounds and the so-called Vigdergauz microstructures that realize the bounds. For three phases, the optimal microstructures are also compared with new rigorous bounds and again it is shown that the method yields designed materials with thermoelastic properties that are close to the bounds. The three-phase design method is illustrated by designing materials having maximum directional thermal expansion (thermal actuators), zero isotropic thermal expansion, and negative isotropic thermal expansion. It is shown that materials with effective negative thermal expansion coefficients can be obtained by mixing two phases with positive thermal expansion coefficients and void.

827 citations


Journal ArticleDOI
TL;DR: This paper focuses its attention on two-dimensional steady-state problems and presents higher order accurate discontinuous finite element solutions on unstructured grids of triangles and shows that, in the presence of curved boundaries, a meaningful high-order accurate solution can be obtained only if a corresponding high- order approximation of the geometry is employed.

780 citations


Journal ArticleDOI
TL;DR: In this article, the plasticity concrete material model in the Lagrangian finite element code DYNA3D was assessed and enhanced, and the main modifications include the implementation of a third independent yield failure surface, removal of the tensile cutoff, and extension of the model in tension.

696 citations


Journal ArticleDOI
TL;DR: It is proved that, in two and more space dimensions, it is impossible to eliminate the so-called pollution effect of the Galerkin FEM.
Abstract: The development of numerical methods for solving the Helmholtz equation, which behaves robustly with respect to the wave number, is a topic of vivid research. It was observed that the solution of the Galerkin finite element method (FEM) differs significantly from the best approximation with increasing wave number. Many attempts have been presented in the literature to eliminate this lack of robustness by various modifications of the classical Galerkin FEM. However, we will prove that, in two and more space dimensions, it is impossible to eliminate this so-called pollution effect. Furthermore, we will present a generalized FEM in one dimension which behaves robustly with respect to the wave number.

631 citations


Journal ArticleDOI
TL;DR: The scaled boundary finite-element method, alias the consistent infinitesimal finite element cell method, is developed in this paper starting from the governing equations of linear elastodynamics and converges to the exact solution in the finite element sense in the circumferential directions.

626 citations



Journal ArticleDOI
TL;DR: A second-order accurate interface tracking method for the solution of incompressible Stokes flow problems with moving interfaces on a uniform Cartesian grid is presented and an implicit quasi-Newton method is developed that allows reasonable time steps to be used.
Abstract: A second-order accurate interface tracking method for the solution of incompressible Stokes flow problems with moving interfaces on a uniform Cartesian grid is presented. The interface may consist of an elastic boundary immersed in the fluid or an interface between two different fluids. The interface is represented by a cubic spline along which the singularly supported elastic or surface tension force can be computed. The Stokes equations are then discretized using the second-order accurate finite difference methods for elliptic equations with singular sources developed in our previous paper [SIAM J. Numer. Anal., 31(1994), pp. 1019--1044]. The resulting velocities are interpolated to the interface to determine the motion of the interface. An implicit quasi-Newton method is developed that allows reasonable time steps to be used.

Journal ArticleDOI
Paul Fischer1
TL;DR: A finite element-based additive Schwarz preconditioner using overlapping subdomains plus a coarse grid projection operator which is applied directly to the pressure on the interior Gauss points can yield as much as a fivefold reduction in simulation time over previously employed methods based upon deflation.

Journal ArticleDOI
TL;DR: In this article, a constitutive model was developed to reproduce the superelastic behavior of shape-memory alloys at finite strains, and the numerical implementation within a finite-element scheme was discussed in detail.

Journal ArticleDOI
TL;DR: In this article, nonlinear dynamic analysis of three-dimensional structural models is used more and more in the assessment of existing structures in zones of high seismic risk and in the development of new structures.
Abstract: In recent years nonlinear dynamic analysis of three-dimensional structural models is used more and more in the assessment of existing structures in zones of high seismic risk and in the development...

Journal ArticleDOI
TL;DR: In this paper, various time-domain finite-element methods for the simulation of transient electromagnetic wave phenomena are discussed, including nodal and edge/facet element basis functions, along with the numerical stability properties of the different methods.
Abstract: Various time-domain finite-element methods for the simulation of transient electromagnetic wave phenomena are discussed. Detailed descriptions of test/trial spaces, explicit and implicit formulations, nodal and edge/facet element basis functions are given, along with the numerical stability properties of the different methods. The advantages and disadvantages of mass lumping are examined. Finally, the various formulations are compared on the basis of their numerical dispersion performance.

Journal ArticleDOI
TL;DR: In this paper, the authors present a method to solve initial and boundary value problems using artificial neural networks, where a trial solution of the differential equation is written as a sum of two parts, the first part satisfies the boundary (or initial) conditions and contains no adjustable parameters.
Abstract: We present a method to solve initial and boundary value problems using artificial neural networks. A trial solution of the differential equation is written as a sum of two parts. The first part satisfies the boundary (or initial) conditions and contains no adjustable parameters. The second part is constructed so as not to affect the boundary conditions. This part involves a feedforward neural network, containing adjustable parameters (the weights). Hence by construction the boundary conditions are satisfied and the network is trained to satisfy the differential equation. The applicability of this approach ranges from single ODE's, to systems of coupled ODE's and also to PDE's. In this article we illustrate the method by solving a variety of model problems and present comparisons with finite elements for several cases of partial differential equations.

Journal ArticleDOI
TL;DR: This paper considers the so-called "inexact Uzawa" algorithm for iteratively solving linear block saddle point problems, and shows that the linear method always converges as long as the preconditioners defining the algorithm are properly scaled.
Abstract: In this paper, we consider the so-called "inexact Uzawa" algorithm for iteratively solving linear block saddle point problems. Such saddle point problems arise, for example, in finite element and finite difference discretizations of Stokes equations, the equations of elasticity, and mixed finite element discretization of second-order problems. We consider both the linear and nonlinear variants of the inexact Uzawa iteration. We show that the linear method always converges as long as the preconditioners defining the algorithm are properly scaled. Bounds for the rate of convergence are provided in terms of the rate of convergence for the preconditioned Uzawa algorithm and the reduction factor corresponding to the preconditioner for the upper left-hand block. In the case of nonlinear iteration, the inexact Uzawa algorithm is shown to converge provided that the nonlinear process approximating the inverse of the upper left-hand block is of sufficient accuracy. Bounds for the nonlinear iteration are given in terms of this accuracy parameter and the rate of convergence of the preconditioned linear Uzawa algorithm. Applications to the Stokes equations and mixed finite element discretization of second-order elliptic problems are discussed and, finally, the results of numerical experiments involving the algorithms are presented.

Journal ArticleDOI
TL;DR: In this paper, the Galerkin finite element solution to the Helmholtz equation has been studied for continuous and discrete spaces with Dirichlet-Robin boundary conditions, and the results on the phase difference between the exact and the Galerikin solution for arbitrary $p$ have been shown.
Abstract: In this paper, which is part II in a series of two, the investigation of the Galerkin finite element solution to the Helmholtz equation is continued. While part I contained results on the $h$ version with piecewise linear approximation, the present part deals with approximation spaces of order $p \ge 1$. As in part I, the results are presented on a one-dimensional model problem with Dirichlet--Robin boundary conditions. In particular, there are proven stability estimates, both with respect to data of higher regularity and data that is bounded in lower norms. The estimates are shown both for the continuous and the discrete spaces under consideration. Further, there is proven a result on the phase difference between the exact and the Galerkin finite element solutions for arbitrary $p$ that had been previously conjectured from numerical experiments. These results and further preparatory statements are then employed to show error estimates for the Galerkin finite element method (FEM). It becomes evident that the error estimate for higher approximation can---with certain assumptions on the data---be written in the same form as the piecewise linear case, namely, as the sum of the error of best approximation plus a pollution term that is of the order of the phase difference. The paper is concluded with a numerical evaluation.

01 Jan 1997
TL;DR: In this article, the problem of using measured modal parameters to detect and locate damage in plate-like structures is investigated, and a method based on the changes in the strain energy of the structure is discussed.
Abstract: In this paper the problem of using measured modal parameters to detect and locate damage in plate-like structures is investigated. Many methods exist for locating damage in a structure given the modal properties before and after damage. Unfortunately, many of these methods require a correlated finite element model or mass normalized mode shapes. If the modal properties are obtained using ambient excitation then the mode shapes will not be mass normalized. In this paper a method based on the changes in the strain energy of the structure will be discussed. This method has been successfully applied to beam-like structures, that is, structures characterized by one-dimensional curvature. In this paper the method will be generalized to plate-like structures that are characterized by two-dimensional curvature. This method only requires the mode shapes of the structure before and after damage. To evaluate the effectiveness of the method it will be applied to simulated data.

Journal ArticleDOI
TL;DR: An expanded mixed finite element approximation of second-order elliptic problems containing a tensor coefficient is presented, and it is shown that rates of convergence are retained for the finite difference method.
Abstract: We present an expanded mixed finite element approximation of second-order elliptic problems containing a tensor coefficient. The mixed method is expanded in the sense that three variables are explicitly approximated, namely, the scalar unknown, the negative of its gradient, and its flux (the tensor coefficient times the negative gradient). The resulting linear system is a saddle point problem. In the case of the lowest order Raviart--Thomas elements on rectangular parallelepipeds, we approximate this expanded mixed method by incorporating certain quadrature rules. This enables us to write the system as a simple, cell-centered finite difference method requiring the solution of a sparse, positive semidefinite linear system for the scalar unknown. For a general tensor coefficient, the sparsity pattern for the scalar unknown is a 9-point stencil in two dimensions and 19 points in three dimensions. Existing theory shows that the expanded mixed method gives optimal order approximations in the $L^2$- and $H^{-s}$-norms (and superconvergence is obtained between the $L^2$-projection of the scalar variable and its approximation). We show that these rates of convergence are retained for the finite difference method. If $h$ denotes the maximal mesh spacing, then the optimal rate is $O(h)$. The superconvergence rate $O(h^{2})$ is obtained for the scalar unknown and rate $O(h^{3/2})$ for its gradient and flux in certain discrete norms; moreover, the full $O(h^{2})$ is obtained in the strict interior of the domain. Computational results illustrate these theoretical results.

Journal ArticleDOI
TL;DR: In this paper, a geometrically non-linear version of the EAS-approach is applied which is based on the enhancement of the Green-Lagrange strains instead of the displacement gradient as originally proposed by Simo and Armero.
Abstract: Well-known finite element concepts like the Assumed Natural Strain (ANS) and the Enhanced Assumed Strain (EAS) techniques are combined to derive efficient and reliable finite elements for continuum based shell formulations. In the present study two aspects are covered: The first aspect focuses on the classical 5-parameter shell formulation with Reissner–Mindlin kinematics. The above-mentioned combinations, already discussed by Andelfinger and Ramm for the linear case of a four-node shell element, are extended to geometrical non-linearities. In addition a nine-node quadrilateral variant is presented. A geometrically non-linear version of the EAS-approach is applied which is based on the enhancement of the Green–Lagrange strains instead of the displacement gradient as originally proposed by Simo and Armero. In the second part elements are derived in a similar way for a higher order, so-called 7-parameter non-linear shell formulation which includes the thickness stretch of the shell (Buchter and Ramm). In order to avoid artificial stiffening caused by the three dimensional displacement field and termed ‘thickness locking’, special provisions for the thickness stretch have to be introduced. © 1997 John Wiley & Sons, Ltd.

Journal ArticleDOI
TL;DR: In this article, the use of finite element methods (FEM) in developing mechanistic variation simulation models for deformable sheet metal parts with complex two or three dimensional free form surfaces was proposed.
Abstract: Traditional variation analysis methods, such as Root Sum Square method and Monte Carlo simulation, are not applicable to sheet metal assemblies because of possible part deformation during the assembly process. This paper proposes the use of finite element methods (FEM) in developing mechanistic variation simulation models for deformable sheet metal parts with complex two or three dimensional free form surfaces. Mechanistic variation simulation provides improved analysis by combining engineering structure models and statistical analysis in predicting the assembly variation. Direct Monte Carlo simulation in FEM is very time consuming, because hundreds or thousands of FEM runs are required to obtain a realistic assembly distribution. An alternative method, based on the Method of Influence Coefficients, is developed to improve the computational efficiency, producing improvements by several orders of magnitude. Simulations from both methods yield almost identical results. An example illustrates the developed methods used for evaluating sheet metal assembly variation. The new approaches provide an improved understanding of sheet metal assembly processes.

Journal ArticleDOI
TL;DR: In this article, an advanced lumped parameter joint model is developed and identified by experimental investigations for an isolated bolted joint, which is implemented in a Finite Element program for calculating the dynamic response of assembled structures incorporating the influence of micro-and macroslip of several bolted joints.
Abstract: The nonlinear transfer behaviour of an assembled structure such as a large lightweight space structure is caused by the nonlinear influence of structural connections. Bolted or riveted joints are the primary source of damping compared to material damping, if no special damping treatment is added to the structure. Simulation of this damping amount is very important in the design phase of a structure. Several well known lumped parameter joint models used in the past to describe the dynamic transfer behaviour of isolated joints by Coulomb friction elements are capable of describing global states of slip and stick only. The present paper investigates the influence of joints by a mixed experimental and numerical strategy. A detailed Finite Element model is established to provide understanding of different slip-stick mechanisms in the contact area. An advanced lumped parameter model is developed and identified by experimental investigations for an isolated bolted joint. This model is implemented in a Finite Element program for calculating the dynamic response of assembled structures incorporating the influence of micro- and macroslip of several bolted joints.

Book
01 May 1997
TL;DR: In this article, the authors bring together concepts behind different branches of work on non-linear finite elements, such as plasticity, beams and rods in two-dimensional framework, shells; hyper-elasticity; rubber; kinematic hardening; yield criteria with volume effects; large rotations; three-dimensional beams and rod; and stability theory.
Abstract: From the Publisher: Taking an engineering rather than a mathematical bias, this work brings together concepts behind different branches of work on non-linear finite elements. Practically-oriented, it acknowledges the increasing role of computers in non-linear analysis by including a series of computer applications. Topics include plasticity; beams and rods in two-dimensional framework; shells; hyper-elasticity; rubber; kinematic hardening; yield criteria with volume effects; large rotations; three-dimensional beams and rods; and stability theory. Also examines advanced solution procedures for non-linear analysis such as line searchers, quasi-Newton and acceleration techniques, arc-length methods, automatic increments and re-starts.

Journal ArticleDOI
TL;DR: In this article, a new numerical method is presented for propagating elastic waves in heterogeneous earth media, based on spectral approximations of the wavefield combined with domain decomposition techniques.
Abstract: A new numerical method is presented for propagating elastic waves in heterogeneous earth media, based on spectral approximations of the wavefield combined with domain decomposition techniques. The flexibility of finite element techniques in dealing with irregular geologic structures is preserved, together with the high accuracy of spectral methods. High computational efficiency can be achieved especially in 3D calculations, where the commonly used finite-difference approaches are limited both in the frequency range and in handling strongly irregular geometries. The treatment of the seismic source, introduced via a moment tensor distribution, is thoroughly discussed together with the aspects associated with its numerical implementation. The numerical results of the present method are successfully compared with analytical and numerical solutions, both in 2D and 3D.

Journal ArticleDOI
J. N. Reddy1
TL;DR: In this article, a locking-free finite element model using the form of the exact solution of the Timoshenko beam theory is developed, which yields exact nodal values for the generalized displacements for constant material and geometric properties of beams.

Journal ArticleDOI
TL;DR: A mathematical description of cardiac anatomy is presented for use with finite element models of the electrical activation and mechanical function of the heart, and parameters defined at the nodes of the finite element mesh are fitted.
Abstract: A mathematical description of cardiac anatomy is presented for use with finite element models of the electrical activation and mechanical function of the heart. The geometry of the heart is given in terms of prolate spheroidal coordinates defined at the nodes of a finite element mesh and interpolated within elements by a combination of linear Lagrange and cubic Hermite basis functions. Cardiac microstructure is assumed to have three axes of symmetry: one aligned with the muscle fiber orientation (the fiber axis); a second set orthogonal to the fiber direction and lying in the newly identified myocardial sheet plane (the sheet axis); and a third set orthogonal to the first two, in the sheet-normal direction. The geometry, fiber-axis direction, and sheet-axis direction of a dog heart are fitted with parameters defined at the nodes of the finite element mesh. The fiber and sheet orientation parameters are defined with respect to the ventricular geometry such that 1) they can be applied to any heart of known dimensions, and 2) they can be used for the same heart at various states of deformation, as is needed, for example, in continuum models of ventricular contraction.

Journal ArticleDOI
TL;DR: In this article, a large class of dissipative materials is described by a time and frequency-dependent viscoelastic constitutive model and the derivation of the numerical model is given.
Abstract: Purely elastic material models have a limited validity. Generally, a certain amount of energy absorbing behaviour can be observed experimentally for nearly any material. A large class of dissipative materials is described by a time- and frequency-dependent viscoelastic constitutive model. Typical representatives of this type are polymeric rubber materials. A linear viscoelastic approach at small and large strains is described in detail and this makes a very efficient numerical formulation possible. The underlying constitutive structure is the generalized Maxwell-element. The derivation of the numerical model is given. It will be shown that the developed isotropic algorithmic material tensor is even valid for the current configuration in the case of large strains. Aspects of evaluating experimental investigations as well as parameter identification are considered. Finally, finite element simulations of time-dependent deformations of rubber structures using mixed elements are presented.