Исследование аналитической модели температурного поля диска газотурбинного двигателя

Размещено на http://www.allbest.ru/

Самарский государственный технический университет

Исследование аналитической модели температурного поля диска газотурбинного двигателя

Иван Александрович Данилушкин - к.т.н., докторант.

Антон Андреевич Московцев - аспирант.

E-mail: idanilushkin@mail.ru

443100, г. Самара, ул. Молодогвардейская, 244

Аннотация

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

Ключевые слова: диск газотурбинного двигателя, термопластическое упрочнение, математическая модель температурного поля диска, объект с распределенными параметрами, функция Грина.

Работа выполнена при финансовой поддержке грантов РФФИ №09-08-00297-а, №10_08_00754_а; ФЦП «Научные и научно-педагогические кадры инновационной России на 2009-2013 годы», заявка НК_66П/11, заявка 2010-1.3.1-230-009/8; АВЦП «Развитие научного потенциала высшей школы», проект №2.1.2/4236.

Усталостное разрушение деталей газотурбинных двигателей, таких как валы, лопатки, диски турбин, является основной причиной выхода их из строя. В конструкции турбины соединение рабочих лопаток с пазом диска является крайне напряженным и ответственным местом. Соединение передает на диск все нагрузки, которые под действием центробежных сил влияют на лопатку. Долговечность паза диска во многом определяется качеством поверхностного слоя, для упрочнения которого могут быть применены различные методы. Одним из таких методов является термопластическое упрочнение, включающее две стадии: нагрев до заданной температуры и последующее спрейерное охлаждение. При этом нагрев должен осуществляться с заданной точностью, определяемой технологией процесса термопластического упрочнения. Учет ряда конструктивных ограничений, с одной стороны, и требований к снижению энергопотребления, с другой, привели к созданию установки термопластического упрочнения елочного паза дисков турбоагрегатов [1, 2]. Применение индукционного нагрева для нагрева края диска позволяет, варьируя форму индуктора, получить желаемое температурное распределение на нагреваемом участке. Однако для решения задач выбора оптимальной формы индуктора и оптимального управления процессом термопластического упрочнения необходимо разработать модель температурного поля, учитывающую пространственную распределенность теплоисточников по поверхности нагреваемого диска.

Полученная в работе [3] математическая модель температурного поля диска в полярных координатах не учитывает потери с торцевых поверхностей диска. Эти потери могут быть учтены в модели двумя способами: 1) путем охвата передаточной функции, полученной в [3], распределенной обратной связью по температуре и нахождением передаточной функции замкнутой системы; 2) за счет добавления слагаемого, описывающего потери с торцевых поверхностей, в исходную постановку задачи моделирования. Ниже рассмотрен второй способ.

Температурное поле диска в полярных координатах описывается уравнением теплопроводности:

,

, , ,(1)

где - температурное поле диска, - температуропроводность, , c - теплопроводность и удельная теплоемкость материала диска, - плотность материала, R - внешний радиус диска, - функция, описывающая распределение удельной мощности теплоисточников. Температурным распределением по толщине диска можно пренебречь в силу его конструктивных параметров. Функция , учитывающая потери тепла за счет конвективного теплообмена между торцевыми поверхностями диска и окружающей средой, принимает вид [4]

,(2)

где - коэффициент теплообмена между средой и диском, - температура среды, H - толщина диска газотурбинного двигателя. Коэффициент 2 появляется для учета обеих торцевых поверхностей диска.

Конструкция индуктора такова, что позволяет рассматривать задачу в симметричной постановке относительно полярной оси. Тогда краевые условия для уравнения (1) записываются следующим образом:

;(3)

;(4)

.(5)

Для границы пренебрегая потерями с боковой поверхности диска газотурбинного двигателя, записываем:

.(6)

Начальное распределение температуры по радиусу диска примем равным нулю: температурный теплоисточник перадаточный

.(7)

Примем также без потери общности . Тогда с учетом (2) уравнение (1) можно переписать в виде:

(8)

Модель (8), дополненная соответствующими краевыми (3)-(6) и начальным условием (7), описывает температурное поле в диске с учетом потерь с торцевых поверхностей диска. Методы структурной теории распределенных систем позволяют получить решение для неоднородного дифференциального уравнения с помощью операции пространственного интегрирования произведения функции распределенных теплоисточников и функции Грина уравнения по области определения [3], [4].

Функция Грина (импульсная переходная функция) [5] есть решение задачи (3) - (8) при воздействии импульсного источника в точке , , задаваемого уравнением вида

.(9)

Найдем функцию Грина объекта с помощью конечных интегральных преобразований. Чтобы избавиться от операции дифференцирования по переменной , используется конечное косинус-преобразование Фурье [6]. Изображение функции вычисляется по формуле

, .(10)

Переход от изображения функции к ее оригиналу осуществляется по формуле

.(11)

Применяя преобразование (10) к дифференциальному уравнению теплопроводности (8) с учетом граничных условий (3), (4), получим

,(12)

.(13)

Чтобы избавиться от операции дифференцирования по переменной , используем конечное интегральное преобразование Ханкеля [6]. Изображение функции находится по формуле

,(14)

где - корень характеристического уравнения

.(15)

Переход от изображения функции к ее оригиналу осуществляется по формуле

.(16)

Применяя к уравнению (12) преобразование Ханкеля (14), получаем

,(17)

.(18)

Применив к уравнению (17) преобразование Лапласа, получим

,(19)

,(20)

тогда выражение (19) принимает вид

.(21)

Выполнив обратное преобразование Лапласа, получаем

.(22)

Граничное условие второго рода (7) определяет вид характеристического уравнения для нахождения собственных чисел конечного интегрального преобразования Ханкеля (15). Собственные числа обращают в нуль интеграл (14) во всех случаях, кроме случая n = 0. Поэтому при n = 0 выражение для обратного преобразования Ханкеля (16) принимает вид

,(23)

т. е. появляется слагаемое, соответствующее , а при n > 0 выражение для обратного преобразования Ханкеля (16) принимает вид

,(24)

где - положительные корни трансцендентного уравнения (15), , .

Применив к (23), (24) обратное косинус-преобразование Фурье (11), получим функцию Грина объекта

. (25)

Изображение функции Грина по Лапласу есть передаточная функция объекта в терминах структурной теории распределенных систем [7].

Функция распределения источников тепла, создаваемых индуктором в зоне нагрева (см. рис. 1), описывается выражением

.(26)

Рис. 1. Область распределения источников тепла

Температурное поле диска для произвольного распределения теплоисточников описывается выражением [7]

.(27)

Подставив (25) и (26) в (27), получаем

.(28)

Выражение (28) было использовано для расчета температурного поля диска со следующими параметрами: м, м, рад, кг/м³, Дж/(кгК), Вт/(мК), Вт/(м²К), МВт/м³.

1 - рад; 2 - рад; 3 - рад; 4 - рад

1 - м; 2 - м; 3 - м; 4 - м; 5 - м

а) б)

1 - рад; 2 - рад; 3 - рад; 4 - рад

1 - м; 2 - м; 3 - м; 4 - м; 5 - м

в) г)

Рис. 2. Температурное поле диска

На рис. 2 представлены результаты расчета в разные моменты времени.

По формуле (28) были рассчитаны температурные поля диска в разные моменты времени. На рис. 2, а; 2, б приведены температурные распределения по радиальной и угловой координатам при c; на рис. 2, в; 2, г - при c. При расчете принято, что , . Из-за учета в расчетах конечного числа собственных функций возникает вычислительная ошибка (в некоторых точках температура диска становится ниже температуры среды), которая уменьшается с увеличением количества учитываемых собственных функций.

Полученное выражение для функции Грина (25), а также выражение (28) для температурного поля диска при непрерывно действующих источниках тепла будут использоваться в дальнейшем для решения задач оптимизации конструкции индуктора и синтеза алгоритма оптимального управления процессом нагрева при термопластическом упрочнении диска газотурбинного двигателя.

Библиографический список

1. Головачев А.Л., Данилушкин А.И., Мишанин Е.А. Система индукционного нагрева для термообработки елочного паза дисков турбоагрегатов // Вестник Саратовского государственного технического университета. - 2006. - №1(10). - C. 108-113.

2. Головачев А.Л. Разработка и исследование индукционных систем для ремонтно-восстановительных технологий роторов газотурбинных двигателей: Автореф. дис. … канд. техн. наук: 05.09.10 / Самара: Самар. гос. техн. ун-т, 2009. - 20 с.

3. Данилушкин И.А., Московцев А.А. Математическая модель температурного поля диска газотурбинного двигателя // Вестник Самарского государственного технического университета. Сер. Технические науки. - Вып. 2(24). - 2009. - С. 218-221.

4. Лыков А.В. Теория теплопроводности. - М.: Высшая школа, 1967. - 599 с.

5. Рапопорт Э.Я. Структурное моделирование объектов и систем управления с распределенными параметрами. - М.: Высшая школа, 2003. - 299 с.

6. Снеддон И. Преобразования Фурье. - М.: Издательство иностранной литературы, 1955. - 668 с.

7. Бутковский А.Г. Характеристики систем с распределенными параметрами. - М.: Наука, 1979. - 224 с.

Размещено на Allbest.ru