Rational photovoltaic batteries mounting on the earth remote sensing microsatellite

Full Text

Космические средства дистанционного зондирования Земли (ДЗЗ) являются одним из основных перспективных направлений развития космических систем. Космические системы ДЗЗ интенсивно развиваются и становятся неотъемлемой частью информационного обеспечения многих направлений развития экономики государств [1]. Современные тенденции в области ДЗЗ, разработка и запуск мини- и микроспутников, активное внедрение спутниковых данных в работу государственных учреждений и высшей школы, развитие перспективных программ космического мониторинга, повышающие требования в части детальности и оперативности космической съемки, ставят задачу более детального исследования отдельных систем КА [2]. Учитывая высокую интенсивность развития и внедрения в работу космических аппаратов для дистанционного зондирования Земли, актуальность приобретает проектирование и разработка систем электроснабжения (СЭС) космических аппаратов данного класса, а также оптимизация их по массе с ограничением по стоимости и надежности. В соответствии с выполняемыми задачами рассматриваемый класс космических аппаратов (КА) имеет следующие особенности: орбиты солнечносинхронные, т. е. направление ПОС (угол между плоскостью орбиты спутника и направлением на Солнце) всегда практически постоянно; из-за небольших размеров КА панели фотоэлектрических батарей (БФ), как правило, неориентируемые; одна из осей КА преимущественно ориентирована в надир, формируя местную вертикаль. На этапе проектирования представляет интерес зависимость потребной массы СЭС, при ограничении стоимости и высокой надежности СЭС, от методов размещения неподвижных панелей солнечных батарей относительно местной вертикали КА для различных параметров орбиты (высота и местное время прохождения восходящего узла). Учитывая высокую надежность современных комплектующих элементов силовой электроники, входящих в СЭС КА, в задачу проектирования входит рациональное расположение батарей фотоэлектрических для стабильного за виток коэффициента использования и мощности при минимальной массе и стоимости. 119 Авиационная и ракетно-космическая техника Для того чтобы оценить потребную массу и стоимость СЭС космического аппарата, проведем анализ рациональной установки БФ для штатной работы СЭС в удельных параметрах мощности фотоэлектрической батареи и емкости химической батареи (БХ). Для этого решения большинство вариантов расположения панелей БФ можно свести к конусу с углом раскрытия от 0° до 180°. При угле раскрытия 0° имеем цилиндрическую поверхность, ось цилиндра направлена на центр Земли. При угле раскрытия 180° имеем плоскость, перпендикулярную направлению на центр Земли [3]. Поэтому в качестве рабочей модели геометрии БФ для дальнейших расчетов был выбран конус, поверхность которого покрыта фотоэлектрическими элементами. Основными варьируемыми параметрами были угол между плоскостью орбиты спутника и направлением на Солнце (ПОС - угол в) и угол раскрытия конуса (2у). Исходя из параметров орбиты космического аппарата (высота и местное время прохождения восходящего узла, которым определяется угол ПОС) и угла установки панелей, был рассчитан коэффициент использования солнечной батареи за виток: Ки 1 Твит Т- f КЕосв(T)dТ Т вит 0 (1) К (Т) = _f Коr.вdS , бок (2) dS - r'dФ'dr sin у : Sо f RdRd ф sin у 5 С учетом неравномерности освещенности эффективная освещённая площадь (условный эквивалент площади, поглощающий ту же энергию при условии нормального падения на него светового потока) составит R +фкр D 7D 7 ^св.эфф ( т) - R f RSR-r'c0s ^(тХ 0 -фкр 1 (6) С учетом выражений (3)-(6) интегральный коэффициент освещенности площади поверхности конуса (К£ осв) КЕ осв (Т) - S0CB. эфф (Т) ‘■Е осв (7) где S&K - площадь боковой поверхности конуса. где Твит - длительность витка, мин; К осв - суммарный коэффициент освещенности всей фотоэлектрической батареи в каждый момент времени на витке; т - единица времени. Суммарный коэффициент освещенности определяется выражением: где Sq^ - площадь боковой поверхности конуса; Косв -мгновенное значение коэффициента освещенности для единичного участка площади на поверхности конуса, он равен cos Q. (рис. 1), где Q. - это угол между направлением на Солнце (вектор S) и вектором нормали к плоскости произвольного единичного элемента на поверхности конуса (векторN); dS - единица площади поверхности конуса: (3) где r - радиальная координата участка dS; ф - угловая координата вектора участка dS относительно оси х. Площадь освещенного участка в каждый момент времени (S^) найдем с помощью выражения R +фкр . , , ...........(4) 0 -фкр где фкр - угол раскрытия для конического сектора освещенной части конуса; R - радиус основания конуса. Для каждого угла а (угол между осью х и вектором направления на Солнце S ) значение фкр находится из условия cos Q = 0. Его можно определить как образующую вдоль поверхности конуса, для которой вектор нормали N перпендикулярен направлению на Солнце: фкр - arccos (tg а' tg у). (5) Рис. 1. Расчетная модель геометрии БФ В результате расчетов были получены зависимости коэффициента использования солнечного излучения от варьируемых параметров: угла между плоскостью орбиты спутника и направлением на Солнце (ПОС - угол в) и угла раскрытия конуса (2у). Для отображения полученных результатов в дальнейшем использовался угол у, т. е. половина полного угла раскрытия. Эти два значения являются определяющими количества солнечной энергии, поступающей на единицу площади. На рис. 2 приведен результат расчета ожидаемого максимального коэффициента использования солнечного излучения для различных значений угла у при высоте орбиты 300 км. Значения у варьируются в пределах 0°-90°. Максимальное значение коэффициента использования солнечного излучения составляет 0,32-0,33 и достигается в крайних положениях плоскости орбиты КА относительно Солнца (0° и 90°) и крайних значениях угла установки панелей (плоскость или цилиндр), при соотношениях в = 0° и у = 90° и наоборот. Минимальное значение Кисп для любых у наблюдается в пределах 40°< в < 60° независимо от высоты и варьируется в пределах 0,22-0,23. 120 Вестник СибГАУ. № 1(53). 2014 Чтобы посчитать необходимую массу БФ и БХ, сначала найдем мощность, генерируемую БФ: NБФ(т) - КТсов(т)• S5 где Е - солнечная постоянная ( : = 1?6Ф ’Вт/м2); ПыР-коэффициент полезного действия фотоэлектрической батареи; Sбок - площадь боковой поверхности конуса (условно принимаем как 1 м2). Тогда среднюю мощность нагрузки за виток можно принять как N 1 БФ ср ' Т 1 вит J NБФ (T)dТ. (9) NБХ (т) - NБФ i - NБФ ср, (10) где NБФ j - мощность БФ в определенный момент времени. 6б -Щ БХ min Uраб • 60 (12) бБх min не учитывает необходимого запаса на деградацию емкости БХ в процессе работы, а также ограниченные возможности реальной БХ воспринимать большие величины зарядного тока. Рекомендуемое значение максимального зарядного тока для литий-ионных аккумуляторов составляет величину 0,5С. Среднее значение зарядного тока на освещенном участке: дБХ min I,. T (13) Циклограмма нагрузки на борту КА редко носит четко детерминированный характер. Как правило, для больших промежутков времени включение потребителей происходит с нечеткой логикой по конкретной задаче на данное время функционирования. Поэтому по закону больших чисел с некоторым допущением можно принять в качестве циклограммы нагрузки постоянное значение (NM ср). Тогда поток мощности через накопитель (Nex) равен разности мощностей генератора и текущего значения мощности нагрузки (которое было принято как константа (NM ср)): где Тосв - длительность освещенного участка на витке. Тогда необходимый коэффициент запаса по емкости Кд, обеспечивающий приемлемые значения зарядного тока, определяем из условия I • 1 ч -*зар.ср 1 2 - 0 5 Кд 'дБХ min C учетом выражения (14) 1 Кд (14) (15) д Т • 0 5С ' 1 осв Соответственно, установленная емкость БХ берется с учетом коэффициента запаса: дБХуст - Kg дБХ min. (16) Массу СЭС космического аппарата составляют масса БХ, масса БФ и масса управляющего и коммутационного оборудования. Удельные характеристики этих элементов были приняты на основе данных [4], где представлены доступные для приобретения образцы. Для определения массы и стоимости БФ за основу был взят солнечный элемент NanoPower Solar P100U-A/B, масса которого составляет 0,059 кг, размеры 82,5 мм*98 мм, установленная мощность 2,3 Вт, КПДБФ 30 %. Принятые расчетные удельные характеристики - 0,03 кг/Вт, 46610 у. е./кг [5]. Для определения массы и стоимости БХ, управляющего и коммутационного оборудования за основу был взята аккумуляторная батарея для микроспутников NanoPower P31U Power Supply. Масса батареи 0,095 кг, общая масса 0,2 кг, установленная емкость 5,2 Ач [6]. у, град Рис. 2. Изменение максимального значения коэффициента использования солнечного излучения для различных значений угла ПОС Проинтегрировав по времени выражение (6), получим зависимость от времени текущей заряженности аккумулятора: ЩБХ (t)-J N БХ (тУ Т. (11) 0 Минимально необходимую для нормальной работы СЭС электрохимическую емкость БХ (А ч) с учетом рабочего напряжения аккумуляторной батареи определим по формуле где Wex max - максимальное значение емкости БХ за виток; Wex min - минимальное значение емкости БХ за виток; ираб - напряжение БХ. Рис. 3. Зависимость изменения емкости химических батарей от массы Дополнительно был также проведен анализ аккумуляторов фирмы SAFT космического назначения. Результаты представлены в виде графика (рис. 3). 121 Авиационная и ракетно-коcмическая техника В результате проведенного анализа была получена величина удельной емкости дуд для единичного аккумулятора 40 Ач/кг. На основе имеющихся данных рассчитана масса аккумуляторной батареи (МБХ): Мбх - дБХ^• n, (17) Ууд где дуд - удельная величина емкости химической батареи; n - количество аккумуляторов. Принятые расчетные удельные характеристики -удельная масса 0,007 кг/Вт, удельная стоимость 16500 у. е./кг. Массу управляющего и коммутационного оборудования определили следующим образом: Муп,ком - NБФ max • mуп,ком. уд , (18) где NБФ max - максимальная мощность БФ за виток; myn,™. уд - удельная масса управляющего и коммутационного оборудования. В дальнейшем все рассчитанные значения массы и стоимости были отнесены к средней мощности нагрузки. Главной задачей расчетов было определить зависимость изменения массы и стоимости СЭС космического аппарата класса микроспутник для дистанционного зондирования Земли от параметров орбиты космического аппарата, в частности, угол между плоскостью орбиты спутника и направлением на Солнце и угла установки панелей [7]. В результате проведенных расчетов были получены следующие зависимости массы от угла установки панелей на высоте орбиты 300 км. Для наглядности полученных результатов при построении графиков использовалось удельное значения массы (кг/Вт). Из полученных графиков (рис. 4) можно сделать вывод, что основной вклад в массу всей системы энергоснабжения космического аппарата делает масса БФ, а наиболее незначительной является масса БХ. При угле между плоскостью орбиты спутника и направлением на Солнце 45° суммарная масса СЭС остается практически постоянной независимо от угла установки панелей и варьируется в пределах 0,1-0,12 кг/Вт. При в = 0° масса максимальна при установке панелей под углом 0° и достигает значения 0,16 кг/Вт, а при 40° < у < 90° изменяется в пределах 0,09-0,1 кг/Вт. Если угол ПОС = 90°, то с точки зрения массы наиболее выгодным является положение солнечных батарей под углом 0° < у < 55°. В этих пределах масса СЭС варьируется от 0,06 до 0,1 кг/Вт, при дальнейшем увеличении угла установки БФ ее масса, а следовательно, и масса всей СЭС резко возрастает. На рис. 5 приведен график изменения суммарной массы СЭС от угла установки панелей (у) для различных значений угла ПОС. Угол в варьировался в пределах 0°-90°. Минимальная масса СЭС (0,06 0,065 кг/Вт) наблюдается при соотношениях крайнего положениях плоскости орбиты относительно Солнца (в = 90°) и крайнего значения угла установки панелей (У = 0°). Значительному увеличение массы соответствуют углы установки солнечных батарей 55°-90° при 75° < в < 90°. Масса системы энергоснабжения микроспутника изменяется в пределах 0,09-0,16 независимо от углов в и у. С учетом стоимости NanoPower Solar P100U-A/B (2750 у. е.) и стоимости NanoPower P31U Power Supply (3300 у. е.) были оценены затраты на систему энергоснабжения для различных значений углов между плоскостью орбиты и направлением на Солнце и углов установки солнечных батарей (рис. 6). Также была получена зависимость изменения стоимости СЭС с учетом стоимости запуска от угла установки панелей для различных значений угла ПОС (рис. 6). Y. град а Y, град б у, град в Рис. 4. Изменение массы СЭС в зависимости от различных значений угла установки: а - при ПОС 0°; б - при ПОС 45°; в - при ПОС 90° 122 Вестник СибГАУ. № 1(53). 2014 Y> I рил Рис. 5. Изменение массы СЭС в зависимости от различных значений угла установки для разных значений угла между плоскостью орбиты спутника и направлением на Солнце Y. град а у,град б Y. град в Рис. 6. Изменение стоимости СЭС в зависимости от различных значений угла установки: а - при ПОС 0°; б - при ПОС 45°; в - при ПОС 90° Y, град Рис. 7. Изменение стоимости СЭС в зависимости от различных значений угла установки для разных значений угла ПОС Полученные графики (рис. 5, 6) повторяют предыдущие результаты и еще раз подтверждают, что стоимость СЭС для 0° < в < 60° варьируется в пределах 4000-8500 у. е. независимо от угла установки панелей. Минимальная (3300-3600 у. е.) и максимальная (13000-14500 у.е.) стоимость достигается при крайних значениях: в = 90° - у = 0° и в = 90° - у = 90° соответственно. Представленная методика дает возможность рассчитать ожидаемое значение максимального коэффициента использования солнечного излучения для различных значений высоты орбиты, угла установки БФ (угол у) и местного времени прохождения восходящего узла, определяемого углом ПОС (угол в). На основе полученных результатов были определены необходимые для нормальной работы СЭС мощность фотоэлектрической батареи и емкость химической батареи, в соответствии с которыми получены графики изменения массы и стоимости СЭС космического аппарата ДЗЗ для различных значений углов у и в. Масса и стоимость СЭС КА были рассчитаны для различных значений угла установки панелей БФ (0-90°) и угла между плоскостью орбиты спутника и направлением на Солнце (0-90°). Полученные данные позволяют сделать вывод, что при угле ПОС, равном 45°, суммарная масса СЭС остается практически постоянной независимо от угла установки панелей и варьируется в пределах 0,1-0,12 кг/Вт. При в = 0° масса максимальна при установке панелей под углом 0° и достигает значения 0,16 кг/Вт, а при 40° < у < 90° изменяется в пределах 0,09-0,1 кг/Вт. Если угол ПОС 90°, то с точки зрения массы наиболее выгодным является положение солнечных батарей под углом 0° < у < 55°. В этих пределах масса СЭС варьируется от 0,06 до 0,1 кг/Вт, при дальнейшем увеличении угла установки БФ значительно увеличивается ее масса, а следовательно, и масса всей СЭС. Графики, отображающие зависимость удельной стоимости СЭС от угла установки панелей для различных значений угла ПОС, повторяют полученные результаты для массы и иллюстрируют, что минимальная стоимость СЭС (3300-3600 у. е.) возможна при условии установки солнечных батарей под углом 0° и угле ПОС 90°.

About the authors

Sergii Victorovich Gubin

National Airspace University Named after N. E. Zhukovskiy “Kharkov Aviation Institute”

Email: gubinsv@d4.khai.edu
Candidate of Engineering Science, associate professor, Deputy head of the Department of Rocket-space engines and power plants of space crafts

Yurii Alekseevich Shepetov

National Airspace University Named after N. E. Zhukovskiy “Kharkov Aviation Institute”

Email: shepetov9@d4.khai.edu
Candidate of Engineering Science, associate professor of Rocket-space engines and power plants of space crafts department, Deputy Laboratory Chair of microsatellite designing

Katerina Olexandrivna Dolzhukova

National Airspace University Named after N. E. Zhukovskiy "Kharkov Aviation Institute"

Email: katerina.dolzhukova@yandex.ua
postgraduate student of rocket-space engines and power plants of space crafts department

References

Supplementary files