Mixed Dimensional Coupling in Finite Element Stress Analysis

Many analysis models utilize finite elements of reduced dimension. However, to capture stress concentrations at local details, it would be desirable to combine the reduced dimensional element types with higher dimensional elements in a single finite element model. It is therefore important in such cases to integrate into the analyses some scheme for coupling the element types that conforms to the governing equations of the problem. In this paper, a novel method that can correctly couple beams to solids, beams to shells and shells to solids for elastic problems is presented. The approach adopted is to equate the work done on either side of the interface between dimensions, and this leads to multi-point constraint equations, thus providing a relationship among nodal degrees of freedom between the differing element types. Example results show that the proposed technique does not introduce any spurious stresses at the dimensional interfaces. ID="A1" Correspondence and offprint requests to: C. G. Armstrong, School of Mechanical and Manufacturing Engineering, The Queen's University of Belfast, Ashby Building, Stranmillis Road, Belfast BT9 5AH, Northern Ireland. E-mail: c.armstrong@qub.ac.uk     

A mesh-geometry based method for coupling 1D and 3D elements

When conducting a finite element analysis (FEA) one way to reduce the total number of degrees of freedom is to use a mixed-dimensional model. Using beam elements to model long and slender components can significantly reduce the total number of elements. Problems arise when trying to connect elements with different dimensions in part due to incompatible degrees of freedom between different types of finite elements. This paper focuses on problems that occur in coupling beams and solids, which means coupling 1D and 3D finite elements. This paper presents a mesh-based solution to these problems only using specific arrangements of classical 1D and 3D finite elements without requiring the use of additional constraint equations. Two alternative solutions are detailed, evaluated and compared in this paper through series of computational experiments. The implementation of both solutions is also presented and involves mesh and geometry processing operations along with an adaptation of boundary representation (BREP) classical data structures.     

MooNMD – a program package based on mapped finite element methods

The basis of mapped finite element methods are reference elements where the components of a local finite element are defined. The local finite element on an arbitrary mesh cell will be given by a map from the reference mesh cell. This paper describes some concepts of the implementation of mapped finite element methods. From the definition of mapped finite elements, only local degrees of freedom are available. These local degrees of freedom have to be assigned to the global degrees of freedom which define the finite element space. We will present an algorithm which computes this assignment. The second part of the paper shows examples of algorithms which are implemented with the help of mapped finite elements. In particular, we explain how the evaluation of integrals and the transfer between arbitrary finite element spaces can be implemented easily and computed efficiently. Communicated by: M.S. Espedal, A. Quarteroni, A. Sequeira     

MooNMD – a program package based on mapped finite element methods

The basis of mapped finite element methods are reference elements where the components of a local finite element are defined. The local finite element on an arbitrary mesh cell will be given by a map from the reference mesh cell. This paper describes some concepts of the implementation of mapped finite element methods. From the definition of mapped finite elements, only local degrees of freedom are available. These local degrees of freedom have to be assigned to the global degrees of freedom which define the finite element space. We will present an algorithm which computes this assignment. The second part of the paper shows examples of algorithms which are implemented with the help of mapped finite elements. In particular, we explain how the evaluation of integrals and the transfer between arbitrary finite element spaces can be implemented easily and computed efficiently.     

A computer algebra based finite element development environment

A finite element development environment based on the technical computing program Mathematica is described. The environment is used to automatically program standard element formulations and develop new elements with novel features. Source code can also be exported in a format compatible with commercial finite element program user-element facilities. The development environment is demonstrated for three mixed Petrov–Galerkin plane stress elements: a standard formulation, an advanced formulation incorporating rotational degrees of freedom and a standard formulation in which the stiffness matrix is integrated analytically, before being exported as ANSYS user elements. The results presented illustrate the accuracy of the standard mixed formulation element and the enhancement of performance when rotational degrees of freedom are added. Further, the analytically integrated element shows that computational requirements can be greatly reduced when analytical integration schemes are used in the formation.     

A simple method to impose rotations and concentrated moments on ANC beams

Recently introduced ANC beam elements furnish a simple formulation that allows to solve nonlinear problems of beams, including those with large displacements and strains, as well as complex nonlinear (inelastic) materials. The success and simplicity of these finite elements is mainly due to the fact that the only nodal degrees of freedom that they employ are displacements, and rotations are thus completely avoided. This in turn makes it very difficult to apply concentrated moments or to impose rotations at specific nodes of a finite element mesh. In this article, we present a simple enhancement to this beam formulation that allows to apply these two types of boundary conditions in a simple manner, making ANC beam elements more versatile for both multibody and structural applications.     

Symmetry considerations for anisotropic shells

The different types of symmetry exhibited by anisotropic shells for various loadings and boudary conditions are identified, and a simple procedure is presented for exploiting these symmetries in the finite element analysis. Examples are given of anisotropic cylindrical and doubly-curved shells where use of symmetry can significantly reduce the number of independent degrees of freedom in their finite element models.     

Isogeometric shell discretizations for flexible multibody dynamics

This work aims at including nonlinear elastic shell models in a multibody framework. We focus our attention to Kirchhoff–Love shells and explore the benefits of an isogeometric approach, the latest development in finite element methods, within a multibody system. Isogeometric analysis extends isoparameteric finite elements to more general functions such as B-splines and NURBS (Non-Uniform Rational B-Splines) and works on exact geometry representations even at the coarsest level of discretizations. Using NURBS as basis functions, high regularity requirements of the shell model, which are difficult to achieve with standard finite elements, are easily fulfilled. A particular advantage is the promise of simplifying the mesh generation step, and mesh refinement is easily performed by eliminating the need for communication with the geometry representation in a CAD (Computer-Aided Design) tool. Target applications are wind turbine blades and twist beam rear suspensions. First numerical examples demonstrate an impressive convergence behavior of the isogeometric approach even for a coarse mesh, while offering substantial savings with respect to the number of degrees of freedom.     

Three- and four-noded planar elements using absolute nodal coordinate formulation

This paper investigates two new types of planar finite elements containing three and four nodes. These elements are the reduced forms of the spatial plate elements employing the absolute nodal coordinate approach. Elements of the first type use translations of nodes and global slopes as nodal coordinates and have 18 and 24 degrees of freedom. The slopes facilitate the prevention of the shear locking effect in bending problems. Furthermore, the slopes accurately describe the deformed shape of the elements. Triangular and quadrilateral elements of the second type use translational degrees of freedom only and, therefore, can be utilized successfully in problems without bending. These simple elements with 6 and 8 degrees of freedom are identical to the elements used in conventional formulation of the finite element method from the kinematical point of view. Similarly to the famous problem called “flying spaghetti” which is used often as a benchmark for beam elements, a kind of “flying lasagna” is simulated for the planar elements. Numerical results of simulations are presented.     

Stabilized finite element formulation for elastic–plastic finite deformations

This paper presents a stabilized finite element formulation for nearly incompressible finite deformations in hyperelastic–plastic solids, such as metals. An updated Lagrangian finite element formulation is developed where mesh dependent terms are added to enhance the stability of the mixed finite element formulation. This formulation circumvents the restriction on the displacement and pressure fields due to the Babuška–Brezzi condition and provides freedom in choosing interpolation functions in the incompressible or nearly incompressible limit, typical in metal forming applications. Moreover, it facilitates the use of low order simplex elements (i.e. P1/P1), reducing the degrees of freedom required for the solution in the incompressible limit when stable elements are necessary. Linearization of the weak form is derived for implementation into a finite element code. Numerical experiments with P1/P1 elements show that the method is effective in incompressible conditions and can be advantageous in metal forming analysis.