Mechanical Engineering
MEM 355 Final Project, Fall 2020

Due: Friday, December 4, 2020

Preliminaries A rigid aircraft has 6 degrees of freedom comprising three translational coordinates and three rotational coordinates. Thus, it has 12 states: 3 space coordinates X,Y,Z and 3 angle coordinates , ,? ? ? roll, pitch and yaw as well as 6 corresponding velocities translational in body coordinates, , ,u v w , and , ,p q r . It is sometimes convenient to replace the translation velocity vector , ,u v w by its polar representation

, ,V ? ? where V is the speed (magnitude of the velocity vector), ? is the angle of attack and ? is the sideslip angle.

If a flat earth is assumed and variation of atmospheric density is ignored then the equations are invariant with respect to the X,Y,Z coordinates and heading ? . So we can ignore them, reducing the number of equations to 8. Furthermore, if the equations are linearized at some trim point they decouple into two smaller sets: a) the longitudinal dynamics consisting of the 4 states , , ,u w q ? (equivalently

, , ,V q? ? ) and b) the lateral dynamics consisting of the remaining 4 states , , ,v p r ? (equivalently, , , ,p r? ? ). Note that the longitudinal dynamics define the motion in the vertical plane.

Part A Consider the F-18 Longitudinal Dynamics:

0.0239 28.3172 0 32.2 3.8114 0.001 0.0003 0.3621 1 0 0.0515 0

0 2.2115 0.2532 0 2.8791 0 0 0 1 0 0 0

e

T

V V d

q qdt ?? ? ?

? ?

? ? ? ?? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ? ? ?? ? ? ? ? ? ? ?= + ? ?? ? ? ? ? ? ? ?? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?

0 1 0 1 1 0 0 0

V

V q ? ?

?

? ? ? ??? ? ? ? ? ?=? ? ? ? ? ?? ? ? ? ? ? ? ?

1. Compute describe the open loop system modes. 2. Compute the 2×2 system transfer matrix , ,e T V? ? ?? :

( ) ( ) ( ) ( ) ( )

11 12

21 22

G s G s G s

G s G s ? ?

= ? ? ? ?

3. Consider the single input – single output system e? ?? . Describe the relevant features of its transfer function, ( )11G s . Obtain a state feedback control gain, K, that places the closed loop poles at ( )0.05 0.09, 1.05 1.1j j? ± ? ± . Obtain an observer feedback gain L that places observer poles in accordance with ‘Rynaski’ design, e.g., place observer poles at ( )1.149, 8, 8.5, 9? ? ? ? .

4. Determine the equivalent compensator, ( )cG s , obtained in Part 3 and compute the closed loop system transfer function. Compute the closed loop transfer function for compensator ( )cG s and plant

( )11G s .

Part B Consider the lateral dynamics of a Boeing 747 at 0 774 fpsV = and 40,000 feet. The control input is the rudder deflection r? :

0.0558 0.9968 0.0802 0.0415 0.00729 0.598 0.115 0.0318 0 0.475

, 3.05 0.388 0.4650 0 0.153 0 0.0805 1 0 0

r

r rd y r

p pdt

? ?

?

? ?

? ?? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ?? ? ? ? ? ? ? ?= + = ? ? ? ? ? ? ? ?? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?

Nomenclature: ? sideslip angle, r yaw rate, p roll rate, ? roll angle. We wish to design a compensator to improve yaw rate control performance. The standard

configuration of such a control system is shown the figure below and includes a ‘washout filter1.’ The rudder actuator dynamics need to be considered as well as shown.

Aircraft

Washout Filter

Rudder ActuatorCompensator

– +

( )cG s 10

10s +

1 3

1 3

s s +

r? rr

Construct a set of state equations that includes the rudder actuator and the washout filter.

1. Choose the following set of closed loop poles: { }0.0051, 0.468, 0.279 0.628 , 0.279 0.628 , 1.106, 9.89j j? ? ? + ? ? ? ?

and find a state feedback controller to achieve those poles. 2. Compare the closed loop and open loop modes. Discuss how they have changed. 3. Design an observer for the system with poles 5 times faster than controller poles:

{ }0.0253, 2.34, 1.39 3.14 , 1.39 3.14 , 5.53, 49.5j j? ? ? + ? ? ? ? 4. Compute the compensator transfer function.

1 Turning an aircraft requires the pilot to coordinate rudder and aileron to establish a desired yaw rate and roll angle. The presence of the washout filter reduces the impact of the rudder on roll dynamics, thereby making the pilot’s job easier.