CC BY
74
УДК 539.3:621.396.67
А. В. Бельков, А. И. Величко, С. В. Пономарев, В. А. Солоненко
МОДЕЛИРОВАНИЕ НАДУВНОГО КОСМИЧЕСКОГО РЕФЛЕКТОРА
Рассматривается компьютерное моделирование крупногабаритного трансформируемого рефлектора КА надувного типа, основанного на постановке задачи механики деформируемого твердого тела.
Ключевые слова: компьютерное моделирование, напряженно-деформированное состояние, надувной рефлектор.
Моделирование крупногабаритных трансформируемых и надувных конструкций рефлекторов является актуальным направлением разработки и создания конструкций систем спутниковой связи и зондирования поверхности Земли, так как экспериментальная отработка подобных конструкций требует больших материальных и временных затрат. Для решения современных задач требуются все более крупногабаритные антенны с высокой точностью формы зеркала рефлектора. Разработка таких конструкций, позволяющих улучшить качество приема и увеличить размеры зон сканирования, является одной из главных целей в области космических программ. Основные требования к конструкциям рефлекторов заключаются в высокой точности формы отражающей поверхности и наведения, высокой температурной стабильности и радиоотражающей способности антенных систем.
Модель надувного рефлектора ARISE, разрабатываемая компанией L'Garde, Inc, представлена на рис. 1 [1].
Надувной купол с
радиоотражающим
покрытием
Солнечные панел
Система наполнения (силовая установка, система коммутаций, система надува)
Надувной обод
Адаптивная система питания
радиочастотное малое зеркало
Рис. 1. Проектная модель надувной антенны ARISE
Задача о напряженно-деформированном состоянии подобных вантово-оболочечных конструкций формулируется с позиций нелинейной теории упругости и в вариационной постановке [2]. Решение задачи осуществлялось методом конечных элементов [3].
Связь деформаций с перемещениями рассматривается в виде
e<xp = (1/2)(Mx,p + «Р.а + «е.аИе,р X
Компоненты тензора напряжений Кирхгофа и компоненты тензора деформаций связаны следующей зависимостью:
ст*' = a jxp
(exp ex p).
Элементы матрицы упругости, зависящие от принадлежности к разнородным элементам конструкции и уровня предварительных напряжений рассчитываются по формуле
ajxxp = aikl (X, ст(0)»).
Компоненты тензора температурных деформаций имеют вид
eXp=JAT8„p.
Для конструкций рефлектора граничные условия в перемещениях задаются на части поверхности S конструкции рефлектора, где имеет место крепление к космическому аппарату
u = u
на внешней части поверхности конструкции Sp в условиях космического пространства
G% (Sj + И, j ) = 0.
При моделировании наземных условий выполняется равенство
s4 (Sj+«,-, j)=p^
где pa - атмосферное давление на поверхности Земли.
Таким образом, при известном температурном поле в элементах конструкции рефлектора задача МДТТ считается поставленной в замкнутой форме, так как количество определяемых функций перемещений соответствует количеству разрешающих уравнений. Для надувных конструкций на внутренних поверхностях надувных элементов, где действует давление газа p, имеет место краевое условие
СТЧ- (Sj + u, j) = pj.
Аналитическое решение подобных задач представляется затруднительным. Для получения приближенных численных решений выгоднее опираться на вариационные постановки:
Jj(a'jxP (eap- eXp ))Seff - P •Sr}dV - J (p2 - p, )n -Sr • dS ,
где P - вектор массовых сил; F - вектор поверхностных сил.
Трехмерная модель надувного космического рефлектора, построенная в программном комплексе ANSYS, представлена на рис. 2. За основу взята модель надувной антенны ARISE (см. рис. 1).
Исходным состоянием является конфигурация рефлектора, вытекающая из технического задания с нулевыми напряжениями в элементах конструкции. После задания требуемых краевых условий в элементах конструкции получается напряженная конфигурация, находящаяся в состоянии статического равновесия и обеспечиваю-
щая с достаточной точностью проектируемую форму отражающей поверхности рефлектора. Перемещения точек конструкции и возникающие в материале напряжения после приложения нагрузки в нашем случае давления показаны на рис. 3 и 4.
Рис. 2. Конечно-элементная 3Б-модель надувного рефлектора
Рис. 3. Отклонения точек конструкции относительно начального ненапряженного состояния
объясняется значительно меньшей толщиной материала (рис. 3, 4). Исходя из этого и учитывая, что наиболее жесткие требования выдвигаются к отражающей поверхности, в дальнейшем будем рассматривать только ее с целью сокращения объемов вычислений.
Компьютерное моделирование показывает, что для оболочечных надувных конструкций имеет место масштабный фактор. Одной стороной его проявления является то, что при увеличении диаметра рефлектора и при уменьшении давления в куполе мы наблюдаем увеличение перемещений точек отражающей поверхности рефлектора относительно идеального параболоида.
Графики, на которых изображены перемещения точек отражающей поверхности в зависимости от радиуса конструкции при одинаковом значении давления в куполе и толщине материала, из которого он изготовлен, приведены на рис. 5.
Рис. 4. Напряжения в элементах конструкции относительно начального ненапряженного состояния
Из проведенных вычислений следует, что основная нагрузка приходится именно на купол рефлектора, что
Рис. 5. Зависимости относительных перемещений точек отражающей поверхности от радиуса конструкции. По оси У - перемещения, отнесенные к радиусу конструкции, по оси X - безразмерный радиус. Давление р = 10 Па; толщина пленки Н = 900 мк; а - радиус конструкции
При надувании подобных конструкций в материале возникают сжимающие напряжения, приводящие к искажению формы отражающей поверхности на краю купола (рис. 6). Чтобы избежать искажения отражающей поверхности по периметру надувного купола, в радиальном направлении прикладывается внешняя растягивающая сила Г, которая в реальной конструкции есть силовое воздействие от окружающего тора.
Зависимость приложенной по периметру в радиальных направлениях растягивающей погонной силы Г, при которой в материале остаются только растягивающие напряжения, отнесенной к толщине материала Н, из которого изготовлен купол, и модулю упругости, от диаметра конструкции изображена на рис. 7.
Зависимость напряжений, возникающих в материале надувного купола, от величины, приложенной растягивающей погонной силы, приводится на рис. 8.
Взяв только нижнюю часть купола и закрепив его по периметру, посмотрим, как меняется геометрия отражающей поверхности, и какие напряжения возникают в материале.
Отклонения точек отражающей поверхности от идеального параболоида после приложения давления (надува конструкции) приводятся на рис. 9. Отсюда видно, что даже при таких граничных условиях имеют место довольно серьезные смещения точек отражающей поверхности. СКО составляет 141 мм.
Рис. 6. Краевые искажения отражающей поверхности надувного рефлектора при внутреннем давлении 10 Па (для наглядности перемещения увеличены в 2 раза)
Рис. 9. Отклонения отражающей поверхности от идеального параболоида до регулировки ее начальной геометрии
Еще одной стороной масштабного фактора является рост напряжений в материале купола при увеличении размеров конструкции, что ведет к еще большей потере формы. Напряжения, возникающие в материале купола после приложения давления, показаны на рис. 10.
¿3 ^ 43
Рис. 7. Зависимость приложенной растягивающей силы для удаления искажений отражающей поверхности (¥!НЕ, где Н = 0,000 9 м) от диаметра конструкции
: и." : ]и I :и - ш и.1 I.
Рис. 8. Зависимость напряжений в надувном куполе от приложенной в радиальном направлении растягивающей силы в двумерном случае
Рис. 10. Напряжения, возникающие в материале отражающей поверхности
Решение этих проблем осуществляется использованием в конструкции более жесткого надувного тора [4], фиксирующего купол с помощью растяжек и позволяющего корректировать форму отражающей поверхности, а также с помощью регулирования начальной формы отражающей поверхности рефлектора, т. е. изначально нам нужна поверхность с геометрией, несколько отличной от вырезки из родительского параболоида. Под регулировкой начальной формы подразумевается определение такой исходной геометрии
э
отражающей поверхности, что при подаче газа в купол рефлектора мы получим рабочую поверхность с требуемыми радиотехническими характеристиками. После проведения регулировки СКО равно 3,5 мм, т. е. находится в области допустимых значений. Отклонения точек отражающей поверхности от идеального параболоида после регулировки отображены на рис. 11. Видно, что центральная область отражающей поверхности имеет отклонения порядка 1 мм, но ближе к краю купола картина ухудшается. Это связано с кривизной поверхности офсетного параболоида и проиллюстрировано на рис. 9. На нем хорошо видна область максимальных отклонений.
ДЕнньеоЕегичинеЕьпагненной регулировки позволяют изготоЕить отражающую поверхность, котораябу-дег удовлетворять всшраqиаraниче:кимтрй:авEHияvl
Гоместивотражаш.уюгоЕфхность впогесилтяже-сти и выполнив расчаыдгядв/х погт&ий «нЕиавв^» и «чаша вниз» (см. рис. 12), мы выяснили, что разница погаж^ий еетсче< между состоянием в невесомости и пслож^иями «чаша ввфх» и «чаша вниз» совпадает с тaчнaoтъюдaтраьffa> знака пaсгезEПяraй Расхождение
- \ ■■■(! - Т1ТП - С*:г" -пи---!-
Рис. 11. Отклонения отражающей поверхности от идеального параболоида после регулировки ее начальной геометрии
составляет менее 2 %, что может быть объяснено кривизной поверхности.
В результате проведенных расчетов приходим к выводу, что хотя на отражающей поверхности и присутствуют области с достаточно большими отклонениями от идеального параболоида, но это не мешает нам добиться того, чтобы значение СКО находилось в требуемых пределах (порядка 3.. .5 мм).
Библиографические ссылки
1. Greschik G., Mikulas M. M., Palisoc A. Torus-less inflated membrane reflector with an exact parabolic center // AIAA. 2004. Vol. 42. №> 12. Р. 2579-2584.
2. Васидзу К. Вариационные методы в теории упругости и пластичности. М. : Мир, 1987.
3. Зенкевич О., Победри Б. Е. Метод конечных элементов в технике. М. : Мир, 1975.
4. Лопатин А. В., Рутковская М. А. Оценка жесткости обода большой космической антенны // Королевские чтения : тез. докл. Всерос. молодеж. науч. конф. Т. 1. Самара : Изд-во СНЦ РАН, 2003. С. 22-23.
1 ,>4111 .1-11Ы
11^1^1 . нкьа
Рис. 12. Разница положений точек отражающей поверхности между положениями в невесомости: «чаша вверх» - шкала a; «чаша вниз» - шкала Ь
A. V Belkov, S. V Ponomarev, A. I. Velichko, V A. Solonenko MODELLING OF THE INFLATABLE SPACE REFLECTOR
Computer modeling of large-sized deployable reflector of the inflatable type for spacecraft, based on statement of problem of mechanics of a deformable firm body is considered in this article.
^^^words: mmputer modeling, intense-deformed condition, inflatable reflector.
© Бельков А. В., Величко А. И., Пономарев С. В., Солоненко В. А., 2010
