Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,074

РАЗРАБОТКА МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ ТЕПЛОПЕРЕДАЧИ В МНОГОСЛОЙНЫХ КОНСТРУКЦИЯХ

Тугов В.В. 1 Акимов И.А. 2
1 Оренбургский государственный университет
2 Оренбургский государственный педагогический университет
Данная работа посвящена разработке математических моделей и изучению теплопередачи, происходящей в многослойных цилиндрических конструкциях для полного цикла (предварительный разогрев, полимеризация и охлаждение) процесса производства композиционных материалов на промышленных установках. Работа продолжает исследования теплофизических процессов в многослойных конструкциях с фазовыми переходами. С целью повышения эффективности синтеза технологических процессов с фазовыми переходами предлагается использование метода изотермических поверхностей, как основы для построения математического и алгоритмического обеспечения управляющих систем. Для решения поставленных модельных задач предлагается использовать методы интегральных преобразований и метод конечных интегральных преобразований Ханкеля, которые дают возможность решать задачи, связанные с получением многослойных структур. Универсальность применяемого математического аппарата позволяет создавать программное обеспечение для проведения имитационного моделирования рассматриваемого процесса на компьютере, а также системы управления производством композиционных материалов на различных режимах работы.
теплопередача
полимеризация
многослойные конструкции
композиционные материалы
фазовый переход
1. Сафонов А.А. Математическое моделирование процесса полимеризации при изготовлении полимерных композиционных материалов / А.А. Сафонов, Б.А. Щеглов // Проблемы машиностроения и надежности машин. – 2011. – № 6. – С. 55–60.
2. Трофимов П.А. Изготовление композиционных материалов для летательных аппаратов / П.А. Трофимов // Вестник ОГПУ. – 2011. – № 3(59). – С. 17–20.
3. Тугов В.В. Исследование влияния погрешностей технологий на эксплуатационные характеристики авиаконструкций из композиционных материалов / В.В. Тугов, И.А. Акимов // Научно-технический вестник Поволжья. – 2016. – № 1. – С. 58–61.
4. Гремячкин В.М. Тепло- и массоперенос при гетерогенных процессах. Учебное пособие / В.М. Гремячкин. – М.: Изд-во МГТУ им. Н.Э. Баумана, 2010. – 175 c.
5. Карташов Э.М. Аналитические методы в теории теплопроводности твердых тел / Э.М. Карташов. – М.: Высшая школа, 1985. – 479 c.
6. Акимов И.А. Исследование массопереноса на втором этапе производства композиционных материалов методом полимеризации в установках автоматического ведения технологического процесса (АВТП) / И.А. Акимов, А.И. Акимов, Т.И. Жумагазеев // Научно-технический вестник Поволжья. – 2015. – № 2. – С. 63–67.
7. Шафеев М.Н. О сходимости метода изотермических поверхностей / М.Н. Шафеев // ВИНИТИ № 6597-В86. – Уфа, 1986. – 13 с.
8. Козлов В.Н. Зависимость механических свойств композиционных материалов от температурного режима полимеризации / В.Н. Козлов, А.И. Акимов, М.А. Фатыхов // Известия высших учебных заведений. Северо-Кавказский регион. Технические науки, 2009. – C. 83–89.

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

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

Полный цикл производства КМ включает в себя три этапа: I этап – режим предварительного разогрева; II этап – режим полимеризации; III этап – режим охлаждения. Для исследования технологических процессов необходимы сведения о теплопередаче на всех этапах производства композиционных материалов.

Цель работы – разработать и исследовать математические модели процессов теплопередачи в многослойных цилиндрических изделиях в установках автоматического ведения технологического процесса (АВТП) на всех этапах производства КМ.

Постановка и решение задачи на первом этапе

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

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

Исследуем изделие в форме полого многослойного цилиндра, внутри которого размещается технологический мешок, уложенного в пресс-форму. К пресс-форме подключается электрообогрев. При этом температура может передаваться только в радиальном направлении, так как по длине изделия z и по φ, то есть по цилиндрическим координатам, температура одинакова. Параметр теплопроводности λi tug02.wmf для каждого слоя одинаков, так как слои состоят из одного и того же КМ, но под воздействием температуры могут незначительно изменяться. Величина теплового потока находится в зависимости от глубины проникновения и от начального нагрева.

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

Исследуемое изделие можно рассматривать как многослойный цилиндр бесконечной длины, так как его длина значительно больше диаметра, то есть tug03.wmf (рис. 1).

При сделанных выше предположениях постановку задачи на первом этапе математически можно описать так [4]:

tug04a.wmf

tug04b.wmf, tug05.wmf,

tug06.wmf, tug07.wmf (1)

начальное условие

tug08.wmf, (2)

где Uk – температура в k-ом слое; r – радиус цилиндра; ak – температуропроводность k-ого слоя изделия; wk – функция источников.

tugov1.wmf

Схема расположения слоев изделия

Запишем в общем виде граничные условия на свободных поверхностях (подбирая соответствующие коэффициенты, можно получить граничные условия I, II или III рода с переменным потенциалом среды и произвольно меняющимся потоком тепла):

при r = R0

tug09.wmf; (3)

при tug10.wmf

tug11.wmf, (4)

где U0 и tug12.wmf – соответственно температуры на внешней и внутренней поверхности изделия;

g1, g2 – потоки тепла с внешней и внутренней поверхностей изделия;

α1, α2, β1, β2, χ1, χ2 – коэффициенты теплообмена между средами.

Кроме того, наложим условия сопряжения между отдельными слоями:

tug13.wmf;

tug14.wmf,

tug15.wmf.

Методом конечных интегральных преобразований Ханкеля [5] найдем решение исходной задачи для каждого слоя в виде

tug16.wmf (5)

Постановка и решение задачи на втором этапе

Процесс полимеризации происходит на втором этапе производства КМ. При этом исследуемое изделие считаем бесконечно полым цилиндром, так как его длина больше диаметра tug17.wmf. Цилиндрические координаты φ и z исключаем из рассмотрения, так как изделие находится внутри пресс-формы и тепло распространяется только в радиальном направлении по изотермическим поверхностям, а по длине и по φ теплопередача будет одна и та же.

Применяя такие предположения, постановка задачи по радиальной схеме в цилиндрической системе координат имеет следующий вид [6]:

tug18.wmf tug19.wmf, tug20.wmf, (6)

tug21.wmf,

tug22.wmf для зоны полимеризации – I, tug23.wmf для противоположной зоны – II, при начальных условиях

tug24.wmf (7)

и при граничных условиях

tug25.wmf; tug26.wmf; (8)

tug27.wmf; tug28.wmf;

tug29.wmf; tug30.wmf,

при tug31.wmf;

tug32.wmf,

где ak, λk, αk – коэффициенты температуропроводности, теплопроводности и теплопередачи; σ – теплота полимеризации; a – безразмерный параметр теплового влияния в k-ом слое; tug33.wmf, tug34.wmf – достаточно гладкие функции своих аргументов, которые подлежат определению через условия (1) и (2); tug35.wmf – температура k-го слоя изделия.

Используем метод изотермических поверхностей [7] для решения поставленной задачи, при этом заменяем истинное распределение температуры tug36.wmf внутри каждой области tug37.wmf и tug38.wmfпри фиксированных положениях границы tug39.wmf нестационарным распределением температуры tug40.wmf, tug41.wmf.

Итоговое решение задачи имеет вид

tug42.wmf (9)

где последовательность {τj} определяется рекуррентной формулой

tug43.wmf, (10)

в которой tug44.wmf;

tug45.wmf, (11)

где

tug46.wmf, (12)

tug47.wmf, tug48.wmf, tug49.wmf.

Постановка и решение задачи на третьем этапе

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

Исследованиями установлено, что в производстве лонжеронов лопастей вертолетов в установках АВТП, скорость остывания готового изделия не должна превышать 2–3 градуса в минуту [8]. Нарушение данного процесса остывания приводит к опасным остаточным термонапряжениям. Поэтому необходимо составить математическую модель естественного остывания изделия в установках АВТП в зависимости от температуры окружающей среды, так как АВТП просто отключается от сети после завершения второго этапа производства КМ. Далее изделие охлаждается от температуры полимеризации (160–165 °С) до температуры окружающей среды (цеха).

Лонжероны лопастей вертолетов имеют форму полого цилиндра. Длина изделия значительно длиннее диаметра tug50.wmf, поэтому изделие можно рассматривать как полый бесконечный цилиндр с внутренним радиусом R0 и с внешним радиусом R.

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

tug51.wmf,

tug52.wmf, tug53.wmf, (13)

с начальным условием

tug54.wmf, (14)

и граничными условиями

tug55.wmf, (15)

tug56.wmf,

где R – радиус внешней поверхности; R0 – радиус внутренней поверхности изделия; tug57.wmf – температура изделия; a – температуропроводность материала; r – переменная в радиальном направлении; τ – время.

Так как tug58.wmf и tug59.wmf, то ищем решение задачи в виде

tug60.wmf. (16)

При этом функция v(r) должна удовлетворять дифференциальному уравнению

tug61.wmf (17)

и граничным условиям

tug62.wmf, tug63.wmf.

Функция tug64.wmf удовлетворяет дифференциальному уравнению

tug65.wmf,

tug66.wmf, tug67.wmf (18)

при начальном условии

tug68.wmf

и граничных условиях

tug69.wmf, tug70.wmf.

Применим классический метод Фурье и окончательно получим распределение температуры в изделии в зависимости от r и τ в виде

tug71.wmf (19)

Собственные числа kn определяются из уравнения

tug72.wmf.

Заключение

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


Библиографическая ссылка

Тугов В.В., Акимов И.А. РАЗРАБОТКА МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ ТЕПЛОПЕРЕДАЧИ В МНОГОСЛОЙНЫХ КОНСТРУКЦИЯХ // Фундаментальные исследования. – 2017. – № 8-2. – С. 320-324;
URL: http://www.fundamental-research.ru/ru/article/view?id=41668 (дата обращения: 21.09.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074