Models of deformation of stiffened orthotropic shells under dynamic loading Текст научной статьи по специальности «Строительство и архитектура»

УДК 539.3, 531.3, 001.891.573

Models of Deformation of Stiffened Orthotropic Shells under Dynamic Loading

Alexey A. Semenov*

Saint Petersburg State University of Architecture and Civil Engineering 2nd Krasnoarmeyskaya, 4, Saint Petersburg, 190005

Russia

Received 23.03.2016, received in revised form 02.08.2016, accepted 08.09.2016 Two models of deformation of reinforced orthotropic shells under dynamic loading are considered in this paper. One such model is in the form of equations of motion and another model is in the form of a system of ordinary differential equations. Mathematical models are based on the hypotheses of the Kirchhoff -Love theory of shells. They take into account the geometric nonlinearity, orthotropic material properties and reinforcement elements. All relations of the models are in general form, and they can be used for a wide range of structures (shallow shells of double curvature, cylindrical, conical, spherical and toroidal shells and panels, etc.). An important feature of the proposed model is the ability to introduce stiffeners both discretely and by the method of constructive anisotropy (MCA) in accordance with their shear and torsional rigidity. The second model is derived by applying the Kantorovich method to the functional of the total energy of deformation of a shell. The resulting initial value problem is easier to solve than the system of equations of motion in partial derivatives.

Keywords: mathematical model, shell, dynamic loading, orthotropy, geometric nonlinearity, the equations of motion, method of constructive anisotropy. DOI: 10.17516/1997-1397-2016-9-4-485-497.

Introduction

Since 40-ies the intensive development of the theory of shells has begun in connection with the technological demand. During the first three decades researchers dealt mainly with static problems in the theory of thin plates and shells. Since the 70s of the last century the emphasis in studies of plates and shells is shifted to dynamic problems. Above all, this is due to demands of the aviation and space technology. However, the study of the dynamic behaviour of structures is also essential for shipbuilding and civil engineering.

The most studied is the behaviour of single-layer isotropic shell structures under dynamic loading. In the last decade, structures fabricated of composite materials are of great interest. For example, such materials as carbon fibre reinforced epoxy resin, graphite-fiber reinforced and boron-fiber reinforced plastics are widely used. The behavior of shells fabricated of these materials has been insufficiently studied [1-12].

The main objectives of the study of shell structures under dynamic loading are to investigate their stability, strength and vibrations as evidenced by review articles and monographs [13-23].

* sw.semenov@gmail.com © Siberian Federal University. All rights reserved

One should mention the extensive review article of E. A.Kogan and A .A. Yurchenko [13] which is devoted to free and forced nonlinear oscillations of multilayered thin elastic plates and shells under periodic loading.

Most studies are devoted to the vibration of isotropic structures [1,4,7, 8,12, 24-30].

As for geometry of structures, the vast majority of studies are devoted to closed cylindrical shells [4,5,7,11,26,27,30-32] because it is connected with the solution of a number of very important practical problems, for example, the occurrence of dynamic periodic load caused by the internal fluid flow. Such problems can be considered as axially symmetric. Studies on conic [3,24] and the spherical [28, 33] closed shells have received little attention. Non-axisymmetric problems for closed shells [34] are also rarely addressed, although they are important in studying structures with form imperfections.

Thus, the study of strength, stability and nonlinear oscillations of reinforced orthotropic shells of revolution under dynamic loading is a vital issue.

The aim of this work is to construct a mathematical model of deformation of reinforced orthotropic shells under dynamic loading.

1. Mathematical models of deformation of reinforced orthotropic shells under dynamic loading

We consider an orthotropic shell supported by the ribs on the internal side of the shell. It is subjected to the influence of an external uniformly distributed lateral load q = q (x, y, t).

To form relations of mathematical model we first consider the expressions for the forces and moments. The unknowns are three displacement functions U = U (x,y,t) ,V = V (x,y,t) ,W = W (x, y, t). Using the Kirchhoff - Love model for orthotropic shells, we have

Nx —

Ei

1 — V12V21 En

N =_2_

Ny

1 — 1^12^21

Nxy — G12

Mx

(h + Fy) Yxy + 2SyX12

E1

(h + Fx) (£x + l^21£y) + Sx (X1 + V21X2) (h + Fy) (£y + V12£x) + Sy (X2 + V12X1)

(h + Fx) Yxy + 2SxX12

Nyx — G12

My

Mx

G

12

1 — V12V21 E2

1 — V12V21 ' h3

h3

Sx (&x + V21£y) + + Jxi (X1 + V21X2) Sy (£y + V12£x) + (h~. + Jy] (X2 + V12X1)

V 12

SyYxy + 2[j2 + Jy) X12

My

G

12

h3

SxYxy + 2 [j2 + Jx 1

X12

(1)

where h is thickness of the shell; E1,E2 are elastic moduli; G12 is the shear modulus; ¡i12,^21 is the Poisson's ratio; Fx, Fy, Sx, Sy, Jx, Jy are ribs characteristics which are described below; £x,£y,Yxy are strains that account for the geometric nonlinearity; x1,X2,X12 are curvature and torsion functions of the form

1 de1 1 dA 1 de2 1 dB 1 de2 1 de1 1 (dA dB \

X1 — AdX+ABdye2> X2 — Bdy + ABdXeu 2X12 — AdX + Bdy — ABydy 1 + dXey'

where 01

1 dW

A^X'

1 dW

B~oy

and A, B are Lame parameters.

Stiffeners can be introduced with the use of the method of constructive anisotropy [35] or discretely. When using the method of constructive anisotropy, we have

hri E (h°rj E hnrÄ rj F E hrj E 1 hri ^ ^ ) ü ;

..,»-- . . ^ , - . . . , ■ ^ . , ab b

3=1 V i=l / 3 = 1 i=l \ 3=1

Sx = > ' — + E (SOri _ E Si3nr3 ^ rj . Syy = j SOn + j ( _ j Sjnr3 \ Ti

n -, i m / 7 3 n -, ii \ m -, 3 n I -, i m

V-^ --t- + ^ h3r3 sr h 3rir3 \ r3. F V^ h3r3 + V^ I hiri sr

¿-T^ ~~a~j~à> F = b — + Ь T"

i=l b 3=1 \ i=l ab J 3 = 1 ^ \ b 3=1

n i m 3 n i3 m 3 n i m

:Sr+е 3_:3 *=E33+E(sr_E

i=l b 3=1 \ i=l ab J 3 = 1 ^ \ b 3 = 1

Jlri + m f J3 r3 ^ 33 r3 . T = m J3 r3 + ^ ( J't- m J3 rir3

^"T+ MT""l^Jä"' Jy = ^+ ^ "T

i=l b 3 = 1 \ i=l ab J 3=1 ^ \ b 3 = 1 äb

^ v^ hiri h3r3 v^ hi3rir3 , . , ~ f a — ai\

fp=e-j-+E -a- _ E -or- a=(a _ «> a b=ьв\—)

where h. is the rib height; indexes i and j indicate the number of rib located parallel to the x and y axis, respectively; n, m are numbers of ribs; hij = min {hl,hß } ; r is the rib width; a,, b are limit values along x, y axes; a1 is the offset from the origin along the x axis;

. hi (h + hi) . hj (h + hj) .. hij (h + hij) . 2 . ( ..2 1 ( .)3

^^-, Sj = —i-, =-i-, J1- = 0.25h h1+0.5h(hM2 + -(h?)3 ,

Jj = Q.25h2hj + 0.5h (hj)2 + 3 (hj)3 , J.j = 0.2bh2hP + 0.5h (tij)2 + - ft*)3 .

The method of constructive anisotropy can be used if the ratio of distance between ribs to shell length is not more than 0.1.

When ribs are introduced discretely, we have

nm

Fx = Fy = Fp = J2 hj 5 (x - xj) + J2 hi $ (y - Vi) ^ * (x - x3) * (v - Vi)'

j=1 i=1 i=1j=1 m n n m

Sx = Sy = j2 Sj 5 (x - Xj) +y, Si5 (y - Vi) -^y, Sij 5 (x - Xj) 5 (y - Vi), j=1 i=1 i=1j=1 m n n m

Jx = Jy =53 Jj5 (x - xj) + y Ji5 (y - Vi) - X) y Jij5 (x - xj) 5 (y - Vi),

j=1 i=1 i=1 j=1

where S (x - xj) and S (y - yi) are differences of two Heaviside functions: S (x - xj) =

U (x - aj) - U (x - bj); S (y - yi) = U (y - ci) - U (y - di), where aj = xj - rj/2, bj = xj + rj/2, ci = yi - ri/2, di = yi + ri/2.

In the presence of reinforcement, rib thickness is defined as

m n n m

H (x, y) = Y hj 5 (x - xj) hiS (y - Vi) - Y^ YI hij S (x - xj) S (y - yi),

j = 1 i=1 i=1 j=1

then the thickness of the entire structure is h + H.

The total deformation energy functional of the shell under dynamic loading is

I = ! 1 (K - Ep) dt, (2)

Jt0

Jx

where K is the kinetic energy of deformation; Ep = n — A is the difference of potential energy of deformation of the system and the work of external forces that corresponds to the functional of the static problem; t is time.

The kinetic energy of deformation of reinforced shell has the form [35, 36]

K = P

Jo J-h/2 \dt) + V+ V)

la! J 0

ABdxdydz,

(3)

where p = — (y,g are the specific gravity of the shell material and the gravity acceleration).

g

Upon integrating (3) with respect to z, we obtain

K = Pp

¡u: {"•+4 f)2 +( & +( w

ABdxdy.

(4)

Next, we consider two models. First model is in the form of equations of motion and second model is in the form of a system of ordinary differential equations.

1.1. Equations of motion

Since the behavior of the system in the time interval [t0,ti] must satisfy the condition [36]

A / (K — n + A) dt = 0,

(5)

we seek the first variation of functional (2), and then we set it equal to zero. First, we find the variation of the kinetic energy

sl'Kdt=pI0 £/ +# ik

8UfdU\ 3VidV\ 3WfdW\ >-A - + 2-A - + 2-A -

2mA{7t + 2 ~dtA\m) + 2^\otl

ABdxdydt.

One should change this expression such that there are no variations of the time derivatives of the unknown functions U(x, y, t), V(x, y, t), W(x, y, t) with respect to t under the integral sign. Upon integrating by parts, we obtain

p ti p ti p a pb

A / Kdt = —p (h + Fp)

J t0 J t0 J ai J 0

o2u d2 u d 2w^7t;

AU + AV + —^77;— AW

dt2

dt2

+P P i6 {(h + Fp)

a10

dU STT 3Vr dW

ltAU + ltAV + ~wAW

dt2

t=ti

t = to

ABdxdydt+ ABdxdy.

The second term in this equation allows us to formulate initial conditions, specifically, either all functions are equal to zero or first derivatives of these functions are equal to zero at the initial moment of time.

The first term in this expression, from which the equations of motion are formed (or rather, the inertial terms of these equations) can be reduced to the form

2

2

t

d 2U ,, „,32V„r ,, „„8 2W

/ti pa pb r

.1 Jo l(h + Fp) AU +(h + Fp) d^TAV + (h + Fp) ^AW j ABdxdy dt.

Since potential energy Ep does not depend on derivatives of the unknown functions with />ti />ti

respect to t, we have A I Epdt = AEpdt. Taking into account (1), functional Ep has the

t0 t0

form [35]

Ep = 1 Ja j0 {Nx£x + Ny£y + \(N*y + ».-) Yxy + MX + MyX2 +

+ (Mxy + Myx) X12 — 2qW} ABdxdy.

Let us find the first variation of this functional and set it equal to zero. Then we carry out transformation in order to exclude derivatives of displacements under the variation sign. As a result, we obtain

dA N

dy 5

ab

ai J0

AV +

dBNx dB dANxy

-X — Ny dx +-

dx

dy

+ ~oyNxy

AU +

BAN,,

+

dy

d (B3W \ d (A dW \ AB k nx + ky Ny) + -^{a^x) + dy \BlyNyj

dBNxy + » 3B__ + »xy dx

dx Ny | +

d few » dy\~dxNxv

+

d (dW„ \ d2 [B\ d2 /A \

-l — N~y) + ox2{amx) + oy? \ BMy)

+ 3x^ dy »xy

d2Mxy d f l dB \

+ 2_—____M —

dxdy dx \A dx yJ

AW^ dxdy+

+

+

a1 i: {

d (l dA \ — ~ \B~dyMx) + ABq_

ANxy AU + ANy AV +

dy \B dy

AdW ^ dW l dA d (A

Ely» + ~dxNxy — BdyMx + d~y[B^

{am,)

+

dM,

xy

dx

AW

-Mxy A

xyK'

fdW\ A r fdW\

\Jw) — BMy AU")

y=b

dx+

y=o

BNxAU + BNxy AV +

BdW

dW

l dB

d fB

A dx Nx + — ATxM' + di(AM'

(am-)

+

dMx

dy

AW

B „ r idW\ _ r fdW\ — -M- a( dx) — Mxy a i - '

A

V dy J

dy = 0.

Neglecting transverse shear and rotational inertia, the equations of motion for orthotropic reinforced shells take the form

(h+f) g (h+f) g

, „ dBNx dB dA dANxy AB = — Ny TT + irN-y + ;

dx dx dy dy

AB =

dANy dBNr

dy

+

dB dA^T a + N-y^~ — ^ N-;

dx dx dy

xy

p (h + F) ?2W AB = AB (k-N- + ky Ny) + KB^N.) + MB^». I +

(BdW \ d (AdW \ dx\A^xNx) + &y\B^dyNy)

(6)

+

d (BW„ \ d (BW„ \ d2 (B\ d2 /A \

d-y{lxNxy) + d~x{-WNxy) + AMx) + dy2 \BM.) d2Mxy d ( l dBr ) d ( l dA

+

+2

xy (l dB \ d ( l dA \

dxy— dx{A -dxMy)— ^{b -d~yMx)+ABq.

System (6) is a system of differential equations in partial derivatives. Solution of this system presents difficulties.

x=a

x=ai

p

p

1.2. Initial value problem for the ODE system

In what follows we obtain the initial value problem for the system of ordinary differential equations using the Kantorovich method [37] to transform three-dimensional functional of the total energy of deformation of the shell into the one-dimensional functional. First, we introduce dimensionless parameters

x — ai y - aA

£ =-, V =T , A , h = hkx

a — ai b bB

k>n — hky,

U = aUA - = bVB

h2

h2

W W = W

F = Fp >p = h'

a2 A2\j (

Ei

■t, P

a4A4q

'a2 A2\j (1 — ni2 Mi) P ' h4Ei '

Upon substituting these parameters into (2), we obtain

A = aA B = bB - = z h ' h ' h

(7)

rh -

I= K — Ep) dt, where K

J to

'¿(¿-pu^OA^^ F 2(1 — H12V21) a4 A4

Eh K, EP = Eh Ep,

and

k=(i+>4 % î+*( dv î+Ai w

A Bd£dn.

This procedure was considered in details in [38]. Functional I can be expressed as

I = S 1 So So $ 1 UVW) — $2 (U, -, W) ABd^dr/dt.

(8)

Here $ i corresponds to K, $2 corresponds to Ep. The dot indicates derivative with respect to time.

According to the Kantorovich method, unknown functions can be represented as

N

N

N

U = J2 Ui(t)X 1(£)Y 1(n), V = J2 Vi(t)X2(£)Y2(n), W = £ Wi(t)X3(£)Y3(n), (9)

i=i

i=i

i=i

where Ui(t),Vi,Wi(t) are unknown functions, and all other functions are known. Examples of such functions for various boundary conditions are shown in Tab. 1 and variables

l = 1VN,p = 1,..., Vn.

Table 1. Approximating functions for various boundary conditions

Type of boundary condition Function At £ = 0, e =1 At n = 0, n =1

Simple support U (£,n) ,V (£,n) ,W (£,n) sin (pn£) sin (lnn)

Clamp U (£,n) ,V (£,n) sin (pn£) sin (lnn)

W (£,n) sin2 (pn£) sin2 (lnn)

It should be noted that approximating functions are provided for the case of an open asymmetrical structure. They can be applied, for example, for conical panels. In the case of symmetrical structure approximating functions are considered in [39].

2

We take into account (9), when differentiating $ 1 and $2.

Let us take variation of functional (8). The condition of its minimum is

SI = 5u I + 5y I + Sw I = 0,

that is, we obtain the system

N t r 1 r 1 Sui = £ ! /

i=1J t0 J0 J0 N rti r 1 r 1 svi = £ / /

i=1Jt0 Jo Jo N ft 1 r 1 r 1 W = £ t /

i=1Jt0 Jo Jo

^ SUt - ^ 5Ul dUi l dUi l

^ * - # sv]

dVi l dVi

<9$ 1 ■ d$2 4- SWi - —2 SWi

dWi

0W,

A BdÇdndt = 0, ÀBdÇdndt = 0, ABdÇdndt = 0, l =1, 2,... ,N.

(10)

The sum sign used in (9) is factored out from the integral signs in relations (10). Moreover, unknown functions with both index i and index l are in $\ and $2. This is shown below in expanded form.

Upon integrating by parts and factoring out SUl, SVl, SWl from each equation, we obtain the equations of motion

N 1 1 £ / ' d io$ 1 dt y dUi

N .1 .1 l=j0 j0 ' d fd$ 1 dt y dV

N .1 .1 i=1 o o ' d id $ 1 dt \vdWWi

+

+

+

d$ 2 dUi

d$ 2 ' dVi

d$ 2 dWi

ABdd = 0, ABdd = 0,

ABd^dn = 0, l = 1, 2,... ,N.

Now we can calculate integrals of known approximating functions with respect to £, n. Here

N 1 1 r 1 N 1 1 r 1 N 1

tUL\tlAB*n UL

=i

=i '

d$ 2 dWi

ABd£dn are the

left-hand sides of the equations of statics. They can be obtained with the use of the Ritz method applied to functional Ep [40]:

oo

N

E

i=1 N

i= o o N

i= o o

d$ 2 dUi

d$2 '

dVi

d$ 2 dWi

N

ABdd = Y, [UiC 1 (i, l) + ViC2(i, l) + WiCs(t, l)] + A 1(l) = 0,

i= N

ABdidn = Y, [UiC6(i, l) + ViCr(i, l) + WiC8(i, l)] + Â2(l) = 0,

(11)

N

ABd^dn = J2 [UiC11(i, l) + ViC12(i, l) + WiC13(i, l)] + Â3(l)-

i=

-CP (l)P = 0, l =1,2,...,N.

Here C1,..., C13, Cp are coefficients, A1 (l), A2(l), A3(l) are nonlinear terms of equations.

The remaining terms of equation (10) are transformed as follows

N

T-^d $ 1 _ 1

= Wt _ A

1=1 i

N

(1 + Fp) 2 ^ UiX1 (i) Y1 (l) ■ X1 (l) Y1 (l)

i=i

M =A

A2

=1 ^{m)

E f1 f1 1

f={Jo Jo dt \K3Ul )

N

(1 + FP)J2 UiX 1(i) Y1 (i)

i=1

X 1(l) Y 1(l),

N

ABdd = 53 [uiB1(i,l)

i=1

After applying similar transformation to other equations of system (10), we obtain

E f f\d i^V

i^Jo Jo [ dt ( dVi )

E f1 f1 i d ()

i=1 Jo Jo dt ydWi J

ABdd [VUiB3 (i, l)

i=1 N

ABdd = 53 WWiB5(l,l).

Here

B1(l,l)_ f f ti (1 + Fp) X 1(i) Y 1(i) X 1(l) Y 1(l) ABd^dn, Jo Jo A

n1 2\2

(1 + Fp) X2 (l) Y2 (l) X2 (l) Y2 (l) ABd^dn,

. A

B5(l,l)_ [ [ 2(1 + Fp) X3(l) Y3(l) X3(l) Y3(l) ABd£dn-

oo

Thus, the equations of motion of a stiffened orthotropic shell has the form

N

]T \viB1(l,l)\ + n (l)_o,

i=1 N

J2 \ViB3(l,l)\ +n (l)_0,

i=1 N

J2wUiB5(l,l) + n3 (l)_0, l _1, 2,...,N,

(12)

where n1(l),n2(l), n3(l) are the left-hand sides of system (11).

Thus, the mathematical model of deformation of reinforced orthotropic shells is reduced to the initial value problem for a system of ordinary differential equations. Let us reduced system of ordinary differential equations (12) to the normal form. We introduce

Then system (12) can be written in the form

N

^ 1+n (l)_0,

J2[Y1iB1(l,l)

i=1 N

J2 \Y2iB3 (l, l)

+ n (l)_0,

i=1

J2N=1 ^3i B5 (l, l) + n (l)_0, l _1, 2,..., N.

Hi _ Y1i, V _ Y2i, Wi _ Y3i.

The matrix of coefficients of system (13) can be schematically represented in the form

B1 0 0

B=l 0 B3 0

1 0 0 B5

The dimension of this matrix is 3N x 3N. Now system (13) can be written as

BY (t)_n(t,X),

where Y(t) _ (Y1, Y2, Y3)T, n _ ( - n1(l), -n2(l), -n3(l))T, X _ (Ui,Vi,Wi)T, l _ 1,...,N,l _1,...,N.

Matrix n is the column matrix of the right-hand sides of system (6). Let us find the inverse of matrix B. Each element of the matrix B is known. Then we obtain

B^By (t) _ B^n (t, X)_ F (t, X).

Thus we have a normal system of ordinary differential equations

Y (t) _ F (t,X), Ui _Y1i, Vi _Y2i, Wi _Y3i, l _1,...,N,

which can be represented in the form X (t) _ Y (t).

This system is solved with the initial conditions at t _ 0

Ui _ Vi _ Wi _0, Ui _ Vi _ Wi _0, l _1,...,N. (14)

Coefficients of the matrix B and coefficients in the expressions for n are obtained with the use of boundary conditions. In what follows we assume that q _ A ■ t.

Upon substituting initial conditions (14) into the column matrix n with 3N elements, N elements of the matrix are not equal to zero. This is due to the presence of a load term that depends on t.

To solve the initial value problem for systems of ordinary differential equations (11), (12) various numerical methods can be used. Since they are stiff systems one can use, for example, an implicit Runge-Kutta scheme.

2. Numerical results

The initial value problem for the ODE system is solved with the use of computing software Maple 2015. The system is numerically solved using an implicit Runge-Kutta scheme.

As an example, we consider a shallow shell of double curvature (A _ B _ 1, kx _ ky _ 1/R). The contour of the shell is simply supported. The shell is subjected to the dynamic external uniformly distributed lateral load q _ q (x,y,t) _ 100t. The dimensional parameters of the shell are a _ b _ 60h, R1 _ R2 _ 225h, h _ 0.09 m. The shell is made of steel with the modulus of elasticity E _ 2.1 ■ 105 MPa and Poisson ratio n _ 0.3. The system is solved in terms of dimensionless parameters and N _ 9.

Fig. 1 shows the relationship between load P and displacement Wc at the center of the shell. An inflection point on the curve "load-displacement" is the point of stability loss under dynamic loading. Such criterion was proposed by A. S. Vol'mir. It is seen from Fig. 1 that after buckling load P _ 232 shell oscillates around the equilibrium state of the static problem.

P

-« 4<> ® « « t 5 » o « »»4 <, O o 0 o o ° o ® O 4 •> <> 8®

J J«»» > " •> 0 fl o t, O O 4 o O ® ® » » O » » o o 4

J (

/

/

0 12 3 4 5 6

Fig. 1. Relationship between load P and displacement Wc at the center of the shell

3. Conclusions

The mathematical model of the deformation of stiffened orthotropic shells under dynamic loading was presented. The model is represented in the form of an initial value problem for a system of ordinary differential equations that is easier to solve than the system of equations of motion in partial derivatives. The proposed model takes into account the geometric nonlinearity, and can be applied to orthotropic materials. Stiffeners can be introduced as discrete elements or with the use of the method of constructive anisotropy.

The work was supported by the Ministry of Education and Science of the Russian Federation (project 3801).

References

[1] H.Abramovich, V.A.Zarutskii, Stability and vibrations of nonclosed circular cylindrical shells reinforced with discrete longitudinal ribs, Int. Appl. Mech., 44(2008), no. 1, 16-22.

[2] V.F.Meish, A.S.Kairov, Vibrations of reinforced cylindrical shells with initial deflections under nonstationary loads, Int. Appl. Mech., 41(2005), no. 1, 42-48.

[3] N.A.Shul'ga, S.Yu.Bogdanov, Forced axisymmetric nonlinear vibrations of reinforced conical shells, Int. Appl. Mech., 39(2003), no. 12, 1447-1451.

[4] S.G.Suleymanova, Free vibrations of a filled cylindrical shell longitudinally strengthened and loaded with axial contracting forces, Proc. of IMM of Azer, XXVII(2007), 135-140.

[5] G.I.Belikov, S.Yu.Kalashnikov, Layout operation of stiffened cylindrical shells of the bridge conduits under free oscillations, Bulleten Volgograd. Gos Univ. Arhitectury i Grazhd. En-giner. Ser.: Civ. Eng. and Arch., (2011), no. 25(44), 14-20 (in Russian).

[6] Yu.P.Dyachenko, E.J.Elenitskiy, D.V.Petrov, Non-stationary problems of the dynamics of stepped section plates and rotation cylindrical shells, Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, (2011), no. 2(23), 278-288 (in Russian).

[7] F.S.Latifov, S.G.Suleymanova, A problem of free vibrations medium-filled cylindrical shells reinforced by a cross system of ribs and loaded with axial contracting forces, Mehanika mashin, mehanizmov i materialov, (2009), no. 1, 59-61 (in Russian).

[8] M.A.Mehtiyev, Nonlinear parametric vibrations of stiffened cylindrical shell with a vis-coelastic filler, Mehanika mashin, mehanizmov i materialov, (2011), no. 3(16), 28-30 (in Russian).

[9] N.A.Nazarov, Oscillations of shallow shells, reinforced by stiffening ribs, Applied Mechanics, 1(1965), no. 3, 24-31.

[10] Yu.I.Nemchinov, Yu.A.Talbatov, Free oscillations of shallow cylindrical shells reinforced by stiffening ribs, Struct. Mehan. i Anal. Konst., (1975), no. 3, 17-22 (in Russian).

[11] A.I.Seyfullayev, K.A.Novruzova, Oscillations of a longitudinally reinforced orthotropic cylindrical shell filled with a viscous fluid, Eastern-European Journal of Enterprise Technologies, (2015), no. 3/7(75), 29-33 (in Russian).

[12] G.V.Tertishniy, Effect of stiffeners on the frequency of free oscillations of a shallow cylindrical panel, Bulleten Kazan. Tehnolog. Universiteta, (2011), no. 19, 217-224 (in Russian).

[13] E.A.Kogan, A.A.Yurchenko, Nonlinear oscillations of a three-layer and multi-layer plates and shells during periodic impacts (survey), Izvestiya MGTU MAMI. Seriya estestv. nayk, 4(2014), no. 1(19), 55-70 (in Russian).

[14] V.D.Kubenko, P.S.Koval'chuk, Nonlinear problems of the dynamics of elastic shells partially filled with a liquid, Int. Appl. Mech., 36(2000), 421-448.

[15] V.D.Kubenko, P.S.Koval'chuk, Nonlinear problems of the vibration of thin shells (review), Int. Appl. Mech., 34(1998), no. 8, 703-728.

[16] M.Amabili, M.P.Paidoussis, Review of studies on geometrically nonlinear vibrations and dynamics of circular cylindrical shells and panels, with and without fluid-structure interaction, Appl. Mech. Rev., 56(2003), 349-381.

[17] F.Moussaoui, R.Benamar, Non-linear vibrations of shell-type structures: a review with bibliography, Journal of Sound and Vibration, 255(2002), 161-184.

[18] M.S.Qatu, R.W.Sullivan, W.Wang, Recent research advances on the dynamic analysis of composite shells: 2000-2009 (Review), Composite Structures, 93(2010), 14-31.

[19] S.K.Sahu, P.K.Datta, Research advances in the dynamic stability behavior of plates and shells: 1987-2005 — part I: conservative systems, Appl. Mech. Rev., 60(2007), no. 2, 60-65.

[20] N.A.Abrosimov, V.G.Bazhenov, Nonlinear problems of composite structural dynamics, Nizhni Novgorod, Izd. NNGU, 2002 (in Russian).

[21] I.Ya.Amiro, V.A.Zarutskii, V.G.Palamarchuk, Dynamics of ribbed shells, Kiev, Naukova dumka, 1983 (in Russian).

[22] Y.Kimihiko, Review of research in Japan on nonlinear oscillations of elastic structures, ISME Int. Journ. C, 39(1996), no. 3, 439-449.

[23] A.N.Blinov, On the Lower Critical Load of the Elastic Cylindrical Shell with Axial Compression, Journal of Siberian Federal University. Mathematics and Physics, 5(2012), no. 3, 359-362 (in Russian).

[24] T.Ueda, Non-linear free vibrations of conical shells, Journal of Sound and Vibration, 64(1979), no. 1, 85-95.

[25] M.Amabili, R.Garziera, R.Muharlyamov, K.Ryabova, Nonlinear oscillations of shallow shells of double curvature, Vestnik Kazanskogo Technolog. Universiteta, 18(2015), no. 6, 158-162 (in Russian).

[26] V.M.Dubrovin, T.A.Butina, Modeling of the dynamic stability of a cylindrical shell under the axial compressive load, Mat. Mod. Chisl. Met., (2015), no. 2(6), 46-57 (in Russian).

[27] G.I.Kolosov, Vibration and Secular Instability of Equilibrium States of Cylindrical Shells under Axial Compression, Cosmonavtika i Raketostroenie, (2012), no. 2(67), 145-150. (in Russian)

[28] V.V.Platonov, The stability of the transversely isotropic spherical shell under the normal dynamic loading, Vestnik St.-Petersburg Univ. Ser.1, (2010), no. 3, 105-110 (in Russian).

[29] L.Yu.Stupishin, A.G.Kolesnikov, I.V.Solomatnikov, Geometric nonlinear shallow shells optimum forms examination for maximum lowest frequencies of small free oscillations, Trudy Yugo-Zapadn. Gos. Univer., (2011), no. 5-2, 313-316 (in Russian).

[30] N.A.Taranuha, G.S.Leyzerovich, Non-linear free flexural oscillations thin circle cylindrical shells, Far Eastern Mathematical Journal, (2000), no. 1, 102-110 (in Russian)

[31] W.W.Bendukov, W.W.Derushev, M.M.Lurie, P.N.Owcharov, About influence of filler on critical parameters of an impulse of pressure at dynamic loss of stability of a cylindrical envelope, Nauchn. bulleten Moskov. Gos. Tehnich. Univer. Grazhdan. Aviatsii, (2005), no. 84(2), 131-137 (in Russian)

[32] A.V.Krysko, V.A.Krysko, N.E.Saveleva, Chaotic vibrations of closed cylindrical shells and plates. Part I, Vestnik SSTU, (2005), no. 3(8), 32-61 (in Russian).

[33] D.Li, A time-mode approach to nonlinear vibrations of orthotropic thin shallow spherical shells, Int. J. Solids Structures, 30(1993), no. 22, 3113-3128.

[34] V.A.Krysko, I.V.Kravtsova, Dynamics and statistics of sectorial shells, Vestnik SSTU, (2004), no. 2(3), 27-36 (in Russian).

[35] V.V.Karpov, O.V.Ignat'ev, A.Yu.Salnikov, Nonlinear mathematical models of deformation of shells of variable thickness and algorithms for their research, Moscow, ASV; SPb, SPbSUACE, 2002 (in Russian).

[36] A.S.Vol'mir, Nonlinear dynamics of plates and shells, Moscow, Nauka, 1972 (in Russian).

[37] L.V.Kantorovich, A direct method for the approximate solution of problem of the minimum of a double integral, Trudy Akad. Nauk SSSR, Matematich. i estest. nauki, (1933), no. 5, 647-652 (in Russian)

[38] V.V.Karpov, A.A.Semenov, Dimensionless parameters in the theory of reinforced shells, PNRPU Mehanika Bulleten, (2015), no. 3, 74-94 (in Russian)

[39] V.V.Karpov, The strength and stability of reinforced shells of revolution. In two parts. Part 1. Models and algorithms of research of the strength and stability of supported shells of revolution, Moscow, Fizmatlit, 2010 (in Russian).

[40] A.A.Semenov, Algorithms for the research of strength and stability of reinforced orthotropic shells, Struct. Meh. Eng. Konst. i Sroit., (2014), no. 1, 49-63 (in Russian).

Математические модели деформирования подкрепленных ортотропных оболочек при динамическом нагружении

Алексей А. Семенов

Санкт-Петербургский государственный архитектурно-строительный университет

2-я Красноармейская, 4, Санкт-Петербург, 190005

Россия

В данной 'работе 'рассматриваются два варианта математической модели деформирования подкрепленных ортотропных оболочек при динамическом нагружении: в виде уравнений движения и в виде системы обыкновенных дифференциальных уравнений. Математические модели строятся на основе гипотез теории оболочек Кирхгофа-Лява, учитывают геометрическую нелинейность и ортотропию материала, а также возможность наличия подкрепления конструкции. Все соотношения моделей даются в общем виде и при указании соответствующих параметров Ляме могут быть использованы для широкого класса различных конструкций (пологих оболочек двоякой кривизны, цилиндрических, конических, сферических и тороидальных оболочек и их панелей и др.). Важной особенностью предложенной модели является возможность введения ребер жесткости как дискретно, так и по методу конструктивной анизотропии с учетом их сдвиговой и крутильной жесткости. Второй вариант математической модели выводится путем применения к функционалу полной энергии деформации оболочки метода Л. В. Канторовича (метод сведения трехмерного функционала к одномерному). Полученная начальная задача решается существенно проще, чем система уравнений движения в частных производных.

Ключевые слова: математическая модель, оболочки, динамическое нагружение, ортотропия, геометрическая нелинейность, уравнения движения, метод конструктивной анизотропии.