книги из ГПНТБ / Рикардс Р.Б. Устойчивость оболочек из композитных материалов
.pdf4.2. |
Устойчивость под действием |
внешнего давления |
123; |
||||
оболочки |
h/R —1/10, 1/25 |
|
|
||||
и 1/100; |
график построен |
|
|
||||
для |
|
величины |
RjL= 1. |
|
|
||
Параметр h-*~1, когда мо |
|
|
|||||
дуль межслойного сдвига |
|
|
|||||
стремится к бесконечнос |
|
|
|||||
ти. Для сравнительно тон |
|
|
|||||
ких оболочек (кривая 3) |
|
|
|||||
влияние поперечных сдви |
|
|
|||||
гов при реальных величи |
|
|
|||||
нах |
модуля |
межслойного |
|
|
|||
сдвига |
мало |
ощутимо. В |
|
|
|||
то же время, у сравни |
|
|
|||||
тельно |
толстых |
оболочек |
|
|
|||
(кривые 1 и 2 ) |
это влия |
Рис. 4.2.2. Кривые изменения |
параметра k~ |
||||
ние |
значительно. В слу |
в зависимости от величины модуля меж |
|||||
чае |
внешнего |
давления |
слойного сдвига. |
|
|||
на |
критическое состояние |
|
|
||||
конструкции и на форму потери устойчивости влияет относи-- тельная длина оболочки. На рис. 4.2.3 представлен график из менения параметра k в зависимости от длины оболочки. Кривые- 1—4 соответствуют параметрам | 4 = ^5 = 0,1; 0,05; 0,02 и 0,01; для всех кривых принято h/R= 1/10. С уменьшением длины увеличи вается разница между результатами, полученными на основе классической теории и найденными с помощью расчетов по фор--
муле (2. 15).
Из рис. 4.2.3 видно, что у коротких оболочек, например,, когда L/R = 0,5, погрешность, вносимая при расчете с использо ванием гипотезы Кирхгофа—Лява, даже в случае большой;
Рис. 4.2.3. Зависимость критического внешнего давления от длины оболочки. Обозначения те же, что на.
рис. 4.2.1.
Глава IV . Устойчивость оболочэк при статических нагрузках |
124 |
■сдвиговой жесткости (кривая ]), значительна. В случае длинных оболочек расчет критических нагрузок по классической теории и теории типа Тимошенко дает близкие результаты. У коротких оболочек большее влияние поперечных сдвигов на критическую нагрузку можно объяснить тем, что с уменьшением длины увели чивается число волн, по которым' выпучивается оболочка, т. е. локальные кривизны больше как в поперечном, так и в продоль ном направлениях. Следовательно, при выпучивании коротких оболочек деформации поперечных сдвигов больше, чем длин ных. Эксперимент, проведенный на оболочках из стеклопластика, показал, что под действием внешнего давления критическая на грузка падает с уменьшением длины по сравнению с теорети ческим расчетом по классической теории. Отметим также, что с уменьшением модуля межслойного сдвига изменяется форма потери устойчивости. Так, например, для оболочки с парамет рами h/R = 1/25; R/L = 1; £4 = £5 = оо; 0,05; 0,01 число волн по ок ружности, при котором значение критического давления (2.15) имеет минимум, /г*= 5,7 и 8 соответственно. С уменьшением мо дуля межслойного сдвига увеличивается число волн, по которым
стремится |
выпучиваться оболочка. Для иллюстрации на |
рис. 4.2.4 |
приведены кривые, которые характеризуют зависи |
мость внешнего давления от параметра волнообразования. Гра фик построен для h/R= 1/25 и R/L= 1; кривые 1—5 соответст вуют параметрам £4 = | 5 = оо, 0,3; 0,1; 0,05 и 0,01. Параметр волнообразования может принимать только дискретные целочис ленные значения вблизи точки минимума (на рисунке обозна
чены кружками). |
устремив £4 = ^5->-0, можно получить |
|||||||||
Отметим, что из (2. 15), |
||||||||||
формулу |
для |
определения критического |
усилия, когда модуль |
|||||||
|
|
|
межслойного сдвига стремится |
|||||||
|
|
|
к нулю. Интересно отметить, |
|||||||
|
|
|
что |
критическое |
давление |
в |
||||
|
|
|
этом |
случае |
не равно нулю, а |
|||||
|
|
|
составляет |
|
примерно |
5—10% |
||||
|
|
|
от |
величины |
классического |
|||||
|
|
|
давления. При этом оболочка |
|||||||
|
|
|
стремится |
выпучиваться с об |
||||||
|
|
|
разованием |
большого |
количе |
|||||
|
|
|
ства волн по окружности. |
|
||||||
|
|
|
Под нагрузкой оболочка из |
|||||||
|
|
|
стеклопластика до потери ус |
|||||||
|
|
|
тойчивости |
|
может |
расслаи |
||||
|
|
|
ваться |
[66], |
следовательно, |
в |
||||
|
|
|
некоторых зонах модуль меж |
|||||||
разования. |
1—5 |
соответствуют |
слойного |
сдвига оказывается |
||||||
54=15=°°; |
0,3; 0,1; 0,05; 0,01. |
равным нулю и нагрузка, кото |
||||||||
4.3. Геометрически нелинейная задача |
12S |
рую может принять оболочка, падает. Эксперименты [66] указы вают, что у сравнительно толстых стеклопластиковых оболочек потеря устойчивости происходит после расслоения или растрес кивания (к этому вопросу вернемся в следующем параграфе при исследовании закритических деформаций оболочки).
Из приведенных расчетов можно сделать вывод о пределах применимости гипотезы Кирхгофа—Лява при расчете на мгно венную критическую нагрузку цилиндрической оболочки под
внешним давлением. При относительной |
толщине |
h / R ^ 0,04, |
||
длине |
L/R^ . 2 и |
при величине параметра |
модуля межслойного |
|
сдвига |
^4 = ^5 —0,1 |
расчет следует вести |
на базе |
уточненных |
теорий. |
|
|
|
|
4.3. ГЕОМЕТРИЧЕСКИ НЕЛИНЕЙНАЯ ЗАДАЧА
ПРИ ВНЕШНЕМ ДАВЛЕНИИ НА ЦИЛИНДРИЧЕСКУЮ ОБОЛОЧКУ
Известно [44], что более полную картину потери устойчивости и послекритического поведения дает решение в геометрически нелинейной постановке. В экспериментальной работе [66], посвя щенной исследованию стеклопластиковых оболочек, установ лено, что после потери устойчивости у оболочек с параметром R/h = 20—150 наблюдается сдвиг слоев. Это свидетельствует о превышении предела прочности на срез. Классическая теория оболочек (гипотеза Кирхгофа—Лява) и уточненная геометри чески линейная теория не дают возможности достаточно убеди тельно объяснить этот факт.
Рассмотрим задачу об устойчивости в «большом» цилин дрической ортотропной стеклопластиковой оболочки под дейст вием внешнего давления*.
Оси х, у, z направим вдоль образующей, по окружности и нормали. Примем, что оси ортотропии материала совпадают с главными направлениями оболочки. Считаем, что оболочка имеет начальные неправильности в виде геометрических откло нений срединной поверхности от идеальной формы. Используем систему уравнений устойчивости (3.5.25) для оболочек боль шого прогиба с начальными несовершенствами, заменив упруго вязкие операторы упругими:
(Vy—Yv(0). Y*~Y*(0)} |
+ |
|
+ НАыЦ\2{уу — yv{0), |
®—“>«»} =0; |
(4-3. 1) |
* Подобная задача с использованием классической теории решена в
{31, 225 и др.].
Глава IV . Устойчивость оболочек при статических нагрузках |
126 |
У= 4 ;1I{’YX—Yx(0), Yy —Yy(0)} + |
(4.3.2) |
|
+ ^55'Пп {Yx —y.t<0). o' —o>(o)} =0; |
||
z ^ X (Y x - Yx(0), YyYy(0)} - |
|
|
|
|
(4.3.3) |
1 |
d2 { w - w (0)) _ |
|
V=V*F+ — |
dx2 |
|
HR |
(4.3.4) |
|
- k {w , w} +Х{Ш(0)>аУ(о)} =0. |
||
Граничным условиям шарнирного опирания оболочки удовлет ворим, представив неизвестные функции w, ух и уу, следующим образом:
жении прогиба w учтено, что вмятины к центру кривизны имеют большую амплитуду. Экспериментальные исследования по фор ме выпучивания, приведенные в главе VI, свидетельствуют о том, что, действительно, волны, по которым происходит потеря устойчивости, имеют большую амплитуду к центру кривизны.
Принимаем, что форма начального волнообразования нахо дится «в резонансе» с волнообразованием оболочки в процессе деформации. Такое допущение несколько увеличивает влияние начальных неправильностей в сторону запаса устойчивости [44], но значительно облегчает решение. Тогда для начальных проги бов и углов поворота нормального волокна
- |
w0 |
|
^ |
w0= —- = *i0sinoasin By+ х2Эsin2 ах; |
. |
||
|
h |
... |
|
•у**0*= *3ocos ax sin $y + xw sin 2oa; |
|
|
|
Yv(0) = x50 sin <xxcos $y. |
|
(4.3.6) |
|
В главе VI приведены результаты экспериментальных изме* рений формы начальных несовершенств цилиндрических стекло пластиковых оболочек и ее развитие при нагружении и во вре мени. Показано, что не всегда форма выпучивания ховпадает с формой начальных несовершенств. г ,
.4-3. Геометрически нелинейная задача |
127 |
Подставим (3.5) и (3.6) в уравнение (3.4). После интегри рования этого уравнения получим выражение для функции на пряжений, в граничных условиях которой учтено, что на обо лочку действует внешнее давление (напряжение 022° = —qR/h):
F=- |
=[ |
х102) ~ |
h2s |
--^гЬрл-Ярг^г' 1 (х2 — х20) ]cos2ax + -r^-npiX
Хп -^ ц ^ !2 —Л'ю2) cos 2 ^y+ — (xlx2 - x l0x20) х
У/2j
Xsin3axsin &у-{—^-[«“2р2- 1(*1 —х10) — |
|
, k 2 |
|
dR |
(4.3.7) |
- (x1x2- x Io^2o)]sin ах sin fiy— ^ j - x 2. |
Здесь безразмерные параметры материала и геометрии оболочки определены по (2. 12), а коэффициенты ki и k2 следующим спо собом:
^1 = 9ц % 2 {n~2K2pi2 --2 U l{-%2- x\xrx+ U -{+ ^ n 2n~2p r 2l r lP~l :
k 2 = ц~1\2~1п~2я2рI2—2 |3£1~112-1 + |б~1+ п2л~2р\-21 г 1ц_!•
Далее решаем систему (3.1) —(3.3) методом Бубнова—Га- леркина, подставив в правые части этих уравнений разложения (3.5) и (3.6), а также выражение (3.7) для функции напряже ний. Согласно этому методу, интегрированием по поверхности оболочки необходимо найти следующие уравнения:
2ЯЯ L
I I Z sin ах sin fiydxdy = 0;
2яН Ь
J J Z sin2 axdxdy = 0 \
оо
2ЯЯ L
II Y cos ах sin fiydxdy = 0;
оо
2ЯЯ L |
2яЯ |
L |
|
||
J |
J X sin ах cos $ydxdy = 0. |
I |
S |
У sin 2аxdxdy = 0; |
(4. 3. 8) |
О |
0 |
0 |
0 |
|
|
Глава IV . Устойчивость оболочек при статических нагрузках |
12а |
После интегрирования получаем систему, состоящую из двух, нелинейных и трех линейных алгебраических уравнений относи тельно пяти неизвестных х4 (i —1, 2,... , 5):
У 1 — a i-'6 3 + О 2 Х 1 Х 2 2 — 03X 1X 2 + Q4X1 + 05X2 + 03X3 + О7Х5 + |
Og — 0; |
У 2 — Ь \Х \ 2Х 2 — Ь 2Х 12 - \ - Ь з Х { + 6 4 X 2 + ^5X4 + Ь 6 = 0 ] |
|
Уз— С\Х\ — с2х3- с3х5 + сА= 0; |
|
У*= d[X2 — d2Xj + d3 = 0; |
|
г/s = / 1-^1 — / 2-^3 — / 3-^5 + / 4 = 0 . |
( 4 . 3 . 9 ) |
Здесь коэффициенты a, b, с, d, f выражаются через геометричес кие и физические параметры материала, величину нагрузки, ам плитуд начальных прогибов и углы поворота нормального во локна. Представим коэффициенты а* для первого уравнения системы:
Oi=—-l (Zin2pl2n - 2 +%2n - 2p r 2n2);
“2“ |
х |
ж |
+ ~ й г ) ; |
|
|
|
«з=4- Ь ^ _2р г 2Р2_! + 4 - «~2р2“'; |
|
|||||
|
О |
|
|
/22 |
|
|
« 4 = |
0 , |
• 0 |
- 4 р 2 - 2 ------— q n ~ 2p |
r 2 - a i X w 2 -\— — ^ 2 l t _ 2 p + 2p 2 _ I ^ 2 0 ; |
|
|
|
Z&2 |
|
|
£ |
О |
|
0,5 = |
1 |
|
-2 |
-1 |
|
|
" о и ~ П |
2Р 2 |
Хх\0 ~ а 2х \Ъх 2Ъ, |
|
|||
|
1 к<1 |
|
|
|
|
|
« 6 = -2 ^ " [liJtp i« _2P2_ 1+ (|з + |
2 5 б ) я -,р г 1р2- 1 ]; |
|
||||
|
|
[Е2П- 2р Г 2Р2- 1 « + ( | з + 2|б) « “ ’рг- 1 ]; |
|
|||
« 8= |
-« б Х з 0 -О 7Х50 + - 4 — (M- 2p2_ 1X l0 X 20 -n _4P2_2JClo). |
(4. 3. 10а) |
||||
|
|
|
|
2к2 |
|
|
Таким же образом получены коэффициенты для второго уравне ния системы:
4.3. Геометрически нелинейная задача |
129 |
Ьз= - ~ п - 2р2 - 1х 1й\ J.K2
Ь^=~%2П-2р Г 2П~2^2~2\
Ьь= - ^ 1 \Щ\Р2~1П-2\
6б= -у^-|2Л'-2р1~2р2-1-Х:К>2 — ^5 X4 3 — Ь^Х2 0 - |
(4.3. 106) |
Коэффициенты di и fi содержат параметры межслойной сдви говой жесткости:
С1=|5Яр1Р2, |
С2 ——^(^,\Л2р12р22 + ^6р22П2) +£5; |
|
||
С з ~ ~ ^ 2 (^3+ |
1б) Л Р 1Р 2 2ГГ, с 4 = С2Х 30 + |
С3Х50 — С\ХГю; |
|
|
^1 = Е 5Я р1р2, |
^ 2 = _T 'E l I1t2p i 2p22 + |
i5 , |
d 3 = d 2X^Q |
й [Х 2о1 |
|
о |
|
|
|
/l=^4«p2, |
/2=-у2'(|з + 5б)яр1Цр22; |
|
|
|
[з= _^2'(^бЯ2р12р22+ ^2р22я2) + |4; |
fi = f2^30“Ь/3^50 |
fiMo* (4. 3. 10в) |
||
В выражениях (3. 10а) — (3. 10в) безразмерные параметры попрежнему определяются формулами (2. 12).
Решение системы (3. 9) можно найти итерационным методом Ньютона—Рафсона на ЭЦВМ, для чего необходимо задавать начальные приближения к решению и использовать способ «по следовательного нагружения» [145], т. е. создать картину, подоб ную нагружению в реальной конструкции.
На рис. 4.3.1 сплошной линией показано нагружение в мо дели, прерывистой — нагружение в реальной оболочке. Сначала находим верхнюю критическую нагрузку qB по формуле (2. 15), затем шагом, например 0,01 qB, изменяем нагрузку, причем за первое начальное приближение неизвестных Xi принимаем их значения при q = 0, т. е. величины начальных амплитуд прогибов и углов поворота нормального волокна. В итоге мы получаем
9 - 1744
Глава IV . Устойчивость оболочек при статических нагрузках |
130 |
Рис. 4.3.1. К определению верхней предель ной и нижней критических нагрузок.
кривую нагрузка—амплитуда прогиба. Сходимость при решении системы (3.9) была хорошей только в докритической стадии. Для закритической стадии (после «хлопка» в модели) по лучить кривую нагрузка—амплитуда прогиба не удалось. Для улучшения сходимости система (3.9), состоящая из пяти урав нений, была преобразована в систему двух уравнений путем ис ключения из первого и второго уравнений системы (3.9) неиз вестных х3, х4 и х5. Из четвертого уравнения системы получаем выражения для неизвестного
х4= |
d\X2 + d3 |
(4.3.11) |
|
|
d2 |
Решая совместно третье и пятое уравнения системы (3.9), полу чаем выражения для неизвестных х3 и х5:
/4С3 —сДз_______Cjf3— C3fi ^
f2c3 — c2f3 / 2 ^ 0 £2 / з
^2/4—£4/2_____Cif2 — fjC2 ^ |
(4.3. |
12) |
||
Czfz — C'Ji |
c2f3 — c3f2 |
|||
|
|
|||
Подставляя (3.11) и (3.12) в первое и второе уравнения сис темы (3.9), получаем систему двух нелинейных алгебраических уравнений относительно двух неизвестных Х\ и х2, сумма кото рых f= x i+ x 2 представляет амплитуду прогибов:
4.3. Геометрически нелинейная задача |
131 |
у 1—а\Х1ъ + а2Хух22 — агх хХ2+ у а 4—
|
а б(с \ ! з ~ сзЫ -\-a 7 {c \ h ~ flC2) |
Ху+ |
|
||
|
/гсз — сз}з |
|
|
|
|
|
|
a7{c2fi~cif2) - = |
|
||
+ |
+ G8 +' аб((4с3 — с4з) |
0; |
|||
|
|
сз!2~с2!з |
|
|
|
|
|
M l |
|
— ~}— — |
0. (4.3.13) |
У2 — ЬуХу2Х2— Ь2Ху2-}-ЬууХу-р (ь^-\------ ----+ |
|||||
|
' |
«о |
' |
do |
|
Сходимость итерационного процесса для системы |
(3. 13) значи |
||||
тельно лучше, чем для (3.9). Используя предложенный способ «последовательных нагружений» в сочетании с методом Нью тона—Рафсона, мы получили кривые нагрузка—амплитуда про гиба для закритического состояния. Характерная диаграмма на грузка—амплитуда прогиба представлена на рис. 4.3.1. На грузку, при которой в модели происходит «хлопок», будем счи
тать предельной |
(точка А на рис. 4.3. 1), |
а нагрузку, при кото |
рой имеет место |
«выхлоп», примем за |
нижнюю критическую |
(точка С на рис. |
4.3.1). Величины qn и qH полностью определя |
|
ются геометрическими и физическими параметрами, а также за висят от начальных несовершенств оболочки, обусловленных технологией изготовления.
Графики на рис. 4.3.1 и 4.3.2 построены при следующих значениях параметров: | 1= | 2=1; £з= £б = 0,15; g4 = E5 = 0,l; р! = = 0,55; р2 = 0,04; п* = 4. На рис. 4.3.2 принята амплитуда началь ного прогиба /o= ^io + -^2o = 0,l; 0,5; 0,9 для кривых 1—3 соответ-
U_1.11_______________ 1 |
J |
||||
0 |
I |
2 |
3 |
4 |
5 |
Рис. 4.3.2. Характерные диаграммы нагрузка—амплитуда прогиба для различных значений начальной ам плитуды. 1—3 соответствуют ампли туде начального прогиба f0=0,l; 0,5; 0,9.
Рис. 4.3.3. Развертка поверх ности цилиндрической оболоч ки (заштрихованы зоны пов реждения).
9*
Глава IV . Устойчивость оболочек при статических нагрузках |
132 |
ственно. В вычислениях принято допущение (2.1.30). |
В этом |
случае величины х30, х40 и х50 выражаются через амплитуду на чальных прогибов Хю и х20 следующим образом:
-Гзо = х40 = лрфг^го, -*5o = ?lP2^io- (4.3.14)
Для приведенного примера значение верхней критической на грузки qB~ 8,2210~2, значение верхней предельной и нижней критической нагрузок для кривой 1 <р= 7,85-10~2 и ^н = 6,98-10-2 соответственно. Из рис. 4.3.2 видно, что в случае больших на чальных искривлений с увеличением нагрузки амплитуда моно тонно возрастает, т. е. отсутствует «хлопок». Разница между верхней и верхней предельной критической нагрузками неболь шая (меньше, чем у изотропной оболочки с такими же парамет рами толщины и длины). Этот факт установлен эксперимен тально и исследован теоретически [66, 205]. Показано, что для анизотропных оболочек (например, из боропластика) влияние начальных неправильностей в виде отклонений срединной по верхности от идеальной формы играет меньшую роль, чем у изо тропных. Тем не менее, в оболочках из армированных пластиков существенную роль может играть неоднородность материала. Это обстоятельство можно учитывать, если отнести случайную неоднородность к отклонению поверхности оболочки от идеаль ной формы и тем самым приближенно оценивать возможность снижения критической нагрузки.
Ранее было отмечено, что в экспериментах на оболочках из стеклопластика после потери устойчивости материал может рас слаиваться. Для толстостенных оболочек расслаивание воз можно и до потери устойчивости или в процессе потери. Макси мальное сдвиговое напряжение, например огз, в оболочке после потери устойчивости (точка В на рис. 4.3. 1) можем определить следующим образом. Подставляя в соотношение для деформа ций поперечных сдвигов
d(w — wо)
823= - |
( У у ~ У у Ф ) ) + |
ду |
|
|
|
|
|
выражения для уу, уут и w, w0 (3.5) и (3.6), |
получаем |
||
823 = [ - |
(*5-*5о) + (X i-x10)ttp2]sinaxcos $У- |
(4-3- 15) |
|
Соответствующие сдвиговые напряжения огз определяются по
закону |
Гука 023= 823623 |
(здесь 6 2з — |
модуль межслойного |
сдвига). |
Максимальными |
напряжениями |
огз характеризуются |
точки оболочки, где sin ax cos рг/= 1, что соответствует х — Ц 2 и
y = nR/m |
(т= 1 ,2 ,...,/г*). Примем А ц = А 22 = 2,5-105 кгс/см2. |
Тогда для |
нашего примера максимальные сдвиговые напряжения |
