Frictional+ViscoElastic interaction between particles. More...

#include <HertzMindlinViscoInteraction.h>

Inheritance diagram for CHertzMindlinViscoInteraction:

Collaboration diagram for CHertzMindlinViscoInteraction:

[legend]

Public Types
typedef CHertzMindlinViscoIGP	ParameterType

typedef double(CHertzMindlinViscoInteraction::*	ScalarFieldFunction) () const

typedef std::pair< bool, double >(CHertzMindlinViscoInteraction::*	CheckedScalarFieldFunction) () const

typedef Vec3(CHertzMindlinViscoInteraction::*	VectorFieldFunction) () const

Public Member Functions
	CHertzMindlinViscoInteraction ()

	CHertzMindlinViscoInteraction (CParticle , CParticle )

	CHertzMindlinViscoInteraction (CParticle , CParticle , const CHertzMindlinViscoIGP &)

virtual	~CHertzMindlinViscoInteraction ()

virtual void	calcForces ()

virtual bool	isPersistent ()

void	setTimeStepSize (double dt)

std::pair< bool, double >	getAbsFrictionalForce () const

std::pair< bool, double >	getAbsFN () const

std::pair< bool, double >	getAbsMuFN () const

std::pair< bool, double >	getSlipVelocity () const

std::pair< bool, double >	getNormalStress () const

std::pair< bool, double >	getMaxFricStress () const

std::pair< bool, double >	getAbsFrictionalStress () const

double	getAbsForceDeficit () const

double	getPotentialEnergy () const

double	getSlipping () const

double	getSticking () const

double	getDissipatedEnergy () const

double	Count () const

Vec3	getForce () const

Vec3	getNormalForce () const

virtual Vec3	getPos () const

std::pair< bool, double >	getMuEff (const Vec3 &, const Vec3 &) const

std::pair< bool, double >	getMuEffXY () const

std::pair< bool, double >	getMuEffXZ () const

Public Member Functions inherited from APairInteraction
	APairInteraction ()

	APairInteraction (CParticle , CParticle )

virtual	~APairInteraction ()

const CParticle *	first () const

const CParticle *	second () const

CParticle *	first ()

CParticle *	second ()

pair< int, int >	getPairID () const

void	setPP (CParticle , CParticle )

void	checkIDs ()

virtual bool	hasTag (int, int) const

virtual Vec3	getPosFirst () const

virtual Vec3	getPosSecond () const

esys::lsm::quintuple< Vec3, double, Vec3, double, Vec3 >	getRaw2Data () const

template<class TmplParticle >
void	setPP (const vector< TmplParticle * > &pp)

virtual void	saveRestartData (std::ostream &)

virtual void	loadRestartData (std::istream &)

Public Member Functions inherited from AInteraction
	AInteraction ()

virtual	~AInteraction ()

bool	initialized () const

virtual void	calcHeatTrans ()

virtual void	calcHeatFrict ()

vector< int >	getAllID () const

int	getID ()

double	Count () const

Static Public Member Functions
static ScalarFieldFunction	getScalarFieldFunction (const string &)

static CheckedScalarFieldFunction	getCheckedScalarFieldFunction (const string &)

static VectorFieldFunction	getVectorFieldFunction (const string &)

static string	getType ()

Protected Attributes
double	m_E
	Young's modulus. More...

double	m_nu
	Poisson ratio. More...

double	m_r0
	equilibrium distance More...

double	m_mu
	coefficient of friction More...

double	m_COR
	coefficient of restitution More...

double	m_dt
	time step More...

Vec3	m_Ffric
	current frictional force More...

Vec3	m_force_deficit
	difference between fric. force & force necessary for slip More...

Vec3	m_cpos
	contact position More...

Vec3	m_normal_force
	current normal force More...

double	m_dn
	caching displacement for E_pot More...

bool	m_is_slipping
	static/dynamic status of the interaction More...

bool	m_is_touching
	contact status of the interaction More...

double	m_E_diss
	dissipated energy More...

Protected Attributes inherited from APairInteraction
CParticle *	m_p1

CParticle *	m_p2

Protected Attributes inherited from AInteraction
vector< int >	m_id
	id's of the particles involved More...

int	m_iid
	interaction id More...

bool	m_init

Friends
class	TML_PackedMessageInterface

std::ostream &	operator<< (std::ostream &, const CHertzMindlinViscoInteraction &)

Detailed Description

Frictional+ViscoElastic interaction between particles.

Author: Qi Shao & Dion Weatherley

Date: 2016-12-14 10:00:00 (wed, 14 dec 2016)

Member Typedef Documentation

◆ CheckedScalarFieldFunction

typedef std::pair<bool,double>(CHertzMindlinViscoInteraction::* CHertzMindlinViscoInteraction::CheckedScalarFieldFunction) () const

◆ ParameterType

typedef CHertzMindlinViscoIGP CHertzMindlinViscoInteraction::ParameterType

◆ ScalarFieldFunction

typedef double(CHertzMindlinViscoInteraction::* CHertzMindlinViscoInteraction::ScalarFieldFunction) () const

◆ VectorFieldFunction

typedef Vec3(CHertzMindlinViscoInteraction::* CHertzMindlinViscoInteraction::VectorFieldFunction) () const

Constructor & Destructor Documentation

◆ CHertzMindlinViscoInteraction() [1/3]

CHertzMindlinViscoInteraction::CHertzMindlinViscoInteraction ( )

◆ CHertzMindlinViscoInteraction() [2/3]

CHertzMindlinViscoInteraction::CHertzMindlinViscoInteraction	(	CParticle *	p1,
		CParticle *	p2
	)

constructor for CHertzMindlinViscoInteraction without friction parameters, only calls the constructor of APairInteraction with the 2 particle pointers

◆ CHertzMindlinViscoInteraction() [3/3]

CHertzMindlinViscoInteraction::CHertzMindlinViscoInteraction	(	CParticle *	p1,
		CParticle *	p2,
		const CHertzMindlinViscoIGP &	param
	)

References CHertzMindlinViscoIGP::dt, CBasicParticle::getPos(), CBasicParticle::getRad(), CHertzMindlinViscoIGP::m_COR, CHertzMindlinViscoIGP::m_E, CHertzMindlinViscoIGP::m_nu, and CHertzMindlinViscoIGP::mu.

Here is the call graph for this function:

◆ ~CHertzMindlinViscoInteraction()

CHertzMindlinViscoInteraction::~CHertzMindlinViscoInteraction ( )

virtual

Member Function Documentation

◆ calcForces()

void CHertzMindlinViscoInteraction::calcForces ( )

virtual

Calculate viscoelastic and frictional forces.

Implements APairInteraction.

References CParticle::applyForce(), CParticle::getMass(), CBasicParticle::getPos(), CBasicParticle::getRad(), CParticle::getVel(), m_COR, m_cpos, m_dn, m_dt, m_E, m_E_diss, m_Ffric, m_force_deficit, m_is_slipping, m_is_touching, m_mu, m_normal_force, m_nu, APairInteraction::m_p1, APairInteraction::m_p2, Vec3::norm(), and Vec3::unit().

Here is the call graph for this function:

◆ Count()

double CHertzMindlinViscoInteraction::Count ( ) const

return 1 if particles are in contact, 0 otherwise

References m_is_touching.

Referenced by getScalarFieldFunction().

Here is the caller graph for this function:

◆ getAbsFN()

pair< bool, double > CHertzMindlinViscoInteraction::getAbsFN ( ) const

get current normal force

References m_is_touching, m_normal_force, and Vec3::norm().

Referenced by getCheckedScalarFieldFunction().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getAbsForceDeficit()

double CHertzMindlinViscoInteraction::getAbsForceDeficit ( ) const

get the force needed to overcome friction and make the interaction slip

References m_force_deficit, and Vec3::norm().

Here is the call graph for this function:

◆ getAbsFrictionalForce()

pair< bool, double > CHertzMindlinViscoInteraction::getAbsFrictionalForce ( ) const

get current frictional/stopping force

References CBasicParticle::getPos(), CBasicParticle::getRad(), m_Ffric, APairInteraction::m_p1, APairInteraction::m_p2, and Vec3::norm().

Referenced by getCheckedScalarFieldFunction().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getAbsFrictionalStress()

pair< bool, double > CHertzMindlinViscoInteraction::getAbsFrictionalStress ( ) const

get current frictional/stopping stress (f_fric/r^2)

References CBasicParticle::getPos(), CBasicParticle::getRad(), m_Ffric, APairInteraction::m_p1, APairInteraction::m_p2, and Vec3::norm().

Referenced by getCheckedScalarFieldFunction().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getAbsMuFN()

pair< bool, double > CHertzMindlinViscoInteraction::getAbsMuFN ( ) const

get max. frictional force, i.e. coeff. of friction * normal force

References CParticle::getMass(), CBasicParticle::getPos(), CBasicParticle::getRad(), CParticle::getVel(), m_COR, m_dn, m_E, APairInteraction::m_p1, APairInteraction::m_p2, Vec3::norm(), and Vec3::unit().

Referenced by getCheckedScalarFieldFunction().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getCheckedScalarFieldFunction()

CHertzMindlinViscoInteraction::CheckedScalarFieldFunction CHertzMindlinViscoInteraction::getCheckedScalarFieldFunction ( const string & name )

static

Get the particle member function which returns a checked scalar field of a given name.

Parameters

name	the name of the field

References getAbsFN(), getAbsFrictionalForce(), getAbsFrictionalStress(), getAbsMuFN(), getMaxFricStress(), getMuEffXY(), getMuEffXZ(), getNormalStress(), getSlipVelocity(), and NULL.

Here is the call graph for this function:

◆ getDissipatedEnergy()

double CHertzMindlinViscoInteraction::getDissipatedEnergy ( ) const

return the amount of energy dissipated during the last time step

References m_E_diss.

Referenced by getScalarFieldFunction().

Here is the caller graph for this function:

◆ getForce()

Vec3 CHertzMindlinViscoInteraction::getForce ( ) const

get net force on particle1 imposed by this interaction. Returns Vec3::ZERO if particles are not in contact.

References m_Ffric, m_is_touching, and m_normal_force.

Referenced by getVectorFieldFunction().

Here is the caller graph for this function:

◆ getMaxFricStress()

pair< bool, double > CHertzMindlinViscoInteraction::getMaxFricStress ( ) const

get max. frictional stress, i.e. coeff. of friction * normal stress

References CParticle::getMass(), CBasicParticle::getPos(), CBasicParticle::getRad(), CParticle::getVel(), m_COR, m_dn, m_E, APairInteraction::m_p1, APairInteraction::m_p2, Vec3::norm(), and Vec3::unit().

Referenced by getCheckedScalarFieldFunction().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getMuEff()

pair< bool, double > CHertzMindlinViscoInteraction::getMuEff	(	const Vec3 &	dir,
		const Vec3 &	norm
	)		const

Calculate effective coefficient of friction for this interaction for a given direction of the applied shear force. If the effective coefficient of friction is infinite, -1 is returned.

Parameters

dir	the direction of the applied shear force

Returns: the effective coefficient of friction if it is finite, -1 otherwise and -2 for no contact

References CBasicParticle::getPos(), CBasicParticle::getRad(), CParticle::getVel(), m_mu, APairInteraction::m_p1, APairInteraction::m_p2, Vec3::norm(), and Vec3::unit().

Referenced by getMuEffXY(), and getMuEffXZ().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getMuEffXY()

std::pair<bool,double> CHertzMindlinViscoInteraction::getMuEffXY ( ) const

inline

References getMuEff().

Referenced by getCheckedScalarFieldFunction().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getMuEffXZ()

std::pair<bool,double> CHertzMindlinViscoInteraction::getMuEffXZ ( ) const

inline

References getMuEff().

Referenced by getCheckedScalarFieldFunction().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getNormalForce()

Vec3 CHertzMindlinViscoInteraction::getNormalForce ( ) const

If the particles are in contact, get normal force, if not in contact return (0,0,0)

References m_is_touching, and m_normal_force.

Referenced by getVectorFieldFunction().

Here is the caller graph for this function:

◆ getNormalStress()

pair< bool, double > CHertzMindlinViscoInteraction::getNormalStress ( ) const

get current normal stress

References CBasicParticle::getRad(), m_is_touching, m_normal_force, APairInteraction::m_p1, APairInteraction::m_p2, and Vec3::norm().

Referenced by getCheckedScalarFieldFunction().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getPos()

virtual Vec3 CHertzMindlinViscoInteraction::getPos ( ) const

inlinevirtual

Implements APairInteraction.

References m_cpos.

◆ getPotentialEnergy()

double CHertzMindlinViscoInteraction::getPotentialEnergy ( ) const

get the potential energy stored in the interaction

References m_E, and m_normal_force.

Referenced by getScalarFieldFunction().

Here is the caller graph for this function:

◆ getScalarFieldFunction()

CHertzMindlinViscoInteraction::ScalarFieldFunction CHertzMindlinViscoInteraction::getScalarFieldFunction ( const string & name )

static

Get the particle member function which returns a scalar field of a given name.

Parameters

name	the name of the field

References Count(), getDissipatedEnergy(), getPotentialEnergy(), getSlipping(), getSticking(), and NULL.

Here is the call graph for this function:

◆ getSlipping()

double CHertzMindlinViscoInteraction::getSlipping ( ) const

Get the static/dynamic status of the interaction. Returns 1 for a contact in dynamic friction, 0 for static or no contact

References m_is_slipping.

Referenced by getScalarFieldFunction().

Here is the caller graph for this function:

◆ getSlipVelocity()

pair< bool, double > CHertzMindlinViscoInteraction::getSlipVelocity ( ) const

get "force deficit", i.e. the force needed to make the contact dynamic

get the slipping velocity, i.e. the absolute value of the tangential part of the relatve particle velocity

References CBasicParticle::getPos(), CBasicParticle::getRad(), CParticle::getVel(), APairInteraction::m_p1, APairInteraction::m_p2, Vec3::norm(), and Vec3::unit().

Referenced by getCheckedScalarFieldFunction().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getSticking()

double CHertzMindlinViscoInteraction::getSticking ( ) const

Get "sticking" contacts, i.e. return 1 if the contact is touching but not slipping, 0 otherwise

References m_is_slipping, and m_is_touching.

Referenced by getScalarFieldFunction().

Here is the caller graph for this function:

◆ getType()

static string CHertzMindlinViscoInteraction::getType ( )

inlinestatic

◆ getVectorFieldFunction()

CHertzMindlinViscoInteraction::VectorFieldFunction CHertzMindlinViscoInteraction::getVectorFieldFunction ( const string & name )

static

Get the particle member function which returns a vector field of a given name.

Parameters

name	the name of the field

References getForce(), getNormalForce(), and NULL.

Here is the call graph for this function:

◆ isPersistent()

bool CHertzMindlinViscoInteraction::isPersistent ( )

virtual

References CBasicParticle::getPos(), CBasicParticle::getRad(), APairInteraction::m_p1, and APairInteraction::m_p2.

Here is the call graph for this function:

◆ setTimeStepSize()

void CHertzMindlinViscoInteraction::setTimeStepSize ( double dt )

References m_dt.

Friends And Related Function Documentation

◆ operator<<

std::ostream& operator<<	(	std::ostream &	,
		const CHertzMindlinViscoInteraction &
	)

friend

◆ TML_PackedMessageInterface

friend class TML_PackedMessageInterface

friend

Member Data Documentation

◆ m_COR

double CHertzMindlinViscoInteraction::m_COR

protected

coefficient of restitution

Referenced by calcForces(), getAbsMuFN(), getMaxFricStress(), TML_PackedMessageInterface::pack(), and TML_PackedMessageInterface::unpack().

◆ m_cpos

Vec3 CHertzMindlinViscoInteraction::m_cpos

protected

contact position

Referenced by calcForces(), and getPos().

◆ m_dn

double CHertzMindlinViscoInteraction::m_dn

protected

caching displacement for E_pot

Referenced by calcForces(), getAbsMuFN(), getMaxFricStress(), TML_PackedMessageInterface::pack(), and TML_PackedMessageInterface::unpack().

◆ m_dt

double CHertzMindlinViscoInteraction::m_dt

protected

time step

Referenced by calcForces(), TML_PackedMessageInterface::pack(), setTimeStepSize(), and TML_PackedMessageInterface::unpack().

◆ m_E

double CHertzMindlinViscoInteraction::m_E

protected

Young's modulus.

Referenced by calcForces(), getAbsMuFN(), getMaxFricStress(), getPotentialEnergy(), TML_PackedMessageInterface::pack(), and TML_PackedMessageInterface::unpack().

◆ m_E_diss

double CHertzMindlinViscoInteraction::m_E_diss

protected

dissipated energy

Referenced by calcForces(), and getDissipatedEnergy().

◆ m_Ffric

Vec3 CHertzMindlinViscoInteraction::m_Ffric

protected

current frictional force

Referenced by calcForces(), getAbsFrictionalForce(), getAbsFrictionalStress(), getForce(), TML_PackedMessageInterface::pack(), and TML_PackedMessageInterface::unpack().

◆ m_force_deficit

Vec3 CHertzMindlinViscoInteraction::m_force_deficit

protected

difference between fric. force & force necessary for slip

Referenced by calcForces(), and getAbsForceDeficit().

◆ m_is_slipping

bool CHertzMindlinViscoInteraction::m_is_slipping

protected

static/dynamic status of the interaction

Referenced by calcForces(), getSlipping(), and getSticking().

◆ m_is_touching

bool CHertzMindlinViscoInteraction::m_is_touching

protected

contact status of the interaction

Referenced by calcForces(), Count(), getAbsFN(), getForce(), getNormalForce(), getNormalStress(), and getSticking().

◆ m_mu

double CHertzMindlinViscoInteraction::m_mu

protected

coefficient of friction

Referenced by calcForces(), getMuEff(), TML_PackedMessageInterface::pack(), and TML_PackedMessageInterface::unpack().

◆ m_normal_force

Vec3 CHertzMindlinViscoInteraction::m_normal_force

protected

current normal force

Referenced by calcForces(), getAbsFN(), getForce(), getNormalForce(), getNormalStress(), and getPotentialEnergy().

◆ m_nu

double CHertzMindlinViscoInteraction::m_nu

protected

Poisson ratio.

Referenced by calcForces(), TML_PackedMessageInterface::pack(), and TML_PackedMessageInterface::unpack().

◆ m_r0

double CHertzMindlinViscoInteraction::m_r0

protected

equilibrium distance

Referenced by TML_PackedMessageInterface::pack(), and TML_PackedMessageInterface::unpack().

The documentation for this class was generated from the following files:

Model/HertzMindlinViscoInteraction.h
Model/HertzMindlinViscoInteraction.cpp

Public Types

Public Member Functions

Static Public Member Functions

Protected Attributes

Friends

Detailed Description

Member Typedef Documentation

◆ CheckedScalarFieldFunction

◆ ParameterType

◆ ScalarFieldFunction

◆ VectorFieldFunction

Constructor & Destructor Documentation

◆ CHertzMindlinViscoInteraction() [1/3]

◆ CHertzMindlinViscoInteraction() [2/3]

◆ CHertzMindlinViscoInteraction() [3/3]

◆ ~CHertzMindlinViscoInteraction()

Member Function Documentation

◆ calcForces()

◆ Count()

◆ getAbsFN()

◆ getAbsForceDeficit()

◆ getAbsFrictionalForce()

◆ getAbsFrictionalStress()

◆ getAbsMuFN()

◆ getCheckedScalarFieldFunction()

◆ getDissipatedEnergy()

◆ getForce()

◆ getMaxFricStress()

◆ getMuEff()

◆ getMuEffXY()

◆ getMuEffXZ()

◆ getNormalForce()

◆ getNormalStress()

◆ getPos()

◆ getPotentialEnergy()

◆ getScalarFieldFunction()

◆ getSlipping()

◆ getSlipVelocity()

◆ getSticking()

◆ getType()

◆ getVectorFieldFunction()

◆ isPersistent()

◆ setTimeStepSize()

Friends And Related Function Documentation

◆ operator<<

◆ TML_PackedMessageInterface

Member Data Documentation

◆ m_COR

◆ m_cpos

◆ m_dn

◆ m_dt

◆ m_E

◆ m_E_diss

◆ m_Ffric

◆ m_force_deficit

◆ m_is_slipping

◆ m_is_touching

◆ m_mu

◆ m_normal_force

◆ m_nu

◆ m_r0