Безразмерные параметры в теории подкрепленных оболочек Текст научной статьи по специальности «Строительство и архитектура»

Карпов В.В., Семенов А.А. Безразмерные параметры в теории подкрепленных оболочек // Вестник ПНИПУ. Механика. - 2015. - № 3. - С. 74-94. DOI: 10.15593/peim.mech/2015.3.07

Karpov V.V., Semenov A.A. Dimensionless parameters in the theory of reinforced shells. PNRPU Mechanics Bulletin. 2015. No. 3. Рр. 74-94. DOI: 10.15593/perm.mech/2015.3.07

ВЕСТНИК ПНИПУ. МЕХАНИКА № 3, 2015

PNRPU MECHANICS BULLETIN

http://vestnik.pstu.ru/mechanics/about/inf/

DOI: 10.15593/perm.mech/2015.3.07 УДК 539.3

БЕЗРАЗМЕРНЫЕ ПАРАМЕТРЫ В ТЕОРИИ ПОДКРЕПЛЕННЫХ ОБОЛОЧЕК В.В. Карпов, А.А. Семенов

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

О СТАТЬЕ

АННОТАЦИЯ

Получена: 23 апреля 2015 г. Принята: 17 июля 2015 г. Опубликована: 30 сентября 2015 г.

Ключевые слова: оболочки,

безразмерные параметры, математическая модель, подкрепленные оболочки, параметры Ляме, ортотропия, оболочки вращения, устойчивость, безразмерная нагрузка, безразмерный функционал

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

© ПНИПУ

© Карпов Владимир Васильевич - доктор технических наук, профессор, e-mail: [email protected] Семенов Алексей Александрович - кандидат технических наук, e-mail: [email protected]

Vladimir V. Karpov - Doctor of Technical Sciences, Professor, e-mail: [email protected] Alexey A. Semenov - Ph.D. in Technical Sciences, e-mail: [email protected]

74

Карпов В.В., Семенов А.А. /Вестник ПНИПУ. Механика 3 (2015) 74-94

DIMENSIONLESS PARAMETERS IN THE THEORY OF REINFORCED SHELLS

V.V. Karpov, A.A. Semenov

Saint Petersburg State University of Architecture and Civil Engineering, St. Petersburg, Russian Federation

ARTICLE INFO

ABSTRACT

Received: 23 April 2015 Accepted: 17 July 2015 Published: 30 September 2015

Keywords:

shell, dimensionless parameters, mathematical model, reinforced shells,

Lame parameters, orthotropy, a shell of revolution, stability, dimensionless load, dimensionless functional

In this paper a variant of the dimensionless parameters for a wide range of shell structures is proposed. For shallow shells of have a rectangular form, the dimensionless parameters are used for a long time, as for the general form of shells; there is no single form of dimensionless ratios for each type of shell, because Lame parameters differ not only in values, but also in dimensions. Therefore, the work shows the dimensionless relations for deformations, stresses, forces, moments and functional of total potential energy of deformation. The relations consider geometric nonlinearity, transverse shifts, orthotropy of material and the introduction of ribs on the structural anisotropy of the method in accordance with their shear and torsional rigidity. The authors show a further approach to solving the strength and stability tasks of different types of shells in the dimensionless parameters. Some methods for solving nonlinear problems of stability do not look quite correct, when using the dimensional parameters (e.g., a technique, based on the method of continuation of solution for the best parameter). In dimensionless parameters, all calculations are beyond doubt. The calculations of some shell structures are provided in dimensionless and dimensional parameters and their consistency is shown. The proposed approach allows one to obtain values for the calculation of a series of such shells, which is more convenient to optimize the choice of design parameters. We have examined some of the differences in critical loads, obtained in dimensionless and dimensional solution of the problem.

© PNRPU

Введение безразмерных параметров позволяет одним расчетом получить напряженно-деформированное состояние целой серии подобных оболочек и выбрать наиболее рациональные размерные параметры оболочек, а для изотропных оболочек еще и материал оболочки.

При расчете пологих изотропных оболочек прямоугольного плана безразмерные параметры использовались в работах [1-8]. Для оболочек вращения в осесимметричной постановке безразмерные параметры использовались в работе [9], что является частным случаем параметров, предлагаемых в данной работе.

В работах [10-17] используются различные безразмерные параметры для каждого конкретного вида оболочек и решаемых задач. Однако целесообразно разработать единые безразмерные параметры, пригодные для различного вида оболочек независимо от класса решаемой задачи.

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

1. Теория и методы

1.1. Основные соотношения математической модели

1.1.1. Соотношения для гладких оболочек

Математическая модель деформирования оболочек состоит из трех групп соотношений:

- геометрические соотношения, связывающие деформации и перемещения;

- физические соотношения, связывающие напряжения и деформации;

75

Karpov V.V., Semenov A.A. / PNRPU Mechanics Bulletin 3 (2015) 74-94

- функционал полной потенциальной энергии деформации оболочки, из условия минимума которого выводятся уравнения равновесия.

Геометрические соотношения в срединной поверхности оболочки с учетом геометрической нелинейности принимают известный вид [18]

1 dU 1 т rdA

S X =

+ Ут7 - kxW +102,

A dx AB dy

1 dV 1 SB 1 2

s y =--+---U----kW + -02

B dy AB dx

2

1 dV 1 dU 1 T TdA 1 SB

У xy =-— + -—- — U~ ~ V~ + 0102,

A dx B dy AB dy AB dx

(1)

01 =-

sw_

A dx

+ kU

02 =-

1 dW , T.

-----+ kyV

B dy y ,

Л

где sx, s - деформации удлинения вдоль координат x, y срединной поверхности; уxy -

xy

деформации сдвига в плоскости XOY; kx = —, k =— - главные кривизны оболочки

R

R,

вдоль осей x и y ; R1, R2 - главные радиусы кривизны вдоль осей x и y ; A, B - параметры Ляме, характеризующие геометрию оболочки; U = U ( x, y), V = V (x, y), W = W ( x, y) -перемещения точек срединной поверхности оболочки вдоль осей x, y, z.

Будем учитывать поперечные сдвиги, тогда

Уxz = kf (z)

* x -01

Уyz = k f (z)

* y-02

(2)

Здесь f (z) - функция, характеризующая распределение напряжений Txz и тyz по толщине оболочки; *x = *x (x, y), *y = *y (x, y) - углы поворота нормали в плоскостях XOZ

и YOZ соответственно; k = 5.

6

Для гладких оболочек [11]

f (z ) = 6

^ z2 ^

v 4 h у

где h - толщина оболочки.

Физические соотношения при линейно-упругом деформировании для ортотропного материала при плоском напряженном состоянии будут иметь вид [19]

а = E S x 21S y + z ( 21Х 2 )

а x 1 1 — ^12^21

а = e2 Sy + ^12Sx + z (2 +^12X1)

а у 1 1 — ^12^21

(3)

76

Карпов В.В., Семенов А.А. /Вестник ПНИПУ. Механика 3 (2015) 74-94

т xy = G12

У xy + 2 ^Xl2

^xz G13Уxz ; ^yz G23y yz *

Здесь E1,E2,-12,-21,G12,G13,G23 - механические характеристики материала, а %1,Х2,Х12 -функции изменения кривизн и кручения,

Х1 =

1 aw x 1 ал

a ax лв ay

wу ; Х2 =

1 aw у 1 ав

в ay лв ax

w.

2Х12 =

1 aw y 1 aw x 1 (ал„, ав„, Л

A dx в ду Ав V ду

—w_ +—w,

ax

Функционал полной потенциальной энергии деформации гладкой ортотропной оболочки можно представить следующим образом:

Ep = . Eh . [jfe +G,s; + 2ци8,Eу + G12У2, + G,k(w, -01 )2

2(1 M-12^21.

^1 0

G13k(wу - 02)2 + ^(x? + G2X2 + 2Ц21Х1Х2 + 4G12X122) - lABdxdy*

q(1 M-12^21 X

(4)

12 V- — ■ —— — E1h

С учетом того что для ортотропного материала Ey21 = E2-12, здесь введены следующие обозначения:

G = E2 G = G12 (1 М-12^21. G = G13 (1 --12-21. G = G23 (1 -Ml2-21X G2 = e , g12 = e , G13 = e , G23 = E

Выражения усилий и моментов для гладких ортотропных оболочек имеют вид [20]

E1

N = ■

1 ^12^21

Mx =

( Л! E2 ( Л!

h 8 x + -218 у I N = 2 ’ / 1 h 8 у + -128 x I

V JJ 1 - —12—21 V JJ

, Nxy = Nyx = G12 hy xy ,

E

1 ^12^21

Mxy = Myx = G12

( h3 Л

V12 J

(Х1 + - 21Х 2)

E

My = Г-21 —12—21

(h3 л

V12 J

(2+-12Х1)

(5)

(h3 л

V12 J

Х12

Q, = Gukh(wx -01), Qy = Gr^kh(w, -02),

где Nx, Ny, Nxy, Nyx - нормальные усилия в направлении осей x, у и сдвиговые усилия в плоскости XOY соответственно; Mx,My,Mxy,Myx - изгибающие моменты в направлении осей x, у и крутящие моменты; Qx, Qy - поперечные силы в плоскостях XOZ и YOZ *

1.1.2. Соотношения для подкрепленных оболочек

Для учета подкрепления конструкции ребрами жесткости при достаточно большом числе ребер будем использовать метод конструктивной анизотропии [21]* Суть этого метода заключается в сведении оболочки дискретно-переменной толщины к оболочке рав-

77

Karpov V.V., Semenov A.A. / PNRPU Mechanics Bulletin 3 (2015) 74-94

новеликой по жесткости оболочке постоянной толщины, что дает возможность учитывать такие существенные факторы, как сдвиговая и крутильная жесткость ребер. Также этот метод может быть модифицирован для расчета оболочек, ослабленных вырезами.

Для оболочки, подкрепленной ребрами жесткости, f (z) принимает вид [21]

f (z ) = -

(h + H)

2>z+2 Iz -2-HJ.

HI II ll HI

H (x y ) = I h s(x - xj )+I h - y )-II h'J s(x - xj )(y - y)

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

где h - высота; индексы i и j указывают номер ребра, расположенного параллельно оси x и y соответственно; n, m - количество ребер; hiJ = min {hi, hJ J ; б(x - xj ) и 8(y - yi)

представляют собой разности двух единичных функций б(x - xj )= U(x - aj )- U(x - bj );

s( y - y )= U (y - ci)-U (y - di) > где aj=xj - rj/2, bj = xj + rj/2, c>=y> - r/2>

di = yi+ri /2; r - ширина ребра. Следовательно, толщина всей конструкции равна h + H. Функционал полной потенциальной энергии деформации будет иметь вид [20]

ЕР = WMVl//{(( + Fx ) + G2 ( + Fy ) +^21 (h+F, + Fy)

\X12\X-2i) a1 0

e x e y +

2G12(2h + Fx + Fy) + Gnk(h + Fx)(, -01 )2 + G23*(h + Fy)(, -02) + 2X1

+ ^21 ( + Sy ) exX2 + ^21 (Sx + Sy ) eyX1 + 2G2SyeyX2 + 2G12 (Sx + Sy ) УxyX12

h

Г U3 Г

+

h

12

-+

X12 + (6)

+g2

12 + Jy

+ 2G

Г h3

12

X2 + M"21 Л

Г h3

V

+ Jx + Jy

X1X2 +

+ Jx + J у

6 x y

У

Х22 - 2 q (1 ^21) w 1 ABdxdy.

E1 J

Здесь Fx, Fy, Sx, Sy, Jx, Jy - площадь поперечного сечения ребра, приходящаяся на

единицу длины сечения; статический момент и момент инерции этого сечения, которые имеют вид [20]

hr

ir m Г hr ”hjr.r. Y

hr.

Гм

hp m hry

\

fx = У-^+I —1 -I—^ —; Fy = I—j+I -I-x 1=1 ь j=1 V a 1=1 ab j a y £ a 1=1 V b % ab

Sir

ir m Г Shy jl S^rr,л

SJrj " Г Siri m,SlJrirj

sx = Z-r+I -У--Z-ir -j; Sy = £-У-+I -J--I

i=1 b = a i=1 ab a j=1 a = b =

a i=1 ab

j=1 V -- i=1 —'

j=1

i=1 v b j=1 ab

(7)

78

Карпов В.В., Семенов А.А. /Вестник ПНИПУ. Механика 3 (2015) 74-94

n Jiy m

J]r, JL J4r,r

J, = Z-r+Z -Z-^ A• Jy =Z-^+Z -г-Z

,,1 b j .11 a i ,1 ab

J]r,

( П

a

Cl

J,ri m J]rirj

b ,=1 ab

где

hi (h + h) hj (h + hj) hj (h + hj) . S2 i, чэ

Si = —^ Sj = —Sij = —Ji = 0,25h2h + 0,5h(hi)+ 3(hij ,

Jj = 0,25h2hj + 0,5h (hj j + 1(hj )3, Jij = 0,25h2hj + 0,5h (h,j )2 + ).

Здесь переменные a, b позволяют ребрам жесткости иметь постоянную ширину

a - a

и определяются как a = aA, b = bB\

V 2 J

Выражения усилий и моментов для подкрепленных ортотропных оболочек имеют вид [20]

N =-

Ei

1 ^12^21

Ny = , El

1 ^12^21

Nxy = G12

Mx =

(h + Fy ) jy + 2SyX12 E1

1 ^12^21 E

My =1-----~

1 ^12^21

(h + Fx ) (( + Л21^y ) + )x ( + Д21X2 )

(h + Fy ) y + ^12S x ) + Sy (2 + ^12X1)

(h + Fx )xy + 2SxX12

(X1 +Д21Х2 )

(2+^12X1)

Nyx = G12

Sx (( + ^y) Sy (y + ^12Bx)

(h3 Л

----+ JX

v 12 x

(h3 Л

12 + y

J

Mxy = G12

(h3 ^ (h3 ^

(N + S' ?- 1 П T 12+Jy, V1^ J X12 , Myx = G12 (N + S' ?- 1 П T + Jx 12 x v1^ J X12

Ox = GL3* (h + Fx X'Fx-01), Qy = G23* (h + Fy )(ш y-02).

1.2. Безразмерные параметры для гладких ортотропных оболочек

Введем безразмерные параметры

U =

x 5 a - a y a1 = —, л = т, a b X = aA bB,

aUA V = bVB!, W = W ш x

h2 h2 h ,

k(=hkx , kЛ=hky ,

ш =________у_

^ h

h

(8)

(9)

79

Karpov V.V., Semenov A.A. / PNRPU Mechanics Bulletin 3 (2015) 74-94

p,12p,21 )a 2 A2 a = a y (1 i ^12^21) a A — —xya A

Eh2 , a y = E2 h2 , —y = G12h2

и l®4 ^1 II 4^ b = bA, z = z, h h

где a, b - линейные размеры оболочки в направлении осей x и y соответственно; X -безразмерный коэффициент; E, ц, z - новая (безразмерная) система координат оболочки.

Важную роль в формулировании безразмерных параметров играют параметры Ляме A, B поверхности оболочки. Значения параметров Ляме для цилиндрических, сферических, конических, тороидальных оболочек приведены в табл. 1.

Параметры Ляме для основных видов оболочек вращения

Таблица 1

Вид оболочки A B R1, м R2, м a b

Пологая, прямоугольная в плане 1 1 r1 r2 метр метр

Цилиндрическая 1 r (Ю r метр радиан

Сферическая r r • sin x r r радиан радиан

Коническая 1 x • sin 9 (Ю x • tan 9 метр радиан

Тороидальная r d + r • sin x r d + r • sin x sin x радиан радиан

Здесь R1, R2 - главные радиусы кривизны оболочки вдоль осей x, y ; d - смещение

образующего сектора от оси вращения для тороидальной оболочки; 9 - угол конусности конической оболочки.

Как видно из табл. 1, параметр b - безразмерный; если a - размерное, то A - безразмерное и наоборот, т.е. aA - размерная величина (м); B = B(x) - размерная величи-

на (м), т.е. bB - размерная величина. Следовательно, если ввести безразмерные параметры Ляме следующим образом: A = aA, B = bBB, то для всех рассматриваемых видов обо-

h h

лочек они будут безразмерными.

Расчет напряженно-деформированного состояния и устойчивости оболочечных конструкций в безразмерных параметрах позволяет получить критическую нагрузку для целой серии подобных оболочек (рис. 1, 2). При этом прогибы, усилия и моменты для различных подобных оболочек находятся масштабированием полученных безразмерных значений этих величин.

Перейдем к безразмерным параметрам (9) в соотношениях (1), (2), заготовив вначале выражения -A—SB, 9Х, 92 в безразмерных параметрах:

1 SB 1 h SB h SB

AB Sx AB ab SE, AaBb SE,;

80

Карпов В.В., Семенов А.А. /Вестник ПНИПУ. Механика 3 (2015) 74-94

Рис. 1. Подобные сферические оболочки

Рис. 2. Подобные конические оболочки

01 =-

с

02 =-

с h dW 1 , h2 -Л

+ —k%.—- U

at d% h q at h dW 1 h2 -Л

h C dW ,

------+k%U

aA "■

s„ =

V

1 dU

+ -k — V bB d^ h ц bB

- KW+102 =

h

J

2

=----X

aA

h JdW

A dx x 2 1 a2A2

dU

d%

+kV дц ц ,

- k% A2W +102

=at =at i5>

' h2 _

= 1dV 1 dB y B dy AB dx

h2 C, 2 dV 1 dB -

u - kW+102 =

2 л 2

a 2 A

1 - 27\2

X2 — +=—U - kц A 2W +-X2 0

dц B d% '' 2

1 dV 1 dU 1 dB

a2 A

h2 _

s y;

2 Л2 Sx;

2 ,2 “У’

a 2 A

^xy A dx B dy AB dx

h2 C X dV .dU . 1 dB- -0 x

V d% dц B d% 12

V + 0!02 =

h2 2

a2 A2

Уxz = kf (Z) Уyz = kf (z)

a2A2 Xjxy;

W x -01

W y-02

h_

aA

h

kf (z)

=oa Uf (]

—1 =

1 dW x 1 h dW x

A dx A Aa2

V x -01

Vy -0: h _ ;

a 2 2 —1;

A2a

h _

= aA У xz; h _

= at У yz;

—2 =

1 dW y

1 dB

B dy AB dx

W x =

h

A2 a 2

dW y 1 dB -

___y. j у

dц B d% x

h - .

a2 2 — 2;

A2 a

1 dW y 1 dW x

2 — =-!A^ + -XA±-_LdB W =^_ 12 A dx B dy AB dx y a2A2

X-

dW y

- + X

dW

1 dB

Л

-X=—W y dц B d% y

h

a2 A2

X2 —12.

81

Karpov V.V., Semenov A.A. / PNRPU Mechanics Bulletin 3 (2015) 74-94

При линейно-упругом деформировании физические соотношения для ортотропных оболочек вращения имеют вид (3), перейдем в них к безразмерным параметрам (9):

_ x =■

Ei

1 Ц-12^21

+ ^21s y + z ( +Д21Х2 )

E h2

1 -д12д21 a 2 A2

+ Д21Sy + z (X1 +Д21X2 )

E1 h2 _ ;

e2

_ y =

E2 h2

1 ^12^21

1 2 a2 x

_ 1 -^12^21 a A

e y x + z (2 +^12X1)

1 -д12д21 a 2 A2

- xy = G12 У xy + 2 zX12

+ ^12ex + z (X2 + ^12X1)

= G

h2

12 2 >2 a 2 A2

e2 h2 _ ;

_ 1 -д12д21 a2 A2 y’

2

У xy + 2 z X12

= G h2 -

12 a2A2 -xy ’

-xz = G13kf (z) -yz = G23kf (z)

* x -01

= G

h2

13 2 ,<2

a 2 A2

kf (z ) A

*y -02

= G

h2

23 a2A2

kf (z )AX

* x -01

* y -02

h2

G13 a2A2 -xz;

= G h2 -

23 a2A2 -yz'

Здесь

f (z )=6 f 4- z 2

Функционал полной потенциальной энергии деформации оболочки запишется следующим образом:

E =

^p п,„4 Л4

Eh

- 4 л\(л \lj|Sx +G2Sy +(21 + G2^12)SxSy + G12Xy.

2a4a4(1 -^12^21)) l y

+ -

(X12 + G2X2 + (21 + G2^12))2 + 4G12X2x122) + G13kA2(x - 01) + (10)

12

G23kA2X2(*y - 02)2 - 2(1 - |д12|д21)PW lABd^dri ■-

+

Eh7 -

E1h______e .

2a4A4(1 -д12д21 ) p

Таким образом, под знаком интеграла остаются только безразмерные слагаемые.

При необходимости можно получить безразмерные значения для усилий и моментов:

N =МА

1 М12М21

( x + ^ 21е y )= 1

( + ^ 21 Sy )

E1h

(1 M12M21))

tvt E2h i

Ny = 1------2------( y + ^12 e x ) =

1 -^21

)■ I

= Ejh h2 (_

^12^21 a A

( +^21ey )= E1hNx;

E2h h2 /_

2 -( + ^12ex ) =

^12^21 a A

2 a2\ у

) = e2 h;

82

Карпов В.В., Семенов А.А. /Вестник ПНИПУ. Механика 3 (2015) 74-94

h,

Nxy = Nyx = Gnhy xy = G12h-TAT ^xy =-=t ЬУху = G12hNxy

a A A

Mx = -

El

1 ^12^21

( h3 л V12 2

(+^21X2) E1h3

Eh h /_ _ ч

1 (X1 + Д21Х2 ) =

My = TT2—

1 M^12^21

( h3 ^

V12 2

12 (1 ^12^21 ))

(X2 +^12X1)

12 (1 -Д12Д21) a2A2 ((1+^21X2 ) = E1h Mx;

E2 h3

12 (1 -Д12Д21) a2A

(X2 + M"12(1 ) = E2h My;

Mxy = Myx = G12 6 X12 = G12 6A2 ^(12 = G12h Mxy;

Q, = Gukh( -01) = G^kt(x -01) = GnhQx;

Oy = G23M(y -02) = G,jtX((y -02) = G-ahQ!.

1.3. Безразмерные параметры для подкрепленных ортотропных оболочек

Для подкрепленных оболочек, помимо предложенных выше, необходимо ввести дополнительные безразмерные параметры

F = — F = S = — S = — J = — J = — h = — x h ’ у h’ x h2’ y h2’ x h3’ y h3’ h ’

f (z z +1)(z 4-H)=f (z)'

- xj -

xj = > rj =

a

r, h, a. - bt

j h = j a ■ = —, b1 =-j

a ’ j h ’ a a

r h, c, — d,

. , h l ci =~, dt =.

' b l ~ h ’ , b , b

■" b ’

Толщина всей конструкции будет равна 1 + H .

Функционал полной потенциальной энергии деформации будет иметь вид

Ер =

E1h7

2a A (1 ^12^21) 51

1 1

Я{(1 + Fx )Sx2 + G2 (1 + Fy )) +^21 (2 + Fx + Fy )SxSy L

(11)

+2 G12 (2 + Fx +Fy )2 +G13kA2 (1 + Fx )(-01 )2 +G23kA2X2 (1 + Fy)(-02 )2 +

+2Sx Sx(1 +^21 ( + Sy )SxX2 +Д21 (x + Sy )Sy(1 + 2G2Sy SyX2 + 2G12 (x + Sy )2Уxy X12 + (12)

+ 1^ 12 + Jx ^|X12 + G2 V12 + JL ^jX2 +^ 21 ^ 6 + Jx + —y ') X1X2 + 2G12 V 6 + —x + Jy | ^^ _

-2(1 - д12д21 l A5dE,dц =

Eh7 -

_______________E .

2a4 A4 (1 -д12д21) p

83

Karpov V.V., Semenov A.A. / PNRPU Mechanics Bulletin 3 (2015) 74-94

Выражения усилий и моментов для подкрепленных ортотропных оболочек запишутся следующим образом:

E

Nx = - E

1 Ц-12^21 Eft h2

(h + Fx ) (( + ^215y ) + ^x ( + Д21Х2 )

1 -д12д21 a2 AA E

Ny = TZ^-

1 M12M21

E2h h2

1 -д12д21 a 2 A2

(1 + Fx ) )( + ^21Sy ) + Sx ( + Д21X2 ) = E1hNx, (( + Fy ) («5y + ^12sx ) + Sy (X2 + ^12X1)

(1 + Fy ) ( + ^125x ) + )y ((2 + Д12Х1 )

= E2 hNy,

Nxy = G12 (h + Fy )xy + 2SyX12 G ш h2 = G12hX 2 a 2 A2 (1 + Fy )y xy + 2 Sy X12

Nyx = G12 (h + Fx )Уxy + 2SxX12 h2 = G12 h\— a 2 A2 (1 + Fx )yxy + 2SxX12

= G12hNxy,

= G12 hNyx ,

Mx = ■

E

Sx (Sx +^21sy ) +

E1h2 h2

1 -д12д21 a 2 A

My=v~E

1 ^12^21

= E2 h2 h2

1 ^12^21

Sx ( +^215y

f h3 1

----+ J

v 12 x

(X1 +Д21Х2 )

2 л2

) + ^ 12 + Jx l(X1 +^21X2)

Sy (Sy +^12Sx )

f h3 ТЛ 12 + y

= F1h2Mx,

(X2 +^12X1)

1 -д12д21 a2 A2

Sy (Sy +^12 sx ) 11! + Jy jfe +^12 X1)

= E2h2 My,

Mxy = G12

Myx = G12

1 f h3 1 1 ^ ~ h2

(N + S' ?- >3 1 12+J 1 <N = G12h Л 2 A 2 a 2 A2

(N + S' ?- 1 1 f h3 1 12+J‘, У X12 =G12 h2x a2A2 _

1

Sy Jxy + 21 12 + Jy I X12

Sx У xy + 21 12 + Jx | X12

= ^tfMxy,

= G12 h22yx,

Ox = G13A (A + F, )( x -01) = Guhk-A (1 + Fx)((x -01) = GuhQx,

Qy = G2,k(h + Fy)((y -02) = G23hk-hA(1 + Fr)((. -02) = G23AQ,..

aA

1.4. Безразмерные параметры для изотропных пологих оболочек при решении задач динамики и термоупругости

При решении уравнений в смешанной форме вводится функция напряжений в срединной поверхности оболочки Ф( x, y). В безразмерных параметрах она имеет вид

- . ч Ф(x, y)

Ф( .

84

Карпов В.В., Семенов А.А. /Вестник ПНИПУ. Механика 3 (2015) 74-94

При решении задач динамики вводится безразмерный параметр времени [22]

Здесь у, g - удельный вес материала оболочки и ускорение силы тяжести.

При решении задач термоупругости вводится безразмерный параметр температуры __ (ха2

[22] T =—— T, где а - коэффициент теплового расширения материала. h

1.5. Некоторые частные случаи

Рассмотрим частный случай для гладких пологих изотропных оболочек, прямоугольных в плане. В этом случае параметры Ляме будут A = 1, B = 1, а безразмерные параметры примут вид

b 2k„

е x y

t = - Л = т,

a b

1 =

a

U =

aU

If'

V =

bV_ h2 ''

b

W =

k =

a_kx

h

k" = h

о x =■

°xa

Eh2

2

_ °ya

о y = —

Eh2

xy

W - ш xa

—, ш x =

h h ’

2 4

x a — _ xy P a q

G12 h2, h4 E ’

- ¥ b

ш = y

y h ''

a = ■

h

z = ■

h

(13)

Следует отметить, что в данном варианте безразмерных параметров, помимо подстановки параметров Ляме и выполнения упрощений, связанных с изотропией материала, по-другому выбраны выражения для главных кривизн оболочки [22].

Выражения для деформаций примут следующий вид:

dU , ш 1 fdW' dx 21 dx

h2 f $U kW 1 -------ktW +-

a

fdW ^ ^ ^2

't

dV

s y =---kyW + -

y dy y 2

f dW ^

dy

=^2 12

a

h2 _

dV

- 1

---kW+-

d^ ц 2

f dW ^

2

d^

2 Sx;

a

=h_ >2.

a2 y;

У xy =

dV dU dW dW h

= ~1

dx dy dx dy a

dV dU dW dW

\ 1-2

Уxz = kf (z)

Уyz = kf (z)

ш x -

dW

dx

ш Ж

y dy _

v dt, d^ dt, дц j

dW

= h .

2 1— xy ;

a

= ~kf (z) a

= - 1k f (z) a

Ш x -x dt _

ш Ж

y d^

h_

=--xz;

a

h >У

=->-yz;

a

X1 =■

h d¥ x h _

______=_____y .

2 ж 2 ^1’

dx a2 dt a2

85

Karpov V.V., Semenov A.A. / PNRPU Mechanics Bulletin 3 (2015) 74-94

_8Vy _ h . 2 8Vy _ h X 2- .

-2 _ ~ _ 2 X ~ _ 2 X —2. 8y a 8n a

2 -12 _

8V v 8Vx h

_^X

8x Cy a2

8V y 8V y-+ ■

Л

8E, 8n

_ "~Г X2i2-

a

При линейно-упругом деформировании физические соотношения для изотропных оболочек вращения будут иметь вид (3) при E1 _ E2 _ E, ц12 _ц21 _ц,

G12 _ G13 _ G23 _

E

2 (1 + ц)

. Перейдем в них к безразмерным параметрам (13):

— _

Gy _■

E 8 x +Це y + z (1 +Ц-2 ) _ Eh2 ex +цх 2 ey + z (-1 +цх 2-2 )

1 -ц2 1 -ц2 a2

E 8y + ^x +z (-2 +ц-1 ) _ Eh2 X28y +^ex +z (X2-2 +ц-1)

1 -ц2 1 -ц2 a2

— -

_ Е— 2;

a

——

_ —-y.

a

E

X _

xy

2 (1+ ц)

E X h2

Y xy + 2 z-12 _ 2 (1 + ц)Х a2 Y xy + 2 z -12

, h2

_eT2 •

h_ h _

Xxz _ G13Y xz _ G13 _ Yxz ; X yz _ G23Y yz _ G23~ Xy yz .

a a

Функционал полной потенциальной энергии деформации пологой оболочки прямоугольного плана с учетом (13) будет иметь вид [22]

—p_ 3(4 Ш8 2+е y+ 2це 'е y+ G 2y+ Gk ( x -0> )2 +

' ( ц ) 0 0

х2 h2

+Gk(Vy -02)2 +12(-2 +—2 + 2Ц-,-2 + 4G—22)- 2W}dxdy _

Eh

5 1 1

Eh

— \ 2

,(1 - 2) 4J/]ex + X%2 + 2цХ28x8y + GX% + Gk(Vx-0,)2 + GkX2(-02)

2 (1 ц ) 0 0 ^

h- (-2 +X4-2 + 2цХ2-1-2 + 4GX2-22)- 2P (1 -ц2 )W J di;rfn _ [Tj

Ep.

„4 P

a

Здесь

G _ G12 (1 -ц )

_ E '

1.6. Дальнейший подход к решению задач прочности и устойчивости изотропных и ортотропных оболочечных конструкций в безразмерных и размерных параметрах

Для минимизации функционала полной потенциальной энергии деформации используем метод Ритца. Для решения задач в размерных параметрах представим искомые функции в виде

86

Карпов В.В., Семенов А.А. /Вестник ПНИПУ. Механика 3 (2015) 74-94

N N N

U(x,y) = jU(I)Z,(I); V(x,y) = j V(I)Z2(I); W(x,у) = jW(I)Z3(I);

I=1 I=1 I=1

NN

Vx (x, y) = j PS (I) Z4( I); V у (x, y) = j PN( I )Z,( I),

I=1 I=1

где U (I), V (I), W (I), PS (I), PN (I) - неизвестные числовые коэффициенты,

a Z1 (I) - Z5 (I) - известные аппроксимирующие функции аргументов x и y , удовлетворяющие заданным краевым условиям на контуре оболочки; N - количество членов разложения.

Для решения задач в безразмерных параметрах искомые функции

_ N _ _ _ N _ _ _ N _ _

U(5,v) = jU(I)Z,(I), V(5,vl) = jV(I)Z2(I); W(5,r,) = jW(I)Z,(I);

I=1 I=1 I=1

_ N ___ _ _ N ____ _

V x (5, v) = j PS (I) Z,(I); V y (5, v) = j PN (I) Z,( I),

I=1 I=1

здесь U (I), V (I), W (I), PS (I), PN (I) - неизвестные безразмерные числовые коэффициенты; Z1 (I) - Z5 (I) - известные аппроксимирующие функции аргументов 5 и v,

удовлетворяющие заданным краевым условиям на контуре оболочки.

Подставив разложения искомых функций в функционал (12) и проведя процедуру метода Ритца, получим систему нелинейных алгебраических уравнений. Полученную систему будем решать с помощью метода продолжения решения по наилучшему параметру [23, 24]. Точность применения метода продолжения решения по наилучшему параметру можно свести к анализу угла между векторами приращений неизвестных параметров с прошлого и текущего шага решения (сумма квадратов которых на каждом шаге дает квадрат длины дуги приращения в многомерном пространстве). Так, если угол между векторами решений на предыдущем и текущем шаге оказался слишком большим, то производится уменьшение параметра длины дуги в 10 раз и откат итерационного процесса на два шага назад (для рассматриваемых в данной работе задач такое ограничение угла было выбрано 0,0177 рад). Если значение угла находится в пределах от 0,0094 до 0,0177 рад, то длина дуги сохраняется такой же, как на предыдущем шаге. Если значение угла не превысило 0,0094 рад, то длина дуги увеличивается в 1,5 раза. Таким образом, можно достичь достаточной точности решения вблизи особых точек и в то же время быстро проходить ровные участки кривой.

Точность метода Ритца зависит от числа разложений искомых функций в ряды. С помощью вычислительного эксперимента было показано, что для рассматриваемых оболочек достаточно принять N = 16 для обеспечения достаточно высокой точности [25].

Такой подход позволяет исследовать прочность и устойчивость оболочек, обходить особые точки кривой «нагрузка - прогиб», получать значения верхней и нижней критических нагрузок, находить точки бифуркации и исследовать закритическое поведение конструкции [23, 24].

87

Karpov V.V., Semenov A.A. / PNRPU Mechanics Bulletin 3 (2015) 74-94

2. Расчеты

2.1. Характеристика рассматриваемых конструкций

Все рассматриваемые далее конструкции имеют шарнирно-неподвижное закрепление контура и находятся под действием внешней статической равномерно-распределенной поперечной нагрузки, направленной по нормали к срединной поверхности.

В дальнейшем будут рассматриваться конструкции, выполненные из следующих изотропных и ортотропных материалов (табл. 2).

Таблица 2

Механические характеристики материалов рассматриваемых конструкций

Характе- ристика Ортотропные Изотропные

Углепластик Стеклопластик Сталь

M60J/Epoxy T-10/УПЭ 22-27

E1, МПа 3,3105 0,294-105 2,1-105

Ml2 0,32 0,123 0,3

E2, МПа 0,59104 1,78-104 2,1-105

G12, МПа 0,39104 0,301-104 0,807-105

G13, МПа 0,39104 0,301-104 0,807-105

G23, МПа 0,39104 0,301-104 0,807-105

Fx+, МПа 1760 508 -

F", МПа -780 -209 -

F2+, МПа 30 246 -

F2 , МПа -168 -117 -

Fu, МПа 39 43 -

ст, МПа - - 1720

Плотн., кг/м3 1500 1800 7800

Используя безразмерные параметры, будем одним расчетом проводить исследование прочности и устойчивости целой серии подобных оболочек. В методе Ритца примем

N = 16.

2.2. Пологие оболочки, квадратные в плане

Рассмотрим вначале подобные квадратные в плане пологие оболочки с размерными и безразмерными параметрами (табл. 3).

Таблица 3

Пример размерных и безразмерных параметров для подобных квадратных

в плане пологих оболочек

Вариант a, м R1, м h, м a R1

1 18 45,3 0,03 600 1510 238,41

2 27 67,95 0,045

3 36 90,6 0,06

4 54 135,9 0,09

88

Карпов В.В., Семенов А.А. /Вестник ПНИПУ. Механика 3 (2015) 74-94

Проведенные расчеты в размерных параметрах для оболочки варианта 4 из углепластика M60J показали, что критическая нагрузка первой потери устойчивости составила Pkr = 5198,543 (размерное значение qkr = 0,01323 МПа), нагрузка потери прочности

Ppr = 12670,4 (размерное значение qpr = 0,0322 МПа). Эти нагрузки будут одинаковыми

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

Максимальный прогиб при потере устойчивости Wmax был равен 0,858, при переводе

в размерные параметры это будет 0,077 м. Значения, полученные переходом к размерным параметрам для остальных вариантов конструкций, показаны в табл. 4.

Таблица 4

Размерные и безразмерные результаты расчета для подобных квадратных

в плане пологих оболочек

Параметр (безр.) Значение Параметр (разм.) При переводе в размерные варианты Сравнение (расчет в разм.)

Вар. 1 Вар. 2 Вар. 3 Вар. 4 Значение Разница, %

Входные параметры конструкции

а 600 а, м 18 27 36 54 54 -

1510 R1, м 45,3 67,95 90,6 135,9 135,9 -

- - h, м 0,03 0,045 0,06 0,09 0,09 -

Результаты расчета

Pkr 5198,543 qkr, МПа 0,01323 0,01318 0,4

Ppr 12670 qpr, МПа 0,0322 0,0308 4,3

Wmax 0,858 Wnax, м 0,026 0,039 0,051 0,077 0,088 14,6

2.3. Панели цилиндрических оболочек

Рассмотрим далее подобные панели цилиндрических оболочек, выполненных из стеклопластика. Расчет будем проводить для конструкции с параметрами

а = 2000, R = 540, угол разворота к.

Критическая нагрузка первой потери устойчивости Pkr = 3,554132 -107 (размерное значение qkr = 0,065307 МПа), нагрузка потери прочности Ppr = 2,194-107 (размерное значение qpr = 0,0403 МПа). Эти нагрузки будут одинаковыми для всех рассмотренных вариантов оболочек. Максимальный прогиб при потере устойчивости Wmax был равен

20,62, при переводе в размерные параметры это будет 0,2062 м.

Были также проведены расчеты этой же оболочки в размерных параметрах (а = 20 м, R = 5,4 м, h = 0,01 м). Значения, полученные при расчете и в безразмерных, и в размерных параметрах показаны в табл. 5.

89

Karpov V.V., Semenov A.A. / PNRPU Mechanics Bulletin 3 (2015) 74-94

Размерные и безразмерные результаты расчета для подобных цилиндрических панелей

Таблица 5

Пара- метр (безр.) Безразм. значение Параметр (разм.) При переводе в размерные варианты Сравнение (расчет в разм.)

Вар. 1 Вар. 2 Вар. 3 Вар. 4 Значение Разница, %

Входные параметры конструкции

a 2000 а, м 5 10 15 20 20 -

- - h, м 0,0025 0,005 0,0075 0,01 0,01 -

Результаты расчета

Pkr 3,554 10' qkr, МПа 0,065307 0,066250 1

Ppr 2,194107 qpr, МПа 0,0403 0,0407 1

Wax 20,62 Wmax, м 0,0515 0,1031 0,1546 0,2062 0,2192 6

2.4. Панели усеченных конических оболочек

Рассмотрим панели усеченных гладких конических оболочек, выполненных из изотропного материала (стали). Размерные параметры и соответствующие им безразмерные параметры приводятся в табл. 6. Чтобы было возможным рассматривать оболочки как подобные, был взят одинаковый угол конусности 9 = 0,78 рад и одинаковый угол разворота

b = п рад. Также соблюдается правило равенства отношений —, — для всех рассматри-

h h

ваемых вариантов. Результаты расчетов прочности и устойчивости для этих конструкций приводятся в табл. 7. Здесь потеря прочности происходит одновременно с потерей устойчивости, поэтому значения соответствующих нагрузок объединены. Сравнение с расчетом в размерных параметрах проводилось по конструкции варианта 2.

Таблица 6

Пример размерных и безразмерных параметров для подобных панелей усеченных конических оболочек

Вариант а, м а1, м h, м а — а1, м а = а / h а1 = а1 / h

1 12,5 2,5 0,005 10 2500 500

2 25 5 0,010 20

3 37,5 7,5 0,015 30

4 50 10 0,020 40

Таблица 7

Размерные и безразмерные результаты расчета для подобных панелей усеченных конических оболочек

Безразм. параметр Безразм. значение Размерный параметр При переводе в размерные варианты Сравнение (расчет в разм.)

Вар. 1 Вар. 2 Вар. 3 Вар. 4 Значение Разница, %

Входные параметры конструкции

а 2500 а, м 12,5 25 37,5 50 25 -

- - h, м 0,005 0,010 0,015 0,020 0,010 -

Результаты расчета

P = P kr pr 5.07 -107 qkr qpr, МПа 0,2730 0,2847 4

W п max 1 34.27 Wmax 1, м 0,1714 0,3427 0,5141 0,6854 0,3515 2

90

Карпов В.В., Семенов А.А. /Вестник ПНИПУ. Механика 3 (2015) 74-94

3. Анализ

Для оболочек, изготовленных из одного и того же материала, с одинаковыми значениями безразмерных параметров A, B, , кц и разными размерными значениями

a, b, h, R1, R2 (подобные оболочки) значения критической нагрузки потери устойчивости (или прочности) P и напряжений ax, ay, Т будут одинаковыми. Соответственно, одинаковыми будут и размерные значения нагрузки и напряжений q, ax, ay, тxy .

Чтобы найти размерные значения прогиба, усилий и моментов, необходимо выполнить следующие действия:

W = hW, Nx = ElhNx, Ny = E2 hNy, Nxy = G12 hNy, Nyx = G12 hNyx,

Mx = Eh 2MX, My = E2 h 2My, Mxy = G12 h 2Mxy, Myx = G12 h 2Myx,

Qx = Gi3hQx, Qy = G2hQy.

Заключение

Введение безразмерных параметров при расчете оболочечных конструкций позволяет получить более обширную информацию о напряженно-деформированном состоянии оболочек и выявить особенности деформирования для целой серии подобных оболочек.

Как видно из проведенных исследований, решения в размерных и безразмерных параметрах несколько отличаются друг от друга. Значения критических нагрузок близки, но кривые «нагрузка - прогиб» не совпадают.

При решении в размерных параметрах искомые функции имеют разный порядок. Так, прогибы W имеют порядок 10-1, а U,V,¥x,¥y - порядок 10_2-10-3. Кроме того,

при действии с малыми числами накапливается инструментальная ошибка. При N = 16 порядок системы алгебраических уравнений равен 81.

При решении задачи в безразмерных параметрах порядок искомых функций выравнивается и инструментальная ошибка накапливается меньше. Однако порядок параметра нагрузки P может существенно отличаться от порядка остальных параметров и для каждой конкретной задачи быть разным. Правильный подбор параметров точности расчета в методе продолжения решения по наилучшему параметру позволяет достаточно точно находить значения критических нагрузок, несмотря на разницу в порядках. Так что расчеты в безразмерных параметрах следует признать более справедливыми.

Работа выполнена в рамках государственного задания Минобрнауки России, проект № 3801.

Библиографический список

1. Амельченко В.В., Неверов И.В., Петров В.В. Решение нелинейных задач теории пологих оболочек путем вариационных итераций / Изв. АН СССР. Механика твердого тела. - 1969. - № 3. -С. 62-68.

2. Петров В.В. Метод последовательных нагружений в нелинейной теории пластинок и оболочек. - Саратов: Изд-во Сарат. гос. ун-та, 1975. - 119 с.

91

Karpov V.V., Semenov A.A. / PNRPU Mechanics Bulletin 3 (2015) 74-94

3. Петров В.В., Филатов В.Н. Расчет гибких пластинок вариационным методом В.З. Власова // Изв. вузов. Строительство и архитектура. - 1970. - № 2.

4. Петров В.В. Расчет гибких пластин и пологих оболочек вариационным методом В.З. Власова // Прикладная механика. - 1966. - Т. II. - Вып. 5.

5. Карпов В.В., Филатов В.Н. Закритические деформации гибких пластин в температурном поле с учетом изменения свойств материала от нагревания // Материалы VII Всесоюз. конф. по теории оболочек и пластин (Днепропетровск, 1969). - М.: Наука, 1970. - С. 276-280.

6. Карпов В.В., Петров В.В. Уточнение решений при использовании шаговых методов в теории гибких пластинок и оболочек // Изв. АН СССР. МТТ. - 1975. - № 5. - С. 189-191.

7. Крысько В.А., Федоров П.Б. Исследование динамической устойчивости гибкой пологой оболочки в зависимости от механических и тепловых характеристик // Прикладная механика. -1984. - Т. XX, № 3. - С. 45-49.

8. Grigorenko Ya.M., Kryukov N.N., Ivanova Yu.I. Solution of two-dimensional problems of the statics of flexible shallow shells by spline approximation // International Applied Mechanics. - 1995. -Vol. 31. - No. 4. - P. 255-260.

9. Валишвили Н.В. Методы расчета оболочек вращения на ЭЦВМ. - М.: Машиностроение, 1976. - 278 с.

10. Elastic Buckling Analysis of Ring and Stringer-stiffened Cylindrical Shells under General Pressure and Axial Compression via the Ritz Method / Arani A. Ghorbanpour, A. Loghman, A.A. Mosallaie Barzoki, R. Kolahchi // Journal of Solid Mechanics. - 2010. - Vol. 2. - No. 4. - P. 332-347.

11. Вольмир А.С. Нелинейная динамика пластинок и оболочек. - М.: Наука, 1972. - 432 с.

12. Вольмир А.С. Устойчивость деформируемых систем. - М.: Наука, 1967. - 984 с.

13. Zenon del Prado, Gonsalves Paulo B., Pai'doussis Michael P. Non-linear vibrations and instabilities of orthotropic cylindrical shells with internal flowing fluid // International Journal of Mechanical Sciences. - 2010. - Vol. 52. - P. 1437-1457. DOI: 10.1016/j.ijmecsci.2010.03.016

14. Янковский А.П. Асимптотический анализ решения нелинейной задачи нестационарной теплопроводности слоистых анизотропных неоднородных оболочек при граничных условиях первого рода на лицевых поверхностях // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. - 2014. -№ 1 (34). - С. 168-185. DOI: 10.14498/vsgtu1281

15. Yasin M. Yaqoob, Kapuria S. An efficient layerwise finite element for shallow composite and

sandwich shells // Composite Structures. - 2013. - Vol. 98. - P. 202-214. DOI:

10.1016/j.compstruct.2012.10.048

16. Large deflection analysis of laminated composite plates resting on nonlinear elastic foundations by the method of discrete singular convolution / A.K. Baltacioglu, O. Civalek, B. Akgoz, F. Demir // International Journal of Pressure Vessels and Piping. - 2011. - Vol. 88. - Iss. 8-9. - P. 290-300. DOI: 10.1016/j.ijpvp.2011.06.004

17. Wu C.-P., Wang Y.-M., Hung Y.-C. Asymptotic finite strip analysis of doubly curved laminated shells // Computational Mechanics. - 2001. - Vol. 27. - P. 107-118. DOI: 10.1007/s004660000218

18. Новожилов В.В. Теория тонких оболочек. - Л.: Судпромиздат, 1962. - 431 с.

19. Амбарцумян С.А. Теория анизотропных оболочек. - М.: Физматлит, 1961. - 384 с.

20. Карпов В.В., Семенов А.А. Математическая модель деформирования подкрепленных ор-тотропных оболочек вращения // Инженерно-строительный журнал. - 2013. - № 5. - С. 100-106. DOI:10.5862/MCE.40.11

21. Карпов В.В. Прочность и устойчивость подкрепленных оболочек вращения: в 2 ч. Ч. 1: Модели и алгоритмы исследования прочности и устойчивости подкрепленных оболочек вращения. - М.: Физматлит, 2010. - 288 с.

22. Карпов В.В., Игнатьев О.В., Сальников А.Ю. Нелинейные математические модели деформирования оболочек переменной толщины и алгоритмы их исследования. - М.: АСВ; СПб.: Изд-во С.-Петерб. гос. архит.-строит. ун-та, 2002. - 420 с.

92

Карпов В.В., Семенов А.А. /Вестник ПНИПУ. Механика 3 (2015) 74-94

23. Семенов А.А. Алгоритмы исследования прочности и устойчивости подкрепленных орто-тропных оболочек // Строительная механика инженерных конструкций и сооружений. - 2014. -№ 1. - С. 49-63.

24. Шалашилин В.И., Кузнецов Е.Б. Метод продолжения решения по параметру и наилучшая параметризация (в прикладной математике и механике). - М.: Эдиториал УРСС, 1999. - 224 с.

25. Карпов В.В. Прочность и устойчивость подкрепленных оболочек вращения: в 2 ч. Ч. 2: Вычислительный эксперимент при статическом механическом воздействии. - М: Физматлит, 2011. - 248 с.

References

1. Amel’chenko V.V., Neverov I.V., Petrov V.V. Reshenie nelineinykh zadach teorii pologikh obolochek putem variatsionnykh iteratsii [Solution of nonlinear problems in the theory of shallow shells by variational iteration]. Izvestiia Akademii naukSSSR. Mekhanika tverdogo tela, 1969, no. 3.

2. Petrov V.V. Metod posledovatel’nykh nagruzhenii v nelineinoi teorii plastinok i obolochek [The method of successive loadings in the nonlinear theory of plates and shells]. Saratovskii gosudarstvennyi universitet, 1975. 119 p.

3. Petrov V.V., Filatov V.N. Raschet gibkikh plastinok variatsionnym metodom V.Z. Vlasova [Calculation of flexible plates by the variational method of V.Z. Vlasov]. Izvestiia vuzov. Stroitel'stvo i arkhitektura, 1970, no. 2.

4. Petrov V.V. Raschet gibkikh plastin i pologikh obolochek variatsionnym metodom V.Z. Vlasova [Calculation of flexible plates and shallow shells by the variational method of V.Z. Vlasov]. Prikladnaia mekhanika, 1966, vol. II, iss. 5.

5. Karpov V.V., Filatov V.N. Zakriticheskie deformatsii gibkikh plastin v temperaturnom pole s uchetom izmeneniia svoistv materiala ot nagrevaniia [Supercritical deformation of flexible plates in the temperature field, taking into account the changes in material properties from heating]. Materialy VII Vsesoiuznoi konferentsii po teorii obolochek i plastin (Dnepropetrovsk, 1969). Moscow: Nauka, 1970. P.276-280.

6. Karpov V.V., Petrov V.V. Utochnenie reshenii pri ispol'zovanii shagovykh metodov v teorii gibkikh plastinok i obolochek [Clarification solutions using stepper methods in the theory of flexible plates and shells]. Izvestiia Akademii nauk SSSR. MTT, 1975, no. 5, pp. 189-191.

7. Krys’ko V.A., Fedorov P.B. Issledovanie dinamicheskoi ustoichivosti gibkoi pologoi obolochki v zavisimosti ot mekhanicheskikh i teplovykh kharakteristik [Investigation of the dynamic stability of a shallow shell flexible depending on the mechanical and thermal characteristics]. Prikladnaia mekhanika, 1984, vol. XX, no. 3, pp. 45-49.

8. Grigorenko Ya.M., Kryukov N.N., Ivanova Yu.I. Solution of two-dimensional problems of the statics of flexible shallow shells by spline approximation. International Applied Mechanics, 1995, vol. 31, no. 4, pp. 255-260.

9. Valishvili N.V. Metody rascheta obolochek vrashcheniia na ETsVM [Methods for calculating the shells of revolution by computer]. Moscow: Mashinostroenie, 1976. 278 p.

10. Ghorbanpour Arani A., Loghman A., Mosallaie Barzoki A.A., Kolahchi R. Elastic Buckling Analysis of Ring and Stringer-stiffened Cylindrical Shells under General Pressure and Axial Compression via the Ritz Method. Journal of Solid Mechanics, 2010, vol. 2, no. 4, pp. 332-347.

11. Vol’mir A.S. Nelineinaia dinamika plastinok i obolochek [Nonlinear dynamics of plates and shells]. Moscow: Nauka, 1972. 432 p.

12. Vol’mir A.S. Ustoichivost’ deformiruemykh system [Stability of deformable systems]. Moscow: Nauka, 1967. 984 p.

13. Zenon del Prado, Gonsalves Paulo B., Paidoussis Michael P. Non-linear vibrations and instabilities of orthotropic cylindrical shells with internal flowing fluid. International Journal of Mechanical Sciences, 2010, vol. 52, pp. 1437-1457. DOI: 10.1016/j.ijmecsci.2010.03.016

93

Karpov V.V., Semenov A.A. / PNRPU Mechanics Bulletin 3 (2015) 74-94

14. Iankovskii A.P. Asimptoticheskii analiz resheniia nelineinoi zadachi nestatsionamoi teploprovodnosti sloistykh anizotropnykh neodnorodnykh obolochek pri granichnykh usloviiakh pervogo roda na litsevykh poverkhnostiakh [Asymptotic analysis of solutions of a nonlinear problem of unsteady heat conduction of layered anisotropic inhomogeneous shells under boundary conditions of the first kind on the front surfaces]. Vestnik Samarskogo gosudarstvennogo tekhnicheskogo universiteta. Fiziko-matematicheskie nauki, 2014, no. 1 (34), pp. 168-185. DOI: 10.14498/vsgtu1281

15. Yasin M. Yaqoob, Kapuria S. An efficient layerwise finite element for shallow composite and sandwich shells. Composite Structures, 2013, vol. 98, pp. 202-214. DOI: 10.1016/j.compstruct.2012.10.048

16. Baltacioglu A.K., Civalek O., Akgoz B., Demir F. Large deflection analysis of laminated composite plates resting on nonlinear elastic foundations by the method of discrete singular convolution. International Journal of Pressure Vessels and Piping, (2011), vol. 88, iss. 8-9, pp. 290-300. DOI: 10.1016/j.ijpvp.2011.06.004

17. Wu C.-P., Wang Y.-M., Hung Y.-C. Asymptotic finite strip analysis of doubly curved laminated shells. Computational Mechanics, 2001, vol. 27, pp. 107-118. DOI: 10.1007/s004660000218

18. Novozhilov V.V. Teoriia tonkikh obolochek [The theory of thin shells]. Leningrad: Sudpromizdat, 1962. 431 p.

19. Ambartsumian S.A. Teoriia anizotropnykh obolochek [The theory of anisotropic shells]. Moscow: Fizmatlit, 1961. 384 p.

20. Karpov V.V., Semenov A.A. Matematicheskaia model' deformirovaniia podkreplennykh ortotropnykh obolochek vrashcheniia [Mathematical model of deformation of orthotropic reinforced shells of revolution]. Inzhenerno-stroitel'nyi zhurnal, 2013, no. 5, pp. 100-106. DOI:10.5862/MCE.40.11

21. Karpov V.V. Prochnost’ i ustoichivost’ podkreplennykh obolochek vrashcheniia. Ch. 1: Modeli i algoritmy issledovaniia prochnosti i ustoichivosti podkreplennykh obolochek vrashcheniia [The strength and stability of reinforced shells of revolution. Part 1: Models and algorithms of research strength and stability of reinforced shells of revolution]. Moscow: Fizmatlit, 2010. 288 p.

22. Karpov V.V., Ignat’ev O.V., Sal’nikov A.Iu. Nelineinye matematicheskie modeli deformirovaniia obolochek peremennoi tolshchiny i algoritmy ikh issledovaniia [Nonlinear mathematical models of deformation of shells of variable thickness and algorithms for their research]. Moscow: ASV; St. Peterburgskii gosudarstvennyi arkhitekturno-stroitelnyi universitet, 2002. 420 p.

23. Semenov A.A. Algoritmy issledovaniia prochnosti i ustoichivosti podkreplennykh ortotropnykh obolochek [Algorithms for the research of strength and stability of reinforced orthotropic shells]. Stroitel’naia mekhanika inzhenernykh konstruktsii i sooruzhenii, 2014, no. 1, pp. 49-63.

24. Shalashilin V.I., Kuznetsov E.B. Metod prodolzheniia resheniia po parametru i nailuchshaia parametrizatsiia (v prikladnoi matematike i mekhanike) [The method of the parameter continuation and best parameterization (in applied mathematics and mechanics)]. Moscow: Editorial URSS, 1999. 224 p.

25. Karpov V.V. Prochnost’ i ustoichivost’ podkreplennykh obolochek vrashcheniia. Ch. 2: Vychislitel’nyi eksperiment pri staticheskom mekhanicheskom vozdeistvii [The strength and stability of reinforced shells of rotation. Part 2: Computer experiment in static mechanical action]. Moscow: Fizmatlit, 2011. 248 p.

94