Published in Comput. Meth. Appl. Mech. Engng., Vol. 194 (21-24), pp. 2406-2443, 2005
doi: 10.1016/j.cma.2004.07.039
ABSTRACT
A family of rotation-free three node triangular shell elements is presented. The simplest element of the family is based on an assumed constant curvature field expressed in terms of the nodal deflections of a patch of four elements and a constant membrane field computed from the standard linear interpolation of the displacements within each triangle. An enhanced version of the element is obtained by using a quadratic interpolation of the geometry in terms of the six patch nodes. This allows to compute an assumed linear membrane strain field which improves the in-plane behaviour of the original element. A simple and economic version of the element using a single integration point is presented. The efficiency of the different rotation-free shell triangles is demonstrated in many examples of application including linear and non linear analysis of shells under static and dynamic loads, the inflation and de-inflation of membranes and a sheet stamping problem.
1 INTRODUCTION
Triangular shell elements are very useful for the solution of large scale shell problems such as those occurring in many practical engineering situations. Typical examples are the analysis of shell roofs under static and dynamic loads, sheet stamping processes, vehicle dynamics and crash-worthiness situations. Many of these problems involve high geometrical and material non linearities and time changing frictional contact conditions. These difficulties are usually increased by the need of discretizing complex geometrical shapes. Here the use of shell triangles and non-structured meshes becomes a critical necessity. Despite recent advances in the field [1]–[6] there are not so many simple shell triangles which are capable of accurately modelling the deformation of a shell structure under arbitrary loading conditions.
A promising line to derive simple shell triangles is to use the nodal displacements as the only unknowns for describing the shell kinematics. This idea goes back to the original attempts to solve thin plate bending problems using finite difference schemes with the deflection as the only nodal variable [7]–[9].
In past years some authors have derived a number of thin plate and shell triangular elements free of rotational degrees of freedom (d.o.f.) based on Kirchhoff's theory [10]–[26]. In essence all methods attempt to express the curvatures field over an element in terms of the displacements of a collection of nodes belonging to a patch of adjacent elements. Oñate and Cervera [14] proposed a general procedure of this kind combining finite element and finite volume concepts for deriving thin plate triangles and quadrilaterals with the deflection as the only nodal variable and presented a simple and competitive rotation-free three d.o.f. triangular element termed BPT (for Basic Plate Triangle). These ideas were extended and formalized in [20] to derive a number of rotation-free thin plate and shell triangles. The basic ingredients of the method are a mixed Hu-Washizu formulation, a standard discretization into three node triangles, a linear finite element interpolation of the displacement field within each triangle and a finite volume type approach for computing constant curvature and bending moment fields within appropriate non-overlapping control domains. The so called cell-centered and cell-vertex triangular domains yield different families of rotation-free plate and shell triangles. Both the BPT plate element and its extension to shell analysis (termed BST for Basic Shell Triangle) can be derived from the cell-centered formulation. Here the control domain is an individual triangle. The constant curvatures field within a triangle is computed in terms of the displacements of the six nodes belonging to the four elements patch formed by the chosen triangle and the three adjacent triangles. The cell-vertex approach yields a different family of rotation-free plate and shell triangles. Details of the derivation of both rotation-free triangular shell element families can be found in [20].
An extension of the BST element to the non linear analysis of shells was implemented in an explicit dynamic code by Oñate et al. [25] using an updated Lagrangian formulation and a hypo-elastic constitutive model. Excellent numerical results were obtained for non linear dynamics of shells involving frictional contact situations and sheet stamping problems [17,18,19,25].
A large strain formulation for the BST element using a total Lagrangian description was presented by Flores and Oñate [23]. A recent extension of this formulation is based on a quadratic interpolation of the geometry of the patch formed by the BST element and the three adjacent triangles [26]. This yields a linear displacement gradient field over the element from which linear membrane strains and constant curvatures can be computed within the BST element.
In this paper the formulation of the BST element is revisited using an assumed strain approach. The content of the paper is the following. First some basic concepts of the formulation of the original BST element using an assumed constant curvature field are given. Next, the basic equations of the non linear thin shell theory chosen based on a total Lagrangian description are presented. Then the non linear formulation of the BST element is presented. This is based on an assumed constant membrane field derived from the linear displacement interpolation and an assumed constant curvature field expressed in terms of the displacements of the nodes of the four element patch using a finite volume type approach. An enhanced version of the BST element is derived using an assumed linear field for the membrane strains and an assumed constant curvature field. Both assumed fields are obtained from the quadratic interpolation of the patch geometry following the ideas presented in [26]. Details of the derivation of the tangent stiffness matrix needed for a quasi-static implicit solution are given for both the BST and EBST elements. An efficient version of the EBST element using one single quadrature point for integration of the tangent matrix is presented. An explicit scheme adequate for dynamic analysis is briefly described.
The efficiency and accuracy of the standard and enhanced versions of the BST element is validated in a number of examples of application including linear and non linear analysis of shells under static and dynamic loads, the inflation and de-inflation of membranes and a sheet stamping problem.
2 FORMULATION OF THE BASIC PLATE TRIANGLE USING AN ASSUMED CONSTANT CURVATURE FIELD
Let us consider a patch of four plate three node triangles (Figure 1). The nodes 1, 2, and 3 in the main central triangle (M) are marked with circles while the external nodes in the patch (nodes 4, 5 and 6) are marked with squares. Mid side points in the central triangle are also marked with smaller squares. Kirchhoff's thin plate theory will be assumed to hold. The deflection is linearly interpolated within each three node triangle in the standard finite element manner as
|
(1) |
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 L_{i}^{e}}
are the linear shape functions (area coordinates) of the three node triangle, 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 w_{i}^{e}}
are nodal deflections and superindex 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}
denotes element values.
|
| Figure 1: Patch of three node triangular elements including the central triangle (M) and three adjacent triangles (1, 2 and 3) |
The curvatures within the central triangle can be expressed in terms of a constant assumed curvatures field as
|
(2) |
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 \kappa }}
is the 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 \hat{\boldsymbol \kappa }}
is the assumed constant curvature field defined as
|
(3) |
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 area of the central triangle in Figure 1.
Integrating by parts Eq.(3) and substituting the resulting expression for 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 \hat{\boldsymbol \kappa }}
into Eq.(2) gives the constant curvature field within the element as
|
(4) |
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 \Gamma _{M}}
is the boundary of the central triangle 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{n}=\left(n_{x},n_{y} \right)}
is the boundary normal. Equation (4) defines the assumed constant curvature field within the central triangle in terms of the deflection gradient along the sides of the triangle. Equation (4) can be found to be equivalent to the standard conservation laws used in finite volume procedures as described in [27,28].
The computation of the line integral in Eq.(4) poses a difficulty as the deflection gradient is discontinuous along the element sides. A simple method to overcome this problem is to compute the deflection gradient at the element sides as the average values of the gradient contributed by the two triangles sharing the side [20,28]. Following this idea the constant curvature field with the element is computed as
|
(5) |
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{w}^{p}=[w_{1},w_{2},w_{3},w_{4},w_{5},w_{6}]^{T}}
is the deflection vector of the six nodes in the patch. In Eq.(5) the sum extends over the three sides of 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}
, 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_{j}^{M}}
are the lengths of the element sides and superindexes 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}
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 j}
refer to the central triangle and to each of the adjacent elements, respectively. The standard sum convention for repeated indexes is used. Note that triangular area coordinates satisfy
|
(6) |
Note also that the constant curvature field is expressed in terms of the six nodes of the four element patch linked to the 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 expression of the 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 3\times{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 \mathbf{B}_{b}}
matrix can be found in [14,20].
The virtual work expression is written as
|
(7) |
where m is the bending moment field related to the curvatures by the standard constitutive equations
|
(8) |
In Eqs.(7) and (8) 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}
is the plate 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 E}
is the Young's modulus, 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 }
is the 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 \delta{\boldsymbol \kappa }}
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 \delta w}
are the virtual curvatures and the virtual deflection, respectively, 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 q}
is a distributed vertical load.
Substituting the approximation for the vertical deflection and the assumed constant curvature field into (7) leads to the standard linear system of equations
|
(9) |
where the stiffness 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{K}}
and the equivalent nodal force 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}}
can be found by assembly of the element contributions given by
|
(10) |
|
(11) |
Note 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{K}^{e}}
is a 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 6\times{6}}
matrix, whereas 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}^{e}}
has the same structure than for the standard linear triangle. The explicit form 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{K}^{e}}
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{f}^{e}}
can be found in [14].
The resulting Basic Plate Triangle (BPT) has one degree of freedom per node and a wider bandwidth than the standard three node triangles as each triangular element is linked to its three neighbours through Eq.(5).
Examples of the good performance of the BPT element for analysis of thin plates can be found in [14,20]. The extension of the BPT element to the analysis of shells yields the Basic Shell Triangle (BST) [20]. Different applications of the BST element to linear and non linear analysis of shells are reported in [14,17–20,23,25,26].
The ideas used to derive the BPT element will now be extended to derive two families of Basic Shell Triangles using a total Lagrangian description.
3 BASIC THIN SHELL EQUATIONS USING A TOTAL LAGRANGIAN FORMULATION
3.1 Shell kinematics
A summary of the most relevant hypothesis related to the kinematic behaviour of a thin shell are presented. Further details may be found in the wide literature dedicated to this field [8,9].
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 a point in the undeformed and the deformed configurations can be respectively written as a function of the coordinates 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 {\boldsymbol \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 arc-length curvilinear principal 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 from the point to the middle surface in the undeformed configuration. 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 from the point to the middle surface measured on the deformed configuration. 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 }
relates the thickness at the present and initial configurations as:
|
(14) |
This approach implies a constant strain in the normal direction. 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 }
will not be considered as an independent variable and will be computed from purely geometrical considerations (isochoric behaviour) via a staggered iterative update. Besides this, the usual plane stress condition of thin shell theory will be adopted.
A convective system is computed at each point as
|
(15) |
|
This can be particularized for the points on the middle surface as
|
The covariant (first fundamental form) metric tensor of the middle surface is
|
(20) |
The Green-Lagrange strain vector of the middle surface points (membrane strains) is defined as
|
(21) |
with
|
(22) |
The curvatures (second fundamental form) of the middle surface are obtained by
|
(23) |
The deformation gradient tensor is
|
(24) |
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) can be written as
|
(25) |
In the derivation of expression (25) 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}}
and the terms 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}}
have been neglected.
Equation (25) 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 not a unit tensor 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
|
(26) |
Note 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 \delta \chi _{ij}=\delta \kappa _{ij}} .
For computational convenience the following approximate expression (which is exact for initially flat surfaces) will be adopted
|
(27) |
This expression 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
|
(28) |
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 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}}
.
The resultant stresses (axial forces and moments) are obtained by integrating across the original thickness the second Piola-Kirchhoff stress 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 { \boldsymbol \sigma }}
using the actual distance to the middle surface for evaluating the bending moments. This gives
|
(29) |
|
(30) |
With these values the virtual work can be written as
|
(31) |
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 \delta \mathbf{u}}
are virtual displacements, 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{\boldsymbol \varepsilon }_{m}}
is the virtual Green-Lagrange membrane strain 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 \delta{\boldsymbol \kappa }}
are the virtual curvatures 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}}
are the surface loads. Other load types can be easily included into (31).
3.2 Constitutive models
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
|
(32) |
Two types of material models are considered here: 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
|
(33) |
A constant linear relationship between the (plane) Hencky stresses and the logarithmic elastic strains is adopted giving
|
(34) |
These constitutive equations are integrated using a standard return algorithm. The following Mises-Hill [29] yield function with non-linear isotropic hardening is chosen
|
(35) |
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 \sigma _{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 e_{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 simple Mises-Hill yield function allows, as a first approximation, to treat rolled thin metal sheets with planar and transversal anisotropy.
For the case of rubbers, the Ogden [30] model extended to the compressible range is considered. The material behaviour is characterized by the strain energy density per unit undeformed volume defined as
|
(36) |
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
|
(37) |
we define now
|
(38) |
which gives
|
(39) |
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 Hencky stresses in the convective coordinate system as
|
(40) |
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.
We define the rotated Hencky and second Piola-Kirchhoff stress 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
|
(43) |
The relationship between the rotated Hencky and Piola-Kirchhoff stresses 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 \left(\alpha , \beta=1,2 \right)}
|
(44) |
The second Piola-Kirchhoff stress tensor can be computed by
|
(45) |
The second Piola-Kirchhoff stress 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 {\boldsymbol \sigma }}
of Eqs.(29–30) can be readily extracted from the 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}} tensor.
4 TOTAL LAGRANGIAN FORMULATION OF THE BASIC SHELL TRIANGLE
4.1 Definition of the element geometry and discretization of the displacement field
The rotation-free BST element has three nodes with three displacement degrees of freedom at each node. As for the BPT element a patch is defined by the central triangle and the three adjacent elements (Figure 1). This four elements patch helps to define the curvature field within the central triangle (the BST element) in terms of the displacements of the six patch nodes.
The node-ordering in the patch is the following (see Figure 1)
- The nodes in the main element (M) are numbered locally as 1, 2 and 3. They are defined counter-clockwise around the positive normal.
- The sides in the main element are numbered locally as 1, 2, and 3. They are defined by the local node opposite to the side.
- The adjacent elements (which are part of the patch) are numbered with the number associated to the common side.
- The extra nodes of the patch are numbered locally as 4, 5 and 6, corresponding to nodes on adjacent elements opposite to sides 1, 2 and 3 respectively.
- The connectivities in the adjacent elements are defined beginning with the extra node as shown in Table 1.
| Element | N1 | N2 | N3 |
| M | 1 | 2 | 3 |
| 1 | 4 | 3 | 2 |
| 2 | 5 | 1 | 3 |
| 3 | 6 | 2 | 1 |
The following local Cartesian coordinate system is defined for the patch. In the main element the unit 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}_{1}} (associated to the local coordinate 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}} ) is directed along side 3 (from node 1 to 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 \mathbf{t}_{3}}
(associated to the coordinate 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 unit normal to the plane, and finally 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}=\mathbf{t}_{3}\times \mathbf{t}_{1}}
(associated to the coordinate 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 _{2}}
).
The coordinates and the displacements are linearly interpolated within each three node triangle in the mesh in the standard manner, i.e.
|
(46) |
|
(47) |
In above 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 \varphi }_{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 \mathbf{u}_{i}}
contain respectively the three coordinates and the three displacements of 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}
.
4.2 Computation of the membrane strains
The Green-Lagrange membrane strains are expressed by substituting the linear displacement interpolation into Eq.(22). This gives
|
(48) |
The membrane strain field is constant within each triangle similarly as in the standard CST element. The variation of the membrane strains is obtained by
|
(49) |
with
|
(50) |
and
|
(51) |
4.3 Computation of bending strains (curvatures)
We will assume the following constant curvature field within each element
|
(52) |
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 \hat{\kappa }_{\alpha \beta }}
is the assumed constant curvature field defined by
|
(53) |
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}^{0}}
is the area (in the original configuration) of the central element in the patch.
Substituting Eq.(53) into (52) and integrating by parts the area integral gives the curvature vector within the element in terms of the following line integral
|
(54) |
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 n_{i}}
are the components (in the local system) of the normals to the element sides in the initial configuration 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}^{0}}
.
For the definition of 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}} , the linear interpolation over the central element is used. In this case the tangent plane components are
|
(55) |
|
(56) |
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^{0}_{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 the side 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}^{M}}
. Equation (56) 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 order to compute the line integral of equation (54) the averaging procedure described in Section 2 is used. Hence along each side of the triangle the average value 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 \varphi }_{^{\prime }\alpha }}
between the main triangle and the adjacent one is taken leading to
|
(57) |
where the sum extends over the three elements adjacent to the central triangle 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} .
Noting 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{t}_{3}\cdot \mathbf{\boldsymbol \varphi }_{^{\prime }\alpha }^{M}=0}
in the main triangle and using (6) it can be found [23]
|
(58) |
This can be seen as the projection of the local derivatives 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 \varphi }}
in the adjacent triangles (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{\boldsymbol \varphi }_{^{\prime }\alpha }^{i}}
where index 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}
denotes values associated to the adjacent elements) over the normal to the main triangle 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}}
. As the triangles have a common 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{t}_{3}\cdot \mathbf{\boldsymbol \varphi }_{^{\prime }s}^{i}=0} , 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{\boldsymbol \varphi }_{^{\prime }s}^{i}}
is the derivative 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 \varphi }}
along the side. Hence only the derivative 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 \varphi }}
along the side 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{\boldsymbol \varphi }_{^{\prime }n}^{i}}
) has non-zero component over 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}} . This gives
|
(59) |
An alternative form to express the curvatures, which is useful when their variations are needed, is to define the vectors
|
(60) |
This gives
|
(61) |
The 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. The variation of the curvatures can be obtained as
|
(62) |
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 {\boldsymbol \varphi }^{\alpha }}
have been included
|
(63) |
with
|
In above expressions superindexes 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 L_{j}^k}
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 \delta \mathbf{u}_{j}^k}
refer to element numbers in the patch whereas subscripts denote node numbers of each element in the patch. As usual the superindex 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}
denotes values in the central triangle (Figure 1). Note that as expected the curvatures (and their variations) in the central element are a function of the nodal displacements of the six nodes in the four elements patch. Note also the isochoric approach
|
(66) |
The derivation of Eq.(62) can be found in [26]. This equation can be rewritten in the form
|
(67) |
where
|
(68) |
is the virtual displacement vector of the patch and
|
(69) |
is the curvature matrix relating the virtual curvatures within the central element and the 18 virtual displacements of the six nodes in the patch.
The form of 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}_{b}}
is given in the Appendix.
5 ENHANCED BASIC SHELL TRIANGLE
An enhanced version of the BST element (termed EBST) has been recently proposed by Flores and Oñate [26]. The main features of the element formulation are the following:
- The geometry of the patch formed by the central element and the three adjacent elements is quadratically interpolated from the position of the six nodes in the patch.
- The membrane strains are assumed to vary linearly within the central triangle and are expressed in terms of the (continuous) values of the deformation gradient at the mid side points of the triangle.
- The assumed constant curvature field within the central triangle is obtained by expression (54) using now twice the values of the (continuous) deformation gradient at the mid side points.
Details of the derivation of the EBST element are given below.
5.1 Definition of the element geometry and computation of membrane strains
As mentioned above a quadratic approximation of the geometry of the four elements patch is chosen using the position of the six nodes in the patch. It is useful to define the patch in the isoparametric space using the nodal positions given in the Table 2 (see also Figure 2).
| 1 | 2 | 3 | 4 | 5 | 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 \xi } | 0 | 1 | 0 | 1 | -1 | 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 \eta } | 0 | 0 | 1 | 1 | 1 | -1 |
The quadratic interpolation is defined by
|
(70) |
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 \zeta=1-\xi-\eta} )
|
(71) |
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 the mid side points of the central element of the patch 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 G_{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 G_{2}}
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 G_{3}}
in Figure 2. This choice has the following advantages.
|
| Figure 2: Patch of elements in the isoparametric space. |
- Gradients at the three mid side points depend only on the nodes belonging to the two elements adjacent to each 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 of two adjacent elements, the same values are obtained, as the coordinates of the same four points are used. This in practice means that the gradients at the mid-side points are independent of the element where they are computed. A side-oriented implementation of the finite element will therefore lead to a unique evaluation of the gradients per side.
The Cartesian derivatives of the shape functions are computed at the original configuration by the standard expression
|
(72) |
where the Jacobian matrix at the original configuration is
|
(73) |
The deformation gradients on the middle surface, associated to an arbitrary spatial Cartesian system and to the material cartesian system defined on the middle surface are related by
|
(74) |
The membrane strains within the central triangle are now obtained using a linear assumed strain field 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 \hat{\boldsymbol \varepsilon }_{m}} . If, for example, Green Lagrange strains are used, i.e.
|
(75) |
with
|
(76) |
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 \varepsilon }_{m}^{i}}
are the membrane strains computed at the three mid side 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_{i}}
(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=1,2,3}
see Figure 2). In Eq.(76) 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}_{1}=(1-2\zeta )}
, etc.
The gradient at each mid side point is computed from the quadratic interpolation (70):
|
(77) |
Substituting Eq.(22) into (76) and using Eq.(20) gives the membrane strain vector as
|
(78) |
and the virtual membrane strains as
|
(79) |
We note that the gradient 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 G_{i}}
depends only on the coordinates of the three nodes of the central triangle and on those of an additional node in the patch, associated 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.
Combining Eqs.(79) and (77) gives
|
(80) |
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 \delta \mathbf{a}^{p}}
is the patch displacement vector (see Eq.(68)) 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{B}_{m}} is the membrane strain matrix. An explicit form of this matrix is given in the Appendix.
Differently from the original BST element the membrane strains within the EBST element are now a function of the displacements of the six patch nodes.
5.2 Computation of curvatures
The constant curvature field assumed for the BST element is chosen again here. The numerical evaluation of the line integral in Eq.(54) results in a sum over the integration points at the element boundary which are, in fact, the same points used for evaluating the gradients when computing the membrane strains. 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 integration 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_{i}} ). In this way the curvatures can be computed by
|
(81) |
In the standard BST element [20,23] 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{\boldsymbol \varphi }_{\prime \alpha }^{i}}
is computed as the average of the linear approximations over the two adjacent elements (see Section 4.3). In the enhanced version, the gradient is evaluated 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 G_{i}}
from the quadratic interpolation
|
(82) |
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 G_{i}} ) where the gradient is computed.
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}}
in Eq.(82) can be seen as a reference direction. If a different direction than that given by Eq.(56) 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.(56) 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 using the 6-node isoparametric interpolation.
The variation of the curvatures can be obtained as
|
(83) |
where the definitions (61) and (63) still hold but with the new 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{h}_{ij}}
given by [26]
|
(84) |
In Eq.(83)
|
(85) |
The expression of the curvature 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}_b}
is given in the Appendix. Details of the derivation of Eq.(83) can be found in [26].
5.3 The EBST1 element
A simplified and yet very effective version of the EBST element can be obtained by using one point quadrature for the computation of all the element integrals. This element is termed EBST1. Note that this only affects the membrane stiffness matrices and it is equivalent to using a assumed constant membrane strain field defined by an average of the metric tensors computed at each side.
Numerical experiments have shown that both the EBST and the EBST1 elements are free of spurious energy modes.
6 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 the bending 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 \boldsymbol \varphi _{^{\prime }n}^{0}}
that does not change during the process.
|
| Figure 3: Local Cartesian system for the treatment of symmetry boundary conditions |
The tangent plane at the boundary (mid-side point) is expressed in terms of two orthogonal unit vectors referred to a local-to-the-boundary Cartesian system (see Figure 3) defined as
|
(86) |
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 \boldsymbol \varphi _{^{\prime }n}^{0}}
is fixed during the process 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 \bar{\boldsymbol \varphi }_{^{\prime }s}}
emerges from 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 in the selected original convective Cartesian system (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},\mathbf{t}_{2} } ) is
|
(87) |
the intersection line (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} ) 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_{i}^{M}} , i.e.
|
(88) |
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}^{i} =\left[n_{1},n_{2}\right]^{T}=\left[\mathbf{n\cdot t}_{1},\mathbf{n\cdot t}_{2}\right]^{T}}
(resolved in the selected original convective Cartesian system) leads to
|
(89) |
where, noting 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 \lambda }
is the determinant of the gradient, 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 \boldsymbol \varphi _{^{\prime }n}^{i}}
can be approximated by
|
(90) |
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}_{\alpha \beta }}
(Eqs. 60 and 84) results in
|
(91) |
For the computation of the curvature variations, 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 now (see [26])
|
(92a) |
|
(92b) |
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 \boldsymbol \varphi _{^{\prime }n}}
has been neglected.
For a simple supported (hinged) side, the problem is not completely defined. The simplest choice is to neglect the contribution to the side rotations from the adjacent element missing in the patch in the evaluation of the curvatures via Eq.(54) [20,23]. This is equivalent to assume that the gradient at the side is equal to the gradient in the central element, i.e.
|
(93) |
More precise changes can be however introduced to account for the different natural boundary conditions. One may assume that the curvature normal to the side is zero, and consider a contribution of the missing side to introduce this constraint. As the change of curvature parallel to the side is also zero along the hinged side, this obviously leads to zero curvatures in both directions. Denoting the contribution to the curvatures of the interior 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 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}
) by
|
It can be easily shown that in order to set the normal curvature to zero the contribution of the simple supported 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} ) should be
|
(94) |
For the case of a triangle with two sides associated to hinged sides, the normal curvatures to both sides must be zero. Denoting 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{n}^{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 \mathbf{n}^{j}}
the normal to the sides, and 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{m}^{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 \mathbf{m}^{j}}
the dual base (associated to the base 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}^{i}-}
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}^{j}}
), the contribution from the hinged 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 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}
) can be written as a function of the contribution of the only interior 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 k} ):
|
(95) |
For a free edge the same approximation can be used but due to Poisson's effect this will lead to some error. The curvature variations of these contributions can be easily computed.
For the membrane formulation of element EBST, the gradient at the mid-side point of the boundary is assumed equal to the gradient of the main triangle.
7 IMPLICIT SOLUTION SCHEME
For a step 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}
the configuration 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{\boldsymbol \varphi }^{n}}
and the plastic strains 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 \varepsilon }_{p}^{n}}
are known. The configuration 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{\boldsymbol \varphi }^{n}}
is obtained by adding the total displacements to the original configuration 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{\boldsymbol \varphi }^{n}=\mathbf{\boldsymbol \varphi }^{0}+\mathbf{u}^{n}}
. The stresses are computed at each triangle using a single sampling (integration) point at the center 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_{L}}
integration points (layers) through the thickness. The plane stress state condition of the classical thin shell theory is assumed, so that for every layer three stress components are computed, (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 \sigma _{11}}
,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 \sigma _{22}} , 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 \sigma _{12}} ) referred to the local Cartesian system.
The computation of the incremental stresses is as follows:
- Evaluate the incremental displacements: 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 \mathbf{u}^{n}=\mathbf{K}_{T}^{n}\mathbf{r}^{n}}
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{K}_{T}}
is the tangent stiffness matrix and r is the residual force vector defined by for each 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/":): \mathbf{r}^e_i =\int \int _A L_i {\boldsymbol t}\, dA - \int \int _{A^\circ } ({\boldsymbol B}_{m_i}^T {\boldsymbol \sigma }_m + {\boldsymbol B}_{b_i}^T {\boldsymbol \sigma }_b)dA (96) - Generate the actual configuration 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{\boldsymbol \varphi }^{n+1}=\mathbf{\boldsymbol \varphi }^{n}+\Delta \mathbf{u}^{n}}
- Compute the 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 }^{n+1}\mathbf{ }}
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 \kappa _{\alpha \beta }^{n+1}}
. Then at each layer 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}
compute the (approximate) right Cauchy-Green tensor (27)
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/":): \mathbf{C}_{k}^{n+1}=\mathbf{a}^{n+1}+z_{k}{\boldsymbol \chi }^{n+1} (97) - Compute the total (32) and elastic (33) deformations at each layer 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}
- Compute the trial elastic stresses (34) at each layer 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}
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/":): \mathbf{T} _{k}^{n+1}=\mathbf{D}\left[ {\boldsymbol \varepsilon }_{e}\right] _{k}^{n+1} (99) - Check the plasticity condition and return to the plasticity surface. If necessary correct the plastic strains 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[{\boldsymbol \varepsilon }_{p}\right] _{k}^{n+1}} at each layer
- Compute the second Piola-Kirchhoff stress 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 \boldsymbol \sigma _k^{n+1}}
and the generalized 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/":): {\boldsymbol \sigma }^{n+1}_{m} =\frac{h^{0}}{N_{L}}\sum _{k=1}^{N_{L}}\boldsymbol \sigma _{k}^{n+1} w_{k} 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/":): {\boldsymbol \sigma }^{n+1}_{b} =\frac{h^{0}}{N_{L}}\sum _{k=1}^{N_{L}}\boldsymbol \sigma _{k}^{n+1}z_{k} w_{k} (100) - Compute the residual force vector. The contribution for the 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}
th element is given 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/":): (\mathbf{r}^{M})^{n+1}=-A_{M}^{0}\left[ \begin{array}{cc} \mathbf{B}_{m}^{T} & \mathbf{B}_{b}^{T}\end{array} \right] ^{n+1}\left[ \begin{array}{c} \boldsymbol \sigma _{m}\\ \boldsymbol \sigma _{b}\end{array} \right] ^{n+1} (101)
The expression of the tangent stiffness matrix for the element is given below. Details of the derivation can be found in [23],[26].
|
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 w_{k}}
is the weight of the through-the-thickness integration point. Recall 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 z_{k}}
is the current distance of the layer to the mid-surface and not the original distance. However, for small strain plasticity this distinction is not important.
This computation of stresses is adequate for an implicit scheme independent of the step size and it is exact for an elastic problem.
7.1 Tangent stiffness matrix
As usual the tangent stiffness matrix is split into material and geometric components. The material tangent stiffness matrix is computed by the integral
|
(102) |
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{B}=\mathbf{B}_{m}+\mathbf{B}_{b}}
includes:
- a membrane contribution 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}_{m}}
given by Eq.(51) or Eq.(80).
- a bending contribution 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}_{b}}
given by Eq.(69) or Eq.(85) which is constant over the element.
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{D}_{ep}}
is the elastic-plastic tangent constitutive matrix integrated through the thickness.
A three point quadrature is used for integrating the stiffness terms 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}_{m}^{T}\mathbf{D}_{ep}\mathbf{B}_{m}}
(recall that for the EBST element the membrane strains vary linearly within the element) whereas one point quadrature is chosen for the rest of 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 \mathbf{K}^{M}}
.
7.2 Geometric tangent stiffness matrix
The geometric stiffness is written as
|
(103) |
where subscripts 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}
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}
denote as usual membrane and bending contributions. For the BST element the membrane part is the same than for the standard constant strain triangle, leading to
|
(104) |
For the EBST element the membrane part is computed as the sum of the contributions of the three sides, i.e.
|
(105) |
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 N_{ij}={\sigma _{m}}_{ij}}
are the axial forces defined in Eq.(29).
The geometric stiffness associated to bending moments is much more involved and can be found in [26]. Numerical experiments have shown that the bending part of the geometric stiffness is not so important and can be disregarded in the iterative process.
Again three and one point quadratures are used for computing the membrane and bending contributions to the geometric stiffness matrix. We note that for elastic-plastic problems a uniform one point quadrature has been chosen for integrating both the membrane and bending stiffness matrices.
8 EXPLICIT SOLUTION SCHEME
For simulations including large non-linearities, such as frictional contact conditions on complex geometries or large instabilities in membranes, convergence is difficult to achieve with implicit schemes. In those cases an explicit solution algorithm is typically most advantageous. This scheme provides the solution for dynamic problems and also for static problems if an adequate damping is chosen.
The dynamic equations of motion to solve are of the form
|
(106) |
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{M}}
is the mass 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}}
is the damping matrix and the dot means the time derivative. The solution is performed using the central difference method. To make the method competitive a diagonal (lumped) 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{M}}
matrix is typically used 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}}
is taken proportional 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 \mathbf{M}}
. As usual, mass lumping is performed by assigning one third of the triangular element mass to each node in the central element.
The explicit solution scheme can be summarized as follows. At each time step 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}
where displacements have been computed:
- Compute the internal forces 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}^{n}} . This follows the same steps (2-8) described for the implicit scheme in the previous section.
- Compute the accelerations at time 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_{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/":): \ddot{\mathbf{u}}^{n} = {\boldsymbol M}_d^{-1} [ \mathbf{r}^{n} - \mathbf{C} \dot{\mathbf{u}}^{n-1/2} ] - Compute the velocities at time 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_{n+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/":): \dot{\mathbf{u}}^{n+1/2} = \dot{\mathbf{u}}^{n-1/2} \ddot{\mathbf{u}}^{n} \delta t - Compute the displacements at time 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_{n+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/":): \mathbf{u}^{n+1} = \mathbf{u}^{n} +\dot{\mathbf{u}}^{n+1/2} \delta t - Update the shell geometry
- Check frictional contact conditions
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 M}_d}
is the diagonal (lumped) mass matrix.
Further details of the implementation of the standard BST element within an explicit solution scheme can be found in [25].
9 EXAMPLES
In this section several examples are presented to show the good performance of the rotation-free shell elements BST, EBST and EBST1. The first five static examples are solved using an implicit code. The rest of the examples are computed using the explicit dynamic scheme. For the explicit scheme the EBST element is always integrated using one integration point per element (EBST1 version) although not indicated.
9.1 Patch tests
The three elements considered (BST, EBST and EBST1) satisfy the membrane patch test defined in Figure 4. A uniform axial tensile stress is obtained in all cases.
|
| Figure 4: Patch test for uniform tensile stress |
|
| Figure 5: 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 5 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 shifted from the center, the results obtained with the EBST and EBST1 elements are not correct. This however does not seems to preclude the excellent performance of these elements, as proved in the rest of the examples analyzed. On the other hand, the BST element gives correct results in all torsion patch tests if natural boundary conditions are imposed in the formulation. If this is not the case, incorrect results are obtained even with structured meshes.
9.2 Cook's membrane problem
This example is used to assess the membrane performance of the EBST and EBST1 elements and to compare it with the standard 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 6.a shows the geometry and the applied load. Figure 6.b plots the vertical displacement of the upper vertex as a function of the number of nodes in the mesh. Results obtained with other isoparametric elements have also been plotted for comparison. They include the constant strain triangle (CST), the bilinear quadrilateral (QUAD4) and the linear strain triangle (LST) [4]. Note that as this is a pure membrane problem the BST and the CST elements give identical results.
|
|
| Figure 6: Cook's membrane problem (a) Geometry (b) Results | |
From the plot shown it can be seen that the enhanced element with three integration points (EBST) 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 than the former. On the other hand, if a one point quadrature is used (EBST1) the convergence in the reported displacement is notably better than for the rest of the elements.
9.3 Cylindrical roof
In this example an effective membrane interpolation is of primary importance. The geometry is a cylindrical roof supported by a rigid diaphragm at both ends and it is loaded by a uniform dead weight (see Figure 7.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 7.a shows orientation B).
Tables 3, 4 and 5 present the normalized vertical displacements at the crown (point A) and at the midpoint of the free side (point B) for the two orientations of the structured meshes and for the non-structured mesh. 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 [31].
|
|
| Figure 7: Cylindrical roof under dead weight. 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 \times 10^{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. (a) Geometry and mesh for orientation B. (b) Displacement of point B for both (structured) mesh orientations | |
| Point-A | Point-B | |||||
| NDOFs | EBST | EBST1 | BST | CBST | EBST1 | BST |
| 16 | 0.65724 | 0.91855 | 0.74161 | 0.40950 | 0.70100 | 1.35230 |
| 56 | 0.53790 | 1.03331 | 0.74006 | 0.54859 | 1.00759 | 0.75590 |
| 208 | 0.89588 | 1.04374 | 0.88491 | 0.91612 | 1.02155 | 0.88269 |
| 800 | 0.99658 | 1.01391 | 0.96521 | 0.99263 | 1.00607 | 0.96393 |
| 3136 | 1.00142 | 1.00385 | 0.99105 | 0.99881 | 1.00102 | 0.98992 |
| Point-A | Point-B | |||||
| NDOFs | EBST | EBST1 | BST | CBST | EBST1 | BST |
| 16 | 0.26029 | 0.83917 | 0.40416 | 0.52601 | 0.86133 | 0.89778 |
| 56 | 0.81274 | 1.10368 | 0.61642 | 0.67898 | 0.93931 | 0.68238 |
| 208 | 0.97651 | 1.04256 | 0.85010 | 0.93704 | 1.00255 | 0.86366 |
| 800 | 1.00085 | 1.01195 | 0.95626 | 0.99194 | 1.00211 | 0.95864 |
| 3136 | 1.00129 | 1.00337 | 0.98879 | 0.99828 | 1.00017 | 0.98848 |
| Point-A | Point-B | |||||
| NDOFs | EBST | EBST1 | BST | EBST | EBST1 | BST |
| 851 | 0.97546 | 0.8581 | 0.97598 | 0.97662 | 1.0027 | 0.97194 |
| 3311 | 0.98729 | 0.9682 | 0.98968 | 0.98476 | 1.0083 | 0.98598 |
| 13536 | 0.99582 | 0.9992 | 1.00057 | 0.99316 | 0.9973 | 0.99596 |
Plots in Figure 7.b show the normalized displacement of point-B for structured meshes as a function of the number of degrees of freedom for each case studied. An excellent convergence for the EBST element can be seen. The version with only one integration point (EBST1) presents a behavior a little more flexible and converges from above for structured meshes. Table 5 shows that both the EBST and the EBST1 elements have an excellent behavior for non structured meshes.
9.4 Open semi-spherical dome with point loads
The main problem of shell finite elements with initially curved geometry is the so called membrane locking. The EBST element 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 8.a shows the discretized geometry (only one quarter of the geometry is considered due to symmetry).
|
|
| Figure 8: Pinched hemispherical shell with a hole, (a)geometry, (b)normalized displacement | |
In Figure 8.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: three membrane locking free elements, namely the original linear BST element, a transverse shear-deformable quadrilateral (QUAD) [32] and an assumed strain quadratic triangle (TRIC) [3]; a transverse shear deformable quadratic triangle (TRIA) (standard displacement formulation for membrane part) [2] that is vulnerable to locking.
From the plotted results it can be seen that the EBST element presents slight membrane locking in bending dominated problems with initially curved geometries. This locking is much less severe than in a standard quadratic triangle. Membrane locking disappears when only one integration point is used (EBST1 element).
9.5 Inflation of a sphere
The example is the inflation of a spherical shell under internal pressure. An incompressible Mooney-Rivlin constitutive material has 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 [33] (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^{0}} ):
|
In this numerical simulation the same geometric and material parameters used in [22] 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^{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 h^{0}=0.02}
. The three meshes of EBST1 element considered to evaluate convergence are shown in Figure 9.a. The value of the actual radius as a function of the internal pressure is plotted in Figure 9.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 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} .
| |
| (a) | |
| |
| (b) | |
| Figure 9: Inflation of sphere of Mooney-Rivlin material. (a) Meshes of EBST1 elements used in the analysis. (b) Radius as a function of the internal pressure. | |
9.6 Clamped spherical dome under impulse pressure loading
The geometry of the dome and the material properties chosen are shown in Figure 10. A uniform pressure load of 600 psi is applied to the upper surface of the dome. The different meshes used in the analysis are shown in Figure 11. One fourth of the dome is considered only due to symmetry. Two different analyses under elastic and elastic-plastic conditions were carried out. The number of thickness layers in eq.(100) is four. Numerical experiments show that this suffices to provide an accurate solution for large elastic-plastic problems [25]. Results are obtained using the explicit scheme.
|
| Figure 10: Spherical dome under impulse pressure. Geometry and material |
|
| Figure 11: Spherical dome under impulse pressure. Meshes used in the analysis. Mesh-1 with 338 elements, Mesh-2 with 1250 elements and Mesh-3 with 2888 elements |
Figure 12 shows results for the time history of the central deflection using different meshes and elastic material properties for both BST and EBST1 elements. Results are almost identical for mesh-2 and mesh-3, showing the excellent convergence properties. The coarsest mesh shows some differences between both elements, but for the finer meshes the results are almost identical. Figure 13 shows similar results but now for an elastic-plastic material. The excellent convergence of the BST and EBST1 elements is again noticeable.
|
| Figure 12: Spherical dome under impulse pressure. History of central deflection for elastic material |
|
| Figure 13: Spherical dome under impulse pressure. History of central deflection for elastic-plastic material |
Results obtained with the present elements compare very well with published results using fine meshes. See for example ABAQUS Explicit example problems manual [34] and WHAMS-3D manual [35], showing plots comparing results using different shell elements.
A summary of results for the central deflection at significant times is given in Tables 6 and 7. Further details on the solution of this problem with the standard BST element can be found in [25].
| Element/mesh | 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/":): t = 0.2 ms | 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/":): t = 0.4 ms | 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/":): t = 0.6 ms | 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/":): t = 0.8 ms |
| BST Coarse | -0.05155 | -0.09130 | 0.04414 | -0.08945 |
| BST Medium | -0.04542 | -0.09177 | 0.03863 | -0.08052 |
| BST Fine | -0.04460 | -0.09022 | 0.03514 | -0.08132 |
| EBST1 Coarse | -0.05088 | -0.08929 | 0.04348 | -0.08708 |
| EBST1 Medium | -0.04527 | -0.09134 | 0.03865 | -0.07979 |
| EBST1 Fine | -0.04453 | -0.09004 | 0.03510 | -0.08099 |
| Element/mesh | 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/":): t = 0.2 ms | 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/":): t = 0.4 ms | 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/":): t = 0.6 ms | 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/":): t = 0.8 ms |
| BST Coarse | -0.05888 | -0.05869 | -0.02938 | -0.06521 |
| BST Medium | -0.05376 | -0.06000 | -0.02564 | -0.06098 |
| BST Fine | -0.05312 | -0.05993 | -0.02464 | -0.06105 |
| EBST1 Coarse | -0.05827 | -0.05478 | -0.02792 | -0.06187 |
| EBST1 Medium | -0.05374 | -0.05884 | -0.02543 | -0.06080 |
| EBST1 Fine | -0.05317 | -0.05935 | -0.02458 | -0.06085 |
9.7 Cylindrical panel under impulse loading
The geometry of the cylinder and the material properties are shown in Figure 14. A prescribed initial normal velocity 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 v_{o}=-5650}
in/sec is applied to the points in the region shown modelling the effect of the detonation of an explosive layer. The panel is assumed clamped along all the boundary. One half of the cylinder is discretized only due to symmetry conditions. Three different meshes 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 6\times{12}}
, 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 12\times{32}}
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 18\times{48}}
triangles were used for the analysis. The deformed configurations for 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 time =1 msec}
are shown for the three meshes in Figure 15.
|
| Figure 14: Cylindrical panel under impulse loading. Geometry and material properties |
|
| Figure 15: Impulsively loaded cylindrical panel. Deformed meshes for 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/":): time =1 msec |
The analysis was performed assuming an elastic-perfect plastic material behaviour (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 \sigma _y = k_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 k_y'=0}
). A study of the convergence of the solution with the number of thickness layers showed again that four layers suffice to capture accurately the non linear material response [25].
A comparison of the results obtained with the BST and EBST1 elements using the coarse mesh and the finer mesh is shown in Figure 16 where experimental results reported in [36] have also been plotted for comparison purposes. Good agreement between the numerical and experimental results is obtained. Figures 16 show the time evolution of the vertical displacement of two reference points along the center line located at 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 y=6.28} in 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 y=9.42} in, respectively. For the finer mesh results between both elements are almost identical. For the coarse mesh it can been seen again that the element BST is more flexible than element EBST1.
|
Figure 16: Cylindrical panel under impulse loading. Time evolution of the displacement of two points along the crown line. Comparison of results obtained with BST and EBST1 elements (mesh 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/":): 6\times{12}
elements and mesh 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/":): 18\times{48}
elements) and experimental values
|
The numerical values of the vertical displacement at the two reference points obtained with the BST and EBST1 elements after a time of 0.4 ms using the 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 16\times{32}}
mesh are compared in Table 8 with a numerical solution obtained by Stolarski et al. [37] using a curved triangular shell element and the 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 16\times{32}} mesh. Experimental results reported in [36] are also given for comparison. It is interesting to note the reasonable agreement of the results for 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 y=6.28}
in. and the discrepancy of present and other published numerical solutions with the experimental value for 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 y=9.42} in.
| Vertical displacement (in.) | ||
| element/mesh | 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/":): y=6.28
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/":): y=9.42
in |
BST (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 6\times 12}
el.) |
-1.310 | -0.679 |
BST (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 18\times 48}
el.) |
-1.181 | -0.587 |
EBST1 (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 6\times 12}
el.) |
-1.147 | -0.575 |
EBST1 (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 18\times 48}
el.) |
-1.171 | -0.584 |
| Stolarski et al. [37] | -1.183 | -0.530 |
| Experimental [36] | -1.280 | -0.700 |
The deformed shapes of the transverse section for 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 y=6.28}
in. and the longitudinal section for 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=0}
obtained with the both elements for the coarse and the fine meshes after 1ms. are compared with the experimental results in Figures 17 and 18. Excellent agreement is observed for the fine mesh for both elements.
|
| Figure 17: Cylindrical panel under impulse loading. Final deformation (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/":): t=1 msec
) of the panel at the cross section 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/":): y=6.28 in Comparison with experimental values. |
|
| Figure 18: Cylindrical panel under impulse loading. Final deformation (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/":): t=1 msec
) of the panel at the crown line (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/":): x=0.00 in ). Comparison with experimental values. |
9.8 Airbag Membranes
Inflation/deflation of a circular airbag
This example has been taken from Ref.[22] 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 [22] it is also discussed the important regularizing properties of the bending energy that, when disregarded, leads to massive wrinkling in the compressed zones.
The airbag geometry is initially circular with an undeformed radius 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 0.35} . The 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=6\times 10^{7}} Pa, 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.3}
and 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 \rho = 2000}
kg/mFailed 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} . A symmetrical solution has been assumed and, hence, only one quarter of the geometry has been modelled. Only the normal displacement 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.0004} m 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} Pa. Pressure is linearly increased 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.15}
sec.
Figure 19 shows the final deformed configurations for a mesh with 10201 nodes and 20000 EBST1 elements. The figure on the left (a) corresponds to an analysis including full bending effects and the right figure (b) is a pure membrane analysis.
We note that when the bending energy is included a more regular final pattern is obtained. Also the final pattern is rather independent of the discretization (note that the solution is non unique due to the strong instabilities). On the other hand, the pure membrane solution shows in the center of the modelled region a wrinkling pattern where the width of the wrinkle is the length of the element.
|
| Figure 19: Inflation of a circular airbag. Deformed configurations for final pressure. (a) bending formulation; (b) membrane formulation. |
Figure 20 shows the results obtained for the de-inflation process. Note that the spherical membrane falls down due to the self weight. The configurations are, of course, non-unique.
|
| Figure 20: Inflation and deflation of a circular air-bag. |
Inflation/deflation of a closed tube
The next problem is the study of the inflating and de-inflating of a tube with a semi-spherical end cap. The tube diameter 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 D=2} m, its total length 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 L=5} m and the 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=5\times 10^{-3}} m. The material has the following properties: 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=4\times 10^{8}} Pa, 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.35 }
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 \rho =3\times 10^{3}}
kg/mFailed 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} . The tube is inflated fast until a pressure 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 10^4} Pa and then is de-inflated under self weight. The analysis is performed with a mesh of 4176 EBST1 elements and 2163 nodes modelling a quarter of the geometry. The evolution of the tube walls during the de-inflating process can be seen in Figure 21. Note that the central part collapses as expected, while the semi-spherical cap remains rather unaltered.
The same analysis is repeated for a longer and thinner tube (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=6} m 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=3\times 10^{-3}} m). The same material than in the previous case was chosen. The evolution of the tube walls is shown in Figure 22. Note that the central part collapses again but in a less smoother manner due to the smaller thickness.
|
| Figure 21: Inflation and deflation of a closed tube. 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/":): L=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/":): D=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/":): h=5\times 10^{-4} . |
|
| Figure 22: Inflation and deflation of a closed tube. 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/":): L=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/":): D=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/":): h=3\times 10^{-4} . |
Inflation of a square airbag
The last example of this kind is the inflation of a square airbag supporting a spherical object. This example resembles a problem studied (numerically and experimentally) in [38], where fluid-structure interaction is the main subject. Here the fluid is not modelled and a uniform pressure is applied over all the internal surfaces. The lower surface part of the airbag is limited by a rigid plane and on the upper part a spherical dummy object is set to study the interaction between the airbag and the object.
The airbag geometry is initially square with an undeformed side length of 0.643m. The constitutive material chosen is a linear isotropic elastic one 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 E=5.88\times 10^8} Pa, 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}
and a density 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 \rho = 1000}
kg/mFailed 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}
. Only one quarter of the geometry has been modelled due to symmetry. The 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.00075} m and the inflation pressure is 250000Pa. 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.15} sec. The spherical object has a 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=0.08} m and a mass of 4.8kg (one quarter) and is subjected to gravity load during all the process.
The mesh has 8192 EBST1 elements and 4225 nodes on each surface of the airbag. Figure 23 shows the deformed configurations for three different times. The sequence on the left of the figure corresponds to an analysis including full bending effects and that on the right is the result of a pure membrane analysis. A standard penalty formulation is used in order to treat the frictionless contact.
|
| Figure 23: Inflation of a square airbag against an spherical object. Deformed configurations for different times. Left figure: results obtained with the full bending formulation. Right figure: results obtained with a pure membrane solution. |
9.9 S-rail sheet stamping
The final problem corresponds to one of the sheet stamping benchmark tests proposed in NUMISHEET'96 [39]. The analysis comprises two parts, namely, simulation of the stamping of a S-rail sheet component and springback computations once the stamping tools are removed. Figure 24 shows the deformed sheet after springback.
|
| Figure 24: Stamping of a S-rail. Final deformation of the sheet after springback obtained in the simulation. The triangular mesh of the deformed sheet is also shown |
The detailed geometry and material data can be found in the proceedings of the conference [39] or in the web [40]. The mesh used for the sheet has 6000 three node triangular elements and 3111 points (Figure 24). The tools are treated as rigid bodies. The meshes used for the sheet and the tools are those provided by the benchmark organizers. The material considered here is a mild steel (IF) with Young Modulus 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.06 GPa}
and Poisson 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.3}
. Mises yield criterion was used for plasticity behaviour with non-linear isotropic hardening defined 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 \sigma _y(e^p) = 545(0.13+e^p)^{0.267} [MPa]} . A uniform friction of 0.15 was used for all the tools. A low (10kN) blank holder force was considered in this simulation.
Figure 25 compares the punch force during the stamping stage obtained with both BST and EBST1 elements for the simulation and experimental values. Also for reference the average values of the simulations presented in the conference are included. Explicit and implicit simulations are considered as different curves. There is a remarkable coincidence between the experimental values and the results obtained with BST and EBST1 elements.
|
| Figure 25: Stamping of a S-rail. Punch force versus punch travel. Average of explicit and implicit results reported at the benchmark are also shown. |
Figure 26 plots the 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 Z}
coordinate along line B"–G" after springback. The top surface of the sheet does not remain plane due to some instabilities due to the low blank holder force used. Results obtained with the simulations compare very well with the experimental values.
|
| Figure 26: Stamping of a S-rail. Z-coordinate along line B–-G after springback. Average of explicit and implicit results reported at the benchmark are also shown. |
10 CONCLUDING REMARKS
We have presented in the paper two alternative formulations for the rotation-free basic shell triangle (BST) using an assumed strain approach. The simplest element of the family is based on an assumed constant curvature field expressed in terms of the nodal deflections of a patch of four elements and a constant membrane field computed from the standard linear interpolation of the displacements within each triangle. An enhanced version of the BST element is obtained by using a quadratic interpolation of the geometry in terms of the six patch nodes. This allows to compute an assumed linear membrane strain field which improves the in-plane behaviour of the original element. A simple and economic version of the new EBST element using a single integration point has been presented. The efficiency of the different rotation-free shell triangles has been demonstrated in many examples of application including linear and non linear analysis of shells under static and dynamic loads, the inflation and de-inflation of membranes and a sheet stamping problem.
The enhanced rotation-free basic shell triangle element with a single integration point (the EBST1 element) has proven to be an excellent candidate for solving practical engineering shell and membrane problems involving complex geometry, dynamics, material non linearity and frictional contact conditions.
ACKNOWLEDGEMENTS
The problems analyzed with the explicit formulation were solved with the computer code STAMPACK [41] where the rotation-free elements here presented have been implemented. The support of the company QUANTECH (www.quantech.es) providing the code STAMPACK is gratefully acknowledged.
REFERENCES
[1] E. Oñate. A review of some finite element families for thick and thin plate and shell analysis. Publication CIMNE N.53, May 1994.
[2] 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.
[3] J.H. Argyris, M. Papadrakakis, C. Apostolopoulou and S. Koutsourelakis. The TRIC element. Theoretical and numerical investigation. Comput. Meth. Appl. Mech. Engrg., 182, 217–245, 2000.
[4] O.C. Zienkiewicz and R.L. Taylor. The finite element method. Solid Mechanics. Vol II, Butterworth-Heinemann, 2000.
[5] H. Stolarski, T. Belytschko and S.-H. Lee. A review of shell finite elements and corotational theories. Computational Mechanics Advances, Vol. 2 N.2, North-Holland, 1995.
[6] E. Ramm and W.A. Wall. Shells in advanced computational environment. In V World Congress on Computational Mechanics, J. Eberhardsteiner, H. Mang and F. Rammerstorfer (Eds.), Vienna, Austria, July 7–12, 2002. http://wccm.tuwien.ac.at.
[7] D. Bushnell and B.O. Almroth, “Finite difference energy method for non linear shell analysis”, J. Computers and Structures, Vol. 1, 361, 1971.
[8] S.P. Timoshenko. Theory of Plates and Shells, McGraw Hill, New York, 1971.
[9] A.C. Ugural. Stresses in Plates and Shells, McGraw Hill, New York, 1981.
[10] R.A. Nay and S. Utku. An alternative to the finite element method. Variational Methods Eng., Vol. 1, 1972.
[11] J.K. Hampshire, B.H.V. Topping and H.C. Chan. Three node triangular elements with one degree of freedom per node. Engng. Comput. Vol. 9, pp. 49–62, 1992.
[12] R. Phaal and C.R. Calladine. A simple class of finite elements for plate and shell problems. I: Elements for beams and thin plates. Int. J. Num. Meth. Engng., Vol. 35, pp. 955–977, 1992.
[13] R. Phaal and C.R. Calladine. A simple class of finite elements for plate and shell problems. II: An element for thin shells with only translational degrees of freedom. Int. J. Num. Meth. Engng., Vol. 35, pp. 979–996, 1992.
[14] E. Oñate and Cervera M. Derivation of thin plate bending elements with one degree of freedom per node. Engineering Computations, Vol. 10, pp 553–561, 1993.
[15] 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.
[16] 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 Technology, Ed. M.J.M. Barata Marques, 18th IDDRG Biennial Congress, Lisbon, 1994.
[17] J. Rojek and E. Oñate. Sheet springback analysis using a simple shell triangle with translational degrees of freedom only. Int. J. of Forming Processes, Vol. 1, No. 3, 275–296, 1998.
[18] J. Rojek, E. Oñate and E. Postek. Application of explicit FE codes to simulation of sheet and bulk forming processes. J. of Materials Processing Technology, Vols. 80-81, 620–627, 1998.
[19] J. Jovicevic and E. Oñate. Analysis of beams and shells using a rotation-free finite element-finite volume formulation, Monograph 43, CIMNE, Barcelona, 1999.
[20] E. Oñate and F. Zárate. Rotation-free plate and shell triangles. Int. J. Num. Meth. Engng., 47, pp. 557–603, 2000.
[21] 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, 2039-2072.
[22] 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.
[23] F.G. Flores 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.
[24] G. Engel, K. Garikipati, T.J.R. Hughes, M.G. Larson, L. Mazzei and R.L. Taylor. Continuous/discontinuous finite element approximation of fourth-order elliptic problems in structural and continuum mechanics with applications to thin beams and plates, and strain gradient elasticity. Comput. Methods Appl. Mech. Engrg., Vol. 191, 3669–3750, 2002.
[25] E. Oñate, P. Cendoya and J. Miquel. Non linear explicit dynamic analysis of shells using the BST rotation-free triangle. Engineering Computations, 19 (6), 662–706, 2002.
[26] F.G. Flores and E. Oñate. Improvements in the membrane behaviour of the three node rotation-free BST shell triangle using an assumed strain approach. Computer Methods in Applied Mechanics and Engineering, in press, 2003.
[27] O.C. Zienkiewicz and E. Oñate. Finite Elements vs. Finite Volumes. Is there really a choice?. Nonlinear Computational Mechanics. State of the Art. (Eds. P. Wriggers and R. Wagner). Springer Verlag, Heidelberg, 1991.
[28] E. Oñate, M. Cervera and O.C. Zienkiewicz. A finite volume format for structural mechanics. Int. J. Num. Meth. Engng., 37, pp. 181–201, 1994.
[29] R. Hill. A Theory of the Yielding and Plastic Flow of Anisotropic Metals. Proc. Royal Society London, Vol. A193, pp. 281, 1948.
[30] 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.
[31] H.C. Huang, Static and Dynamic Analysis of Plates and Shells, page 40, Springer-Verlag, Berlin, 1989.
[32] E.N. Dvorkin and K.J. Bathe. A continuum mechanics based four node shell element for general non-linear analysis. Eng. Comp., 1, 77–88, 1984.
[33] A. Needleman. Inflation of spherical rubber ballons. Int. J. of Solids and Structures, 13, 409–421, 1977.
[34] Hibbit, Karlson and Sorensen Inc. ABAQUS, version 5.8, Pawtucket, USA, 1998.
[35] WHAMS-3D. An explicit 3D finite element program. KBS2 Inc., Willow Springs, Illinois 60480, USA.
[36] H.A. Balmer and E.A. Witmer. Theoretical experimental correlation of large dynamic and permanent deformation of impulsively loaded simple structures. Air force flight Dynamic Lab. Rep. FDQ-TDR-64-108, Wright-Patterson AFB, Ohio, USA, 1964.
[37] H. Stolarski, T. Belytschko and N. Carpenter. A simple triangular curved shell element. Eng. Comput., Vol. 1, 210–218, 1984.
[38] P.O. Marklund and L. Nilsson. Simulation of airbag inflation processes using a coupled fluid structure approach. Computational Mechanics, 29, 289–297, 2002.
[39] NUMISHEET'96, Third International Conference and Workshop on Numerical Simulation of 3D Sheet Forming Processes, NUMISHEET'96, E.H. Lee, G.L. Kinzel and R.H. Wagoner (Eds.), Dearbon-Michigan, USA, 1996.
[40] http://rclsgi.eng.ohio-state.edu/%Elee-j-k/numisheet96/
[41] STAMPACK. A General Finite Element System for Sheet Stamping and Forming Problems, Quantech ATZ, Barcelona, Spain, 2003 (www.quantech.es).
APPENDIX
A.1Curvature matrix for the BST element
|
with
|
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/":): \mathbf{B}_{b}^{T}=
| 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^{M}_{2,1} L^{2}_{2,1} +L^{M}_{3,1} L^{3}_{3,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/":): +L^{M}_{2,2} L^{2}_{2,2} +L^{M}_{3,2} L^{3}_{3,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/":): +L^{M}_{2,2} L^{2}_{2,1} +L^{M}_{2,1} L^{2}_{2,2} +L^{M}_{3,2} L^{3}_{3,1} +L^{M}_{3,1} L^{3}_{3,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 L^{M}_{1,1} L^{1}_{3,1} +L^{M}_{3,1} L^{3}_{2,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/":): L^{M}_{1,2} L^{1}_{3,2} +L^{M}_{3,2} L^{3}_{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/":): L^{M}_{1,2} L^{1}_{3,1} +L^{M}_{1,1} L^{1}_{3,2} +L^{M}_{3,2} L^{3}_{2,1} +L^{M}_{3,1} L^{3}_{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 L^{M}_{1,1} L^{1}_{2,1} +L^{M}_{2,1} L^{2}_{3,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/":): L^{M}_{1,2} L^{1}_{2,2} +L^{M}_{2,2} L^{2}_{3,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/":): L^{M}_{1,2} L^{1}_{2,1} +L^{M}_{1,1} L^{1}_{j,3} +L^{M}_{2,2} L^{2}_{3,1} +L^{M}_{2,1} L^{2}_{3,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 L^{M}_{1,1} L^{1}_{1,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/":): L^{M}_{1,2} L^{1}_{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/":): L^{M}_{1,2} L^{1}_{1,1} +L^{M}_{1,1} L^{1}_{1,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 L^{M}_{2,1} L^{2}_{1,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/":): L^{M}_{2,2} L^{2}_{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/":): L^{M}_{2,2} L^{2}_{1,1} +L^{M}_{2,1} L^{2}_{1,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 L^{M}_{3,1} L^{3}_{1,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/":): L^{M}_{3,2} L^{3}_{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/":): L^{M}_{3,2} L^{3}_{1,1} +L^{M}_{3,1} L^{3}_{1,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 -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 L^{M}_{1,1}\rho _{11}^{1}+L^{M}_{1,2}\rho _{11}^{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/":): L^{M}_{1,1}\rho _{22}^{1}+L^{M}_{i,2}\rho _{22}^{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/":): L^{M}_{1,1}\rho _{12}^{1}+L^{M}_{1,2}\rho _{12}^{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 L^{M}_{2,1}\rho _{11}^{1}+L^{M}_{2,2}\rho _{11}^{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/":): L^{M}_{2,1}\rho _{22}^{1}+L^{M}_{2,2}\rho _{22}^{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/":): L^{M}_{2,1}\rho _{12}^{1}+L^{M}_{2,2}\rho _{12}^{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 L^{M}_{3,1}\rho _{11}^{1}+L^{M}_{3,2}\rho _{11}^{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/":): L^{M}_{3,1}\rho _{22}^{1}+L^{M}_{3,2}\rho _{22}^{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/":): L^{M}_{3,1}\rho _{12}^{1}+L^{M}_{3,2}\rho _{12}^{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 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/":): 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/":): 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 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/":): 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/":): 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 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/":): 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/":): 0
|
A.2 Membrane strain matrix and curvature matrix for the EBST element
A.2.1 Membrane strain 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/":): \mathbf{B}_{m}^{T}=\frac{1}{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 N^{1}_{1,1}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }1} + N^{2}_{1,1}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }1} + N^{3}_{1,1}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }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/":): N^{1}_{1,2}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }2} + N^{2}_{1,2}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }2} + N^{3}_{1,2}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }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 N^{1}_{2,1}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }1} + N^{2}_{2,1}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }1} + N^{3}_{2,1}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }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/":): N^{1}_{2,2}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }2} + N^{2}_{2,2}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }2} + N^{3}_{2,2}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }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 N^{1}_{3,1}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }1} + N^{2}_{3,1}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }1} + N^{3}_{3,1}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }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/":): N^{1}_{3,2}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }2} + N^{2}_{3,2}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }2} + N^{3}_{3,2}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }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 N^{1}_{4,1}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }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/":): N^{1}_{4,2}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }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 N^{2}_{5,1}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }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/":): N^{2}_{5,2}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }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 N^{3}_{6,1}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }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/":): N^{3}_{6,2}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }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 N^{1}_{1,2}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }1} +N^{1}_{1,1}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }2} +N^{2}_{1,2}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }1} +N^{2}_{1,1}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }2} +N^{3}_{1,2}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }1} +N^{3}_{1,1}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }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 N^{1}_{2,2}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }1} +N^{1}_{2,1}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }2} + N^{2}_{2,2}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }1} +N^{2}_{2,1}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }2} + N^{3}_{2,2}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }1} +N^{3}_{2,1}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }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 N^{1}_{3,2}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }1} +N^{1}_{3,1}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }2} +N^{2}_{3,2}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }1} +N^{2}_{3,1}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }2} +N^{3}_{3,2}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }1} +N^{3}_{3,1}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }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 N^{1}_{4,2}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }1} +N^{1}_{4,1}\mathbf{\boldsymbol \varphi }^{1}_{^{\prime }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 N^{2}_{5,2}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }1} +N^{2}_{5,1}\mathbf{\boldsymbol \varphi }^{2}_{^{\prime }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 N^{3}_{6,2}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }1} +N^{3}_{6,1}\mathbf{\boldsymbol \varphi }^{3}_{^{\prime }2} }
|
A.2.2 Curvature 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/":): \mathbf{B}_{b}^{T}=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 L_{1,1}\left(N_{1,1}\right)_{G_{1}} +L_{2,1}\left(N_{1,1}\right)_{G_{2}} +L_{3,1}\left(N_{1,1}\right)_{G_{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/":): L_{1,2}\left(N_{1,2}\right)_{G_{1}} +L_{2,2}\left(N_{1,2}\right)_{G_{2}} +L_{3,2}\left(N_{1,2}\right)_{G_{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 L_{1,1}\left(N_{2,1}\right)_{G_{1}} +L_{2,1}\left(N_{2,1}\right)_{G_{2}} +L_{3,1}\left(N_{2,1}\right)_{G_{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/":): L_{1,2}\left(N_{2,2}\right)_{G_{1}} +L_{2,2}\left(N_{2,2}\right)_{G_{2}} +L_{3,2}\left(N_{2,2}\right)_{G_{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 L_{1,1}\left(N_{3,1}\right)_{G_{1}} +L_{2,1}\left(N_{3,1}\right)_{G_{2}} +L_{3,1}\left(N_{3,1}\right)_{G_{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/":): L_{1,2}\left(N_{3,2}\right)_{G_{1}} +L_{2,2}\left(N_{3,2}\right)_{G_{2}} +L_{3,2}\left(N_{3,2}\right)_{G_{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 L_{1,1}\left(N_{4,1}\right)_{G_{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/":): L_{1,2}\left(N_{4,2}\right)_{G_{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 L_{2,1}\left(N_{5,1}\right)_{G_{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/":): L_{2,2}\left(N_{5,2}\right)_{G_{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 L_{3,1}\left(N_{6,1}\right)_{G_{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/":): L_{3,2}\left(N_{6,2}\right)_{G_{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 L_{1,2}\left(N_{1,1}\right)_{G_{1}}+L_{1,1}\left(N_{1,2}\right)_{G_{1}} +L_{2,2}\left(N_{1,1}\right)_{G_{2}}+L_{2,1}\left(N_{1,2}\right)_{G_{2}} +L_{3,2}\left(N_{1,1}\right)_{G_{3}}+L_{3,1}\left(N_{1,2}\right)_{G_{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 L_{1,2}\left(N_{2,1}\right)_{G_{1}}+L_{1,1}\left(N_{2,2}\right)_{G_{1}} +L_{2,2}\left(N_{2,1}\right)_{G_{2}}+L_{2,1}\left(N_{2,2}\right)_{G_{2}} +L_{3,2}\left(N_{2,1}\right)_{G_{3}}+L_{3,1}\left(N_{2,2}\right)_{G_{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 L_{1,2}\left(N_{3,1}\right)_{G_{1}}+L_{1,1}\left(N_{j,3}\right)_{G_{1}} +L_{2,2}\left(N_{3,1}\right)_{G_{2}}+L_{2,1}\left(N_{j,3}\right)_{G_{2}} +L_{3,2}\left(N_{3,1}\right)_{G_{3}}+L_{3,1}\left(N_{j,3}\right)_{G_{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 L_{1,2}\left(N_{4,1}\right)_{G_{1}}+L_{1,1}\left(N_{4,3}\right)_{G_{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 L_{2,2}\left(N_{5,1}\right)_{G_{2}}+L_{2,1}\left(N_{5,3}\right)_{G_{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 L_{3,2}\left(N_{6,1}\right)_{G_{3}}+L_{3,1}\left(N_{6,3}\right)_{G_{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/":): -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(L_{1,1}\rho _{11}^{1}+L_{1,2}\rho _{11}^{2}\right)} | 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/":): \left(L_{1,1}\rho _{22}^{1}+L_{i,2}\rho _{22}^{2}\right) | 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/":): \left(L_{1,1}\rho _{12}^{1}+L_{1,2}\rho _{12}^{2}\right) |
| 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(L_{2,1}\rho _{11}^{1}+L_{2,2}\rho _{11}^{2}\right)} | 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/":): \left(L_{2,1}\rho _{22}^{1}+L_{2,2}\rho _{22}^{2}\right) | 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/":): \left(L_{2,1}\rho _{12}^{1}+L_{2,2}\rho _{12}^{2}\right) |
| 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(L_{3,1}\rho _{11}^{1}+L_{3,2}\rho _{11}^{2}\right)} | 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/":): \left(L_{3,1}\rho _{22}^{1}+L_{3,2}\rho _{22}^{2}\right) | 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/":): \left(L_{3,1}\rho _{12}^{1}+L_{3,2}\rho _{12}^{2}\right) |
| 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} | 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/":): 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/":): 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 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/":): 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/":): 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 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/":): 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/":): 0
|
In this last expression 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,j} =L_{i,j}^{M}}
Document information
Published on 01/01/2005
DOI: 10.1016/j.cma.2004.07.039
Licence: CC BY-NC-SA license
Share this document
claim authorship
Are you one of the authors of this document?