src/Systems/vacuum.cxx

   1 // vacuum.cxx - a vacuum pump connected to the aircraft engine.
   2 // Written by David Megginson, started 2002.
   3 //
   4 // This file is in the Public Domain and comes with no warranty.
   5
   6 #ifdef HAVE_CONFIG_H
   7 #  include <config.h>
   8 #endif
   9
  10 #include "vacuum.hxx"
  11
  12 #include <cstring>
  13
  14 #include <Main/fg_props.hxx>
  15
  16
  17 VacuumSystem::VacuumSystem ( SGPropertyNode *node )
  18     :
  19     _name(node->getStringValue("name", "vacuum")),
  20     _num(node->getIntValue("number", 0)),
  21     _scale(node->getDoubleValue("scale", 1.0))
  22 {
  23     for ( int i = 0; i < node->nChildren(); ++i ) {
  24         SGPropertyNode *child = node->getChild(i);
  25         if (!strcmp(child->getName(), "rpm"))
  26             _rpms.push_back(child->getStringValue());
  27     }
  28 }
  29
  30 VacuumSystem::~VacuumSystem ()
  31 {
  32 }
  33
  34 void
  35 VacuumSystem::init()
  36 {
  37     unsigned int i;
  38     std::string branch;
  39     branch = "/systems/" + _name;
  40
  41     SGPropertyNode *node = fgGetNode(branch.c_str(), _num, true );
  42     _serviceable_node = node->getChild("serviceable", 0, true);
  43     for ( i = 0; i < _rpms.size(); i++ ) {
  44       SGPropertyNode_ptr _rpm_node = fgGetNode(_rpms[i].c_str(), true);
  45       _rpm_nodes.push_back( _rpm_node );
  46     }
  47     _pressure_node = fgGetNode("/environment/pressure-inhg", true);
  48     _suction_node = node->getChild("suction-inhg", 0, true);
  49 }
  50
  51 void
  52 VacuumSystem::bind ()
  53 {
  54 }
  55
  56 void
  57 VacuumSystem::unbind ()
  58 {
  59 }
  60
  61 void
  62 VacuumSystem::update (double dt)
  63 {
  64     // Model taken from steam.cxx
  65
  66     double suction;
  67     unsigned int i;
  68
  69     if (!_serviceable_node->getBoolValue()) {
  70         suction = 0.0;
  71     } else {
  72         // select the source with the max rpm
  73         double rpm = 0.0;
  74         for ( i = 0; i < _rpm_nodes.size(); i++ ) {
  75           double tmp = _rpm_nodes[i]->getDoubleValue() * _scale;
  76           if ( tmp > rpm ) {
  77             rpm = tmp;
  78           }
  79         }
  80         double pressure = _pressure_node->getDoubleValue();
  81         // This magic formula yields about 4 inhg at 700 rpm
  82         suction = pressure * rpm / (rpm + 4875.0);
  83
  84         // simple regulator model that clamps smoothly to about 5 inhg
  85         // over a normal rpm range
  86         double max = (rpm > 0 ? 5.39 - 1.0 / ( rpm * 0.00111 ) : 0);
  87         if ( suction < 0.0 ) suction = 0.0;
  88         if ( suction > max ) suction = max;
  89     }
  90     _suction_node->setDoubleValue(suction);
  91 }
  92
  93 // end of vacuum.cxx