def side_force_coeff_test_cases():
constants = SkywalkerX8Constants()
b = constants.wing_span
state1 = State()
state1.vx = 20.0
state1.vy = 0.0
state1.vz = 0.0
state2 = State()
state2.vx = 28.6362532829
state2.vy = 1.0
state2.vz = 0.0
state2.ang_rate_x = 5 * np.pi / 180
state2.ang_rate_z = 5 * np.pi / 180
wind = np.zeros(6)
airspeed = np.sqrt(np.sum(calc_airspeed(state2, wind)**2))
zero_input = ControlInput()
aileron_input = ControlInput()
aileron_input.aileron_deflection = 2.0 * np.pi / 180.0
return [
(state1, wind, 0.0, zero_input, 2.99968641720902E-08),
(state1, wind, 1.0, zero_input, -7.94977329537982E-06),
(state2, wind, 0.0, aileron_input,
-0.008604744865183 + -0.085 * b / (2 * airspeed) * state2.ang_rate_x +
0.005 * b / (2 * airspeed) * state2.ang_rate_z +
0.0433 * aileron_input.aileron_deflection),
(state2, wind, 1.0, aileron_input,
-0.007089388672593 + -0.133 * b / (2 * airspeed) * state2.ang_rate_x +
0.002 * b / (2 * airspeed) * state2.ang_rate_z +
0.0433 * aileron_input.aileron_deflection)
]
def yaw_moment_coeff_test_cases():
constants = SkywalkerX8Constants()
b = constants.wing_span
state1 = State()
state1.vx = 20.0
state1.vy = 0.0
state1.vz = 0.0
state2 = State()
state2.vx = 28.6362532829
state2.vy = 1.0
state2.vz = 0.0
state2.ang_rate_x = 5 * np.pi / 180
state2.ang_rate_z = 5 * np.pi / 180
wind = np.zeros(6)
airspeed = np.sqrt(np.sum(calc_airspeed(state2, wind)**2))
zero_input = ControlInput()
aileron_input = ControlInput()
aileron_input.aileron_deflection = 2.0 * np.pi / 180.0
return [
(state1, wind, 0.0, zero_input, 4.9176697574439E-06),
(state1, wind, 1.0, zero_input, 1.96093394589053E-05),
(state2, wind, 0.0, aileron_input,
0.000825947539055 + 0.027 * b / (2 * airspeed) * state2.ang_rate_x +
-0.022 * b / (2 * airspeed) * state2.ang_rate_z -
0.00339 * aileron_input.aileron_deflection),
(state2, wind, 1.0, aileron_input,
0.001052911121301 + 0.017 * b / (2 * airspeed) * state2.ang_rate_x +
-0.049 * b / (2 * airspeed) * state2.ang_rate_z -
0.00339 * aileron_input.aileron_deflection)
]
def roll_moment_coeff_test_cases():
constants = SkywalkerX8Constants()
b = constants.wing_span
state1 = State()
state1.vx = 20.0
state1.vy = 0.0
state1.vz = 0.0
state2 = State()
state2.vx = 28.6362532829
state2.vy = 1.0
state2.vz = 0.0
state2.ang_rate_x = 5 * np.pi / 180
state2.ang_rate_z = 5 * np.pi / 180
wind = np.zeros(6)
airspeed = np.sqrt(np.sum(calc_airspeed(state2, wind)**2))
zero_input = ControlInput()
aileron_input = ControlInput()
aileron_input.aileron_deflection = 2.0 * np.pi / 180.0
return [
(state1, wind, 0.0, zero_input, -8.40821757613653E-05),
(state1, wind, 1.0, zero_input, -7.34515369827804E-05),
(state2, wind, 0.0, aileron_input,
-0.00380800071177 + -0.409 * b / (2 * airspeed) * state2.ang_rate_x +
0.039 * b / (2 * airspeed) * state2.ang_rate_z +
0.12 * aileron_input.aileron_deflection),
(state2, wind, 1.0, aileron_input,
-0.003067251004494 + -0.407 * b / (2 * airspeed) * state2.ang_rate_x +
0.158 * b / (2 * airspeed) * state2.ang_rate_z +
0.12 * aileron_input.aileron_deflection)
]
def test_aileron_deflection():
control_input = ControlInput()
aileron_deflection = 0.2
control_input.aileron_deflection = aileron_deflection
assert control_input.aileron_deflection == pytest.approx(
aileron_deflection)
assert np.allclose(control_input.control_input,
[0.0, aileron_deflection, 0.0, 0.0])
def test_dynamics_forces():
control_input = ControlInput()
prop = SimpleTestAircraftNoMoments(control_input)
t = 0
for i in range(-50, 101, 50):
control_input.throttle = 0.8
control_input.elevator_deflection = i
control_input.aileron_deflection = i
control_input.rudder_deflection = i
state = State()
state.vx = 20.0
state.vy = 1
state.vz = 0
params = {"prop": prop, "wind": no_wind()}
update = dynamics_kinetmatics_update(t=t,
x=state.state,
u=control_input.control_input,
params=params)
V_a = np.sqrt(np.sum(calc_airspeed(state, params['wind'].get(0.0))**2))
forces_aero_wind_frame = np.array([
-np.abs(control_input.elevator_deflection),
control_input.aileron_deflection, -control_input.rudder_deflection
])
forces_aero_body_frame = wind2body(forces_aero_wind_frame, state,
params['wind'].get(0))
force_propulsion = np.array([(2 * control_input.throttle)**2 - V_a**2,
0, 0])
force_gravity = inertial2body(
np.array([0, 0, prop.mass() * GRAVITY_CONST]), state)
forces_body = forces_aero_body_frame + force_propulsion + force_gravity
vx_update_expect = (1 / prop.mass()) * forces_body[0]
vy_update_expect = (1 / prop.mass()) * forces_body[1]
vz_update_expect = (1 / prop.mass()) * forces_body[2]
# No moments
ang_rate_x_update_expect = 0
ang_rate_y_update_expect = 0
ang_rate_z_update_expect = 0
assert np.allclose(vx_update_expect, update[6])
assert np.allclose(vy_update_expect, update[7])
assert np.allclose(vz_update_expect, update[8])
assert np.allclose(ang_rate_x_update_expect, update[9])
assert np.allclose(ang_rate_y_update_expect, update[10])
assert np.allclose(ang_rate_z_update_expect, update[11])
def test_dynamics_moments():
control_input = ControlInput()
t = 0
for i in range(-50, 51, 50):
control_input.throttle = i
control_input.elevator_deflection = i
control_input.aileron_deflection = i
control_input.rudder_deflection = i
prop = SimpleTestAircraftNoForces(control_input)
state = State()
state.vx = 20.0
state.vy = 1
state.vz = 0
state.ang_rate_x = 0.157079633
state.ang_rate_y = 0.157079633
state.ang_rate_z = 0.157079633
params = {"prop": prop, "wind": no_wind()}
update = dynamics_kinetmatics_update(t=t,
x=state.state,
u=control_input.control_input,
params=params)
moments_aero = np.array([
control_input.elevator_deflection,
control_input.aileron_deflection, control_input.rudder_deflection
])
omega = np.array(
[state.ang_rate_x, state.ang_rate_y, state.ang_rate_z])
coreolis_term = prop.inv_inertia_matrix().dot(
np.cross(omega,
prop.inertia_matrix().dot(omega)))
ang_rate_x_update_expect = (2 / 3) * moments_aero[0] - (
1 / 3) * moments_aero[2] - coreolis_term[0]
ang_rate_y_update_expect = (1 / 2) * moments_aero[1] - coreolis_term[1]
ang_rate_z_update_expect = (2 / 3) * moments_aero[2] - (
1 / 3) * moments_aero[0] - coreolis_term[2]
assert np.allclose(ang_rate_x_update_expect, update[9])
assert np.allclose(ang_rate_y_update_expect, update[10])
assert np.allclose(ang_rate_z_update_expect, update[11])