Automated control system of heat fluxes at thermovacuum testing of spacecrafts

Texto integral

Современная стратегия экспериментальной отработки космических аппаратов (КА) основывается преимущественно на наземной отработке систем и частей КА. При этом имитация штатных условий эксплуатации производится в наиболее полном возможном объеме [1]. Тепловакуумные испытания (ТВИ) являются важной составляющей предполетной отработки КА [2]. Во время ТВИ имитируются такие условия космоса, как вакуум, солнечный поток, резкие перепады температуры по поверхности КА. Современные системы для проведения ТВИ имеют возможность использования сотен устройств подвода мощности (УПМ) и датчиков температуры, что делает необходимым автоматизацию задачи получения определенного температурного режима на испытываемом КА. Представленная система предназначена для проведения различных тепловакуумных испытаний с количеством контрольных точек до 750 и управляющих воздействий до 380 (каналов). Работы по разработке автоматизированной системы велись в рамках создания горизонтальной вакуумной установки ГВУ-600 для ОАО «ИСС» имени академика М. Ф. Решетнева». Система осуществляет автоматическое управление ходом испытаний в крупноразмерной горизонтальной 137 Авиационная и ракетно-коамическая техника вакуумной камере ГВУ-600 с объёмом более 600 м3 (рис. 1), обеспечивает анализ нештатных ситуаций и автоматическое управление алгоритмом испытаний, контролирует температуру различных точек изделия в диапазоне от -150 до +150 °С, позволяет проводить в автоматическом режиме один из видов тепловакуумных испытаний: термобалансные или термостатические. Термобалансные испытания подразумевают подачу постоянной мощности независимо от температуры либо релейное управление мощностью в заданном температурном диапазоне. Термостатические (термоциклические) тепловакуумные испытания заключаются в последовательном выводе испытываемого КА на определенные температурные режимы (ТР) и поддержание заданных ТР в течение требуемого времени. Под температурным режимом подразумевается такое состояние объекта, при котором температура, измеряемая в контролируемых точках, существенно не меняется со временем. Рис. 1. Крупноразмерная горизонтальная вакуумная камера ГВУ-600 с объёмом более 600 м3 Задача получения заданного температурного режима. Сложность проведения испытаний заключается в том, что любой источник тепла в общем случае оказывает влияние на любую точку объекта, вызывая изменение температуры на ней. Поэтому практически невозможно перейти к конечному числу одноконтурных регуляторов. В общем случае существует N управляемых УПМ и T датчиков, фиксирующих температуру в интересующих точках испытываемого изделия. Влияние изменения подводимой мощности на изменение температуры в каждой точке неизвестно. F(nb n2, ..., nN) = (tb t2, ..., tT) - неизвестная зависимость ТР от значений подводимых мощностей. Входными данными являются: исходный ТР (t{0], t20], ..., tT0]), исходное значение подводимых мощностей (n[0], n20], ..., n[N), требуемый ТР (?/L], ?2[L], ..., tTL) и точность 5t получения ТР. Искомым является вектор значений подводимых мощностей («l], n2L], ..., n[N), при котором получается ТР (t{L], t2L], ..., tTL]), отличающийся от требуемого не более чем на величину 5t в каждой контролируемой точке: F(,,[L], „24 ..., »iL]) = (<IL], 4L], ..., <TL]), Vi e{1, 2.....N} |(/L] - t[L]| < dt. Алгоритм управления тепловыми потоками. Решение задачи основано на нахождении приращения F в текущей точке, при изменении мощности на каждом УПМ на известную величину An,, i e {1, 2, ., N}, и вычислении вектора коэффициентов (k1, k2, ..., kN) в предположении, что F локально линейна, а точнее, что F(n1 +k1An1, «2 +k2An2, nN +kNAnN)-F(nl,n2, .., nN) = =k1 •(F(n +Anb И2, .., %)-F(,«2, nN))+ (1) +kN •(F(n1, n2, ..., nN +A«n)-F(n1, n2, ..., nN)). Реализация алгоритма состоит в последовательном изменении мощности на каждом УПМ на известную величину и ожидании выхода системы на ТР. Для ожидания ТР используется три параметра - время ожидания реакции системы, время ожидания баланса и допуск на баланс. Время реакции системы определяет интервал, в течение которого после изменения мощности не контролируется выход на ТР. Допуск на баланс задает минимальное отклонение температуры в любой контролируемой точке за время, меньшее времени ожидания баланса, при котором констатируется температурная нестабильность системы. Ожидание ТР заканчивается, если за время ожидания баланса температура во всех точках изменилась на величину, не превышающую допуска на баланс. После получения ТР для текущего УПМ мощность на нем возвращается в исходное состояние и меняется мощность следующего. В результате этого процесса получается набор ТР: F(np] + A1, n20].....^Н", г21].....T) ... (2) F(, (].....«N'+An .....Г), t [i] где tj - температура, зарегистрированная на j-м датчике при изменении мощности на i-м УПМ. Из (1) и (2) получаем систему из T уравнений с N неизвестными: 138 Вестник СибГАУ. № 1(53). 2014 (/[1] _ ЛЦ t[2] _ 7[L] h h h h t[1] _ 7[L] t[2] _ 7[L] 2 2 2 2 t[1] _ /M t[2] _ 7_[i] V‘ r t[N] _ 7[L] Л t[N] _ 7[L] l2 l2 [N] _ r[L] 2 V kN J (t[0] _ 7[L] \ t[0] _ J[L] l2 l2 t[0] _~t[L] LT ir- V r (3) r J После успешного решения системы (3) [3] значение подводимых мощностей устанавливается в (n[L], 4L].....n[LN )- - (wj[0] + k1An1,4°] + k2Дп2, ..., n[N +kNAnN) и ожидается выход системы на ТР (/{ [L] t[L] 2 [ L ] )• В общем случае (исходя из практического предположения, что Т > N) система (3) переопределена, т. е. любое решение имеет невязку, которая связана с нелинейностью F, неточностью определения выхода системы на ТР, а также с погрешностью измерения температуры, нелинейностью УПМ и неконтролируемым подводом мощности. Еще одним фактором, влияющим на наличие и величину невязки, является теоретическая достижимость требуемого ТР. Система датчиков и УПМ может быть сконфигурирована так, что запрашиваемый ТР не может быть получен ни при каких значениях мощностей. Невязку, вызванную нелинейностью F, можно уменьшать путем последовательного применения вышеописанного алгоритма, добиваясь на каждом этапе минимизации невязки и беря в качестве исходного состояния для следующего этапа состояние, полученное на текущем. Результаты. Нами проведены опытные испытания алгоритма управления. В результате получены данные, подтверждающие работоспособность алгоритма. В эксперименте было задействовано семь каналов мощности (рис. 2) и пять контрольных точек измерения температуры (рис. 3). В соответствии с регламентом описанного выше алгоритма сначала было проведено последовательное изменение мощности на каждом УПМ на известную величину с ожиданием выхода системы на ТР. После получения ТР для текущего УПМ мощность на нем возвращалась в исходное состояние и менялась мощность следующего (правая часть графиков на рис. 2 и 3 до времени 16.00). Затем в соответствии с (3) устанавливались вычисленные значения подводимых мощностей и ожидался выход системы на ТР (правая часть графиков на рис. 2 и 3 после времени 16.00). Как видно из этих рисунков, выход на заданный ТР успешно осуществлен в автоматическом режиме. [Г г -?пп 1 2 3 4 5 / i' \ / . \ R 1- V 1 у ® Г / V 1 Vu А j Л 11 .1 J . , . * 1 ш JI /11 i\\ Л 07:00:00 10:00:00 13:00:00 16:00:00 19:00:00 22:00:00 01:00:00 04:00:00 07:00:00 21112012 21 112012 21.11.2012 21.11.2012 21.11.2012 21.11.2012 22.11.2012 22.11.2012 22.11.2012 Рис. 2. Графики зависимости мощности от времени для семи УПМ, участвующих в эксперименте по проверке действия алгоритма 07:00:00 10:00:00 13:00:00 16:00:00 19:00:00 22:00:00 01:00:00 04:00:00 07:00:00 21.11.2012 21.11.2012 21.11.2012 21.11.2012 21.11.2012 21.11.2012 22.11.2012 22.11.2012 22.11.2012 Рис. 3. Графики зависимости температуры от времени для пяти контролируемых точек, используемых в эксперименте по проверке действия алгоритма 139 Авиационная и ракетно-коамическая техника В работе представлена система управления тепловыми потоками при проведении тепловакуумных (термобалансных или термостатических) испытаний космических аппаратов с количеством контрольных точек до 750 и управляющих воздействий до 380 (каналов). Она обеспечивает анализ нештатных ситуаций и автоматическое управление алгоритмом испытаний, контролирует температуру различных точек изделия в диапазоне от -150 до +150 °С. Приведен разработанный алгоритм автоматического вывода испытываемого объекта на заданное температурное распределение. Результаты опытных испытаний алгоритма свидетельствуют о его работоспособности.

Sobre autores

Andrei Shevlyakov

Technological Design Institute of Scientific Instrument Engineering, Siberian Branch of Russian Academy of Sciences

Email: ghost2613@mail.ru
senior programmer

Ignat Vykhristyuk

Technological Design Institute of Scientific Instrument Engineering of Siberian Branch of Russian Academy of Sciences

Email: uic@ngs.ru
junior scientific worker

Aleksander Verkhoglyad

Technological Design Institute of Scientific Instrument Engineering of Siberian Branch of Russian Academy of Sciences

Email: verhog@tdisie.nsc.ru
head of laboratory

Vladimir Khalimanovich

JSC “Information Satellite System” named after academician M. F. Reshetnev”

Email: office@iss-reshetnev.ru
Candidate of Phisical and Mathematical Sciences, Professor, Director of the Industrial Center for Large-Sized Foldable Mechanical Structures

Valery Hristich

JSC “Information Satellite System” named after academician M. F. Reshetnev”

Email: hristich@iss-reshetnev.ru
associate professor, chief specialist

Bibliografia

Arquivos suplementares