книги / Нелинейные задачи динамики цилиндрических композитных оболочек
..pdfВ подавляющем большинстве решений нелинейных задач тео рии оболочек используются упрощенные выражения для -a*î, т**:
„ |
1 |
dvi* |
, |
1 |
dAi |
|
|
Xi=xг=—:--- х-----1 т~к-----> |
|
|
|||||
|
At |
даi |
|
AiAj |
да,• |
|
(2.13) |
|
|
T*i=Ti, |
|
|
|
||
соответствующие модели типа Тимошенко, и |
|
|
|||||
|
1 |
двг |
|
1 |
àA< |
' |
|
|
Ai |
даi |
|
AiAj |
даj |
(2.14) |
|
|
1 |
àd, |
| |
1 |
дЛ, |
|
|
|
01 |
|
|||||
|
At |
дои |
|
AiAj |
daj |
|
соответствующие модели Кирхгофа—Лява. Возможность пренеб режения нелинейными членами в указанных формулах можно
обосновать, введя ограничение снизу на показатели изменяемости деформированного состояния оболочки. Так, при условиях
/ |
|
hj^kiAj |
duï°\ |
|
(2.15) |
|
|
1 dvi° |
1 |
0ТН0С11тельиая |
|||
iсогласно определению ^u= —ô ^7» taj= ^ô |
|
|||||
поправка, вносимая нелинейными членами в |
величину |
(2.5), |
||||
|
fv°\ |
fv02 \ |
величину т*=т*1+т*2— |
|||
имеет порядок 01-^-I +OJ-^ô1+0(0г), в |
||||||
O(0i) +O(02), |
а в величину х*{(2.12) —О |
|
Таким образом, |
точность формул (2.13), (2.14) понижается с увеличением углов поворота и с уменьшением показателей изменяемости %ц, %ц. От метим, что ограничения типа (2.15) обычно вводятся в техниче
ской теории (теории пологих оболочек).
Уравнения обобщенного закона Гука в пренебрежении квадра тами компонент тензора деформации ец и поперечным обжатием для упругой ортотропной оболочки записываются в форме
Gii=BiiSii+BijEjj\ |
= |
Gi3=Bi3&i3. |
(2.16) |
Уравнения равновесия оболочки в обобщенных усилиях и мо ментах имеют вид (1.42), (1.57):
9 |
<л т. Ï , |
д |
М Г* |
ï I дА{ T* |
dAi T* |
I |
~ д Т {А‘Т и) +~д^ ■{AiT |
+^ Г |
~дТ |
i,+ |
|||
|
+AiAikiriz+AiAi(X\+P\) =0; |
(2.17) |
||||
д |
(А2Г 1г)+- |
д |
(А\Т*2ъ) А\А%(k\T*i\+/г27’*22) + |
|||
да\ |
|
да2 |
|
|
|
|
|
|
+Л,Л2(Л'*з+Р*з) =0; |
|
(2.18) |
д |
пА- |
ÔA; |
ЛГн- |
дсц СAjM*u) + дщ |
оо.) |
даi |
|
- AiAjT** +ЛИ;(Р*,+Ф*3-) =0. |
(2-19) |
||
Здесь, согласно (1.71)—(1.74) |
|
|
|
Га =сн ( e,+^-0i2) +С„ (e;+ÿ.0/) |
+С,ли®(»,°+е,); |
||
Г*,,=С66(Yi+Yi+ |
+Ca»j°(Oi°+0j); |
(2.20) |
|
î’*i3=Ci3(üi°+0i) +C«.0i( et+-^-0i2) +С,д (е;+у6/) + |
|||
+C660j(Yi+Yj+9i0j); |
|
(2.21) |
|
M"u=Dii [x*i+Aj (Ei+~2-0i2)] +Dij [x*j+Éj( |
0j2)] ; |
||
АГу=£>66[t*i+T*j+k,(Yi+Y;+e,0j)] ; |
(2.22) |
||
Г* =CMBi+Vi°) +С,з(0,+о,°) (Yj-o,°0j) +Du(kt+ks) X |
|||
X (0i+î)i°) (xi+Xj) +Dj3(k(+kj) {Oj+vj'>)Ti-kjvi’>(Djjy.j+DtiKi) + |
|||
+kjD^vp(Xi+Xj) + {ki—kj)DSsQj{Xi+Xj). |
(2.23) |
Величины Cki, Dm —мембранные и изгибные жесткости ортотропного однородного либо многослойного (с симметричной укладкой слоев по толщине) пакета, вычисляемые по формулам [21, 24,
180].
Выражения (2.22), (2.23) можно упростить, следуя подходу, используемому в технической теории (теории пологих оболочек). Пренебрегая в (2.22) мембранными деформациями и используя
для х%-, т*г формулы (2.13), имеем
М**.=DuXi+Dijxy, M*ij=Dee Ы+xj)л |
(2.24) |
Кроме того, основываясь на малости величин yj. Ui°0j по сравне нию с единицей и учитывая введенные выше ограничения на по казатели изменяемости деформированного состояния оболочки,
можно существенно упростить выражение для Т**а (2.23), отбро сив в нем нелинейные слагаемые. Таким образом,
7*;-Сй(Л,+1»,°). |
(2.25) |
Входящие в уравнения (2.17)—(2.19) поверхностные нагрузки Х*и Х*з, массовые силы Р**, Р*3, моменты поверхностных и массо
вых сил Ф*и F*u согласно формулам (1.46), (1.43), |
(1.61), (1.60), |
|
при учете малости е,-, у*. 0«2> |
по сравнению с единицей имеют |
|
вид |
|
(2.26) |
X*i= (Xi+-Xr) + (Хз+-Хз”)и£°; |
х*3= (X3+-Xz~) + №+-*г)01+ №+-Яг)е2; |
(2.27) |
|
h/2 |
/1/2 |
|
P*i= J [1 -\-az(ki-t-kj)]Fida3, |
Р*з= J [\-\~a^(ki+k2)]Fzdaz; |
|
-h/2 |
-h/2 |
(2.28) |
ф*,=А{ (^-ь+^г) +4“ |
(^ +^> + (^+-Х,-)Тг] }; |
F*i= J {(Fj-Fzv?)[\+az{ki+kj)]+atFj{%i+Kj)}abdaz.
-h/2 |
(2.30) |
В результате подстановки формул (2.6)—(2.8), (2.13), (2.20), (2.21), (2.24), (2.25), (2,26)—(2.30) в систему (2.17)—(2.19) по следняя может быть записана относительно пяти неизвестных
ФУНКЦИЙUi°, W3°, Vf0.
При использовании кинематической модели Кирхгофа—Лява имеем три уравнения равновесия (2.17), (2.18). Выражения для обобщенных усилий в приближении технической теории оболочек записываются в виде (1.78), (1.83):
Гц=с„ (si+Y0!2) +Cij(e,+Y.0j2 ) ;
T*ü=Сбб (y*+Уз+ ; |
(2.31) |
|
|
T*iz=— \ —-— {Aj[Duv.i-\- (Z)i2+2D66)Kj]}— |
|
AiAj ucci |
|
- ~a!a~ да1 КД|г+2Рбб)«.+Р»кЛ+Сцв,( Ei+ y 6**) + |
|
+Cjjfl< (ej+Y0i2 ) +C660j(Yi+7i+9<9i) |
(2.32) |
Поверхностные нагрузки X*u Х*з, согласно (1.48), при учете
малости ей уи 9г2 по сравнению с единицей, задаются форму лами
Х*,= (Xi+-Xr) - (Х3+-Хг) 9г |
(2.33) |
и (2.27). Для массовых сил Р*и Р*з сохраняются выражения (2.28). Для моментов поверхностных сил при условиях (1.13), которые уже использовались ранее, вместо (2.29) можно приме
нять упрощенную формулу
ф%=-| №++*г)-
Моменты массовых сил Р\-, согласно (1.65), в рассматриваемом приближении имеют вид
/1/2
F*i= J {(Fi+F&j) [1+а3(ki+kj)]+a3Fj(xf+%,•)}azdaz.
-W |
(2.35) |
Путем последовательных подстановок (2.31), (2.32), (2.27), (2.28), (2.33)—(2.35), (2.6), (2.14) в систему (2.17), (2.18) по
следняя может быть записана относительно трех неизвестных функций иД 1/3°.
2.2. НЕЛИНЕЙНЫЕ УРАВНЕНИЯ ТЕХНИЧЕСКОЙ ТЕОРИИ ОРТОТРОПНЫХ КРУГОВЫХ ЦИЛИНДРИЧЕСКИХ ОБОЛОЧЕК
Рассмотрим частный случай —замкнутую круговую ортотропную цилиндрическую оболочку радиуса R. Введем координаты: а\=х —вдоль образующей, а2=у —по окружности, а3=z —по нормали к срединной поверхности (координата г положительна в направлении внешней нормали). При этом Ai=A2=\, &i=0,
В дальнейшем будем использовать обозначения: Ui°=u,
U2°=V, Ü3°=W, Vi°=yx, u2°=yy.
Уравнения равновесия оболочки в рамках нелинейной техни ческой теории с учетом деформаций поперечных сдвигов, полу
чаемые из (2.17)—(2.19), имеют вид |
|
||||
|
<эг„ .+J IiL +^.l+p*I=0; |
|
|||
|
дх |
ду |
|
|
|
|
дТ*22 |
дТ*12 |
+Х*2+Р*2=0; |
|
|
<97*13 |
ду |
дх |
|
|
|
^7*23 |
1 |
(2.36) |
|||
~df~+~dy1 |
Ж Г22+^*з+Р*з=°; |
||||
0Л1*“ + -^ rL-TVt+F*î+<!>tl=0; |
|
||||
дх |
ày |
|
|
|
|
дМ*22 |
дМ'ц |
-7**+7*1+Ф*1-0. |
|
||
ày |
дх |
|
Обобщенные усилия и моменты определяются по формулам (2.20),
(2.21), (2.24), (2.25), в которых, |
согласно |
(2.6), |
(2.8), (2.13), |
||||
ei=- |
du |
dv |
w |
Yi=. |
du |
У2=- |
du |
дх |
62=“ду |
R |
дх |
ày |
Величины |
Ei, 82, Yi» уг, 0i, 02 |
определяются по |
формулам (2.37). |
|||||
Кроме того, согласно (2.14), (2.12), |
|
|
|
|||||
_ |
^'01 |
_ __ д02 |
|
_ _ д02 |
_ _ |
d0i |
||
Ул~ |
дх |
' *2_ |
ду |
|
Tl |
дх |
Т2 |
ду |
Подставив далее (2.37) в |
(2.42), (2.43), находим |
|
(2.43) |
|||||
|
|
|||||||
- |
du 1 ( dw Y Л &v w 1 ( dw Y |
|
||||||
£,, = ôF+T |
du |
e22=â7+X T \ dy ) |
|
|||||
|
|
~ |
dv |
dw |
dw |
|
|
|
|
|
812 |
dy |
dx * dx |
dy |
d2w |
|
|
|
|
d2w |
|
d2w |
|
|
||
|
|
X2=~ |
^ |
' : |
T,=T2=_ |
dxdy |
|
При рассмотрении оболочки, имеющей начальный прогиб ш°,
согласно подходу, |
изложенному в |
1.7, формулы |
(2.44) |
заменя |
||||||||
ются на |
- |
да |
|
1 [7 dw |
V |
( dwQ\21 |
|
|
||||
|
|
’ |
|
|||||||||
|
8,1 “ дх* 2 |
I \ |
дх |
/ |
\ ox |
/ J |
|
|||||
|
dv |
|
|
|
|
1 Г/ dw \2 |
/ |
du |
\21 |
^ |
||
622=-^--t— |
R |
|
|
2 K |
<9ÿ ) |
( |
/ -I: |
|||||
|
ày |
du |
dv |
dw dw |
|
|
<?*/ |
|
|
|||
|
~ |
|
dw° dw° |
|
|
|||||||
|
6,2 |
dy* dx* dx dy |
dx |
dy |
|
|
||||||
|
|
d2{w—w°) |
|
|
d2(w —w°) |
|
||||||
*> = ------ 13---- » У>2~----------- |
|
|||||||||||
|
|
|
dx2 |
|
d2(w —w°) |
of |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
dxdy |
|
|
|
|
||
Используя (2.2) и |
(2.45), для деформаций |
|
получаем следующие |
|||||||||
выражения: |
|
|
|
|
|
|
|
|
|
|
|
|
du |
1 |
Г( dw |
Y |
( dw° \21 |
|
d2(w —wQ) |
|
|||||
В”= Ж +Т |
K - â r) |
- ( л г) |
|
|
-----; |
|||||||
dv |
w—wQ |
1 Г/ |
dw Y |
I dw° Y1 |
|
|||||||
e22= d f+~R ~~+Т1\~дГ/ '\~д<Г1 J" |
|
|||||||||||
|
|
|
|
|
d2{w—w°) |
|
|
|
(2.46) |
|||
|
|
|
|
~z |
àf |
|
: |
|
|
|
||
du |
dv |
dw |
|
|
d2(w—w°) |
|||||||
dw |
_ dw° |
dwQ |
|
|||||||||
e,n=~dÿ*~dx*~dx |
dy |
dx |
|
dy |
|
2 |
dxdy |
|
Исходя из приведенных уравнений и формул нетрудно вывести уравнения равновесия оболочки в перемещениях. Подставляя по
следовательно (2.41), (2.45) в первые |
три |
уравнения |
системы |
||||||||
(2.36), получаем |
д2и |
dw d2w |
dw° d2w° |
\ |
|
|
|||||
|
/ |
|
|
||||||||
|
' |
dx2 |
dx dx2 |
dx |
dx2 |
/ |
|
|
|||
|
Г â2v |
1 d(w—w°) |
dw |
ô2w |
|
||||||
dw° |
12L dxdij |
R |
dx |
à2v |
ày |
dxdy |
|
||||
d2w° |
1 |
/ à2u |
|
d2w |
dw |
|
|||||
~~ dy |
àxdy J |
66' dy2 |
dxdy |
dxdy dy |
|
||||||
d2w° |
|
dw° |
t dw' d2w |
âw° |
à2w° |
\ |
|
||||
~~àxâÿ |
|
dÿ~+~dx |
àf |
dx |
dy2 |
|
l + |
|
|||
/ |
d2u |
|
dw |
+**i+P*,=0; |
d2w® \ |
|
|
(2.47) |
|||
|
à2w |
|
dw° |
|
|
|
|||||
Cl2 \ ~dxày |
dx |
dxdy |
dw |
dx |
dxdy / |
+Сгг X |
*] |
||||
P d2v |
| ^ |
ô[w —ш°) ^ |
à2w |
|
àwQ |
à2w0 |
|||||
X L"cty2"” |
R |
|
дУ |
~дУ |
àÿ2 |
|
dy |
|
dy2"J |
+ |
|
„ / |
à2u |
, |
d2w |
dw |
à2w° |
|
dw° |
|
|||
+Сб6\ |
àxdy |
|
dx2 + |
dx2 |
dy ~~~~d~x2 |
|
dïT + |
|
|||
dw |
d2w |
dw° |
d2w° \ |
|
|
|
|
(2.48) |
|||
+ -1Г ~Ш Г --дГ -Щ Г ) |
|
|
|
|
dw° - - Н ày
dw Г ди , |
1 / dw |
V |
ду 1 дх |
дх ■) |
|
dw Г ди до |
dw |
dw |
~дГ L~àf+ дх +- дх |
ày |
dw° |
dw° ]} =X*3+/>*3+^-(F*2+0*2) +^ _ (f.|+®., |
дх |
ày |
|
(2.49) |
Уравнения равновесия в смешанной форме можно получить
непосредственно из (1.88), (1.91). Учитывая, что, согласно (1.86), для цилиндрической оболочки
Тп = |
d2F |
Т22~ |
d2F |
d2F |
(2.50) |
ду2 |
дх2 |
Т12=Т21= —дхду |
и используя формулы для изменения кривизн и кручения (2.45), приходим к известной системе двух уравнений:
D u * 1” - * ' |
+2(Д„+2Д66) |
дх2ду2 |
d4'*~W°] + |
|||||||
|
дх4 |
|
|
|
|
ду4 |
||||
|
d2F |
|
d2w |
|
d2F |
d2w |
d2F :+2 |
d2w |
d2F |
|
1 дх2 |
|
дх2 |
|
ду2 |
ду2 |
дх2 |
дхду |
дхду |
||
|
=**3+Р*3+- ^ |
|
(^ +cr2)+^ _ (F.I+(I). |
(2.51) |
||||||
22 |
d*F |
|
|
|
|
d*F |
+Лц |
d*F |
d2(w —îü°) |
|
дх* —h (2Л12+Лбб) |
|
дх2ду2 |
ду* |
дх2 |
||||||
/ |
d2w |
\ 2 |
/ |
d2w° |
\ 2 |
d2w |
d2w |
d2w° |
d2w* |
|
+ \ |
дхди |
/ |
' |
дхду |
I |
дх2 |
ду2 |
дх2 |
ду2 |
|
|
|
|
|
|
|
|
|
|
|
(2.52) |
При выводе (2.52) учтено, что в начальном напряженном состоя нии уравнение совместности деформаций имеет вид
К =0.
На основе систем (2.47)—(2.49) и (2.51), (2.52) построены решения рассматриваемых в последующих главах нелинейных за дач деформирования несовершенных ортогропных цилиндрических
оболочек.
2.3. ЛИНЕАРИЗОВАННЫЕ УРАВНЕНИЯ ДЛЯ ОРТОТРОПНОЙ КРУГОВОЙ ЦИЛИНДРИЧЕСКОЙ ОБОЛОЧКИ
Вывод линеаризованных уравнений устойчивости тонкостенных ортотропных оболочек изложен в 1.8. Для частного случая — круговой цилиндрической оболочки, согласно (1.98), (1.100),имеем
следующие уравнения:
дТ'п |
дГп |
+P\+X'i =0\ |
дТ'а |
дТ'а |
Г'м |
+Р'2+Х'2=0; |
|
дх |
ày |
дТ'хз ^ |
|
ду |
дх |
R |
|
|
|
дт'п |
.- ^ - + Р ' з+Х'3=0; |
(2.53) |
|||
|
|
дх |
ày |
А |
|
|
|
|
|
дМ'п |
, дМ'21 —Г^з+Ф^+^г—0; |
|
|||
|
|
дх |
ду |
|
|
|
|
|
|
дМ*22 |
дМ'а |
— +Ф,1+F'1—0. |
|
||
|
|
ду |
дх |
|
|
|
|
Здесь Т'ц, Т'ц, T'i3 —«дополнительные» усилия, возникающие в оболочке при отклонении от докритического равновесного состоя ния. В рамках технической теории они имеют вид (1.102):
T'i^Cne'i +C^e^; |
Г'-^Сгге'г+С^е'ь |
Т'\2=Т/2\= С^{у'\-\гуГ2)\ |
|
r/i3=Ci3(0/i +Y/*); Г/2з= С23 (О^+У у)\ |
(2-54) |
||
М\I=D\\y!\ -\-D\2Kr2\ |
М'22=022^2+Di2vJ\; |
М\2=M'2i =D66(x'i+т'г). |
Величины X'i, Х'3 представляют собой «дополнительные» по верхностные силы, воздействующие на оболочку только при откло нении ее от докритического равновесного состояния. В соответ ствии с (1.99) для цилиндрической оболочки они записываются
вформе
Х\=-^~ (r„V2) + - |- (TfflVs) +?1(e'l+e'a)+bd‘W,2+‘M'*)> дх оу
Х'2= - ~ |
(Тю0® '^ -—- (7'hVi)+— |
6,2+^2(e/i + e'2) + |
||
|
|
+ôcUriY,i +<73Y'I/); |
(2-55) |
|
Х ' з |
= ( М |
оу |
' г ) “ -^г—e'i+^3(6'i+e'2) + |
|
дх |
|
а |
|
+бс(^10,1+^20/2) • Моменты поверхностных сил Ф'г выражаются в виде (1.101).
Входящие в (2.54), (2.55) величины е'*, у'и 0'г-, у/и т'г в рамках технической теории определяются по формулам (2.37), (2.38).
Если предположить, что единственными массовыми силами являются силы инерции, то для цилиндрической
ветствии с (1.104), имеем
P'v--- рл( |
д2и' |
i_ |
h2 |
à2yx |
|
|
dt2 |
Г |
12R |
dt‘ |
|
|
|
Р'2--- p/t( |
d2v' |
J_ |
h2 |
à^.j |
b |
(2.56) |
dt2 |
T |
12R |
dt2 |
|||
P'3=-p/i- d2w' |
|
|
|
|||
/i3 |
àWv |
dt2 |
d2v' |
|
|
|
1 |
|
|
||||
F'i = -p “ПГ |
dt2 |
|
dt2 |
b |
(2.57) |
|
Л3 |
/ à2y'x |
|
d2u' |
|
||
R |
). |
|
||||
F'i = - рЧГ1 |
dt2 |
dt2 |
|
где p —плотность материала оболочки.
При использовании модели Кирхгофа—Лява, исключая Т\з и Т\з из последних двух уравнений системы (2.53), получаем
|
дГп |
, |
дТ'п +Р\+Х\ =0; |
||
дТ'22 |
дх |
|
ày |
i'X2 \ |
|
дТ\ |
|
дМ'22 |
дМ'| |
) +Р'2+Х'2+ |
|
ày |
дх■ Ч г ( |
ày |
>Х |
> |
|
|
+ ^-(®'.+^i)=0; |
(2.58) |
|||
|
д2М\х |
д2М'2j |
а2(М'12+м'21) =Р'з+*'з+ |
||
|
дх2 |
ду2~ |
дхду |
|
|
|
+ -4 - (Ф'>+га)+ -4- (ф',+г,). |
||||
|
дх |
|
|
ду |
|
Взаключение рассмотрим выражения для деформаций, усилий
имоментов более общие, чем используемые в технической теории оболочек. Из формул (1.2) и (1.3) в линейном случае вытекает
(штрихи в дальнейшем опускаем):
Bii- |
8г +азХг |
еij=- |
Vi+a3Tf |
+ yj+а3т. |
(2.59) |
1+аз&« |
1+а3А:,- |
1+а3£; |
Разлагая эти выражения в ряды по степеням а3и сохраняя все члены до первой степени включительно, для цилиндрической обо
лочки находим
8ц=81+азХЬ 822=82+азХ2; Bi2=Y +a3*. |
(2.60) |