Chapter Four
Chapter Four
Chapter Four
CHAPTER FOUR
FINITE ELEMENT FORMULATION
FOR STIFFENED SHELL
4.1 Introduction
Various shell theories have been formulated by adopting simplified
hypotheses, which render them to specific cases.
The finite element method has been proven to be an efficient tool in
analyzing shell structures with arbitrary shapes, load conditions, and
support conditions.
Three distinct approaches to the finite element representation of
shell structures have been employed:
(a)
(b)
(c)
55
finite elements with different degrees of complexity have been formulated
for deep and shallow shells [8,20,35,36]. The inadequate membrane
representation was eliminated, but the simplicity of the element was lost.
It was necessary to introduce additional degrees of freedom namely first
and second order derivatives, which lead to difficulties at folds, and
boundaries.
An alternative simple element, which was easy to apply, was first
presented in [2]. This element was derived from a three dimensional
element by applying a so-called degeneration process. The classical
Kirchhoff hypothesis was no longer needed in this formulation since both
in-plane and shear strains were accounted for. Although this type of
element seemed very promising when it was introduced, it was only a
short time until the weakness of this model was detected. As the thickness
of the element was reduced, the results obtained worsened due to locking.
An improvement of the model was achieved by reduced and selective
integration to avoid the locking [39,28].
The above difficulties led to improve the behavior of the element
against locking and spurious mode by interpolating the strains at given
sampling points to the Gauss sampling points [5,19] which appear to be
efficient in representing the shell structure.
Finite element formulation for the stiffened shell is employed by
1.curved elements formulated on the basis of shell theories
2.elements derived from three dimensional elements.
For the shallow shell formulations, two types of elements were
used, the shallow shell element and the bar element for the formulation of
the stiffened shell structure.
56
formulation is employed to represent the thick shell the situation is
restricted to specific shapes [7]. It is also important to indicate the
complex formulation and implementation of the shallow shell element for
shell structures.
It has been proved that the degenerated shell formulation is a
simple and general method to represent the stiffened and unstiffened shell
structures. When the degenerated element is used there are two methods
to represent the stiffeners. The first is by using the curved beam element
and the second is by representing the stiffeners by shell elements. When
the analysis is performed better results are obtained when the stiffeners
are represented by the curved beam elements [3,33,38]. This is due to the
crude approximation for the stiffeners when modeled by shell elements
which is in fact nothing more than using five degrees of freedom for a
case which has only two degrees at any section. Accurate results have
been obtained for both the mode and the magnitude when the stiffeners
are represented by a beam element in the finite element analysis [38].
In this chapter, a finite element formulation for stiffened shell is
presented. The modeling of the stiffened shell by the degenerated shell
element connected to the curved beam is adopted in the present work.
Compatibility must exit between the beam and the degenerated shell
element to predict the correct behavior of the stiffened shell [22].
Thannon [38] developed a nonlinear analysis based on the linear work of
Jirousek [22]and was applied to the analysis of reinforced concrete.
Layered approach for the degenerated shell element and for the curved
beam is adopted in the present work. The reinforcement is represented by
the smeared approach and by the embedded bar formulation. The
assumed strain formulation for the degenerated shell element is adopted
in the present work.
57
The consistent mass matrix for both the degenerated shell element
and for the curved beam is adopted. The displacement, velocity and
acceleration compatibility between the curved beam and the degenerated
shell element must be existed. Suitable transformation for stiffness matrix
( or B-matrix) was introduced by Jirousek[22] for linear analysis and then
developed by Thannon[38] through layers adopted in this work and it is
applied to the mass matrix. The damping matrix is obtained by
proportioning of the mass and stiffness matrices which is known as
Rayleigh damping.
4.2 Degenerated Isoparametric Elements
The degenerated shell element has been originally introduced by
Ahmed [2]. In this element Mindlin type theory is employed. The basic
idea of this element is:
(a)
58
definition of independent rotations and displacements degrees of freedom
permits transverse shear deformation to be taken into account, since
rotation are not tied to the slope of the midsurface. This approach is
equivalent to using a general shell theory and reduced to the hypotheses
of Reissner and Mindlin when applied to plates.
4.3 Coordinate System
In the formulation of the degenerated shell elements different coordinate
systems are shown in Figure (4.1):
4.3.1 Global Coordinate Set
This is a Cartesian coordinate system, freely chosen, in which the
geometry is defined in space, as well as in the global stiffness matrix and
applied force vector.
4.3.2 Curvilinear Coordinate System ( , , )
In this system
approximately normal to the middle surface and relates the points through
the thickness to the middle surface.
4.3.3Local Cartesian Coordinate System ( x , y , z )
This local coordinates system is defined at the sampling Gauss points,
wherein the stresses and strains are calculated. At such point, the z
direction is taken to be normal the surface
cross product of the
x x
y y
z z
(4.1)
direction as:
Figure (5.1) Coordinate system (a) nodal and curvilinear (b) local axes
system.
59
x
y
x
z
(4.2)
(4.3)
(4.4)
V3 k
V3 k y k
z
k
v3k
top
xk
yk
z
k
(4.5)
bot
v3xk
3k v3yk
V3 k z
v3k
The vector
v1k
(4.6)
is perpendicular to
v3k
z plane, so that:
v1k j v3 k / j v3 k
Or if
v3k
is in the y-direction,
v1k v3 k i / v3 k i
where
i and
(4.7)
(4.8)
60
(4.9)
v 2 k v1k v 3 k / v1k v3 k
y
N k ( , ) 2 y k
z
z
k
Here
N k ( , )
top
xk
(1 )
N k ( , )
yk
2
zk
(4.10)
bot
y N k yk
z
z
k
N k
mid
V3 k
2
(4.11)
61
Figure(4.1) Coordinate system: (a) nodal and curvilinear system (b) local
axes system
62
4.5 Displacement Field
The displacement through the element is defined by the three Cartesian
components of the midsurface node displacement k and two rotations of
the nodal vector
v3k
v2k
and
v1k
(of unite
k
k
1k
k
2
k 1
w
wk k 1
v2k
(4.12)
When the shell element is combined with a beam element, the local
( k , ) rotations must be expressed in terms of the global rotations (
k , k , k )
k
k
Tk k
k
k
where
(4.13)
Tk
geometric consideration as
v
z
v3k
v vy
3k
Tk
T
2k
T
1k
v3zk
v3yk
0
v3xk
v3xk
0
(4.14)
63
The five significant strain components are:
x
y
xy
xz
yz
u
x
v
y
u v
y x
u w
z x
v w
z y
(4.15)
in the local
coordinates system. The local derivatives are obtained from the global
derivatives of displacement u , v , and w by the following operation,
u
x
u
y
u
v
x
v
y
v
z
w
T
y
u
x
u
y
u
z
v
x
v
y
v
z
w
x
w
y
w
(4.16)
where is a transformation matrix
The derivatives of the displacements with respect to the global
coordinates are given by
u
x
u
y
u
v
x
v
y
v
z
w
1
J
w
w
w
(4.17)
64
x
x
J
x
z
z
z
(4.18)
xz
yz
(4.19)
where 0 may represent any initial strain due, for instance, to thermal
expansion and,
1
1
0
0
0
0
0
0
1
0
0
E 0 0
2
D
1
1 2 0 0
0
0
2r
1
0
0
0 0
2r
(4.20)
65
obtained when reduced integration is used [39,28] but the problem
appears again for very thin situation [29].
(ii) 9-Node Lagrangian
The performance of the 9-node element appears to be superior to the
Serendipity element when reduced or selective integration is used [29].
However, a serious problem may arise with this element when reduced
integration is employed. The stiffness matrix exhibits rank deficiency
which originates the appearance of spurious mechanisms, i.e., zero energy
modes[6,27]. Theses communicable mechanisms may result in very
erratic solutions, especially when the restraints imposed by the boundary
conditions are minimal. This drawback is not exhibited by the 8-node
Serendipity element with reduced integration.
(iii)Heterosis Element
This 9-node element was initially developed for plates [15]. The
Serendipity shape function is employed for the vertical deflection w and
Lagrangian shape functions for the two rotations. This heterosis concept
has been extended to the shell analysis using the Serendipity shape
function for the three translation degrees of freedom [18]. This element
exhibits improved characteristics when compared with the previous
quadratic elements. Its stiffness matrix possesses the correct rank when
selectively integrated and the element does not lock even in very thin
situation. Adopting uniform reduced integration, the element behaves
better than the 9-node Lagrangian element, with the appearance of
spurious mechanisms being improbable for practical application.
4.9 Numerical Integration of Shell Element
Although the shape functions of the finite element are simple
polynomials, the expressions for the stiffness become rather complex for
curved element due to the mapping between the natural coordinates
( , , )
( x, y , z ) .
66
element, a numerical integration is therefore necessary. Numerical
integration is widely used for isoparametric shell element.
Two Gauss points through the thickness for linear analysis are
sufficient to capture the bending behavior, because of a linear assumption
of strain through the thickness is assumed. High order Gaussian
quadrature rules (for instance 5-7 Gauss points) have been advocated for
nonlinear analysis, but a simple and general procedure to discretize and
integrate through the thickness is offered by the layered model.
In the shell plane (surface
67
Explanation based on mathematical consideration is also presented.
This explanation showed that identical terms are obtained for the
analytical and reduced approaches, but extra terms appear when the full
integration is used [17,23,27].
It was observed, however, that the reduced integrated element has
some deficiencies. The 8-node element of the Serendipity family to which
the reduced integration technique is applied may produce the locking
phenomena in problems below a certain thickness ration, even though its
performance is improved [29]. The 9-node Lagrangian and Heterosis
elements do not exhibit this drawback. Also for the 8-node element the
results may decline if the elements are distorted [29]. The 9-node
Lagrangian element displays undesirable zero or low energy modes when
the reduced integration technique is used [6,27]. This problem appears
also for the Heterosis element as shown below in the table which shows
the number of such zero energy modes for a quadratic isoparametric
element.
Although Lagrange shape functions may yield better membrane
behavior than the Serendipity pattern in highly distorted elements, the
Heterosis is preferred for thin shell application in view of its lower
number of mechanism.
Table(4.1) Number of zero energy modes
8-node S Heterosis
9-node L
No. of membrane mechanism 1(*)
1(*)
2+1(*)
No. of bending mechanism
1(*)
2+1(*)
3+1(*)
((*) not communicable in mesh of two or more elements)
xy
Excessive shear
0.6
1/ 3
Correctunder
shearlinear
Figure (4.2) Linear strain distribution in quadratic element
moment condition
68
4.9.2 Selective Integration
The concept of the selective integration was first employed by Doherty,
Wislon and Taylor [10], to obtain improved bending behavior in simple 4node elements. One-point Gauss quadrature was used on shear strain
terms, and 2x2 Gauss points quadrature was used to integrate the
remaining terms. Selective integration was first employed for the
degenerated shell element by Pawsey and Clough [28].
In the selective integration technique, bending (or bending and
membrane) energies are integrated using the normal rule, while the shear
and membrane (or just shear) terms are computed using the reduced
integration. A considerable improvement in the results by this technique
is noticed for plate and shell problems [16,17,27].
The selective integration, however, does not remove the entire
aforementioned difficulties, i.e. the selectively integrated element may
still lead to zero energy spurious modes in Lagrangian family elements.
Several modifications to the reduced or selectively integrated elements
were suggested to avoid spurious modes. The Heterosis element neither
locks nor posses unwanted zero energy modes in plate bending [15].
Using the selective integration technique, the number of mechanisms
becomes as follows:
Table(4.2) Number of zero energy modes
Heterosis
No. of membrane mechanism
1(*)
No. of bending mechanism
----((*) not communicable in mesh of two or more elements)
9-node L
2+1(*)
1
69
The assumed strain method was implemented to the 9-node and the 8node elements.
4.10.1 Assumed Shear Strain
The elimination of the shear locking is achieved by interpolating the
shear strain in the natural coordinate system. It was found that the
assumed shear strain
and
(4.21)
c1 c 2 c3 c 4 c5 2 c6 2
(4.22)
1
3
and
on which three points must be taken, Figure (4.3). One of the three
is on
0. .
component between two adjacent elements, the other two points should
be b 1 and
b 1
at lines
and
a .
Now
points
b 1, 0 and b 1 at
are three
. Therefore
Pi ( )Q j ( ) ij
(4.23)
i 1 j 1
3
ij
Pi ( )Q j ( )
(4.24)
i 1 j 1
where
is linear in
P1 ( z )
z z
2b b
, P2 ( z ) 1
, P3 ( z )
z z
2b b
(4.25)
70
Q1 ( z )
1
z
1
2
a
Q2 ( z )
1
z
1
2
a
(4.26)
For the 8-node element the assumed shear strains should be taken
as polynomials of at least the following degree
b1 b2 b3 b4 b5 2
(4.27)
c1 c 2 c3 c 4 c5 2
(4.28)
Ri ( , ) i
(4.29)
i 1
n
S i ( , ) i
(4.30)
i 1
R1
1
1 1 R5
4
a
4
(4.31)
R2
1
1 1 R5
4
a
4
(4.32)
R3
1
1 1 R5
4
a
4
(4.33)
R4
1
1 1 R5
4
a
4
(4.34)
R5 1
a
2
1
(4.35)
and
S i ( , ) Ri ( , )
(4.36)
71
coincides with the mid-surface of the shell element) can help to eliminate
the membrane locking behavior. To ensure that a special coordinate
system is uniquely defined, symmetry between the natural coordinate
system and the orthogonal coordinate system existed [19].
In the orthogonal curvilinear coordinate system, the in-plane
displacements can be expressed as follows:
u r u r 0 t r
(4.37)
u s u so t s
(4.38)
72
where
u ro
and
u so
u ro
and
are quadratic in
u so
and as well
(4.39)
em ss u s 0, s
(4.40)
em rs
1
(u r 0 , s u s 0 ,r )
2
(4.41)
Therefore, for
e m rr
3
e m ss
and
em rs
em rr Pi ( )Q j ( )e ij m rr
(4.42)
3
2
1 (a)
e m rs P ( ) i Q( ) j (e ij m rr / 2)
2
i 1 j 1
(4.43)
i 1 j 1
em ss P ( ) i Q j ( )e ij m ss
(4.44)
3
2
1 (b)
e m rs Pi ( )Q j ( )(e ij m rs / 2)
2
i 1 j 1
(4.45)
i 1 j 1
where
Pi (z )
4.26). Thus
and
Qi (z )
e m rr
and
e m ss
1 (b )
em rs
2
-direction and
similar to that used for the assumed shear strain, Figure (4.5). Finally, the
in-plane shear strain is quadratic in both directions, that is
e m rs
1 (a)
1
e m rs em( b)rs
2
2
(4.46)
73
for the assumed shear strain mentioned in section (4.10.1). Then, for
e m rr
and
1 (a)
em rs
2
em rr Ri ( , )emi rr
i 1
(4.47)
5
1 (a)
em rs Ri ( , )(emi rs / 2)
2
i 1
(4.48)
and for
e m ss
and
1 (b )
em rs
2
em ss S i ( , )emi ss
i 1
5
1 (b)
em rs S i ( , )(emi rs / 2)
2
i 1
(4.50)
(4.49)
74
4.11 Comparison Between the Assumed Strain Element and
Quadratic Shell Elements
Huang [19] performed patch test, locking test, and convergence test and
compared the results with the Heterosis and 9-node Lagrangian elements
with selective integration.
4.11.1 Patch Test
The patch test was performed with meshes of arbitrary quadrilateral
elements.
For bending patch test the assumed strain element and the 9-node
Lagrangian element passed the test but the Heterosis element did not. For
twisting patch test the assumed strain element passed the test while both
the 9-node Lagrangian element and the Heterosis element did not. For
shear patch test all elements passed this test.
4.11.2 Convergence Test
This test is a comparison between the analytical solution and the finite
element solution. Convergence studies are carried out on a uniformly
loaded, simply supported, square plate and also on a clamped circular
plate.
Comparison with the exact Reissner solution on a clamped circular
plate showed that the assumed strain element and the Heterosis element
provided a good idealization whereas the 9-node Lagrangian element
exhibited significant oscillation in the transverse displacement caused by
spurious zero energy mode [15].
Comparison was also carried out for a thin circular plate and the
results showed that for the assumed strain element a much better rate of
convergence was obtained.
For a thin skew plate of Morleys 30 0, the result of the finite
element solution was compared with the exact thin solution. The result of
75
the assumed strain element and the 9-node Lagrangian element gave good
agreement and better than the Heterosis element.
For Razzaques 600 -skew thin plate which was solved by the finite
difference solution, good convergence was obtained for all types of
elements.
More details are given in reference [19].
4.11.3 Locking Test
The element is tested to shear and membrane locking by imposing pure
bending and the results showed that no locking occurred.
4.12 Eccentric Curvilinear Beam Element
The best physical and geometric representation of the stiffened shell is by
modeling both the shell and the beam in their real behavior. The shell
element and the beam element were connected, for the analysis of
stiffened shell, displacement compatibility between and the beam element
and the shell element must be exist. This compatibility is achieved by
relating the displacement of the beam element to the axes at the
connected region between the shell and the beam (joint axes) to the
middle surface of the shell element (reference axes) [22].
4.12.1 Reference, Joint and Centroidal Axes
Figure (4.6) illustrates a typical 3-noded curvilinear beam element, with
six degrees of freedom comprising three global displacements ui, vi, wi
and three global rotations L, L,L. The points G1, G2, and G3 are the
nodes on the element centroidal axis and is a curvilinear coordinate
which varies between (-1,+1).
xG
xGL
3
y
N L ( ) y GL
G
L 1
z G
z GL
(4.51)
76
When introducing such an element to a thick shell element, two
additional axes should be introduced, the reference axis and the joint axis.
These are shown in Figure (4.6). These additional axes guarantee the
compatibility between the shell element and the beam elements.
The position of the points G L and E L with respect to the reference
axis can be specified by using vectors
e L , g L .
Consequently nodal
(4.52)
(4.53)
g ( ) N L ( ) g L ,
b( ) N L ( )bL
L 1
L 1
and
j b
b
(4.55)
where
a b g
and
c a b
(4.54)
c
k
c
77
The matrix of the orthogonal transformation from the local to the
global axes will therefore be given as:
( ) i( )
j ( )
k ( )
(4.56)
M Z
EA
GAY
GAZ
;D
GJ
EI Z
EI Y
u G , x
vG , x
wG , x
G , x
G , x
G , x
(4.57)
where, . , denotes local reference frame and , x for differentiation with
respect to x
A: cross sectional area for each layer
Ay Az
E: Modulus of elasticity
G: Shear modulus
J: Torsional constant
I y , I z :
The displacements and the rotations in the local and global systems
are related through [22]:
u 'G
'
vG
w'G
'
uG G
G
'
T vG ; G T G
G
wG G '
(4.58)
78
Neglecting the possible small difference between the local axes x
and the tangent to the centroidal axis it follows that:
d
1 d
dx t d
(4.59)
xGL
dN
L
t( )
y GL
L 1 d
z GL
3
(4.60)
BGL
0
N L .k T
N L j T
1 dN L T
t d
(4.61)
where
Ln
Ln
T
BGL
DBGL td
(4.62)
(4.63)
79
ARL
z L
y L
z L
0
x L
y L
x L
0
(4.64)
where
BRL BGL TL
(4.65)
I
TRL
0
(4.66)
0
ARL
80
81
82
4.14 Dynamic Equilibrium Equation
The dynamic equilibrium equations are obtained by using the principle of
virtual work [4,32] which states that for any arbitrary kinematically
consistent set of displacements, the virtual work done must be equal to
that done by external forces irrespective of the material behavior, as
dv
st
tds u
b u cu dv
(4.67)
and
b, u
are
cu
d , d
and
u , u and
by
d
can
u
the expressions
u N i ( , , )d i N d , u N d
(4.68)
i 1
m
u N i ( , , )d i N d
(4.69)
i 1
m
u N i ( , , )di N d
(4.70)
i 1
Bi d i Bd
(4.71)
i 1
d T Md Cd p ( d ) d T f
(4.72)
83
in which the mass matrix M , the damping matrix C , the internal force
vector
p (d )
f e have
the following
element contributions
Me
Ndv
(4.73)
ve
Ce
cN
dv
Ndv
(4.74)
ve
pe
(4.75)
ve
fe
tdS
se
where
Nbdv
(4.76)
ve
s e and v e
Md Cd p ( d ) f e
For linear elastic situation, the stresses are related to the strain
as follows:
D DBd
(4.78)
p (d ) can
be written as
p (d ) K s d
Ks
(4.79)
DBdv
ve
Me
The resulting element matrix is a full matrix and the global mass matrix
has the same structure as the stiffness matrix. A very convenient
simplification is achieved by lumping the element mass to nodes. This is
done in order to avoid the factorization process and to solve the
associated problem by explicit algorithm. Several lumping schemes are
proposed and in general depending on intuitive reasoning [14].
84
The sub-matrix of the element mass matrix linking node ( i ) and (
j
N j dv
(4.80)
ve
which does not lead to a diagonal matrix, when the adopted shape
functions are identical to those used in the evaluation of the element
stiffness matrix.
In the present work a consistent mass matrix is adopted for both the
degenerated shell element and the beam element.
Acceleration compatibility between the degenerated shell element
and the beam element is achieved by the following transformation:
T
M RL TRL
M GLTRL
where
TRL :
M GL is
(4.81)
(4.82)
d 0
is adopted so
that[14,32]:
cij d mij
where d
(4.83)
2 r r
in which r and r are the damping ratio and the circular frequency for
the r th mode.
85
4.17 Layered Model
In reinforced concrete shell problems, the nonlinear behavior of
compressive concrete, concrete cracking and reinforcement response need
a convenient representation across the thickness. Layers of a different
thickness can be employed, Figure (4.9). The natural coordinate
varies from
which
x
( y)
h/2
h / 2
x
( y)
dz
h n i
x i
2 i 1 ( y )
(4.84)
-bending moments,
M
x
( y)
( xy )
h/2
h / 2
zdz
x
( y)
xy
h2
4
i i
i
x
i 1
( y)
xy
(4.85)
-shear force
Q
x
( y)
h/2
h / 2
xz
( yz )
dz
h n i
xz i
2 i 1 ( yz )
(4.86)
86
1 1 1
Ke
DB J ( , , ) ddd
1 1 1
(4.87)
1 1
h j
1 1
where
DB J ( , , )
2h j
d d
j 1
J ( , ,
(4.88)
jth
87
element). This means less number of elements are required to represent a
particular structure.
In the present study simple implementation is formulated for
reinforcement with arbitrary location in the degenerated shell element.
Perfect bond is assumed between the reinforcement and the
concrete, which means that the strain of concrete and reinforcement has
the same values at the common points.
The basic idea is only to introduce the correct direction cosines
relating the uniaxial strain in steel in the local coordinate system of the
bar to the global coordinate system. Then the same strain-displacement
matrix of the shell element will be used with the corresponding direction
cosines of the embedded bar.
4.18.1 Geometric Definition of Bar Element
Figure (4.10) shows the representation of the reinforcement by a bar
element in the degenerated three dimensional Ahmed element. In the
present study 3-noded curved element is used to represent the geometry
of the reinforcement. These elements are used to calculate the direction
cosines at the Guass points and to calculate the length of the bar using a
numerical integration method in order to evaluate the stiffness matrix and
stresses.
It is only sufficient to give the natural coordinate of the bar element
within the parent element (shell element). Then the Cartesian coordinates
of the bar element are calculated using equation (4.10). By using these
coordinates and the shape functions below, one can calculate the Jacobian
matrix, which is used to compute the length of the bar at the Gauss points
as follows:
The global coordinate of the bar is related to the natural coordinate
using the three-noded elements:
88
x i 1 N i x i
n
(4.89)
y i 1 N i y i
n
(4.90)
z i 1 N i z i
n
where
(4.91)
Ni
N1
1
( 1)
2
(4.92)
N 2 (1 2 )
N3
where
(4.93)
1
( 1)
2
(4.94)
Jb
(4.95)
(4.96)
Embedded bar
2
Mid surface
89
u
T u
x
b
x
y
u
v
x
v
y
v
z
w
x
w
b
y
w
y
z
b
Jb
where
Jb
(4.98)
Or
b T
x
x
y
x
z
x
(4.99)
90
the B-matrix is calculated.
(4.100)
E
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
(4.94)
(4.101)
Db Bb dv
(4.102)
91
References
1. Abdul-Razzk, A. A. (1996). "Nonlinear finite element analysis of
fibrous reinforced concrete structural members." Ph.D. Thesis, Mosul
University, Iraq, 234 pps.
2. Ahmed, S., Irons, B. M., and Zeinkiewicz, O. C. (1970). " Analysis
of thick and thin shell structure by curved finite elements."
International Journal for Numerical Methods in Engineering, Vol. 2,
pp. 419-451.
3. Arnesen, A. (1979). "Analysis of reinforced concrete shells
considering material and geometric nonlinearities." Division of
Structural mechanics, Norwegian Institute of Technology, University
of Trondheim, Norway, Report No. 79-1.
4. Bathe, K. J. (1996). "Finite element procedure." Prentice-Hall Inc.
5. Bathe, K. J., Iosilevich, A., and Chapelle, D. (2000)."An evaluation
of the MITC shell elements." Computers and Structures, Vol.75, pp. 130.
6. Bicanic, N. (1979). "Spurious modes in two dimensional
isoparametric elements." International Journal for Numerical Methods
in Engineering, Vol. 14, pp. 1545-1557.
7. Bouabdalla, M. S., Batoz, J. L. (1996). "Formulation and evaluation
of a finite element model for the linear analysis of stiffened composite
cylindrical panels." Finite Element in Analysis and Design, Vol. 21,
pp. 265-189.
8. Barut, A., Madenci, E., Tessler, A., and Starnes, J. H. (2000). "A
new stiffened shell element for geometrically nonlinear analysis of
composite laminates." Computers and Structures, 77, pp. 11-40.
92
9. Chang T. Y., Taniguchi H., Chen W.F (1987)."Nonlinear finite
element analysis of reinforced concrete panel. " ASCE, Vol. 113, No.1.
10. Doherty, W. P., Wilson, E. L., and Taylor; R. L. (1969). "Stress
analysis of axi-symmetric solids using higher order quadrilateral finite
elements. " Structural Engineering Laboratory, Report. No. SESM 693, University of California, Berkeley.
11. Elwi, A. E., Hrudy, T. M. (1989)."Finite element model for curved
embedded reinforcement." ASCE, Vol. 115, No.4.
12. Fezans, G., and Verchery, G. (1982)."Some results on the behavior
of degenerated shell element." Nuclear Engineering Design, Vol. 70,
pp. 27-35.
13. Gallagher; R. H. (1976). "Problem and progress in thin shell finite
elements analysis, in finite elements for thin shells and curved
members." D. G. Ashwell and R.H Gallagher (eds.), John Wiley &
Sons, London.
14. Hinton, E. (1988). "Numerical methods and software for dynamic
analysis of plates and shells." Pineridge Press, Swansea UK.
15. Hughes, J. R., and Cohen, M. (1979). "The heterosis finite element
plate bending. " Computers and structures, Vol.9, pp. 445-450.
16. Hughes, J. R., Cohen, M., and Haroun, M. (1978). "Reduced and
selective integration technique in finite element analysis of plates."
Nuclear Engineering Design, Vol.46, pp. 203 222.
17. Hughes, T. J. R., Taylor, R. L., and Kanoknuchulchai, W. (1977). "A
simple and efficient finite element for plate bending. " International
Journal for Numerical Methods in Engineering, Vol. 11, pp. 15291543.
93
18. Hughes, J. R. and Liu, W. K. (1981). "Nonlinear finite element
analysis of shells. Part 1, three dimensional shells." Computer
Methods in Applied Mechanics and Engineering, Vol.26, pp. 331-362.
19. Huang, H. C. (1989). "Static and dynamic analyses of plates and
shells." Springer-Verlag, Berlin Heidelberg.
20. Idelsuhn, S. (1981). "On the use of deep, shallow or flat shell finite
elements for the analysis of thin shell structure." Computer Methods in
Applied Mechanics and Engineering, Vol.26, pp. 321-330.
21. Jiang, J., and Olson, M. D. (1991). "Nonlinear dynamic analysis of
blast loaded cylindrical shell structures." Computers and Structures,
Vol.41, No.1, pp. 41-52.
22. Jirousek, J. (1981). "A family of variable section curved beam and
thick
shell
or
membrane
stiffening
isoparametric
element."
94
26. Morris, I. R., and Dawe, D. J. (1980). "Free vibration of curved
plate assemblies with diaphragm ends." Journal of Sound and
Vibration, Vol.73, No.1, pp. 1-17.
27. Parich, H. (1979). "A critical survey of the 9-node degenerated shell
element with special emphasis on thin shell application and reduced
integration."
Computer
Method
in
Applied
Mechanics
and
95
34. Ranjbaran, A. (1991). "Embedding of reinforcement in reinforced
concrete elements implemented in DENA." Computers and Structures,
Vol. 40, No. 4, pp. 925-930.
35. Sinha, G. and Mukhopadhyay, M. (1997). "Static, free and forced
vibration analysis of arbitrary non-uniform shells with tapered
stiffeners. Computers and structures." Vol.62, No.5, pp. 919-933.
36. Sinha, G. and Mukhopadhyay, M. (1995)."Transient dynamic
response of arbitrary stiffened shells by the finite element methods
Journal of Vibration and Acoustics." Vol. 117, pp. 11-16.
37. Stolartski, H. and Blytschko, T. (1982). "Membrane locking and
reduced integration for curved elements." Journal of Applied
Mechanics, Vol.49, pp. 172-176.
38. Thannon, A. Y.