#include <constltp_impl.h>

Inheritance diagram for CubicViscoElasticGenericConstitutiveLaw< Vec3, Mat3x3 >:

Collaboration diagram for CubicViscoElasticGenericConstitutiveLaw< Vec3, Mat3x3 >:

Public Member Functions
	CubicViscoElasticGenericConstitutiveLaw (const TplDriveCaller< Vec3 > *pDC, const Vec3 &PStress, const Vec3 &Stiff1, const Vec3 &Stiff2, const Vec3 &Stiff3, const Mat3x3 &StiffPrime)

virtual	~CubicViscoElasticGenericConstitutiveLaw (void)

virtual ConstitutiveLaw3D *	pCopy (void) const

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

virtual void	Update (const Vec3 &Eps, const Vec3 &EpsPrime=Zero3)

Public Member Functions inherited from ElasticConstitutiveLaw< Vec3, Mat3x3 >
	ElasticConstitutiveLaw (const TplDriveCaller< Vec3 > *pDC, const Vec3 &PStress)

virtual	~ElasticConstitutiveLaw (void)

ConstLawType::Type	GetConstLawType (void) const

void	SetValue (DataManager pDM, VectorHandler &X, VectorHandler &XP, SimulationEntity::Hints ph=0)

virtual Hint *	ParseHint (DataManager pDM, const char s) const

Public Member Functions inherited from ConstitutiveLaw< Vec3, Mat3x3 >
	ConstitutiveLaw (void)

virtual	~ConstitutiveLaw (void)

virtual void	AfterConvergence (const Vec3 &Eps, const Vec3 &EpsPrime=mb_zero< Vec3 >())

virtual const Vec3 &	GetEpsilon (void) const

virtual const Vec3 &	GetEpsilonPrime (void) const

virtual const Vec3 &	GetF (void) const

virtual const Mat3x3 &	GetFDE (void) const

virtual const Mat3x3 &	GetFDEPrime (void) const

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 i) 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	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 TplDriveOwner< Vec3 >
	TplDriveOwner (const TplDriveCaller< Vec3 > *pDC=0)

virtual	~TplDriveOwner (void)

void	Set (const TplDriveCaller< Vec3 > *pDC)

TplDriveCaller< Vec3 > *	pGetDriveCaller (void) const

Vec3	Get (const doublereal &dVar) const

Vec3	Get (void) const

virtual bool	bIsDifferentiable (void) const

virtual Vec3	GetP (void) const

Private Attributes
Vec3	Stiff1

Vec3	Stiff2

Vec3	Stiff3

Additional Inherited Members
Public Types inherited from ConstitutiveLaw< Vec3, Mat3x3 >
typedef ConstitutiveLaw< Vec3, Mat3x3 >::ErrNotAvailable	Err

Public Types inherited from SimulationEntity
typedef std::vector< Hint * >	Hints

Protected Member Functions inherited from ElasticConstitutiveLaw< Vec3, Mat3x3 >
virtual std::ostream &	Restart_int (std::ostream &out) const

Protected Attributes inherited from ElasticConstitutiveLaw< Vec3, Mat3x3 >
Vec3	PreStress

Protected Attributes inherited from ConstitutiveLaw< Vec3, Mat3x3 >
Vec3	Epsilon

Vec3	EpsilonPrime

Vec3	F

Mat3x3	FDE

Mat3x3	FDEPrime

Protected Attributes inherited from WithLabel
unsigned int	uLabel

std::string	sName

Protected Attributes inherited from TplDriveOwner< Vec3 >
TplDriveCaller< Vec3 > *	pTplDriveCaller

Detailed Description

template<>
class CubicViscoElasticGenericConstitutiveLaw< Vec3, Mat3x3 >

Definition at line 1703 of file constltp_impl.h.

Constructor & Destructor Documentation

CubicViscoElasticGenericConstitutiveLaw< Vec3, Mat3x3 >::CubicViscoElasticGenericConstitutiveLaw	(	const TplDriveCaller< Vec3 > *	pDC,
		const Vec3 &	PStress,
		const Vec3 &	Stiff1,
		const Vec3 &	Stiff2,
		const Vec3 &	Stiff3,
		const Mat3x3 &	StiffPrime
	)

inline

Definition at line 1711 of file constltp_impl.h.

References ConstitutiveLaw< Vec3, Mat3x3 >::FDEPrime.

         : ElasticConstitutiveLaw3D(pDC, PStress),
                 Stiff1(Stiff1), Stiff2(Stiff2), Stiff3(Stiff3)
         {
                 ConstitutiveLaw3D::FDEPrime = StiffPrime;
         };

virtual CubicViscoElasticGenericConstitutiveLaw< Vec3, Mat3x3 >::~CubicViscoElasticGenericConstitutiveLaw ( void )

inlinevirtual

Definition at line 1721 of file constltp_impl.h.

References NO_OP.

                                                                {
                 NO_OP;
         };

Member Function Documentation

virtual ConstitutiveLaw3D* CubicViscoElasticGenericConstitutiveLaw< Vec3, Mat3x3 >::pCopy ( void ) const

inlinevirtual

Implements ConstitutiveLaw< Vec3, Mat3x3 >.

Definition at line 1725 of file constltp_impl.h.

References ConstitutiveLaw< Vec3, Mat3x3 >::FDEPrime, CubicViscoElasticGenericConstitutiveLaw< T, Tder >::pCopy(), TplDriveOwner< Vec3 >::pGetDriveCaller(), ElasticConstitutiveLaw< Vec3, Mat3x3 >::PreStress, and SAFENEWWITHCONSTRUCTOR.

                                                      {
                 ConstitutiveLaw3D* pCL = 0;
 
                 typedef CubicViscoElasticGenericConstitutiveLaw<Vec3, Mat3x3> cl;
                 SAFENEWWITHCONSTRUCTOR(pCL,
                                 cl,
                                 cl(ElasticConstitutiveLaw3D::pGetDriveCaller()->pCopy(),
                                         ElasticConstitutiveLaw3D::PreStress,
                                         Stiff1, Stiff2, Stiff3, ConstitutiveLaw3D::FDEPrime));
                 return pCL;
         };

Here is the call graph for this function:

virtual std::ostream& CubicViscoElasticGenericConstitutiveLaw< Vec3, Mat3x3 >::Restart ( std::ostream & out ) const

inlinevirtual

Reimplemented from ConstitutiveLaw< Vec3, Mat3x3 >.

Definition at line 1737 of file constltp_impl.h.

References ConstitutiveLaw< Vec3, Mat3x3 >::FDEPrime, ElasticConstitutiveLaw< Vec3, Mat3x3 >::Restart_int(), and Write().

                                                            {
                 out << "cubic elastic generic, ",
                         Write(out, Stiff1, ", ") << ", ",
                         Write(out, Stiff2, ", ") << ", ",
                         Write(out, Stiff3, ", ") << ", ",
                         Write(out, ConstitutiveLaw3D::FDEPrime, ", ") << ", ";
                 return ElasticConstitutiveLaw3D::Restart_int(out);
         };

Here is the call graph for this function:

virtual void CubicViscoElasticGenericConstitutiveLaw< Vec3, Mat3x3 >::Update	(	const Vec3 &	Eps,
		const Vec3 &	EpsPrime = `Zero3`
	)

inlinevirtual

Implements ConstitutiveLaw< Vec3, Mat3x3 >.

Definition at line 1746 of file constltp_impl.h.

References ConstitutiveLaw< Vec3, Mat3x3 >::Epsilon, ConstitutiveLaw< Vec3, Mat3x3 >::EpsilonPrime, ConstitutiveLaw< Vec3, Mat3x3 >::F, grad::fabs(), ConstitutiveLaw< Vec3, Mat3x3 >::FDE, ConstitutiveLaw< Vec3, Mat3x3 >::FDEPrime, TplDriveOwner< Vec3 >::Get(), ElasticConstitutiveLaw< Vec3, Mat3x3 >::PreStress, and Mat3x3::Put().

                                                                            {
                 ConstitutiveLaw3D::Epsilon = Eps;
                 ConstitutiveLaw3D::EpsilonPrime = EpsPrime;
                 Vec3 v1 = Eps - ElasticConstitutiveLaw3D::Get();
                 ConstitutiveLaw3D::F = ElasticConstitutiveLaw3D::PreStress + ConstitutiveLaw3D::FDEPrime*EpsPrime;
 
 #if defined(MBDYN_X_WORKAROUND_GCC_3_2) || defined(MBDYN_X_WORKAROUND_GCC_3_3)
                 Vec3 FTmp;
 #endif // MBDYN_X_WORKAROUND_GCC_3_2 || MBDYN_X_WORKAROUND_GCC_3_3
 
                 for (int iCnt = 1; iCnt <= 3; iCnt++) {
                         doublereal e1 = v1(iCnt);
                         doublereal f1 = fabs(e1);
                         doublereal e2 = e1*e1;
                         doublereal f2 = f1*e1;
                         doublereal e3 = e2*e1;
                         
 #if defined(MBDYN_X_WORKAROUND_GCC_3_2) || defined(MBDYN_X_WORKAROUND_GCC_3_3)
                         ConstitutiveLaw3D::FDE.Put(iCnt, iCnt,
                                 Stiff1(iCnt)
                                 + 2.*Stiff2(iCnt)*f1
                                 + 3.*Stiff3(iCnt)*e2);
                         FTmp(iCnt) = Stiff1(iCnt)*e1
                                 + Stiff2(iCnt)*f2
                                 + Stiff3(iCnt)*e3;
 #else // ! MBDYN_X_WORKAROUND_GCC_3_2 && ! MBDYN_X_WORKAROUND_GCC_3_3
                         ConstitutiveLaw3D::FDE(iCnt, iCnt) = Stiff1(iCnt)
                                 + 2.*Stiff2(iCnt)*f1 + 3.*Stiff3(iCnt)*e2;
                         ConstitutiveLaw3D::F(iCnt) += Stiff1(iCnt)*e1
                                 + Stiff2(iCnt)*f2 + Stiff3(iCnt)*e3;
 #endif // ! MBDYN_X_WORKAROUND_GCC_3_2 && ! MBDYN_X_WORKAROUND_GCC_3_3
                 }
 
 #if defined(MBDYN_X_WORKAROUND_GCC_3_2) || defined(MBDYN_X_WORKAROUND_GCC_3_3)
                 ConstitutiveLaw3D::F += FTmp;
 #endif // ! MBDYN_X_WORKAROUND_GCC_3_2 && ! MBDYN_X_WORKAROUND_GCC_3_3
         };