Тепловой процесс при термообработке гуммированных изделий выражает изменение состояния тепловой системы, которое осуществляется теплопроводностью, конвекцией или излучением. Его можно описать детерминированными или стохастическими соотношениями. Классическая теория теплопередачи основана на детерминированном выражении закономерностей теплопередачи посредством основного уравнения теплопроводности с соответствующими граничными и начальными условиями, при этом ни одна из переменных и ни один из физических параметров в уравнениях не носят случайного характера. Детерминированная математическая модель в абсолютном большинстве случаев является исходной моделью для тепловых задач, возникающих при гуммировании [3, 5].
При нестационарном процессе, неоднородных граничных условиях, а также при граничных условиях IV рода невозможно получить решение задачи аналитическими методами [6]. Применение же численных методов значительно расширяет классы решаемых задач теплообмена при изготовлении гуммированных объектов. При дискретизации на первом этапе выбираются узловые точки в пространственной и временной областях. На втором этапе составляется система алгебраических уравнений относительно значений искомых температур в этих узловых точках. На третьем - проводится решение системы, и находятся значения температур. Рассмотрим гуммированный объект произвольной формы в двумерной системе координат с площадью S и поверхностью L [6]:
(1)
с нулевыми начальными условиями:
Т(х, у, 0) = 0, (2)
при граничных условиях третьего рода:
(3)
где λ - теплопроводность; T - температура; α - коэффициент теплоотдачи на границе; qυ, qs - объемная и поверхностная плотности мощности источников теплоты. Метод конечных элементов (МКЭ) основан на определении температурного поля путем приближенного решения соответствующей вариационной задачи. Вариационная формулировка рассматриваемой краевой дифференциальной задачи (1)-(3) состоит в следующем. Задача решения уравнения (1) с граничными условиями (3) эквивалентна задаче определения функции Т(х, y), минимизирующей функционал I[T(х, у)] вида
(4)
В МКЭ приближение для искомой функции Т(х,у) отыскивается в виде
(5)
где ат - неизвестные постоянные коэффициенты, a fm(x,y) - известные функции пространственных координат. Разобьем гуммированный объект на N треугольных элементов и введем М узлов во всех вершинах треугольников. Присвоим сквозную нумерацию всем элементам (n = 1, ..., N) и всем узлам (m = 1, ..., М). Искомыми величинами в МКЭ являются приближенные значения температуры tm в узлах m = 1, ..., М. Распределение температуры в каждом элементе запишется в виде суммы, в которую входят три функции формы элемента, умноженные на приближенные значения температуры в его трех узловых точках. Распределение температуры в n-м треугольном элементе t(n)(x,y) имеет вид:
(6)
где ,
,
- линейные функции координат ?,yх,y, равные единице в узлах i, j или k соответственно и равные нулю в двух других узлах. Таким образом, для функции формы
должны выполняться равенства
(7)
Аналогичные соотношения имеют место для и
. Используя условия (7), можно выразить коэффициенты
,
,
(m = i, j, k) через координаты узлов i, j, k. Здесь и далее индекс элемента (n) будем опускать. Из условий (7) следует:
;
;
, (8)
где S - площадь треугольника, значение которой вычисляется через координаты узлов:
(9)
Градиент температуры в каждом элементе имеет постоянное значение, и производные по координатам определяются соотношениями:
(10)
Система уравнений для определения температур в узлах составляется на основании условий минимума функционала (4). Этот функционал можно представить в виде суммы интегралов по всем элементам:
(11)
(12)
Условия минимума функционала можно с учетом (11) переписать в виде
m = 1, ..., M. (13)
Вычислим сначала интеграл по площади элемента:
(14)
Если стороны элемента принадлежат границе области, то в функционале следует учесть интеграл по этим сторонам. Пусть граничной является сторона Lij. Тогда
(15)
Здесь учтено, что Fk(x,y) = 0 на стороне Lij. Если граничной является сторона Lik, то выражение для интеграла запишется аналогично, но вместо Lij подставляем Lik, а вместо tj - tk. Если же к границе прилегает сторона Ljk, то рассматриваемый интеграл равен нулю, так как функция формы Fi(х, у) для узла i равна нулю на стороне Ljk и, следовательно, распределение температуры на этой стороне не зависит от ti. Таким образом, окончательно получаем для производной от функционала I(n) по температуре ti следующее выражение:
(16)
Аналогично можно получить выражения для производных от I(n) по температурам tj и tk. Проанализируем теперь с учетом (16) структуру системы (13). Видно, что производная представляет собой сумму произведений неизвестных температур ti,tj,tk на постоянные известные коэффициенты, зависящие от координат узлов и параметров задачи, а также постоянных известных членов, не зависящих от искомых температур. Приравнивая левые части уравнений (13) к нулю, получаем линейную систему разностных уравнений относительно неизвестных температур узловых точек. Запишем полученные выше соотношения для частных производных функционала в матричном виде GT = Ф, где G - глобальная матрица теплопроводности размером М×М; Т - вектор-столбец искомых значений температур в М узлах; Ф - глобальный вектор-столбец тепловых потоков.
(17)
Особенностью объектов с распределенными параметрами (ОРП) является то, что они имеют пространственную протяженность и их состояние характеризуется одной или несколькими величинами, зависящими не только от времени, но и от точки области пространства, в которой формируются свойства ОРП [1]. Часто требуется, чтобы функционал
, γ ≥ 1, (18)
где T*(τ) - заданная программа изменения температуры T(l, τ), достигал своего минимального значения [5]. Задача заключается в том, чтобы создать такую систему управления заданием регулятора температуры, в которой отклонение средней температуры гуммированного изделия, выходящего из оборудования,
0 ≤ y ≤ l, 0 ≤ τ ≤ T
при у = l было наименьшим в смысле (1). При этом гарантируется достижение заданных свойств в гуммированном изделии. Тем не менее необходимо создание автоматизированной системы управления технологическим процессом на базе адаптивных математических моделей [3]. Процесс изменения температурных полей в вулканизуемом изделии описывается уравнением нестационарной теплопроводности типа Фурье, а оценка степени вулканизации и качества вулканизуемого изделия интегральным нелинейным уравнением Аррениусовского типа [4, 5]. Следовательно, данный процесс как объект управления принадлежит к ОРП, описываемым нелинейными дифференциальными уравнениями в частных производных параболического типа. Градиентные методы решения задач оптимального управления ОРП, описываемые параболическими уравнениями, рассмотрены в [4, 5]. В соответствии с этими методами требуется, регулируя температуру внешней среды изделия управлением u(t), сделать распределение температуры в изделии Т(х,τ,и) к заданному моменту времени τ, равному заданному распределению температуры Тз(х). Температура внешней среды для однородного нагреваемого стержня 0 < x < l
(19)
где a ≤ u(t) ≤ b - числа, выражающие крайние допустимые значения температуры; un+1(t) - управление в n + 1 приближении; - вспомогательное приближение, определяемое как
(20)
где T, α, ν, l - положительные постоянные коэффициенты, a φ(t) - решение краевой задачи.
(21)
(22)
(23)
(24)
Из (3) для множества
a ≤ u(t) ≤ b 0 ≤ t ≤ τ
следует:
(25)
В целях повышения помехоустойчивости алгоритма формулу (25) запишем в виде
(26)
где Δ - заданная постоянная величина. Коэффициент αn > 0, входящий в формулу (19), находится из условия
min qn(α) = qn(αn),
где
(27)
(28)
Откуда видно, что qп(α) достигает экстремума в точке
(29)
это значит, что квадратный трехчлен (28) достигает своего минимума на отрезке 0 ≤ α ≤ 1 в точке при
и в точке
при
, то есть можно записать
Полученное αn подставляется в (19), откуда следует (n + l)-е приближение для решения задачи. Если αn = 0 или
,
квадратный трехчлен (28) вырождается, итерационный процесс прекращается, и управление un(t) = u* будет неполным оптимальным решением задачи. Суммарные потери F(t0, tk) от начального момента t0 до конечного tk определяются как
(30)
где f(τ) - текущие потери производства в момент τ[3, 5]. Если f(τ) можно связать с оперативной информацией об объекте
состоящей из полученных к моменту t измерений текущих и прошлых значений сигналов ω, u, и v, то задача сводится к построению алгоритма, преобразующего J(τ) в сигнал оптимального управления uоп(τ), доставляющего минимум функционалу (30), то есть
(31)
С учетом цифровой реализации критерия оптимизации , где Т - фиксированный период. При этом класс допустимых сигналов управления ограничен семейством ступенчатых функций
(32)
где un - варьируемые переменные. Оптимальный сигнал un вычисляется как
(33)
где Sn(un,1) - прогноз будущих потерь в промежутке tn+l, tk. При данной J(tn) и данном un, при условии, что все будущие значения сигнала (12) будут строиться в моменты времени tn+1 = (n + i)T (i = 1, 2, ..., n - k) аналогично (13), то есть на основании полученной в соответствующие моменты времени информации J(tn+i) будут реализованы оптимальные значения . Если подынтегральная функция в (30) задана в виде
(34)
где f1, f2 - выпуклые, неотрицательные определенные функции от u и v, то оценить суммарные потери непосредственно в виде функции от условного прогноза можно формулой
(35)
К достоинствам описанных методов следует отнести возможность их реализации на базе средств современной вычислительной техники. Можно сделать вывод о том, что наиболее подходящими для программной интерпретации на ЭВМ являются модели, реализуемые при помощи численных методов, среди которых можно выделить МКЭ. Наиболее перспективными путями повышения качества изделий и сокращения продолжительности процесса являются автоматическая корректировка режимов вулканизации и автоматическое изменение параметров граничных условий. Задачу синтеза адаптивных дискретных устройств управления режимами работы ОРП целесообразно решать на основе методов динамической оптимизации условного прогнозирования и математического моделирования процессов в реальном и ускоренном масштабах времени.
Рецензенты:
-
Аншелес В.Р., д.т.н., профессор, профессор Череповецкого государственного университета, г. Череповец.
-
Синицын Н.Н., д.т.н., профессор, зав. кафедрой «Промышленная теплоэнергетика» Череповецкого государственного университета, г. Череповец.
Работа поступила в редакцию 20.02.2012.
Библиографическая ссылка
Осипов С.Ю., Потапов Н.В., Скоробогатова А.Ю., Осипов Ю.Р. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ТЕПЛООБМЕНА КАК ОСНОВА АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ КИНЕТИКОЙ НЕИЗОТЕРМИЧЕСКОЙ ВУЛКАНИЗАЦИИ ЭЛАСТОМЕРНЫХ ПОКРЫТИЙ // Фундаментальные исследования. – 2012. – № 3-3. – С. 647-651;URL: https://fundamental-research.ru/ru/article/view?id=29775 (дата обращения: 17.08.2024).