13.2.3: Aircraft motion exercise
- Page ID
- 78427
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)In order to obtain the 4D trajectory, the flight will be divided into three phases: climb, cruise, and descent.
1. The climb phase will be assumed to be a symmetric flight into the vertical plane non considering any wind and assuming the heading angle to be zero. Therefore, the ODE system to be used is as follows:
\[m \dot{V} = T - D - m \cdot g \cdot \sin \gamma, \nonumber \]
\[m V \dot{\gamma} = L - m \cdot g \cdot \cos \gamma,\nonumber \]
\[\dot{x}_e = V \cdot \cos \gamma,\nonumber \]
\[\dot{h}_e = V \cdot \sin \gamma,\nonumber \]
\[\dot{m} = -T \cdot \eta.\nonumber \]
where according to BADA, \(\eta = ( \tfrac{C_{f1}}{1000*60} ) \cdot (1 + \tfrac{V}{C_{f2}} )\), with \(V\) in knots.
In order to integrate the system, one should:
1.1 Set the initial conditions for all the state variables, initial and final time. This conditions must be selected according to typical values of aircraft performance.
1.2 Set the control variables \((T(t), C_L (t))\),7 for instance to the following values:
* \(T = 0.8 \cdot T_{\max}\), where \(T_{\max} = C_{tc1} \cdot (1 - \tfrac{h_e}{C_{tc2}} + C_{tc3} \cdot h_e^2)\), \(h_e\) in feet,
* \(C_L = C_{L_{opt}}\).
1.3 Use a suitable numerical method to solve the resulting system.
2. The cruise phase will be assumed to be a symmetric flight into the horizontal plane not considering any wind. Therefore, the ODE system to be used is as follows:
\[m \dot{V} = T - D,\nonumber \]
\[m V \dot{\chi} = L \sin \mu, \nonumber \]
\[\dot{x}_e = V \cos \chi,\nonumber \]
\[\dot{y}_e = V \sin \chi, \nonumber \]
\[\dot{m} = -T \eta,\nonumber \]
being \(L = \tfrac{mg}{\cos \mu}\).
In order to solve the system, one should:
2.1 Set the initial conditions for all the state variables, initial and final time. Initial conditions and initial time will coincide with the final conditions of the previous phase. Set the final time according to typical values of aircraft performance.
2.2 Set the control variables \((T(t), \mu (t))\) to the following values:
* \(T = 0.5 * T_{\max}\),
* \(\mu = 0\).
2.3 Use a suitable numerical method to solve the resulting system.
3. The landing phase will be assumed to be gliding performance not considering any wind. Therefore, the ODE system to be used is as follows:
\[m \dot{V} = -D - mg \sin \gamma, \nonumber \]
\[m V \dot{\gamma} = L - mg \cos \gamma, \nonumber \]
\[\dot{x}_e = V \cos \gamma,\nonumber \]
\[\dot{h}_e = V \sin \gamma, \nonumber \]
3.1 Set the initial conditions for all the state variables, initial and final time. Initial conditions and initial time will coincide with the final conditions of the previous phase. Set the final time according to typical values of aircraft performance.
3.2 Set the control variable \((C_L(t))\) to the following value:
* \(C_L = C_{L_{opt}}.\)
3.3 Use a suitable numerical method to solve the resulting system.
7. Notice that \(C_L\) acts as control.