GEOMETRICALY NONLINEAR MODEL OF AN TRANSFORMED RIM OF A BIG SPACE ANTENNA WITH FLEXIBLE COMPOSITE ELEMENTS

Full Text

В настоящее время активно развиваются космические телекоммуникационные услуги широкополосного доступа. В связи с этим актуальна тенденция роста диаметра космических антенн, которые используются для передачи больших объемов данных. У космических антенн есть большой диаметр в операционной конфигурации. Во время вывода на орбиту внутренний объем под обтекателем ракеты-носителя ограничен. Поэтому имеются существенные ограничения, накладываемые на любые большие бортовые антенны космических аппаратов: антенна должна быть сжато упакована, чтобы поместить ее под обтекатель ракеты-носителя, и должна быть развернута на орбите с высокой точностью. Несколько структурных моделей трансформируемых космических антенн могут быть найдены в литературе [1-5]. Одна из таких моделей -антенна, состоящая из обода и приложенной к нему тонкой радиоотражающей мембраны. Концепция такой космической антенны была представлена ранее в работе [5]. Развертывание обода на орбите может быть выполнено несколькими способами. Самым перспективным является трансформация, когда развертывание обода антенны производится запасенной потенциальной энергией деформации его упругих гибких элементов. Применение композиционных материалов в дизайне такой оправы является самым эффективным. Рассмотрена новая модель трансформируемого составного обода большой космической антенны с использованием гибких элементов. При моделировании использована геометрически нелинейная модель, представленная в работе авторов [6]. При этом использованы статические уравнения, описывающие геометрически нелинейную деформацию элемента -гибкой полосы - во время разворачивания обода антенны. В целях удобства проектирования была аналитически решена обратная задача. В этой задаче были исследованы и найдены усилия, смещения и потенциальная энергия, возникающие в гибкой полосе при ее деформации от прямолинейного состояния до сильно криволинейного. Полученное аналитическое решение позволит определять оптимальную форму изогнутой полосы и оценивать зарезервированную энергию. Эти результаты позволяют подойти к проектированию основных упругих элементов обода антенны. Модель трансформируемого обода антенны. Основной элемент трансформируемого обода антенны - криволинейная композитная полоса, которая в ненагруженном состоянии имеет форму, показанную на рис. 1. Рис. 1. Изогнутая полоса Такую изогнутую полосу можно получить формовкой на шаблоне, повторяя ее контур. Две изогнутых полосы формируют основной упругий элемент обода (рис. 2). Развернутый обод антенны состоит из двух кругов изогнутых углепластиковых полос и труб, соединяющих их (рис. 3). Рис. 2. Гибкий упругий элемент Через полосы на внутреннем круге продеты гибкие кабели. Кабели натягивают, прежде чем поместить антенну под обтекатель ракеты-носителя. Точки контакта кривых полос сдвигаются по радиусу к центру антенны (рис. 4). * Исследование выполнено при поддержке Министерства образования и науки РФ, соглашение 14.В37.21.0405 «Моделирование композитных элементов крупногабаритных трансформируемых механических систем космических аппаратов связи и телекоммуникаций». 75 Математика, механика, информатика Рис. 3. Развернутый обод антенны Рис. 4. Деформирование полос при развертывании обода антенны Рис. 5. Свернутый обод антенны Рис. 6. Развернутая антенна с вантовой системой В свернутой конфигурации (рис. 5) кривые полосы выправляются и резервируют энергию вследствие структурной деформации. Развертывание антенны на орбите (рис. 6) осуществляется освобождением зарезервированной энергии. Кабель размотан, и происходит упругое движение элементов по радиусу от центра антенны до внешнего круга (см. рис. 4). Начальными параметрами для проектирования являются внутренний диаметр обода в развернутом состоянии и внешний диаметр обода в свернутой конфигурации. Первый размер определяет диаметр апертуры, второй - диаметр пространства для хранения свернутой антенны под обтекателем ракеты-носителя. Геометрически нелинейная модель гибкой полосы - элемента обода антенны. Рассмотрим задачу определения формы изогнутой полосы. Форма может быть найдена из решения обратной задачи, когда изначально прямоугольная полоса изгибается к необходимому положению (рис. 7). Упругие параметры прямоугольной полосы: B = Ex hb, D — ExI — Exb hL 12 где h - толщина полосы; b - ширина полосы; Ex - упругий модуль материала полосы в направлении оси х. В нашей модели прямоугольную полосу будем рассматривать как гибкий упругий нерастяжимый стержень, консольно закрепленный. Концы стержня могут свободно смещаться вдоль заданных направляющих. Длина стержня L. 1 * \ и \ \ \ 1 R -► ^\0 А* Рис. 7. Изогнутая полоса. Система координат Деформация стержня имеет место в секторе с центральным углом ф. Величина угла определяется как ф = п/n, где n - число упругих элементов, состоящих из двух изогнутых полос (см. рис. 1); R - заданный внутренний радиус развернутого обода антенны. Рассмотрим далее изгиб стержня для заданных величин 76 Вестник Сибирского государственного аэрокосмического университета имени академика М. Ф. Решетнева R, L и ф (см. рис. 7). Найдем форму изогнутого стержня, отклонение конца стержня и запасенную энергию. Выберем декартову систему координат так, что изначально прямой стержень расположен вдоль оси х, защемлен на левом конце в начале координат. На правом конце рассмотрим два варианта его соединения со вторым стержнем (полосой): защемление и шарнирное соединение. Стержень подвергается деформации изгиба под действием неизвестной равнодействующей силы реакции со стороны второго стержня Р, которая сосредоточена на правом конце стержня и из соображений симметрии задачи имеет постоянное направление под углом ф к центру обода. Px и Py - декартовы компоненты силы Р. Обозначим через l - криволинейную координату -длину вдоль изогнутой линии стержня, 0(l) - угол наклона касательной в текущей точке к линии стержня; тогда d0/dl - кривизна кривой линии стержня. Зная зависимость 0(l), можно определить форму изогнутого стержня в параметрическом задании, с параметром криволинейной длины l, выражениями l l x (l ) = J cos 9(l )dl, y (l ) = J sin 9(l )dl. 0 0 Приведем далее систему нелинейных уравнений равновесия стержня, следуя работе [6]. Поскольку отсутствует внешнее давление (распределенная нагрузка), то уравнение равновесия элемента стержня имеет следующий вид: dF ~dl ■ = 0 , F = const, (1) ■ = -Fn = Fx sin 9 - Fy cos 9, dl n x y где M - изгибающий момент. Закон Гука: d9 M = Ddl (2) (3) Из уравнений (1)-(3) имеем одно нелинейное уравнение для определения зависимости 0(l): d 29 D—— + P sin ф sin 9+ P cos ф cos 9= 0 , dl2 или d29 P . ——+— sin (9 + ф) = 0 . dl2 D Преобразуем это уравнение равновесия стержня в уравнение типа нелинейного маятника. Введем безразмерную длину t = l/L, изменяющуюся от 0 до 1, и сделаем замену у = 0 + ф. Введем обозначение для неизвестного параметра q2 = PL2/D. Тогда получим уравнение для угла у типа уравнения нелинейного маятника: d У 2 • —— + q sin у = 0. dt2 Решение этого уравнения хорошо известно: (4) Y(t) = 2arcsin[k sn(qt + Fb k)], dy(t)/dt = 2kq cn(qt + Fb k), 0(t) = -ф + 2arcsin[k sn(qt + Fb k)], (5) где функции sn и cn - эллиптические синус и косинус Якоби. Модуль эллиптических функций k и параметр F, играют роль констант интегрирования, и их связь с параметром q и углом ф определяется далее из граничных условий для каждого конкретного случая изгиба стержня. Введем обозначение для аргумента эллиптических функций u = qt + Fj. Используя выражение (5), получим cos y(t) = 1 - 2k2 sn2 u, sin y(t) = 2k sn u dn u. Интегрируя геометрические соотношения dx/dl = = cos0, dy/dl = sin0, получим координаты произвольной точки стержня в параметрическом виде, определяющие изогнутую форму стержня: t t — = J cos 9dt = J cos (у-ф) = L 0 0 t t = cos ф J cos ydt + sin фJ sin ydt = X0 cos ф + Y0 sin ф, 0 0 t t где F - упругая сила. Пусть Fx = - P sin ф, Fy = P cos ф -ее компоненты в декартовой системе координат, Fn, FT - компоненты в касательной системе координат. Уравнение моментов: dM = cos — = J sin 9dt = J sin (у-ф) = L 0 0 t t ф J sin ydt - sin фJ cos ydt = Y0 cos ф - X0 sin ф. (6) Использованные здесь обозначения X0 и Y0, полученные с помощью (4), имеют вид t t X0 = J cos ydt = J(l - 2k2 sn2 u )t = 0 0 = -1 + — [ Е (am u )- Е (am Fj ), q Y0 =J sin ydt = J 2k sn u dn udt = — [cn FJ - cn u ], (7) 0 0 q где £(am u) - неполный эллиптический интеграл второго рода от эллиптической амплитуды Якоби. Выражения (6) и (7) задают форму профиля изогнутого стержня в параметрическом виде с параметром t - приведенной криволинейной длиной. Определим теперь запасенную потенциальную энергию изогнутого стержня: U = — \M2dl = — J {—1 dt. 2D 2 L I dt i 77 Математика, механика, информатика Введем безразмерный параметр полной потенциальной энергии изогнутого стержня f U —- 3D L f, т. е- f — 6 j (,jt dt. (8) Таким образом, зная решение 0(t), можно точно вычислить запасенную энергию. Рассмотрим случай защемления в соединении двух стержней в элементе обода антенны. Определим связь между параметрами k, F1, q и ф из граничных условий: 1) на левом конце стержня 0(0) = 0; 2) на правом конце стержня 0(1) = п/n = ф. Применим первое условие, используя решение (4) уравнения (5). Получим sn F1 = sin^^/k, следовательно sin ф/ 2' F — F arcsin Соответственно, модуль k изменяется в пределах sin^^2) < k < 1. Применив второе условие, аналогично получим, что sn(q + F1) = sin (ф)/2. (9) Из-за периодичности эллиптического синуса sn уравнение (9) имеет бесконечное количество корней, каждое из них будет определять отдельную ветвь (моду) решений с разным числом перегибов стержня. Исходя из разумного предположения, что изогнутый стержень в реальной задаче будет иметь одну точку перегиба, выберем второй корень уравнения (9), который даст соотношение q = 2K(k) - F2 - F1, F2 — FI arcsin sin ф (10) где K(k) и F(ф, k) - соответственно полный и неполный эллиптические интегралы первого рода. Поскольку правый конец стержня смещается по линии под углом ф, то исходя из геометрии задачи (см. рис. 7), получим дополнительное соотношение У1 (r + x1) — Шф, (11) где у1 и х1 - координаты правого конца изогнутого стержня, а у1 - прогиб стержня. Используя выражения (6) и (7), найдем координаты правого конца стержня: х1 — sin ф71 + cos фХ1, y — sin фХ1 - cos ф^ , 2 Х1 —-1 + —(2E(k)-E2 -E1), q 2k q 1 _sin ф/2 + I - sin (ф) (12) E1 — EI arcsin sin ф/ 2 E2 — E arcsin sin ф k у у 2 ^ I k где E(k) и E^, k) - соответственно полный и неполный эллиптические интегралы второго рода. Далее рассмотрим число упругих элементов n = 18, 24, 32, 48. Внутренний радиус обода антенны зададим r = 2, 3, 4, 5, ..., где введен безразмерный радиус r = R/L. Система алгебраических уравнений (10)-(12) позволяет найти неизвестные параметры решения q и k. Последовательно задавая известные параметры ф = п/n, n и r, подставляя выражение (10) для параметра q(k) в (12), а выражения x1(k) и y1(k) из (12) в геометрическое соотношение (11), получим одно нелинейное алгебраическое уравнение для определения зависимости модуля k от n и r. Это уравнение можно легко решить любым численным методом. Определив модуль k, можно вычислить параметр q, используя выражение (10). Далее, используя параметрические выражения (6), (7) с найденными значениями параметров, можно построить изогнутые формы стержня. Координаты точек стержня в параметрическом задании x(t) и y(t) определяются выражениями (6) и (7) подстановкой значений q и F1 из (10). Таким образом, каждому значению параметров n и r будет соответствовать своя форма изгиба стержня, определяемая одним вычисляемым параметром - модулем k, который находится численным решением алгебраического уравнения (10). Минимальному значению k = sinф/2 соответствует минимальное значение силы P и неизогнутая форма стержня. Когда сила Р ^ да и модуль k ^ 1, мы имеем максимально возможную сильно искривленную форму стержня. Значение прогиба правого конца изогнутого стержня y1(k) можно вычислить с помощью выражений (12), зная параметры решения q и k. Ранее в работе [5] был введен безразмерный параметр z: z — У1 (1+r )гяф В работе [5] были найдены значения этого параметра, а также параметра энергии f по линейной теории. Найденные в [5] значения являются оценкой сверху. Сделав предельный переход k ^ 1 в уравнениях (10)-(12), получим ограничение на область значений параметров n и r: п 1 <- : n r (13) или приближенно nr < n. Это соотношение важно при проектировании обода антенны для выбора разрешенных значений параметров n и r. Последовательные формы изогнутого гибкого элемента обода при r = 2, 3, 4 для случая n = 18, построенные по аналитическим выражениям (6) и (7), приведены на рис. 8, а для случая n = 32 - на рис. 9. Вычисленные значения модуля k путем решения алгебраической системы (10)-(12) приведены в табл. 1, вычисленные по геометрически нелинейной модели с использованием выражения (12) значения безразмерного параметра прогиба z, приведены в табл. 2, значения безразмерного параметра запасенной энергии f одного стержня, вычисленные с помощью выражения (8) для различных n и r, - в табл. 3. В этих таблицах прочерки поставлены там, где наступает ограничение (13). 78 Вестник Сибирского государственного аэрокосмического университета имени академика М. Ф. Решетнева Полная сохраненная энергия обода антенны равна сумме энергий деформации всех кривых полос: ТЛолн = 2nf Большую запасенную энергию обода антенны обеспечивает разбиение на возможно большее число элементов. Таким образом, развита новая модель трансформируемого композитного обода большой космической антенны. В свернутой конфигурации изогнутые элементы обода выправляются и резервируют энергию вследствие структурной деформации. Открытие антенны на орбите осуществляется за счет выпуска зарезервированной энергии. Выполненный анализ с точным учетом геометрической нелинейности позволяет определять форму упругих полос, которые формируют обод антенны на случай любых сильных изгибов элементов. Полученные аналитические решения позволяют оценивать форму и величину прогибов элементов, запасенную энергию кривой полосы, и, впоследствии, вычислять полную сохраненную энергию, требуемую для разворачивания антенны. Такой анализ необходим на ранней стадии дизайна обода антенны. Все аналитические решения легко реализуются в пакете Maple. Таблица 1 Значения модуля k(n, r) n r 2 3 4 5 6 7 8 9 10 18 0,460 0,598 0,755 - - - - - - 24 0,349 0,451 0,559 0,678 0,841 - - - - 32 0,263 0,340 0,418 0,500 0,587 0,686 0,821 - - 48 0,176 0,227 0,278 0,331 0,384 0,439 0,497 0,558 0,623 Таблица 2 Значения безразмерного параметра прогиба z(n, r) n r 2 3 4 5 6 7 8 9 10 18 0,946 0,923 0,890 - - - - - - 24 0,970 0,958 0,945 0,927 0,892 - - - - 32 0,983 0,977 0,970 0,963 0,954 0,942 0,920 - - 48 0,993 0,990 0,987 0,984 0,981 0,978 0,974 0,970 0,965 Таблица 3 Значения параметра запасенной энергииf(n, r) n r 2 3 4 5 6 7 8 9 10 18 0,412 0,889 1,822 - - - - - - 24 0,226 0,462 0,830 1,432 2,861 - - - - 32 0,125 0,250 0,431 0,682 1,039 1,588 2,768 - - 48 0,055 0,108 0,182 0,278 0,400 0,553 0,745 0,991 1,316 о -0.5 -I -1.5 Рис. 8. Последовательные формы изогнутого гибкого элемента обода с защемлением при r = 2, 3, 4 для случая n = 18 y/L Рис. 9. Последовательные формы изогнутого гибкого элемента обода с защемлением при r = 2, 3, 4, 5, 6, 7, 8 для случая n = 32

About the authors

A. V. Lopatin

Email: lopatin@krasmail.ru

Yu. V. Zakharov

K. G. Okhotkin

Email: okg2000@mail.ru

V. V Vilyanen

Email: vilyanen_valya@mail.ru

A. V. Pashkovskiy

Email: linch07@mail.ru

References

Supplementary files