CC BY
101
2008
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА Математика и механика
№ 1(2)
МЕХАНИКА
УДК 539.3:621.396.67
А.В. Бельков, С.В. Пономарев
МОДЕЛИРОВАНИЕ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ НАДУВНОЙ КОНСТРУКЦИИ КОСМИЧЕСКОГО РЕФЛЕКТОРА
Рассматривается моделирование напряженно-деформированного состояния надувной конструкции космического рефлектора.
Ключевые слова: моделирование, надувная космическая антенна, рефлектор, напряженно-деформированное состояние, метод конечных элементов.
Развитие систем спутниковой связи делает необходимым разработку и создание новых космических рефлекторов. За рубежом активно развивается направление надувных конструкций, отверждаемых в космосе. На рис. 1 показана экспериментальная надувная рефлекторная антенна диаметром 14 м с куполом, под низким давлением, изготовленным из лавсановой пленки толщиной 6,5^10-6 м, торусом под высоким давлением и тремя стойками под высоким давлением, разработанная L’Garde, Inc. [1, 2]. Соединение купола с торусом при помощи растяжек позволило обеспечить требуемые характеристики отражающей поверхности рефлектора.
Рис. 1. Летная надувная рефлекторная антенна (разработана L’Garde, Inc.)
Следует отметить работы по исследованию подобных конструкций [3, 4]. Задача о колебании обода (торуса) большой космической антенны рассматривается в работах А.В. Лопатина и М.А. Рутковской [4].
На рис. 2 изображена проектная модель 50-метровой надувной антенны, разрабатываемая L’Garde, Inc. [5]. Данная экспериментальная модель была взята за основу (в плане конструкции) для моделирования космического надувного рефлектора.
Надувной купол с радиоотражающим покрытием
Надувной обод
Солнечные панели
Система наполнения (силовая установка, система коммуникаций система надува)
Адаптивная система питания
Радиочастотное малое зеркало
Рис. 2. Проектная модель надувной антенны ARISE (разработана L’Garde, Inc.)
Математическая постановка задачи
Моделирование подобных конструкций требует геометрически нелинейной постановки задачи механики деформируемого твердого тела (МДТТ) с учетом температурных деформаций.
Связь деформаций с перемещениями рассматривается в виде
еар = (1/2)(иар + Щ,а + ие,аие,р). (1)
Для моделирования механического поведения в напряженном состоянии, которое возникает при раскрытии рефлектора, использованы упрощенные зависимости между напряжениями и деформациями. Компоненты тензора напряжений
Кирхгофа и компоненты тензора деформаций связаны зависимостью
а5 =а**( - е’ав), (2)
где а!]И = а]'к1(Х,о(-0')]') - элементы матрицы упругости, зависящие от принадлежности к разнородным элементам конструкции и уровня предварительных напряжений; = -ЭЛ75а|3 - компоненты тензора температурных деформаций; О - коэф-
фициент линейного расширения; АТ = Т - Т0.
Уравнения равновесия записываются в виде
[(8* + + Р, = 0, (3)
где Р, - компоненты вектора массовых сил, 8,у - символ Кронекера.
Краевые условия для надувных конструкций рефлекторов
Для конструкций рефлектора граничные условия в перемещениях задаются на части поверхности Би конструкции рефлектора, где имеет место крепление к космическому аппарату:
и = и; (4)
на внешней части поверхности конструкции Бр в условиях космического пространства
0>еп](8у + иу) = 0. (5)
При моделировании наземных условий
0е И/(8,у + Иу) = ра,
где ра - атмосферное давление на поверхности Земли.
Таким образом, при известном температурном поле в элементах конструкции рефлектора задача МДТТ считается поставленной в замкнутой форме, так как количество определяемых функций перемещений соответствует количеству разрешающих уравнений.
Для надувных конструкций на внутренних поверхностях надувных элементов, где действует давление газа р, имеет место краевое условие
ое и,(8,у + и,*) = р. (6)
В приведенной постановке исходным состоянием является конфигурация рефлектора, вытекающая из технического задания с нулевыми напряжениями в элементах конструкции. После задания требуемых краевых условий в элементах конструкции получается напряженная конфигурация, находящаяся в состоянии статического равновесия и обеспечивающая с достаточной точностью проектируемую форму отражающей поверхности рефлектора. Аналитическое решение подобных задач представляется затруднительным. Для получения приближенных численных решений выгоднее опираться на вариационные постановки. Приведенной дифференциальной постановке соответствует эквивалентная вариационная постановка в виде принципа виртуальной работы:
I{(а ))8е,у - Р 8г}с1У- / ^ 8г =0, (7)
где 8г - вариация вектора перемещений; Р - вектор массовых сил; ^ - вектор поверхностных сил. На части поверхности Su заданы нулевые перемещения, что со-
ответствует закреплению конструкции.
Моделирование
Для решения поставленной задачи использован метод конечных элементов [6]. Так как основным фактором для зеркальных антенн является форма отражающей поверхности, то основное внимание направлено на центральную купольную часть надувного рефлектора. Построенная в программном комплексе ANSYS конечно элементная трехмерная модель купола надувного космического рефлектора дает возможность оценить перемещения точек отражающей поверхности при надувании купола рефлектора. На рис. 3 приведены графики, на которых изображены перемещения точек отражающей поверхности в зависимости от радиуса конструкции при одинаковом значении давления в куполе и толщине материала, из которого он изготовлен.
г/а
Рис. 3. Зависимости относительных перемещений точек отражающей поверхности от радиуса конструкции. По оси У - перемещения, отнесенные к радиусу конструкции, по оси X - безразмерный радиус. Давление р = 10 Па, толщина пленки к = 0,0009 м, а - радиус конструкции
При надувании подобных конструкций в материале возникают сжимающие напряжения, приводящие к искажению формы отражающей поверхности на краю купола, что показано на рис. 4. Чтобы избежать искажения отражающей поверхности по периметру надувного купола в радиальном направлении прикладывается внешняя растягивающая сила Г, которая в реальной конструкции есть силовое воздействие от окружающего тора.
Для оценки величины данной силы использована упрощенная конечноэлементная двумерная осесимметричная модель. Чтоб показать адекватность данной модели был выполнен расчет с одинаковыми параметрами, такими как толщина пленки, диаметр конструкции, давление, для двух случаев - двумерного и трехмерного.
Рис. 4. Краевые искажения отражающей поверхности надувного рефлектора при внутреннем давлении 10 Па (для наглядности перемещения увеличены в 2 раза)
На рис. 5 изображена та же модель что и на рис. 4, но с приложенной по периметру в радиальном направлении силой.
Рис. 5. Отражающая поверхность рефлектора с приложенными по периметру внешними растягивающими силами
На рис. 6 изображена зависимость приложенной по периметру в радиальных направлениях растягивающей погонной силы Г, при которой в материале остаются только растягивающие напряжения, отнесенной к толщине материала Н, из которого изготовлен купол, и модулю упругости, от диаметра конструкции.
ПНЕ 0,25
0,2
0,15
0,1
0,05
0
0 5 10 15 20 25 30 35 40 45 Д м
Рис. 6. Зависимость приложенной растягивающей силы для удаления искажений отражающей поверхности (ПкЕ, где к = 0,0009 м) от диаметра конструкции
На рис. 7 приводится зависимость напряжений, возникающих в материале надувного купола, от величины, приложенной растягивающей погонной силы.
ае/£
0
-0,5-10 3 -1-10 3 -1,5-10 3 -2-10 3 -2,5-10 3
1-10 5 1-10 4 1-10 3 1-10 2 1-10 1 Е/ИБ
Рис.7. Зависимость напряжений в надувном куполе от приложенной в радиальном направлении растягивающей силы в двумерном случае
Результаты моделирования
Результаты компьютерного моделирования показывают, что для оболочечных надувных конструкций имеет место масштабный фактор. Одной стороной его проявления является то, что при увеличении диаметра рефлектора мы наблюдаем увеличение перемещений точек купола рефлектора (рис. 3), приводящее к увеличению отклонений от заданной формы. Другой стороной масштабного фактора является рост напряжений в материале купола при увеличении размеров конструкции, что ведет к большей потере формы особенно по краю купола (рис. 4). Ввиду таких изменений формы отражающей поверхности приходим к выводу, что необходим учет этих факторов при раскрое и создании купола для обеспечения требуемых геометрических параметров. С этой же целью введена внешняя радиальная погонная сила, моделирующая силовое воздействие окружающего тора. Проведенный анализ необходимой величины этой силы для формирования нужных растягивающих напряжений позволяет получить требуемую форму отражающей поверхности.
ЛИТЕРАТУРА
1. FreelandR.E. Significance of the inflatable antenna experiment technology // AIAA-98-2004.
2. FreelandR.E., Bilyeu G.D., Veal G.R., Steiner M.D., Carson D.E. Large inflatable deployable antenna flight experiment results // IAF-97-1.3.01.
3. Greschik G., Mikulas M.M., Palisoc A. Torus-less inflated membrane reflector with an exact parabolic center // AIAA. 2004. V. 42. No. 12. P. 2579 - 2584.
4. Лопатин А.В., Рутковская М.А. Оценка жесткости обода большой космической антенны // Королевские чтения: Тез. докл. Всерос. молодеж. науч. конф. Т.1. Самара: Изд-во СНЦ РАН, 2003. С. 22 - 23.
5. Salama M., Kuo C.P., Garba J., Wada B. Thomas M. On-orbit shape correction of inflatable structures // Jet Propulsion Laboratory California Institute of Technology Pasadena, CA. 91109.
6. Зенкевич О. Метод конечных элементов в технике / Под ред. Б. Е. Победри. М.: Мир, 1975. 541 с.
Принята в печать 31.03.08.
