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

Статистическая механика композитных материалов

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

Соотношения (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]

я

 

Ьш = Р [(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 — с вероятно­ стью 1v. Величина 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=

 

 

 

 

 

 

 

 

 

 

 

 

о

о

**

 

 

о 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 =

________________ .

 

 

 

 

 

 

 

С п и С3333—

C i из

 

 

 

 

 

Ci 133 С2222 — Ci 122 С2233

 

(3.35)

 

 

«13 =

 

Спи C2222 — C1122

 

 

 

 

 

 

 

 

 

 

Е3 =

 

 

 

 

 

 

 

 

 

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 »

где

— Ош 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