CC BY
92
UDC 629.78:531.395
Doi: 10.31772/2587-6066-2019-20-1-68-73
For citation: Zhang Zikun, Zimin V. N., Krylov A. V., Churilin S. A. [The definite questions of simulation of transformable space structures dynamics]. Siberian Journal of Science and Technology. 2019, Vol. 20, No. 1, P. 68-73. Doi: 10.31772/2587-6066-2019-20-1-68-73
Для цитирования: Джан Ц., Зимин В. Н., Крылов А. В., Чурилин С. А. Некоторые аспекты моделирования динамики трансформируемых космических конструкций // Сибирский журнал науки и технологий. 2019. Т. 20, № 1. С. 68-73. Doi: 10.31772/2587-6066-2019-20-1-68-73
THE DEFINITE QUESTIONS OF SIMULATION OF TRANSFORMABLE SPACE STRUCTURES DYNAMICS
Zikun Zhang, V. N. Zimin*, A. V. Krylov, S. A. Churilin
Bauman Moscow State Technical University 5, 2-nd Baumanskaya St., building 1, Moscow, 105005, Russian Federation E-mail: zimin@bmstu.ru
This paper describes large transformable space structures with various configuration in the folded transport position and in the open working one. As an example, simulation of transformable space structures dynamics is shown for the antenna circuit foldable load-bearing frame with diameter of 5 m. For investigation of the foldable frame deployment dynamics, a design scheme presented by a system of rigid bodies connected with each other by hinges is accepted as it is simple, but at the same time it considers features of the structure well enough. For performing stress analysis of the foldable frame elements during deployment, the frame shape at the certain time point of deployment, when relative velocities of adjacent elements are ultimate, is chosen. As a results of calculation using MSC.Adams software, positions, velocities and accelerations of the centres of mass of the foldable frame elements as well as the angular velocities and accelerations of the elements for each time step of the deployment are obtained. To perform stress analysis of the fold-able load-bearing frame, finite element model of the frame is developed using MSC.Patran/Nastran software. As a results of investigation of stressed and deformed states of antenna circuit foldable frame elements both without taking into account damping and with consideration of damping, stresses arising in the foldable frame elements at the certain time points during deployment are found.
Keywords: transformable structures, deployable structures, dynamics, modeling, simulation, movement, configuration, model.
НЕКОТОРЫЕ АСПЕКТЫ МОДЕЛИРОВАНИЯ ДИНАМИКИ ТРАНСФОРМИРУЕМЫХ
КОСМИЧЕСКИХ КОНСТРУКЦИЙ
Ц. Джан, В. Н. Зимин*, А. В. Крылов, С. А. Чурилин
Московский государственный технический университет имени Н. Э. Баумана Российская Федерация, 105005, Москва, 2-я Бауманская ул., 5, стр. 1 E-mail: zimin@bmstu.ru
Рассматриваются крупногабаритные трансформируемые космические конструкции, имеющие разные конфигурации в транспортном состоянии и рабочем положении. Моделирование динамики трансформируемых космических конструкций рассмотрено на примере складного несущего силового каркаса антенного контура диаметром 5 м. Для анализа динамики раскрытия принимается простая, достаточно хорошо учитывающая особенности конструкции расчетная схема в виде системы абсолютно твердых тел, связанных между собой шарнирными соединениями. В качестве расчетной схемы для определения напряженно-деформированного состояния элементов каркаса при раскрытии принимается его форма в определенный момент времени, когда относительные скорости соседних элементов максимальны. В результате расчета в программном комплексе MSC.Adams получены координаты, скорости и ускорения центров масс элементов каркаса, а также их угловые скорости и ускорения. Для определения напряженно-деформированного состояния каркаса в программном комплексе MSC.Patran/Nastran построена его конечно-элементная модель. В результате расчета напряженно-деформированного состояния элементов каркаса без учета и с учетом демпфирования получены значения возникающих в них напряжений в определенные моменты времени.
Ключевые слова: трансформируемые конструкции, динамика, моделирование, движение, конфигурация, модель.
Introduction. The intensive development of space technology poses the task of creating of fundamentally new large space structures [1; 2]. A necessary stage of their development is a preliminary investigation of the characteristics of both free and controlled such structures' motion by means of mathematical modeling of their dynamics.
Mathematical modeling of dynamics of such structures is a subject of wide research. The main methods of analysis of the mechanical behavior of large space structures are described in [3]. Transformable large space structures have different configurations in the transport state and in the operating position in orbit. As a rule, deployment of the developed transformable structures, is unique process for each considered system, however, in some cases it is possible to find a common way of modeling the dynamics of their deployment. A model that is simple but well enough considering the special features of structures consisting of tens, hundreds and even thousands elements interconnected by hinges is accepted as a design scheme [4-9]. System of differential equations of transformable structure elements motion is usually known as a model of dynamics of such structure, and a numerical solution of these equations is usually understood as the mathematical modeling of structure dynamics.
Calculation model for the analysis of the deployment of the circuit antenna frame. Among the large space structures, a special place is occupied by the multilink transformable circuit antennas, the load-bearing frame of which provides the minimal dimensions of the structure in the folded transport position and required rigidity of the structure in the unfolded working one (fig. 1). The frame consists of one-type elements connected by elastic hinges. The elements can be made of both metals and composite materials (fig. 2). The hinges contain one-sided stops to provide of deployment of the reflector into the opened working position, making desired shape of the reflecting surface its maintenance during operating life. Deployment of considered structures occurs automatically due to potential energy accumulated in the elastic elements of the hinges when the reflector had been being folded into the transport state. Computer modeling of steps of the multi-link load-bearing frame deployment gives us an opportunity to consider various schemes of this process and evaluate correctness of the technical solutions incorporated in the structure's scheme at the project level.
Modeling of deployment of the load-bearing frame with a diameter of 5 m using MSC.Adams and MSC. Patran/Nastran software is under consideration [10].
Fig. 1. The spacecraft with circular transformable antennas with diameters of 20 m
Рис. 1. Космический корабль с трансформируемыми кольцевыми антеннами диаметром 20 м
a b
Fig. 2. The elements of transformable space structures: a - shaped rectangular panels; b - thin-walled quill rods of unidirectional carbon fiber reinforced polymer
with longitudinal fibers orientation
Рис. 2. Элементы трансформируемых космических конструкций: а - прямоугольные профилированные панели; б - тонкостенные трубчатые стержни из однонаправленного углепластика с продольным расположением волокон
The folding load-bearing frame consists of two packets of shaped rectangular panels 550x 230x 6 mm (13 panels per package) one side of which is pivotally connected with the spacecraft, and another side the same way is connected with short closing panel.
As a design scheme for studying of dynamics of the antenna folding frame deployment using MSC.Adams, a system of rigid bodies interconnected by hinges is accepted. At a certain relative position of the adjacent panels during the deployment, restrictions limiting mutual angular displacement of panels are applied to provide holding of the antenna frame in the operating unfolded position. The panels of the foldable antenna frame are connected with each other by external and internal hinges modeled by means of "axial joint" elements (Revolute Joint in terms of Adams software). Torsion springs in the hinges are modeled by means of elastic and dissipative elements with linear dependencies of the moments from the opening angles and the relative angular velocities of the panels (Torsion Spring in terms of Adams software), the stiffness coefficient, viscous resistance coefficient and preliminary angles of twist of the springs for which are given. The external and internal hinges have different preliminary angles of twist, which have been chosen the way to provide a circle shape of the folding load-bearing frame in the operating unfolded position.
The panels are modeled by rigid rectangular parallelepipeds with dimensions of 550 x 230 x 6 mm. The density of the equivalent material of the panels in Adams is chosen so that the mass of the simulated non-profiled panels is equal to the mass of the real profiled panels. The moments of inertia of the panels in Adams are taken equal to the moments of inertia of the real profiled panels.
Free ends of the panels placing in the base of the left and right packets are considered to be hinged.
The power characteristic of each of the elastic and dis-sipative elements modeling the torsion spring in the hinges [11-13] is determined by the following ratio
Mt (9,., 9i) = ct -9,.) cp,
where c, is the stiffness coefficient of the i -th spring; | is the coefficient of viscous resistance of the i -th spring; 9,. is the current opening angle of adjacent panels; 9^ is
the preliminary angle of twist of the i -th spring.
At a certain relative position of the adjacent panels when they are opening, and the opening angle c reaches a certain value pstop corresponding to operating position
of panels, restrictions limiting the mutual angular displacement of the panels are applied to panels by the certain making of their hinges. Physically the restrictions are made as various stops, which are modeled by massless elastic and damping elements with nonlinear dependency of the moment from the opening angle and can be written as
M ( • ) i if 9i< 9stop'
Mstop1 (Pi, CPi) = i , , • .. . ,
" [-CStopi (9i -Pstopi) Istopj Cpi, f Pi ^Pstop1
where cstop is the stiffness coefficient of the i-th elastic element; |stop is the coefficient of viscous resistance of the i-th damping element; pst is value of the opening
angle of adjacent panels, when the placing on the stop is getting; 9 is the relative angular velocity of the adjacent panels.
When moving, adjacent panels can rotate towards each other and touch each other. The model contains stops preventing a contact between the panels. They are presented by massless elastic elements with nonlinear dependency of the moment from the opening angle
^ , X f 0 if Pi
Mcont_(Pi) =i ( ) f ^
" [-Ccont1 (9i -Pcont,), if Pi ^Pcont1,
where cconti is the stiffness coefficient of the i-th elastic element; 9cont is value of the i-th opening angle, when
the adjacent elements are in contact.
As a result of calculations, dependencies of coordinates, speeds and accelerations of centres of mass of the load-bearing frame panels from time, as well as dependencies their angular speeds and accelerations from time were received.
Analysis of stressed state of the frame elements. The stressed state of the elements of the load-bearing frame during its deployment is determined by the impact loads arising when the adjacent frame panels are placing on the stops. The impact loads are obtained from the analysis of the frame deployment dynamics using the MSC. Adams software.
Simulation shows that the panels of the load-bearing frame are placed on the stops at the different time points. At each time point, various groups of panels in different places of the frame reach their stops. Several time points are chosen when relative velocities of the adjacent panels from the certain group are ultimate.
To investigate the stressed state of the load-bearing frame panels during deployment, in MSC. Patran/Nastran software the frame shape at the selected moments is adopted as a design scheme. At each considered moment it is supposed that panels are placed on their stops and the structure behaves as an elastic rod with specified characteristics. This approach gives the safety margin, since the mobility of frame panels relative to each other leads to reduction of real stresses due to the loss of kinetic energy in the joints.
Knowing value of the angular velocity of the hinged first panel of one of the packets and knowing values of the angular velocities of the other panels of the packet, the velocity field at the certain points of each panel of the packet is determined by the formula v n = v + © x r ,
where vp is the speed of the poles, for which each joint
on the antenna frame from the attachment base point to the short closing panel is consistently accepted; © is the angular velocity vector of the corresponding panel; r = rn - rp , rp is the radius vector of the pole, rn is the
radius vector of the point, where the velocity is determined. The radius vectors are specified in the inertial frames of reference whose origins are located at the places of hinge fastening for each frame packet on a rigid base. Velocities of the opposite packets panels are taken equal in magnitude and opposite in direction. Thus, the opening of the frame occurs symmetrically.
The angular velocities of the panels as well as the velocities of their centres of mass (hereinafter used to test calculation of velocities of the selected model points) obtained for each time step of the deployment calculation are exported from Adams to a certain file as a table. To find velocities of the model selected points from the dependency given above for each selected moment of time a C++ program was written.
To perform stress analysis of the frame panels, finite element model of the frame was built using MSC.Patran/Nastran software. Analysis of the stressed and deformed state of antenna frame was carried out by means of Direct Transient Response analysis (SOL 109, study of transient process) both without taking into account damping and with consideration of damping.
Each panel is divided into ten "beam" finite elements. In total, the model consists of 270 such elements. All finite elements including elements at the panel boundaries are rigidly connected to each other. The panels placed at the base of the left and right packet are considered to be clamped by their free ends.
As a dynamic model of the folding frame the finite element model is accepted
M u + D u + Ku = P,
where M, D, K are mass, damping and stiffness matrices respectively, P is the external loads vector, i is the vector of nodal displacements, u is the vector of nodal velocities. In our case P = 0 .
Calculated values of velocities in the selected points of the frame model are taken as initial conditions for calculation of the transient process using MSC. Nastran software, that is, when t = 0: u (0) = u0, u (0) = u 0, where
u0, u 0 are vectors of the initial nodal displacements and velocities.
In accordance with the recommendations presented in the MSC. Nastran user manuals, transient response
a
analysis time integration step should be selected taking the conditions as at least ten steps per response period time for the maximal frequency of interest. Then, if we take into account first ten frequencies of natural oscillations of the antenna frame in the plane of deployment, in our case the time integration step will be equal to
T 1 1 At =-=-« 0,01sec.
10 10 • f10 10 • 9 Hz
At the same time, it is recommended to take the time integration step finding within At = 10 4...10 3 sec . interval, so to obtain more accurate stress values, the time integration step during the calculations was taken equal to At = 103 sec.
Structural damping in Direct Transient Response analysis is taken into account by usage of the equivalent viscous damping. For the sinusoidal displacement response with constant amplitude, the structural damping force is constant and the viscous damping force is proportional to the frequency. For such a response, the two damping forces are equal at some frequency. Having done calculations without damping it was clear that the frequency of the structure oscillations under the impact had closed value to the second natural oscillations frequency in the plane of the opening of the frame. So the second natural oscillations frequency was chosen as a frequency of equivalence of structural damping and viscous one.
As a results of investigation of stressed and deformed states of antenna circuit foldable frame elements both without taking into account damping and with consideration of damping, stresses arising in the foldable frame elements at the certain time points during deployment are found.
Fig. 3 shows the stressed state of the antenna frame at the time point of t = 5.925 sec., when small stresses occur, looking all considered cases. The values of the stresses occurring in the panels are shown the same.
b
Fig. 3. Stressed state of the antenna frame at time point of t = 5.925 sec.: a - without consideration of damping; b - with consideration of damping
Рис. 3. Напряжения в момент времени t = 5,925 с.: a - без демпфирования; б - с демпфированием
Conclusion. It is possible to reduce or even exclude impact loads during the process of uncontrolled deployment by reducing the initial potential energy of the springs, i.e. by changing the elastic characteristics of the spring elements included in the structure. At the same time, the angular velocities at the time of placing of adjacent structural elements on the stops will accordingly decrease. This will lead to decreasing of the velocities determining the magnitude of the impact impulse acting on the considered elastic transformable structure. However, this way is not always possible, since besides the smooth and reliable deployment of the transformable structures, it is necessary to guarantee their following operation in orbit. Therefore, in order to provide a controlled transformation of the system from the transport folded position to the operating opened one, it is necessary either to develop the structure with an additional cable synchronization system, or to replace the springs by operating drives. This will significantly complicate the structure and increase its mass. It is possible to create large transformable space structures in orbit, using for their deployment actuators made of materials with a shape memory effect [14; 15].
There are several types of actuators based on materials with shape memory effect. One of them can be called a one-dimensional direct-acting actuator, the another one -an actuator with aftereffect, which uses a spring or differential drive making an active retarding force and damping for creating force of actuator itself.
References
1. Lopatin A. V., Rutkovskaya M. A. [The review of designs of modern transformed space antennas. Part 1]. VestnikSibGAU. 2007, No. 2(15), P. 51-57 (In Russ.).
2. Lopatin A. V., Rutkovskaya M. A. [The review of designs of modern transformed space antennas. Part 2]. Vestnik SibGAU. 2007, No. 3(16), P. 78-81 (In Russ.).
3. Banichuk N. V., Karpov I. I., Klimov D. M. et al. Mekhanica bol'shikh kosmicheskikh konstruktsiy [Mechanics of large space structures]. Moscow, Faktorial Publ., 1997, 302 p.
4. Wittenburg J. Dinamika system tviordykh tel [Dynamics of multibody systems]. Moscow, Mir Publ., 1980, 292 p.
5. Boykov V. G. [Simulation of the mechanical systems dynamics using EULER software]. SAPR i Graphika. 1998, No. 1, P. 38-48 (In Russ.).
6. Boykov V. G. [Software complex for the automated dynamic analysis of multicomponent mechanical systems EULER]. SAPR i Graphika. 2000, No. 9, P. 17-20 (In Russ.).
7. Bakulin D. V., Borzikh S. V., Ososov N. S., Shchiblev Yu. N. [Modeling dynamics of opening large-sized solar batteries]. Izvestiya RAN. Matematicheskoe modelirovanie. 2004, Vol. 16, No. 6, P. 88-92 (In Russ.).
8. Borzikh S. V., Ilyasova I. R. [Simulation and experimental development of deployment of large-sized multi-link solar batteries of spacecrafts]. Engineering Journal: Science and Innovation. 2012, No. 8 (8), P. 60-68 (In Russ.).
9. Kuznetsova A. O. [Research of the dynamics of the motion of uncovered mechanical systems with elastic ties]. Vestnik SibGAU. 2005, No. 3, P. 135-138 (In Russ.).
10. Schesniak S., Romanov A., Khitrov I. et al. [Designing and calculation of large-sized deployable structures using MSC. Software simulation software technology]. CADmaster. 2009, No. 2-3 (47-48), P. 28-38 (In Russ.).
11. Krylov A. V., Churilin S. A. [Simulation of Deployment of Multi-Link Closed Space Structures]. Engineering Journal: Science and Innovation. 2012, No. 8 (8), P. 80-91 (In Russ.).
12. Krylov A. V., Churilin S. A. [Methods of deployable multilink structures elements stressed state investigation]. Russian physics journal. 2013, Vol. 56, No. 7-3, P. 170-172 (In Russ.).
13. Zimin V. N., Krylov A. V., Meshkovskiy V. Ye. et al. [Features of the Calculation Deployment Large Transformable Structures of Different Configurations]. Science and Education of the Bauman MSTU. Electronic Journal. 2014, No. 10, P. 179-191 (In Russ.). Doi: 10.7463/1014.0728802.
14. Molodtsov G. A., Bitkin B. E., Simonov V. F. et al. Formostabil'nyye i intellektual'nyye konstruktsii iz kompozitsionnykh materialov [Form-stable and intelligent structures made of composite materials]. Moscow, Mashinostroyenie Publ., 2000, 352 p.
15. Liand C., Rogers C. A. Design of shape memory alloy actuators. Journal of intelligent material systems and structures. 1997, Vol. 8, P. 303-313.
Библиографические ссылки
1. Лопатин А. В., Рутковская М. А. Обзор конструкций современных трансформируемых космических антенн. Ч. 1 // Вестник СибГАУ. 2007. № 2(15). С. 51-57.
2. Лопатин А. В., Рутковская М. А. Обзор конструкций современных трансформируемых космических антенн. Ч. 2 // Вестник СибГАУ. 2007. № 3(16). С. 78-81.
3. Механика больших космических конструкций / Н. В. Баничук, И. И. Карпов, Д. М. Климов [и др.]. М. : Факториал, 1997. 302 с.
4. Виттенбург Й. Динамика систем твердых тел. М. : Мир, 1980. 292 с.
5. Бойков В. Г. Моделирование динамики механических систем в программном комплексе EULER // САПР и графика. 1998. № 1. С. 38-48.
6. Бойков В. Г. Программный комплекс автоматизированного динамического анализа многокомпонентных механических систем EULER // САПР и графика. 2000. № 9. С. 17-20.
7. Бакулин Д. В., Борзых С. В., Ососов Н. С., Щиблев Ю. Н. Моделирование процесса раскрытия крупногабаритных солнечных батарей // Известия РАН. Математическое моделирование. 2004. № 6. С. 88-92.
8. Борзых С. В., Ильясова И. Р. Моделирование и экспериментальная отработка процесса раскрытия крупногабаритных многозвенных солнечных батарей космических аппаратов // Инженерный журнал: наука и инновации. 2012. № 8. С. 60-68.
9. Кузнецова А. О. Исследование динамики движения раскрывающихся механических систем
с упругими связями // Вестник СибГАУ. 2005. № 3. С. 135-138.
10. Проектирование и расчет крупногабаритных раскрывающихся конструкций с помощью программных комплексов MSC. Software / С. С. Щесняк, А. В. Романов, И. В. Хитров [и др.] // CAD master. 2009. № 2-3 (47-48). С. 28-38.
11. Крылов А. В., Чурилин С. А. Моделирование развертывания многозвенных замкнутых космических конструкций // Инженерный журнал: наука и инновации. 2012. № 8 (8). C. 80-91.
12. Крылов А. В., Чурилин С. А. Методика определения напряженно-деформированного состояния элементов трансформируемых многозвенных конструкций // Известия вузов. Физика. 2013. Т. 56, № 7/3. С. 170-172.
13. Особенности расчета раскрытия крупногабаритных трансформируемых конструкций различных конфигураций / В. Н. Зимин, А. В. Крылов, В. Е. Мешковский [и др.] // Наука и образование. 2014. № 10. С. 179-191. Doi: 10.7463/1014.0728802.
14. Формостабильные и интеллектуальные конструкции из композиционных материалов / Г. А. Молодцов, Б. Е. Биткин, В. Ф. Симонов [и др.]. М. : Машиностроение, 2000. 352 с.
15. Liand C., Rogers C. A. Design of shape memory alloy actuators // Journal of intelligent material systems and structures. 1997. Vol. 8. P. 303-313.
© Zhang Zikun, Zimin V. N., Krylov A. V., Churilin S. A., 2019
Zhang Zikun - postgraduate student; Bauman Moscow State Technical University. E-mail: 89738415@qq.com.
Zimin Vladimir Nikolaevich - first pro-rector, scientific work pro-rector, Dr. Sc., professor; Bauman Moscow State Technical University. E-mail: zimin@bmstu.ru
Krylov Aleksey Vladimirovich - Ph. D., associate professor; deputy head of the Spacecrafts and Carrier Rockets Department; Bauman Moscow State Technical University. E-mail: kav1982@bmstu.ru
Churilin Sergey Aleksandrovich - principal engineer; Bauman Moscow State Technical University. E-mail: churilin@bmstu.ru
Джан Цыкунь - аспирант; Московский государственный технический университет имени Н. Э. Баумана. E-mail: 89738415@qq.com.
Зимин Владимир Николаевич - доктор технический наук, профессор, первый проректор - проректор по научной работе; Московский государственный технический университет имени Н.Э. Баумана. E-mail: zimin@bmstu.ru.
Крылов Алексей Владимирович - кандидат технических наук, доцент, заместитель заведующего кафедрой космических аппаратов и ракет-носителей; Московский государственный технический университет имени Н. Э. Баумана. E-mail: kav1982@bmstu.ru
Чурилин Сергей Александрович - ведущий инженер; Московский государственный технический университет имени Н. Э. Баумана. E-mail: churilin@bmstu.ru.
