POLYTECHNI C of BARI FACULTY of ENGI NEERI NG COURSE of DEGREE in - - PowerPoint PPT Presentation

polytechni c of bari
SMART_READER_LITE
LIVE PREVIEW

POLYTECHNI C of BARI FACULTY of ENGI NEERI NG COURSE of DEGREE in - - PowerPoint PPT Presentation

POLYTECHNI C of BARI FACULTY of ENGI NEERI NG COURSE of DEGREE in ELECTRONI C ENGI NEERI NG GRADUATION THESIS Bistatic SAR for Earth observation: application of the Dynamic I nversion technique to the orbital and attitude control


slide-1
SLIDE 1

POLYTECHNI C of BARI

FACULTY of ENGI NEERI NG

COURSE of DEGREE

in ELECTRONI C ENGI NEERI NG

slide-2
SLIDE 2

GRADUATION THESIS

Bistatic SAR for Earth

  • bservation: application of

the “Dynamic I nversion” technique to the orbital and attitude control

slide-3
SLIDE 3

Chairman :

  • Prof. Luciano GUERRI ERO

(Polytechnic of Bari) Co-chairman :

  • Eng. Leonardo MAZZI NI

(Alenia Spazio – Rome) Student:

Andrea TEDESCO

slide-4
SLIDE 4

Thesis Aim

To verify the feasibility and, in which case, to plan a low cost orbital and attitude control system, even using some magnetic- torquer.

slide-5
SLIDE 5

“BI SSAT “ Mission

(BI static Sar SATellite)

The mission deals with a SAR of bistatic type (bistatic SAR - Synthetic Aperture Radar), equipped with receiving

  • nly antenna

(microwave system), for Earth

  • bservation.
slide-6
SLIDE 6

Mission Features

The BISSAT will fly in formation tandem with one of the satellites from the COSMO/SkyMed

  • constellation. The

proposed system does will not require any modification of existing COSMO design and

  • perations.

The proposed system will carry some meaningful

advantages which:

  • Low cost of the mission,

thanks also to the drastic reduction of the weight and the dimensions of the satellite (uses only receiving antenna);

  • Added value to the

COSMO mission.

slide-7
SLIDE 7

Scientific Purposes

The bistatic mission allows the simultaneous acquisition of two set of data of the same area under different angles of view, that allows to obtain:

  • A digital elevation model (DEM) of the

ground of observed area, which can be useful:

  • 1. to map the topography of the territory

and to characterize eventual crumbling (liable to slide down) motions;

  • 2. to correct the data of the satellite master.
slide-8
SLIDE 8
  • Low rugged images (lakes, seas,…),

allowing:

  • 1. The implementation of new models for

the location of algae or superficial pollution.

  • Stereoscopiche images;
  • Techniques of interferometric
  • bservation in real time in the hypothesis
  • f being able to realize compatible

baseline with the maintenance of the interferometric coherence.

slide-9
SLIDE 9

COSMO/ SkyMed

  • COSMO (COnstellation
  • f Small satellite in the

Mediterranean basin Observation) is an ASI space program for Earth observation, for scopes: 1. Management of the environmental risks; 2. Commercial applications; 3. Military applications.

  • I t is a constellation of 7

satellites:

  • 4 equipped with SAR;
  • 3 with optical instruments.

I t offers: 1. High and highest resolution images ; 2. Fast system’s response time; 3. Ability of acquisition in every moment; 4. Ability to acquire data of immense areas with single passage; 5. Ability of along-track stereoscopic vision and interferometry; 6. Total accessibility also through mobile reception stations.

slide-10
SLIDE 10

Choice of the Type of Formation

Type of formation that allows parallel orbits of the two satellites to introduce a difference between the right ascensions of the ascending nodes () and a difference between the transit times of the same ascending nodes (t), in

  • rder to avoid every

eventual danger of collision between them near the poles.

slide-11
SLIDE 11

Orbital requirement of the two Satellites

COSMO

  • Altitude:

hC = 619 km

  • Inclination:

iC = 98.5°

  • Right ascension of the

ascending node:

C = 0° (reference orbit)

  • Off-nadir sighting angle
  • f the antenna:

 = 35°

BI SSAT

  • Altitude: hB= 619 km
  • Inclination: iB= 98.5°
  • Right ascension of the

ascending node :

B = 1.5°

  • Angular delay:

B = 0.8°

  • Minimum safety range

from COSMO: 100 km

slide-12
SLIDE 12

Minimal and Maximum Distance Produced

The previous orbital elements produce the following minimal and maximum distances:

Dmin = 97.7348 km Dmax = 276.604 km

slide-13
SLIDE 13

Distance and Slant-Range

Orbital Period: 97.084 min

(the same one for both satellites inasmuch as it is considered the same altitude)

  • Periodic course of the distance

between BISSAT and the COSMO

  • SLC = 774.5718 km

(constant for all orbital period since it always aims with the same angle- shot)

SLB(min) = 690.216 km SLB(Max) = 895.488 km

(because of its continuous pursuit of the trace on the earth)

slide-14
SLIDE 14

Propagation of the two Satellities

slide-15
SLIDE 15
slide-16
SLIDE 16
slide-17
SLIDE 17

Nominal Attitude Angles

  • Bissat, in order to follow the

trace to COSMO earth, needs of two angles:

  • “Azimuth” angle:

B(min) = 23.924 ° B(max) = 43.261 °

  • “Elevation” angle:

B(min) = 24.977 ° B(Max) = 43.559 °

slide-18
SLIDE 18

DYNAMI C I NVERSI ON

The Dynamic I nversion is a technique developed for the plan of control systems for high performances, proposed especially for non-linear dynamic systems.

This is a methodology for the use of continuous- time dynamics in order to supply an estimate of the variable roots in the time of functions dependents from the time.

I n other words, the basis idea is to define a prior the wished functions of time on the outputs of the controlled system, in order to use,after, an inversion of the dynamic system to determine the inputs (m1, m2, m3) that generate the assigned outputs (, , , ’, ’, ’).

slide-19
SLIDE 19

The optimal solution, called reference, will be

  • btained through a feedforward branch.

Moreover, a feedback will become necessary to stabilize the dynamic around the same reference solution. Such feedback branch limits the sensors noise effects and effects of the non-modelling disturbance pair.

slide-20
SLIDE 20

Advantages of the D.I .

The technique of the Dynamic I nversion, from the point of view of continuous-time systems, allows to describe wished solutions that respect determined constraints and to describe the structure of such solutions.

The use of total stabilizing feedback in non-linear field requires the demonstration, preliminary matter, of the stability and can converge on non optimal solutions.

The Dynamic I nversion allows to realize a ring of feedforward inner the one of feedback, in order to this one become simpler, obtaining, normally, greatly advanced performances.

slide-21
SLIDE 21

D.I . Technique

Given a esprimibile system by means of:

) 1 ( ) ( ) , , ( T x B t x x f x   

  

I f  a constraint on its solutions:

) 2 ( , 1 ) , ( n i con t t x

i

    

) 3 ( ) , (           

 

t x x t x

i

) 4 ( 2 ) , (

2 2

                    

      

t x t x x x x x x x t x

j i j i i

And they are:

slide-22
SLIDE 22

Then, necessary and sufficient condition so that the constraint (2) is satisfied t, is that

  • ne dynamics such is had that the eq. (4) it is

verified t and that the eq. (3) and (2) they are satisfied in the initial conditions (t= 0).

t t x

i

   ) , (

t t x t x t x

i t i t i

      

    

) , ( ) , ( ) , (

slide-23
SLIDE 23

2 } ) ( ) , , ( { ) , , ( ) (

2 2

                                  

    

t x t x x x x x T x B t x x f x t x x G T x B x

j i j i

Where the last one is obtained replacing (1) in the (4):

 

t t x x G T x B x t x x n i con t x

t t

                             

   

) , , ( ) ( , , 1 ) , ( 

slide-24
SLIDE 24

Attitude Determination

  • The plan for the attitude control is guided from

the idea to make so that it is independent from the yaw angle (), that is, it is free of rotate round z axis (yaw axis) of the body reference of the BI SSAT, which is directed towards the Earth.

  • The problem has been dealed with considering

the equations that govern dynamics and the cinematic of the “rigid body”, simplifying therefore the modelling of the real system, even if can be present flexible elements (solar panels, antennas,…).

  • The rotation sequence of the Euler angles choice

is 3-1-2 (-- “Type 1”).

slide-25
SLIDE 25

I mplementation of the D.I .

slide-26
SLIDE 26

Final Result

In the hypotheses made previously (B = 1.5°; B =

0.8°) and for values of 0 = 0° and ’0 = 0 [°/ sec],

I have obtained following graphs relative to the Euler angles and the terrestrial magnetic field (necessary in

  • rder to use the magnetic-torquer):
slide-27
SLIDE 27

In order that, I shall use the torqrods, I estimated also the orders of the control torques and the necessary magnetic moments to produce demanded attitude.

slide-28
SLIDE 28
  • I searched “periodic solutions”, that is those

solutions for which the values of the Euler angles

and their derivatives, turn out to be equal, after the completion of one orbit, to the respective initial values (0= finale; 0= finale; 0= finale; ).

  • They are attractive since guarantee the same

behavior, of the satellite, for every orbit and therefore it is sure that the maximum magnetic moment made on one orbit is still the maximum after a sure number of orbits.

Periodic Solution

slide-29
SLIDE 29

Family of Periodic Solutions

  • It is possible in exact way to determine the periodic
  • solutions. They are those in correspondence of the

intersections of the curve of unitary value of ’final/’0 with that one of the number of the turns completed from the BISSAT.

  • -18°  0  0°

(101)

  • 0.111* vB  ’0  0.511* vB

(1 5)

slide-30
SLIDE 30
slide-31
SLIDE 31
  • In absolute, the periodic solution is not said to be the

best one, also introducing the mentioned advantages

  • already. In fact specific problems can exist in which

to consider a family of aperiodic solutions, with magnetic moment particularly low, is particularly convenient.

 I found approximate periodic solution, somewhat

interesting for values of:

0 = -10° ’0 = 0.311* vB

(vB = BISSAT angular speed)

slide-32
SLIDE 32

Max| m| = 3.547 [A·m2] min| m| = 1.882 [A·m2] B = 1.5° B = 0.8 [°/ sec] 0 = -10°

’0 = 0.311* vB

slide-33
SLIDE 33

Dimensioning and choice of the TORQRODS

Magnetic torquers are generally planar coils of uniform

wire (sometimes copper) rigidly placed along the attitude axes of the spacecraft. When a voltage is applied across a coil winding, a current is created, setting up a magnetic dipole. This dipole interacts with the Earth’s magnetic field, causing the coil to attempt to align its own magnetic field in a direction

  • pposite to that of the Earth’s.
slide-34
SLIDE 34

For the approximated periodic solution:

Max| m| = 3.547 [A·m2]

+ 50% of

  • ver-sizing

M = 5.32 [A·m2]

The chosen model is:

“TR10CFN”

that amply satisfies our requirements.

Resistance Scale Mass (2) Length
  • Dia. (1)
No. at 25°C (3) Factor (4) lbs in in
  • f
(ohms) (Am 2/mA) (kg) (cm) (cm) Coils (nom) (max) (max) (max) (max) Catalog Linear Saturation Residual Voltage Power Voltage Power Number Moment Moment Moment (V) (W) (V) (W) (min) (min) (max) (nom) (nom) (nom) (nom) 0.22 5 0,347222222 (0.10) (0.12) (1.27) 0.09 16.00 0.07 (0.4) (40.7) (1.8) 1.00 15.00 0.07 (0.5) (38.1) (1.8) 1.50 22.25 0,451388889 (0.7) (56.6) (1.7) 2.01 19.05 0.09 (1.0) (49.6) (2.3) 0,177083333 25.01.00 1.01 (1.7) (63.8) (2.8) 3.09 25.00.00 1.01 (1.8) (63.5) (2.8) 9.08 16.02 2.00 (4.5) (41.2) (5.1) 7.02 34.00.00 1.04 (3.3) (86.4) (3.6) 4.06 33.02.00 1.00 (2.1) (84.4) (2.1) 0,34375 29.42.00 1.04 (3.5) (75.0) (3.6) 8.00 38.31.00 1.26 (3.63) (0.98) (3.2) 10.09 37.01.00 1.06 (5.0) (94.3) (4.1) 5.52 35.43.00 1.00 (2.5) (90.0) (2.5) 6.01 36.00.00 1.00 (2.8) (92.0) (2.6) 7.07 30.25.00 1.03 (3.5) (76.9) (3.3) 7.07 30.25.00 1.03 (3.5) (76.9) (3.3) 11.35 36.04.00 1.04 (5.2) (92.5) (3.6) 9.06 36.00.00 1.03 (4.4) (91.5) (3.1) 12.25 48.00.00 1.02 (5.6) (122.0) (3.1) 17.00 50.25.00 1.03 (7.8) (127.7) (3.3) 27.00.00 44.00.00 2.01 (12.3) (111.8) (5.4) 30.00.00 55.01.00 1.06 (30.6) (140.0) (4.1)

x

Moments (Am 2) Linear Moment Saturation Moment Notes (Absolute Values)

TR1UPN

1.01 1.13 0.01 5.06 0.18 7.00 0.28 176 0.029 1 Caseless design for low
  • weight. (5)

TR10CFN

12 15 0.01 10.05 0.07 14.00 1.03 150 0.17 1 g case and 2 mounting blocks. TR10CFR 12 15 0.01 15.04 0.09 20.04 1.06 265 0.21 2 g case, no connector, (pigtail leads), TR30CFN 30 35 0.03 14.01 1.01 17.02 1.06 179 0.40 1 g case and 2 mounting blocks. TR30CFR 35 40 0.02 23.05 4.02 26.02.00 5.04 132 0.20 2 g case, 3 mounting blocks. TR60CFR 60 75 0.03 9.03 2.01 12.04 3.08 40 0.26 2 g case, 3 mounting blocks. TR65CAR 65 80 0.04 9.02 2.02 12.03 3.09 39 0.28 2 Aluminum
  • case. (5)
TR100CFN 110 130 0.02 10.09 5.09 13.05 8.09 20 0.21 1 g case, 2 mounting blocks. TR100CFR 100 130 1.00 8.04 0.07 12.02 1.04 106 1.20 2 g case, 3 mounting blocks. TR100UPR 110 135 1.00 20.00 2.05 25.06.00 4.01 162 0,063888889 2 Caseless design for low
  • weight. (5)
TR100CAR 110 130 1.00 14.06 1.07 18.05 2.07 126 0,065972222 2 Aluminum
  • case. (5)
TR130CFR 130 170 1.00 18.05 1.02 26 25 276 0,109722222

2x

g case, 3 mounting blocks. TR140CFR 140 170 1.00 9.05 0.08 12.02 1.04 110 0,088888889 2 g case, 3 mounting blocks. TR150UPR 150 180 1.00 24.00.00 4.05 30.02.00 7.03 124 0,056944444 2 Caseless design for low
  • weight. (5)
TR160UPR 160 190 1.02 19.02 3.01 24.02.00 5.00 118 1.00 2 Caseless design for low
  • weight. (5)
TR180CAR 180 220 1.00 22.05 7.05 29.05.00 12.08 68 0.55 2 Aluminum
  • case. (5)
TR180CFR 180 220 1.00 22.05 7.05 29.09.00 12.08 68 0.55

2

g case, 3 mounting blocks. TR230CFR 230 270 1.05 15.08 3.08 18.09 5.04 66 0,067361111

2

Caseless design for low
  • weight. (5)
TR230UPR 230 290 1.05 25.02.00 5.00 32.08.00 8.05 127 1.02 2 Caseless design for low
  • weight. (5)
TR300UPR 300 370 2.05 15.09 3.03 21.09 6.03 76 1.47 2 Caseless design for low
  • weight. (5)
TR460UPR 460 550 3.00 20.02 4.08 25.09.00 7.09 85 2.00 2 Caseless design for low
  • weight. (5)
TR500CFR 500 600 2.00 31.00.00 5.06 39.01.00 8.07 175 2.09 2 g case and 3 mounting blocks. TR790UPR 790 950 4.00 26.06.00 8.08 34.00.00 14.05 80 2.05 2 Caseless design for low
  • weight. (5)
slide-35
SLIDE 35

Conclusions

  • I t is possible to realize a system that

advantageously turns out economic, as regards traditional realizations, since they can be utilized:

  • One sun sensor and
  • One magnetometer,

in order to determine the attitude;

  • Three magnetic-torquers as

actuators, for the restoration of the attitude, one for every body axis.

slide-36
SLIDE 36

End of Presentation