CC BY
3Труды МАИ. 2021. № 121 Trudy MAI, 2021, no. 121
Научная статья УДК 534.1
DOI: 10.34759/trd-2021-121-05
МЕТОД РАСЧЕТА УПРУГИХ КОЛЕБАНИЙ ЦИКЛИЧЕСКИ СИММЕТРИЧНОЙ КОНСТРУКЦИИ
Татьяна Витальевна Гришанина1н, Евгения Егоровна Гусева2
1,2Московский авиационный институт (национальный исследовательский университет), МАИ, Москва, Россия
[email protected] 2evgenia-kesha@mail. ru
Аннотация. Представлен новый подход для составления уравнений колебаний упругих конструкций, имеющих циклически симметричную структуру. Упругие элементы конструкции, работающие на растяжение, кручение и изгиб, соединенные между собой, рассматриваются как конечные элементы. Перемещения и углы поворотов в узлах, которые расположены симметрично относительно центральной оси конструкции, выражаются через шесть обобщенных координат узла. Для такого представления используются тригонометрические функции ортогональные на дискретном множестве равноотстоящих точек. Таким образом, уравнения колебания распадаются на отдельные группы уравнений, имеющих размерность шестого порядка.
Ключевые слова: циклически симметричные системы, стержневые системы, антенна, уравнения колебаний, несвязанные уравнения
Для цитирования: Гришанина Т.В., Гусева Е.Е. Метод расчета упругих колебаний циклически симметричной конструкции // Труды МАИ. 2021. № 121 DOI: 10.34759/trd-2021-121-05
METHOD FOR CALCULATING ELASTIC VIBRATIONS OF A CYCLICALLY SYMMETRIC STRUCTURE
Tatiana V. Grishanina13, Evgeniia E. Guseva2
1,2Moscow Aviation Institute (National Research University), Moscow, Russia 4. [email protected] 2evgenia-kesha@mail. ru
Abstract. The article presents a new approach to the oscillation equations composing of elastic structures with cyclically symmetric structure. A cyclically symmetric elastic system is under consideration. At the к = 0, 1, 2, ..., N - 1 nodes of the system, located on a circle at equidistant points with angular coordinates of 0k, the identic elastic rods of constant cross-section are being connected, working in tension-compression, torsion and bending-shear in and out of the plane of the system. To compose the oscillations equations of the system, both displacement and rotation angles components, symmetrical with respect to the radial plane passing through the k-th node, are being represented as cosine
expansions in the circumferential direction, while skew symmetric ones are represented as sine expansions with wave numbers of n = 0, 1, 2, ..., N/2.
Expressions of potential and kinetic energies for all elements of the system are being composed. With account for the cos n0k and sin n0k functions orthogonality conditions for different n on a system of equidistant points, the terms of these expressions for different n from the set of n = 0, 1, 2, ..., < N/2 are being obtained uncoupled among themselves. As the result, the Lagrange equations in generalized coordinates for a cyclically symmetric system with 6N degrees of freedom disintegrate into separate groups of six equations for each number n being accounted for.
Thus, the solution of the problem of a cyclically symmetric system oscillations of еру 6N order is being reduced to solving a number of problems for uncoupled subsystems of equations of the sixth order, representing separately harmonics n = 0, 1, 2, ..., < N/2. Keywords: cyclically symmetric systems, rod systems, antenna, oscillation equations, uncoupled equations
For citation: Grishanina T.V., Guseva E.E. Method for calculating elastic vibrations of a cyclically symmetric structure. Trudy MAI, 2021, no. 121. DOI: 10.34759/trd-2021-121-05
Введение
Антенна является одним из основных элементов бортового радиокомплекса
космического аппарата, предназначенного для обеспечения связи с наземными и
другими объектами. Космические аппараты обычно оснащены малыми антеннами,
которые, как правило, имеют монолитную структуру. Большие антенны, которые
3
развертываются непосредственно в космосе, представляют собой стержневую систему, на которые натягивается отражающая пленка [1-9]. Схема такой антенны, вынесенной на трехгранной штанге стержневой конструкции, приведена на рис. 1.
В настоящее время составные конструкции разворачивающихся антенн рассчитывают, используя программные комплексы, основанные на методе конечных элементов [10-12]. Для расчета конструкций, имеющих регулярную структуру, можно использовать редуцированные модели с использованием условий цикличности [13, 14], а также - условий ортогональности синусов и косинусов на дискретном множестве равноотстоящих точек [15, 16].
Расчетная модель циклических систем
Рассмотрим метод расчета колебаний циклически симметричных систем, представляющих конструкцию развертывающейся антенны. Данная конструкция состоит из следующих основных элементов (см. рис. 2): 1 - рама из упругих стержней (в развернутом состоянии стержни в местах соединения друг с другом (узлах) жестко фиксируются, образуя правильный многоугольник); 2 - центральная стойка, расположенная в центре рамы, перпендикулярно её плоскости; 3 - растяжки, связывающие узлы рамы с крайними точками стойки; 4 - пленочные элементы, расположенные в плоскости рамы.
Рис. 1. Схема антенны Рис. 2. Циклически симметричная конструкция
Каждый узел такой антенны при пространственных колебаниях имеет шесть степеней свободы, представляющих перемещения в трех ортогональных направлениях и поворотам относительно осей, совпадающих с этими направлениями. Если система имеет N узлов, то общее число степеней свободы (без учета растяжек и стойки) будет равно 6N. Ортогональность функций sinndk и cos ndk на дискретном множестве N равностоящих точек при n =0, 1, ..., < N/2 [16] позволяет выразить перемещения и углы поворота в узлах через эти функции. Тогда для циклически симметричной системы из уравнений Лагранжа при 0 < n < N /2 получим шесть уравнений, определяющих свободные колебания этой системы, для каждого n. Коэффициенты жесткостей и масс элементов можно получить из выражений потенциальной и кинетической энергии системы, записанных в обобщенных координатах.
шк "к
к - 1
¥к+1
к+1
к-1
и
к ик
и
к +1
Рк+1
а)
б)
Рис. 3 Перемещения в узлах рамы
Перемещения и углы поворота узлов целесообразно рассматривать относительно двух систем координат: системы, связанной непосредственно со стержнем; системы, связанной с центром симметрии.
Рассмотрим к-ый стержень рамы, соединяющий узлы к и к + 1, (рис. 3). Его положение определяется углом вк, который отсчитывается от плоскости симметрии:
ж
9к = 90 + 2ка, а = — (к = 0,1,..., N -1); в0 = 0.
Кинематические условия сопряжения в к-ом узле:
и>к = ик БШа + ик СОБа, = ¥к СОБа + Рк БШа, Рк = Рк СОБа - \ук БШа, и^ = ик соБа - Бта, wk = ш, = $.
(1)
ук ~ уук--
В (к+1)-ом узле:
"кк+1= »к+1 СОБа - ик+1 siпа, ¡^+1= ¥к+1 СОБа - Рк+1 эта
pkk+l=Рk+lcosа + Wk+lsiпа, ик+1= uk+lcosа + "k+lsiпа, шк+1= шк+и $к+1=$к+1.
к
N
Ортогональность функций sinn6k и cosn6k при 0 < n <— на дискретном множестве N равностоящих вк точек позволяет принять:
(3)
ük = qicosn9k, (pk = q2cosnOk> wk = ^3cosnOk> Wk = q4sinnOk, Uk = q5 sinnOk, 3k = qe sinnOk,
где q- = q (/) (/ = 1, 2, ..., 6) - обобщенные координаты в узле k = 0.
N
Далее ограничимся рассмотрением случая когда 1 < n < — .
Потенциальная энергия рамы равна сумме потенциальных энергий её элементов, которые состоят из энергии растяжения, энергии кручения и энергии изгиба-сдвига стержня в двух плоскостях [1, 17, 18]. Поперечные сечения стержней в пределах их длин будем считать постоянными. 1) Растяжение
1 Т7Т7 N-1
П = 1 Er^Ui« -U)2, (4)
2 1 k=0
где E - модуль упругости, F - площадь поперечного сечения стержня, l - длина стержня.
Используя выражения (1) и (2) с учетом (3) выражение (4) преобразуется к
виду:
EF N-i
П1 =-^ (A cos nOk + B sin nOk )2 .
l k=0
Здесь A = 2 cos na(q5 cos a sin na + qx sin a cos na), B = -2 sin na(q5 cos a sin na + q sin a cos na).
Функции sinn6k и cosn6k ортогональны при различных n = 0, 1, 2, ... < N/2. Кроме того, выполняются условия
N-1
N
N-1
N
N-1
Zsin2 пвк = —, Zcos2 пвк = —, Zsinпвк cosпвк = 0.
к=0
2
к=0
2
к=0
1 EF N 2 2 Таким образом: П1 =---(A + Bz)
или
^ 1 EF N л/. . -
П1 = —2 cosa sin na + q1 sin a cos na)
(5)
2) Кручение
1 GJ n-1
JJ _ 1 GJ кр V"1/ ,„к
2
кр £ 2
l к=0
(6)
где GJ^ - крутильная жесткость стержня.
После аналогичных преобразований координат выражение (6) принимает вид:
1 GJ N
2 = ——^ 4(q2 cosasinna-q4 sinacosna)
П =
2
(7)
3) Изгиб-сдвиг в плоскости рамы [17, 19]
N-1
Пз= - Z (b (oh- ик )2 - 2сз {vi- vi )(3h+ $) + d3 (3¡к+1-$ )2 + 2с313кк+Л) (8)
2
к=0
где
к к 1 + Ък
Ы = 12EI—T; с, = 6E/-T; d = EI-; к =
Ъ lЪ Ъ l2 Ъ /
, 12 EI Л
1 + -г-
V l2 GFc ,
параметры EI, GFc определяются в плоскости рамы.
После перехода к центральным координатам, с учетом ортогональности получим следующее выражение для потенциальной энергии изгиба в плоскости рамы:
л- 1 N ЛГ1 , ч2
= cos na + q1 cosasin na) +
+ c3(q5 sinacosna + q1 cosasinna)q6 cosna+ (9)
+ d3q6 sin2 na +1 c3/q6 (cos2 na - sin2 na)].
4) Изгиб-сдвиг из плоскости рамы Аналогично (8) получаем:
N -1
k 2 , ,..k\/,,,k , ,,,k\ , j /,.,k ,.,k\2 . о „ 7,,,k ,,,k~
IN -i
П4 = - X (64 (wV wk )2 + 2C4KV W )(wL+ W) + ¿4GWV W )2 + 2C4IWÍ+1WÍ ), (10)
2 k=0
где
^ ^ - +
64 = 12EI-T; c4 = 6Е1—Г; = EI-; ^
4 l3 4 l2 4 l
' 12 EI Л 1 + ^-
v l2 ^Fc ,
-1
параметры EI, GFC определяются для плоскости, перпендикулярной плоскости рамы.
1 N 2
П4 = - у№ q3Sin na-8c4 q3sin na(q2Sinasin na + q4cosacos na) +
+ 4d4(q2 sinacos na- q4 cosasin na) + 2c4l (q4 cos a(cos na- sin na) - (11) - q2 sin2 na(cos2 na - sin2 na) + 4q2 q4 sin a cos a sin nacos na)].
Потенциальная энергия всей рамы
Прама = П1 + П2 + П3 + П4
Её можно записать в виде:
j 6 6
Прама = X X kij,рамаq iq j , (12)
о ^ ^ У
2 i=1 j=1
где kij рама - коэффициенты обобщенных жесткостей рамы (kijpрама = kfl,рама).
EF oí 11
В частности: k11 = N-j-2sin a cos na + 2Nb3 cos asín na, k12 = 0, и т.д.
Аналогично формуле (4) записывается потенциальная энергия растяжения для нижних и верхних растяжек, неподвижно закрепленных на крайних сечениях стойки
1 EF N-1
n = 1 E1 раст Y (u 2 + u 2 ) раст. ~ i /Ли k ,верх k ,ниж. / •
2 'раст. k=0
Учитывая кинематические условия сопряжения верхней и нижней растяжек с к-ым узлом
Uk,верх = Wk sin (- Vk cos( и щ= W cos ( + Vk sin (,
потенциальная энергия всех растяжек примет вид
1 EF n-1 i EF kt
Праст. = 1 //cos2 nffk (q2 + Í12) = 1N (q2 + q2).
2 'раст. k=0 2 'раст. 2
Таким образом, потенциальная энергия всей циклически симметричной конструкции запишется в виде:
1 6 6
П = 1 //klñlq} , (13)
2 i=1 j=1
где kij = kij,рама + kij,раст. ( kij = kji ).
Кинетическая энергия системы в общем случае записывается аналогично (13)
1 6 6
т=1 //mmj, (14)
2 i=1 j=1
где mij - коэффициенты обобщенных масс (miJ- = ).
В первом приближении конструкцию можно представить в виде сосредоточенных в узлах масс, соединенных невесомыми стержнями. Кинетическую энергию растяжек и инерцию вращения можно не учитывать. Тогда кинетическая энергия конструкции записывается в виде:
T =11![m(uk2 + »I + ^)] =1 mN (q12 + qз2 + q52), 2 k=0 22
где т - масса, сосредоточенная в узле.
Уравнения свободных колебаний системы получаем как уравнения Лагранжа в
обобщенных координатах [17, 20] по отдельности для каждой гармоники п = 0, 1, 2,
... < N/2:
6
^(ш1Д] + ^) = 0, (/ = 1,2,..., 6). (15)
} =1
Используя эти уравнения, можно получить собственные частоты и формы колебаний циклически симметричной системы.
Выводы
Приведенный метод преобразования циклически симметричной системы к основным обобщенным координатам, в которых уравнения для каждого п при 0 < п < N /2 распадаются на отдельные группы, может быть также использован и в общих случаях, например для упругого диска с лопатками и др.
Представленный в статье метод позволяет составить удобную программу для компьютера, с помощью которой можно определять соответственные частоты и формы колебаний циклических систем.
Список источников
1. Усюкин В.И. Строительная механика конструкций космической техники. - М.: Машиностроение, 1988. - 392 с.
2. Лопатин А.В., Рутковская М.А. Обзор конструкций современных трансформируемых космических антенн // Вестник Сибирского государственного университета науки и технологи имени Академика М.Ф. Решетнева. 2007. № 2. С. 51-57.
3. Лопатин А.В., Рутковская М.А. Обзор Конструкций современных трансформируемых космических антенн // Вестник Сибирского государственного университета науки и технологи имени Академика М.Ф. Решетнева. 2007. № 3. С. 78-81.
4. Сгадова Н.А., Струлев И.М. Анализ формы отражающей поверхности параболической антенны деформированной под действием весовой нагрузки // Труды МАИ. 2010. № 38. URL: http://trudymai.ru/published.php?ID= 14168
5. Русских С.В., Шклярчук Ф.Н. Динамика раскрытия космической зонтичной антенны, состоящей из многозвенных гибких радиальных стержней // Космонавтика и ракетостроение. 2020. № 2 (113). С. 86-98.
6. Русских С.В., Шклярчук Ф.Н. Расчет формообразования космической зонтичной антенны при сильном изгибе радиальных стержней, связанных по параллелям растяжными тросами // Известия Российской академии наук. Механика твердого тела. 2021. № 5. С. 99-112. DOI: 10.31857/S0572329921050093
7. Куи М.Х. Геометрические модели внешней компоновки солнечных антенн космических аппаратов // Труды МАИ. 2015. № 82. URL: http://trudymai.ru/published.php?ID=58836
8. Русских С.В., Нагорнов А.Ю., Шавня Р.А. Метод численного решения нелинейной задачи формообразования зонтичных антенн, состоящих из гибких радиальных стержней, соединенных растяжимыми тросовыми элементами // 17 -я Международная конференция «Авиация и космонавтика - 2018»: тезисы докладов. (Москва, 19-23 ноября 2018). - М.: Люксор, 2018. С. 341-342.
9. Гряник М.В., Ломан В.И. Развертываемые зеркальные антенны зонтичного типа. -М.: Радио и связь, 1987. - 72 с.
10. Образцов И.Ф., Савельев Л.М., Хазанов Х.С. Метод конечных элементов в задачах строительной механики летательных аппаратов. - М.: Высшая школа, 1985. - 390 с.
11. Баничук Н.В., Карпов И.И., Климов Д.М. и др. Механика больших космических конструкций. - М.: Изд-во «Факториал», 1997. - 302 с.
12. Цурков В.И. Динамические задачи большой размерности. - М.: Изд-во «Наука», 1998. - 288 с.
13. Бидерман В.Л. Теория механических колебаний. - М.: Высшая школа, 1980. -408 с.
14. Шклярчук Ф.Н., Бобылев Д.С. Колебания упругих систем, образованных из последовательно соединенных однотипных модулей // Вопросы прочности тонкостенных конструкций: сборник статей. - М.: Изд-во МАИ, 1989. С. 52-56.
15. Власов В.З. Избранные труды. Т I. - М.: Изд-во АН СССР, 1962. - 528 с.
13
16. Власов В.З. Избранные труды. Т III. - М.: Изд-во АН СССР, 1964. - 472 с.
17. Гришанина Т.В., Шклярчук Ф.Н. Колебания упругих систем. - М.: Изд-во МАИ, 2013. - 100 с.
18. Гришанина Т.В. Задачи по теории колебаний упругих систем. - М.: Изд-во МАИ, 1998. - 48 с.
19. Гнездилов В.А., Гришанина Т.В., Нагорнов А.Ю. Деформация плоской статически неопределимой стержневой системы при потере устойчивости стержней // Труды МАИ. 2017. № 95. URL: http://trudymai.ru/published.php?ID=84435
20. Бабаков И.М. Теория колебаний. - М.: Наука, 1968. - 560 с.
References
1. Usyukin V.I. Stroitel'naya mekhanika konstruktsii kosmicheskoi tekhniki (Construction mechanics of space technology structures), Moscow, Mashinostroenie, 1988, 392 p.
2. Lopatin A.V., Rutkovskaya M.A. Vestnik Sibirskogo gosudarstvennogo universiteta nauki i tekhnologi imeni AkademikaM.F. Reshetneva, 2007, no. 2, pp. 51-57.
3. Lopatin A.V., Rutkovskaya M.A. Vestnik Sibirskogo gosudarstvennogo universiteta nauki i tekhnologi imeni Akademika M.F. Reshetneva, 2007, no. 3, pp. 78-81.
4. Sgadova N.A., Strulev I.M. Trudy MAI, 2010, no. 38. URL: http://trudymai.ru/eng/published.php?ID=14168
5. Russkikh S.V., Shklyarchuk F.N. Kosmonavtika i raketostroenie, 2020, no. 2 (113), pp. 86-98.
6. Russkikh S.V., Shklyarchuk F.N. Izvestiya Rossiiskoi akademii nauk. Mekhanika
tverdogo tela, 2021, no. 5, pp. 99-112. DOI: 10.31857/S0572329921050093
14
7. Kui M.Kh. Trudy MAI, 2015, no. 82. URL: http://trudymai.ru/eng/published.php?ID=58836
8. Russkikh S.V., Nagornov A.Yu., Shavnya R.A. 17-ya Mezhdunarodnaya konferentsiya «Aviatsiya i kosmonavtika - 2018», Moscow, Lyuksor, 2018, pp. 341-342.
9. Gryanik M.V., Loman V.I. Razvertyvaemye zerkal'nye antenny zontichnogo tipa (Deployable umbrella-type mirror antennas), Moscow, Radio i svyaz', 1987, 72 p.
10. Obraztsov I.F., Savel'ev L.M., Khazanov Kh.S. Metod konechnykh elementov v zadachakh stroitel'noi mekhaniki letatel'nykh apparatov (The finite element method in the problems of structural mechanics of aircraft), Moscow, Vysshaya shkola, 1985, 390 p
11. Banichuk N.V., Karpov I.I., Klimov D.M. et al. Mekhanika bol'shikh kosmicheskikh konstruktsii (Mechanics of large space structures), Izd-vo «Faktorial», 1997, 302 p.
12. Tsurkov V.I. Dinamicheskie zadachi bol'shoi razmernosti (Dynamic problems of large dimension), Moscow, Izd-vo «Nauka», 1998, 288 p.
13. Biderman V.L. Teoriya mekhanicheskikh kolebanii (Theory of mechanical vibrations), Moscow, Vysshaya shkola, 1980, 408 p.
14. Shklyarchuk F.N., Bobylev D.S. Voprosy prochnosti tonkostennykh konstruktsii (Issues of strength of thin-walled structures, collection of articles), Moscow, Izd-vo MAI, 1989, pp. 52-56.
15. Vlasov V.Z. Izbrannye Trudy (Selected works). Vol. I. Moscow, Izd-vo AN SSSR, 1962, 528 p.
16. Vlasov V.Z. Izbrannye trudy (Selected works). Vol. III. Moscow, Izd-vo AN SSSR, 1964, 472 p.
17. Grishanina T.V., Shklyarchuk F.N. Kolebaniya uprugikh system system (Vibrations of elastic systems), Moscow, Izd-vo MAI, 2013, 100 p.
18. Grishanina T.V. Zadachi po teorii kolebanii uprugikh system (Problems in the theory of oscillations of elastic systems), Moscow, Izd-vo MAI, 1998, 48 p.
19. Gnezdilov V.A., Grishanina T.V., Nagornov A.Yu. Trudy MAI, 2017, no. 95. URL: http://trudymai.ru/eng/published.php?ID=84435
20. Babakov I.M. Teoriya kolebanii (Theory of vibrations), Moscow, Nauka, 1968, 560 p.
Статья поступила в редакцию 25.10.2021; одобрена после рецензирования 30.10.2021; принята к публикации21.12.2021.
The article was submitted on 25.10.2021; approved after reviewing on 30.10.2021; accepted for publication on 21.12.2021.
