Статистическая механика композитных материалов
..pdfСоотношения (3.30) имеют место для стеклопластиков при Р > 0,22. Окончательно формулы для макроскопиче ских модулей упругости однонаправленных материалов запишем в виде
|
q(l* + |
G*)2 |
|
С°Ш1 = 1 + 2т + Р { \ - Р ) |
|
|
|
|
\ — q(l* + G*){\— 2P) + |
||
+ |
gl*2 |
|
|
\ — gG*(l— 2P) |
|
|
|
C\ 122 = l + P (1 |
q (/* + G*)2 |
||
— P) |
|
|
|
|
1 — q(l* + G*) (1 — 2P) |
||
_________ gG*2 |
1 |
(3.31) |
|
|
\ — g G * { \ — 2P) |
J |
|
|
|
||
, , |
P (l - P) ql* (l* + G*) |
||
U 1 1 33 = / - 4 - -------------------------------------------------------, |
|||
|
l — q(l*+ G*){\— 2P) |
||
G*2323 — |
P { l —P)fG** |
|
|
ТП -\~ 1 — /G* (1 — 2P) |
’ |
Сзззз = / + 2m +
P(1 — P) ql*2
1 _ 9 (/*_PG*)(1— 2P)
На рис. 18 приведены результаты вычисления по фор мулам (3.31) и (3.27) макроскопических модулей упру
гости однонаправленного |
стеклопластика |
(Е1 = 7- |
1010 Н/м2; £И = 0,3-1010 Н/м*; |
^ = 0,2; дп = 0,35). Для |
модуля Е 1 (перпендикулярно волокнам) расчет по этим формулам дает завышенные (по сравнению с эксперимен тальными) значения. Отчасти это обусловлено тем, что в расчете не учитываются такие неизбежные дефекты ма териалов, как пористость и отклонения волокон, от задан ного направления армирования. Влияние этих факторов обсуждается ниже (п. 3 гл. 3 и п. 2 гл. 5). Область эффек тивного применения формул (3.31) ограничена компози тами, составленными из компонентов, модули Юнга ко торых отличаются не более чем на порядок.
Для исследования зависимости макроскопических модулей упругости однонаправленного композита от фор мы сечения и взаимного расположения волокон арматуры
ill
применим решение статистической краевой задачи теории упругости в реализациях (см. п. 6 гл. 2). В корреляцион ном приближении^ получаем формулу, аналогичную
(3.13). Интегралы I*s), |
входящие в эту формулу соглас |
но (3.11), вычисляются |
по контурам L.<?, ограничиваю |
щим сечение волокон: |
|
"l[s) = J |
def G(x/? х') n(sWLs, |
где n(s) — орт внешней нормали к контуру Ls; G(x, х') — |
тензор Грина (2.60).
Пусть волокна арматуры имеют круглое сечение^По мещаем точку 7И(х) в центр круга (радиус круга может быть произвольным). Тогда n(0)= (fi) и I(0)= (Iijmn), при чем величины Iijmn определяются соотношениями (3.25), т. е. поправки, полученные путем решения задачи в реа
лизациях, в частном случае совпадают с hjjmn из (3.26). В общем случае они содержат слагаемые, зависящие от формы сечений и взаимного расположения волокон. Ана лиз формул (3.13) показывает, что влияние отмеченных
факторов на макроскопические модули упругости |
мате |
|||
риала невелико. |
армирующего' вещества |
|||
Влияние |
неоднородности |
|||
можно оценить по формулам |
(3.14), подставив 1^тп из |
|||
соотношений (3!25). Для построения последующих |
при |
|||
ближений |
используется алгоритм, |
изложенный в п. 1 |
||
гл. 3 применительно к хаотически |
армированному |
ма |
териалу.
Однонаправленный материал с анизотропной армату рой. Более сложную проблему представляет вычисление модулей упругости однонаправленного материала, арми рующие элементы в котором анизотропны. Как следует из результатов п. 5 гл. 1, даже при кубической или гекса гональной. симметрии армирующих элементов тензор средних модулей упругости С= С1Р + СП(1—Р) не явля ется изотропным (в простейшем случае С — трансвер сально-изотропный тензор). Поэтому при вычислении поправок необходимо использовать тензор Грина для соответствующей анизотропной среды (с модулями С). Аналогичная задача возникает при вычислении макроско пических модулей упругости композитных материалов с текстурой вращения произвольного вида [82]. Средние
112
модули упругости таких материалов вычисляются на основании результатов п. 5 гл. 1; они образуют транс версально-изотропный тензор. Представив тензор Грина для неограниченной трансверсально-изотропной среды
[93] в виде G(x, х') = ” Н(ф, -&), после интегрирования по
г имеем [82] |
2я |
я |
|
||
Ьш = Р [(1 — Р) С*С*' + D<!>] J |
f Н (ф, ft) <Шср. |
0J О
Интегралы по Ф, входящие в эту формулу, не выража ются, через элементарные функции, поэтому дальнейшие вычисления производятся численным интегрированием на ЭВМ.
К другим эффектам, учитываемым с помощью мето дов статистической механики, относится влияние на макроскопические свойства однонаправленного компози та пористости (см. п. 2 гл. 5), искривлений [120, 121], длины волокон [122, 123], их разориентации [124] и не линейного поведения связующего [113—117].
3. МОДУЛИ УПРУГОСТИ слоистых МАТЕРИАЛОВ
Метод вычисления. Рассмотрим композитный матери ал, состоящий из анизотропных слоев с детерминирован
ными модулями упругости Cijmm заданными в системе
координат (хь *2 , *з), связанной со слоем (рис. 19). Полагаем, что толщина слоев мала по сравнению с тол щиной тела. Если каждый слой изготовлен из однона правленно армированного композитного материала, то толщину слоя считаем большой по сравнению с характер ным размером армирующих элементов (например, диа метром волокон). Это обеспечивает детерминированность
модулей Cijmtl. При случайном расположении слоев мо дули упругости материала, отнесенные к окрестности точ ки радиусом порядка толщины слоя, в системе координат *2, *3, связанной с телом, являются случайными функ
циями и задаются формулой преобразования (1.23). Для вычисления макроскопических модулей упругости
воспользуемся результатами И. М. Лифшица и Л. Н. Ро-
8. Зак. 674 |
113 |
зенцвейга [53] применительно к одномерному случаю. Пусть плоская плита, неограниченная в плоскости х2х3, имеет конечную толщину в направлении Х\ и состоит из большого числа однородных тонких слоев, причем край ние слои одинаковы. На обе поверхности плиты действу ют равномерно распределенные силы.
Обозначим 0 цтп — компоненты |
случайных |
модулей |
|
упругости слоя в системе координат |
(х и х2, х3), |
связан |
|
ной с телом; Д — определитель, составленный |
из вели |
||
чин Qijmn' A = |20ifij— |
Аа — алгебраическое до |
полнение определителя |
Д; Z)= |<Дt- ,/Д>|; |
D\j — алге |
браическое дополнение определителя D. |
проекции на |
|
Используя уравнения |
равновесия в |
ось. Хи уравнения совместности деформаций и физические уравнения, записываем компоненты тензора макроскопи ческих модулей упругости сплошного композитного ма
териала |
[53]: |
|
|
|
|
|
C' n / = _ ^ D '(1+6li6l' ): |
|
|||
|
Сим = |
|
; |
(3-32) |
|
c ;4rs = |
^ < ( 2 0 |
p, lv- 6 lv0 P9ll) ^ v _ \ / - 0 16rs |
|||
- |
< 2 0 Jlglv - |
6lv0 Pell) 0 1Brs A |
v . ) + < @pqra > , |
||
где i, j, |
a, 0, y, |
6=1, |
2, 3; p, q, |
r, s = 2, |
3. |
В случае ортотропных слоев вывод формул для макро- |
|||||
скопических модулей |
упругости |
слоистых |
композитов |
||
может быть упрощен [1.6, 125]. |
|
|
114
По формулам (3.32) можно вычислить макроскопиче ские модули упругости слоистого материала произволь ной структуры при заданном расположении слоев по на правлениям (распределении косинусов Я1;). Если слои из однонаправленно армированного материала, их модули
упругости Сijnui вычисляются по формулам п. 2 гл. 3.
В общем случае слоистой структуры матрица косину сов Xij имеет вид
/ 1 |
0 |
0 |
|
(Хи) = [ |
0 |
C0S(P |
— sin ф |
\ |
0 |
sin ф |
cos.cp/ |
где ср — случайный угол поворота слоя относительно си стемы-координат, связанной с телом (хи х2, х3). Для не которых конкретных видов структуры плотности распре деления величин Хц приведены в работах [125—127].
Ортогональное армирование. Пусть однонаправленные
слои стеклопластика с модулями упругости С//тя (ось х направлена вдоль волокон) расположены параллель но плоскости х2х3 так, что направление волокон (х3) в
слое параллельно х2 с вероятностью v и х3 — с вероятно стью 1—v. Величина v численно равна степени армиро вания — объемной доле слоев, волокна в которых парал лельны х2. Так как направленность волокон слоя вдоль х2— случайное событие,-поставим ему в соответ
ствие индикатор |
Я(Я=1, если х 3 ||х2, и Я = 0, если х3 ||х3). |
||||
Матрица направляющих косинусов слоев имеет вид |
|||||
|
|
1 |
0 |
0 |
|
|
(Я..) = |
0 |
X |
Я—1 |
|
|
\ |
0 |
1—Х X |
||
Получаем: 0 Ш 1 |
== Сцц; |
61212 = |
С1212А,2 + С*iз1з (1 Я)2; |
||
01313 — ^*131 з ^ 2 + |
£ 4 2 1 2 (1 — |
А,)2; |
Д = 4 С |
1111 © 121201з1з’» |
|
^4©1212©1313 |
0 |
|
0 |
||
(*„) = |
о |
2Сц1101з1з |
0 |
||
|
о* |
|
0 |
|
2Сц ц 0 1212. |
8* |
115 |
|
Сип |
0 |
О |
|
|
|
|
|
|
D = |
О |
/ 1 А |
о |
|
|
|
^ 201212 / |
|
|
|
О |
О |
/ _ ! |
А |
|
|
|
\ 2©1313 |
/ |
Далее, учитывая, что для задания упругих свойств ортотропного материала нужно знать 9 независимых по стоянных, по формулам (3.32) находим:
|
Спи = Сип; |
СU22 = |
< ©1122 ) '» |
||||||||
Cl 133 = < ©1133 ) ; |
|
С2323= |
< ©2323 > » |
||||||||
' 1212 |
/ _ |
1 |
А-1 |
|
М313 = / - |
|
1 |
А-1 |
|||
|
\ 0 f |
/ |
|
|
\ |
|
0, |
/ |
|||
|
' |
|
|
' |
|
||||||
|
'1212 |
|
|
|
^1313 |
(3.33) |
|||||
|
|
|
|
|
|
|
|
|
|
||
С2233= ( ©2233 ) |
^1111 ( ( ©1122©1133 ) |
|
|
( ®1122 ) ( ©1133 ) )» |
|||||||
С2222= < ©2222 ) — С'п\\ ( ( ®1 122 > — < ©1122 ) 2)'у |
|||||||||||
Сзззз = |
< ©зззз ) — Cj7n ( < ©1133 > — < ©изз ) 2)- |
||||||||||
Заметим, что |
< |
1 |
> - ]ф<® \2\2>- |
В Ы Ч И С Л Я Я |
|||||||
01212 |
|||||||||||
|
|
|
в |
формулах |
(3.38), оконча |
||||||
математические ожидания |
|||||||||||
тельно получаем следующие значения постоянных: |
|||||||||||
Cl И 1= |
Спи; |
Cl 122 = |
Cl 122 v-\-C\\ зз (1 |
— v)\ |
|||||||
Ci i33 = |
Ci i3 3 1 ; -(- Ci i22 (1 |
— v)\ |
C2 3 2 3 '= |
Ci3 i3 *, |
|||||||
|
Cl212 = |
Ci212 Cj3i3 |
|
|
|
|
|||||
|
Cl313^ + |
Cl 2 12 (1—v) |
|
|
|||||||
|
|
|
|
(3.34) |
|||||||
|
|
|
C1212 Ciзi3 |
|
|
|
|||||
|
Ci313 = |
|
|
» |
|||||||
|
' 1313 |
— |
---------- тг?---- ------ : |
||||||||
|
|
|
Ci2 1 2 ^ + Ci зiз (1 —v) |
|
|
||||||
C2222 = Cl 111 v 4" |
C3333 (1 — v) — ^2233» |
116
C2233 |
Сзззз = C3 3 3 3 D + |
С п и |
(1 — у) — h2233> |
|||||||||
— Cl 133 + Л23', |
|
^2233 = |
CJj J |
J |
(Cl 122- |
С ,,зз)М1— V)' |
||||||
Как видно из |
формул |
|
(3.34), модули С2222, С3333 и |
|||||||||
С°22зз |
отличаются от средних значений, найденных в ра |
|||||||||||
боте |
[126], |
однако вычисления |
|
показывают, |
что для |
|||||||
стеклопластиков поправки Л2233 малы. |
|
|
||||||||||
Технические постоянные находятся по формулам: |
||||||||||||
|
|
|
Е i= |
|
|
|
D° |
|
|
|
|
|
|
|
|
|
о |
о |
** |
|
|
о 2 |
|
|
|
|
|
|
|
G 2 2 2 2 G 3 3 3 3 — |
С2233 |
|
|
|||||
|
|
|
|
О |
о |
|
|
|
о |
о |
|
|
|
|
|
C 1 1 2 2 G 3 3 3 3 — G 1 1 3 3 C 2 2 3 3 |
|
|
|||||||
|
|
« 1 2 |
= |
|
<1111 ^3333 |
— с 1133 |
|
|
||||
|
|
|
|
|
|
|
||||||
|
|
|
Е2 = |
________1Г________ . |
|
|
||||||
|
|
|
|
|
С п и С3333— |
C i из |
|
|
||||
|
|
|
Ci 133 С2222 — Ci 122 С2233 |
|
(3.35) |
|||||||
|
|
«13 = |
|
Спи C2222 — C1122 |
|
|||||||
|
|
|
|
|
|
|
||||||
|
|
Е3 = |
|
|
D° |
|
|
|
|
|
||
|
|
C1111 C2222 |
|
Cn 22 |
|
|
||||||
|
|
|
|
|
|
|
||||||
|
|
|
Cl 111 C-2233 -- C] 122 Cl 133 |
|
|
|||||||
|
|
n 32 |
~~ |
|
С ип C3333 |
— C]233 |
|
|
||||
|
|
|
|
|
|
|
||||||
|
О |
|
О |
|
О |
|
О |
|
|
о |
о |
|
|
G12 |
= С 1 2 1 2 ; |
|
G i 3 = C i 3 i 3 ; |
G23 = |
0 2 3 2 3 » |
||||||
где |
D° — Ош 1C2222G3333 -f- 2 Ci 122G1133C2233 — |
|||||||||||
|
||||||||||||
|
— |
С ц ц С 2 23 3 |
— G 2222G П 3 3 — C 3 3 3 3 C 1 122- |
|||||||||
В качестве иллюстрации |
|
на |
рис. 20, а, б приведены |
|||||||||
зависимости модуля |
|
Юнга |
|
и |
коэффициента Пуассона |
|||||||
ортотропного |
стеклопластика (для главных |
направле |
ний) от степени армирования v. Экспериментальные зна
чения модуля Юнга (о) и коэффициента Пуассона |
(о) с |
|
9 5 %-ныМИ доверительными |
границами указаны |
для |
АГ-4 С согласно [128], а для |
27-63С — согласно |
[129]. |
117
Рис. 20. Зависимость модуля Юнга (1) и коэффициента Пуассона (2) ортотропных стеклопластиков АГ-4С (а) и 27-63С (б) от степени армирования без отклонении слоев (сплошные линшг) и при среднем квадратическом угле отклонения 5ф=0,1 (пунктирные линии)
Отличие экспериментальных значений от вычисленных можно объяснить отклонениями направлений слоев от заданных при укладке в пресс-форму и прессовании.
Влияние случайных отклонений. Пусть отклонения слоев от заданных направлений малы, поэтому можно принять
costp = |
(( 1 |
с вероятностью о; |
|
|
( ср |
с вероятностью (1 — у); |
|
sin ф = |
( ф с вероятностью у; |
|
|
{ |
с вероятностью (1 --и); |
|
|
|
I 1 |
|
|
< cos2 ф > = у + |
(1 |
< sin2 ф > = 1 — у + |
vDф; |
|
|
< cos4 ф > = у; |
|
< sin4 Ф > = 1 —v, |
< sin2 фcos2 ф > = Dy, |
DV |
( A P — дисперсия ф). В этом случае, вычисляя по форму лам (3.33), получаем следующие зависимости модулей упругости ортотропного материала от Dv:
С\ 111 = Cl 11 Г. С2222 = Cl п ! v -)- Сзззз (1—V)--
— 2 [С] 111 v + Сзззз (1 у) — Cj 1зз — 2Сiзiз] Цр —
118
К 2и (1 - |
2Цр) + |
2d*C;i33 [i> + |
(1 - |
2v) Оф]+ |
|||||||
с 1111 |
|
|
|
|
|
|
|
|
|
|
|
+ С*пзз — |
C \ \ 22}; |
d* = |
Ci 122 — |
С ц 3з; |
|
||||||
C1212 = |
|
C1212 С т з |
|
|
|
|
|
||||
C1212 + f*[v + (1—2v) Dy] ’ |
|
|
|||||||||
|
|
|
|
||||||||
C1313 |
= |
|
C1212 Сшз |
|
|
|
|
|
|||
C'\3i3 — f*[v + (1 — 2v)Dq>] ’ |
|
|
|||||||||
|
|
|
|
||||||||
f* = C\3\3 — C1212; |
Ci 122 = |
Сц33 + |
d* [y+ (l —2y) Dj; |
||||||||
Ci i 33 — C i i33 — d* [v + |
(1 — 2u) Оф]; |
(3.36) |
|||||||||
C2233 = Ci i33 -f (Ci 111 + |
C3333 — 2Ci i33) |
||||||||||
+ |
|
||||||||||
H— —;------[C1122 C11 з3— |
C 1122 C1133 — d*2D^]\ |
|
|||||||||
C i 111 |
|
|
|
|
|
|
|
|
|
|
|
C2323 = |
C1313 -f- (Ci 111 -f- C3333 — 2C1133 — |
|
|||||||||
— 2C1з13) Оф; |
C3333 = C3^33 ti + |
C jin (l — u) — |
|
||||||||
— 2 [C3333 у + Ci 111 (1 — v) — C1133 — 2C1313] Цр-(- |
|
||||||||||
-\— |
----- {Ci 133 — C1122 + |
2d*C1122 X |
|
||||||||
Ci 111 |
|
|
|
|
|
|
|
|
|||
X [V + |
( 1 |
- 2v) D J |
- d**v( 1 - |
2Ц р)}. |
|
|
|||||
При Дф=0, как нетрудно убедиться, |
получаем соотно |
||||||||||
шения (3.34). |
|
|
|
|
отклонений |
слоев |
на |
||||
На рис. 20 показано влияние |
|||||||||||
упругие постоянные ортотропных стеклопластиков. |
|
||||||||||
Полагая в формулах (3.36) и=1, находим соотноше |
|||||||||||
ния для однонаправленного |
(в среднем) |
материала, |
ко |
||||||||
торые позволяют |
учесть влияние |
отклонений слоев на |
|||||||||
упругие свойства. |
Это влияние для стеклопластика, со |
||||||||||
стоящего из компонентов со свойствами: Е 1= 7 |
1010 Н/м2; |
||||||||||
£ п = 0,3- Ю10 Н/м2*; |
/г1 = 0,2; |
дп = 0,3 при Р = 0,46 и пори |
стости 5%, показано на рис. 21.
Равномерное армирование в плоскости (текстура вра щения). Этот случай соответствует равномерному рас-
119
Рис. 21. Зависимость постоян ных упругости однонаправленного стеклопластика слоистой структу ры от среднего квадратического угла отклонения слоев stp:
£ 3; 2—п2 з ' 3~ Е \ ; 4~ °2 Ъ
0,05 0,10 0,15 Sep
пределению угла <р по направлениям в плоскости армиро вания. При вычислении моментов распределения, входя щих в формулы (3.32), воспользуемся результатами п. *5 гл. 1. Тогда для макроскопических модулей упругости рассматриваемого композита получаем [12]:
'in i = СшГ, СП22 = Ci 1зз = — (Cj 122 + Сцзз);
Ci212 = Сiзiз = 2 Сi2 1гС1313 (Ci212 + С1 2 12) г\
С2233 = — (Ci 111 + С3333) ----С1313 +
+-7-Ci 133 — (Ci 122 —Сцзз) 2 (8 С1111)"1; 4
С°2323 = "Г" С1з1з -1— — (СJill + С3333 — 2Сцзз);
2 о
С2222 — Сзззз = С2233 + 2С2323.
I ^ налогичнь,й случай рассматривался в работах [125 130], причем в работе [130] найдены лишь средние зна
чения упругих податливостей. |
|
v |
Переход к техническим постоянным |
осуществляется |
|
ПО формулам (3.35). |
^ |
иилиеии |
В отношении упругих свойств равномерное в плоско-
120