]> git.mxchange.org Git - flightgear.git/blobdiff - src/FDM/JSBSim/models/propulsion/FGNozzle.cpp
Andreas Gaeb: fix #222 (JSBSIm reset problems)
[flightgear.git] / src / FDM / JSBSim / models / propulsion / FGNozzle.cpp
index bb5d91955a9e016c1a850a313210c3b3159cccc4..f1d0fb12e1911e07bc223f71fbe70d41656c6bf2 100644 (file)
@@ -5,23 +5,23 @@
  Date started: 08/24/00
  Purpose:      Encapsulates the nozzle object
 
- ------------- Copyright (C) 2000  Jon S. Berndt (jsb@hal-pc.org) -------------
+ ------------- Copyright (C) 2000  Jon S. Berndt (jon@jsbsim.org) -------------
 
  This program is free software; you can redistribute it and/or modify it under
- the terms of the GNU General Public License as published by the Free Software
+ the terms of the GNU Lesser General Public License as published by the Free Software
  Foundation; either version 2 of the License, or (at your option) any later
  version.
 
  This program is distributed in the hope that it will be useful, but WITHOUT
  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
- FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
+ FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public License for more
  details.
 
- You should have received a copy of the GNU General Public License along with
+ You should have received a copy of the GNU Lesser General Public License along with
  this program; if not, write to the Free Software Foundation, Inc., 59 Temple
  Place - Suite 330, Boston, MA  02111-1307, USA.
 
- Further information about the GNU General Public License can also be found on
+ Further information about the GNU Lesser General Public License can also be found on
  the world wide web at http://www.gnu.org.
 
 FUNCTIONAL DESCRIPTION
@@ -35,14 +35,19 @@ HISTORY
 INCLUDES
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/
 
+#include <iostream>
 #include <sstream>
+#include <cstdlib>
 
 #include "FGNozzle.h"
-#include <models/FGAtmosphere.h>
+#include "models/FGAtmosphere.h"
+#include "input_output/FGXMLElement.h"
+
+using namespace std;
 
 namespace JSBSim {
 
-static const char *IdSrc = "$Id$";
+static const char *IdSrc = "$Id: FGNozzle.cpp,v 1.13 2009/10/26 03:49:58 jberndt Exp $";
 static const char *IdHdr = ID_NOZZLE;
 
 /*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -50,44 +55,26 @@ CLASS IMPLEMENTATION
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/
 
 
-FGNozzle::FGNozzle(FGFDMExec* FDMExec, Element* nozzle_element, int num) : FGThruster(FDMExec)
+FGNozzle::FGNozzle(FGFDMExec* FDMExec, Element* nozzle_element, int num)
+                    : FGThruster(FDMExec, nozzle_element, num)
 {
-
-  if (nozzle_element->FindElement("pe"))
-    PE = nozzle_element->FindElementValueAsNumberConvertTo("pe", "PSF");
+  if (nozzle_element->FindElement("area"))
+    Area = nozzle_element->FindElementValueAsNumberConvertTo("area", "FT2");
   else {
-    cerr << "Fatal Error: Nozzle exit pressure must be given in nozzle config file." << endl;
+    cerr << "Fatal Error: Nozzle exit area must be given in nozzle config file." << endl;
     exit(-1);
   }
-  if (nozzle_element->FindElement("expr"))
-    ExpR = nozzle_element->FindElementValueAsNumber("expr");
-  else {
-    cerr << "Fatal Error: Nozzle expansion ratio must be given in nozzle config file." << endl;
-    exit(-1);
-  }
-  if (nozzle_element->FindElement("nzl_eff"))
-    nzlEff = nozzle_element->FindElementValueAsNumber("nzl_eff");
-  else {
-    cerr << "Fatal Error: Nozzle efficiency must be given in nozzle config file." << endl;
-    exit(-1);
-  }
-  if (nozzle_element->FindElement("diam"))
-    Diameter = nozzle_element->FindElementValueAsNumberConvertTo("diam", "FT");
+/*
+  if (nozzle_element->FindElement("pe"))
+    PE = nozzle_element->FindElementValueAsNumberConvertTo("pe", "PSF");
   else {
-    cerr << "Fatal Error: Nozzle diameter must be given in nozzle config file." << endl;
+    cerr << "Fatal Error: Nozzle exit pressure must be given in nozzle config file." << endl;
     exit(-1);
   }
-
+*/
   Thrust = 0;
-  ReverserAngle = 0.0;
   Type = ttNozzle;
-  Area2 = (Diameter*Diameter/4.0)*M_PI;
-  AreaT = Area2/ExpR;
-
-//  char property_name[80];
-//  snprintf(property_name, 80, "propulsion/c-thrust[%u]", EngineNum);
-//  PropertyManager->Tie( property_name, &ThrustCoeff );
-
+  
   Debug(0);
 }
 
@@ -95,40 +82,28 @@ FGNozzle::FGNozzle(FGFDMExec* FDMExec, Element* nozzle_element, int num) : FGThr
 
 FGNozzle::~FGNozzle()
 {
-//  char property_name[80];
-//  snprintf(property_name, 80, "propulsion/c-thrust[%u]", EngineNum);
-//  PropertyManager->Untie( property_name );
-
   Debug(1);
 }
 
 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
-double FGNozzle::Calculate(double CfPc)
+double FGNozzle::Calculate(double vacThrust)
 {
   double pAtm = fdmex->GetAtmosphere()->GetPressure();
-  Thrust = max((double)0.0, (CfPc * AreaT + (PE - pAtm)*Area2) * nzlEff);
-  vFn(1) = Thrust * cos(ReverserAngle);
+  Thrust = max((double)0.0, vacThrust - pAtm*Area);
 
-  ThrustCoeff = max((double)0.0, CfPc / ((pAtm - PE) * Area2));
+  vFn(1) = Thrust * cos(ReverserAngle);
 
   return Thrust;
 }
 
 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
-double FGNozzle::GetPowerRequired(void)
-{
-  return PE;
-}
-
-//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
 string FGNozzle::GetThrusterLabels(int id, string delimeter)
 {
   std::ostringstream buf;
 
-  buf << Name << "_Thrust[" << id << ']';
+  buf << Name << " Thrust (engine " << id << " in lbs)";
 
   return buf.str();
 }
@@ -170,10 +145,7 @@ void FGNozzle::Debug(int from)
   if (debug_lvl & 1) { // Standard console startup message output
     if (from == 0) { // Constructor
       cout << "      Nozzle Name: " << Name << endl;
-      cout << "      Nozzle Exit Pressure = " << PE << endl;
-      cout << "      Nozzle Expansion Ratio = " << ExpR << endl;
-      cout << "      Nozzle Efficiency = " << nzlEff << endl;
-      cout << "      Nozzle Diameter = " << Diameter << endl;
+      cout << "      Nozzle Exit Area = " << Area << endl;
     }
   }
   if (debug_lvl & 2 ) { // Instantiation/Destruction notification