#include <ginaccltp.h>

Inheritance diagram for GiNaCElasticConstitutiveLaw< doublereal, doublereal >:

Collaboration diagram for GiNaCElasticConstitutiveLaw< doublereal, doublereal >:

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

virtual	~GiNaCElasticConstitutiveLaw (void)

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

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

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

Public Member Functions inherited from SymbolicElasticConstitutiveLaw< doublereal, doublereal >
	SymbolicElasticConstitutiveLaw (const TplDriveCaller< doublereal > *pDC, const doublereal &PStress, std::vector< std::string > &epsilon, std::vector< std::string > &expression)

virtual	~SymbolicElasticConstitutiveLaw (void)

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
GiNaC::symbol	gEps

GiNaC::ex	gExpr

GiNaC::ex	gExprDEps

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 GiNaCElasticConstitutiveLaw< doublereal, doublereal >

Definition at line 242 of file ginaccltp.h.

Constructor & Destructor Documentation

GiNaCElasticConstitutiveLaw< doublereal, doublereal >::GiNaCElasticConstitutiveLaw	(	const TplDriveCaller< doublereal > *	pDC,
		const doublereal &	PStress,
		std::vector< std::string > &	epsilon,
		std::vector< std::string > &	expression
	)

Definition at line 262 of file ginaccltp.h.

References GiNaCElasticConstitutiveLaw< T, Tder >::gEps, GiNaCElasticConstitutiveLaw< T, Tder >::gExpr, GiNaCElasticConstitutiveLaw< T, Tder >::gExprDEps, and GiNaCElasticConstitutiveLaw< T, Tder >::Update().

 : SymbolicElasticConstitutiveLaw<doublereal, doublereal>(pDC, PStress, epsilon, expression),
 gEps(epsilon[0])
 { 
         ConstitutiveLaw<doublereal, doublereal>::FDE = 0.;
 
         GiNaC::lst l; 
 
         l.append(gEps);
 
         try {
                 gExpr = GiNaC::ex(expression[0], l);
 
         } catch (std::exception& e) {
                 silent_cerr("GiNaCElasticConstitutiveLaw<doublereal, doublereal>: expression parsing failed: " << e.what() << std::endl);
                 throw e;
         }
 
         try {
                 gExprDEps = gExpr.diff(gEps);
 
         } catch (std::exception& e) {
                 silent_cerr("GiNaCElasticConstitutiveLaw<doublereal, doublereal>: expression differentiation wrt/ Eps failed: " << e.what() << std::endl);
                 throw e;
         }
 
         silent_cout("\tGiNaCElasticConstitutiveLaw:" << std::endl
                 << "\t\tEps:              \"" << gEps << "\"" << std::endl
                 << "\t\tConstitutive law: \"" << gExpr << "\"" << std::endl
                 << "\t\tDer/Eps:          \"" << gExprDEps << "\"" << std::endl);
 
         // try and evaluate the constitutive law
         try {
                 Update(ElasticConstitutiveLaw<doublereal, doublereal>::Epsilon, ElasticConstitutiveLaw<doublereal, doublereal>::Epsilon);
         }
         catch (std::exception& e) {
                 silent_cerr("GiNaCElasticConstitutiveLaw<doublereal, doublereal>::GiNaCElasticConstitutiveLaw: Update() failed (" << e.what() << ")" << std::endl);
                 throw e;
         }
 }

Here is the call graph for this function:

GiNaCElasticConstitutiveLaw< doublereal, doublereal >::~GiNaCElasticConstitutiveLaw ( void )

virtual

Definition at line 307 of file ginaccltp.h.

References NO_OP.

 {
         NO_OP;
 };

Member Function Documentation

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

virtual

Implements ConstitutiveLaw< doublereal, doublereal >.

Definition at line 313 of file ginaccltp.h.

References GiNaCElasticConstitutiveLaw< T, Tder >::gEps, GiNaCElasticConstitutiveLaw< T, Tder >::gExpr, GiNaCElasticConstitutiveLaw< T, Tder >::pCopy(), and SAFENEWWITHCONSTRUCTOR.

 {
         ConstitutiveLaw<doublereal, doublereal>* pCL = 0;
 
         std::ostringstream      eps;
         std::ostringstream      expr;
 
         eps << gEps;
         expr << gExpr;
 
         std::vector<std::string> epsilon;
         std::vector<std::string> expression;
 
         epsilon.push_back(eps.str());
         expression.push_back(expr.str());
 
         typedef GiNaCElasticConstitutiveLaw<doublereal, doublereal> cl;
         SAFENEWWITHCONSTRUCTOR(pCL, 
                 cl, 
                 cl(ElasticConstitutiveLaw<doublereal, doublereal>::pGetDriveCaller()->pCopy(), 
                         ElasticConstitutiveLaw<doublereal, doublereal>::PreStress,
                         epsilon, expression));
       
         return pCL;
 }

Here is the call graph for this function:

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

virtual

Reimplemented from ConstitutiveLaw< doublereal, doublereal >.

Definition at line 340 of file ginaccltp.h.

References GiNaCElasticConstitutiveLaw< T, Tder >::gEps, GiNaCElasticConstitutiveLaw< T, Tder >::gExpr, and ElasticConstitutiveLaw< T, Tder >::Restart_int().

 {
         out << "symbolic elastic isotropic, epsilon, \"" << gEps
                 << "\", expression, \"" << gExpr << "\"";
 
         return Restart_int(out);
 }

Here is the call graph for this function:

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

virtual

Implements ConstitutiveLaw< doublereal, doublereal >.

Definition at line 349 of file ginaccltp.h.

References GiNaCElasticConstitutiveLaw< T, Tder >::gEps, TplDriveOwner< T >::Get(), GiNaCElasticConstitutiveLaw< T, Tder >::gExpr, and GiNaCElasticConstitutiveLaw< T, Tder >::gExprDEps.

 {
         GiNaC::lst l;
 
         ElasticConstitutiveLaw<doublereal, doublereal>::Epsilon = Eps;
 
         doublereal e = Eps - ElasticConstitutiveLaw<doublereal, doublereal>::Get();
 
         l.append(gEps == e);
 
         GiNaC::ex f_expr = gExpr.subs(l);
         ConstitutiveLaw<doublereal, doublereal>::F = 
                 ElasticConstitutiveLaw<doublereal, doublereal>::PreStress
                 + GiNaC::ex_to<GiNaC::numeric>(f_expr).to_double();
 
         GiNaC::ex f_derEps = gExprDEps.subs(l);
         ConstitutiveLaw<doublereal, doublereal>::FDE = GiNaC::ex_to<GiNaC::numeric>(f_derEps).to_double();
 }