DEVELOPMENT OF A THERMAL MATHEMATICAL MODEL FOR THE DIAGNOSIS OF POLYMER COMPOSITE MATERIALS USED IN AIRCRAFT MANUFACTURING AT NON-STATIONARY TEMPERATURES

Full Text

ВВЕДЕНИЕ В авиации полимерные композиционные материалы (ПКМ) очень часто эксплуатируются в условиях нестационарных температур. Поэтому возникает вопрос о проверки соответствия материалов условиям эксплуатации. Основные преимущества углеродных ПКМ состоят в высокой теплостойкости, малой плотности, стойкости к тепловому удару и облучению. Они длительное время могут работать при температурах до 773 °К в окислительной среде и 3273 °К - в инертной среде и вакууме. При этом их прочность с ростом температуры повышается в 1,5 - 2 раза. ОПИСАНИЕ ЗАДАЧИ Для проведения всесторонней диагностики свойств ПКМ, требуется рассмотреть тепловую модель, описывающую поведение ПКМ в условиях нестационарного нагрева. По своим структурам ПКМ большинство материалов состоят из некоторого числа слоев, поэтому интерес представляет поведение температуры как внутри слоя, так и внутри самого материала. Основными параметрами образца ПКМ являются e - диэлектрическая проницаемость и tgd - тангенс угла потерь. ПУТИ РЕШЕНИЯ ЗАДАЧИ Исследуемый материал подвергается одностороннему нестационарному нагреву, вследствие которого образуется неравномерное распределение температуры нагрева диагностируемого образца по толщине. Запишем уравнение теплопроводности с целью получения выражения описывающего распределение температуры нагрева по толщине. Для упрощения описательного процесса рассмотрим один слой, который будет являться базовым элементом для многослойной модели. Слои между собой будут сшиваться через граничные условия температур на их поверхности. Границы раздела слоем могут быть реальными или условными. Будем считать, что температура в каждой точке образца ПКМ в момент времени описывается функцией . В начальный момент времени все части тела имеют одинаковую температуру (температуру тела в нормальных условиях). Для решения поставленной задачи составим теплофизическое уравнение для образца ПКМ, считая, что образец-тело однородное по толщине и не имеет на поверхности изломов, трещин, ребер и т.д. В случае наличие различных неоднородностей температурный рельеф будет отличен в местах неоднородностей. Уравнение распространения тепла в произвольном объеме V, ограниченный гладкой поверхностью S, можно представить в виде нескольких слагаемых [1, 2]: 1) Изменение количества тепла в объеме V за промежуток времени , (1) где k - коэффициент внутренней теплопровод-ности, n - внутренняя нормаль к поверхности S, Т - температура нагрева тела. 2) На изменение температуры объема на DT за промежуток времени Dt надо потратить количество тепла , (2) где - теплоемкость вещества, - плотность вещества. 3) Функцией обозначим плотность тепловых источников находящихся в теле (количество выделяемого тепла в единицу времени в единицу объема). Количество поглощенной или выделяемой энергии за промежуток времени равно , (3) Полный баланс энергии для заданного объема будет представлен в виде: Q2=Q1+Q3. (4) Применив для второго интеграла теорему Остроградского, в целом получили уравнение теплопроводности для неоднородного изотропного тела. . (5) Будем считать, что тело однородное и параметры , и . (6) где ; . По теории теплопроводности для нахождения температуры внутри тела в любой момент времени необходимо знать начальное условие и граничное условие (тепловой режим на границе тела). Граничные условия могут задаваться разными способами в зависимости от того, что известно. Граничные условия на поверхности тела обычно задается тепловым потоком ; (7) или в каждой точке на поверхности тела задается температура - известная функция точки на поверхности и времени. Начальные условия это распределение температуры по толщине тела. . (8) Будем считать: нагрев производится СВЧ энергией, для связывания в единую систему электродинамические модели ПКМ и тепловую модель; тепловой поток производящий нагрев распространяется вдоль одной координаты (координаты z см. рис. 1), перпендикулярной поверхности тела, вглубь тела. Вследствие этого будем считать, что распространения тепла по двум другим координатам не происходит и температура нагрева тела имеет распределение только по координате z. В силу сделанных утверждений можно перейти от решения трехмерной задачи теплопроводности к решению одномерной. Для решения поставленной задачи потребуется уточнение характеристик условий нагрева и физического процесса нагрева. Для простоты описания будем считать: - излучение монохроматическое (источник нагрева -электромагнитная волна); - состояние тела равновесное, перенос тепла осуществляется за счет теплопроводности; - теплоотдачей всех видов пренебрегаем, в силу быстротечности процесса нагрева; - для исследуемого материала берется модель неограниченной пластины (рис.1) [3, 4]. Решение уравнения теплопроводности будем производить в духе электродинамики. Нагревательная энергия, которая поглощается образцом ПКМ будет равна: , (9) где - напряженность электрического поля, производящего нагрев импедансного тела; - импеданс исследуемого тела, определяемый как ,(10) - коэффициент затухания электромагнитной волны в импедансном теле , (11) где Теперь составим одномерное уравнение теплопроводности для ПКМ образца, считая, что температура нагрева (избыточная температура) распространяется только вдоль одной координаты (координаты Z). , (12) Т - избыточная температура, ; - коэффициент температуропроводности, ; - коэффициент теплопроводности; - теплоемкость; - плотность вещества; - Время облучения (нагрева) материала; - координата распространения электромагнитной волны, ; - толщина ПКМ образца. Это типичное уравнение теплопроводности, в котором первый член в правой части отвечает за количество тепла прошедшее путем теплопередачи через границу перпендикулярную плоскости ХY. Второй член правой части отвечает за количество тепла поглощенное материалом за время облучения. Для конкретного решения уравнения теплопроводности (1) необходимо уточнить начальные и граничные условия. Начальное условие: задается начальная температура по всей поверхности материала, в данном случае или . (13) Граничные условия: без учета теплоотдачи нагреваемого тела применяются граничные условия второго рода , (14) Решением уравнения (13) в общем виде [1, 2] имеет вид (12) и будет искаться через фундаментальную функцию: , (15) где , (16) (17) Выражение (16) представляет собой распределение начальной температуры тела по его толщине. Численное решение выражения (16) приводится ниже. . (18) Выражение (18) показывает распределение начальной температуры тела по толщине непосредственно перед началом исследуемого нагрева. Распределение избыточной или нагреваемой температуры по толщине имеет вид (17) ниже приведено численное решение: . (19) Выражение (19) качественно показывает дискретный прирост температуры нагрева по толщине тела и по времени нагрева. Задавшись дискретами по толщине и по времени мы получаем временное распределение приращенной температуры по толщине образца ПКМ. Точность вычисления температурного распределения будет зависеть от точности выбранного численного метода, и от шага дискретов по толщине и времени. Полученные результаты справедливы для поглощающих материалов, для полупрозрачных материалов математическая модель будет представлена как поглощающей, рассеивающей и испускающей. По этому в математической модели для полупрозрачных материалов необходимо рассматривать процессы рассеивания, испускания и переноса энергии [5]. РАСЧЕТНЫЙ АЛГОРИТМ ТЕМПЕРАТУРНЫХ ЗАВИСИМОСТЕЙ И ПКМ Как отмечалось раннее, при одностороннем нагреве диэлектрического материала происходит распределение температуры по толщине исследуемого материала. Из-за повышения температуры тела происходит изменение и . Поэтому при измерении температурных зависимостей и необходимо учитывать их изменение. Изменение электродинамических параметров импедансного тела и от первоначальных (в нормальных условиях) будет представлено как: - относительную диэлектрическую проницаемость в виде линейной аппроксимации от температуры нагрева: ; - тангенс угла диэлектрических потерь в виде параболической аппроксимации: , где и - параметры ПКМ в нормальных условиях. Для теоретического исследования температурных зависимостей и необходимо решить систему уравнений, состоящую из уравнений электродинамики и уравнений теплофизики. Решение уравнения теплофизики дает распределение температуры нагрева образца ПКМ по толщине при заданных параметрах нагревательного излучения и длительности воздействия на тело, тем самым производя моделирование теплового режима. Решение уравнений электродинамики позволяет определить комплексную диэлектрическую проницаемость в условиях нагрева по комплексным коэффициентам отражения или прохождения. Алгоритм последовательность действий вычисления температурных зависимостей приведен ниже: Входные данные: 1. Измеренные значения: Задание: погрешности : интервалов: 2. Вычислить температуру 3. Задание аппроксимаций ; 4. Вычислить 5. Оценка плохо - вернуться на корректировку п.3 6. Оценка плохо - вернуться на корректировку п.2 7. Вывод Описание алгоритма вычисления. Во-первых, производиться измерение комплексной диэлектрической проницаемости образца ПКМ и его температуры в обычных условия. Следующим шагом производится задание параметров нагревательного излучения или элемента в времени воздействия, аппаратура нагрева. В процессе нагрева образца ПКМ, через заданные временные интервалам производятся измерения температуры нагрева образца ПКМ (поверхности и внутри тела на заданных глубинах), и комплексных коэффициентов отражения и прохождения измерительного излучения. Контрольные измерения температуры позволяют производить корректировку точности вычисления распределения температуры по толщине исследуемого тела. Для чистоты диагностики измерение температуры проводится бесконтактным методом (пирометром) Рис. 1. Модель неограниченной пластины

About the authors

German Vyacheslavovich Dmitrienko

Ulyanovsk State Technical University, Separate Structural Unit “Institute of Aviation Technologies and Management”

Email: dmitrienko.german@yandex.ru
Doctor of Engineering Sciences, Associate Professor, Professor at the Department of Aircraft Engineering Ulyanovsk

Ekaterina Nikolaevna Zguralskaya

Ulyanovsk State Technical University, Separate Structural Unit “Institute of Aviation Technologies and Management”

Senior Lecturer at the Department of Aircraft Engineering Ulyanovsk

Georgy Leonidovich Rivin

Ulyanovsk State Technical University, Separate Structural Unit “Institute of Aviation Technologies and Management”

Email: avia@ulstu.ru
Candidate of Engineering Sciences, Associate Professor at the Department of Aircraft Engineering Ulyanovsk

Alexander Aleksandrovich Fedorov

Ulyanovsk State Technical University, Separate Structural Unit “Institute of Aviation Technologies and Management”

Candidate of Technical Sciences, Associate Professor at the Department of Aircraft Engineering Ulyanovsk

References

Supplementary files