JSBSim Flight Dynamics Model 1.2.2 (22 Mar 2025)
An Open Source Flight Dynamics and Control Software Library in C++
Loading...
Searching...
No Matches
FGTurbine.cpp
1/*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2
3 Module: FGTurbine.cpp
4 Author: David Culp
5 Date started: 03/11/2003
6 Purpose: This module models a turbine engine.
7
8 ------------- Copyright (C) 2003 David Culp (daveculp@cox.net) ---------
9
10 This program is free software; you can redistribute it and/or modify it under
11 the terms of the GNU Lesser General Public License as published by the Free Software
12 Foundation; either version 2 of the License, or (at your option) any later
13 version.
14
15 This program is distributed in the hope that it will be useful, but WITHOUT
16 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
17 FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more
18 details.
19
20 You should have received a copy of the GNU Lesser General Public License along with
21 this program; if not, write to the Free Software Foundation, Inc., 59 Temple
22 Place - Suite 330, Boston, MA 02111-1307, USA.
23
24 Further information about the GNU Lesser General Public License can also be found on
25 the world wide web at http://www.gnu.org.
26
27FUNCTIONAL DESCRIPTION
28--------------------------------------------------------------------------------
29
30This class descends from the FGEngine class and models a turbine engine based
31on parameters given in the engine config file for this class
32
33HISTORY
34--------------------------------------------------------------------------------
3503/11/2003 DPC Created
3609/08/2003 DPC Changed Calculate() and added engine phases
37
38%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
39INCLUDES
40%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/
41
42#include <iostream>
43#include <sstream>
44
45#include "FGFDMExec.h"
46#include "math/FGFunction.h"
47#include "math/FGRealValue.h"
48#include "FGTurbine.h"
49#include "FGThruster.h"
50#include "input_output/FGXMLElement.h"
51#include "input_output/string_utilities.h"
52
53using namespace std;
54
55namespace JSBSim {
56
57/*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
58CLASS IMPLEMENTATION
59%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/
60
61FGTurbine::FGTurbine(FGFDMExec* exec, Element *el, int engine_number, struct Inputs& input)
62 : FGEngine(engine_number, input), FDMExec(exec)
63{
64 Type = etTurbine;
65
66 MilThrust = MaxThrust = 10000.0;
67 TSFC = std::make_unique<FGSimplifiedTSFC>(this, 0.8);
68 ATSFC = std::make_unique<FGRealValue>(1.7);
69 IdleN1 = 30.0;
70 IdleN2 = 60.0;
71 MaxN1 = MaxN2 = 100.0;
72 Augmented = AugMethod = Injected = 0;
73 BypassRatio = BleedDemand = 0.0;
74 IdleThrustLookup = MilThrustLookup = MaxThrustLookup = InjectionLookup = nullptr;
75 N1_spinup = 1.0; N2_spinup = 3.0; IgnitionN1 = 5.21; IgnitionN2 = 25.18; N1_start_rate = 1.4; N2_start_rate = 2.0;
76 N1_spindown = 2.0; N2_spindown = 2.0;
77 InjectionTime = 30.0;
78 InjectionTimer = InjWaterNorm = 0.0;
79 EPR = 1.0;
80 disableWindmill = false;
81
82 Load(exec, el);
83 Debug(0);
84}
85
86//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
87
88void FGTurbine::ResetToIC(void)
89{
90 FGEngine::ResetToIC();
91
92 N1 = N2 = InjN1increment = InjN2increment = 0.0;
93 N2norm = 0.0;
94 correctedTSFC = TSFC->GetValue();
95 AugmentCmd = InjWaterNorm = 0.0;
96 InletPosition = NozzlePosition = 1.0;
97 Stalled = Seized = Overtemp = Fire = Augmentation = Injection = Reversed = false;
98 Cutoff = true;
99 phase = tpOff;
100 EGT_degC = in.TAT_c;
101 OilTemp_degK = in.TAT_c + 273.0;
102}
103
104//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
105// The main purpose of Calculate() is to determine what phase the engine should
106// be in, then call the corresponding function.
107
108void FGTurbine::Calculate(void)
109{
110 double thrust;
111
112 RunPreFunctions();
113
114 ThrottlePos = in.ThrottlePos[EngineNumber];
115
116 if (ThrottlePos > 1.0) {
117 AugmentCmd = ThrottlePos - 1.0;
118 ThrottlePos -= AugmentCmd;
119 } else {
120 AugmentCmd = 0.0;
121 }
122
123 // When trimming is finished check if user wants engine OFF or RUNNING
124 if ((phase == tpTrim) && (in.TotalDeltaT > 0)) {
125 if (Running && !Starved) {
126 phase = tpRun;
127 N1_factor = MaxN1 - IdleN1;
128 N2_factor = MaxN2 - IdleN2;
129 N2 = IdleN2 + ThrottlePos * N2_factor;
130 N1 = IdleN1 + ThrottlePos * N1_factor;
131 OilTemp_degK = 366.0;
132 Cutoff = false;
133 } else {
134 phase = tpOff;
135 Cutoff = true;
136 EGT_degC = in.TAT_c;
137 }
138 }
139
140 if (!Running && Cutoff && Starter) {
141 if (phase == tpOff) phase = tpSpinUp;
142 }
143
144 // start
145 if ((Starter == true) || (in.qbar > 30.0)) {
146 if (!Running && !Cutoff && (N2 > 15.0)) phase = tpStart;
147 }
148
149 if (Cutoff && (phase != tpSpinUp)) phase = tpOff;
150 if (in.TotalDeltaT == 0) phase = tpTrim;
151 if (Starved) phase = tpOff;
152 if (Stalled) phase = tpStall;
153 if (Seized) phase = tpSeize;
154
155 switch (phase) {
156 case tpOff: thrust = Off(); break;
157 case tpRun: thrust = Run(); break;
158 case tpSpinUp: thrust = SpinUp(); break;
159 case tpStart: thrust = Start(); break;
160 case tpStall: thrust = Stall(); break;
161 case tpSeize: thrust = Seize(); break;
162 case tpTrim: thrust = Trim(); break;
163 default: thrust = Off();
164 }
165
166 Thruster->Calculate(thrust); // allow thruster to modify thrust (i.e. reversing)
167
168 RunPostFunctions();
169}
170
171//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
172
173double FGTurbine::Off(void)
174{
175 Running = false;
176 FuelFlow_pph = Seek(&FuelFlow_pph, 0, 1000.0, 10000.0);
177 // some engines have inlets that close when they are off. So, if a flag is true disable windmilling
178 if (disableWindmill == false) {
179 // Need a small non-zero increment for acceleration otherwise acceleration will be 0 if N1 = 0
180 N1 = Seek(&N1, in.qbar/10.0, N1/2.0 + 0.1, N1/N1_spindown);
181 N2 = Seek(&N2, in.qbar/15.0, N2/2.0 + 0.1, N2/N2_spindown);
182 } else {
183 N1 = Seek(&N1, 0, N1/2.0, N1/N1_spindown);
184 N2 = Seek(&N2, 0, N2/2.0, N2/N2_spindown);
185 }
186 EGT_degC = Seek(&EGT_degC, in.TAT_c, 11.7, 7.3);
187 OilTemp_degK = Seek(&OilTemp_degK, in.TAT_c + 273.0, 0.2, 0.2);
188 OilPressure_psi = N2 * 0.62;
189 NozzlePosition = Seek(&NozzlePosition, 1.0, 0.8, 0.8);
190 EPR = Seek(&EPR, 1.0, 0.2, 0.2);
191 Augmentation = false;
192 return 0.0;
193}
194
195//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
196
197double FGTurbine::Run()
198{
199 double idlethrust, milthrust, thrust;
200
201 idlethrust = MilThrust * IdleThrustLookup->GetValue();
202 milthrust = (MilThrust - idlethrust) * MilThrustLookup->GetValue();
203
204 Running = true;
205 Starter = false;
206
207 N1_factor = MaxN1 - IdleN1;
208 N2_factor = MaxN2 - IdleN2;
209 if ((Injected == 1) && Injection && (InjWaterNorm > 0)) {
210 N1_factor += InjN1increment;
211 N2_factor += InjN2increment;
212 }
213 N2 = Seek(&N2, IdleN2 + ThrottlePos * N2_factor,
214 N2SpoolUp->GetValue(), N2SpoolDown->GetValue());
215 N1 = Seek(&N1, IdleN1 + ThrottlePos * N1_factor,
216 N1SpoolUp->GetValue(), N1SpoolDown->GetValue());
217 N2norm = (N2 - IdleN2) / N2_factor;
218 thrust = idlethrust + (milthrust * N2norm * N2norm);
219 EGT_degC = in.TAT_c + 363.1 + ThrottlePos * 357.1;
220 OilPressure_psi = N2 * 0.62;
221 OilTemp_degK = Seek(&OilTemp_degK, 366.0, 1.2, 0.1);
222
223 if (!Augmentation) {
224 correctedTSFC = TSFC->GetValue();
225 FuelFlow_pph = Seek(&FuelFlow_pph, thrust * correctedTSFC, 1000.0, 10000.0);
226 if (FuelFlow_pph < IdleFF) FuelFlow_pph = IdleFF;
227 NozzlePosition = Seek(&NozzlePosition, 1.0 - N2norm, 0.8, 0.8);
228 thrust = thrust * (1.0 - BleedDemand);
229 EPR = 1.0 + thrust/MilThrust;
230 }
231
232 if (AugMethod == 1) {
233 if ((ThrottlePos > 0.99) && (N2 > 97.0)) {Augmentation = true;}
234 else {Augmentation = false;}
235 }
236
237 if ((Augmented == 1) && Augmentation && (AugMethod < 2)) {
238 thrust = MaxThrustLookup->GetValue() * MaxThrust;
239 FuelFlow_pph = Seek(&FuelFlow_pph, thrust * ATSFC->GetValue(), 5000.0, 10000.0);
240 NozzlePosition = Seek(&NozzlePosition, 1.0, 0.8, 0.8);
241 }
242
243 if (AugMethod == 2) {
244 if (AugmentCmd > 0.0) {
245 Augmentation = true;
246 double tdiff = (MaxThrust * MaxThrustLookup->GetValue()) - thrust;
247 thrust += (tdiff * std::min(AugmentCmd, 1.0));
248 FuelFlow_pph = Seek(&FuelFlow_pph, thrust * ATSFC->GetValue(), 5000.0, 10000.0);
249 NozzlePosition = Seek(&NozzlePosition, 1.0, 0.8, 0.8);
250 } else {
251 Augmentation = false;
252 }
253 }
254
255 if ((Injected == 1) && Injection && (InjWaterNorm > 0.0)) {
256 InjectionTimer += in.TotalDeltaT;
257 if (InjectionTimer < InjectionTime) {
258 thrust = thrust * InjectionLookup->GetValue();
259 InjWaterNorm = 1.0 - (InjectionTimer/InjectionTime);
260 } else {
261 Injection = false;
262 InjWaterNorm = 0.0;
263 }
264 }
265
266 if (Cutoff) phase = tpOff;
267 if (Starved) phase = tpOff;
268
269 return thrust;
270}
271
272//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
273
274double FGTurbine::SpinUp(void)
275{
276 Running = false;
277 FuelFlow_pph = 0.0;
278 N2 = Seek(&N2, IgnitionN2, N2_spinup, N2/2.0);
279 N1 = Seek(&N1, IgnitionN1, N1_spinup, N1/2.0);
280 EGT_degC = Seek(&EGT_degC, in.TAT_c, 11.7, 7.3);
281 OilPressure_psi = N2 * 0.62;
282 OilTemp_degK = Seek(&OilTemp_degK, in.TAT_c + 273.0, 0.2, 0.2);
283 EPR = 1.0;
284 NozzlePosition = 1.0;
285 if (Starter == false) phase = tpOff;
286 return 0.0;
287}
288
289//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
290
291double FGTurbine::Start(void)
292{
293 if ((N2 > 15.0) && !Starved) { // minimum 15% N2 needed for start
294 Cranking = true; // provided for sound effects signal
295 if (N2 < IdleN2) {
296 N2 = Seek(&N2, IdleN2, N2_start_rate, N2/2.0);
297 N1 = Seek(&N1, IdleN1, N1_start_rate, N1/2.0);
298 EGT_degC = Seek(&EGT_degC, in.TAT_c + 363.1, 21.3, 7.3);
299 FuelFlow_pph = IdleFF * N2 / IdleN2;
300 OilPressure_psi = N2 * 0.62;
301 if ((Starter == false) && (in.qbar < 30.0)) phase = tpOff; // aborted start
302 }
303 else {
304 phase = tpRun;
305 Running = true;
306 Starter = false;
307 Cranking = false;
308 }
309 }
310 else { // no start if N2 < 15%
311 phase = tpOff;
312 Starter = false;
313 }
314
315 return 0.0;
316}
317
318//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
319
320double FGTurbine::Stall(void)
321{
322 EGT_degC = in.TAT_c + 903.14;
323 FuelFlow_pph = IdleFF;
324 N1 = Seek(&N1, in.qbar/10.0, 0, N1/10.0);
325 N2 = Seek(&N2, in.qbar/15.0, 0, N2/10.0);
326 if (ThrottlePos < 0.01) {
327 phase = tpRun; // clear the stall with throttle to idle
328 Stalled = false;
329 }
330 return 0.0;
331}
332
333//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
334
335double FGTurbine::Seize(void)
336{
337 N2 = 0.0;
338 N1 = Seek(&N1, in.qbar/20.0, 0, N1/15.0);
339 FuelFlow_pph = Cutoff ? 0.0 : IdleFF;
340 OilPressure_psi = 0.0;
341 OilTemp_degK = Seek(&OilTemp_degK, in.TAT_c + 273.0, 0, 0.2);
342 Running = false;
343 return 0.0;
344}
345
346//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
347
348double FGTurbine::Trim()
349{
350 double idlethrust = MilThrust * IdleThrustLookup->GetValue();
351 double milthrust = (MilThrust - idlethrust) * MilThrustLookup->GetValue();
352 double N2 = IdleN2 + ThrottlePos * N2_factor;
353 double N2norm = (N2 - IdleN2) / N2_factor;
354 double thrust = (idlethrust + (milthrust * N2norm * N2norm))
355 * (1.0 - BleedDemand);
356
357 if (AugMethod == 1) {
358 if ((ThrottlePos > 0.99) && (N2 > 97.0)) {Augmentation = true;}
359 else {Augmentation = false;}
360 }
361
362 if ((Augmented == 1) && Augmentation && (AugMethod < 2)) {
363 thrust = MaxThrust * MaxThrustLookup->GetValue();
364 }
365
366 if (AugMethod == 2) {
367 if (AugmentCmd > 0.0) {
368 double tdiff = (MaxThrust * MaxThrustLookup->GetValue()) - thrust;
369 thrust += (tdiff * std::min(AugmentCmd, 1.0));
370 }
371 }
372
373 if ((Injected == 1) && Injection) {
374 thrust = thrust * InjectionLookup->GetValue();
375 }
376
377 return thrust;
378}
379
380//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
381
382double FGTurbine::CalcFuelNeed(void)
383{
384 FuelFlowRate = FuelFlow_pph / 3600.0; // Calculates flow in lbs/sec from lbs/hr
385 FuelExpended = FuelFlowRate * in.TotalDeltaT; // Calculates fuel expended in this time step
386 if (!Starved) FuelUsedLbs += FuelExpended;
387 return FuelExpended;
388}
389
390//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
391
392double FGTurbine::GetPowerAvailable(void) {
393 if( ThrottlePos <= 0.77 )
394 return 64.94*ThrottlePos;
395 else
396 return 217.38*ThrottlePos - 117.38;
397}
398
399//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
400
401double FGTurbine::Seek(double *var, double target, double accel, double decel) {
402 double v = *var;
403 if (v > target) {
404 v -= in.TotalDeltaT * decel;
405 if (v < target) v = target;
406 } else if (v < target) {
407 v += in.TotalDeltaT * accel;
408 if (v > target) v = target;
409 }
410 return v;
411}
412
413//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
414
415bool FGTurbine::Load(FGFDMExec* exec, Element *el)
416{
417 Element* function_element = el->FindElement("function");
418
419 while(function_element) {
420 string name = function_element->GetAttributeValue("name");
421 if (name == "IdleThrust" || name == "MilThrust" || name == "AugThrust"
422 || name == "Injection" || name == "N1SpoolUp" || name == "N1SpoolDown"
423 || name == "N2SpoolUp" || name == "N2SpoolDown")
424 function_element->SetAttributeValue("name", string("propulsion/engine[#]/") + name);
425
426 function_element = el->FindNextElement("function");
427 }
428
429 FGEngine::Load(exec, el);
430
431 ResetToIC();
432
433 if (el->FindElement("milthrust"))
434 MilThrust = el->FindElementValueAsNumberConvertTo("milthrust","LBS");
435 if (el->FindElement("maxthrust"))
436 MaxThrust = el->FindElementValueAsNumberConvertTo("maxthrust","LBS");
437 if (el->FindElement("bypassratio"))
438 BypassRatio = el->FindElementValueAsNumber("bypassratio");
439 if (el->FindElement("bleed"))
440 BleedDemand = el->FindElementValueAsNumber("bleed");
441 if (el->FindElement("ignitionn1"))
442 IgnitionN1 = el->FindElementValueAsNumber("ignitionn1");
443 if (el->FindElement("ignitionn2"))
444 IgnitionN2 = el->FindElementValueAsNumber("ignitionn2");
445 if (el->FindElement("idlen1"))
446 IdleN1 = el->FindElementValueAsNumber("idlen1");
447 if (el->FindElement("idlen2"))
448 IdleN2 = el->FindElementValueAsNumber("idlen2");
449 if (el->FindElement("maxn1"))
450 MaxN1 = el->FindElementValueAsNumber("maxn1");
451 if (el->FindElement("maxn2"))
452 MaxN2 = el->FindElementValueAsNumber("maxn2");
453 if (el->FindElement("n1spinup"))
454 N1_spinup = el->FindElementValueAsNumber("n1spinup");
455 if (el->FindElement("n2spinup"))
456 N2_spinup = el->FindElementValueAsNumber("n2spinup");
457 if (el->FindElement("n1startrate"))
458 N1_start_rate = el->FindElementValueAsNumber("n1startrate");
459 if (el->FindElement("n2startrate"))
460 N2_start_rate = el->FindElementValueAsNumber("n2startrate");
461 if (el->FindElement("n1spindown"))
462 N1_spindown = el->FindElementValueAsNumber("n1spindown");
463 if (el->FindElement("n2spindown"))
464 N2_spindown = el->FindElementValueAsNumber("n2spindown");
465 if (el->FindElement("augmented"))
466 Augmented = (int)el->FindElementValueAsNumber("augmented");
467 if (el->FindElement("augmethod"))
468 AugMethod = (int)el->FindElementValueAsNumber("augmethod");
469 if (el->FindElement("injected"))
470 Injected = (int)el->FindElementValueAsNumber("injected");
471 if (el->FindElement("injection-time")){
472 InjectionTime = el->FindElementValueAsNumber("injection-time");
473 InjWaterNorm =1.0;
474 }
475 if (el->FindElement("injection-N1-inc"))
476 InjN1increment = el->FindElementValueAsNumber("injection-N1-inc");
477 if (el->FindElement("injection-N2-inc"))
478 InjN2increment = el->FindElementValueAsNumber("injection-N2-inc");
479 if (el->FindElement("disable-windmill"))
480 disableWindmill = el->FindElementValueAsBoolean("disable-windmill");
481 string property_prefix = CreateIndexedPropertyName("propulsion/engine", EngineNumber);
482
483 IdleThrustLookup = GetPreFunction(property_prefix+"/IdleThrust");
484 MilThrustLookup = GetPreFunction(property_prefix+"/MilThrust");
485 MaxThrustLookup = GetPreFunction(property_prefix+"/AugThrust");
486 InjectionLookup = GetPreFunction(property_prefix+"/Injection");
487
488 JSBSim::Element* tsfcElement = el->FindElement("tsfc");
489 if (tsfcElement) {
490 string value = tsfcElement->GetDataLine();
491 try {
492 TSFC = std::make_unique<FGSimplifiedTSFC>(this, atof_locale_c(value));
493 } catch (InvalidNumber&) {
494 TSFC = std::make_unique<FGFunction>(FDMExec, tsfcElement, to_string(EngineNumber));
495 }
496 }
497
498 JSBSim::Element* atsfcElement = el->FindElement("atsfc");
499 if (atsfcElement) {
500 string value = atsfcElement->GetDataLine();
501 try {
502 ATSFC = std::make_unique<FGRealValue>(atof_locale_c(value));
503 } catch (InvalidNumber&) {
504 ATSFC = std::make_unique<FGFunction>(FDMExec, atsfcElement, to_string((int)EngineNumber));
505 }
506 }
507
508 // Pre-calculations and initializations
509 N1SpoolUp = GetPreFunction(property_prefix+"/N1SpoolUp");
510 if (!N1SpoolUp)
511 N1SpoolUp = std::make_shared<FGSpoolUp>(this, BypassRatio, 1.0);
512
513 N1SpoolDown = GetPreFunction(property_prefix+"/N1SpoolDown");
514 if (!N1SpoolDown)
515 N1SpoolDown = std::make_shared<FGSpoolUp>(this, BypassRatio, 2.4);
516
517 N2SpoolUp = GetPreFunction(property_prefix+"/N2SpoolUp");
518 if (!N2SpoolUp)
519 N2SpoolUp = std::make_shared<FGSpoolUp>(this, BypassRatio, 1.0);
520
521 N2SpoolDown = GetPreFunction(property_prefix+"/N2SpoolDown");
522 if (!N2SpoolDown)
523 N2SpoolDown = std::make_shared<FGSpoolUp>(this, BypassRatio, 3.0);
524
525 N1_factor = MaxN1 - IdleN1;
526 N2_factor = MaxN2 - IdleN2;
527 OilTemp_degK = in.TAT_c + 273.0;
528 IdleFF = pow(MilThrust, 0.2) * 107.0; // just an estimate
529
530 bindmodel(exec->GetPropertyManager().get());
531 return true;
532}
533
534//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
535
536string FGTurbine::GetEngineLabels(const string& delimiter)
537{
538 std::ostringstream buf;
539
540 buf << Name << "_N1[" << EngineNumber << "]" << delimiter
541 << Name << "_N2[" << EngineNumber << "]" << delimiter
542 << Thruster->GetThrusterLabels(EngineNumber, delimiter);
543
544 return buf.str();
545}
546
547//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
548
549string FGTurbine::GetEngineValues(const string& delimiter)
550{
551 std::ostringstream buf;
552
553 buf << N1 << delimiter
554 << N2 << delimiter
555 << Thruster->GetThrusterValues(EngineNumber, delimiter);
556
557 return buf.str();
558}
559
560//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
561
562void FGTurbine::bindmodel(FGPropertyManager* PropertyManager)
563{
564 string property_name, base_property_name;
565 base_property_name = CreateIndexedPropertyName("propulsion/engine", EngineNumber);
566 property_name = base_property_name + "/n1";
567 PropertyManager->Tie( property_name.c_str(), &N1);
568 property_name = base_property_name + "/n2";
569 PropertyManager->Tie( property_name.c_str(), &N2);
570 property_name = base_property_name + "/injection_cmd";
571 PropertyManager->Tie( property_name.c_str(), this,
572 &FGTurbine::GetInjection, &FGTurbine::SetInjection);
573 property_name = base_property_name + "/seized";
574 PropertyManager->Tie( property_name.c_str(), &Seized);
575 property_name = base_property_name + "/stalled";
576 PropertyManager->Tie( property_name.c_str(), &Stalled);
577 property_name = base_property_name + "/bleed-factor";
578 PropertyManager->Tie( property_name.c_str(), this, &FGTurbine::GetBleedDemand, &FGTurbine::SetBleedDemand);
579 property_name = base_property_name + "/MaxN1";
580 PropertyManager->Tie( property_name.c_str(), this,
581 &FGTurbine::GetMaxN1, &FGTurbine::SetMaxN1);
582 property_name = base_property_name + "/MaxN2";
583 PropertyManager->Tie( property_name.c_str(), this,
584 &FGTurbine::GetMaxN2, &FGTurbine::SetMaxN2);
585 property_name = base_property_name + "/InjectionTimer";
586 PropertyManager->Tie( property_name.c_str(), this,
587 &FGTurbine::GetInjectionTimer, &FGTurbine::SetInjectionTimer);
588 property_name = base_property_name + "/InjWaterNorm";
589 PropertyManager->Tie( property_name.c_str(), this,
590 &FGTurbine::GetInjWaterNorm, &FGTurbine::SetInjWaterNorm);
591 property_name = base_property_name + "/InjN1increment";
592 PropertyManager->Tie( property_name.c_str(), this,
593 &FGTurbine::GetInjN1increment, &FGTurbine::SetInjN1increment);
594 property_name = base_property_name + "/InjN2increment";
595 PropertyManager->Tie( property_name.c_str(), this,
596 &FGTurbine::GetInjN2increment, &FGTurbine::SetInjN2increment);
597 property_name = base_property_name + "/atsfc";
598 PropertyManager->Tie(property_name.c_str(), ATSFC.get(), &FGParameter::GetValue);
599 property_name = base_property_name + "/tsfc";
600 PropertyManager->Tie(property_name.c_str(), &correctedTSFC);
601 auto node = PropertyManager->GetNode(property_name.c_str(), false);
602 node->setAttribute(SGPropertyNode::WRITE, false);
603}
604
605//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
606
607int FGTurbine::InitRunning(void)
608{
609 FDMExec->SuspendIntegration();
610 Cutoff=false;
611 Running=true;
612 N1_factor = MaxN1 - IdleN1;
613 N2_factor = MaxN2 - IdleN2;
614 N2 = IdleN2 + ThrottlePos * N2_factor;
615 N1 = IdleN1 + ThrottlePos * N1_factor;
616 Calculate();
617 FDMExec->ResumeIntegration();
618 return phase=tpRun;
619}
620
621//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
622// The bitmasked value choices are as follows:
623// unset: In this case (the default) JSBSim would only print
624// out the normally expected messages, essentially echoing
625// the config files as they are read. If the environment
626// variable is not set, debug_lvl is set to 1 internally
627// 0: This requests JSBSim not to output any messages
628// whatsoever.
629// 1: This value explicity requests the normal JSBSim
630// startup messages
631// 2: This value asks for a message to be printed out when
632// a class is instantiated
633// 4: When this value is set, a message is displayed when a
634// FGModel object executes its Run() method
635// 8: When this value is set, various runtime state variables
636// are printed out periodically
637// 16: When set various parameters are sanity checked and
638// a message is printed out when they go out of bounds
639
640void FGTurbine::Debug(int from)
641{
642 if (debug_lvl <= 0) return;
643
644 if (debug_lvl & 1) { // Standard console startup message output
645 if (from == 0) { // Constructor
646
647 }
648 if (from == 2) { // called from Load()
649 cout << "\n Engine Name: " << Name << endl;
650 cout << " MilThrust: " << MilThrust << endl;
651 cout << " MaxThrust: " << MaxThrust << endl;
652 cout << " BypassRatio: " << BypassRatio << endl;
653 cout << " TSFC: " << TSFC->GetValue() << endl;
654 cout << " ATSFC: " << ATSFC->GetValue() << endl;
655 cout << " IdleN1: " << IdleN1 << endl;
656 cout << " IdleN2: " << IdleN2 << endl;
657 cout << " MaxN1: " << MaxN1 << endl;
658 cout << " MaxN2: " << MaxN2 << endl;
659 cout << " Augmented: " << Augmented << endl;
660 cout << " AugMethod: " << AugMethod << endl;
661 cout << " Injected: " << Injected << endl;
662 cout << " MinThrottle: " << MinThrottle << endl;
663
664 cout << endl;
665 }
666 }
667 if (debug_lvl & 2 ) { // Instantiation/Destruction notification
668 if (from == 0) cout << "Instantiated: FGTurbine" << endl;
669 if (from == 1) cout << "Destroyed: FGTurbine" << endl;
670 }
671 if (debug_lvl & 4 ) { // Run() method entry print for FGModel-derived objects
672 }
673 if (debug_lvl & 8 ) { // Runtime state variables
674 }
675 if (debug_lvl & 16) { // Sanity checking
676 }
677 if (debug_lvl & 64) {
678 if (from == 0) { // Constructor
679 }
680 }
681}
682}
JSBSim::Element::FindElement
Element * FindElement(const std::string &el="")
Searches for a specified element.
Definition FGXMLElement.cpp:396
JSBSim::Element::SetAttributeValue
bool SetAttributeValue(const std::string &key, const std::string &value)
Modifies an attribute.
Definition FGXMLElement.cpp:275
JSBSim::Element::FindElementValueAsBoolean
bool FindElementValueAsBoolean(const std::string &el="")
Searches for the named element and returns the data belonging to it as a bool.
Definition FGXMLElement.cpp:453
JSBSim::Element::GetAttributeValue
std::string GetAttributeValue(const std::string &key)
Retrieves an attribute.
Definition FGXMLElement.cpp:267
JSBSim::Element::GetDataLine
std::string GetDataLine(unsigned int i=0)
Gets a line of data belonging to an element.
Definition FGXMLElement.cpp:340
JSBSim::Element::FindElementValueAsNumberConvertTo
double FindElementValueAsNumberConvertTo(const std::string &el, const std::string &target_units)
Searches for the named element and converts and returns the data belonging to it.
Definition FGXMLElement.cpp:487
JSBSim::Element::FindNextElement
Element * FindNextElement(const std::string &el="")
Searches for the next element as specified.
Definition FGXMLElement.cpp:414
JSBSim::Element::FindElementValueAsNumber
double FindElementValueAsNumber(const std::string &el="")
Searches for the named element and returns the data belonging to it as a number.
Definition FGXMLElement.cpp:436
JSBSim::FGEngine
Base class for all engines.
Definition FGEngine.h:104
JSBSim::FGEngine::ResetToIC
virtual void ResetToIC(void)
Resets the Engine parameters to the initial conditions.
Definition FGEngine.cpp:77
JSBSim::FGFDMExec
Encapsulates the JSBSim simulation executive.
Definition FGFDMExec.h:184
JSBSim::FGFDMExec::SuspendIntegration
void SuspendIntegration(void)
Suspends the simulation and sets the delta T to zero.
Definition FGFDMExec.h:555
JSBSim::FGFDMExec::GetPropertyManager
std::shared_ptr< FGPropertyManager > GetPropertyManager(void) const
Returns a pointer to the property manager object.
Definition FGFDMExec.h:421
JSBSim::FGFDMExec::ResumeIntegration
void ResumeIntegration(void)
Resumes the simulation by resetting delta T to the correct value.
Definition FGFDMExec.h:558
JSBSim::FGModelFunctions::GetPreFunction
std::shared_ptr< FGFunction > GetPreFunction(const std::string &name)
Get one of the "pre" function.
Definition FGModelFunctions.cpp:131
JSBSim::FGTurbine::Seek
double Seek(double *var, double target, double accel, double decel)
A lag filter.
Definition FGTurbine.cpp:401
JSBSim::FGTurbine::Calculate
void Calculate(void)
Calculates the thrust of the engine, and other engine functions.
Definition FGTurbine.cpp:108
JSBSim::FGTurbine::ResetToIC
void ResetToIC(void)
Resets the Engine parameters to the initial conditions.
Definition FGTurbine.cpp:88
JSBSim::FGTurbine::CalcFuelNeed
double CalcFuelNeed(void)
The fuel need is calculated based on power levels and flow rate for that power level.
Definition FGTurbine.cpp:382
JSBSim::FGTurbine::FGTurbine
FGTurbine(FGFDMExec *Executive, Element *el, int engine_number, struct Inputs &input)
Constructor.
Definition FGTurbine.cpp:61