Follower loads for high order finite elements

Finite element modelling of hydrostatic compaction where the applied pressure acts normal to the deformed surface requires a geometric nonlinear formulation and follower load terms [1, 5, 7]. These concepts are applied to high order [6] (p-FEM) elements with hierarchic shape functions. Applying the blending function method allows to precisely describe curved boundaries on coarse meshes. High order elements exhibit good performance even for high aspect ratios and strong distortion and therefore allow an efficient discretization of thin-walled structures. Since high order finite elements are less prone to locking effects a pure displacement-based formulation can be chosen. After introducing the basic concept of the p-version the application of follower loads to geometrically nonlinear high order elements is presented. For the numerical solution the displacement based formulation is linearized yielding the basis for a Newton-Raphson iteration. The accuracy and performance of the high order finite element scheme is demonstrated by a numerical example. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Isogeometrical approach for cable type structures allowing large sliding contact

A special combination of two finite elements approaches is considered in order to derive a formulation suitable for cable type structures. First, a finite beam element allowing arbitrarily large rotations together with a high order approximation technique including a C¹-continuous approximation is used to model cable behavior. Afterwards, a covariant approach for a curve-to-curve contact is exploited to incorporate contact between a pair of cable elements. This leads to a cable type element with contact. In numerical examples both “soft” and “hard” cables are represented. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Thermo-mechanical modeling of turbine blades taking into account structural defects

In turbine blades of aero-engines typical defects are cracks due to high mechanical and thermal loads. The extended finite element method (XFEM) is used for simulations of fracture mechanics problems with cracks. Discontinuities in the displacement and temperature field are allowed and the crack opening displacement and crack tip stress field are reproduced accurately. Since crack closure and non-physical penetration of the crack surfaces may occur under certain load conditions, it becomes necessary to enforce the non-penetration condition for crack surfaces. This contact formulation is assumed to be frictionless. The node-to-segment approach proposed in [3] is extended to ten-node tetrahedral elements with quadratic shape functions. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Different weak FE‐formulations for 2‐D frictionless Large Deformation Contact based on the Mortar Method

The application of the mortar method in contact mechanics is motivated by the limited use of well known elements, for example the node‐to‐segment (NTS) element, see [1]. The NTS element contains a strong projection of the displacements from one contacting body to the next. Coupling of this element type with higher order shape functions leads to a loss of accuracy of contact stresses. In contrast to this, the mortar element has the advantage of a weak projection. Therefore, consistent coupling with continuum elements of higher order is possible. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical aspects of the XFEM — The XFEM p–Version

The XFEM is known to approximate the displacements and stresses around a crack tip in a very efficient way. But as we will present in this paper we have to deal with a phenomenon coming along with this method that compels us to use higher order shape functions for those elements that are enriched by the crack tip functions. For the computation of the stress–intensity–factors we are using a J–integral over a circular domain Ω. The accuracy of the results depend on • the radius of Ω • the number of elements used in the XFEM computation • the number of nodes which were enriched by the crack tip functions (number of layers) and • the shape functions which were used for the standard FE term For more information about the XFEM we refer to [1]. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

p-version finite elements and finite cells for finite strain elastoplastic problems

In this paper, the p-version finite element method and its fictitious domain extension, the finite cell method, are extended to the finite strain J₂ plasticity. High-order shape functions are used for the finite element approximation of volume-preserving plastic dominated deformations. The accuracy and efficiency of p-version elements and cells in the finite plastic strain range are evaluated by the computation of two benchmark problems. It is shown that they provide locking free behavior and simplified meshing. These results are verified in comparison with the results of h-version elements in F-bar formulation. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Weak Galerkin finite element methods for the biharmonic equation on polytopal meshes

A new weak Galerkin (WG) finite element method is introduced and analyzed in this article for the biharmonic equation in its primary form. This method is highly robust and flexible in the element construction by using discontinuous piecewise polynomials on general finite element partitions consisting of polygons or polyhedra of arbitrary shape. The resulting WG finite element formulation is symmetric, positive definite, and parameter‐free. Optimal order error estimates in a discrete H² norm is established for the corresponding WG finite element solutions. Error estimates in the usual L² norm are also derived, yielding a suboptimal order of convergence for the lowest order element and an optimal order of convergence for all high order of elements. Numerical results are presented to confirm the theory of convergence under suitable regularity assumptions. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1003–1029, 2014     

A combined nodal continuous-discontinuous finite element formulation for the Maxwell problem

Continuous Galerkin formulations are appealing due to their low computational cost, whereas discontinuous Galerkin formulation facilitate adaptative mesh refinement and are more accurate in regions with jumps of physical parameters. Since many electromagnetic problems involve materials with different physical properties, this last point is very important. For this reason, in this article we have developed a combined cG-dG formulation for Maxwell’s problem that allows arbitrary finite element spaces with functions continuous in patches of finite elements and discontinuous on the interfaces of these patches. In particular, the second formulation we propose comes from a novel continuous Galerkin formulation that reduces the amount of stabilization introduced in the numerical system. In all cases, we have performed stability and convergence analyses of the methods. The outcome of this work is a new approach that keeps the low CPU cost of recent nodal continuous formulations with the ability to deal with coefficient jumps and adaptivity of discontinuous ones. All these methods have been tested using a problem with singular solution and another one with different materials, in order to prove that in fact the resulting formulations can properly deal with these problems.     

hp-Adaptive Finite Element Methods for Variational Inequalities

In this work, we combine an hp–adaptive strategy with a posteriori error estimates for variational inequalities, which are given by contact problems. The a posteriori error estimates are obtained using a general approach based on the saddle point formulation of contact problems and making use of a yposteriori error estimates for variational equations. Error estimates are presented for obstacle problems and Signorini problems with friction. Numerical experiments confirm the reliability of the error estimates for finite elements of higher order. The use of the hp–adaptive strategy leads to meshes with the same characteristics as geometric meshes and to exponential convergence. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A spatial shear deformable beam finite element based on the absolute nodal coordinate formulation

In the present paper a three-dimensional beam finite element undergoing large deformations is proposed. Since the definition of the proposed finite element is based on the absolute nodal coordinate formulation (ANCF), no rotational coordinates occur in the formulation. In the current approach, the orientation of the cross section is parameterized by means of slope vectors. Since those are no unit vectors, the cross-section can deform, similar to existing thick beam and shell elements. The nodal displacements and the directional derivatives of the displacements are chosen as nodal coordinates, but in contrast to standard ANCF elements, the proposed formulation is based on the two transversal slope vectors per node only. Different approaches for the virtual work of elastic forces are presented: a continuum mechanics based formulation, as well as a structural mechanics based formulation, which is in accordance with classical nonlinear beam finite elements. Since different interpolation functions as in standard ANCF elements are used, a much better convergence rate (up to order four) can be obtained. Therefore, the present element has high potential for application in geometrically nonlinear problems. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modelling of Wave Propagation using Spectral Finite Elements

For the analysis of wave propagation at high frequencies, the spectral finite element method (SFEM) is under investigation. In contrast to the conventional finite element method high-order shape functions are used. They are composed of Lagrange polynomials with nodes at the Gauß-Lobatto-Legendre points. The Gauß-Lobatto-Legendre integration scheme is applied in order to obtain a diagonal mass matrix. So, the resulting system equations can be solved efficiently. In the numerical examples, spectral finite elements with shape functions of different order are applied to a plane strain problem. The numerical examples cover structures without and with stiffness discontinuities. It is shown that the results agree well with analytical solutions. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A weak Galerkin mixed finite element method for the Helmholtz equation with large wave numbers

In this article, a new weak Galerkin mixed finite element method is introduced and analyzed for the Helmholtz equation with large wave numbers. The stability and well‐posedness of the method are established for any wave number k without mesh size constraint. Allowing the use of discontinuous approximating functions makes weak Galerkin mixed method highly flexible in term of little restrictions on approximations and meshes. In the weak Galerkin mixed finite element formulation, approximation functions can be piecewise polynomials with different degrees on different elements and meshes can consist elements with different shapes. Suboptimal order error estimates in both discrete H¹ and L² norms are established for the weak Galerkin mixed finite element solutions. Numerical examples are tested to support the theory.     

Mixed finite element methods for the Signorini problem with friction

In this article, we propose and study different mixed variational methods in order to approximate the Signorini problem with friction using finite elements. The discretized normal and tangential constraints at the contact interface are expressed by using either continuous piecewise linear or piecewise constant Lagrange multipliers in the saddle?point formulation. A priori error estimates are established and several numerical examples corresponding to the different choices of the discretized normal and tangential constraints are carried out. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Comparison of a mixed least-squares formulation using different approximation spaces

The main goal of the present work is the comparison of the performance of a least-squares mixed finite element formulation where the solution variables (displacements and stresses) are interpolated using different approximation spaces. Basis for the formulation is a weak form resulting from the minimization of a least-squares functional, compare e.g. [1]. As suitable functions for standard interpolation polynomials of Lagrangian type are chosen. For the conforming discretization of the Sobolev space vector-valued Raviart-Thomas interpolation functions, see also [2], are used. The resulting elements are named as P_mP_k and RT_mP_k. Here m (stresses) and k (displacements) denote the approximation order of the particular interpolation function. For the comparison we consider a two-dimensional cantilever beam under plain strain conditions and small strain assumptions. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A finite element formulation based on the theory of a Cosserat point

The theory of Cosserat points is the basis of a 3D finite element formulation allowing for large deformations in structural mechanics, that recently was presented by [1]. First attempts have revealed, that this formulation is free of showing undesired locking or hourglassing-phenomena. It additionally shows excellent behaviour for any type of incompressible material, for large deformations and sensitive structures such as plates or shells. Within the theory of Cosserat points, the position vectors X and x , are described through director vectors D _i and d _i by use of trilinear shape functions Nⁱ for an 8-node brick element. The special choice of shape functions Nⁱ allows for director vectors with which the deformation can be split into a homogeneous and an inhomogeneous part. This split enables the use of stiffnesses that correspond to different deformation modes. Analytical solutions to the inhomogeneous deformation modes are incorporated in the formulation and avoid the undesired phenomena. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical experiments on optimal shape parameters for radial basis functions

A numerical investigation on a technique for choosing an optimal shape parameter is proposed. Radial basis functions (RBFs) and their derivatives are used as interpolants in the asymmetric collocation radial basis method, for solving systems of partial differential equations. The shape parameter c in RBFs plays a major role in obtaining high quality solutions for boundary value problems. As c is a user defined value, inexperienced users may compromise the quality of the solution, often a problem of this meshless method. Here we propose a statistical technique to choose the shape parameter in radial basis functions. We use a cross‐validation technique suggested by Rippa 6 for interpolation problems to find a cost function Cost(c) that ideally has the same behavior as an error function. If that is the case, the parameter c that minimizes the cost function will be an optimal shape parameter, in the sense that it minimizes the error function. The form of the cost and error functions are analized for several examples, and for most cases the two functions have a similar behavior. The technique produced very accurate results, even with a small number of points and irregular grids. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

H-p clouds—an h-p meshless method

A new methodology to build discrete models of boundary-value problems is presented. The h-p cloud method is applicable to arbitrary domains and employs only a scattered set of nodes to build approximate solutions to BVPs. This new method uses radial basis functions of varying size of supports and with polynomial-reproducing properties of arbitrary order. The approximating properties of the h-p cloud functions are investigated in this article and several theorems concerning these properties are presented. Moving least squares interpolants are used to build a partition of unity on the domain of interest. These functions are then used to construct, at a very low cost, trial and test functions for Galerkin approximations. The method exhibits a very high rate of convergence and has a greater flexibility than traditional h-p finite element methods. Several numerical experiments in 1-D and 2-D are also presented. © 1996 John Wiley & Sons, Inc.     

Frequency-dependent vibration analysis of symmetric cross-ply laminated plate of Levy-type by spectral element and finite strip procedures

This research describes spectral finite element formulation for vibration analysis of rectangular symmetric cross-ply laminated composite plates of Levy-type based on classical lamination plate theory (CLPT). Formulation based on SFEM includes partial differential equations of motion, spectral displacement field, dynamic shape functions, and spectral element stiffness matrix (SESM). In this paper, vibration analysis of composite plate is investigated in two sections: free vibrations and forced vibrations. In free vibrations, natural frequencies are calculated for different Young’s moduli ratios and boundary conditions. In forced vibrations, plate vibrations are investigated under high-frequency concentrated impulsive loads. The resulting responses due to spectral element formulation are compared with those of (time-domain) finite element and analytical formulations, whenever available. The results demonstrate the superiority of SFEM with respect to FEM, in reducing computational burden, simultaneously increasing numerical accuracy, specifically for excitations of high-frequency content. By reducing the time duration of impulsive loads, and consequently increasing the modal contribution of higher modes in the transient response of plate, the accuracy of FEM responses decreases substantially accompanied with a high volume of computations, while the accuracy of the SFEM response results is very high and simultaneously, with a low computational burden. Practically, SFEM follows very closely exact analytical solutions.     

A smart displacement based (SDB) beam element with distributed plasticity

The standard displacement based inelastic beam element suffers of approximations related to the inability of the cubic polynomial interpolation functions to properly describe the displacement response of the beam when exhibiting inelastic behaviour. The increase of the number of finite elements, or the use of higher order functions with additional internal degrees of freedom, are common remedies suggested to improve the approximation leading to an unavoidable reduction of the computational efficiency. Alternatively, it has been shown that the development of force based finite elements, based on the adoption of exact force shape functions, lead to more accurate results, although requiring different and more complicated iterative solution strategies. Within this scenario, this paper proposes a new inelastic beam element, within the context of the displacement based approach, based on variable displacement shape functions, whose analytic expressions are related to the plastic deformation evolution in the beam element. The adaptive generalised displacement shape are obtained by identifying, at each step, an equivalent tangent beam, characterised by abrupt variations of flexural stiffness, as a suitable representation of the current inelastic state of the beam. The presented approach leads to the formulation of a Smart Displacement Based (SDB) beam element whose accuracy appears to be comparable to those obtained through a force based approach but requiring a reduced implementation effort and a more straightforward approach. The term ‘smart’ aims at emphasizing the ability of the element to upgrade the displacement field according to the current inelastic state.     

Exponential Finite Elements for a Phase Field Fracture Model

Sharp interface material models can be related to phase field models by introducing an order parameter, whose value is assigned to the different phases of a material. The elastic material law is coupled to the evolution equation of the order parameter and cracking is addressed as a phase transition problem instead of a moving boundary value problem. A regularization parameter ϵ controls the width of the diffuse cracks represented by the order parameter and the underlying sharp interface model can be recovered from the phase field model by the limit ϵ → 0. However, in numerical simulations using standard finite elements with linear shape functions, the minimum value of ϵ is restricted by the grid size and therefore the discretization of the crack field requires extensive mesh refinement for small values of ϵ. In this work, we construct special 2d shape functions which take into account the exponential character of the crack field and its dependence on the parameter ϵ. Especially in simulations with small values of ϵ and a rather coarse mesh, the elements with exponential shape functions perform significantly better than standard linear elements. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)