книги / Прикладные задачи устойчивости многослойных композитных оболочек
..pdf7.3. Расчёт критического внешнего давления |
191 |
Обычно при потере устойчивости от внешнего давления выполня ется условие гг2 1. В этом случае из соотношения (7.15) после ми нимизации по параметру п получаются следующие расчётною формулв1 для определения критического колвцевого усилия:
fjnKp _ rjnnc |
}не) |
(7.16) |
К 2 - Т Т = 7 ^ Г |
\ / в Л О Г 3 |
ip.
Коэффициент у в соответствии с [8] вычисляется по формуле:
7 = 3 —2,1-г\ |
О < — < 0,8. |
(7.17) |
|
|
Г2 |
Г2 |
|
При малой конусности |
> 0,8^ |
полагают 7 = |
1 . |
|
у г2 |
|
|
Полученные расчётною формулв1 применимы, если выполняется
условие |
|
|
1,75тг |
А/= |
(7.18) |
К 2 < 0,5 Т2Ж; Т2Ж = 7 |
В ф \ . |
крЛ-х.
Если заполнителе достаточно жёсткий на поперечною сдвиги: К 2 > > 0,5 Тж, то можно использоватв модель прямолинейного элемента (типа сдвиговой модели Тимошенко). Тогда разрешающее соотношение
имеет вид |
п |
В 2 |
В\ Af |
|
|||
Т2 = |
,2 Д2 |
гР2m ax |
(7.19) |
1 |
+ и>рп2 |
Пи |
После минимизации этого выражения по параметру п можно по лучить следующую формулу для определения критического кольцевого усилия [72]:
грКр _ грж ' |
Q |
Т * \ |
|
' - |
l e |
t - |
, 7а д |
Зная значения критических усилий, можно найти критическое внешнее давление. Для оболочек со слабым на поперечные сдвиги заполнителем (К2 < 0,5 Тж):
|
рКР=рпс + _ К ^ . |
рнс = 7 1 |
^ . |
4/^1 { D r ) 3 |
(721) |
|||
|
|
|
Рmax |
I p ^ax |
|
|
||
В случае достаточно жёсткого на поперечные сдвиги заполнителя |
||||||||
( к 2 > \ |
Т» критическое внешнее давление определяется формулой |
|||||||
ркр = р» |
, _ л т ? |
|
1,75л |
(7.22) |
||||
|
|
B XDI |
||||||
|
|
|
16 |
к2 |
^ |
£Ртах |
|
|
|
|
|
|
|
е 3/2 |
|
|
|
При |
действии |
внешнего |
давления |
в |
расчётные |
теоретиче |
||
ские формулы |
вводится |
поправочный |
эмпирический коэффициент |
|||||
куст = 0,8 -7 0,9 |
[23, |
27]. |
|
|
|
|
|
192Гл. 7. Устойчивость многослойных композитных конических оболочек
7.4.Устойчивость при действии крутящего момента
исдвигающих усилий
При исследовании устойчивости от действия крутящего момента Мкр или сдвиговое усилий S ограничимся слабоконическими пологи ми трёхслойными оболочками. Тогда можно восполвзоваться зависимо стями для цилиндрических оболочек. Особенноствю потери устойчиво сти подобных оболочек является то, что в случае слабв1х на поперечные сдвиги заполнителей при потере устойчивости образуется много волн как в кольцевом, так и в осевом направлениях [66]. В то же время до статочно жёсткие на поперечные сдвиги оболочки теряют устойчивость по классической схеме, образуя одну полуволну в осевом направлении. Поэтому необходимо рассмотреть отдельно каждый из этих типов потери устойчивости.
В случае трёхслойных оболочек со слабым на поперечные сдви ги заполнителем по аналогии с цилиндрическими оболочками можно получить следующие зависимости для расчёта критических усилий сдвига [66]:
S KP= y / ( K i + T - p ) К >; г кр k opTrI Z ^ j 1 - ^
ГГ = - ^ |
J B 2 [)»>■-, |
т г = — |
A/ S A ; |
ртах = ^ т т . (7.23) |
Р т а х |
V |
Р т а х |
V |
COS 0 |
Если трёхслойные или многослойные оболочки имеют высокую жёсткость на поперечные сдвиги, то по аналогии с цилиндрическими оболочками для расчёта критических усилий сдвига можно получить
[ 66] : |
|
rp:ж |
|
|
|
|
|
|
,S'Kp |
1 - 0 ,3 4 ^ - |
|
|
|
л 2 |
|
S УК |
3,3 |
1,75тг |
(7.24) |
|
|||
|
|
^ Р т а х |
|
Расчётные зависимости (7.23) для оболочек с малыми жесткостями
на поперечные сдвиги применяются, когда выполняется условие |
|
лГ кЖ 2 < 0,33 5 Ж. |
(7.25) |
Расчётные зависимости (7.24) для оболочек с большой жёсткостью
на поперечные сдвиги применяются, когда выполняется условие |
|
||
/ ^ г ^ > 0 , 3 |
3 5 ж. |
(7.26) |
|
В пограничной области |
~ |
0,33 £ ж) следует провести рас |
чёты по формулам (7.23), (7.24) и выбрать минимальный результат.
7.4. Действие крутящего момента и сдвигающих усилий |
193 |
|||||||
Критические |
значения крутящего |
момента М кр и эквивалентной |
||||||
перерезывающей силы QKp определяются зависимостями: |
|
|||||||
y/Ki K2 < 0,33 ^ |
: |
|
|
|
|
|
|
|
|
№> = 2np2mi^ /( K i + T-P) K2; |
|
|
|||||
|
QKp = npminV(Ki + T r ) K2 ; |
|
(7.27) |
|||||
y/Ki K2 > 0,33 ^ |
: |
|
|
|
|
|
|
|
|
П^Ж.\ |
|
|
6,6п pmmn |
VbD 5 ; |
|
||
М кр = М ж |
1 - 0,34 Т2 |
\ ; М ж= |
|
|||||
|
K2) ; |
|
|
71/2 |
|
|
||
QKp |
T2?\ |
; |
Q |
= |
3,3п pm/4 |
^ в D . |
(7.28) |
|
K2 |
; |
|||||||
|
|
|
71/2 |
|
|
Соотношения (7.27), (7.28) получены с помощью формул (7.23), (7.24) и с учетом зависимостей
М = 2nR2S ; Q = nRS.
7 С.Н. Сухинин
Г л а в а 8
УСТОЙЧИВОСТЬ многослойных
СФЕРИЧЕСКИХ ОБОЛОЧЕК ПРИ ДЕЙСТВИИ ВНЕШНЕГО ДАВЛЕНИЯ
8.1. Разрешающие соотношения
Рассмотрим устойчивость сферического сегмента с углом полураствора $. Значению д = 7г/2 соответствует полусфера, значению "д = д — полная сфера. Разрешающие соотношения на основе модели лома ной линии применительно к исследованию устойчивости трёхслойных и многослойных сферических оболочек получаются из зависимостей (1.1)-т-(1.5), если положить в них вместо радиусов кривизн R \, R 2 радиус сферической оболочки R = R\ = R 2. При S = 0 получаются разрешающие соотношения для анализа устойчивости сферического сегмента в случае действия наружного давления [71]:
|
T = ( D f + |
nz?f) |
|
|
A2 |
B 2 |
1 + Ф2 |
|
|
|
R2 |
|
1+ ф2 + F2 (ф) |
A2 |
’ |
|
|||
|
|
|
|
|
|||||
F\ |
(ф) = 1 + а-1Ф2 + Ф\ФА\ |
|
F2 (ф) = 1 + сх2ф2 + /32фА\ |
|
|||||
|
Ш7Г |
|
|
|
pR |
£ 2 |
|
|
В2 |
ф = п / A; A |
T = T\ T2 |
~ Y ; A |
D \ ’ ih P |
Вф |
|||||
|
|
^ A 2 , 0 |
\ |
B 2 |
0 |
|
|
/0 N |
|
|
ai = —---- b 2 V2 |
a2 =—-----2 г/2- |
|
|
(8.1) |
||||
|
|
|
|
|
B \2 |
|
|
|
|
Здесь D™c |
— суммарная изгибная |
жёсткость |
несущих |
слоёв, |
D\ = |
= 77“с + D* — изгибная жёсткость трёхслойного пакета; D"c + П D* —
эффективная |
изгибная |
жёсткость трёхслойного пакета; В 2 = В2х |
х (1 —^1^2); |
В2, i/\, щ |
— жёсткость пакета с лёгким заполнителем |
в широтном направлении и соответствующие коэффициенты попереч ной деформации; р — наружное давление; од, а2, фи Ф2 — параметры анизотропии трёхслойной конструкции. Величина П(А,ф), определяю щая влияние поперечных сдвигов (0 < Q ф 1), имеет вид (1.5).
Критическое значение усилия Ткр и соответствующие значения Акр, фкр параметров волнообразования получаются из (8.1) минимизацией по А и ф.
Из разрешающего уравнения (8.1), построенного на основе общей модели ломаной линии, как частные случаи получаются зависимости,
8.2. Устойчивость изотропных оболочек (классическая модель) |
195 |
|||
соответствующие моделям более низкого уровня. Так, если Q.D* |
||||
3> D™, то применима модель |
прямой линии, |
и |
соотношение |
(8.1) |
примет вид: |
|
|
\+ф2 |
|
А |
А2ВД) |
_ |
|
|
1 + Л2 (ип +си2Ф2) (1 + Ф2) R2 |
2АЩ{ф)' |
|
||
А |
А |
|
|
,0 г,\ |
W1 - |
Ш2 |
|
|
(8-2) |
Здесь А> А — жёсткости ортотропной оболочки на изгибы в двух
направлениях; шь |
^ |
- |
безразмерные податливости оболочки на по |
|||
перечные |
сдвиги; |
К\, |
А |
— жёсткости |
пакета на |
поперечные сдви |
ги [41]. |
|
|
|
|
|
|
В том |
случае, |
когда |
жёсткость на |
поперечные |
сдвиги заполни |
теля велика, т. е. податливости сщ, малы и выполняется условие А2(сщ Т и>2'ф2) -С 1, справедлива классическая модель оболочек с неиз менной нормалью. Тогда из соотношений (8.2) получается зависимость для классических сферических оболочек:
= А |
х*Щ) |
1+Ф |
R 2 |
1 + А |
2АЩ{ф) |
8.2. Устойчивость изотропных оболочек (классическая модель)
Рассмотрим сначала решение задачи об устойчивости изотропных сферических оболочек, построенное на основе классической модели оболочек с неизменной нормалью. Количественный критерий приме нимости классической модели будет дан ниже. В случае изотропных оболочек соотношение (8.3) приводится к виду:
Т = - ^ х + — ; х = А2(1 + ф 2) = А2 + п2; |
|
|
Rz |
х |
|
|
В = В { 1 - г /2). |
(8.4) |
Минимизируя выражение (8.4) по параметру х, получим классиче ский результат:
TKp = P ^ R = 2 v ^ ; Л2р + п2р = А2, = |
(8.5) |
Как видно, в случае изотропных сферических оболочек форма волнообразования при потере устойчивости (параметры Акр и пкр) однозначно не определена: параметры А и п могут быть любыми, лишь бы выполнялось условие А2р + п2р = А^.. Таким образом, одному значению критической нагрузки Ткр соответствует множество форм волнообразования. Это позволяет предположить [4, 69], что результаты
7*
196 Гл. 8. Устойчивость многослойных сферических оболочек
расчёта критического значения наружного давления в общем случае могут значительно отличаться от экспериментальных данных, что и на блюдается в изотропных однородных сферических оболочках.
В случае однородных по толщине h оболочек с характеристиками упругости Е н г жёсткости определяются соотношениями
В = Е h , D = 12(1 _ v2y
а формулы для расчёта критических параметров примут вид [8, 81]:
Ркр — |
Eh2 |
1,21 Е |
|
R 2 |
|||
|
Д2’ |
||
|
/ M l - * ' 2) |
|
К Р + пкР = >L = 3,3 |
R |
(8 |
.6) |
|
8.3.Изотропные трёхслойные оболочки
слёгким заполнителем
Рассмотрим устойчивость при действии наружного давления трёх слойных изотропных сферических оболочек. В этом случае разрешаю щее соотношение (8.1) приводится к виду:
Т = |
D* |
х |
В |
|
Dnc + |
д 2 + Т : |
|
||
|
1 + U!X |
|
||
D* |
х = А2 |
(1 + ф2) = А2 + п2. |
(8.7) |
|
LU |
Д Д 2’
Здесь б?3, S — соответственно модуль упругости заполнителя при поперечных сдвигах и его толщина; Н — расстояние между нейтраль ными поверхностями несущих слоёв.
Если заполнитель достаточно слабый на поперечные сдвиги (шх з> 1, К < 0,5 Т*), то из соотношений (8.7) прямой минимизацией
можно получить:
|
Ткр = T0HV l - Ь 2 + К; |
|
|
Л2р + гг2р = А2С\ / 1 - Ь 2; Ь = Щ ^< |
1; |
|
1о |
|
Анс2 |
B R 2 |
(8.8) |
£ )чс ’ |
При выводе формул (8.8) принималось во внимание характерное для трёхслойных конструкций соотношение D* к D.
Если жёсткость К заполнителя на поперечные сдвиги достаточно ве лика (UJX <с 1, К ф 0,5Т^), то анализ устойчивости можно проводить
8.3. Изотропные трёхслойные оболочки с лёгким заполнителем |
197 |
на основе модели прямой линии, В этом случае соотношение (8.7) примет вид:
D |
В |
D |
D |
Ъ 2 Х + — ’ |
х = А2(1 + ф2) = А2 + п2. (8.9) |
||
Rz |
х |
|
1 + и>х’ |
Здесь величина шх <^i 1, и, как показал численный анализ, она слабо влияет на форму волнообразования. Учитывая это, проминимизируем выражение (8.9) по параметру х , считая величину D постоянной. Тогда
ГГ.Ж |
/уж |
(8.10) |
Т кр= /-g— ; |
xKP = \ i V Y T T 0 - £0 = 4 , . |
|
V 1 + £о |
|
|
Если пренебречь квадратом малой величины £о, то расчётные фор мулы (8.10) примут вид:
Т кр = |
7 ;>. |
_о_ |
|
|
4К |
|
Т'ж. |
B R 2 |
|
l + ifi- |
|
А кр + « к р |
(8.П) |
|
|
4К |
|
Как видно из (8.8), (8.10), (8.11), критические усилия определяются тремя величинами: Тцс — критическое усилие, соответствующее обо лочке с раздельно работающими несущими слоями; Тж — критическое усилие трёхслойной оболочки с абсолютно жёстким (К —*• оо) на по перечные сдвиги заполнителем; К — жёсткость заполнителя на попе речные сдвиги — критическое усилие при так называемой сдвиговой форме потери устойчивости (b = 1). Этими тремя величинами (Т™, Тж, К) определяются обобщённые жёсткости трёхслойной изотропной сферической оболочки в задаче об устойчивости.
Соотношение между обобщёнными жесткостями оболочки опреде ляет две характерные области потери устойчивости. Так , в области слабых на поперечные сдвиги заполнителей (К < 0,5Т^, b < 1), как следует из соотношений (8.8), критическое усилие пропорционально жёсткости К заполнителя и слабо зависит от жёсткости Т™ несу щих слоёв, поскольку для силовых трёхслойных конструкций харак терна зависимость К Т™. Если заполнитель достаточно жёсткий (К > 0,5 Т*, b > 1), то, как следует из соотношений (8.10), (8.11), критическое усилие определяется обобщённой жёсткостью Тж, а попе речные сдвиги играют поправочную роль и слабо влияют на критиче ское усилие. Поправка от сдвигов определяется величиной
|
Т0Ж |
_ В |
А с |
Ж |
( 8 . 12) |
|
~ Ж Ж ' |
Если Дс <с 1, то можно применять классическую модель оболочек с неизменной нормалью и пользоваться формулами (8.5).
198 Гл. 8. Устойчивость многослойных сферических оболочек
Зависимости (8.10), (8.11) пригодны также для расчёта на устойчи вость многослойных сферических оболочек, сопротивляющихся по мо дели прямой линии либо по классической модели. Для этого доста точно при вычислении величины Г * подставлять соответствующие жёсткости оболочки на растяжение и изгиб; кроме того, необходимо принимать во внимание жёсткость К многослойного пакета на попе речные сдвиги [32, 41].
8.4. Устойчивость многослойных ортотропных оболочек (классическая модель)
Рассмотрим ортотропные оболочки сначала на основе классической модели. Критические характеристики таких оболочек определяются минимизацией выражения (8.3) по параметрам волнообразования А и 1ф. Проводя минимизацию по параметру А, найдём [71]:
Т = ^ = к ( ф ) ^ ^ Т Ш = к ( ф ) Т 0; То = | 0 ^ А ;
к(ф) |
Р\{Ф) Л2 _ \ 2 ! + V > 2 |
^2 _ |
1В 2Д 2 |
(8.13) |
|
В Д Г |
ж |
V А |
|||
|
|
Величина Го соответствует критическому наружному давлению при осесимметричной форме (ф = 0) потери устойчивости. При этом кри тическое значение параметра А становится равным Аж (Акр = Аж).
В общем случае волнообразование неосесимметрично (ф Ф 0), и для определения критических параметров Ткр, Акр, фкр необходимо прове сти минимизацию функции к(ф) по ф. Проводя минимизацию, получим формулы для расчёта критической нагрузки и параметров волнообра зования при потере устойчивости многослойных ортотропных сфериче ских оболочек при действии наружного давления (од < а 2):
Т кр = Р к р И |
кортГ), |
к,орт |
В Д кр) _ |
/ а 1+ 2 i/3 ~ ^ |
|||
|
2 |
|
|
Р 2 ( Ф к р ) |
у «2 + 2\ f f i |
|
|
|
2 |
|
I + ф^р |
|
А2 |
2 |
|
|
Акр |
|
|
|
кр' |
|
|
|
|
|
|
^кр |
|
||
|
|
В Д к р ) Г 2 ( ^ к р ) |
|
|
|
|
|
/2 |
\ J (A-/?2)2+ ( a 2- a i) { P \ o t 2 - l k o L i ) - (А - |
/32) |
I |
,01/1. |
|||
фкр = |
- -------------------- |
д--------- |
д--------------------------- |
|
|
« - 7=-. |
(8.14) |
|
|
P\0t2-P20t\ |
|
|
у/Р |
|
|
Как |
следует из |
(8.14), |
параметры |
волнообразования Акр |
и пкр |
в рассматриваемом |
случае однозначно определены. Если закономер |
ности, замеченные |
для цилиндрических оболочек [4, 69], верны |
и для сферических, |
то однозначная определённость волнообразования |
8.4. Многослойные ортотропные оболочки (классическая модель) 199
при потере устойчивости указывает на возможность удовлетворитель ного согласования теоретических и экспериментальных результатов. Как известно, в случае изотропных оболочек, когда форма волнооб разования теоретическим путём однозначно не определена, различие между и теоретическими и экспериментальными значениями критиче ских усилий достигает Зч-4 раз [8, 18]. Систематизированные экспе риментальные результаты по устойчивости ортотропных сферических оболочек нам неизвестны.
При подстановке в формулы (8.13), (8.14) соответствующих жёстко стей можно рассчитать критическое внешнее давление для вафельных сферических оболочек [40].
Для однородных по толщине h ортотропных сферических оболочек зависимости (8.14) принимают вид:
укр _ РкрД _ £ |
Ф',кр |
ор т Т0; |
Л2 |
= А2 |
____1+ \J~fi_____ |
А1 = 3,46, |
I Еч( 1 —V\V2) |
R |
|||
^кр |
''я |
а 1 + 2^/Д ) («2 + 2\[Ji) |
|
ЖI |
h' |
|||
|
|
|
|
|||||
|
|
|
|
|
|
|||
|
_ |
л/ЁфЁч |
h2 |
_ |
/2G]2( 1+ yJV\V2 ) |
(8.15) |
||
|
~ |
/3 (1 -1 /1 1 ^ ) |
R ' |
°рт |
V |
|
||
|
|
|
|
Здесь Ei, Е2, G12, v\, щ — характеристики упругости ортотропного материала оболочки.
Если анизотропия оболочек такова, что а\ > «2, то оболочки теряют устойчивость по осесимметричной форме. Это следует из соотношения (8.13), поскольку в этом случае величина к(ф) достигает минимума при фКр = 0 (/сорт = О- Критические параметры при этом определяются соотношениями:
Ткр= ^ = Т о = | У в д ; Акр= А 0; ^ кр= 0 (пкр = 0). (8.16)
Расчёт критических нагрузок в случае а\ = ач также проводится по формуле (8.16). Однако форму волнообразования в этом случае нельзя определить однозначно. Анализ зависимости (8.3) показывает, что критические параметры волнообразования при он = ач определя ются формулой
АкР + у/Р-«кР = А^ |
B 4 R 2 |
(8.17) |
|
|
“ о Г |
Таким образом, форма волнообразования при потере устойчивости от наружного давления в случае а\ = ач может быть любой, лишь бы выполнялось условие (8.17).
200Гл. 8. Устойчивость многослойных сферических оболочек
8.5.Устойчивость трёхслойных ортотропных оболочек
При исследовании устойчивости трёхслойных сферических оболо чек примем наиболее общую модель — модель ломаной линии. При этом критические параметры получаются минимизацией выражения (8.1) по Л и ф. Применяя приемы асимптотического анализа [72], можно получить простые формулы для расчёта критических нагрузок и параметров волнообразования [71].
Так, в области заполнителей, слабо сопротивляющихся поперечным сдвигам, когда
К i + К 2ф2 |
R 2 |
Л2 Fi(tp) |
(8.18) |
А |
функцию П влияния поперечных сдвигов можно представить в виде
R 2 К х + К 2ф2 ( |
R 2 К х+ К 2ф2 |
(8.19) |
||
А ' АЩ ( ф ) |
А ~ |
' AЩ ( ф ) |
||
|
Подставляя это значение П в формулу (8.1), получим:
A ^ I A ) . А |
1 А + А < А 2. |
|
R 2 Х+ F2 ' х + 1 + ф2 ’ |
|
|
2А |
1 + 1 У |
(8.20) |
В* В 2 1 |
1 + ф2 |
|
|
|
Проминимизируем это выражение по параметру х и найдём:
'НС |
1 - |
2А А А ; |
+ |
А + А А . |
Т = к(ф) Тс |
|
|
||
О |
|
1 + ф2 |
|
1 + А |
\V
т0™ = д v 'A A
Л2 |
= Л2 |
1 +Ф1р |
1 |
( 2А |
1 + -^ ^2Р |
|
|
( . ) |
|||||
^кр |
''нс |
|
|
'Тук |
1 |
8 21 |
|
|
^ ( А р ) т2(а р) |
^ |
Г " |
+ Ф1р |
|
|
|
|
|
|||
|
|
|
|
|
|
Критическое |
значение |
параметра ф находится минимизацией по |
ф величины Т |
в формуле |
(8.21). Рассмотрим случай изотропного |