MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies > Class Template Reference

Contains contact force properties for contacts between particles with two different species. More...

#include <MixedSpecies.h>

Inheritance diagram for MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >:

Public Member Functions
	MixedSpecies ()
	The default constructor. More...
	MixedSpecies (const MixedSpecies &s)
	The default copy constructor. More...
	MixedSpecies (const Species< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies > &s)
	Creates a mixed species with the same force properties as a Species. More...
virtual	~MixedSpecies ()
	The default destructor. More...
MixedSpecies < NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies > *	copy () const final
	Creates a deep copy of the MixedSpecies from which it is called. More...
void	copyInto (BaseSpecies *bs) const final
	Copies the content of this into the species bs, if they are of the same type. More...
void	read (std::istream &is) final
	Reads the species properties from an input stream. More...
void	write (std::ostream &os) const final
	Writes the MixedSpecies properties to an output stream. More...
std::string	getName () const final
	Returns the name of the MixedSpecies as it is used in the restart file. More...
BaseInteraction *	getNewInteraction (BaseInteractable *const P, BaseInteractable *const I, unsigned timeStamp) const final
	When a contact between two particles is determined, an Interaction object is created, as the type of Interaction depends on the MixedSpecies type. More...
BaseInteraction *	getEmptyInteraction () const final
void	deleteEmptyInteraction (BaseInteraction *interaction) const final
bool	getUseAngularDOFs () const final
	Returns true if torques have to be calculated. More...
void	mixAll (BaseSpecies *const S, BaseSpecies *const T) final
	sets the MixedSpecies properties by mixing the properties of two particle species More...
Public Member Functions inherited from BaseSpecies
	BaseSpecies ()
	The default constructor. More...
	BaseSpecies (BaseNormalForce *normalForce, BaseFrictionForce *frictionForce_, BaseAdhesiveForce *adhesiveForce)
	BaseSpecies (const BaseSpecies &p)
	The copy constructor. More...
	~BaseSpecies ()
	The default destructor. More...
void	setHandler (SpeciesHandler *handler)
	Sets the pointer to the handler to which this species belongs. More...
SpeciesHandler *	getHandler () const
	Returns the pointer to the handler to which this species belongs. More...
Mdouble	getInteractionDistance () const
	returns the largest separation distance at which adhesive short-range forces can occur. More...
BaseNormalForce *	getNormalForce () const
BaseFrictionForce *	getFrictionForce () const
BaseAdhesiveForce *	getAdhesiveForce () const
void	setInteractionDistance (Mdouble interactionDistance)
	Sets BaseSpecies::interactionDistance_. More...
Public Member Functions inherited from BaseObject
	BaseObject ()=default
	Default constructor. More...
	BaseObject (const BaseObject &p)=default
	Copy constructor, copies all the objects BaseObject contains. More...
virtual	~BaseObject ()=default
	virtual destructor More...
virtual void	moveInHandler (unsigned int index)
	Except that it is virtual, it does the same thing as setIndex() does. More...
void	setIndex (unsigned int index)
	Allows one to assign an index to an object in the handler/container. More...
void	setId (unsigned long id)
	Assigns a unique identifier to each object in the handler (container) which remains constant even after the object is deleted from the container/handler. More...
unsigned int	getIndex () const
	Returns the index of the object in the handler. More...
unsigned int	getId () const
	Returns the unique identifier of any particular object. More...
void	setGroupId (unsigned groupId)
unsigned	getGroupId () const

Additional Inherited Members
Static Public Member Functions inherited from BaseSpecies
static Mdouble	average (Mdouble a, Mdouble b)
	Returns the harmonic mean of two variables. More...
static Mdouble	averageInf (Mdouble a, Mdouble b)
	Returns the harmonic mean of two variables, returning inf if either is inf. More...
Protected Attributes inherited from BaseSpecies
BaseNormalForce *	normalForce_
	A pointer to the normal force parameters This pointer is used by the Interaction's to get a pointer to the species The pointer is set in the constructors of SPecies and MixedSpecies. More...
BaseFrictionForce *	frictionForce_
	A pointer to the friction force parameters This pointer is used by the Interaction's to get a pointer to the species The pointer is set in the constructors of SPecies and MixedSpecies. More...
BaseAdhesiveForce *	adhesiveForce_
	A pointer to the adhesive force parameters This pointer is used by the Interaction's to get a pointer to the species The pointer is set in the constructors of SPecies and MixedSpecies. More...

Detailed Description

template<class NormalForceSpecies, class FrictionForceSpecies = EmptyFrictionSpecies, class AdhesiveForceSpecies = EmptyAdhesiveSpecies>
class MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >

Contains contact force properties for contacts between particles with two different species.

See Species for details.

Definition at line 42 of file MixedSpecies.h.

Constructor & Destructor Documentation

template<class NormalForceSpecies , class FrictionForceSpecies , class AdhesiveForceSpecies >

MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >::MixedSpecies

(

)

The default constructor.

Definition at line 102 of file MixedSpecies.h.

References DEBUG, and logger.

        : BaseSpecies(this,this,this), NormalForceSpecies(), FrictionForceSpecies(), AdhesiveForceSpecies()
{
    logger(DEBUG, "MixedSpecies::MixedSpecies() finished");
}

template<class NormalForceSpecies , class FrictionForceSpecies , class AdhesiveForceSpecies >

MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >::MixedSpecies

(

const MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies > &

s

)

The default copy constructor.

Definition at line 109 of file MixedSpecies.h.

References BaseSpecies::adhesiveForce_, DEBUG, BaseSpecies::frictionForce_, logger, BaseSpecies::normalForce_, and BaseForce::setBaseSpecies().

        : BaseSpecies(s), NormalForceSpecies(s), FrictionForceSpecies(s), AdhesiveForceSpecies(s)
{
    normalForce_ = this;
    frictionForce_ = this;
    adhesiveForce_ = this;
    normalForce_->setBaseSpecies(this);
    frictionForce_->setBaseSpecies(this);
    adhesiveForce_->setBaseSpecies(this);
    logger(DEBUG, "MixedSpecies::MixedSpecies(const MixedSpecies &p) finished");
}

template<class NormalForceSpecies , class FrictionForceSpecies , class AdhesiveForceSpecies >

MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >::MixedSpecies

(

const Species< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies > &

s

)

Creates a mixed species with the same force properties as a Species.

Definition at line 122 of file MixedSpecies.h.

References BaseSpecies::adhesiveForce_, DEBUG, BaseSpecies::frictionForce_, logger, BaseSpecies::normalForce_, and BaseForce::setBaseSpecies().

        : BaseSpecies(s), NormalForceSpecies(s), FrictionForceSpecies(s), AdhesiveForceSpecies(s)
{
    normalForce_ = this;
    frictionForce_ = this;
    adhesiveForce_ = this;
    normalForce_->setBaseSpecies(this);
    frictionForce_->setBaseSpecies(this);
    adhesiveForce_->setBaseSpecies(this);
    logger(DEBUG, "MixedSpecies::MixedSpecies(const MixedSpecies &p) finished");
}

template<class NormalForceSpecies , class FrictionForceSpecies , class AdhesiveForceSpecies >

MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >::~MixedSpecies ( )

virtual

The default destructor.

Definition at line 136 of file MixedSpecies.h.

References DEBUG, and logger.

{
    logger(DEBUG, "MixedSpecies::~MixedSpecies() finished");
}

Member Function Documentation

template<class NormalForceSpecies , class FrictionForceSpecies , class AdhesiveForceSpecies >

MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies > * MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >::copy ( ) const

finalvirtual

Creates a deep copy of the MixedSpecies from which it is called.

MixedSpecies copy method. It calls to copy constructor of this MixedSpecies, useful for polymorphism.

Implements BaseSpecies.

Definition at line 145 of file MixedSpecies.h.

Referenced by SpeciesHandler::readAndAddObject().

{
    return new MixedSpecies(*this);
}

template<class NormalForceSpecies , class FrictionForceSpecies , class AdhesiveForceSpecies >

void MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >::copyInto ( BaseSpecies *  bs ) const

finalvirtual

Copies the content of this into the species bs, if they are of the same type.

Useful for polymorphism: speciesHandler.getObject(i)->copyInto(bs); creates a deep copy (i.e. also copies properties of the derived species), whereas bs = speciesHandler.getObject(i); would only create a shallow copy.

Implements BaseSpecies.

Definition at line 159 of file MixedSpecies.h.

References getName(), logger, and WARN.

{
    if (bs == nullptr)
    {
        logger(WARN, "Error in %::copyInto: cannot copy into a nullptr");
        return;
    }
    auto* s = dynamic_cast<MixedSpecies<NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies>*>(bs);
    if (s == nullptr)
    {
        logger(WARN, "Error in %::copyInto: copying of % failed", getName(), s->getName());
        return;
    }
    *s = *this;
}

template<class NormalForceSpecies , class FrictionForceSpecies , class AdhesiveForceSpecies >

void MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >::deleteEmptyInteraction ( BaseInteraction *  interaction ) const

finalvirtual

Implements BaseSpecies.

Definition at line 263 of file MixedSpecies.h.

{
    Interaction<typename NormalForceSpecies::InteractionType, typename FrictionForceSpecies::InteractionType, typename AdhesiveForceSpecies::InteractionType>* interactionDestroyer;
    interactionDestroyer = dynamic_cast<Interaction<typename NormalForceSpecies::InteractionType, typename FrictionForceSpecies::InteractionType, typename AdhesiveForceSpecies::InteractionType>*>(interaction);
    delete interactionDestroyer;
}

template<class NormalForceSpecies , class FrictionForceSpecies , class AdhesiveForceSpecies >

BaseInteraction * MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >::getEmptyInteraction ( ) const

finalvirtual

Implements BaseSpecies.

Definition at line 257 of file MixedSpecies.h.

{
    return new Interaction<typename NormalForceSpecies::InteractionType, typename FrictionForceSpecies::InteractionType, typename AdhesiveForceSpecies::InteractionType>();
}

template<class NormalForceSpecies , class FrictionForceSpecies , class AdhesiveForceSpecies >

std::string MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >::getName ( ) const

finalvirtual

Returns the name of the MixedSpecies as it is used in the restart file.

Returns the name of the MixedSpecies as it is used in the restart file. The name of the species is a concatenation of the names of the three force components, e.g.

MixedSpecies<LinearViscoelasticNormalSpecies,SlidingFrictionSpecies,ReversibleAdhesiveSpecies> species; std::cout << species.getName();

will output "LinearViscoelasticSlidingFrictionReversibleAdhesiveMixedSpecies". The EmptyFrictionSpecies and the EmptyAdhesiveSpecies return empty strings, such that

MixedSpecies<LinearViscoelasticNormalSpecies> species; std::cout << species.getName();

will output "LinearViscoelasticMixedSpecies".

Returns

The name of the MixedSpecies.

Implements BaseObject.

Definition at line 231 of file MixedSpecies.h.

{
    return NormalForceSpecies::getBaseName()
           + FrictionForceSpecies::getBaseName()
           + AdhesiveForceSpecies::getBaseName() + "MixedSpecies";
}

template<class NormalForceSpecies , class FrictionForceSpecies , class AdhesiveForceSpecies >

BaseInteraction * MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >::getNewInteraction ( BaseInteractable *const  P, BaseInteractable *const  I, unsigned  timeStamp  ) const

finalvirtual

When a contact between two particles is determined, an Interaction object is created, as the type of Interaction depends on the MixedSpecies type.

The input parameters of this function are directly passed into the constructor for the new interaction. See Interaction for details.

Parameters

[in]	P	first of the two objects that interact
[in]	I	second of the two objects that interact
[in]	timeStamp	current value of DPMBase::time_

Returns

pointer to the newly created Interaction.

Implements BaseSpecies.

Definition at line 247 of file MixedSpecies.h.

{
    // JMFT: memory is allocated here, so it will need to be freed later
    return new Interaction<typename NormalForceSpecies::InteractionType, typename FrictionForceSpecies::InteractionType, typename AdhesiveForceSpecies::InteractionType>(
            P, I, timeStamp);
}

template<class NormalForceSpecies , class FrictionForceSpecies , class AdhesiveForceSpecies >

bool MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >::getUseAngularDOFs ( ) const

finalvirtual

Returns true if torques have to be calculated.

Returns true for any FrictionForceSpecies except EmptyFrictionSpecies, because for spherical particles, torques are only caused by tangential forces. See SpeciesHandler::useAngularDOFs for more details

Returns

true iff torques have to be calculated

Implements BaseSpecies.

Definition at line 278 of file MixedSpecies.h.

{
    return FrictionForceSpecies::getUseAngularDOFs();
}

template<class NormalForceSpecies , class FrictionForceSpecies , class AdhesiveForceSpecies >

void MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >::mixAll ( BaseSpecies *const S  , BaseSpecies *const T    )

finalvirtual

sets the MixedSpecies properties by mixing the properties of two particle species

Uses the harmonic mean for most properties. Calls the mix function for each of the force species from which MixedSpecies is derived.

Parameters

[in]	S	the first of two species whose properties are mixed to create the new species
[in]	T	the second of two species whose properties are mixed to create the new species

Implements BaseSpecies.

Definition at line 290 of file MixedSpecies.h.

References BaseNormalForce::getConstantRestitution(), BaseObject::getName(), BaseSpecies::getNormalForce(), and logger.

{
    logger.assert_always(T!= nullptr && S!= nullptr,"Arguments of mixAll cannot be null pointers");

    logger.assert_always(S->getNormalForce()->getConstantRestitution() == T->getNormalForce()->getConstantRestitution(), "mixing two LinearPlasticViscoelasticNormalSpecies, but only one has constantRestitution");
    NormalForceSpecies::setConstantRestitution(S->getNormalForce()->getConstantRestitution());

    const auto TN = dynamic_cast<NormalForceSpecies*> (T);
    const auto TF = dynamic_cast<FrictionForceSpecies*> (T);
    const auto TA = dynamic_cast<AdhesiveForceSpecies*> (T);
    logger.assert_always(TN!= nullptr && TF!= nullptr && TA!= nullptr,
            "Cannot mix two species of different type (% and %)",S->getName(),T->getName());

    const auto SN = dynamic_cast<NormalForceSpecies*> (S);
    const auto SF = dynamic_cast<FrictionForceSpecies*> (S);
    const auto SA = dynamic_cast<AdhesiveForceSpecies*> (S);
    logger.assert_always(SN!= nullptr && SF!= nullptr && SA!= nullptr,
            "Cannot mix two species of different type (% and %)",S->getName(),T->getName());

    NormalForceSpecies::mix(SN,TN);
    FrictionForceSpecies::mix(SF,TF);
    AdhesiveForceSpecies::mix(SA,TA);
}

template<class NormalForceSpecies , class FrictionForceSpecies , class AdhesiveForceSpecies >

void MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >::read ( std::istream &  is )

finalvirtual

Reads the species properties from an input stream.

Called by SpeciesHandler::readAndAddObject

Parameters

[in]

is

input stream (typically the restart file)

Reimplemented from BaseSpecies.

Definition at line 201 of file MixedSpecies.h.

References BaseSpecies::read(), BaseObject::setId(), and BaseObject::setIndex().

{
    //note: name is already read by SpeciesHandler::readAndAddObject
    std::string dummy;
    unsigned int id, index;
    is >> dummy >> id;
    is >> dummy >> index;
    BaseObject::setId(id);
    BaseObject::setIndex(index);
    BaseSpecies::read(is);
    NormalForceSpecies::read(is);
    FrictionForceSpecies::read(is);
    AdhesiveForceSpecies::read(is);
}

template<class NormalForceSpecies , class FrictionForceSpecies , class AdhesiveForceSpecies >

void MixedSpecies< NormalForceSpecies, FrictionForceSpecies, AdhesiveForceSpecies >::write ( std::ostream &  os ) const

finalvirtual

Writes the MixedSpecies properties to an output stream.

It prints human readable MixedSpecies information to the output stream, typically to Files::restartFile::fstream_. The basic species information is written in ParticleSpecies::write; then the three force types write additional information to the stream.

Parameters

[out]

os

output stream (typically the restart file)

Reimplemented from BaseSpecies.

Definition at line 184 of file MixedSpecies.h.

References BaseObject::getId(), BaseObject::getIndex(), getName(), and BaseSpecies::write().

{
    os << getName();
    os << " idA " << BaseObject::getId();
    os << " idB " << BaseObject::getIndex();
    BaseSpecies::write(os);
    NormalForceSpecies::write(os);
    FrictionForceSpecies::write(os);
    AdhesiveForceSpecies::write(os);
}

---

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

• Species/MixedSpecies.h