AerodynamicModal Class Reference

#include <aeromodal.h>

Inheritance diagram for AerodynamicModal:

Collaboration diagram for AerodynamicModal:

Public Types
enum	RigidF_t { NO_RIGID = 0, RIGID = 6 }
Public Types inherited from Elem
enum	Type {   UNKNOWN = -1, AIRPROPERTIES = 0, INDUCEDVELOCITY, AUTOMATICSTRUCTURAL,   GRAVITY, BODY, JOINT, JOINT_REGULARIZATION,   BEAM, PLATE, FORCE, INERTIA,   ELECTRICBULK, ELECTRIC, THERMAL, HYDRAULIC,   BULK, LOADABLE, DRIVEN, EXTERNAL,   AEROMODAL, AERODYNAMIC, GENEL, SOCKETSTREAM_OUTPUT,   RTAI_OUTPUT = SOCKETSTREAM_OUTPUT, LASTELEMTYPE }
Public Types inherited from SimulationEntity
typedef std::vector< Hint * >	Hints
Public Types inherited from ToBeOutput
enum	{ OUTPUT = 0x1U, OUTPUT_MASK = 0xFU, OUTPUT_PRIVATE = 0x10U, OUTPUT_PRIVATE_MASK = ~OUTPUT_MASK }
Public Types inherited from AerodynamicElem
enum	Type {   UNKNOWN = -1, INDUCEDVELOCITY = 0, AEROMODAL, AERODYNAMICBODY,   AERODYNAMICBEAM, AERODYNAMICEXTERNAL, AERODYNAMICEXTERNALMODAL, AERODYNAMICLOADABLE,   AIRCRAFTINSTRUMENTS, GENERICFORCE, LASTAEROTYPE }

Public Member Functions
	AerodynamicModal (unsigned int uLabel, const StructNode *pN, const Modal *pMJ, const Mat3x3 &RaTmp, const DofOwner *pDO, doublereal Cd, const int NModal, const int NAero, RigidF_t rgF, const int Gust, const doublereal Vff, SpMapMatrixHandler *pAMat, FullMatrixHandler *pBMat, FullMatrixHandler *pCMat, FullMatrixHandler *pD0Mat, FullMatrixHandler *pD1Mat, FullMatrixHandler *pD2Mat, flag fout)
	~AerodynamicModal (void)
std::ostream &	Restart (std::ostream &out) const
unsigned int	iGetNumDof (void) const
DofOrder::Order	GetDofType (unsigned int i) const
Elem::Type	GetElemType (void) const
void	WorkSpaceDim (integer *piNumRows, integer *piNumCols) const
VariableSubMatrixHandler &	AssJac (VariableSubMatrixHandler &WorkMat, doublereal dCoef, const VectorHandler &, const VectorHandler &)
SubVectorHandler &	AssRes (SubVectorHandler &WorkVec, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
void	Output (OutputHandler &OH) const
unsigned int	iGetInitialNumDof (void) const
void	InitialWorkSpaceDim (integer *piNumRows, integer *piNumCols) const
VariableSubMatrixHandler &	InitialAssJac (VariableSubMatrixHandler &WorkMat, const VectorHandler &XCurr)
SubVectorHandler &	InitialAssRes (SubVectorHandler &WorkVec, const VectorHandler &XCurr)
AerodynamicElem::Type	GetAerodynamicElemType (void) const
virtual void	GetConnectedNodes (std::vector< const Node * > &connectedNodes) const
Public Member Functions inherited from Elem
	Elem (unsigned int uL, flag fOut)
virtual	~Elem (void)
virtual std::ostream &	DescribeDof (std::ostream &out, const char *prefix="", bool bInitial=false) const
virtual void	DescribeDof (std::vector< std::string > &desc, bool bInitial=false, int i=-1) const
virtual std::ostream &	DescribeEq (std::ostream &out, const char *prefix="", bool bInitial=false) const
virtual void	DescribeEq (std::vector< std::string > &desc, bool bInitial=false, int i=-1) const
virtual void	AssMats (VariableSubMatrixHandler &WorkMatA, VariableSubMatrixHandler &WorkMatB, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual bool	bInverseDynamics (void) const
void	SetInverseDynamicsFlags (unsigned uIDF)
unsigned	GetInverseDynamicsFlags (void) const
bool	bIsErgonomy (void) const
bool	bIsRightHandSide (void) const
virtual VariableSubMatrixHandler &	AssJac (VariableSubMatrixHandler &WorkMat, const VectorHandler &XCurr)
virtual SubVectorHandler &	AssRes (SubVectorHandler &WorkVec, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr, const VectorHandler &XPrimePrimeCurr, InverseDynamics::Order iOrder=InverseDynamics::INVERSE_DYNAMICS)
virtual int	GetNumConnectedNodes (void) const
Public Member Functions inherited from WithLabel
	WithLabel (unsigned int uL=0, const std::string &sN="")
virtual	~WithLabel (void)
void	PutLabel (unsigned int uL)
void	PutName (const std::string &sN)
unsigned int	GetLabel (void) const
const std::string &	GetName (void) const
Public Member Functions inherited from SimulationEntity
	SimulationEntity (void)
virtual	~SimulationEntity (void)
virtual bool	bIsValidIndex (unsigned int i) const
virtual DofOrder::Order	GetEqType (unsigned int i) const
virtual void	SetValue (DataManager *pDM, VectorHandler &X, VectorHandler &XP, SimulationEntity::Hints *h=0)
virtual Hint *	ParseHint (DataManager *pDM, const char *s) const
virtual void	BeforePredict (VectorHandler &, VectorHandler &, VectorHandler &, VectorHandler &) const
virtual void	AfterPredict (VectorHandler &X, VectorHandler &XP)
virtual void	Update (const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual void	DerivativesUpdate (const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual void	Update (const VectorHandler &XCurr, InverseDynamics::Order iOrder)
virtual void	AfterConvergence (const VectorHandler &X, const VectorHandler &XP)
virtual void	AfterConvergence (const VectorHandler &X, const VectorHandler &XP, const VectorHandler &XPP)
virtual unsigned int	iGetNumPrivData (void) const
virtual unsigned int	iGetPrivDataIdx (const char *s) const
virtual doublereal	dGetPrivData (unsigned int i) const
virtual std::ostream &	OutputAppend (std::ostream &out) const
virtual void	ReadInitialState (MBDynParser &HP)
Public Member Functions inherited from ToBeOutput
	ToBeOutput (flag fOut=fDefaultOut)
virtual	~ToBeOutput (void)
virtual void	OutputPrepare (OutputHandler &OH)
virtual void	Output (OutputHandler &OH, const VectorHandler &X, const VectorHandler &XP) const
virtual flag	fToBeOutput (void) const
virtual bool	bToBeOutput (void) const
virtual void	SetOutputFlag (flag f=flag(1))
Public Member Functions inherited from AerodynamicElem
	AerodynamicElem (unsigned int uL, const DofOwner *pDO, flag fOut)
virtual	~AerodynamicElem (void)
virtual bool	NeedsAirProperties (void) const
virtual const InducedVelocity *	pGetInducedVelocity (void) const
Public Member Functions inherited from ElemWithDofs
	ElemWithDofs (unsigned int uL, const DofOwner *pDO, flag fOut)
virtual	~ElemWithDofs (void)
Public Member Functions inherited from DofOwnerOwner
	DofOwnerOwner (const DofOwner *pDO)
virtual	~DofOwnerOwner ()
virtual const DofOwner *	pGetDofOwner (void) const
virtual integer	iGetFirstIndex (void) const
virtual void	SetInitialValue (VectorHandler &X)
Public Member Functions inherited from AirPropOwner
	AirPropOwner (void)
virtual	~AirPropOwner (void)
virtual void	PutAirProperties (const AirProperties *pAP)
virtual flag	fGetAirVelocity (Vec3 &Velocity, const Vec3 &X) const
virtual doublereal	dGetAirDensity (const Vec3 &X) const
virtual doublereal	dGetAirPressure (const Vec3 &X) const
virtual doublereal	dGetAirTemperature (const Vec3 &X) const
virtual doublereal	dGetSoundSpeed (const Vec3 &X) const
virtual bool	GetAirProps (const Vec3 &X, doublereal &rho, doublereal &c, doublereal &p, doublereal &T) const
Public Member Functions inherited from InitialAssemblyElem
	InitialAssemblyElem (unsigned int uL, flag fOut)
virtual	~InitialAssemblyElem (void)
Public Member Functions inherited from SubjectToInitialAssembly
	SubjectToInitialAssembly (void)
virtual	~SubjectToInitialAssembly (void)

Protected Member Functions
void	AssVec (SubVectorHandler &WorkVec)

Protected Attributes
const StructNode *	pModalNode
const Modal *	pModalJoint
Vec3	P0
Mat3x3	R0
const Mat3x3	Ra
doublereal	Chord
unsigned int	NStModes
unsigned int	NAeroStates
unsigned int	NGust
SpMapMatrixHandler *	pA
FullMatrixHandler *	pB
FullMatrixHandler *	pC
FullMatrixHandler *	pD0
FullMatrixHandler *	pD1
FullMatrixHandler *	pD2
MyVectorHandler *	pq
MyVectorHandler *	pqPrime
MyVectorHandler *	pqSec
MyVectorHandler *	pxa
MyVectorHandler *	pxaPrime
MyVectorHandler *	pgs
MyVectorHandler *	pgsPrime
const doublereal	gustVff
const doublereal	gustXi
RigidF_t	RigidF
Protected Attributes inherited from WithLabel
unsigned int	uLabel
std::string	sName
Protected Attributes inherited from ToBeOutput
flag	fOutput
Protected Attributes inherited from AirPropOwner
const AirProperties *	pAirProperties

Detailed Description

Definition at line 69 of file aeromodal.h.

Member Enumeration Documentation

enum AerodynamicModal::RigidF_t

Enumerator
NO_RIGID
RIGID

Definition at line 107 of file aeromodal.h.

                      {
                NO_RIGID = 0,
                RIGID = 6
        };

Constructor & Destructor Documentation

AerodynamicModal::AerodynamicModal	(	unsigned int	uLabel,
		const StructNode *	pN,
		const Modal *	pMJ,
		const Mat3x3 &	RaTmp,
		const DofOwner *	pDO,
		doublereal	Cd,
		const int	NModal,
		const int	NAero,
		RigidF_t	rgF,
		const int	Gust,
		const doublereal	Vff,
		SpMapMatrixHandler *	pAMat,
		FullMatrixHandler *	pBMat,
		FullMatrixHandler *	pCMat,
		FullMatrixHandler *	pD0Mat,
		FullMatrixHandler *	pD1Mat,
		FullMatrixHandler *	pD2Mat,
		flag	fout
	)

Definition at line 43 of file aeromodal.cc.

References ASSERT, DEBUGCOUTFNAME, StructDispNode::GetNodeType(), StructNode::GetRCurr(), StructDispNode::GetXCurr(), Mat3x3::MulTV(), NAeroStates, NGust, NStModes, P0, pgs, pgsPrime, pModalNode, pq, pqPrime, pqSec, pxa, pxaPrime, R0, Ra, RigidF, SAFENEWWITHCONSTRUCTOR, and Node::STRUCTURAL.

: Elem(uLabel, fout),
AerodynamicElem(uLabel, pDO, fout),
InitialAssemblyElem(uLabel, fout),
pModalNode(pN), pModalJoint(pMJ),
Ra(RaTmp),
Chord(Cd),
NStModes(NModal),
NAeroStates(NAero),
NGust(Gust),
pA(pAMat), pB(pBMat), pC(pCMat),
pD0(pD0Mat), pD1(pD1Mat), pD2(pD2Mat),
pq(0), pqPrime(0), pqSec(0),
pxa(0), pxaPrime(0),
pgs(0), pgsPrime(0),
gustVff(Vff), gustXi(0.707),
RigidF(rgF)
{
        DEBUGCOUTFNAME("AerodynamicModal::AerodynamicModal");

        R0 = pModalNode->GetRCurr()*Ra;
        P0 = R0.MulTV(pModalNode->GetXCurr());

        SAFENEWWITHCONSTRUCTOR(pq, MyVectorHandler,
                MyVectorHandler(NStModes+RigidF));
        SAFENEWWITHCONSTRUCTOR(pqPrime, MyVectorHandler,
                MyVectorHandler(NStModes+RigidF));
        SAFENEWWITHCONSTRUCTOR(pqSec, MyVectorHandler,
                MyVectorHandler(NStModes+RigidF));
        SAFENEWWITHCONSTRUCTOR(pxa, MyVectorHandler,
                MyVectorHandler(NAeroStates));
        SAFENEWWITHCONSTRUCTOR(pxaPrime, MyVectorHandler,
                MyVectorHandler(NAeroStates));
        SAFENEWWITHCONSTRUCTOR(pgs, MyVectorHandler,
                MyVectorHandler(2*NGust));
        SAFENEWWITHCONSTRUCTOR(pgsPrime, MyVectorHandler,
                MyVectorHandler(2*NGust));

        ASSERT(pModalNode != 0);
        ASSERT(pModalNode->GetNodeType() == Node::STRUCTURAL);
}

Here is the call graph for this function:

AerodynamicModal::~AerodynamicModal

(

void

)

Definition at line 102 of file aeromodal.cc.

References DEBUGCOUTFNAME, pA, pB, pC, pD0, pD1, pD2, pq, pqPrime, pqSec, pxa, pxaPrime, and SAFEDELETE.

{
        DEBUGCOUTFNAME("AerodynamicModal::~AerodynamicModal");
        if (pq != 0) {
                SAFEDELETE(pq);
        }
        if (pqPrime != 0) {
                SAFEDELETE(pqPrime);
        }
        if (pqSec != 0) {
                SAFEDELETE(pqSec);
        }
        if (pxa != 0) {
                SAFEDELETE(pxa);
        }
        if (pxaPrime != 0) {
                SAFEDELETE(pxaPrime);
        }
        if (pA != 0) {
                SAFEDELETE(pA);
        }
        if (pB != 0) {
                SAFEDELETE(pB);
        }
        if (pC != 0) {
                SAFEDELETE(pC);
        }
        if (pD0 != 0) {
                SAFEDELETE(pD0);
        }
        if (pD1 != 0) {
                SAFEDELETE(pD1);
        }
        if (pD2 != 0) {
                SAFEDELETE(pD2);
        }
}

Member Function Documentation

VariableSubMatrixHandler & AerodynamicModal::AssJac ( VariableSubMatrixHandler &  WorkMat, doublereal  dCoef, const VectorHandler &  XCurr, const VectorHandler &  XPrimeCurr  )

virtual

Implements Elem.

Definition at line 149 of file aeromodal.cc.

References Chord, DEBUGCOUT, AirPropOwner::fGetAirVelocity(), AirPropOwner::GetAirProps(), StructNode::GetRCurr(), StructDispNode::GetVCurr(), StructDispNode::GetXCurr(), gustVff, gustXi, DofOwnerOwner::iGetFirstIndex(), Modal::iGetModalIndex(), FullSubMatrixHandler::IncCoef(), NAeroStates, NGust, NStModes, pA, pB, pC, pD0, pD1, pD2, pModalJoint, pModalNode, FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), Ra, FullSubMatrixHandler::ResizeReset(), RigidF, VariableSubMatrixHandler::SetFull(), WorkSpaceDim(), and Zero3.

Referenced by InitialAssJac().

{
        DEBUGCOUT("Entering AerodynamicModal::AssJac()" << std::endl);

        FullSubMatrixHandler& WM = WorkMat.SetFull();

        integer iNumRows = 0;
        integer iNumCols = 0;
        WorkSpaceDim(&iNumRows, &iNumCols);
        WM.ResizeReset(iNumRows, iNumCols);

        integer iModalIndex = pModalJoint->iGetModalIndex();
        for (unsigned int iCnt = 1; iCnt <= NStModes; iCnt++) {
                WM.PutRowIndex(iCnt, iModalIndex+NStModes+iCnt);
                WM.PutColIndex(iCnt, iModalIndex+iCnt);
                WM.PutColIndex(NStModes+iCnt, iModalIndex+NStModes+iCnt);
        }

        integer iAeroIndex = iGetFirstIndex();
        for (unsigned int iCnt = 1; iCnt <= NAeroStates+2*NGust; iCnt++) {
                WM.PutRowIndex(NStModes+iCnt, iAeroIndex+iCnt);
                WM.PutColIndex(2*NStModes+iCnt, iAeroIndex+iCnt);
        }

        /* Dati del nodo rigido */
        const Vec3& X0(pModalNode->GetXCurr());
        const Vec3& V0(pModalNode->GetVCurr());
        const Mat3x3& Rn(pModalNode->GetRCurr());
        Mat3x3 RR(Rn*Ra);
        Vec3 Vr(V0);

        doublereal rho, vs, p, T;
        GetAirProps(X0, rho, vs, p, T);         /* p, T are not used yet */

        Vec3 VTmp(Zero3);
        if (fGetAirVelocity(VTmp, X0)) {
                Vr -= VTmp;
        }

        /* velocità nel riferimento nodale aerodinamico */
        VTmp = RR.MulTV(Vr);
        doublereal nV = std::abs(VTmp(1));
        /* doublereal CV=Chord/(2*nV); */
        doublereal qd  = 0.5*rho*nV*nV;
        doublereal qd1 = 0.25*rho*nV*Chord; /* qd*CV */
        doublereal qd2 = 0.125*rho*Chord*Chord; /* qd*CV*CV */

        /* parte deformabile :
         *
         * |                                ||aP|
         * |                                ||  |
         * |   cKae     Mae+cCae     -cqC   ||bP|
         * |                                ||  |
         * |-c(1/CV)B     0      I-c(1/CV)A ||xa|
         *
         * con
         * Mae=-qd CV^2 D2
         * Cae=-qd CV D1
         * Kae=-qd D0
         */
        for (unsigned int i = 1; i <= NStModes; i++) {
                for (unsigned int j = 1; j <= NStModes; j++) {
                        WM.IncCoef(i, j,
                                -dCoef*qd*pD0->operator()(RigidF + i, RigidF + j));
                        WM.IncCoef(i, j + NStModes,
                                -qd2*pD2->operator()(RigidF + i, RigidF + j)
                                -dCoef*qd1*pD1->operator()(RigidF + i, RigidF + j));
                }
        }

        for (unsigned int j = 1; j <= NAeroStates; j++) {
                for (unsigned int i = 1; i <= NStModes; i++) {
                        WM.IncCoef(i, j + 2*NStModes,
                                -dCoef*qd*pC->operator()(RigidF + i, j));
                        WM.IncCoef(j + NStModes, i,
                                -dCoef*(2*nV/Chord)*pB->operator()(j, RigidF + i));
                }
        }

        for (unsigned int j = 1; j <= NAeroStates; j++) {
                for (unsigned int i = 1; i <= NAeroStates; i++) {
                        WM.IncCoef(i + NStModes, j + 2*NStModes,
                                1.*(i == j) - dCoef*(2*nV/Chord)*pA->operator()(i, j));
                }
        }

        if (NGust) {
                for (unsigned int i = 1; i <= NStModes; i++) {
                        WM.IncCoef(i, 2*NStModes + NAeroStates + 1,
                                -dCoef*qd*pD0->operator()(RigidF + i, RigidF + NStModes + 1));
                        WM.IncCoef(i, 2*NStModes + NAeroStates + 3,
                                -dCoef*qd*pD0->operator()(RigidF + i, RigidF + NStModes + 2));
                        WM.IncCoef(i, 2*NStModes + NAeroStates + 2,
                                -qd2*pD2->operator()(RigidF + i, RigidF + NStModes + 1)
                                -dCoef*qd1*pD1->operator()(RigidF + i, RigidF + NStModes + 1));
                        WM.IncCoef(i, 2*NStModes + NAeroStates + 4,
                                -qd2*pD2->operator()(RigidF + i, RigidF + NStModes + 2)
                                -dCoef*qd1*pD1->operator()(RigidF + i, RigidF + NStModes + 2));
                }

                for (unsigned int i = 1; i <= NAeroStates; i++) {
                        WM.IncCoef(i + NStModes, 2*NStModes + NAeroStates + 1,
                                -dCoef*(2*nV/Chord)*pB->operator()(i, RigidF + NStModes + 1));
                        WM.IncCoef(i + NStModes, 2*NStModes + NAeroStates + 3,
                                -dCoef*(2*nV/Chord)*pB->operator()(i, RigidF + NStModes + 2));
                }

                for (unsigned int i = 0; i < NGust; i++) {
                        WM.IncCoef(i*2 + 1 + NStModes+NAeroStates,
                                i*2 + 1 + 2*NStModes + NAeroStates, 1.);
                        WM.IncCoef(i*2 + 1 + NStModes+NAeroStates,
                                i*2 + 2 + 2*NStModes + NAeroStates, -dCoef);
                        WM.IncCoef(i*2 + 2 + NStModes+NAeroStates,
                                i*2 + 1 + 2*NStModes+NAeroStates,
                                dCoef*gustVff*gustVff);
                        WM.IncCoef(i*2 + 2 + NStModes + NAeroStates,
                                i*2 + 2 + 2*NStModes + NAeroStates,
                                1. + dCoef*2*gustXi*gustVff);
                }
        }

        return WorkMat;
}

Here is the call graph for this function:

SubVectorHandler & AerodynamicModal::AssRes ( SubVectorHandler &  WorkVec, doublereal  dCoef, const VectorHandler &  XCurr, const VectorHandler &  XPrimeCurr  )

virtual

Implements Elem.

Definition at line 278 of file aeromodal.cc.

References AssVec(), DEBUGCOUTFNAME, StructNode::GetRCurr(), StructDispNode::GetVCurr(), StructNode::GetWCurr(), StructNode::GetWPCurr(), StructDispNode::GetXCurr(), StructDispNode::GetXPPCurr(), DofOwnerOwner::iGetFirstIndex(), StructDispNode::iGetFirstMomentumIndex(), Modal::iGetModalIndex(), NAeroStates, NGust, NStModes, P0, pgs, pgsPrime, pModalJoint, pModalNode, pq, pqPrime, pqSec, MyVectorHandler::Put(), MyVectorHandler::PutCoef(), SubVectorHandler::PutRowIndex(), pxa, pxaPrime, R0, Ra, VectorHandler::ResizeReset(), RigidF, and RotManip::VecRot().

{
        DEBUGCOUTFNAME("AerodynamicModal::AssRes");
        WorkVec.ResizeReset(RigidF + NStModes + NAeroStates + 2*NGust);

        if (RigidF) {
                integer iFirstIndex = pModalNode->iGetFirstMomentumIndex();
                for (int iCnt = 1; iCnt <= RigidF; iCnt++) {
                        WorkVec.PutRowIndex(iCnt, iFirstIndex + iCnt);
                }

                const Vec3& X0(pModalNode->GetXCurr());
                const Mat3x3& Rn(pModalNode->GetRCurr());

                Mat3x3 RR(Rn*Ra);

                // q
                pq->Put(1, RR.MulTV(X0) - P0);
                pq->Put(4, RotManip::VecRot(RR.MulMT(R0))); // nota: wrappa; se serve si può eliminare

                // dot{q}
                const Vec3& V0(pModalNode->GetVCurr());
                const Vec3& W0(pModalNode->GetWCurr());
                pqPrime->Put(1, RR.MulTV(V0 + X0.Cross(W0)));
                pqPrime->Put(4, RR.MulTV(W0));

                // ddot{q}
                const Vec3& XPP0(pModalNode->GetXPPCurr());
                const Vec3& WP0(pModalNode->GetWPCurr());
                pqSec->Put(1, RR.MulTV(XPP0));  // verificare: non mancano i termini di trasporto ecc?
                pqSec->Put(4, RR.MulTV(WP0));
        }

        integer iModalIndex = pModalJoint->iGetModalIndex();

        for (unsigned int iCnt = 1; iCnt <= NStModes; iCnt++) {
                WorkVec.PutRowIndex(RigidF + iCnt, iModalIndex + NStModes + iCnt);
        }

        integer iAeroIndex = iGetFirstIndex();
        for (unsigned int iCnt = 1; iCnt <= NAeroStates + 2*NGust; iCnt++) {
                WorkVec.PutRowIndex(RigidF + NStModes+iCnt, iAeroIndex + iCnt);
        }

        /* Recupera i vettori {a} e {aP} e {aS} (deformate modali) */
        // FIXME: get this info from modal joint?
        for (unsigned int  iCnt = 1; iCnt <= NStModes; iCnt++) {
                pq->PutCoef(iCnt + RigidF, XCurr(iModalIndex + iCnt));
                pqPrime->PutCoef(iCnt + RigidF, XPrimeCurr(iModalIndex + iCnt));
                pqSec->PutCoef(iCnt + RigidF, XPrimeCurr(iModalIndex + NStModes + iCnt));
        }

        /* Recupera i vettori {xa} e {xaP}  */
        for (unsigned int  iCnt = 1; iCnt <= NAeroStates; iCnt++) {
                pxa->PutCoef(iCnt, XCurr(iAeroIndex + iCnt));
                pxaPrime->PutCoef(iCnt, XPrimeCurr(iAeroIndex + iCnt));
        }

        for (unsigned int  iCnt = 1; iCnt <= 2*NGust; iCnt++) {
                pgs->PutCoef(iCnt, XCurr(iAeroIndex + NAeroStates + iCnt));
                pgsPrime->PutCoef(iCnt, XPrimeCurr(iAeroIndex + NAeroStates + iCnt));
        }

        AssVec(WorkVec);

        return WorkVec;
}

Here is the call graph for this function:

void AerodynamicModal::AssVec ( SubVectorHandler &  WorkVec )

protected

Definition at line 430 of file aeromodal.cc.

References VectorHandler::Add(), Chord, DEBUGCOUTFNAME, AirPropOwner::fGetAirVelocity(), AirPropOwner::GetAirProps(), StructNode::GetRCurr(), StructDispNode::GetVCurr(), StructDispNode::GetXCurr(), gustVff, gustXi, VectorHandler::IncCoef(), MyVectorHandler::IncCoef(), MatrixHandler::MatVecIncMul(), NAeroStates, NGust, NStModes, pA, pB, pC, pD0, pD1, pD2, pgs, pgsPrime, pModalNode, pq, pqPrime, pqSec, Vec3::Put(), pxa, pxaPrime, Ra, MyVectorHandler::Reset(), RigidF, and Zero3.

Referenced by AssRes(), and InitialAssRes().

{
        DEBUGCOUTFNAME("AerodynamicModal::AssVec");

        /* Dati del nodo rigido */
        const Vec3& X0(pModalNode->GetXCurr());
        Vec3 V0(pModalNode->GetVCurr());
        const Mat3x3& Rn(pModalNode->GetRCurr());
        Mat3x3 RR(Rn*Ra);
        Vec3 Vr(V0);

        Vr(1) = 0.0; // perturbazione di velocita', verificare

        doublereal rho, vs, p, T;
        GetAirProps(X0, rho, vs, p, T); /* p, T are not used yet */

        Vec3 Vair(Zero3);
        if (fGetAirVelocity(Vair, X0)) {
                Vr -= Vair;
        }
        /* velocita' nel riferimento nodale aerodinamico */
        V0 = RR.MulTV(Vr);
        doublereal nV = std::abs(V0(1));
        MyVectorHandler TmpA(NAeroStates);
        TmpA.Reset();
        MyVectorHandler Tmp(RigidF + NStModes);
        Tmp.Reset();
        MyVectorHandler TmpP(RigidF + NStModes);
        TmpP.Reset();
        MyVectorHandler TmpS(RigidF + NStModes);
        TmpS.Reset();

        // FIXME: use matrix/vector product?
        for (unsigned int i = 1; i <= NAeroStates; i++) {
                for (unsigned int j = 1; j <= NStModes + RigidF; j++) {
                        TmpA.IncCoef(i, pB->operator()(i, j)*(pq->operator()(j)));
                }
        }

        pA->MatVecIncMul(TmpA, *pxa);

        /* doublereal CV = Chord/(2*nV); */
        doublereal qd = 0.5*rho*nV*nV;
        doublereal qd1 = 0.25*rho*nV*Chord; /* qd*CV */
        doublereal qd2 = 0.125*rho*Chord*Chord; /* qd*CV*CV */

        for (unsigned int i = 1; i <= NAeroStates; i++) {
                WorkVec.IncCoef(RigidF + NStModes + i, -pxaPrime->operator()(i) + (2*nV/Chord)*(TmpA(i)));
        }

        for (unsigned int i = 1; i <= NStModes + RigidF; i++) {
                for (unsigned int j = 1; j <= NStModes + RigidF; j++) {
                        Tmp.IncCoef(i, pD0->operator()(i, j)*(pq->operator()(j)));
                        TmpP.IncCoef(i, pD1->operator()(i, j)*(pqPrime->operator()(j)));
                        TmpS.IncCoef(i, pD2->operator()(i, j)*(pqSec->operator()(j)));
                }
        }

        for (unsigned int i = 1; i <= NStModes + RigidF; i++) {
                for (unsigned int j = 1; j <= NAeroStates; j++) {
                        Tmp.IncCoef(i, pC->operator()(i, j)*(pxa->operator()(j)));
                }
        }

        if (RigidF) {
                Vec3 F(Zero3);
                Vec3 M(Zero3);

                for (unsigned int i = 1; i <= 3; i++) {
                        F.Put(i, qd*Tmp(i) + qd1*TmpP(i) + qd2*TmpS(i));
                        M.Put(i, qd*Tmp(i + 3) + qd1*TmpP(i + 3) + qd2*TmpS(i + 3));
                }

                // std::cout << F << std::endl;
                F = RR*F;
                M = RR*M;

#if 0
                // tolti perche' altrimenti inizia un Dutch Roll
                F(1) = 0.;
                F(2) = 0.;
                M(1) = 0.;
                M(3) = 0.;
#endif

                WorkVec.Add(1, F);
                WorkVec.Add(4, M);
        }

        for (unsigned int i = 1 + RigidF; i <= NStModes + RigidF; i++) {
                WorkVec.IncCoef(i, qd*Tmp(i) + qd1*TmpP(i) + qd2*TmpS(i));
        }

        if (NGust) {
#if 0 /* unused */
                doublereal Vyg = Vair(2)/nV;
                doublereal Vzg = Vair(3)/nV;
#endif
                for (unsigned int i = 1; i <= NAeroStates; i++) {
                        WorkVec.IncCoef(RigidF + NStModes + i,
                                (2*nV/Chord)*pB->operator()(i, RigidF + NStModes + 1)*pgs->operator()(1)
                                + (2*nV/Chord)*pB->operator()(i, RigidF + NStModes + 2)*pgs->operator()(3));
                }

                if (RigidF) {
                        Vec3 F(Zero3);
                        Vec3 M(Zero3);

                        for (unsigned int i = 1; i <= 3; i++) {
                                F.Put(i, qd*pD0->operator()(i, RigidF + NStModes + 1)*pgs->operator()(1)
                                        + qd*pD0->operator()(i, RigidF + NStModes + 2)*pgs->operator()(3)
                                        + qd1*pD1->operator()(i, RigidF + NStModes + 1)*pgs->operator()(2)
                                        + qd1*pD1->operator()(i, RigidF + NStModes + 2)*pgs->operator()(4)
                                        + qd2*pD2->operator()(i, RigidF + NStModes + 1)*pgsPrime->operator()(2)
                                        + qd2*pD2->operator()(i, RigidF + NStModes + 2)*pgsPrime->operator()(4));
                                M.Put(i, qd*pD0->operator()(i + 3, RigidF + NStModes + 1)*pgs->operator()(1)
                                        + qd*pD0->operator()(i + 3, RigidF + NStModes + 2)*pgs->operator()(3)
                                        + qd1*pD1->operator()(i + 3, RigidF + NStModes + 1)*pgs->operator()(2)
                                        + qd1*pD1->operator()(i + 3, RigidF + NStModes + 2)*pgs->operator()(4)
                                        + qd2*pD2->operator()(i + 3, RigidF + NStModes + 1)*pgsPrime->operator()(2)
                                        + qd2*pD2->operator()(i + 3, RigidF + NStModes + 2)*pgsPrime->operator()(4));
                        }

                        // std::cout << F << std::endl;
                        F = RR*F;
                        M = RR*M;
                        WorkVec.Add(1, F);
                        WorkVec.Add(4, M);
                }

                for (unsigned int i = 1 + RigidF; i <= NStModes + RigidF; i++) {
                        WorkVec.IncCoef(i, qd*pD0->operator()(i, RigidF + NStModes + 1)*pgs->operator()(1)
                                + qd*pD0->operator()(i, RigidF + NStModes + 2)*pgs->operator()(3)
                                + qd1*pD1->operator()(i, RigidF + NStModes + 1)*pgs->operator()(2)
                                + qd1*pD1->operator()(i, RigidF + NStModes + 2)*pgs->operator()(4)
                                + qd2*pD2->operator()(i, RigidF + NStModes + 1)*pgsPrime->operator()(2)
                                + qd2*pD2->operator()(i, RigidF + NStModes + 2)*pgsPrime->operator()(4));
                }

                for (unsigned int i = 0; i < NGust; i++) {
                        WorkVec.IncCoef(RigidF + i*2 + 1 + NStModes + NAeroStates,
                                -pgsPrime->operator()(1 + i*2) + pgs->operator()(2 + i*2));
                        if (nV != 0) {
                                WorkVec.IncCoef(RigidF + i*2 + 2 + NStModes + NAeroStates,
                                        -pgsPrime->operator()(2 + i*2)
                                        - gustVff*gustVff*pgs->operator()(1 + i*2)
                                        - 2*gustXi*gustVff*pgs->operator()(2 + i*2)
                                        + gustVff*gustVff*Vair(2 + i)/nV);

                        } else {
                                WorkVec.IncCoef(RigidF + i*2 + 2 + NStModes + NAeroStates,
                                        -pgsPrime->operator()(2 + i*2)
                                        - gustVff*gustVff*pgs->operator()(1 + i*2)
                                        - 2*gustXi*gustVff*pgs->operator()(2 + i*2));
                        }
                }

                // std::cout << X0 << std::endl;
                // std::cout << WorkVec(3) << std::endl;
        }
}

Here is the call graph for this function:

AerodynamicElem::Type AerodynamicModal::GetAerodynamicElemType ( void  ) const

inlinevirtual

Implements AerodynamicElem.

Definition at line 213 of file aeromodal.h.

References AerodynamicElem::AEROMODAL.

                                                               {
                return AerodynamicElem::AEROMODAL;
        };

virtual void AerodynamicModal::GetConnectedNodes ( std::vector< const Node * > &  connectedNodes ) const

inlinevirtual

Reimplemented from Elem.

Definition at line 221 of file aeromodal.h.

References pModalNode.

                                                                         {
                connectedNodes.resize(1);
                connectedNodes[0] = pModalNode;
        };

DofOrder::Order AerodynamicModal::GetDofType ( unsigned int  i ) const

inlinevirtual

Reimplemented from Elem.

Definition at line 149 of file aeromodal.h.

References ASSERT, DofOrder::DIFFERENTIAL, and iGetNumDof().

                                                       {
                /* gradi di liberta' differenziali (eq. modali) */
                ASSERT(i < iGetNumDof());
                return DofOrder::DIFFERENTIAL;
        };

Here is the call graph for this function:

Elem::Type AerodynamicModal::GetElemType ( void  ) const

inlinevirtual

Implements Elem.

Definition at line 156 of file aeromodal.h.

References Elem::AEROMODAL.

                                         {
                return Elem::AEROMODAL;
        };

unsigned int AerodynamicModal::iGetInitialNumDof ( void  ) const

inlinevirtual

Implements SubjectToInitialAssembly.

Definition at line 187 of file aeromodal.h.

References NAeroStates, and NGust.

                                                   {
                return NAeroStates + NGust*2;
        };

unsigned int AerodynamicModal::iGetNumDof ( void  ) const

inlinevirtual

Reimplemented from Elem.

Definition at line 144 of file aeromodal.h.

References NAeroStates, and NGust.

Referenced by GetDofType().

                                            {
                return NAeroStates + NGust*2;
        };

VariableSubMatrixHandler& AerodynamicModal::InitialAssJac ( VariableSubMatrixHandler &  WorkMat, const VectorHandler &  XCurr  )

inlinevirtual

Implements SubjectToInitialAssembly.

Definition at line 200 of file aeromodal.h.

References AssJac(), and DEBUGCOUTFNAME.

        {
                DEBUGCOUTFNAME("AerodynamicModal::InitialAssJac");
                AssJac(WorkMat, 0, XCurr, XCurr);
                return WorkMat;
        };

Here is the call graph for this function:

SubVectorHandler & AerodynamicModal::InitialAssRes ( SubVectorHandler &  WorkVec, const VectorHandler &  XCurr  )

virtual

Implements SubjectToInitialAssembly.

Definition at line 350 of file aeromodal.cc.

References AssVec(), grad::atan(), DEBUGCOUTFNAME, StructNode::GetRCurr(), StructDispNode::GetVCurr(), StructNode::GetWCurr(), StructDispNode::GetXCurr(), DofOwnerOwner::iGetFirstIndex(), Modal::iGetModalIndex(), NAeroStates, NGust, NStModes, P0, CGR_Rot::Param, pgs, pgsPrime, pModalJoint, pModalNode, pq, pqPrime, pqSec, MyVectorHandler::Put(), MyVectorHandler::PutCoef(), SubVectorHandler::PutRowIndex(), pxa, pxaPrime, R0, Ra, VectorHandler::ResizeReset(), RigidF, and Zero3.

{
        DEBUGCOUTFNAME("AerodynamicModal::AssRes");
        const Vec3& X0(pModalNode->GetXCurr());
        const Mat3x3& Rn(pModalNode->GetRCurr());
        R0 = Rn*Ra;
        P0 = R0.MulTV(X0);

        WorkVec.ResizeReset(RigidF + NStModes + NAeroStates + 2*NGust);

        if (RigidF) {
                integer iFirstIndex = pModalNode->iGetFirstIndex();
                for (int iCnt = 1; iCnt <= RigidF; iCnt++) {
                        WorkVec.PutRowIndex(iCnt, iFirstIndex + RigidF + iCnt);
                }

                const Vec3& X0(pModalNode->GetXCurr());
                const Mat3x3& Rn(pModalNode->GetRCurr());
                Mat3x3 RR(Rn*Ra);

                // q
                pq->Put(1, RR.MulTV(X0) - P0);

                Vec3 g(CGR_Rot::Param, RR.MulMT(R0));
                doublereal d(g.Norm());
                if (d > std::numeric_limits<doublereal>::epsilon()) {
                        pq->Put(4, g*(2./d*atan(d/2.)));

                } else {
                        pq->Put(4, Zero3);
                }

                // dot{q}
                const Vec3& V0(pModalNode->GetVCurr());
                const Vec3& W0(pModalNode->GetWCurr());
                pqPrime->Put(1, RR.MulTV(V0));
                pqPrime->Put(4, RR.MulTV(W0));

                // ddot{q}
                pqPrime->Put(1, Zero3);
                pqPrime->Put(4, Zero3);
        }

        integer iModalIndex = pModalJoint->iGetModalIndex();

        for (unsigned int iCnt = 1; iCnt <= NStModes; iCnt++) {
                WorkVec.PutRowIndex(RigidF + iCnt, iModalIndex + NStModes + iCnt);
        }

        integer iAeroIndex = iGetFirstIndex();
        for (unsigned int iCnt = 1; iCnt <= NAeroStates + 2*NGust; iCnt++) {
                WorkVec.PutRowIndex(RigidF + NStModes + iCnt, iAeroIndex + iCnt);
        }

        /* Recupera i vettori {a} e {aP} e {aS} (deformate modali) */
        for (unsigned int iCnt = 1; iCnt <= NStModes; iCnt++) {
                pq->PutCoef(iCnt + RigidF, XCurr(iModalIndex + iCnt));
                pqPrime->PutCoef(iCnt + RigidF, 0);
                pqSec->PutCoef(iCnt + RigidF, 0);
        }

        /* Recupera i vettori {xa} e {xaP}  */
        for (unsigned int iCnt = 1; iCnt <= NAeroStates; iCnt++) {
                pxa->PutCoef(iCnt, XCurr(iAeroIndex + iCnt));
                pxaPrime->PutCoef(iCnt, 0);
        }

        for (unsigned int iCnt = 1; iCnt <= 2*NGust; iCnt++) {
                pgs->PutCoef(iCnt, XCurr(iAeroIndex + NAeroStates + iCnt));
                pgsPrime->PutCoef(iCnt, 0.);
        }

        AssVec(WorkVec);

        return WorkVec;
}

Here is the call graph for this function:

void AerodynamicModal::InitialWorkSpaceDim ( integer *  piNumRows, integer *  piNumCols  ) const

inlinevirtual

Implements SubjectToInitialAssembly.

Definition at line 193 of file aeromodal.h.

References NAeroStates, NGust, and NStModes.

                                                                          {
                *piNumRows = NAeroStates + NStModes + NGust*2;
                *piNumCols = NAeroStates + 2*NStModes + NGust*2;
        };

void AerodynamicModal::Output ( OutputHandler &  OH ) const

virtual

Reimplemented from ToBeOutput.

Definition at line 595 of file aeromodal.cc.

References OutputHandler::AeroModals(), ToBeOutput::bToBeOutput(), WithLabel::GetLabel(), NAeroStates, NGust, pxa, and pxaPrime.

{
        if (bToBeOutput()) {
                OH.AeroModals() << std::setw(8) << GetLabel() << " ";
                for (unsigned int iCnt = 1; iCnt <= NAeroStates; iCnt++) {
                        OH.AeroModals() << " " << pxa->operator()(iCnt);
                }
                for (unsigned int iCnt = 1; iCnt <= NAeroStates + 2*NGust; iCnt++) {
                        OH.AeroModals() << " " << pxaPrime->operator()(iCnt);
                }
                OH.AeroModals() << std::endl;
        }
}

Here is the call graph for this function:

std::ostream & AerodynamicModal::Restart ( std::ostream &  out ) const

virtual

Implements Elem.

Definition at line 142 of file aeromodal.cc.

{
        return out << "  /* aerodynamic modal: not implemented yet */" << std::endl;
}

void AerodynamicModal::WorkSpaceDim ( integer *  piNumRows, integer *  piNumCols  ) const

inlinevirtual

Implements Elem.

Definition at line 163 of file aeromodal.h.

References NAeroStates, NGust, and NStModes.

Referenced by AssJac().

                                                                        {
                *piNumRows = NAeroStates + NStModes + NGust*2;
                *piNumCols = NAeroStates + 2*NStModes + NGust*2;
        };

Member Data Documentation

doublereal AerodynamicModal::Chord

protected

Definition at line 81 of file aeromodal.h.

Referenced by AssJac(), and AssVec().

const doublereal AerodynamicModal::gustVff

protected

Definition at line 103 of file aeromodal.h.

Referenced by AssJac(), and AssVec().

const doublereal AerodynamicModal::gustXi

protected

Definition at line 104 of file aeromodal.h.

Referenced by AssJac(), and AssVec().

unsigned int AerodynamicModal::NAeroStates

protected

Definition at line 83 of file aeromodal.h.

Referenced by AerodynamicModal(), AssJac(), AssRes(), AssVec(), iGetInitialNumDof(), iGetNumDof(), InitialAssRes(), InitialWorkSpaceDim(), Output(), and WorkSpaceDim().

unsigned int AerodynamicModal::NGust

protected

Definition at line 84 of file aeromodal.h.

Referenced by AerodynamicModal(), AssJac(), AssRes(), AssVec(), iGetInitialNumDof(), iGetNumDof(), InitialAssRes(), InitialWorkSpaceDim(), Output(), and WorkSpaceDim().

unsigned int AerodynamicModal::NStModes

protected

Definition at line 82 of file aeromodal.h.

Referenced by AerodynamicModal(), AssJac(), AssRes(), AssVec(), InitialAssRes(), InitialWorkSpaceDim(), and WorkSpaceDim().

Vec3 AerodynamicModal::P0

protected

Definition at line 77 of file aeromodal.h.

Referenced by AerodynamicModal(), AssRes(), and InitialAssRes().

SpMapMatrixHandler* AerodynamicModal::pA

protected

Definition at line 86 of file aeromodal.h.

Referenced by AssJac(), AssVec(), and ~AerodynamicModal().

FullMatrixHandler* AerodynamicModal::pB

protected

Definition at line 87 of file aeromodal.h.

Referenced by AssJac(), AssVec(), and ~AerodynamicModal().

FullMatrixHandler* AerodynamicModal::pC

protected

Definition at line 88 of file aeromodal.h.

Referenced by AssJac(), AssVec(), and ~AerodynamicModal().

FullMatrixHandler* AerodynamicModal::pD0

protected

Definition at line 89 of file aeromodal.h.

Referenced by AssJac(), AssVec(), and ~AerodynamicModal().

FullMatrixHandler* AerodynamicModal::pD1

protected

Definition at line 90 of file aeromodal.h.

Referenced by AssJac(), AssVec(), and ~AerodynamicModal().

FullMatrixHandler* AerodynamicModal::pD2

protected

Definition at line 91 of file aeromodal.h.

Referenced by AssJac(), AssVec(), and ~AerodynamicModal().

MyVectorHandler* AerodynamicModal::pgs

protected

Definition at line 100 of file aeromodal.h.

Referenced by AerodynamicModal(), AssRes(), AssVec(), and InitialAssRes().

MyVectorHandler* AerodynamicModal::pgsPrime

protected

Definition at line 101 of file aeromodal.h.

Referenced by AerodynamicModal(), AssRes(), AssVec(), and InitialAssRes().

const Modal* AerodynamicModal::pModalJoint

protected

Definition at line 76 of file aeromodal.h.

Referenced by AssJac(), AssRes(), and InitialAssRes().

const StructNode* AerodynamicModal::pModalNode

protected

Definition at line 75 of file aeromodal.h.

Referenced by AerodynamicModal(), AssJac(), AssRes(), AssVec(), GetConnectedNodes(), and InitialAssRes().

MyVectorHandler* AerodynamicModal::pq

protected

Definition at line 93 of file aeromodal.h.

Referenced by AerodynamicModal(), AssRes(), AssVec(), InitialAssRes(), and ~AerodynamicModal().

MyVectorHandler* AerodynamicModal::pqPrime

protected

Definition at line 94 of file aeromodal.h.

Referenced by AerodynamicModal(), AssRes(), AssVec(), InitialAssRes(), and ~AerodynamicModal().

MyVectorHandler* AerodynamicModal::pqSec

protected

Definition at line 95 of file aeromodal.h.

Referenced by AerodynamicModal(), AssRes(), AssVec(), InitialAssRes(), and ~AerodynamicModal().

MyVectorHandler* AerodynamicModal::pxa

protected

Definition at line 97 of file aeromodal.h.

Referenced by AerodynamicModal(), AssRes(), AssVec(), InitialAssRes(), Output(), and ~AerodynamicModal().

MyVectorHandler* AerodynamicModal::pxaPrime

protected

Definition at line 98 of file aeromodal.h.

Referenced by AerodynamicModal(), AssRes(), AssVec(), InitialAssRes(), Output(), and ~AerodynamicModal().

Mat3x3 AerodynamicModal::R0

protected

Definition at line 78 of file aeromodal.h.

Referenced by AerodynamicModal(), AssRes(), and InitialAssRes().

const Mat3x3 AerodynamicModal::Ra

protected

Definition at line 79 of file aeromodal.h.

Referenced by AerodynamicModal(), AssJac(), AssRes(), AssVec(), and InitialAssRes().

RigidF_t AerodynamicModal::RigidF

protected

Definition at line 113 of file aeromodal.h.

Referenced by AerodynamicModal(), AssJac(), AssRes(), AssVec(), and InitialAssRes().

---

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

• aeromodal.h
• aeromodal.cc