Scribd Logo
Upload

Welcome to Scribd!

  • 18 views

    Jet Engine

    Uploaded by

    Ebubechim Chiburikem Ogunka-Nnoka
    The document analyzes the restoring force and numerical solution of an oscillating system. Equations are derived for the restoring force and natural frequency. Plots compare the numerical solution to the exact solution over time.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd

    Jet Engine

    Uploaded by

    Ebubechim Chiburikem Ogunka-Nnoka
    18 views2 pages
    The document analyzes the restoring force and numerical solution of an oscillating system. Equations are derived for the restoring force and natural frequency. Plots compare the numerical solution to the exact solution over time.

    Copyright

    © © All Rights Reserved

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    The document analyzes the restoring force and numerical solution of an oscillating system. Equations are derived for the restoring force and natural frequency. Plots compare the numerical solution to the exact solution over time.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    18 views2 pages

    Jet Engine

    Uploaded by

    Ebubechim Chiburikem Ogunka-Nnoka
    The document analyzes the restoring force and numerical solution of an oscillating system. Equations are derived for the restoring force and natural frequency. Plots compare the numerical solution to the exact solution over time.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    of 2

    I n [ ] : = j = 240; m = 40; l = 2.7; k = 4; g = 9.81; A = 90.

    0 * Pi / 180;
    fun1[x_] := δ * x + β * Sin[x]
    k m*g*l
    δ= ;β= ;
    j + m * l2  j + m * l2 
    β (4 - 4 * Cos[A] - A * Sin[A])
    α2 = δ + ;
    A2
    2 * β (A * Sin[A] + (2 * Cos[A]) - 2)
    α4 = ;
    A4
    α4 * A 2
    φ= α2 + α4 * A2  ; m = ; Tp = 4 * EllipticK[m] / φ; ωn = 2 * Pi / Tp;
    2 * α2 + α4 * A2 
    fun4[x_] := α2 * x + α4 * x3 ; (*EQUIVALENT SOLUTION*)
    1
    freq = ;
    Tp

    (*PLOT RESTORING FORCE PLOT*)


    Show[{Plot[fun1[x], {x, - A, A}, PlotStyle  { Blue }],
    Plot[fun4[x], {x, - A, A}, PlotStyle  { Red, Dashing[Large]}]},
    AxesLabel  {"u", "f(u),g(u)"}, PlotRange  All]

    (*NUMERICAL SOLUTION *)
    sol = NDSolve[{x ''[t] + fun1[x[t]]  0 , x[0]  A, x '[0]  0} , x, {t, 0, Tp} ]
    solFun = x /. sol[ [1]]

    (*DISPLACEMENT PLOT OF NUMERICAL VS CUBICATION*)


    Show[
    Plot[{A * JacobiCN[φ * t, m]}, {t, 0, Tp}],
    Plot[Evaluate[x[t] /. sol], {t, 0, Tp}, PlotStyle  { Black , Dashed } ],

    (*VELOCITY PLOT OF NUMERICAL VS CUBICATION*)


    Plot[{- A * φ * JacobiSN[φ * t, m] * JacobiDN[φ * t, m]}, {t, 0, Tp}],
    Plot[Evaluate[x '[t] /. sol], {t, 0, Tp}, PlotStyle  { Blue, Dashed }],
    AxesLabel  {"t", "u,du/dt"}, PlotRange  All]
    2 JET ENGINE.nb

    O u t [ ] =
    f(u),g(u)

    u
    -1.5 -1.0 -0.5 0.5 1.0 1.5

    -1

    -2

    O u t [ ] =

    Domain: 0., 5.23


    x  InterpolatingFunction 
    Output: scalar

    O u t [ ] =

    Domain: 0., 5.23


    InterpolatingFunction 
    Output: scalar

    O u t [ ] =
    u,du/dt

    t
    1 2 3 4 5

    -1

    -2

    You might also like