книги / Строительная механика, динамика и устойчивость композитных конструкций
..pdfКомпоненты тензора линейных деформаций в произвольной точке оболочки определяют выражениями
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