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

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

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

Компоненты тензора линейных деформаций в произвольной точке оболочки определяют выражениями

x ex

z x

 

 

 

 

 

z y

 

y ey

 

e

z

0

 

 

 

 

 

 

 

 

 

 

 

 

(1.12.4)

xy exy z xy

 

 

 

w

 

xz

x

 

 

 

 

 

 

 

1 w

 

 

 

 

 

 

уz

 

 

 

R

 

 

 

 

 

Тензор (1.12.4) – поле возможных деформаций первого родав выражении энергетического критерия устойчивости ( Э),

для построения выражений для компонент тензора учтем линейные слагаемые в тензоре мембранных деформаций:

 

 

 

 

 

 

1

с 2

 

1

2

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

2

 

 

 

 

 

2

xz

 

 

 

 

 

 

 

 

1

 

w

2

 

1

 

2

 

 

 

 

 

 

 

 

 

 

y

2

 

 

 

 

 

 

 

 

2

yz ,

(1.12.5)

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

xy

w

w

xz yz .

 

 

 

 

 

 

 

 

 

 

 

 

 

x R

 

 

 

 

 

 

Запишем выражение для приращения полной энергии цилиндрической оболочки:

Э Т е dV 0 Т

 

dV .

(1.12.6)

 

V

V

 

Подстановка (1.12.4) в первое слагаемое позволит произвести интегрирование по толщине оболочки, т.е. по координате z (от – h/2 до h/2, где h – толщина оболочки)

2 е

 

2

 

2

 

2

Э

11 x

22 y

2 12

x y 66 xy

 

 

S 0

 

 

 

 

 

 

71

 

 

А e

 

2

А

 

 

e

y

2 2А e e

y

А e 2

 

 

 

 

 

11

x

 

 

 

22

 

 

 

 

12

x

 

66 xy

 

 

 

 

 

55

 

xy

2

 

 

44

 

yz

2

 

 

Rdxdy

 

 

 

 

 

 

А

 

 

 

 

 

 

 

А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 l

0

 

1

 

2

 

 

 

0 1

2

 

 

0

 

 

 

 

 

(1.12.7)

 

 

Tx

 

2

xz Ty

 

2

yz

Txy

xz yz Rdxdy.

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В выражении Э

 

Dij

компоненты тензора изгибных жест-

костей

( D C h3

/ 12 );

 

A

 

 

 

тензор

мембранных

жесткостей

 

ij

ij

 

 

 

 

 

 

 

 

ij

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( A C h );

T 0 ,T 0 ,T 0

 

тензор мембранных усилий для первона-

ij

ij

x

y

xy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

чального (невозмущенного) состояния цилиндрической оболочки.

( Э) 0

(1.12.8)

Для замкнутой цилиндрической оболочки на ее торцах должны быть заданы граничные условия:

перемещение U или осевая елка Tx ;

перемещение v или угловое усилие Sxy ;

перемещение w или поперечное усилие Qx ;

угол поворота нормали или изгибающий момент Mx . Если на оболочку действуют:

осевое усилие (сжимающее)

закручивающий момент

внешнее давление

тотензорымембранныхусилийопределяются следующимобразом

Tx0 N , 2 Rh

Ty0 pR p2 RL , 2L

72

T 0

 

M

.

(1.12.9)

2 R2

xy

 

 

 

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

 

 

 

 

 

 

 

 

Э U ,V , w

 

 

 

 

 

 

 

 

u

 

u

 

v

 

v

 

w w 2 w 2 w

 

2 w

 

(1.12.10)

 

x

 

 

 

x

 

 

 

x

 

 

 

x2

 

2

 

x

 

 

 

F w,

 

,

 

,

 

,

 

,

 

,

 

,

 

,

 

,

 

, Rdxd .

 

Составим систему уравнений Остроградского – Эйлера:

 

F

 

 

F

 

 

 

F

 

 

2

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

x ( w/ x)

 

y ( w/ y)

x

2

 

 

(

2

w/ x

2

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

F

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0;

 

 

 

 

 

y2

( 2 w/ y2 )

x y( 2 w/ x y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.12.11)

 

 

 

 

 

F

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

( u/ x)

 

y

 

( u/ )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

x ( v/ x)

y ( v/ )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так,

для

осесимметричной

формы потери

 

устойчивости

(u = v 0) при действии продольной нагрузки

d 4 w Еx h

xdx4 R2

w T 0

d 2 w

0.

(1.12.12)

x

dx2

 

 

При проведении практических расчетов устойчивости оболочечных конструкций обычно используют вариационную подстановку ( Э) = 0 для получения разрешающих соотношений, например, в форме МКЭ.

1.13. Нелинейные задачи теории устойчивости

Линейная теория устойчивости, в рамках которой получен энергетический критерий бифуркационной устойчивости и линеаризованные уравнения устойчивости, не позволяет исследовать поведение механической системы после потери устойчивости.

73

Причины этого являются:

1.Конечные перемещения (большие деформации).

2.Неупругое поведение материала (физическая нелинейность) при наличии больших деформаций

3.Динамические эффекты, сопровождающие переход из неустойчивого состояния равновесия в смежное устойчивое.

Кроме того, для элементов конструкций могут не выполняться и основные гипотезы теории устойчивости:

1.Первоначальное состояние системы может не описываться соотношениями линейной теории упругости (линейной теории стержней, пластин и оболочек).

2.Первоначальная форма может иметь отклонения от идеальной (стержень не прямой, пластина не плоская, кольцо не круглое).

3.Действующие на систему нагрузки не являются «мертвыми», неизменяемыми по величине и направлению (например, нестационарное воздействие).

Все это приводит к усложнению задачи устойчивости, определяющие соотношение при этом получаются нелинейными.

Таким образом, при анализе закритического поведения конструкций, а также в некоторых случаях и при анализе условий потери устойчивости механических систем необходимо рассматривать нелинейные соотношения устойчивости (что, безусловно, является более сложной задачей).

1.14.Закритическое поведение систем после потери устойчивости

Вариационные принципы механики:

1.Полная энергия стремится к минимуму (метод Лагранжа).

2.Дополнительная энергия стремится к минимуму (метод Кастильяно).

3.Мощность стремится к минимуму (метод Журдена). Полная энергия представляет некоторый функционал, функ-

ция под интегралом не является квадратичной относительно неизвестной функции (поля перемещений).

74

Э FdV.

(1.14.1)

V

 

Для функционалов такого вида несправедлив энергетический бифуркационный критерий устойчивости

( Э) .

(1.14.2)

Однако из вариационного исчисления известно, что условие( Э) является необходимым условием экстремума, а характер экстремума (минимум или максимум) определяется знаком второй вариации

2Э ( Э) 0 min

(1.14.3)

и, следовательно, предельное состояние для механической системы описывается уравнением

2Э=

(1.14.4)

которое и является критерием устойчивости.

Легко убедиться, что бифуркационный критерий частный случай приведенного соотношения. Разложение в ряд Тейлора приращения полной энергии дает

Э Э

1

 

2Э ...,

(1.14.5)

2!

 

 

 

следовательно, Э = Э, если в разложение Тейлора ограничиться одним слагаемым, и тогда

2Э ( Э) ( Э).

(1.14.6)

Учитывая, что F в выражении для полной энергии содержит достаточно большое число неизвестных функций, то в выраже-

нии 2Э будет достаточно большое число слагаемых вида:

 

 

2

F

 

 

2Э ...

 

p q ... dV.

(1.14.7)

p q

V

 

 

причем, если неизвестных функций n, то слагаемых n2 .

75

1.15. Динамическая устойчивость

Если J интеграл действия, то необходимое условие экстремума функционала (принцип Гамильтона)

2 J 0,

(1.15.1)

а условие устойчивости движения

 

2 J 0.

(1.15.2)

1.16. Учет начальных несовершенств

Поле перемещений, как стационарное, так возмущенное, представляется в виде:

U u0 u ,

(1.16.1)

где U полное поле перемещений; u0 начальное поле переме-

щений; u отклонение от первоначального поля перемещений. Наличие некоторых функций {u0} в предположении, что не-

совершенства (их форма) известны, дают дополнительные слагаемые в выражения для 2Э, что может повлиять на знак. При этом изменение может происходить как в сторону увеличения, так и в сторону снижения критической нагрузки.

Кроме того, следует иметь в виду, что начальные прогибы несовершенства имеют статистическую природу, и задача оценки устойчивости конструкции в этом случае является вероятностной.

76

Глава 2. ФИЗИКО-МЕХАНИЧЕСКИЕ СВОЙСТВА КОМПОЗИЦИОННЫХ МАТЕРИАЛОВ

Композиционные материалы благодаря высоким удельным физико-механическим свойствам нашли широкое применение в ответственных высоконагруженных элементах конструкций, например: в ракетной технике – корпуса, обтекатели, сопла двигателей в авиации – панели фюзеляжа, оперения, плоскостей

ввоенной технике – защитные шлемы, жилеты, на транспорте – шины, пневмоамортизаторы, корпуса спортивных автомобилей,

всудостроении – корпуса кораблей, яхт, шлюпок, в строительстве – панели, перекрытия зданий и сооружений, в технических видах спорта, например в бобслее, – экипировка спортсменов.

Прогнозирование эффективных физико-механических свойств композиционных материалов и рациональное (оптимальное) проектирование композиционных структур, например, на основе варьирования таких структурных параметров, как взаимное расположение

ифизико-механические свойства элементов структуры, составляют одну из центральных задач механики композитов, что позволяет создаватьматериалыс заранее заданными свойствами.

Вволокнистых композиционных материалах элементы структуры-волокна диаметром 5÷200 мкм, физико-механические свойства некоторых наиболее часто используемых видов волокон приведены в табл. 2.1.

Строение тонкостенной конструкции из волокнистого композиционного материала можно представить в виде многоуровневой структурной модели, например: армирующие волокна в связующем или матрице, ткань, лента, слой, пакет слоев.

1. Статический разброс прочности волокон. Функция на-

дежности (рис. 2.1) вида

 

 

 

 

 

 

R( ) exp

 

 

 

(2.1)

 

 

s

 

 

 

 

 

 

 

 

 

77

есть распределение Вейбулла, где s – параметр уровня, α – параметр формы; вероятность разрушения

 

 

P( ) 1 R( ).

 

 

(2.2)

 

 

 

 

 

 

Таблица 2 . 1

Физико-механические свойства волокон

 

 

 

 

 

 

 

 

 

 

Диа-

Плот-

Модуль

 

Предел

Удельные

Волокно

метр d,

ность ρ,

Юнга

 

проч-

свойства ГПа м³/кг

 

мкм

кг/м³

Е, ГПа

 

ности

Е/ρ

σ/ρ

 

 

σ, ГПа

Стеклянное

5÷20

2500

89

 

3,5

0,036

0,0014

Углеводородное ВМ

1950

400

 

2,1

0,205

0,0011

Углеводородное ВП

 

1760

260

 

3,3

0,148

0,0019

Борное

100÷20

2500

400

 

3,5

0,16

0,0014

Органическое

1450

120

 

2,8

0,083

0,0019

0

 

Стальное

7800

210

 

2,9

0,027

0,0004

 

 

Рис. 2.1. Функция надежности

Среднее значение и коэффициент вариации Кv быть рассчитаны по формулам:

s Г 1 1/ ,

 

 

 

 

 

2

1/2

 

 

Г 1

 

 

 

 

 

 

 

 

Kv

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1 ;

 

Г

2

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

могут

(2.3)

значение коэффициента а, например, для стекловолокна 3...4 и для углеволокна 6...9.

78

2. Масштабный эффект прочности волокон. Пусть на волокна, функция надежности

 

 

 

 

 

R( ) exp

 

 

 

.

 

 

s

 

 

 

 

 

 

 

l0 дли-

(2.4)

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

 

 

l

 

 

 

Ri ( ) R( ) exp

 

 

 

,

(2.5)

 

l

l0

 

s

 

 

 

 

 

 

 

 

 

 

 

 

l0

 

 

 

 

 

 

 

 

 

 

 

где l0 базовая длина; s, a – параметры волокна на базовой длине; закон Малмейстера при β≈1/2

 

l

 

 

 

 

 

 

R( ) exp

 

 

 

 

 

 

.

(2.6)

l

 

 

 

 

 

 

 

s

 

 

 

 

0

 

 

 

 

 

 

3. Матрицы. Физико-механические свойства различных матриц приведены в табл. 2.2; в круглых скобках – соответствующие удельные значения Е/ρ [ГПа м3 /кг] и σр /ρ [МПам3 /кг].

2.1. Упругие характеристики анизотропных материалов

Обобщенный закон Гука имеет вид

ij aijkl k ,

(2.1.1)

ij Aijkl kl ,

(2.1.2)

где i,j,k,l = 1,2,3; а – тензор податливости, А – тензор жесткости. Соотношения (2.1.1), (2.1.2) могут быть записаны в матрич-

ной форме

j aij j ,

(2.1.3)

i Aij j .

(2.1.4)

79

Таблица 2 . 2 Физико-механические свойства матриц

 

 

Плот-

 

Пределы

Тепло-

 

 

Модуль

прочности

стабиль

 

 

ность

 

Матрицы

,

Юнга

на растяжение

биль-

 

 

£, ГПа

и сжатие, МПа

ность

 

 

кг/м3

 

 

 

S+

S

Т, °С

Термореактивные:

 

 

 

 

 

1.

Полиэфирные

1200–

2,8–3,8

30–70

80–150

50–80

1350

 

 

 

 

 

 

2.

Фенолформальдегидные

1200–

7–11

40–70

100–125

140–180

 

 

1300

 

 

 

 

3.

Эпоксидные

1200–

2,4–4,2

35–100

90–160

130–150

1300

 

 

 

 

 

 

4.

Кремнийорганические

1350–

6,8–10

25–50

60–100

250–280

1400

5.

Полиамидные

1410–

3,2–5

90–95

250–280

250–320

1430

 

 

 

 

 

 

Термопластичные:

 

 

 

 

 

1.

Нейлоновые

1140

2,8

83

65

2.

Сополимер, этиленовые

1700

1,4

45

74

Углеродная

1900

30 (0,016)

32 (0,017)

62(0,033)

3000

Металлические:

 

 

 

 

 

1.

Алюминиевые

2700

72 (0,027)

70 (0,026)

780

2.

Титановые

4710

117 (0,025)

200 (0,042)

~900

3.

Стальные

7800

210 (0,027)

420 (0,054)

~1000

с использованием соответствий между компонентами, например, тензора и матрицы напряжений

11

12

13

 

1

6

5

 

 

 

 

22

23

 

 

2

4

 

,

(2.1.5)

 

 

 

 

 

 

33

 

 

 

3

 

 

 

sim

 

 

sim

 

 

 

 

где индексы i, j 1,6 .

Упругие характеристики однонаправленного волокнистого композиционного материала. Будем использовать следующие предположения: а) композит обладает трансверсально-изотропной симметрией упругих свойств, б) матрица и волокна в композите – одно-

80