The geometric non-linear shell elements which are currently available in XFINAS is based on the Updated Lagrangian method. The formulation of these elements uses Mindlin-Reissner theory, which the formulation of the linear and geometrical stiffness is exactly defined by incorporating membrane forces, bending moment and transverse shear resultant force. In order to remove rigid body rotations, the co-rotational method is used. The material is assumed to be isotropic, laminated composite and elasto-plastic material. Plasticity is handled by applying the von Mises yield condition and Prandl-Reuss flow rule to discrete points through the shell wall. By using six degrees of freedom per node, the present element can model stiffened plate and shell structures.

The transverse shear stiffness in the laminated composite materials is defined by an equilibrium approach instead of using the shear correction factor. Four macroscopic progressive failure criteria for fiber-reinforced composite under loading were implemented. These criteria may be used to check for first ply failure of composite structures or in determining ultimate loads for complete laminate failure and getting risk assessment.


An assumed strain quasi-conforming shell element

A new Quasi-conforming formulations of 4-node stress resultant shell element ( XSHELL41 ) is developed and implemented in the XFINAS for the solution of stability problems of stiffened plates and shells. The shear-locking behaviour is eliminated by using the Quasi-conforming method. The stiffness matrices for the present elements are explicitly expressed and the stresses are taken accurately at the nodal points. Compared to elements using Gauss integration, where the stresses are most accurate at the integration points, the extrapolation procedure needed for post processing is eliminated in the present shell element.

 

< Fig. 2 Mid-Surface Geometry and Local Coordinates of 4 Node Shell >

 


An assumed natural strain shell element

A 4-nodes assumed strain shell element ( XSHELL42 ) and 8 nodes assumed strain shell element ( XSHELL82 ) are developed and implemented in XFINAS for the solution of stability problems of stiffened plates and shells. The element is free of both membrane and shear locking behaviour by using the assumed strain method such that the element performs very well in the thin shells.

 

< Fig. 3 Mid-Surface Geometry and Local Coordinates of 8 Node Shell >

 


Example Linear analysis of Pinched cylinder (Fig. 4) with rigid end diaphragms using 4-node quasi-conforming shell element is carried out. This example is the most severe benchmark test for 4-node shell element. In addition large deformation analysis of the pinched elastoplastic cylinder with strain hardening using the same examples is carried out.

 

Mesh

S4R5
(ABAQUS)

MITC4
(ADINA)

QPH
(DYNA3D)

XSHELL41
(XFINAS)

4x4

-

0.37

0.37

0.629

6X6

0.602

-

-

-

8x8

-

0.74

0.74

0.929

11X11

0.875

-

-

-

16x16

-

0.93

0.93

0.997

21X21

0.974

-

-

-

< Fig. 4 Pinched Cylinder >

< Table Normalized results of pinched cylinder >

 

< Fig. 5 Load-deflection curve of pinched cylinder >

 

 

 

<point a in Fig. 5>

<point b in Fig. 5>

 

 

 

 

point c in Fig. 5

point d in Fig. 5

 

 

 

point e in Fig. 5

 

< Fig. 6 Deformed shape of elastoplastic cylinder >

 

 

Example Large deformation of elasto-plastic analysis of Agelidis's imperfect stiffened shell under compression is carried out using 8-node shell element.

 

 

 

< Fig. 7 Stiffened shell under compression >

< Fig. 8 Load-deflection curve >

 

< Fig. 9 Preprocessing of stiffened panel >