Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Нелинейные задачи динамики цилиндрических композитных оболочек

..pdf
Скачиваний:
4
Добавлен:
12.11.2023
Размер:
9.34 Mб
Скачать

В подавляющем большинстве решений нелинейных задач тео­ рии оболочек используются упрощенные выражения для -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)