Геометрическая модель раскрывающейся крупногабаритной космической конструкции ферменного типа Текст научной статьи по специальности «Физика»

УДК 629.783:514.85

В. Е. Мешковский

ГЕОМЕТРИЧЕСКАЯ МОДЕЛЬ РАСКРЫВАЮЩЕЙСЯ КРУПНОГАБАРИТНОЙ КОСМИЧЕСКОЙ КОНСТРУКЦИИ ФЕРМЕННОГО ТИПА

Рассмотрена геометрическая модель, содержащая информацию о форме и размерах раскрывающейся крупногабаритной космической конструкции ферменного типа, ее элементах и их взаимном расположении. В соответствующую математическую модель входят соотношения, определяющие координаты центров узловых шарниров лицевого и тыльного поясов, положения геометрических осей цилиндрических шарниров складывающихся и диагональных стержней. Предложенная модель может использоваться при проектировании и разработке математических моделей для исследования собственных динамических характеристик конструкций рассматриваемого типа и при исследовании процесса их раскрытия. E-mail: sm11@sm.bmstu.ru

Ключевые слова: крупногабаритные космические конструкции, ферменные конструкции, раскрывающиеся антенны, геометрическая модель.

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

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

Складывающиеся стержни тыльного пояса

х

Диагональные стержни

Рис. 1. Каркас рефлектора антенны

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

На конструкцию космической антенны накладывается ряд ограничений, в частности ограничения на массовые и габаритные характеристики. Конструкция антенны должна: иметь минимальную массу, ограниченные габариты в сложенном состоянии для возможности ее размещения под обтекателем ракеты-носителя и большую площадь рабочей поверхности в раскрытом положении; быть достаточно жесткой в рабочем состоянии и обеспечивать надежность раскрытия из транспортного состояния в рабочее.

Разработка конструкции антенны, удовлетворяющей указанным требованиям, как и разработка любого технического объекта, основывается на математическом моделировании [1].

На начальном этапе использована упрощенная геометрическая модель, предполагающая пошаговое определение координат точек, соответствующих центрам узловых шарниров лицевого и тыльного поясов (далее центры узловых шарниров), основанная на методике расчета, предложенной ОКБ МЭИ — разработчиком космических антенн ферменного типа.

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

Рис. 2. Отображение лицевого пояса рефлектора на плоскость Оиу

Т = {М е М | Мр, р е N Л [1, ^уз.л]}, где М, N - множества действи-

тельных

и

натуральных

чисел,

N

уз.л

= 2N0 (2Ni + 1) — Ni (Ni — 1) + 1 — число узловых шарниров лицевого пояса.

Если в плоскости Ouv (рис. 2) определить множество D =

= {(u, v) Е Z2 | (и > —N0 A u ^ N0 A v > N1 A v ^ N1 A v ^ —u + +N0 A v > —u — N0)}, где Z — множество целых чисел, то F = S(D), где S : D ^ R.

Представим F в виде F = (у F(s)) U ( U Fjr)), (i, j Е N A [1,6]),

где F(s) = S(D(s)), Fj(r) = S(Djr)),

D(s) = D П {(u, v) e Z2 | x(+2u + xk+iv > 0 A u + xk-iv > 0}

d+2u + xk+iv > 0 A xk1)u + ^'(1)

— множество внутренних точек i-го сектора с углом 60о,

j =D П {(u, v) e Z2 | |xj-i|u+|xj+i|v=0 A xj+iu+j-j)v > 0}

(i) 4

(1)

(1)

(1)

(2)

— множество точек, принадлежащих ] -му лучу, представляющему собой границу двух смежных секторов,

хП1} = sign(sin(nп/3)), хП2) = sign(cos(nп/3)).

Определим множество Тк = Т^ и Тг) и Т^, где к е N Л [1, 6], к0 = к — 1, если к > 2, и к0 = 6, если к = 1, которое будем назы-

а

вать множеством точек k-го сектора лицевого пояса. Форма и размер лицевого пояса характеризуется числом пролетов длиной L каркаса рефлектора: N0 — вдоль оси Ou и N1 — вдоль оси Ov.

На первом шаге расчета определяем координаты центров узловых шарниров на основной кривой, являющейся образующей для поверхности вращения и расположенной в координатной плоскости Oxz (плоскость Q) глобальной системы координат Oxyz с началом в вершине поверхности вращения (точка O) лицевого пояса (см. рис. 1). Длины всех складывающихся стержней, расположенных вдоль этой кривой, равны L, а координаты y всех рассматриваемых на этом шаге центров равны нулю. Для центра в начале координат x = y = z = 0. Далее последовательно находим координаты центров x!Q), ziQ), n Е N Л [1,N0], из системы нелинейных алгебраических уравнений

(xn xn-1) + (zn zn-1) L = 0

Ф^) -znQ) = 0,

где x0Q) = z0Q) = 0; Ф(x) — функция кривой, являющейся образующей поверхности вращения.

Первое уравнение определяет расстояние между центрами

í (Q) (Q)) í (Q) (Q) ) т

[xn ,zn) и (xn+1 ,zn+1), равное L, второе — расположение центра на основной кривой.

В частных случаях функция Ф^) имеет следующий вид: Ф^) = Sx2

для параболоида вращения, где S = — (f — фокусное расстояние), и

Ф^) = 0 для пластины.

На втором шаге, используя оператор R^=n/3, определяем координаты центров узловых шарниров, лежащих на кривой в плоскости P, проходящей через ось симметрии поверхности вращения и расположенной под углом 60° к плоскости Q, в которой лежит основная кривая (см. рис. 2):

xnP) = xnQ) cosф, y(P = xnQ) sinф, znP) = znQ).

Под оператором R будем понимать отображение Y = R (X) с помощью линейного ортогонального оператора с матрицей преобразования

(cos ф — sin ф 0 sin ф cos ф 0 0 0 1

соответствующей вращению относительно оси симметрии (ось Oz) поверхности лицевого пояса на угол ф.

На третьем шаге определяем координаты центров, лежащих внутри первого сектора (к = 1). Рассматриваем ячейку М1М2М3М4 (см. рис. 2). При известных координатах центров М2(г — 1, ]) и М4(г, ] — 1) находим координаты центра М3(г, ]) в глобальной системе координат Ожуг из следующей системы нелинейных алгебраических уравнений:

(xi,j - Хм)2 + (yi,j - Ум)2 + (Zi,j - ZM)2 - h2 = 0,

где Хм = 0,5(xi,j-i + Xi-i,j); ум = 0,5(ум-1 + yi_i,j); zm = 0,5(zi,j-i + + Zi— 1 ,j ); A Xi,j — 1 Xi— 1 ,j ; B yi,j-1 yi— 1 ,j ; C Zi,j— 1 Zi— 1,j ; D = -(AXm + ByM + CZm); h2 = L2 - 0,25(A2 + B2 + C2).

В системе (1) первое уравнение представляет собой уравнение плоскости, проходящей через середину отрезка М2М4 перпендикулярно ему, т.е. нормальный вектор N этой плоскости имеет координаты А, В, С. Второе уравнение системы (1) определяет принадлежность точки М3 поверхности вращения Ф(х2 + у2), а третье — расстояние от точки М3 до точки М (здесь и далее полужирным шрифтом обозначаются векторные величины).

Таким образом, рассматриваемая система (1) получена на основе одного из основных условий (равенство расстояний от точек М2 и М4 до точки М3 (см. рис. 2)), заложенных в построение рефлектора и необходимых для приведения каркаса рефлектора в транспортное плотноупакованное состояние и обеспечения его раскрытия на орбите.

На заключительном шаге с использованием операторов Яфф=кп/3, к € [2, 6], находим координаты центров узловых шарниров, расположенных в остальных пяти секторах, т.е. множества Тк =

= Яф=кп/3(Т1).

После того как определены координаты центров узловых шарниров лицевого пояса, ведется расчет координат центров узловых шарниров тыльного пояса (множество В).

Так как длины Ь1 всех диагональных стержней равны друг другу, то координаты центров тыльного пояса находятся из условия их равноудаленности от соответствующих трех центров лицевого пояса (рис. 3). В этом случае центр Q лежит на прямой перпендику-

лярной плоскости треугольника М1М2М3 и проходящей через центр 01 описанной окружности. Исходя из условий построения лицевого пояса длины отрезков М1М3 и М2М3 равны Ь.

Единичный вектор нормали к плоскости треугольника М1М2М3 определяется так:

AXi,j + Byij + CZi,j + D = 0, $(x2,j + y'j) - Zi,j = 0,

(1)

n0 = cos a i + cos в j + cos y k,

(2)

Рис. 3. Элементарная ячейка (тетраэдр) рефлектора

A

B

C

где cos а = ±—; cos в = ; cos 7 = ±—; N = у/A 2 + B2 + C2;

A = (У2 - У1 )(Z3 - Zi) - (уз - yi)(Z2 - Zi); B = (X3 - Xi)(z2 - zi) --(X2 - Xi)(Z3 - Zi); C = (X2 - Xi)(уз - yi) - (хз - Xi)(y2 - yi); i, j, k — орты осей системы координат Oxyz; х, y^, z^, i = 1, 2, 3, — координаты точек Mi ,M2 ,M3.

Знак в выражении (2) для направляющих косинусов выбирается из условия расположения центров узловых шарниров тыльного пояса с внешней стороны образующей поверхности, т.е. должно выполняться условие

(no ) < 0,

где (...,____) означает скалярное произведение векторов.

В этом случае координаты узла Q тыльного пояса определяются по зависимостям (см. рис. 3)

xq = хо1 + Hi cos а, yQ = уо1 + Hi cos в, zq = zo1 + Hi cos 7,

где XOl, yOl, zOl — координаты точки Oi; Hi=y/L2 - R2; R = 0,5L2/H; H = V/(X3 - Xp)2 + (уз - yp)2 + (z3 - zp)2; xp = 0,5(xi + X2); yp = = 0,5(yi + У2); zp = 0,5(zi + z2).

В данном случае треугольник Mi M2 M3 — равнобедренный, что определяет соответствующие выражения для высоты H и радиуса R описанной окружности.

Координаты точки Oi находятся по зависимостям

хз + Axp Уз + Ayp z3 + Azp

Хо1 = "TTX", Уо1 = т+т, zoi = "ГТА",

где A = R/(H - R).

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

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

В полярной целочисленной системе координат каждому центру ставится в соответствие пара чисел (m, n), где m £ N определяет значение полярного угла, а n £ N — полярного радиуса. В качестве примера на рис. 4 показан закон изменения m и n для N0 = 6, Ni = 3. Так, для n = 1 число m принимает значения 1, 2, ... 6. Аналогично для n = 2 и 3. Для n = 6 число m может принимать значения от 1 до 36. Однако пары чисел (m, n), соответствующие областям off, не принадлежат множеству D. Поэтому для n = 6 значения m составляют от 1 до 4, от 16 до 22, от 34 до 36. Аналогично и для n = 4 и 5. Здесь необходимо отметить следующее. Если в плане лицевой пояс представляет собой правильный шестиугольник, т.е. Ni = N0, то область off представляет собой пустое множество, т.е. в этом случае все пары чисел (m, n) (см. рис. 4) принадлежат множеству D.

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

v „

^o = arctg-—— + C 2u + v

где C = 0, если 2u + v> 0 Л v> 0; С = п, если 2u + v < 0; С = 2п,

Рис. 4. Полярная целочисленная система координат

3

если 2u + v > 0 Л v < 0; С = п/2, если 2u + v = 0 Л v > 0; С = - п

2

если 2u + v = 0 Л v< 0.

Далее определяем номер сектора

k = E ^о) + 1

где E(...) — целая часть числа, и соответствующие полярные координаты:

|u| + |v| при k = 1, 4, n = |v| при k = 2, 5,

|u| при k = 3, 6;

m = E

3n

n

+1

(3)

Для найденных значений п и т находим номер центра:

_ Г 3п(п — 1) + т, п ^ N1,

Пц = \ 3Жх(N1 + 1) + 2(2Жх + 1)(п — N1 — 1) + по, п > N1,

где п0 = т, если к = 1; п0 = т — (3п — 2^ — 1), если к = 3, 4; п0 = т — 2(3п — 2Nl — 1), если к = 6.

Для к = 2, 5 номер пц определяется по первой зависимости в выражениях (3), так как для этих секторов п ^ N. Приведенные зависимости учитывают, что значения целочисленных полярных координат

могут изменяться в следующих интервалах:

n G [1,N0], m G [1, 6n] при n < Nb

m G [1, N1 + 1], если k = 1 Л n > Ni; m G [3n - N1 + 1, 3n + N1 + 1], если k = 3, 4 Л n > N1; m G [6n — N1 + 1, 6n], если k = 6 Л n > N1.

Данные интервалы изменения полярных координат соответствуют лицевому поясу, контур которого в системе координат Ouv подобен контуру, изображенному на рис. 2.

Алгоритм вычисления прямоугольных декартовых координат x, y, z точки M некоторого узла лицевого пояса, лежащей в секторе k и имеющей косоугольные координаты u и v, состоит из следующих шагов.

Шаг 1 — определяются полярные координаты r, < (нецелочисленные):

r = у/£2 + П2, < = arctg | + C1,

где £ = u + v cos a; n = v sin а, а = п/3, C1 = 0, если £ > 0 Л п > 0; п,

3

если £ < 0; 2п, если £ > 0 Л п < 0; п/2, если £ = 0 Л п > 0; - п, если

2

£ = 0 Л п < 0.

Шаг 2 — определяются косоугольные координаты точки M(1), лежащей в первом секторе и являющейся прообразом точки M, т.е. M = R,=W3(M(1)):

V(1) = Eo (¡Ira) , U(1) = Eo(£(1) — v(1) cos^

где £(1) = rcos<p(1); n(1) = rsin<(1); <(1) = < — (k — 1)a, Eo(...) — функция округления числа.

Шаг 3 — определяются декартовы координаты x, y, z точки M:

X = r(1) cos <(k), y = r(1) sin <(fc), z = Z(1),

где r(1) = \/x21) + У21)» <(k) = arctg~ +(k — 1)a; x(1), У(1), z(1) —

x(1)

прямоугольные декартовые координаты точки M(1).

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

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

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

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

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

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

Вводятся четыре системы координат:

— глобальная система Оху^ с базисом {е} = (г ] к) (см. рис. 1);

— система Охрургр (рис. 6) с базисом {ер} = (гр ] кр), координатная плоскость Охргр которой проходит через ось О^ и центр узлового шарнира Мр;

— система Мр£пС с базисом {д } = (тр Ьр пр), центр которой (точка Мр) совпадает с центром узлового шарнира, координатная плоскость Мрсовпадает с плоскостью корпуса узлового шарнира, а ось Мр£ направлена по касательной к кривой гр = f (хр), являющейся линией пересечения координатной плоскости Охргр с образующей поверхностью (см. рис. 6);

— система Мр£дщСд (рис.7) с базисом ^} = (трд Ьрд щд), координатная плоскость Мр£дпд которой совпадает с плоскостью корпуса узлового шарнира, а ось Мр£д проходит через середину (точка Ард) геометрической оси.

В данном случае р Е N Л [1, ^з.л], а д Е N Л [1, 6]. Далее связь между указанными базисами определяется с помощью матриц направляющих косинусов {«ро)}, {«р1»}, {ард»}:

{ер}т = {аро)}{е}т, {др}т = {а^}{ер}т, {дрд}т = {а^}{др}т.

Для углов , 0р, врд = + — (д — 1) на рис. 6,7 матрицы {«р0»},

3

{ар1»}, {«рд»} имеют следующий вид:

rj

Рис. 7. Положение осей цилиндрических шарниров складывающихся стержней в узловом шарнире

К0)} =

cos sin 0 — sin cos 0

, {«¡1}

0

0

1

cos #p 0 sin #p

— sin 0 cos #p

{apq }

cos epq sin ß,

0

— sin epq cos epq 0

pq

0

0

1

где tg ^р = Ур/хр.

Поскольку вектор тр в точке Мр направлен по касательной к кривой

% = /(хр), то ^р = бг^^

dx p

Зная угол ^>p, можно определить радиус-векторы точек A

■pq

r Apq = r1 T

pq

— h1 n

pq

(4)

и направляющие векторы Ьр? осей цилиндрических шарниров, проходящих через точки (см. рис. 7). В формуле (4) п и Л4 соответствуют координатам точки в системе координат Мрп<?С? (рис. 8).

После определения положения осей вращения складывающихся стержней (определено положение систем осей МрС?) могут быть найдены аналогичные геометрические параметры для осей цилиндрических шарниров трех диагональных стержней, сходящихся в данном узловом шарнире (см. рис. 8).

Из конструкции узловых шарниров следует, что если направление координатной оси Мр^ для д =1 определено, то направляющие век-

0

1

0

Рис. 8. Положение осей цилиндрических шарниров диагональных стержней в узловом шарнире

торы осей вращения диагональных стержней

Si = cos YiTp1 + sin 7¿bp1, i G N Л [1, 3],

2п 2п

где Y1 = —3-, 72 = 0 Y3 = у.

Координаты точек, соответствующих серединам осей вращения диагональных стержней, определяются векторами (см. рис. 8)

Гвг = r2Si — h2«p1. (5)

В формуле (5) r2 и h2 сооответствуют координатам точки B в системе координат Мр£р1пр1Ср1 (см. рис. 8). Положение геометрических осей диагональных стержней показано на рис. 7 штриховыми линиями.

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

Поверхность тыльного пояса задается дискретным множеством точек B = {Q G R | Qr, r G N Л [1, NУз.т]}, где N^ = N1(4No — N1) — число узловых шарниров тыльного пояса. Множество B является отображением множества F, т.е. T : F ^ B, и оно основано на равенстве расстояний точки Qr G B тыльного пояса от соответствующих трех точек лицевого поса c координатами Mpi (u, v), Mp2 (u, v — 1), Mp3 (u + + 1,v) G F. Эти четыре точки являются вершинами тетраэдра. В рассмотренном выше алгоритме определения координат центров узловых шарниров тыльного пояса предполагалось, что предварительно необходимо полностью определить множество F.

Однако множество B можно определить и другим способом. Как указывалось выше, точки лицевого пояса находят пошаговым расчетом: сначала определяют точки первого сектора F1, а затем с помощью оператора R — точки остальных пяти секторов: Fk = R^=kn/3(F1), где k G [2, 6]. После того как определены точки лицевого пояса в первом секторе F1, используя отображение T, определяют соответствующие точки тыльного пояса B1 = T(F1), а затем посредством оператора находят и остальные точки как лицевого, так и тыльного поясов. Из этого алгоритма следует, что точки тыльного пояса, как и точки лицевого пояса, расположены на поверхности вращения, не имеющей какого-либо аналитического представления в отличие от поверхности, на которой расположены центры узловых шарниров лицевого пояса (для антенны — это параболоид вращения).

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

принадлежащих кривой, являющейся образующей для поверхности вращения, которая задается рядом дискретных точек. С поверхностью связывается некоторая цилиндрическая система координат. Затем все точки поверхности тыльного пояса отображаются на полуплоскость, проходящую через ось аппликат (ось вращения поверхности) Ог и полярную ось Ог. В полученном множестве проводится упорядочение точек с целью исключения повторяющихся, т.е. точек, имеющих одинаковые аппликаты и полярные радиусы, и расположения их в порядке возрастания полярных радиусов. В результате получаем дискретное множество точек М = {(г, г) € М2 | гк+1 > гк > 0, к € N Л [1, Км]}, определяющих положение центров узловых шарниров тыльного пояса для положительных значений полярного радиуса. Множество М позволяет численно определить значения первых производных с использованием В-сплайнов [2].

Чтобы уменьшить погрешности определения значений первых производных на концах отрезка аппроксимации, размер этого отрезка выбираем следующим образом. В цилиндрической системе координат полагаем г > 0. Однако далее, рассматривая систему координат Огг, будем считать, что г может принимать и отрицательные значения. Тогда можно зеркально отобразить множества М относительно оси Ог, т.е. в этом случае дискретный ряд точек будет задан на отрезке [—гКм, +гКм ]. Это позволяет более точно определить значения производных в точках, расположенных около начала координат. Далее для повышения точности определения значений первых производных в точках, соответствующих центрам узловых шарниров, расположенных на периферии конструкции рефлектора, вычисление множеств точек Т и В, соответствующих центрам узловых шарниров лицевого и тыльного поясов, проводится для увеличенного значения N0). Результаты вычисления показали, что достаточно увеличить это значение на единицу.

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

Для определения указанных углов можно поступить следующим образом. Положение осей цилиндрических шарниров в корпусах узловых шарнирах фиксировано, т.е. они располагаются под углом 60° по отношению друг к другу (см. рис. 7). Для обеспечения раскрытия рефлектора оси цилиндрических шарниров двух смежных узловых шарниров, т.е. шарниров, связанных одним складывающимся стержнем, должны быть параллельны. Однако из-за фиксированных значений углов в каждом отдельном узловом шарнире это требование

не выполняется. Тем не менее благодаря податливости складывающихся стержней непараллельность осей цилиндрических шарниров не препятствует полному раскрытию рефлектора. Итак, в развернутом состоянии узловые шарниры могут быть несколько повернуты относительно друг друга. Чтобы определить положение узловых шарниров, будем считать, что сумма квадратов расстояний между точками Apq и Ars двух смежных узловых шарниров с центрами Mp и Mr для каждого пояса имеет минимальное значение:

Ф[^] = ЙП Y1 1 rApq- rArs (W) |2, (6)

^скл.ст

где p = ... Wn}; p,r, G N Л [1, n]; q, s G N Л [1, 6]; ^кл.ст — число

складывающихся стержней рассматриваемого пояса.

В зависимости (6) для лицевого пояса n = ^уз.л, Дскл ст = Дскл ст.л = = 2N0(6N1 + 1) — N1(3N1 — 1), а для тыльного пояса n = ^узт,

^скл.ст ^скл.ст.т ^скл.ст.л 3(2N0 + N1).

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

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

СПИСОК ЛИТЕРАТУРЫ

1. Зарубин В. С. Математическое моделирование в технике: Учеб.для вузов / Под ред. В.С.Зарубина, А.П.Крищенко. - М.: Изд-во МГТУ им.Н.Э.Баумана, 2001. - 496 с. (Сер. Математика в техническом университете; Вып. XXI, заключительный).

2. Стечкин С. Б., Субботин Ю. Н. Сплайны в вычислительной математике. - М.: Наука, 1976. - 248 с.

Статья поступила в редакцию 13.02.2009

Виталий Евгеньевич Мешковский родился в 1946 г., в 1971г. окончил МВТУ им. Н.Э. Баумана. Старший преподаватель кафедры "Космические аппараты и ракеты-носители", по совместительству зав. отделом СМ1-1 НИИСМ МГТУ им. Н.Э. Баумана. Автор 37 научных публикаций в области проектирования и динамики раскрывающихся крупногабаритных космических конструкций.

V.Ye. Meshkovskii (b. 1946) graduated from the Bauman Moscow Higher Technical School in 1971. Senior teacher of "Spacecrafts and Rocket Boosters" of the Bauman Moscow State Technical University. Author of 37 publications in the field of design and dynamics of unwrapping large-sized space constructions.