Сопряженное механическое и электродинамическое моделирование трансформируемых космических рефлекторов

Решая данные системы нелинейных неравенств, можно определить диапазон углов атаки и относительных высот движения, в которых крыло обладает статической устойчивостью, а следовательно, обеспечена безопасная эксплуатация аппарата.

Таким образом, решены линейные стационарные задачи безотрывного обтекания для прямоугольного крыла, треугольного крыла с углом поперечного V> 0 и V < 0, получены численные значения аэродинамических коэффициентов с (а, к) и тг (а, к), относительной координаты аэродинамического фокуса X/ (а, к). Построены функциональные зависимости су=/(а, к), тг = /(а, к), X/=/(а, к). Используя критерии Иродова, для оценки устойчивости несущей поверхности вблизи экрана получены многочлены, позволяющие определить диапазон углов атаки и относительных высот движения. Кроме этого, сделаны выводы, что для крыльев треугольной формы с углами поперечного V < 0 или V > 0 наблюдается меньшее смещение аэродинамического фокуса при изменении высоты движения, однако отсутствие шайб на крыльях с углом

поперечного V > 0 приводит к их худшей несущей способности.

Библиографический список

1. Трехмерное отрывное обтекание тел произвольной формы / С. М. Белоцерковский [и др.] ; под ред. С. М. Бе-лоцерковского. М. : ЦАГИ, 2000. 265 с.

2. Численные методы в сингулярных интегральных уравнениях и их применение в аэродинамике, теории упругости, электродинамике / С. М. Белоцерковский, И. К. Лифанов. М. : Наука, 1985. 253 с.

3. Белавин, Н. И. Экранопланы / Н. И. Белавин. СПб. : Судостроение, 1977. 232 с.

4. Жуков, В. И. Особенности аэродинамики, устойчивости и управляемости экраноплана / В. И. Жуков. М. : ЦАГИ, 1997. 81 с.

5. Формалев, В. Ф. Численные методы / В. Ф. Форма-лев, Д. Л. Ревезников. 2-е изд., испр., доп. М. : Физматлит, 2006. 400 с.

M. I. Antipin, I. N. Gusev

THE SERFACE ANALYSIS OF HOWERCRAFT

Bearing surfaces W. I. G.craft the any form are investigated. Values of aerodynamic factors c (a, h) and mz (a, h), relative coordinate of aerodynamic focus x(a, h), distribution of aerodynamic loading on a bearing surface for three forms of bearing surfaces are received numerical. Functional dependences cy = f(a, h), m^ = f(a, h), xf = f(a, h) are constructed.

УДК 539.3:621.396.67

А. С. Евдокимов, С. В. Пономарев, Н. А. Тестоедов, Д. Б. Усманов

СОПРЯЖЕННОЕ МЕХАНИЧЕСКОЕ И ЭЛЕКТРОДИНАМИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТРАНСФОРМИРУЕМЫХ КОСМИЧЕСКИХ РЕФЛЕКТОРОВ

Рассмотрена комплексная методика компьютерного моделирования перспективных рефлекторов, основанная на механике деформированного твердого тела и радиофизике.

Моделирование крупногабаритных трансформируемых рефлекторов является актуальным направлением разработки и создания конструкций систем спутниковой связи и зондирования поверхности Земли, так как экспериментальная отработка подобных конструкций требует больших материальных и временных затрат. Спутник ^игауа с ободной конструкцией рефлектора [ 1] показан на рис. 1.

Основными конструктивными элементами ободных космических рефлекторов являются ферменный обод, обеспечивающий заданный профиль отражающей поверхности и ориентацию рефлектора, сама отражающая поверхность, а также вантовая система. Ферменный обод диаметром 12 м представляет собой стержневую конструкцию, собранную из жестких углепластиковых элементов (рис. 2).

Основные требования к конструкциям рефлекторов заключаются в высокой точности формы отражающей поверхности и наведения, высокой температурной стабильности и радиоотражающей способности антенных систем.

Варианты методик расчета радиотехнических характеристик крупногабаритных рефлекторов рассматривались в работах J. Ruze [1], М. ^ Thomson [2], монографии М. В. Гряника [3], диссертации Д. Б. Усманова [4]. Однако учет искажений отражающей поверхности производился на основе экспериментальных измерений или сильно упрощенных допущений о деформациях отражающей поверхности.

Целью комплексного моделирования является повышение точности за счет использования для расчета радиотехнических характеристик равновесной формы от-

1О5

ражающеи поверхности, полученной в результате компьютерного моделирования напряженно-деформированного состояния (НДС) рефлектора с позиций механики деформируемого твердого тела при отработочных нагрузках на Земле и функционировании в космосе.

Рис. 1. Спутник ТЬитауа

Фронтальная

сеть

СетепоХотно

Растяжки

Ферменный

барабан

Эм, / Эх, Эм2/ Эл, Эмз/ Эл:,

0

0

Эм, / Эх. 0

Эм,/ Эх3

В =

0

0

Эм./Эх2 0

Эм./ Эх3

Э/Эх,

0

0

0

0

дм-з/ Эх. 0

Эм. / Эх3

0

дм,/ дх2 0

д м,/ Эх, дм,/ Эх3 0

0

Эм. /дх2 0

дм 2/ дх, дм2 / Эх3 0

0

Эм. /Эх. 0

Эм. / дх, Эм. /Эх3 0

0

0

дм,/ Эх3 0

дм,/ дх 2 дм,/ дх,

0 0

0 0

дм2/ дхз Эм./ Эх.

0 0

Эм./Эх2 Эм./Эх.

Эм./ Эх, дм./Эх,

0 0 д/дх2 0 д/дх3

д/дх2 0 д/дх, д/дх3 0

0 д/дх3 0 д/дх2 д/дх,

Б Е V

'1 -V V V 0 0 0

V 1 -V V 0 0 0

V V 1 -V 0 0 0

0 0 0 (1 - .V)/. 0 0

0 0 0 0 (1 - .V)/. 0

0 0 0 0 0 (1 - .V)/.

Е и п - модуль упругости и коэффициент Пуассона материала конечного элемента соответственно.

Система уравнений равновесия имеет вид

ВТ у + р = с и где р - вектор объемной силы.

Добавление начальных и граничных условий делает постановку задачи с позиций механики деформируемого твердого тела полной. В результате такой постановки получается нелинейная задача.

Основой для отыскания деформированной формы отражающей поверхности и остальной части рефлектора является принцип возможных перемещений для упругодинамических задач [5]:

|[с йтщ + уде - рдц^У = 0,

ТыХьная сеть

Рис. 2. Основные компоненты конструкции рефлектора

Постановка задачи о НДС рефлектора. Пусть изменение положения точки элемента конструкции при деформировании обозначается вектором перемещений й(І) = [й1, й2, й3]т. Связь между тензором деформаций и вектором перемещений представима в виде е = Вй + Aq /2, где

q = [Эм1/Эх1 ди2/дх1 ди3/дх1 ди1/дх2 ди2/дх2 ди3/дх2 ди1/дх3 ди2/дх3 ди3/дх3]т;

е = [ е11е22е33е12е23е31]Т-

Тензор напряжений имеет компоненты о = [о11 о22 о33 о12 о23 о31 ]т. Физические соотношения между напряжениями и деформациями, учитывающие температурные деформации и начальные напряжения, имеют следующий вид: о = Б(е - ет) + о0, где ет = [а1АТ а2Ат а3АТ 0 0 0]т - вектор температурных деформаций; s0 - вектор начальных напряжений; Б - матрица материальных констант:

где 5 - вариации соответствующих величин; р - вектор массовых сил; и - вектор ускорения.

Для получения численного решения задачи о напряженно-деформированном состоянии рефлектора использовался метод конечных элементов [6]. Пространственная геометрическая модель рефлектора разбивалась на конечные элементы, в которых аппроксимация перемещений и выбиралась как линейная зависимость от узловых перемещений и = Ти, где ¥ - матрица, образованная базисными функциями; и - вектор узловых перемещений конечного элемента.

Тогда выражение деформаций через перемещения и напряжений через деформации в векторном виде примут вид

е = В0и + Л0/2,

5е = (В0 + ЛО)5П ,

0 = Ои, о = Б(е - еТ) + о0,

Во = В¥ .

Граничные условия учитывались заданием нулевых перемещений узлов конечных элементов в местах крепления рефлектора к штанге, связывающего его с космическим аппаратом. Влияние температурных нагрузок на напряженно-деформированное состояние конструкции моделировалось через зависимость напряжений от деформаций внутри каждого конечного элемента. Решение подобной нелинейной задачи позволяет получить равновесную форму отражающей поверхности и напряженно-деформированное состояние при заданных геометрических параметрах рефлектора, физико-механических характеристиках материалов элементов конструкции, начально-краевых условиях и нагрузках.

х

Для моделирования напряженно-деформированного состояния на поверхности Земли рассмотрены три положения рефлектора: «чаша вниз», «чаша вверх», «вертикальная вывеска». Граничные условия соответствовали полному закреплению ободной конструкции в узлах связи со штангой от космического аппарата. Объемная нагрузка соответствовала ускорению свободного падения у поверхности Земли. В плоскости раскрыва рефлектора задавались начальные напряжения о0 отражающей поверхности, соответствующие рабочим натяжениям сетепо-лотна. В качестве обобщенной меры отклонения отражающей поверхности рефлектора в равновесном состоянии использовалось среднеквадратичное значение отклонений (СКО) полученной расчетной поверхности в узлах конечноэлементной сетки от поверхности соответствующего параболоида.

При расчете форм и частот собственных колебаний было рассмотрено закрепление рефлектора в точках соединения со штангой от космического аппарата. Такие граничные условия соответствуют консольному закреплению ободной фермы.

Решение механических задач проводилось методом конечных элементов, реализованным в программном комплексе ANSYS. Средствами ANSYS была построена геометрическая и конечно-элементная модели космического рефлектора с фермой натяжения. Методика решения нелинейной задачи основывалась на ранее разработанной методике для зонтичных рефлекторов [7].

Электродинамическое моделирование. Определение основных радиотехнических характеристик антенн связано с получением выражения для электромагнитного поля в дальней зоне, когда источниками поля являются заданные сторонние токи j на отражающей поверхности рефлектора.

Система уравнений Максвелла имеет вид [8] rot H = -1шЕ + j,

rot Е = /юцН.

С помощью соотношений H = rot А,

Е = /Ю|х4 +—grad divA Ю£

вводится векторный потенциал A. Уравнения Максвелла для векторного потенциала в декартовой системе координат дают векторное уравнение Гельмгольца:

ДА + к2 A = -j,

где k = ю(е • ц)1/2 - волновое число.

В каждой точке поверхности зеркала рефлектора возникает поверхностный ток, фаза, амплитуда и направление которого определяется соотношением

js = 2 [«оH ]

где j - вектор плотности поверхностного тока в данной точке зеркала; H - вектор напряженности первичного магнитного поля облучателя в этой точке; п0 - орт нормали к поверхности зеркала в этой же точке.

Соотношение (1) является точным только в случае падения плоской волны на плоскую бесконечно большую поверхность. Напряженность первичного магнитного поля облучателя определяется формулой

г

где С - коэффициент, не зависящий от г и характеризующий направленные свойства облучателя; г - расстояние от фазового центра облучателя до точки, в которой определяется поле. Для построения картины токов, возникающих на отражающей поверхности рефлектора под влиянием поля облучателя, необходимо знать распределение вектора С в пространстве, т. е. векторную диаграмму направленности облучателя. Распределение тока на зеркале определяется по формуле (1), которая в скалярной форме приобретает вид

}а = 2(пуИг - пгИу),

4 = 2(пИ. - пИ,),

/„ = 2(пхИу - пуИх).

Зная распределение тока на поверхности зеркала, можно определить направленные свойства параболической антенны. Для этого необходимо проинтегрировать по всей поверхности зеркала выражение для напряженности поля, создаваемого элементом поверхности зеркала, рассматривая его как элементарный электрический вибратор.

Поле излучения параболоида можно представить в виде

£ = -е-т [у,е-Ж 2гХ * !х

Е =-12ь ^ Ух е~*г [ ] е~*Ж,

нер 2гХ [ Лу

где Б - поверхность параболоида, используемая в качестве антенны; Е - напряженность поля, созданного токами основной поляризации; Епер - напряженность поля, созданного токами перекрестной поляризации. Угол q - определяется относительно луча, идущего прямо от облучателя до точки приема.

В качестве численного метода использовалось сочетание метода моментов и метода физической оптики. Отражающая поверхность разбивалась на элементарные площадки треугольной формы. Наибольшее количество узлов в некоторых расчетах достигало величины 1,5 • 106из соображений, что длина элемента не должна превышать 10-1 длины волны. Для расчета радиотехнических характеристик применялся пакет FEKO [9] - система 3D электромагнитного моделирования.

Результатыы механического анализа. Для положения рефлектора «чаша вверх» результаты показаны на рис. 3 и 4.

Распределение перемещений в направлении, перпендикулярном плоскости раскрыва, показывает, что отражающая поверхность рефлектора в положениях чашей «вверх» имеет СКО меньше, чем чашей «вниз». Полученные результаты согласуются с результатами [4].

Характер распределения интенсивности напряжения отражающей поверхности, позволяющий говорить о достаточной равномерности натяжения сетеполотна и, соответственно, возможных незначительных отклонениях величины коэффициента отражения радиоволн от оптимального значения, показан на рис. 4.

Динамические расчеты используются для определения собственных частот и форм колебаний, которые являются важными характеристиками, учитываемыми при проектировании конструкции в целях учета в условиях динамического силового воздействия (рис. 5).

•кем. ЗОІАГТІОИ

Рис. 3. Распределение перемещений (м) по оси 2 при воздействии силы тяжести

МООМ. ІС4ЛГПШ

4М» 0848 МШ М20І М1?1

---------ц;и---------------ШИ----------------ШІІ--------------->028?-------------инь

Рис. 4. Интенсивность напряжения отражающей поверхности

МОСЛЬ юипісм

НЛЇ 28 2006

.008)14 -.008И4 -.0018*7 .001

__________-.«««_________________.КПП____________________.ТОО»-От______________________________________.00»»_.РРИИ

а

МСФАІ. ІОЮТІСЖ

-.011*78 .00*717 .0)1*)) .0И»|7

-.гит -.ті їм .кпим .імім* .пит

б

Рис. 5. Вид перемещений (м): а - для 1-ой собственной частоты; б - для 2-ой собственной частоты

По результатам расчетов построен график зависимости собственных частот рефлектора от натяжения фронтальной сети (рис. 6). Полученные результаты по значениям частот собственных колебаний хорошо согласуются с расчетными данными рефлектора AstroMesh [10].

Рис. 6. График частот собственных колебаний рефлектора от номера собственной формы

Результаты электродинамического анализа. Одной из основных характеристик антенн является диаграмма направленности - зависимость излучаемого поля от положения точки наблюдения. Трехмерный вид основного лепестка диаграммы направленности (ДН), полученный расчетным методом, представлен на рис. 7.

Рис. 7. Геометрия модели параболического рефлектора

Сечения остронаправленных ДН удобнее и точнее изображать в прямоугольной системе координат, поскольку угловой масштаб здесь может быть выбран произвольно в соответствии с шириной ДН.

Диаграммы направленности для идеального параболоида и отражающей поверхности с расчетной равновесной формой для ободного рефлектора показаны на рис. 8. Граничные условия соответствовали закреплению рефлектора в точках соединения с штангой от космического аппарата. Нагружение конструкции производилось температурным полем, которое привело к соответствующему деформированному состоянию отражающей поверхности рефлектора и изменению положения оси рефлектора. Вследствие чего, во-первых, уменьшился главный лепесток, во-вторых, смещена вся диаграмма и, в-третьих, имеются изменения боковых лепестков. Результаты качественно соответствуют результатам, опубликованным в работах М. В. Гряника, В. И. Ломона [3] и других авторов. При этом наибольшие отклонения отражающей поверхности от идеального параболоида достигали 2 мм.

Диаграммы направленности деформированной отражающей поверхности рефлектора при различном СКО показаны на рис. 9. При расчете радиотехнических харак-

теристик в качестве облучателя использовалась модель рупора.

r \ \

T"T \ \

T \ \ \

\ \ L.\

D.4 1 f \ \

/7 1 \ \

i 1 \\

t1 Jj \ \ \ b

тгаяяштттт

Рис. 8. Диаграммы направленности в вертикальной

плоскости:------для идеальной параболической поверхности;

------расчетной равновесной формы отражающей

поверхности рефлектора

Рис. 9. Диаграммы направленности рефлектора

Последние результаты демонстрируют, что резко изменяется уровень бокового излучения, вследствие чего существенно уменьшается коэффициент направленного действия.

В полученных результатах ширина диаграммы направленности антенны для идеальной параболической поверхности равна 1,8° при частоте 1 ГГц и 0° 10,5’ при частоте 10,4 ГГц.

На основе подходов механики деформируемого твердого тела и электродинамики реализована комплексная методика компьютерного моделирования перспективных трансформируемых космических рефлекторов, позволяющая более точно учитывать форму и напряженность отражающей сетчатой поверхности, сократить объем экспериментальных работ при создании оптимальных конструкций ободных рефлекторов по заданным ДН и прогнозировать эффективность функционирования трансформируемых космических антенн в условиях космического пространства.

Библиографический список

1. Ruze, J. Antenna Tolerance Theory - A Review / J. Ruze // Proceedings of the IEEE. 1966. Vol. 54 (April). Р. 633-640.

2. Thomson, M. W. Astromesh deployable reflectors for Ku- and Ka-band commercial satellites / M. W. Thomson. AIAA. 2002. Р. 2032.

3. Гряник, М. В. Развертываемые зеркальные антенны зонтичного типа / М. В. Гряник, В. И. Ломан. М. : Радио и связь, 1987. 72 с.

4. Усманов, Д. Б. Моделирование напряженно-деформированного состояния крупногабаритного трансформируемого рефлектора : дис. ... канд. физ.-мат. наук / Д. Б. Усманов. Томск, 2006. 179 с.

5. Сахаров, А. С. Метод конечных элементов в механике твердых тел / А. С. Сахаров, И. Альтенбах. Киев : Вища школа. Головное изд-во, 1982. 480 с.

6. Зенкевич, О. Метод конечных элементов в технике / О. Зенкевич ; под ред. Б. Е. Победри. М. : Мир, 1975. 541 с.

7. Бутов, В. Г. Моделирование температурных деформаций рефлекторов космических аппаратов / В. Г. Бутов [и др.] // Физика. 2004. № 10. Приложение. С. 10-18.

8. Дмитриев, В. И. Численные методы решения задач синтеза излучающих систем / В. И. Дмитриев, Н. И. Березина. М. : Изд-во Моск. ун-та, 1986. 112 с.

9. FEKO User’s Manual Suite 4.2., 2004.

10. TRW-built AstroMesh reflector deployed aboard Thuraya spacecraft [Электронный ресурс] / TRW, Inc. Электрон. дан. Режим доступа: http://www.trw.com. Загл. с экрана.

A. S. Evdokimov, S. V. Ponomarev, N. A. Testoedov, D. B. Usmanov

THE COMPLEX MECHANICAL AND ELECTRODYNAMIC MODELING OF TRANSFORMED SPACE REFLECTORS

The complex technique of computer modeling of the perspective reflectors based on mechanics of the deformed firm body and electrodynamics is considered.