%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/
#include "FGAerodynamics.h"
+#include "FGFactorGroup.h"
+#include "FGCoefficient.h"
static const char *IdSrc = "$Id$";
static const char *IdHdr = ID_AERODYNAMICS;
-extern short debug_lvl;
-
/*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
CLASS IMPLEMENTATION
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/
-FGAerodynamics::FGAerodynamics(FGFDMExec* FDMExec) : FGModel(FDMExec),
- vMoments(3),
- vForces(3),
- vFs(3),
- vLastFs(3),
- vDXYZcg(3)
+FGAerodynamics::FGAerodynamics(FGFDMExec* FDMExec) : FGModel(FDMExec)
{
Name = "FGAerodynamics";
Coeff = new CoeffArray[6];
- if (debug_lvl & 2) cout << "Instantiated: FGAerodynamics" << endl;
+ Debug(0);
}
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
}
delete[] Coeff;
- if (debug_lvl & 2) cout << "Destroyed: FGAerodynamics" << endl;
+ Debug(1);
}
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
bool FGAerodynamics::Run(void)
{
- float alpha, beta;
unsigned int axis_ctr,ctr;
if (!FGModel::Run()) {
- alpha = Translation->Getalpha();
- beta = Translation->Getbeta();
-
vLastFs = vFs;
vFs.InitMatrix();
vFs(axis_ctr+1) += Coeff[axis_ctr][ctr]->TotalValue();
}
}
+ //correct signs of drag and lift to wind axes convention
+ //positive forward, right, down
+ vFs(1)*=-1; vFs(3)*=-1;
+ //cout << "Aircraft::vFs: " << vFs << endl;
+ vForces = State->GetTs2b()*vFs;
- vForces = State->GetTs2b(alpha, beta)*vFs;
-
- // see http://home.earthlink.net/~apeden/jsbsim_moments_due_to_forces.txt
- // for details on this
+ vDXYZcg(eX) = -(Aircraft->GetXYZrp(eX)
+ - MassBalance->GetXYZcg(eX))*inchtoft;
+ vDXYZcg(eY) = (Aircraft->GetXYZrp(eY)
+ - MassBalance->GetXYZcg(eY))*inchtoft;
+ vDXYZcg(eZ) = -(Aircraft->GetXYZrp(eZ)
+ - MassBalance->GetXYZcg(eZ))*inchtoft;
- vDXYZcg(eX) = -(Aircraft->GetXYZrp(eX) - MassBalance->GetXYZcg(eX))/12.0;
- vDXYZcg(eY) = (Aircraft->GetXYZrp(eY) - MassBalance->GetXYZcg(eY))/12.0;
- vDXYZcg(eZ) = -(Aircraft->GetXYZrp(eZ) - MassBalance->GetXYZcg(eZ))/12.0;
-
- vMoments(eL) = vForces(eZ)*vDXYZcg(eY) - vForces(eY)*vDXYZcg(eZ);
- vMoments(eM) = vForces(eX)*vDXYZcg(eZ) - vForces(eZ)*vDXYZcg(eX);
- vMoments(eN) = vForces(eY)*vDXYZcg(eX) - vForces(eX)*vDXYZcg(eY);
+ vMoments = vDXYZcg*vForces; // M = r X F
for (axis_ctr = 0; axis_ctr < 3; axis_ctr++) {
for (ctr = 0; ctr < Coeff[axis_ctr+3].size(); ctr++) {
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-bool FGAerodynamics::LoadAerodynamics(FGConfigFile* AC_cfg)
+bool FGAerodynamics::Load(FGConfigFile* AC_cfg)
{
string token, axis;
axis = AC_cfg->GetValue("NAME");
AC_cfg->GetNextConfigLine();
while ((token = AC_cfg->GetValue()) != string("/AXIS")) {
- ca.push_back(new FGCoefficient(FDMExec, AC_cfg));
- if (debug_lvl > 0) DisplayCoeffFactors(ca.back()->Getmultipliers());
+ if( token == "COEFFICIENT" ) {
+ ca.push_back( new FGCoefficient(FDMExec) );
+ ca.back()->Load(AC_cfg);
+ cm[ca.back()->Getname()]=ca.back();
+ } else if ( token == "GROUP" ) {
+ ca.push_back( new FGFactorGroup(FDMExec) );
+ ca.back()->Load(AC_cfg);
+ cm[ca.back()->Getname()]=ca.back();
+ }
}
Coeff[AxisIdx[axis]] = ca;
AC_cfg->GetNextConfigLine();
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-void FGAerodynamics::DisplayCoeffFactors(vector <eParam> multipliers)
-{
- unsigned int i;
-
- cout << " Non-Dimensionalized by: ";
-
- for (i=0; i<multipliers.size();i++) cout << State->paramdef[multipliers[i]];
-
- cout << endl;
-}
-
-//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
string FGAerodynamics::GetCoefficientStrings(void)
{
string CoeffStrings = "";
} else {
CoeffStrings += ", ";
}
- CoeffStrings += Coeff[axis][sd]->Getname();
+ CoeffStrings += Coeff[axis][sd]->GetCoefficientStrings();
}
}
-
return CoeffStrings;
}
string FGAerodynamics::GetCoefficientValues(void)
{
string SDValues = "";
- char buffer[10];
bool firstime = true;
for (unsigned int axis = 0; axis < 6; axis++) {
} else {
SDValues += ", ";
}
- sprintf(buffer, "%9.6f", Coeff[axis][sd]->GetSD());
- SDValues += string(buffer);
+ SDValues += Coeff[axis][sd]->GetCoefficientValues();
}
}
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-float FGAerodynamics::GetNlf(void)
+double FGAerodynamics::GetLoD(void)
{
- if (fabs(Position->GetGamma()) < 1.57) {
- return (vFs(eZ)/(MassBalance->GetWeight()*cos(Position->GetGamma())));
- } else {
- return 0.0;
- }
-}
-
-//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-float FGAerodynamics::GetLoD(void)
-{
- float LoD;
-
if (vFs(1) != 0.00) return vFs(3)/vFs(1);
else return 0.00;
}
-
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-void FGAerodynamics::Debug(void)
+// The bitmasked value choices are as follows:
+// unset: In this case (the default) JSBSim would only print
+// out the normally expected messages, essentially echoing
+// the config files as they are read. If the environment
+// variable is not set, debug_lvl is set to 1 internally
+// 0: This requests JSBSim not to output any messages
+// whatsoever.
+// 1: This value explicity requests the normal JSBSim
+// startup messages
+// 2: This value asks for a message to be printed out when
+// a class is instantiated
+// 4: When this value is set, a message is displayed when a
+// FGModel object executes its Run() method
+// 8: When this value is set, various runtime state variables
+// are printed out periodically
+// 16: When set various parameters are sanity checked and
+// a message is printed out when they go out of bounds
+
+void FGAerodynamics::Debug(int from)
{
- //TODO: Add your source code here
+ if (debug_lvl <= 0) return;
+
+ if (debug_lvl & 1) { // Standard console startup message output
+ if (from == 0) { // Constructor
+
+ }
+ }
+ if (debug_lvl & 2 ) { // Instantiation/Destruction notification
+ if (from == 0) cout << "Instantiated: FGAerodynamics" << endl;
+ if (from == 1) cout << "Destroyed: FGAerodynamics" << endl;
+ }
+ if (debug_lvl & 4 ) { // Run() method entry print for FGModel-derived objects
+ }
+ if (debug_lvl & 8 ) { // Runtime state variables
+ }
+ if (debug_lvl & 16) { // Sanity checking
+ }
+ if (debug_lvl & 64) {
+ if (from == 0) { // Constructor
+ cout << IdSrc << endl;
+ cout << IdHdr << endl;
+ }
+ }
}