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

книги / Строительная механика и расчеты композитных конструкций на прочность

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

Вычислим деформации и изменения кривизны. Согласно полубезмоментной теории е;. = уху = 0; аех = х.ч. = 0, поэтому необхо­

димо определить только ех и агу

 

 

 

 

ди

Э2Ф

 

_

 

ех = — = —-у =

тх 81п пу;

 

 

дх

дх2

'

 

 

= _ Э ^ + _1_Эу _ _ Л Э^Ф__1_Э2Ф _

(2.176)

}

Эу2

Я Эу

Эу4

/? ду2

 

п А^ л ^ з ш

Я

х з т й ,

= -Ф 1Ш- ^ п 2[п2 -1)зт/йхз1

Угол со;. поворота нормали вокруг оси Ох

 

 

 

Эи>

1

, Э3Ф

1 ЭФ _

 

 

 

ду + Я У

ду'

Я ду

(2.177)

 

 

 

 

 

 

= Ф |* .(лл 3 - ^ я ) и п |Г

 

= Ф„„, Д*-п (л2 -1) зш тх соз пу.

При формулировке задачи устойчивости будем учитывать, что имеем дело с нормальным внешним давлением р. Распределен­ ные силы ос таются нормальными к деформированной поверхно­ сти (рис. 2.61), поэтому появляется проекция —рту на начальную касательную. Работа касательных сил на перемещениях 5у будет

5А = —р(йуЬу.

Поскольку согласно полубезмоментной теории цилиндричес­ ких оболочек из внутренних силовых факторов учитывают толь­

ко осевую силу ТУ* и изгибающий момент

соотношения упру­

гости запишем только для них:

 

N. = Вхех; М х = Б22я у,

(2.178)

б

 

У Л = -

Рис. 2.61. Деформирование поперечного сечения при потере ус­ тойчивости:

а — до потери устойчивости; б — после потери устойчивости (соответствующие линии у и уу параллельны)

где Вх — коэффициент мембранной жесткости (погонная жест­

кость в направлении образующей), Вх = Вп - В22/В 22; 3 22 — ми­ нимальная погонная изгибная жесткость в окружном направле­

нии, />22 = ^22-С Ъ /В 22 ■

Запишем вариационную формулировку задачи устойчивости для рассматриваемого случая:

2яЯ/

I

I(№х8ех + М у8 х у+АГу<оу8<оу + рыу8ч}с1х<1у = 0, (2.179)

о

о

где Ыу — начальное окружное погонное усилие от внешнего дав­

ления,

N ° = - р Я ; 8еЛ., 8®я

8о>я

8у принимаются в следующем

виде:

 

 

 

 

 

 

8ех = - 8Ф,„пт 2 з т т х з т п у ,

 

8ху = - 8Ф т„ -^32 2 - 1 )$ т т х 5т пу,

 

 

 

 

 

(2.180)

 

8соу = 8Фтл -р-п(п2 - 1)$ т т х соз пу,

 

-

Э5Ф

г .

1 .

_

 

8у= - = -

8Ф„т

—Л8Ш

тх созпу.

 

 

ду

 

К

 

Поскольку

 

 

 

 

2яй I

 

2пЯ I

 

|

| $ т 2 етхзт2 пус1хс1у = |

$$ т 2 тх соз2 пу сЬсНу =

о о

 

о о

2

после подстановки (2.176)—(2.178), (2.180) в (2.179) получим

Вхп г и П п - ~ п А(и2 - 1)2 =

=- 1)2 + Р ^ п 2(п7 - \) = р ± { п 2- \ ) п \

Отсюда находим выражение, из которого определяется значе­ ние критического давления,

ВхтАЯг М ”2-о

л4 2 - 1)

Л1

Из (2.181) видно, что наименьшее давление по числу т полу­ волн в осевом направлении соответствует т = 1. Тогда вместо (2.181) запишем

(2.182)

Для длинных оболочек (Л //«: 1) первым слагаемым в (2.182) можно пренебречь, и при п = 2 получим

(2.183)

Критическое давление, при котором происходит потеря ус­ тойчивости, определяется из (2.183) минимизацией значения р по числу и волн (рис. 2.62). При большом числе волн п » 2 можно считать п- - п г — 1 = т| » 1. Параметр г) можно рассмат­ ривать как непрерывный и выполнить для (2.182) минимизацию аналитически.

Для этого представим структуру зависимости (2.182) в виде

(2.184)

Тогда из условия минимума /;(г|)получим

(2.185)

Р

2

4

6...

л

Рис. 2.62. Определение критического давления

Подстановка г\ (2.185) в (2.184) дает значение критического давления

3аЬ3'4 Ь(3а)1/А

=3(3а)3/А + Ь*»

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

Р кР

 

(2.186)

’ /л Тл

 

при этом, как следует из (2.185), число волн можно оценить как

(2.187)

Пример 2.15

►Дать оценку критического давления для цилинарической обо­ лочки, изготовленной из изотропного материала и из КМ, имею­ щего структуру [±<р]. Для КМ воспользоваться нитяной моделью,

а) Изотропный материал

 

3 = ЕИ\

Д22 = -

 

 

 

 

" 1 2 (1 - ц г)'

Для длинной оболочки / »

Я согласно (2.183)

р

=

4(1- р 2)

-0*27-

Рхр Я2

[Я ]

{Я )

Для оболочки средней длины согласно (2.186)

.

г Г * У * ’1*

____ .

0,92е Щ ) *

Р“ 5,5 1 Ш

Е ( ' ) ( л )

 

[12(1 -(Л1)]■Ч Т * ‘

б) Композиционный материал [±<р] Воспользовавшись нитяной моделью, вычислим

Е,кг . а

Для длинной оболочки / »

Я согласно (2.183)

22

_ (И \3 1 . 4

= -

/ = г ' Ы

4 ЯП *•

Для оболочки средней длины по (2.186)

 

_ с

Я (И ]5П 5,

1 Я ^

1

/ Ы

— р -С 0 5 ф 3 1 П 'с р =

12 3

 

п / и

(2.188)

»0,853Е,у1 — I

созф5т 3 ф.

Примерная зависимость значения критического давления от угла укладки показана на рис. 2.63.

<р,° К

30 0.092

45 0,213

60 0,277

75 0,199

о

30

60

90

Рис. 2.63. Зависимость значения критического давления от угла ф укладки

Из условия

— (с О З ф З Ш 3 ф ) = - 5 Н 1 4 ф + СОЗ2 (ЗЗН 12 ф ) =

= з ш 2 ф ( 3 с 0 5 2 ф - 5Ш 2 ф ) = 0

получим 1§3ф = 3 и при ф = 60е значение критического давления будет максимальным. ■

Пример 2.16

►Для длинной цилиндрической оболочки, нагруженной вне­ шним давлением, определить толщину А, при которой потеря устойчивости и разрушение происходят одновременно. Структу­ ра пакета (907(±45°)/90°], И%= А, А45 = 2А (рис. 2.64). Воспользо­ ваться нитяной моделью.

Рис. 2.64. Структура многослойного пакета

Согласно (2.165) критическое давление Ркр для длинной цилин­ дрической оболочки

Щ 2 К3

Определим изгибную жесткость Л22 для данного пакета (см. рис. 2.63)

Тогда

(2.188)

Мембранные жесткости данного пакета имеют виц:

= Е}к

Воспользуемся соотношениями упругости:

= Яп ех + Вп в} — 0;

Му = В12ех + В27еу = РкрЯ.

Отсюда определяем*

*п/В и

или

Ркр к 1

А 2 '

Наиболее нагруженными будут слои кольцевой подмотки = 90°:

! = Е\ (ех С052 ф + еу ып2 ф) = ркр -

(2.190)

где а, — предел прочности ОКМ при сжатии вдоль направления армирования. Подставив (2.189) в (2.190), получим

 

- т

’ 29 Л I _ „

29

 

 

' - Е' Ш

т * 2 - г , д а

4 -

 

Тогда

 

 

 

 

 

Е| 729

 

 

где

можно принять за характеристику материала.

 

Пример 2.17

 

 

 

►Для оболочки средней длины опреде­

2

 

лить толщину А, обеспечивающую устойчи­

 

 

 

вость оболочки, нагруженной внешним дав­

90° Ж

А/3

лением р. Структура многослойного пакета

А/3

[±4570790°], А45 = А0 = А*, = А/3 (рис. 2.65).

Для расчета воспользоваться нитяной моде­

± 4 5 °

А/3

лью.

 

 

 

 

 

Примем в качестве координатной повер­

0

 

Рис. 2.65. Структура

хности внутреннюю поверхность оболочки

многослойного

па­

(см. рис.

2.64).

 

 

 

 

Подсчитаем необходимые мембранные, смешанные и изгибные жесткости:

кета

ви = в„ . в 1а! ( 1 + 1) = г 1а^ ;

Ъ - Ъ - Й - М з т И у т ) - ^ ° - 0335;

Д , - Д , - А . е 4

Согласно (2.186) запишем Л*2 /"Л**

(1Г“^ ,9;* - К ^ юГ

2 .7.2. Кручение

Рассмотрим задачу об устойчивости длинной тонкой много­ слойной цилиндрической оболочки при кручении (рис. 2 .66).

Для решения воспользуемся гипотезами полубезмоментной теории цилиндрической оболочки (2.173). Введем функцию пе­ ремещений Ф так, чтобы приближенно описывать винтовые склад­ ки (2.174). Функцию Ф аппроксимируем (см. рис. 2.66)

Ф = Фтп$ т {тх + пу),

т = тп/1; п = п/Я .

Необходимые для решения деформационные соотношения будут иметь вид:

=~ф тпт 51П (,пх+«3').

Э2в>

1 Эу

Э4Ф

1 Э2Ф

(2.191)

г Эу2

Я ду

Эу 4

Я ду2

 

= ~Ф/пл -р-и2 2 - 1^ $ т(тх + пу).

Рис. 2.66. Определение по терн устойчивости цилинд рической оболочки при кру чении

Углы поворота нормали

(2.192)

Соотношения упругости

(2.193)

Возможные значения деформации бед., изменения кривизны бае^ и углы поворота нормали аппроксимируем аналогичным (2.191), (2.192) образом:

(2.194)

Запишем вариационную формулировку задачи устойчивости полубезмоментной цилиндрической оболочки при кручении:

2яЯ /

О О

где

 

 

 

 

I

(2.196)

 

8т|^. = 8(01,00^) =

 

+ 8ш>,сол..]

 

 

 

Подстановка (2.191)—(2.194), (2.196) в (2.195) дает

8Фщл Фшп

+ Д а ^ б « 4 («2 -

!)2 -

 

( " 2 " 0 '" ] = °-

Отсюда получаем

 

 

 

 

 

ВхЕ?тг

В>22п(пг ~\)

 

(2.197)

 

~ 2 п г ( я 2 “ 1 ) +

2 В ? т

 

 

 

 

Представим структуру формулы (2.197) в виде

 

 

 

1

 

(2.198)

 

= ат

+ Ь—,

 

 

 

 

т

 

 

 

ВХКУ

.

5 22я (л2 -1)

 

 

 

" 2 н3 ( « - - 1) ’

 

2 Л3

 

Для длинных оболочек можно ожидать, что число полуволн в осевом направлении т достаточно большое число. Поэтому в (2.198) выполним минимизацию, считая т непрерывным пара­ метром:

К

Ьт

Отсюда получим

Ь14

т

( 3 * Г

3аЬ*

Ь{3а)ил

3(3а)*4+

\ ( ъ - Г ь * =

Ь'!'