Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 276304, 16 pages

http://dx.doi.org/10.1155/2013/276304

Research Article

An 8-Node Shell Element for Nonlinear Analysis of Shells Using the Refined Combination of Membrane and Shear Interpolation Functions

Woo-Young Jung¹ and Sung-Cheon Han²

¹Department of Civil Engineering, Gangneung-Wonju National University, 7 Jukheon, Gangneung 210-702, Republic of Korea

²Department of Civil & Railroad Engineering, Daewon University College, 599 Shinwol, Jecheon 390-702, Republic of Korea

Received 6 April 2013; Accepted 19 June 2013

Academic Editor: Hung Nguyen-Xuan

Copyright © 2013 Woo-Young Jung and Sung-Cheon Han. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

An improved 8-node shell finite element applicable for the geometrically linear and nonlinear analyses of plates and shells is presented. Based on previous first-order shear deformation theory, the finite element model is further improved by the combined use of assumed natural strains and different sets of collocation points for the interpolation of the different strain components. The influence of the shell element with various conditions such as locations, number of enhanced membranes, and shear interpolation is also identified. By using assumed natural strain method with proper interpolation functions, the present shell element generates neither membrane nor shear locking behavior even when full integration is used in the formulation. Furthermore, to characterize the efficiency of these modifications of the 8-node shell finite elements, numerical studies are carried out for the geometrically linear and non-linear analysis of plates and shells. In comparison to some other shell elements, numerical examples for the methodology indicate that the modified element described locking-free behavior and better performance. More specifically, the numerical examples of annular plate presented herein show good validity, efficiency, and accuracy to the developed nonlinear shell element.

1. Introduction

The 8-node isoparametric serendipity shell finite elements were suffering from locking effects due to smaller thickness. To avoid this deficiency, reduced and selective integration techniques in the finite element method have been proposed; however, spurious zero-energy kinematic modes occur and may disturb the finite element response in a mesh.

Hinton and Huang [1] developed the 8-node element model, but the element showed poor accuracy rather than the 9-node quadrilateral element. Lakshminarayana and Kailash [2] presented an 8-node shell element with locking free. In order to describe the locking problem, they used appropriately chosen interpolation functions based on Hinton and Huang’s concept.

Bathe and Dvorkin [3] proposed an 8-node shell element-MITC8 to avoid membrane and shear locking problem. The strain tensor was expressed in terms of the covariant components and contravariant base vectors. The performance of this element was quite satisfying and suggested as the promising results in very complex shell structures. Bucalem and Bathe [4] had improved in previous publications the MITC8 shell elements [3], and the results provided conservative and unconservative performance in the finite element method.

Kim and Park [5] and Kim et al. [6] presented an 8-node shell finite element. In an 8-node shell element, the persistence of locking problems was found to continue through numerical experiments on the standard test problem of MacNeal and Harder [7]. In recent year, Han et al. [8] studied the new combination of sampling points for the assumed natural strain and concluded that while it performed quite effectively in some cases, in a few analyses of very thin-walled structures, the accuracy of the element showed less than the predicted results.

Han et al. [8] used the membrane strain sampling points by Bathe et al. [9] and Han et al. [10]. The shear strain sampling points by Hwang [11] are selected in Han et al. [8].

In present study, the interpolation functions by Polit et al. [12] for the sampling points (Han et al., [8]) are used and the strain component of center point is replaced by mean of the components at two points (Bathe and Dvorkin, [3]).

The primary goal of this paper is to propose an improvement of the most useful curved quadrilateral shell finite element, which is clearly, from a practical point of view, the 8-node element and then in order to improve the 8-node ANS shell element, a new combination of sampling points and shaper functions are adopted for the analysis of plates and shells. Also, in order of validation of the present shell element models, the numerical examples are investigated and compared with those solutions from the previous literatures.

2. Improvement of Shell Element

2.1. Kinematics of Shell

Figure 1 shows the geometry and current kinematics of an 8-noded shell element with six degrees of freedom.

Figure 1: Geometry of 8-node shell element with six degrees of freedom.

Kinematic equations for the first-order shear deformation theory including extension of the normal line can be obtained from the 3D equations of the theory of elasticity by using a well-known first-order approximation of a vector function with respect to the coordinate (Rikards et al, [14]). Further well-known geometric relations for the shell in normal coordinates are also used.

Let us assume that vector characterizes the position of an arbitrary point of the shell in the initial reference state (see point in Figure 2) and vector is the position of the same point (point ) in the deformed state. The position of a point at the mid-surface of the shell (point ) in the initial state is characterized by vector , and position of the same point in the deformed state (point ) is characterized by vector . Normal curvilinear coordinates at the mid-surface of the shell in the initial state are defined by the right-handed triad of the base vectors . Here, a unit vector is normal to the mid-surface of shell. Therefore (see Figure 2),

Figure 2: Kinematics of the first-order shear deformation theory.

Similarly, curvilinear coordinates in the deformed state are defined by the triad of vectors . In the deformed state, the vector may not be perpendicular to the mid-surface of shell.

Vector function can be expanded in a Taylor’s series with respect to coordinate normal to the mid-surface of shell as follows:

Here, is a second-order tensor (see Luri’e, [18]), which characterizes the gradient of strains with respect to coordinate , where is a Hamiltonian operator and is a dyadic product of the tensors.

Further, the notations for displacement at the mid-surface and displacement of an arbitrary point of shell are introduced. The following expressions can be written (see Figure 2) as

From expressions (1)–(3), the representation of the displacement of an arbitrary point of the shell for the first-order approximation can be obtained from

where is vector of rotation at the mid-surface of shell

where

The incremental form of the displacement field may be written in terms of the nodal incremental vector as

where and is the rotational matrix of the normal vector at node

where is the transformation matrix between the initial shell normal and the deformed normal can be written as the result of a sequence of finite rotations , , and as follows:

The incremental membrane and bending and transverse shear strains can be separated into linear and nonlinear parts (see Han et al., [10]). The linear strain parts can be expressed as:

The shell theory outlined above is the so-called first-order shear deformation theory with six degrees of freedom.

2.2. Strain Energy and Stiffness of Shell

The strain energy of the shell represented as a three-dimensional body is given by the expression, where in curvilinear coordinates, the contravariant stress tensor is contracted with the covariant strain tensor

where is a volume element of the shell.

The Hooke’s law for each layer can be written as

Here, is the component of the elasticity tensor in the mid-surface metrics. Substituting (12) into (11) as the strain energy can be expressed as

The plane stress condition is enforced. For single layer plates and shells, these equations are quite straightforward. The membrane forces, the bending moments, and the transverse shear forces can be obtained by integrating the relevant stresses through the thickness using the equivalent constitutive equations:

where , and .

The membrane-bending and shear strain rigidity have and matrices, respectively

where is the elastic constitutive coefficient

where is the shear correction factor and is the shear modulus. When the entire cross section is isotropic elastic, then the rigidity matrix reduces to the following form:

where and are the elastic constitutive coefficients of membrane-bending and transverse shear parts.

2.3. Various Enhanced Strain Interpolation

Many researchers addressed 8-node shell finite elements (Hinton and Huang [1]; Lakshminarayana and Kailash [2]; MacNeal and Harder [19]; Bucalem and Bathe [4]; Kim et al. [5, 6]; Han et al. [8]). In this study, for the improved efficient 8 node shell element, the usual 8-nodes of Lagrangian displacement interpolations are employed, and the various combinations of assumed natural strain interpolation functions are used. Figure 3 lists various patterns of sampling points that can be employed for membrane and in-plane shear and transverse shear strain interpolations for the improved 8-node shell element. Based on Figure 3, the pattern is used for membrane ( and ), and the pattern is used for membrane (), as well as transverse shear (). The pattern and pattern are used for in-plane and transverse shear, respectively.

(a) Pattern (following Bathe et al. [9]) 

(b) Pattern (following Han et al. [10])

(c) Pattern (following Han et al. [10])

(d) Patterns , (following Hwang [11])

(e) Patterns (following Polit et al. [12])

Figure 3: Sampling points of the 8-node ANS shell element.

Han et al. [8] used the transverse shear interpolation functions by Huang [11] in the patterns. The interpolation functions by Polit et al. [12] are used in the and patterns. In the and patterns, the strain component of center point is replaced by the mean of the components at points and (Bathe and Dvorkin [3]). The six cases of the combinations of various sampling points are used in the analysis.

3. Incremental Equation of Equilibrium

The generalized Hook’s law at large strain does not represent an approximate material behavior description because stress-strain relation is non-linear. From the practical point of view, Hook’s law is only applicable to small strain, for which constitutive tensor is constant coefficient.

The following incremental equilibrium equation is obtained

where superscript which is generally used as the current configuration is ignored in (18) and superscript is the adjust incremented configuration; is the external virtual work in .

The total tangent stiffness comprises the material stiffness and the geometric stiffness. The linear part of the Green strain tensor is used to derive the material stiffness matrix, and the non-linear part of the Green strain tensor is used to derive the geometric stiffness matrix.

3.1. Linear Element Stiffness Matrix

If the strain displacement (8) is substituted into (18), the linearized element material stiffness matrix is obtained as follows:

Finally, the element stiffness matrix has size on the reference surface of shell element

where the submatrix of is shown in Han et al. [10].

3.2. Geometric Stiffness Matrix

In order to obtain an accurate geometric stiffness matrix, the stresses should be evaluated accurately. The accuracy of the computation of stresses for formulation of geometric stiffness matrix is maintained by obtaining the same interpolated strains in the computation of linear stiffness matrix. The stresses are computed at the integration points based on these strains. Substituting the non-linear part of strain into (18), the following geometric stiffness matrix () is obtained:

The geometric stiffness matrix in the natural coordinate is analytically integrated through the thickness. By the transformation of the natural to the global frame, the element geometric stiffness matrix is obtained on the global frame with submatrix as follows:

where, the sub-matrix of is shown in Han et al. [10]. Then the final assembled incremental non-linear equilibrium equation can be written as

where, and are the external and internal forces, respectively.

The equilibrium equation must be satisfied throughout the complete history of loading, and the non-linear processing will be stopped only when the out of balance forces are negligible within a certain convergence limit. If it is necessary to extend the stability analysis beyond the limit point, that is, in the so-called post-buckling range, appropriate solution procedures must be applied. One approach is to use the arc-length control method in conjunction with the Newton-Raphson method to extend the stability analysis beyond the limit point, by Crisfield [21].

3.3. Torsional Effect

The element stiffness matrix may be written in a matrix form using the equivalent constitutive equations. Finally, the element stiffness matrix has size on the reference surface of shell element. The torsional stiffness term was formed as described by Kanok-Nukulchai [22] and added to the stiffness term.

In this study, based on the procedure proposed by Kanok-Nukulchai [22], the drilling degree of freedom will be tied to the in-plane twist by a penalty functional through additional artificial strain energy as

where is a parameter to be determined (the value of 0.1 suggested); is the shear modulus; is the volume of the element; and is the volume element. After integration throughout the thickness, (24) can be written as

If is chosen to be large relative to the factor (which appears in the bending energy), (25) will play the role of penalty function and result in the desired constraint:

At the Gauss points, A two-by-two Gauss integration scheme is applied for the evaluation of the torsional stiffness in order to avoid the overconstrained situation. To derive a torsional stiffness from (25), the local variables are expressed in terms of global nodal variables, , by shape functions. This gives (25) in the form (Kanok-Nukulchai [22])

The torsional stiffness term is added as described by Kanok-Nukulchai [22].

4. Numerical Results of the 8-Node Shell Element

In this section, we describe numerical tests of improved 8-node shell elements listed in Table 1. They are based on the various combinations of assumed natural strain interpolation functions. Several numerical examples are solved to conduct the performance of the shell element in linear and non-linear applications. The patch test, distortion test, and other various numerical tests of the present shell elements are carried out and validated using FEAP program (Zienkiewicz and Taylor [23, 24]). All the results of the shell element showed very good agreement with references. Several examples are given to demonstrate the efficiency and accuracy of the present shell element.

Table 1: List of shell elements used for comparison.

Name
Description (Patterns referred to Figure 2)

Present ()
Sampling points for membrane (Pattern ), in-plane shear (Pattern ), and transverse shear (Pattern )

Present ()
Sampling points for membrane (Pattern ), in-plane shear (Pattern ), and transverse shear (Pattern )

Present ()
Sampling points for membrane (Pattern ), in-plane shear (Pattern ), and transverse shear (Pattern )

Present ()
Sampling points for membrane (Pattern ), in-plane shear (Pattern ), and transverse shear (Pattern )

Present ()
Sampling points for membrane (Pattern ), in-plane shear (Pattern ), and transverse shear (Pattern )

Present ()
Sampling points for membrane (Pattern ), in-plane shear (Pattern ), and transverse shear (Pattern )

QUAD8
8-node shell element (MacNeal and Harder [7])

QUAD8*
8-node shell element (STRAND 7 [15])

QUAD8**
8-node ANS shell element (Lakshminarayana and Kailash [2])

XSHELL41
4-node quasiconforming shell element (Kim et al. [16])

XSHELL-8-ANS
8-node assumed natural strain shell element (XFINAS [17])

MITC8
8-node shell element using a mixed interpolation of tensorial components (Bathe and Dvorkin [3])

4.1. Patch Test

In the study, the basic patch tests proposed by Simo et al. [13] were performed, and the results illustrated in Figure 4. A patch of five elements is used in this study. Also, Boundary and loading conditions are illustrated in Figures 5 and 6, respectively. Tables 2, 3, and 4 present the normalized solutions of nodal displacements on the right edges. The normalized solutions are presented in the nondimensional form using the following equation:

Table 2: Results of patch test under bending (reference solution: ).

Combinations of various sampling points

Normalized solutions
1.000
1.000
1.000
1.000
1.000
1.000

Table 3: Results of patch test under out-of-plane shear (reference solution: ).

Combinations of various sampling points

Normalized solutions
0.053
0.989
0.989
0.989
0.989
1.000

Table 4: Results of patch test under in-plane tension (reference solution: ).

Combinations of various sampling points

Normalized solutions
1.000
0.965
0.992
0.992
0.992
0.992

Figure 4: Mesh for patch test (Simo et al., [13]). Length of the square ; Young’s modulus ; Poisson’s ratio ; and thickness .

(a) Bending test  

(b) Transverse shear test

(c) In-plane tension test

Figure 5: Boundary displacement conditions for patch tests.

(a) Bending test

(b) Transverse shear test

(c) In-plane tension test

Figure 6: Loading conditions for patch tests.

Patch test results indicate that various sampling points for assumed natural strain method in Tables 2 to 4 can represent fields of constant moment, constant in-plane tension, and transverse shear forces. In some cases, the present shell element can pass the patch test when patterns , , , and are used in all cases. In other cases, the shell element cannot pass the patch test when pattern is used in pure transverse shear case. Therefore, we found that it is very important to apply the best combination of sampling points for 8-node shell elements.

4.2. Bending of Rectangular Plate

The simply supported and clamped rectangular plate problem under uniform and central point loadings is applied to test shear locking by changing aspect ratios. Two aspect ratios of and 5 were considered, and a quarter was modeled due to symmetry. To test the general applicability, the plate is analyzed using both a rectangular mesh and a distorted mesh. Reference solutions (MacNeal and Harder [7]) are given in Tables 5–10.

Table 5: Results of rectangular plate (; reference solution: ).

Element size
Normalized solutions

QUAD8
QUAD8**

1.009
1.000
1.000
1.000
1.016

Table 6: Results of rectangular plate (; reference solution: ).

Element size
Normalized solutions

QUAD8
QUAD8**

1.057
1.001
1.001
0.997
1.090

Table 7: Results of rectangular plate (; reference solution: ).

Element size
Normalized solutions

QUAD8
QUAD8**

1.009
1.000
1.000
1.000
1.010

Table 8: Results of rectangular plate (; reference solution: ).

Element size
Normalized solutions

QUAD8
QUAD8**

0.845
0.845
0.897
0.975
0.867

1.029
0.991
0.991
—
—

Table 9: Results of simply supported plate with rectangular and distorted meshes (reference solutions: and ).

Mesh types
Element
Uniform load
Concentrated load

Rectangular mesh

1.009
1.022

1.000
0.999

Distorted mesh

1.010
1.016

1.000
0.999

Table 10: Results of clamped plate with rectangular and distorted meshes (reference solutions: and ).

Mesh types
Elements
Uniform load
Concentrated load

Rectangular mesh
Choi et al*
0.965
0.990

XSHELL41
0.997
0.992

1.052
1.057

1.002
1.001

Distorted mesh
Choi et al.*
0.976
1.011

XSHELL41
1.023
0.995

1.049
1.039

0.999
0.998

Results of Choi et al. [20] and XSHELL41 are referenced from Kim et al. [16].

4.2.1. Rectangular Mesh

(a) For the Case of  .

The reference vertical deflection at the center of the simply supported plate (shown in Figure 7) under uniform load is 4.062, and the clamped plate under concentrated load is 5.60. In Tables 5 and 6, the numerical results obtained by using different types of existing elements are listed as follows.

(1)Simply supported plate—uniform load.(2)Clamped plate—concentrated load.

Figure 7: Simply supported and clamped rectangular plate.

(b) For the Case of  .

The reference deflection at the center of the simply supported plate under uniform load is 12.97, and the clamped plate under concentrated load is 7.23. In Tables 7 and 8, the numerical results obtained by using different types of elements are listed as follows.

(1)Simply supported plate—uniform load.(2)Clamped plate—concentrated load.

Locking problems occurred when the pattern is applied to transverse shear strain in thin plates (Han et al. [8]). This particular example showed the limit of 8-node shell element for the application of pattern in transverse shear strain. A specifically designed application must be developed in order to remove these errors in the 8-node shell element because other patterns, and , still have the errors. Therefore, in this study, the improved patterns are used and the case of shows more correct results.

4.2.2. Distorted Mesh

In order to test the general applicability of the present formulation for bending, a thin square plate, , is analyzed using both a rectangular mesh and a distorted mesh. A clamped boundary condition is chosen because it is considered to be more severe compared with simple supports. Two types of loading are considered, concentrated load and distributed load. In Figure 8, the problem geometry and material properties are provided. A mesh is used on one-quarter of the plate. The results are compared in Table 9 to those provided by Timoshenko and Woinowosky-Krieger [25]. The accuracy of the results obtained for both cases is maintained and compared with Choi et al. [20] and Kim et al. [16].

Figure 8: Distorted mesh of simply supported and clamped rectangular plate.

The reference vertical deflection at the center of the plate is the simply supported plate under uniform load which is 4.062 and the clamped plate under concentrated load which is 5.60. In Tables 9 and 10, the numerical results obtained by using different types of existing elements are listed as follows.

(1)Simply supported plate—uniform and concentrated load.(2)Clamped plate—uniform and concentrated load.

The cases of and are compared for the test of distorted mesh. It is also found that the locking phenomenon does not happen and the solutions with rectangular mesh are of similar accuracy.

4.3. Curved Beam Problem

The cantilever beam is clamped at one end and loaded by a unit force at the other. The force is applied either along the in-plane axis or along the out-of-plane axis as seen in Figure 9. In the curved cantilever beam, combination of the principal deformation modes is evoked by a single in-plane or out-of-plane shear load at the free end. Note also that the element shape is quite rectangular, which will test the effect of slight irregularity. Theoretical displacement of free end under in-plane load is proposed 0.08734 by MacNeal and Harder [7]. Recently, a different solution of 0.08854 was given by Young. [26]. Reference solutions (Young [26] and MacNeal and Harder [7]) are given in Tables 11-12.

Table 11: Results of curved beam () under in-plane shear (reference solution: ).

Element size
Normalized solutions

QUAD8*
QUAD8**

0.991
0.991
0.991
0.869
0.942

Table 12: Results of curved beam () under out-of-plane shear (reference solution: ).

Element size
Normalized solutions

QUAD8*
QUAD8**

0.970
0.972
0.972
0.958
0.962

Figure 9: Curved Cantilever Beam.

Results presented in Tables 11 and 12 show the outstanding performance of the proposed element.

4.4. Twisted Beam Problem

The initial geometry of the beam is twisted, but the initial strain is equal to zero. The twisted beam problem (shown in Figure 10), which was introduced by MacNeal and Harder [7], was proposed to test the effect of element warping. The performance of elements of thickness 0.32 was investigated under in-plane and out-of-plane shear loads. The warp of each element is . Numerical results in Tables 13 and 14 are listed with displacements of free end normalized to the reference solution (MacNeal and Harder [7]).

Table 13: Results of twisted beam () under in-plane shear (reference solution: ).

Element size
Normalized solutions

QUAD8
QUAD8**

0.997
0.997
0.993
0.998
0.998

Table 14: Results of twisted beam () under out-of-plane shear (reference solution: ).

Element size
Normalized solutions

QUAD8
QUAD8**

1.001
1.004
0.997
0.998
0.998

Figure 10: Twisted cantilever beam.

Results are presented in Tables 13 and 14, showing more correct results than the result of references in spite of the mesh.

4.5. Pinched Hemispherical Shell with 18° Hole

The hemispherical shell with an 18° hole is loaded by two pairs of equal but opposite external forces; see Figure 11. The shell undergoes an almost inextensional deformation, and a membrane locking can destroy the solution. For the shown FE mesh, the elements are flat and trapezoidal. Only a quarter of the shell is modeled because of the double symmetry.

Figure 11: Definition of loading points and initial configuration of hemispherical shell.

For this problem, the reference deflection at point A in the -direction is 0.093 (Simo et al. [13]). The geometry and material properties are shown in Figure 11. The numerical results are given in Table 15.

Table 15: Results of hemispherical shell with 18° hole (reference solution: ).

Element size
Normalized solutions

QUAD8
QUAD8**

0.542
0.543
0.516
0.832
0.387

0.915
0.915
0.902
1.003
0.773

0.984
0.984
0.978
C*
0.950

The grade C means that the QUAD8 gives accuracy 20% > error > 10% (MacNeal and Harder [7]).

This indicates that the results obtained by our elements are more accurate than the solution by the QUAD8** element. Note that the error of mesh for the QUAD8 element is greater than 10%.

4.6. Pinched Cylinder with End Diaphragm

This problem, shown in Figure 12, is yet another challenging test for the ability of the element in modeling inextensional bending modes and complex membrane states. Most of the 4-node shell elements do not show good behavior in this test, but the present 8-node solid element and others passed this test. The vertical deflection of point A was investigated. The reference value is (Belytschko et al. [27]). The results with different meshes and elements are presented in Table 16. A good performance of the current elements was again observed.

Table 16: Results of pinched cylinder with end diaphragm (reference solution: ).

Element size
Normalized solutions

XSHELL-8-ANS
MITC8

0.687
0.720
0.781
0.738
0.833 ()

1.006
1.024
0.964
0.945
0.990

Figure 12: Problem of pinched cylinder with end diaphragms.

The cases of and show that the results are more conservative than the reference solution. However, the solution in the case of is more accurate than that in the reference element XSHELL-8-ANS.

4.7. Nonlinear Analysis of Slit Annular Plate

To validate the present shell element for geometrically nonlinear analysis, a numerical example of annular plate under end shear force is solved. We examine an annular clamped plate subjected to a distributed transverse shear force (see Figure 13). This good benchmark problem was considered by Buechter and Ramm [28], Brank et al. [29], and Sansour and Kollmann [30], among others, and Sze et al. [31], presented ABAQUS results for the problem. In particular, we use the ABAQUS results presented in a recent paper of Sze et al. [31] for popular benchmark problems of nonlinear shell analysis (because of the tabulated data). In Figure 13, the geometry and elastic material properties for the isotropic case are given.

Figure 13: The slit annular plate loaded with various end forces .

The line force is applied at one end of the slit, while the other end of the slit is fully clamped. In the present example, the plate is modeled by polar coordinates. The problem, however, can be solved by using the standard Cartesian coordinates as well. A regular mesh of elements is considered in the present analysis. Solutions obtained by the mesh of are too stiff and suffer from locking. The arc-length method is used to obtain the load-deflection curve in these examples. The shear load versus displacement curves for three characteristic directions are depicted in Figure 14. Solutions obtained with the present formulation are in good agreement with those obtained by Sze et al. [31]. Figure 15 shows the deformed configuration of the isotropic annular plate for .

Figure 14: Load-deflection curves for the slit annular plate lifted by line force .

Figure 15: The deformed shape of slip annular plate at .

4.8. Nonlinear Analysis of Hinged Shell

This is a widely used example for testing geometrically nonlinear shell elements. In this study, a quadrant of the shell is modeled with mesh sizes, as shown in Figure 16. Three different thickness values are considered; that is, , 6.35, and 12.7 mm, respectively. The problem has been considered in Sabir and Lock [32], Horrigmoe and Bergan [33], Rhiu and Lee [34], Sze et al. [31], and Arciniega and Reddy [35]. The arc-length control method is used to trace the equilibrium path, and a tolerance coefficient of is used. The total length is  mm, curvature of the circular edges is 1/2,540 mm, and the circumferential length is 508 mm. The material properties are Young’s modulus , and Poisson’s ratio .

Figure 16: Hinged cylindrical shell; curved edges are free, and straight edges are hinged .

For the case of 12.7 mm, the curve of the central deflection versus load is given in Figure 17, showing the good correlation between the present solution and the existing solutions.

Figure 17: Load-deflection curves for the hinged shell ().

For the case of 6.35 mm, the solution derived with four-node shell elements by Sze et al. [31] is used as the major reference. The load-deflection curve is shown in Figure 18. It is seen that the solution in this study is very close to the solution of Sze et al. [31] in Figures 17 and 18.

Figure 18: Load-deflection curves for the hinged shell ().

In the third case of the present analysis, even a further reduction of the panel thickness was presumed taking 5 mm. The load-deflection curve of very thin panel was presented in Figure 19. One can remark that the decrease of the thickness of the shell has reduced its stiffness; nevertheless, the range of rotations has not exceeded the limits of small rotations.

Figure 19: Load-deflection curves for the hinged shell ().

4.9. Nonlinear Analysis of Open-Ended Cylindrical Shell

The large deformation analysis of a stretched cylinder, as shown in Figure 20, is carried out with free ends. It is subjected to a pair of concentrated forces. One octant of the cylinder is analyzed with mesh sizes. Material properties are and . The cylinder length is 10.35, the radius is 4.953, and the thickness is 0.094.

Figure 20: Problem definition and an initial configuration. .

The radial displacements at points A, B, C, and D of the shell are shown in Figure 21. It is noticed that the load-deflection curves of the shell have two different regions: the first is dominated by bending stiffness with large displacements; the second, at load level of , is characterized by a very stiff response of the shell. The present converged solutions are compared with the results obtained by Sze et al. [31], who used the commercial code ABAQUS. The present results agree well with those obtained by Sze et al. [31] for the three curves considered. The deformed shape of the cylindrical shell, at load level of , is presented in Figure 22.

Figure 21: Load-deflection curves for the open-ended cylindrical shell.

Figure 22: The deformed shape of open-ended shell at .

4.10. Nonlinear Analysis of Pinched Cylindrical Shell

Here, the cylindrical shell which is clamped at one end and subjected to a pair of concentrated loads at the other free end is studied. This example provides a severe test for finite element formulations. This example was considered in Brank et al. [29], Park et al. [36], Sansour and Kollmann [30], and Arciniega and Reddy [35] among others.

The two forces act in opposite directions as shown in Figure 23. Material properties are Young’s modulus, , Poisson ratio, . The cylinder length is , and the radius with thickness . One-quarter of the structure is analyzed with element mesh. In Figure 24, the radial displacements at points A and B of the shell are shown and compared with the reference solution of Sze et al. [31]. Figure 24 shows that the load-deflection curve goes beyond the highest physically possible deflection of the loading point. The main purpose of the presented example was to test the behavior of the numerical model, although the authors are aware that the result of a real structure can differ from the presented solution due to a possible contact between two deformed regions of the shell in the postbuckling region. The deformed shape of the semicylindrical shell subjected to an end pinching load, at load level of 0.35, is presented in Figure 25.

Figure 23: Problem definition and an initial configuration .

Figure 24: Load-deflection curves for the pinched cylindrical shell.

Figure 25: The deformed shape of pinched cylindrical shell at 0.35.

5. Concluding Remarks

In this study we have improved an 8-node ANS finite element model for linear and non-linear analysis of plate and shell structures. The accuracy and robustness of the improved formulation and its implementation are illustrated via a variety of shell geometries and deformations of shells.

The objective of this paper is to present some results using the 8-node shell element when the sampling points for the strain components are changed. The 8-node ANS shell element using the sampling points earlier proposed (Bathe and Dvorkin [3], Lakshminarayana and Kailash [2], Bucalem and Bathe [4], Kim et al. [6], and Han et al., [8]) is considerably a powerful element and has shown good convergence behavior in many problem solutions. However, the research paper of predicted strain components show that the proper combination of sampling points is not predicted with excellent accuracy or reveals persistence of locking problems.

In order to improve the 8-node ANS shell element, a new combination of sampling points is adopted for linear analysis of isotropic and laminated composite structures. The present assumed strain method completely removes both membrane and shear locking behavior even when full integration is used in the formulation. A new combination of the sampling points for the in-plane normal, in-plane shear, and transverse shear strain components results in significantly better results. Optimality in the convergence behavior is retained, and all strain components are predicted with reasonable accuracy. Other combinations of sampling points are also considered but do not result into a further improvement in the predictive capability of the element. From several numerical examples, the new combination 8-node ANS shell element shows better performance compared with other combination shell elements. The present combination of sampling points () could be easily implemented into finite element code and used for the practical purpose. In particular, numerical examples for annular plate presented herein clearly show the validity of the present approach and the accuracy of the developed nonlinear shell element. Future work will be useful to extend this work to dynamic and buckling analysis of isotropic and laminated composite shell structures.

Acknowledgment

This work was supported by the National Research Foundation of Korea grant funded (NRF) by the Korean Government (MEST) (no. 2012-0008762).

References

E. Hinton and H. C. Huang, “A family of quadrilateral Mindlin plate elements with substitute shear strain fields,” Computers and Structures, vol. 23, no. 3, pp. 409–431, 1986.

H. V. Lakshminarayana and K. Kailash, “A shear deformable curved shell element of quadrilateral shape,” Computers and Structures, vol. 33, no. 4, pp. 987–1001, 1989.

K.-J. Bathe and E. N. Dvorkin, “Formulation of general shell elements—the use of mixed interpolation of tensorial components,” International Journal for Numerical Methods in Engineering, vol. 22, no. 3, pp. 697–722, 1986.

M. L. Bucalem and K.-J. Bathe, “Higher-order MITC general shell elements,” International Journal for Numerical Methods in Engineering, vol. 36, no. 21, pp. 3729–3754, 1993.

K. D. Kim and T. H. Park, “An 8-node assumed strain element with explicit integration for isotropic and laminated composite shells,” Structural Engineering and Mechanics, vol. 13, no. 4, pp. 387–410, 2002.

K. D. Kim, G. R. Lomboy, and S. C. Han, “A co-rotational 8-node assumed strain shell element for postbuckling analysis of laminated composite plates and shells,” Computational Mechanics, vol. 30, no. 4, pp. 330–342, 2003.

R. H. MacNeal and R. L. Harder, “A proposed standard set of problems to test finite element accuracy,” Finite Elements in Analysis and Design, vol. 1, no. 1, pp. 3–20, 1985.

S.-C. Han, W. Kanok-Nukulchai, and W.-H. Lee, “A refined finite element for first-order plate and shell analysis,” Structural Engineering and Mechanics, vol. 40, no. 2, pp. 191–213, 2011.

K.-J. Bathe, P.-S. Lee, and J.-F. Hiller, “Towards improving the MITC9 shell element,” Computers and Structures, vol. 81, no. 8–11, pp. 477–489, 2003.

S. C. Han, K. D. Kim, and W. Kanok-Nukulchai, “An element-based 9-node resultant shell element for large deformation analysis of laminated composite plates and shells,” Structural Engineering and Mechanics, vol. 18, no. 6, pp. 807–829, 2004.

H.-C. Huang, Static and Dynamic Analyses of Plates and Shells, Springer, London, UK, 1989.

O. Polit, M. Touratier, and P. Lory, “A new eight-node quadrilateral shear-bending plate finite element,” International Journal for Numerical Methods in Engineering, vol. 37, no. 3, pp. 387–411, 1994.

J. C. Simo, D. D. Fox, and M. S. Rifai, “On a stress resultant geometrically exact shell model. II. The linear theory; computational aspects,” Computer Methods in Applied Mechanics and Engineering, vol. 73, no. 1, pp. 53–92, 1989.

R. Rikards, A. Chate, and O. Ozolinsh, “Analysis for buckling and vibrations of composite stiffened shells and plates,” Composite Structures, vol. 51, no. 4, pp. 361–370, 2001.

STRAND 7 Verification Manual, G+D Computing Pty Ltd, Sydney, 2000

K. D. Kim, G. R. Lomboy, and G. Z. Voyiadjis, “A 4-node assumed strain quasi-conforming shell element with 6 degrees of freedom,” International Journal for Numerical Methods in Engineering, vol. 58, no. 14, pp. 2177–2200, 2003.

XFINAS Validation Manual, AIT, Bangkok, 2008

A. I. Luri’e, Theory of Elasticity, Nauka Publishing House, Moscow, Russia, 1970.

R. H. MacNeal and R. L. Harder, “Eight nodes or nine?” International Journal for Numerical Methods in Engineering, vol. 33, no. 5, pp. 1049–1058, 1992.

C.-K. Choi, P.-S. Lee, and Y.-M. Park, “Defect-free 4-node flat shell element: NMS-4F element,” Structural Engineering and Mechanics, vol. 8, no. 2, pp. 207–231, 1999.

M. A. Crisfield, “A fast incremental/iterative solution procedure that handles “snap-through”,” Computers and Structures, vol. 13, no. 1–3, pp. 55–62, 1981.

W. Kanok-Nukulchai, “Simple and efficient finite element for general shell analysis,” International Journal for Numerical Methods in Engineering, vol. 14, no. 2, pp. 179–200, 1979.

O. C. Zienkiewicz and R. L. Taylor, The Finite Element Method, McGraw-Hill, London, UK, 1989.

O. C. Zienkiewicz and R. L. Taylor, The Finite Element Method, Butterworth-Heinemann, London, UK, 2000.

S. P. Timoshenko and S. Woinowosky-Krieger, Theory of Plates and Shells, McGraw-Hill, New York, NY, USA, 1959.

W. C. Young, Roark’s Formulas for Stress and Strain, McGraw-Hill, New York, NY, USA, 6th edition, 1989.

T. Belytschko, H. Stolarski, W. K. Liu, N. Carpenter, and J. S.-J. Ong, “Stress projection for membrane and shear locking in shell finite elements,” Computer Methods in Applied Mechanics and Engineering, vol. 51, no. 1–3, pp. 221–258, 1985.

N. Buechter and E. Ramm, “Shell theory versus degeneration. A comparison in large rotation finite element analysis,” International Journal for Numerical Methods in Engineering, vol. 34, no. 1, pp. 39–59, 1992.

B. Brank, F. B. Damjanić, and D. Perić, “On implementation of a nonlinear four node shell finite element for thin multilayered elastic shells,” Computational Mechanics, vol. 16, no. 5, pp. 341–359, 1995.

C. Sansour and F. G. Kollmann, “Families of 4-node and 9-node finite elements for a finite deformation shell theory. An assessment of hybrid stress, hybrid strain and enhanced strain elements,” Computational Mechanics, vol. 24, no. 6, pp. 435–447, 2000.

K. Y. Sze, X. H. Liu, and S. H. Lo, “Popular benchmark problems for geometric nonlinear analysis of shells,” Finite Elements in Analysis and Design, vol. 40, no. 11, pp. 1551–1569, 2004.

A. B. Sabir and A. C. Lock, “The application of finite elements to the large deflection geometrically non-linear behaviour of cylindrical shells,” in Variational Methods in Engineering, C. A. Brebia and H. Tottenham, Eds., vol. 7/66-7/75, Southampton University Press, 1973.

G. Horrigmoe and P. G. Bergan, “Nonlinear analysis of free-form shells by flat finite elements,” Computer Methods in Applied Mechanics and Engineering, vol. 16, no. 1, pp. 11–35, 1978.

J. J. Rhiu and S. W. Lee, “A nine-node finite element for analysis of geometrically non-linear shells,” International Journal for Numerical Methods in Engineering, vol. 24, pp. 581–604, 1988.

R. A. Arciniega and J. N. Reddy, “Tensor-based finite element formulation for geometrically nonlinear analysis of shell structures,” Computer Methods in Applied Mechanics and Engineering, vol. 196, no. 4–6, pp. 1048–1073, 2007.

H. C. Park, C. Cho, and S. W. Lee, “Efficient assumed strain element model with six DOF per node for geometrically non-linear shells,” International Journal for Numerical Methods in Engineering, vol. 38, no. 24, pp. 4101–4122, 1995.