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

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

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

0г (1.6) пренебрегается тангенциальным перемещением щ°. Соот­ ветственно этому формулы (1.80) принимают вид

____ 1___/

'

1

ди3°

\

1

dAj

dus0

 

Хг

Ai

даi

Ai

dai

*

AiAj2

даj

daj

( 1...84)

_

1

d

/_1

ди3°

\

1

dAj

duzQ

 

Ai

dai

'

Aj

daj

'

A^Aj

daj

dai

 

Во-вторых, в первых двух уравнениях равновесия (1.42) пренебре­ гается перерезывающими усилиями 7’*г-3, что дает

 

 

 

 

 

 

 

dAj

 

 

 

 

+AiAj(X*i+P*i)=0

dai

■?jj+

 

 

 

 

 

(1.85)

(здесь учтено, согласно (1.78), что Т*ц= Тц, T*ij=Tij).

и, следуя

Если в уравнениях

(1.85)

положить X*i=0, P*i=0

[94], ввести функцию усилий F такую, что

 

 

 

 

1

д /

1

dF

\

1

dAj

dF

 

 

и Aj

daj

\ Aj

daj

' + AMj

àai

dai

_

1

d2F

1

 

dAj

dF

1

dAj

dF

13

AiAj

daidaj+ AiAj2

dai

daj

Ai2Aj

daj

dai

то в приближении технической теории (1.85) удовлетворяются тождественно. Выражение для Г\-3 принимает вид

T*iz—Tiz+TuQi+TijQj=-

—--------{Aj\Pii%i-\- (Z)i2+2Z)6e)«i]}—

1

dAi

AiAj

dai

Г 1

à

!

1

dF

\

 

 

~

' d^T 1 (Dii+2D№)**+DH*A + 6i L17

daj

\

Aj

daj

1 +

1

a n

Г

1

à2F

J __

dAj

dF

 

AMj

àat àat l +6i

AiAj

ôaida,

AiAj2

dat

daj

 

+ J —

dAi

Éf_1

+ г ;+Ф*5-.

 

 

(1.87)

 

Ai2Aj

daj

àat J

 

 

 

 

В результате система уравнений равновесия сводится к одному уравнению:

àïT {л 7 ^ № " * ‘+

 

|

\

- x

S

r [(Z)l2+

Г d

(

/ 1 dF

1

àAi

dF 1

+2066)х,+ДиК2]+0, [ _

— —

) +-j^- ôai

ôai J +

+02 [-

1

d2F

 

1

дА2 dF

\_dM__dF_

 

Ai

daida.2

‘ AiA2 dai

da2 + Ai2 da2 àai ■и

 

 

+ - + { 4 ---—{A,[û22x2+ (Oi2+2D66)hi]>-

(1-88)

 

 

да2

l A2

o<X2

 

 

 

Г à I

1

àF \

 

 

1 dAi

1'(°»+20^

+

ы

 

~ ъ Ж

+e4 âsr (

 

/ +

^ ü 'U ^ F l

Г

 

1

iPF

1

àAi

dF [

A22

da.2

da2 J ^ 1*-

 

A2

da2<5ai

Л2Л1

<5^2

<5ai

\ ] _

1

A4*

dF

11

Г

 

d

i

1 _c?F_ ,

<5aidn,

<5da,i

+ Л/1"2

oaidcci

c/da2

JJ/

ÆlLЛ‘*- ' о'агda2

^'

Л2AÏ da2* Л1A-

/ J

L

dai

' Л1 dai

'

A2

da2

da2 J

dai

 

+«■,)] +

 

 

 

+-+[A, ( f .+Ф*,)] +A,A2(X*3+P*3) =0, da2

содержащему неизвестные функции м3° и F. Второе уравнение мо­ жет быть получено из уравнения Гаусса (3.32), приведенного

в [196]:

д

i

 

д

/л -

ч

 

dЛ1*

1 „

^ei2

dai

 

1 г—

(Л2е22) - h— Ё12_ —Л1—

l Ai L

 

 

 

 

da2

2

da2

д

Г 1 Г

д

,л -

ч

 

dЛ2-

1 „

 

<^ei2

 

 

 

 

 

 

--г—e^-__ Лл

dai

+ d^T1 а2 1в5-и,в,|)‘

dai

2

 

=x4ii42(T2—Х1Х2—Л2У1 —/г1X2),

- - <5А2 11 , "в1,л Г И +

- dЛ1 11 822 da2 -1*

(1.89)

где eü=ei+j0i2; ei2=Yi+y2+0i02. Выражая эти величины сог­ ласно (1.78) через усилия,

g

_Cj}Tii—CijTij

Тх2

(1.90)

СцС22—С,22

Сбб

 

подставляя их в (1.89) и используя в преобразованиях первые два уравнения равновесия, приходим к уравнению

 

—12

 

[

* (. > ü k r . ) -

( ^22‘

 

2

à

1’

1

/ L dai '

A, ôa,

22 /

 

1

/

A« \

da.2 '

Л2 а^ТГи) J + И

h-AI2- - ^ - ) X

Г à

,/

1

ôA, _

\

d

/ 1 <?A2

\ 1

X l-^r-

\ Л2

 

1

(а Г ^ Г

11)J +

Lоаг

<Ja2 Г“ )

dai

 

д

 

 

 

 

(1.91)

+ dot]Н ^ ^

^

22Г“ + ( Л|2+~

^

Г1|1 } +

+^ { ^ ^

[

л,,7',,+( л,2+4

1 ) 7'22]} =

где

=А\А2 (т2—kjx2 —k2x\—k\y,2),

 

С»

 

 

Си

 

1

Ац=.

Au—

 

СцС22~С122

С11С22—Ci22

Л66 =

 

 

 

 

 

 

Сб6 (1.92)

Подстановка в (1.91) выражений (1.86) и (1.84) дает уравне­ ние относительно функций ы3°, F, однако записать его не удается в столь же компактной форме, как для изотропного материала. Отметим, что в последнем случае Л22—12—бб/2=0, Ац—А \2—

-Л66/2=0, 2412+Лбб/2=^ц=Л22=^, и введением оператора

1 Г д ГА2 dF 1 | д ГAj dF j \ ^ ^~ А\А2 I да\ LА\ да\ J да,2 I- Л2 да2 -! /

(1.91) приводится к виду

(1.93)

V2VV=Eh(т2—Х1Х2—k2y.i—k\%2).

Для «весьма пологой» ортотропной оболочки, когда Л) и Л2 при дифференцировании ведут себя как постоянные [21], из

(1.84), (1.86) следует:

 

 

1

д2и3°

_____1

d2uz°

 

Kt

Ai2

дсц2

*

 

Л1Л2 daida2

 

 

1

д2р .

т - т

21

1

d2F

(1.94)

 

Aj2

daf

 

*

12

Л1Л2

d<xida2

 

a уравнение (1.88) принимает вид

 

 

 

(1.95)

где

 

LiUsq+VrF—£(u3°, F) =Z,

 

 

 

 

 

 

 

 

 

D22

d4

£,= Du

d4

2(Z)12+2D66)

___

 

 

dai4

 

 

Л,2Л22

 

dai2da22

Л24 da24

 

VR=-

£2

d2

+

âîi

d2

 

 

 

dai2

Л22 da22

 

C(fufz)=' A\2A2

I d2U

d2f2

 

<?2/i

 

\ dai2

d2f2

dai2 da22

 

 

 

 

 

d2U

 

 

 

 

 

 

 

 

да\да2

да\да2

 

 

 

Z=X \+P\+~JL(F\+<b'2)+ ~ -? —(F>l+<S>*l).

А\ да\

Л2 о&2

■Соответственно уравнение (1.91) в рассматриваемом случае запи­ сывается в форме

£2F-VBu30+y£(«30,tt30)=0, (1.96)

где

_ Л22

дА

2Л12+Л66

д4

Ли

д4

2“ "Л7

da,4+

А?А22

да\2да.22

+ Л24

да24 '

Уравнения (1.95), (1.96) были получены в [24].

1.7. УЧЕТ НАЧАЛЬНЫХ НЕСОВЕРШЕНСТВ ФОРМЫ ОБОЛОЧКИ

В процессе изготовления оболочек или монтажа конструкций, •содержащих оболочечные элементы, неизбежно возникают общие или локальные технологические несовершенства, приводящие к от­ клонению формы срединной поверхности от «идеальной». Наличие таких несовершенств с общих позиций может быть описано по­ средством введения начальных деформаций и, как следствие, на­ чальных напряжений. Нелинейная теория оболочек с учетом на­ чального напряженно-деформированного состояния построена в работах X. М. Муштари [194, 195]. Она изложена также в моно­

графии [196]. При этом начальные деформации гц°, в\2° опреде­ ляются через начальные перемещения иг°, н3° по формулам теории

пологих оболочек.

В дальнейшем будем рассматривать частный случай, когда на­ чальные несовершенства носят чисто изгибный характер, т. е. Ui°=0 вдоль всей поверхности оболочки. Отклонение формы сре­ динной поверхности от «идеальной» полностью задается полем нормального перемещения w° (си, аг). Следовательно, начальные деформации определяются в рамках технической теории по фор­

мулам

а

1

 

 

 

1 / 1

dw° \ 2

 

 

eiio=eio+-I e ,o '= ^ + T ( ^ 7- ^ r )1

 

 

а

 

 

 

 

1

dw°

dw°

 

 

S,20=Yi0+Y2,+ei%°= A.Aj

da.

Sa,

 

0

1

d

(

1

dw°

 

1

dAi

dw°

Ki°-

Ai

даi

'

Ai

даi

AiAj2

daj

dctj ; ( 97)

.

1

d

(

1

dw°

 

1

dAi

àw°

T

A,-

àai

\

At

даj

Ai2Aj

don

âoи

Будем считать, что начальные «искажения» формы оболочки не вызывают начальных напряжений. Поэтому при начальных ус­

ловиях на перемещения Ui°\ =0, из°|*=о=ш° и отсутствии в момент /=0 внешних воздействий на оболочку все напряжения при /=0 должны быть равны нулю. Уравнения равновесия при

этом должны удовлетворяться тождественно.

Таким образом, для учета начальных геометрических^несовер­

шенств следует заменить га на (е™—е**0), ei2 — иа (&12—S120),. Хг — на (Хг-Хг°), Т{ —нa (ъ-Тг°) либо оставить их в прежнем

виде, понимая под ег*, В12, и*, тг «дополнительные» величины, воз­ никающие под действием нагрузки.

1.8. ЛИНЕАРИЗАЦИЯ УРАВНЕНИИ РАВНОВЕСИЯ

Проблема вывода линеаризованных уравнений устойчивости оболочек рассматривалась многими авторами. Краткий обзор при­ мененных подходов и библиография по этому вопросу приведены в работе В. В. Болотина [70], указывающего на наличие боль­ шого числа вариантов записи «параметрических» («нагрузочных») членов, обусловленного разнообразием путей получения уравне­ ний устойчивости. В этой связи особо отметим работы [148, 196], где линеаризация осуществляется исходя из нелинейных уравне­ ний равновесия оболочки. Аналогичный подход использован при выводе уравнений устойчивости в [134, 135].

Для получения линеаризованных уравнений будем исходить из уравнений равновесия (1.42) и выражений для обобщенных усилий (1.41). Рассмотрим частный случай, когда в докритическом состоянии возникают только нормальные усилия Тц°, что харак­ терно для нагружения оболочки продольными усилиями и внеш­ ним давлением. Деформациями в докритическом состоянии будем

пренебрегать. Обозначим через Т'и,

Т\3 «дополнительные»

усилия, возникающие в оболочке при

отклонении от докритиче-

ского равновесного состояния. Следуя методу, изложенному в [134, 135], заменим в (1.41) Тц на Тц°+Т'ц, подставим получен­ ные выражения для обобщенных усилий в (1.42), вычтем урав­ нения равновесия докритического состояния и проведем линеари­

зацию. В конечном результате получим следующие уравнения устойчивости:

9 {А,Ги)+- 9

дА ■

 

Т'ц+

•(A,T'jt) Н—rr—T'ij—■dAj

даi

&ccj

dccj

àtxi

(1.98)

 

+AtAikiT'i3+AiAiP'i+X'i^0-,

 

(AjT'i3)

(AiT’n)-AiAiikiT'ii+kjT'ij) A-AiAjP'а+ЛГ',-0.

где

х' ^ ~ к ( А т ) +~ t~ w m + -ÏJ- w -

- дА •W

+Лi i4 Ггг0е'г+Л*4j[?* (e'<+e'j) +ôc (fty'i+q3Vi0')];

 

(1,99)

^'3=~dâ—

—даp HiîW’Oy —i4iy42(^i7’n06/2+/î27,220e/ri)+

 

+^Иг[?з(е/|+£'2)+ôc (<7iO'iЧ-^О’г)];

qi=Xi+—Xr;q3=Xz+—Xz~; ôc=l и Опри следящей н консерва­ тивной нагрузках соответственно. Уравнения (1.98), как видно, от­

личаются от обычных уравнений равновесия линейной теории обо­ лочек только наличием «нагрузочных» членов X'i, Х'з. Формулы (1.99) являются частным случаем выражений, полученных в [134, 135, 148]. Величины Х'и Х'з, Р'и Р'ъ представляют собой «допол­ нительные» поверхностные и массовые силы, воздействующие на оболочку только при отклонении ее от докритического равновес­ ного состояния.

Линейные уравнения для моментов согласно (1.57) имеют вид

д

dAj

М'п-

даг {AiM'u) +~ è r {AiM,l<)+" й г т tj

даг

(1.100)

-AiAiT'iz+AiAjWi+F'i) =0,

 

где <D'j и F'j — «дополнительные» моменты поверхностных и мас­ совых сил, возникающие при отклонении от докритического рав­

новесного состояния. Ограничиваясь

в (1.60), (1.61) линейными

членами, находим

 

71/2

h

F'<=J [1+а3(ki+kj)]a3da3\

Ф^=-^-{Х,++Х:Г)-

-*/*

(1.101)

Усилия T'a, T'a, T'iz и моменты М'и, М'ц выражаются через «до­ полнительные» деформации по формулам линейной теории обо­

лочек:

Ги=Сиг'i+Cije'j;

Т'ц=С66(y'i+y'j);

(1102)

T'i3=Ciз(0'г+Vi0');

M'ii=DiiK'i+£>,•jh'j;

M'ij=D6e{T'i+Vj).

 

При использовании гипотез Кирхгофа—Лява усилия Г',з опреде­

ляются из уравнений (1.100).

Для решения задач устойчивости оболочек на основе уравне­ ний, записанных в смешанной форме, необходимо положить Р,'=0,

X'i=0, выразить усилия Т'ц, Т'ц через функцию F по формулам (1.86) и осуществить аналогичные проведенным в 1.6 преобразо­ вания, отбрасывая нелинейные члены. В результате (1.88) прини­ мает вид

 

даtai

I Ai dai

{Л2[0,1к',+ (О,2+2Д66)*'2]>-

Ai

оА2

[(^12+21)бб)/С/1+^22^,2] J-

д

Г 1

д

dai

да2

I А2

оа2 X

 

 

 

 

1

dAi

 

 

X {Ai[D22y/2+ (D12+2D66)x'i]} — A2

da2 [(^12+2Z)66)^2+

 

 

д

( *

dF

)+ —

dA2

dF

 

 

да2 V а2

dai >

Ai

dai

dai ] -

 

 

f 1

dF

) |

1

àAl

dF

 

 

 

1 А\

dai

f

A2

da2

da2 i+

dai

[Л2(Г2+Ф'2)] +

da2 [Ai (F'i +O'i)] +AlA2P',+X'3=0.

 

 

 

 

 

 

 

 

(1.103).

В уравнении совместности деформаций (1.91) следует опустить нелинейные слагаемые в правой части.

При решении задач динамической устойчивости усилия Тц° являются функциями времени. Учитывая, что силы инерции F'i=

= —p

d2u'i

д2и'ъ

i

*

 

Fг——p -^2

(P —

плотность материала оболочки;:

u'it и'г — «дополнительные» перемещения), по формулам (1.1),. (1.43), (1.101) находим

Г d2Uj0'

h2

(ki+ki)

d2Vi°'

];

P'z=-çh

dW'

1 РН di2

+ 12

Г

 

dt2

 

д2и?' j

dt2 ;

F'i= -

h3

 

+ {ki+kj)

(1.104)

12

L

Ôt2

dt2

Приведенные в данной главе уравнения и формулы позволяют решать нелинейные задачи деформирования, задачи статической и динамической устойчивости тонкостенных ортотропных оболочек как с учетом деформаций поперечных сдвигов, так и на основе кинематической модели Кирхгофа—Лява.

ГЛАВА 2

У р а в н е н и я н е л и н е й н о й т е о р и и о р т о т р о п н ы х ц и л и н д р и ч е с к и х о б о л о ч е к

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

жит все необходимые исходные формулы и уравнения, используемые далее в монографии.

Теория тонкостенных цилиндрических оболочек разрабатыва­ лась многими авторами. Обстоятельные исследования по линейной теории представлены в монографиях [21, 24, 94, 123, 143, 209, 254, 322, 441]. Нелинейная теория изотропных цилиндрических оболочек изложена в [196]. Уравнения нелинейной теории для случая анизотропного материала рассматривались в монографиях [92, 224]. Не останавливаясь на историческом аспекте развития теории цилиндрических оболочек, укажем лишь на обзорныестатьи [10, 13, 107, 115, 134, 136, 138, 154, 199, 203, 204, 291], содержа­ щие обширную библиографию по этой проблеме.

2.1. ОСНОВНЫЕ УРАВНЕНИЯ НЕЛИНЕЙНОЙ ТЕОРИИ

ТОНКОСТЕННЫХ ОРТОТРОПНЫХ ОБОЛОЧЕК

Пусть аь сх2, а3 — координаты точки в произвольной криво­ линейной ортогональной системе координат. Совместим коорди­

натную поверхность а3=0 со срединной поверхностью оболочки. Координатные линии ai, a2 направим вдоль линий главных кри­ визн этой поверхности. Соответствующие орты обозначим еь е2. Примем, что направление орта е3 совпадает с направлением внеш­ ней нормали к срединной поверхности оболочки.

Обозначим через Wi(ai, аг, а3), ы2(ai, а2, а3), и3(аь а2, а3) проекции вектора перемещения произвольной точки оболочки на

направления ei, е2, е3 и примем, что они являются линейными функциями координаты а3 [208]:

Wi =Wi°+a3üi°; M3=w3°+a3u30; i=l,2.

(2.1)

Здесь U]°(ai, а2), «г0(ai, а2), ц30(аь а2) — проекции вектора пере­ мещения произвольной точки {ai, a2} срединной поверхности обо­ лочки на направления еь е2, е3; i>i°(ai, a2), a2°(ai, a2) — функции, характеризующие повороты сечений, проходящих через нормаль

ккоординатной поверхности, при деформировании оболочки. Компоненты тензора деформаций, согласно [210], определя­

ются формулами (1.2), (1.3). Подстановкой в них (2.1), в пре­ небрежении величинами порядка a3£i по сравнению с единицей, получены выражения для деформаций оболочки (1.12). Дальней­ шие их упрощения основывались на следующих допущениях:

1) показатели изменяемости вдоль координатных линий ai, a2 характеристик деформированного состояния и коэффициентов пер­ вой квадратичной формы срединной поверхности оболочки Л* ог­ раничены сверху и удовлетворяют условиям (1.13);

2) деформации удлинений и сдвигов — величины более высо­ кого порядка малости, чем углы поворота нормальных сечений, в соответствии с чем отбрасываются их квадраты, взаимные про­ изведения и произведения на остальные неизвестные функции.

В результате были получены формулы

(1.15),

согласно ко­

торым

 

 

 

 

гц=гц+а.гУ*i\

ei2=ej24-a3 (t*i

2);

(2.2)

где

ei3=üi°+6i,

 

(2.3)

I

А

 

 

A

 

(2.4)

i—вН—2"®i2' ®l2=7l^"Y2+-ei02Ï

я*,-=и<+0<Х(з; r*i=ti+0jXi3;

(2.5)

1

дщ°

1

dAi

Uj°-{-kiUz0',

 

Ai

даi

AiAj

daj

dAi

(2.6)

Уi=

1

 

 

1

dai

 

 

Ai

AiAj

daj

 

 

Ôi=

1

duj

--kiUi°\

 

 

Ai

dai

 

Ai

dvi° t

1

dAi Vj°+kiV3°;

(2.7)

dai

AiAj

daj

dAj

 

 

1

dvj°

 

1

(2.8)

 

Ai

dai

AiAj

daj

 

 

 

 

1

dv3°

—kiVi0.

 

 

 

Ai

dai

 

 

 

В этих формулах индекс /=2 при i= 1 и /= 1при i=2. Деформа­ ция езз, согласно (1.11), выражается в виде

833=f30+4-(oi0’+»20!+°30!)- (2,9)

В главе 1 был рассмотрен вариант нелинейной теории ортотропных оболочек, основанный на кинематических допущениях Кирхгофа—Лява. При этом в соответствии с трактовкой В. В. Но­ вожилова [208] предполагалось, что при определении изменения направления орта ез величиной поперечных сдвигов можно пре­ небречь. Это приводит к условию (1.16), из которого, с учетом малости и3° по сравнению с единицей, получается соотношение

(1.17):

Vi°~-Qt. (2.10) Кроме того, было получено выражение для функции v3° (1.27)

ü3°« —(*1181 +0282) —r-0i2(l +ai) — O22 (1+П2) » (2.11)

содержащее комбинации коэффициентов Пуассона ортотропного материала (1.24). При рассмотрении материалов с сильной ани­

зотропией, для которых лчз<1, V23<1, выражение (1.27) мало

отличается от

предложенной

В. В. Новожиловым [208] фор­

мулы

(1.19).

(2.10),

(2.11),

ограничиваясь линейными относи­

С

учетом

тельно удлинений и

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

членами, выражения для у*и т** (2.5) преобразуем к виду

 

 

1

<?0г

1 dAi