#include <constltp_impl.h>

Inheritance diagram for DoubleLinearViscoElasticConstitutiveLaw< Vec3, Mat3x3 >:

Collaboration diagram for DoubleLinearViscoElasticConstitutiveLaw< Vec3, Mat3x3 >:

Public Member Functions
	DoubleLinearViscoElasticConstitutiveLaw (const TplDriveCaller< Vec3 > *pDC, const Vec3 &PStress, doublereal dStiff, doublereal dUppLimStrain, doublereal dLowLimStrain, doublereal dSecondStiff, doublereal dStiffPrime, doublereal dSecondStiffPrime)

virtual	~DoubleLinearViscoElasticConstitutiveLaw (void)

ConstLawType::Type	GetConstLawType (void) const

virtual ConstitutiveLaw< Vec3, Mat3x3 > *	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
doublereal	dStiffness

doublereal	dUpperLimitStrain

doublereal	dLowerLimitStrain

doublereal	dSecondStiffness

doublereal	dStiffnessPrime

doublereal	dSecondStiffnessPrime

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 DoubleLinearViscoElasticConstitutiveLaw< Vec3, Mat3x3 >

Definition at line 1927 of file constltp_impl.h.

Constructor & Destructor Documentation

DoubleLinearViscoElasticConstitutiveLaw< Vec3, Mat3x3 >::DoubleLinearViscoElasticConstitutiveLaw	(	const TplDriveCaller< Vec3 > *	pDC,
		const Vec3 &	PStress,
		doublereal	dStiff,
		doublereal	dUppLimStrain,
		doublereal	dLowLimStrain,
		doublereal	dSecondStiff,
		doublereal	dStiffPrime,
		doublereal	dSecondStiffPrime
	)

inline

Definition at line 1938 of file constltp_impl.h.

References ConstitutiveLaw< T, Tder >::FDE, ConstitutiveLaw< T, Tder >::FDEPrime, Mat3x3DEye_Manip::Manipulate(), and Mat3x3DEye.

      : ElasticConstitutiveLaw3D(pDC, PStress),
      dStiffness(dStiff),
      dUpperLimitStrain(dUppLimStrain),
      dLowerLimitStrain(dLowLimStrain),
      dSecondStiffness(dSecondStiff),
      dStiffnessPrime(dStiffPrime),
      dSecondStiffnessPrime(dSecondStiffPrime)
    {
       Mat3x3DEye.Manipulate(FDE, dStiffness);
       Mat3x3DEye.Manipulate(FDEPrime, dStiffnessPrime);
    };

Here is the call graph for this function:

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

inlinevirtual

Definition at line 1958 of file constltp_impl.h.

References NO_OP.

                                                           {
       NO_OP;
    };

Member Function Documentation

ConstLawType::Type DoubleLinearViscoElasticConstitutiveLaw< Vec3, Mat3x3 >::GetConstLawType ( void ) const

inlinevirtual

Implements ConstitutiveLaw< Vec3, Mat3x3 >.

Definition at line 1962 of file constltp_impl.h.

References ConstLawType::VISCOELASTIC.

                                                      {
                 return ConstLawType::VISCOELASTIC;
         };

virtual ConstitutiveLaw<Vec3, Mat3x3>* DoubleLinearViscoElasticConstitutiveLaw< Vec3, Mat3x3 >::pCopy ( void ) const

inlinevirtual

Implements ConstitutiveLaw< Vec3, Mat3x3 >.

Definition at line 1966 of file constltp_impl.h.

References DoubleLinearViscoElasticConstitutiveLaw< T, Tder >::pCopy(), TplDriveOwner< T >::pGetDriveCaller(), ElasticConstitutiveLaw< T, Tder >::PreStress, and SAFENEWWITHCONSTRUCTOR.

                                                             {
       ConstitutiveLaw<Vec3, Mat3x3>* pCL = 0;
 
       typedef DoubleLinearViscoElasticConstitutiveLaw<Vec3, Mat3x3> cl;
       SAFENEWWITHCONSTRUCTOR(pCL,
                             cl,
                             cl(pGetDriveCaller()->pCopy(),
                                PreStress,
                                dStiffness,
                                dUpperLimitStrain,
                                dLowerLimitStrain,
                                dSecondStiffness,
                                dStiffnessPrime,
                                dSecondStiffnessPrime));
 
       return pCL;
    };

Here is the call graph for this function:

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

inlinevirtual

Reimplemented from ConstitutiveLaw< Vec3, Mat3x3 >.

Definition at line 1984 of file constltp_impl.h.

References ElasticConstitutiveLaw< T, Tder >::Restart_int().

                                                       {
       out << "double linear viscoelastic, "
         << dStiffness << ", "
         << dUpperLimitStrain << ", "
         << dLowerLimitStrain << ", "
         << dSecondStiffness << ", "
         << dStiffnessPrime << ", "
         << dSecondStiffnessPrime << ", ";
       return Restart_int(out);
    };

Here is the call graph for this function:

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

inlinevirtual

Implements ConstitutiveLaw< Vec3, Mat3x3 >.

Definition at line 1995 of file constltp_impl.h.

References Vec3::dGet(), ConstitutiveLaw< T, Tder >::Epsilon, ConstitutiveLaw< T, Tder >::EpsilonPrime, ConstitutiveLaw< T, Tder >::F, ConstitutiveLaw< T, Tder >::FDE, ConstitutiveLaw< T, Tder >::FDEPrime, TplDriveOwner< T >::Get(), and ElasticConstitutiveLaw< T, Tder >::PreStress.

                                                                       {
       Epsilon = Eps;
       EpsilonPrime = EpsPrime;
 
       Vec3 PreStrain = Get();
       Vec3 CurrStrain = Epsilon-PreStrain;
       doublereal dCurrStrain = CurrStrain.dGet(3);
 
       if (dCurrStrain <= dUpperLimitStrain && dCurrStrain >= dLowerLimitStrain) {
          FDE.Put(3, 3, dStiffness);
          FDEPrime.Put(3, 3, dStiffnessPrime);
          F = PreStress+CurrStrain*dStiffness+EpsilonPrime*dStiffnessPrime;
       } else {
          FDE.Put(3, 3, dSecondStiffness);
          FDEPrime.Put(3, 3, dSecondStiffnessPrime);
 
          if (dCurrStrain > dUpperLimitStrain) {
             F = PreStress
               +Vec3(CurrStrain(1)*dStiffness + EpsilonPrime(1)*dStiffnessPrime,
                     CurrStrain(2)*dStiffness + EpsilonPrime(2)*dStiffnessPrime,
                     dUpperLimitStrain*dStiffness
                         +(dCurrStrain-dUpperLimitStrain)*dSecondStiffness
                         +EpsilonPrime(3)*dSecondStiffnessPrime);
          } else /* if (dCurrStrain < dLowerLimitStrain) */ {
             F = PreStress
               +Vec3(CurrStrain.dGet(1)*dStiffness + EpsilonPrime(1)*dStiffnessPrime,
                     CurrStrain.dGet(2)*dStiffness + EpsilonPrime(2)*dStiffnessPrime,
                     dLowerLimitStrain*dStiffness
                         +(dCurrStrain-dLowerLimitStrain)*dSecondStiffness
                         +EpsilonPrime(3)*dSecondStiffnessPrime);
          }
       }
    };