Прогнозирование влажностно-упругих и термоупругих свойств древесно-минерального композита

ДЕРЕВООБРАБОТКА

ПРОГНОЗИРОВАНИЕ ВЛАЖНОСТНО-УПРУГИХ И ТЕРМОУПРУГИХ СВОЙСТВ ДРЕВЕСНО-МИНЕРАЛЬНОГО КОМПОЗИТА

В.И. ЗАПРУДНОВ, проф., зав. каф. геодезии и строительного делаМГУЛ, д-р техн. наук, А.С. ЩЕРБАКОВ, проф., зав. каф. безопасности жизнедеятельности МГУЛ, д-р техн. наук

В ряде работ, опубликованных нами ранее [1-4], мы рассматривали методы прогнозирования жесткостных и прочностных характеристик древесно-минерального композита (арболита), основанные на применении модели стохастических неоднородных сред.

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

При эксплуатации конструкций из древесно-минерального композита за счет его гигроскопичности под влиянием климатических факторов происходят значительные колебания влажности материала, вызывающие относительные деформации порядка 0,1 - 0,5 %. Это приводит к появлению напряжений, особенно в многослойных панелях, снижающих долговечность и прочность сооружений из древесно-минерального композита.

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

zaprudnov@mgul.ac.ru; scherb@mgul.ac.ru

Щ = nj

Закон связи между относительным изменением влажности AW и деформациями имеет вид

j=VW

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

j = Щ9,

где a - тензор коэффициентов температур-

Jk

ного расширения.

В изотропном случае

ajk = aj

Учитывая соотношения (2) и (3), зависимости между тензором напряжений oj.k, тензором деформаций s приращением тем-

k

пературы 9 и относительным изменением влажности AW представим в форме

°jk = kjkngSmg - в*9 -

где

РЛл YjAW = kA

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

pjk =p5j*, p = (3k + 2p)a.

Y* =§ jk ^, Y =

Тензоры макроскопических постоянных определяются соотношениям^

°jk = Kkng Sng - e*jk 9 - Y*jkAW,

где черта сверху означает осреднение по макрообъему.

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

Ы = ХУ A)-e,k <9-Y*k

Тем самым задача определения макроскопических тензоров kjkng, рjk, уjk оказывается эквивалентной определению закона

ЛЕСНОЙ ВЕСТНИК 3/2012

125

ДЕРЕВООБРАБОТКА

связи между математическими ожиданиями

Ы-Ы • w .< aw •

Подставляя (5) в уравнения равновесия без учета массовых сил, получаем

Дополним уравнение (10) условием на бесконечности U0 (да) = 0, тогда для определения деформаций арболита получим уравнение

«У=1 ,(- *(2>)«

KAn^-M-YngAW)*.,

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

«(1) = (е) + K(U) * (A^ -р0- yAW)(2) (12) интегральный оператор К действует по правилу

кS* Ф =J Uj) g (x - x<2V2’ dv'2 +

+JUn(jk)(x<1> -x<2>y2>cos(n,xg2VS'<2>).

Предположим, что древесно-минеральный композит составлен из двух компонентов: связующего цемента и органического заполнителя. Осредним уравнение (12) при условии, что в точке x(2) находится материал заполнителя. Для этого умножим (12) на условную плотность распределения

f

и выполним интегрирование по A, е, Р, у. Здесь мы предполагаем также, что древесно-минеральный композит находится под воздействием равномерной температуры и равномерно увлажнен, т.е.

(в) = 0 = const, ^AW) = AW = const .(14) После интегрирования получаем интегральное уравнение для определения условных моментных функций тензора деформаций

е f = (е) + K(1,2) * (A0fе ffPff + A0Vf pfc -

A0 Pf-V0 Pfc-r'AWf-yc ШрГс .

Используемые здесь обозначения эквивалентны введенным в уравнение

sc = (г) + K(1,2) * (A0VcPcc + A0 f е fcPf ) .

Для условных вероятностей Pf, P^ имеют место, очевидно, соотношения

Pff + Pfc = 1, 5 fPfc = 5cPf, Pfc (0) = 0;

Pfc И = 5f, Pff (0) = 1, Pff M = 5 f . (16) Используя равенства (16) и ограничиваясь приближением однородности деформи-

рованного состояния связующего и заполнителя есо « ec = const, еfc «е* = const, из (15) получаем

еf =

Осредняя выражение (5), находим соотношение, связывающее средние напряжения древесно-минерального композита со средними деформациями, средними деформациями заполнителя, температурой и влажностью

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

Уравнение (17) получено без каких-либо предположений относительно вида микроструктуры древесно-минерально-

го композита. Ранее при интегрировании в ec =(«) + K(1,2) * Pcc (A0cec - A0f еf ) было принято, что плотность Pcc зависит лишь от расстояния между точками. Это соответствует хаотическому расположению частиц заполнителя в пространстве.

Проводим построение функции Pff для случая, когда связующее древесно-минерального композита содержит расположенные, но ориентированные в одном направлении волокна заполнителя, имеющие форму эллипсоидов вращения. Для определенности предположим, что волокна заполнителя ориентированы вдоль оси Ox3 , тогда древесно-минеральный композит будет статистически изотропен в плоскости x x а значит,

Pff

Вследствие статистической изотропии свойств в плоскости xx , x2 достаточно построить функцию Pff в произвольной плоскости, проходящей через ось x3. Проведем в этой плоскости некоторую прямую x3 =mr. Предположим, что длины отрезков прямой, образованных точками пересечения прямой с поверхностями раздела между связующим и заполнителем, распределены по экспоненциальным законам

126

ЛЕСНОЙ ВЕСТНИК 3/2012

ДЕРЕВООБРАБОТКА

Л/

Pc (T ) =

1 к Л

T /

—exp

ю / V ю / ,

1 к

— exp -

Юс V Юс у

где Tf - длина отрезка прямой, проходящего через волокно заполнителя; тс - длина отрезка, проходящего через связующее.

Функцию PfX) на этой прямой можно рассматривать как условность перехода некоторой системы из /-состояния в /-состояние. При этом из соотношений (20) следует, что условная плотность P//(x) должна удовлетворять дифференциальному уравнению

dPff

d Ixl

f

+

1

- + -

1

Л

Vю f

ю„

Pff(|x) = —

ю

и условиям (16). Решение уравнения (21) имеет вид

к и Л

Pf (|xС = - +- exp

-сю

f у

где ю / - средняя длина хорды пересечения прямой x3 =mr с эллипсоидом вращения, которым моделируется волокно заполнителя.

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

п2 a1a2y[\+ m

2

ю / ='

b(a2 + m2al)1/2 ’ ^ ^

где a, а2 - размеры полуосей эллипсоида вращения соответственно в поперечном и продольном направлениях.

Подставляя в (23) m= x3/r и принимая во внимание выражение (22) для условной плотности P

Pf (x) = -/ + - exP

окончательно получаем

к „ п—Т\

8’ 2 ”2

r

x

е 2 \ 2 2 -c п \ ai a2

Используя выражение (24), можно выполнить интегрирование в уравнении (17) и найти явный вид тензора.

Пусть отношение продольного и поперечного размеров волокон заполнителя к = a2/ a1 > 1, т.е. древесно-минеральный композит содержит волокна заполнителя в форме вытянутых эллипсоидов вращения. В этом случае получаем пять независимых компонентов трансверсально изотропного тензора K * Pff :

I =

3 ^ к + м)

к 2

1 1б(м)(к + 2|м)[ (к2 -1)2 (к2 -1)5/2

P

12 =

2{ м

к1

- + -

к

2 1\3/2

к2 - 1

P

к2 (1 + 2к 2 ) - к(1 - 4к2 ) 1б(м)(к + 2|м)|_ (к2 -1)2 (к2 -1)5/2

2 к

P

= 3

Т3 =

14 =

4

3к2

(к2 -1)2 (к2 -1)5/2

P

2к2

+

3к

2(м)(к + 2|м)[(к2 -1)2 (к2 -1)5/2

P

+

+

I5 =

1 м к + м

1

+

к

P

к2 -1 (к2 -1)3/2

3к 2 к

2(м)(к + 2^[ (к2 -1)2 (к2 -1)5/2

P

+

+

4< М>

2 - к2

к

P

к2 -1 (к2 -1)3/2

P = 1w(k - V k -1), где в (25) обозначено

Т = к * P • Т =___K * P •

^ _Л1И1 12~г Л1122 А/’

-

Т = — K * P • Т = — к * P •

73 _ - Л1133 rffi А4~ - Л3333 rffi cc 1

Т = — K * P

*5 - - Л1313 // .

Если заполнитель имеет форму сплюснутых эллипсоидов вращения к = a2 / a1 < 1, формулы (25) принимают вид

к 2

т = 3( к + м) Т1 =

+

1б

2\5/2

б

+

+

Т 2 =

2< м> к + м

к2

+

к

б

1 - к (1 - к2)3/2

к2

1б( м)(к + 2м

+

2\5/2

б

1

1

1

ЛЕСНОЙ ВЕСТНИК 3/2012

127

ДЕРЕВООБРАБОТКА

I =

14 =

4

К + м

3k2

+ k(1 + 4k2) Q

2 ^ м) (К + 2м i Г i

(i-k2)2 (i -k2) 2+k2 3k

2\5/2

2\5/2

Q

k2

м)_1 - k2 (1 - k2 )3/2

Q

I5 =

К + м

3k2

+ k(1 + 2k2 ) Q

2

25/2

4

2 - k2

k

23/2

1 - k2

Q = arcsin-\A - k2

Q

Если k = a2/ a1 = 1, т.е. включения имеют сферическую форму, из (25) и (27) путем предельного перехода получаем

Л = 14 -

15

i

Используя формулы (25-28), из (17) находим средние деформации заполнителя ef и подставляем их значения в (18). В результате проделанных операций получаем следующие соотношения между средними напряжениями, деформациями, температурой и влажностью:

(° jk) = ж! (е pp) 5* + 2ж2 (е Jk) +

+жЦе^) 5jk -P1Q -rlAW,

(°33 ) = ж4 (еpp ) + (^ - ) (е3^ - Р2б - Y*2AW

Здесь макроскопические влажностнотермоупругие постоянные древесно-мине-

* * *

рального композита жi, р j, у j определяются формулами

ж* =Кс +4 f [а(К f -Кс)+а2(К f+2м f -Кс- 2мс)]; ж2 = мс +4 fP (м f-мс); ж3 = 4 f (ai- 2a2- pi )(м f-мс);

ж4 = Кс + 2м с + 4 f [Ь1 (К f -Кс ) +

+b2(Кf + 2мf -Кс - 2мс)];

ж55 = 2[мс +4 fp2 (м f-мс)];

Pi =(в)+4 f [d1 (К f +м f -Кс -мс )+d2 (К f -Кс )]; в2 =(в)+4 f [d1 (Кf -Кс )+d2 (К f + 2м f -Кс-2мс )]; Y* = М + 4f [d3 (Кf +мf -Кс -мс ) + d4 (Кf -Кс )];

y2 = (y)+4 f d (К f-Кс)+

+d4

где

a = Д {1 - К'(2/3 +14) - 2м'/4 - 2/21 --2m2 (I -13) + 2[Кт2 + 2м' (l2 + m2)]S};

а2 = Д [(К - /2 )(213 +14 ) +

+2(м -m2)I3 + 2(м l2 -К m2)S];

2

b1 =д[(К -l2)I+2(м -m2)I3 + 2(м l2-К m2)S];

b2 = д{1 - 2К 1 - 2м (I -13 ) -12 (2I3 +14 ) -

-2m214 + 2[2К м2 + м (l2 + 2m2 )S]} ;

Pi =1 -2”’,(1‘-14 di = 44с(Pf-Рс)(I-2n'S);

P2 =

1 - 2м

1 - 2m015. j = 1

"; " = д

д

----pi; d2 = т4с

1 - 2м h

Д = 1 -К (21 - 213 +14) --2м (I -13 +14) + 2м (3К + 2м )S ;

2

d = Д 4с (Y f-Y с)(I - 2м S);

Д

d4 = Д4с (y f -Yc )(2I3 +14- 2p's );

Д

I = /1 + 12 + /3 ; S = 114 - 2I32 - I314 ;

к'=

м = (4с-4 f )(м f-мс); /2 = К 2 -(К);

т2 =мс -

Технические упругие постоянные арболита: модули Юнга в продольном направлении E*’, в поперечном направлении £*, коэффициент Пуассона v* в плоскости, параллельной направлению волокон, модуль сдвига в той же плоскости ми модуль сдвига в поперечной плоскости м*, а также коэффициенты температурных и

* ^ i i i i

влажностных деформаций а1, а2, п1, П2 можно выразить через коэффициенты

Е*' = с; Е*' =

z

4ж2с

2ж2ж4 + с

i' i'

; v = —; м = ж;

z

м = ж2 ; а1 =

Pi i л i i л i — л i i

1ж4 - в2ж1 ^ в2z - 2в1ж1

; а2 =

П1 =

c c

Yiж4 -Y2жi , ni = y2z - 2Yiжi , ; П 2 _ ;

2 r\/ i i i\

с = zж4 - 2ж1 ; z = 2

На рис. 1 - 3 представлены кривые изменения макроскопических упругих постоянных древесно-минерального композита (арболита) соответственно продольного и поперечного модулей Юнга E*33 и E* а также мо-

128

ЛЕСНОЙ ВЕСТНИК 3/2012

ДЕРЕВООБРАБОТКА

О 0,1 0,4 0,6 0

Рис. в

Рис. 7

дуля сдвига ц 12 плоскости, перпендикулярной направлению волокон - от концентрации связующего £ Из этих рисунков видно, что увеличение соотношения k приводит к возрастанию про-и уменьшению

5 Рис. 6

и коэффициентом влажностных деформаций Пс = 0,36.

Соответствующие величины для заполнителя равны

6 1

af = 3-10е

; % = 5,7.

дольного модуля Юнга E 33 поперечных модулей E11 и ц*12.

На рис. 4 - 7 изображены изменения макроскопических коэффициентов температурного и влажностного расширения в продольном и поперечном направлениях по отношению к оси ориентации волокон.

При расчетах принималось, что связующее обладает коэффициентом температурного расширения

61

a = 15-10е

градС

градС

Библиографический список

1. Щербаков, А.С. Арболит / А.С. Щербаков, Л.П. Хо-рошун, B.C. Подчуфаров. - М.: Лесная пром-сть, 1979. - 200 с.

2. Хорошун, Л.П. Вычисление упругих свойств арболита / Л.П. Хорошун, А.С. Щербаков // Научные тр. Московского лесотехнического института.

- 1976. - Вып. 93. - С. 161-166.

3. Щербаков, А.С. Некоторые вопросы теории прочности и деформативности арболита / А.С. Щербаков // Научные труды МЛТИ. - 1960. - Вып. 127.

- С. 5-20.

4. Щербаков, А.С. Вопросы теории упругости и деформативности композиционных материалов / А.С. Щербаков и др. // Научные труды МЛТИ.

- 1969. - Вып. 222. - С. 131-145.

5. Щербаков, А.С. Методы расчета и прочности и деформативности арболита / А.С. Щербаков // Научные труды МЛТИ. - 1969. - Вып. 216. -С. 5-11.

ЛЕСНОЙ ВЕСТНИК 3/2012

129