Models a gas cell. More...

#include <FGGasCell.h>

Inheritance diagram for FGGasCell:

Collaboration diagram for FGGasCell:

Classes
struct	Inputs

Public Member Functions
	FGGasCell (FGFDMExec exec, Element el, unsigned int num, const struct Inputs &input)
	Constructor. More...

void	Calculate (double dt)
	Runs the gas cell model; called by BuoyantForces.

int	GetIndex (void) const
	Get the index of this gas cell. More...

const FGMatrix33 &	GetInertia (void) const
	Get the moments of inertia of the gas cell (including any ballonets) More...

double	GetMass (void) const
	Get the current mass of the gas cell (including any ballonets) More...

const FGColumnVector3 &	GetMassMoment (void) const
	Get the moment due to mass of the gas cell (including any ballonets) More...

double	GetPressure (void) const
	Get the current gas pressure inside the gas cell. More...

double	GetTemperature (void) const
	Get the current gas temperature inside the gas cell. More...

const FGColumnVector3 &	GetXYZ (void) const
	Get the center of gravity location of the gas cell (including any ballonets) More...

double	GetXYZ (int idx) const
	Get the center of gravity location of the gas cell (including any ballonets) More...

Public Member Functions inherited from FGForce
	FGForce (FGFDMExec *FDMExec)
	Constructor.

	FGForce (const FGForce &force)

	~FGForce ()
	Destructor.

const FGColumnVector3 &	GetActingLocation (void) const

double	GetActingLocationX (void) const

double	GetActingLocationY (void) const

double	GetActingLocationZ (void) const

const FGColumnVector3 &	GetAnglesToBody (void) const

double	GetAnglesToBody (int axis) const

virtual const FGColumnVector3 &	GetBodyForces (void)

double	GetBodyXForce (void) const

double	GetBodyYForce (void) const

double	GetBodyZForce (void) const

const FGColumnVector3 &	GetLocation (void) const

double	GetLocationX (void) const

double	GetLocationY (void) const

double	GetLocationZ (void) const

const FGColumnVector3 &	GetMoments (void) const

double	GetPitch (void) const

TransformType	GetTransformType (void) const

double	GetYaw (void) const

void	SetActingLocation (double x, double y, double z)
	Acting point of application. More...

void	SetActingLocation (const FGColumnVector3 &vv)

double	SetActingLocationX (double x)

double	SetActingLocationY (double y)

double	SetActingLocationZ (double z)

void	SetAnglesToBody (double broll, double bpitch, double byaw)

void	SetAnglesToBody (const FGColumnVector3 &vv)

void	SetLocation (double x, double y, double z)

void	SetLocation (const FGColumnVector3 &vv)

void	SetLocationX (double x)

void	SetLocationY (double y)

void	SetLocationZ (double z)

void	SetPitch (double pitch)

void	SetTransformType (TransformType ii)

void	SetYaw (double yaw)

const FGMatrix33 &	Transform (void) const

void	UpdateCustomTransformMatrix (void)

Public Member Functions inherited from FGJSBBase
	FGJSBBase ()
	Constructor for FGJSBBase.

virtual	~FGJSBBase ()
	Destructor for FGJSBBase.

void	disableHighLighting (void)
	Disables highlighting in the console output.

std::string	GetVersion (void)
	Returns the version number of JSBSim. More...

void	PutMessage (const Message &msg)
	Places a Message structure on the Message queue. More...

void	PutMessage (const std::string &text)
	Creates a message with the given text and places it on the queue. More...

void	PutMessage (const std::string &text, bool bVal)
	Creates a message with the given text and boolean value and places it on the queue. More...

void	PutMessage (const std::string &text, int iVal)
	Creates a message with the given text and integer value and places it on the queue. More...

void	PutMessage (const std::string &text, double dVal)
	Creates a message with the given text and double value and places it on the queue. More...

int	SomeMessages (void)
	Reads the message on the queue (but does not delete it). More...

void	ProcessMessage (void)
	Reads the message on the queue and removes it from the queue. More...

Message *	ProcessNextMessage (void)
	Reads the next message on the queue and removes it from the queue. More...

Public Attributes
const struct Inputs &	in

Additional Inherited Members
Public Types inherited from FGForce
enum	TransformType { tNone, tWindBody, tLocalBody, tCustom }

Public Types inherited from FGJSBBase
enum	{ eL = 1, eM, eN }
	Moments L, M, N.

enum	{ eP = 1, eQ, eR }
	Rates P, Q, R.

enum	{ eU = 1, eV, eW }
	Velocities U, V, W.

enum	{ eX = 1, eY, eZ }
	Positions X, Y, Z.

enum	{ ePhi = 1, eTht, ePsi }
	Euler angles Phi, Theta, Psi.

enum	{ eDrag = 1, eSide, eLift }
	Stability axis forces, Drag, Side force, Lift.

enum	{ eRoll = 1, ePitch, eYaw }
	Local frame orientation Roll, Pitch, Yaw.

enum	{ eNorth = 1, eEast, eDown }
	Local frame position North, East, Down.

enum	{ eLat = 1, eLong, eRad }
	Locations Radius, Latitude, Longitude.

enum	{ inNone = 0, inDegrees, inRadians, inMeters, inFeet }
	Conversion specifiers.

Static Public Member Functions inherited from FGJSBBase
static double	CelsiusToFahrenheit (double celsius)
	Converts from degrees Celsius to degrees Fahrenheit. More...

static double	CelsiusToKelvin (double celsius)
	Converts from degrees Celsius to degrees Kelvin. More...

static double	CelsiusToRankine (double celsius)
	Converts from degrees Celsius to degrees Rankine. More...

static double	Constrain (double min, double value, double max)
	Constrain a value between a minimum and a maximum value.

static bool	EqualToRoundoff (double a, double b)
	Finite precision comparison. More...

static bool	EqualToRoundoff (float a, float b)
	Finite precision comparison. More...

static bool	EqualToRoundoff (float a, double b)
	Finite precision comparison. More...

static bool	EqualToRoundoff (double a, float b)
	Finite precision comparison. More...

static double	FahrenheitToCelsius (double fahrenheit)
	Converts from degrees Fahrenheit to degrees Celsius. More...

static double	FeetToMeters (double measure)
	Converts from feet to meters. More...

static double	GaussianRandomNumber (void)

static double	KelvinToCelsius (double kelvin)
	Converts from degrees Kelvin to degrees Celsius. More...

static double	KelvinToFahrenheit (double kelvin)
	Converts from degrees Kelvin to degrees Fahrenheit. More...

static double	KelvinToRankine (double kelvin)
	Converts from degrees Kelvin to degrees Rankine. More...

static double	MachFromVcalibrated (double vcas, double p, double psl, double rhosl)
	Calculate the Mach number from the calibrated airspeed. More...

static double	PitotTotalPressure (double mach, double p)
	Compute the total pressure in front of the Pitot tube. More...

static double	RankineToCelsius (double rankine)
	Converts from degrees Rankine to degrees Celsius. More...

static double	RankineToKelvin (double rankine)
	Converts from degrees Rankine to degrees Kelvin. More...

static double	sign (double num)

static double	VcalibratedFromMach (double mach, double p, double psl, double rhosl)
	Calculate the calibrated airspeed from the Mach number. More...

Static Public Attributes inherited from FGJSBBase
static short	debug_lvl = 1

static char	highint [5] = {27, '[', '1', 'm', '\0' }
	highlights text

static char	halfint [5] = {27, '[', '2', 'm', '\0' }
	low intensity text

static char	normint [6] = {27, '[', '2', '2', 'm', '\0' }
	normal intensity text

static char	reset [5] = {27, '[', '0', 'm', '\0' }
	resets text properties

static char	underon [5] = {27, '[', '4', 'm', '\0' }
	underlines text

static char	underoff [6] = {27, '[', '2', '4', 'm', '\0' }
	underline off

static char	fgblue [6] = {27, '[', '3', '4', 'm', '\0' }
	blue text

static char	fgcyan [6] = {27, '[', '3', '6', 'm', '\0' }
	cyan text

static char	fgred [6] = {27, '[', '3', '1', 'm', '\0' }
	red text

static char	fggreen [6] = {27, '[', '3', '2', 'm', '\0' }
	green text

static char	fgdef [6] = {27, '[', '3', '9', 'm', '\0' }
	default text

Protected Member Functions inherited from FGJSBBase
void	Debug (int)

Static Protected Member Functions inherited from FGJSBBase
static std::string	CreateIndexedPropertyName (const std::string &Property, int index)

Protected Attributes inherited from FGForce
FGFDMExec *	fdmex

FGMatrix33	mT

TransformType	ttype

FGColumnVector3	vActingXYZn

FGColumnVector3	vFn

FGColumnVector3	vH

FGColumnVector3	vMn

FGColumnVector3	vOrient

FGColumnVector3	vXYZn

Static Protected Attributes inherited from FGJSBBase
static const double	degtorad = 0.017453292519943295769236907684886

static const double	fpstokts = 1.0/ktstofps

static const double	fttom = 0.3048

static int	gaussian_random_number_phase = 0

static const double	hptoftlbssec = 550.0

static const double	in3tom3 = 1.638706E-5

static const double	inchtoft = 0.08333333

static const double	inhgtopa = 3386.38

static const std::string	JSBSim_version = "1.0 " __DATE__ " " __TIME__

static const double	kgtolb = 2.20462

static const double	kgtoslug = 0.06852168

static const double	ktstofps = 1.68781

static const double	lbtoslug = 1.0/slugtolb

static Message	localMsg

static const double	m3toft3 = 1.0/(fttomfttomfttom)

static double	Mair = 28.9645

static unsigned int	messageId = 0

static std::queue< Message >	Messages

static const std::string	needed_cfg_version = "2.0"

static const double	psftoinhg = 0.014138

static const double	psftopa = 47.88

static const double	radtodeg = 57.295779513082320876798154814105

static double	Reng = 1716.56

static double	Rstar = 1545.348

static const double	SHRatio = 1.40

static const double	slugtolb = 32.174049

Detailed Description

Models a gas cell.

Author: Anders Gidenstam

Configuration File Format:

<buoyant_forces>
  <gas_cell type="{HYDROGEN | HELIUM | AIR}">
    <location unit="{M | IN}">
      <x> {number} </x>
      <y> {number} </y>
      <z> {number} </z>
    </location>
    <x_width unit="{M | IN}"> {number} </x_width>
    <y_radius unit="{M | IN}"> {number} </y_radius>
    <z_radius unit="{M | IN}"> {number} </z_radius>
    <max_overpressure unit="{PA | PSI}"> {number} </max_overpressure>
    <valve_coefficient unit="{M4*SEC/KG | FT4*SEC/SLUG}"> {number} </valve_coefficient>
    <fullness> {number} </fullness>  
    <heat>
      {heat transfer coefficients} [lbs ft / sec]
    </heat>
    <ballonet>
      <location unit="{M | IN}">
        <x> {number} </x>
        <y> {number} </y>
        <z> {number} </z>
      </location>
      <x_width unit="{M | IN}"> {number} </x_width>
      <y_radius unit="{M | IN}"> {number} </y_radius>
      <z_radius unit="{M | IN}"> {number} </z_radius>
      <max_overpressure unit="{PA | PSI}"> {number} </max_overpressure>
      <valve_coefficient unit="{M4*SEC/KG | FT4*SEC/SLUG}"> {number} </valve_coefficient>
      <fullness> {number} </fullness>  
      <heat>
       {heat transfer coefficients} [lb ft / (sec Rankine)]
      </heat>
      <blower_input>
       {input air flow function} [ft^3 / sec]
      </blower_input>
    </ballonet>
  </gas_cell>
</buoyant_forces>

Definition of the gas cell configuration file parameters:

type - One of HYDROGEN, HELIUM or AIR.
location - Location of cell center in the aircraft's structural frame. Currently this is were the forces of the cell is applied.
{x|y|z}_radius - Radius along in the respective direction (both ends).
{x|y|z}_width - Width in the respective direction. NOTE: A 'x', 'y', 'z'-radius/width combination must be specified.
fullness - Initial fullness of the cell, normally [0,1], values >1 initialize the cell at pressure.
max_overpressure - Maximum cell overpressure (excess is automatically valved off).
valve_coefficient - Capacity of the manual valve. The valve is considered to be located at the top of the cell. The valve coefficient determine the flow out of the cell according to: dVolume/dt = ValveCoefficient * DeltaPressure.
heat - Zero or more FGFunction:s describing the heat flow from the atmosphere into the gas cell. Unit: [lb ft / (sec Rankine)]. If there are no heat transfer functions at all the gas cell temperature will equal that of the surrounding atmosphere. A constant function returning 0 results in adiabatic behaviour.
ballonet - Zero or more ballonets, i.e. air bags inside the gas cell. Ballonets are used to maintain the volume of the gas cell and keep its internal pressure higher than that of the surrounding environment.
- location - Location of ballonet center in the aircraft's structural frame.
- {x|y|z}_radius - Radius along in the respective direction (both ends).
- {x|y|z}_width - Width in the respective direction.
- max_overpressure - Maximum ballonet overpressure (excess is automatically valved off).
- valve_coefficient - Capacity of the exit valve between the ballonet and the atmosphere. The valve coefficient determine the flow out of the cell according to: dVolume/dt = ValveCoefficient * DeltaPressure.
- heat - Zero or more FGFunction:s describing the heat flow from the enclosing gas cell into the ballonet. Unit: [lb ft / (sec Rankine)]
- blower_input - One FGFunction describing the air flow into the ballonet. Unit: ft³ / sec

Definition at line 172 of file FGGasCell.h.

Constructor & Destructor Documentation

◆ FGGasCell()

FGGasCell	(	FGFDMExec *	exec,
		Element *	el,
		unsigned int	num,
		const struct Inputs &	input
	)

Constructor.

Parameters

exec	Executive a pointer to the parent executive object
el	Pointer to configuration file XML node
num	Gas cell index number.

Definition at line 65 of file FGGasCell.cpp.

   : FGForce(exec), in(input)
 {
   string token;
   Element* element;
 
   FGPropertyManager* PropertyManager = exec->GetPropertyManager();
   MassBalance = exec->GetMassBalance();
 
   gasCellJ = FGMatrix33();
   gasCellM = FGColumnVector3();
 
   Buoyancy = MaxVolume = MaxOverpressure = Temperature = Pressure =
     Contents = Volume = dVolumeIdeal = 0.0;
   Xradius = Yradius = Zradius = Xwidth = Ywidth = Zwidth = 0.0;
   ValveCoefficient = ValveOpen = 0.0;
   CellNum = num;
 
   // NOTE: In the local system X points north, Y points east and Z points down.
   SetTransformType(FGForce::tLocalBody);
 
   type = el->GetAttributeValue("type");
   if      (type == "HYDROGEN") Type = ttHYDROGEN;
   else if (type == "HELIUM")   Type = ttHELIUM;
   else if (type == "AIR")      Type = ttAIR;
   else                         Type = ttUNKNOWN;
 
   element = el->FindElement("location");
   if (element) {
     vXYZ = element->FindElementTripletConvertTo("IN");
   } else {
     cerr << "Fatal Error: No location found for this gas cell." << endl;
     exit(-1);
   }
   if ((el->FindElement("x_radius") || el->FindElement("x_width")) &&
       (el->FindElement("y_radius") || el->FindElement("y_width")) &&
       (el->FindElement("z_radius") || el->FindElement("z_width"))) {
 
     if (el->FindElement("x_radius")) {
       Xradius = el->FindElementValueAsNumberConvertTo("x_radius", "FT");
     }
     if (el->FindElement("y_radius")) {
       Yradius = el->FindElementValueAsNumberConvertTo("y_radius", "FT");
     }
     if (el->FindElement("z_radius")) {
       Zradius = el->FindElementValueAsNumberConvertTo("z_radius", "FT");
     }
 
     if (el->FindElement("x_width")) {
       Xwidth = el->FindElementValueAsNumberConvertTo("x_width", "FT");
     }
     if (el->FindElement("y_width")) {
       Ywidth = el->FindElementValueAsNumberConvertTo("y_width", "FT");
     }
     if (el->FindElement("z_width")) {
       Zwidth = el->FindElementValueAsNumberConvertTo("z_width", "FT");
     }
 
     // The volume is a (potentially) extruded ellipsoid.
     // However, currently only a few combinations of radius and width are
     // fully supported.
     if ((Xradius != 0.0) && (Yradius != 0.0) && (Zradius != 0.0) &&
         (Xwidth  == 0.0) && (Ywidth  == 0.0) && (Zwidth  == 0.0)) {
       // Ellipsoid volume.
       MaxVolume = 4.0  * M_PI * Xradius * Yradius * Zradius / 3.0;
     } else if  ((Xradius == 0.0) && (Yradius != 0.0) && (Zradius != 0.0) &&
                 (Xwidth  != 0.0) && (Ywidth  == 0.0) && (Zwidth  == 0.0)) {
       // Cylindrical volume.
       MaxVolume = M_PI * Yradius * Zradius * Xwidth;
     } else {
       cerr << "Warning: Unsupported gas cell shape." << endl;
       MaxVolume = 
         (4.0  * M_PI * Xradius * Yradius * Zradius / 3.0 +
          M_PI * Yradius * Zradius * Xwidth +
          M_PI * Xradius * Zradius * Ywidth +
          M_PI * Xradius * Yradius * Zwidth +
          2.0  * Xradius * Ywidth * Zwidth +
          2.0  * Yradius * Xwidth * Zwidth +
          2.0  * Zradius * Xwidth * Ywidth +
          Xwidth * Ywidth * Zwidth);
     }
   } else {
     cerr << "Fatal Error: Gas cell shape must be given." << endl;
     exit(-1);
   }
   if (el->FindElement("max_overpressure")) {
     MaxOverpressure = el->FindElementValueAsNumberConvertTo("max_overpressure",
                                                             "LBS/FT2");
   }
   if (el->FindElement("fullness")) {
     const double Fullness = el->FindElementValueAsNumber("fullness");
     if (0 <= Fullness) { 
       Volume = Fullness * MaxVolume; 
     } else {
       cerr << "Warning: Invalid initial gas cell fullness value." << endl;
     }
   }  
   if (el->FindElement("valve_coefficient")) {
     ValveCoefficient =
       el->FindElementValueAsNumberConvertTo("valve_coefficient",
                                             "FT4*SEC/SLUG");
     ValveCoefficient = max(ValveCoefficient, 0.0);
   }
 
   // Initialize state
   SetLocation(vXYZ);
 
   if (Temperature == 0.0) {
     Temperature = in.Temperature;
   }
   if (Pressure == 0.0) {
     Pressure = in.Pressure;
   }
   if (Volume != 0.0) {
     // Calculate initial gas content.
     Contents = Pressure * Volume / (R * Temperature);
     
     // Clip to max allowed value.
     const double IdealPressure = Contents * R * Temperature / MaxVolume;
     if (IdealPressure > Pressure + MaxOverpressure) {
       Contents = (Pressure + MaxOverpressure) * MaxVolume / (R * Temperature);
       Pressure = Pressure + MaxOverpressure;
     } else {
       Pressure = max(IdealPressure, Pressure);
     }
   } else {
     // Calculate initial gas content.
     Contents = Pressure * MaxVolume / (R * Temperature);
   }
 
   Volume = Contents * R * Temperature / Pressure;
   Mass = Contents * M_gas();
 
   // Bind relevant properties
   string property_name, base_property_name;
 
   base_property_name = CreateIndexedPropertyName("buoyant_forces/gas-cell", CellNum);
 
   property_name = base_property_name + "/max_volume-ft3";
   PropertyManager->Tie( property_name.c_str(), &MaxVolume, false );
   PropertyManager->GetNode()->SetWritable( property_name, false );
   property_name = base_property_name + "/temp-R";
   PropertyManager->Tie( property_name.c_str(), &Temperature, false );
   property_name = base_property_name + "/pressure-psf";
   PropertyManager->Tie( property_name.c_str(), &Pressure, false );
   property_name = base_property_name + "/volume-ft3";
   PropertyManager->Tie( property_name.c_str(), &Volume, false );
   property_name = base_property_name + "/buoyancy-lbs";
   PropertyManager->Tie( property_name.c_str(), &Buoyancy, false );
   property_name = base_property_name + "/contents-mol";
   PropertyManager->Tie( property_name.c_str(), &Contents, false );
   property_name = base_property_name + "/valve_open";
   PropertyManager->Tie( property_name.c_str(), &ValveOpen, false );
 
   Debug(0);
 
   // Read heat transfer coefficients
   if (Element* heat = el->FindElement("heat")) {
     Element* function_element = heat->FindElement("function");
     while (function_element) {
       HeatTransferCoeff.push_back(new FGFunction(PropertyManager,
                                                  function_element));
       function_element = heat->FindNextElement("function");
     }
   }
 
   // Load ballonets if there are any
   if (Element* ballonet_element = el->FindElement("ballonet")) {
     while (ballonet_element) {
       Ballonet.push_back(new FGBallonet(exec,
                                         ballonet_element,
                                         Ballonet.size(),
                                         this, in));
       ballonet_element = el->FindNextElement("ballonet");
     }
   }
 
 }

Here is the call graph for this function:

Member Function Documentation

◆ GetIndex()

int GetIndex ( void ) const

inline

Get the index of this gas cell.

Returns: gas cell index.

Definition at line 196 of file FGGasCell.h.

196 {return CellNum;}

◆ GetInertia()

const FGMatrix33& GetInertia ( void ) const

inline

Get the moments of inertia of the gas cell (including any ballonets)

Returns: moments of inertia matrix in the body frame in slug ft².

Definition at line 215 of file FGGasCell.h.

215 {return gasCellJ;}

◆ GetMass()

double GetMass ( void ) const

inline

Get the current mass of the gas cell (including any ballonets)

Returns: gas mass in slug.

Definition at line 210 of file FGGasCell.h.

210 {return Mass;}

Here is the caller graph for this function:

◆ GetMassMoment()

const FGColumnVector3& GetMassMoment ( void ) const

inline

Get the moment due to mass of the gas cell (including any ballonets)

Note that the buoyancy of the gas cell is handled separately by the FGForce part and not included here.

Returns: moment vector in the structural frame in lbs in.

Definition at line 222 of file FGGasCell.h.

222 {return gasCellM;}

◆ GetPressure()

double GetPressure ( void ) const

inline

Get the current gas pressure inside the gas cell.

Returns: gas pressure in lbs / ft².

Definition at line 230 of file FGGasCell.h.

230 {return Pressure;}

◆ GetTemperature()

double GetTemperature ( void ) const

inline

Get the current gas temperature inside the gas cell.

Returns: gas temperature in Rankine.

Definition at line 226 of file FGGasCell.h.

226 {return Temperature;}

◆ GetXYZ() [1/2]

const FGColumnVector3& GetXYZ ( void ) const

inline

Get the center of gravity location of the gas cell (including any ballonets)

Returns: CoG location in the structural frame in inches.

Definition at line 201 of file FGGasCell.h.

201 {return vXYZ;}

Here is the caller graph for this function:

◆ GetXYZ() [2/2]

double GetXYZ ( int idx ) const

inline

Get the center of gravity location of the gas cell (including any ballonets)

Returns: CoG location in the structural frame in inches.

Definition at line 206 of file FGGasCell.h.

206 {return vXYZ(idx);}

The documentation for this class was generated from the following files:

Classes

Public Member Functions

Public Attributes

Additional Inherited Members

Detailed Description

Configuration File Format:

Constructor & Destructor Documentation

◆ FGGasCell()

Member Function Documentation

◆ GetIndex()

◆ GetInertia()

◆ GetMass()

◆ GetMassMoment()

◆ GetPressure()

◆ GetTemperature()

◆ GetXYZ() [1/2]

◆ GetXYZ() [2/2]