Управление температурой быстрых термических процессов

Данная статья является продолжением темы быстрых термических процессов, которой я занимаюсь в рамках диссертационной работы. Логично, что после задачи идентификации, которая рассматривалась в предыдущем топике, стоит постараться решить задачу управления процессами такого рода.Установка представляет собой высокотемпературную водоохлаждаемую камеру, в которой происходит осаждение элементов из металлоорганических соединений на сапфировые подложки. Рост производится при температурах от 550 до 1000 градусов Цельсия и атмосферном давлении. Полученные таким образом полупроводники обладают надежностью и высоким КПД, что позволяет использовать их в космических технологиях.

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

c26338cde63446bf987ef4b04fd61717.jpgОдним из основных элементов управляющей цепи является высокочастотный генератор с индуктивным элементом.

be20773d5c844b808f197b01207f63e5.jpg

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

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

5b9ddef397464f3da3755e8339166d7d.jpg

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

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

a43e649d369cd2b8384725daf9c7658f.png

где T — измеряемая температура, P — управляющая мощность индуктора, B — коэффициент при управляющем воздействие, A1 — коэффициент характеризующий лучистое распространение энергии, A2 — коэффициент характеризующий конвекцию, С — скорость нагрева точки измерения от окружающей среды. В ходе экспериментов выявилось, что коэффициент A1 зависит в основном от типа материала, из которого изготовлен приемник (подложкодержатель), A2 растет с увеличением температуры и скорости продувки химических реагентов через лайнер (кварцевый кожух, который локализует зону химической реакции). Получается, для качественного результата идентификации параметров A1 и C нужно измерять температуру газа на входе в реактор и на выходе из него, а так же, не помешало бы, знать температуру лайнера. Но это пока остается в планах на будущее.

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

Для 550 градусов: 029dd295b245df2c4ee145fa50f9365a.png Для 1000 градусов: 6ebe1c2cf293b79f179f0f389fe9e9fc.png Все изменения значений параметров имеют физическую природу. В частности, при увеличении температуры графита меняются его электромагнитные характеристики. Мы видим снижение эффективности преобразования электромагнитной энергии в тепло на 17%. Оставшиеся коэффициенты сами по себе зависят от температуры.

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

ea78a244fe6d4e26a2b2b22489db6c38.png

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

Далее приведены некоторые результаты работы для 550 и 1000 градусов Цельсия. В качестве управления подавалось ступенчатое воздействие с амплитудой 20 градусов.

93750bc238474dedb1a210fdafaddb48.jpga1eb5a49b8b3489fa3f5f6cb13783665.jpg

Стоит отметить, что пропорциональный коэффициент регулятора завышен, так как видно перерегулирование в 2 градуса в режиме удержания 550 градусов. Математическое моделирование и экспериментальные данные имеют следующий график невязки для режима 1000 градусов. Выглядит весьма неплохо, хотя трудно говорить о качестве модели замкнутой системы.

8917bbc3ea3b4b8ab88dbe6ad56e06ab.jpg

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

958acc8cfb8845aba1d3027acfbd94e0.jpg

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

© Habrahabr.ru