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

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

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

Энерготехнологический комплекс «нагрев - обработка давлением» предназначен для производства изделий широкой номенклатуры (валки, профили) из стандартных металлических заготовок, как правило, цилиндрической формы. Технология прессования через матрицу на гидравлических прессах предварительно нагретых заготовок, в свою очередь, является наиболее широко распространенной для получения различных видов профилей из алюминиевых сплавов.

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

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

, (1)

согласно которому температура в фильере матрицы (), т.е. максимальная температура в зоне деформирования, не должна превышать допустимого предела , что гарантирует отсутствие на изделиях трещин основного вида необратимого брака изготавливаемых полуфабрикатов [1, 2].

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

. (2)

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

(3)

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

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

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

Процедура непосредственной -параметризации управляющих воздействий проводится на основе метода последовательной параметризации [4].

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

Рассмотрим в качестве объекта управления нестационарный процесс прессования цилиндрической заготовки, описываемый двумерным уравнением теплопроводности Фурье в относительных единицах вида (4) с краевыми условиями (5)-(6) и внутренними источниками тепла пластического формоизменения :

(4)

(5)

, (6)

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

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

. (7)

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

. (8)

Соответствующее заданному значению целевое множество может быть описано в форме требования достижения с допустимой погрешностью предельного значения температуры в каждый момент процесса прессования :

.(9)

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

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

.(10)

Выбирая начальное состояние объекта (4)-(6) в качестве управляющего воздействия, осуществляется переход к следующей задаче оптимального управления по начальным условиям.

Задача 1. Для объекта, описываемого неоднородным уравнением Фурье (4) с краевыми условиями (5)-(7), необходимо определить такое оптимальное распределение температуры по длине заготовки перед процессом прессования, рассматриваемое в качестве начального условия, которое обеспечивает минимальное значение критерия (10) в каждый момент процесса прессования на временном интервале в условиях выполнения ограничения (1).

Особенность рассматриваемой задачи заключается в формировании функции максимума в (9), (10) не в пространственной, а во временной области изменения управляемой величины.

Редукция к задаче математического программирования

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

В роли рассматриваемых параметров , могут выступать непосредственно параметры функций, описывающих искомое одномерное температурное распределение .

Для решения задачи 1 используем, в пренебрежении радиальной неравномерностью начального температурного состояния, -параметризованное представление управляющего воздействия в форме кусочно-линейной функции, повторяющей требуемый характер распределения температуры по длине слитка при его градиентном нагреве, который обеспечивает достаточно малую величину :

(11)

Зависимость (11) определяет набор температурных профилей, представленных на рис. 1. В качестве оптимизируемых параметров выступают значения температур , на концах y=0, y=1 заготовки по ее длине, температура в точке и сама координата .

Рис. 1. Температурные профили нагреваемой заготовки перед прессованием

В результате параметризованное управление (11) характеризуется вектором , : ; ; ; , где для S=2 принимается , , а при S=3 выполняются соотношения , .

При выборе в виде упорядоченной определенным образом последовательности S параметров минимально достижимые значения в классе управлений, однозначно характеризуемых величиной S, монотонно убывают с возрастанием и образуют последовательность неравенств вида:

.(12)

Решение уравнения объекта (4) с -параметризованным начальным условием и соответствующими граничными условиями позволяет получить температуру в фильере матрицы и значение критерия оптимальности (10) в виде явных функций, соответственно и , вектора . В результате без каких-либо погрешностей в рамках используемых моделей осуществляется точная редукция задачи (10) к виду:

. (13)

Решение задачи (13) для каждого , где по определению:

, (14)

может быть получено по общей схеме альтернансного метода.

Метод базируется на специальных альтернансных свойствах вектора оптимальных решений, являющихся аналогом известных условий экстремума в теории нелинейных чебышевских приближений [3]. Особенностью задачи (14) является ее формулировка в терминах ошибок равномерного приближения управляемой величины к не в пространственной, а во временной области ее определения в связи с минимизацией функции максимума (13) во временной области изменения температуры в фильере матрицы. Указанное обстоятельство требует проведения дополнительных исследований альтернансных свойств по форме оптимальной кривой изменения во времени температуры пластической зоны в процессе прессования с постоянной скоростью.

Альтернансные свойства оптимального процесса изменения во времени температуры пластической зоны.

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

Согласно этим свойствам, равные величине одинаковые значения максимальных отклонений:

,

достигаются в некоторые моменты времени на временном интервале , общее число которых должно быть равным числу искомых неизвестных в (13).

Это означает, что альтернансные свойства решения задачи (13) позволяют записать систему соотношений:

; (15)

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

Таким образом, в (15):

. (16)

Если дополнить равенства (15) условиями существования экстремума максимума в моменты времени , которые заранее не определены, то мы приходим к следующей системе соотношений:

(17)

замкнутой в условиях (16) относительно всех неизвестных величин ; и при , где - точки экстремума на отрезке , общее число которых .

Решение задачи трехпараметрической оптимизации.

Задача 1 решалась с использованием численной модели процесса прессования металла, построенной на основе известных зависимостей для поля скоростей деформации [4, 5].

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

Данному случаю соответствует трехпараметрическое представление температурного распределения при (рис. 1, в). При этом оптимизируемыми параметрами являются температуры и на торцах заготовки, а также температура в точке с заранее фиксированной координатой y01.

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

В соответствии с правилом (16) , т.е. максимальное отклонение температуры в очке матрицы от заданной должно достигаться в процессе прессования четыре раза, а число экстремумов кривой не должно превышать [4, 5].

Рис. 2. Температурные профили нагреваемой заготовки перед прессованием и общий вид зависимости температуры пластической зоны от времени в процессе прессования: а - , б -, в -

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

гидравлический прессование алюминиевый сплав

Установленные альтернансные свойства температурных зависимостей в процессе прессования с заданной скоростью заготовок цилиндрической формы позволяют найти решение задачи 1 с помощью специальных систем уравнений (17).

Результаты решения задачи оптимизации процесса прессования с постоянной скоростью VП = 450 мм/мин цилиндрической заготовки диаметром D=250 мм и длиной L=430 мм из сплава Д16 представлены на рис. 2 и приведены в таблице в абсолютных единицах для температуры в фильере матрицы оптимизируемых параметров профиля и точности соответственно.

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

Литература

Рапопорт Э.Я. Оптимизация процессов индукционного нагрева металла. М.: Металлургия, 1993. - 279 с.

Перлин И.Л., Гайтбарг Л.Х. Теория прессования металлов. - М.: Металлургия, 1975. - 448 с.

Рапопорт Э.Я. Альтернансный метод в прикладных задачах оптимизации. М.: Наука, 2000. - 336 с.

Плешивцева Ю.Э. Последовательная параметризация управляющих воздействий и полубесконечная оптимизация алгоритмов управления технологическими объектами с распределенными параметрами: Автореф. дис. … д-ра техн. наук: 05.13.06. / Самара: Самар. гос. техн. ун-т, 2009. - 50 с.

Афиногентов А.А. Моделирование и оптимальное управление технологическим комплексом «нагрев - обработка металла давлением»: Автореф. дис. … канд. техн. наук: 05.13.06. / Самара: Самар. гос. техн. ун-т, 2008. - 22 с.

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