3 11 |
СВО БО ДН Ы Е К О Л ЕБА Н И Я |
351 |
Затруднительно также отыскание общего решения уравнения (1. 23). Поэтому ограничиваемся классом решений, являющихся решениями следующего дифференциального уравнения:
и удовлетворяющих граничным условиям для и>, т. е. условиям непрерывности и однозначности w на сфере.
Из (1.23), согласно (1.24), для частот собственных колебаний трансверсально изотропной сферической оболочки получим
2 _ gE п __р\ с'г (X— I )2+ |
1 + х*х |
|
(1.25) |
<0 |
ТоД 2 { |
’ |
( l + A * X ) ( |
X - l + v ) |
|
|
|
|
Формулой (1.25) собственные частоты оболочки представлены |
с помощью пока еще неизвестного параметра X. |
|
|
Значения параметра X определим из уравнения (1.24), |
которое |
подробно изучено. |
Решение |
уравнения |
(1.24) |
ищется |
в виде |
|
|
|
sin к В, |
|
|
Условие однозначности |
функции W (а., (3) на сфере |
требует, |
чтобы к было нулем |
или целым числом (& = 0 , 1 , 2 ,.. .). |
|
Подстановка значения W в (1.24) с учетом обозначения cos а = |
= х приводит к известному уравнению |
|
|
|
|
- * ’ > - £ ] + ( * - А |
) ^ |
0 |
|
|
(к = |
0 , |
1, 2, ...), |
|
|
|
собственными функциями которого являются присоединенные функции Лежандра.
Нетривиальное ограниченное решение этого уравнения суще
ствует лишь при собственных значениях |
|
Х„ = |
п(п + |
1) |
(ге = |
0 , 1, 2, ...), |
(1.26) |
которым соответствуют следующие собственные функции: |
|
^ |
(*) = |
( ! - * T |
* £ EP. (* ) |
(1-27) |
|
(* = |
0 , |
1, 2, |
...), |
|
где полином Лежандра |
|
|
|
|
Ря(х) = 2яп ! ^ К * 2— !)*]■
352 |
К О Л ЕБА Н И Я И УСТОЙЧИВОСТЬ |
О БО ЛО ЧЕК |
[ГЛ . I l l |
|
Возвращаясь |
к |
уравнению |
(1.24), |
находим |
систему |
его |
решений: |
|
|
|
|
|
|
|
|
к— 0: |
W0(а, |
(3) = |
Рпcos а; |
|
|
|
к= |
1: |
W_! (а, |
Р) = |
P^l) (cos а) sin |3, |
|
|
|
|
|
Wj (а, |
Р) = |
(cos а) cos Р; |
|
|
|
к — п: |
1¥_я (a, ft) = |
PM (cos а) sin п р, |
|
|
|
|
|
W„ (а, |
Р) = |
Р№ (cos а) cos п р. |
|
|
|
Напомним, что полиномы Лежандра Р п(а:) в интервале |
( — 1, |
+ 1 ), т. е. в интервале изменения а (0, я), имеют п нулей. Присоеди ненные функции Лежандра P ik) (х) имеют п—к нулей.
Как известно, sin п р и cos пр обращаются в нуль на 2п мери дианах, а Р (пк) (х) (в интервале 0 ^ а ^ тс) — на п—к широтах. Следовательно, этими меридианами и широтами (узловыми линиями) оболочка разбивается на участки, внутри которых Wn не меняет свой знак. Это значит, что параметр \ определяет вид формы коле баний оболочки. Чем больше параметр X, тем меньше размеры по луволн.
3. Осесимметричные свободные колебания анизотропной кру говой цилиндрической оболочки. Осесимметричные колебания круговой цилиндрической оболочки, в каждой точке которой име ется лишь одна плоскость упругой симметрии, параллельная сре динной поверхности оболочки, согласно (1.3.36) и исходным по ложениям настоящего параграфа описываются следующим урав нением:
д*ш |
I о |
|
. |
| Д 4 |
т0 , |
(1.28) |
ds |
4- 2т |
-— г -\-nw -\- -г.-----— h |
' |
' |
1 |
Du |
g |
|
где
§ 1] |
СВО БО ДН Ы Е КО Л ЕБА Н И Я |
353 |
Решение уравнения (1.28) ищем в виде |
|
w{s, |
t) = Акs |
i n s |
i n <s>kt, |
(1.30) |
где A k — неизвестные |
постоянные, |
<ой — частоты |
колебаний. |
Принимая (1.30), полностью |
удовлетворяем граничным усло |
виям, а удовлетворяя |
уравнению (1.28), получим |
для частот |
колебаний следующую |
формулу: |
|
|
Полученная здесь формула для определения частот колебаний рассматриваемой оболочки позволяет ставить задачи оптимального проектирования, конечно для рас
сматриваемого частного случая. Изме |
|
няя материал оболочки или, что очень |
|
важно, ориентацию главных направ |
|
лений упругости материала оболочки, |
|
можно |
проектировать |
оболочки |
|
с требуемыми |
частотами |
свободных |
|
колебаний. |
|
|
|
|
|
Пусть, |
например, оболочка изго |
|
товлена из ортотропного материала, |
|
главные направления упругости ко |
|
торого в |
каждой |
точке |
оболочки |
|
с главными геометрическими направ |
|
лениями оболочки |
(s и 0) могут со |
|
ставлять произвольный угол. |
В ча |
|
стности, |
пусть |
в |
каждой |
точке |
|
оболочки главное направление упру |
Рис. 61. |
гости составляет с главным геоме |
|
трическим |
направлением |
s |
угол ? |
|
(рис. 61). Пусть, далее, оболочка составлена из однородного материала (Са:=hBik); тогда для коэффициентов /тайга, входящих
|
в расчетную формулу |
(1.31), будем иметь |
|
|
т __о ^16(Д11^26 — |
|
|
|
Вп{ВпВвв-Щв) ' |
(1.32) |
|
__12Л- (ВцВ%2 — В \г) 8 еs ~Ь ^.В^В^ф ^ъ — Дц^?б — ^22^16 |
|
|
|
/г- |
В п (ВцВвд — В j’e) |
|
Имея значения упругих постоянных материала оболочки в глав ных направлениях упругости, с помощью формул (1.1.42) илй (II.9.38) легко определить значения коэффициентов Bik, входя щих в (1.32).
Имея значения всех В(к, с помощью формул (1.32) определим значения /тайга при различных значениях угла <ри далее с помощью
354 |
КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК |
[ГЛ. III |
формулы (1.31) найдем значения частот колебаний оболочки при различной ориентации ортотропного материала в теле оболочки.
Рассматривая эти результаты, замечаем, что, изменяя лишь ориентацию материала в теле оболочки, мы можем существенно изменить динамические характеристики оболочки. Безусловно, полученные результаты относятся кчастному случаю осесимметрич ных колебаний оболочки при определенных граничных условиях, при определенной геометрии оболочки и при определенно задан ном ортотропном материале. Однако, несмотря на это, сделанные выше выводы имеют достаточно общий характер и могут считаться корректными.
§ 2. Некоторые вопросы статической устойчивости анизотропных оболочек
За последние годы вопросы устойчивости оболочек привле кают всеобщее внимание исследователей. Имеется громадное ко личество работ, посвященных этому весьма актуальному вопросу современной теории оболочек.
В настоящей главе приводится лишь ничтожная часть этих исследований, представляющая интерес с точки зрения теории анизотропных слоистых оболочек.
Общая система дифференциальных уравнений устойчивости анизотропных слоистых оболочек достаточно сложна и громоздка и не всегда может быть использована для решения многочислен ных задач, представляющих интерес с точки зрения приложений. Однако для выяснения многих вопросов теории и для решения кон кретных задач устойчивости общая система уравнений может быть существенно упрощена.
Упрощенная система дифференциальных уравнений устойчи вости анизотропных оболочек может быть получена на основании уравнений теории пологих оболочек (см. гл. I, §§ 5 и 14).
Полагая, что основное напряженное состояние оболочки до потери устойчивости является безмоментным, и ограничиваясь точностью классической теории пологих оболочек, получим для компонент интенсивности фиктивной поверхностной нагрузки
известные представления: |
|
|
|
|
|
Z* — — Т\ч — Т\жг — S \ |
X* = 0, |
У* = 0. |
(2.1) |
Заменяя в соответствующих уравнениях |
(например, |
(1.5.14), |
(1.5.21), |
(1.5.27), (1.5.46), |
(1.14.12), |
(1.14.30), (1.14.41) |
и др.) |
Z , X , Y |
значениями Z * , |
X * , Y *, |
получим уравнения локаль |
ной устойчивости соответствующей анизотропной слоистой обо лочки.
Картина несколько изменяется при рассмотрении задач устой чивости по уточненным теориям. В этом случае представления (2.1) не всегда приемлемы, ибо возникают вопросы, связанные
§ 2J СТАТИЧЕСКАЯ УСТОЙЧИВОСТЬ АНИЗОТРОПНЫХ ОБОЛОЧЕК |
355 |
с учетом искривления и нарушения нормальности поперечного элемента оболочки. Эти вопросы требуют дополнительных ис следований, более корректных, чем это сделано до сих пор неко торыми авторами.
Здесь, сознавая незавершенность исследований в этой области, при рассмотрении задач устйчивости оболочек с учетом попереч ных сдвигов все-таки будем пользоваться представлениями (2.1),
которые в первом приближении позволяют выяснить влияние по
перечных сдвигов на устой |
|
чивость |
анизотропной обо |
|
лочки. |
|
|
|
1. Устойчивость |
пологой |
|
ортотропной цилиндрической |
|
панели. Рассмотрим |
задачу |
|
статической устойчивости по |
|
логой ортотропной |
круго |
|
вой цилиндрической панели |
|
(i?1=co, |
R 2= R ), |
сжатой |
Рис. 62. |
вдоль образующих равномер |
|
но распределенной нагрузкой |
|
с интенсивностью р, приложенной по торцам оболочки (рис. 62).
Пусть оболочка изготовлена так, что главные направления упругости материала ее совпадают с координатными направле ниями а, [3(А = \ , 2?=1). Пусть, далее, оболочка шарнирно оперта по всему контуру (а = а , а=0; $=Ь, [3=0).
Очевидно, что в рассматриваемой оболочке в безмоментном
состоянии появляются следующие внутренние силы: |
|
|
Ц = — Р, Л>= 0, S° = 0 . |
|
|
(2.2) |
Тогда, согласно (2. 1), получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.3) |
Разрешающую |
систему |
получим |
из |
системы |
уравнений |
(1.7.56). Принимая |
(1.7) и |
вместо Z подставляя |
|
Z*, |
получим |
1 d*F , АЗ г <?2 |
<?2 |
|
| |
, |
|
№ ( |
д* . |
|
i T ^ '+ T 2' L ^ z'11+ |
|
|
|
|
|
То |
|
|
|
+ «55 | г) |
- т ] * + |
р t 1 - |
-JJ (% A i + |
|
|
-f- аиЬ22) -(- а44а55щ |
(L nL 22 |
£f2)J |
|
= О, |
(2.4) |
(L aL B - Щ |
|
|
(amLn + |
aHL22) + |
|
|
|
|
■ |
hi |
,, |
, |
LJll)J 0a2 |
0, |
|
|
~ T aU «55 100 '^11^22 |
|
где для операторов L ik имеем формулы |
(1.9). |
|
|
|
356 |
КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК |
[ГЛ. III |
Решение системы (2.4) представим в виде |
|
F(n, |
Р) = |
sin |
sin |
К = |
|
f |
(a, |
p) = |
f mnsinXIasinX2p, |
Х2 = - ^ , |
|
где m a n — целые числа, которые показывают числа полуволн
по направлениям а и р соответственно.
Принимая (2.5), удовлетворим условиям шарнирного опира ю т по всему контуру.
Подставляя значения F и Ч*- из (2.5) в систему уравнений (2.4), получим относительно коэффициентов Fmn и Чгтп однород
ную систему двух алегебраических уравнений. Приравнивая нулю определитель этой системы, будем иметь для определения крити ческой силы следующую формулу:
£ - = т я г l b [ЦЬи+ 2X1X2612+ х^ 2+
|
+ |
(аиЦ + |
orJ ‘D (Ьпь2, - ft?2) ] + |
^ |
, (2.6) |
где, |
как и раньше, |
|
|
|
|
ди„= |
1. + |
Tf7 (as6fcn + |
аМ |
+ flaOssт (ьн6и — ьи). |
|
Ьц — Х»Яц -f- Х|566, fe22 = |
Х|В22 -{- Х?506, Ь12 = |
Х(Х2 (В 12 -f- Soe)- |
Формула (2.6) получена на основании уточненной теории, учитывающей явления поперечных сдвигов. В частном случае классической теории, т. е. когда пренебрегают явлениями, свя занными с поперечными сдвигами, полагая a44=l/G 23= 0 , я55= = 1/G13= 0, получим
р ;„= - щ - (цьп + |
2х м 2+ чъп) |
|
^Рввыо______Х( |
(2.8) |
|
|
|
Д2 |
ЬпЬТ1— Щг |
|
В случае трансверсально изотропной оболочки, когда в каж |
дой точке плоскость изотропии материала |
|
оболочки |
параллельна |
срединной поверхности оболочки, в силу (1.15) получим |
|
. |
D |
(Xf + |
Xj)2 |
, |
E h |
X? |
|
(2.9) |
~ |
X2 |
1 + h*a? (Х| + Ц ) |
' |
Д2 |
(Х| + |
XI)2 * |
|
|
|
где, как всегда, |
|
|
... |
l b 2 |
|
Е |
|
|
(2.10) |
D |
Eh? |
|
|
|
|
|
|
|
|
|
|
|
|
|
12 (1 — V2) » |
П ~ ' 10 о2 |
(1 — v2) G ' " |
|
|
Рассматривая формулы (2.6)— (2.10), |
|
замечаем, |
что критиче |
ская сила вообще и величины, входящие в зти формулы и представ-
§ 2] СТАТИЧЕСКАЯ УСТОЙЧИВОСТЬ АНИЗОТРОПНЫХ ОВОЛОЧЕК |
357 |
ляющие явления, связанные с учетом поперечных сдвигов, в ча стности зависят от формы волнообразования (т. е.от чисел т и п )
срединной поверхности панели при потери устойчивости. С уве личением чисел та=1, 2, 3,. . . и п = 1, 2, 3,. . . влияние членов,
представляющих поперечные сдвиги, увеличивается.
Далее, следует указать, что с увеличением параметра h*
увеличивается расхождение между значениями критических сил, найденными по классической теории (k * = 0 ) и по уточненной тео
рии. Причем, как и следовало ожидать, критическая сила, най денная по уточненной теории, всегда меньше соответствующей критической силы, опрделенной по классической теории.
Наконец, укажем, что, минимизируя р*пп по т и п, из (2.6), (2.8) и (2.9) получим значения т и п, при которых р*тп получат
свои минимальные значения, т. е. найдем верхние критические значения напряжения, при котором оболочка теряет устойчивость, 2. Две задачи устойчивости замкнутой трансверсально изо тропной цилиндрической оболочки. Рассмотрим замкнутую круго вую цилиндрическую оболочку, изготовленную из трансверсально изотропного материала так, что в каждой точке оболочки плоскость изотропии параллельна срединной поверхности оболочки. Система
координат выбрана так, что i?1=co, |
R 2= R , |
А ~ 1, 5 =1 . Пусть, |
далее, |
рассматриваемая оболочка шарнирно |
оперта по |
торцам |
(а = 0, |
а=1). |
|
|
|
|
а. О б о л о ч к а н а х о д и т с я п о д с о в м е с т н ы м |
д е й с т в и е м о с е в о г о |
и |
р а в н о м е р н о |
р а с |
п р е д е л е н н о г о п о п е р е ч н о г о д а в л е н и й . Н а г р у з к а , п р и х о д я щ а я с я на е д и н и ц у д л и н ы д у г и п о п е р е ч н о г о с е ч е н и я о б о л о ч к и , р ав н а р, а и н т е н с и в н о с т ь р а в н о м е р н о г о п о п е р е ч
н о г о д а в л е н и я |
д. |
Для начального безмоментного состоя |
ния имеем |
|
|
|
Т\ = |
— р, T% = — gR, S° = 0, |
(2.11) |
тогда для Z* получим |
|
— p d?w |
|
Z* = |
( 2. 12) |
|
|
da2 |
|
Следовательно, для решения задачи статической устойчиво сти рассматриваемой оболочки, согласно (1.8.23), получим сле дующую систему дифференциальных уравнений:
Д2|,_____ V |
^ |
I |
1 |
d3w |
* |
|
R |
доз Т |
Л |
дадр2 |
|
2 |
+ ч |
d3w |
1 |
<?3ц> |
(2.13) |
|
Л |
да^д$ |
~ Л " Тр* ’ |
|
|
DA*w + § - (1 - Л *Д )^ + (1 - h * A )l* (p ^
358 КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК [ГЛ. III
где
D : |
Eh3 |
|
Eh2 |
|
Д _ |
02 |
(2.14) |
12(1 — м2)» |
**= ■ 10 (1 — v2) G' * — 0а2"1"<$2* |
Решение системы (2.13) |
представим следующим образом: |
|
|
и = |
. |
тка . |
иЙ |
, |
|
|
|
A |
cos — |
sin |
|
|
|
|
V = |
г, |
. mm |
cos |
ий |
|
(2.15) |
|
В |
sin —— |
, |
|
|
|
10= |
г, |
. «пса . |
иЗ |
, |
|
|
|
t |
Sin —— Sin |
-£■ |
|
|
где А , В , |
С — некоторые постоянные, I — длина |
оболочки, т, |
п — целые |
числа, |
характеризующие волнообразование |
средин |
ной поверхности |
оболочки. |
|
|
|
Принимая (2.15), удовлетворяем условиям шарнирного опира- |
ния по торцам. |
|
|
|
|
Подставляя значения и (а , (3), |
v (а, В), w ( а, |
В) |
из (2.15) |
в систему |
уравнений (2.13), после |
некоторых преобразований |
получим следующую однородную систему алгебраических урав
нений относительно неизвестных постоянных А , |
В, С: |
А ( X 2 |
п2) = |
С Х (v X 2 — п 2) , |
|
|
В (X 2 + |
и2) = |
— С [(2 + |
v ) X 2 + |
п21п, |
(2.16) |
(Х2 + и2)2 |
( |
Х+ |
рХ2 + |
|
9Д п 2 _ ^ _ ^ |
С \С21 + h* (X2+ |
И«) |
1 |
(Х2-f « 2)2 |
E h |
|
где |
|
|
hi |
|
|
|
m n R |
|
|
|
E № |
|
Cz = 1 2 |
( 1 |
— v * ) Д 2 * |
h* = 10(1 — v2) RiG " (2.17) |
Очевидно, что тривиальное решение системы |
(2.16) А = В = |
—С= 0 нас не будет интересовать. |
|
|
|
Из условия существования нетривиального решения, из треть его уравнения однородной системы (2.16), получим соотношение для определения критической комбинации нагрузок р и q:
рХ2-f-q R n i |
_ .9. (X2-Н n i)i |
,____ X* |
(2 18) |
E h |
~ |
1 -f h* (Х2 + «2) |
(Х2 -f « 2)2 |
|
В случае классической теории |
(h *= 0) формула (2.18) |
пере |
пишется следующим образом: |
|
|
|
рХ2-|-дД«2 |
о /*, tt I |
«КО I |
Х4 |
|
E h |
|
: С2(X2 - f П2)2 ■ |
(Х2 -f. «2)2 • |
(2.19) |
Уравнение (2.18) отличается от соответствующего уравнения классической теории (2.19) лишь наличием члена с множителем/г*,
§ 2] СТАТИ ЧЕСКАЯ УСТОЙЧИВОСТЬ АН И ЗО ТРО П Н Ы Х ОБО ЛО ЧЕК |
359 |
т. е. члена, учитывающего явления поперечных сдвигов. Поэтому дальнейшее исследование уравнения (2.18) ничем не отличается от анализа классического уравнения (2.19).
Не вдаваясь в подробности, укажем, что с помощью уравнения (2.18) нетрудно определить те комбинации чисел т и п , при ко торых нагрузки р и q получат свои критические комбинации в слу
чае учета поперечных сдвигов. |
|
б. О б о л о ч к а н а х о д и т с я п о д |
д е й с т в и е м |
к р у т я щ и х м о м е н т о в М 0, п р и л о ж е н н ы х п о
к о н ц а м о б о л о ч к и . |
Такая |
нагрузка в оболочке вызывает |
следующие внутренние силы безмоментного состояния: |
|
|
Z* = 0, |
Г» = 0, S° = MJ2*Д2. |
|
(2.20) |
тогда для Z* получим |
MQ d2w |
|
|
|
|
|
(2. 21) |
|
z*= ui?2 dctdp • |
|
В этом случае первые два уравнения системы |
(1.8.23) |
остаются |
неизменными, |
а третье уравнение |
принимает |
следующий вид: |
ДД*ш + |
г * , . |
d *w |
|
|
|
§ .(1 |
- 5 - ( 1 - ^ д) д гт а - = ° - |
<2-22> |
|
|
Граничные условия свободного опирания удовлетворяются |
приближенно, |
а именно: |
|
|
|
|
при а = 0, |
a= l |
W — 0, |
|
|
(2.23) |
т. е. удовлетворяются лишь условия, накладываемые на нормаль ное перемещение.
Такая форма удовлетворения граничных условий позволяет решить задачу с помощью лишь уравнения (2.22), без учета пер вых двух уравнений системы (1.8.23).
Решение уравнения (2.22), удовлетворяющее условиям (2.23), можно представить в виде
w = ^ C { sin |
(2.24) |
»=i |
|
тогда условия (2.23) сводятся к следующим уравнениям: |
<?1 + <?2 = |
0, |
Ci sin (r*x - п | - )+ С2sin ( ц2£- — п |-) = |
(2.25) |
0. |
Условия (2.25) имеют место при любом р и отличных от нуля значениях С{, лишь если
360 КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК [ГЛ. Ill
т. |
е. когда |
2mnR |
|
|
|
|
|
1*2 1*1— |
(т = |
1, |
2, ...). |
(2.26) |
|
I |
|
|
|
|
|
Подставляя (2.24) в (2.22), получим |
|
|
|
2 |
с -( щ г ^ г ,яж ^ + ^ rt [i + £ |
W ■+" * ) ] - |
|
|
М, |
|
|
|
|
|
|
s - [ i |
|
|
|
" Ю = ° - |
< 2 -2 7 > |
Отсюда для определения критического значения крутящего мо |
мента получим следующее выражение: |
|
|
|
|
|
((4 + |
п2)2 |
|
IА |
(2.28) |
|
«ft- [1 + h* (|4 + пЩ |
П(ft2+ » 2)2!• |
|
|
где с* и h* определяются по формулам (2.17). Здесь рД£=1, 2) — |
некоторые параметры, которые удовлетворяют условию |
(2.26). |
Кроме того, параметры ft должны быть такими, чтобы выражение |
(2.28) для критического значения крутящего момента имело одно |
и то же значение при £=1 |
и £—2. |
|
|
|
|
На основании сказанного выше нетрудно каким-либо извест |
ным способом определить критическое |
значение крутящего мо |
мента с учетом явлений, связанных с поперечным сдвигом.
3. Об устойчивости трансверсально изотропной сферической оболочки. Пусть замкнутая сферическая оболочка находится под действием равномерно распределенного по поверхности нормаль ного внешнего давления g=const.
Под действием равномерно распределенного нормального дав ления q в оболочке в безмоментном состоянии появляются следую
щие внутренние силы: |
|
|
= |
Т\ = — Ц -, S0 = 0. |
(2.29) |
Тогда в силу (1.8.26) и (2.29) из (2.1) получим |
|
Z* = |
- £ - ( * + 2 ) w, |
(2.30) |
где А — безразмерный (А = В , B = R sin а) оператор Лапласа,
который в принятых здесь географических координатах имеет вид
Д = ж т й ( 81па1 ) + Ж т й -