Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.ssau.ru/files/science/conferences/itnt2015/itnt_2015_2.pdf Дата изменения: Tue Aug 4 10:53:46 2015 Дата индексирования: Mon Apr 11 03:33:08 2016 Кодировка:

RESEARCH OF EQUILIBRIUM STATE OF SPACE TETHER SSTEM CONSIDERING ATMOSPHERIC D RAG
D , L A, Cq N P U
The equilibrium equations of tether and subsatellite will be given with the consideration of tether mass, tether deformation and atmospheric drag. Space tether system equilibrium state will be obtained through solving the equilibrium equations of tether and subsatellite. In the part of simulation and analysis, two cases will be simulated. One case is that tether length is fixed and subsatellite mass is variable and the other is subsatellite mass is fixed and tether length is variable. The results of equilibrium state, tether tension and tether angle of deviation will be analysed and contrasted. Simulation results provide a reference for choosing the parameters of tether system in future space applications.

Key words atmospheric dragspace tether systemequilibrium state

1. Introduction B ' urgent understanding of atmospheric environment, meanwhile in order to adapt to modern atmospheric science research and satisfy the needs of national defense, aerospace, communication , etc, the study of atmospheric sounding technology has always been a research hotspot for national researchers. Due to the existence of atmospheric drag for atmosphere of 100 to 200 kilometers high, it is hard for traditional satellite to keep running in this height range. So it cannot obtain data effectively. However, atmosphere of this height range can be detected effectively by using space tether system equilibrium state. Owing to different payload masses and the height of payload, gravity gradient force and atmospheric drag effects on payload are different. So tether critical lengths of space tether system with various probes are also different. In order to make sure that space tether system with atmospheric probe can work regular at equilibrium state and avoid tether fracture, it is necessary to analyze and calculate the equilibrium of space tether system. Berryman J etc. [1] discuss the possibility of space tether system running on LEO with more than one payloads equilibrium state existence. Natarajan[2] studies on two-body space tether system stabilization on equilibrium state and analyzes three equilibrium states by using linear time-dependent dynamics model. Pettazzi, L. etc.[3] put forward a distributed navigation technology based on equilibrium configuration, which can make several satellites arrive at designated spots in space. Bensong Yu[4] mentions equilibrium state modeling ignoring various external perturbation factors. This paper analyses space tether system equilibrium state of different influence factors. Atmospheric drag and tether flexibility are taking into consideration in dynamics modeling. 2. Space Tether System Equilibrium State Analysis Assumptions as follows are made during dynamics modeling. (1) spacecraft runs in circular orbit and its mass is much bigger than satellite and tether under spacecraft. (2) tether flexibility is considered. (3) tether cross section is rounded. (4) spacecraft motion is limited in orbital plane. (5) payload shape is sphere. Payload radius can be ignored compared to tether length.

10


2.1 Tether Equilibrium State Analysis As tether flexibility is considered, continuous medium model is adopted. Equation(1) can be obtained by force analysis on tether microelement.
(S )S
2r T ( S S , t ) T ( S , t ) qS t 2

1
r

where T (S S , t ) T (S , t ) is tether force increment, S is tether curve coordinate, tether microelement radius vector, t is time, q is distributed load. Equation(1) divides S and make S 0 . Tether motion equation can be obtained.
(S )
2 r T q t 2 S

is

2

Considering flexible tether cannot bear lateral load, tether force direction is along tangential.
T T

1 r S

3

where is unit tangent vector and

r . S

Tether equilibrium state condition can be obtained from tether motion equation.
dT q 0 dS

4

This equation is applied to both unstretchable tether and stretchable tether. In order to make reaerch results universal, stretchable tether is studied in this paper. There are two kinds of distributed load q [5]: mass force such as gravity and surface force such as aerodynamic force. Gravity is described as:
qg

g T
(T )

5 is tether radius rate of

where g is acceleration of gravity. is tether linear density. change. Inertia force is described as:
q

T





ae

6

where ae 2 r is acceleration vector. is space tether system orbit angular velocity.
r x2 y

' . Aerodynamic force is described as:
2

1 qR C a dTVV sin 2

7

where V is velocity tether microelement relative to atmosphere, C is aerodynamic drag coefficient, a is atmospheric density. is tether microelement angle of attack, d T is tether diameter. Research on space tether system equilibrium state is usually based on orbit motion coordinate system[5]. Interaction of tether microelement with distributed load is as Fig.1.

11


where is tether microelement tangential direction, n is tether microelement normal direction, q is distributed load interaction with tether microelement. O ' . S1 is spacecraft.
dx sin dS dy cos dS

8 9

where x y are tether microelement coordinates.
Y
Y S1

S1

k

TS

e

n O

q
G

S2

X
O



R
X

k

Fig. 1 Interaction of tether microelement with distributed load

Fig. 2 Forces operating on the payload

Gravity, inertia force and aerodynamic force acting on tether microelement are projected on microelement tangent direction and normal direction. Equation(10) can be obtained combing Equation(4).
dT ( g 2 r ) sin( ) qR dS (T ) T


10
n

d (2 r g ) cos( ) qR dS (T )
n

11

1 1 where qR CdT aV 2 | sin | cos , qR CdT aV 2 sin 2 , sin x / r , cos y / r .


2

2

Equation(8),(9),(10) and (11) are tether microelement equilibrium state equations. 2.2 Payload Equilibrium State Analysis In order to obtain initial value condition to solve tether equilibrium state equations, payload needs to be force analyzed and its equilibrium state equations need to be built. Payload force situation is as Fig.2, where S1 is spacecraft. S 2 is payload. O ' core. TS is tether force. G is payload gravity. R is payload aerodynamic force. e is payload centrifugal force. Payload is treated as mass point. Payload equilibrium state condition is as follows:
TS G R e 0

12

Payload equilibrium state equations are as follows:
(G e ) sin k R cos k TS cosk 0

13

12


(G e ) cos k R sin k TS sin k 0

14

where G

2 mK , e m2 rk , R CS a (rk )V Sm . rk is payload center of mass radius rk2 2

vector. m is payload mass. K is earth gravitational constant. is spacecraft orbit angular velocity. CS is payload front drag coefficient. Sm is payload cross sectional area. V is payload linear velocity relative to atmosphere. Equations(13),(14) are payload equilibrium state equations. 3. Equilibrium State Calculation and Simulation Simulation parameters are as follows: K 3.986 105 km3 /s2 , the earth mean radius:
R3 6371.02km ,

spacecraft

orbital

altitude: H 0 300km ,

tether

linear

density:

0.1667kg/km , tether radius: r0 3 104 m tether Young modulus: E 2.5 1010 N/m2 ,

tether drag coefficient: C 2.2 , payload front drag coefficient: CS 2 . 3.1 Equilibrium State Calculation and Analysis when Tether Length is fixed Tether length is selected as S 150km .

Fig. 3 Equilibrium position when tether length

Fig. 4 Relationship between tether tension and natural coordinates when tether length S 150km
300 S=150km S=176km S=180km S=181km
Y/k m

S 150km
7 6 5 4

250

/°

200

3 2 1 0

m=50kg m=35kg m=16.0035kg m=15kg 0 50 100 150

150

100

0

10

S/k m

20 X/km

30

40

Fig. 5 Relationship between tether angle of deviation and natural coordinates when tether length S 150km

Fig. 6 Equilibrium position when subsatellite mass m 50kg

13


It can be seen from the simulation result Fig.3 that the horizontal distance between subsatellite and spacecraft increases along with subsatellite mass decreases when the system reaches equilibrium position. From Fig.4 it can be seen that tether tension of the junction point of subsatellite gradually decreased along with subsatellite mass decreases. Tether tension variation trend which is from the junction point of subsatellite to the junction point of spacecraft changes from increasing to decreasing. When subsatellite mass m 16.0035kg , tether tension of the junction point of spacecraft reaches zero. When subsatellite mass m 16.0035kg , tether tension appears negative value. It means tether is loosed. Fig. 5 shows variation trend of tether angle of deviation. It can be seen that when subsatellite mass gets smaller, it is more influenced by atmospheric drag and may cause tether loosed. 3.2 Equilibrium State Calculation and Analysis when Subsatellite Mass is fixed Subsatellite mass is selected as m 50kg
40
20 S=150km S=176km

30

15

S=180km S=181km

T/N

/°

20

S=150km S=176km S=180km S=181km

10

10

5

0 0 50 100 S/km 150 200
0 0 50 100 S/km 150 200

Fig. 7 Relationship between tether tension and natural coordinates when tether length m 50kg

Fig. 8 Relationship between tether angle of deviation and natural coordinates when tether length m 50kg

It can be seen from Fig.6 that distance between subsatellite and spacecraft increases along with tether length increases when the system reaches equilibrium position. From Fig.7 it can be seen that tether tension of the junction point of subsatellite gradually increased along with tether length increases. Tether tension variation trend which is from the junction point of subsatellite to the junction point of spacecraft changes from increasing to decreasing. When tether length S 181km , tether tension appears negative value. It means tether is loosed. Fig. 8 shows variation trend of tether angle of deviation. It can be seen that when tether length gets longer, it is more influenced by atmospheric drag and may cause tether loosed. 4. Conclusion Tether and payload equilibrium state equations are built in this paper by analyzing space tether system equilibrium state influenced by atmospheric drag in this paper. Space tether system equilibrium state is obtained by calculating Tether and payload equilibrium state equations. It is obtained by simulation that when tether length is fixed, tether will get loosed if subsatellite mass larger than its critical value. When subsatellite mass is fixed, tether will get loosed if tether length longer than its critical length.
References [1] Berryman J., Schaub, H.. Static Equilibrium Configurations in GEO Coulomb Spacecraft Formations [C]. AAS/ AIAA Space Flight Mechanics Meeting. San Diego: American Astronautical Society, 2005. [2] Natarajan, A.. A Study of Dynamics and Stability of Two-Craft Coulomb Tether Formations [D]. Blacksburg: Virginia Polytechnic Inst. and State Univ, 2007. [3] Pettazzi, L., Izzo, D.. Self-Assembly of Large Structures in Space Using Intersatellite Coulomb Forces [C]. 56th International Astronautica Congress. Fukuoka: International Astronautical Federation / AIAA ( IAF ), 2005. [4] Yu Bensong. Dynamics and Control of Flexible Tethered Satellite in Complex Space Environment [D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2011. [5] Zabolotnov Yuriy. Intoduction to Dynamics and Control in Space Tether System [M]. Beijing: Science Press, 2013:105-115.

14