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

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

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

МЕТОД ПОДСТАНОВКИ 0 РЕШЕНИИ ДВУМЕРНОЙ ЗАДАЧИ НАГРЕВА ТЕЛ С ПОМОЩЬЮ ПОДВИЖНЫХ ИСТОЧНИКОВ ТЕПЛА

Ключевые слова: система, распределенные параметры, подвижное управление, метод подстановки, температурное поле, поверхность.

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

В опубликованных работах предложен метод подстановки для расчета распределенных управле-

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

Суть этих методов состоит в следующем.

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

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

Данная работа посвящена обобщению метода подстановки для расчета распределенных управлений [1], успешно применяемого для одномерных систем, на двумерные системы нагрева подвижными источниками тепла.

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

Состояние объекта определяется функцией —(х, 1), х = (хр х2) е В, 1 1 10, описывающей распределение температуры в области В. Функция —(х, 1) удовлетворяет уравнению (вообще говоря, нелинейному) теплопроводности:

срд- - цд\ + Ф(0 = /(х, 1),

х = (Хр х2) е В, 1 1 0, (1.1)

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

-(х, 0) = -0(х), х е В. (1.2)

Здесь оператор Ь[ —], действующий по пространственной переменной х на функцию — = —(х, 1) определяется как

а с, р, X и а — теплотехнические параметры материала объекта, соответственно теплоемкость, плотность, коэффициент теплопроводности и коэффициент теплообмена поверхности и внешней среды, зависящие в общем случае от координат точки х и температуры; ф(—) — функция, описывающая теплоотвод с поверхности объекта во внешнюю среду;

д— — производная по направлению внешней

нормали к границе Г области В. Функция /(х, 1) описывает воздействие от подвижного источника, зависящее от параметров источника, таких как интенсивность, форма, т. е. распределение его интенсивности в пространстве и времени, и некоторых других параметров, например, дисперсии.

При построении аппроксимирующей функции будем рассматривать установившийся режим, мно-гоцикловый режим движения источника, при котором состояние объекта —(х, 1) по прошествии достаточно большого числа периодов становится функцией, близкой к периодической по времени при любом фиксированном х е В с тем же периодом Т, что и период движения источника. В работах [3, 4] показано, что в установившемся режиме

д—(х. 0 =о,и уравнение объекта (1.1) упрощается

и принимает вид

-Д—] + ф(—) = /(х, 1), х = (х1, х2) е В,

Обозначим желаемое состояние объекта — *(х), а аппроксимирующую функцию —А (х). Подставив ее в уравнение (1.1), получим требуемое распределенное воздействие:

К аппроксимирующей функции —А (х) предъявляются следующие требования:

1) удовлетворять граничным и начальным условиям (1.2);

2) быть дважды кусочно-дифференцируемой;

3) аппроксимировать желаемое состояние — *(х) в области Вх с требуемой точностью;

4) функция (1.5), полученная из уравнения (1.1), должна быть кусочно-непрерывной и удовлетворять в области В условиям: /*(х) 1 0, |/*(х)с1х < Цтяу

(Цтяу — максимальная интенсивность (мощность) источника) и быть равной нулю вне области В.

Требования 1 и 2 необходимы для того, чтобы

функция —А (х) удовлетворяла уравнению (1.4) с граничными и начальными условиями (1.2). Это и делает возможным расчет воздействия /*(х) по формуле (1.5).

Основную трудность в методе подстановки составляет построение аппроксимирующей функции

—А (х), удовлетворяющей требованиям 1—4. Одна

из особенностей состоит в том, что функция —А (х) должна быть построена в области В, в то время как желаемое состояние — *(х) определено в области В\ с В. Поэтому функция —А (х) строится отдельно в области Вх и в области В\ВХ с последующим сопряжением этих функций.

Другая особенность состоит в том, что способ построения функции —А (х) в области В существенно зависит от способа задания состояния — *(х), х е Вг В частности, оно может быть задано в явном виде (например, аналитически) или в виде сеточной функции, определенной в области В1 со значениями в каждом узле и т. п. Практически состояние —*(х) часто удобно задавать в виде кусочно-постоянной функции. В каждом случае требуется соответствующий способ построения аппроксимирующей функции. Эти способы обсудим далее.

2. АЛГОРИТМ РАСЧЕТА АППРОКСИМИРУЮЩЕЙ ФУНКЦИИ

Так как состояние — *(х) задано в области Вр а функция —А (х) определена в области В(ВХ с В), то

Рис. 1. Схема разбиения области Ь

построить функцию —А (х) означает аппроксимировать состояние — *(х) в области В1 и продолжить аппроксимированную функцию на область В, выполняя при этом требования 1—4.

Конфигурация области В1 может быть различной и зависит от условий задачи, от вида траектории (кольцевой, спиралевидной, пилообразной, построчной) перемещения источника теплового воздействия. Область В удобно разбить на три не-пересекающиеся области (рис. 1). На рисунке обозначено: дВ1 — граница области Вр Г^ — линия, отстоящая на заданное расстояние ^ от границы дВр О2 — область, ограниченная линиями

дВ1 и Г^; О3 = ВДО2; О1 = В\ВГ Функция —А (х)

строится отдельно для областей О1 и О3, а область О2 используется для «сшивки» функций, построенных в областях О1 и О3.

Тогда аппроксимирующую функцию можно

представить в виде —А (х) = £ —г(х)хг(х), где фун-

кция —.(х) определена в области О. %г(х) =

Функция —1(х), х е О1, определяется как решение краевой задачи (1.1), (1.2) в области О1, с граничным условием вида —1(х) х е8П^ = И» х едВ, где функция ц(х), х е дВ1, может быть выбрана достаточно произвольно, но так, чтобы не нарушались перечисленные в § 1 условия 1—4.

Способ построения функции —3(х), х е О3, зависит от способа задания состояния — *(х), х е В1.

Если функция — *(х) задана в явном виде (например, аналитически) и удовлетворяет требованиям 2 и 3, то функция —3(х) совпадает с ней.

Если функция — *(х) задана приближенно, например, в виде сеточной функции в области В1 со значениями в каждом узле, то для построения функции —(х) во всей области применяются методы аппроксимации (сплайнами, методом наименьших квадратов, тригонометрическими полиномами, поверхностями Кунса [5, 6] и др.).

Часто состояние — *(х) задается в виде кусочно-постоянной функции следующим образом.

Рис. 2. Схема расчета сопряжений кусочно-постоянной функции:

а — пример расположения прямоугольных областей; б — график перехода температуры от начального значения к требуемому

В области В1 выделяются прямоугольные области £от, т = 1, 2, . М (рис. 2, а). Местонахождение, взаимное расположение, а также число этих областей может быть любым, что определяется условиями задачи: в каком месте объекта и до какого значения нужно довести температуру поверхности. Эти области могут как соприкасаться, так и быть достаточно разрежены. В каждой из этих областей функция — *(х) постоянна и равна Ст. В других точках области В1 функция — *(х) равна С = т1Пт=1,2. ^.

Тогда функция —А (х) строится следующим образом. Назначаются некоторые окрестности каждой области (например, шириной ^от), т. е. определяется область «сшивки» О2. Далее производится сглаживание функции — *(х) в — окрестностях каждой области так, чтобы обеспечить

непрерывный переход функции —А (х) от значения С к значению Ст при выполнении требований 2 и 3 (рис. 2, б).

Аппроксимирующую функцию представим в виде:

|Cm(x1> x2) пРи x Є

X при X Є т = 1, 2, . М.

Здесь через ¿т обозначены ^-окрестности областей а си(х19 Х2) = ст1 > (хх) ст2> (Х2). Функции С.> (х), ] = 1, 2 определяются по формулам

📎📎📎📎📎📎📎📎📎📎

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

Похожие темы научных работ по математике , автор научной работы — Финягина Валерия Ивановна

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