Revision as of 15:33, 22 October 2018

Abstract

In this paper an assumed strain approach is presented in order to improve the membrane behaviour of a thin shell triangular element. The so called Basic Shell Triangle (BST) has three nodes with only translational degrees of freedom and is based on a Total Lagrangian Formulation. As in the original BST element the curvatures are computed resorting to the surrounding elements (patch of four elements). Membrane strains are now also computed from the same patch of elements which leads to a non-conforming membrane behaviour. Despite this non-conformity the element passes the patch test. Large strain plasticity is considered using a logarithmic strain-stress pair. A plane stress behaviour with an additive decomposition of elastic and plastic strains is assumed. A hyperplastic law is considered for the elastic part while for the plastic part an anisotropic quadratic (Hill) yield function with non-linear isotropic hardening is adopted. The element, termed CBST, has been implemented in an explicit (hydro-)code adequate to simulate sheet-stamping processes and in an implicit static/dynamic code. Several examples are given showing the good performance of the enhnaced rotation-free shell triangle.

1 Introduction

Development of efficient and robust shell finite elements for modelling large-scale industrial applications has been a field of constant research in the last thirty years. The continuining efforts of many scientists has focused in the direction of deriving simple elements (preferably triangles) applicable to large scale computation, typical of practical shell problems. Some recent successful shell triangles are reported in [1,2]. For a comprehensive review of shell elements see [3]. From the theoretical point of view Kirchhoff hypothesis are enough for the simulation of the essential aspects of most shells problems, including structures in civil, mechanical and naval engineering, sheet metal forming, crash-worthiness analysis and others. Due the well-known difficulties associated to Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle C^{1}}

continuity, “thick shell” approaches that need Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle C^{0}}
 continuity only have been extensively used, although they are computationally  more expensive due to the number of degrees of freedom necessary to obtain thin   solutions with similar precision. The development of shell elements without  transversal shear strains then collides with the need for CFailed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle ^{1}}
continuity, which is not  simple to satisfy in general conditions and has led to the development of  non-conforming approaches. Among the many options of this kind, the use of shell elements with only  displacements as degrees of freedom (rotation-free) is very attractive from  the computational point of view.

Few shell elements exists in the literature with only translations as degrees of freedom. They are generally based on non-conforming approaches. A state of the art can be found in Ref. [4]. Recently Cirak y Ortiz ^[5] have developed a conforming thin shell element, but it departs in some aspects from the standard finite element formulation. The use of rotation-free shell elements in large scale simulations of sheet metal forming processes is growing ^[6]-[9], specially in explicit integration codes.

In Ref. [10] we presented a thin shell triangular element with only displacements as degrees of freedom, based on a total lagrangian formulation. The element is an extension of the rotation-free Basic Shell Triangle (BST) originally developed by Oñate y Zárate ^[4] using an updated lagrangian formulation The element has three nodes and a linear approximation of the displacement field. The membrane part is identical to the standard constant strain triangle. The bending part is also constant and is computed from an integration over the element boundary using information from the deflection gradients of the adjacent elements. This leads to a non-conforming element for the bending part, but it converges to the right solution and is robust. The BST element has a very good performance in bending for structured meshes, but it is less accurate for non-structured grids. On the other hand, for membrane dominated problems, as is the case for most sheet metal forming simulations, it requires fine discretizations associated to the constant strain triangle behaviour.

In this paper an extension of the rotation-free BST element is presented. The displacement gradient terms needed for the membrane and bending strains are computed from the nodal displacements of a patch of element using a new procedure. This allows to improve the membrane behaviour and obtain a smoother curvature field. The outline of this paper is as follows. In sections 2 and 3 we introduce the kinematics of the shell and the constitutive models used in the numerical simulations. Sections 4 to 6 describe the finite element approximation, including the new proposal for the gradient evaluation from a patch of elements, the computation of the metric tensor and the curvatures and the derivation of the necessary element expressions for the computer implementation. Sections 7 and 8 sumarize the numerical experiments performed in the linear and non-linear range. Numerical results demonstrate the excellent performance of the enhanced three node rotation-free shell triangle.

2 Shell kinematics

A summary of the most relevant hypothesis related to the kinematic behaviour of the shell are presented. Further details may be found in the wide literature dedicated to this field [3].

Consider a shell with undeformed middle surface occupying the domain Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Omega ^{0}}

in Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle R^{3}}
with a boundary Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Gamma ^{0}}

. At each point of the middle surface a thickness Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h^{0}}

is defined. The positions Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{x}^{0}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{x}}
of any point in the undeformed and the deformed configurations  can be respectively written as a function of the position of the middle  surface Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{\varphi }}
and the normal Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{t}_{3}}
at the point as

where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \xi _{1},\xi _{2}}

are curvilinear local coordinates defined over the  middle surface of the shell, and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \zeta }
is the distance in the undeformed  configuration of the point to the middle surface. The product Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \zeta \lambda }
 is the distance of the point to the middle surface measured on the deformed  configuration. This implies a constant strain in the normal direction  associated to the parameter Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \lambda }
relating the thickness at the present and initial   configurations, i.e.

A convective system is defined at each point as

This can be particularized for the points on the middle surface as

This allows to introduce the covariant (first fundamental form) and contravariant metric tensors of the middle surface as

and the curvatures (second fundamental form) of the middle surface as

The deformation gradient can be written as

The product Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{F}^{T}\mathbf{F=U}^{2}=\mathbf{C}}

(where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{U}}
is the right stretch tensor, and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{C}}
the right Cauchy-Green  deformation tensor) is

where the derivatives of the thickness ratio Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \lambda _{^{\prime }a}}

have been  neglected. Neglecting also the term associated to Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \zeta ^{2}}
and introducing  the definition of the covariant metric tensor Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle a_{\alpha \beta }}
at the middle  surface and the curvatures Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle k_{\alpha \beta }}
gives

Above expression shows that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{U}^{2}}

  is non-zero at the original  configuration for curved surfaces (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \kappa _{ij}^{0}\neq{0}}

). The changes of curvature of the middle surface are computed by

For computational convenience the following approximate expression (which is exact for initially flat surfaces) will be adopted

Above expression of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{U}^{2}}

is useful  to compute different Lagrangian strain  measures. An advantage of these measures is that they are associated to  material fibres, what makes it easy to take into account  material anisotropy. It is also useful to compute  the eigen decomposition of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{U}}
as

where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \lambda _{\alpha }}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{r}_{\alpha }}
are the eigenvalues and  eigenvectors of the right stretch tensor Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{U}}

.

In order to treat plasticity at finite strains an adequate stress-strain pair must be used. The Hencky measures will be adopted here. The (logarithmic) strains are defined as

Consistently, the Hencky stress tensor Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{T}}

 will be  used as the stress measure. Using a total lagrangian formulation it is convenient to work with  the the second Piola-Kirchhoff stresses (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{S}}

) for the evaluation of the residual forces. We define the rotated tensors as

where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{R}_{L}}

is the rotation tensor obtained from the eigenvectors  of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{U}}
given by

The relationship between the Hencky and Piola-Kirchhoff stresses is

The second Piola-Kirchhoff stress tensor can be finally computed by

The generalized stresses (forces and moments) are obtained by integrating across the original thickness the second Piola-Kirchhoff stress tensor using the actual distance to the middle surface for the evaluation of the bending moments, i.e.

With these values the weak form of the equilibrium equations can be written as

where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol K}}

is the virtual curvature vector and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{E}}
is the virtual  Green-Lagrange strain vector of the middle surface, with

3 Constitutive models

In this section, a brief description of the constitutive models used in the numerical examples presented in Sections 7 and 8 is given. Two types of materials are described: an elastic-plastic material associated to thin rolled metal sheets and a hyper-elastic material for rubbers.

In the case of metals, where the elastic strains are small, the use of a logarithmic strain measure reasonably allows to adopt an additive decomposition of elastic and plastic components as

A constant linear relationship between (plane) stresses and elastic strains is also adopted giving

These constitutive equations are integrated using a standard return algorithm. The following Mises-Hill^[11] yield function with non-linear isotropic hardening is adopted

where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle F, G, H}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N}
define the non-isotropic shape of the yield  surface and the parameters Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \kappa ,}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \varepsilon _{0}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle n}
define its size  as a function of the effective plastic strain Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle e^{p}}

.

The Mises-Hill yield function has the advantage of simplicity and it allows, as a first approximation, to treat rolled thin metal sheets with planar and transversal anisotropy.

For the case of rubbers, the Ogden^[12] model extended to the compressible range is considered. The material behaviour is characterized by the strain energy density per unit undeformed volume of the form

where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle K}

is the bulk modulus of the material, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle J}
is the determinant of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{U}}

, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mu _{i}}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha _{i}}
are material parameters, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mu _{i}\,,\,\alpha _{i}}
are real numbers such that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mu _{i}\alpha _{i}>0}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle  (\forall i=1,N)}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N}
is a positive integer.

The stress measures associated to the principal logarithmic strains are denoted by Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \beta _{i}} . They can be computed noting that

we define now

which gives

The values of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \beta _i} , expressed in the principal strains directions, allow to evaluate the stresses in the convective coordinate system as

We note that the Hencky stress tensor Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{T}}

can be easily  particularized for the plane stress case.

4 Interpolation functions and gradient evaluation

The starting point of the enhanced rotation-free of BST element is to discretize the shell surface with a standard 3-node triangular mesh. The difference with a standard finite element method is that for the computation of strains within an element, the configuration of the three adjacent triangular elements is used. Then at each triangle a patch of four elements, formed by the central triangle and the three adjacent ones is considered (see Figure 1.a). This allows to define a quadratic interpolation of the geometry from the position of the 6 nodes. In the isoparametric space, we keep the vertices of the 3-node master element (standard linear triangle) with the positions

• node 1: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \left( \xi ,\eta \right) =\left( 0,0\right) }

• node 2: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \left( \xi ,\eta \right) =\left( 1,0\right) }

• node 3: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \left( \xi ,\eta \right) =\left( 0,1\right) }

and the three additional nodes that complete the patch with the positions

• node 4: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \left( \xi ,\eta \right) =\left( 1,1\right) }

• node 5: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \left( \xi ,\eta \right) =\left( -1,1\right) }

• node 6: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \left( \xi ,\eta \right) =\left( 1,-1\right) }

Figure 1: Patch of four elements (a)in spacial coordinates (b)in natural coordinates

It is possible to define now a set of (non-standard) quadratic shape functions over over the six node element as:

This interpolation allows to compute the displacement gradients at selected points in order to use an assumed strain approach. The computation of the gradients is performed at three points over the boundary of the central element of the patch. These points are located at the mid-point of each side, denoted by G1, G2 and G3 in Figure 1.b. This election has the following characteristics.

• Gradients at these points depend only on the nodes pertaining to the two elements adjacent to the side. This can be easily verified by sampling the derivatives of the shape functions at each mid-side point.

• When gradients are computed at the common mid-side point from two adjacent elements, the same values are obtained, as the coordinates of the same four points are used. That is, the gradients at the mid-side points are independent of the element where they are computed. A side-oriented implementation of the finite element could lead to a unique evaluation of the gradients per side.

Next, some details of the implementation chosen are presented. We denote by Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{t}_{1}}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{t}_{2}}
the orthogonal unit vectors over the  tangent plane (undeformed configuration) associated to a conveniently selected  local cartesian system. The cartesian derivatives of the shape functions are computed at the original configuration by the standard expression

where the Jacobian matrix at the original configuration is

With the previous definitions the deformation gradient on the middle surface, associated to an arbitrary spatial cartesian system and to the material cartesian system defined on the middle surface, is

The covariant metric tensor can be computed from above expression as

The component of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol g}

can be used to compute any convenient membrane strain measures. For example the  Green-Lagrange strain tensor is obtained by

In the usual FEM matrix notation

The virtual strains are obtained by the variation of above expression as

In sides of elements where the adjacent element does not exist (i.e. at the boundaries of the shell), the gradient on that side is made equal to the gradient computed over the 3 nodes of the central element of the patch. From the definition of the gradients at the three mid-side points of the triangle, the element formulation can be interpreted as an assumed strain approach, where the metric tensor is interpolated over the element using the values computed on the sides by

5 Computation of curvatures

Curvatures are assumed to be constant over the element. They are computed by performing an average over the central element as

where (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A^{M}}

is the original area of the element).

Integrating by parts Eq.(43), the following integral over the element boundary Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Gamma ^{M}}

is  obtained:

The three distinct components of the curvature tensor are

The numerical evaluation of the boundary integral in Eq.(44) results in a sum over the integration points on the element boundary which are, in fact, the same points used for the computation of the gradients in the membrane formulations, i.e.

For the definition of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{t}_{3}} , the linear interpolation over the central element is used. In this case the tangent plane components are (with Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle a_{i}}

 and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle b_{i}}
being the usual cartesian projections of the sides when area coordinates  are used)

From these expressions it is also possible to compute in the original configuration the element area Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A^{M}} , the outer normals Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \left( n_{1} ,n_{2}\right) ^{i}}

at each side and its lengths Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle l_{i}}

. Eq.(47) also allows to evaluate the thickness ratio Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \lambda }

in the deformed configuration and the actual normal Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{t}_{3}}

. In Eq.(46), Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{N}^{I}}

are the linear shape functions of the three node linear triangle [3].

As one integration point is used over each side, it is not necessary to distinguish between sides (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle I} ) and integrations points (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle G} ) any more. In this way the curvatures can be computed by:

In the original BST element^[4,10] the gradient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{\varphi }_{\alpha }^{I}}

was computed as the average of the linear approximations over  the two adjacent elements. In the new version, the gradient is  computed at each side Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle I}
from the quadratic interpolation

Note again than at each side the gradients depend only on the positions of the three nodes of the central triangle and of an extra node (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle I+3} ), associated precisely to the side (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle I} ) where the gradient is computed.

An alternative form to express the curvatures, which is useful when their variations are needed, is to define:

This gives

This last expression allows to interpret the curvatures as the projections of the vectors Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{h}_{ij}}

over the normal of the central element. Direction  tFailed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle _{3}}
can be seen as a reference direction.  If a different  direction  is chosen, at an angle Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \theta }
with the former, this has an  influence of order Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \theta ^{2}}
in the projection. This justifies Eq.(47) for the  definition of tFailed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle _{3}}
as a function exclusively of the three nodes of  the central triangle,  instead of the 6-node  isoparametric interpolation.

5.1 Variation of the curvatures

From Eq.(52) the variation of the curvatures is obtained as

The first term is crucial in the evaluation of the curvatures variation. The variations of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{h}_{ij}}

are:

where

The variations of tFailed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle _{3}}

can be computed as shown in the  original BST element [10].

leading to

where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{\tilde{\varphi }}_{^{\prime }\alpha }}

are the  contravariant base vectors in the central triangle

We then have

Substituting the last expression in Eq.(53), results

where the projections of the vectors Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{h}_{ij}}

over the  contravariant base vectors Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{\tilde{\varphi }}_{^{\prime }\alpha }}
have  been included

These projections together with Eq.(52) allows to write

5.2 Boundary conditions

Elements at the domain boundary, where an adjacent element does not exist, deserve a special attention. The treatment of essential boundary conditions associated to translational constraints is straightforward, as they are the natural degrees of freedom of the element. The conditions associated to the normal vector are crucial in this formulation. For clamped sides or symmetry planes, the normal vector Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{t}_{3}}

must be kept fixed (clamped case), or  constrained to move in the plane of symmetry (symmetry case). The former case can be seen as a  special case of the latter, so we will consider symmetry planes only. This  restriction can be imposed through the definition of the tangent plane at the  boundary, including the normal to the plane of symmetry Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{\varphi  }_{^{\prime }n}^{\left( 0\right) }}
that does not change during the process.

Figure 2: Local cartesian system for the treatment of boundary conditions

The tangent plane at the boundary is expressed in terms of two orthogonal unit vectors referred to as local-to-the-boundary cartesian system (see Figure 2) and defined as

where vector Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{\varphi }_{^{\prime }n}^{\left( 0\right) }}

is fixed while  direction Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{\bar{\varphi }}_{^{\prime }s}}
emerges as the intersection of  the symmetry plane with the plane defined by the central element (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M}

). The plane (gradient) defined by the central element is

the intersection of this plane with the plane of symmetry can be written in terms of the position of the nodes that define the side (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle J}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle K}

) and the original length of the side Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle L_{0}} , i.e.

That together with the outer normal to the side Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{n =}\left[ \nu _{1},\nu _{2}\right]^T } (resolved in the selected original convective cartesian system) leads to

where the normal component of the gradient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{\varphi }_{^{\prime }n}}

is

In this way the contribution of the gradient at side Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle I}

to vectors Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{h}_{ij}}
results in

The only difference with Eq.(54) necessary for the computation of the curvature variations, is that now the contribution from the gradient at side Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle I}

is

where the influence of variations in the length of vector Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{\varphi }_{^{\prime }n}}

has been neglected.

In the case of natural boundary conditions (either free or simple supported) the gradient at the sides is supposed to be equal to the gradient in the central element, similarly as in the membrane approach, i.e.

Using Eq.(71), vectors Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{h}_{ij}}

and their variations can be easily  computed. A second possibility is to make use of the natural boundary  condition constraining the normal curvature to a zero  value. In simple supported  sides where the curvature along the side is zero, this leads to zero values for  both bending moments.

6 Stiffness matrix evaluation

When a predictor-corrector scheme is used to trace the movement of the shell, the derivative (stiffness matrix) of the weak form of the momentum equations (23) is needed. As usual, material and geometric components are computed separately. The material part does not offer difficulties and it is computed from the integral

where matrix Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{B}={\boldsymbol B}_m+{\boldsymbol B}_\phi }

 includes:

a) a membrane part Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol B}_m}

computed at each mid side point Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle I}
from

Typically one integration point is used for computing the terms in Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol B}_m} .

b) a bending part Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol B}_\phi }

which results from Eq.(60) that is  constant over the element.

In Eq.(72), C is the elasticity matrix integrated in the thickness. If material non-linearity is considered, a layer wise numerical integration across the thickness is required This leads to the tangent or algorithmic elastic-plastic constitutive matrix Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{C}_{ep}} .

6.1 Geometric stiffness (membrane part)

The geometric stiffness due to membrane forces results from computing

This can be written as the sum of the contributions of the three sides, i.e.

6.2 Geometric stiffness (bending part)

The geometric stiffness associated to bending moments is more involved. Its expression stems from

Recalling expressions (51) and (52) it can be written

where

The first two terms lead to symmetric components. The second one Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \left( \delta \mathbf{t}_{3}\cdot \Delta \mathbf{h}_{ij}\right) }

 can be expressed as

Then, a first component of the geometric bending stiffness can be written as

The last term in (78), can be obtained observing that (sum in Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha }

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \beta }

)

then

with

where the covariant metric tensor at the middle surface Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle a^{\alpha \beta }}

has  been used.

Denoting the sum

the term

results in

The last expression (87) has components only in the nodes of the central element, which stems from the definition used for tFailed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle _{3}} .

Numerical experiments have shown that the bending part of the geometric stiffness is not so important and can be disregarded in the iterative process.

7 Numerical examples in the linear range

In this section, a summary of the results obtained in the assessment of the proposed element in the linear range are presented. In this first part of the numerical experiments, results have been obtained with a static/dynamic code that uses a implicit solution of the discretized equilibrium equations. In the examples, the original BST element^[10] is compared with the enhanced version here proposed, denoted by CBST when a numerical integration is performed with three points, and CBST1 when only one integration point is used. Note that one integration quadrature is equivalent to averaging the metric tensors computed at each side.

7.1 Membrane patch test

Figure 3: Membrane patch test

The main objective of the present formulation is to obtain a (non conforming) membrane element with a similar performance than the linear strain triangle, that satisfies the patch test. To assess this, a square domain subjected to nodal forces associated to a constant stress state (in both directions and in shear) is considered. Figure 3 shows the patch of elements and the necessary nodal forces to obtain a uniform tensile stress in direction Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x_{1} } . Note that the nodal forces used correspond to the constant strain triangle, not to the linear one. For the distorted mesh shown, correct results are obtained using either 1 or 3 integration points.

7.2 Bending patch test (torsion)

Figure 4: Patch test for uniform torsion

The element bending formulation does not allow to apply external bending moments (there are not rotational DOFs). Hence it is not possible to analyse a patch of elements under loads leading to a uniform bending moment. A uniform torsion can be considered if a point load is applied at the corner of a rectangular plate with two consecutive free sides and two simple supported sides. Figure 4 shows three patches leading to correct results both in displacements and stresses. All three patches are structured meshes. When the central node in the third patch is displaced from the center, the results obtained with the CBST element are not correct. The original element BST gives correct results in all cases, if natural boundary conditions are imposed in the formulation. If this is not the case, incorrent results are obtained even with structured meshes.

7.3 Cook's membrane problem

This example is used to assess the membrane performance of the CBST element and to compare it with the linear triangle (constant strain) and the quadratic triangle (linear strain). This example involves important shear energy and was proposed to assess the distortion capability of elements. Figure 5.a shows the geometry and the applied load. Figure 5.b plots the vertical displacement of the upper vertex as a function of the number of nodes in the mesh. In this figure results obtained with other isoparametric elements have been also plotted for comparison. They include the constain strain triangle (CST), the bilinear quadrilateral (QUAD4) and the linear strain triangle (LST). Note that in this membrane problem both the BST and the CST elements give identical results.

Figure 5: Cook membrane problem (a) Geometry (b) Results

From the plot shown it can be seen that the new element with three integration points (CBST) gives values slightly better that the constant strain triangle for the most coarse mesh (only two elements). However, when the mesh is refined, a performance similar to the linear strain triangle is obtained, that is dramatically superior that the former. On the other hand, if a one point quadrature is used (CBST1) the convergence in the reported displacement is notably better that for the rest of the elements.

7.4 Cylindrical roof

In this example an effective membrane interpolation is of primary importance. Hence this is good test to assess the new element. The geometry is a cylindrical roof supported by a rigid diaphragm at both ends and it is loaded by a uniform dead weight (see 6.a.). Only one quarter of the structure is modelled due to symmetry conditions. Unstructured and structured meshes are considered. In the latter case two orientations are possible (Figure 6 shows orientation B).

Tables 1, 2 and 3 present the normalized vertical displacements at the crown (point A) and in the midpoint of the free side (point B) for the two orientations of structured meshes and non-structured meshes. Values used for normalization are Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle u_{A}=0.5407}

y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle u_{B}=-3.610}
that are quoted in reference [13].

Figure 6: Cylindrical roof under dead weight. E=3 EFailed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): ^{6} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \nu=0.0 , Thickness =3.0, shell weight =0.625 per unit area.

Figure 6.b shows the normalized displacement of point-B over the structured meshes as a function of the number of degrees of freedom. An excelent convergence the CBST element can be seen. The version with only one integration point (CBST1) presents a behavior a little more flexible and converges from above. Table 3 shows that for non-structured meshes the element converges to the correct value but more slowly.

7.5 Open semi-espherical dome with point loads

The main problem of finite elements with initially curved geometry is the so called membrane locking. The CBST element proposed has a quadratic interpolation of the geometry, then it may suffer from this problem. To assess this we resort to an example of inextensional bending. This is an hemispherical shell of radius Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle r=10}

and  thickness Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h=0.04}
with an 18Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle ^{o}}
hole in the pole and free at all  boundaries,  subjected to two inward and two outward forces 90Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle ^{o}}
apart.  Material properties are Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle E=6.825\times{10}^{7}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \nu=0.3}

. Figure 7.a shows the discretized geometry (only one quarter of the geometry is considered due to symmetry).

Figure 7: Pinched hemispherical shell with a hole, (a)geometry, (b)normalized displacement

In Figure 7.b the displacements of the points under the loads have been plotted versus the number of nodes used in the discretization. Due to the orientation of the meshes chosen, the displacement of the point under the inward load is not the same as the displacement under the outward load, so in the figure an average (the absolute values) has been used. Results obtained with other elements have been included for comparison: two membrane locking free elements namely the original linear BST and a transverse shear-deformable quadrilateral (QUAD); a transverse shear deformable quadratic triangle (TRIA)^[1] that is vulnerable to locking and an assumed strain version of this triangle (TRIA-1)^[1] that alleviates membrane locking but does not eliminate it.

From the plotted results it can be seen that the CBST element presents membrane locking in bending dominated problems with initially curved geometries. This locking is much less severe than in a standard quadratic triangle and even than for the assumed strain version used for comparison. When only one integration point is used (CBST1 element) membrane locking disappears.

8 Non linear numerical examples

Results for examples with geometric and material non-linearities are presented next using the CBST1 element. Due to the features of the modelized problems, with strong nonlinearities associated to instabilities and frictional contact conditions, a code with explicit integration of the dynamic equilibrium equations has been used^[14]. This code allows to obtain pseudo-static solutions throught dynamic relaxation.

8.1 Inflation of a sphere

The example is the inflation of a spherical shell under internal pressure. An incompressible Mooney-Rivlin constitutive material have been considered. The Ogden parameters are Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N=2} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha _{1}=2} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mu _{1}=40} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha _{2}=-2} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mu _{2}=-20} . Due to the simple geometry an analytical solution exists^[15] (with Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \gamma =R/R^{\left( 0\right) }} ):

In this numerical simulation the same geometric and material parameters used in Ref. [16] have been adopted: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle R^{\left( 0\right) }=1}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h^{\left( 0\right) }=0.02}

. The three meshes considered to evaluate convergence are shown in Figure 8.a. The value of the actual radius as a function of the internal pressure is plotted in Figure 8.b for the different meshes and is also compared with the analytical solution. It can be seen that with a few degrees of freedom it is possible to obtain an excellent agreement for the range of strains considered. The final value corresponds to a thickness radius ratio of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h/R=0.00024} .

Figure 8: Inflation of sphere of Mooney-Rivlin material. (a) Meshes used in the analysis (b) Change of radius as a function of the internal pressure.

8.2 Inflation of an air-bag

This example has also been taken from Ref.[16] where it is shown that the final configuration is mesh dependent due to the strong instabilities leading to a non-uniqueness of the solution. In [16] it is also discussed the important regularizing properties of the bending energy, that when disregarded leads to massive wrinkling in the compressed zones.

The air bag geometry is initially square with an undeformed diagonal of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 1.20} . The constitutive material is a linear isotropic elastic one with modulus of elasticity Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle E=5.88\times{10}^{8}}

and Poisson's ratio Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \nu=0.4}

. Only one quarter of the geometry has been modelled due to symmetry. Only the normal to the original plane is constrained along the boundaries. The thickness considered is Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h=0.001}

and the inflation pressure is Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 5000}

. Using a density Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \delta=1000} , pressure is linearly incremented from 0 to the final value in Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle t=0.1} .

With comparative purposes and also to backup the comments in Ref. [16] two analyses have been performed, a purely membranal one and one including bending effects. Figure 9 shows the final deformed configurations for three meshes with 289, 1089 and 4225 nodes. The top row corresponds to a full analysis including bending and the central row is a pure membrane analysis. The bottom row is also an analysis including bending where the mesh orientation has been changed.

Figure 9: Inflation of a square air-bag . Deformed configurations for three different meshes with 800, 3136 and 12416 degrees of freedom.

The top and bottom lines show the final shapes change according to the degree of discretization and mesh orientation due to instabilities and non uniqueness of the solution. The central row shows the noteworthy increment of wrinkles as the mesh is refined. Finally in Figure 10 presents the convergence of the maximum displacement of the central point as the mesh is refined.

Figure 10: Inflation of a square air-bag. Convergence of the maximum displacement versus the number of degrees of freedom in the model.

8.3 Deep drawing of a rolled sheet

The element performance in problems involving large strains and anisotropic plastic behaviour is assesed in a benchmark proposed in a recent NUMISHEET^[17]> meeting. This is the deep drawing of a circular mild steel sheet with a spherical punch of radius Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 50 mm}

(Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle R_{p}}

). The original radius of the sheet (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle R_{s}} ) is Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 100 mm}

 (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \frac{R_{L}}{R_{P}}=\beta=2}

). The drawing depth is (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle DD} ) is 85 mm (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \frac{DD}{R_{P}}=1.7} ) and the force over the blank holder is Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 80 kN} . The constitutive material of the rolled sheet is characterized by three tensile tests along three different directions respect to rolling direction (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle RD} ) (Table 4). Note that the modellization requires to treat the material as transversally anisotropic, as this improves the deep drawability and if not considered leads to an early localized necking of the material. The yield function originally proposed by Hill^[11] has been chosen for both associative and non-associative potential functions. For elastic-plastic problems a numerical integration along the thickness direction is necessary. This is accomplished here using four integration points.

Only one quarter of the geometry has been considered due to symmetry. Tools are modelized as rigid surfaces, the punch has been defined by 1439 points and 2730 triangles. For the die 744 points and 490 quadrilaterals have been used. Finally the blank holder has been discretized with 155 points and 120 quadrilaterals. Figure 11.a shows a perspective of the tools and Figure 11.a the final deformed sheet. The sheet has been modelled with 6370 3-node triangles and 3284 nodes (9274 DOFs).

Figure 11: Deep drawing of a circular sheet. (a)geometry of the tools, (b)final deformed shape of the sheet.

The reported results are associated to three different meridiands: A at 90Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle ^{0}}

of the rolling direction, B: at 45Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle ^{0}}
of the rolling  direction and C: in the rolling direction. The numerical results are  compared with a set of experimental values reported by Thyssen Krupp Stahl AG  (who proposed the benchmark and supplied the sheet samples)[17]. It must be noted  that there was a great dispersion in the experimental data send to NUMISHEET,  so it is difficult to extract conclusions from a unique comparison. In any case the  experimental data supplied  seems to be consistent enought to have an  idea if the numerical model gives reasonable results.

The first element of comparison associated to the planar anisotropy, are the different displacements (draw in) of the points in the boundary of the sheet. Table 5 shows the values obtained for the three reference meridians chosen.

The values for the non-associative model are much like the experimental ones, those obtained with the associative model are a little lower. The largest draw-in is at the meridian at 45Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle ^{o}}

of rolling direction ( B). In the numerial simulation the draw-in at meridian  C is larger than in meridian  A, while in the experimental results   the latter are larger than the former. It should be reminded that other  numerical results presented at NUMISHEET and most of the numerical simulations  showed the same tendency of present values.

The experimental data supplied are the three logarithmic principal strains associated to the three meridians defined above. Here one strain at each meridian has been choosen for comparison. Figure 12.a shows the meridian strain (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle E_{1}} ) along meridian A (direction transversal to rolling); Figure 12.b shows the circunferencial strain (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle E_{2}} ) along meridian B and in Figure 12.c the thickness strain (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle E_{3}} ) along meridian C (rolling direction) has been plotted.

The experimental data show a saw-tooth profile (specially for the thickness strain) that is difficult to accept. These values are therefore not reliable point-wise but as a whole. Based on this fact and on the above mentioned differences between experimental values, it can be said that the present simulations agree quite well with the experimental values, specially when the non-associative plasticity model was used.

In the central part of the sheet the experimental data gives strains lower than the numerical simulation. It can be said that experimentally the sheet has better drawability. In the external zone there are less differences for in plane strains but for thickness strains they are larger. Experimentally thickness strains are more or less uniform between Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle -0.1}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle -0.2}
what  differs from present numerical simulations and also from other numerical  simulations and experimental data presented at NUMISHEET [17].

Finally Figure 12 shows resistance to punch as a funtion of punch travel. In the final part simulation forces are smaller than experimental values. This may be due to an incorrect definition of friction with the tools.

Figure 12: Results from NUMISHEET benchmark. (a) Meridional strain Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): E_{1} along transversal direction to rolling (meridian A); (b) hoop strain Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): E_{2} along direction at 45Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): ^{0} of rolling (meridian B); (c) Thickness strain Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): E_{3} along rolling direction (meridian C); (d) Punch force versus Punch travel;

From the present results it is difficult to draw definite conclusions about what should be improved in the numerical model due to all the factors involved, including the constitutive model, the thin shell theory, the contact treatment, the reliability of the experimental data of both the material characterization and the deep drawing, etc..

9 Conclusions

A membrane and bending non-conforming rotation-free shell finite element has been presented. The element pass the membrane patch test and the numerical experiments performed do not show problems. From the bending point of view the element is a little stiffer than the original BST element, but presents a smoother continuity of displacement gradients. Convergence rate in membrane dominated problems is similar to the linear strain triangle. Using three integration points the elements is vulnerable to membrane locking for initially curved surfaces. This locking effect dissapears when only one integration point is used, which the usual case for elastic-plastic problems and for any large scale simulation.

Acknowledgments

The first author is a member of the scientific staff of the Science Research Council of Argentina (CONICET). The support provided by grants of CONICET and Agencia Córdoba Ciencia S.E. is gratefully acknowledged.

BIBLIOGRAPHY

[1] F.G. Flores, E. Oñate and F. Zárate “New assumed strain triangles for non-linear shell analysis”, Computational Mechanics, 17, 107-114, 1995.

[2] J.H. Argyris, M. Papadrakakis, C. Aportolopoulun and S. Koutsourelakis. The TRIC element. Theoretical and numerical investigation. Comput. Meth. Appl. Mech. Engrg., 182, 217–245, 2000.

[3] O.C. Zienkiewicz and R.L. Taylor. The finite element method. Solid Mechanics. Vol II, Butterworth-Heinemann, 2000.

[4] E. Oñate and F. Zárate, “Rotation-free plate and shell triangles”, Int. J. Num. Meth. Engng, 47, pp. 557–603, 2000.

[5] F. Cirak and M. Ortiz, Subdivision surfaces: A new paradigm for thin-shell finite element analysis, Int. J. Num. Meths in Engng, vol. 47, 2000, págs. 2039-2072.

[6] D.Y. Yang, D.W. Jung, L.S. Song, D.J. Yoo and J.H. Lee, “Comparative investigation into implicit, explicit and iterative implic/explicit schemes for simulation of sheet metal forming processes” NUMISHEET'93, A. Makinouchi, E. Nakamachi, E. Oñate and R.H. Wagoner (Eds.), RIKEN, pp. 35–42, Tokyo, 1993.

[7] M. Brunet and F. Sabourin, “Prediction of necking and wrinkles with a simplified shell element in sheet forming”, Int. Conf. of Metal Forming Simulation in Industry, Vol. II, pp. 27–48, B. Kröplin (Ed.), 1994.

[8] G. Rio, B. Tathi and H. Laurent, “A new efficient finite element model of shell with only three degrees of freedom per node. Applications to industrial deep drawing test”, in Recent Developments in Sheet Metal Forming Technoloy, Ed. M.J.M. Barata Marques, 18th IDDRG Biennial Congress, Lisbon, 1994.

[9] Rojek J., Oñate E. and Postek E ., Application of explicit FE codes to simulation of sheet and bulk metal forming processes Journal of the Materials Processing Thecnology, 41, págs. 620-627, 1998.

[10] F.G. Flores F.G. and E. Oñate A basic thin shell triangle with only translational DOFs for large strain plasticity, Int. J. Num. Meths in Engng, vol. 51, pp 57-83, 2001.

[11] R. Hill, “A Theory of the Yielding and Plastic Flow of Anisotropic Metals”, Proc. Royal Society London, vol. A193, pp. 281, 1948.

[12] R.W. Ogden “Large deformation isotropic elasticity: on the correlation of theory and experiments for incompressible rubberlike solids”, Proceedings of the Royal Society of London, vol. A326, pp. 565-584, 1972.

[13] H.C. Huang, Static and Dynamic Analysis of Plates and Shells, page 40, Springer-Verlag, Berlin, 1989.

[14] STAMPACK, A General Finite Element System for Sheet Stamping and Forming Problems, Quantech ATZ, Barcelona, Spain, 2003 (www.quantech.es).

[15] A. Needleman, Inflation of spherical rubber ballons, Int. J. of Solids and Structures, 13, 409–421, 1977.

[16] F. Cirak and M. Ortiz, Fully Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle C^{1}} -conforming subdivision elements for finite deformations thin-shell analysis , Int. J. Num. Meths in Engng, vol. 51, 2001, 813-833.

[17] NUMISHEET'99, Fourth International Conference and Workshop on Numerical Simulation of 3D Sheet Forming Processes, NUMISHEET'99, J.C. Gelin and P. Picart (Eds.), BURS Edition, Besancon, France, 1999.