SciELO - Scientific Electronic Library Online

vol.60 número2A comparative study on three different construction methods of stiffened plates-strength behaviour and ductility characteristicsInfluence of the residual stresses on the strength of steel columns considering advanced analysis with distributed plasticity índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados




Links relacionados


Rem: Revista Escola de Minas

versión impresa ISSN 0370-4467versión On-line ISSN 1807-0353

Rem: Rev. Esc. Minas v.60 n.2 Ouro Preto abr./jun. 2007 



Elastic–plastic analysis of metallic shells at finite strains


Análise elasto–plástica de cascas metálicas por deformações finitas



Eduardo Moraes Barreto CampelloI; Paulo de Mattos PimentaII; Peter WriggersIII

IPolytechnic School at the University of São Paulo, P.O. Box 61548, São Paulo, 05424–970, Brazil. E–mail:
IIPolytechnic School at the University of São Paulo, P.O. Box 61548, São Paulo, 05424–970, Brazil. E–mail:
IIIInstitute of Mechanics and Computational Mechanics, University of Hanover, Appelstrasse 9A, Hanover, 30167, Germany. E–mail: wriggers@ibnm.uni–




The geometrically–exact finite–strain variable–thickness shell model of [1] is reviewed in this paper and extended to the case of metallic elastoplastic shells. Isotropic elasticity and von Mises yield criterion with isotropic hardening are considered. The model is implemented within a triangular finite element and is briefly assessed by means of two numerical examples.

Keywords: Elastoplasticity, Elastic–plastic shells, finite strains, finite rotations, variable–thickness, triangular shell finite element.


Nesse trabalho, a teoria de cascas geometricamente exata, que considera a variação da espessura apresentada em [1], é estendida para o caso de materiais metálicos. Elasticidade finita isótropa, nas deformações logarítmicas, e o critério de plastificação de von Mises são adotados. A teoria é implementada, computacionalmente, com a ajuda do elemento finito triangular T6–3i de [1], cujo desempenho é avaliado por meio de dois exemplos numéricos.

Palavras–chave: Elastoplasticidade, cascas elastoplásticas, deformações finitas, rotações finitas, variação da espessura, elemento finito triangular de casca.



1. Introduction

The main purpose of this work is to extend the geometrically–exact finite–strain variable–thickness shell formulation of [1] to the case of metallic elastoplastic shells. The shell motion is described by a director theory with a standard Reissner–Mindlin kinematical assumption, and the Euler–Rodrigues formula is utilized in order to represent finite rotations in a fully exact manner.

Our approach to elastoplasticity combines the multiplicative decomposition of the deformation gradient with the concept of an instantaneous natural configuration. The elastic domain is defined on the principal axes of the right Cauchy–Green elastic strain tensor, through a logarithmic neo–Hookean material law. A von Mises yield function with isotropic hardening is adopted for constraining the admissible states, and as a result, the classical radial return method (see e.g. [2,3]) may be straightforwardly employed for integration of the stresses.

The model is implemented over the six–node triangular finite element of [1] and is briefly assessed by means of two numerical examples.


2. Shell theory

A flat reference configuration is assumed for the shell mid–surface at the outset. A unit orthogonal system is defined (see Figure 1) so that points in this configuration are described by the vector field

where describes the position of points at the middle surface and is the shell reference director. Z Î H [–hb, ht] is the thickness coordinate.



Let now {e1, e2, e3} be a local orthogonal system on the current configuration, with as depicted in Figure 1. Here is the rotation tensor and q the rotation vector. On this configuration the position of the shell material points is given by

x = z + se3 (2)

where represents the position of points at the middle surface and

describes the position of points along the current director e3. Note that a quadratic distribution for s was assumed. This is the simplest assumption that avoids the so–called Poisson's locking. The scalar parameters and are transversal degrees–of–freedom that depicts the constant and the linear parts of the strains related to the thickness change.

The rotation tensor Q may be expressed by the well–known Euler–Rodrigues formula

Q = I + h1 (q)Q + h2 (q)Q 2 (4)

in which q = ||q|| is the true rotation angle and Q = Skew (q) is the skew–symmetric tensor whose axial vector is q . Here, the following trigonometric functions are introduced for use along the text

The displacement vector for points on the reference mid–surface is u = zz . Components of u and q together with scalars p and q define the 8 degrees–of–freedom of this shell theory:

d = [ u q p q ]T (6)

The deformation gradient F is evaluated from differentiation of (2) with respect to x, and after some algebra one has

with the notation (·),a = ¶) / ¶xa and (·)'= ¶ (·) /¶z for derivatives. Tensor G in (8)4 is

G = I + h2 (q)Q + h3 (q)Q 2 (9)

Linearization of (7) with respect to d yields the virtual deformation gradient:

where dW = dQQT is a skew–symmetric tensor representing the director virtual spin. Virtual strains of (10) are obtained directly from (8)1–2, rendering


Here, Z,a = Skew (z,a) and

G,a = h2 (q) Q,a + h3 (q) QQ,a + Q,aQ ) + h4 (q) ( q . q,a)Q + h5 (q) ( q . q,a)Q2 (13)


with Q,a = Skew (q,a) .

Let now the first Piola–Kirchhoff stress tensor be expressed by P = Q P r , where is its back–rotated counterpart. Vectors are the back–rotated stress vectors acting on a cross–section normal to ei . With the aid of (10), it is not difficult to show that the shell internal virtual work may be written as

in which W Ì Â2 is the shell domain on the reference middle plane and

are 18x1–vectors grouping the cross–sectional resultants and the corresponding strains. In this case,

The external virtual work may be expressed by

in which are the surface and body loadings respectively. Vector groups the generalized external forces, and is given by


From the virtual work theorem one arrives at the following local equilibrium equations

where are the sectional forces at the current configuration. The Fréchet derivative of d W = d Wintd Wext with respect to d leads to the tangent bilinear form


Submatrices Ga of G can be found in [1], while submatrices of D are computed by the chain rule with the aid of the following quantities


3. Tangent Moduli

Let S be the second Piola–Kirchhoff stress tensor and E the Green–Lagrange strain tensor. The fourth–order tensor of tangent moduli associated with the pair {S, E } is defined by

The back–rotated first Piola–Kirchhoff stress tensor is related to S through P r = F rS. If one introduces the following third– and second–order tangent tensors

then it is possible to write

Expressions (26)1–2 can be used in (16), while (26)3–5 can be used in (23) for the computation of D.


4. Constitutive equation

4.1 Logarithmic neo–Hookean material

Let li (i = 1,2,3) be the principal values of U = (FT F)1/2 and V = (FFT)1/2 . Let vi (i = 1, 2, 3) be a principal triad of V. Let ei = ln li be the principal logarithmic strains, and let us define the following strain invariants: J = ln J = e1 + e2 + e3 , where J = l1l2l3, and , with (notice that .). A neo–Hookean isotropic elastic material can be defined by the specific strain energy function below

where k and µ are material parameters. From this expression, the Kirchhoff–Trefftz stress tensor reads

in which Si are the principal values of S. The second Piola–Kirchhoff stress tensor is then given by

where is a fourth–order tensor that performs a congruent transformation on S. The related tensor of tangent moduli , as first derived in [2] (see also [4,5]), can be written as follows

4.2 Von Mises plasticity

Our basic kinematical assumption in describing an elastoplastic deformation is the multiplicative decomposition of the deformation gradient, in the form F = Fe Fp . Here Fe represents the elastic part of the deformation andFp the natural configuration. We assume that Fp = Fp (t) (t = time) describes the plastic deformation, with Fp (0) being the initial natural configuration set as Fp(0) = I .

The elastic and the plastic deformation gradients have the polar decompositions Fe = ReUe = VeRe and Fp = RpUp, where Ue = (FeTFe)1/2 is the right elastic stretching tensor, Ve = ReUeReT = (FeFeT)1/2 is the left elastic stretching tensor, Re is the elastic rotation tensor, Up is the right plastic stretching tensor and Rp is the plastic rotation tensor.

We relate the Kirchhoff–Trefftz stress tensor to the elastic part of the deformation through the isotropic neo–Hookean law of (27). A von Mises yield function with isotropic hardening constrains the elastic domain as below

In this case, Se = ReTSRe is the Green–Naghdi stress tensor and is its deviatoric part, while is the deviatoric part of S. Still in (32), a is the equivalent plastic strain. From the normality law we have

where the superposed dot denotes time differentiation.

Let us consider now the material state at instant ti, at which Fi,and ai are supposed to be known. One thus has

Performing the polar decomposition , the stress tensor at this instant may be computed by where . At instant ti +1 only the deformation gradient Fi +1 is assumed to be known, and the following trial elastic state is defined assuming that the plastic part of the deformation has not changed

With , the stress tensor in this state is obtained via . If Ft = F (St) < 0, then the trial state is an admissible one, i.e. Da and

However, if Ft > 0 , there is an incremental plastic stretching DUp such that

The incremental stress integration algorithm is then formulated by the following constrained minimization problem, which must be solved for

For the von Mises yield function (32) the analytical solution of (38) is given by

which is the well–known radial return method. Hence, Da is the solution of the equation

(in the case of linear hardening an explicit solution is attained). The radial return algorithm is summarized below:

1) Trial step:

2) Radial return algorithm:


5. Numerical examples

The present model is implemented over the six–node flat triangular element introduced by the authors in [1]. The displacements u are approximated with quadratic functions, while a linear interpolation scheme is adopted for the rotation vector q and for the thickness parameters p and q. No special techniques such as ANS, EAS or reduced–integration with hourglass control are necessary and the element is a pure displacement–based one. Either shear or membrane locking are not observed.

5.1 Bending of square plate

A quadratic plate under constant dead load is considered here similarly to [6,7]. The vertical displacements at the edges of the plate are restricted to zero, and symmetry conditions are adopted. Problem data and analysis results are shown in Figure 2.

5.2 Pinched hemisphere

A relatively thick hemispherical shell (see [8,9]) is loaded by two pairs of collinear forces, as in Figure 3. Symmetry conditions are used so that only one quadrant of the structure needs to be modeled. Geometry and material data are also shown in Figure 3, together with the results obtained.


6. Conclusions

The geometrically–exact finite–strain shell model of [1] was extended here to the case of elastic–plastic shells. Good performance of the formulation can be observed via numerical assessment. We believe that simple but reliable triangular elements combined with stable hardening–plasticity models are an excellent tool for the nonlinear elastoplastic analysis of shell structures.


7. Acknowledgements

Fellowship funding from FAPESP (Fundação de Amparo à Pesquisa do Estado de São Paulo), CNPq (Conselho Nacional de Pesquisa) and DFG (Deutsche Forschungsgemeinschat), as well as the material support from IBNM (Institut für Baumechanik und Numerische Mechanik) are gratefully acknowledged by the first two authors in this work.


8. References

[1] PIMENTA, P.M., CAMPELLO, E. M. B. , WRIGGERS, P. A fully nonlinear multi–parameter shell model with thickness variation and a triangular shell finite element. Comput. Mech., 34, p.181–193, 2004.         [ Links ]

[2] PIMENTA, P. M. Finite–deformation soil plasticity on principal axes, Proceedings of the 3rd COMPLAS, CIMNE, Barcelona, 1992.         [ Links ]

[3] SIMO, J. C., HUGHES, T. J. R. Computational Inelasticity, Springer–Verlag, New York, 1998.         [ Links ]

[4] PIMENTA, P. M. A stress integration algorithm for the analysis of elastic–plastic solids by the finite element method. I. Small strain analysis. Technical Report BT/PEF–8913, Department of Structural and Foundation Engineering, Polytechnic School at the University of São Paulo, 1989.         [ Links ]

[5] PIMENTA, P. M. A stress integration algorithm for the analysis of elastic–plastic solids by the finite element method. II. Large deformation analysis. Technical Report BT/PEF–8914, Department of Structural and Foundation Engineering, Polytechnic School at the University of São Paulo, 1989.         [ Links ]

[6] BÜCHTER, N., RAMM, E. E., ROEHL, D. "3–D extension of nonlinear shell formulation based on the EAS concept. Int. J. Numer. Meth. Engrg., 37, p. 2551–2568, 1994.         [ Links ]

[7] RAMM, E. E., MATZENMILLER, A. Consistent linearization in elastoplastic shell analysis. Eng. Comp., 5, p.289–299, 1988.        [ Links ]

[8] EBERLEIN, R., WRIGGERS, P. FE concepts for finite elastoplastic strains and isotropic stress response in shells: theoretical and computational analysis. Comput. Meth. Appl. Mech. Engrg., 171, p. 243–279, 1999.         [ Links ]

[9] SIMO, J. C., KENNEDY, J. G. On a stress resultant geometrically exact shell model. Part V. Nonlinear plasticity: formulation and integration algorithms. Comp. Meth. Appl. Mech. Engrg., 96, p. 133–171, 1992.        [ Links ]



Artigo recebido em 27/11/2006 e aprovado em 29/01/2007.

Creative Commons License Todo el contenido de esta revista, excepto dónde está identificado, está bajo una Licencia Creative Commons