МЕХАНИКА

УДК 531.8

Р. М. Киракосян

Об уточненных уравнениях устойчивости анизотропных пластин

(Представлено академиком С. А. Амбарцумяном 22/VII 2004)

   1. Система уравнений устойчивости трехмерной теории упругости перед последним упрощением имеет вид ([1], с. 155, система v.13):
   
x
й
к
л
sўxx + eўxxsxx0 + ж
з
и
1
2
eўxy - wўz ц
ч
ш
sxy0 + ж
з
и
1
2
eўxz + wўy ц
ч
ш
sxz0 щ
ъ
ы
+
+
y
й
к
л
sўxy + eўxxsxy0 + ж
з
и
1
2
eўxy - wўz ц
ч
ш
syy0 + ж
з
и
1
2
eўxz +y ц
ч
ш
syz0 щ
ъ
ы
+
+
z
й
к
л
sўxz + eўxxsxz0 + ж
з
и
1
2
eўxy - wўz ц
ч
ш
syz0 + ж
з
и
1
2
eўxz +y ц
ч
ш
szz0 щ
ъ
ы
= 0, Cycl.
(1.1)

   Здесь sij0 - напряжения исходного состояния равновесия, sijў, eijў и wijў - дополнительные напряжения, деформации и углы поворотов, возникающие вследствие потери устойчивости.
   Для eijў и wijў имеем [1]:
eўxx = uў
x
, eўyy = vў
y
, eўzz = w
z
, eўxy = uў
y
+ vў
x
, eўyz = vў
z
+ w
y
, eўxz = uў
z
+ w
x
,
wўx = 1
2
ж
з
и
w
y
- vў
z
ц
ч
ш
,  wўy = 1
2
ж
з
и
uў
z
- w
x
ц
ч
ш
,  wўz = 1
2
ж
з
и
vў
x
- uў
y
ц
ч
ш
,
(1.2)

где uў, vў, w - дополнительные перемещения по осям x,y,z соответственно.
   Пренебрегая членами, умноженными на eўij, по сравнению с членами, умноженными на wjў, система (1.1) приводится к окончательному виду ([1], с. 157, система (v. 16)):

x
(sўxx - wўzsxy0 +ysxz0) +
y
(sўxy - wўzsyy0 +ysyz0) +
+
z
(sўxz - wўzsyz0 +yszz0) = 0,  Cycl.
(1.3)

   В  результате  этого  упрощения  члены  типа   [1/2]eўij ± wўk,  точные  значения  которых  равны   [(uўi)/(xj)],  заменяются  их  приближенными  значениями [1/2]([(uўi)/(xj)]-[(uўj)/(xi)]). Например, члены
1
2
eўzx - wўy = 1
2
ж
з
и
uў
z
+ w
x
ц
ч
ш
- 1
2
ж
з
и
uў
z
- w
x
ц
ч
ш
= w
x
,
1
2
eўzy +x = 1
2
ж
з
и
w
y
+ vў
z
ц
ч
ш
+ 1
2
ж
з
и
w
y
- vў
z
ц
ч
ш
= w
y
(1.4)
заменяются значениями
-wўy = 1
2
ж
з
и
w
x
- uў
z
ц
ч
ш
  и    wўx = 1
2
ж
з
и
w
y
- vў
z
ц
ч
ш

соответственно. В силу этого при составлении уравнений устойчивости пластин на основе системы (1.3) в выражении фиктивной поперечной нагрузки Z2 появляются лишние члены, имеющие характер поперечных сдвигов [2]. Поэтому в теории пластин более предпочтительной следует считать систему (1.1), применение которой не связано с особо серьезными трудностями и не приводит к появлению отмеченных членов.
   2. Рассмотрим прямоугольную пластину постоянной толщины h0 и размеров в плане a0,b0. Материал пластины в каждой точке имеет лишь одну плоскость упругой симметрии, параллельную срединной плоскости. Пластину отнесем к правой системе декартовых координат Oxyz, совместив плоскость Oxy со срединной плоскостью пластины. Пусть в пластине создано однородное напряженное состояние
sxx0,  syy0,  sxy0,       sxz0 = syz0 = szz0 = 0.
(2.1)

   Рассмотрим задачу устойчивости пластины с учетом влияний поперечных сдвигов, нормального напряжения zz и изменения первоначальных размеров h0, a0, b0. При составлении уравнений устойчивости будем исходить из сравнительно более точной системы (1.1).
   Пользуясь законом Гука анизотропного тела и геометрически линейными соотношениями [3], для толщины и размеров в плане в момент наступления потери устойчивости с учетом (2.1) будем иметь:
h = h0(1 + ezz0) = h0(1 + a13sxx0 + a23syy0 + a36sxy0),
a = a0(1 + exx0) = a0(1 + a11sxx0 + a12syy0 + a16sxy0),
b = b0(1 + eyy0) = b0(1 + a12sxx0 + a22syy0 + a26sxy0).
(2.2)

Здесь aij - упругие постоянные материала.
   Влияния поперечных сдвигов будем учитывать на основе уточненной теории пластин [3]. В рамках этой теории
sўxz = 1
2
ж
з
и
h2
4
- z2 ц
ч
ш
j(x,y),       sўyz = 1
2
ж
з
и
h2
4
- z2 ц
ч
ш
y(x,y),
(2.3)

где j и y - функции, характеризующие распределение поперечных сдвигов, возникающих вследствие потери устойчивости пластины.
   С учетом (2.3), закона Гука и геометрических соотношений имеем:
eўxx = -z 2w
x2
+ z
2
ж
з
и
h2
4
- z2
3
ц
ч
ш
ж
з
и
a55 j
x
+ a45 y
x
ц
ч
ш
,    1
2
eўxz - wўy = w
x
,
eўyy = -z 2w
y2
+ z
2
ж
з
и
h2
4
- z2
3
ц
ч
ш
ж
з
и
a44 y
y
+ a45 j
y
ц
ч
ш
,    1
2
eўyz +x = w
y
,
1
2
eўxy - wўz = uў
y
= -z 2w
xy
+ z
2
ж
з
и
h2
4
- z2
3
ц
ч
ш
ж
з
и
a55 j
y
+ a45 y
y
ц
ч
ш
,
1
2
eўxy + wўz = vў
x
= -z 2w
xy
+ z
2
ж
з
и
h2
4
- z2
3
ц
ч
ш
ж
з
и
a44 y
y
+ a45 j
y
ц
ч
ш
.
(2.4)

   Умножив первые два уравнения системы (1.1) на zdz, а третье уравнение - на dz и проинтегрировав полученные уравнения по толщине пластины, с учетом (1.1), (2.3), (2.4) и нулевых поверхностных условий находим:
Mx
x
+ Mxy
y
= Nx + h3
12
м
н
о
sxx0 й
к
л
3w
x3
- h2
10
ж
з
и
a55 2j
x2
+ a45 2y
x2
ц
ч
ш
щ
ъ
ы
+ syy0 й
к
л
3w
xy2
-
- h2
10
ж
з
и
a55 2j
y2
+ a45 2y
y2
ц
ч
ш
щ
ъ
ы
+ 2sxy0 й
к
л
3w
x2y
- h2
10
ж
з
и
a55 2j
xy
+ a45 2y
xy
ц
ч
ш
щ
ъ
ы
ь
э
ю
,
My
y
+ Mxy
x
= Ny + h3
12
м
н
о
syy0 й
к
л
3w
y3
- h2
10
ж
з
и
a44 2y
y2
+ a45 2j
y2
ц
ч
ш
щ
ъ
ы
+ sxx0 й
к
л
3w
x2y
-
- h2
10
ж
з
и
a44 2y
x2
+ a45 2j
x2
ц
ч
ш
щ
ъ
ы
+ 2sxy0 й
к
л
3w
xy2
- h2
10
ж
з
и
a44 2y
xy
+ a45 2j
xy
ц
ч
ш
щ
ъ
ы
ь
э
ю
,
Nx
x
+ Ny
y
= -h ж
з
и
sxx0 2 w
x2
+ syy0 2 w
y2
+ 2sxy0 2 w
xy
ц
ч
ш
.
(2.5)

   Через Mx, My, Mxy и Nx, Ny обозначены моменты и поперечные силы пластины соответственно.
   Сравним эти уравнения с уравнениями изгиба пластины [3]

Mx
x
+ Mxy
y
= Nx - hX1,    My
y
+ Mxy
x
= Ny - hY1,    Nx
x
+ Ny
y
= -Z2,
(2.6)
где
X1 = (X+ - X-)/2,    Y1 = (Y+ - Y-)/2,    Z2 = Z+ + Z-,
(2.7)

X±, Y±, Z± - интенсивности поверхностной нагрузки по осям x,y,z на поверхностях z = ±h/2 соответственно. Это сравнение приводит к равенствам:

X1 = - h2
12
м
н
о
sxx0 й
к
л
3w
x3
- h2
10
ж
з
и
a55 2j
x2
+ a45 2y
x2
ц
ч
ш
щ
ъ
ы
+ syy0 й
к
л
3w
xy2
-
- h2
10
ж
з
и
a55 2j
y2
+ a45 2y
y2
ц
ч
ш
щ
ъ
ы
+ 2sxy0 й
к
л
3w
x2y
- h2
10
ж
з
и
a55 2j
xy
+ a45 2y
xy
ц
ч
ш
щ
ъ
ы
ь
э
ю
,
(2.8)
Y1 = - h3
12
м
н
о
syy0 й
к
л
3w
y3
- h2
10
ж
з
и
a44 2y
y2
+ a45 2j
y2
ц
ч
ш
щ
ъ
ы
+ sxx0 й
к
л
3w
x2y
-
- h2
10
ж
з
и
a44 2y
x2
+ a45 2j
x2
ц
ч
ш
щ
ъ
ы
+ 2sxy0 й
к
л
3w
xy2
- h2
10
ж
з
и
a44 2y
xy
+ a45 2j
xy
ц
ч
ш
щ
ъ
ы
ь
э
ю
,
(2.9)
Z2 = Tx0 2 w
x2
+ Ty0 2 w
y2
+ 2S0 2 w
xy
.
(2.10)
   Здесь Tx0, Ty0, S0 - тангенциальные усилия исходного состояния пластины:
Tx0 = sxx0h,       Ty0 = syy0h,       S0 = sxy0h.
(2.11)

   Смысл (2.8)-(2.10) заключается в том, что уравнения устойчивости пластины можно формально рассматривать как уравнения задачи изгиба при действии фиктивной поверхностной нагрузки с компонентами X1, Y1, Z2. Однако смысл X1 и Y1 следует понимать с особой оговоркой. А именно, они не входят в выражения усилий и моментов, но фигурируют в первых двух уравнениях (2.6).
   Имея в виду формулы [3],
Nx = h3
12
j,    Ny = h3
12
y.
(2.12)

   Из третьего уравнения (2.5) после соответствующих дифференцирований заключаем, что члены выражений (2.8) и (2.9) с участием вторых производных функций j и y пренебрежительно малы по сравнению с теми членами этих выражений, в которых фигурируют третьи производные прогиба пластины w. Для наглядности рассмотрим одномерный случай, когда y = 0, а j и w зависят только от координаты x.
   Тогда
dj
dx
= - 12
h2
sxx0 d2j
dx2
Ю d2j
dx2
= - 12
h2
sxx0 d3j
dx3
.
(2.13)
   В силу этого
X1 = - h2
12
sxx0 ж
з
и
1 + 6
5
a55sxx0 ц
ч
ш
d3w
dx3
.
(2.14)

   Поскольку a55 = 1/G13, где G13 - модуль сдвига материала в плоскости 0xz, то член 6a55sxx0/5 в несколько тысяч раз меньше единицы и им можно свободно пренебречь. Таким образом, в выражениях (2.8) и (2.9) следует оставить только члены с участием производных прогиба. Выражения компонент фиктивной поверхностной нагрузки примут следующий окончательный вид:

X1 = - h
12
ж
з
и
Tx0 3w
x3
+ Ty0 3w
xy2
+ 2S0 3w
x2y
ц
ч
ш
= - h
12
Z2
x
,
Y1 = - h
12
ж
з
и
Ty0 3w
y3
+ Tx0 3w
x2y
+ 2S0 3w
xy2
ц
ч
ш
= - h
12
Z2
y
,
Z2 = Tx0 2w
x2
+ Ty0   2w
y2
+ 2S0 2w
xy
.
(2.15)

   Выражение фиктивной поперечной нагрузки Z2 совпадает с аналогичным выражением классической теории пластин [4].
   Подставив выражения усилий и моментов [3] в (2.5), с учетом (2.15) получим окончательные уравнения устойчивости анизотропных пластин при учете влияний поперечных сдвигов, нормального напряжения sz и изменения размеров:
(B11 + sxx0) 3w
x3
+ (B12 + 2B66 + syy0) 3w
xy2
+ (3B16 + 2sxy0) 3w
x2y
+ B26 3w
y3
-
- h2
10
й
к
л
(a55B11 + a45B16 + A1) 2j
x2
+ (2a55B16 + a45B12 + a45B66 + A3) 2j
xy
+
+ (a55B66 + a45B26) 2j
y2
+ (a45B11 + a44B16) 2y
x2
+ (2a45B16 + a44B12 + a44B66 +
+ A1) 2y
xy
+ (a45B66 + a44B26 + A3) 2y
y2
щ
ъ
ы
+ j = 0,
(B22 + syy0) 3w
y3
+ (B12 + 2B66 + sxx0) 3w
x2y
+ (3B26 + 2sxy0) 3w
xy2
+ B16 3w
x3
-
- h2
10
й
к
л
(a44B22 + a45B26 + A2) 2y
y2
+ (2a44B26 + a45B12 + a45B66 + A3) 2y
xy
+
+ (a44B66 + a45B16) 2y
x2
+ (a45B22 + a55B26) 2j
y2
+ (2a45B26 + a55B12 + a55B66 +
+ A2) 2j
xy
+ (a45B66 + a55B16 + A3) 2j
x2
щ
ъ
ы
+ y = 0,
j
x
+ y
y
+ 12
h2
ж
з
и
sxx0 2w
x2
+ syy0 2w
y2
+ 2S0 2w
xy
ц
ч
ш
= 0.
(2.16)

   С целью оценки влияния тангенциальных компонент фиктивной поверхностной нагрузки на значения критических напряжений пластины рассмотрим одномерный случай. Для критического значения сжимающего напряжения sxx0 = -s0 шарнирно-опертой по краям бесконечной полосы в рамках классической теории пластин получаем:
sкл0 = B11p2h2
12a2
,    s0 = B11p2h2
12a2 + p2h2
.
(2.17)

Здесь a - ширина полосы, sкл0 соответствует случаю неучета, а s0 - случаю учета X1. Как и следовало ожидать, s0 < sкл0, т.е. учет тангенциальной фиктивной нагрузки приводит к уменьшению величины критического напряжения. Поправка в процентах составляет:
D = sкл0 - s0
sкл0
100% = 100p2h2
12a2 + p2h2
.
(2.18)
   В случаях h/a = 1/5 и 1/4 имеем D » 3 и 5% соответственно.
   Для сравнения отметим, что поправка только от поперечного сдвига будет:
D = sкл0 - s0
sкл0
100% = 100a55B11p2h2
10a2 + a55B11p2h2
.
(2.19)

   При a55B11 = 3, h/a = 1/5 и 1/4 эта поправка составляет 10.6 и 15.6% соответственно. Следовательно, влияние тангенциальной фиктивной нагрузки в рассматриваемом случае примерно в 3 раза меньше влияния поперечного сдвига. Не исключено, что картина может существенно измениться в случае анизотропных пластин конечных размеров.
   В работе [2] уравнения устойчивости ортотропных пластин при учете влияния поперечных сдвигов выведены исходя не из системы (1.1), а из сранительно грубой системы (1.3). В этой работе тангенциальные компоненты фиктивной нагрузки X1,Y1 не учитываются, а для поперечной компоненты Z2 в обозначениях [3] получается выражение
Z2 = Tx0 ж
з
и
2w
x2
+ h2
16
a55 j
x
ц
ч
ш
+ Ty0 ж
з
и
2w
y2
+ h2
16
a44 y
y
ц
ч
ш
+
       + S0 й
к
л
2 2w
xy
+ h2
16
ж
з
и
a55 j
y
+ a44 y
x
ц
ч
ш
щ
ъ
ы
.
(2.20)

   В отличие от (2.15) выражение (2.20) содержит слагаемые с производными функций j и y, что автор считает результатом более последовательного подхода. Однако, как было показано выше, при использовании более точной системы (1.1) такие члены не появляются. Следовательно, появление этих членов является сдедствием того, что при выводе третьего уравнения устойчивости (2.16) величины w/x и w/y заменены их приближенными значениями (w/x - ¶uў/z)/2 и (w/y - ¶vў/z)/2 соответственно. Добавим, что при получении этого уравнения допущена арифметическая неточность, по причине которой в выражении (2.20) вместо множителя -h2/24 фигурирует множитель h2/16. Эта неточность и привела к заметному понижению значения критической нагрузки пластины в рассматриваемой конкретной задаче [2]. Отметим, что замена неточного множителя на его точное значение приводит к малому увеличению критической нагрузки [2] по сравнению с [5].
   3. Краевые условия задачи устойчивости пластин при учете изменения размеров формулируются как обычно [3], только с той разницей, что в этих условиях вместо первоначальных известных размеров a0 и b0 фигурируют неизвестные размеры a и b, соответствующие моменту наступления потери устойчивости. Это приведет к дополнительным осложнениям при решении конкретных задач.

     Институт механики НАН РА

Литература

     1. Новожилов В.В.  Основы нелинейной теории упругости. М.-Л. Гос. изд. технико-теоретич. литературы. 1948. 212 с.
     2. Томашевский В.Т.  - Труды VI Всесоюзной конференции по теории оболочек и пластин. Баку. 1966. М. Наука. 1966. С. 753-761.
     3. Амбарцумян С.А.  Теория анизотропных пластин. М. Наука. 1987. 360 с.
     4. Вольмир А.С.  Устойчивость упругих систем. М. Физматгиз. 1963. 880 с.
     5. Амбарцумян С.А., Хачатрян А.А.  - Изв. АН СССР. ОТН. Механика и машиностроение. 1960. N1. С.113-123.