CC BY
23
Buckling analysis
THEORETICAL APPROACH FOR THE GEOMETRICALLY NONLINEAR
BUCKLING ANALYSIS OF SINUSOIDAL VELAROIDAL SHELLS
GIL-OULBE MATHIEU, PhD, Assoc. Prof., SHAMBINA S.L., PhD, Assoc. Prof. DAU TYEKOLO, PhD, Assis. Prof. RUDN University, Russia 117198, Moscow, Miklukho-Maklaya street, 6;
E-mail: [email protected], [email protected], tiek.d@ hotmail.com
The main subjects of the article are theoretical strength and buckling analyses of a sinusoidal velaroidal shell subjected to its self-weight and uniformly distributed load in geometric nonlinearity. A short history of finite element applications to shell buckling is given. The finite element method is used in its matrix formulation. The elastic stiffness matrixes in the local coordinate system of the membrane element are defined in their general form. The out-of-plane geometric stiffness matrix for the plate along the same lines is derived.
KEY WORDS: sinusoidal velaroidal surfaces, nonlinear buckling analysis, linear elastic stiffness matrices, normal stiffness, geometric stiffness matrix, chain rule differentiation.
Velaroidal surface is a surface of translation on the flat rectangular plan with a generating curve of variable curvature [1],[2]. Thus, the surface is limited by four mutually orthogonal contour straight lines (kx = ky = 0) lying in the same plane.
A sinusoidal velaroid generates by two families of half waves of the sinusoids lying in mutually perpendicular planes and facing by convexities into the same side [3]. Each set of sinusoids has the identical period. Sinusoidal velaroid is limited by a flat rectangular contour.
The history of finite element applications to shell buckling is extensive going back to the work of Clough and Johnson (1968). The natural mode contribution of Argyris et al. (1977) was a major addition to shell theory. It was recently modified to include elastoplastic effects (Argyris et al. 2000). Horrigmoe and Bergan (1978) used the co-rotational method for nonlinear analysis and Bathe and Ho (1980), Hsiao (1987), Mohan and Kapania (1997), Peng and Crisfield (1992) improved element performance along those lines. The 3-D elasticity "degenerate" element of Ahmad et al. (1970) was followed by, among others, Bathe and Balourchi (1980), Hughes and Lui (1981), Dvorkin and Bathe (1984), and Buechter and Ramm (1992). In an excellent review, Ibrahimbegovic (1997) addresses the various approaches and the complex issues involved.
Here the derivation of the geometric stiffness matrix is somewhat different but consistent with the approach used throughout this text. The linear equilibrium equations for a flat triangular shell element in its local coordinates system are first perturbed to yield the in-plane geometric stiffness ma- Fig.l. A sinusoidal velaroidal shell trix. Then out-of-plane considerations that
involve the effect of rigid body rotations on member forces yield an out-of-plane geometric stiffness matrix. The shell element that was chosen for that purpose combines the constant stress triangle (CST) flat triangular membrane element (Zienkeiwicz (1977)) and of the discrete Kirchhoff theory (DKT) flat triangular plate element (Batoz et al. (1980)).
Let's consider a sinusoidal velaroidal shell [4],[5] with the inner radius r0 = 0, the outer variables radii from 10m to 20m and the number of waves n= 8 (Fig. 1).
Structural Mechanics of Engineering Constructions and Buildings, 2017, № 3 The Geometric Stiffness Matrix of Triangular Element Shells.
The local geometric stiffness matrix of the shell element is split into three distinct matrices:
[jre-\shell _ [jre-\mem I [jre~\plate . vjre'ishell
L^GiTOTAL — LagJ/P "T LaGJ/P "T" LaGJ0P , (1)
where the first, second and third terms on the R.H.S. of Eq.1. represent the in-plane geometric stiffness matrix of the membrane, the in-plane geometric stiffness matrix of the plate and the out-of-plane geometric stiffness matrix of the shell element respectively. The total, 'tangential' stiffness matrix for use in nonlinear analysis will include, in addition, the linear elastic stiffness matrices of a plane stress triangular element (membrane) and that of a triangular plate element.
The geometrically nonlinear triangular shell element has eighteen local degrees of freedom (DOF's): 3 displacements and 3 rotations at each node. The membrane element contributes to nine displacement DOF's only. The basic three noded constant stress triangular flat element has only six local displacement DOF's that are shown in Figure 2. The out-of-plane contribution (the normal stiffness) of the membrane element to the basic local shell element is a displacement DOF in the direction normal to the plane of the element.
Fig. 2. Triangular membrane element Fig.3. Triangular plate element
The flexural element contributes to eighteen local DOF's. The basic three noded plate triangular flat element has only nine DOF's that are shown in Figure 3. The inplane contribution to the basic local element adds two displacement DOF's in the plane of the element. The out-of-plane contribution adds a rotation DOF in the direction normal to the plane of the element.
In-Plane Contribution of the Triangular Membrane Element.
The elastic stiffness matrix in the local coordinate system of the membrane element has the general form:
[K§]mem = JA BmemDBmem t
dxdy. (2)
The local in-plane contribution of the membrane element to the geometric stiffness matrix is the gradient of the nodal force vector:
\ise-\mem _ r7Fe
[AG]/P _ v "mem. (3)
In-Plane Contribution of the Triangular Plate Bending Element.
The general form of the elastic stiffness matrix of the DKT triangular plate element may be written as:
[Ke]Plate = 2A jl J^ ^ Bplated^df1, (4)
where ^ and r/ are the usual area coordinates and Bplate, which appears explicitly in Batoz et al. (1980) as a function of r/ is slightly more complex than the Bmem of the membrane. Here again the local in-plane contribution of the plate element to the geometric stiffness matrix will be derived as the gradient for fixed M, of the element nodal force vector which is given as:
i i-v
Opiate - I 2A J J 12BlplateDBviated^dqjqe -
0 0
1 l-V
where
1 l-V
-2 aJJ B^;iatel1^DBplateqe)d^di1-2Aj J B1plateMd^dr]
I 0 0 0
M - {Mxx,Myy,Mxy}T,
Fe — Fe Fe\T■ rplate lri ,rj ,rm) ;
Fy {-^Vz, M, Mj-y^ v i, j,
qe — {wx 6xl 6yl w2 6x2 6y2 w3 6x3 0y3). The in-plane geometric stiffness matrix may be written symbolically as:
(5)
(6)
(7)
(8)
' (An)T' «7* (A \plate ' (/iimJip
plate - UFe - _ v "plate ~ (Aim)Plate V Jmsjp
(A \plate L (nmiJlp (a \plate (/imiJip (A \plate mm)ip -
d Ffz ^ ^Vz 0"
dxs dys
where (A \pla.te -(/irsJip ~ dMrx dMrx dxs dys 0
dMry dMry 0
L dxs dys
(9)
(10)
The expressions for the individual terms of the geometric stiffness matrix were obtained in closed form using symbolic algebra.
When applied to the plate, chain rule differentiation with respect to the coordinates yields
dF;iate- £ [¿f24// (Zplate)fixedMdi;d^dxr
r=i,],m \ 0 0
1 1 -1
+ Wr\2A SS (BPlate)fixedMd^dr1jdyr
0 0
1 1-V
+
X [3^(2^/ S BPl*teMfixedd^dll)dXr
■=i,j,m r\ 0 0 /
1 1 -v
+ ^(2Af J BplateMfixed dZdv)dyr
0 0
(11)
where dxr =ur; dyr =vr; r=i,j,m.
The first expression on the R.H.S. of Eq. l1 returns the elastic stiffness matrix whereas the second expression becomes the geometric stiffness matrix. It is left to clarify what is meant by Mfixed in Eq. 11. The moment vector M, of Eq. 6 contains three components that are functions of ^ and ^ and defined as:
mxx(S,v) - (M;
1
XX
M^H + (M& - Mlxx)V + ML,
37
(12)
Myyv) = (M^y - MiyyX + (M™ - Miyy)^ + Miyy,
Mxy(t,ri) = (Mïy - Mlxy)t + {M™ - Mixy)r] + Mixy, where
[Mlxx Miyy Mlxy}T = M(t = 0,^ = 0),
{m4 Mjy}T = M(t = 1,^ = 0),
(M- M™ Mfyf = M(^ = 0,^ = 1).
(13)
(14)
(15)
(16) (17)
When the components of M (Eqs. 12-14) are inserted into Eq. 11 it is the values of these components at the nodes (Eqs. 15-17) that are held fixed.
Out-of-Plane Contribution to the Shell Geometric Stiffness Matrix.
This section will derive the out-of-plane geometric stiffness matrix for the plate along the same lines and subsequently present a combined out-of-plane contribution. The derivation starts with the change in a vector G, due to a small rotation that is given by Goldstein (1950) as
dG = a)XG, (18)
where w is the rigid body rotation vector due to changes in the geometry. In terms of joint displacements with respect to the local coordinate system, two components 6 and -cp. of the rigid body rotation are obtained from Figure 4. The third component is Fig. 4. A triangular finite element chosen for the plate arbitrarily, as the local z-rotation in its coordinate system
of node i. Recall that for the membrane this component is included in the in-plane contribution:
; <*y' +-t(5i)z> ^z' = 9iz, (19)
where (Sr)z' is the displacement in the local (primed) z-direction of node r. At each node, the forces, moments and rotations may be written as
Fr = Frzk'; Mr = Mrxi' + Mryj' + Mrzk'; w = Mxd' + uyij' + Mz>k'
for r = i, j, m and the changes in the force and moment vectors are given in more detail as
(20) (21)
(22) (23)
0 F 1 rz 0 (Ux>
dFr = -Fr x = -F 1 rz 0 0 K'
0 0 0 [wz>
0 Mrz Mry
= —Mr x w = -Mrz 0 Mrx
Mry -Mrx 0
¡y'
The out-of-plane stiffness contribution of the plate is now obtained as:
(dFe)plate = [KG]
plate OP
i£ = -Fexw =
where
Fe =
(Pi 1
Mi
Fj M, F 1 m vmJ s Fr =
0 Mrz -Mry
-Mrz 0 Mrx
. Mry -Mrx 0
Now w can be described in terms of the displacement vector as:
di
u = ASe = [Ai Aj Am] ^ dj
, dr
(24)
where
(At) =
0 0
e — c
ae 1
0 0 —-
e
000 000
L0 0 00 0 1-
0 0 c ae 0 0 0
M = 0 0 1 Q 0 0 0
0 0 0 0 0 0
0 0 1/a 0 0 0
(Am) = 0 0 0 0 0 0
0 0 0 0 0 0
(25)
(26)
(27)
and [di] =
(Si)x
(Si)y
(Si)
(0i)x (0i)y ■(0i)z
=
[dm] =
(Sm)y (Sm)z
(^m)x (^m)y
.(O
(28)
Carrying out Eq. 22 for the plate and adding the out-of-plane contribution to the geometric stiffness matrix, in the local coordinate system more explicitly (adjusted to 18 DOF's) for the membrane, results in the following out-of-plane shell geometric stiffness matrix with respect to the local coordinate system:
[jreTshell [KG]OP
(Ai)! (Ai)2 (Aj\ (Aj)2
(Am)i (Am)2
(Aih (Am)3 J
(29)
where
r)l —
0 0 ar 0 0 0
0 0 br 0 0 0
0 0 ar 0 0 0
0 0 cr 0 0 -Mry
0 0 dr 0 0 Mrx
_0 0 er 0 0 0
(30); (Ar)2 =
0 0 fr 0 0 0"
0 0 gr 0 0 0
0 0 Pr 0 0 0
0 0 hr 0 0 0
0 0 ir 0 0 0
_0 0 Jr 0 0 0-
(31)
(Ar)3 =
0 0 0 0 0 0"
0 0 ky 0 0 0
0 0 Ar 0 0 0
0 0 0 0 0 0
0 0 lr 0 0 0
_0 0 mr 0 0 0-
(32)
ил e~c (JvV . n __(p) c (JvV . _(£V)y. _ _ Гг V ae e ' Pr (Гг V ae e ' Лг a ' "r e
e _ с Mrz e _ с Mry e _ с Mrx
ae e ae a ae e
/t v4 С » ^^^ f > * С
r ~ < fl'r *rz ~; "г ~ < ^r '^rz ~;
с Frz Mrz
jr MyZ ; ^^ ; If
ae a a
Conclusions: The theoretical approach of the buckling and strength analyses of the sinusoidal velaroidal shell is worked out using the finite element method. The inplane and the out-of-plane shell geometric stiffness matrices with respect to the local coordinate system are obtained. This result gives a possibility of further numerical strength and buckling analyses of the sinusoidal velaroidal shells.
R e f e r e n c e s
1. Mihailescu M., Horvath I. (1977). Velaroidal shells for covering universal industrial halls, Acta Techn. Acad. Sci. Hung., 85 (1-2), pp. 135-145.
2. Krivoshapko S.N., Mamieva I.A. (2012) Outstanding space structures of the last 20 years, Montazhnye i spetsialnye raboty v stroitelstve, № 12, pp.8- 14.
3. Krivoshapko S.N., Gil-oulbe Mathieu. (2013). Geometry and strength of a shell of velaroidal type on annulus plan with two families of sinusoids, International Journal of Soft Computing and Engineering (IJSCE) ISSN: 2231-2307, Volume-3, Issue-3, July 2013.
4. Krivoshapko S.N., Shambina S.L. (2009). Investigation of surfaces of velaroidal type with two families of sinusoids on annular plan, Structural Mechanics of Engineering Constructions and Buildings, № 4, pp. 9-12.
5. Gogoberidze YaA.(1950). Pokrytie "Darbazi", Tbilisi: «Tekhnika da Shroma», 278 p.
Л и т е р а т у р а
1. Mihailescu M., Horvath I. (1977). Velaroidal shells for covering universal industrial halls// Acta Techn. Acad. Sci. Hung., 85 (1-2), pp. 135-145.
2. Кривошапко С.Н., Мамиева И.А. Выдающиеся пространственные сооружения последних 20 лет // Монтажные и специальные работы в строительстве. - 2012, № 12. - С. 8-14.
3. Krivoshapko S.N., Gil-oulbe Mathieu. (2013). Geometry and strength of a shell of velaroidal type on annulus plan with two families of sinusoids// International Journal of Soft Computing and Engineering (IJSCE) ISSN: 2231-2307, Volume-3, Issue-3, July 2013.
4. Кривошапко С.Н., Шамбина С.Л. Исследование поверхностей велароидального типа с двумя семействами синусоид на кольцевом плане// Строительная механика инженерных конструкций и сооружений. - 2009, № 4. - С. 9-12.
5. ГогоберидзеЯ.Ф. Перекрытия «Дарбази». - Тбилиси: Техника да Шрома, 1950. - 278 с.
К РАСЧЕТУ НА УСТОЙЧИВОСТЬ В ГЕОМЕТРИЧЕСКИ НЕЛИНЕЙНОЙ ПОСТАНОВКЕ ОБОЛОЧЕК В ВИДЕ СИНУСОИДАЛЬНОГО ВЕЛАРОИДА
ЖИЛЬ-УЛБЕ МАТЬЕ, С.Л. ШАМБИНА, ДАУ ТЬЕКОЛО
Российский университет дружбы народов, Москва
В статье рассматривается расчет на прочность и устойчивость в геометрической нелинейной постановке для синусоидальных велароидальных оболочек под действием собственного веса и равномерно распределенной нагрузки. Приводится краткая история применения метода конечного элемента для изучения потери устойчивости оболочки. Метод конечного элемента используется в матричной формулировке. Эластичная матрица жесткости в локальной системе координат мембранного элемента определяется в общем виде. Получена матрица геометрической жесткости вне плоскости для пластины вдоль тех же линий.
КЛЮЧЕВЫЕ СЛОВА: синусоидальные велароидальные оболочки, расчет на устойчивость в нелинейной постановке, линейная матрица упругой жесткости, нормальная жесткость, геометрическая матрица жесткости.
