Design and Implementation of Robust Decentralized Control Laws for ...

0 downloads 0 Views 2MB Size Report
Design and Implementation of. Robust Decentralized. Control. Laws for the ACES Structure at Marshall. Space Flight Center. Emmanuel. G. Collins,. Jr., Douglas.
NASA

Contractor

Design Robust Laws

Report

and Implementation of Decentralized Control for the ACES

at Marshall

Space

Emmanuel

Jr.,

and

G. Collins,

David

Harris

under

Florida

for Research Contract

Center NAS1-18872

National Aeronautics and Space Administration Office of Management Scientific and Technical Information Division

1990

Flight

Douglas

J. Phillips,

C. Hyland

Aerospace

Melbourne,

Langley

Structure

Corporation

Government

Prepared

4310

Systems

Division

Center

Abstract

Many mance

large space system

requirements

it is imperative based

vibration Space

that

Flight

especially

the benefits This

Investigator control

report

The

the study

accuracy.

of active

The

results The

structure

of line-of-sight

of this

be practically

conducted experiment

testbed

is the

is dynamically control

control

to satisfy

for these concepts

control

an experiment

structure.

ACES

active vibration In order

vibration

describes

Program.

for a flexible

Center.

allows

will require

such as line-of-sight

experiments.

CSI Guest

concepts

traceable

to become

successfully

as part

perfor-

operational

demonstrated

by Harris

ACES

critical

in groundof the NASA

demonstrate

structure

at NASA

to future

space

active Marshall

systems

and

issues.

iii PRECED]F,]G PAGE

BLAI',]Rix;OT FILMED

I_'_.L-L_INIENTIONALLY

8tANg

Table

1.0

Introduction

2.0

Description

3.0

Modeling

5.0

................................

of

the

ACES

Procedure

3.1

Choosing

3.2

Model

4.0

of Contents

the

Control

Structure

for

the

Modeling

Development

Design

1-1

ACES

the

the

ACES

4.1

Design

Process

for the

Detector

4.2

Design

Process

for the

Base

4.3

Design

Process

for the

LMED

Performance

of

the

Realization

Structure

..................

to IMC

Force

to Pulse

Disturbances

5.2

Response

Due

to to an

RCS

5.3

Response

Due

to a Crew

5.4

Some

Gimbal

4-1

Loops

to Accelerometer

4-1 ............

Loops

4-2

..........

4-3

.................

5-1

....................

Disturbances

Disturbances

5-1

..................

5-2

....................

Implementation

Results

5-2 .............

5-2

7.0

References

.................................

7-1

8.0

Appendix

.................................

8-1

BLANK

Conclusions

3-3

Closing

PAGE

and

.......

6.0

PRECEDING

Remarks

Algorithm

................

Controller

Due

on the

Loops

to AGS

Response

Remarks

3-1 3-1

Eigensystem

5.1

Final

...............

......................

Gyro

Integrated

2-1

Structure

Procedure Using

for

....................

NOT

FILMED

.....................

6-1

_L..._.._NIEN]IOttALL'_

BL_K

List of Tables 5.1.1

Det-X Response Improvement to a BET-X As the Controllers Were Combined

Pulse

................

5-3

5.1.2

Det-Y Response Improvement to a BET-Y As the Controllers Were Combined

Pulse

................

5-4

List 1.1

The Maximum Entropy/Optimal Design Equations

The ACES

Test Article

of Figures

Projection

..................

1-3

...........................

2-5

AGS-X to BGYRO-X: FE Model Bode Plot

Magnitude vs. FRF

Comparison

of

................

3-4

3.1.2

AGS-Y to BGYRO-Y: FE Model Bode Plot

Magnitude vs. FRF

Comparison

of

................

3-5

3.1.3

AGS-X to BGYRO-X: Phase Comparison FE Model Bode Plot vs. FRF

3.1.4

IMC-X to DET-Y: Magnitude Comparison FE Model Bode Plot vs. FRF

of

.................

3-7

3.1.5

IMC-Y to DET-X: Magnitude Comparison FE Model Bode Plot vs. FRF

of

.................

3-8

3.1.6

Pattern of Strong and Weak in the System Dynamics

Iterative

Procedure

Final Steps in Model AGS-X to BGYRO-X

for System

of

Interactions

3-6

....................

Identification

Development and AGS-Y

..................

3-9

..................

3-10

for ..................... to BGYRO-Y Loops

3-11

vii PRECEDING

PAGE

BLANK

NOT

FILMED

I__NRNTIONALL_

BLANK

3.2.3

Final Steps in Model Development for IMC-X to DET-Y and IMC-Y to DET-X

.................... Loops

3.2.4

AGS-X to BGYRO-X: Magnitude ERA Model Bode Plot vs. FRF

Comparison

of

.......

3.2.5

AGS-Y to BGYRO-Y: Magnitude ERA Model Bode Plot vs. FRF

Comparison

of

................

3.2.6

AGS-X to BGYRO-X: Phase Comparison ERA Model Bode Plot vs. FRF

3.2.7

.........

3-13

3-14

..................

3-15

IMC-X to DET-Y: Magnitude Comparison ERA Model Bode Plot vs. FRF

of ..................

3-16

3.2.8

IMC-Y to DET-X: ERA Model Bode

of ..................

3-17

4.1.1

Disturbance Tracking to DET Loops

..................

4-5

4.1.2

DET-X Response to a BET-X Pulse: ..................... Open Loop vs. Closed Loop with the IMC to DET

Controller

DET-Y Response to a BET-Y Pulse: ..................... Open Loop vs. Closed Loop with the IMC to DET

Controller

4.1.3

of

3-12

Magnitude Comparison Plot vs. FRF

Problem

for the IMC

4-6

4-7

4.2.1

Influence of Maximum Entropy Region for the AGS to BGYRO

Design on Phase Loops

4.2.2

Influence of Maximum Entropy Region for the AGS to BGYRO

Design on Magnitude Loops

4.2.3

Maximum Entropy Design Robustified the Notches for the High Frequency Modes

4.2.4

DET-X Response to a BET-X Pulse: ..................... Open Loop vs. Closed-Loop with the AGS to BGYRO

..o

VIII

in the Performance

in the Performance

...................

.......

4-8

.....

4-9

4-10

4-11 Controller

4.2.5

4.2.6

BGYRO-Y

Response

Open

vs.

DET-Y Open

4.2.7

4.3.3

4.3.4

Loop

Response Loop

vs.

to a BET-Y

Response

Open

Loop

vs.

Block

Diagram

for LMED

Feedback

Response vs.

Closed

ACCEL-1X

Response

Open

vs.

Closed

ACCEL-2X

Response

Open

vs.

Loop

Closed

4.3.5

BGYRO-X Open Loop

4.3.6

ACCEL-1Y

Response

Open

vs.

4.3.7

5.1.3

5.1.5

5.1.6

5.1.7

CIosed

Response

Open

vs.

Loop

DET-X

Closed

Response Loop

vs.

Response

Open

vs.

Response

Open

vs.

Closed

ACCEL-2X

Response

Open

vs.

DET-Y Open

Loop

Response Loop

with

vs.

Closed

Loop

with

Loop

with

with

Loop

with

Loop

with

with

Loop

with

the

the

4-16

Controller

the

LMED

4-20

4-21

Controller

..................... Integrated

5-5 Controller

...................

Integrated

5-8

Controller

..................... Integrated

5-7

Controller

...................

Integrated

5-6

Controller

...................

Integrated

ix

4-19

Controller

...................

Pulse: the

4-18

Controller

...................

Pulse: the

4-17

Controller

LMED

Pulse: the

Pulse: with

...................

...................

Pulse:

to a BET-X Loop

4-15

...................

Pulse: the

to a BET-X Loop

4-14 Controller

...................

LMED

Pulse: with

Controller

Pulse: ................... the LMED Controller

to a BET-Y Loop

to BGYRO

LMED

Pulse: the

4-13

to BGYRO

LMED

Pulse: the

4-12 Controller

...................

AGS

Pulse: the

to BGYRO

..................... AGS

Loops

to a BET-Y Loop

...................

AGS

Pulse: the

to a BET-X

to a BET-Y Closed

the

to a BET-X

to a BET-X

Closed

ACCEL-1X Loop

Loop

to a BET-X Closed

BGYRO-Y Loop

to a BET-X

Response to a BET-Y vs. Closed Loop with

ACCEL-2Y

Open 5.1.4

Loop

to a BET-Y with

Open

Pulse: the

Pulse: with

Closed-Loop

BGYRO-Y

Loop

with

Closed-Loop

BGYRO-X

Loop

to a BET-X

Closed-Loop

5-9 Controller

5.1.8

5.1.9

5.1.10

BGYRO-X

Response

Open

vs.

Loop

to a BET-Y

Closed

ACCEL-1Y

Response

Open

vs.

Loop

ACCEL-2Y Open Loop

RCS

Loop

the

to a BET-Y

Closed

Loop

Response vs. Closed

Disturbance

Pulse:

with

...................

Integrated

Pulse:

with

the

...................

Integrated

Profile

BGYRO-Y

Response

Open

Loop

vs.

Crew

Disturbance

DET-X Open

Response Loop

vs.

..........................

to a BET-X

Closed

Loop

Profile

with

RCS the

Loop

Crew with

the

5-12

5-13

Disturbance: Integrated

................ Controller

..............

5-14

5-15

Controller

..........................

to a BET-X Closed

5-11

Controller

to a BET-Y Pulse: ................... Loop with the Integrated Controller

DET-X Response to a BET-X RCS Disturbance: Open Loop vs. Closed Loop with the Integrated

5.2.3

5-10

Controller

5-16

Disturbance: Integrated

x

............... Controller

5-17

1. INTRODUCTION Many

large space system

performance

requirements

roughness.

In order vibration

discusses

an experiment

results

Figure

primary

reduction

of the

cantly

the beam

reduced-order

and

law implementation

The

Section

3, the strong

involving

actuators

eight

The jection

primary Approach

simultaneous controller

Entropy/Optimal high performance, (Figure

1.1) consist

equations MEOP

interaction

two of the feedback

paper

Program.

space

missions,

testbed,

shown

in

that the

provide

limitations

substantial

controllers

to signifi-

on controller

simplicity

reduction

in control

on throughput

of even

design

at NASA

by designing

design

[12-13].

by using

As will be described

in the dynamics

decentralized

was facilitated

of the ACES

control

more

fully in

structure

laws for selected

the

allowed

system

loops

Optimal

Pro-

and eight sensors.

for Uncertain

Systems

for control

robustness

Projection robust

versus (MEOP)

control

of four coupled

control

design

(OPUS)

of four fundamental

was used to develop approach

benefits

This

Investigator

Complexity

for control

patterns

in this experiment

[1-10].

The

issues in control system

law design matrix

equations

design

actuator

A subset

[2-6] which

for flexible

OPUS

design:

performance.

methodology

when the plant is known perfectly

projection

the

(MSFC).

was also placed

architectures).

used

developed

chosen

and

Center

required

due to stringent

models

methodology

trade-off order

interest

system

and weak

the control

emphasis

that

The

controllers

of this objective

controller

Algorithm

us to approach

to design

surface

to future

structure.

Flight

on rms

processors.

of the

Realization

was

Space

critical

experiments.

traceable

for a flexible

Particular

decentralized

qualified

development

Eigensystem

vibration.

it is imperative

directly

to satisfy

constraints

in ground-based

Marshall

satisfaction

is of paramount

space

control

at NASA

The

operational

features

of this experiment

LOS errors.

and

as part of the NASA CSI Guest

has

vibration

control systems

accuracy

demonstrated

which

structure

(LOS)

to become

by Harris

active

attenuate

state-of-the-art

concepts

experiment,

objective

active vibration

as line-of-sight

conducted

2.1, is the ACES

will require

be practically

demonstrate

The

(i.e.,

control

of this

successfully

such

for these

of active

The

concepts

is the process

sizing, sensor

of OPUS

was developed

structures.

The

is the

MEOP

design

specialize

to the standard

and a full order

controller

is desired.

[7-8] was used to develop

vibration

control

loops. 1-1

suppression.

laws for a tracking

accuracy, Maximum

particularly

which

laws for active

allows for the

to allow equations

LQG

In this experiment

The discrete-time problem

Riccati

associated

optimal with

The paperis organizedas follows.Section2 describesthe basicACES configuration.Section 3 then discusses the finite elementmodelprovidedby MSFCand the developmentof the models actually usedfor control design. This sectionalsomotivatesthe decentralizedapproachchosen for the control design. Next, Section4 describescontrol designfor eight systemloops which wereselectedby analysisof test data collectedfrom the structure. Section5 then presentsthe experimentalresults. It is seenthat very significantperformanceimprovementis achievedwhen the eight feedbackloopsareclosed.Finally,Section6 presentsclosingremarksandconclusions.

1-2

0=

A_Q

+ QA T + Vt -

QZ'Q

+ r±Qy_,Qr

T

P

P

+ Z

a2AiQAT

+ Z

i:I

0 = ATp

a2AiQAT

i=I

+ PA T + R1 - PZP

+ rWp_Pr±

P

P

+ _ _ ATPA,+ _ 4 ATPA, i=i

0 = (A_ -

_P)Q

+ Q(A_

0 = (A_ - Q2)TP

_"= c2P(_P)

Figure coupled

1.1 The MEOP by a projection

flexible

structures.

-

i=1

Z'P)

+ Q_,Q

-

+ P(A_ - QF_,) + P£P P = rank

QP

T

- rTpZPr±

rank

_) :

#,

(.)# denotes the group generalized

design matrix

rank

r±QZ'Qr

:

nc

inverse.

equations consist of four modified Riccati and Lyapunov equations, r and allow high performance, robust control law developement for

1-3

2.

DESCI_IPTION The

ACES

ployable,

built

(JPL). The

form

The

approximately

of the beam

Inc.

It was supplied

which

is triangular

properties.

When

continuously

the beam

basic

test

article

is a de-

is a spare

Voyager

and is very lightly

damped.

section.

its full length. each having

the Astromast

article

Labora-

in cross

along

test

by the Jet Propulsion

5 pounds)

into 91 sections

fully deployed,

The

to MSFC

(about

beam

divide

The

lightweight

and extend

its shape,

MSFC.

45 feet in length.

is extremely

configuration

the Astromast appendages

(Figure 2.1) consists

tip and the pointing creates

the

frequency

Research,

at NASA

Three

The

cross members,

equal

exhibits

longerons

length

and mass

a longitudinal

twist

of

260 degrees.

The ACES

Overall,

is located

is a symmetric

elastic

STRUCTURE.*

approximately

Astromast

give the beam

and similar

the

"nested"

structure

modes

configuration study

beam,

Astromast

ACES

testbed

by ASTRO

the corners

which

THE

experimental

lightweight

Astromast tory

OF

is very

(more

modal

and

40 modes traceable

and counterweight

arms at the Astromast

frequencies

flexible

than

is dynamically

gimbal

of an antenna

characteristic

lightly

under

damped.

10 Hz).

to future

space

base.

The addition

of Large

It contains

Space many

As illustrated

systems

legs appended of structural

Structures closely

by Figure

and is particularly

to

(LSS).

spaced,

2.1, the

low

ACES

responsive

to the

of LOS issues. The precise

voltage

input

Variable servo

motion

to the BET

Differential

controllers

adjust

of the Base Excitation servo control

Transformer

compare

the position

system. whose

the commanded

input

of the BET.

The closed-loop

limitations

of the hydraulic

to be position

commands

to the BET.

two pointing

Motion

gimbals,

Figure

1.1 shows

control

design

gimbals

the

Compensation

location

are each

* This description

positioned

detector

of each

of the

the laser beam on the

of the ACES

movements

are monitored

LVDT

consists

and associated components

in the center

a commanded by a Linear

signals

and

the disturbances

of a 5-mW electronics, of the

primarily

from

and two power

The detector

to increase

[11]

The

the

within

are chosen

two 12-inch

IMC system.

of the detector. appendage

laser,

The

automatically

allows any type of BET movement

In this experiment

System

2-1

by supplying

to the

controller

is taken

is obtained

are fed back to the servo controllers.

voltage

end of a flexible

testbed

BET

outputs

system.

(IMC)

a four quadrant

is to position

The

(LVDT)

the frequency

The Image

Table (BET)

mirrors, supplies.

goal of the and pointing

difficulty

of the

control problem. The the controller

design

lack of information (i.e., there

about

the

appendage

is no accelerometer

motion

also adds

or gyro at the location

complexity

of the

gimbals

to

or the

detector). In addition Gimbal

to the two IMC gimbals,

System

applications, provides

(AGS),

which

torque

control

of 27 amps servo

which

outputs

can generate

azimuth

torquer

however,

torquers,

be set manually

of azimuth

freedom.

of-10

that

designed

gimbal

in the

azimuth.

AGS

receives

Astromast.

The

to +10 volts.

causes

of generating

13.8 ft-lbs

allows

the

5 degrees test

range

torque

from

system)

the

in the

a current

limit

the AGS

command

the

signals.

used in the SSC labora-

range

± 30 degrees.

of _ 5 degrees.

at any position to any

system

to the current,

of approximately

to be rotated

gimbal

Because

proportional

over an angular

of rotation

article

measure.

to produce

pointing

commands

represents

and servo amplifiers

over an angular

The

This saturation

torque

the Advanced

accuracy

acquisition

be designed

of torque

i

for high

data

an applied

must

also include

as a protective

with the power supply

to allow

This

system

servo amplifier

in the COSMEC

is capable

actuators

a third

of the

into the AGS

37.5 ft-lbs

control

on an HP 9000 via the COSMEC

a current

used

gimbal

with

base

over the range

The AGS gimbal tory,

at the

is built

algorithms

two-axis

augmented

(implemented

inputs

amplifier

control

has been

actuation

algorithm

form of analog

a precision,

the available

about

The It can,

the 360 degrees

position

desired

without

remounting. Linear apply

Momentum

forces

having

and measure

orthogonal

(LVDT's). These

The

axes,

were

modes.

Y axes of the

inertial

to the structure

and

(LMEDs)

packages

and

the Astromast

a colocated

Each

frame

the

accelerometer

sensor/actuator

LMED

package

Variable

at intermediate

ability

is at rest, shown

colocated

two Linear

are positioned

to maximize

reference

provide

accelerations.

two accelerometers,

selected

When

Devices

the resulting

two LMED

locations

structural

Exchange

of these

in Figure measures

2.1.

points

package

The

along

Transformers the

is aligned applies

which

two LMEDs

to control

LMED

the resulting

contains

Differential

devices

each LMED

pairs

Astromast.

the

dominant

with the X and a horizonal

acceleration

force

at the actuator

location. The LMED applied travels

is linear permanent

to the structure along

coil which

a single

extends

as a reaction shaft

inside

magnet

motor

against

on a pair of linear the magnet

assembly

whose

magnet

the acceleration bearings. from

2-2

functions

of proof mass.

The armature one end.

as a proof mass.

The

The

magnet

of the motor

magnet

assembly

Force is assembly

is of a hollow moves

along

the shaft on each

with respect end

by a bracket

accelerometer the proof

which

is mounted

mass with

In addition and LVDT's gyros

to the coil which is fixed to the LMED

in line with

respect

rate

package

small

angular

digital

converter

period

of approximately

card

rates The

dynamic

for warmup, the

power

included

on board

The signals

time

reduces

the instrument

package.

One channel

of each pair

to that

and negative

the

As in the

cases

COSMEC

system.

As mentioned

to 12-bit ATM

requires

rate

Rate

and the

of

accelerometers

include

Gyros.

signals binary

three-axis

rate

However,

since

and

They

were not

words

during

are designed

of an ATM

gyro packages

1.5 amps

ASTROMAST provide

requires

are

rate

gyro

by the

analog-to-

require

a warmup

warmup

0.9 amp

measured

resolution

1.2 amps

per package.

are different

require

and then

1.25

The

by two identical

finer than

0.0001

approximately

at 28 Volts

from the ATM

of the accelerometer carries

by a decrease

previously

position

DC.

20 minutes After

accelerometer

g and a

warmup

electronics

are

Two channels

are

package.

acceleration

of the other

A linear

devices.

output

of 25 to 30 Hz. They

package

of the synchronization

acceleration

the

implementation

(ATM)

analog

The

tip of the

from the accelerometers

carries

Mount

package

to about

of freedom

identical

Each

each

for each degree

channel

mass.

at the tip and base.

measurement

The accelerometers

required

other

is constrained

to measure

devices

for controller

is converted

of ± 3 g with a bandwidth

requirement

proof

IMC System

measurement

The

system.

and

packages.

which

magnet

at 28 volts DC.

at the base

during

to the

is utilized

accelerometers

Telescope

signal

COSMEC

both

accelerometer range

the available

precisely.

40 minutes.

accelerations

three-axis

very

force

with the

only the remaining

analog

of the

amps after stabilization, The

associated

at the base are Apollo

is ± 45 volts.

An LVDT

at the tip are not available

we will describe

The rate gyros to measure

detector

centering

The

assembly.

as well as three-axis

gyros

used for evaluation,

the shaft.

with the LMEDs,

at the tip and base

the three-axis

a small

to the LMED

to the two-axis

associated

provides

package.

a 2.4-kHz

information. channel,

the computer

these

square Zero

positive

in frequency

instruments,

package,

system

2-3

i.e., six channels

acceleration acceleration

signal,

is represented by an increase

by a hardware

of an HP 9000 digital

and the

by a signal in frequency,

to the synchronization

are monitored

consists

per accelerometer

wave synchronization

as compared

signals

rate gyros.

channel. card

in the

computer

inter-

facedwith the COSMECInput/Output system.The HP 9000performsthe controlalgorithm,data storage,real-timeplotting, andthe strapdownalgorithm(describedin the next section).The HP 9000is a 32-bitmachinewith an 18-MHzclockrate. It includesan HPIB interfacecard,two 16-bit parallelinterfacecards,512kbytesof extra memory,anda floppydiscdrive. The benchmarktest timesfor processingthe presentcontrol andstrapdownalgorithms,plotting, andstorageare .010 to .013millisecondsper sample. The COSMECis a highlymodifiedAIM-65microcomputersystemusedforI/O processing.The primary purposesof the COSMECare to processthe sensorinputs, to provideforce and torque commandsfor the actuators,and to off-loadcontrol and sensordata to the computersystem. Currently, the COSMECperformsthesetaskswith 25 sensorinputs and nine actuator outputs, whilemaintaininga 50-Hzsamplingrate. Thecycletime for COSMECoperationis approximately 5 milliseconds. In our controldesignandimplementationweused8 controlinputsand8 measurement outputs. The inputswerethe X andY torquesof the IMC gimbals,the X andY torquesof the AGSgimbals andthe X and Y forcesof the two LMED packages.The measurements consistedof the X and Y detector(DET) positionoutputs,the X andY basegyro (BGYRO)rate outputsandthe X andY outputs of the LMED accelerometers.

2-4

/

o

=,.2.

N_

© g m_ _._

_.-

2-5

_

3.

the

MODELING

PROCEDURE

The first critical

step in control

process

that

procedure

that

the system 3.1

identification

Initial

designer

has

model

development

for flexible

model

design

and validation.

models

Realization

(FEM)

structures

by using

Algorithm.

that

the

three data

FEM

Below, we describe

a system

We then

identification

present

details

data

collected

(i) modify

an appropriate

which is in some sense a hybrid

based on the finite element

by comparing

is inadequate

options:

using

is usually

is evaluated

from input-output

essentially test

STRUCTURE

development

control

the Eigensystem

Procedure

shows

input-output

our

Modeling

gleaned

comparison

ACES

of

procedure.

finite element

information

this

the

the

model

The initial with

upon

THE

design is model

led us to of develop is based

Choosing

FOR

from the actual

FEM,

system

of the FEM

its time and frequency

for control the

system

identification developed

physical

design,

(ii) develop

and models

approach.

apparatus.

then

a new

algorithm

responses

the

model

If

control based

on

or (iii) develop

using system

a

identification

techniques. In this noise,

experiment

sine-sweeps

system

modes

to limitiations inputs

the

inputs

or delta-functions.

were

on the length

system

finite element DET-Y,

loops

model.

(ii) IMC-Y

As evidenced ent frequency

contain

1.4 Hz mode

the FEM

predicts

trends

than

FRF's

that

yielded

which

of Figures

lower frequency to generate

real behavior

that

the broad-band in the

response

(iv) AGS-Y

the FEM

finite

bode

TO

loop.

for the AGS-X

3-1

Also notice

to BGYRO-X

(FRF's)

plots

of the

(i) IMC-X

to

BGYRO-Y. differ-

For example,

Figure

model

that

time

significantly

the open loop LOS performance.

to BGYRO-Y

due

random

functions

predicted

to some of the modes

the

allowable

loops:

the structure. element

about

this is primarily

for four system and

random

information

dynamics

frequency

testing

loop the

peaks corresponding

for the AGS-Y

most

broad-band

with the corresponding

3.1.1-3.1.5,

by actually

influenced

the

and the fact that

to BGYRO-X,

to BGYRO-X most

allowed

show this comparison

(iii) AGS-X

to be either

We conjecture

were compared

those obtained

AGS-X

positive

inputs

was used

3.1.1-3.1.5

do not show the large magnitude shows analogous

data

and these

for the

chosen

the dominant

by the comparisons

that

the

test

to DET-X,

3.1.1 shows

were

and delta-functions.

excite

Figures

responses

system

of the time-histories

The input-output

of selected

The

the sine-sweeps

did not significantly

windows.

to the

past

Bode

plot does

The

FEM's

not also

8 Hz. Figure

3.1.2

as seen by Figure

3.1.3

loop, while the

FRF reveals

that between2 Hz and4Hz the phaselags-90° by as muchas 25° (evenwhen the computational delaydueto the 50 Hz samplerate is not takeninto account).This phaselag is probablydueto actuatorandsensordynamics. The FRF's ofthe IMC-X to DET-Y andIMC-Y to DET-X loops,shownrespectivelyin Figures 3.1.4and3.1.5,revealthat theseloopsareinfluencedverylittle by the flexiblemodesof the structure. It followsthat the IMC gimbalsarenot capableof controllingflexiblemodesto improveLOS performance.Thus,if oneconsidersthe four actuatorinputs(IMC-X, IMC-Y, AGS-XandAGS-Y) andthe four sensoroutputs(DET-X, DET-Y, BGYRO-X,BGYRO-Y), it is not necessary to feed backthe BGYRO outputs about

the behavior

outputs

of the flexible

do not contain

much

provided

by the BGYRO's.

the DET

outputs

by Figure above DET-Y,

were

and

above

the

loops

comparable

analysis FRF's

four

IMC-Y

In summary,

and accounting

should

account

space

of their

of the LMED's

data

cannot

outputs

contain

information

In addition,

gimbals

that

be improved

the decentralized AGS-Y

the DET

is not already by feeding

to the AGS gimbals.

to BGYRO-X,

performance finite

back

As illustrated

structure

described

to BGYRO-Y,

models

IMC-X

and actuator

and actuator

was in transition.

useful

system

actual

dynamics

sensors

loops. and

control some

to

test data.

described

four dominant

structure.

Also,

of the trends design

design.

Thus,

seen

studies,

the

in the

the

FEM

it was necessary

Modifying

to

the finite element

would have been a very time-consuming

to use the Eigensystem

Realization

It is important

actuators,

models

Algorithm

to note that generated

(ERA)

since all test

from

this data

dynamics.

input-output

In addition,

showed

involving

for preliminary

control

a model based upon

and actuators

structure

a centralized

and was thus adequate

the control

masses.

control

model

of the four dominant

obtaining

for the four sensors

with

element

and thus we decided

for the sensor

proof

that

for high performance

by using

We had difficulty itations

control.

for the AGS

that within

AGS-X

primarily

cannot

performance

also revealed

the

for the sensor

process

were collected

information

with a decentralized

the FEM or develop

state

IMC's

or the BGYRO

achievable

by the test

model

the

of the test data revealed

although

as a model

and expensive

useful

loops:

performance

that

modify

data

dominant

which

the achievable

of test data

to the

either

to develop

Thus,

analysis

revealed

was inadequate

(if any)

since the BGYRO's

to DET-X.

achievable

generated

modes

to the IMC gimbals

3.1.6, analysis

there

to the IMC gimbals

early

In particular,

data

from the

in the project MSFC

3-2

in joint

LMED's

the internal consultation

due to the stroke control

lim-

configuration

with the guest

inves-

tigatorsdecidedto havethe springsremovedfrom the LMED's andreplacedby internal position loops.The uncertaintyregardingthe dynamicsof the LMED's ledus to delaycontroldesignactivities involving thesedevices.Ultimately,it wasdecidedto simplyfeedbackthe colocatedLMED accelerometers to the corresponding LMED forceaxisandto determinesimpledynamicsfor these controllersby usingcrudemodelsof the loops,developedby informationprovidedto us by MSFC andour knowledgeof proofmassdevices.Aswill be seenin Section5,thesecolocatedLMED loops did providesignificantperformanceimprovement. 3.2

Model

Development

The collection usually steps

procedure

in developing

design

17th and DET-Y

the

for the

19th order and IMC-Y

implementation

filters

AGS-X

were effective

loops.

to BGYRO-X models.

model.

because

Figures

for control

first-order

all-pass

for the computational

delay.

AGS-Y

The control

discrete-time to design

maturity

to BGYRO-Y

design

models

design

design

models.

The

were

control

were respectively

Ideally,

directly

using

the AGS to BGYRO

in the theory

is

the final

filters

for the loops

the controller

We chose to design

of a greater

models

3.2.2 and 3.2.3 show

Note that

and

it is better

Algorithm

high fidelity

3.2.1.

loops to account

algorithm

models

IMC-X

for digitial discrete-

feedback

and software

to

loops

development

setting. show

from

the ERA

by Figure

were both 4th order

of the system

3.2.4-3.2.8

generated

show that

loops

of a control

for the continuous-time

FRF's

to BGYRO

Realization

to obtain

for the four major

to DET-X

using continuous-time

Figures

of test data

continuous-time

time representations

Eigensystem

as illustrated

models

into the AGS

models

the

and manipulation

an iterative

incorporated

Using

comparisons

the test data.

models

closely

in emulating

The resemble

of the

ERA

magnitude

plots

the FRF's.

the computational

3-3

models

of the

of Figures

As illustrated

delay

four

system

3.2.4,

3.2.5,

by Figure

in the system.

loops

with

the

3.2.7, and 3.2.8

3.2.6, the all-pass

AGSoX 10 0

TO BGYRO-X

....................

10-1

FIE MODEL

10-2 Z 10-3

10-4

......

AGS-X

10o

_02

TO BGYRO-X

FRF 10"i

ud

I

.j'

10 "4

W"

il',,t t l,( ti

i

E -i I

I t,-

10.5 I 10-2

,

10-1

10o FREQ

Figure

3.1.1

BGYRO-X (which sponding

A comparison

loop most

shows

influences

to some

of the

that

of the the

LOS higher

finite

,

FRF

data

element

performance) frequency

and model

and

does

modes.

3-4

h

10 i

10 2

IN HZ

finite

element

neglects not

show

the

Bode contribution

the

large

plot

for

the

of the magnitude

AGS-X 1.4

to

Hz mode

peaks

corre-

10 0 _

AGS-Y

10-I

FINITE

TO BGYRO-Y

ELEMENT

MODEL

10-4

10-5_ 10-2

.............. 10-:

...............101 FREQ IN HZ

AGS-Y

TO BGYRO-Y

10 o FRF lff t

10-2

10-3

10-4 ¸

......

10-2

'10_l FREQ

Figure

3.1.2

BGYRO-Y (which sponding

A loop

most

comparison shows

influences

to some

of the

that

of the the

LOS higher

finite

......

101

FRF

data

element

and model

performance)

and

frequency

modes.

does

3-5

'102

IN HZ

finite

element

neglects not

show

the

Bode contribution

the

large

plot

for of the

magnitude

the

AGS-Y

to

1.7 Hz

mode

peaks

corre-

AGS-X 2OO

,

TO BGYRO-X ,

v

,,v.

,

i

,i,.,

150 100

m

50 0 -50

a.

-I00

-150

FE MODEL

"2_2

FREQ

AGS-X

......

_02

......

102

IN HZ

TO BGYRO-X

200 150 100 u2 t_

5O 0 -50 -100

......................................................

-150

-201! .2

10-1

10 ! FREQ

Figure behavior

3.1.3 while

For the

the AGS-X FRF reveals

to BGYRO-X that between

loop 2 Hz

25 ° .

3-6

the and

IN HZ

finite 4 Hz

element loop predicts positive real the phase lags -90 ° by as much as

IMC-X 10 3

TO

DET-Y

..................

10 2 gd

FE

MODEL

10 x

ERE

_

x_

FREQ

Figure

3.1.4

A comparison

of the FRF data

loop shows that the finite element and the dominance of the IMC-X much

higher

and finite

IN HZ

element

Bode plot for the IMC-X

to DET-Y

model correctly predicts the small influence of the flexible modes mode but predicts much lower damping in the gimbal mode and

loop gain.

3-7

IMC-Y 10 3

TO

DET-X

i

i

t

t

!

t

A

I

I

I

L

I

J

i

A

I

10 2

101

10 0

10-1 1(_ 2

i

I

I

A

I

I

i

I

I

10-1

i

I

I

I

I

I

I

10 o FREQ

IN HZ

Figure 3.1.5 A comparison of the FRF data and finite element Bode plot for the IMC-Y to DET-X loop shows that the finite element model over estimates the influence of the flexible modes.

3-8

l

101

0

._

E

e_

Q

o_

o_

0

U

3-9

r,D

C

0 Wla_ C _0

"U

c.D 0

0

u_

0




O'9 cW

rru.l

O_

C_

3-10

(I)

AGS-X

TO

13 ERA

BRATE-X

STATES

-1 STATE

LOST

+4 STATES

FOR

+1 STATE

FOR

17 th ORDER

(IV)

AGS-Y

IN CONVERTING HIGHER ALL-PASS

FROM

FREQUENCY

DISCRETE-TIME UNMODELED

TO EMULATE

CONTINUOUS-TIME

TO MODES

COMPUTATIONAL

DESIGN

CONTINUOUS-TIME

DELAY

MODEL

TO BRATE-Y

17 ERA

STATES

-1 STATE

LOST

-2 STATES +4 STATES +1 STATE 19 th ORDER

FOR FOR FOR

IN CONVERTING DELETED HIGHER ALL-PASS

HIGH

FROM

DISCRETE-TIME

FREQUENCY

FREQUENCY

CONTINUOUS-TIME

MODES

COMPUTATIONAL

DESIGN

Figure 3.2.2 The final steps in developing control and AGS-Y to BGYRO-Y loops yielded respectively

3-11

CONTINUOUS-TIME

MODE

UNMODELED

TO EMULATE

TO

DELAY

MODEL

design models for the AGS-X to BGYRO-X 17th and 19th order continuous-time models.

(I)

IMC-X 6 ERA

TO DET-Y STATES

-4 SPURIOUS +1 DELAY

STATES STATE

+1 FILTER 4th ORDER

(II)

IMC-Y USED

(DISCRETE-TIME)

DISTURBANCE

STATE

DISCRETE-TIME

DESIGN

MODEL

TO DET-X IMC-X

TO

DET-Y

DESIGN

MODEL

(THE OPEN LOOP GAIN AND THE DAMPING OF THE DOMINANT MODE WAS MODIFIED, HOWEVER)

Figure IMC-Y

3.2.3 The final to DET-X loops

steps in developing control design models yielded 4th order discrete-time models.

3-12

for the

IMC-X

to DET-Y

and

AGS-X

TO BGYRO-X

10 o

10-I

ERA

MODEL

\

10-3

10 4

10-5 10-2

............... 10-1

i_ FR_Q

AGS-X

............... 101

io_

IN HZ

TO BGYRO-X

10o FR_ 10-I

t

10.2

10-3

10-4

10-5 10-2

10-1

10o FREQ

Figure generated

3.2.4 from

The

ERA

test

data.

model

for

the

AGS-X

x

IN HZ

to

3-13

BGYRO-X

loop

closely

resembles

the

FRF

AGS-Y 100

TO BGYRO-Y

....................

10-I

ERA

MODEL

t_ t

0.2

_

10.3

10 4

1@5 10-a

,

10-1

10 o FREQ

AGS-Y

10 o -

l

,

i

101

L

,

10 a

IN HZ

TO

BGYRO-Y

10-1

10.2

10-3

10-4

10-5 10. 2

...... FREQ

Figure generated

3.2.5 from

The

ERA

test

data.

model

for

the

AGS-Y

_0 2

IN HZ

to

3-14

BGYRO-Y

loop

closely

resembles

the

FRF

AGS-X 2OO

,

TO

,

,

,

BGYRO-X

,,,

,

,

,

,

,,,

150 100 50 0 U3

/ i

-50 J -100 ERA

-150

MODEL

10-1

°2Op(y2

10 o FREQ

AGS-X

101

lo 2

IN HZ

TO BGYRO-X

200 150

50

io

1

-50

-_oo_

............................................ '_

"150 t -2

010_

2

10.1

i

,

_

,

10 o

,

,_1

FREQ

Figure

3.2.6.

As

shown

modeled

the

computational

here delay

for

the

AGS-X

by using

to

all-pass

3-15

IN HZ

BGYRO-X filters.

loop,

the

ERA

models

effectively

IMC-X

TO

DET-Y

101 i

i

i

i

i

i

i

i

i

i

i

/,?-_

.J

ERA

MODEI_ .............................

"

FRF

10°

10-1 10-2 FREQ

Figure 3.2.7 The ERA from test data.

model

for the IMC-X

IN I--IZ

to DET-Y

3-16

loop closely

resembles

the FRF

generated

IMC-Y

TO

DET-X

10 z

Z

101 _

10 o 10-2

I

I

t

I

f

I

I

I

I

I

[

L

ERA

model

for the

I

i

I

I

I

t

t

I

l

L

I

10 0 FREQ

Figure 3.2.8 The from test data.

I

10-1

IMC-Y

3-17

loop

closely

1

]

101

IN HZ

to DET-X

I

resembles

the

FRF

generated

4.

CONTROL

DESIGN

Once

settled

we had

structure

of the

a three

step

IMC gimbal

(ii)

design

of the AGS gimbal

(iii)

design

of the LMED

Below,

we give details

scribing

the

resultant

THE

design

integrated

delayed

Section

5. The sample

4.1

Design The

Process

design

essentially

processes

identical.

as a disturbance white

noise process

as the

positive

parameter

the low frequency synthesized

IMC-X

a filter

(i.e., modes

[7-8] by minimizing J(e)

the

along

to DET-Y

loop

with

data

controller

Loops

with

ACES

the

paths

IMC-Y

1/(z-

to be tracked

1-{-e), which

zero.

This filter accounted

below

1.5 Hz).

the quadratic

Optimal

loops

closed)

is

was 50 Hz.

to DET-X problem

loop were

was formulated by filtering

a discrete-time

for the system

Projection

de-

performance

w_ was modeled

approximates

data

the

all feedback

feedback

and

experimental

describing

4.1.1, for each loop the control

The disturbance

e approaches

modes

step

Experimental

to IMC

by Figure

problem.

wl through

for each

(i.e.,

Detector

As illustrated

for the

loops.

rate for each of the three

for the

tracking

procedures

controller

for the

design

and

accelerometer

improvement.

of the

control

loops,

to base gyro loops,

performance

architecture,

process:

force to colocated

of the

STRUCTURE

controller

to detector

improvement until

ACES

on a decentralized

was essentially

(_) design

FOR

control

a

integrator

biases

and

laws H(z)

also were

cost function

= lim E[qW(e)q(e)-{-

puW(e)u(e)],

p > 0.

/c---* oo

The

controllers

trollers

implemented

contained

integrators

both loops the design controller

gains

models

were the limiting which

were 4th order

are given in Appendix

The resultant

performance responses

to the x and

y axes

of the Base

eliminate

system

biases

as e --_ 0. Thus,

to effectively

eliminate

while the controllers

the

the implemented line-of-sight

implemented

con-

biases.

For

were 3rd order.

The

A.

improvement

loop and closed-loop

the

were able

controllers

is illustrated

by Figures

of the x and y axes of the detectors Excitation

and also improved

Table.

Notice

that

respectively the

the LOS performance

tracking. 4-1

4.1.2-4.1.3

feedback

which

show open-

to pulse

commands

loops

by providing

were

able

low frequency

to

4.2

Design

Process

The design quently

processes

to be referred

design

Entropy

of varying

only the modes

The Maximum

controllers

which

yielded

the ME approach Figure

varies

over

loop transfer

and were unstable region

robustness

Figure full-order

4.2.2

than

implication

controllers.

is that

of the compensator

two highest by Figure

authority

frequency 4.2.5,

implying several

when implemented.

the

in the

those

In practice,

The higher

interval,

ME design

shown was able

were then

In our control

designs

LOS perfor-

in developing

implemented.

on the

The phase

the

times.

stable

The utility

of

of a full-order

LQG compensator

plot of the expect,

the ME designs

Thus

phase

of the

Nyquist

As one would

However

region.

yield

Notice

corresponding

these

became

the ME designs

notched in the

design that

designs

positive

provided

robust

controllers actually

on the

were real in

the needed

performance

that

provide

insight

reduced

the high frequency

modes

models

to robustify

4-2

of Figures

the controller

magnitudes

reduced

order

into the choice

of the

controllers.

had high gain, i.e., the

3.2.4 and notches.

of a

An-

order that

shape

robustness.

are effectively

aid in synthesizing

ERA

magnitude

the ME compensator

thus providing

ME designs

and is a numerical

modes

that

compensators,

the ME designs

controllers

continuous-

controllers

to be crucial

design

of ME uncertainty

performance

the full-order

that

region.

influence

of the LQG

the modes

the open-loop

when

3 Hz).

rate feedback.

in the performance describes

improvement

(i.e., less than

toward

These

dominate

proved

influenced

4.2.1-4.2.3.

the origin

tending

compensator

are smoother

order

encircles

most

for the

mode.

prewarping.

of ME uncertainty

region

that

for the Y-loop

and robustness.

design

performance

modes

3, the

model

was used to synthesize

with frequency

robustness

influence

this frequency

function

the performance stability

the

the

(subse-

in Section

while the design

Similarly,

synthesis

authorities

by Figures

in the performance

widely

nonrobust

describes

model

3 Hz since these modes

(ME)

significant

modes.

to BGYRO-Y

As mentioned

For the X-loop

(MEOP)

control

less than

is illustrated

4.2.1

compensator

continuous-time

transformation

Entropy

loop and the AGS-Y

1.7 Hz and 2.3 Hz bending

Projection

orders,

Loops

were very similar.

model.

were

by using the bilinear

we penalized mance.

continuous-time

Optimal

Gimbal

to BGYRO-X

was a 17th order

LOS performance

time controllers discretized

to AGS

were 1.4 Hz and 2.4 Hz bending

most influenced Maximum

Gyro

to as the X and Y loops)

was a 19th order

LOS performance

Base

for the AGS-X

model for the X-loop

Y-loop

other

for the

3.2.5.

That

As illustrated

is, the controller

notcheswereincreasedin both width anddepth. The controllerswhich yieldedthe best performancewhenimplementedwerea 4th order controller for the X-loop and a 6th ordercontrollerfor the Y-loop. The controllergainsarepresented in AppendixA. The resultantperformanceimprovementis shownin Figures4.2.4--4.2.7 which showopenand closedloop responses of the detectorsandbasegyrosto pulsecommandsto the x andy axesof the BaseExcitation Table.Noticethat significantperformanceimprovementwasachievedin both the detectorandbasegyroresponses 4.3

Design

Process

of the

In this subsection to as LMED-1 The

control

and hereafter

while

design

the corresponding

LMED

the

two-axis

was based LMED

Force

the two-axis LMED

proof-mass

axis.

(ii) LMED-1Y

LMED-2Y

to ACCEL-2Y.

It was assumed

assumed

dynamics

From LMED

4.3.1

controller

it follows

and along

that

device

closest

to the

Thus

closest

to the base will be referred

tip will be referred

that H(s)

the open

loops is shown

the transfer

for control

(iii) LMED-2X

of each

in each loop. in Figure

function

from

a given axis) to the force applied

outputs

design

to ACCEL-2X,

loop dynamics

can be utilized

to as LMED-2.

accelerometer

the four loops utilized

to ACCEL-1Y,

for each of the feedback

Figure

location

the same

LMED

Loops

back each of the four colocated

to ACCEL-1X,

so that

device

on feeding

LMED-1X

was identical

to Accelerometer

of the

A block

to

were (i) and

(iv)

four loops

diagram

of the

4.3.1.

the

beam

(by the LMED

velocity along

(at the

the same

given axis) is

/_5¢ = H(s)s 2 + Ds + k._kd _p s 2 + D--s+

The

design

frequency

goal was to choose (say around

10 Hz) in order objectives

stroke

limitations

low frequency, filter.

could of the

low frequency

damping

LMED

stroke

proof

positive

above

Thus,

transfer

real to some

by simply mass

even this controller a first order

the

to the beam

be accomplished

is not implementable.

Unfortunately,

such that

1 Hz) and remains

to provide

design

H(s)

devices,

modes

choosing

H(s)

this controller,

the stroke

high pass filter was then 4-3

significantly

in this frequency

H(s) was initially

caused

function

chosen limitations

cascaded

= 1/s 2. which

is positive higher

real at low

frequency

band.

In theory

However, has very

with the second

the

due to the high

to be a second-order to be violated.

(say

gain

at

low-pass To limit the

order low pass

filter. The resultantcontrollerwasthus of the form ks H(s) The

low and

high pass portions

transformation control

gains

The Figure

with

It is seen that

more

the LMED

clearly

closed

in Figures

seen if Figures

of the

4.3.5-4.3.7 disturbance.

AGS gimbals,

the LMED frequency

in the

controllers

and was then

beam

performance

4.3.3

and 4.3.4 which

frequency

which When

vibration

that

harmonics

integrated

damping

to a BET-X

significantly ACCEL-1Y

with the feedback

the LOS performance

in the detector

responses.

4-4

loops.

The

to the higher

pulse

even

and ACCEL-2X.

of the responses

reduced.

Similar

and ACCEL-2Y

especially

frequency

is demonstrated

magnitude

controllers

4.3.2.-4.3.7.

pulse disturbance.

of ACCEL-1X

both the peak were

the bilinear

of the

in Figures

to a BET-X

show the responses

reveal

in each

is demonstrated

in providing

improvement

by using

A.

of BGYRO-Y

aided

show the BGYRO-Y,

loops improved

implemented

vibration

especially

responses

separately

are given in Appendix

loop

higher

pulse

of the higher

controller

were discretized

and closed loop responses

loop accelerometer

influences

BET-Y

prewarping

loop attenuation

The

The closed the

frequency

4.3.2 shows the open

harmonics.

of the controller

of the discretized

closed

= (s + _)(s _ + 2_,_,,s + w_,)"

involving by reducing

results

responses the IMC

and are to a and

the influence

..,._

l wt

r'

FILTER

z -

--_=_

H (z)

- CONTROLLERS

J(e)

Figure

4.1.1

The

loop was formulated

_=_

--

control

lim

SYSTEM

1+ •

LOW

as a disturbance

w, _- white

DESIGNED

E[qW(e)q(e)

problem

for both tracking

ACCOUNTS

1

,,@--=-=

WERE

white noise

problem.

4-5

AND

FREQUENCY

MODES

noise

BY

MINIMIZING

+ 9uW(e)u(e)]

the IMC-X

BIASES

FOR

to DET-Y

p > O.

loop and the IMC-Y

to DET-X

OPEN

LOOP

RESPONSE

5.0

STFIRT CONTR, }L

71

^10

TO BET-X

m

PULSE

1/

ON

4.0 m-

3.0

_,..'_t/'_,_-

y

I

,IIll Igl

Lr--

1.0

'dIlUI '

Z -.0 B B m

_EC

m

,

l

-5-_.

0 '

I

• :7;

,

,

.6

I

,3

z

Ii

12

At

I_

II

18

:'

TIME (SECS)

RESPONSE

S.0 2X10

WITH

IMC

TO DET

FEEDBACK

LOOPS

iS

-I m

4.0

m m

X E

...... lii"

"wl,._

z 1,0 rm

m

I XI ) -5-_.

Figure addition

0 '

I

I

• :2',

4.1.2. The IMC to bias correction

I

.6

I

I

,

,

_.4

,3

to DET feedback loops to improve the DET-X 4-6

were able to provide response to a BET-X

I

,

,

t7.7

low frequency tracking pulse disturbance.

,

:3,D

in

OPEN

5.0

LOOP

RESPONSE

TO BET-Y

PULSE

m

-×i0

1:21

-I .0

-e. o

E-'.-:-_-:. 0 LO

'_ , I_"

_ii*',,

-4.

1 I

I

I

TIME

RESPONSE

WITH

-t

IMC TO DET

START' COI4TR')L

-1 m

LOOPS

I I ,"16/09 Oil P

1:21 Irr

IlL)1,

-

, t! _.i_,,

tl

. 71 m m

-5.

FEEDBACK

.0

-2.0

-4

:3,0

(SECS)

5.0 E)< i 0

i

e

.¢j

.

_,.yl •/'l

I'_. i',_.. Yt

!

!1,

i 5EC

'

l

.3

,

I

_

.6

[

I

_3

t

!

I "-

t

Xl ) i

1

l. 8

2.

7'.4

I

L

2..7

:3.E_

TIME (SECS)

Figure 4.1.3. The IMC to DET feedback in addition to bias correction to substantially

loops were able to provide low frequency tracking improve the DET-Y response to a BET-Y pulse

disturbance.

4-7

COMPENSATOR

PHASE

IN THE

REGION

PERFORMANCE

SOLID-= W/OUT MAXIMUM DASHED=-WITH

2o8

ENTROPY

MAXIMI]M

ENTROPY

lo8

uu r_

i:i .........

"4

...... '::::--LLLII.I. ""

........... l..i" ,,

-lo8

-20O

I

i

i

i

I

1

2

3

4

5

HZ

Figure positive

4.2.1. Maximum Entropy design real in the performance region.

rendered

4-8

the

compensators

for the

AGS

to BGYRO

loops

COMPENSATOR

MAGNITUDE

i

SHAPE IN THE PERFORMANCE

'l

i

REGION

;

)

SOL/D-= W/OUT MAXIMUM DASHED=-WITH

ENTROPY

MAXIMUM

ENTROPY

102 • / .........

-4 //

J"

,\

/" / _'\

// \

('

,

/ / \'\

,,

"\

10i

/

I

I

I

1

1

2

3

4

5

6

HZ

4.2.2.

For the

AGS

to BGYRO

loops

Maximum

Entropy

design

magnitudes in the performance region, thus providing performance that the robust controllers were effectively reduced-order controllers.

4-9

smoothed robustness

out

the

and

compensator also indicating

COMPENSATOR

NOTCH

103 SOLID= W/OUT MAXIMUM DASHED=

WITH MAXIMUM

ENTROPY ENTROPY

lo 2 \

/

'\

/

\

/

/ /

\ 10 !

10 s

4 I

2

;

8

1o

I

|

x2 I

1',

1;

HZ

4.2.3. high

For frequency

the

AGS modes

to BGYRO by increasing

loops their

Maximum width

Entropy and

4-10

depth.

design

robustified

the

notches

for the

OPEN LOOP RESPONSE 5.0

TO

START CONTR

71 m

4.0 p,m

m

3.0

I/llllil illl liil

2.0

1.0

Ill

-.0

1

.0

TIME (SECS)

RESPONSE

WITH AGS TO BGYRO

5.0 --

m

LOOPS

1/'1(9/98 kJL ON P

START CONTR

-'1

ZXIO

FEEDBACK

In

4.0 L--

X

x-

_2.8

&t

"lfil/Iv

_

1,0

SF'C _.5

o C"t

_-

_

•2_i

,

,

.:7

I

t

.8

,5

t

I .2

_

I ._

i

J

,

1.8

i

t

2.1

l

I

F'.4

TIME (SECS)

Figure 4.2.4. The reduction beam vibrationof the AGS substantially improve the DET-X

response to a BET-X

4-11

to BGYRO

feedback loops was able to

pulse disturbance.

OPEN

LOOP

RESPONSE

TO BET-X

PULSE

0"02 / 0.0151 0.01

_, 0.oo5

+° e_.0.005 -0.01 -0.015

-0.020

1_0

I

t

15

20

25

30

TiME (SECS)

RESPONSE

WITH

AGS

TO BGYRO

FEEDBACK

LOOPS

0"02 I 0.015' 0.01

;>.

II

!

0.005

6 -0.005 -0.01 -0.015 _

-0.020

10

L

15

2'0

9_5

30

TIME (SECS) Figure 4.2.5. For a BET-X pulse a comparison response with the AGS to BGYRO feedback loops of the beam vibration.

4-12

of the open loop BGYRO-Y response to the closed reveals significant closed loop damping

OPEN

LOOP

RESPONSE

TO

BET-Y

0

i- I

-I

-2

1

0 _=. ,.=.

i !,_

ll!!jit

m

El

-4

.0

V

tl'_u1

p

-5

, i

• :3

!3, • 0

.6

TIME (SECS)

RESPONSE

WITH

AGS

TO

BGYRO

5.0 -XIO -I

FEEDBACK

START i CON'rR,

m

LOOPS 9

1/

5:

ON

)L

.0

=- li, l_!

-2.0

-'

l,lll,,l/.ll,

0

%;%(, "%._,,-.

l',jlg',._j' _,l' dr'J"b[ _r"

--4 . E1 I 5Ei-j I

-5,_





.Ei '

.J

'

.9

11 ,

I

' 1._

1.2

1.8

I

2.

;:'.4

to BGYRO

feedback

c

:3,0

,

TIME (SEem)

Figure

4.2.6.

to substantially

The

reduction

improve

the

in beam DET-Y

vibration response

of the to a BET-Y

4-13

AGS

pulse

disturbance.

loops

was

able

OPEN 0.02 0 • 015'

LOOP

RESPONSE

TO BET-Y

PULSE

i

0.01 0.005

6 0

o -0.005 -0.01

i

lo

_5

i

I

TIME

RESPONSE

WITH

io

2_

3O

(SECS)

AGS TO BGYRO

FEEDBACK

LOOPS

0.02 0.015 0.01 0.005

6

o -0.005 -0.01 -0.015

"0"020

5

i'0

15

2'0

25

30

TIME(SECS) Figure 4.2.7. For a BET-Y pulse a comparison response with the AGS to BGYRO feedback loops of the beam vibration.

4-14

of the open loop BGYRO-X response to the closed reveals significant closed loop damping

fc

kmkd x s/p

p= beam

acceleration

x s/p= relative

velocity

of proof mass

m = mass of proof mass f c = force

applied

k s = position D = inherent k

loop stiffness viscous

rp = position

with

4.3.1. control

The

damping

of the LMED

mkd = motor force constants

H(s) = compensator

Figure

to structure

LMED

designs

transfer

function

command.

assumed

that

this block

law H(s).

4-15

diagram

described

each feedback

loop

OPEN

LOOP

0.02

RESPONSE

,

TO BET-X

1

PULSE ,

0.015 0.01 0.005

6

o

c3 _-O.OOfi

4).01 -0.015

-0.02{

3O

TIME (SECS)

RESPONSE

WITH

COLOCATED

LMED

FEEDBACK

LOOPS

0.02 0.015 0`01 0.005

_-0`005

-0.01 -0.015

-o.02_

}

1'o

io

i_

23

_o

TIME (SECS) Figure

4.3.2.

response with harmonics.

For a BET-X the

LMED

loops

pulse closed

a comparison reveals

some

4-16

of the closed

open

loop

BGYRO-Y

loop damping

response

of the higher

to the

frequency

OPEN

0.6

LOOP

RESPONSE

TO BET-X

PULSE

0.4

0.2

m

0

< -0.2

-0.4

-0.6 0

3O TIME

RESPONSE

0.6

WITH

(SEGS)

COLOCATED

LMED

FEEDBACK

LOOPS

0.4

O.2 X

o

cJ