ViscoElasticAxialJoint Class Reference

#include <vehj4.h>

Inheritance diagram for ViscoElasticAxialJoint:

Collaboration diagram for ViscoElasticAxialJoint:

Public Member Functions
	ViscoElasticAxialJoint (unsigned int uL, const DofOwner *pDO, const ConstitutiveLaw1D *pCL, const StructNode *pN1, const StructNode *pN2, const Mat3x3 &tilde_R1h, const Mat3x3 &tilde_R2h, flag fOut)
	~ViscoElasticAxialJoint (void)
virtual void	AfterConvergence (const VectorHandler &X, const VectorHandler &XP)
virtual ConstLawType::Type	GetConstLawType (void) const
virtual VariableSubMatrixHandler &	AssJac (VariableSubMatrixHandler &WorkMat, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual void	AssMats (VariableSubMatrixHandler &WorkMatA, VariableSubMatrixHandler &WorkMatB, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual SubVectorHandler &	AssRes (SubVectorHandler &WorkVec, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
SubVectorHandler &	AssRes (SubVectorHandler &WorkVec, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr, const VectorHandler &XPrimePrimeCurr, InverseDynamics::Order iOrder=InverseDynamics::INVERSE_DYNAMICS)
virtual void	InitialWorkSpaceDim (integer *piNumRows, integer *piNumCols) const
virtual VariableSubMatrixHandler &	InitialAssJac (VariableSubMatrixHandler &WorkMat, const VectorHandler &XCurr)
virtual SubVectorHandler &	InitialAssRes (SubVectorHandler &WorkVec, const VectorHandler &XCurr)
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 DofOrder::Order	GetDofType (unsigned int) 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 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	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	AfterConvergence (const VectorHandler &X, const VectorHandler &XP, const VectorHandler &XPP)
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 DeformableAxialJoint
	DeformableAxialJoint (unsigned int uL, const DofOwner *pDO, const ConstitutiveLaw1D *pCL, const StructNode *pN1, const StructNode *pN2, const Mat3x3 &tilde_R1h, const Mat3x3 &tilde_R2h, flag fOut)
virtual	~DeformableAxialJoint (void)
virtual Joint::Type	GetJointType (void) const
virtual std::ostream &	Restart (std::ostream &out) const
virtual void	Output (OutputHandler &OH) const
virtual void	AfterPredict (VectorHandler &X, VectorHandler &XP)
void	SetValue (DataManager *pDM, VectorHandler &X, VectorHandler &XP, SimulationEntity::Hints *ph=0)
virtual void	SetInitialValue (VectorHandler &)
virtual bool	bInverseDynamics (void) const
virtual Hint *	ParseHint (DataManager *pDM, const char *s) const
virtual unsigned int	iGetNumDof (void) const
virtual void	WorkSpaceDim (integer *piNumRows, integer *piNumCols) const
virtual unsigned int	iGetInitialNumDof (void) const
virtual void	GetConnectedNodes (std::vector< const Node * > &connectedNodes) const
virtual unsigned int	iGetNumPrivData (void) const
virtual unsigned int	iGetPrivDataIdx (const char *s) const
virtual doublereal	dGetPrivData (unsigned int i) const
Public Member Functions inherited from Joint
	Joint (unsigned int uL, const DofOwner *pD, flag fOut)
virtual	~Joint (void)
virtual Elem::Type	GetElemType (void) const
std::ostream &	Output (std::ostream &out, const char *sJointName, unsigned int uLabel, const Vec3 &FLocal, const Vec3 &MLocal, const Vec3 &FGlobal, const Vec3 &MGlobal) const
virtual void	Update (const VectorHandler &XCurr, InverseDynamics::Order iOrder=InverseDynamics::INVERSE_DYNAMICS)
bool	bIsPrescribedMotion (void) const
bool	bIsTorque (void) const
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 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
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)
Public Member Functions inherited from ConstitutiveLawOwner< T, Tder >
	ConstitutiveLawOwner (const ConstitutiveLaw< T, Tder > *pCL)
virtual	~ConstitutiveLawOwner (void)
ConstitutiveLaw< T, Tder > *	pGetConstLaw (void) const
void	Update (const T &Eps, const T &EpsPrime=mb_zero< T >())
void	AfterConvergence (const T &Eps, const T &EpsPrime=mb_zero< T >())
const T &	GetF (void) const
const Tder &	GetFDE (void) const
const Tder &	GetFDEPrime (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 i) const
virtual std::ostream &	OutputAppend (std::ostream &out) const

Protected Member Functions
virtual void	AssMats (FullSubMatrixHandler &WMA, FullSubMatrixHandler &WMB, doublereal dCoef)
virtual void	AfterPredict (void)
virtual void	AssVec (SubVectorHandler &WorkVec)
Protected Member Functions inherited from DeformableAxialJoint
virtual void	AssMatM (FullSubMatrixHandler &WMA, doublereal dCoef)
void	AssMatMDE (FullSubMatrixHandler &WMA, doublereal dCoef)
virtual void	AssMatMDEPrime (FullSubMatrixHandler &WMA, FullSubMatrixHandler &WMB, doublereal dCoef)
Protected Member Functions inherited from Joint
virtual void	OutputPrepare_int (const std::string &type, OutputHandler &OH, std::string &name)
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
doublereal	dThetaRef
doublereal	dThetaCurr
doublereal	dOmega
Protected Attributes inherited from WithLabel
unsigned int	uLabel
std::string	sName
Protected Attributes inherited from ToBeOutput
flag	fOutput
Protected Attributes inherited from DeformableAxialJoint
const StructNode *	pNode1
const StructNode *	pNode2
Mat3x3	tilde_R1h
Mat3x3	tilde_R2h
bool	bFirstRes
doublereal	dTol
doublereal	dM
doublereal	dMDE
doublereal	dMDEPrime
Vec3	M
Mat3x3	MDE
Mat3x3	MDEPrime
Protected Attributes inherited from GravityOwner
Gravity *	pGravity
Protected Attributes inherited from ConstitutiveLawOwner< T, Tder >
ConstitutiveLaw< T, Tder > *	pConstLaw

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 }
Public Types inherited from Joint
enum	Type {   UNKNOWN = -1, DISTANCE = 0, DISTANCEWITHOFFSET, CLAMP,   SPHERICALHINGE, PIN, UNIVERSALHINGE, UNIVERSALROTATION,   UNIVERSALPIN, PLANEHINGE, PLANEROTATION, PLANEPIN,   AXIALROTATION, PLANEDISP, PLANEDISPPIN, INPLANE,   INPLANECONTACT, J_INLINE, ROD, RODBEZIER,   DEFORMABLEHINGE, DEFORMABLEDISPJOINT, DEFORMABLEJOINT, DEFORMABLEAXIALJOINT,   VISCOUSBODY, LINEARVELOCITY, ANGULARVELOCITY, LINEARACCELERATION,   ANGULARACCELERATION, PRISMATIC, DRIVEHINGE, DRIVEDISP,   DRIVEDISPPIN, IMPOSEDORIENTATION, IMPOSEDDISP, IMPOSEDDISPPIN,   IMPOSEDKINEMATICS, BEAMSLIDER, BRAKE, GIMBAL,   POINT_SURFACE_CONTACT, TOTALJOINT, TOTALPINJOINT, TOTALEQUATION,   TOTALREACTION, MODAL, SCREWJOINT, LASTJOINTTYPE }

Detailed Description

Definition at line 388 of file vehj4.h.

Constructor & Destructor Documentation

ViscoElasticAxialJoint::ViscoElasticAxialJoint	(	unsigned int	uL,
		const DofOwner *	pDO,
		const ConstitutiveLaw1D *	pCL,
		const StructNode *	pN1,
		const StructNode *	pN2,
		const Mat3x3 &	tilde_R1h,
		const Mat3x3 &	tilde_R2h,
		flag	fOut
	)

Definition at line 934 of file vehj4.cc.

References AfterPredict().

: Elem(uL, fOut),
DeformableAxialJoint(uL, pDO, pCL, pN1, pN2, tilde_R1h, tilde_R2h, fOut),
dThetaRef(0.)
{
        // force update of MDE/MDEPrime as needed
        AfterPredict();
}

Here is the call graph for this function:

ViscoElasticAxialJoint::~ViscoElasticAxialJoint

(

void

)

Definition at line 950 of file vehj4.cc.

References NO_OP.

{
        NO_OP;
}

Member Function Documentation

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

virtual

Reimplemented from SimulationEntity.

Definition at line 956 of file vehj4.cc.

References ConstitutiveLawOwner< T, Tder >::AfterConvergence(), dOmega, and dThetaCurr.

{
        ConstitutiveLaw1DOwner::AfterConvergence(dThetaCurr, dOmega);
}

Here is the call graph for this function:

void ViscoElasticAxialJoint::AfterPredict ( void  )

protectedvirtual

Implements DeformableAxialJoint.

Definition at line 963 of file vehj4.cc.

References DeformableAxialJoint::bFirstRes, dOmega, dThetaCurr, dThetaRef, ConstitutiveLawOwner< T, Tder >::GetFDE(), ConstitutiveLawOwner< T, Tder >::GetFDEPrime(), StructNode::GetRRef(), StructNode::GetWRef(), DeformableAxialJoint::MDE, DeformableAxialJoint::MDEPrime, DeformableAxialJoint::pNode1, DeformableAxialJoint::pNode2, Mat3x3::Tens(), DeformableAxialJoint::tilde_R1h, DeformableAxialJoint::tilde_R2h, ConstitutiveLawOwner< T, Tder >::Update(), and RotManip::VecRot().

Referenced by ViscoElasticAxialJoint().

{
        /* Computes strains, updates constitutive law and generates
         * MDE and MDEPrime */

        Mat3x3 R1h(pNode1->GetRRef()*tilde_R1h);
        Mat3x3 R2h(pNode2->GetRRef()*tilde_R2h);

        /* Current strain in material reference frame (node 1) */
        dThetaCurr = dThetaRef = RotManip::VecRot(R1h.MulTM(R2h))(3);

        /* Relative angular velocity */
        Vec3 e1z(R1h.GetVec(3));
        dOmega = e1z.Dot(pNode2->GetWRef() - pNode1->GetWRef());

        /* Updates constitutive law */
        ConstitutiveLaw1DOwner::Update(dThetaRef, dOmega);

        /* don't repeat the above operations during AssRes */
        bFirstRes = true;

        /* Tangent matrices are updated and projected in the global
         * reference frame; they won't change during the solution
         * of the current time step, according to the updated-updated
         * approach */
        MDE = e1z.Tens(e1z*ConstitutiveLaw1DOwner::GetFDE());
        MDEPrime = e1z.Tens(e1z*GetFDEPrime());
}

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 994 of file vehj4.cc.

References AssMats(), StructDispNode::iGetFirstMomentumIndex(), StructDispNode::iGetFirstPositionIndex(), DeformableAxialJoint::pNode1, DeformableAxialJoint::pNode2, FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), FullSubMatrixHandler::ResizeReset(), VariableSubMatrixHandler::SetFull(), and Elem::WorkSpaceDim().

{
        FullSubMatrixHandler& WM = WorkMat.SetFull();

        /* Dimensiona e resetta la matrice di lavoro */
        integer iNumRows = 0;
        integer iNumCols = 0;
        WorkSpaceDim(&iNumRows, &iNumCols);
        WM.ResizeReset(iNumRows, iNumCols);

        /* Recupera gli indici */
        integer iNode1FirstPosIndex = pNode1->iGetFirstPositionIndex() + 3;
        integer iNode1FirstMomIndex = pNode1->iGetFirstMomentumIndex() + 3;
        integer iNode2FirstPosIndex = pNode2->iGetFirstPositionIndex() + 3;
        integer iNode2FirstMomIndex = pNode2->iGetFirstMomentumIndex() + 3;

        /* Setta gli indici della matrice */
        for (int iCnt = 1; iCnt <= 3; iCnt++) {
                WM.PutRowIndex(iCnt, iNode1FirstMomIndex + iCnt);
                WM.PutColIndex(iCnt, iNode1FirstPosIndex + iCnt);
                WM.PutRowIndex(3 + iCnt, iNode2FirstMomIndex + iCnt);
                WM.PutColIndex(3 + iCnt, iNode2FirstPosIndex + iCnt);
        }

        AssMats(WM, WM, dCoef);

        return WorkMat;
}

Here is the call graph for this function:

void ViscoElasticAxialJoint::AssMats ( FullSubMatrixHandler &  WMA, FullSubMatrixHandler &  WMB, doublereal  dCoef  )

protectedvirtual

Definition at line 1067 of file vehj4.cc.

References DeformableAxialJoint::AssMatM(), DeformableAxialJoint::AssMatMDE(), and DeformableAxialJoint::AssMatMDEPrime().

Referenced by AssJac(), and AssMats().

{
        AssMatM(WMA, dCoef);
        AssMatMDE(WMA, dCoef);
        AssMatMDEPrime(WMA, WMB, dCoef);
}

Here is the call graph for this function:

void ViscoElasticAxialJoint::AssMats ( VariableSubMatrixHandler &  WorkMatA, VariableSubMatrixHandler &  WorkMatB, const VectorHandler &  XCurr, const VectorHandler &  XPrimeCurr  )

virtual

Reimplemented from Elem.

Definition at line 1028 of file vehj4.cc.

References AssMats(), StructDispNode::iGetFirstMomentumIndex(), StructDispNode::iGetFirstPositionIndex(), DeformableAxialJoint::pNode1, DeformableAxialJoint::pNode2, FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), FullSubMatrixHandler::ResizeReset(), VariableSubMatrixHandler::SetFull(), and Elem::WorkSpaceDim().

{
        FullSubMatrixHandler& WMA = WorkMatA.SetFull();
        FullSubMatrixHandler& WMB = WorkMatB.SetFull();

        /* Dimensiona e resetta la matrice di lavoro */
        integer iNumRows = 0;
        integer iNumCols = 0;
        WorkSpaceDim(&iNumRows, &iNumCols);
        WMA.ResizeReset(iNumRows, iNumCols);
        WMB.ResizeReset(iNumRows, iNumCols);

        /* Recupera gli indici */
        integer iNode1FirstPosIndex = pNode1->iGetFirstPositionIndex() + 3;
        integer iNode1FirstMomIndex = pNode1->iGetFirstMomentumIndex() + 3;
        integer iNode2FirstPosIndex = pNode2->iGetFirstPositionIndex() + 3;
        integer iNode2FirstMomIndex = pNode2->iGetFirstMomentumIndex() + 3;

        /* Setta gli indici della matrice */
        for (int iCnt = 1; iCnt <= 3; iCnt++) {
                WMA.PutRowIndex(iCnt, iNode1FirstMomIndex + iCnt);
                WMA.PutColIndex(iCnt, iNode1FirstPosIndex + iCnt);
                WMA.PutRowIndex(3 + iCnt, iNode2FirstMomIndex + iCnt);
                WMA.PutColIndex(3 + iCnt, iNode2FirstPosIndex + iCnt);

                WMB.PutRowIndex(iCnt, iNode1FirstMomIndex + iCnt);
                WMB.PutColIndex(iCnt, iNode1FirstPosIndex + iCnt);
                WMB.PutRowIndex(3 + iCnt, iNode2FirstMomIndex + iCnt);
                WMB.PutColIndex(3 + iCnt, iNode2FirstPosIndex + iCnt);
        }

        AssMats(WMA, WMB, 1.);
}

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 1078 of file vehj4.cc.

References AssVec(), StructDispNode::iGetFirstMomentumIndex(), DeformableAxialJoint::pNode1, DeformableAxialJoint::pNode2, SubVectorHandler::PutRowIndex(), VectorHandler::ResizeReset(), and Elem::WorkSpaceDim().

{
        /* Dimensiona e resetta la matrice di lavoro */
        integer iNumRows = 0;
        integer iNumCols = 0;
        WorkSpaceDim(&iNumRows, &iNumCols);
        WorkVec.ResizeReset(iNumRows);

        /* Recupera gli indici */
        integer iNode1FirstMomIndex = pNode1->iGetFirstMomentumIndex() + 3;
        integer iNode2FirstMomIndex = pNode2->iGetFirstMomentumIndex() + 3;

        /* Setta gli indici della matrice */
        for (int iCnt = 1; iCnt <= 3; iCnt++) {
                WorkVec.PutRowIndex(iCnt, iNode1FirstMomIndex + iCnt);
                WorkVec.PutRowIndex(3 + iCnt, iNode2FirstMomIndex + iCnt);
        }

        AssVec(WorkVec);

        return WorkVec;
}

Here is the call graph for this function:

SubVectorHandler & ViscoElasticAxialJoint::AssRes ( SubVectorHandler &  WorkVec, const VectorHandler &  XCurr, const VectorHandler &  XPrimeCurr, const VectorHandler &  XPrimePrimeCurr, InverseDynamics::Order  iOrder = InverseDynamics::INVERSE_DYNAMICS  )

virtual

Reimplemented from Elem.

Definition at line 1106 of file vehj4.cc.

References ASSERT, AssVec(), DeformableAxialJoint::bFirstRes, StructDispNode::iGetFirstPositionIndex(), InverseDynamics::INVERSE_DYNAMICS, DeformableAxialJoint::pNode1, DeformableAxialJoint::pNode2, SubVectorHandler::PutRowIndex(), VectorHandler::ResizeReset(), and Elem::WorkSpaceDim().

{
        ASSERT(iOrder == InverseDynamics::INVERSE_DYNAMICS);

        bFirstRes = false;

        /* Dimensiona e resetta la matrice di lavoro */
        integer iNumRows = 0;
        integer iNumCols = 0;
        WorkSpaceDim(&iNumRows, &iNumCols);
        WorkVec.ResizeReset(iNumRows);

        /* Recupera gli indici */
        integer iNode1FirstMomIndex = pNode1->iGetFirstPositionIndex() + 3;
        integer iNode2FirstMomIndex = pNode2->iGetFirstPositionIndex() + 3;

        /* Setta gli indici della matrice */
        for (int iCnt = 1; iCnt <= 3; iCnt++) {
                WorkVec.PutRowIndex(iCnt, iNode1FirstMomIndex + iCnt);
                WorkVec.PutRowIndex(3 + iCnt, iNode2FirstMomIndex + iCnt);
        }

        AssVec(WorkVec);

        return WorkVec;
}

Here is the call graph for this function:

void ViscoElasticAxialJoint::AssVec ( SubVectorHandler &  WorkVec )

protectedvirtual

Definition at line 1139 of file vehj4.cc.

References VectorHandler::Add(), DeformableAxialJoint::bFirstRes, dOmega, dThetaCurr, ConstitutiveLawOwner< T, Tder >::GetF(), StructNode::GetRCurr(), StructNode::GetWCurr(), DeformableAxialJoint::M, DeformableAxialJoint::pNode1, DeformableAxialJoint::pNode2, VectorHandler::Sub(), DeformableAxialJoint::tilde_R1h, DeformableAxialJoint::tilde_R2h, ConstitutiveLawOwner< T, Tder >::Update(), and RotManip::VecRot().

Referenced by AssRes().

{
        Mat3x3 R1h(pNode1->GetRCurr()*tilde_R1h);
        Vec3 e1z(R1h.GetVec(3));

        if (bFirstRes) {
                bFirstRes = false;

        } else {
                Mat3x3 R2h(pNode2->GetRCurr()*tilde_R2h);

                /* orientazione intermedia */
                dThetaCurr = RotManip::VecRot(R1h.MulTM(R2h))(3);

                /* velocita' relativa nel riferimento intermedio */
                dOmega = e1z.Dot(pNode2->GetWCurr() - pNode1->GetWCurr());

                /* aggiorna il legame costitutivo */
                ConstitutiveLaw1DOwner::Update(dThetaCurr, dOmega);
        }

        M = e1z*ConstitutiveLaw1DOwner::GetF();

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

Here is the call graph for this function:

virtual ConstLawType::Type ViscoElasticAxialJoint::GetConstLawType ( void  ) const

inlinevirtual

Implements DeformableAxialJoint.

Definition at line 419 of file vehj4.h.

References ConstLawType::VISCOELASTIC.

                                                             {
                return ConstLawType::VISCOELASTIC;
        };

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

virtual

Implements SubjectToInitialAssembly.

Definition at line 1168 of file vehj4.cc.

References FullSubMatrixHandler::Add(), ConstitutiveLawOwner< T, Tder >::GetF(), ConstitutiveLawOwner< T, Tder >::GetFDE(), ConstitutiveLawOwner< T, Tder >::GetFDEPrime(), StructNode::GetRCurr(), Mat3x3::GetVec(), StructNode::GetWCurr(), StructDispNode::iGetFirstPositionIndex(), InitialWorkSpaceDim(), MatCross, DeformableAxialJoint::MDE, DeformableAxialJoint::MDEPrime, DeformableAxialJoint::pNode1, DeformableAxialJoint::pNode2, FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), FullSubMatrixHandler::ResizeReset(), VariableSubMatrixHandler::SetFull(), FullSubMatrixHandler::Sub(), Mat3x3::Tens(), and DeformableAxialJoint::tilde_R1h.

{
        FullSubMatrixHandler& WM = WorkMat.SetFull();

        /* Dimensiona e resetta la matrice di lavoro */
        integer iNumRows = 0;
        integer iNumCols = 0;
        InitialWorkSpaceDim(&iNumRows, &iNumCols);
        WM.ResizeReset(iNumRows, iNumCols);

        /* Recupera gli indici */
        integer iNode1FirstPosIndex = pNode1->iGetFirstPositionIndex() + 3;
        integer iNode1FirstVelIndex = iNode1FirstPosIndex + 6;
        integer iNode2FirstPosIndex = pNode2->iGetFirstPositionIndex() + 3;
        integer iNode2FirstVelIndex = iNode2FirstPosIndex + 6;

        /* Setta gli indici della matrice */
        for (int iCnt = 1; iCnt <= 3; iCnt++) {
                WM.PutRowIndex(iCnt, iNode1FirstPosIndex + iCnt);
                WM.PutColIndex(iCnt, iNode1FirstPosIndex + iCnt);
                WM.PutColIndex(3 + iCnt, iNode1FirstVelIndex + iCnt);

                WM.PutRowIndex(3 + iCnt, iNode2FirstPosIndex + iCnt);
                WM.PutColIndex(6 + iCnt, iNode2FirstPosIndex + iCnt);
                WM.PutColIndex(9 + iCnt, iNode2FirstVelIndex + iCnt);
        }

        Vec3 e1z(pNode1->GetRCurr()*tilde_R1h.GetVec(3));
        const Vec3& W2(pNode2->GetWCurr());

        MDE = e1z.Tens(e1z*ConstitutiveLaw1DOwner::GetFDE());
        MDEPrime = e1z.Tens(e1z*ConstitutiveLaw1DOwner::GetFDEPrime());

        Mat3x3 Tmp(MDE + MDEPrime*Mat3x3(MatCross, W2));
        WM.Add(4, 7, Tmp);
        WM.Sub(1, 7, Tmp);

        Tmp += Mat3x3(MatCross, e1z*ConstitutiveLaw1DOwner::GetF());
        WM.Add(1, 1, Tmp);
        WM.Sub(4, 1, Tmp);

        WM.Add(1, 4, MDEPrime);
        WM.Add(4, 10, MDEPrime);

        WM.Sub(1, 10, MDEPrime);
        WM.Sub(4, 4, MDEPrime);

        return WorkMat;
}

Here is the call graph for this function:

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

virtual

Implements SubjectToInitialAssembly.

Definition at line 1221 of file vehj4.cc.

References VectorHandler::Add(), DeformableAxialJoint::bFirstRes, dOmega, dThetaCurr, ConstitutiveLawOwner< T, Tder >::GetF(), StructNode::GetRCurr(), StructNode::GetWCurr(), StructDispNode::iGetFirstPositionIndex(), InitialWorkSpaceDim(), DeformableAxialJoint::M, DeformableAxialJoint::pNode1, DeformableAxialJoint::pNode2, SubVectorHandler::PutRowIndex(), VectorHandler::ResizeReset(), VectorHandler::Sub(), DeformableAxialJoint::tilde_R1h, DeformableAxialJoint::tilde_R2h, ConstitutiveLawOwner< T, Tder >::Update(), and RotManip::VecRot().

{
        /* Dimensiona e resetta la matrice di lavoro */
        integer iNumRows = 0;
        integer iNumCols = 0;
        InitialWorkSpaceDim(&iNumRows, &iNumCols);
        WorkVec.ResizeReset(iNumRows);

        /* Recupera gli indici */
        integer iNode1FirstPosIndex = pNode1->iGetFirstPositionIndex() + 3;
        integer iNode2FirstPosIndex = pNode2->iGetFirstPositionIndex() + 3;

        /* Setta gli indici della matrice */
        for (int iCnt = 1; iCnt <= 3; iCnt++) {
                WorkVec.PutRowIndex(iCnt, iNode1FirstPosIndex + iCnt);
                WorkVec.PutRowIndex(3 + iCnt, iNode2FirstPosIndex + iCnt);
        }

        Mat3x3 R1h(pNode1->GetRCurr()*tilde_R1h);
        Vec3 e1z(R1h.GetVec(3));

        if (bFirstRes) {
                bFirstRes = false;

        } else {
                Mat3x3 R2h(pNode2->GetRCurr()*tilde_R2h);

                dThetaCurr = RotManip::VecRot(R1h.MulTM(R2h))(3);
                dOmega = e1z.Dot(pNode2->GetWCurr() - pNode1->GetWCurr());

                ConstitutiveLaw1DOwner::Update(dThetaCurr, dOmega);
        }

        M = e1z*ConstitutiveLaw1DOwner::GetF();

        WorkVec.Add(1, M);
        WorkVec.Sub(4, M);

        return WorkVec;
}

Here is the call graph for this function:

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

inlinevirtual

Implements SubjectToInitialAssembly.

Definition at line 466 of file vehj4.h.

Referenced by InitialAssJac(), and InitialAssRes().

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

Member Data Documentation

doublereal ViscoElasticAxialJoint::dOmega

protected

Definition at line 394 of file vehj4.h.

Referenced by AfterConvergence(), AfterPredict(), AssVec(), and InitialAssRes().

doublereal ViscoElasticAxialJoint::dThetaCurr

protected

Definition at line 392 of file vehj4.h.

Referenced by AfterConvergence(), AfterPredict(), AssVec(), and InitialAssRes().

doublereal ViscoElasticAxialJoint::dThetaRef

protected

Definition at line 391 of file vehj4.h.

Referenced by AfterPredict().

---

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

• vehj4.h
• vehj4.cc