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

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

Осипов Ю.Р. 1 Демакова А.О. 1 Шилова А.О. 1 Ночевалова Ю.А. 1 Осипов С.Ю. 2 Богданов Д.А. 1
1 ФГБОУ ВО «Вологодский государственный университет»
2 ФГБОУ ВО «Тверской государственный технический университет»
Рассмотрены исследования температурных режимов процесса восстановления гуммировочных покрытий с целью повышения их качества. Процесс протекает в условиях действующих возмущений стохастической природы, которые являются факторами неопределенности и влияют на ход вулканизации. Разработана математическая модель тепловых режимов вулканизации. Стохастическая задача решена с применением вероятностных критериев. Построены кинетические кривые теплового режима прогрева покрытия, полученные экспериментальным путем и при решении вероятностной задачи по предложенной математической модели.
. теплообмен
температура
гуммировочное покрытие
стохастическая задача
неопределенность
вероятность
случайная величина
управление
1. Боднер В.А. Оптимизация терминальных стохастических систем / В.А. Боднер, Н.Е. Роднищев, Е.П. Юриков. – М.: Машиностроение, 1987. – 208 с.
2. Кафаров В.В. Развитие идей перспективного стохастического программирования для задач химической технологии / В.В. Кафаров, В.И. Бодров, В.Г. Матвейкин // ДАН СССР. – 1989. – Т. 308. – № 4. – С. 918–921.
3. Осипов С.Ю. Технологические режимы при восстановлении гуммировочных покрытий с учетом технико-экономических критериев [Текст]: монография / С.Ю. Осипов, Ю.Р. Осипов, С.А. Шлыков; М-во образования и науки Российской Федерации, Вологодский гос. ун-т. – Вологда: ВоГУ, 2016. – 163 с.
4. Рудобашта С.П. Диффузия в химико-технологических процессах / С.П. Рудобашта, Э.М. Карташов. – М.: КолосС, 2010. – 478 c.
5. Юдин Д.Б. Задачи и методы стохастического программирования / Д.Б. Юдин. – М., 1979.

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

В связи с тем, что, как это отмечалось в работе [3], оптимальная стратегия интенсификации процесса вулканизации не инвариантна вариациям возмущающих воздействий, для обеспечения заданного качества вулканизации необходимо учитывать факторы неопределенности. Задача интенсификации данного процесса заключается в повышении качества восстановительного ремонта при минимуме энергетических затрат и определяется температурным режимом его проведения. Использование стохастической модели процесса вулканизации приводит к необходимости постановки и решения задач интенсификации на основе теории оптимального управления случайными процессами.

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

Во многих задачах управления в условиях неполной информации, связанных с повторяющимися ситуациями, нет необходимости в том, чтобы ограничения удовлетворялись при каждой реализации случая. Затраты на накопление информации или другие затраты, обеспечивающие исключение невязок в условии задачи, могут превышать достигаемый при этом эффект. Часто конкретное содержание задачи требует лишь, чтобы вероятность попадания решения в допустимую область превышала заданное число α > 0. В тех случаях, когда возможные невязки в отдельных ограничениях вызывают различный ущерб, имеет смысл дифференцированно подходить к разным условиям. Чтобы уравновесить ущерб, определяемый невязками в разных условиях задачи, следует ограничить снизу вероятность выполнения каждого из них различными числами α > 0. Обычно α1 > 1/2. Подобные постановки задач стохастического программирования называют моделями с вероятностными ограничениями.

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

Результаты оценки адекватности полученной модели [2, 4] доказали ее пригодность в целях моделирования и интенсификации процесса термообработки при ремонте эластомерных обкладок гуммированных аппаратов. В результате получим

osip01.wmf (1)

osip02.wmf osip03.wmf osip04.wmf osip05.wmf – вектор управляющих воздействий; osip06.wmf osip07.wmf osip08.wmf osip09.wmf.

С соответствующими начальными условиями:

osip10.wmf.

Критерий оптимальности имеет вид

osip11.wmf

osip12a.wmf

osip12аа.wmf (2)

где osip13.wmf osip14.wmf – соответственно расчетная и измеряемая в «характерных» точках температура, при условиях связи, представленных в виде уравнения теплопроводности при конкретных граничных условиях и ограничениях на конструктивные параметры и управляющие воздействия

osip15.wmf,

и на мощность нагревателя osip16.wmf osip17.wmf.

Для решения задачи воспользуемся методами оптимизации систем со случайными параметрами [1, 5]. При этом задача (1–2) может быть формализована следующим образом:

osip18.wmf

osip19.wmf

osip20.wmf

osip21.wmf

osip22.wmf (3)

где osip23.wmf, osip24.wmf

osip25.wmf, osip26.wmf представляет собой часть компонентов вектора фазовых координат, не используемых в (2). Определим необходимые условия оптимальности для рассматриваемой задачи, используя методы математической теории экстремальных задач [1].

Вариации δv, где v = {u,а}, в точке v будем называть возможными, если достаточно малое перемещение из точки v∈Θ в направлении вариаций не выводит за пределы множества Θ, т.е. существует такое малое положительное e > 0, (v + εδv)∈Θ.

Функционал (3) является выпуклым и ограниченным, что следует из [3] и для того, чтобы он достигал минимума на Θ в точке v* = {u*,a*}, необходимо и достаточно, чтобы в этой точке не было подходящих вариаций [1].

Решение задачи начинается с определения некоторого исходного un(τ), an, pn, n = 0, 1, 2,...; где n – номер текущей итерации, р = (х0, b).

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

Задача Коши для уравнения (1) решается с заданным приближением un, an. Определяется исходное приближение osip27.wmf и вычисляются значения функционалов М[osip28.wmf] и функций ψs(u), osip29.wmf при каждом фиксированном τ∈[τ0, τ1].

Далее формируется множество индексов

osip30.wmf osip31.wmf

osip32.wmf (4)

osip33.wmfosip34.wmf

osip35.wmf osip36.wmf

osip37.wmf osip38.wmf (5)

где osip39.wmf, множества osip40.wmf совпадают при en = 0, если osip41.wmf, где e – параметры, вводимые для обеспечения сходимости итерационного алгоритма.

Согласно методу возможных направлений [3, 4], введем искусственную переменную s, вектор с положительными компонентами q, и определим конус возможных вариаций:

osip42.wmf

osip43.wmf osip44.wmf (6)

Выберем из всех допустимых вариаций osip45.wmf, osip46.wmf удовлетворяющих (6) вариации, имеющие наибольшую скорость убывания функционала (2). Вводя нормализацию osip47.wmf, osip48.wmf||, приходим к необходимости решения задачи вида

osip49.wmf

osip50.wmf (7)

osip51.wmf

где osip52.wmf – непрерывные случайные функции, определенные следующим образом:

osip53.wmf

Произвольный вектор (osip54.wmf, osip55.wmf, σn), удовлетворяющий (7), при σ ≥ 0 определяет возможные подходящие вариации, так как

osip56.wmf

osip57.wmf osip58.wmf (8)

– есть конус возможных вариаций [ ],

osip59.wmf

Следовательно, за счет максимизации sn выбираются такие подходящие вариации osip60.wmf osip61.wmf в направлении которых функционал (2) убывает, а функция М[Н(t, х, u, а, b, L)] переменного u почти для всех τ∈[τ0, τ1] возрастает достаточно быстро. Основное значение нормализации – это обеспечение единственности решения вспомогательной задачи выбора подходящих вариаций. В направлении наискорейшего спуска osip62.wmf задача о выборе возможных вариаций (7) сводится к отысканию конечномерного вектора (an, dan), обеспечивающего

osip63.wmf osip64.wmf

osip65.wmf osip66.wmf (9)

osip67.wmf

Здесь коэффициенты при неизвестных apn – известные величины, вычисляемые при un, an, pn вдоль n-х траекторий по формулам

osip68.wmf (10)

Случайные величины и вектор-функции, входящие в правые части (9), (10), моделируются по их заданным вероятностным характеристикам. Наиболее простую вычислительную схему имеют алгоритмы, построенные на основе метода Монте-Карло (статистических испытаний). Для каждой выборки значений реализаций случайных векторов и вектор-функций при помощи численных методов интегрируется система (10), причем численное интегрирование проводится многократно [3]. В результате получается некоторое количество решений (10) для различных значений реализаций векторов и вектор-функций.

Задача выбора подходящих вариаций (9) является задачей линейного программирования с двухсторонними ограничениями на переменные. При ее решении возникают затруднения, связанные с тем, что задача линейного программирования (8) имеет бесконечное число ограничений. Поэтому, как известно [1, 4], для того, чтобы преодолеть эту трудность, интервал времени [t0, t1] покрывается счетной сеткой с шагом t, (узлы сетки tv = nt(n = 0,.., t1 = (Nt)] и улучшение управляющей функции производится при каждом фиксированном τv∈τ0,τ1]. Таким образом, решение (9) сводится к приближенному.

osip1.wmf

Кинетические кривые теплового режима прогрева покрытия марки 1751 (1,5) + 2566 (1,5 + 1,5 + 1,5) при δст = 4 мм до и по результатам решения вероятностной задачи

Величина шага hn, на который следует сместиться в направлении подходящих вариаций osip69.wmf, вычисляется как наименьший, положительный корень уравнений

osip70.wmf (11)

osip71.wmf osip72.wmf (12)

Частные производные и функции вычисляются в точке (un(t) + hnosip73.wmf(τ)| τ∈τ’’, an + hnosip74.wmf); hn – единственная переменная; un(t),an, osip75.wmf, osip76.wmf – известные векторы; τ» = ∪sτs. Если какое-либо из уравнений (12) не имеет положительных корней, то его наименьший положительный корень полагается равным ∞.

Определим ΔI0 = I0(u n+1) – I0(un). Поскольку решение задачи имеет в большинстве случаев быструю сходимость по функционалу и несколько медленнее сходится по управлению, то находится также ΔН = M[H(u n+1) – H(un)] для каждого [τ0, τ1]. Производим оценку ΔI0, ΔН. Если |ΔI0| < ε0, |ΔН| < ε0, J(vn, εn) = J(vn), после этого итерационный процесс прекращается. В противном случае переходим к определению un(τ), an, pn и продолжаем вычисления используя новое приближение

osip77.wmf

osip78.wmf

Таким образом, метод возможных направлений предусматривает выбор начального приближения, выбор подходящих вариаций, в направлении которых функционал I0(v) убывает, а функция М[Н] переменного u почти для всех τ∈[τ0, τ1] возрастает достаточно быстро, обеспечение условий сходимости процесса к оптимальному решению, определение длины шага и оценку приближений.

Результаты решения задач интенсификации для различных марок ремонтируемых гуммировочных покрытий [3, 4] представлены на рисунке, где под буквой «а» показаны законы изменения токов нагревательных элементов плит, под буквой «б» – изменение во времени средней по поверхности температуры вулканизируемого участка ремонтируемого покрытия. Кривая 1 описывает изменения средней по поверхности температуры при эксперименте по результатам реализации детерминированной задачи процесса вулканизации, кривая 2 характеризует изменение температуры, полученное при реализации решения вероятностной задачи.

На рисунке «а» представлено управляющее воздействие процессом вулканизации при ремонте покрытия марки 1751 (1,5) + 2566 (1,5 + 1,5 + 1,5) при δст = 4 мм, которое представляет собой кусочно-постоянную функцию. График имеет две ярко выраженные площадки (основной нагрев и остывание в индукторе), что не противоречит физико-химическим основам процесса вулканизации. Первая площадка, находящаяся на уровне 4,2 А, имеет продолжительность 2400 с и соответствует процессу разогрева, вторая – соответствует процессу выдержки.

При сравнении результатов можно сделать обоснованный вывод, что введение в математическую модель формализованных возмущений стохастической природы позволяет существенно снижать максимальные значения токов (до 5–7 %) и, как следствие, обеспечивать меньшие затраты электроэнергии.

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


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

Осипов Ю.Р., Демакова А.О., Шилова А.О., Ночевалова Ю.А., Осипов С.Ю., Богданов Д.А. ЗАДАЧА ИНТЕНСИФИКАЦИИ ТЕПЛООБМЕНА ПРИ ВОССТАНОВЛЕНИИ ЭЛАСТОМЕРНЫХ ПОКРЫТИЙ НА ОСНОВЕ ТЕОРИИ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ // Фундаментальные исследования. – 2016. – № 12-1. – С. 96-101;
URL: https://fundamental-research.ru/ru/article/view?id=41053 (дата обращения: 29.03.2024).

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

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