from simupy.block_diagram import BlockDiagram import simupy_flight import numpy as np from nesc_testcase_helper import plot_nesc_comparisons, int_opts, benchmark from nesc_testcase_helper import ft_per_m, kg_per_slug planet = simupy_flight.Planet( gravity=simupy_flight.earth_J2_gravity, winds=simupy_flight.get_constant_winds(), atmosphere=simupy_flight.atmosphere_1976, planetodetics=simupy_flight.Planetodetic( a=simupy_flight.earth_equitorial_radius, omega_p=simupy_flight.earth_rotation_rate, f=simupy_flight.earth_f ) ) Ixx = 0.001894220*kg_per_slug/(ft_per_m**2) #slug-ft2 Iyy = 0.006211019*kg_per_slug/(ft_per_m**2) #slug-ft2 Izz = 0.007194665*kg_per_slug/(ft_per_m**2) #slug-ft2 Ixy = 0.0*kg_per_slug/(ft_per_m**2) #slug-ft2 Iyz = 0.0*kg_per_slug/(ft_per_m**2) #slug-ft2 Izx = 0.0*kg_per_slug/(ft_per_m**2) #slug-ft2 m = 0.155404754*kg_per_slug #slug x = 0. y = 0. z = 0. S_A = 0.22222/(ft_per_m**2)
from simupy.block_diagram import BlockDiagram import simupy_flight import numpy as np from nesc_testcase_helper import plot_nesc_comparisons, int_opts, benchmark from nesc_testcase_helper import ft_per_m, kg_per_slug planet = simupy_flight.Planet(gravity=simupy_flight.earth_J2_gravity, winds=simupy_flight.get_constant_winds(wy=20 / ft_per_m), atmosphere=simupy_flight.atmosphere_1976, planetodetics=simupy_flight.Planetodetic( a=simupy_flight.earth_equitorial_radius, omega_p=simupy_flight.earth_rotation_rate, f=simupy_flight.earth_f)) Ixx = 3.6 * kg_per_slug / (ft_per_m**2) #slug-ft2 Iyy = 3.6 * kg_per_slug / (ft_per_m**2) #slug-ft2 Izz = 3.6 * kg_per_slug / (ft_per_m**2) #slug-ft2 Ixy = 0.0 * kg_per_slug / (ft_per_m**2) #slug-ft2 Iyz = 0.0 * kg_per_slug / (ft_per_m**2) #slug-ft2 Izx = 0.0 * kg_per_slug / (ft_per_m**2) #slug-ft2 m = 1.0 * kg_per_slug #slug x = 0. y = 0. z = 0. S_A = 0.1963495 / (ft_per_m**2) b_l = 1.0 c_l = 1.0