Inheritance diagram for NLSFViscoElasticConstitutiveLaw< doublereal, doublereal >:

Collaboration diagram for NLSFViscoElasticConstitutiveLaw< doublereal, doublereal >:

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

virtual	~NLSFViscoElasticConstitutiveLaw (void)

ConstLawType::Type	GetConstLawType (void) const

virtual ConstitutiveLaw < doublereal, doublereal > *	pCopy (void) const

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

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

Public Member Functions inherited from ElasticConstitutiveLaw< doublereal, doublereal >
	ElasticConstitutiveLaw (const TplDriveCaller< doublereal > *pDC, const doublereal &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< doublereal, doublereal >
	ConstitutiveLaw (void)

virtual	~ConstitutiveLaw (void)

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

virtual const doublereal &	GetEpsilon (void) const

virtual const doublereal &	GetEpsilonPrime (void) const

virtual const doublereal &	GetF (void) const

virtual const doublereal &	GetFDE (void) const

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

virtual	~TplDriveOwner (void)

void	Set (const TplDriveCaller< doublereal > *pDC)

TplDriveCaller< doublereal > *	pGetDriveCaller (void) const

doublereal	Get (const doublereal &dVar) const

doublereal	Get (void) const

virtual bool	bIsDifferentiable (void) const

virtual doublereal	GetP (void) const

Private Attributes
ConstLawType::Type	CLType

doublereal	FDE0

doublereal	FDEPrime0

std::vector< const DifferentiableScalarFunction * >	FDEsc

std::vector< const DifferentiableScalarFunction * >	FDEPrimesc

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

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

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

Protected Attributes inherited from ElasticConstitutiveLaw< doublereal, doublereal >
doublereal	PreStress

Protected Attributes inherited from ConstitutiveLaw< doublereal, doublereal >
doublereal	Epsilon

doublereal	EpsilonPrime

doublereal	F

doublereal	FDE

doublereal	FDEPrime

Protected Attributes inherited from WithLabel
unsigned int	uLabel

std::string	sName

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

Detailed Description

template<>
class NLSFViscoElasticConstitutiveLaw< doublereal, doublereal >

Definition at line 178 of file constltp_nlsf.cc.

Constructor & Destructor Documentation

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

inline

Definition at line 188 of file constltp_nlsf.cc.

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

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

virtual NLSFViscoElasticConstitutiveLaw< doublereal, doublereal >::~NLSFViscoElasticConstitutiveLaw ( void )

inlinevirtual

Definition at line 204 of file constltp_nlsf.cc.

References NO_OP.

                                                        {
                 NO_OP;
         };

Member Function Documentation

ConstLawType::Type NLSFViscoElasticConstitutiveLaw< doublereal, doublereal >::GetConstLawType ( void ) const

inlinevirtual

Implements ConstitutiveLaw< doublereal, doublereal >.

Definition at line 208 of file constltp_nlsf.cc.

References NLSFViscoElasticConstitutiveLaw< T, Tder >::CLType.

                                                      {
                 return CLType;
         };

virtual ConstitutiveLaw<doublereal, doublereal>* NLSFViscoElasticConstitutiveLaw< doublereal, doublereal >::pCopy ( void ) const

inlinevirtual

Implements ConstitutiveLaw< doublereal, doublereal >.

Definition at line 212 of file constltp_nlsf.cc.

References NLSFViscoElasticConstitutiveLaw< T, Tder >::CLType, NLSFViscoElasticConstitutiveLaw< T, Tder >::FDE0, NLSFViscoElasticConstitutiveLaw< T, Tder >::FDEPrime0, NLSFViscoElasticConstitutiveLaw< T, Tder >::FDEPrimesc, NLSFViscoElasticConstitutiveLaw< T, Tder >::FDEsc, NLSFViscoElasticConstitutiveLaw< T, Tder >::pCopy(), and SAFENEWWITHCONSTRUCTOR.

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

Here is the call graph for this function:

virtual std::ostream& NLSFViscoElasticConstitutiveLaw< doublereal, doublereal >::Restart ( std::ostream & out ) const

inlinevirtual

Reimplemented from ConstitutiveLaw< doublereal, doublereal >.

Definition at line 223 of file constltp_nlsf.cc.

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

                                                            {
                 silent_cerr("NLSFViscoElasticConstitutiveLaw: 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:

virtual void NLSFViscoElasticConstitutiveLaw< doublereal, doublereal >::Update	(	const doublereal &	Eps,
		const doublereal &	EpsPrime = `0.`
	)

inlinevirtual

Implements ConstitutiveLaw< doublereal, doublereal >.

Definition at line 258 of file constltp_nlsf.cc.

References NLSFViscoElasticConstitutiveLaw< T, Tder >::CLType, ConstLawType::ELASTIC, NLSFViscoElasticConstitutiveLaw< T, Tder >::FDE0, NLSFViscoElasticConstitutiveLaw< T, Tder >::FDEPrime0, NLSFViscoElasticConstitutiveLaw< T, Tder >::FDEPrimesc, NLSFViscoElasticConstitutiveLaw< T, Tder >::FDEsc, TplDriveOwner< T >::Get(), and ConstLawType::VISCOUS.

                                                                                     {
                 ConstitutiveLaw<doublereal, doublereal>::Epsilon = Eps;
                 ConstitutiveLaw<doublereal, doublereal>::EpsilonPrime = EpsPrime;
 
                 doublereal dEps = ElasticConstitutiveLaw<doublereal, doublereal>::Epsilon
                         - ElasticConstitutiveLaw<doublereal, doublereal>::Get();
                 doublereal dEpsPrime = EpsPrime;
 
                 ConstitutiveLaw<doublereal, doublereal>::F = ElasticConstitutiveLaw<doublereal, doublereal>::PreStress;
                 if (CLType & ConstLawType::ELASTIC) {
                         ConstitutiveLaw<doublereal, doublereal>::F += FDE0*dEps;
                         ConstitutiveLaw<doublereal, doublereal>::FDE = FDE0;
                 }
                 if (CLType & ConstLawType::VISCOUS) {
                         ConstitutiveLaw<doublereal, doublereal>::F += FDEPrime0*dEpsPrime;
                         ConstitutiveLaw<doublereal, doublereal>::FDEPrime = FDEPrime0;
                 }
 
                 /*
                  * 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 f = 0.;
 
                 if ((CLType & ConstLawType::ELASTIC) && FDEsc[0]) {
                         f += (*FDEsc[0])(dEps);
                         doublereal df1DE = FDEsc[0]->ComputeDiff(dEps);
                         ConstitutiveLaw<doublereal, doublereal>::FDE += df1DE;
                 }
 
                 if ((CLType & ConstLawType::VISCOUS) && FDEPrimesc[0]) {
                         f += (*FDEPrimesc[0])(dEpsPrime);
                         doublereal df2DEPrime = FDEPrimesc[0]->ComputeDiff(dEpsPrime);
                         ConstitutiveLaw<doublereal, doublereal>::FDEPrime += df2DEPrime;
                 }
 
                 ConstitutiveLaw<doublereal, doublereal>::F += f;
         };