CC BY
60
Ташкинов М.А. Моделирование упругого поведения многокомпонентных композиционных материалов с использованием приближенных решений стохастических краевых задач // Вестник ПНИПУ. Механика. - 2015. -№ 3. - С. 165-181. DOI: 10.15593/perm.mech/2015.3.12
Tashkinov M.A. Modeling of elastic behavior of multicomponent composite materials basing on the approximate solution of stochastic boundary value problems. PNRPU Mechanics Bulletin. 2015. No. 3. Рр. 165-181. DOI: 10.15593/perm.mech/2015.3.12
ВЕСТНИК ПНИПУ. МЕХАНИКА № 3, 2015
PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/inf/
DOI: 10.15593/perm.mech/2015.3.12 УДК 539.3
МОДЕЛИРОВАНИЕ УПРУГОГО ПОВЕДЕНИЯ МНОГОКОМПОНЕНТНЫХ КОМПОЗИЦИОННЫХ МАТЕРИАЛОВ С ИСПОЛЬЗОВАНИЕМ ПРИБЛИЖЕННЫХ РЕШЕНИЙ СТОХАСТИЧЕСКИХ КРАЕВЫХ ЗАДАЧ
М.А. Ташкинов
Пермский национальный исследовательский политехнический университет, Пермь, Россия
О СТАТЬЕ
АННОТАЦИЯ
Получена: 1 августа 2015 г. Принята: 5 сентября 2015 г. Опубликована: 30 сентября 2015 г.
Ключевые слова:
моментные функции, случайная микроструктура, статистические характеристики, композиты с металлической матрицей, многокомпонентные представительные объемы, геометрическая модель, локальные поля напряжений и деформаций, стохастическая краевая задача, метод последовательных приближений, теория упругости
Статья посвящена методам статистической механики, их разработке и применению для изучения микроструктурного поведения многокомпонентных композитов. Объектом исследования являются структурно-неоднородные материалы, состоящие их более чем двух компонентов.
Целью работы является разработка аналитического инструментария для анализа полей микроструктурных напряжений и деформаций в многокомпонентных средах, позволяющего учитывать геометрические и физико-механические свойства компонентов, на основе вычисления статистических характеристик локальных полей напряжений и деформаций в компонентах.
Исследование поведения компонентов микроструктуры композитов основано на концепции представительного объема материалов. Предполагается, что компоненты являются однородными и изотропными. Информация о внутренней морфологии представительных объемов формализуется с помощью моментных функций различных порядков. В качестве характеристик процессов деформирования выступают статистические моменты (статистики) полей напряжений и деформаций в компонентах материала. Аналитические выражения для статистических характеристик локальных полей напряжений и деформаций получены с использованием решения краевой задачи теории упругости в стохастической постановке. Краевая задача решена с помощью метода функций Грина для упругой среды.
Разработанная аналитическая модель позволяет учитывать как геометрические параметры микроструктуры, так и физико-механические свойства компонентов. Впервые получены выражения для моментов первого и второго порядка локальных полей напряжений для многокомпонентных материалов. Исследованы частные случаи композитов с титановой (Ti) матрицей, армированных случайно расположенными частицами карбида кремния (SiC). Выполнен анализ влияния микроструктурных параметров на поведение каждой из фаз в отдельности. Представлены численные результаты, полученные для статистик полей напряжений и деформаций.
© ПНИПУ
© Ташкинов Михаил Анатольевич - кандидат физико-математических наук, старший научный сотрудник, e-mail: m.tashkinov@pstu.ru
Mikhail A. Tashkinov - Ph.D. in Physical and Mathematical Sciences, Senior Researcher, e-mail: m.tashkinov@pstu.ru
165
Tashkinov M.A. /PNRPUMechanics Bulletin 3 (2015) 165-181
MODELING OF ELASTIC BEHAVIOR OF MULTICOMPONENT COMPOSITE MATERIALS BASED ON APPROXIMATE SOLUTION OF STOCHASTIC BOUNDARY VALUE PROBLEMS
M.A. Tashkinov
Perm National Research Polytechnic University, Perm, Russian Federation
ARTICLE INFO
ABSTRACT
Received: 1 August 2015 Accepted: 5 September 2015 Published: 30 September 2015
Keywords: correlation functions, random microstructure, statistical characteristics, metal matrix composites, multiphase representative volume elements, geometry model, local stress and strain fields, stochastic boundary value problem, successive approximations, theory of elasticity
This work is devoted to development and application of the statistical mechanics approaches to studying microstructural behavior of multi-component composites. The object of the study is heterogeneous materials consisting of more than two components.
The aim is to develop analytical tools for the analysis of the microstructural stress and strain fields in multicomponent media taking into account geometrical and mechanical properties of components and basied on calculation of the statistical characteristics of the local fields of stress and strain in components.
Research of the microstructure behavior of components of composites is based on the concept of representative volume elements. It is assumed that the components are homogeneous and isotropic. Information about the internal morphology of the representative volume is formalized by means of correlation functions of various orders. Characteristics of the deformation processes are the statistical moments (statistics) of stress and strain fields in the components of the material. Analytical expressions for the statistical characteristics of the local stress and strain fields are obtained using the solution of the boundary problem of elasticity theory in a stochastic formulation. The boundary value problem is solved using the Green's functions method for the elastic medium.
The developed analytical model takes into account both the geometrical parameters of the microstructure and mechanical properties of components. For the first time expressions for the moments of the first and second order of local stress fields for the multicomponent materials were derived. Case studies of composites with titanium (Ti) matrix reinforced with randomly distributed particles of silicon carbide (SiC) were investigated. Analysis of the influence of microstructural parameters on the behavior of each phase separately was performed. Numerical results for the statistics of stress and strain fields are presented.
© PNRPU
Введение
Появление и широкое распространение новых композиционных материалов связано с необходимостью решения сложных прикладных задач, продиктованных потребностями высокотехнологичных отраслей промышленности, таких как аэрокосмическая, автомобильная, электронная и др. Создание аналитических инструментов, позволяющих предоставить рекомендации для оптимального проектирования композитов, наиболее соответствующих условиям эксплуатации конкретных приложений, является актуальной проблемой. В отдельный класс можно выделить многокомпонентные композиты, состоящие их двух и более фаз, отличающихся по своим исходным свойствам. Развитие подходов к многомасштабному моделированию механической микроструктуры таких материалов дает возможность прогнозирования эффективного поведения композиционного материала с учетом микроструктурных особенностей, таких как размер частиц фаз, их объемная доля, ориентация, дисперсия и кластеризация [1-6].
Метод, используемый в данной работе, основан на статистических подходах и теории случайных функций. Его суть заключается в том, что поведение отдельных компонентов микроструктуры при нагружении представительного объема оценивается с помощью статистических характеристик локальных полей напряжений и деформаций, таких как моменты первого и второго порядка [7-9], которые вычисляются с помощью решения стохастических краевых задач.
166
Ташкинов М.А. /Вестник ПНИПУ. Механика 3 (2015) 165-181
Стохастические краевые задачи могут быть сформулированы в реализациях случайной структуры [10] и в моментных функциях [9, 11-13], когда усреднение числа реализаций заменяется усреднением по объему. Решение краевых задач в последнем случае сводится к интегрально-дифференциальному уравнению, содержащему функцию Грина [13]. Уравнение может быть решено различными способами, каждый из которых основан на определенном упрощающем предположении о статистических свойствах материальных характеристик полей композита, что позволяет упростить бесконечную цепочку стохастических уравнений, образованных путем последовательного статистического усреднения определяющих соотношений для компонентов в неоднородных средах [14]. Начало изучению деформирования композитов со случайными структурами в рамках статистического подхода положено фундаментальными работами И.М. Лифшица и Л.Н. Розенцвейга, в которых приведена постановка стохастической краевой задачи теории упругости для поликристаллических сред [15]. Впоследствии этот подход получил развитие в работах В.В. Болотина [15], В.А. Ломакина [12], Т.Д. Шермергора [13], С.Д. Волкова и соавторов [11], Г. А. Ванина [10], а также М. Берана [14], Ю.В. Соколкина, А. А. Ташкинова [9, 17], А.А. Панькова [18], Л.П. Хорошуна [19], Л.А. Сараева [20] и др.
Несмотря на большое количество вариаций статистических подходов, существуют различия в получаемых различными методами результатах для эффективных свойств композитов и в большей мере для статистических характеристик полей напряжений и деформаций в компонентах композитов. Это обусловлено статистической нелинейностью исследуемой задачи и построением приближенных решений, которые в различной степени соответствуют физической модели композита, в частности, его структуре. Отсюда необходимость в получении более точного приближенного решения задачи в рамках существующего подхода, а также обобщения подхода для решения более широкого круга задач.
Данная работа посвящена методологии получения статистических характеристик для многокомпонентных матричных композитов со случайно распределенными включениями на основе решения стохастических краевых задач теории упругости для неоднородных сред.
1. Методология описания характеристик случайной микроструктуры
Внутренняя микроструктура композитов зачастую представляет собой случайную гетерогенную многофазную среду. Согласно статистическому подходу для определения морфологии таких материалов на микромасштабном уровне могут быть введены индикаторные функции XC (r) [8, 9, 21, 22]. Значения этих функций зависят от положения радиуса-вектора в представительном объеме. В частности, XC (г) = 1, если радиус-вектор находится в компоненте C, и XC (г) = 0 в других случаях. Тогда тензор модулей упругости для представительного объема определяется как сумма произведений тензора модулей упругости Cjki для каждой компоненты и соответствующих индикаторных функций X(z) (r). В результате он принимает форму кусочно-постоянной функции координат:
C,« {?) = SX(z,(r=)C®. (1)
z=1
167
Tashkinov M.A. /PNRPUMechanics Bulletin 3 (2015) 165-181
Гипотеза эргодичности случайных функций Хс (r) позволяет рассматривать усреднение по ансамблю реализаций представительных объемов равным усреднению по объему [11, 14]. Другими словами, (Х^z)(r)^ = (Х^zдля каждого r и г1. Принимая это во
внимание, тензор Cijkl (r) и индикаторные функции Хс (r) для представительных объемов можно разложить на два слагаемых - среднее значение и флуктуацию:
Cmnkl (r) = {с„ш (r)) + Cm,nkl (r), (2)
M?) = <Xc(r)) + XC(f). (3)
где флуктуации определяют степень, с которой значение функции отличается от усредненного. Другими словами, они показывают, как случайная структура будет отличаться от периодической структуры с теми же самыми геометрическими параметрами [24]. Усреднение уравнения (1) дает постоянный изотропный тензор, который зависит от объемных долей компонент:
('Cjj(r}) = £PMj), (4)
z=1
где P(z) =(X(z)(r- объемная доля фазы z в представительном объеме.
Кроме констант, количественно определяющих структурные геометрические основы, инструменты статистического подхода предлагают механизмы оценки пространственного взаимодействия между микроструктурными компонентами [8]. Так, одним из наиболее широко используемых механизмов являются моментные функции, которые характеризуют связь микрочастиц внутри представительного объема [9, 13, 22]. Таким образом, мо-ментная функция n-го порядка может быть определена как среднее скалярного произведения n случайных величин в различных положениях радиусов-векторов r,r1,...,rn.
Выражение для моментной функции n-го порядка флуктуаций случайной индикаторной функции в компоненте C можно записать в общем виде:
(Х'с (Й)ХС (?2)...XC (rn )) = {(XC (й) - Pc )(XC (r2) - Pc ) ... (XC (rn ) - Pc )), (5)
где pC - объемная доля.
Моментные функции чувствительны к таким параметрам, как распределение, ориентация и форма микроструктурных составляющих. Теоретически микроструктурная морфология композитов может быть однозначно определена бесконечным числом момент-ных функций. Эти функции используются в качестве геометрических характеристик в задачах получения статистик высоких порядков для микроскопических полей напряжений и деформаций [8, 9, 13, 16, 22, 24, 27, 28].
Методология моментных функций требует некоторых начальных экспериментальных данных относительно внутренней структуры представительного объема композита. Обычно для получения микроструктурных характеристик используются экспериментальные методы визуализации [23], такие как микрокомпьютерная томография. Другой способ - создание моделей с помощью некоторых заранее определенных геометрических параметров.
168
Ташкинов М.А. /Вестник ПНИПУ. Механика 3 (2015) 165-181
Для того чтобы получить значения моментных функций микроструктуры (5), необходимо знать поля значений индикаторных функций (4). Они могут быть вычислены путем выполнения алгоритма, который разбивает представительный объем сеткой с фиксированным шагом и проверяет наличие фаз в каждом ее узле. Точность функций при этом зависит от шага сетки.
В терминах статистического подхода локальные поля, такие как поля перемещений ит, деформаций s^- и напряжений aij- внутри представительного объема со случайной
микроструктурой, как правило, считаются случайными функциями. Они зависят от радиуса-вектора и могут быть представлены аналогично формулам (3) и (4):
(Г) = (um (Г )) + U'm (Г) , (6)
(r ) = (sy (r )) + s'y (r ), (7)
(Г ) = (ау. (r )) + a'j (r ). (8)
Многоточечные статистики, рассчитанные для полей напряжений и деформаций в представительном объеме и в каждой фазе отдельно, могут быть использованы для анализа влияния микроструктурных параметров на поведение гетерогенных материалов. Средние (по представительному объему) значения в формулах (6)-(8) являются константами, поэтому статистики строятся для флуктуаций. Распространенные модели для расчета эффективных характеристик материала, как правило, используют моменты локальных полей первого и второго порядка [7, 8]. Моменты первого порядка представляют собой средние значения, в то время как моменты второго порядка также называют дисперсиями. Моменты полей деформаций могут быть определены с использованием флуктуаций перемещений
и соотношения Коши:
V (г ) = 2 (j (г ) + иj. (г ))
. Флуктуации напряжений
можно пред-
ставить с учетом уравнения состояния а- (r ) = Ci]kl (r )skl (r ) в следующем виде:
a'j (r ) = a'j (r ) ~{aij) = C'ykl (r ) (s kl (r Й ~(C'ykl (f)s'kl (r ^ + (Cjkl (r ^ skl (r ). (9)
Формулы для моментов первого и второго порядка деформаций и напряжений в компонентах были впервые получены С. Волковым и выражаются через смешанные моменты, объединяющие флуктуации индикаторных функций и локальные поля деформаций [11, 25]:
— =+TjC(?$ ^C ^ К ^)}, (10)
j =Ы +^1^ (гк(г)), I11»
R(f Kp(r)C = R(r)s'af (r)) + -(s,,)c{sap )C +
(12)
+ ^ ^ ((C (( ((p (( + e,j{ГС (r Kp (Г ^ + eap(X’C (г )sj (г ^),
169
Tashkinov M.A. /PNRPUMechanics Bulletin 3 (2015) 165-181
fa](A)<P(A))c = fa'] (F)<P(A)) + faj){°*v)-{°«)c Kp )C +
(13)
( (r)a,y (F)a«P (Ч ( + faj) fa (r Kp (r )) + {°av) fac (ЧК; (A)) ) •
Численный расчет статистик локальных полей связан с нахождением флуктуаций перемещений, которые могут быть найдены из решения краевой задачи в стохастической постановке [9, 25].
2. Решение стохастической краевой задачи в упругой постановке
Стохастическая краевая задача содержит следующие уравнения:
Gy, j(r) = 0, (14)
ч(Ч)=12 (ui, y«)+uji«)), (15)
Vj (4 ) = cyki (4 )eki(r X (16)
ui(4)l?еГг = j, (17)
где (14) - уравнение равновесия; j обозначает частную производную
д/
'дх
; (15) - соот-
ношение Коши; (16) - уравнение состояния. Адгезия между компонентами считается идеальной.
Для получения решения для флуктуации перемещений в виде рекуррентного интегрально-дифференциального уравнения был использован метод функций Грина [13, 25]:
д<(yJ(r ) = с dGim (г, t\)
дх.
[■
дх,
Cmnkl (Г1 )ekl + Cmnkl (A )
дик11} (A) дх.
’1n
dVx,
(18)
где x - порядок аппроксимации; Cmnkl (r) = i^(z)(r)С^Пи • Последовательная процедура,
z=1
когда каждое предыдущее приближение подставляется в правую часть, позволяет получить более точные решения [25, 26].
Когда размеры представительного объема много меньше характерных размеров твердого тела, а также в случае изотропного тензора (ci]mn (r), в качестве функции Грина
Gim (r, A1) может быть использован тензор Кельвина-Сомильяны [11, 25]:
Gmk (r, a ) = Ат-^Ч + B-(Ги Pm )(r Pk)
'-p\
-p|3
(19)
где А и B - константы, зависящие от свойств материала микроструктурных фаз. С использованием первого приближения уравнения (18)
ди'(1)(г)
дх,
= e
kl
i CAkl j Gim , (A, A )Ч)(А)] dV1
z=1
(20)
170
Ташкинов М.А. /Вестник ПНИПУ. Механика 3 (2015) 165-181
флуктуации деформаций из уравнения (7), которые используются в формулах для моментов, выражаются с помощью уравнения (20):
f .. \
C1 f ( (( ) + Gm>, (f,f)){V{!,
v-
i(f) = -2(<J(f)+uj,(f)) =2 ekS
Решение (20) содержит флуктуации индикаторной функций под интегралом. Формулы для моментов (-0)—(-3) также включают эти флуктуации. Таким образом, в аналитических выражениях выделяются в явном вид структурные моментные функции (5). Например:
(к;y)(fК(f))“> = 2еи CG'k f (G„,,J (f, f) + Gi„,, (f, f)))dV- +
ек$(г, f-)
\ V
+ S cmi,к (f,г-)+g
z=l, ИФ y V
jm.
G r )№) JV !(r,r- >)——)X--------------LdV-
u-hn
(22)
где K(2)(f,f) = {K’(y —(f )K(y)(fi)j - моментная функция второго порядка.
Можно показать, что второе слагаемое в выражении (22), содержащее произведение флуктуаций индикаторных функций различных компонент ^К(y)(f)K(z)(f)^, обращается
в ноль. Так, по определению, флуктуация индикаторной функции в компоненте C выражается как
KC (f) = (Kc (f)) -Kc (f),
(23)
где (Xc (f)} является константой и отражает объемную долю компоненты в представительном объеме. Тогда осредненное произведение флуктуаций индикаторных функций для двух разных компонент вырождается в константу:
(К(x)(f )К(y) (f ^ = ((P(x) - К(x) (f )) (P(y) - K(y)(f ))) = Const, (24)
что обращает производную в подынтегральном выражении в уравнении (22) в ноль. Следовательно,
(К(y)(f )si (f ^ = 2 eklC(myil f {G,m,( ((, f ) + Gjm,, (f, f ))^^V-2 V CX-n
(25)
Аналогично
(к y )(f )sj(f )s«p(f ^=-4 ekie^AcmniicSAx
fr/ w -.d3 K(3,(f, f-, f2)
xf f (G,m.i((f ) + Gjm_i(f,f ))((((2 ) + Gp,.„ (f, f2)) ~^dVd-,
(26)
V- V,
л
где KiC)(f, f-) = (kC(f )KC (/D) - моментная функция второго порядка; (f, f-, f,) =
= {K (f )KC ((Dk;, ft)) - моментная функция третьего порядка.
Ш
Tashkinov M.A. /PNRPUMechanics Bulletin 3 (2015) 165-181
С учетом равенства (9) смешанные моменты, входящие в выражения для моментов первого и второго порядка полей напряжений в компонентах и содержащие флуктуации полей напряжений, выражаются через флуктуации полей деформаций:
= e.
(z )
'ijkl
(27)
(К y )(r )oJ(f)) =
- £ p, r) P( z C
z=l, zФ y
- (K( y)(r )К( y)(r )e'a (r ^ Ciy + (K( y)(r )г’к1(r ^ (Crn) ■
Момент (oj (r )o,p (r)j равен дисперсии напряжений в представительном объеме как
однородной среде:
fa (fKP (r))
= ekiei)h
(28)
v d(2)c(z)C(z) - V P P C(y)C(z)
/ j^(z) ^ iikl^raRAh / j P(y)P(zy^ijkl^apфh
4z=1 z=1, z Ф y
+ea(C,w,) £ Cj,( К z )(r )eih (f)\ + eth(Cm) £ C<«„( X(, )(f)s’t,1 (r))-
z=1 z =1
+ (C,„) (CaW,}{ (r )eih (r)) - e„ £ CjC ■\h {К’, ,(r)K(:, (r)e;„ (r)) -
z=1
-eth £ cjcк; z,(? )K ..,;? к, (f)}-{cm) £ Cй„( к; ; ,(? %„(f к ;f)) -
z=1 z=1
-(C,^) £ cjzi!{ к; z (f)<l, (n)+£ Clgc^h (К z ,(m; ,(r)e;„ (r)4, (n),
z=1 z=1
где = {КС (r )X'C (r)) = pC (1 - pC) - центральный момент второго порядка (дисперсия)
для индикаторной функции компоненты C ■
В силу свойств индикаторных функций принимать значения, равные либо нулю, либо единице, связь моментов, содержащих произведение индикаторных функций от одного и того же аргумента, выражается следующим образом:
(КС (г )Кс (r )г'ы (г)) = (1 - 2 Pc )(КС (r )г'ы (г)),
(к (r )К'с (r )К'с (r Ki (r)) = (1 - 3DC2) ) (КС (r Ум (r)),
(К'с (r)Кс(rУы(rKh(r)) = (1 - 2Pc ) (Кс(rУм(r)sU(r)) + Dc} {г'ы(rKh(r^,
(К'с(r )К'с(r )К'с(r Ум(r Kh(r ^ = (1 - 3Dc})( Кс(r Ум(r Kh(r ^ + Dc)( г'и(r Kh(r ^,
где DCn =(1 - pC У pC +(-pC У (1 - PC ) - центральный момент n -го порядка для индикаторной функции компоненты C ■
(К( y )(r )o’j (r ^ =
(29)
= e
kl
DQCk - £ P;yPzCk +(Ky,(?)6l(?))l(CjkJ-(1 - 2p)C)>\
z=1, z Ф y
172
Ташкинов М.А. /Вестник ПНИПУ. Механика 3 (2015) 165-181
Для дисперсий напряжений в представительном объеме
f n n Л
К (Г Кц (Г) = v* I - I Py) P( , filled h
у z=1 z=1, z Ф y
+4,(cw) I c^( ц z У к>, (П)+%{f0} ICUV,( Ц z>(f У» (f))+
z=1 z=1
+{С#в)(СаИв)( (f ** (f)) - e„, I СЩС,UV, (1 - 2 P( z) )(Ц z ,(r Ц,, (r)) -
z=1
I esc,uv, (1 - 2 p( .., x^; z ук, (о) - (c„) I c&( ц .■ у к, (f к, (f>)-
z=1 z=1
-{ Qv.) I fzl (Ц z)((f )Ek, (?)} +
z=1
+I cgc,UV, ( - 2 P( z) (z ( (f ))+d<?>( e;, (?&, (f ))).
(30)
Аналогичным образом выводится выражение для (X'c(r )aj(r )a,p(r )) .
Полученные аналитические выражения для статистик деформаций и напряжений зависят от свойств материала (через константы, входящие в функцию Грина) и формализованной информации о микроструктурной геометрии (посредством моментных функций различных порядков). Порядок необходимых моментных функций зависит от точности решения (18). С возрастанием его приближения выражения, содержащие интегралы в уравнениях (22) и (23), становятся более громоздкими и включают моментные функции более высокого порядка.
3. Численная реализация методики на примере композитов с металлической матрицей
В качестве частного случая были рассмотрены композиты с металлической матрицей с матрицей из титана (Ti), армированные включениями из карбида кремния (SiC). По сравнению с традиционными материалами и композитами с полимерной матрицей, используемыми для многих прикладных задач, такие композиты проявляют себя лучшим образом при высоких температурах, большую прочность, повышенную устойчивость к износу и высокоцикловой усталости, имеют более высокое отношение жесткости к весу, значительную электрическую и тепловую проводимость.
Микроструктурные составляющие представительного объема композитов были смоделированы как случайно расположенные, непересекающиеся сферы различных размеров. Были рассмотрены три фазы внутри представительного объема: Ti-матрица, армирующие SiC-частицы и пустоты (поры) в разных пропорциях. Такие модели соответствуют данным изображений, полученных из микроструктурных экспериментов, например в [6] и [29].
Геометрические характеристики моделируемых представительных объемов отображены в табл. 1. Поскольку предложенная методика не зависит от масштабных факторов, все величины, характеризующие размеры, представлены в условных масштабных единицах. В качестве представительного объема был рассмотрен куб со стороной 200.
173
Tashkinov M.A. /PNRPUMechanics Bulletin 3 (2015) 165-181
Таблица 1
Геометрические характеристики многокомпонентных представительных объемов
Предст. объем Объемная доля, p, % Количество включений, N Минимальный радиус Максимальный радиус
Фаза #1 (SiC) Фаза #2 (Поры) Фаза #1 (SiC) Фаза #2 (Поры) Фаза #1 (SiC) Фаза #2 (Поры)
1 10 1 222 8 8 8 12
2 15 1 327 8 8 8 12
3 20 1 445 8 8 8 12
4 25 1 603 8 8 8 12
6 10 4 351 8 8 8 12
7 10 8 511 8 8 8 12
8 10 12 685 8 8 8 12
Первые четыре структуры отличаются объемной долей первой фазы (SiC). Структуры 6-8 имеют одинаковую долю SiC-включений, но различные концентрации пор (от 4 до 16 %). Все величины представлены в условных масштабных единицах. ЭБ-графические изображения этих структур представлены на рис. 1.
а
б
в
д
е
г
с
Л1Г‘
t//C
Рис. 1. Трехкомпонентные представительные объемы материалов с различными объемными долями фаз: а - 10-1 %; б - 15-1 %; в - 20-1 %; г - 25-1 %; д - 10-4 %; е - 10-8 %;
ж - 10-12 %
Заданные механические свойства SiC-включений и Ti-матрицы представлены в табл. 2. Поры считались заполненными вакуумом.
Таблица 2
Механические свойства фаз для случая Ti + SiC-композитов
Включения (SiC) Матрица (Ti)
Упругий модуль EC = 410 ГПа EM = 116 ГПа
Коэффициент Пуассона vC = 0,14 v m = 0,32
174
Ташкинов М.А. /Вестник ПНИПУ. Механика 3 (2015) 165-181
Нагружение представительного объема соответствовало простому сдвигу и было определено путем задания компонент макромасштабного тензора деформаций etj в граничных условиях (17):
' 0 10—6 0 ^
e - 10—6 0 0 . (31)
0 V 0 0 7
В случае когда матрица и включения являются упругими, однородными и изотропными, тензор упругих модулей для каждой компоненты является постоянным и может быть представлен в следующем виде:
Cj - Xw8„8* + ци (8*8, +8„8Д ), (32)
где X(z), ц(z) - постоянные Ламе компоненты z,
X( z)
E(z) v( z) (33)
(1 + v( z)) — 2v( z)),
E(z) (z) - E
(34)
' 2 (1 + v(z)),
где E(z) - упругий модуль; v( z) - коэффициент Пуассона.
Коэффициенты A и B функции Грина (19) в случае упругих компонентов связаны с константами Ламе:
A - (X) + (3ц) в - + M
8тс(ц)(Х + 2ц) 8тс(ц)(Л + 2ц)
(35)
где (X)-ZX(z)p(z)> (^-S^(z)P(z).
z-1 z-1
Методология, описанная в п. 1, позволяет вычислять дискретные значения момент-ных функций, которые не могут быть непосредственно использованы в аналитических выражениях для статистических характеристик. Переход от дискретных значений к непрерывным функциям осуществляется с помощью аппроксимации. Аппроксимация функций может быть проведена путем построения нелинейной аналитической модели со структурой, которая соответствует некоторому заданному выражению для последовательно вычисленных дискретных значений моментных функций. С помощью вычислительных экспериментов ранее было показано [25, 30-32], что степенное выражение с двумя периодическими функциями и тремя независимыми коэффициентами в наибольшей степени подходит для аппроксимации (разница составляет не более 1-1,5 %):
fKn) (rn)- exp(—qAr)
+ c3 sin
(36)
Здесь Ar - расстояние между радиусами-векторами; ci - коэффициенты аппроксимации. Степень a в соотношении (36) равна 2 для функций второго порядка и 1 - для функций
175
Tashkinov M.A. /PNRPUMechanics Bulletin 3 (2015) 165-181
более высоких порядков. Предполагается, что все компоненты являются изотропными, что позволяет принимать все расстояния между радиусами-векторами одинаковыми [24]:
Ат = |/- - Я,| = r - г,| = |r2 - ^1 = - = |rn-i - rn\. (37)
Для каждого частного случая были вычислены моменты первого и второго порядка полей напряжений и деформаций в компонентах. Подробности об используемых при этом численных методиках содержатся в работах [30, 31].
На рис. 2 отображены изменения статистических характеристик как для титановой матрицы, так и для включений в зависимости от объемной доли частиц карбида кремния.
Графики демонстрируют близкую к линейной зависимость для моментов первого порядка - средних значений напряжения и деформации, в то время как дисперсии полей во включениях значительно выше, чем в матрице. С увеличением объемной доли средние деформации как в матрице, так и во включениях непропорционально увеличиваются, что соответствует предсказываемому поведению. Дисперсии напряжений и деформаций в компонентах имеют тенденцию к схождению с увеличением объемной доли армирующих включений.
0,12 0,14 0,16 0,18 0Д0 0,22 0,24
а
0,12 0,14 0,16 0,18 0Д0 0,22 0Д4
в
0,12 0,14 0,16 0,18 0Д0 0,22 0Д4
г
б
Рис. 2. Значения моментов первого (а, в) и второго (б, г) порядков для представительных объемов Ti+SiC-композитов в зависимости от объемной доли карбида кремния: включения (SiC),
‘ матрица (Ti)
176
Ташкинов М.А. /Вестник ПНИПУ. Механика 3 (2015) 165-181
Рис. 3. Моменты первого (а, в) и второго (б, г) порядков для Ti+SiC-композитов в зависимости от объемной доли пор: включения (SiC), - ; - матрица (Ti)
б
а
в
г
Рис. 3 позволяет проследить влияние микроструктурных дефектов, которые были смоделированы, как поры. Как видно, увеличение пористости практически не оказывает влияния на матрицу. В то же время для статистик первого и второго порядка наблюдаются нелинейные зависимости от объемной доли пор, что свидетельствует о том, что не только количество дефектов, но и их распределение и позиционирование в представительном объеме оказывает влияние на напряжения и деформации в компонентах.
Заключение
Данная работа демонстрирует комплексную методику для получения статистических характеристик полей микроструктурных напряжений и деформаций в многокомпонентных материалах с помощью статистического подхода. Описание геометрии случайной микроструктуры проводилось с использованием методологии моментных функций. Аналитические выражения для статистик были получены с помощью стохастического решения краевой задачи в виде перемещений и учитывают внутреннюю информацию о морфологии посредством моментных функций различных порядков. Точность численных вычислений для статистик с использованием представленной модели может быть улучшена путем использования последовательных приближений решений уравнения (18) более высоких порядков.
Одним из возможных применений для самих статистик локальных микроструктурных полей является создание статистических моделей разрушения в компонентах, где они используются для обнаружения начала разрушения, а также для создания на их основе новых критериев прочности, которые могут адекватно отражать процессы деформации и разрушения.
177
Tashkinov M.A. /PNRPUMechanics Bulletin 3 (2015) 165-181
Исследование возможных вариаций механических и геометрических свойств микроструктуры позволяют оценить их влияние на внутренние поля напряжений и деформаций. Эти результаты могут быть использованы при проектировании новых композитов с предопределенными свойствами и оптимизированной микроструктурой, которые могут быть пригодны для широкой области возможных применений.
Работа выполнена при поддержке гранта Президента Российской Федерации для государственной поддержки молодых российских ученых - кандидатов наук (МК-5172.2015.1).
Библиографический список
1. Liu K.C., Ghoshal А. Validity of random microstructures simulation in fiber-reinforced composite materials // Composites. Part B: Engineering. - 2014. - Vol. 57. - Р. 56-70. DOI: 10.1016/j.compositesb.2013.08.006
2. Rasool A., Bohm H.J. Effects of particle shape on the macroscopic and microscopic linear behaviors of particle reinforced composites // International Journal of Engineering Science. - 2012. - Vol. 58. -Р. 21-34. DOI: 10.1016/j.ijengsci.2012.03.022
3. Mishnaevsky Jr.L., Derrien K., Baptiste D. Effect of microstructure of particle reinforced composites on the damage evolution: probabilistic and numerical analysis // Composites Science and Technolog. - 2004. - Vol. 64. - Iss. 12. - Р. 1805-1818. DOI: 10.1016/j.compscitech.2004.01.013
4. Influence of waviness and curliness of fibres on mechanical properties of composites / A.Y. Matveeva, S.V. Pyrlin, M.M.D. Ramos, H.J. Bohm, F.W.J. van Hattum // Computational Materials Science. - 2014. - Vol. 87. - Р. 1-11. DOI: 10.1016/j.commatsci.2014.01.061
5. Yu M., Zhu P., Ma Y. Effects of particle clustering on the tensile properties and failure mechanisms of hollow spheres filled syntactic foams: A numerical investigation by microstructure based modeling // Materials & Design. - 2013. - Vol. 47. - Р. 80-89. DOI: 10.1016/j.matdes.2012.12.004
6. Effects of random particle dispersion and size on the indentation behavior of SiC particle reinforced metal matrix composites / R. Ekici, M. Kemal Apalak, M. Yildirim, F. Nair // Materials & Design. -2010. - Vol. 31. - Iss. 6. - Р. 2818-2833. DOI: 10.1016/j.matdes.2010.01.001.
7. Ghossein E., Levesque M. Homogenization models for predicting local field statistics in ellipsoidal particles reinforced composites: Comparisons and validations // International Journal of Solids and Structures. - 2015. - Vol. 58. - Р. 91-105. DOI: 10.1016/j.ijsolstr.2014.12.021.
8. Buryachenko V. Micromechanics of heterogeneous materials. - New York: Springer, 2007. - 686 p.
9. Вильдеман В.Э., Соколкин Ю.В., Ташкинов А.А. Механика неупругого деформирования и разрушения композиционных материалов. - М.: Наука, 1997. - 288 с.
10. Ванин Г.А. Микромеханика композиционных материалов. - Киев: Наукова думка, 1985. -302 с.
11. Волков С.Д., Ставров В.П. Статистическая механика композитных материалов. - Минск: Изд-во Белорус. гос. ун-та, 1978. - 208 с.
12. Ломакин В.А. Статистические задачи механики твердых деформируемых тел. - М.: Наука, 1970. - 139 с.
13. Шермергор Т.Д. Теория упругости микронеоднородных сред. - М.: Наука, 1976. - 400 с.
14. Beran M. J. Statistical continuum theories. - New-York: Wiley. Intersci. Publ., 1968. - 493 p.
15. Лифшиц И.М., Розенцвейг Л.Н. К теории упругих свойств поликристаллов // Журн. эксперимент. и теорет. физики. - 1946. - № 16. - Вып. 11. - С. 967-980.
16. Болотин В.В., Москаленко В.К К расчету макроскопических постоянных сильно изотропных композиционных материалов // Изв. АН СССР. Механика твердого тела. - 1967. - № 3. -С.106-111.
178
Ташкинов М.А. /Вестник ПНИПУ. Механика 3 (2015) 165-181
17. Соколкин Ю.В., Ташкинов А.А. Механика деформирования и разрушения структурнонеоднородных тел. - М.: Наука, 1984. - 116 с.
18. Паньков А.А. Статистическая механика пьезокомпозитов. - Пермь: Изд-во Перм. гос. техн. ун-та, 2009. - 480 с.
19. Хорошун Л.П. Методы случайных функций в задачах о макроскопических свойствах микронеоднородных сред // Прикл. механика. - 1978. - Т. 14. - Вып. 2. - С. 3-17.
20. Сараев Л.А., Глущенков В.С. Неупругие свойства многокомпонентных композитов со случайной структурой. - Самара: Изд-во Самар. ун-та, 2004. - 163 с.
21. Tashkinov M.A., Vildeman V.E., Mikhailova N.V. Method of successive approximations in a stochastic boundary-value problem in the elasticity theory of structurally heterogeneous media // Composites: Mechanics, Computations, Applications. - 2011. - Vol. 2. - No. 1. - P. 21-37.
22. Torquato S. Random heterogenous materials, microstructure and macroscopic properties. -Springer, 2001. - 701 p.
23. Lee H., Gillman A.S., Matous K. Computing overall elastic constants of polydisperse particulate composites from microtomographic data // Journal of the Mechanics and Physics of Solids. - 2011. -Vol. 59. - Iss. 9. - P. 1838-1857. DOI: 10.1016/j.jmps.2011.05.010
24. Torquato S. Modeling of physical properties of composite materials // International Journal of Solids and Structures. - 2000. - Vol. 37. - P. 411-422.
25. Tashkinov M.A. Methods of Stochastic Mechanics for Characterization of Deformation in Randomly Reinforced Composite Materials // Mechanics of Advanced Materials / Eds. V.V. Silberschmidt, V.P. Matveenko. - Springer, 2015. - P. 43-78. DOI: 10.1007/978-3-319-17118-0_3
26. Tashkinov M. Statistical characteristics of structural stochastic stress and strain fields in polydisperse heterogeneous solid media // Computational Materials Science. - 2014. - Vol. 94. - P. 4450. DOI: 10.1016/j.commatsci.2014.01.050
27. Tashkinov M. Probabilistic description of stochastic processes of structural failure in advanced polydisperse composites // Proceedings of 11th. World Congress on Computational Mechanics (WCCM XI) / Eds. Eugenio Onate, Xavier Oliver and Antonio Huerta. - Barcelona: CIMNE, 2014.
28. Tashkinov M., Mikhailova N., Wildemann V. Characterisation of Micro-Scale Mechanical Behaviour of Composites Using Stochastic Boundary Value Problem Solutions // Proceedings of 20th International Conference on Composite Materials. - Copenhagen: Denmark, 2015. - P. 3315-3322, available at: http://iccm20.org/fullpapers/file?f=6OwZadRnYW (accessed 22 August 2015).
29. Hybrid composites - Metallic and ceramic reinforcements influence on mechanical and wear behavior / G. Miranda, M. Buciumeanu, S. Madeira, O. Carvalho, D. Soares, F.S. Silva // Composites. Part B: Engineering. - 2015. - Vol. 74. - P. 153-165. DOI: 10.1016/j.compositesb.2015.01.007
30. Ташкинов М. Статистические характеристики полей напряжений и деформаций в компонентах композитов со сферическими включениями при различных видах макрооднородного напряженно-деформированного состояния // Решение инженерных задач на высокопроизводительном вычислительном комплексе Пермского национального исследовательского политехнического университета: моногр. / под ред. В.Я. Модорского. - Пермь: Изд-во Перм. нац. исслед. политехн. ун-та, 2014. - C. 172-192.
31. Ташкинов М.А. Применение системы WOLFRAM MATHEMATICA для вычисления многомерных интегралов в задачах стохастической механики // Высокопроизводительные параллельные вычисления на кластерных системах (HPC 2014): материалы XIV Междунар. конф. - Пермь: Изд-во Перм. нац. исслед. политехн. ун-та, 2014. - С. 432-439.
32. Ташкинов М.А., Михайлова Н.В. Многоточечные приближения высших порядков в краевой задаче упругости полидисперсных композитов со случайной структурой // Вестник Нижего-род. ун-та им. Н.И. Лобачевского. - 2011. - № 4-4. - С. 1799-1800.
179
Tashkinov M.A. /PNRPUMechanics Bulletin 3 (2015) 165-181
References
1. Liu K.C., Ghoshal А. Validity of random microstructures simulation in fiber-reinforced composite materials. Composites Part B: Engineering, 2014, vol. 57, pp. 56-70. DOI:
10.1016/j.compositesb.2013.08.006
2. Rasool A., Bohm H.J. Effects of particle shape on the macroscopic and microscopic linear behaviors of particle reinforced composites. International Journal of Engineering Science, 2012, vol. 58, pp. 21-34. DOI: 10.1016/j.ijengsci.2012.03.022
3. Mishnaevsky Jr.L., Derrien K., Baptiste D. Effect of microstructure of particle reinforced composites on the damage evolution: probabilistic and numerical analysis. Composites Science and Technolog., 2004, vol. 64, iss. 12, pp. 1805-1818. DOI: 10.1016/j.compscitech.2004.01.013
4. Matveeva A.Y., Pyrlin S.V., Ramos M.M.D., Bohm H.J., van Hattum F.W.J. Influence of waviness and curliness of fibres on mechanical properties of composites. Computational Materials Science, 2014, vol. 87, pp. 1-11. DOI: 10.1016/j.commatsci.2014.01.061
5. Yu M., Zhu P., Ma Y. Effects of particle clustering on the tensile properties and failure mechanisms of hollow spheres filled syntactic foams: A numerical investigation by microstructure based modeling. Materials & Design, 2013, vol. 47, pp. 80-89. DOI: 10.1016/j.matdes.2012.12.004
6. Ekici R., Kemal Apalak M., Yildirim M., Nair F. Effects of random particle dispersion and size on the indentation behavior of SiC particle reinforced metal matrix composites. Materials & Design, 2010, vol. 31, iss. 6, pp. 2818-2833. DOI: 10.1016/j.matdes.2010.01.001.
7. Ghossein E., Levesque M. Homogenization models for predicting local field statistics in ellipsoidal particles reinforced composites: Comparisons and validations. International Journal of Solids and Structures, 2015, vol. 58, pp. 91-105. DOI: 10.1016/j.ijsolstr.2014.12.021
8. Buryachenko V. Micromechanics of heterogeneous materials. New York: Springer, 2007. 686 p.
9. Vil’deman V.E., Sokolkin Yu.V., Tashkinov A.A. Mekhanika neuprugogo deformirovaniia i raz-rusheniia kompozitsionnykh materialov [Mechanics of nonlinear deformation and failure of composite materials]. Moscow: Nauka, 1997. 288 p.
10. Vanin G.A. Mikromekhanika kompozitsionnykh materialov [Micromechanics of composite materials]. Kiev: Naukova dumka, 1985. 302 p.
11. Volkov S.D., Stavrov V.P. Statisticheskaia mekhanika kompozitnykh materialov [Statistical mechanics of composite materials]. Minsk: Belorusskii gosudarstvennyi universitet, 1978. 208 p.
12. Lomakin V.A. Statisticheskie zadachi mekhaniki tverdykh deformiruemykh tel [Statistical Problems of the Mechanics of Solid Deformable Bodies]. Moscow: Nauka, 1970. 139 p.
13. Shermergor T.D. Teoriya uprugosti mikroneodnorodnyh sred [Elasticity theory of microinhomogeneous media]. Moscow: Nauka, 1976. 400 p.
14. Beran M.J. Statistical continuum theories. New-York: Wiley. Intersci. Publ., 1968. 493 p.
15. Lifshitz I.M., Rosenzweig L.N. K teorii uprugikh svoistv polikristallov [On the theory of the elastic properties of polycrystals]. Journal of Experiment and Theor. Physics, 1946, vol. 16, pp. 967-980.
16. Bolotin V.V., Moskalenko V.K. K raschetu makroskopicheskikh postoiannykh sil’no izotropnykh kompozitsionnykh materialov [Calculation of macroscopic permanent strongly isotropic composites]. Izvestiia Akademii naukSSSR. Mekhanika tverdogo tela, 1967, no. 3, pp. 106-111.
17. SokolkinYu.V., Tashkinov A.A. Mekhanika deformirovaniia i razrusheniia strukturno-neodnorodnykh tel [Mechanics of deformation and failure of structurally heterogeneous bodies]. Moscow: Nauka, 1984. 116 p.
18. Pan’kov A.A. Statisticheskaia mehanika piezokompozitov [Statistical mechanics of piezocomposites]. Permskii gosudarstvennyi tekhnicheskii universitet, 2009. 480 p.
19. Khoroshun L.P. Metody sluchainykh funkcii v zadachakh o makroskopicheskikh svoistvakh mikroneodnorodnykh sred [Methods of random functions in the problems of macroscopic properties of microinhomogeneous media]. Prikladnaia mekhanika, 1978, vol. 14, no. 2, pp. 3-17.
180
Ташкинов М.А. /Вестник ПНИПУ. Механика 3 (2015) 165-181
20. Saraev L.A., Glushchenkov V.S. Neuprugie svoistva mnogokomponentnykh kompozitov so sluchainoi strukturoi [Inelastic properties of multi-component composites with random structure]. Samarskii universitet, 2004. 163 p.
21. Tashkinov M.A., Vildeman V.E., Mikhailova N.V. Method of successive approximations in a stochastic boundary-value problem in the elasticity theory of structurally heterogeneous media. Composites: Mechanics, Computations, Applications, 2011, vol. 2, no. 1, pp. 21-37.
22. Torquato S. Random heterogenous materials, microstructure and macroscopic properties. Springer, 2001. 701 p.
23. Lee H., Gillman A.S., Matous K. Computing overall elastic constants of polydisperse particulate composites from microtomographic data. Journal of the Mechanics and Physics of Solids, 2011, vol. 59, no. 9, pp. 1838-1857. DOI: 10.1016/j.jmps.2011.05.010
24. Torquato S. Modeling of physical properties of composite materials. International Journal of Solids and Structures, 2000, vol. 37, pp. 411-422.
25. Tashkinov M.A. Methods of Stochastic Mechanics for Characterization of Deformation in Randomly Reinforced Composite Materials. Mechanics of Advanced Material. Eds. V.V. Silberschmidt, V.P. Matveenko. Springer, 2015, pp. 43-78. DOI: 10.1007/978-3-319-17118-0_3
26. Tashkinov M. Statistical characteristics of structural stochastic stress and strain fields in polydisperse heterogeneous solid media. Computational Materials Science, 2014, vol. 94, pp. 44-50. DOI: 10.1016/j.commatsci.2014.01.050
27. Tashkinov M. Probabilistic description of stochastic processes of structural failure in advanced polydisperse composites. Proceedings of 11th. World Congress on Computational Mechanics (WCCM XI). Eds. Eugenio Onate, Xavier Oliver and Antonio Huerta. Barcelona: CIMNE, 2014.
28. Tashkinov M., Mikhailova N., Wildemann V. Characterisation of Micro-Scale Mechanical Behaviour of Composites Using Stochastic Boundary Value Problem Solutions. Proceedings of 20th International Conference on Composite Materials. Copenhagen, Denmark, 2015, 19-24 July, pp. 3315-3322, available at: http://iccm20.org/fullpapers/file?f=6OwZadRnYW (accessed 22 August 2015).
29. Miranda G., Buciumeanu M., Madeira S., Carvalho O., Soares D., Silva F.S. Hybrid composites -Metallic and ceramic reinforcements influence on mechanical and wear behavior. Composites. Part B: Engineering, 2015, vol. 74, pp. 153-165. DOI: 10.1016/j.compositesb.2015.01.007
30. Tashkinov M. Statisticheskie kharakteristiki polei napriazhenii i deformatsii v komponentakh kompozitov so sfericheskimi vkliucheniiami pri razlichnykh vidakh makroodnorodnogo napriazhenno-deformirovannogo sostoianiia [Statistical characteristics of stress and strain fields in the components of the composite with spherical inclusions subjected to various macroscopic stress-strain state]. Reshenie inzhenernykh zadach na vysokoproizvoditel’nom vychislitel’nom komplekse Permskogo natsional’nogo issledovatel’skogo politekhnicheskogo universiteta. Ed. by V.Ia. Modorskii. Permskii natsionalnyi issledovatelskii politekhnicheskii universitet, 2014, pp. 172-192.
31. Tashkinov M.A. Primenenie sistemy WOLFRAM MATHEMATICA dlia vychisleniia mnogomernykh integralov v zadachakh stokhasticheskoi mekhaniki [Application of WOLFRAM MATHEMATICA for the calculation of multidimensional integrals in problems of stochastic mechanics]. Materialy XIV Mezhdunarodnoi konferentsii “Vysokoproizvoditel’nye parallel’nye vychisleniia na klasternykh sistemakh (HPC 2014)”. Permskii natsionalnyi issledovatelskii politekhnicheskii universitet, 2014, pp.432-439.
32. Tashkinov M.A., Mikhailova N.V. Mnogotochechnye priblizheniia vysshikh poriadkov v kraevoi zadache uprugosti polidispersnykh kompozitov so sluchainoi strukturoi [Multipoint higher-order approximation in the boundary value problem of elasticity polydisperse composites with random structure]. VestnikNizhegorodskogo universiteta imeni N.I. Lobachevskogo, 2011, vol. 4, no. 4, p. 1799-1800.
181
