Rotational control of a satellite


by Erik_at_DTU
Tags: control, rotational, satellite
Erik_at_DTU
Erik_at_DTU is offline
#1
Mar11-06, 12:36 PM
P: 28
Hey,

We're working on a mathematic project, which is about rotational control of a satellite. It's in some ways connected to physics,because we're supposed to analyse the dynamics of a satellite. First a one-axis, then a three axis model.

Anyway, most of the tasks are about differential equations, and it doesn't involve any advanced physics.

ONE-AXIS MODEL - problem formulation
A simplified model is given by:
[tex] I \dot{ \omega } (t) = N_{ctrl} + N_{ext} \qquad (1) [/tex]

where [tex]\omega(t)[/tex] is the angular speed as a function of time, [tex]\dot{\omega}(t)[/tex] is the derived function. [tex]N_{ctrl}[/tex] is the control torque generated by a momentum wheel (reaction wheel), and [tex]N_{ext}[/tex] is the disturbance. [tex]I[/tex] is the inertiemomentum for the satellite. Finally we set the angular position to be [tex]g(t)[/tex], so we have:
[tex] \dot{g}(t)=\omega(t) \qquad (2) [/tex]

The disturbance [tex]N_{ext}[/tex] on the satellite is a torque caused by air resistance. It is modelled by the following expression:
[tex]N_{ext}=N_0(1+\sin^2(2 \pi \frac {t}{T_0})) \qquad (3) [/tex]

Where [tex]T_0[/tex] is the periode and [tex]N_0[/tex] is a constant. This can be rewritten by knowing:
[tex]\sin^2(2 \pi \frac{t}{T_0})=\frac{1}{2} (1-\cos(4 \pi \frac{t}{T_0})) \qquad (4) [/tex]

So [tex]N_{ext}[/tex] is a constant influence superimposed with a periodical influence with the periode [tex]2 T_0[/tex]

The reaction wheels influence on the system is described by a constant [tex]N_{nom}[/tex] multiplied with the control signal u :
[tex]N_{ctrl}=N_{nom}u \qquad (5) [/tex]

We have the following data:

[tex]I=14.33 [/tex] [tex]kgm^2 [/tex]
[tex]T_0=90[/tex] [tex]min = 5400[/tex] s
[tex]N_0=2 \cdot 10^{-4}[/tex] [tex]Nm [/tex]
[tex]N_{nom}=0.010[/tex] [tex]Nm[/tex]

Now we are ready to start on the tasks. There will be around 10 tasks, but we won't ask any questions before running into problems, so some of the tasks (like the first one) will just be uploaded with a solution. So if we make any mistakes, please make a point of it :)

Task 1
Find the complete solution to the homogene differential equation [tex]I \cdot \ddot{g}(t)=N_{ctrl}+N_{ext}[/tex], when [tex]N_{ctrl}=0[/tex] and [tex]N_{ext}=0[/tex]. Then find solutions to the following initial conditions:

[tex]1: \qquad \omega(0)=0, \qquad g(0)= \frac{1}{2} [/tex]
[tex]2: \qquad \omega(0)=0.01, \qquad g(0)= \frac{1}{2}[/tex]

SOLUTION
Complete solution: [tex]g(t)=C_1 \cdot t+C_2[/tex], where [tex]C_1[/tex], [tex]C2[/tex] are constants.

[tex]1: \qquad \frac{1}{2}[/tex]
[tex]2:\qquad \frac{1}{100} \cdot t+\frac{1}{2}[/tex]

Task 2
The satellite can be regulated with an regulator, given by:
[tex]u=k_1 \cdot \omega(t)+k_2 \cdot g(t) \qquad (6)[/tex]

Where [tex]k_1[/tex] and [tex]k_2[/tex] are regulator parameters. Set up a condition (state) model (a linear differential equation of second order) for the closed loop system given by equations (1), (2) and (6).

Show that it's possible to write the model as:
[tex]I \cdot \ddot{g}(t)-N_{nom} \cdot k_1 \cdot \dot{g}(t)-N_{nom} \cdot k_2 \cdot g(t)=\frac{3}{2} \cdot N_0-\frac{1}{2} \cdot N_0 \cdot \cos(4\pi \frac{t}{T_0}) \qquad (7)[/tex].

SOLUTION
Can somebody help us with figuring out what a condition model is?

Now our answer is just:
[tex] I \cdot \ddot{ g } (t) = N_{ctrl} + N_{ext} [/tex]

Which if we move [tex]N_{ctrl}[/tex] over to the left side, and use equations (3)-(5), would give us equation (7).

Any suggestions?

More tasks to come (hopefully with solutions :))
Phys.Org News Partner Science news on Phys.org
SensaBubble: It's a bubble, but not as we know it (w/ video)
The hemihelix: Scientists discover a new shape using rubber bands (w/ video)
Microbes provide insights into evolution of human language

Register to reply

Related Discussions
dimmer or router speed control as resistance heater level control Electrical Engineering 11
Physics help: rotational mechanics satellite Introductory Physics Homework 4
Rotational spectra - thermal population of rotational levels Advanced Physics Homework 0
Rotational control of a satellite - mainly diff. eq. Introductory Physics Homework 0
Rotational equilibrium, subtending, and rotational acceleration Introductory Physics Homework 1