CC BY
119
Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. — 2008. — № 2 (17). — С. 161—170
Математическое моделирование
УДК 539.3:621.396.67
КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ПЕРСПЕКТИВНЫХ КОСМИЧЕСКИХ РЕФЛЕКТОРОВ
А. В. Бельков1, А. С. Евдокимов1, С. В. Пономарев1, Д. Б. Усманов2
1 НИИ прикладной математики и механики при Томском государственном университете, 634050, г. Томск, пр-т Ленина, 36.
2 ОАО «Информационные спутниковые системы им. ак. М. Ф. Решетнева», 662972, г. Железногорск, ул. Ленина, 52.
E-mails: khirurgasibmail.com, eas1985amail.ru, psvh3psy.tsu.ru
Рассматривается комплексная методика компьютерного моделирования обод-ных рефлекторов, основанная на связанных постановках задач механики деформированного твердого тела и электродинамики, а также компьютерного моделирования напряжённо-деформированного состояния надувной конструкции космического рефлектора.
Ключевые слова: напряжённо-деформированное состояние, метод конечных элементов, трансформируемый сетчатый рефлектор, диаграмма направленности, надувная антенна, коэффициент направленного действия.
Введение. Для современной спутниковой связи требуются крупногабаритные антенны с высокой точностью формы зеркала рефлектора. Использование трикотажного металлического сетеполотна в качестве отражающей поверхности позволяет получить улучшенные удельные массовые характеристики для больших габаритов конструкций рефлекторов. Рассматривается трансформируемый космический рефлектор с силовой схемой, включающей ферменный обод, вантовую систему и отражающую поверхность. Ферменный обод представляет собой стержневую конструкцию, собранную из жёстких углепластиковых элементов.
Варианты методик расчёта радиотехнических характеристик крупногабаритных рефлекторов рассматривались в работах J. Ruze (см., например, [1]). В монографии М. В. Гряника [2] приводится расчёт для зонтичных конструкций, однако учёт искажений отражающей поверхности производился на основе экспериментальных измерений.
1. Механическое моделирование. Особенность подобных конструкций и необходимость учёта температурных деформаций при эксплуатации в кос-
Бельков Алексей Викторович — младший научный сотрудник НИИ прикладной математики и механики при Томском государственном университете.
Евдокимов Александр Семёнович — младший научный сотрудник НИИ прикладной математики и механики при Томском государственном университете.
Пономарев Сергей Васильевич — зав. лаб. НИИ прикладной математики и механики при Томском государственном университете, к.ф.-м.н.
Усманов Давид Бисенович — главный специалист ОАО «Информационные спутниковые системы им. ак. М. Ф. Решетнева,».
мосе требуют общей геометрически нелинейной постановки задачи с учётом предварительной напряжённости элементов конструкции и температурных деформаций. Связь деформаций с перемещениями в элементах конструкции рефлектора рассматривается в виде
eaß = 1 (Ua,ß + Uß,a + Ue,aU0ß).
Как показывают результаты экспериментов, материалы, используемые при создании элементов конструкций крупногабаритных рефлекторов, являются анизотропными и физически нелинейными. Однако для моделирования механического поведения в напряжённом состоянии, которое возникает в раскрытом рефлекторе, можно использовать упрощенные зависимости между напряжениями и деформациями.
Компоненты тензора напряжений Кирхгофа и компоненты тензора деформаций связаны зависимостью
= aijae(eaß — eTaß),
где aijkl = atjkl(X, a((0",%j) —элементы матрицы упругости, зависящие от принадлежности к разнородным элементам конструкции и уровня предварительных напряжений; е^ = §ATöaß — компоненты тензора температурных деформаций; § — коэффициент линейного расширения; AT = T — To.
Уравнения равновесия для ободной конструкции и для надувной конструкции рассматриваются в виде
[(öij + uhj)(aej + a(0)&j)] + Рг = рщ, [(0г] + Ui,3 )a°j ]fi + Рг = 0,
где Рг —компоненты вектора массовых сил; р — удельная масса материала; щ — компонента ускорения; — компоненты тензора предварительных на-
пряжений в элементах конструкции. Задание требуемых напряжений в элементах вантовой конструкции необходимо для создания жёсткости, обеспечивающей с достаточной точностью форму отражающей поверхности рефлектора. При этом важно, чтобы форма ответственных элементов напряжённой вантовой конструкции мало отличалась от проектной для реализации требуемых радиотехнических характеристик антенны.
Добавление начальных и граничных условий делает постановку задачи с позиций механики деформируемого твердого тела полной. В результате такой постановки получается нелинейная задача. Граничные условия в перемещениях задаются на части поверхности Su конструкции рефлектора, где имеет место крепление к космическому аппарату:
иг = иг,
на большей части поверхности конструкции Sa, свободной от нагрузок:
aej nj + a(0)kj nk u$,j = 0.
Начальные условия:
иг(Х, 0) = иг(0); Uг(X, 0) = Щг(0).
Таким образом, задача считается поставленной в замкнутой форме, так как количество определяемых функций перемещений соответствует количеству разрешающих уравнений.
В приведённой постановке исходным состоянием является конфигурация рефлектора, вытекающая из технического задания, с нулевыми напряжениями в элементах конструкции. После задания требуемых напряжений в элементах конструкции получается напряжённая начальная конфигурация, находящаяся в состоянии статического равновесия с требуемой жёсткостью, обеспечивающей с достаточной точностью проектируемую форму отражающей поверхности рефлектора.
Приведённой дифференциальной постановке соответствует эквивалентная вариационная постановка в виде принципа виртуальной работы [3]:
(а(0)^ + (еав - е0в- Р¿г + с1У = 0,
где 5т — вариация вектора перемещений; Р — вектор массовых сил; Ь — время. На части поверхности Би заданы нулевые перемещения, что соответствует закреплению конструкции. Аналитическое решение подобных задач не представляется возможным.
Для получения численного решения задачи о напряжённо-деформированном состоянии рефлектора использовался метод конечных элементов [4]. Решение проводилось методом конечных элементов, реализованным в программном комплексе ЛЫВУБ. Средствами ЛЫБУВ была построена геометрическая и конечно-элементная модели космического рефлектора с фермой натяжения (рис. 1).
Влияние температурных нагрузок на напряжённо-деформированное состояние конструкции моделировалось через зависимость напряжений от деформаций внутри каждого конечного элемента. Решение подобной нелинейной задачи позволяет получить равновесную форму отражающей поверхности и напряжённо-деформированное состояние. Методика решения нелинейной задачи основывалась на ранее разработанной методике для зонтичных рефлекторов [5].
В качестве обобщённой меры отклонения отражающей поверхности рефлектора в равновесном состоянии использовалось среднеквадратическое значение отклонений (СКО) полученной расчётной поверхности в узлах конечно-элементной сетки от поверхности соответствующего параболоида.
2. Электродинамическое моделирование. Определение основных радиотехнических характеристик антенн связано с получением выражения для электромагнитного поля в дальней зоне, когда источниками поля являются заданные сторонние токи ] на отражающей поверхности рефлектора [7].
Рис. 1. Конечно-элементная модель рефлектора
Система уравнений Максвелла имеет вид
гоЬИ = —шеЕ + 3, ю1Е = ш^И.
Токовый метод определения направленных свойств антенны базируется на известном распределении поверхностных токов на внутренней поверхности зеркала. Вектор плотности тока с учётом коэффициента отражения в данной точке поверхности зеркала можно определить с учётом ориентации векторов И в падающей и отраженной волнах по формуле
3 а = [поИ ]Я(а),
где 3 3 — вектор плотности поверхностного тока в данной точке зеркала; И — вектор напряжённости магнитного поля, создаваемого падающей волной облучателя в данной точке на поверхности зеркала; по — орт нормали к поверхности зеркала в этой же точке; Я(а) — коэффициент отражения, зависящий от напряжения а, реализующегося в отражающей поверхности из сетеполот-на.
Напряжённость первичного магнитного поля облучателя определяется по формуле
И = В е-,
r
где B — величина, не зависящая от r и характеризующая направленные свойства облучателя; r — расстояние от фазового центра облучателя до точки, в которой определяется поле. Для построения картины токов, возникающих на отражающей поверхности рефлектора под влиянием поля облучателя, необходимо знать распределение вектора B в пространстве, то есть векторную диаграмму направленности (ДН) облучателя.
Зная распределение тока на поверхности зеркала, можно определить направленные свойства параболической антенны. Для этого необходимо проинтегрировать по всей поверхности зеркала выражение для напряжённости поля, создаваемого элементом поверхности зеркала, рассматривая его как элементарный электрический вибратор.
Поле излучения параболоида можно представить в виде
E = -ЩТр1 е-гвГ1 jsx exp(-q)dS,
где Sa — поверхность параболоида, используемая в качестве антенны; E — напряжённость поля, созданного токами основной поляризации. Угол q определяется относительно луча, идущего прямо от облучателя до точки приёма.
Диаграмма направленности параболической антенны с круглым симметричным относительно фокальной оси раскрывом в секторе направлений, близких к направлению максимального излучения, определяется соотношением ф° 2п • ф
/ / sin / \
F (9) = ! I jsx cos^ exp (-ОД sin 9 tg ф cos(á - p)J dpd-ф.
0 0
Для расчёта радиотехнических характеристик применялся пакет FEKO — система 3D электромагнитного моделирования. Идеальная и деформированная модели рефлектора были импортированы из программного комплекса ANSYS в программный комплекс FEKO при помощи файла ASCII.
3. Результаты механического и электродинамического анализов. В результате механического моделирования было получено напряжённое деформированное состояния рефлектора, позволяющее говорить о достаточной равномерности натяжения сетеполотна и, соответственно, о возможных незначительных отклонениях величины коэффициента отражения радиоволн от оптимального значения (рис. 2). Полученные результаты согласуются с результатами работы [6].
Рис. 2. Интенсивность напряжения отражающей поверхности, Па
Одной из основных характеристик антенн является диаграмма направленности - зависимость излучаемого поля от положения точки наблюдения. Сечения остронаправленных ДН удобнее и точнее изображать в прямоугольной системе координат, поскольку угловой масштаб здесь может быть выбран произвольно в соответствии с шириной ДН.
На рис. 3 приведены диаграммы направленности идеальной параболической поверхности и расчётной равновесной формы отражающей поверхности для ободного рефлектора. Граничные условия соответствовали закреплению рефлектора в точках соединения со штангой от космического аппарата. На-гружение конструкции производилось температурным полем, которое привело к соответствующему деформированному состоянию отражающей поверхности рефлектора и изменению положения оси рефлектора. В полученных результатах ширина диаграммы направленности антенны для идеальной параболической поверхности равна 0° 12' 00'' при частоте 8,0 ГГц.
1
0,8
0,6
0,4
0,2
0
-1 -0,8 -0,6 -0,4 -0,2 0 0,2 0,4 0,6 0,8 1
Рис. 3. Диаграммы направленности для идеальной параболической поверхности и расчётной равновесной формы отражающей поверхности рефлектора с учётом СКО: 1 — идеальный параболоид; 2 — СКО = 0,5 мм; 3 — СКО = 2,0 мм; 4 — СКО = 6,5 мм
В результате вычислений можно сделать вывод, что при увеличении СКО от параболического профиля диаграмма направленности существенно отличается от идеальной диаграммы.
Во-первых, увеличивается ширина диаграммы направленности от 0° 12' 00'' до 0° 18' 00''; во-вторых, возрастают боковые лепестки; в-третьих, происходит смещение главного лепестка относительно центра; в-четвёртых, существенно падает коэффициент усиления с 35 Дб до 25 Дб.
Расчёт диаграммы направленности и других характеристик антенны, основанный на приближённом определении распределения токов на её рабочей поверхности, обеспечивает достаточно точные результаты в пределах главного лепестка антенны и прилегающих к нему боковых лепестков.
В связи с тем, что за рубежом активно развивается направление надувных конструкций, отверждаемых в космосе, была поставлена задача о моделировании таких рефлекторов. Для моделирования в плане конструкции была выбрана проектная модель 50-ти метровой надувной антенны, разрабатываемой L'Garde, ¡по., с куполом, под низким давлением, изготовленным из лавсановой пленки толщиной 6,5 ■ 10-6 м, торусом под высоким давлением и тремя стойками под высоким давлением. Соединение купола с торусом при помощи растяжек позволило обеспечить требуемые характеристики отражающей поверхности рефлектора.
Следует отметить работы по исследованию подобных конструкций. Задача о колебании обода (торуса) большой космической антенны рассматривается в работах А. В. Лопатина и М.А. Рутковской [8].
4. Краевые условия для надувных конструкций рефлекторов. Для конструкций рефлектора граничные условия в перемещениях задаются на части поверхности Би конструкции рефлектора, где имеет место крепление к космическому аппарату:
иг — иг;
на внешней части поверхности конструкции Бр в условиях космического пространства
ав3 п] + и^) = 0.
При моделировании наземных условий использовались соотношения
авз п (5гз + иг,^) = ра,
где ра — атмосферное давление на поверхности Земли.
Таким образом, при известном температурном поле в элементах конструкции рефлектора задача механики деформируемого твёрдого тела считается поставленной в замкнутой форме, так как количество определяемых функций перемещений соответствует количеству разрешающих уравнений.
Для надувных конструкций на внутренних поверхностях надувных элементов, где действует давление газа р, имеет место краевое условие
а<)3 пз (5Ц + иг,з) = Рз.
В приведённой постановке исходным состоянием является конфигурация рефлектора с нулевыми напряжениями в элементах конструкции. После задания требуемых краевых условий в элементах конструкции получается напряжённая конфигурация, находящаяся в состоянии статического равновесия и обеспечивающая с достаточной точностью проектируемую форму отражающей поверхности рефлектора. Аналитическое решение подобных задач представляется затруднительным. Для получения приближенных численных решений выгоднее опираться на вариационные постановки. Приведенной дифференциальной постановке соответствует эквивалентная вариационная постановка в виде принципа виртуальной работы:
{агзав(еав — ета13))5егз — Р5г}йУ — ^ ¥5г = 0,
где 5г — вариация вектора перемещений; Р — вектор массовых сил; ¥ — вектор поверхностных сил. На части поверхности Би заданы нулевые перемещения, что соответствует закреплению конструкции.
5. Моделирование. Для решения поставленной задачи использован метод конечных элементов [4]. Так как основным фактором для зеркальных антенн является форма отражающей поверхности, то основное внимание направлено на центральную купольную часть надувного рефлектора. Построенная в программном комплексе ЛЫБУВ конечно-элементная трёхмерная модель купола надувного космического рефлектора даёт возможность оценить перемещения точек отражающей поверхности при надувании купола рефлектора. На рис. 4 приведены графики, на которых изображены перемещения точек отражающей поверхности в зависимости от радиуса конструкции при одинаковом
иг/а 0,02
0,015
0,01
0,005
0
Й-
-X-
-о -О--о-
-ф--6-
-о
-0-
-х-
-О--о-
-X-
-о -0-о-
-X-
-о
-О-
50 м
30 м
10 м
-о
1 м
0
0,2
0,4
0,6
0,8
г/а
Рис. 4. Зависимости относительных перемещений пх точек отражающей поверхности от радиуса г конструкции: ось У — перемещения, отнесённые к радиусу конструкции; ось X — безразмерный радиус; давление р =10 Па; толщина плёнки — Н = 0,0009 м; а —
радиус конструкции
-.001451- --.410Е-03 ЯИЩШЯ НЩ ..002-111
'•ШУИ <&&вю 'ШШ5 -ШЩ
Рис. 5. Краевые искажения отражающей поверхности надувного рефлектора при внутреннем давлении 10 Па (для наглядности перемещения увеличены в 2 раза)
-0,001
-0,002
Р/(кЕ)
Рис. 6. Зависимость напряжений в надувном куполе от приложенной в радиальном направлении растягивающей силы в двумерном случае
значении давления в куполе и одинаковой толщине материала, из которого он изготовлен.
При надувании подобных конструкций в материале возникают сжимающие напряжения, приводящие к искажению формы отражающей поверхности на краю купола, что показано на рис. 5. Чтобы избежать искажения отражающей поверхности по периметру надувного купола, в радиальном направлении прикладывается внешняя растягивающая сила F, которая в реальной конструкции есть силовое воздействие от окружающего тора.
При введении данной силы видно, что напряжения в материале купола стали уменьшаться, что и проиллюстрировано на рис. 6. Здесь напряжения отнесены к модулю упругости E с целью обезразмеривания величины.
6. Результаты моделирования. Результаты компьютерного моделирования показывают, что для оболочечных надувных конструкций имеет место масштабный фактор. Одной стороной его проявления является то, что при увеличении диаметра рефлектора мы наблюдаем увеличение перемещений точек купола рефлектора (рис. 4), приводящее к увеличению отклонений от заданной формы. Другой стороной масштабного фактора является рост напряжений в материале купола при увеличении размеров конструкции, что ведет к большей потере формы, особенно по краю купола (рис. 5). Ввиду таких изменений формы отражающей поверхности приходим к выводу, что необходим учёт этих факторов при раскрое и создании купола для обеспечения требуемых геометрических параметров. С этой же целью введена внешняя радиальная погонная сила, моделирующая силовое воздействие окружающего тора. Проведенный анализ необходимой величины этой силы для формирования нужных растягивающих напряжений позволяет получить требуемую форму отражающей поверхности.
Заключение. На основе связанной задачи механики деформируемого твердого тела и электродинамики реализована комплексная методика компьютерного моделирования ободных трансформируемых космических рефлекторов, позволяющая более точно учитывать форму и напряжённость отражающей сетчатой поверхности. Компьютерное моделирование позволяет сократить объем экспериментальных работ при создании оптимальных конструкций космических рефлекторов и прогнозировать эффективность функционирования их в условиях космического пространства.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Ruze, J. Antenna Tolerance Theory — A Review [Text] / J. Ruze // Proc. of the IEEE. — 1966. — Vol. 54. — P. 633-640.
2. Гряник, М. В. Развертываемые зеркальные антенны зонтичного типа [Текст] / М. В. Гря-ник, В. И. Ломан. — М.: Радио и связь, 1987. — 72 с.
3. Сахаров, А. С. Метод конечных элементов в механике твердых тел [Текст] / А. С. Сахаров, И. Альтенбах. — Киев: Вища школа, 1982. — 480 с.
4. Зенкевич, О. Метод конечных элементов в технике [Текст] / О. Зенкевич.—М.: Мир, 1975. — 541 с.
5. Моделирование температурных деформаций рефлекторов космических аппаратов [Текст] / В. Г. Бутов, С. В. Пономарев [и др.] // Физика. — 2004. — № 10. (Приложение)— С. 10-18.
6. Усманов, Д. Б. Моделирование напряжённо-деформированного состояния крупногабаритного трансформируемого рефлектора: Диссерт. ... к.ф.-м.н. —Томск, 2006. — 179 с.
7. Дмитриев, В. И. Численные методы решения задач синтеза излучающих систем [Текст] / В. И. Дмитриев, Н. И. Березина. —М.: Изд-во Моск. ун-та, 1986. — 112 с.
8. Лопатин, А. В. Оценка жёсткости обода большой космической антенны [Текст] / А. В.
A. V. Bel'kov, A. S. Evdokimov, S. V. Ponomarev, D.B. Usmanov
Лопатин, М. А. Рутковская / В сб.: Королевские чтения. Тез. докл. Всерос. молодеж. науч. конф. - Самара: СНЦ РАН, 2003. - Т. 1. — С. 22-23.
Поступила в редакцию 29/У1/2008; в окончательном варианте — 17/Х/2008.
MSC: 74K25, 74K35, 74S05
COMPUTER MODELING OF ADVANCED SPACE REFLECTORS
A. V. Bel'kov1, A.S. Evdokimov1, S. V. Ponomarev1, D.B. Usmanov2
1 Scientific Research Institute of Applied Mathematics and Mechanics by Tomsk State University,
634050, Tomsk, Lenina prosp., 36.
2 JSC "Information Satelline Systems Reshetnev Company", 662972, Zheleznogorsk, Lenin str., 52.
E-mails: khirurgasibmail.com, eas1985amail.ru, [email protected]
Complex technique of computer modeling for truss-tension reflectors, based on the interconnection of tasks of deformed firm body mechanics and electrodynamics, and also computer modeling of intensively-deformed condition of an expandable construction of a space reflector.
Key words: computer modelling, reflector, method of final elements, intense-deformed condition, pattern, gain.
Original article submitted 29/VI/2008; revision submitted 17/X/2008.
Bel'kov Alexey Viktorovich, Younger Scientific Researcher of the Scientific Research Institute of Applied Mathematics and Mechanics by Tomsk State University. Evdokimov Alexandr Semenovich, Younger Scientific Researcher of the Scientific Research Institute of Applied Mathematics and Mechanics by Tomsk State University. Ponomarev Sergey Vasil'evich, Ph. D. (Phis. & Math.), Scientific Research Institute of Applied Mathematics and Mechanics by Tomsk State University.
Usmanov David Bisenovich, Main Expert of the JSC "Information Satelline Systems Reshetnev Company".
