CC BY
56
НАУЧНОЕ ИЗДАНИЕ МЕТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 - 48211. Государственная регистрация №0421200025. КБМ 1994-0408
электронный научно-технический журнал
Моделирование нестационарного внутреннего теломассопереноса в теплозащитных кострукциях на основе трехмерного конечно-элементного анализа
# 10, октябрь 2013 Б01: 10.7463/1013.0606069
Димитриенко Ю. И., Коряков М. Н., Балакшин А. И.
УДК 539.3
Россия, МГТУ им. Н.Э. Баумана [email protected] [email protected] [email protected]
Полимерные композитные материалы до сих пор являются одним из основных типов материалов, применяемых в теплозащитных системах высокоскоростных летательных аппаратов, движущихся в плотных слоях атмосферы Земли, а также других планет (Марса, Венеры). В работах [1-13] была разработана математическая модель процессов внутреннего тепломассопереноса в полимерных теплозащитных материалах при нестационарном нагреве, учитывающая процессы термодеструкции, газообразования и фильтрации газов в порах материала. Модель была применена для расчета одномерных процессов теплопереноса в элементах конструкций с относительно простой геометрической формой: пластинах, цилиндрических оболочках [11-13]. Вопросы
моделирования внутреннего тепломассопереноса в композиционных материалах при высоких температурах рассматривались также в современных работах зарубежных авторов [14-26]. В настоящей работе предложен метод численного конечно-элементного расчета процессов внутреннего теплопереноса в трехмерных элементах теплозащитных конструкций из композитов.
Система трехмерных уравнений внутреннего тепломассопереноса в теплозащитных конструкциях. Под действием высоких температур, обычно свыше 300оС, в полимерных композиционных материалах происходят внутренние физикохимические превращения - процессы термодеструкции. В результате композит при нагреве представляет собой многофазную систему, в работах [1, 2] была предложена модель 4-х фазной системы, в которой 1-я фаза — наполнитель (волокна, армирующие
частицы), 2-я фаза — полимерная матрица; 3-я фаза - твердая пиролитическая фаза, образующаяся из полимерной при высоких температурах, 4-я фаза - газ, также образующийся вследствие термодеструкции исходной фазы 2 и находящийся в порах.
Обозначим объемные концентрации фаз как щ., г = 1...4. Эти концентрации фаз (рх, г = 1,4 постоянно меняются при нагреве, что приводит к изменению характеристик всего материала: коэффициентов теплопроводности, плотности, теплоемкости. Система
уравнений изменения концентраций фаз, а также уравнение энергии всей 4-х фазной системы, образует систему уравнений внутреннего теплопереноса в композите, и записывается в тензорной форме следующим образом [1]:
где р , щ , сг , г = 1,4 - плотности , объёмные доли и удельные теплоёмкости фаз, в -температура, общая для всех фаз, t - время, х - радиус-вектор, V - область, занятая композитной конструкцией, V - набла -оператор [27], J - скорости массопереноса (термодеструкции), АН - теплота термодеструкции, для которых имеют место соотношения:
Здесь р - плотность всего композита, p - поровое давление газообразных продуктов термодеструкции, с - удельная теплоёмкость , ц - коэффициент вязкости газа, E -энергия активации терморазложения, R - универсальная газовая постоянная, Г -коэффициент газификации. Плотности всех твердых фаз полагаются постоянными: Р = const, i = 1,2,3. Тензор коэффициентов газопроницаемости K и тензор теплопроводности X композита являются функциями от объемных долей фаз и температуры:
dt 4 '
дй
рс— = V-(X-V#) + ^4c4V#-K-Vp - J Ah; x eV
d^P = V- ( K-Vp ) + Jr;
dt
4У4 _
(1)
pc = ^Lpc% , 1 = Z^-, р = Yр, % = 021 -<?2)(1 -r)—,
(2)
K = K E exp(n^4/3 ^ X = X(^1, ^2, ^3, й),
(3)
где К0 , пк - константы, а Е - метрический тензор. Конкретные выражения для функций X = X (щ1, щ2, щ3, в) приведены в [1 ].
На нагреваемой части внешней границы Е1 области V заданы тепловой поток и давление, на части поверхности Е2 заданы температура и давление, а на оставшейся части
Е2 заданы условия герметичности и теплоизоляции:
х єЕ1: —п 'Л -Ув = qe = а(ве — в) — єав4 , р = ре (4)
X £^2 : в = ве р = ре , (5)
х єЕ3 ,: п • Ур = 0; п- X- У0 = 0, (6)
где в температура поверхности композитной конструкции, ве - температура внешней поверхности пограничного слоя потока, обтекающего конструкцию, р = ре - давление газового потока на конструкцию, £6 - интегральный коэффициент излучения поверхности конструкции, <г - коэффициент Стефана-Больцмана. Унос материала конструкции не учитывался. Вся поверхность композита состоит из 3-х частей: Е = Е1 +Е2 +Е3.
Начальные условия к системе (1)-(3) имеют вид
/ = 0: Щ2 = Щ2 , Р = P0, в = в0. (7)
В итоге получаем постановку задачи внутреннего тепломассопереноса (1)-(7) относительно 3-х неизвестных скалярных функций: р,в, щ2, зависящих от х, ^.
Вариационная постановка задачи внутреннего тепломассопереноса. Дадим вариационную постановку задачи внутреннего тепломассопереноса (1)-(7). Введем вариации 3-х неизвестных функций: 8(р2, 8р и 8в, которые удовлетворяют граничным
условиям 1-го рода на внешних поверхностях:
х є£1: 8р = 0 , (8)
х єЕ2 : 80 = 0 8р = 0, (9)
Домножим первое уравнение системы (1) на 8(р2, второе на 8р, а третье -80, и проинтегрируем получившиеся уравнения по области V, в результате получим:
|р28(р2 ^ йУ = -18ср2МУ, (10)
V д У
\8р др4р4 =| 8 рУ • (К-Ур)йУ + j8pJГdV (11)
У д у
лл
І рс8в—йУ = \бв7- (X • V в)йУ + | Зв(р4е4у в- (К • Чр)с1У -18вЗ ШУ, (12)
Преобразуем в этих уравнениях интегралы, содержащие ковариантные производные, следующим образом:
\PpV-(К-Vp)dV = |У-(р?К-Vp)dV-\VPp-К-VрйУ =
V V V
= Г п-рК-VpdЕ-fVp- К-^8р<йУ = -fVp- К-^8р<йУ , (13)
Ех +^2 +Е3 V V
$Р№- (X - V6)dV = Г п - (Рех - Vв)dV - Г VP6 - (X - V6)dV =
V V V
= Г PеqedV - Г VPе - (X - Vе)dV,
Ех V
здесь использована формула Гаусса-Остроградского, а также граничные условия (4)-(6),
(8), (9). Подставляя формулы (13) в (11) и (12), получаем искомую вариационную
остановку задачи внутреннего тепломассопереноса
|рР«>2 ^V = -\P,f2_JdУ, (14)
V д V
р дР4^1 = -^ - К- VРpdV + Г Р pJГdV (15)
V д V
Л Л
Г pcРе—dV = Г РеqedV - Г VРе - (X - Vе)dV +1 Реср4е^е - (К - Vp)dV -1 РеJ£М¥, (16)
Производная по времени в (15) преобразуется следующим образом:
др4^>4 _ ц др цр дв
(17)
ді Яв ді Яв2 ді
Конечно-элементная аппроксимация задачи тепломассопереноса. Согласно общей концепции метода конечных элементов [28-30], область У решения вариационных уравнений (14)-(16) разобьем на систему большого числа подобластей (конечных элементов) Уе , в каждом из которых решение системы (14)-(16) ищем в виде аппроксимаций по узловым функциям температуры 3, порового давления у и концентрации полимерной фазы и (значения функций в узлах КЭ):
(3.Л
в={фв}’ 3}=к ... Фв)
1хп Пх1
Р = {ФР Г-М ’ ^ ={Ф}Т• {и} , (18)
}Т-
1хп пх1 1хп пх1
где {Фе},{Фp},{Ф} - координатные столбцы функций форм КЭ [29], все имеющие
одинаковую структуру {ф ^, £2, Ь3, Ь4 ) и выражающиеся через барицентрические
координаты КЭ Li. Тогда вариации неизвестных функций имеют вид:
Ре = {Фе}TР{3}, Рp = {Фр}TР{у}, Рср2 = {Ф}Р{и}. (19)
Градиенты от неизвестных функций при конечно-элементной аппроксимации
вычисляются следующим образом:
Vp = {Цф } -{у} = [ Вг ]•{у}, (20)
тхп пх1
где {L} - координатный столбец дифференциальных операторов, m - размерность
шх\
пространства, а \Bp} ]={L}&{фp} - матрица, полученная внешним умножением
p J
mхn mх1 1хп
координатных столбцов. Аналогично вычисляем
vе = Щф,}-{5} = [вв}{5}, (21)
тхп Пх1
V
mхn mх1 1хп
где ^{ф^у.
Подставляя выражения (18)-(21) в систему (14)-(16), имеем:
|р2У{и}Т • {Ф} ® {Ф}Т д{u}dV = -|{Ф}Т У{и^^, (22)
ТЛ ді ТҐ
у {ф мф }т цв=$ -в д {фp ^ ^ д д -
-| { у} • [ в, ]Т • [ K ] • [ вp ] - {у у} аУ+1зг [фг }т {у у} №
(23)
І рс {У3 } {Ф,} в {ф, }Т ] ^ к’ = -/ {3 }• [В, Г • Ы • [в в ] • {У3 № -1Ч, {Фв }Т {У3 № +
У І ді д У е2
|У3}Т {Фв^сІУ (24)
+
У
где обозначена функция
¥ = (Р4 c4V , - К ^ - J АИ = у4с4 {3}-[Бе] [К НВгИ.у}-J АИ (25)
Вынося вариации неизвестных функций за скобки, преобразуем эту систему к виду
[M, -[S ]д-дг+[K-]•< y} -[ Rp]•(">=0
[М«]'д{Г+[*»Ш + №1 = 0 (26)
dt
[M ]+ [ R] • (и) = 0
dt
где обозначены матрицы и координатные столбцы:
[М ]=/К}®К)Г R,dV ■ [Me] = S&r }®-»Ф - YpcdV , (27)
[M] = р\{Ф}»{Ф}Т dV, [S]=1{Ф-}® {ф,}Т R,dV,
\кг ]=№- ]• K \b, ]iV , [K^fiBj •[№,№ ,
V V
[R] = J0 f exp (--^1{Ф }®{Ф}Т dV , [ R- ] = rJo f exp ["--§){Ф- }®{Ф}Т dV
V V R,J nxn V V R,J nxn
(R,) = f q, {Ф,}d 1-\{Ф,}^У .
I2 V
Метод численного решения системы уравнений. Полученная система (26) представляет собой нелинейную систему обыкновенных дифференциальных уравнений
относительно узловых функций {у), (5) и (и) . Все матрицы [Mp ] -, \ s ],[ к„ ],
\Me], \Ke], [M] и столбцы правых частей (R) , (Rp), (R,) зависят от неизвестных
функций. Для решения этой системы применим пошаговый метод линеаризации с использованием неявной схемы, в результате приходим к следующей последовательности системы линейных уравнений:
([Mp ]m + At [Kp ]m) - {у)m+1 - [S]m - {3)m+1 - [Rp ]m - (uu)m+1 = [Mp ]m - {y)m - [S]m -{3)m,
([Me]m + At\K,]m)-(ЗГ1 =[M,]mW -(R,)m, (28)
([M]m + At[R]m) - (u)m+1 = [M]m - (u)m
где индекс m - означает номер шага по времени. Система (28) является линейной
относительно неизвестных функций (y)m+1, (5)m+1 и (u)m+1, для ее решения применяем
методы решения больших систем линейных алгебраических уравнений.
Линеаризованная система (28) может быть записана в обобщенном матричном виде
[ к ]т -иг1 — {/}
(29)
где введена обобщенная матрица жесткости
[К ]т —
[М, ]т +м [кг ]т -[Х]т -[Л, ]т
0 [М.Г +ы [кв]' 0
0 0 [М ]т + АГ[ Л]т)
(30)
обобщенный столбец неизвестных ^}т+ — {{у}т+ {$}т+ {п}т+ } и {/}т - столбец правых частей.
Для всех функций будем рассматривать линейную аппроксимацию с использованием 4-х узлового КЭ, для него барицентрические координаты Ь, имеют вид [29]:
Ь = ^—(а, + Ъх + с. у + diz), г — 1,4
I у-тг ^ г г г ^
ЬУ
(31)
У = 1
Ь
1 х у ^
1 Х2 У 2 ^
1 Х3 Уз z3
х2 У 2 1 У 2 1 х2 Z2 1 х2 У 2
а1 = х3 У3 zз , Ъ = - 1 У3 Z3 , С1 = 1 х3 Z3 , — — 1 х3 У3
х4 У 4 Z4 1 У4 Z4 1 х4 Z4 1 х4 У 4
1 Х4 У4 ^
(остальные коэффициенты вычисляются аналогично, по циклическим перестановкам
узлов с заменой знаков на противоположный: (х2, х3, х4),( х3, х4, х1),( х4, х1, х2),( х1, х2, х3)).
Матрица [ В, ] вычисляется следующим образом:
Ва —
дЬ1 дЬ2 дЬ3 дЬ4
дх дх дх дх
дЬ1 дЬ2 дЬ3 дЬ4
ду ду ду ду
дЬ1 дЬ2 дЬ3 дЬ4
1_
ЬУ
Ъ1 Ъ2 СО Ъ4
С1 С2 С3 С4
С1 2 ^3 С3
(32)
дz дz дz дz
Тогда матрицы [К, ] и [Кв] могут быть вычислены без расчета интегралов по КЭ:
[К,] — У [В ,Г ‘К'[В ,] , [К®] = У [Вв] [Я] [Вв]. Для вычисления матриц [М,],
[Ме],[М ],[ £ ], [Л] воспользуемся формулами [29]
а\Ъ\С\С\
Л/ ь°А ьз Ь4 су
(а + Ъ + с + С + 3)\
ЬУ
(33)
г У с 2 У г У
/ ЬСУ — -, / Ь СУ — —, / ЬЬСУ ——
I 4 I 10 I 1 20
т V Р
Т огда, вычисляя переменные — , —2
Яв Яв
и аналогичные другие, участвующие в
выражениях (27) в матрицах, как средние по 4-м узловым значениям, получаем:
[М ] — РрУ I', [М ] ——I', [^]— РУ I',
1 J 20 [ Р ] 20Яв 20Яв2
где обозначены следующие величины
I'
2 1 1 1
1 2 1 1
1 1 2 1
1 1 1 2
Е' —
1
1
1
1
р С — Т £ Р°С С
4 г—1
[Мв ] — -рСУ I', [ Я] — — I', (34)
1 в1 20 20
1/ в — £1/4,9,. , 4 — ^У,
^ г в2 41—1А2’
г —1
V ЯВ у
(35)
Алгоритм конечно-элементного метода решения задачи внутреннего тепломассопереноса. Общая схема алгоритма конечно-элементного метода решения задачи внутреннего тепломассопереноса состоит из следующих этапов:
1 .Разбивка области на КЭ.
2.Нумерация: границ области; узлов на границе; КЭ; узлов КЭ; узловых степеней свободы.
З.Задание свойств материалов: констант рр,рЪ,Г,Гг,Ае,Ае^,С? и других.
4.Задание граничных условий.
5. Задание начальных условий.
Ь. Организация итерационного процесса решения задачи, на каждой итерации:
6.1. вычисляются локальные матрицы «жесткости» и матрицы «масс» конечных элементов, для этого используются значения температуры и давления для конечного элемента из предыдущей итерации (на 1 - начальные значения);
6.2 формируются на их основе глобальная матрица «жёсткости» и векторы «узловых нагрузок»;
6.4. формируется разрешающая СЛАУ с учетом граничных условий;
6.5. осуществляется решение СЛАУ методом локального оптимума.
7. Проверяется условие окончания процесса нагрева конструкции.
8. Осуществляется вывод и обработка результатов расчётов.
Хотя все матрицы (27) являются симметричными, но несимметричной является обобщенная матрица линеаризованной системы (30), поэтому для решения системы (29) использовались методы бисопряженного градиента.
Г енерирование трехмерной КЭ сетки осуществлялось с помощью функций свободно распространяемой библиотеки tetgen. Программная реализация разработанного конечноэлементного метода решения задачи осуществлялась в программной среде УБ-2010. Был написан специализированный программный модуль БетИМТСот, реализующий конечно-элементный метод решения задачи.
Тестовые расчеты. Для тестирования разработанного программного обеспечения использовалось решение одномерной задачи нестационарной теплопроводности об одностороннем нагреве цилиндра
дв_ — дв
дt ~а дг2’ "2’ ' (36)
в I „ = вт, в I „ = в„ Т в I. = вп .
— а~, Я1 < г < ^ 0 < t < Тт.
со следующими значениями параметров: а — X/ рС — 0.3 -104 м2 / с, в0 — вЯ1 — 300^, вЛ2 — 1000^, Я1 — 0,1м, Л1 — 0,2м, Ттах — 50с . Число КЭ в сетке, построенной для тестовых расчетов составляло 5095, число узлов в сетке - 1137. Тестирование проводилось путем сравнения расчетов на основе метода МКЭ и расчетов, полученных с помощью неявной конечно-разностной схемы, для решения которой использовался метод прогонки. На рисунке 1 показаны распределения температуры по толщине цилиндра, полученные с помощью обоих методов. Как видно, расчетные кривые, полученные по обоим методам достаточно хорошо совпадают, волнистый характер графиков температуры, полученных с помощью МКЭ, обусловлен первым порядком аппроксимации температуры в конечном элементе и достаточно грубой КЭ-сеткой. При использовании более мелких сеток характер кривых температуры становится более гладким, однако при этом существенно возрастает время расчета задачи, что связано с трехмерностью задачи.
,1000
е,к
900 800 700 600 . 500 400 300
001 0 02 0.03 0 04 0 05 0.06 0.07 0.08 0.09
Г, м
Рис.1 Результаты тестовых расчетов температуры в цилиндре (гладкие кривые - расчет по конечно-разностному методу, волнистые - МКЭ-расчет), цифры - у кривых - время
нагрева 1 ( с).
Численное конечно-элементное моделирование тепломассопереноса в теплозащитной конструкции. МКЭ расчеты проводились для трехмерной цилиндрической конструкции из полимерного композиционного материала при внешнем нагреве. Геометрические параметры конструкции были следующими: длина цилиндра Ь=0,15 м, внутренний радиус Ю=0,12 м, толщина цилиндра Ь=0,02 м. КЭ сетка состояла из 49160 элементов и 11848 узлов, элемент использованной КЭ сетки для конструкции показан рис. 2. Характеристики полимерного композиционного материала (был выбран стеклопластик на эпоксифенольном связующем с углеродными волокнами) взяты из работы [1], температура внешнего потока ве = 1000К, давление внешнего потока ре = 1атм .
Результаты расчетов показаны на рисунках 3...5. На рисунке 3 приведены распределения температурного поля в(г, ^) по толщине цилиндра (г - радиус) для различных моментов времени 1 прогрева. С увеличением времени нагрева конструкции происходит возрастание температуры как на внешней нагреваемой поверхности, так и внутри самого материала. На рисунке 4 показаны распределения концентрации полимерной фазы композита (р2 по толщине конструкции. Поскольку на внешней поверхности и в подповерхностном слое конструкции температура композита превышает температуру начала интенсивной термодеструкции полимерной матрицы (примерно в«600К), то в зоне максимального нагрева реализуется процесс термодеструкции матрицы и содержание полимерной фазы снижается. С увеличением времени нагрева зона уменьшенного содержания полимерной фазы ср2 продвигается внутрь конструкции.
На рисунке 5 показано распределение порового давления р газообразных продуктов термодеструкции матрицы для разных моментов времени нагрева. Ввиду того, что в подповерхностном слое одновременно с уменьшением доли полимерной фазы увеличивается пористость композита, то в этой зоне резко возрастает коэффициент газопроницаемости композита. Вследствие этого график распределения порового давления имеет характерный профиль с наличием локального экстремума внутри конструкции. Этот максимум по мере прогрева материала смещается внутрь конструкции. Величина порового давления весьма значительна, в проведенных расчетах максимум давления составил 60 атм, при таком значении давлении возможно расслоение композитов, обусловленное повышенным газовыделением вследствие термодеструкции. Аналогичные результаты были получены и протестированы в предыдущих работах, с использованием конечно-разностных одномерных и двумерных схем. Результаты настоящей работы, позволяют значительно расширить возможности численного моделирования нестационарного прогрева теплозащитных конструкций за счет исследования сложных трехмерных конструкций и условий нагрева.
Выводы. Предложен метод численного конечно-элементного расчета процессов внутреннего теплопереноса с учетом процессов терморазложения и фильтрации газов в порах применительно к трехмерным элементам теплозащитных конструкций из композитов. Разработан специализированный программный модуль, реализующий численный конечно-элементный метод решения задачи внутреннего тепломассопереноса в трехмерных элементах конструкций из композитных термодеструктирующих материалов. Проведено тестирование разработанного метода и программного модуля, показавшее его работоспособность и применимость для решения задач исследования нестационарных режимов нагрева теплозащитных конструкций.
Рис. 2. Элемент КЭ сетки, использованной для моделирования тепломассопереноса в
цилиндре
е-е0,к 400
300
200
100
1 п 'у/И\
Ш)
0.1 0.105 0.11 0.115 г,м
Рис. 3. Распределение температуры в(г, ^) по толщине цилиндра.
ф2 0,25
0,2
0,15
ОД
0,05
10Д
0.1 0.105 0.11 0.115 г,м
Рис. 4. Распределение концентрации полимерной фазы ср2 по толщине цилиндра.
Ftt/4
,3
Ц 2y
\
о —--------
0.1 0.105 0.11 0.115 0.12
г, м
Рис. 5. Распределение порового давления p по толщине цилиндра.
Исследование выполнено при поддержке Министерства образования и науки Российской Федерации (Соглашения №14.B37.21.0448, 14.B37.21.1869, НИР № 1.5433.2011), грантов Президента РФ (МК-6421-2012-9, МК-765.2012.8, МК-6573.2013.3) и гранта РФФИ №1208-00998.
Список литературы
1. Димитриенко Ю.И. Механика композиционных материалов при высоких температурах. М.: Машиностроение, 1997. 368 c.
2. Dimitrienko Yu.I. Thermal stresses and heat mass-transfer in ablating composite materials // Int. Journ. of Heat Mass Transfer. 1995. Vol. 38, no. 1. P. 139-146.
3. Dimitrienko Yu.I. Thermal Stresses in Ablative Composite Thin-Walled Structures under Intensive Heat Flows // Int. Journal of Engineering Science. 1997. Vol. 35, no. 1. P. 15-31.
4. Dimitrienko Yu.I. Heat- Mass-Transport and Thermal Stresses in Porous Charring Materials // Journal of Transport in Porous Media. 1997. Vol. 27, no. 2. P. 143-170.
5. Dimitrienko Yu.I. Modelling of Mechanical Properties of Composite Materials under High Temperatures. Part 1. Matrix and Fibres // Int. Journal of Applied Composite Materials. 1997. Vol. 4, no. 4. P. 219-237.
6. Dimitrienko Yu.I. Modelling of Mechanical Properties of Composite Materials under High Temperatures. Part 2. Unidirectional Composites // Int. Journal of Applied Composite Materials. 1997. Vol. 4, no. 4. P. 239-261.
7. Dimitrienko Yu.I. Thermomechanical Behaviour of Composite Materials and Structures under High Temperatures. Part 1. Materials Composites. Part A // Applied Science and Manufacturing. 1997. Vol. 28A. P. 453-461.
8. Dimitrienko Yu.I. Thermomechanical Behaviour of Composite Materials and Structures under High Temperatures. Part 2. Structures. Composites. Part A // Applied Science and Manufacturing. 1997. Vol. 28A. P. 463-471.
9. Dimitrienko Yu.I. Effect of Finite Deformations on Heat-Mass Transfer in Elastomer Ablating Materials // Int. Journ. of Heat Mass Transfer. 1997. Vol. 40, no. 3. P. 699-709.
10. Dimitrienko Yu.I. Internal Heat-Mass-Transfer and Stresses in Thin-Walled Structures of Ablating Materials // Int. Journ. of Heat Mass Transfer. 1997. Vol. 40, no. 7. P. 1701-1711.
11. Димитриенко Ю.И., Минин В.В., Сыздыков Е.К. Численное моделирование процессов тепломассопереноса и кинетики напряжений в термодеструктирующих композитных оболочках // Вычислительные технологии. 2012. Т. 17, № 2. С. 44-60.
12. Димитриенко Ю.И., Минин В.В., Сыздыков Е.К. Моделирование внутреннего тепломассопереноса и термонапряжений в композитных оболочках при локальном нагреве // Математическое моделирование. 2011. Т. 23, № 9. С. 14-32.
13. Димитриенко Ю.И., Минин В.В., Сыздыков Е.К. Моделирование термомеханических процессов в композитных оболочках при локальном нагреве излучением // Механика композиционных материалов и конструкций. 2011. Т. 17, № 1. С. 71-91.
14. McGurn Matthew, DesJardin P., Dodd A. Numerical simulation of expansion and charring of carbon-epoxy laminates in fire environments // International Journal of Heat and Mass Transfer. 2012. V. 55. P. 272-281.
15. Mouritz A.P., Mathys Z. Post-fire mechanical properties of marine polymer composites // Composite Structures. 1999. V. 47. P. 643-653.
16. Mouritz A.P., Mathys Z. Post-fire mechanical properties of glass-reinforced polyester composites // Composites Science and Technology. 2001. V. 61. P. 475-490.
17. Gibson A.G., Wright P.N.H., Wu Y.S., Mouritz A.P., Mathys Z., Gardiner C.P. The integrity of polymer composites during and after fire // Journal of Composite Materials. 2004. V. 38. P. 1283-1307.
18. Kardomateas G.A., Liu L., Birman V., Holmes J.W., Simitses G.J. Thermal buckling of a heat-exposed, axially restrained composite column // Composites Part A (Applied Science and Manufacturing). 2006. V. 37. P. 972-980.
19. Henderson J.B., Wiebelt J.A., Tant M.R. A model for the thermal response of polymer composite materials with experimental verification // Journal of Composite Materials. 1985. V.
19. P. 579-595.
20. Henderson J.B., Wiecek T.E. A mathematical model to predict the thermal response of decomposing expanding polymer composites // Journal of Composite Materials. 1987. V. 21. P. 373-393.
21. Sullivan R.M., Salamon N.J. A finite element method for the thermochemical decomposition of polymeric materials-1. theory // International Journal of Engineering Science. 1992. V. 30. P. 431-441.
22. Sullivan R.M., Salamon N.J. A finite element method for the thermochemical decomposition of polymeric materials-2. carbon phenolic composites // International Journal of Engineering Science. 1992. V. 30. P. 939-951.
23. Sullivan R.M. A coupled solution method for predicting the thermostructural response of decomposing, expanding polymeric composites // Journal of Composite Materials. 1993. V. 27.
P. 408-434.
24. Pering G.A., Farrell P.V., Springer G.S. Degradation of tensile and shear properties of composites exposed to fire or high temperature // Journal of Composite Materials. 1980. V. 14.
P. 54-68.
25. McManus H.L.N., Springer G.S. High temperature thermomechanical behavior of carbon-phenolic and carbon-carbon composites. 1. Analysis // Journal of Composite Materials. 1992. V. 26. P. 206-229.
26. McManus H.L.N., Springer G.S. High temperature thermomechanical behavior of carbon-phenolic and carbon-carbon composites. 2. Results // Journal of Composite Materials. 1992. V. 26. P. 230-255.
27. Димитриенко Ю.И. Тензорное исчисление. М.: Высшая школа, 2001. 575 с.
28. Попов Б.Г. Расчет многослойных конструкций вариационно-матричными методами. М.: Изд-во МГТУ им. Н.Э. Баумана, 1993. 294 с.
29. Зенкевич О. Метод конечных элементов в технике. М.: Мир,1975. 542 с.
30. Алфутов Н.А., Зиновьев П.А., Попов Б.Г. Расчет многослойных пластин и оболочек из композиционных материалов. М.: Машиностроение, 1984. 264 с.
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Modeling of non-stationary internal heat-mass-transfer in heat-shielding structures by the three-dimensional finite-element method
# 10, October 2013
DOI: 10.7463/1013.0606069
Dimitrienko Yu.I., Koryakov M.N., Balakshin A.I.
Bauman Moscow State Technical University, 105005, Moscow, Russian Federation
In this article the authors propose a finite-element method for calculating non-stationary heat-mass-transfer 3D problems in heat-shielding structures made of composite materials with consideration for thermo-decomposition and filtration of gases in pinholes. Specialized software that implements the numerical finite-element method for solving the problem of internal heat-mass-transfer in 3D elements of structures made of composite thermo-destructive materials was developed. In this paper the authors also present results of testing the developed method and software that illustrate working capacity of this method and its applicability to solving investigation problems of non-stationary heating modes of heat-shielding structures.
Publications with keywords: composites, numerical modelling, heat-shield structures, high temperatures, thermal decomposition, gas generation, pore pressure, finite element method Publications with words: composites, numerical modelling, heat-shield structures, high temperatures, thermal decomposition, gas generation, pore pressure, finite element method
References
1. Dimitrienko Yu.I. Mekhanika kompozitsionnykh materialovpri vysokikh temperaturakh [Mechanics of composite materials at high temperatures]. Moscow, Mashinostroenie, 1997. 368 p.
2. Dimitrienko Yu.I. Thermal stresses and heat mass-transfer in ablating composite materials. Int. Journ. of Heat Mass Transfer, 1995, vol. 38, no. 1, pp. 139-146.
3. Dimitrienko Yu.I. Thermal Stresses in Ablative Composite Thin-Walled Structures under Intensive Heat Flows. Int. Journal of Engineering Science, 1997, vol. 35, no. 1, pp. 15-31.
4. Dimitrienko Yu.I. Heat- Mass-Transport and Thermal Stresses in Porous Charring Materials. Journal of Transport in Porous Media, 1997, vol. 27, no. 2, pp. 143-170.
5. Dimitrienko Yu.I. Modelling of Mechanical Properties of Composite Materials under High Temperatures. Part 1. Matrix and Fibres. Int. Journal of Applied Composite Materials,
1997, vol. 4, no. 4, pp. 219-237.
6. Dimitrienko Yu.I. Modelling of Mechanical Properties of Composite Materials under High Temperatures. Part 2. Unidirectional Composites. Int. Journal of Applied Composite Materials, 1997, vol. 4, no. 4, pp. 239-261.
7. Dimitrienko Yu.I. Thermomechanical Behaviour of Composite Materials and Structures under High Temperatures. Part 1. Materials Composites. Part A. Applied Science and Manufacturing, 1997, vol. 28A, pp. 453-461.
8. Dimitrienko Yu.I. Thermomechanical Behaviour of Composite Materials and Structures under High Temperatures. Part 2. Structures. Composites. Part A. Applied Science and Manufacturing, 1997, vol. 28A, pp. 463-471.
9. Dimitrienko Yu.I. Effect of Finite Deformations on Heat-Mass Transfer in Elastomer Ablating Materials. Int. Journ. of Heat Mass Transfer, 1997, vol. 40, no. 3, pp. 699-709.
10. Dimitrienko Yu.I. Internal Heat-Mass-Transfer and Stresses in Thin-Walled Structures of Ablating Materials. Int. Journ. of Heat Mass Transfer, 1997, vol. 40, no. 7, pp. 1701-1711.
11. Dimitrienko Yu.I., Minin V.V., Syzdykov E.K. Chislennoe modelirovanie protsessov teplomassoperenosa i kinetiki napryazheniy v termodestruktiruyushchikh kompozitnykh obolochkakh [Numerical modeling of heat-mass-transfer and stress kinetics in thermodestructing composite shells]. Vychislitel'nye tekhnologii [Computational Technologies], 2012, vol. 17,
no. 2, pp. 44-60.
12. Dimitrienko Yu.I., Minin V.V., Syzdykov E.K. Modelirovanie vnutrennego teplomassoperenosa i termonapryazheniy v kompozitnykh obolochkakh pri lokal'nom nagreve [Modelling of internal heat and mass transfer and stress-strain state in composite shells under local heating]. Matematicheskoe modelirovanie, 2011, vol. 23, no. 9, pp. 14-32.
13. Dimitrienko Yu.I., Minin V.V., Syzdykov E.K. Modelirovanie termomekhanicheskikh protsessov v kompozitnykh obolochkakh pri lokal'nom nagreve izlucheniem [Mathematical modelling of thermomechanical processes in composite shells under local radiant heating]. Mekhanika kompozitsionnykh materialov i konstruktsiy [Journal on Composite Mechanics and Design], 2011, vol. 17, no. 1, pp. 71-91.
14. McGurn Matthew, DesJardin P., Dodd A. Numerical simulation of expansion and charring of carbon-epoxy laminates in fire environments. International Journal of Heat and Mass Transfer, 2012, vol. 55, pp. 272-281.
15. Mouritz A.P., Mathys Z. Post-fire mechanical properties of marine polymer composites. Composite Structures, 1999, vol. 47, pp. 643-653.
16. Mouritz A.P., Mathys Z. Post-fire mechanical properties of glass-reinforced polyester composites. Composites Science and Technology, 2001, vol. 61, pp. 475-490.
17. Gibson A.G., Wright P.N.H., Wu Y.S., Mouritz A.P., Mathys Z., Gardiner C.P. The integrity of polymer composites during and after fire. Journal of Composite Materials, 2004, vol. 38, pp. 1283-1307.
18. Kardomateas G.A., Liu L., Birman V., Holmes J.W., Simitses G.J. Thermal buckling of a heat-exposed, axially restrained composite column. Composites Part A (Applied Science and Manufacturing), 2006, vol. З7, pp. 972-980.
19. Henderson J.B., Wiebelt J.A., Tant M.R. A model for the thermal response of polymer composite materials with experimental verification. Journal of Composite Materials, 198S, vol.
19. pp. S79-S9S.
20. Henderson J.B., Wiecek T.E. A mathematical model to predict the thermal response of decomposing expanding polymer composites. Journal of Composite Materials, 1987, vol. 21, pp. З7З-З9З.
21. Sullivan R.M., Salamon N.J. A finite element method for the thermochemical decomposition of polymeric materials-1. Theory. International Journal of Engineering Science, 1992, vol. З0, pp. 4З1-441.
22. Sullivan R.M., Salamon N.J. A finite element method for the thermochemical decomposition of polymeric materials-2. Carbon phenolic composites. International Journal of Engineering Science, 1992, vol. З0, pp. 9З9-951.
23. Sullivan R.M. A coupled solution method for predicting the thermostructural response of decomposing, expanding polymeric composites. Journal of Composite Materials, 199З, vol. 27, pp. 408-4З4.
24. Pering G.A., Farrell P.V., Springer G.S. Degradation of tensile and shear properties of composites exposed to fire or high temperature. Journal of Composite Materials, 1980, vol. 14, pp. S4-68.
25. McManus H.L.N., Springer G.S. High temperature thermomechanical behavior of carbon-phenolic and carbon-carbon composites, 1. Analysis. Journal of Composite Materials, 1992, vol. 26, pp. 206-229.
26. McManus H.L.N., Springer G.S. High temperature thermomechanical behavior of carbon-phenolic and carbon-carbon composites, 2. Results. Journal of Composite Materials, 1992, vol. 26, pp. 2З0-255.
27. Dimitrienko Yu.I. Tenzornoe ischislenie [Tensor calculus]. Moscow, Vysshaya shkola, 2001. 575 p.
28. Popov B.G. Raschet mnogosloynykh konstruktsiy variatsionno-matrichnymi metodami [Calculation of multilayer structures by variational matrix methods]. Moscow, Bauman MSTU Publ., 199З. 294 p.
29. Zienkiewicz O.C. The Finite Element Method in Engineering Science. McGraw-Hill Publishing Co., London, 1971. 521 p. (Russ. ed.: Zienkiewicz O.C. Metodkonechnykh elementov v tekhnike. Moscow, Mir, 1975. 542 p.).
30. Alfutov N.A., Zinov'ev P.A., Popov B.G. Raschet mnogosloynykhplastin i obolochek iz kompozitsionnykh materialov [Calculation of multilayered plates and shells made of composite materials]. Moscow, Mashinostroenie, 1984. 264 p.
