HBeam Class Reference

#include <hbeam.h>

Inheritance diagram for HBeam:

Collaboration diagram for HBeam:

Classes
class	ErrGeneric

Public Member Functions
	HBeam (unsigned int uL, const StructNode *pN1, const StructNode *pN2, const Vec3 &F1, const Vec3 &F2, const Mat3x3 &R1, const Mat3x3 &R2, const ConstitutiveLaw6D *pd, flag fOut)
virtual	~HBeam (void)
virtual Beam::Type	GetBeamType (void) const
virtual Elem::Type	GetElemType (void) const
virtual std::ostream &	Restart (std::ostream &out) const
virtual void	AfterConvergence (const VectorHandler &X, const VectorHandler &XP)
virtual void	WorkSpaceDim (integer *piNumRows, integer *piNumCols) const
void	SetValue (DataManager *pDM, VectorHandler &, VectorHandler &, SimulationEntity::Hints *ph=0)
virtual void	AfterPredict (VectorHandler &, VectorHandler &)
virtual SubVectorHandler &	AssRes (SubVectorHandler &WorkVec, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual VariableSubMatrixHandler &	AssJac (VariableSubMatrixHandler &WorkMat, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual void	Output (OutputHandler &OH) const
virtual unsigned int	iGetInitialNumDof (void) const
virtual void	InitialWorkSpaceDim (integer *piNumRows, integer *piNumCols) const
virtual void	SetInitialValue (VectorHandler &)
virtual VariableSubMatrixHandler &	InitialAssJac (VariableSubMatrixHandler &WorkMat, const VectorHandler &XCurr)
virtual SubVectorHandler &	InitialAssRes (SubVectorHandler &WorkVec, const VectorHandler &XCurr)
virtual unsigned int	iGetNumPrivData (void) const
virtual unsigned int	iGetPrivDataIdx (const char *s) const
virtual doublereal	dGetPrivData (unsigned int i) const
virtual const StructNode *	pGetNode (unsigned int i) 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 unsigned int	iGetNumDof (void) const
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 DofOrder::Order	GetDofType (unsigned int) 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 Hint *	ParseHint (DataManager *pDM, const char *s) const
virtual void	BeforePredict (VectorHandler &, VectorHandler &, VectorHandler &, VectorHandler &) const
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, const VectorHandler &XPP)
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 ElemGravityOwner
	ElemGravityOwner (unsigned int uL, flag fOut)
virtual	~ElemGravityOwner (void)
virtual doublereal	dGetM (void) const
Vec3	GetS (void) const
Mat3x3	GetJ (void) const
Vec3	GetB (void) const
Vec3	GetG (void) const
Public Member Functions inherited from GravityOwner
	GravityOwner (void)
virtual	~GravityOwner (void)
void	PutGravity (const Gravity *pG)
virtual bool	bGetGravity (const Vec3 &X, Vec3 &Acc) 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 Types
enum	NodeName { NODE1 = 0, NODE2 = 1, NUMNODES = 2 }
enum	Deformations { STRAIN = 0, CURVAT = 1, NUMDEFORM = 2 }

Protected Member Functions
virtual void	AssStiffnessMat (FullSubMatrixHandler &WMA, FullSubMatrixHandler &WMB, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual void	AssStiffnessVec (SubVectorHandler &WorkVec, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual void	AddInternalForces (Vec6 &)
virtual void	DsDxi (void)
virtual void	Omega0 (void)
void	SetBeamType (Beam::Type T)
virtual std::ostream &	Restart_ (std::ostream &out) const
Protected Member Functions inherited from ElemGravityOwner
virtual Vec3	GetS_int (void) const
virtual Mat3x3	GetJ_int (void) const
virtual Vec3	GetB_int (void) const
virtual Vec3	GetG_int (void) const

Protected Attributes
const StructNode *	pNode [NUMNODES]
Vec3	f [NUMNODES]
Vec3	fRef [NUMNODES]
Mat3x3	Rn [NUMNODES]
Mat3x3	R
Mat3x3	RRef
ConstitutiveLaw6DOwner *	pD
Mat6x6	DRef
Vec3	Omega
Vec3	OmegaRef
Vec6	Az
Vec6	AzRef
Vec6	AzLoc
Vec6	DefLoc
Vec6	DefLocRef
Vec3	p
Vec3	g
Vec3	L0
Vec3	L
Vec3	LRef
doublereal	xi
doublereal	dsdxi
doublereal	dxids
bool	bFirstRes
Vec3	Rho
Vec3	Rho0
Mat3x3	pdp [NUMNODES]
Mat3x3	pdP [NUMNODES]
Mat3x3	Ldp [NUMNODES]
Mat3x3	LdP [NUMNODES]
Mat3x3	RdP [NUMNODES]
Mat3x3	RhodP [NUMNODES]
Protected Attributes inherited from WithLabel
unsigned int	uLabel
std::string	sName
Protected Attributes inherited from ToBeOutput
flag	fOutput
Protected Attributes inherited from GravityOwner
Gravity *	pGravity

Private Attributes
Beam::Type	BeamT

Friends
class	AerodynamicBeam

Additional Inherited Members
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 }

Detailed Description

Definition at line 57 of file hbeam.h.

Member Enumeration Documentation

enum HBeam::Deformations

protected

Enumerator
STRAIN
CURVAT
NUMDEFORM

Definition at line 72 of file hbeam.h.

{ STRAIN = 0, CURVAT = 1, NUMDEFORM = 2 };

enum HBeam::NodeName

protected

Enumerator
NODE1
NODE2
NUMNODES

Definition at line 71 of file hbeam.h.

{ NODE1 = 0, NODE2 = 1, NUMNODES = 2 };

Constructor & Destructor Documentation

HBeam::HBeam	(	unsigned int	uL,
		const StructNode *	pN1,
		const StructNode *	pN2,
		const Vec3 &	F1,
		const Vec3 &	F2,
		const Mat3x3 &	R1,
		const Mat3x3 &	R2,
		const ConstitutiveLaw6D *	pd,
		flag	fOut
	)

Definition at line 67 of file hbeam.cc.

References ASSERT, Az, AzLoc, AzRef, DefLoc, DefLocRef, DsDxi(), f, g, StructDispNode::GetNodeType(), L, L0, NODE1, NODE2, Omega, p, pD, pNode, Rn, SAFENEWWITHCONSTRUCTOR, Node::STRUCTURAL, Zero3, and Zero6.

: Elem(uL, fOut),
ElemGravityOwner(uL, fOut),
InitialAssemblyElem(uL, fOut),
bFirstRes(true)
{
        /* Validazione dati */
        ASSERT(pN1 != NULL);
        ASSERT(pN1->GetNodeType() == Node::STRUCTURAL);
        ASSERT(pN2 != NULL);
        ASSERT(pN2->GetNodeType() == Node::STRUCTURAL);

        pNode[NODE1] = pN1;
        pNode[NODE2] = pN2;
        f[NODE1] = F1;
        f[NODE2] = F2;

        Rn[NODE1] = R1;
        Rn[NODE2] = R2;

        /*
        RRef = R = interpolazione di R1 e R2 ;
         */

        pD = NULL;
        SAFENEWWITHCONSTRUCTOR(pD,
                        ConstitutiveLaw6DOwner,
                        ConstitutiveLaw6DOwner(pd));

        Omega = Zero3;
        Az = Zero6;
        AzRef = Zero6;
        AzLoc = Zero6;
        DefLoc = Zero6;
        DefLocRef = Zero6;
        p = Zero3;
        g = Zero3;
        L0 = Zero3;
        L = Zero3;

        DsDxi();
}

Here is the call graph for this function:

HBeam::~HBeam ( void  )

virtual

Definition at line 118 of file hbeam.cc.

References ASSERT, pD, and SAFEDELETE.

{
        /* Distrugge il legame costitutivo */
        ASSERT(pD != NULL);
        if (pD != NULL) {
                SAFEDELETE(pD);
        }
}

Member Function Documentation

virtual void HBeam::AddInternalForces ( Vec6 &  )

inlineprotectedvirtual

Definition at line 159 of file hbeam.h.

References NO_OP.

Referenced by AfterPredict(), and AssStiffnessVec().

                               {
        NO_OP;
    };

void HBeam::AfterConvergence ( const VectorHandler &  X, const VectorHandler &  XP  )

virtual

Reimplemented from SimulationEntity.

Definition at line 347 of file hbeam.cc.

References ConstitutiveLawOwner< T, Tder >::AfterConvergence(), DefLoc, and pD.

{
        pD->AfterConvergence(DefLoc);
}

Here is the call graph for this function:

void HBeam::AfterPredict ( VectorHandler &  , VectorHandler &    )

virtual

Reimplemented from SimulationEntity.

Definition at line 577 of file hbeam.cc.

References AddInternalForces(), Az, AzLoc, AzRef, bFirstRes, DefLoc, DefLocRef, DRef, f, g, ConstitutiveLawOwner< T, Tder >::GetF(), ConstitutiveLawOwner< T, Tder >::GetFDE(), StructNode::GetgRef(), StructNode::GetRRef(), StructDispNode::GetXCurr(), L, L0, LRef, CGR_Rot::MatG, CGR_Rot::MatR, MultRMRt(), MultRV(), NODE1, NODE2, NUMNODES, p, pD, pNode, R, RRef, and ConstitutiveLawOwner< T, Tder >::Update().

{
        /*
         * Calcola le deformazioni, aggiorna il legame costitutivo
         * e crea la FDE
         */
#if 0
        /* Recupera i dati dei nodi */
        Vec3   gNod[NUMNODES];
        Vec3   xTmp[NUMNODES];

        for (unsigned int i = 0; i < NUMNODES; i++) {
                gNod[i] = pNode[i]->GetgRef();
                xTmp[i] = pNode[i]->GetXCurr()+pNode[i]->GetRRef()*f[i];
        }

        Mat3x3 RDelta;
        Vec3 gGrad;

        /* Aggiorna le grandezze della trave nel punto di valutazione */

        /* Posizione */
        p = InterpState(xTmp[NODE1], xTmp[NODE2]);

        /* Matrici di rotazione */
        g = InterpState(gNod[NODE1], gNod[NODE2]);
        RDelta = Mat3x3(MatR, g);
        R = RRef = RDelta*R;

        /* Derivate della posizione */
        L = LRef = InterpDeriv(xTmp[NODE1], xTmp[NODE2]);

        /* Derivate dei parametri di rotazione */
        gGrad = InterpDeriv(gNod[NODE1], gNod[NODE2]);

        /* Per le deformazioni nel sistema del materiale */
        Mat3x3 RTmp(R.Transpose());

        /*
         * Calcola le deformazioni nel sistema locale
         * nei punti di valutazione
         */
        DefLoc = DefLocRef = Vec6(RTmp*L-L0,
                        RTmp*(Mat3x3(MatG, g)*gGrad)+DefLoc.GetVec2());

        /* Calcola le azioni interne */
        pD->Update(DefLoc);
        AzLoc = pD->GetF();

        /* corregge le azioni interne locali (piezo, ecc) */
        AddInternalForces(AzLoc);

        /* Porta le azioni interne nel sistema globale */
        Az = AzRef = MultRV(AzLoc, R);

        /* Aggiorna il legame costitutivo di riferimento */
        DRef = MultRMRt(pD->GetFDE(), RRef);

        bFirstRes = true;
#endif
}

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 497 of file hbeam.cc.

References AssStiffnessMat(), DEBUGCOUTFNAME, StructDispNode::iGetFirstMomentumIndex(), StructDispNode::iGetFirstPositionIndex(), NODE1, NODE2, pNode, FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), FullSubMatrixHandler::ResizeReset(), and VariableSubMatrixHandler::SetFull().

{
        DEBUGCOUTFNAME("HBeam::AssJac => AssStiffnessMat");

        integer iNode1FirstMomIndex = pNode[NODE1]->iGetFirstMomentumIndex();
        integer iNode1FirstPosIndex = pNode[NODE1]->iGetFirstPositionIndex();
        integer iNode2FirstMomIndex = pNode[NODE2]->iGetFirstMomentumIndex();
        integer iNode2FirstPosIndex = pNode[NODE2]->iGetFirstPositionIndex();

        FullSubMatrixHandler& WM = WorkMat.SetFull();

        /* Dimensiona la matrice, la azzera e pone gli indici corretti */
        WM.ResizeReset(12, 12);

        for (int iCnt = 1; iCnt <= 6; iCnt++) {
                WM.PutRowIndex(iCnt, iNode1FirstMomIndex + iCnt);
                WM.PutColIndex(iCnt, iNode1FirstPosIndex + iCnt);
                WM.PutRowIndex(6 + iCnt, iNode2FirstMomIndex + iCnt);
                WM.PutColIndex(6 + iCnt, iNode2FirstPosIndex + iCnt);
        }

        AssStiffnessMat(WM, WM, dCoef, XCurr, XPrimeCurr);

        return WorkMat;
}

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 529 of file hbeam.cc.

References AssStiffnessVec(), DEBUGCOUTFNAME, StructDispNode::iGetFirstMomentumIndex(), NODE1, NODE2, pNode, SubVectorHandler::PutRowIndex(), and VectorHandler::ResizeReset().

{
        DEBUGCOUTFNAME("HBeam::AssRes => AssStiffnessVec");

        integer iNode1FirstMomIndex = pNode[NODE1]->iGetFirstMomentumIndex();
        integer iNode2FirstMomIndex = pNode[NODE2]->iGetFirstMomentumIndex();

        /* Dimensiona il vettore, lo azzera e pone gli indici corretti */
        WorkVec.ResizeReset(12);

        for (unsigned int iCnt = 1; iCnt <= 6; iCnt++) {
                WorkVec.PutRowIndex(iCnt, iNode1FirstMomIndex + iCnt);
                WorkVec.PutRowIndex(6 + iCnt, iNode2FirstMomIndex + iCnt);
        }

        AssStiffnessVec(WorkVec, dCoef, XCurr, XPrimeCurr);

        return WorkVec;
}

Here is the call graph for this function:

void HBeam::AssStiffnessMat ( FullSubMatrixHandler &  WMA, FullSubMatrixHandler &  WMB, doublereal  dCoef, const VectorHandler &  XCurr, const VectorHandler &  XPrimeCurr  )

protectedvirtual

Definition at line 354 of file hbeam.cc.

References FullSubMatrixHandler::Add(), AzRef, ComputeFullInterpolation(), Vec3::Cross(), grad::Cross(), DEBUGCOUTFNAME, DRef, dxids, f, ConstitutiveLawOwner< T, Tder >::GetFDE(), StructNode::GetRCurr(), Vec6::GetVec1(), Vec6::GetVec2(), StructDispNode::GetXCurr(), L, Ldp, LdP, MatCross, MultRMRt(), NODE1, NODE2, NUMNODES, p, pD, pdp, pdP, pNode, R, RdP, Rho, RhodP, Rn, FullSubMatrixHandler::Sub(), Mat6x6::SubMat12(), xi, and Zero3x3.

Referenced by AssJac(), and InitialAssJac().

{
        DEBUGCOUTFNAME("HBeam::AssStiffnessMat");

        /*
         * La matrice arriva gia' dimensionata
         * e con gli indici di righe e colonne a posto
         */

        /* Recupera i dati dei nodi */
        Vec3 yTmp[NUMNODES];
        Vec3 fTmp[NUMNODES];
        Mat3x3 RTmp[NUMNODES];
        for (unsigned int i = 0; i < NUMNODES; i++) {
                RTmp[i] = pNode[i]->GetRCurr()*Rn[i];
                fTmp[i] = pNode[i]->GetRCurr()*f[i];
                yTmp[i] = pNode[i]->GetXCurr() + fTmp[i];
        }
        ComputeFullInterpolation(yTmp, RTmp, fTmp,
                        xi, dxids,
                        p, R,
                        RdP, pdP, pdp,
                        L, Rho,
                        RhodP, LdP, Ldp);
        /* Legame costitutivo (viene generato sempre) */
        DRef = MultRMRt(pD->GetFDE(), R);

        /* Derivate delle deformazioni rispetto alle incognite nodali */
        Mat6x6 AzTmp[NUMNODES];

        for (unsigned int i = 0; i < NUMNODES; i++) {
                /* Delta - deformazioni */
                AzTmp[i] = Mat6x6(Ldp[i]*dCoef,
                                Zero3x3,
                                (LdP[i] + L.Cross(RdP[i]))*dCoef,
                                (RhodP[i] + Rho.Cross(RdP[i]))*dCoef);

                /* Delta - azioni interne */
                AzTmp[i] = DRef*AzTmp[i];

                /* Correggo per la rotazione da locale a globale */
                AzTmp[i].SubMat12((AzRef.GetVec1()*dCoef).Cross(RdP[i]));
                AzTmp[i].SubMat22((AzRef.GetVec2()*dCoef).Cross(RdP[i]));
        }

        Vec3 bTmp[2];

        bTmp[0] = p - pNode[NODE1]->GetXCurr();
        bTmp[1] = p - pNode[NODE2]->GetXCurr();

        for (unsigned int i = 0; i < NUMNODES; i++) {
                /* Equazione all'indietro: */
                WMA.Sub(1, 6*i + 1, AzTmp[i].GetMat11());
                WMA.Sub(1, 6*i + 4,
                        AzTmp[i].GetMat12()
                        - (AzRef.GetVec1()*dCoef).Cross(pdP[i]));

                WMA.Sub(4, 6*i + 1,
                        AzTmp[i].GetMat21()
                        - (AzRef.GetVec1()*dCoef).Cross(pdp[i])
                        + bTmp[0].Cross(AzTmp[i].GetMat11()));
                WMA.Sub(4, 6*i + 4,
                        AzTmp[i].GetMat22()
                        + bTmp[0].Cross(AzTmp[i].GetMat12()));

                /* Equazione in avanti: */
                WMA.Add(7, 6*i + 1, AzTmp[i].GetMat11());
                WMA.Add(7, 6*i + 4,
                        AzTmp[i].GetMat12()
                        - (AzRef.GetVec1()*dCoef).Cross(pdP[i]));

                WMA.Add(10, 6*i + 1,
                        AzTmp[i].GetMat21()
                        - (AzRef.GetVec1()*dCoef).Cross(pdp[i])
                        + bTmp[1].Cross(AzTmp[i].GetMat11()));
                WMA.Add(10, 6*i + 4,
                        AzTmp[i].GetMat22()
                        + bTmp[1].Cross(AzTmp[i].GetMat12()));
        }

        /* correzione alle equazioni */
        Mat3x3 FTmp(MatCross, AzRef.GetVec1()*dCoef);
        WMA.Sub(4, 1, FTmp);
        WMA.Add(10, 7, FTmp);
};

Here is the call graph for this function:

void HBeam::AssStiffnessVec ( SubVectorHandler &  WorkVec, doublereal  dCoef, const VectorHandler &  XCurr, const VectorHandler &  XPrimeCurr  )

protectedvirtual

Definition at line 447 of file hbeam.cc.

References VectorHandler::Add(), AddInternalForces(), Az, AzLoc, ComputeInterpolation(), grad::Cross(), DEBUGCOUTFNAME, DefLoc, dxids, f, ConstitutiveLawOwner< T, Tder >::GetF(), StructNode::GetRCurr(), Vec6::GetVec1(), Vec6::GetVec2(), StructDispNode::GetXCurr(), L, L0, MultRV(), NODE1, NODE2, NUMNODES, p, pD, pNode, R, Rho, Rho0, Rn, VectorHandler::Sub(), Mat3x3::Transpose(), ConstitutiveLawOwner< T, Tder >::Update(), and xi.

Referenced by AssRes(), and InitialAssRes().

{
        DEBUGCOUTFNAME("HBeam::AssStiffnessVec");

        /*
         * Riceve il vettore gia' dimensionato e con gli indici a posto
         * per scrivere il residuo delle equazioni di equilibrio dei tre nodi
         */

        /* Recupera i dati dei nodi */
        Vec3 yTmp[NUMNODES];
        Vec3 fTmp[NUMNODES];
        Mat3x3 RTmp[NUMNODES];
        for (unsigned int i = 0; i < NUMNODES; i++) {
                RTmp[i] = pNode[i]->GetRCurr()*Rn[i];
                fTmp[i] = pNode[i]->GetRCurr()*f[i];
                yTmp[i] = pNode[i]->GetXCurr() + fTmp[i];
        }

        /* Interpolazione generica */
        ComputeInterpolation(yTmp, RTmp, fTmp,
                        xi, dxids,
                        p, R,
                        L, Rho);

        Mat3x3 RT(R.Transpose());
        DefLoc = Vec6(RT*L - L0, RT*Rho - Rho0);

        /* Calcola le azioni interne */
        pD->Update(DefLoc);
        AzLoc = pD->GetF();

        /* corregge le azioni interne locali (piezo, ecc) */
        AddInternalForces(AzLoc);

        /* Porta le azioni interne nel sistema globale */
        Az = MultRV(AzLoc, R);

        WorkVec.Add(1, Az.GetVec1());
        WorkVec.Add(4, (p - pNode[NODE1]->GetXCurr()).Cross(Az.GetVec1()) + Az.GetVec2());
        WorkVec.Sub(7, Az.GetVec1());
        WorkVec.Sub(10, Az.GetVec2() + (p - pNode[NODE2]->GetXCurr()).Cross(Az.GetVec1()));
}

Here is the call graph for this function:

doublereal HBeam::dGetPrivData ( unsigned int  i ) const

virtual

Reimplemented from SimulationEntity.

Definition at line 197 of file hbeam.cc.

References ASSERT, AzLoc, DefLoc, Vec6::dGet(), WithLabel::GetLabel(), and MBDYN_EXCEPT_ARGS.

{
        ASSERT(i > 0 && i <= 6);

        switch (i) {
        case 1:
        case 4:
        case 5:
        case 6:
                return DefLoc.dGet(i);

        case 7:
        case 10:
        case 11:
        case 12:
                return AzLoc.dGet(i);

        case 2:
        case 3:
                silent_cerr("HBeam(" << GetLabel() << "): "
                        "not allowed to return shear strain" << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);

        case 8:
        case 9:
                silent_cerr("HBeam(" << GetLabel() << "): "
                        "not allowed to return shear force" << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);

        default:
                silent_cerr("HBeam(" << GetLabel() << "): "
                        "illegal private data " << i << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
}

Here is the call graph for this function:

void HBeam::DsDxi ( void  )

protectedvirtual

Definition at line 234 of file hbeam.cc.

References ComputeInterpolation(), Vec3::Dot(), dsdxi, dxids, f, WithLabel::GetLabel(), StructNode::GetRCurr(), StructDispNode::GetXCurr(), L, L0, LRef, MBDYN_EXCEPT_ARGS, NUMNODES, p, pNode, R, Rho, Rho0, Rn, RRef, grad::sqrt(), Mat3x3::Transpose(), and xi.

Referenced by HBeam().

{
        /* Calcola il ds/dxi e le deformazioni iniziali */
        Mat3x3 RTmp[NUMNODES];
        Vec3 yTmp[NUMNODES];
        Vec3 fTmp[NUMNODES];
        for (unsigned int i = 0; i < NUMNODES; i++) {
                RTmp[i] = pNode[i]->GetRCurr()*Rn[i];
                fTmp[i] = pNode[i]->GetRCurr()*f[i];
                yTmp[i] = pNode[i]->GetXCurr() + fTmp[i];
        }

        xi = 0.5;
        dsdxi = 1.0;
        /* Calcolo i wder ... */
        ComputeInterpolation(yTmp, RTmp, fTmp,
                        xi, dsdxi,
                        p, R,
                        L, Rho);

        doublereal d = L.Dot();
        if (d > std::numeric_limits<doublereal>::epsilon()) {
                d = std::sqrt(d);
        } else {
                silent_cerr("HBeam(" << GetLabel() << ") "
                        "has singular metric; aborting..." << std::endl);

                throw ErrNullNorm(MBDYN_EXCEPT_ARGS);
        }

        dsdxi = dsdxi*d;
        dxids = 1./dsdxi;

        /* Calcolo le caratteristiche iniziali ... */
        ComputeInterpolation(yTmp, RTmp, fTmp,
                        xi, dxids,
                        p, R,
                        L, Rho);

        /* Grandezze iniziali e di riferimento */
        /* FIXME: fare un temporaneo per i trasposti ... */
        RRef = R;
        Mat3x3 RT(R.Transpose());
        Rho0 = RT*Rho;
        LRef = L;
        L0 = RT*L;
}

Here is the call graph for this function:

virtual Beam::Type HBeam::GetBeamType ( void  ) const

inlinevirtual

Definition at line 193 of file hbeam.h.

References Beam::ELASTIC.

                                             {
        return Beam::ELASTIC;
    };

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

inlinevirtual

Reimplemented from Elem.

Definition at line 301 of file hbeam.h.

References NUMNODES, and pNode.

                                                                     {
        connectedNodes.resize(NUMNODES);
        for (int i = 0; i < NUMNODES; i++) {
            connectedNodes[i] = pNode[i];
        }
    };

virtual Elem::Type HBeam::GetElemType ( void  ) const

inlinevirtual

Implements Elem.

Definition at line 198 of file hbeam.h.

References Elem::BEAM.

                                             {
        return Elem::BEAM;
    };

virtual unsigned int HBeam::iGetInitialNumDof ( void  ) const

inlinevirtual

Implements SubjectToInitialAssembly.

Definition at line 259 of file hbeam.h.

                                                       {
        return 0;
    };

unsigned int HBeam::iGetNumPrivData ( void  ) const

virtual

Reimplemented from SimulationEntity.

Definition at line 129 of file hbeam.cc.

{
        return 12;
}

unsigned int HBeam::iGetPrivDataIdx ( const char *  s ) const

virtual

Reimplemented from SimulationEntity.

Definition at line 135 of file hbeam.cc.

References ASSERT.

{
        ASSERT(s != NULL);

        /*
         * {ex|k{x|y|z}}
         */

        unsigned int idx = 0;

        switch (s[0]) {
        case 'F':
                idx += 6;
        case 'e':
                switch (s[1]) {
                case 'x':
                        idx += 1;
                        break;

                case 'y':
                case 'z':
                        return 0;

                default:
                        return 0;
                }
                break;

        case 'M':
                idx += 6;
        case 'k':
                idx += 3;
                switch (s[1]) {
                case 'x':
                        idx += 1;
                        break;

                case 'y':
                        idx += 2;
                        break;

                case 'z':
                        idx += 3;
                        break;

                default:
                        return 0;
                }
                break;

        default:
                return 0;
        }

        if (s[2] != '\0') {
                return 0;
        }

        return idx;
}

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

virtual

Implements SubjectToInitialAssembly.

Definition at line 656 of file hbeam.cc.

References AssStiffnessMat(), DEBUGCOUTFNAME, StructDispNode::iGetFirstPositionIndex(), NODE1, NODE2, pNode, FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), FullSubMatrixHandler::ResizeReset(), and VariableSubMatrixHandler::SetFull().

{
        DEBUGCOUTFNAME("HBeam::InitialAssJac => AssStiffnessMat");

        /* Dimensiona la matrice, la azzera e pone gli indici corretti */
        FullSubMatrixHandler& WM = WorkMat.SetFull();
        WM.ResizeReset(12, 12);

        integer iNode1FirstPosIndex = pNode[NODE1]->iGetFirstPositionIndex();
        integer iNode2FirstPosIndex = pNode[NODE2]->iGetFirstPositionIndex();

        for (int iCnt = 1; iCnt <= 6; iCnt++) {
                WM.PutRowIndex(iCnt, iNode1FirstPosIndex+iCnt);
                WM.PutColIndex(iCnt, iNode1FirstPosIndex+iCnt);
                WM.PutRowIndex(6+iCnt, iNode2FirstPosIndex+iCnt);
                WM.PutColIndex(6+iCnt, iNode2FirstPosIndex+iCnt);
        }

        AssStiffnessMat(WM, WM, 1., XCurr, XCurr);
        return WorkMat;
}

Here is the call graph for this function:

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

virtual

Implements SubjectToInitialAssembly.

Definition at line 682 of file hbeam.cc.

References AssStiffnessVec(), DEBUGCOUTFNAME, StructDispNode::iGetFirstPositionIndex(), NODE1, NODE2, pNode, SubVectorHandler::PutRowIndex(), and VectorHandler::ResizeReset().

{
        DEBUGCOUTFNAME("HBeam::InitialAssRes => AssStiffnessVec");

        /* Dimensiona il vettore, lo azzera e pone gli indici corretti */
        WorkVec.ResizeReset(12);

        integer iNode1FirstPosIndex = pNode[NODE1]->iGetFirstPositionIndex();
        integer iNode2FirstPosIndex = pNode[NODE2]->iGetFirstPositionIndex();

        for (int iCnt = 1; iCnt <= 6; iCnt++) {
                WorkVec.PutRowIndex(iCnt, iNode1FirstPosIndex+iCnt);
                WorkVec.PutRowIndex(6+iCnt, iNode2FirstPosIndex+iCnt);
        }

        AssStiffnessVec(WorkVec, 1., XCurr, XCurr);
        return WorkVec;
}

Here is the call graph for this function:

virtual void HBeam::InitialWorkSpaceDim ( integer *  piNumRows, integer *  piNumCols  ) const

inlinevirtual

Implements SubjectToInitialAssembly.

Definition at line 268 of file hbeam.h.

                                                  {
        *piNumRows = 12;
        *piNumCols = 12;
    };

void HBeam::Omega0 ( void  )

protectedvirtual

Definition at line 285 of file hbeam.cc.

References dN2, StructNode::GetRCurr(), StructNode::GetWCurr(), CGR_Rot::MatG, CGR_Rot::MatGm1, MatR2gparam(), NODE1, NODE2, NUMNODES, Omega, pNode, and grad::Transpose().

{
#if 0
        /* Modo consistente: */
        Mat3x3 RNod[NUMNODES];
        Vec3 w[NUMNODES];
        for (unsigned int i = 0; i < NUMNODES; i++) {
                RNod[i] = pNode[i]->GetRCurr();
                w[i] = pNode[i]->GetWCurr();
        }

        /*
         * Calcolo i parametri di rotazione della rotazione relativa
         * tra inizio e fine e li dimezzo nell'ipotesi che siano limitati
         */
        Vec3 gTmp(MatR2gparam(RNod[NODE2].Transpose()*RNod[NODE1]));

        /*
         * Le derivate dei parametri di rotazione si ricavano da omega
         */
        Vec3 g1P(Mat3x3(MatGm1, gTmp*(-.5))*w[NODE1]);
        Vec3 g2P(Mat3x3(MatGm1, gTmp*.5)*w[NODE2]);

        Vec3 gPTmp(g1P*dN2[NODE1]+g2P*dN2[NODE2]);
        Omega = Mat3x3(MatG, gTmp)*gPTmp;

#if 0
        /* Modo brutale: interpolo le velocita' dei nodi */
        Vec3 w[NUMNODES];
        for (unsigned int i = 0; i < NUMNODES; i++) {
                w[i] = pNode[i]->GetWCurr();
        }
        Omega[i] = w[NODE1]*dN2[NODE1]+w[NODE2]*dN2[NODE2];
#endif /* 0 */
#endif
}

Here is the call graph for this function:

void HBeam::Output ( OutputHandler &  OH ) const

virtual

Reimplemented from ToBeOutput.

Definition at line 645 of file hbeam.cc.

References AzLoc, OutputHandler::Beams(), ToBeOutput::bToBeOutput(), WithLabel::GetLabel(), Vec6::GetVec1(), and Vec6::GetVec2().

{
        if (bToBeOutput()) {
                OH.Beams() << std::setw(8) << GetLabel() << " "
                        << AzLoc.GetVec1() << " " << AzLoc.GetVec2() << std::endl;
        }
}

Here is the call graph for this function:

const StructNode * HBeam::pGetNode ( unsigned int  i ) const

virtual

Definition at line 704 of file hbeam.cc.

References ASSERT, MBDYN_EXCEPT_ARGS, and pNode.

{
        ASSERT(i >= 1 && i <= 2);
        switch (i) {
        case 1:
        case 2:
                return pNode[i-1];
        default:
                throw HBeam::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
}

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

virtual

Implements Elem.

Definition at line 325 of file hbeam.cc.

References Restart_().

{
        return Restart_(out)<< ';' << std::endl;
}

Here is the call graph for this function:

std::ostream & HBeam::Restart_ ( std::ostream &  out ) const

protectedvirtual

Definition at line 331 of file hbeam.cc.

References f, WithLabel::GetLabel(), Mat3x3::GetVec(), NUMNODES, pD, ConstitutiveLawOwner< T, Tder >::pGetConstLaw(), pNode, R, Write(), and Vec3::Write().

Referenced by Restart().

{
        out << "  beam2: " << GetLabel();
        for (unsigned int i = 0; i < NUMNODES; i++) {
                out << ", " << pNode[i]->GetLabel() << ", reference, node, ",
                f[i].Write(out, ", ");
        }
        out << ", reference, global,"
                << "1, ", (R.GetVec(1)).Write(out, ", ") << ", "
                << "2, ", (R.GetVec(2)).Write(out, ", ") << ", ",
        pD->pGetConstLaw()->Restart(out);

        return out;
}

Here is the call graph for this function:

void HBeam::SetBeamType ( Beam::Type  T )

inlineprotected

Definition at line 173 of file hbeam.h.

References BeamT.

                                 {
        BeamT = T;
    };

virtual void HBeam::SetInitialValue ( VectorHandler &  )

inlinevirtual

Definition at line 275 of file hbeam.h.

References NO_OP.

                                                    {
        NO_OP;
    };

void HBeam::SetValue ( DataManager *  pDM, VectorHandler &  , VectorHandler &  , SimulationEntity::Hints *  ph = 0  )

virtual

Implements ElemGravityOwner.

Definition at line 555 of file hbeam.cc.

References Az, AzRef, bFirstRes, DefLoc, DefLocRef, DRef, ConstitutiveLawOwner< T, Tder >::GetFDE(), L, LRef, MultRMRt(), pD, R, and RRef.

{
        /* Aggiorna le grandezze della trave nei punti di valutazione */
        (Mat3x3&)RRef = R;
        (Vec3&)LRef = L;
        (Vec6&)DefLocRef = DefLoc;
        (Vec6&)AzRef = Az;

        /*
         * Aggiorna il legame costitutivo di riferimento
         * (la deformazione e' gia' stata aggiornata dall'ultimo residuo)
         */
        (Mat6x6&)DRef = MultRMRt(pD->GetFDE(), RRef);

        bFirstRes = true;
}

Here is the call graph for this function:

virtual void HBeam::WorkSpaceDim ( integer *  piNumRows, integer *  piNumCols  ) const

inlinevirtual

Implements Elem.

Definition at line 213 of file hbeam.h.

                                                                            {
        *piNumRows = 12;
        *piNumCols = 12;
    };

Friends And Related Function Documentation

friend class AerodynamicBeam

friend

Definition at line 59 of file hbeam.h.

Member Data Documentation

Vec6 HBeam::Az

protected

Definition at line 100 of file hbeam.h.

Referenced by AfterPredict(), AssStiffnessVec(), HBeam(), and SetValue().

Vec6 HBeam::AzLoc

protected

Definition at line 102 of file hbeam.h.

Referenced by AfterPredict(), AssStiffnessVec(), dGetPrivData(), HBeam(), and Output().

Vec6 HBeam::AzRef

protected

Definition at line 101 of file hbeam.h.

Referenced by AfterPredict(), AssStiffnessMat(), HBeam(), and SetValue().

Beam::Type HBeam::BeamT

private

Definition at line 68 of file hbeam.h.

Referenced by SetBeamType().

bool HBeam::bFirstRes

protected

Definition at line 118 of file hbeam.h.

Referenced by AfterPredict(), and SetValue().

Vec6 HBeam::DefLoc

protected

Definition at line 103 of file hbeam.h.

Referenced by AfterConvergence(), AfterPredict(), AssStiffnessVec(), dGetPrivData(), HBeam(), and SetValue().

Vec6 HBeam::DefLocRef

protected

Definition at line 104 of file hbeam.h.

Referenced by AfterPredict(), HBeam(), and SetValue().

Mat6x6 HBeam::DRef

protected

Definition at line 93 of file hbeam.h.

Referenced by AfterPredict(), AssStiffnessMat(), and SetValue().

doublereal HBeam::dsdxi

protected

Definition at line 114 of file hbeam.h.

Referenced by DsDxi().

doublereal HBeam::dxids

protected

Definition at line 115 of file hbeam.h.

Referenced by AssStiffnessMat(), AssStiffnessVec(), and DsDxi().

Vec3 HBeam::f[NUMNODES]

protected

Definition at line 78 of file hbeam.h.

Referenced by AfterPredict(), AssStiffnessMat(), AssStiffnessVec(), DsDxi(), HBeam(), and Restart_().

Vec3 HBeam::fRef[NUMNODES]

protected

Definition at line 79 of file hbeam.h.

Vec3 HBeam::g

protected

Definition at line 107 of file hbeam.h.

Referenced by AfterPredict(), and HBeam().

Vec3 HBeam::L

protected

Definition at line 110 of file hbeam.h.

Referenced by AfterPredict(), AssStiffnessMat(), AssStiffnessVec(), DsDxi(), HBeam(), and SetValue().

Vec3 HBeam::L0

protected

Definition at line 109 of file hbeam.h.

Referenced by AfterPredict(), AssStiffnessVec(), DsDxi(), and HBeam().

Mat3x3 HBeam::Ldp[NUMNODES]

protected

Definition at line 131 of file hbeam.h.

Referenced by AssStiffnessMat().

Mat3x3 HBeam::LdP[NUMNODES]

protected

Definition at line 132 of file hbeam.h.

Referenced by AssStiffnessMat().

Vec3 HBeam::LRef

protected

Definition at line 111 of file hbeam.h.

Referenced by AfterPredict(), DsDxi(), and SetValue().

Vec3 HBeam::Omega

protected

Definition at line 96 of file hbeam.h.

Referenced by HBeam(), and Omega0().

Vec3 HBeam::OmegaRef

protected

Definition at line 97 of file hbeam.h.

Vec3 HBeam::p

protected

Definition at line 106 of file hbeam.h.

Referenced by AfterPredict(), AssStiffnessMat(), AssStiffnessVec(), DsDxi(), and HBeam().

ConstitutiveLaw6DOwner* HBeam::pD

protected

Definition at line 90 of file hbeam.h.

Referenced by AfterConvergence(), AfterPredict(), AssStiffnessMat(), AssStiffnessVec(), HBeam(), Restart_(), SetValue(), and ~HBeam().

Mat3x3 HBeam::pdp[NUMNODES]

protected

Definition at line 129 of file hbeam.h.

Referenced by AssStiffnessMat().

Mat3x3 HBeam::pdP[NUMNODES]

protected

Definition at line 130 of file hbeam.h.

Referenced by AssStiffnessMat().

const StructNode* HBeam::pNode[NUMNODES]

protected

Definition at line 75 of file hbeam.h.

Referenced by AfterPredict(), AssJac(), AssRes(), AssStiffnessMat(), AssStiffnessVec(), DsDxi(), GetConnectedNodes(), HBeam(), InitialAssJac(), InitialAssRes(), Omega0(), pGetNode(), and Restart_().

Mat3x3 HBeam::R

protected

Definition at line 86 of file hbeam.h.

Referenced by AfterPredict(), AssStiffnessMat(), AssStiffnessVec(), DsDxi(), Restart_(), and SetValue().

Mat3x3 HBeam::RdP[NUMNODES]

protected

Definition at line 133 of file hbeam.h.

Referenced by AssStiffnessMat().

Vec3 HBeam::Rho

protected

Definition at line 125 of file hbeam.h.

Referenced by AssStiffnessMat(), AssStiffnessVec(), and DsDxi().

Vec3 HBeam::Rho0

protected

Definition at line 126 of file hbeam.h.

Referenced by AssStiffnessVec(), and DsDxi().

Mat3x3 HBeam::RhodP[NUMNODES]

protected

Definition at line 134 of file hbeam.h.

Referenced by AssStiffnessMat().

Mat3x3 HBeam::Rn[NUMNODES]

protected

Definition at line 82 of file hbeam.h.

Referenced by AssStiffnessMat(), AssStiffnessVec(), DsDxi(), and HBeam().

Mat3x3 HBeam::RRef

protected

Definition at line 87 of file hbeam.h.

Referenced by AfterPredict(), DsDxi(), and SetValue().

doublereal HBeam::xi

protected

Definition at line 113 of file hbeam.h.

Referenced by AssStiffnessMat(), AssStiffnessVec(), and DsDxi().

---

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

• hbeam.h
• hbeam.cc