aeroelem.cc File Reference

#include "mbconfig.h"
#include "dataman.h"
#include "aerodata.h"
#include "aeroelem.h"
#include "shapefnc.h"
#include "aerodc81.h"
#include "c81data.h"
#include "Rot.hh"
#include <sstream>

Include dependency graph for aeroelem.cc:

Go to the source code of this file.

Functions
static bool	ReadInducedVelocity (DataManager *pDM, MBDynParser &HP, unsigned uLabel, const char *sElemType, InducedVelocity *&pIndVel, bool &bPassive)
void	ReadAerodynamicCustomOutput (DataManager *pDM, MBDynParser &HP, unsigned int uLabel, unsigned &uFlags, OrientationDescription &od)
void	ReadOptionalAerodynamicCustomOutput (DataManager *pDM, MBDynParser &HP, unsigned int uLabel, unsigned &uFlags, OrientationDescription &od)
Elem *	ReadAerodynamicBody (DataManager *pDM, MBDynParser &HP, const DofOwner *pDO, unsigned int uLabel)
Elem *	ReadAerodynamicBeam (DataManager *pDM, MBDynParser &HP, const DofOwner *pDO, unsigned int uLabel)
Elem *	ReadAerodynamicBeam2 (DataManager *pDM, MBDynParser &HP, const DofOwner *pDO, unsigned int uLabel)

Variables
static const bool	bDefaultUseJacobian = false
static const char *	elemnames []
static const doublereal	d13 = 1./sqrt(3.)
static const doublereal	pdsi3 [] = { -1., -d13, d13 }
static const doublereal	pdsf3 [] = { -d13, d13, 1. }
static const doublereal	pdsi2 [] = { -1., 0. }
static const doublereal	pdsf2 [] = { 0., 1. }

Function Documentation

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

Definition at line 2268 of file aeroelem.cc.

References AerodynamicOutput::AEROD_OUT_NODE, AerodynamicOutput::AEROD_OUT_PGAUSS, AerodynamicOutput::AEROD_OUT_STD, Elem::AERODYNAMIC, ASSERT, bDefaultUseJacobian, Elem::BEAM, DEBUGCOUTFNAME, DEBUGLCOUT, Shape::dGet(), DataManager::fReadOutput(), HighParser::GetInt(), WithLabel::GetLabel(), IncludeParser::GetLineData(), DataManager::GetLogFile(), MBDynParser::GetPosRel(), MBDynParser::GetRotRel(), HighParser::GetYesNoOrBool(), AeroData::iGetNumDof(), HighParser::IsArg(), HighParser::IsKeyWord(), MBDYN_EXCEPT_ARGS, MYDEBUG_INPUT, AerodynamicOutput::OUTPUT_NONE, DataManager::pFindElem(), Beam::pGetNode(), ReadAeroData(), ReadInducedVelocity(), ReadOptionalAerodynamicCustomOutput(), SAFENEWWITHCONSTRUCTOR, UNKNOWN_ORIENTATION_DESCRIPTION, and Write().

Referenced by DataManager::ReadOneElem().

{
        DEBUGCOUTFNAME("ReadAerodynamicBeam");

        /* Trave */
        unsigned int uBeam = (unsigned int)HP.GetInt();

        DEBUGLCOUT(MYDEBUG_INPUT, "Linked to beam: " << uBeam << std::endl);

        /* verifica di esistenza della trave */
        Elem* p = pDM->pFindElem(Elem::BEAM, uBeam);
        if (p == 0) {
                silent_cerr("Beam3(" << uBeam << ") not defined "
                                "at line " << HP.GetLineData()
                                << std::endl);
                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
        Beam *pBeam = dynamic_cast<Beam *>(p);
        if (pBeam == 0) {
                silent_cerr("Beam(" << uBeam << ") is not a Beam3 "
                                "at line " << HP.GetLineData()
                                << std::endl);
                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
        ASSERT(pBeam != 0);

        /* Eventuale rotore */
        InducedVelocity* pIndVel = 0;
        bool bPassive(false);
        (void)ReadInducedVelocity(pDM, HP, uLabel, "AerodynamicBeam3",
                pIndVel, bPassive);

        /* Nodo 1: */

        /* Offset del corpo aerodinamico rispetto al nodo */
        const StructNode* pNode1 = pBeam->pGetNode(1);

        ReferenceFrame RF(pNode1);
        Vec3 f1(HP.GetPosRel(RF));
        DEBUGLCOUT(MYDEBUG_INPUT, "Node 1 offset: " << f1 << std::endl);

        Mat3x3 Ra1(HP.GetRotRel(RF));
        DEBUGLCOUT(MYDEBUG_INPUT,
                   "Node 1 rotation matrix: " << std::endl << Ra1 << std::endl);

        /* Nodo 2: */

        /* Offset del corpo aerodinamico rispetto al nodo */
        const StructNode* pNode2 = pBeam->pGetNode(2);

        RF = ReferenceFrame(pNode2);
        Vec3 f2(HP.GetPosRel(RF));
        DEBUGLCOUT(MYDEBUG_INPUT, "Node 2 offset: " << f2 << std::endl);

        Mat3x3 Ra2(HP.GetRotRel(RF));
        DEBUGLCOUT(MYDEBUG_INPUT,
                   "Node 2 rotation matrix: " << std::endl << Ra2 << std::endl);

        /* Nodo 3: */

        /* Offset del corpo aerodinamico rispetto al nodo */
        const StructNode* pNode3 = pBeam->pGetNode(3);

        RF = ReferenceFrame(pNode3);
        Vec3 f3(HP.GetPosRel(RF));
        DEBUGLCOUT(MYDEBUG_INPUT, "Node 3 offset: " << f3 << std::endl);

        Mat3x3 Ra3(HP.GetRotRel(RF));
        DEBUGLCOUT(MYDEBUG_INPUT,
                "Node 3 rotation matrix: " << std::endl << Ra3 << std::endl);

        Shape* pChord = 0;
        Shape* pForce = 0;
        Shape* pVelocity = 0;
        Shape* pTwist = 0;
        Shape* pTipLoss = 0;

        integer iNumber = 0;
        DriveCaller* pDC = 0;
        AeroData* aerodata = 0;

        ReadAeroData(pDM, HP, 3,
                &pChord, &pForce, &pVelocity, &pTwist, &pTipLoss,
                &iNumber, &pDC, &aerodata);

        bool bUseJacobian(false);
        if (HP.IsKeyWord("jacobian")) {
                bUseJacobian = HP.GetYesNoOrBool(bDefaultUseJacobian);
        }

        if (aerodata->iGetNumDof() > 0 && !bUseJacobian) {
                silent_cerr("AerodynamicBeam3(" << uLabel << "): "
                        "aerodynamic model needs \"jacobian, yes\" at line " << HP.GetLineData()
                        << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        OrientationDescription od = UNKNOWN_ORIENTATION_DESCRIPTION;
        unsigned uFlags = AerodynamicOutput::OUTPUT_NONE;
        ReadOptionalAerodynamicCustomOutput(pDM, HP, uLabel, uFlags, od);

        flag fOut = pDM->fReadOutput(HP, Elem::AERODYNAMIC);
        if (HP.IsArg()) {
                if (HP.IsKeyWord("std")) {
                        fOut |= AerodynamicOutput::AEROD_OUT_STD;
                } else if (HP.IsKeyWord("gauss")) {
                        fOut |= AerodynamicOutput::AEROD_OUT_PGAUSS;
                } else if (HP.IsKeyWord("node")) {
                        fOut |= AerodynamicOutput::AEROD_OUT_NODE;
                } else {
                        silent_cerr("AerodynamicBeam3(" << uLabel << "): "
                                "unknown output mode at line " << HP.GetLineData()
                                << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

        } else if (fOut) {
                fOut |= AerodynamicOutput::AEROD_OUT_STD;
        }
        fOut |= uFlags;

        Elem* pEl = 0;

        SAFENEWWITHCONSTRUCTOR(pEl,
                AerodynamicBeam,
                AerodynamicBeam(uLabel, pDO, pBeam, pIndVel, bPassive,
                        f1, f2, f3, Ra1, Ra2, Ra3,
                        pChord, pForce, pVelocity, pTwist, pTipLoss,
                        iNumber, aerodata, pDC, bUseJacobian, od, fOut));

        /* Se non c'e' il punto e virgola finale */
        if (HP.IsArg()) {
                silent_cerr("semicolon expected at line "
                        << HP.GetLineData() << std::endl);
                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        std::ostream& out = pDM->GetLogFile();
        out << "aero3: " << uLabel;

        Vec3 ra1 = Ra1.GetVec(1);
        Vec3 ra3 = Ra1.GetVec(3);
        doublereal dC = pChord->dGet(-1.);
        doublereal dP = pForce->dGet(-1.);
        out
                << " " << pNode1->GetLabel()
                << " ", (f1 + ra1*(dP - dC*3./4.)).Write(out, " ")
                << " ", (f1 + ra1*(dP + dC/4.)).Write(out, " ");

        ra1 = Ra2.GetVec(1);
        ra3 = Ra2.GetVec(3);
        dC = pChord->dGet(0.);
        dP = pForce->dGet(0.);
        out
                << " " << pNode2->GetLabel()
                << " ", (f2 + ra1*(dP - dC*3./4.)).Write(out, " ")
                << " ", (f2 + ra1*(dP + dC/4.)).Write(out, " ");

        ra1 = Ra3.GetVec(1);
        ra3 = Ra3.GetVec(3);
        dC = pChord->dGet(1.);
        dP = pForce->dGet(1.);
        out
                << " " << pNode3->GetLabel()
                << " ", (f3 + ra1*(dP - dC*3./4.)).Write(out, " ")
                << " ", (f3 + ra1*(dP + dC/4.)).Write(out, " ")
                << std::endl;

        return pEl;
} /* End of ReadAerodynamicBeam() */

Here is the call graph for this function:

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

Definition at line 3167 of file aeroelem.cc.

References AerodynamicOutput::AEROD_OUT_NODE, AerodynamicOutput::AEROD_OUT_PGAUSS, AerodynamicOutput::AEROD_OUT_STD, Elem::AERODYNAMIC, bDefaultUseJacobian, Elem::BEAM, DEBUGCOUTFNAME, DEBUGLCOUT, Shape::dGet(), DataManager::fReadOutput(), HighParser::GetInt(), WithLabel::GetLabel(), IncludeParser::GetLineData(), DataManager::GetLogFile(), MBDynParser::GetPosRel(), MBDynParser::GetRotRel(), HighParser::GetYesNoOrBool(), AeroData::iGetNumDof(), HighParser::IsArg(), HighParser::IsKeyWord(), MBDYN_EXCEPT_ARGS, MYDEBUG_INPUT, AerodynamicOutput::OUTPUT_NONE, DataManager::pFindElem(), Beam2::pGetNode(), ReadAeroData(), ReadInducedVelocity(), ReadOptionalAerodynamicCustomOutput(), SAFENEWWITHCONSTRUCTOR, UNKNOWN_ORIENTATION_DESCRIPTION, and Write().

Referenced by DataManager::ReadOneElem().

{
        DEBUGCOUTFNAME("ReadAerodynamicBeam2");

        /* Trave */
        unsigned int uBeam = (unsigned int)HP.GetInt();

        DEBUGLCOUT(MYDEBUG_INPUT, "Linked to beam: " << uBeam << std::endl);

        /* verifica di esistenza della trave */
        Elem *p = pDM->pFindElem(Elem::BEAM, uBeam);
        if (p == 0) {
                silent_cerr("Beam2(" << uBeam << ") not defined "
                                "at line " << HP.GetLineData()
                                << std::endl);
                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
        Beam2* pBeam = dynamic_cast<Beam2 *>(p);
        if (pBeam == 0) {
                silent_cerr("Beam(" << uBeam << ") is not a Beam2 "
                                "at line " << HP.GetLineData()
                                << std::endl);
                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        /* Eventuale rotore */
        InducedVelocity* pIndVel = 0;
        bool bPassive(false);
        (void)ReadInducedVelocity(pDM, HP, uLabel, "AerodynamicBeam2",
                pIndVel, bPassive);

        /* Nodo 1: */

        /* Offset del corpo aerodinamico rispetto al nodo */
        const StructNode* pNode1 = pBeam->pGetNode(1);

        ReferenceFrame RF(pNode1);
        Vec3 f1(HP.GetPosRel(RF));
        DEBUGLCOUT(MYDEBUG_INPUT, "Node 1 offset: " << f1 << std::endl);

        Mat3x3 Ra1(HP.GetRotRel(RF));
        DEBUGLCOUT(MYDEBUG_INPUT,
                   "Node 1 rotation matrix: " << std::endl << Ra1 << std::endl);

        /* Nodo 2: */

        /* Offset del corpo aerodinamico rispetto al nodo */
        const StructNode* pNode2 = pBeam->pGetNode(2);

        RF = ReferenceFrame(pNode2);
        Vec3 f2(HP.GetPosRel(RF));
        DEBUGLCOUT(MYDEBUG_INPUT, "Node 2 offset: " << f2 << std::endl);

        Mat3x3 Ra2(HP.GetRotRel(RF));
        DEBUGLCOUT(MYDEBUG_INPUT,
                "Node 2 rotation matrix: " << std::endl << Ra2 << std::endl);

        Shape* pChord = 0;
        Shape* pForce = 0;
        Shape* pVelocity = 0;
        Shape* pTwist = 0;
        Shape* pTipLoss = 0;

        integer iNumber = 0;
        DriveCaller* pDC = 0;
        AeroData* aerodata = 0;

        ReadAeroData(pDM, HP, 2,
                &pChord, &pForce, &pVelocity, &pTwist, &pTipLoss,
                &iNumber, &pDC, &aerodata);

        bool bUseJacobian(false);
        if (HP.IsKeyWord("jacobian")) {
                bUseJacobian = HP.GetYesNoOrBool(bDefaultUseJacobian);
        }

        if (aerodata->iGetNumDof() > 0 && !bUseJacobian) {
                silent_cerr("AerodynamicBeam2(" << uLabel << "): "
                        "aerodynamic model needs \"jacobian, yes\" at line " << HP.GetLineData()
                        << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        OrientationDescription od = UNKNOWN_ORIENTATION_DESCRIPTION;
        unsigned uFlags = AerodynamicOutput::OUTPUT_NONE;
        ReadOptionalAerodynamicCustomOutput(pDM, HP, uLabel, uFlags, od);

        flag fOut = pDM->fReadOutput(HP, Elem::AERODYNAMIC);
        if (HP.IsArg()) {
                if (HP.IsKeyWord("std")) {
                        fOut |= AerodynamicOutput::AEROD_OUT_STD;
                } else if (HP.IsKeyWord("gauss")) {
                        fOut |= AerodynamicOutput::AEROD_OUT_PGAUSS;
                } else if (HP.IsKeyWord("node")) {
                        fOut |= AerodynamicOutput::AEROD_OUT_NODE;
                } else {
                        silent_cerr("AerodynamicBeam2(" << uLabel << "): "
                                "unknown output mode at line " << HP.GetLineData()
                                << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

        } else if (fOut) {
                fOut |= AerodynamicOutput::AEROD_OUT_STD;
        }
        fOut |= uFlags;

        Elem* pEl = 0;

        SAFENEWWITHCONSTRUCTOR(pEl,
                AerodynamicBeam2,
                AerodynamicBeam2(uLabel, pDO, pBeam, pIndVel, bPassive,
                        f1, f2, Ra1, Ra2,
                        pChord, pForce, pVelocity, pTwist, pTipLoss,
                        iNumber, aerodata, pDC, bUseJacobian, od, fOut));

        /* Se non c'e' il punto e virgola finale */
        if (HP.IsArg()) {
                silent_cerr("semicolon expected at line "
                        << HP.GetLineData() << std::endl);
                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        std::ostream& out = pDM->GetLogFile();
        out << "aero2: " << uLabel;

        Vec3 ra1 = Ra1.GetVec(1);
        Vec3 ra3 = Ra1.GetVec(3);
        doublereal dC = pChord->dGet(-1.);
        doublereal dP = pForce->dGet(-1.);
        out
                << " " << pNode1->GetLabel()
                << " ", (f1 + ra1*(dP - dC*3./4.)).Write(out, " ")
                << " ", (f1 + ra1*(dP + dC/4.)).Write(out, " ");

        ra1 = Ra2.GetVec(1);
        ra3 = Ra2.GetVec(3);
        dC = pChord->dGet(1.);
        dP = pForce->dGet(1.);
        out
                << " " << pNode2->GetLabel()
                << " ", (f2 + ra1*(dP - dC*3./4.)).Write(out, " ")
                << " ", (f2 + ra1*(dP + dC/4.)).Write(out, " ")
                << std::endl;

        return pEl;
} /* End of ReadAerodynamicBeam2() */

Here is the call graph for this function:

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

Definition at line 1361 of file aeroelem.cc.

References AerodynamicOutput::AEROD_OUT_NODE, AerodynamicOutput::AEROD_OUT_PGAUSS, AerodynamicOutput::AEROD_OUT_STD, Elem::AERODYNAMIC, bDefaultUseJacobian, grad::cos(), DEBUGCOUTFNAME, DEBUGLCOUT, Shape::dGet(), DataManager::fReadOutput(), WithLabel::GetLabel(), IncludeParser::GetLineData(), DataManager::GetLogFile(), MBDynParser::GetPosRel(), HighParser::GetReal(), MBDynParser::GetRotRel(), HighParser::GetYesNoOrBool(), AeroData::iGetNumDof(), HighParser::IsArg(), HighParser::IsKeyWord(), MBDYN_EXCEPT_ARGS, MYDEBUG_INPUT, AerodynamicOutput::OUTPUT_NONE, ReadAeroData(), ReadInducedVelocity(), DataManager::ReadNode(), ReadOptionalAerodynamicCustomOutput(), SAFENEWWITHCONSTRUCTOR, grad::sin(), Node::STRUCTURAL, UNKNOWN_ORIENTATION_DESCRIPTION, and Write().

Referenced by DataManager::ReadOneElem().

{
        DEBUGCOUTFNAME("ReadAerodynamicBody");

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

        InducedVelocity* pIndVel = 0;
        bool bPassive(false);
        (void)ReadInducedVelocity(pDM, HP, uLabel, "AerodynamicBody",
                pIndVel, bPassive);

        ReferenceFrame RF(pNode);
        Vec3 f(HP.GetPosRel(RF));

        DEBUGLCOUT(MYDEBUG_INPUT, "Offset: " << f << std::endl);

        Mat3x3 Ra(HP.GetRotRel(RF));

        doublereal dSpan = HP.GetReal();
        DEBUGLCOUT(MYDEBUG_INPUT, "Span: " << dSpan << std::endl);

        Shape* pChord = 0;
        Shape* pForce = 0;
        Shape* pVelocity = 0;
        Shape* pTwist = 0;
        Shape* pTipLoss = 0;

        integer iNumber = 0;
        DriveCaller* pDC = 0;
        AeroData* aerodata = 0;

        ReadAeroData(pDM, HP, 1,
                &pChord, &pForce, &pVelocity, &pTwist, &pTipLoss,
                &iNumber, &pDC, &aerodata);

        bool bUseJacobian(false);
        if (HP.IsKeyWord("jacobian")) {
                bUseJacobian = HP.GetYesNoOrBool(bDefaultUseJacobian);
        }

        if (aerodata->iGetNumDof() > 0 && !bUseJacobian) {
                silent_cerr("AerodynamicBody(" << uLabel << "): "
                        "aerodynamic model needs \"jacobian, yes\" at line " << HP.GetLineData()
                        << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        OrientationDescription od = UNKNOWN_ORIENTATION_DESCRIPTION;
        unsigned uFlags = AerodynamicOutput::OUTPUT_NONE;
        ReadOptionalAerodynamicCustomOutput(pDM, HP, uLabel, uFlags, od);

        flag fOut = pDM->fReadOutput(HP, Elem::AERODYNAMIC);
        if (HP.IsArg()) {
                if (HP.IsKeyWord("std")) {
                        fOut |= AerodynamicOutput::AEROD_OUT_STD;
                } else if (HP.IsKeyWord("gauss")) {
                        fOut |= AerodynamicOutput::AEROD_OUT_PGAUSS;
                } else if (HP.IsKeyWord("node")) {
                        fOut |= AerodynamicOutput::AEROD_OUT_NODE;
                } else {
                        silent_cerr("AerodynamicBody(" << uLabel << "): "
                                "unknown output mode at line " << HP.GetLineData()
                                << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

        } else if (fOut) {
                fOut |= AerodynamicOutput::AEROD_OUT_STD;
        }
        fOut |= uFlags;

        Elem* pEl = 0;
        SAFENEWWITHCONSTRUCTOR(pEl,
                AerodynamicBody,
                AerodynamicBody(uLabel, pDO, pNode, pIndVel, bPassive,
                        f, dSpan, Ra,
                        pChord, pForce, pVelocity, pTwist, pTipLoss,
                        iNumber, aerodata, pDC, bUseJacobian, od, fOut));

        /* Se non c'e' il punto e virgola finale */
        if (HP.IsArg()) {
                silent_cerr("semicolon expected at line "
                        << HP.GetLineData() << std::endl);
                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        Vec3 Ra3 = Ra.GetVec(3);
        doublereal dCm1 = pChord->dGet(-1.);
        doublereal dPm1 = pForce->dGet(-1.);
        doublereal dTm1 = pTwist->dGet(-1.);
        Vec3 Ram1 = Ra*Vec3(std::cos(dTm1), std::sin(dTm1), 0.);

        doublereal dCp1 = pChord->dGet(1.);
        doublereal dPp1 = pForce->dGet(1.);
        doublereal dTp1 = pTwist->dGet(1.);
        Vec3 Rap1 = Ra*Vec3(std::cos(dTp1), std::sin(dTp1), 0.);

        std::ostream& out = pDM->GetLogFile();
        out << "aero0: " << uLabel
                << " " << pNode->GetLabel()
                << " ", (f - Ra3*(dSpan/2.) + Ram1*(dPm1 - dCm1*3./4.)).Write(out, " ")
                << " ", (f - Ra3*(dSpan/2.) + Ram1*(dPm1 + dCm1/4.)).Write(out, " ")
                << " ", (f + Ra3*(dSpan/2.) + Rap1*(dPp1 - dCp1*3./4.)).Write(out, " ")
                << " ", (f + Ra3*(dSpan/2.) + Rap1*(dPp1 + dCp1/4.)).Write(out, " ")
                << std::endl;

        return pEl;
} /* End of ReadAerodynamicBody() */

Here is the call graph for this function:

void ReadAerodynamicCustomOutput	(	DataManager *	pDM,
		MBDynParser &	HP,
		unsigned int	uLabel,
		unsigned &	uFlags,
		OrientationDescription &	od
	)

Definition at line 1295 of file aeroelem.cc.

References IncludeParser::GetLineData(), HighParser::IsArg(), HighParser::IsKeyWord(), MBDYN_EXCEPT_ARGS, AerodynamicOutput::OUTPUT_GP_ALL, AerodynamicOutput::OUTPUT_GP_CONFIGURATION, AerodynamicOutput::OUTPUT_GP_F, AerodynamicOutput::OUTPUT_GP_FORCES, AerodynamicOutput::OUTPUT_GP_M, AerodynamicOutput::OUTPUT_GP_R, AerodynamicOutput::OUTPUT_GP_V, AerodynamicOutput::OUTPUT_GP_W, AerodynamicOutput::OUTPUT_GP_X, AerodynamicOutput::OUTPUT_NONE, and ReadOptionalOrientationDescription().

Referenced by DataManager::ReadControl(), and ReadOptionalAerodynamicCustomOutput().

{
        uFlags = AerodynamicOutput::OUTPUT_NONE;

        while (HP.IsArg()) {
                unsigned uFlag;

                if (HP.IsKeyWord("position")) {
                        uFlag = AerodynamicOutput::OUTPUT_GP_X;

                } else if (HP.IsKeyWord("orientation")) {
                        uFlag = AerodynamicOutput::OUTPUT_GP_R;

                } else if (HP.IsKeyWord("velocity")) {
                        uFlag = AerodynamicOutput::OUTPUT_GP_V;

                } else if (HP.IsKeyWord("angular" "velocity")) {
                        uFlag = AerodynamicOutput::OUTPUT_GP_W;

                } else if (HP.IsKeyWord("configuration")) {
                        uFlag = AerodynamicOutput::OUTPUT_GP_CONFIGURATION;

                } else if (HP.IsKeyWord("force")) {
                        uFlag = AerodynamicOutput::OUTPUT_GP_F;

                } else if (HP.IsKeyWord("moment")) {
                        uFlag = AerodynamicOutput::OUTPUT_GP_M;

                } else if (HP.IsKeyWord("forces")) {
                        uFlag = AerodynamicOutput::OUTPUT_GP_FORCES;

                } else if (HP.IsKeyWord("all")) {
                        uFlag = AerodynamicOutput::OUTPUT_GP_ALL;

                } else {
                        break;
                }

                if (uFlags & uFlag) {
                        silent_cerr("AerodynamicElement(" << uLabel << "): "
                                "duplicate custom output "
                                "at line " << HP.GetLineData()
                                << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

                if (uFlag & AerodynamicOutput::OUTPUT_GP_R) {
                        od = ReadOptionalOrientationDescription(pDM, HP);
                }

                uFlags |= uFlag;
        }
}

Here is the call graph for this function:

static bool ReadInducedVelocity ( DataManager *  pDM, MBDynParser &  HP, unsigned  uLabel, const char *  sElemType, InducedVelocity *&  pIndVel, bool &  bPassive  )

static

Definition at line 1216 of file aeroelem.cc.

References ASSERT, DEBUGLCOUT, HighParser::GetInt(), IncludeParser::GetLineData(), Elem::INDUCEDVELOCITY, HighParser::IsKeyWord(), Elem::LOADABLE, MBDYN_EXCEPT_ARGS, MYDEBUG_INPUT, and DataManager::pFindElem().

Referenced by ReadAerodynamicBeam(), ReadAerodynamicBeam2(), and ReadAerodynamicBody().

{
        bool bReadIV(false);
        bool bReadUDIV(false);
        if (HP.IsKeyWord("rotor")) {
                silent_cerr(sElemType << "(" << uLabel << "): "
                        "\"rotor\" keyword is deprecated; "
                        "use \"induced velocity\" instead "
                        "at line " << HP.GetLineData()
                        << std::endl);

                bReadIV = true;

        } else if (HP.IsKeyWord("induced" "velocity")) {
                bReadIV = true;

        } else if (HP.IsKeyWord("user" "defined" "induced" "velocity")) {
                bReadIV = true;
                bReadUDIV = true;
        }

        if (bReadIV) {
                unsigned int uIV = (unsigned int)HP.GetInt();
                DEBUGLCOUT(MYDEBUG_INPUT,
                        "Linked to InducedVelocity(" << uIV << ")" << std::endl);

                bPassive = false;
                if (HP.IsKeyWord("passive")) {
                        bPassive = true;
                }

                /*
                 * verifica di esistenza del rotore
                 * NOTA: ovviamente il rotore deve essere definito
                 * prima dell'elemento aerodinamico
                 */
                Elem* p;

                if (bReadUDIV) {
                        p = pDM->pFindElem(Elem::LOADABLE, uIV);
                        if (p == 0) {
                                silent_cerr(sElemType << "(" << uLabel << "): "
                                        "user-defined InducedVelocity(" << uIV << ") not defined "
                                        "at line " << HP.GetLineData()
                                        << std::endl);
                                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }

                } else {
                        p = pDM->pFindElem(Elem::INDUCEDVELOCITY, uIV);
                        if (p == 0) {
                                // try a user-defined one?
                                p = pDM->pFindElem(Elem::LOADABLE, uIV);
                                if (p == 0 || !dynamic_cast<InducedVelocity *>(p)) {
                                        silent_cerr(sElemType << "(" << uLabel << "): "
                                                "InducedVelocity(" << uIV << ") not defined "
                                                "at line " << HP.GetLineData()
                                                << std::endl);

                                        throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
                                }

                                silent_cerr(sElemType << "(" << uLabel << "): "
                                        "InducedVelocity(" << uIV << ") not defined; using user-defined InducedVelocity(" << uIV << ") "
                                        "at line " << HP.GetLineData()
                                        << std::endl);
                        }
                }

                pIndVel = dynamic_cast<InducedVelocity *>(p);
                ASSERT(pIndVel != 0);
        }

        return (pIndVel != 0);
}

Here is the call graph for this function:

void ReadOptionalAerodynamicCustomOutput	(	DataManager *	pDM,
		MBDynParser &	HP,
		unsigned int	uLabel,
		unsigned &	uFlags,
		OrientationDescription &	od
	)

Definition at line 1351 of file aeroelem.cc.

References Elem::AERODYNAMIC, DataManager::GetOutput(), HighParser::IsKeyWord(), and ReadAerodynamicCustomOutput().

Referenced by ReadAerodynamicBeam(), ReadAerodynamicBeam2(), and ReadAerodynamicBody().

{
        pDM->GetOutput(Elem::AERODYNAMIC, uFlags, od);
        if (HP.IsKeyWord("custom" "output")) {
                ReadAerodynamicCustomOutput(pDM, HP, uLabel, uFlags, od);
        }
}

Here is the call graph for this function:

Variable Documentation

const bool bDefaultUseJacobian = false

static

Definition at line 159 of file aeroelem.cc.

Referenced by ReadAerodynamicBeam(), ReadAerodynamicBeam2(), and ReadAerodynamicBody().

const doublereal d13 = 1./sqrt(3.)

static

Definition at line 1556 of file aeroelem.cc.

const char* elemnames[]

static

Initial value:

= {
        "AerodynamicBody",
        "AerodynamicBeam2",
        "AerodynamicBeam3",
        0
}

Definition at line 258 of file aeroelem.cc.

Referenced by Aerodynamic2DElem< iNN >::DescribeDof(), Aerodynamic2DElem< iNN >::DescribeEq(), and Aerodynamic2DElem< iNN >::OutputPrepare().

const doublereal pdsf2[] = { 0., 1. }

static

Definition at line 2524 of file aeroelem.cc.

Referenced by AerodynamicBeam2::AssJac(), and AerodynamicBeam2::AssVec().

const doublereal pdsf3[] = { -d13, d13, 1. }

static

Definition at line 1558 of file aeroelem.cc.

Referenced by AerodynamicBeam::AssJac(), and AerodynamicBeam::AssVec().

const doublereal pdsi2[] = { -1., 0. }

static

Definition at line 2523 of file aeroelem.cc.

Referenced by AerodynamicBeam2::AssJac(), and AerodynamicBeam2::AssVec().

const doublereal pdsi3[] = { -1., -d13, d13 }

static

Definition at line 1557 of file aeroelem.cc.

Referenced by AerodynamicBeam::AssJac(), and AerodynamicBeam::AssVec().