Showing papers on "Finite element method published in 2001"

A 99 line topology optimization code written in Matlab

Abstract: The paper presents a compact Matlab implementation of a topology optimization code for compliance minimization of statically loaded structures. The total number of Matlab input lines is 99 including optimizer and Finite Element subroutine. The 99 lines are divided into 36 lines for the main program, 12 lines for the Optimality Criteria based optimizer, 16 lines for a mesh-independency filter and 35 lines for the finite element code. In fact, excluding comment lines and lines associated with output and finite element analysis, it is shown that only 49 Matlab input lines are required for solving a well-posed topology optimization problem. By adding three additional lines, the program can solve problems with multiple load cases. The code is intended for educational purposes. The complete Matlab code is given in the Appendix and can be down-loaded from the web-site http://www.topopt.dtu.dk.

An optimal control approach to a posteriori error estimation in finite element methods

Abstract: This article surveys a general approach to error control and adaptive mesh design in Galerkin finite element methods that is based on duality principles as used in optimal control. Most of the existing work on a posteriori error analysis deals with error estimation in global norms like the ‘energy norm’ or the L2 norm, involving usually unknown ‘stability constants’. However, in most applications, the error in a global norm does not provide useful bounds for the errors in the quantities of real physical interest. Further, their sensitivity to local error sources is not properly represented by global stability constants. These deficiencies are overcome by employing duality techniques, as is common in a priori error analysis of finite element methods, and replacing the global stability constants by computationally obtained local sensitivity factors. Combining this with Galerkin orthogonality, a posteriori estimates can be derived directly for the error in the target quantity. In these estimates local residuals of the computed solution are multiplied by weights which measure the dependence of the error on the local residuals. Those, in turn, can be controlled by locally refining or coarsening the computational mesh. The weights are obtained by approximately solving a linear adjoint problem. The resulting a posteriori error estimates provide the basis of a feedback process for successively constructing economical meshes and corresponding error bounds tailored to the particular goal of the computation. This approach, called the ‘dual-weighted-residual method’, is introduced initially within an abstract functional analytic setting, and is then developed in detail for several model situations featuring the characteristic properties of elliptic, parabolic and hyperbolic problems. After having discussed the basic properties of duality-based adaptivity, we demonstrate the potential of this approach by presenting a selection of results obtained for practical test cases. These include problems from viscous fluid flow, chemically reactive flow, elasto-plasticity, radiative transfer, and optimal control. Throughout the paper, open theoretical and practical problems are stated together with references to the relevant literature.

Optimal scaling for various Metropolis-Hastings algorithms

Abstract: We review and extend results related to optimal scaling of Metropolis–Hastings algorithms. We present various theoretical results for the high-dimensional limit. We also present simulation studies which confirm the theoretical results in finite-dimensional contexts.

Finite element interface models for the delamination analysis of laminated composites: mechanical and computational issues

Abstract: The finite element analysis of delamination in laminated composites is addressed using interface elements and an interface damage law. The principles of linear elastic fracture mechanics are indirectly used by equating, in the case of single-mode delamination, the area underneath the traction/relative displacement curve to the critical energy release rate of the mode under examination. For mixed-mode delamination an interaction model is used which can fulfil various fracture criteria proposed in the literature. It is then shown that the model can be recast in the framework of a more general damage mechanics theory. Numerical results are presented for the analyses of a double cantilever beam specimen and for a problem involving multiple delamination for which comparisons are made with experimental results. Issues related with the numerical solution of the non-linear problem of the delamination are discussed, such as the influence of the interface strength on the convergence properties and the final results, the optimal choice of the iterative matrix in the predictor and the number of integration points in the interface elements. Copyright © 2001 John Wiley & Sons, Ltd.

Arbitrary discontinuities in finite elements

Abstract: A technique for modelling arbitrary discontinuities in finite elements is presented. Both discontinuities in the function and its derivatives are considered. Methods for intersecting and branching discontinuities are given. In all cases, the discontinuous approximation is constructed in terms of a signed distance functions, so level sets can be used to update the position of the discontinuities. A standard displacement Galerkin method is used for developing the discrete equations. Examples of the following applications are given: crack growth, a journal bearing, a non-bonded circular inclusion and a jointed rock mass. Copyright © 2001 John Wiley & Sons, Ltd.

Filters in topology optimization

Abstract: In this article, a modified (‘filtered’) version of the minimum compliance topology optimization problem is studied. The direct dependence of the material properties on its pointwise density is replaced by a regularization of the density field by the mean of a convolution operator. In this setting it is possible to establish the existence of solutions. Moreover, convergence of an approximation by means of finite elements can be obtained. This is illustrated through some numerical experiments. The ‘filtering’ technique is also shown to cope with two important numerical problems in topology optimization, checkerboards and mesh dependent designs. Copyright © 2001 John Wiley & Sons, Ltd.

A new method for modelling cohesive cracks using finite elements

Abstract: A model which allows the introduction of displacements jumps to conventional finite elements is developed. The path of the discontinuity is completely independent of the mesh structure. Unlike so-called ‘embedded discontinuity’ models, which are based on incompatible strain modes, there is no restriction on the type of underlying solid finite element that can be used and displacement jumps are continuous across element boundaries. Using finite element shape functions as partitions of unity, the displacement jump across a crack is represented by extra degrees of freedom at existing nodes. To model fracture in quasi-brittle heterogeneous materials, a cohesive crack model is used. Numerical simulations illustrate the ability of the method to objectively simulate fracture with unstructured meshes. Copyright © 2001 John Wiley & Sons, Ltd.

Application of finite element analysis in implant dentistry: A review of the literature

Abstract: Finite element analysis (FEA) has been used extensively to predict the biomechanical performance of various dental implant designs as well as the effect of clinical factors on implant success. By understanding the basic theory, method, application, and limitations of FEA in implant dentistry, the clinician will be better equipped to interpret results of FEA studies and extrapolate these results to clinical situations. This article reviews the current status of FEA applications in implant dentistry and discusses findings from FEA studies in relation to the bone-implant interface, the implant-prosthesis connection, and multiple-implant prostheses.

The finite element method and its reliability

Abstract: Preface 1. Introduction 2. Mathematical formulation of the model problem 3. The finite element method 4. Local behaviour in the finite element method 5. A-posteriori estimation of the error 6. Guaranteed a-posteriori error estimation, and a-posteriori estimation of the pollution error Appendix Index

A point interpolation method for two-dimensional solids

Abstract: A point interpolation method (PIM) is presented for stress analysis for two-dimensional solids. In the PIM, the problem domain is represented by properly scattered points. A technique is proposed to construct polynomial interpolants with delta function property based only on a group of arbitrarily distributed points. The PIM equations are then derived using variational principles. In the PIM, the essential boundary conditions can be implemented with ease as in the conventional finite element methods. The present PIM has been coded in FORTRAN. The validity and efficiency of the present PIM formulation are demonstrated through example problems. It is found that the present PIM is very easy to implement, and very flexible for obtained displacements and stresses of desired accuracy in solids. As the elements are not used for meshing the problem domain, the present PIM opens new avenues to develop adaptive analysis codes for stress analysis in solids and structures. Copyright © 2001 John Wiley & Sons, Ltd.

A point interpolation method for two-dimensional solids

Abstract: A Point Interpolation Method (PIM) is presented for stress analysis for two-dimensional solids. In the PIM, the problem domain is represented by properly scattered points. A technique is proposed to construct polynomial interpolants with delta function property based only on a group of arbitrarily distributed points. The PIM equations are then derived using variational principles. In the PIM, the essential boundary conditions can be implemented with ease as in the conventional Finite Element Methods. The present PIM has been coded in FORTRAN. The validity and efficiency of the present PIM formulation are demonstrated through example problems. It is found that the present PIM is very easy to implement, and very flexible for obtained displacements and stresses of desired accuracy in solids. As the elements are not used for meshing the problem domain, the present PIM opens new avenue to develop adaptive analysis codes for stress analysis in solids and structures.

Cross-sectional mapping of residual stresses by measuring the surface contour after a cut

Abstract: A powerful new method for residual stress measurement is presented. A part is cut in two, and the contour, or profile, of the resulting new surface is measured to determine the displacements caused by release of the residual stresses. Analytically, for example using a finite element model, the opposite of the measured contour is applied to the surface as a displacement boundary condition. By Bueckner's superposition principle, this calculation gives the original residual stresses normal to the plane of the cut. This contour method is more powerful than other relaxation methods because it can determine an arbitrary cross-sectional area map of residual stress, yet more simple because the stresses can be determined directly from the data without a tedious inversion technique. The new method is verified with a numerical simulation, then experimentally validated on a steel beam with a known residual stress profile.

Higher-Order Numerical Methods for Transient Wave Equations

Abstract: Solving efficiently the wave equations involved in modeling acoustic, elastic or electromagnetic wave propagation remains a challenge both for research and industry. To attack the problems coming from the propagative character of the solution, the author constructs higher-order numerical methods to reduce the size of the meshes, and consequently the time and space stepping, dramatically improving storage and computing times. This book surveys higher-order finite difference methods and develops various mass-lumped finite (also called spectral) element methods for the transient wave equations, and presents the most efficient methods, respecting both accuracy and stability for each sort of problem. A central role is played by the notion of the dispersion relation for analyzing the methods. The last chapter is devoted to unbounded domains which are modeled using perfectly matched layer (PML) techniques. Numerical examples are given.

Theoretical Numerical Analysis: A Functional Analysis Framework

Abstract: Preface 1 Linear Spaces 2 Linear Operators on Normed Spaces 3 Approximation Theory 4 Nonlinear Equations and Their Solution by Iteration 5 Finite Difference Method 6 Sobolev Spaces 7 Variational Formulations of Elliptic Boundary Value Problems 8 The Galerkin Method and Its Variants 9 Finite Element Analysis 10 Elliptic Variational Inequalities and Their Numerical Approximations 11 Numerical Solution of Fredholm Integral Equations of the Second Kind 12 Boundary Integral Equations References Index.

Dynamic real-time deformations using space & time adaptive sampling

Abstract: This paper presents a robust, adaptive method for animating dynamic visco-elastic deformable objects that provides a guaranteed frame rate. Our approach uses a novel automatic space and time adaptive level of detail technique, in combination with a large-displacement (Green) strain tensor formulation. The body is partitioned in a non-nested multiresolution hierarchy of tetrahedral meshes. The local resolution is determined by a quality condition that indicates where and when the resolution is too coarse. As the object moves and deforms, the sampling is refined to concentrate the computational load into the regions that deform the most. Our model consists of a continuous differential equation that is solved using a local explicit finite element method. We demonstrate that our adaptive Green strain tensor formulation suppresses unwanted artifacts in the dynamic behavior, compared to adaptive mass-spring and other adaptive approaches. In particular, damped elastic vibration modes are shown to be nearly unchanged for several levels of refinement. Results are presented in the context of a virtual reality system. The user interacts in real-time with the dynamic object through the control of a rigid tool, attached to a haptic device driven with forces derived from the method.

Aims, scope, methods, history and unified terminology of computer-aided topology optimization in structural mechanics

Abstract: Topology optimization of structures and composite continua has two main subfields: Layout Optimization (LO) deals with grid-like structures having very low volume fractions and Generalized Shape Optimization (GSO) is concerned with higher volume fractions, optimizing simultaneously the topology and shape of internal boundaries of porous or composite continua. The solutions for both problem classes can be exact/analytical or discretized/FE-based. This review article discusses FE-based generalized shape optimization, which can be classified with respect to the types of topologies involved, namely Isotropic-Solid/Empty (ISE), Anisotropic-Solid/Empty (ASE), and Isotropic-Solid/Empty/Porous (ISEP) topologies. Considering in detail the most important class of (i.e. ISE) topologies, the computational efficiency of various solution strategies, such as SIMP (Solid Isotropic Microstructure with Penalization), OMP (Optimal Microstructure with Penalization) and NOM (NonOptimal Microstructures) are compared. The SIMP method was proposed under the terms "direct approach" or "artificial density approach" by Bendsoe over a decade ago; it was derived independently, used extensively and promoted by the author's research group since 1990. The term "SIMP" was introducted by the author in 1992. After being out of favour with most other research schools until recently, SIMP is becoming generally accepted in topology optimization as a technique of considerable advantages. It seems, therefore, useful to review in greater detail the origins, theoretical background, history, range of validity and major advantages of this method.

Classical and Computational Solid Mechanics

Abstract: Tensor analysis stress tensor analysis of strain conservation laws elastic and plastic behaviour of materials linearized theory of elasticity solutions of problems in linearized theory of elasticity by potentials two-dimensional problems in linearized theory of elasticity variational calculus, energy theorems, Saint-Venant's principle Hamilton's principle, wave propagation, applications of generalized co-ordinates elasticity and thermodynamics irreversible thermodynamics and viscoelasticity thermoelasticity viscoelasticity large deformation incremental approach to solving some nonlinear problems finite element methods mixed and hybrid formulations finite element methods for plates and shells finite element modelling of nonlinear elasticity, viscoelasticity, plasticity, viscoplasticity and creep.

Discontinuous hp -Finite Element Methods for Advection-Diffusion-Reaction Problems

Abstract: We consider the hp-version of the discontinuous Galerkin finite element method (DGFEM) for second-order partial differential equations with nonnegative characteristic form. This class of equations includes second-order elliptic and parabolic equations, advection-reaction equations, as well as problems of mixed hyperbolic-elliptic-parabolic type. Our main concern is the error analysis of the method in the absence of streamline-diffusion stabilization. In the hyperbolic case, an hp-optimal error bound is derived; here, we consider only advection-reaction problems which satisfy a certain (standard) positivity condition. In the self-adjoint elliptic case, an error bound that is h-optimal and p-suboptimal by $\frac{1}{2}$ a power of p is obtained. These estimates are then combined to deduce an error bound in the general case. For elementwise analytic solutions the method exhibits exponential rates of convergence under p-refinement. The theoretical results are illustrated by numerical experiments.

High order finite difference and finite volume WENO schemes and discontinuous Galerkin methods for CFD

Abstract: In recent years high order numerical methods have been widely used in computational fluid dynamics (CFD), to effectively resolve complex flow features using meshes which are reasonable for today''s computers. In this paper we review and compare three types of high order methods being used in CFD, namely the weighted essentially non-oscillatory (WENO) finite difference methods, the WENO finite volume methods, and the discontinuous Galerkin (DG) finite element methods. We summarize the main features of these methods, from a practical user''s point of view, indicate their applicability and relative strength, and show a few selected numerical examples to demonstrate their performance on illustrative model CFD problems.

Computational Mechanics of the Heart

Abstract: Finite elasticity theory combined with finite element analysis provides the framework for analysing ventricular mechanics during the filling phase of the cardiac cycle, when cardiac cells are not actively contracting. The orthotropic properties of the passive tissue are described here by a “pole–zero” constitutive law, whose parameters are derived in part from a model of the underlying distributions of collagen fibres. These distributions are based on our observations of the fibrous-sheet laminar architecture of myocardial tissue. We illustrate the use of high order (cubic Hermite) basis functions in solving the Galerkin finite element stress equilibrium equations based on this orthotropic constitutive law and for incorporating the observed regional distributions of fibre and sheet orientations. Pressure–volume relations and 3D principal strains predicted by the model are compared with experimental observations. A model of active tissue properties, based on isolated muscle experiments, is also introduced in order to predict transmural distributions of 3D principal strains at the end of the contraction phase of the cardiac cycle. We end by offering a critique of the current model of ventricular mechanics and propose new challenges for future modellers.