Computes the normal forces due to linear plastic visco elastic interaction.
91 NormalForceInteraction::computeNormalForce();
92 Mdouble radius = 2.0 * NormalForceInteraction::getEffectiveRadius();
93 Mdouble contactArea =
constants::pi * radius * std::max(0.0,NormalForceInteraction::getOverlap());
94 const SpeciesType* species =
static_cast<const SpeciesType*
>(NormalForceInteraction::getBaseSpecies()->getNormalForce());
98 if (pParticle && iParticle)
101 Mdouble heatTransfer = species->getThermalConductivity()
103 * contactArea / NormalForceInteraction::getDistance();
105 Mdouble mdT = heatTransfer / species->getHeatCapacity()
106 * NormalForceInteraction::getHandler()->getDPMBase()->getTimeStep();
double Mdouble
Definition: GeneralDefine.h:34
ThermalSpecies< typename NormalForceInteraction::SpeciesType > SpeciesType
Definition: ThermalInteraction.h:39
Definition: ThermalParticle.h:36
void addTemperature(Mdouble temperature)
Definition: ThermalParticle.h:112
Mdouble getTemperature() const
Definition: ThermalParticle.h:101
const Mdouble pi
Definition: ExtendedMath.h:45