книги / Прямые и обратные задачи механики упругих композитных пластин и оболочек вращения
..pdf22 Гл. 1. Структурные модели композиционного материала
об однородном НДС в структурных компонентах — основное
предположение для этих моделей. |
|
|
||
Более |
сложные |
модели упругого |
поведения |
композицион |
ных материалов, учитывающие неоднородность |
полей струк |
|||
турных |
напряжений |
и деформаций, |
построены |
в работах |
Н. С. Бахвалова, Г. А. Ванина, Т. Исикавы, Р. Кристенсена, С. Кобаяси,
Б.Е. Победри, Ю. В. Соколкина, 3. Хашина, Р. Хилла, Л. П. Хорошуна,
Т.Д. Шермергора, С. Штрикмана и др. Такие модели требуют значи тельного объема вычислений уже на стадии определения эффективных
характеристик КМ, поэтому их использование при практических расчетах НДС многослойных конструкций с переменными параметрами армирования слоев не нашло пока широкого применения. Тем не менее, эти модели могут эффективно использоваться для оценки применимости используемых в расчетной практике приближенных моделей и их уточнения, когда это необходимо.
1.2. Определяющие соотношения полиармированного слоя
Связь между осредненными напряжениями аар, таз и деформация ми еар, 7аз в полиармированном слое (обобщенный закон Гука) имеет вид
&OLOL —Q'OLOL^'OLOL
аа0 = flQ3eQQ + |
арге(3(3 + а зз • 2 еар — аз@ 0; |
(1 .1) |
7аЗ |
Qaa^a3 Т" 9Q/3^”/33I |
|
где 0 — приращение температуры. Соотношения (1.1) называют также соотношениями термоупругости или, при отсутствии температурных слагаемых, просто соотношениями упругости.
1.2.1. |
Структурная модель армированного волокнами компози |
та [32] — |
модель В. В. Болотина (МБ) — стала основанием многих |
исследований в данном направлении и широко используется при моде лировании композитных конструкций. Модель основывается на пред положениях об однородности НДС в изотропных упругих волокнах и во всем объеме изотропной идеально упругой матрицы, совместности деформаций волокон и связующего в направлении армирования, равен ства напряжений в волокнах и матрице по остальным компонентам.
Для вычисления эффективных упругих коэффициентов однонаправ ленного волокнистого композита используются методы осреднения по Фойгту и по Рейссу, дающие следующие формулы:
Е[ —и>аЕа 4- сосЕс\ |
Е 2 |
ЕаЕс |
||
UJa.Ec + |
UJcEa |
|||
|
|
|||
1>12 = CJaUa + CJcUc, |
G = |
GgGc |
||
UJgGc + |
( 1.2) |
|||
|
|
UJcGa ’ |
1.2. Определяющие соотношения полиармированного слоя |
23 |
где отброшены члены, имеющие порядок квадрата от коэффициентов Пуассона. Здесь Е \, Е 2 — эффективные модули упругости в направ лении армирования и поперек его, G — эффективный модуль сдвига, и\2 — эффективный коэффициент Пуассона в плоскости слоя; Е, и, ш с индексами с и а — модули упругости, коэффициенты Пуассона и объ емные содержания связующего и армирующих волокон соответственно, при этом шс + ша = 1.
Из соображений симметрии тензора податливости имеем
1^21 = V\2E2E\
Формулы для эффективных коэффициентов теплового расширения имеют вид
«1 = и)аа а + ojca c, а 2 |
(jJaQLaEa -|- LJQOCEC |
(1.3) |
|
и)аЕа ~\~ШсЕс |
|||
|
|
В описании модели указано, что среди формул для эффективных характеристик наихудшие результаты получаются при использовании осреднения по Рейссу, в частности для G. Приведены оценки для G, получаемые применением вариационного метода, при этом показано, что нижняя оценка
G = ( 1 + |
U>a)Gg + |
LOcGc |
(1.4) |
||
|
(1 + |
U)a)Gc + |
LOcGa |
|
|
дает более точное приближение, |
чем |
выражение (1.2). Здесь, выше |
|||
и далее модули сдвига связующего и арматуры равны |
|
||||
/-1 _ |
Е |
|
|
Е |
|
J-yc |
|
|
|
|
|
с ~ |
2(1- И / с ) ’ |
a ~ 2 (\+ v ay |
|
Коэффициенты тензора эффективных жесткостей для однонаправ ленно армированного слоя в случае плоского напряженного состояния будут иметь вид
^ аааа = , |
- |
* |
» -^1122= « |
2' |
> -^1212 = G. |
(1-5) |
1 |
V\2V2\ |
1 - |
^12^21 |
|
|
Невыписанные величины получаются по правилу симметрии либо по лагаются равными нулю. Здесь и далее считаем а, (3 = 1, 2 и а ^ (3.
При описании коэффициентов соответствующих тензоров для поли армированного слоя проводится осреднение по толщине слоя, при этом осреднение по прослойкам связующего отделяется от осреднения по армированным слоям, в итоге имеем
CLZ E C + U), [^11111!{Ф) + ^22221р{Ф)+
+ 2(Ап22 + 2А т 2)1<\{ф)11{'ф)]р{ф)йф, (1-6)
24 Гл. 1. Структурные модели композиционного материала
azucE c . |
[-41122 + |
(4 ц ц + ^2222“ |
|
0-а.р —------- т + ^2 |
|
||
1- Vr |
О |
|
|
|
|
|
|
|
- |
2(ЛИ22 + ЪАт 2))1\{Ф)11{Ф)]р{Ф)<1ф, |
(1.7) |
7Г |
|
|
|
fl33 —O’zGc + Шг [Л1212 + (^4m i + ^2222“ |
|
||
О |
|
|
1 8 |
|
- |
2 (Л ц 22 + 2Ат2))1\(Ф)1\{Ф)\р{Ф)<1ф, |
|
|
( . ) |
AQ3 = UJZ [A n n la i^) - Ачччч^'ф )-
0 |
|
- ( А п п |
+ 2А ,ш )(11(ф)-11(ф))]11(ф)12(ф)р(ф)<1ф, (1.9) |
где u z и az = 1 —u z |
— удельные толщины слоев, содержащих арма |
туру, и прослоек связующего; 1\{ф) = COS(X/J); ki^P) = sin(0); р(ф) — плотность распределения волокон по углам, нормированная на отрезке [0,7г]. Плотность является обобщенной функцией ф: она может содер жать особенности типа дельта-функции.
В работе [2] для плоского напряженного состояния на основе исход ных предположений модели получены формулы, связывающие осредненные напряжения и деформации с напряжениями и деформациями в компонентах армированного слоя: в волокнах арматуры (добавочный
индекс а) и в связующем (добавочный индекс cl): |
|
||
_ |
Еса\\ + |
U>g(ucEg —VaEc)(T22 |
|
< 7 (ll)c l “ |
|
ujaEa + шсЕс |
’ |
_ Еасг\1 +ujc{vcEa —vaEc)a22 |
( 1. 10) |
||
< 7 (1 1 )a _ |
|
ШаЕа+ШсЕс |
|
|
|
||
|
Л |
|
Л |
0"(22)cl — cr(22)a |
= C’’22, |
CT(a/3)cl — CT(a/3)a |
— <7 a /3- |
Величины a aa, аар соответствуют системе координат, связанной с на правлением армирования в рассматриваемом слое, и определяются по формулам
<7ц = |
1\{Ф)(Ти + 1\{Ф)о22 - |
21\{Ф)к{Ф)°\2, |
|
<?22 = |
1\{Ф )<711 + /? (0 )< 7 2 2 + |
2 / i ( 0 ) / 2 ( 0 V l 2 , |
( 1. 11) |
^12 = 1\{Ф)к{Ф){<711 - <Т22) + Й ( 0 ) - 1\{Ф))<712-
1.2. Определяющие соотношения полиармированного слоя |
25 |
Там же выписана более полная формула для эффективного модуля Е%, в которой учитываются слагаемые порядка квадрата от коэффициентов Пуассона:
Е 2 |
ЕаЕсЕ\ |
( 1. 12) |
Е\ {uJaEc “Ь UJcEd^j UJaUJci^UaEc Z^c-^a) |
Наряду с данной моделью рассмотрим базирующуюся на ней модель армированного слоя [129] (МБГ), которая предназначается авторами для описания свойств резинокордных материалов. В качестве величин, характеризующих удельное содержания компонент, принимаются диа метр D круглого волокна и плотность укладки — характерное число нитей п на ширину I сечения представительного элемента плоскостью, перпендикулярной к направлению укладки волокон. Описанные харак теристики позволяют вычислить объемное содержание арматуры в слое
и а = 7г£>2/4 • n/(lh), |
(М 3) |
где h — толщина слоя. Для определения эффективного |
модуля Е 2 |
используется формула (1.12). Для учета сдвигов в слое и поперечных
сдвигов используются формулы: |
(1.4) для |
Л 1212 = G 12 и Л 1313 = G 13, |
(1.2) для Д.2323 = @ 2 3 • В этой |
модели |
нет описания температур |
ных компонент. При восстановлении напряжений в волокнах их напряженно-деформированное состояние рассматривается как одномер-
ное: |
Л |
|
|
СГ(11 ) а = Е ае п , СГ(12)о = СГ(22)а = ЦаЗ)а = 0, |
(1-14) |
а |
напряжения в связующем не восстанавливаются. |
Здесь |
бц = е\\1\{ф) + 2ei2/i(VOMVO + е22^ (^ ) — продольные деформации
вволокне.
1.2.2.Модели КМ с одномерными и двумерными волокнами (мод ли Ю. В. Немировского) являются развитием идей, заложенных в МБ. Пренебрежение «многомерностью» НДС волокна и рассмотрение его как нити позволило построить существенно более простые, с точки зрения математики, модели КМ с одномерными волокнами (МОВ, МОВУ) [217, 218]. Уравнения этих моделей, в отличие от МБ, линейны относительно таких характеристик КМ, как объемные содержания свя зующего и армирующих волокон, и просты, поэтому их удобно приме
нять при решении задач рационального проектирования композитных конструкций. Модель с двумерными волокнами (МДВ) [220] является, в свою очередь, уточнением МБ с более четкими и точными мате матическими формулировками и выводами формул для характеристик армированного материала.
Ниже будут выписаны соотношения для структурной модели с двухмерными волокнами. Уравнения для перечисленных выше моде лей с одномерными волокнами, а также нитяной модели (НМ) [260] вытекают из них как частные случаи. Материал армированного слоя оболочки полагается идеально упругим, с плоскостью симметрии,
26 Гл. 1. Структурные модели композиционного материала
параллельной отсчетной поверхности оболочки. Приведем основные предположения для рассматриваемых моделей [220].
1.Полиармированный слой (рис. 1.1) представляет собой изотроп ное упругое однородное связующее с внедренной в него регуляр ной сетью однонаправленных упругих изотропных волокон.
2.Число армирующих волокон достаточно велико, так что армиро ванный слой можно считать квазиоднородным.
3.Внедренные в связующее армирующие волокна воспринимают как растягивающие, так и сжимающие усилия.
4.В процессе деформирования удлинения и сдвиги остаются малы ми по сравнению с единицей.
5.Элементы композиции соединены идеально: отсутствует про скальзывание между армирующими элементами и связующим.
6.Расстояние между армирующими элементами достаточно велико
по сравнению с их поперечными размерами и в то же вре мя достаточно мало по сравнению с рассматриваемым элемен том оболочки, поэтому локальными эффектами вблизи волокон и нерегулярностью деформации между смежными волокнами бу дем пренебрегать.
7.а. Армированный слой представляет собой сеть одномерных волокон; нагрузка воспринимается только армирующими волокнами, а напряжения в связующем не учитываются (рис. 1.2, а).
б.Армированный слой — сеть одномерных волокон; попереч ные сдвиговые напряжения воспринимаются только слоя ми связующего, а армированные слои являются абсолютно жесткими на сдвиг (рис. 1.2,6, в).
в.Армирующие волокна имеют прямоугольное поперечное се чение (рис. 1.2, г), а их НДС рассматривается как двумерное.
Предположения |
1 — 6 |
общие для |
всех |
рассматриваемых моде |
лей. Предположение |
7. а |
соответствует |
нитяной модели, предположе |
|
ния 7 б, в — моделям с |
одномерными |
и с |
двумерными волокнами |
|
соответственно. |
|
|
|
|
1.2. Определяющие соотношения полиармированного слоя |
27 |
Компоненты тензоров эффективных жесткостей и эффективных тепловых жесткостей выпишем в виде [9, 220]
N N
0>аа —&гЕс |
^ |
' Аа а , |
|
Ла/З —0.z^cEc Л ^ ^ |
|
|||||||
|
|
|
п=1 |
|
|
|
|
|
|
п=1 |
|
|
|
|
N |
|
|
|
N |
|
|
|
|
N |
|
Яа@ — ~£Г + |
^ |
9а/3> |
ааЗ = |
^ |
ааЗ> |
а33 = |
azGc + ^ а ^ , |
|||||
|
|
п=1 |
|
|
п=1 |
|
|
|
|
п= 1 |
||
|
Л |
|
|
|
N |
|
|
|
|
N |
|
|
|
_ CLZ & CE C . \ ' п |
®3© |
_ \ ' |
|
||||||||
|
^а© — «_ |
“Г / уДд©> |
/ |
^3©> |
|
|||||||
|
|
|
|
|
71= 1 |
|
|
|
71=1 |
|
||
|
|
|
/4 |
| |
о / дп |
| |
о Л га |
\/2 |
/2 |
| лтг |
/4 |
|
а2а = ^ H I l L |
+ 2(^П22 + 2Л^ 12) ^ 71^71 + ^2222^71- |
|||||||||||
аа/3 = |
(^П 11 + -^2222 —4Л 1212) ^ 71^71 + ^П22^атг + ^тг)> |
|||||||||||
ааЗ = [^ПП^атг —-^2222^71 —(^П22 + ^Л ^гХ ^тг —^п)] ^атХ/Зп. |
||||||||||||
а33 = |
И м и + -^2222 —2Л П22)^1тг^2тг + ^?212(^тг —^2TI)2’ |
|||||||||||
а а© = |
^ lW a n + |
^220^/Зтг> |
|
а 3© = |
(^П © “ |
^22©)^п^2тг. |
||||||
^ п п = ^ |
|
|
|
+ 5Е сЕ 2 (и>п » : + |
u j y c) 2X n l ] , |
|||||||
Л7212 = 5u nz E cE |
^ ~ \ |
Л? 122 |
= 5и>” Е сЕ : ( ш п и2 + Ш |
» 2* " 1, |
||||||||
Л 2222 = ^ 2 ЕсЕаХп 1> |
9аа = 9lVcm + Я22^п’ 9а/3 = (?П —Q22)lanlpn, |
|||||||||||
я?1 |
|
= |
2 <Ч” |
( 1 + */с ) ( 1 + |
Я22 = |
<4 % |
|
|||||
|
Xn = |
w "[l - |
K ) 2]S C + |
а£[1 - |
(»/с)2]Я», |
|
Т}п = Шп ( 1 + I'с)Е™ + |
( 1 + I 'а )Е с, |
28 Гл. 1. Структурные модели композиционного материала
|
|
Сп = 2 К ( 1 + К ) Е с + ^ ( 1 + 1' Ж ] . |
|
||||
|
|
|
— |
az |
- |
у |
|
|
|
|
&z — |
~ Го I &Z — |
1 |
|
|
|
|
|
1 - |
(ис) |
|
|
|
|
|
|
N |
|
|
|
|
|
|
u z = Y ^ u z> |
l\n — cos фпу |
hn = |
sirn/'n, |
(1.15) |
|
|
|
n = 1 |
|
|
|
|
|
Здесь |
E ц, |
— |
модуль |
Юнга и коэффициент Пуассона матери |
|||
ала п-го |
семейства |
арматуры; и п , и™ — интенсивности |
армирования |
в поверхности и в направлении толщины оболочки для 7г-го семейства арматуры; и™ = 1 —и п для МДВ и МОВУ, а для МОВ и НМ и™ = 0; фп — угол укладки 7г-го семейства арматуры.
Соотношения МДВ [220] получаются при <5=1; моделям с одно мерными волокнами [218] соответствует <5= 0; формулы для НМ [260] получаем при <5= 0, Ес = 0.
После того как определены средние напряжения и деформации представительного элемента можно найти напряжения и деформации в структурных элементах КМ, а именно в связующем материале и ар матуре.
Для нитяной модели напряжения и деформации арматуры 7г-го слоя будут иметь вид [260]:
а (П)а - Е а е™и)а' е(1\)а ~ е 11^ln + ^ X i U n k n + б2 2 ^2п- |
(1 -1 6 ) |
Для модели с одномерными волокнами к соотношениям (1.16) до бавляются выражения для напряжений в прослойках связующего [235]:
Е с
<Т(аа)с2 - { |
^ |
'(бас* |
^ с ^ р р )у |
*Да/3)с2 — 2 Gc£<xf3y |
|
|
|
|
(1.17) |
Д<*3)с2 |
= |
|
/У(<*3)с2 |
'YCLZ/ ^ Z - |
Для уточненной модели с одномерными волокнами к напряжениям и деформациям, определяемым соотношениями (1.16), (1.17), нужно добавить выражения для напряжений и деформаций связующего в ар мирующем слое
^(пЭы = ^'се(п)с1> e(ii)ci = en*in + Zenhuhn + е22^п- |
(1-18) |
Невыписанные компоненты напряжений для волокон и связующего для МОВ и НМ полагаются равными нулю.
Для МДВ напряжения в арматуре, в связующем армирующего слоя и в связующем между армирующими слоями будут иметь вид [9]:
Е с
б(а/3)с2 |
— б<*/3> *Д аа)с2 |
(ис) |
(б (аа)с2 “Ь ^/с б (^ ^ )с2)> |
||
|
|
1 |
|
|
|
*Да/3)с2 |
2G -c6(a/3)c2i |
Д аЗ)с2 |
7"аЗ> |
Т(<*3)с2 |
G c Д а 3)с2> |
е(П)с1 — е0 !)“ _ |
Л7 |
п |
_ . п |
_ Л' |
|
в11’ |
^(1з)с1 |
— 7(°13)а — 713 |
1.3. Критерии прочности и начального разрушения |
29 |
(1.19)
а
( 1.20)
Здесь Л, р = 1,2 и используется правило суммирования по немым индексам.
1.3. Критерии прочности и начального разрушения полиармированного слоя
Проблемы прочности КМ разрабатывались многими авторами и по лучили в литературе широкое освещение. Как и при построении моде лей КМ, при установлении критериев прочности можно выделить два основных похода — феноменологический и структурный. В рамках пер вого из них КМ рассматривается как квазиоднородная упругая среда, для которой постулируется критерий прочности. Параметры, входящие в его математическую формулировку, определяются из эксперименталь ных данных. Среди феноменологических критериев прочности важное место занимает тензорно-полиномиальный критерий [205], который обобщает практически все известные феноменологические критерии. Однако следует отметить, что даже для относительно простых ви дов напряженного состояния требуется реализовать весьма трудоем кие программы экспериментов и математической обработки получен ных данных. Другой недостаток таких критериев — их формулировка в терминах средних напряжений, что не позволяет выявить механизм возникновения начального разрушения и исследовать направление его дальнейшего развития.
Структурный подход свободен от указанных недостатков. Это на правление базируется на изучении напряжений в элементах субструк туры, для каждого из которых принимается тот или иной критерий
30 Гл. /. Структурные модели композиционного материала
прочности. После определения средних характеристик НДС в кон струкции напряжения в элементах КМ восстанавливаются с помощью уравнений структурной модели. Таким путем вычисляются разруша ющие интенсивности внешних нагрузок всех элементов композита, и наименьшая из них принимается за нагрузку начального разрушения. Этот подход позволяет выявить эффективность работы связующего и армирующих элементов, указать рациональные по прочности пара метры армирования.
Рассмотрим структурный критерий прочности волокнистого КМ [9], который будет использован ниже в конкретных расчетах. Наряду с до пущениями, сформулированными в предыдущем разделе, примем сле дующие два постулата:
—адгезионная прочность связующего не ниже когезионной;
—материалы связующего и армирующих элементов подчиняются условиям прочности вида
^(<Да0). Да3).СГ+.СГ_) = 1, |
(1.21) |
где <Т+, а~ — пределы прочности материала на растяжение и на сжатие соответственно.
В качестве условий прочности можно принять, например, критерий Баландина [120], тогда
F = FB — [o'и + а %1 + ^ЗЗ ~ Оц 022 —ОцСГзз —022033 +
+ Ъо\2 + Зг12з + Зг23 - (а+ - о ~ ){о ц + сг22 + СГ33)] (сг+сг- ) - 1 . (1.22)
Если пределы прочности на растяжение и сжатие совпадают о + = о~ —= а*, то из (1.22) получаем критерий Мизеса
FM{o(al3), T(Q3), O-*) = 1. |
(1.23) |
При использовании критерия Мизеса удобно использовать приве денные интенсивности напряжений в матрице и волокнах n -го семей ства арматуры
bsr = max ISUp FM {Ofa0)cV r (Q3)cl >o£),
^Slip F M (^(Q|3)C2> ДаЗ)с2>0’*) >>
К = ^/sup FM{o?a(3)a, r ("a3)a, a*n) , |
(1.24) |
где V — область занятая слоем. Начальное разрушение в матрице или волокнах произойдет в случае, когда bsc, bs™ > = 1 соответственно.
Основываясь на приведенных зависимостях, линейная задача опре деления нагрузки начального разрушения армированного слоя реша
1.4. Сравнительный анализ расчетных характеристик |
31 |
ется следующим образом. Пусть рассматриваемая оболочка нагружена системой внешних нагрузок, интенсивности которых пропорциональны скалярному параметру Р. В силу линейности дифференциальных урав нений и граничных условий соответствующей краевой задачи статики оболочечных конструкций напряжения в волокнах и связующем будут пропорциональны Р [9]. Тогда нагрузки начального разрушения обрат нопропорциональны интенсивностям напряжений в соответствующих элементах и равны для связующего, волокон n -го семейства и слоя в целом
Р : = Р Ы~ с '. Р Г = Р-(Ь~8Па Г ' , P ' = min{Pc*,P„*"}. (1.25)
где величины bs_C} bsa соответствуют НДС при значении скалярного параметра Р = Р.
Для определения нагрузки начального разрушения в нелинейном случае можно использовать следующий итерационный процесс.
1.Задается начальное значение нагрузки Р '.
2.Решается задача определения НДС оболочки (пластины) при значении скалярного параметра Р = Р '.
3.Вычисляется bs* = max{6s£, bsc} и Р* по приведенным выше формулам.
4. Если |1 - 6s* | > е > 0, то задается новое значение величины Р ' =
= Р* и осуществляется возврат к пункту 2.
5.Искомое значение нагрузки начального разрушения равно Р*.
1.4. Сравнительный анализ расчетных характеристик композиционных материалов с экспериментальными данными
Вопрос о сравнении эффективных упругих и прочностных характе ристик армированных композитов, рассчитанных на основе различных моделей армированной среды, с данными, полученными из экспери ментов, весьма актуален. В литературе в основном представлены экс периментальные данные по продольному и поперечному растяжениям и сдвиговой деформации в плоскости слоя однонаправленно армиро ванного образца.
В табл. 1.1, 1.2 приведены эффективные модули Е ь Е% и G n = G |
|
для однонаправленных композиционных материалов, полученные экс |
|
периментально [185, 345], по формулам (1.2) и МДВ. |
Здесь 5 — |
относительная разность между величинами, полученными |
эксперимен |
тально и по моделям [9, 32].
На |
рис. 1.3 приведены результаты расчетов эффективных модулей |
E i, Е 2 |
однонаправленно армированного бороалюминиевого слоя в за |
висимости от удельного содержания арматуры с применением МДВ (сплошные кривые) и МОВ при az = 1 —ша (пунктирные кривые) и при az — 1 (штриховые кривые). Символами А на графике обозначены соот-