F-8 aircraft (JModelica)

From mintOC
Revision as of 13:47, 15 March 2016 by MadeleineSchroeter (Talk | contribs)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

This page contains the model formulation of the MIOCP F-8 aircraft in JModelica format.

JModelica

The model in JModelica code.

// The F-8 aircraft control problem is based on a very simple aircraft model.
// (c) Sebastian Sager, 2005-2009
// More info on http://mintoc.de/index.php/F-8_aircraft
 
 optimization F8_aircraft_opt (objective = tf,
                         startTime = 0,
                         finalTime = 1)
 
    // Parameters
    parameter Real ksi=0.05236;
    parameter Real tf(free=true,min=0.1);
 
    // The states
    Real x0(start=0.4655, fixed=true); //angle of attack
    Real x1(start=0, fixed=true); //pitch angle
    Real x2(start=0, fixed=true); //pitch rate
 
    // The control signal
    input Real  w(min=0,max=1);
 
  equation
 
    der(x0)= 1*tf*(-0.877*x0 + x2 - 0.088*x0*x2 + 0.47*x0^2 - 0.019*x1^2 -x0^2*x2 + 3.846*x0^3-(0.215*ksi-0.28*x0^2*ksi -0.47*x0*ksi^2 - 0.63*ksi^3)*w-(-0.215*ksi + 0.28*x0^2 -0.47*x0*ksi^2 + 0.63*ksi^3)*(1-w));
    der(x1) = 1*tf*x2;
    der(x2) = 1*tf*(-4.208*x0 - 0.396*x2 - 0.47*x0^2 - 3.564*x0^3+20.967*ksi - 6.265*x0^2*ksi + 46*x0*ksi^2 -61.4*ksi^3-(20.967*ksi - 6.265*x0^2*ksi - 61.4*ksi^3)*2*w);
 
  constraint
    x0(finalTime)=0;
    x1(finalTime)=0;
    x2(finalTime)=0;
 
  end F8_aircraft_opt;