Inheritance diagram for ContContactCL:

Collaboration diagram for ContContactCL:

Public Types
enum	Type { CC_FLORES_ET_AL, CC_HUNT_CROSSLEY, CC_LANKARANI_NIKRAVESH }

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

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

Public Member Functions
	ContContactCL (const TplDriveCaller< doublereal > *pTplDC, doublereal dSign, ContContactCL::Type type, doublereal dRest, doublereal dK, doublereal dExp, doublereal dInitialEpsPrimeTol)

virtual	~ContContactCL (void)

virtual 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)

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

virtual void	AfterConvergence (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 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 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

Protected Attributes
const doublereal	m_dSign

const doublereal	m_dInitialEpsPrimeTol

const doublereal	m_dRest

ContContactCL::Type	m_type

const doublereal	m_dK

const doublereal	m_dExp

bool	m_bActive

bool	m_bToggling

doublereal	m_dInitialEpsPrime

doublereal	m_dDissCoef

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

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

Detailed Description

Definition at line 49 of file module-cont-contact.cc.

Member Enumeration Documentation

enum ContContactCL::Type

Enumerator
CC_FLORES_ET_AL
CC_HUNT_CROSSLEY
CC_LANKARANI_NIKRAVESH

Definition at line 52 of file module-cont-contact.cc.

                   {
                 CC_FLORES_ET_AL,
                 CC_HUNT_CROSSLEY,
                 CC_LANKARANI_NIKRAVESH
         };

Constructor & Destructor Documentation

ContContactCL::ContContactCL	(	const TplDriveCaller< doublereal > *	pTplDC,
		doublereal	dSign,
		ContContactCL::Type	type,
		doublereal	dRest,
		doublereal	dK,
		doublereal	dExp,
		doublereal	dInitialEpsPrimeTol
	)

inline

Definition at line 78 of file module-cont-contact.cc.

References NO_OP.

Referenced by pCopy().

         : ElasticConstitutiveLaw<doublereal, doublereal>(pTplDC, 0.),
         m_dSign(dSign), m_dInitialEpsPrimeTol(dInitialEpsPrimeTol),
         m_dRest(dRest), m_type(type), m_dK(dK), m_dExp(dExp),
         m_bActive(false), m_bToggling(false),
         m_dInitialEpsPrime(0.)
         {
                 NO_OP;
         };

virtual ContContactCL::~ContContactCL ( void )

inlinevirtual

Definition at line 90 of file module-cont-contact.cc.

References NO_OP.

                                      {
                 NO_OP;
         };

Member Function Documentation

virtual void ContContactCL::AfterConvergence	(	const doublereal &	Eps,
		const doublereal &	EpsPrime = `0.`
	)

inlinevirtual

Reimplemented from ConstitutiveLaw< doublereal, doublereal >.

Definition at line 201 of file module-cont-contact.cc.

References m_bActive, m_bToggling, m_dInitialEpsPrime, and m_dSign.

                                                                                               {
 
 #if 0
                 if (m_bActive) {
                 doublereal eps2show = m_dSign*ConstitutiveLaw<doublereal, doublereal>::Epsilon ;
                 doublereal F2show = m_dSign*ConstitutiveLaw<doublereal, doublereal>::F;
                 silent_cout( "  "        << eps2show  << "      "       << F2show << std::endl);
                 }       
 #endif          
 
                 if (m_bToggling) {
                 
 #if 0
                         silent_cout(">> ContContactCL::AfterConvergence() "
                                 "m_bToggling=" << (m_bToggling ? "true" : "false") << " "
                                 "m_bActive=" << (m_bActive ? "true" : "false") << " "
                                 "m_dInitialEpsPrime=" << m_dInitialEpsPrime << std::endl);
 #endif
 
                         if (m_bActive) {
                                 m_bActive = false;
                                 m_dInitialEpsPrime = 0.;
 
                         } else {
                                 m_bActive = true;
                         }
                         m_bToggling = false;
 #if 0
                         silent_cout("<< ContContactCL::AfterConvergence() "
                                 "m_bToggling=" << (m_bToggling ? "true" : "false") << " "
                                 "m_bActive=" << (m_bActive ? "true" : "false") << " "
                                 "m_dInitialEpsPrime=" << m_dInitialEpsPrime << std::endl);
 #endif
                 }
         };

virtual ConstLawType::Type ContContactCL::GetConstLawType ( void ) const

inlinevirtual

Implements ConstitutiveLaw< doublereal, doublereal >.

Definition at line 94 of file module-cont-contact.cc.

References ConstLawType::VISCOELASTIC.

                                                              {
                 return ConstLawType::VISCOELASTIC;
         };

virtual std::ostream& ContContactCL::OutputAppend ( std::ostream & out ) const

inlinevirtual

Reimplemented from SimulationEntity.

Definition at line 196 of file module-cont-contact.cc.

References m_dDissCoef, and m_dInitialEpsPrime.

                                                                 {
                 return out << " " << m_dInitialEpsPrime << " " << m_dDissCoef;
         };

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

inlinevirtual

Implements ConstitutiveLaw< doublereal, doublereal >.

Definition at line 98 of file module-cont-contact.cc.

References ContContactCL(), m_dExp, m_dInitialEpsPrimeTol, m_dK, m_dRest, m_dSign, m_type, TplDriveOwner< doublereal >::pGetDriveCaller(), and SAFENEWWITHCONSTRUCTOR.

                                                                            {
                 ConstitutiveLaw<doublereal, doublereal>* pCL = 0;
 
                 // pass parameters to copy constructor
                 SAFENEWWITHCONSTRUCTOR(pCL, ContContactCL,
                         ContContactCL(pGetDriveCaller()->pCopy(),
                                 m_dSign, m_type, m_dRest, m_dK, m_dExp, m_dInitialEpsPrimeTol));
 
                 return pCL;
         };

Here is the call graph for this function:

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

inlinevirtual

Reimplemented from ConstitutiveLaw< doublereal, doublereal >.

Definition at line 109 of file module-cont-contact.cc.

References CC_FLORES_ET_AL, CC_HUNT_CROSSLEY, CC_LANKARANI_NIKRAVESH, m_dExp, m_dInitialEpsPrimeTol, m_dK, m_dRest, m_dSign, m_type, and ElasticConstitutiveLaw< T, Tder >::Restart_int().

                                                            {
                 out << "continuous contact"
                         << ", sign, " << m_dSign
                         << ", formulation, ";
 
                 switch (m_type) {
                 case CC_FLORES_ET_AL:
                         out << "flores";
                         break;
 
                 case CC_HUNT_CROSSLEY:
                         out << "hunt crossley";
                         break;
 
                 case CC_LANKARANI_NIKRAVESH:
                         out << "lankarani nikravesh";
                         break;
                 }
 
                 out
                         << ", restitution, " << m_dRest
                         << ", kappa, " << m_dK
                         << ", exp, " << m_dExp
                         << ", EpsPrimeTol, " << m_dInitialEpsPrimeTol
                         << ", ", ElasticConstitutiveLaw<doublereal, doublereal>::Restart_int(out);
                 return out;
         };

Here is the call graph for this function:

virtual void ContContactCL::Update	(	const doublereal &	Eps,
		const doublereal &	EpsPrime
	)

inlinevirtual

Implements ConstitutiveLaw< doublereal, doublereal >.

Definition at line 137 of file module-cont-contact.cc.

References CC_FLORES_ET_AL, CC_HUNT_CROSSLEY, CC_LANKARANI_NIKRAVESH, TplDriveOwner< T >::Get(), m_bActive, m_bToggling, m_dDissCoef, m_dExp, m_dInitialEpsPrime, m_dInitialEpsPrimeTol, m_dK, m_dRest, m_dSign, m_type, and grad::pow().

                                                                                {
                 ConstitutiveLaw<doublereal, doublereal>::Epsilon = Eps - ElasticConstitutiveLaw<doublereal, doublereal>::Get();
                 ConstitutiveLaw<doublereal, doublereal>::EpsilonPrime = EpsPrime;
 
                 doublereal x = m_dSign*ConstitutiveLaw<doublereal, doublereal>::Epsilon;
                 if (x < 0.) {
                         if (m_bActive) {
                                 if (!m_bToggling) {
                                         m_bToggling = true;
                                 }
                         }
 
                         ConstitutiveLaw<doublereal, doublereal>::F = 0.;
                         ConstitutiveLaw<doublereal, doublereal>::FDE = 0.;
                         ConstitutiveLaw<doublereal, doublereal>::FDEPrime = 0.;
 
                 } else {
                         doublereal xp = m_dSign*ConstitutiveLaw<doublereal, doublereal>::EpsilonPrime;
 
                         if (!m_bActive && !m_bToggling) {
                                 m_bToggling = true;
                                 m_dInitialEpsPrime = std::max(xp, m_dInitialEpsPrimeTol);
 
                                 switch (m_type) {
                                 case CC_FLORES_ET_AL:
                                         m_dDissCoef = (std::abs(xp) > m_dInitialEpsPrimeTol) * 8./5. * m_dK * (1. - m_dRest) / m_dRest / m_dInitialEpsPrime;
                                         break;
 
                                 case CC_HUNT_CROSSLEY:
                                         m_dDissCoef = (std::abs(xp) > m_dInitialEpsPrimeTol) * 3./2. * m_dK * (1. - m_dRest) / m_dInitialEpsPrime;
                                         break;
 
                                 case CC_LANKARANI_NIKRAVESH:
                                         m_dDissCoef = (std::abs(xp) > m_dInitialEpsPrimeTol) * 3./4. * m_dK * (1. - std::pow(m_dRest, 2)) / m_dInitialEpsPrime;
                                         break;
                                 }
                         }
 
 //                      doublereal xn   = std::pow(x, m_dExp);
 //                      doublereal xnm1 = std::pow(x, m_dExp - 1.);     
 
 //                      doublereal xp_t = ( ( (m_dDissCoef*xp < 0.)&&(abs(m_dDissCoef*xp) > m_dK)) ? abs(m_dK/m_dDissCoef/xp) : 1. );   // Threshold
 
 //                      ConstitutiveLaw<doublereal, doublereal>::F = m_dSign*(m_dK*xn + m_dDissCoef *xn* xp * xp_t);    // FIXME
 //                      ConstitutiveLaw<doublereal, doublereal>::FDE = m_dSign*(m_dExp*m_dK*xnm1 + m_dDissCoef *m_dExp*xnm1* xp_t);     //FIXME
 //                      ConstitutiveLaw<doublereal, doublereal>::FDEPrime = m_dSign*(m_dDissCoef * xn) * xp_t ;
                         
 
                         doublereal xn   = std::pow(x, m_dExp);
                         doublereal xnm1 = std::pow(x, m_dExp - 1.);                     
                         
                         ConstitutiveLaw<doublereal, doublereal>::F = m_dSign*xn*(m_dK + m_dDissCoef*xp);        // FIXME
                         ConstitutiveLaw<doublereal, doublereal>::FDE = m_dSign*xnm1*(m_dExp*m_dK + m_dDissCoef*m_dExp); 
                         ConstitutiveLaw<doublereal, doublereal>::FDEPrime = m_dSign*xn*m_dDissCoef;
                         
 //                      silent_cout("hello "  << std::endl);
                 }
         };