ScanGate document
Trang 1VNU JOURNAL OF SCIENCE, Mathematics - Physics T.XIX, No3 - 2003
ON THE ELASTOPLASTIC STABILITY PROBLEM OF
THE CYLINDRICAL PANELS SUBJECTED TO THE COMPLEX LOADING WITH THE SIMPLY SUPPORTED
AND CLAMPED BOUNDARY CONSTRAINTS Dao Van Dung
Department of Mathematics, College of Science, VNU
Abstract In this paper, an elastoplastic stability problem of the cylindrical panels under the action of the compression force along the generatrix and external pressure has been investigated By the Bubnov-Galerkin method, we have established the expression for determining the critical loads The sufficient condition of extremum for a long cylindrical panels was considered Some numerical results have been also given and discussed 1 Formulation of the stability problem and fundamental equations
Let us consider a round cylindrical panel of thickness h and radius of the middle surface equal to R We choose a orthogonal coordinate system Ozxyz so that the plane Oxy coincides with the middle surface and the axis Ox lies along the generatrix of cylindrical panel while y = RO, with 6,-the angle circular arc and z in the direction of the normal to the middle surface Denote the sides of cylindrical panel by a and b respectively to the axis Ox and Oy
Suppose that the cylindrical panel is simultaneously subjected to the comp force of intensity p(t) along the generatrix and external pressure of intensity q (/) increasing monotonously and depending arbitrarily on any loading parameter ¢ We have to find the
critical values t = t,, Pe = p(ts), Ge = q(t.) at which an instability of the structures
appears In order to investigate the proposed we will use the criterion of bifurcation of equilibrium states and dont take into account the unloading in the cylindrical panel Afterhere we will present the fundamental equations of stability problem
lon
1.1 Pre-buckling process
Trang 2The components of the strain velocity tensor determined according to the theory of elasto-
plastic proc are of the form [1] "Nà ~b+ s4) - Us - 5) tay = (củ + g8) = (9 (9= gỡ), (1.2) = —(211 +222); 612 = €13 = 23 = 0, where 23% be des 1 1 PP + 44 ~ 5P ~ 304 ) P-pg+ge ` ‘The are-length of the strain trajectory is respectively calculated by the formula ds 2 dt mì + nể» +é},)? = F(s,t) (1.3)
ciated with the relations (1.2), (1.3) and boundary con- and strain state at any point M in the structure at any So the equilibrium equations ass
ditions entirely define the stre moment of pre-buckling process 1.2 Post-buckling process
The system of stability equations of the thin cylindrical panel established in [5] is written in the form
» ow es tow “ise aw - 9 7 Pow 7 Pow 1 Las UAT TSAR TS Hy + TN (° 2m 74 Øy? Hôz2
Hy ap ae, N ow
3 3, 5 =0, 5
Agee + Oa age + hôn †T gui TỦ (1.5)
where the coefficients a,, 3; (i = 1,3, 5) are calculated as follows ma1-G(-§)s on =2-5 (1-5) 8 ¬`= ø =1+1(7~1) OP, 6) f=32+(5~I)ÈP 04T), i= 147 (3 - th, For solving the stability problem of cylindrical panel, we consider two types of kinematic constraints following * The cylindrical panel is simply supported at the four edges x = 0, r = a, y= 0, y=b
Trang 32 The solving method for the simply supported cylindrical panel at four edges We find the increment of deflection 6w in the form
ne, 2wru
ôn = ` > mn Amn sin —— sin (2:1) 2.1 it is easy to see that this solution satisfies the kinematic boundary conditions
Substituting the expression of dw into (1.5) we receive the particular solution y as follows Moe max , 2n g= > › Brn sin 2 - sin > (2.2) m=l n=l where
Bon = (ME) helo (SE)+04(%) CE) HACE) 0
It is seen that the system of functions
(m,n =1,2, ,M)
na, 2n 6Umn = sin —— sin C50 a b
is linearly independent Therefore we can apply the Bubnov-Galerkin method for estab- lishing an expression of critical forces
First of all, substitute the expressions of dw and y from (2.1), (2.2) into (1.4),
afthward multiply both sides of the just received equation by éwiy = sin = sin Zin and integrate that equation following x and y Finally we get
ab 5 5 3
25980 L How + How + 9 OP 5w + Piw 1p Dea +84 Ox? dy? os Oy? ˆ hÀN Paạz T4 Ø2 — H0x2 0
37 sin
wn sin Yardy =0 (i,j = 1,2, -,M) (24)
For taking this integral, it needs to use the result
ab i i: j
Nam _ an ine nmy _ 2yjmy { 0 with m#i; nA Jj
sin sin — sin sin = 4 ab
b b with m =i, n=j
After series of calculations, the relation (2.4) gives us
Trang 4Because of the condition on the existence of non-trivial solution i.e Amn # 0 then we receive the expression for determining critical loads 2Á) +(22) (CE) 0) tm() TE) +) CE) anb 3b Putting X = (=) Y¥ =n?, i= —; the relation (2.6) can be rewritten in the other J (2.6) 2na h form 4 AN? (aX tay +22) (8.x ++ 2 Tu cha li 3) 22c x)( 5 a) 2i (2.7) ức 5 Y(pX +4)(8X + + Ở) — Tang Minh ininizing 7?, it means ax tn ĐỂ — 0, ĐỂ — 0, that yield 1 By =O that yields ĐỀN = ———.—n`^ (28) >»e(y+Ÿ)fex ta +Ÿ) (5=)6x+»+#)=(5-)fnz+a+# 2q ‘i F Ps\ _ , Fe ca] tử te+Ÿ)(8X +4 +) ~u (2.9) x Substituting the expressions (2.8) and (2.9) into (2.7) we obtain , AN? Be
2= _ re(p+ 4)" {aux toast Shox ++ 2} (2.10)
where X is found from the equation (2.9)
Applying the loading parameter method [1], and solving simultaneously the equation
(1.3) and (2.10) we can find the critical values t., ps = p(ts) de = @(ts)
Trang 5Substituting the values of as, G5 and X = X, into (2.11) we obtain 4N202 = (a 7 TH = 12 { N Gp-a)? pH zg BN } (2.12) Pie; (-Ủp Nrai † pm
3 The solving method for the simple supported cylindrical panel at y = 0, y = b simultaneously clamped at the sides x = 0, x =a
The kinematic boundary constraints of stability problem are satisfied completely by choosing M M bw= D> D> Cmn (1 = 008 SE) dn “rẽ (3.1) m=1n=1 Using the equation (1.5) and the expression of dw we can find the particular solution # in the form M M 2mmz 2m ved g= > > Dyn COs sin te (3.2) m=1n=l nạ) In(ŒP)'+a(9JŒP)+»(PJT an ‘i Qmrx\ 2n7t
It is possible to prove that the system of functions dWnn = ( — cos 27") sin SAY is linearly independent Then we can use the Bubnov-Galerkin By the same method presented in the above part we change the equation (1.4) into a relation as follows b / a 280 ow lượn C2 9 ow fi củng Gat + Margy =.nG: (p Oa? ` Tay 9Ø Zine 2Jm “PRzazjú~e — b drdụ=0 (¡,j = 1,2, , Äf) (3.4)
For taking this integral above all we substitute éw and ¢ represented by (3.1) and (3.2)
into (3.4), afterwards integrate that received expression We will obtain a system of linear
algebraic equations with the unknowns C,; which is written in the matrix form la;][Œu¿] =
i,j =1,2, ,M (3.5)
Because of the condition on the existence of non-trivial solution i.e Cj; # 0 then the determinant of the coefficients of C;; must be equal to zero
Trang 6Associating this expression with (1.3) we can determine the critical values t., p = p(t.),
qe = Gta)
Note that a development of the determinant (3.6) in general case is complicated mathematically therefore we will take the solution in the first approximation
In this case we choose 6w and ¿ in the form 2mza+ 3um ðtt = Cnn ( = cos wat din a b 8 2m7\2_, = ¬ ome, i, SRY, (th) ch ụ?
Substituting dw and ¿ into (3.4), integrating that relation and taking into account the condition Cian # 0, leads us fos 2) om 2) CE) to 2S) — pe EY eZ) + Bs 9 /2mm ủ 2mz 3m2 —— /2nz 2nzy\4)-! — =0 gì + pe (Oa) AE) ea) At) +a(AF) “pee 2) l b
Trang 7For finding the critical value ¢, of loading parameter t, we need to solve simultane-
ously the equation (1.3) and (3.12)
After determining t we can obtain the critical forces as follows Px = p(te), Ge = a(ta) Now consider the case of a long cylindrical panel Based on {2] leads us 12Nz°?a;0 bẦN = TC cóc C THỦ ý a ccc 3.15 €=1 7<l, 7 (pn + 3q)55 — Cuế Co Tere (3.13) ee a oe es a aay OP Ổ: The minimization of the expression i* in (3.13), i.e On = 0, yields » = ; =1 Moreover 9.2 N 6 oy - Nas >0 (3.14) On? In=n (30 p+ 2B)’ So the sufficient condition of extremum is satisfied Taking into account a5, (5, 7., the relation (3.13) becomes _ ee 1/N) =a?) , SEN Aba l)p“m+õ|† san" 2 _— 122W? = Remarks
1) If the cylindrical panel has a very small curvature i.e KR — +00; g = 0 and
m= 1,n=1 then the expression (2.7) coincides with the result of {1, 5, 7]
2) If b = 27 that means the cylindrical panel becomes a closed round cylindrical shell, then the expressions (2.10), (2.12), (3.12), (3.15) return respectively to the previous well-known results
4 Numerical calevlations and discussion
A numerical analysis is considered on the long cylindrical panel made of the steel 30XICA with an elastic modulus 3G = 2.6- 105 MPa, an yield point ¢, = 400 MPa and
the material function ¢(s) presented in [1]
The relations for determining the critical loads are given in the form:
* Formulae (2.12) and (1.3) for the part a) of the examples * Formulae (3.15) and (1.3) for the part b) of the examples
The numerical results are realized by the program of MATLAB Example 1 The complex loading law is given in the form
Trang 9Example 2 Suppose that the complex loading law is of the form
p= p(t)=po+pil®; po=2MPa, pi =0.1MPa
q=4()=qo+4@if; qo= 1MPa, qị =0.1MPa
Trang 10‘The above received results leads us to some conclusions
1 By using the Bnbnov-Galerkin we have solved the elastoplastic stability problem of the cylindrical panels with two types of kinematic boundary constraints
2 We have shown, for long cylindrical panel, the sufficient conditions of extremum 3 The critical loads of the simply supported cylindrical panels are always stnaller than critical loads of the clamped cylindrical panels (see tables 1, 2, 3, 4 and figures 1, 2) 4 The more the cylindrical panel is thin the more the value of critical stress intensity
7, is small (see table 1, 2, 3, 4)
‘This paper is completed with financial support from the National Basic Research Program in Natural Sciences
References
Dao Huy Bich, Theory of clastoplastic processes, Vietnam National University Pub-
lishing House, Hanoi 1999 (in Vietnamese)
Volmir A S Stability of deformable systems, Moscow 1963 (in Russian) Ulo Tepik, Bifurcation analysis of elastic-plastic cylindrical shells, Int Journal of
Non-linear Mech, 34(1999), 299-311,
4, W 'T Koiter, Buckling and postbuckling behaviour of a cylindrical panel under axial compression, Nat Luchtvaart labort Rep., No 476(1956), Amsterdam Dao Van Dung, On the stability problem of cylindrical panels according to the
theory of elastoplastic processes Proceeding of the Seventh National Congress on Mechanics, Hanoi, 18+20 December 2002, pp 141-150 (in Vietnamese)
3 Dao Van Dung, Solving method for stability problem of elastoplastic cylindrical
shells with compressible material subjected to complex loading processes, Vietnam
Journal of Mechanics, NCST of Vol 23, No 2(2001), pp 69-86