コード例 #1
0
def main(output='flutter-velocity-example.pdf'):
    """
    >>> import io
    >>> fobj = io.BytesIO()
    >>> main(fobj)
    >>> assert len(fobj.getvalue()) > 0
    >>> assert fobj.getvalue()[:4] == b'%PDF'

    """
    zs = np.linspace(0, 20000, 200)
    ps, ts, ss = model_atmosphere(zs)
    rhos = (ps/1000) / (0.2869 * (ts + 273.1))
    vs_t = flutter_velocity_transonic(ps, ss, 14.02, 7, 8.59, 0.5)
    vs_s1 = flutter_velocity_supersonic(rhos, 5026, 11671, 0.0923, 0.0859, 0.0284, 0.0040, 4.3)
    vs_s2 = flutter_velocity_supersonic(rhos, 5026, 11671, 0.0923, 0.0859, 0.0284, 0.0040, 3)
    vs_s3 = flutter_velocity_supersonic(rhos, 5026, 11671, 0.0923, 0.0859, 0.0284, 0.0040, 2)

    plt.figure()
    #plt.plot(zs * 1e-3, vs_t, 'r', label="transonic flutter velocity")
    plt.plot(zs*1e-3, vs_t/343.2, 'c', label="Theodorson")
    plt.plot(zs*1e-3, vs_s1/343.2, 'g', label="Mach 4.3")
    plt.plot(zs*1e-3, vs_s2/343.2, 'r', label="Mach 3")
    plt.plot(zs*1e-3, vs_s3/343.2, 'b', label="Mach 2")

    plt.title('Flutter velocity vs altitude, Core/Booster fin')
    plt.xlabel('Altitude [km]')
    plt.ylabel('Flutter Velocity [Mach]')
    plt.legend()
    plt.savefig(output, format='PDF')
コード例 #2
0
ファイル: plot_fin_flutter.py プロジェクト: lwc24/m3cfd
def main(output="flutter-velocity-example.pdf"):
    """
    >>> import io
    >>> fobj = io.BytesIO()
    >>> main(fobj)
    >>> assert len(fobj.getvalue()) > 0
    >>> assert fobj.getvalue()[:4] == b'%PDF'

    """
    zs = np.linspace(0, 20000, 200)
    ps, ts, ss = model_atmosphere(zs)
    rhos = (ps / 1000) / (0.2869 * (ts + 273.1))
    vs_t = flutter_velocity_transonic(ps, ss, 20, 10, 10, 0.2)
    vs_s = flutter_velocity_supersonic(rhos, 380, 104, 1, 0.3, 0.0313, 0.0855, 3)

    plt.figure()
    plt.plot(zs * 1e-3, vs_t, "r", label="transonic flutter velocity")
    plt.plot(zs * 1e-3, vs_s, "g", label="supersonic flutter velocity, Mach 3")
    plt.title("Flutter velocity vs altitude")
    plt.xlabel("Altitude [km]")
    plt.ylabel("Flutter velocity [m/s]")
    plt.legend()
    plt.savefig(output, format="PDF")
コード例 #3
0
def main(output='flutter-velocity-example.pdf'):
    """
    >>> import io
    >>> fobj = io.BytesIO()
    >>> main(fobj)
    >>> assert len(fobj.getvalue()) > 0
    >>> assert fobj.getvalue()[:4] == b'%PDF'

    """
    zs = np.linspace(0, 20000, 200)
    ps, ts, ss = model_atmosphere(zs)
    rhos = (ps / 1000) / (0.2869 * (ts + 273.1))
    vs_t = flutter_velocity_transonic(ps, ss, 20, 10, 10, 0.2)
    vs_s = flutter_velocity_supersonic(rhos, 380, 104, 1, 0.3, 0.0313, 0.0855,
                                       3)

    plt.figure()
    plt.plot(zs * 1e-3, vs_t, 'r', label="transonic flutter velocity")
    plt.plot(zs * 1e-3, vs_s, 'g', label="supersonic flutter velocity, Mach 3")
    plt.title('Flutter velocity vs altitude')
    plt.xlabel('Altitude [km]')
    plt.ylabel('Flutter velocity [m/s]')
    plt.legend()
    plt.savefig(output, format='PDF')