Inheritance diagram for NLPViscoElasticConstitutiveLaw< T, Tder >:

Collaboration diagram for NLPViscoElasticConstitutiveLaw< T, Tder >:

Public Member Functions
	NLPViscoElasticConstitutiveLaw (const TplDriveCaller< T > pDC, const T &PStress, const ConstLawType::Type &cltype, const Tder &fde0, const std::vector< const DifferentiableScalarFunction > &fdesc, const Tder &fdeprime0, const std::vector< const DifferentiableScalarFunction * > &fdeprimesc)

virtual	~NLPViscoElasticConstitutiveLaw (void)

ConstLawType::Type	GetConstLawType (void) const

virtual ConstitutiveLaw< T, Tder > *	pCopy (void) const

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

virtual void	Update (const T &Eps, const T &EpsPrime=0.)

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

virtual	~ElasticConstitutiveLaw (void)

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< T, Tder >
	ConstitutiveLaw (void)

virtual	~ConstitutiveLaw (void)

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

virtual const T &	GetEpsilon (void) const

virtual const T &	GetEpsilonPrime (void) const

virtual const T &	GetF (void) const

virtual const Tder &	GetFDE (void) const

virtual const Tder &	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< T >
	TplDriveOwner (const TplDriveCaller< T > *pDC=0)

virtual	~TplDriveOwner (void)

void	Set (const TplDriveCaller< T > *pDC)

TplDriveCaller< T > *	pGetDriveCaller (void) const

T	Get (const doublereal &dVar) const

T	Get (void) const

virtual bool	bIsDifferentiable (void) const

virtual T	GetP (void) const

Private Attributes
ConstLawType::Type	CLType

Tder	FDE0

Tder	FDEPrime0

std::vector< const DifferentiableScalarFunction * >	FDEsc

std::vector< const DifferentiableScalarFunction * >	FDEPrimesc

Additional Inherited Members
Public Types inherited from ConstitutiveLaw< T, Tder >
typedef ConstitutiveLaw< T, Tder >::ErrNotAvailable	Err

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

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

Protected Attributes inherited from ElasticConstitutiveLaw< T, Tder >
T	PreStress

Protected Attributes inherited from ConstitutiveLaw< T, Tder >
T	Epsilon

T	EpsilonPrime

T	F

Tder	FDE

Tder	FDEPrime

Protected Attributes inherited from WithLabel
unsigned int	uLabel

std::string	sName

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

Detailed Description

template<class T, class Tder>
class NLPViscoElasticConstitutiveLaw< T, Tder >

Definition at line 46 of file constltp_nlp.cc.

Constructor & Destructor Documentation

template<class T , class Tder >

NLPViscoElasticConstitutiveLaw< T, Tder >::NLPViscoElasticConstitutiveLaw	(	const TplDriveCaller< T > *	pDC,
		const T &	PStress,
		const ConstLawType::Type &	cltype,
		const Tder &	fde0,
		const std::vector< const DifferentiableScalarFunction * > &	fdesc,
		const Tder &	fdeprime0,
		const std::vector< const DifferentiableScalarFunction * > &	fdeprimesc
	)

inline

Definition at line 56 of file constltp_nlp.cc.

References NLPViscoElasticConstitutiveLaw< T, Tder >::FDE0, and NLPViscoElasticConstitutiveLaw< T, Tder >::FDEPrime0.

         : ElasticConstitutiveLaw<T, Tder>(pDC, PStress),
         CLType(cltype),
         FDE0(fde0), FDEPrime0(fdeprime0),
         FDEsc(fdesc), FDEPrimesc(fdeprimesc)
         {
                 ConstitutiveLaw<T, Tder>::FDE = FDE0;
                 ConstitutiveLaw<T, Tder>::FDEPrime = FDEPrime0;
         };

template<class T , class Tder >

virtual NLPViscoElasticConstitutiveLaw< T, Tder >::~NLPViscoElasticConstitutiveLaw ( void )

inlinevirtual

Definition at line 72 of file constltp_nlp.cc.

References NO_OP.

                                                       {
                 NO_OP;
         };

Member Function Documentation

template<class T , class Tder >

ConstLawType::Type NLPViscoElasticConstitutiveLaw< T, Tder >::GetConstLawType ( void ) const

inlinevirtual

Reimplemented from ElasticConstitutiveLaw< T, Tder >.

Definition at line 76 of file constltp_nlp.cc.

References NLPViscoElasticConstitutiveLaw< T, Tder >::CLType.

                                                      {
                 return CLType;
         };

template<class T , class Tder >

virtual ConstitutiveLaw<T, Tder>* NLPViscoElasticConstitutiveLaw< T, Tder >::pCopy ( void ) const

inlinevirtual

Implements ConstitutiveLaw< T, Tder >.

Definition at line 80 of file constltp_nlp.cc.

References NLPViscoElasticConstitutiveLaw< T, Tder >::CLType, NLPViscoElasticConstitutiveLaw< T, Tder >::FDE0, NLPViscoElasticConstitutiveLaw< T, Tder >::FDEPrime0, NLPViscoElasticConstitutiveLaw< T, Tder >::FDEPrimesc, NLPViscoElasticConstitutiveLaw< T, Tder >::FDEsc, and SAFENEWWITHCONSTRUCTOR.

Referenced by NLPViscoElasticConstitutiveLaw< doublereal, doublereal >::pCopy().

                                                             {
                 ConstitutiveLaw<T, Tder>* pCL = NULL;
 
                 typedef NLPViscoElasticConstitutiveLaw<T, Tder> cl;
                 SAFENEWWITHCONSTRUCTOR(pCL, cl,
                         cl(ElasticConstitutiveLaw<T, Tder>::pGetDriveCaller()->pCopy(),
                                 ElasticConstitutiveLaw<T, Tder>::PreStress,
                                 CLType, FDE0, FDEsc, FDEPrime0, FDEPrimesc));
                 return pCL;
         };

template<class T , class Tder >

virtual std::ostream& NLPViscoElasticConstitutiveLaw< T, Tder >::Restart ( std::ostream & out ) const

inlinevirtual

Reimplemented from ConstitutiveLaw< T, Tder >.

Definition at line 91 of file constltp_nlp.cc.

References NLPViscoElasticConstitutiveLaw< T, Tder >::FDE0, NLPViscoElasticConstitutiveLaw< T, Tder >::FDEPrime0, NLPViscoElasticConstitutiveLaw< T, Tder >::FDEPrimesc, NLPViscoElasticConstitutiveLaw< T, Tder >::FDEsc, MBDYN_EXCEPT_ARGS, and ElasticConstitutiveLaw< T, Tder >::Restart_int().

                                                            {
                 silent_cerr("NLPViscoElasticConstitutiveLaw: Restart not implemented"
                         << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
 
 #if 0 /* ScalarFunction restart not implemented yet */
                 out << "nlp viscoelastic, ",
                         FDE0.Write(out, ", ");
 
                 for (unsigned i = 0; i < FDEsc.size(); i++) {
                         out << ", ";
                         if (FDEsc[i]) {
                                 FDEsc[i]->Restart(out);
 
                         } else {
                                 out << "null";
                         }
                 }
 
                 out << ", ", FDEPrime0.Write(out, ", ");
 
                 for (unsigned i = 0; i < FDEPrimesc.size(); i++) {
                         out << ", ";
                         if (FDEsc[i]) {
                                 FDEPrimesc[i]->Restart(out);
 
                         } else {
                                 out << "null";
                         }
                 }
 
                 return ElasticConstitutiveLaw<T, Tder>::Restart_int(out);
 #endif
         };

Here is the call graph for this function:

template<class T , class Tder >

virtual void NLPViscoElasticConstitutiveLaw< T, Tder >::Update	(	const T &	Eps,
		const T &	EpsPrime = `0.`
	)

inlinevirtual

Implements ConstitutiveLaw< T, Tder >.

Definition at line 126 of file constltp_nlp.cc.

References NLPViscoElasticConstitutiveLaw< T, Tder >::CLType, ConstLawType::ELASTIC, ConstitutiveLaw< T, Tder >::F, ConstitutiveLaw< T, Tder >::FDE, NLPViscoElasticConstitutiveLaw< T, Tder >::FDE0, ConstitutiveLaw< T, Tder >::FDEPrime, NLPViscoElasticConstitutiveLaw< T, Tder >::FDEPrime0, NLPViscoElasticConstitutiveLaw< T, Tder >::FDEPrimesc, NLPViscoElasticConstitutiveLaw< T, Tder >::FDEsc, TplDriveOwner< T >::Get(), and ConstLawType::VISCOUS.

                                                                   {
                 ConstitutiveLaw<T, Tder>::Epsilon = Eps;
                 ConstitutiveLaw<T, Tder>::EpsilonPrime = EpsPrime;
 
                 T E = ElasticConstitutiveLaw<T, Tder>::Epsilon
                         - ElasticConstitutiveLaw<T, Tder>::Get();
 
                 ConstitutiveLaw<T, Tder>::F = ElasticConstitutiveLaw<T, Tder>::PreStress;
                 if (CLType & ConstLawType::ELASTIC) {
                         ConstitutiveLaw<T, Tder>::F += FDE0*E;
                         ConstitutiveLaw<T, Tder>::FDE = FDE0;
                 }
                 if (CLType & ConstLawType::VISCOUS) {
                         ConstitutiveLaw<T, Tder>::F += FDEPrime0*EpsPrime;
                         ConstitutiveLaw<T, Tder>::FDEPrime = FDEPrime0;
                 }
 
                 for (unsigned i = 0; i < FDEsc.size(); i++) {
 
                         /*
                          * each diagonal coefficient is:
                          *
                          * f = (f0' + f'(eps))*eps + (f0'' + f''(eps))*epsPrime
                          *
                          * so the Jacobian matrix contributions are
                          *
                          * f/eps = f0' + f' + f'/eps * eps + f''/eps * epsPrime
                          * f/epsPrime = f''
                          */
 
                         doublereal dEps = E(i + 1);
 
                         doublereal f = 0.;
                         doublereal fde = 0.;
 
                         if ((CLType & ConstLawType::ELASTIC) && FDEsc[i]) {
                                 doublereal df1 = (*FDEsc[i])(dEps);
                                 doublereal df1DE = FDEsc[i]->ComputeDiff(dEps);
                                 f += df1*dEps;
                                 fde += df1 + df1DE*dEps;
                         }
 
                         if ((CLType & ConstLawType::VISCOUS) && FDEPrimesc[i]) {
                                 doublereal df2 = (*FDEPrimesc[i])(dEps);
                                 doublereal df2DE = FDEPrimesc[i]->ComputeDiff(dEps);
                                 doublereal dEpsPrime = EpsPrime(i + 1);
                                 f += df2*dEpsPrime;
                                 fde += df2DE*dEpsPrime;
                                 ConstitutiveLaw<T, Tder>::FDEPrime(i + 1, i + 1) += df2;
                         }
 
                         ConstitutiveLaw<T, Tder>::F(i + 1) += f;
                         ConstitutiveLaw<T, Tder>::FDE(i + 1, i + 1) += fde;
                 }
         };