книги / Устойчивость и колебания трехслойных оболочек
..pdfкий смысл— он является коэффициентом приведения длины. Дей ствительно, согласно выражениям для х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Г |
|
||
|
{ |
\ |
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). Эти уравнения в рамках принятых |
гипотез учитывают |
полную инерционную силу и представляют собой довольно слож ную систему параболического типа. Заметим, что если считать несущие слои мембранами (Л = ^2—0), то система становится ги перболической. В общем случае эта система уравнений может быть сведена к одному уравнению восьмого порядка с четными производными, коэффициенты которого могут быть проанализи рованы. Однако здесь мы ограничимся учетом только главного инерционного члена, имея в виду, что влияние остальных инер-
ционных членов на первые частоты незначительно. Для этого в
(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
Формула (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