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

книги / Устойчивость и колебания трехслойных оболочек

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

кий смысл— он является коэффициентом приведения длины. Дей­ ствительно, согласно выражениям для х2 и k по (I. 174) крити­ ческую силу можем определить так:

N

кр

Djt2

»А2Я2

к2л2

(1.175)

к 2

 

 

 

 

 

 

 

где = /Я-1.

Остается для каждого конкретного случая опирания стержня определить минимальное значение Я, соответствующее нетри­ виальному решению уравнения (1.171).

N ,

 

л <

■ N

N ,

 

J L -

 

W /

UV U

: w , СCCL

-----)

\

 

 

{

\

W.

 

 

77

тг,

 

f 77

г

I

 

 

1=0 х=хшх*=о

(1

х=хБ=х*=о

^ х а=хш-ш Ч х-кхс=х Ч ш=о х=хг=хш=о

Х -к Х ^ -Ш ^ Ч )

Х~кхЧт=Хш=0

Х~кХп=Х1=Хш=0 Х-кХБ=^ х шч

Х-кХг=ХБ-киХ^=Хш=0

0,5

г.

1 + Щ к ^ \ 1 +

+1,ЬЗл1Щ1+к)\

 

 

Рис. 8. Критическая сила NHр продольного сжатого стержня при пяти видах краевых условий на концах стержня (£ = 0 и £ = я )

Запишем общее решение уравнения (1. 171) в форме

х (5):= Аг

+ А2 C O S

+ Л3

+

А

 

 

V

-|-Л4 Ch V$-|-Л5$-|-Лв.

(1.176)

Удовлетворяя однородным

краевым

условиям

при £ = 0 и £= л,

приходим сначала к однородной системе линейных уравнений относительно Л,, для существования нетривиального решения ко­ торой требуем далее равенства нулю ее определителя, в резуль­ тате получаем трансцендентное уравнение относительно пара­ метра Я[у2 и х2 через Я2 выражаются посредством (1.174)]. Мини­ мальное значение Я будет соответствовать критической силе.

В табл, на рис. 8 приведены пять случаев опирания стержня. Для каждого случая даны краевые условия и значения ЯтьЗа­ метим, что в двух последних случаях формулы для Ят ш практи­ чески пригодны при ограничениях Ф С 1. Кроме того, разни-

41

ца между первым и пятым случаем состоит в наличии у послед­ него бесконечно жестких диафрагм

(* „ .„ = 1+ 0 ,6 4 * j / - Н - [1 + 1 ,4 3 V kb (1 + * ) ] ) .

Критическая сила, приведенная на рис. 8, равна

В случае е ф е 0 на торцах

M = N (e во)

(1. 177)

и уравнение (1.168) при неоднородных краевых условиях (1.177) имеет единственное решение для произвольных N за исключени­ ем N—Nxp. Так, для свободно опертого стержня прогиб в сере­ дине пролета равен

/6»

Г ( 1 + /6X2) v 4

( 1 —

/6v2) Х-*

-I - 1

(1.178)

Х2 + V2

1

,

.

1

 

 

 

c o s у

Х я

c h

\ я

J

 

 

 

 

 

2

 

при стремлении А—►1 он интенсивно возрастает, становясь не­ ограниченным. В частности при # = 0 наибольший момент и наи­ больший прогиб равны (х= 1/2)

-^гаах—N (e 0О)

cos у Хя

(1.179)

® m ax=(*-<?o)

[ c o s - y Хя ]

Формально выражения (1.179) совпадают с соответствующи­ ми выражениями для однородных стержней, однако в силу со­ отношения

А *=— - — ,

(1.180)

1 — /6x2

V

7

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

N.кр

<

1 ,

(1.181)

1 + /6Х2

 

 

отсюда

 

 

 

1 + * (1 — (О2 о

8 >

(1.182)

Х2 =

 

 

 

42

что и доказывает предыдущее утверждение, поскольку в случае однородного стержня в (1.179) будет фигурировать не %, а р.

Аналогичная ситуация складывается и при наличии попереч­ ной нагрузки д (х ), только теперь отлична от нуля правая часть уравнения (1. 132)

X 4 ' - . — - ? - -

X lv — — Х" =

------.

(1.183)

А»

А»

D nm

 

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

 

 

 

 

 

 

Д/3(1 +АХ2)

 

(1 + jfeX2) tg аи

 

w (т)- 16£>(1 + 2AvX2 + &А2Х4)

[

и3

 

 

 

(1 — Av2)th у v

 

 

 

 

 

 

 

 

 

и3

-]

 

 

(1.

184)

м

ш

= ___ p jj± ± m

 

 

I (1 + Ш 2 ) ji i L

 

 

V 2 /

 

4 (1 +2А»Х2 + »А2Х4) | V

'

и

 

поэтому

tg

и -> оо при Я-»-1 (и = я'к/2; ti=jtv/2).

 

Полагая 0 = 0

и v = oo, из (1. 145) находим

 

 

 

 

 

 

I

\

Pl*{1+ АХ2) г (1

ах2)tgц — и .

 

 

 

 

[

2

J

16D

[

и*

J ’

 

 

 

м

[

1 \

+ m

tg в

 

(1.

185)

 

 

 

 

 

 

 

 

[

2

/

4

и

 

 

 

Формулы, соответствующие однородному стержню, получим,

положив в

(1.185)

k = 0 .

 

 

 

 

 

 

 

11.

КОЛЕБАНИЯ ТРЕХСЛОЙНОГО СТЕРЖНЯ

 

Вернемся к рассмотрению динамических уравнений (1.51) —

(1 .53). Эти уравнения в рамках принятых

гипотез учитывают

полную инерционную силу и представляют собой довольно слож­ ную систему параболического типа. Заметим, что если считать несущие слои мембранами (Л = ^20), то система становится ги­ перболической. В общем случае эта система уравнений может быть сведена к одному уравнению восьмого порядка с четными производными, коэффициенты которого могут быть проанализи­ рованы. Однако здесь мы ограничимся учетом только главного инерционного члена, имея в виду, что влияние остальных инер-

ционных членов на первые частоты незначительно. Для этого в

(1 .51)— (1.53) следует положить Е = 0, D = 0, в результате приходим к следующей системе уравнений:

dN_ - В * —

= 0 ;

 

(1.186)

дх

dt2

 

 

 

дН

д2М

оя

д2w

(1. 187)

Q3= Q ;

----- -—

В *

------

q=Q.

дх

дх2

 

dt2

 

 

Вводя функцию у и в соответствии с (1.65) функцию перемеще­ ний х, приходим к двум независимым динамическим уравнениям

 

 

d2v

1

д2*/

 

(1.188)

 

 

'dl2 ~'~a2

'dt2

 

 

 

 

 

D

ш

д2 \д*Х

\ Eh

J l ( l _

 

 

о -

Р

дх2 ) дх*

1 a2

dt2 \

■ р dt2 I х-

у

Здесь введена скорость распространения продольных волн

e = j / - f - .

(Ы 9 0 )

Рассмотрим свободные поперечные колебания стержня.. Преобразуем уравнение (1.190), используя безразмерную ко­

ординату 1=пх/1 и записывая функцию х (1. 0 в форме

Х(£, t ) = - ^ - X (S)expШ .

(1.191)

71

Вводя в полученное уравнение безразмерные параметры

Ь 2Я 2 _

12/40)2

(1.192)

р/2 ’ Ш _ Л2а2я40 ’

получаем обыкновенное дифференциальное уравнение относи­ тельно X (|) (штрих означает дифференцирование по £)

A-vi.

а>2

 

а>2

(1.193)

X Iv - S - X " - \ ---- 5^ х = о .

kb

»

'

kb

 

Положив далее

 

 

 

 

 

* ( £ )= e x p / s £ ,

 

(1.194)

приходим к кубическому характеристическому уравнению, соот­ ветствующему (1.193)

s3

sa

о>*2s

<o*2

(1.195)

b

о.

 

kb

 

 

44

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

 

 

 

 

«1 =

— X2; sa=p.2; s3= v2.

 

 

(1.196)

Теперь

 

 

 

 

 

 

 

 

 

 

 

 

 

H«+V*= 1+ **— ;

|» V = —

,

 

(1.197)

 

 

 

r '

 

k b

 

 

W kb

 

 

K

и для

определения

корней р,2,

 

v2 имеем квадратное уравнение

 

 

 

 

 

1 + Щ 2

 

0*2

 

 

(1.198)

 

 

 

 

 

kb

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если учесть, что

 

 

 

1 +

 

 

 

 

 

 

 

 

<о*2= Х 4

 

 

(1.199)

 

 

 

 

 

1 + 6X2

 

 

то по

(1.198)

 

 

 

 

 

 

 

 

 

 

а

 

1 + щ г Г

/ ~ 1 + * Х 2 ( 1 — 3 » ) + # Ш

1 ,

 

^

~

2kb

[

V

1 +

6 X 2 ( 1 + &) +

 

|

 

 

 

1 + k m f J . /

1 + 6 X 2 ( 1 3 » ) + 6 2 » X 4 ~|

(1. 200)

 

3

 

 

 

V

~ ~

2kb

L

 

 

1 +

6 X 2 ( 1 + » ) +

* 2 » X 4

J

'

Таким

образом

и здесь корень X можно принять в качестве

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

 

 

Q«o

 

К I

( 1. 201)

 

 

 

 

где /о — приведенная длина стержня, равная

 

 

 

 

10— Х~Ч.

 

(1.202)

Согласно

(1.181)

при Ф<0,25, р2> 0 ,

v2> 0

поэтому общее

решение уравнения

(1. 173)

можно записать в форме

*

(5)= Л,

+ Лаcos X?+

Л5

+

 

 

А

 

V

 

 

- ^ „ c h v S - M s - ^ + A c h m .

(1-203)

 

 

 

и-

 

 

Остается, удовлетворив однородным краевым условиям, сос­ тавить характеристическое уравнение для параметра X. Краевые условия остаются теми же, что и в случае устойчивости стержня, но, кроме свободно опертого стержня Я = 1, значение параметра X существенным образом зависит от значений Л и в меньшей сте-

пени от -ft, при этом приходится отыскивать корни весьма гро^ моздких трансцендентных уравнений.

В заключение заметим, что при •0—+0

V-'

Х2

оо,

(1.204)

1 + Ш

 

£9

 

и в общем решении (1.204) следует положить Л5= Л в = 0. Здесь трансцендентные уравнения значительно проще и их первые кор­ ни могут служить хорошим приближением к корням более об­ щих уравнений.

12. ПОСЛЕКРИТИЧЕСКОЕ ПОВЕДЕНИЕ ТРЕХСЛОЙНОГО СТЕРЖНЯ

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

Пусть s натуральный параметр нейтральной линии изогнуто­ го стержня, отстоящей на расстоянии hi2 С13 от средней линии заполнителя, тогда уравнения изгиба стержня будут (1.60) — (1.61)

dH

ds

Q3;

(1.205)

d?M

N — w = 0.

ds2

ds*

Предполагая деформации стержня малыми, согласно (1.47) имеем

Н= О у\ -У — da

1.1— 9 dx

(1.206)

M = D \у —

L dx Q3— Ohby.

Здесь x — кривизна стержня, вычисленная по точной формуле

X

(1.207)

46

Из уравнения, связывающего Н и Q3, имеем

 

 

у

_d *_ _ / 1 __

— V -

(1.208)

 

Р ds

V

Р

ds2 /

'

Это уравнение будет тождественно удовлетворено,

если поло­

жить

 

 

A2

d2

 

 

 

 

 

 

 

 

 

. = Л ~ £ - £ Л Ч Г ;

 

 

 

 

а = -1 — »

Ш_ <№

 

(1.209)

 

 

 

 

Отсюда

 

 

р

ds

 

 

 

 

 

 

 

 

(\ — — — W ==

d$w

 

( 1. 210)

 

ds*

 

V

Р

<*«*/

 

 

 

теперь для полного момента М имеем выражение через 'У

M = D ( 1 — — W .

(1. 2 1 1 )

\Р ds2/

Используя второе уравнение равновесия (1.206), получим уравнение, связывающее V и w,

D ( 1

1 \

d2lt'

■N - ^ - w = 0.

( 1. 212)

I

р ds*)

ds2

ds*

 

( 1

h* d*\

 

зуя ( 1. 2 11 ), приходим к нелинейному уравнению относительно w

\

vh2 d2 \ 42

D 1

рds*)dts*

^ \

= 0

(1.213)

р rfs2 / ds2

v •

или, разлагая обратную величину радикала в степенной ряд, по­ лучим его в форме

\

* L \ J L [ * 1 L \ \ + — ( — f 4 .Л / ' ^ 4+ . . Д +

р ds*! ds*

\ ds*

[

2

\ ds J

~

8 \ ds J 1

J 1

 

+

Лф

- * 1

1

М

^

=

0.

(1.214)

 

 

V

P

ds*)

ds*

 

 

 

Найдем решение уравнения (1.215), считая, что стержень свободно оперт. Если сила N ненамного превосходит критичес-

кую силу NKр, то естественно ожидать, что упругая кривая бу­ дет иметь вид, близкий к кривой, полученной при исследовании линейной задачи, поэтому положим

w ( s ) = f sin - у - ,

(1.215)

и для отыскания постоянного параметра / используем метод Бубнова. Как известно, метод Бубнова решения дифференци­ альных уравнений состойт в том, что искомая функция аппрок­ симируется рядом координатных функций с произвольными па­ раметрами, которые определяются из требования ортогонально­ сти результата подстановки этого ряда в уравнение к каждой из координатных функций.

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

Ограничиваясь в уравнении (1.214) одним нелинейным чле­ ном, после подстановки выражения для w в форме (1.215) най­ дем

D

T f

{ '

 

(1 + И ) sln ^

+ ° т г 4 / * Х

 

X

(1 +

9Л&) sin i p

— 7 V / -| -(l+ *)s in ^ H -^ JP.

(1.216)

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

k = я2Л2/p/2.

 

 

Вычисляя интеграл

 

 

 

 

 

 

 

 

 

i

 

sin j- d s — 0 ,

 

 

 

^

F

 

найдем

 

 

 

 

 

 

 

 

 

 

 

(1+5 /2)(1+^)/_iV(1+^)/==0-

(1,217)

Вводя Эйлерову критическую силу

 

 

 

 

 

 

N

 

= D

я 2

l + k b

 

 

 

 

 

^

 

 

 

 

-'V KP----

-

1 + k

 

 

 

 

 

 

 

/2

 

 

запишем

(1.217) в форме

 

 

 

 

 

 

 

 

*кр

 

/ а^

'

(1.218)

 

 

 

 

8/2

1 )

 

отсюда

 

 

 

 

 

 

 

 

 

 

 

_

2 / 21

'

N

 

1.

(1.219)

 

 

/

 

 

 

 

 

 

 

я

 

Nкр

 

 

48

Рис. 9. Малый элемент трехслой­ ной оболочки произвольного вида:
/—первый слой; 2—второй слой: 3—тре­ тий слой (заполнитель)

Формула (1.219) по форме совпадает с соответствующей формулой теории однородных стержней, однако следует иметь в виду, что Л/'кр зависит от k и Ф, т. е. от сдвига заполнителя и структуры стержня. Эта формула дает для однородных стерж­

ней достаточно точные результаты

вплоть

до f = 0,2

I, т. е. до

значения силы Af= 1,045 NK-p.

 

 

 

Г л а в а

2

 

 

ОБЩАЯ ТЕОРИЯ ТОНКИХ УПРУГИХ ПОЛОГИХ ОБОЛОЧЕК

ПРИ КОНЕЧНЫХ ПРОГИБАХ

 

В этой главе вариационным методом

получены

основные

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

коэффициент Пуассона v. Теория, не содержащая последнего допу­ щения, при предпосылках, ука­ занных выше, изложена в работах {12, 13, 14].

Будем считать оболочку пологой, различием радиусов кри­ визны слоев пренебрегаем. Принимая за исходную срединную поверхность заполнителя, отнесем ее, учитывая пологость обо­ лочки, к декартовой системе координат х г, хг. Положительную нормальную координату z будем отсчитывать в сторону внешней нормали к исходной поверхности. Называя несущий слой, рас­ положенный со стороны внешней нормали, первым слоем, слой

со стороны

внутренней

нормали — вторым, а заполнитель —

третьим (рис.

9), введем

обозначения:

49

Лл(А3= 2 с)-то л щ и н а слоев (* = 1 ,2 ,3 ) ; h — толщина стенки оболочки;

Е к и V* — модуль упругости и коэффициент Пуассона слоев ( * = 1 ,2 ,3 ) ;

О — модуль поперечного сдвига заполнителя; *//*/= 1 ,2)— кривизны и кручение координатных линий;

/г — размер оболочки в направлении осей JC,( / = 1 ,2);

w— нормальное перемещение точек исходной поверхнос­ ти (прогиб);

u f — тангенциальное перемещение точек заполнителя в направлении оси xl (i= 1 , 2);

a t ( d f/d z ) — закон распределения поперечных сдвигов по толщи­ не заполнителя в направлении оси *,(/== 1 , 2)

[<*i= <*f to , ■*■«), / = / ( * ) ]

Ъц— символ Кронеккера: 8» = 1 , 8^ = 0 при i ^ j . Приведенный коэффициент Пуассона примем в виде

Введем безразмерные жесткостные характеристики и безраз­ мерные толщины слоев

 

 

- l

 

 

Ь =

1 — V*2

fifeAft |

tk= hkh~l

(* = 1 ,2 ,3 ) . (2.2)

1 -v *»

 

 

 

Очевидно, для yk и /* имеем равенства

 

 

ft-1

ft-1

 

(2.3)

 

 

 

Для более компактной записи формул удобно ввести осреднен ный модуль упругости

£ = (1 — V2) А- 1

i h i 1 -v*2

(2.4)

Тогда

 

 

 

£ft*ft

Eh

(2.5)

l — v*2

1 — V2

 

Непосредственно задается, dfldz, поэтому определим

функцию

f(z ) равенством

 

 

 

 

^ / ( 2) r f z = - i - * 3aX

(2 . 6)

50