(adsbygoogle = window.adsbygoogle || []).push({}); Differential equations!

:rofl: 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 signalu:

[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

**Physics Forums - The Fusion of Science and Community**

# Rotational control of a satellite

Have something to add?

- Similar discussions for: Rotational control of a satellite

Loading...

**Physics Forums - The Fusion of Science and Community**