src/FDM/JSBSim/models/propulsion/FGTurboProp.cpp

   1 /*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   2
   3  Module:       FGTurboProp.cpp
   4  Author:       Jiri "Javky" Javurek
   5                based on SimTurbine and Turbine engine from David Culp
   6  Date started: 05/14/2004
   7  Purpose:      This module models a turbo propeller engine.
   8
   9  ------------- Copyright (C) 2004  (javky@email.cz) ---------
  10
  11  This program is free software; you can redistribute it and/or modify it under
  12  the terms of the GNU Lesser General Public License as published by the Free Software
  13  Foundation; either version 2 of the License, or (at your option) any later
  14  version.
  15
  16  This program is distributed in the hope that it will be useful, but WITHOUT
  17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
  18  FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public License for more
  19  details.
  20
  21  You should have received a copy of the GNU Lesser General Public License along with
  22  this program; if not, write to the Free Software Foundation, Inc., 59 Temple
  23  Place - Suite 330, Boston, MA  02111-1307, USA.
  24
  25  Further information about the GNU Lesser General Public License can also be found on
  26  the world wide web at http://www.gnu.org.
  27
  28 FUNCTIONAL DESCRIPTION
  29 --------------------------------------------------------------------------------
  30
  31 This class descends from the FGEngine class and models a Turbo propeller engine
  32 based on parameters given in the engine config file for this class
  33
  34 HISTORY
  35 --------------------------------------------------------------------------------
  36 05/14/2004  Created
  37 02/08/2011  T. Kreitler, added rotor support
  38
  39 //JVK (mark)
  40
  41 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  42 INCLUDES
  43 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/
  44
  45 #include <iostream>
  46 #include <sstream>
  47 #include "FGTurboProp.h"
  48 #include "FGPropeller.h"
  49 #include "FGRotor.h"
  50 #include "models/FGPropulsion.h"
  51 #include "models/FGAuxiliary.h"
  52
  53 using namespace std;
  54
  55 namespace JSBSim {
  56
  57 static const char *IdSrc = "$Id: FGTurboProp.cpp,v 1.19 2011/03/10 01:35:25 dpculp Exp $";
  58 static const char *IdHdr = ID_TURBOPROP;
  59
  60 /*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  61 CLASS IMPLEMENTATION
  62 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/
  63
  64 FGTurboProp::FGTurboProp(FGFDMExec* exec, Element *el, int engine_number)
  65   : FGEngine(exec, el, engine_number),
  66     ITT_N1(NULL), EnginePowerRPM_N1(NULL), EnginePowerVC(NULL)
  67 {
  68   SetDefaults();
  69   thrusterType = Thruster->GetType();
  70
  71   Load(exec, el);
  72   bindmodel();
  73   Debug(0);
  74 }
  75
  76 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  77
  78 FGTurboProp::~FGTurboProp()
  79 {
  80   delete ITT_N1;
  81   delete EnginePowerRPM_N1;
  82   delete EnginePowerVC;
  83   Debug(1);
  84 }
  85
  86 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  87
  88 bool FGTurboProp::Load(FGFDMExec* exec, Element *el)
  89 {
  90   IdleFF=-1;
  91   MaxStartingTime = 999999; //very big timeout -> infinite
  92   Ielu_max_torque=-1;
  93
  94 // ToDo: Need to make sure units are properly accounted for below.
  95
  96   if (el->FindElement("milthrust"))
  97     MilThrust = el->FindElementValueAsNumberConvertTo("milthrust","LBS");
  98   if (el->FindElement("idlen1"))
  99     IdleN1 = el->FindElementValueAsNumber("idlen1");
 100   if (el->FindElement("idlen2"))
 101     IdleN2 = el->FindElementValueAsNumber("idlen2");
 102   if (el->FindElement("maxn1"))
 103     MaxN1 = el->FindElementValueAsNumber("maxn1");
 104   if (el->FindElement("maxn2"))
 105     MaxN2 = el->FindElementValueAsNumber("maxn2");
 106   if (el->FindElement("betarangeend"))
 107     BetaRangeThrottleEnd = el->FindElementValueAsNumber("betarangeend")/100.0;
 108   BetaRangeThrottleEnd = Constrain(0.0, BetaRangeThrottleEnd, 0.99999);
 109   if (el->FindElement("reversemaxpower"))
 110     ReverseMaxPower = el->FindElementValueAsNumber("reversemaxpower")/100.0;
 111
 112   if (el->FindElement("maxpower"))
 113     MaxPower = el->FindElementValueAsNumber("maxpower");
 114   if (el->FindElement("idlefuelflow"))
 115     IdleFF = el->FindElementValueAsNumber("idlefuelflow");
 116   if (el->FindElement("psfc"))
 117     PSFC = el->FindElementValueAsNumber("psfc");
 118   if (el->FindElement("n1idle_max_delay"))
 119     Idle_Max_Delay = el->FindElementValueAsNumber("n1idle_max_delay");
 120   if (el->FindElement("maxstartingtime"))
 121     MaxStartingTime = el->FindElementValueAsNumber("maxstartingtime");
 122   if (el->FindElement("startern1"))
 123     StarterN1 = el->FindElementValueAsNumber("startern1");
 124   if (el->FindElement("ielumaxtorque"))
 125     Ielu_max_torque = el->FindElementValueAsNumber("ielumaxtorque");
 126   if (el->FindElement("itt_delay"))
 127     ITT_Delay = el->FindElementValueAsNumber("itt_delay");
 128
 129   Element *table_element;
 130   string name;
 131   FGPropertyManager* PropertyManager = exec->GetPropertyManager();
 132
 133   while (true) {
 134     table_element = el->FindNextElement("table");
 135     if (!table_element) break;
 136     name = table_element->GetAttributeValue("name");
 137     if (name == "EnginePowerVC") {
 138       EnginePowerVC = new FGTable(PropertyManager, table_element);
 139     } else if (name == "EnginePowerRPM_N1") {
 140       EnginePowerRPM_N1 = new FGTable(PropertyManager, table_element);
 141     } else if (name == "ITT_N1") {
 142       ITT_N1 = new FGTable(PropertyManager, table_element);
 143     } else {
 144       cerr << "Unknown table type: " << name << " in turbine definition." <<
 145       endl;
 146     }
 147   }
 148
 149   // Pre-calculations and initializations
 150
 151   delay=1;
 152   N1_factor = MaxN1 - IdleN1;
 153   N2_factor = MaxN2 - IdleN2;
 154   OilTemp_degK = Auxiliary->GetTAT_C() + 273.0;
 155   if (IdleFF==-1) IdleFF = pow(MilThrust, 0.2) * 107.0;  // just an estimate
 156
 157   // cout << "ENG POWER:" << EnginePowerRPM_N1->GetValue(1200,90) << endl;
 158
 159   return true;
 160 }
 161
 162 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 163 // The main purpose of Calculate() is to determine what phase the engine should
 164 // be in, then call the corresponding function.
 165
 166 void FGTurboProp::Calculate(void)
 167 {
 168   RunPreFunctions();
 169
 170   TAT = Auxiliary->GetTAT_C();
 171   dt = FDMExec->GetDeltaT() * Propulsion->GetRate();
 172
 173   Throttle = FCS->GetThrottlePos(EngineNumber);
 174
 175   RPM = Thruster->GetRPM() * Thruster->GetGearRatio();
 176   if (thrusterType == FGThruster::ttPropeller) {
 177     ((FGPropeller*)Thruster)->SetAdvance(FCS->GetPropAdvance(EngineNumber));
 178     ((FGPropeller*)Thruster)->SetFeather(FCS->GetPropFeather(EngineNumber));
 179     ((FGPropeller*)Thruster)->SetReverse(Reversed);
 180     if (Reversed) {
 181       ((FGPropeller*)Thruster)->SetReverseCoef(Throttle);
 182     } else {
 183       ((FGPropeller*)Thruster)->SetReverseCoef(0.0);
 184     }
 185
 186     if (Reversed) {
 187       if (Throttle < BetaRangeThrottleEnd) {
 188           Throttle = 0.0;  // idle when in Beta-range
 189       } else {
 190         // when reversed:
 191         Throttle = (Throttle-BetaRangeThrottleEnd)/(1-BetaRangeThrottleEnd) * ReverseMaxPower;
 192       }
 193     }
 194   }
 195
 196   // When trimming is finished check if user wants engine OFF or RUNNING
 197   if ((phase == tpTrim) && (dt > 0)) {
 198     if (Running && !Starved) {
 199       phase = tpRun;
 200       N2 = IdleN2;
 201       N1 = IdleN1;
 202       OilTemp_degK = 366.0;
 203       Cutoff = false;
 204     } else {
 205       phase = tpOff;
 206       Cutoff = true;
 207       Eng_ITT_degC = TAT;
 208       Eng_Temperature = TAT;
 209       OilTemp_degK = TAT+273.15;
 210     }
 211   }
 212
 213   if (!Running && Starter) {
 214     if (phase == tpOff) {
 215       phase = tpSpinUp;
 216       if (StartTime < 0) StartTime=0;
 217     }
 218   }
 219   if (!Running && !Cutoff && (N1 > 15.0)) {
 220     phase = tpStart;
 221     StartTime = -1;
 222   }
 223   if (Cutoff && (phase != tpSpinUp)) phase = tpOff;
 224   if (dt == 0) phase = tpTrim;
 225   if (Starved) phase = tpOff;
 226   if (Condition >= 10) {
 227     phase = tpOff;
 228     StartTime=-1;
 229   }
 230
 231   // limiter intervention wanted?
 232   if (Ielu_max_torque > 0.0) {
 233     double torque = 0.0;
 234
 235     if (thrusterType == FGThruster::ttPropeller) {
 236       torque = ((FGPropeller*)(Thruster))->GetTorque();
 237     } else if (thrusterType == FGThruster::ttRotor) {
 238       torque = ((FGRotor*)(Thruster))->GetTorque();
 239     }
 240
 241     if (Condition < 1) {
 242       if ( abs(torque) > Ielu_max_torque && Throttle >= OldThrottle ) {
 243         Throttle = OldThrottle - 0.1 * dt; //IELU down
 244         Ielu_intervent = true;
 245       } else if ( Ielu_intervent && Throttle >= OldThrottle) {
 246         Throttle = OldThrottle + 0.05 * dt; //IELU up
 247         Ielu_intervent = true;
 248       } else {
 249         Ielu_intervent = false;
 250       }
 251     } else {
 252       Ielu_intervent = false;
 253     }
 254     OldThrottle = Throttle;
 255   }
 256
 257   switch (phase) {
 258     case tpOff:    HP = Off(); break;
 259     case tpRun:    HP = Run(); break;
 260     case tpSpinUp: HP = SpinUp(); break;
 261     case tpStart:  HP = Start(); break;
 262     default: HP = 0;
 263   }
 264
 265   Thruster->Calculate(HP * hptoftlbssec);
 266
 267   RunPostFunctions();
 268 }
 269
 270 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 271
 272 double FGTurboProp::Off(void)
 273 {
 274   double qbar = Auxiliary->Getqbar();
 275   Running = false; EngStarting = false;
 276
 277   FuelFlow_pph = Seek(&FuelFlow_pph, 0, 800.0, 800.0);
 278
 279   //allow the air turn with generator
 280   N1 = ExpSeek(&N1, qbar/15.0, Idle_Max_Delay*2.5, Idle_Max_Delay * 5);
 281
 282   OilTemp_degK = ExpSeek(&OilTemp_degK,273.15 + TAT, 400 , 400);
 283
 284   Eng_Temperature = ExpSeek(&Eng_Temperature,TAT,300,400);
 285   double ITT_goal = ITT_N1->GetValue(N1,0.1) + ((N1>20) ? 0.0 : (20-N1)/20.0 * Eng_Temperature);
 286   Eng_ITT_degC  = ExpSeek(&Eng_ITT_degC,ITT_goal,ITT_Delay,ITT_Delay*1.2);
 287
 288   OilPressure_psi = (N1/100.0*0.25+(0.1-(OilTemp_degK-273.15)*0.1/80.0)*N1/100.0) / 7692.0e-6; //from MPa to psi
 289
 290   ConsumeFuel(); // for possible setting Starved = false when fuel tank
 291                  // is refilled (fuel crossfeed etc.)
 292
 293   if (RPM>5) return -0.012; // friction in engine when propeller spining (estimate)
 294   return 0.0;
 295 }
 296
 297 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 298
 299 double FGTurboProp::Run(void)
 300 {
 301   double thrust = 0.0, EngPower_HP, eff_coef;
 302   Running = true; Starter = false; EngStarting = false;
 303
 304 //---
 305   double old_N1 = N1;
 306   N1 = ExpSeek(&N1, IdleN1 + Throttle * N1_factor, Idle_Max_Delay, Idle_Max_Delay * 2.4);
 307
 308   EngPower_HP = EnginePowerRPM_N1->GetValue(RPM,N1);
 309   EngPower_HP *= EnginePowerVC->GetValue();
 310   if (EngPower_HP > MaxPower) EngPower_HP = MaxPower;
 311
 312   eff_coef = 9.333 - (N1)/12; // 430%Fuel at 60%N1
 313   FuelFlow_pph = PSFC * EngPower_HP * eff_coef;
 314
 315   Eng_Temperature = ExpSeek(&Eng_Temperature,Eng_ITT_degC,300,400);
 316   double ITT_goal = ITT_N1->GetValue((N1-old_N1)*300+N1,1);
 317   Eng_ITT_degC  = ExpSeek(&Eng_ITT_degC,ITT_goal,ITT_Delay,ITT_Delay*1.2);
 318
 319   OilPressure_psi = (N1/100.0*0.25+(0.1-(OilTemp_degK-273.15)*0.1/80.0)*N1/100.0) / 7692.0e-6; //from MPa to psi
 320 //---
 321   EPR = 1.0 + thrust/MilThrust;
 322
 323   OilTemp_degK = Seek(&OilTemp_degK, 353.15, 0.4-N1*0.001, 0.04);
 324
 325   ConsumeFuel();
 326
 327   if (Cutoff) phase = tpOff;
 328   if (Starved) phase = tpOff;
 329
 330   return EngPower_HP;
 331 }
 332
 333 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 334
 335 double FGTurboProp::SpinUp(void)
 336 {
 337   double EngPower_HP;
 338   Running = false; EngStarting = true;
 339   FuelFlow_pph = 0.0;
 340
 341   if (!GeneratorPower) {
 342     EngStarting=false;
 343     phase=tpOff;
 344     StartTime = -1;
 345     return 0.0;
 346   }
 347
 348   N1 = ExpSeek(&N1, StarterN1, Idle_Max_Delay * 6, Idle_Max_Delay * 2.4);
 349
 350   Eng_Temperature = ExpSeek(&Eng_Temperature,TAT,300,400);
 351   double ITT_goal = ITT_N1->GetValue(N1,0.1) + ((N1>20) ? 0.0 : (20-N1)/20.0 * Eng_Temperature);
 352   Eng_ITT_degC  = ExpSeek(&Eng_ITT_degC,ITT_goal,ITT_Delay,ITT_Delay*1.2);
 353
 354   OilTemp_degK = ExpSeek(&OilTemp_degK,273.15 + TAT, 400 , 400);
 355
 356   OilPressure_psi = (N1/100.0*0.25+(0.1-(OilTemp_degK-273.15)*0.1/80.0)*N1/100.0) / 7692.0e-6; //from MPa to psi
 357   NozzlePosition = 1.0;
 358
 359   EngPower_HP = EnginePowerRPM_N1->GetValue(RPM,N1);
 360   EngPower_HP *= EnginePowerVC->GetValue();
 361   if (EngPower_HP > MaxPower) EngPower_HP = MaxPower;
 362
 363   if (StartTime>=0) StartTime+=dt;
 364   if (StartTime > MaxStartingTime && MaxStartingTime > 0) { //start failed due timeout
 365     phase = tpOff;
 366     StartTime = -1;
 367   }
 368
 369   ConsumeFuel(); // for possible setting Starved = false when fuel tank
 370                  // is refilled (fuel crossfeed etc.)
 371
 372   return EngPower_HP;
 373 }
 374
 375 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 376
 377 double FGTurboProp::Start(void)
 378 {
 379   double EngPower_HP = 0.0;
 380   double eff_coef;
 381
 382   EngStarting = false;
 383   if ((N1 > 15.0) && !Starved) {       // minimum 15% N2 needed for start
 384     double old_N1 = N1;
 385     Cranking = true;                   // provided for sound effects signal
 386     if (N1 < IdleN1) {
 387       EngPower_HP = EnginePowerRPM_N1->GetValue(RPM,N1);
 388       EngPower_HP *= EnginePowerVC->GetValue();
 389       if (EngPower_HP > MaxPower) EngPower_HP = MaxPower;
 390       N1 = ExpSeek(&N1, IdleN1*1.1, Idle_Max_Delay*4, Idle_Max_Delay * 2.4);
 391       eff_coef = 9.333 - (N1)/12; // 430%Fuel at 60%N1
 392       FuelFlow_pph = PSFC * EngPower_HP * eff_coef;
 393       Eng_Temperature = ExpSeek(&Eng_Temperature,Eng_ITT_degC,300,400);
 394       double ITT_goal = ITT_N1->GetValue((N1-old_N1)*300+N1,1);
 395       Eng_ITT_degC  = ExpSeek(&Eng_ITT_degC,ITT_goal,ITT_Delay,ITT_Delay*1.2);
 396
 397       OilPressure_psi = (N1/100.0*0.25+(0.1-(OilTemp_degK-273.15)*0.1/80.0)*N1/100.0) / 7692.0e-6; //from MPa to psi
 398       OilTemp_degK = Seek(&OilTemp_degK, 353.15, 0.4-N1*0.001, 0.04);
 399
 400     } else {
 401       phase = tpRun;
 402       Running = true;
 403       Starter = false;
 404       Cranking = false;
 405       FuelFlow_pph = 0;
 406     }
 407   } else {                 // no start if N2 < 15% or Starved
 408     phase = tpOff;
 409     Starter = false;
 410   }
 411
 412   ConsumeFuel();
 413
 414   return EngPower_HP;
 415 }
 416
 417
 418 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 419
 420 double FGTurboProp::CalcFuelNeed(void)
 421 {
 422   double dT = FDMExec->GetDeltaT() * Propulsion->GetRate();
 423   FuelFlowRate = FuelFlow_pph / 3600.0;
 424   FuelExpended = FuelFlowRate * dT;
 425   return FuelExpended;
 426 }
 427
 428 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 429
 430 double FGTurboProp::Seek(double *var, double target, double accel, double decel)
 431 {
 432   double v = *var;
 433   if (v > target) {
 434     v -= dt * decel;
 435     if (v < target) v = target;
 436   } else if (v < target) {
 437     v += dt * accel;
 438     if (v > target) v = target;
 439   }
 440   return v;
 441 }
 442
 443 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 444
 445 double FGTurboProp::ExpSeek(double *var, double target, double accel_tau, double decel_tau)
 446 {
 447 // exponential delay instead of the linear delay used in Seek
 448   double v = *var;
 449   if (v > target) {
 450     v = (v - target) * exp ( -dt / decel_tau) + target;
 451   } else if (v < target) {
 452     v = (target - v) * (1 - exp ( -dt / accel_tau)) + v;
 453   }
 454   return v;
 455 }
 456
 457 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 458
 459 void FGTurboProp::SetDefaults(void)
 460 {
 461 //  Name = "Not defined";
 462   N1 = N2 = 0.0;
 463   HP = 0.0;
 464   Type = etTurboprop;
 465   MilThrust = 10000.0;
 466   IdleN1 = 30.0;
 467   IdleN2 = 60.0;
 468   MaxN1 = 100.0;
 469   MaxN2 = 100.0;
 470   Throttle = 0.0;
 471   InletPosition = 1.0;
 472   NozzlePosition = 1.0;
 473   Reversed = false;
 474   Cutoff = true;
 475   phase = tpOff;
 476   Stalled = false;
 477   Seized = false;
 478   Overtemp = false;
 479   Fire = false;
 480   Eng_ITT_degC = 0.0;
 481
 482   GeneratorPower=true;
 483   Condition = 0;
 484   Ielu_intervent=false;
 485
 486   Idle_Max_Delay = 1.0;
 487
 488   Throttle = OldThrottle = 0.0;
 489   ITT_Delay = 0.05;
 490   ReverseMaxPower = 0.0;
 491   BetaRangeThrottleEnd = 0.0;
 492 }
 493
 494 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 495
 496
 497 string FGTurboProp::GetEngineLabels(const string& delimiter)
 498 {
 499   std::ostringstream buf;
 500
 501   buf << Name << "_N1[" << EngineNumber << "]" << delimiter
 502       << Name << "_N2[" << EngineNumber << "]" << delimiter
 503       << Name << "_PwrAvail[" << EngineNumber << "]" << delimiter
 504       << Thruster->GetThrusterLabels(EngineNumber, delimiter);
 505
 506   return buf.str();
 507 }
 508
 509 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 510
 511 string FGTurboProp::GetEngineValues(const string& delimiter)
 512 {
 513   std::ostringstream buf;
 514
 515   buf << N1 << delimiter
 516       << N2 << delimiter
 517       << HP << delimiter
 518       << Thruster->GetThrusterValues(EngineNumber,delimiter);
 519
 520   return buf.str();
 521 }
 522
 523 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 524
 525 int FGTurboProp::InitRunning(void)
 526 {
 527   FDMExec->SuspendIntegration();
 528   Cutoff=false;
 529   Running=true;
 530   N2=16.0;
 531   Calculate();
 532   FDMExec->ResumeIntegration();
 533   return phase==tpRun;
 534 }
 535
 536 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 537
 538 void FGTurboProp::bindmodel()
 539 {
 540   string property_name, base_property_name;
 541   base_property_name = CreateIndexedPropertyName("propulsion/engine", EngineNumber);
 542   property_name = base_property_name + "/n1";
 543   PropertyManager->Tie( property_name.c_str(), &N1);
 544   // property_name = base_property_name + "/n2";
 545   // PropertyManager->Tie( property_name.c_str(), &N2);
 546   property_name = base_property_name + "/reverser";
 547   PropertyManager->Tie( property_name.c_str(), &Reversed);
 548   property_name = base_property_name + "/power-hp";
 549   PropertyManager->Tie( property_name.c_str(), &HP);
 550   property_name = base_property_name + "/itt-c";
 551   PropertyManager->Tie( property_name.c_str(), &Eng_ITT_degC);
 552   property_name = base_property_name + "/engtemp-c";
 553   PropertyManager->Tie( property_name.c_str(), &Eng_Temperature);
 554   property_name = base_property_name + "/ielu_intervent";
 555   PropertyManager->Tie( property_name.c_str(), &Ielu_intervent);
 556 }
 557
 558 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 559 //    The bitmasked value choices are as follows:
 560 //    unset: In this case (the default) JSBSim would only print
 561 //       out the normally expected messages, essentially echoing
 562 //       the config files as they are read. If the environment
 563 //       variable is not set, debug_lvl is set to 1 internally
 564 //    0: This requests JSBSim not to output any messages
 565 //       whatsoever.
 566 //    1: This value explicity requests the normal JSBSim
 567 //       startup messages
 568 //    2: This value asks for a message to be printed out when
 569 //       a class is instantiated
 570 //    4: When this value is set, a message is displayed when a
 571 //       FGModel object executes its Run() method
 572 //    8: When this value is set, various runtime state variables
 573 //       are printed out periodically
 574 //    16: When set various parameters are sanity checked and
 575 //       a message is printed out when they go out of bounds
 576
 577 void FGTurboProp::Debug(int from)
 578 {
 579   if (debug_lvl <= 0) return;
 580
 581   if (debug_lvl & 1) { // Standard console startup message output
 582     if (from == 0) { // Constructor
 583
 584     }
 585     if (from == 2) { // called from Load()
 586       cout << "\n ****MUJ MOTOR TURBOPROP****\n";
 587       cout << "\n    Engine Name: "         << Name << endl;
 588       cout << "      MilThrust:   "         << MilThrust << endl;
 589       cout << "      IdleN1:      "         << IdleN1 << endl;
 590       cout << "      MaxN1:       "         << MaxN1 << endl;
 591
 592       cout << endl;
 593     }
 594   }
 595   if (debug_lvl & 2 ) { // Instantiation/Destruction notification
 596     if (from == 0) cout << "Instantiated: FGTurboProp" << endl;
 597     if (from == 1) cout << "Destroyed:    FGTurboProp" << endl;
 598   }
 599   if (debug_lvl & 4 ) { // Run() method entry print for FGModel-derived objects
 600   }
 601   if (debug_lvl & 8 ) { // Runtime state variables
 602   }
 603   if (debug_lvl & 16) { // Sanity checking
 604   }
 605   if (debug_lvl & 64) {
 606     if (from == 0) { // Constructor
 607       cout << IdSrc << endl;
 608       cout << IdHdr << endl;
 609     }
 610   }
 611 }
 612 }