#include <ginaccltp.h>

Inheritance diagram for GiNaCViscoElasticConstitutiveLaw< T, Tder >:

Collaboration diagram for GiNaCViscoElasticConstitutiveLaw< T, Tder >:

Public Member Functions
	GiNaCViscoElasticConstitutiveLaw (const TplDriveCaller< T > *pDC, const T &PStress, std::vector< std::string > &epsilon, std::vector< std::string > &epsilonPrime, std::vector< std::string > &expression)

virtual	~GiNaCViscoElasticConstitutiveLaw (void)

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

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

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

Public Member Functions inherited from SymbolicViscoElasticConstitutiveLaw< T, Tder >
	SymbolicViscoElasticConstitutiveLaw (const TplDriveCaller< T > *pDC, const T &PStress, std::vector< std::string > &epsilon, std::vector< std::string > &epsilonPrime, std::vector< std::string > &expression)

virtual	~SymbolicViscoElasticConstitutiveLaw (void)

virtual ConstLawType::Type	GetConstLawType (void) const

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
unsigned	dim

std::vector< GiNaC::symbol * >	gEps

std::vector< GiNaC::symbol * >	gEpsPrime

std::vector< GiNaC::ex >	gExpr

std::vector< std::vector < GiNaC::ex > >	gExprDEps

std::vector< std::vector < GiNaC::ex > >	gExprDEpsPrime

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 GiNaCViscoElasticConstitutiveLaw< T, Tder >

Definition at line 691 of file ginaccltp.h.

Constructor & Destructor Documentation

template<class T , class Tder >

GiNaCViscoElasticConstitutiveLaw< T, Tder >::GiNaCViscoElasticConstitutiveLaw	(	const TplDriveCaller< T > *	pDC,
		const T &	PStress,
		std::vector< std::string > &	epsilon,
		std::vector< std::string > &	epsilonPrime,
		std::vector< std::string > &	expression
	)

Definition at line 724 of file ginaccltp.h.

References GiNaCViscoElasticConstitutiveLaw< T, Tder >::dim, GiNaCViscoElasticConstitutiveLaw< T, Tder >::gEps, GiNaCViscoElasticConstitutiveLaw< T, Tder >::gEpsPrime, GiNaCViscoElasticConstitutiveLaw< T, Tder >::gExpr, GiNaCViscoElasticConstitutiveLaw< T, Tder >::gExprDEps, GiNaCViscoElasticConstitutiveLaw< T, Tder >::gExprDEpsPrime, MBDYN_EXCEPT_ARGS, and GiNaCViscoElasticConstitutiveLaw< T, Tder >::Update().

 : SymbolicViscoElasticConstitutiveLaw<T, Tder>(pDC, PStress, epsilon, epsilonPrime, expression)
 { 
         if (typeid(T) == typeid(Vec3)) {
                 dim = 3;
 
         } else if (typeid(T) == typeid(Vec6)) {
                 dim = 6;
 
         } else {
                 throw (typename ConstitutiveLaw<T, Tder>::ErrNotAvailable(MBDYN_EXCEPT_ARGS)); 
         }
 
         gEps.resize(dim);
         gEpsPrime.resize(dim);
         gExpr.resize(dim);
 
         gExprDEps.resize(dim);
         gExprDEpsPrime.resize(dim);
         for (unsigned row = 0; row < dim; row++) {
                 gExprDEps[row].resize(dim);
                 gExprDEpsPrime[row].resize(dim);
         }
 
         ConstitutiveLaw<T, Tder>::FDE = mb_zero<Tder>();
         ConstitutiveLaw<T, Tder>::FDEPrime = mb_zero<Tder>();
 
         GiNaC::lst l; 
 
         for (unsigned row = 0; row < dim; row++) {
                 gEps[row] = new GiNaC::symbol(epsilon[row]);
                 l.append(*gEps[row]);
                 gEpsPrime[row] = new GiNaC::symbol(epsilonPrime[row]);
                 l.append(*gEpsPrime[row]);
         }
 
         for (unsigned row = 0; row < dim; row++) {
                 try {
                         gExpr[row] = GiNaC::ex(expression[row], l);
 
                 } catch (std::exception& e) {
                         silent_cerr("GiNaCViscoElasticConstitutiveLaw<T, Tder>: expression #" << row << " parsing "
                                 "failed: " << e.what() << std::endl);
                         throw e;
                 }
 
                 for (unsigned col = 0; col < dim; col++) {
                         try {
                                 gExprDEps[row][col] = gExpr[row].diff(*gEps[col]);
 
                         } catch (std::exception& e) {
                                 silent_cerr("GiNaCViscoElasticConstitutiveLaw<T, Tder>: expression #" << row << " differentiation "
                                         "wrt/ Eps #" << col << "failed: "
                                         << e.what() << std::endl);
                                 throw e;
                         }
 
                         try {
                                 gExprDEpsPrime[row][col] = gExpr[row].diff(*gEpsPrime[col]);
 
                         } catch (std::exception& e) {
                                 silent_cerr("GiNaCViscoElasticConstitutiveLaw<T, Tder>: expression #" << row << " differentiation "
                                         "wrt/ EpsPrime #" << col << "failed: "
                                         << e.what() << std::endl);
                                 throw e;
                         }
                 }
         }
 
         silent_cout("\tGiNaCViscoElasticConstitutiveLaw:" << std::endl);
         for (unsigned row = 0; row < dim; row++) {
                 silent_cout("\t\tEps[" << row << "]:              \"" << *gEps[row] << "\"" << std::endl);
         }
         for (unsigned row = 0; row < dim; row++) {
                 silent_cout("\t\tEpsPrime[" << row << "]:              \"" << *gEpsPrime[row] << "\"" << std::endl);
         }
         for (unsigned row = 0; row < dim; row++) {
                 silent_cout("\t\tConstitutive law[" << row << "]: \"" << gExpr[row] << "\"" << std::endl);
         }
         for (unsigned row = 0; row < dim; row++) {
                 for (unsigned col = 0; col < dim; col++) {
                         silent_cout("\t\tDer[" << row << "]/Eps[" << row << "][" << col << "]:          \"" << gExprDEps[row][col] << "\"" << std::endl);
                 }
         }
         for (unsigned row = 0; row < dim; row++) {
                 for (unsigned col = 0; col < dim; col++) {
                         silent_cout("\t\tDer[" << row << "]/EpsPrime[" << row << "][" << col << "]:          \"" << gExprDEpsPrime[row][col] << "\"" << std::endl);
                 }
         }
 
         // try and evaluate the constitutive law
         try {
                 Update(ElasticConstitutiveLaw<T, Tder>::Epsilon, ElasticConstitutiveLaw<T, Tder>::Epsilon);
         }
         catch (std::exception& e) {
                 silent_cerr("GiNaCViscoElasticConstitutiveLaw<T, Tder>::GiNaCViscoElasticConstitutiveLaw: Update() failed (" << e.what() << ")" << std::endl);
                 throw e;
         }
 }

Here is the call graph for this function:

template<class T , class Tder >

GiNaCViscoElasticConstitutiveLaw< T, Tder >::~GiNaCViscoElasticConstitutiveLaw ( void )

virtual

Definition at line 830 of file ginaccltp.h.

 {
         for (unsigned row = 0; row < dim; row++) {
                 delete gEps[row];
                 delete gEpsPrime[row];
         }
 };

Member Function Documentation

template<class T , class Tder >

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

virtual

Implements ConstitutiveLaw< T, Tder >.

Definition at line 839 of file ginaccltp.h.

References SAFENEWWITHCONSTRUCTOR.

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

 {
         ConstitutiveLaw<T, Tder>* pCL = 0;
 
         std::vector<std::string> epsilon(dim);
         std::vector<std::string> epsilonPrime(dim);
         std::vector<std::string> expression(dim);
 
         for (unsigned row = 0; row < dim; row++) {
                 std::ostringstream      eps;
                 std::ostringstream      epsPrime;
                 std::ostringstream      expr;
 
                 eps << *gEps[row];
                 epsPrime << *gEpsPrime[row];
                 expr << gExpr[row];
 
                 epsilon[row] = eps.str();
                 epsilonPrime[row] = epsPrime.str();
                 expression[row] = expr.str();
         }
 
         typedef GiNaCViscoElasticConstitutiveLaw<T, Tder> cl;
         SAFENEWWITHCONSTRUCTOR(pCL, 
                 cl, 
                 cl(ElasticConstitutiveLaw<T, Tder>::pGetDriveCaller()->pCopy(), 
                         ElasticConstitutiveLaw<T, Tder>::PreStress,
                         epsilon, epsilonPrime, expression));
       
         return pCL;
 }

template<class T , class Tder >

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

virtual

Reimplemented from ConstitutiveLaw< T, Tder >.

Definition at line 872 of file ginaccltp.h.

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

 {
         out << "symbolic viscoelastic, epsilon";
         for (unsigned row = 0; row < dim; row++) {
                 out << ", \"" << *gEps[row] << "\"";
         }
         out << "\", epsilon prime";
         for (unsigned row = 0; row < dim; row++) {
                 out << ", \"" << *gEpsPrime[row] << "\"";
         }
         out << "\", expression";
         for (unsigned row = 0; row < dim; row++) {
                 out << ", \"" << gExpr[row] << "\"";
         }
 
         return ElasticConstitutiveLaw<T, Tder>::Restart_int(out);
 }

Here is the call graph for this function:

template<class T , class Tder >

void GiNaCViscoElasticConstitutiveLaw< T, Tder >::Update	(	const T &	Eps,
		const T &	EpsPrime = `mb_zero<T>()`
	)

virtual

Implements ConstitutiveLaw< T, Tder >.

Definition at line 891 of file ginaccltp.h.

References ConstitutiveLaw< T, Tder >::F, ConstitutiveLaw< T, Tder >::FDE, ConstitutiveLaw< T, Tder >::FDEPrime, and TplDriveOwner< T >::Get().

Referenced by GiNaCViscoElasticConstitutiveLaw< T, Tder >::GiNaCViscoElasticConstitutiveLaw(), and GiNaCViscoElasticConstitutiveLaw< doublereal, doublereal >::GiNaCViscoElasticConstitutiveLaw().

 {
         GiNaC::lst l;
 
         ElasticConstitutiveLaw<T, Tder>::Epsilon = Eps;
         ConstitutiveLaw<T, Tder>::EpsilonPrime = EpsPrime;
 
         T e = Eps - ElasticConstitutiveLaw<T, Tder>::Get();
 
         for (unsigned row = 0; row < dim; row++) {
                 l.append(*gEps[row] == e(row + 1));
                 l.append(*gEpsPrime[row] == EpsPrime(row + 1));
         }
 
         ConstitutiveLaw<T, Tder>::F = ElasticConstitutiveLaw<T, Tder>::PreStress;
 
         for (unsigned row = 0; row < dim; row++) {
                 GiNaC::ex f_expr = gExpr[row].subs(l);
 
                 ConstitutiveLaw<T, Tder>::F(row + 1)
                         += GiNaC::ex_to<GiNaC::numeric>(f_expr).to_double();
 
                 for (unsigned col = 0; col < dim; col++) {
                         GiNaC::ex f_derEps = gExprDEps[row][col].subs(l);
 
                         ConstitutiveLaw<T, Tder>::FDE(row + 1, col + 1)
                                 = GiNaC::ex_to<GiNaC::numeric>(f_derEps).to_double();
 
                         GiNaC::ex f_derEpsPrime = gExprDEpsPrime[row][col].subs(l);
 
                         ConstitutiveLaw<T, Tder>::FDEPrime(row + 1, col + 1)
                                 = GiNaC::ex_to<GiNaC::numeric>(f_derEpsPrime).to_double();
                 }
         }
 }