Сужение «Вилки» Фойгта-Рейсса в теории упругих структурно неоднородных в среднем изотропных композиционных тел без применения вариационных принципов Текст научной статьи по специальности «Физика»

ЭЛЕКТРОННЫЙ НАУЧНЫЙ ЖУРНАЛ «ДРтот. СЕРИЯ: ЕСТЕСТВЕННЫЕ И ТЕХНИЧЕСКИЕ НАУКИ»

№ 3 2014

УДК 539.3

СУЖЕНИЕ «ВИЛКИ» ФОЙГТА-РЕЙССА В ТЕОРИИ УПРУГИХ СТРУКТУРНО НЕОДНОРОДНЫХ В СРЕДНЕМ

ИЗОТРОПНЫХ КОМПОЗИЦИОННЫХ ТЕЛ БЕЗ ПРИМЕНЕНИЯ

ВАРИАЦИОННЫХ ПРИНЦИПОВ

Тарасюк Иван Александрович

студент

Кравчук Александр Степанович

д-р физ.-мат. наук Белорусский государственный университет, Минск (Беларусь)

аи^ог@аргюг1-]оигпа!. ги

Аннотация. В теории композиционных тел сложился определенный стереотип - для сужения «вилки» Фойгта-Рейсса необходимо использовать вариационный принцип и получить более узкую и, соответственно, точную «вилку» Хашина-Штрикмана. Данная последовательность действий в той или иной форме воспроизводится в большинстве научных работ посвященных механике композиционных материалов. При этом следует подчеркнуть, что экспериментального подтверждения верности получаемого сужения («вилки» Хашина-Штрикмана) до настоящего времени выполнено не было. В данной работе предлагается другой теоретический подход к сужению «вилки» модулей упругости композиционного материала без использования вариационных принципов. Он основан на вычислении эффективных матриц жесткости и податливости, определяющих связи усредненных по правилу смеси напряжений и усредненных по правилу смеси деформаций, вычисленных по гипотезам Фойгта и Рейсса в элементарном объеме композиционного материала. На последнем шаге вычисляется математическое ожидание коэффициентов полученных матриц, определяющих суженую «вилку». В отличие от устоявшегося в научной литературе мнения установлено, что только модули упругости, стоящие на диагоналях матриц податливости, полученных по Фойгту и Рейссу, всегда являются верхней и нижней границей реальных модулей упругости композитов. Остальные коэффициенты матриц могут (в смысле верхней и нижней оценки) меняться местами, что зависит от конкретной реализации упругих характеристик компонент композиционного материала.

Ключевые слова: гомпозиционный изотропный материал; гипотеза Фойгта; гипотеза Рейсса; сужение «вилки» модулей упругости композиционного материала.

CALCULATION OF VOIGT-REUSS «RANGE» IN THE THEORY OF ELASTICITY OF STRUCTURAL HETEROGENEOUS, AVERAGE ISOTROPIC, COMPOSITE BODIES WITHOUT APPLICATION OF VARIATIONAL PRINCIPLES

Tarasyuk Ivan Alexandrovich

student

Kravchuk Alexander Stepanovich

doctor of physical and mathematical sciences Belarus State University, Minsk (Belarus)

Abstract. In the theory of composite bodies formed stereotypical - it is necessary to use a variational principle to narrow the «range» of Voigt-Reuss and obtain, respectively, a more precise «range» of Hashin-Shtrikman. This sequence of operations appears in one form or another in most scientific papers dealing with the mechanics of composite materials. It should be emphasized that the experimental confirmation of validity of resulting narrowing («range» of Hashin-Shtrikman) hasn’t been met so far. In this paper we propose a different approach to narrowing «range» of the elastic modules of the composite material without the use of variational principles. It is based on the calculation of effective flexibility and stiffness matrices defining the relations between averaged by mixture rule stress and averaged mixture rule deformations. It is calculated by Voigt and Reuss hypotheses in a volume element of the composite material. The last step is to calculate the expectation coefficients of obtained matrices which are determined restriction «range». In contrast to the well-established in the scientific literature opinion it is found that only the elastic modules which stand on the diagonal flexibility matrices obtained by Voigt and Reuss can be upper and lower limits of the real elastic moduli of composites. The remaining coefficients of the matrices may change places (in the sense of the upper and lower bound), that depends on the particular implementation of the elastic characteristics of component of the composite material.

Key words: composite isotropic material; Voigt hypothesis; Reuss hypothesis; narrowing «range» of the elastic moduli of the composite material.

Введение. В последние годы существенное развитие получили исследования физико-механических свойств материалов, температурных полей и напряженно-деформированного состояния в элементах конструкций с учетом их структуры. Результаты этих исследований находят практическое применение при расчете на прочность элементов конструкций из неоднородных материалов [1].

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

Постановка задачи. Для решения задачи определения эффективных модулей рассматривается элемент композиционного материала, на границе которого задаются воздействия, имитирующие воздействия, возникающие в испытательных машинах при проведении серии опытов (чистое растяжение, кручение, всестороннее сжатие), для определения полного набора модулей однородного изотропного материала [2]. Макроточкой называется элементарный макрообъем, размеры которого значительно превосходят характерные размеры неоднородностей, однако существенно меньше размеров тела [2].

Таким образом, для решения задачи определения эффективных деформационных характеристик элемента композиционного материала в качестве макроточки примем прямоугольный параллелепипед объемом АххАухА . Размеры параллелепипеда значительно превосходят

характерные размеры неоднородностей, но пренебрежимо малы в сравнении с характерными размерами твердого тела.

Будем считать, что параллелепипед состоит из п изотропных компонент, и предполагается, что значения концентраций ук компонент

композиционного материала известны для всего материала, и они же являются объемными долями компонент для любой макроточки твердого тела. Рассмотрим напряженно-деформированное состояние отдельно взятого элементарного параллелепипеда с координатами (х, у, z), имитирующего напряженно-деформированное состояние твердого тела из композиционного материала в данной точке.

Применение гипотезы Фойгта для изотропного в среднем композиционного материала [3]. Используем гипотезу Фойгта о том, что при простейшем нагружении в макроточке тела из композиционного материала имеет место однородная деформация, т.е. для любой компоненты параллелепипеда верны соотношения: ех1 = ех, еу4 = еу, = е2,

£ху,1 =£ху, £уг,1 =€уг, £2х,1 =£2х . Запишем закон Гука для / -ой компоненты элементарного параллелепипеда:

Для нахождения эффективных характеристик по Фойгту определим

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

(1)

(2)

(3)

(4)

°уі > ,і > аху,і > &У2,і > &2Х,і (1)-(4):

п / у\

(&х)р = £/1 ■ & =(,ЛЛ V, ^/.А 1 £х (1-М р )+<П р (£у +£ )),

I =1

(1+И рX1- 2А р)

Ур (1 + А> рI1- 2<А р)

&у)р 2 'Сту,' |1,/,Л »1_У,Л \чр1-\-/р'

I =1

(еу (1 -А> р )+А р (£х +£ )),

п / У\

{&)р = 2/,■& = (1+МЪр-2(у)) £ (1-М р )+И р (£х + £у)),

&х^1р (Р)р £ху, (^)р {р)р '£yz , (Szx)р (р)р 'ezx,

(5)

(6)

(7)

(8)

где усредненные значения (Е)р и (г)р в (5)-(8) определяются выраже-

ниями:

Ер =

п п \ { п

ЕГ, ■ Е у , • Е у

1+ А, ^ 1 - 2у,

V I=1

I=1

V I=1

, • У,

(1 + М)(1- 2п)

(9)

И р =

п

2

V I=1

(1+ И )(1 - 2а)

п

I

V I=1

,• Е

(1+а )(1 - 2а)

(10)

Е

р 1+ А

р

р

(11)

Применение гипотезы Рейсса для изотропного в среднем композиционного материала [4]. Используем гипотезу Рейсса о том, что при простейшем нагружении в композите имеет место однородное напряженное состояние, т.е. для любой компоненты параллелепипеда

верны соотношения: Ох^ =Ъх , Ъу I- =Ъу , &zI = &2 , &ху,! = & , &yz,/ = & yz ,

°2х1 =°2х. Запишем закон Гука для I -ой компоненты элементарного параллелепипеда:

Єх,і = Е- (<7х - У і {(Ту + (У2 )), (1 2)

Еі

Єуі = Е {ту - уі {тх +Тz )) - (13)

Еі

Є2,і = Е {Т -Уі {тх +Ту ))- (14)

Еі

1 + Уі 1 + Уі 1+ Уі /и с\

Єхуі = —^~ Тху - Єуя,і = —^~ Туя - Є2хі = Т2х ■ (1 5)

Еі Еі Еі

Для нахождения эффективных характеристик по Рейссу определим

средне значения деформаций {е^я - (є^к - (Є^я - (є^)е - (еу^я - {е)я

по всем компонентам композиционного материала в элементарном параллелепипеде с помощью статистического усреднения дискретных случайных величин е^- - Єу ,і - е2 - Єхузі - Єу2- е^ (12)-(15):

п

(£х)я = ' Є =Щ~ Т -У Я '{ТУ + Т ))- (16)

і=1

Я

п

£у)к=Т/< '£уі =Е; Ту-уяТ+т )) - (17)

і=1 ' /я

П

е*)я = У/' Є =Е~ {Т -УУ я ■ {тх +ту ))- (18)

м ^>я

еу), =Т- - (е)„ = Ту- - (Єх),, = Т- - (19)

'х"я Ху"я (С)я’ "х'я &я

где усредненные значения (Е)я и У)я в (16)-(19) определяются выражениями:

г

• , Е ,

V 1=1

А

я

п

V I У У

V г=1 У

г

V I V У

V 1=1 У

я

1+ И

(20)

(21)

(22)

я

Основное неравенство, связывающее усредненные модули упругости по Фойгту и по Рейссу. Покажем, что для эффективных модулей Юнга (Е)р и (е)к выполняется неравенство:

Ер > Ея. (23)

Отметим, что выполняется следующее соотношение:

п п

гЕ

» Ег+АаЕг-аГ2>1(1+А )(1 - И)

(Е)р -]гЕ, =-*=---------------------^-;-—-^->0. (24)

/'=1

Е

=1

,Е I

(1+ А )(1 - 2а)

Ввиду того, что знаменатель в выражении (24) строго положительный, достаточно показать, что числитель является неотрицательным. Выполним следующие преобразования:

п п

п п

пит - Ее е

гЕ

. . 1 '-2а '^^(1+А)(1-2а)

I=1 I=1 I=1 I=1

пп

=ЕЕда

с

2-(п+А)-4Аа/ 2-(а+а/)-А +А) Л

г=1 /=г+1

(1+А)(1-2и)(1+А)(1-2а-) (1+а)(1-2а)(1+а,)(1-И)

пп

=ЕЕггЕЕ

=1 / = +1

<А -А/)2 '

(1+А)(1- 2а)(1+А/ )(1-2а )

Из (25) следует, что (24) выполнено.

Докажем, что выполняется еще одно неравенство:

> 0.

(26)

Ввиду того, что знаменатель в неравенстве (26) строго положительный, достаточно показать, что числитель является неотрицательным. Выполним следующие преобразования:

Из (27) следует, что (26) также выполнено. Используя (24) и (26), можно составить следующую цепь неравенств:

Очевидно, что (28) доказывает неравенство (23). Аналогично можно показать, что (р) не превосходит (р). Покажем, что для эффективных

коэффициентов Пуассона (у)р и (у)к нет связывающего их неравенства.

Для этого приведем примеры, демонстрирующие, что Пр может быть и

больше, и меньше (у)к. Рассмотрим разность:

(27)

п

(28)

7 =1

(29)

Не нарушая общности, будем считать, что у1 >у2 > ... >уп. Тогда для того, чтобы выражение (29) было положительным, достаточно при 7 = 1, п -1 выполнения условия:

Е: >

' \

(1 + у )(1 - 2у)

(1 + у7+, )(1 - 2у+1)

Е,

и, наоборот, для того, чтобы выражение (29) было неположительным,

достаточно при 7 = 1, п -1 выполнения условия:

Е7 <

' 1

(1 + у )(1 - 2у)

(1 + у+, )(1 - 2у+1)

Можно показать, что для отношений , и ~т^ также нет связы-

Ер

вающего их неравенства.

Таким образом, установлено, что только модули упругости (9), (20) и (11), (22), стоящие на диагоналях матриц податливости, полученных по Фойгту и Рейссу, могут являться верхней и нижней границей реальных модулей упругости композитов. Остальные коэффициенты матриц могут (в смысле верхней и нижней оценки) меняться местами, что зависит от конкретной реализации упругих характеристик компонент композиционного материала.

Вычисление «вилки» эффективных характеристик. Перейдем к определению «вилки» эффективных характеристик. Для этого запишем элементы матрицы жесткости J [а], используя усреднение с неопределенным коэффициентом а, где 0 <а< 1:

(ех )а = а(е)р + (1 -аХ ех)я = а

■(

»* 1Еа

СХ -у Я \?у +Ъ)1 =

(30)

(Е(4

(ап 1ау + а2

еУ} а=а1\£>)р + (1 ~а^ £")- =

(1 -а)( Е\ +а Е)

у/л г

я

(ЕАЕя

ау --

(1 -аХ Щр у я +а Щя у 1

1

{Е(а))с

------------------1яУ/р (а + а )

(1 -а)< Е)р +а Е)к 1 х г)

= (31)

а у-у(а))а(ах + а

а \ у

е )а = ае)р + (1"а)( е)я =

(1_аХЩР+аЩя („ (1_аХщгуя+ащяу 1

ШЕя

а-----------------------------------------' ■1’У',я'~', 1я'1р (а,+а у

V

(1-аХ Щр +а Щя

1— к -уа(ах+ау

£„,) = а(е,)_ + (1 -а)^ е

ХУ I а \хУ/Р

а

ХУ/Я

<Е(а)}

а аХУ

^а=ае^+(1 -ахеуу!я ща

а

а ау?

ех ) а = ае*х)Р + (1 - а)( е«)я = 1 ]У(Х}а ах

\Е (а]>

= (32)

(33)

(34)

(35)

где усредненные значения (Е(а))а и (у[а)а в (30)-(35) определяются выражениями:

(Е (а))а

ЩрЕя

(1 -а)( щр + а( Е)

(36)

'Я

а

(1 -а)( Щр + а( Щя 10

(37)

Следует отметить, что для композиционного материала, изотропного в среднем, с эффективными характеристиками (36) и (37) выполняется соотношение, связывающее модуль сдвига с модулем Юнга и коэффициентом Пуассона:

_ (£(«)Г

(О (а))*

*

(38)

Запишем элементы матрицы податливости Е (а), используя усреднение с неопределенным коэффициентом а, где 0 < а < 1:

<*) а_а*Х)Р + _

йЕР

ЮЖ)єип)+МДе+є ))+го-ш

к

х

ЛеМПяПШе+є )1_

е+має\Єу+є2

(39)

*)а^Р +М_

аЕ

(1-й) Е

^(^Мр)(і-2Пр) Ы^к!1-

) Ек ,(

Щ()Г +

є-єу +єУ х

а

Г аЕр . о-йЕк Т 1Г йЕРПР , (1-ойЕкПк )

ЛИПЛИПр) МПДМУ — .(^1 МаИ + П НПЛИПр) МПИП)] Шє +єІ і (1-а)Ек V ,(„ + є є)х

(цпагИпагг({)' 1 (1 ГІЇ*\ _1 йЕР

а .Ъ ( а Л [(1+Мр)1-2Пр ..Г йЕр , а-йЕк Ї МПк)М>ЦГ є є є)” 1Г аЕРПР , (1-аЕІяПя )

"1(1+(Пр)1-2Пр) МПДМУ _і+паєНпає)є 1 па) )+П 11 ІМПЛМПр) МПМП)] Жк+єуІ

) _ (40)

є

Ю а=а%)Р+(*-°а а»)к(1+мк)+(1-аХ (1+Мр ))х

где усредненные значения {Б{а)£ и П(а)е в (39)-(44) определяются выражениями:

т£=

= ШР (1-2М«)+(1-аХБк (1-П ЫБт (1+Мк)+(1-аХБк МПт)) (45) <Бр МП к)(!-2П к )+(1-аХ Бк МП т )(1-2П т) ’

(у(а))‘ =

а Бт V Р (1 + V К)(! - 2 (п К)+(1 - аХ Б к V к (1 + V , )(1 - 2 {V Р) (46)

■

Следует отметить, что для композиционного материала, изотропного в среднем, с эффективными коэффициентами (45) и (46) также выполняется соотношение, связывающее модуль сдвига с модулем Юнга и коэффициентом Пуассона:

О (а)У = ^Е . (47)

1 + <г(а))£

Можно показать, что (Е(а))є > (Е(а))*, (О(а))£ > (О(а))*, а также что

/ ( ^ і ( (у(а)У (у(а))а

для (у(а)) , у(а)) и -------—, ------— нет связывающих их неравенств

(Е(а)Г (Е(а¥

для любого коэффициента а, где 0<а< 1, используя рассуждения, описанные в предыдущем пункте.

Определим уменьшенные (по сравнению с «вилкой» Фойгта-Рейсса) «вилки» модулей упругости К ЕЕ е1 коэффициентов Пуассона

[(у)*;(у)£] и модулей сдвига _(О)*;(О)£

в среднем изотропного компози-

ционного материала как среднее значение соответствующих выражений (36)-(38) и (45)-(47):

Е * = |( Е (а))*Са, Е£ = |( Е (а))£ сіа, (48)

і

У = ^(у(а))*Са, У)£ = |(у(а))£Са, (49)

і і

О * = |( О (а))*Са, (О)£ = |( О(а))£Са. (50)

о о

После преобразования выражений (48)-(50), представим эффективные характеристики в виде:

і

і

о

о

о

о

Е * = (ЕЛЕ*1п К ЕЛЕ*)

Е ' ЁТ-Ё»

= п*-Е*п . ЕЛЕ,п-п>*Мх-,^-,„ (52)

п ' (Е„-Е* + (Е, -(Е),)2 ' (52)

Е. ЕрЕ*1пПЕЛЕ*) (53)

Е '------ОНО;----------------• 1531

Е"' 2 « +Е*)+

нр н* (п, -П* )2 (н р -н) - 2НгН* 1п(Н ,/И*))

(Н, - И * )3 :

н г п Р -н * п * И, н * (п,-п * )1п (И г/ И *)

(54)

И, -И*

(И, -И*)

(55)

2

Е"' 2«о)¥ + Е«). (56)

Е Е

где нр '(1+п Д1 - 2 п,)и н * '(ГЩ^Ю ■

На рисунках 1-4 представлены графики «вилки» Фойгта-Рейсса для модулей Юнга (Е^. (е)* (9). (20) (Рис. 1). коэффициентов Пуассона

пР. п* (10). (21) (Рис. 2) и модулей сдвига (о)р. (о)* (11). (22) (Рис. 3.

4). а также сужений этой «вилки». полученных по предлагаемой методике

усреднением эффективных модулей Юнга (Е)а. (Е)е (51). (54) (Рис. 1).

коэффициентов Пуассона (у)а. (у)£ (52). (55) (Рис. 2). и модулей сдвига

Оа. Е)£ (53). (56) (Рис. 3. 4).

Еь ГПа

Рис. 1. Зависимость усредненных модулей Юнга от состава смеси

(я = 0.4, у2 = 0.6, Е2 = 210 ГПа, п1 = 0.27, п2 = 0.3): непрерывная линия - (Е)р, штрихпунктирная линия - (е)к ,

пунктирная линия - (е)точечная линия - (е)а

ш.

Рис. 2. Зависимость усредненных коэффициентов Пуассона от состава смеси = 0.4, у2 = 0.6, Е = 190 ГПа, Е2 = 210 ГПа, У2 = 0.3 ): непрерывная линия - (у) ^, штрихпунктирная линия - (у) к,

пунктирная линия - (у)£, точечная линия - (у)а

Р У ГПа 160

100

120

140

160

180

200

Еь ГПа

Рис. 3. Зависимость усредненных модулей сдвига от модулей Юнга состава смеси (у1 = 0.4, ^2 = 0.6, Е2 = 210ГПа, п1 = 0.27 , п2 = 0.3): непрерывная линия - (р)р, штрихпунктирная линия - (о)к,

пунктирная линия - (р)е, точечная линия - (р)а

Р У ГПа

Рис. 4. Зависимость усредненных модулей сдвига от коэффициентов Пуассона состава смеси = 0.4, у2 = 0.6,

Е1 = 190 ГПа, Е2 = 210 ГПа, У2 = 0.3 ): непрерывная линия - (о^,

штрихпунктирная линия - (о)к, пунктирная линия - (о)е,

точечная линия - (о) °

1 1

£I(а)а и податливости Е = ^Е(а)<Ла, т.е. оценим близость к нулю

О

где I - единичная матрица, ||л|| = тах^Ц-,] - норма матрицы

Произведем оценку взаимообратности матриц жесткости

11

I =

0 0 выражения:

I • Е - Ц, (57)

О

, - норма матрицы л.

> > ]=1

Отметим, получить аналитическую оценку близости к нулю (57) невозможно, однако при подстановке конкретных значений модулей упругости и коэффициентов Пуассона компонент композиционного материала погрешность взаимообратности матриц жесткости I и податливости Е составляет не более 10%.

Выводы. Предложен метод сужения «вилки» Фойгта-Рейсса, определяющей разброс модуля упругости, коэффициента Пуассона и модуля сдвига композиционного материала без использования вариационных принципов.

Установлено, что только модули упругости (9), (20) и (11), (22), стоящие на диагоналях матриц податливости, полученных по Фойгту и Рейссу, могут являться верхней и нижней границей реальных модулей упругости композитов. Остальные коэффициенты матриц могут (в смысле верхней и нижней оценки) меняться местами, что зависит от конкретной реализации упругих характеристик компонент композиционного материала.

1. Победря Б.Е. Механика композиционных материалов. М.: Изд-во МГУ, 1984. 336 с.

2. Кравчук А.С., Майборода В.П., Уржумцев Ю.С. Механика полимерных композиционных материалов. Экспериментальные и численные методы. М.: Наука, 1985. 304 с.

3. Voigt W. Lehrbuch der Kristallphysik. Berlin: Teubner. 1928. 962 s.

4. Reuss A. Berechung der Fliessgrenze von Mischkristallen auf Grund der Plastizitatsbedingung // Z. Angew. Math. Und Mech. 1929. V 9. № 1. S. 49-58.