Computes the normal forces due to linear plastic visco elastic interaction.
71 NormalForceInteraction::computeNormalForce();
72 Mdouble radius = 2.0 * NormalForceInteraction::getEffectiveRadius();
74 const SpeciesType* species =
static_cast<const SpeciesType*
>(NormalForceInteraction::getBaseSpecies()->getNormalForce());
78 if (pParticle && iParticle)
81 Mdouble heatTransfer = species->getThermalConductivity()
83 * contactArea / NormalForceInteraction::getDistance();
85 Mdouble mdT = heatTransfer / species->getHeatCapacity()
86 * NormalForceInteraction::getHandler()->getDPMBase()->getTimeStep();
ThermalSpecies< typename NormalForceInteraction::SpeciesType > SpeciesType
Definition: ThermalInteraction.h:18
Definition: ThermalParticle.h:15
void addTemperature(Mdouble temperature)
Definition: ThermalParticle.h:91
Mdouble getTemperature() const
Definition: ThermalParticle.h:80
#define max(a, b)
Definition: datatypes.h:23
const Mdouble pi
Definition: ExtendedMath.h:23