Концентрация микродефектов вблизи трещины разрушения в полимерах и композитах на их основе Текст научной статьи по специальности «Физика»

УДК 22.251

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

© А. А. Валишин МГТУ им. Н.Э. Баумана, Москва, 105005, Россия

Статья является продолжением предыдущих работ автора [1—3], в которых рассмотрено формирование зоны вынужденной эластичности перед фронтом трещины разрушения в аморфных стеклообразных полимерах, кинетика разрушения слабых узлов несущего молекулярного каркаса, образование и накопление локальных микродефектов, названных дырками, их упругое взаимодействие. В настоящей статье показано, что взаимодействие дырок вызывает появление вокруг каждой дырки «атмосферы», состоящей из более мелких дырок, что фронт трещины является источником собственного упругого поля и дырки диффундируют навстречу фронту трещины, а также рассчитаны диффузионные потоки дырок.

Ключевые слова: микродефекты, трещина разрушения, полимеры, композиты.

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

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

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

Рассмотрим поведение дырок вблизи фронта трещины разрушения.

Фронт трещины, т. е. край ее клюва, как и изолированная дырка, является источником собственного упругого поля. Это поле нетрудно рассчитать, если представить фронт трещины как линию, вдоль которой непрерывно распределены элементарные источники поля — центры дилатации. Тогда тензор деформаций упругого поля фронта трещины примет вид

= £ J& - 3 тЧ *(--» )<*>■ (1)

—H ^ '

Здесь радиус-вектор r' = (--) отнесен к локальной системе координат, связанной с каждой точкой фронта трещины. Чтобы выполнить интегрирование, необходимо приведение к единой системе, привязанной к фронту. В формуле (1) vo — объем, в котором локализован элементарный акт разрушения (см. подробнее [3]); 8ik — символ Кронекера; H — полудлина фронта трещины; h (zo) — функция распределения центров дилатации на фронте трещины, так что h (zo )Az o — число таких центров на малом участке фронта Azo. Расписывая покомпонентно тензор деформаций, получаем

S11 = ^Г ff"Г — 3-г\ h (zo)dzo;

—H l r3 r5 J v ;

H

612 =■

3 Vo r -y

—H r

¡ ~-yh(zo)dzo;

813 =— ^j-^ h(zo )dzo,

I -.3 „5

н - н 3 ^ Н У ( - г0 ) и

823 =--Г- I —^--И(г0 И0;

4л _н г

н

н ( 1 -з / \2 ^

У0 (• 1 3 ( - г0) ,/ \ ,

833 = — I Г3--г1- ^) 0.

V

Примем для простоты, что центры дилатации распределяются по фронту трещины равномерно: И (г0 ) = 1/2 Н. Тогда, вычисляя интегралы, получаем

811 =-822 =- 4лН р4 '

8 £У. (2)

812 = . „ 4 ; (2) 4лН р4

812 = 823 = 833 = 0.

^2 2

х + у — расстояние от точки наблюдения до фронта

трещины. Из формул (2) ясно, что полученный тензор деформаций не зависит от координаты г, отсчитываемой вдоль фронта, т. е. собственное упругое поле фронта трещины обладает в направлении этой координаты трансляционной симметрией. Формулы (2) удобнее записать в полярных координатах в плоскости (х, у):

2 2 V х - У

811 =-822 =- 4лН р2

У0 Бт 2ф.

У0 СОБ 2ф

812 = 0 ТТ 2 ;

8лН р2

813 = 823 = 833 = 0

(ф — полярный угол). След этого тензора деформаций равен нулю, т. е. упругое поле фронта трещины чисто сдвиговое. Это поле обратно пропорционально квадрату расстояния от фронта, в то время как

поле изолированной дырки убывает как г [3], т. е. поле фронта трещины распространяется на значительно большие расстояния.

Теперь нетрудно найти энергию взаимодействия дырки, попавшей в поле фронта трещины, с этим полем, а также силу, действующую на дырку в этом поле:

3AGV2 1

ивз = - 0

вз 4

32л H р

^4cos2 2ф + sin2 2ф j;

F = &AGV0 cos ф; (3.

* я-л-2 тт2 5 ' (3)

H р „ _ SAgvq

sin ф

Г =---'

* 8л2Н2 р5 '

Fz = 0.

Здесь $ — объем дырки (см. [3]); AG — возмущение объемного модуля упругости при элементарном акте разрушения (см. [3]).

Из формул (3) следует, что в упругом поле фронта трещины на дырку действует сила притяжения, т. е. дырка, зародившаяся в этом поле или попавшая в него в процессе перемещения, не покидает поля. Таким образом, упругое взаимодействие дырки с фронтом трещины приводит к образованию «облака» дырок вблизи фронта.

От рассмотрения отдельных дырок перейдем к анализу их коллективного поведения. В реальной эластической зоне число дырок может быть значительным, поэтому оправданно их усредненное описание. Дырки — результат флуктуаций теплового движения в слабых местах эластической зоны (в слабых узлах несущего каркаса), это стабильные, долгоживущие образования, время существования которых определяется их размерами. Мелкие дырки могут «залечиваться», в то время как крупные, раз возникнув, не исчезают. С течением времени число дырок увеличивается до некоторого критического уровня — тогда наступает коллапс эластической зоны и она теряет устойчивость. В объеме эластической зоны дырки распределяются очень неравномерно: они группируются в скопления около крупных дырок и около фронта трещины.

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

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

В [1] установлено, что по окончании формирования эластической зоны, т. е. по завершении вынужденной высокоэластической ползучести, в ней устанавливается однородное напряженно-деформированное состояние. Это состояние характеризуется напряжением, равным пределу вынужденной эластичности, и постоянной деформацией. На фоне такого состояния возникают локальные возмущения, создаваемые появившимися дырками. Было показано, что сфера влияния дырки имеет радиус не более двух ее диаметров, т. е. возмущение локализуется вблизи дырки. Пока дырок мало, возмущенные участки разбросаны беспорядочными «пятнами» в объеме эластической зоны и не перекрываются. По мере роста числа дырок и их объединения в скопления число и размеры «пятен» увеличиваются и перед коллапсом вся эластическая зона покрывается «пятнами» упругого возмущения, а напряженно-деформированное состояние зоны становится очень сложным.

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

число дырок в малом объеме АV эластической зоны в произвольный момент. Интеграл по объему эластической зоны выражает полное число дырок в зоне:

|п (М, г= р V)Ф(г). (4)

Формула (4) связывает локальную концентрацию дырок п (М, г)

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

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

До (п, Т)= КТ 1пп + у(Т), (5)

где у(Т) — функция только температуры. Через концентрацию п

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

Д = + Мвз>

или

д = кт 1п п+у(т )-мав22. (6)

Применим к обеим частям формулы (6) операцию пространственного градиента V:

= +^вз. (7)

Второе слагаемое здесь с точностью до знака — это сила, действующая на дырку со стороны поля эластической зоны. Поскольку деформационное состояние зоны однородно, то Vuвз = 0. Это означает, что однородное поле эластической зоны никак не воздействует на дырки. Первое слагаемое — это тоже сила (с точностью до знака), специфическая сила, которая действует на частицы растворенного вещества со стороны частиц растворителя в процессе их теплового движения и называется диффузионной силой. В нашем случае это сила, действующая на дырки со стороны теплового движения атомов среды, их окружающей, т. е. некоторая средняя сила, определяемая «толчками» атомов «растворителя» и определяющая броуновское движение дырки. Диффузионная сила проявляется в растворе как осмотическое давление, которое пропорционально температуре и описывается законом Вант-Гоффа [13].

Равновесие дырок по отношению к их диффузионному перемещению определяется условием независимости химического потенци-

ала | от пространственных координат (однородность химического потенциала). В этом случае V| = 0, тогда из формулы (7) следует, что наряду со вторым слагаемым VuM равно нулю и первое слагаемое V|0, т. е. диффузионная сила также отсутствует. Отсюда V| = 0, т. е., как ясно из формулы (7), химический потенциал ансамбля дырок не зависит от пространственных координат. Это означает, что дырки распределены в объеме эластической зоны в среднем равномерно. В то же время число дырок увеличивается со временем. Из всего этого следует, что на начальном этапе образование дырок происходит равновесно по отношению к пространственному распределению дырок при сохранении пространственной однородности.

Рассмотрим поведение мелких дырок в поле крупной дырки, которая создает вокруг себя упругое поле [3]. Энергия взаимодействия дырки с этим полем

ивз = "Т^ MGvo2 (5^0

8л r

где r — расстояние от центра крупной дырки, куда помещено начало координат, до точки нахождения мелкой дырки. На мелкую дырку в этом поле действует сила

F = "VuM = "-^2 MGv2 (5«0 )2 Г-. (8)

4л r

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

Мелкие дырки имеют очень незначительный радиус влияния, и потому их можно считать невзаимодействующими, даже если их достаточно много. Тогда такую систему дырок можно рассматривать в приближении слабого твердого раствора. Химический потенциал «облака» мелких дырок

|(n, T, r) = KTlnn + 4 + v(T), (9)

r

где обозначено

a = —32 SAGv2 (5n0 ).

8 л

Равновесие «облака» мелких дырок определяется условием |(r ) = const. Тогда из (9) получаем равновесную плотность дырок в «облаке»:

а -6 -г

п(г) = п0 (Т)ект . (10)

Следовательно, крупная дырка, вовлекая в сферу своего влияния мелкие дырки, а также инициируя их дополнительное зарождение, создает вокруг себя некую «атмосферу», в которой с течением времени устанавливается равновесное распределение плотности и которая по мере удаления от центра притяжения становится все более «разреженной» в соответствии с формулой (10). Аналогичная «атмосфера» устанавливается и вблизи фронта трещины, только ее плотность убывает с расстоянием медленнее.

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

Если равновесие в «атмосфере» мелких дырок еще не установилось, то отличие от равновесного состояния характеризуется градиентом химического потенциала = V(Яо + Vuвз. Этот градиент вызывает диффузионный поток мелких дырок к крупным. Поток пропорционален градиенту ^, т. е. средней силе, действующей на дырку: сумме диффузионной и внешней сил, и направлен противоположно градиенту ^, т. е. в сторону результирующей силы. В равновесии, когда = 0, эти силы, будучи направлены в разные стороны, уравновешиваются.

Диффузионный поток определяется формулой

3 = пУ = -у^0 + ивз),

где V — скорость поступательного перемещения мелких дырок; у — коэффициент пропорциональности. Видно, что поток состоит из двух частей: диффузионного броуновского потока

31 =-У^0

и вынужденного потока

32 =-У^вз.

и V.

В связи с этим и скорость дырок V разбивается на две скорости V .

Диффузионный поток

3 = nV1 = ^ц0 =-у^° Vn = -DVn, (11)

дп

где коэффициент диффузии О определяется соотношением

д|- КТ

О = у—- = у-.

дп п

Здесь использована формула (5). Вынужденный поток

/ = пУ2 = -уУмВз =ур = -2а у 4-, (12)

г

где обозначено:

а = АСу2(5по)2. 8л

Здесь использована формула (8) для силы, действующей на мелкую дырку в поле крупной. Отсюда вынужденная скорость

V = = -1 У^ =-Кз,

п П

где константа Ъ = у/п называется вынужденной подвижностью дырки. Результирующий поток дырок

/ = / + / =-О Уп - 2а у Г = -О Уп + у ^. (13)

г

Из формул (11), (12) ясно, что потоки / и ./2, а также скорости

VI и V направлены в противоположные стороны. Диффузионный

поток /1 направлен противоположно градиенту химического потенциала, т. е. в сторону убывания концентрации дырок от притягивающего центра, а вынужденный поток /2, наоборот, направлен к этому

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

Сравнив выражения для коэффициента диффузии и подвижности, получим О = Ъ КТ. Это аналог известной в теории броуновского движения формулы Эйнштейна, связывающей два кинетических коэффициента — коэффициент диффузии и подвижность дырок.

Величину, обратную подвижности, называют коэффициентом внутреннего трения:

, 1 кт

к = — =-.

Ь Б

Пространственное и временное распределение дырок в неравновесной «атмосфере» описывается уравнением неразрывности

д п ^

— + 3 = g(г, X), (14)

дХ

где g (г, X) — функция источника дырок. Смысл ее состоит в том, что величина g (г, X)ДVЛX определяет число дырок, зародившихся

в малом объеме Л V за малое время ДХ. Интегрируя эту функцию по объему эластической зоны, получаем число дырок, зародившихся в зоне за время ДХ:

{ g (г, X) dV = ^ V )ф( X )Лх ,

где функция ф(х) определена в работе [2]. Интегрируя еще раз, теперь по времени, получим полное число дырок во всей эластической зоне в момент X:

X

(г, X^Х = р V)Ф(Х).

0

Подставляя в уравнение (14) поток 3 из формулы (13), после преобразований получаем

дп 10аБ -8 ч

— = БДп-п-г 8 + g(г, X).

дХ КТ У }

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

Итак, около каждой дырки и около фронта трещины имеется своя «атмосфера», состоящая из более мелких дырок различной плотности и протяженности. В целом в эластической зоне устанавливается некая иерархия: дырка со своей «атмосферой» входит в состав «атмосферы» еще более крупной дырки, эта система является частью еще большей системы и т. д. В поле фронта трещины и в поле крупных дырок или их скоплений возникают диффузионные потоки дырок. Требуется рассмотреть заключительный этап эволюции зоны вынужденной эластичности, заканчивающийся ее коллапсом. Ключевой вопрос здесь — какова критическая концентрация дырок, т. е. в какой момент наступает коллапс. Но для этого потребуется другой подход.

ЛИТЕРАТУРА

[1] Валишин А.А., Степанова Т.С. Особенности квазихрупкого разрушения полимеров и композитов на их основе. Инженерный журнал: наука и инновации, 2012, вып. 2 (2). doi: 10.18698/2308-6033-2012-2-52

[2] Валишин А.А., Степанова Т.С. Кинетика зарождения локальных микродефектов при квазихрупком разрушении полимеров и композитов на их основе. Инженерный журнал: наука и инновации, 2013, вып. 9 (21). doi: 10.18698/2308-6033-2013-9-1119.

[3] Валишин А.А., Миронова Т.С. Силовые упругие поля локальных микродефектов в напряженных полимерах и композитах на их основе. Инженерный журнал: наука и инновации, 2014, вып. 8 (32). doi: 10.18698/2308-60332014-8-1241.

[4] Димитриенко Ю.И., Сборщиков С.В., Соколов А.П., Шпакова Ю.В. Численное моделирование процессов разрушения тканевых композитов. Вычислительная механика сплошной среды, 2013, т. 6, № 4, с.389-402. DOI: 10.7242/1999-6691/2013.6.4.43.

[5] Димитриенко Ю.И., Соколов А.П. Многомасштабное моделирование упругих композиционных материалов. Математическое моделирование, 2012, т. 24, № 5, с. 3-20.

[6] Dimitrienko Y.I., Sokolov A.P. Elastic Properties of Composite Materials. Mathematical Models and Computer Simulations, 2010, vol. 2, no. 1, pp. 116-130.

[7] Dimitrienko Yu.I. Thermal Stresses and Heat Mass-Transfer in Ablating Composite Materials. Int. J. of Heat Mass Transfer, 1995, vol. 38, no. 1, pp. 139-146.

[8] Dimitrienko Yu.I. Thermal Stresses in Ablative Composite Thin-Walled Structures under Intensive Heat Flows. Int. J. of Engineering Science, 1997, vol. 35, no. 1, pp. 15-31.

[9] Looyehl M.R.E., Samanta A., Jihan S. Modeling of Reinforced Polymer Composites Subject to Thermo-Mechanical Loading. Int. J. for Numerical Methods in Engineering, 2005, v. 63, no. 6, pp. 898-925.

[10] MeManns H.N., Springer G.S. Hugh Temperature Thermomechanical Behavior of Carbon-Phenolic Composites: I Analysis, II Results. J. Composite Materials, 1992, vol. 26, pp. 206-255.

[11] Baia Yu, Valleea Till, Keller T. Modeling of Thermal Responses for FRP Composites under Elevated and High Temperatures. Composites Science and Technology, 2008, v. 68, no. 1, pp. 47-56.

[12] Димитриенко Ю.И. Нелинейная механика сплошной среды. Москва: Физматлит, 2009, 624 с.

[13] Гуревич Л.Э. Основы физической кинетики. Москва, ГИТТЛ, 1940, с. 14.

Статья поступила в редакцию 27.06. 2015

Ссылку на эту статью просим оформлять следующим образом:

Валишин А.А. Концентрация микродефектов вблизи трещины разрушения в полимерах и композитах на их основе. Инженерный журнал: наука и инновации. 2015, вып. 6.

URL: http://engjournal.ru/catalog/msm/pmcm/1409.html

Валишин Анатолий Анатольевич — д-р физ.-мат. наук, профессор кафедры «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана, профессор кафедры «Высшая и прикладная математика» Московского государственного университета тонких химических технологий им. М.В. Ломоносова. e-mail: enf@mail.ru

A.A. BanuMun

The concentration of microdefects near the crack fracture in polymers and composites based on them

© A.A. Valishin

The article continues and develops previous studies [1—3], which described the formation of zones offorced elasticity in front of the crack fracture in amorphous glassy polymers, the kinetics of destruction of weak nodes carrying the molecular skeleton formation and accumulation of local micro-defects, called holes, their elastic interaction. In this paper, it is shown that the interaction of holes leads to the fact that each hole is surrounded by "atmosphere" of smaller holes. It is also shown that the crack front is the source of its own elastic field. It is shown that the holes diffuse toward the front of the crack. Diffusive flux of holes was designed.

Keywords: microdefects, crack destruction, polymers, composites.

REFERENS

[1] Valishin A.A., Stepanova T.S. Inzhenernyi zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovation, 2012, issue no. 2 (2). DOI: 10.18698/2308-6033-2012-2-52.

[2] Valishin A.A., Stepanova T.S. Inzhenernyi zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovation, 2013, iss. 9 (21). DOI: 10.18698/2308-6033-2013-9-1119.

[3] Valishin A.A., Mironov T.S. Inzhenernyi zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovation, 2014, issue 8 (32). DOI: 10.18698/2308-6033-2014-8-1241.

[4] Dimitrienko Yu.I., Sborschikov S.V., Sokolov A.P., Shpakova Yu.V. Vychislitel'naya mekhanika sploshnoi sredy — Computational Continuum Mechanics, 2013, vol. 6, no. 4. pp. 389-402. DOI: 10.7242/19996691/2013.6.4.43.

[5] Dimitrienko Yu. I. ,. Sokolov A.P. Matematicheskoe Modelirovanie — Mathematical Models and Computer Simulations, 2012, vol. 24, no. 5, pp. 3-20.

[6] Dimitrienko Y. I., Sokolov A.P. Elastic Properties of Composite Materials.

Mathematical Models and Computer Simulations, 2010, vol. 2, no. 1, pp. 116130.

[7] Dimitrienko Yu.I. Thermal Stresses and Heat Mass-Transfer in Ablating Composite Materials. Int. J. of Heat Mass Transfer, 1995, vol. 38, no. 1, pp. 139-146.

[8] Dimitrienko Yu.I. Thermal Stresses in Ablative Composite Thin-Walled Structures under Intensive Heat Flows. Int. J. of Engineering Science, 1997, vol. 35, no. 1, pp. 15-31.

[9] Looyehl M.R.E., Samanta A., Jihan S. Modeling of Reinforced Polymer Composites Subject to Thermo-Mechanical Loading. Int. J. for Numerical Methods in Engineering, 2005, v. 63, no. 6, pp. 898-925.

[10] MeManns H.N., Springer G.S. Hugh Temperature Thermomechanical Behavior of Carbon-Phenolic Composites: I Analysis, II Results. J. Composite Materials, 1992, vol. 26, pp. 206-255.

[11] Baia Yu, Valleea Till, Keller T. Modeling of Thermal Responses for FRP Composites under Elevated and High Temperatures. Composites Science and Technology, 2008, v. 68, no. 1, pp. 47-56.

[12] Dimitrienko Yu.I. Nelineinaya mekhanika sploshnoi sredy [Nonlinear Continuum Mechanics]. Moscow, Fizmatlit Publ., 2009, 624 p.

[13] Gurevich L.E. Osnovy fizicheskoi kinetiki [Fundamentals of physical kinetics]. Moscow, GITTL , 1940, pp. 14.

Valishin A.A., Dr. Sci. (Phys.-Math.), professor of the Computational Mathematics and Mathematical Physics Department at the Bauman Moscow State Technical University; professor of the Higher and Applied Mathematics Department at Lomonosov Moscow State University of Fine Chemical Technologies. e-mail: enf@mail.ru