ExtRigidForceEDGE Class Reference

#include <modaledge.h>

Inheritance diagram for ExtRigidForceEDGE:

Collaboration diagram for ExtRigidForceEDGE:

Public Member Functions
	ExtRigidForceEDGE (DataManager *pDM)
unsigned	Recv (ExtFileHandlerBase *pEFH, unsigned uFlags, unsigned &uLabel, Vec3 &f, Vec3 &m, std::vector< doublereal > &a)
void	Send (ExtFileHandlerBase *pEFH, unsigned uFlags, unsigned uLabel, const Vec3 &x, const Mat3x3 &R, const Vec3 &v, const Vec3 &w, const std::vector< doublereal > &q, const std::vector< doublereal > &qP)
Public Member Functions inherited from ExtForceEDGE
	ExtForceEDGE (DataManager *pDM)
virtual	~ExtForceEDGE (void)
bool	Prepare (ExtFileHandlerBase *pEFH, unsigned uLabel, bool bRigid, unsigned uModes)
Public Member Functions inherited from ExtModalForceBase
virtual	~ExtModalForceBase (void)

Additional Inherited Members
Public Types inherited from ExtModalForceBase
enum	BitMask {   EMF_NONE = 0x0U, EMF_RIGID = 0x1U, EMF_MODAL = 0x2U, EMF_ALL = (EMF_RIGID | EMF_MODAL),   EMF_RIGID_DETECT = 0x10U, EMF_MODAL_DETECT = 0x20U, EMF_RIGID_DETECT_MASK = (EMF_RIGID_DETECT | EMF_RIGID), EMF_MODAL_DETECT_MASK = (EMF_MODAL_DETECT | EMF_MODAL),   EMF_DETECT_MASK = (EMF_RIGID_DETECT | EMF_MODAL_DETECT), EMF_ERR = 0x10000000U }
Protected Attributes inherited from ExtForceEDGE
const AirProperties *	pAP

Detailed Description

Definition at line 58 of file modaledge.h.

Constructor & Destructor Documentation

ExtRigidForceEDGE::ExtRigidForceEDGE

(

DataManager *

pDM

)

Definition at line 71 of file modaledge.cc.

References NO_OP.

: ExtForceEDGE(pDM)
{
        NO_OP;
}

Member Function Documentation

unsigned ExtRigidForceEDGE::Recv ( ExtFileHandlerBase *  pEFH, unsigned  uFlags, unsigned &  uLabel, Vec3 &  f, Vec3 &  m, std::vector< doublereal > &  a  )

virtual

Implements ExtModalForceBase.

Definition at line 88 of file modaledge.cc.

References ASSERT, buf, ExtModalForceBase::EMF_ALL, ExtModalForceBase::EMF_NONE, ExtModalForceBase::EMF_RIGID, ExtModalForceBase::EMF_RIGID_DETECT_MASK, ExtFileHandlerBase::GetInStream(), MBDYN_EXCEPT_ARGS, mbedge_eat_field(), mbedge_eat_sep(), mbedge_goto_eol(), and STRLENOF.

{
        std::istream& fin = *pEFH->GetInStream();
        unsigned uOutFlags = ExtModalForceBase::EMF_NONE;

        ASSERT((uFlags & ExtModalForceBase::EMF_RIGID) != 0);
        ASSERT((uFlags & ~ExtModalForceBase::EMF_ALL) == 0);

        // fake
        uLabel = 0;

        // cycle
        unsigned lineno = 0;
        while (true) {
                char buf[BUFSIZ], *p;
                if (mbedge_goto_eol(fin, buf, sizeof(buf))) {
                        if (!fin) {
                                break;
                        }
                }

                lineno++;

                if (buf[0] == '*') {
                        continue;
                }

                if (strncasecmp(buf, "body_forces", STRLENOF("body_forces")) == 0) {
                        size_t buflen = sizeof(buf) - STRLENOF("body_forces");
                        p = &buf[0] + STRLENOF("body_forces");

                        p = mbedge_eat_sep(p, buflen);
                        if (p == 0) {
                                silent_cerr("ExtRigidForceEDGE: unable to skip separator "
                                        "at line=" << lineno << ", \"" << buf[sizeof(buf) - buflen] << "\"" << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }

                        p = mbedge_eat_field(p, buflen, "R");
                        if (p == 0) {
                                silent_cerr("ExtRigidForceEDGE: unable to skip field \"R\" "
                                        "at line=" << lineno << ", \"" << buf[sizeof(buf) - buflen] << "\"" << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }

                        p = mbedge_eat_sep(p, buflen);
                        if (p == 0) {
                                silent_cerr("ExtRigidForceEDGE: unable to skip separator "
                                        "at line=" << lineno << ", \"" << buf[sizeof(buf) - buflen] << "\"" << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }

                        p = mbedge_eat_field(p, buflen, "1");
                        if (p == 0) {
                                silent_cerr("ExtRigidForceEDGE: unable to skip field \"1\" "
                                        "at line=" << lineno << ", \"" << buf[sizeof(buf) - buflen] << "\"" << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }

                        p = mbedge_eat_sep(p, buflen);
                        if (p == 0) {
                                silent_cerr("ExtRigidForceEDGE: unable to skip separator "
                                        "at line=" << lineno << ", \"" << buf[sizeof(buf) - buflen] << "\"" << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }

                        p = mbedge_eat_field(p, buflen, "6");
                        if (p == 0) {
                                silent_cerr("ExtRigidForceEDGE: unable to skip field \"6\" "
                                        "at line=" << lineno << ", \"" << buf[sizeof(buf) - buflen] << "\"" << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }

                        p = mbedge_eat_sep(p, buflen);
                        if (p == 0) {
                                silent_cerr("ExtRigidForceEDGE: unable to skip separator "
                                        "at line=" << lineno << ", \"" << buf[sizeof(buf) - buflen] << "\"" << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }

                        p = mbedge_eat_field(p, buflen, "0");
                        if (p == 0) {
                                silent_cerr("ExtRigidForceEDGE: unable to skip field \"0\" "
                                        "at line=" << lineno << ", \"" << buf[sizeof(buf) - buflen] << "\"" << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }

                        p = mbedge_eat_sep(p, buflen);
                        if (p == 0) {
                                silent_cerr("ExtRigidForceEDGE: unable to skip separator "
                                        "at line=" << lineno << ", \"" << buf[sizeof(buf) - buflen] << "\"" << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }

                        if (p[0] != '\0' && p[0] != '\n') {
                                silent_cerr("ExtRigidForceEDGE: no line terminator "
                                        "at line=" << lineno << ", \"" << p << "\"" << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }

                        if (uFlags & ExtModalForceBase::EMF_RIGID_DETECT_MASK) {
                                doublereal b[6];
                                for (int i = 0; i < 6; i++) {
                                        fin >> b[i];
                                }
                                f = Vec3(&b[0]);
                                m = Vec3(&b[3]);

                                uOutFlags |= (uFlags & ExtModalForceBase::EMF_RIGID_DETECT_MASK);
                        }

                        // skip to eol
                        if (mbedge_goto_eol(fin, buf, sizeof(buf))) {
                                silent_cerr("ExtRigidForceEDGE: unable to skip to end of line "
                                        "at line=" << lineno << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }
                        continue;
                }

                silent_cerr("ExtRigidForceEDGE: unexpected line=" << lineno << ", "
                        "\"" << buf << "\"" << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        return uOutFlags;
}

Here is the call graph for this function:

void ExtRigidForceEDGE::Send ( ExtFileHandlerBase *  pEFH, unsigned  uFlags, unsigned  uLabel, const Vec3 &  x, const Mat3x3 &  R, const Vec3 &  v, const Vec3 &  w, const std::vector< doublereal > &  q, const std::vector< doublereal > &  qP  )

virtual

Implements ExtModalForceBase.

Definition at line 245 of file modaledge.cc.

References ASSERT, ExtModalForceBase::EMF_ALL, ExtModalForceBase::EMF_RIGID, ExtFileHandlerBase::GetOutStream(), AirProperties::GetVelocity(), Mat3x3::MulTV(), ExtForceEDGE::pAP, and R.

{
        std::ostream& fout = *pEFH->GetOutStream();
        ASSERT((uFlags & ExtModalForceBase::EMF_RIGID) != 0);
        ASSERT((uFlags & ~ExtModalForceBase::EMF_ALL) == 0);

        Vec3 vv = v;
        Vec3 v_infty;

        if (pAP && pAP->GetVelocity(x, v_infty)) {
                vv -= v_infty;
        }

        Vec3 vB = R.MulTV(vv);
        Vec3 wB = R.MulTV(w);

        fout << "* MBDyn to EDGE rigid body dynamics\n"
                "body_dynamics,N,0,0,3\n"
                "* Body linear velocity in body axes\n"
                "VRELV,R,1,3,0\n"
                << vB(1) << " " << vB(2) << " " << vB(3) << "\n"
                "* Body angular velocity in body axes\n"
                "VRELM,R,1,3,0\n"
                << wB(1) << " " << wB(2) << " " << wB(3) << "\n"
                "* Body reference frame cosines (listed by columns)\n"
                "OMGMAN,R,3,3,0\n"
                // NOTE: this is the transpose of the node's R
                << R(1, 1) << " " << R(2, 1) << " " << R(3, 1) << "\n"
                << R(1, 2) << " " << R(2, 2) << " " << R(3, 2) << "\n"
                << R(1, 3) << " " << R(2, 3) << " " << R(3, 3) << "\n";
}

Here is the call graph for this function:

---

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

• modaledge.h
• modaledge.cc