genfm.h File Reference

#include "aerodyn.h"

Include dependency graph for genfm.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
struct	GenericAerodynamicData
struct	GenericAerodynamicData::GenericAerodynamicCoef
class	GenericAerodynamicForce

Functions
Elem *	ReadGenericAerodynamicForce (DataManager *pDM, MBDynParser &HP, const DofOwner *pDO, unsigned int uLabel)

Function Documentation

Elem* ReadGenericAerodynamicForce	(	DataManager *	pDM,
		MBDynParser &	HP,
		const DofOwner *	pDO,
		unsigned int	uLabel
	)

Definition at line 802 of file genfm.cc.

References Elem::AERODYNAMIC, GenericAerodynamicData::bAlphaFirst, Eye3, DataManager::fReadOutput(), IncludeParser::GetFileName(), IncludeParser::GetLineData(), MBDynParser::GetPosRel(), HighParser::GetReal(), MBDynParser::GetRotRel(), HighParser::GetYesNo(), HighParser::IsKeyWord(), M_PI, MBDYN_EXCEPT_ARGS, ReadGenericAerodynamicData(), DataManager::ReadNode(), SAFENEWWITHCONSTRUCTOR, Node::STRUCTURAL, and Zero3.

{
        /* Nodo */
        const StructNode* pNode = pDM->ReadNode<const StructNode, Node::STRUCTURAL>(HP);

        /* The offset is in the reference frame of the node */
        ReferenceFrame RF(pNode);
        Vec3 f(Zero3);
        if (HP.IsKeyWord("position")) {
                f = HP.GetPosRel(RF);
        }

        /* the orientation is in flight mechanics (FIXME?) reference frame:
         * X forward, Y to the right, Z down */
        Mat3x3 Ra(Eye3);
        if (HP.IsKeyWord("orientation")) {
                Ra = HP.GetRotRel(RF);
        }

        /* 1. by default, which means that coefficients are only normalized
         * by the dynamic pressure */
        doublereal dRefSurface = 1.;
        doublereal dRefLength = 1.;
        bool bAlphaFirst(false);

        if (HP.IsKeyWord("reference" "surface")) {
                dRefSurface = HP.GetReal();
        }

        if (HP.IsKeyWord("reference" "length")) {
                dRefLength = HP.GetReal();
        }

        if (HP.IsKeyWord("alpha" "first")) {
                if (!HP.GetYesNo(bAlphaFirst)) {
                        silent_cerr("GenericAerodynamicForce(" << uLabel << "): "
                                "invalid \"alpha first\" value at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }
        }

        /* TODO: allow to reference previously loaded data */
        GenericAerodynamicData *pData = 0;
        if (HP.IsKeyWord("file")) {
                doublereal dScaleAngle = 1.;
                if (HP.IsKeyWord("angle" "units")) {
                        if (HP.IsKeyWord("radians")) {
                                dScaleAngle = 1.;

                        } else if (HP.IsKeyWord("degrees")) {
                                dScaleAngle = M_PI/180.;

                        } else {
                                silent_cerr("GenericAerodynamicForce(" << uLabel << "): "
                                        "unknown angle unit at line " << HP.GetLineData() << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }

                } else if (HP.IsKeyWord("scale" "angles")) {
                        doublereal d = HP.GetReal(dScaleAngle);
                        if (d <= 0.) {
                                silent_cerr("GenericAerodynamicForce(" << uLabel << "): "
                                        "invalid angle scale factor " << d
                                        << " at line " << HP.GetLineData() << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }
                        dScaleAngle = d;
                }

                doublereal dScaleLength = 1.;
                if (HP.IsKeyWord("scale" "lengths")) {
                        doublereal d = HP.GetReal(dScaleLength);
                        if (d <= 0.) {
                                silent_cerr("GenericAerodynamicForce(" << uLabel << "): "
                                        "invalid length scale factor " << d
                                        << " at line " << HP.GetLineData() << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }
                        dScaleLength = d;
                }

                std::string fname(HP.GetFileName());
                pData = ReadGenericAerodynamicData(fname, dScaleAngle, dScaleLength, bAlphaFirst);

        } else if (HP.IsKeyWord("reference")) {
                silent_cerr("GenericAerodynamicForce(" << uLabel << "): "
                        "references to generic aerodynamic data not implemented yet "
                        "at line " << HP.GetLineData() << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);

        } else {
                silent_cerr("GenericAerodynamicForce(" << uLabel << "): "
                        "keyword \"file\" or \"reference\" expected "
                        "at line " << HP.GetLineData() << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        flag fOut = pDM->fReadOutput(HP, Elem::AERODYNAMIC);

        Elem *pEl = 0;
        SAFENEWWITHCONSTRUCTOR(pEl, GenericAerodynamicForce,
                GenericAerodynamicForce(uLabel, pDO,
                        pNode, f, Ra,
                        dRefSurface, dRefLength, bAlphaFirst,
                        pData, fOut));

        return pEl;
}

Here is the call graph for this function: