TwoParticleClassicalWilletCollisionSelfTest.cpp File Reference

#include "Mercury2D.h"
#include "Walls/InfiniteWall.h"
#include "Species/Species.h"
#include "Species/HertzianViscoelasticMindlinLiquidBridgeClassicalWilletSpecies.h"

Classes
class	TwoParticleClassicalWilletCollision
	In this file two particles are symmetrically placed in a bi-axial box are allowed to jump around under gravity. It tests walls gravity and symmetry. More...

Functions
int	main (int argc UNUSED, char *argv[] UNUSED)

Function Documentation

main()

int main	(	int argc	UNUSED,
		char *argv[]	UNUSED
	)

{
 
    TwoParticleClassicalWilletCollision twoParticleClassicalWilletCollisionProblem;
    twoParticleClassicalWilletCollisionProblem.setName("TwoParticleClassicalWilletCollisionSelfTest");
 
    HertzianViscoelasticMindlinLiquidBridgeClassicalWilletSpecies species;
    species.setDensity(2000);
    species.setDissipation(0.16);
    species.setSlidingDissipation(0.16);
    Mdouble particleElasticModulus = 5e6; // taking 100 times that of soft material
    Mdouble particlePoissonRatio = 0.35;
    Mdouble particleSlidingFrictionCoefficient = 0.50;
    //Mdouble particleRadius = 0.0011;
    Mdouble particleShearModulus = particleElasticModulus/ 2 /(1+particlePoissonRatio);
    Mdouble particleEffectiveShearModulus = particleShearModulus/(2-particlePoissonRatio);
    Mdouble particleEffectiveElasticModulus = particleElasticModulus/(1-particlePoissonRatio*particlePoissonRatio);
    species.setEffectiveElasticModulus(particleEffectiveElasticModulus/2);
    species.setEffectiveShearModulus(particleEffectiveShearModulus/2);
    species.setSlidingFrictionCoefficient(particleSlidingFrictionCoefficient);    
    species.setLiquidBridgeVolume(8e-11);
    species.setContactAngle(0);
    species.setSurfaceTension(0.070);
 
    Mdouble radius = 0.0005;
    Mdouble rayleighTime = helpers::getRayleighTime(radius,species.getEffectiveShearModulus(), particlePoissonRatio,species.getDensity());
 
    twoParticleClassicalWilletCollisionProblem.setTimeStep(rayleighTime*0.1);
    twoParticleClassicalWilletCollisionProblem.speciesHandler.copyAndAddObject(species);
 
    twoParticleClassicalWilletCollisionProblem.setTimeMax(0.25);
    twoParticleClassicalWilletCollisionProblem.setSaveCount(twoParticleClassicalWilletCollisionProblem.getTimeMax()/twoParticleClassicalWilletCollisionProblem.getTimeStep());
    twoParticleClassicalWilletCollisionProblem.fStatFile.setFileType(FileType::MULTIPLE_FILES_PADDED);
    twoParticleClassicalWilletCollisionProblem.solve();
}

References BaseHandler< T >::copyAndAddObject(), e(), DPMBase::fStatFile, ParticleSpecies::getDensity(), helpers::getRayleighTime(), DPMBase::getTimeMax(), DPMBase::getTimeStep(), MULTIPLE_FILES_PADDED, UniformPSDSelfTest::radius, ParticleSpecies::setDensity(), File::setFileType(), DPMBase::setName(), DPMBase::setSaveCount(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), DPMBase::solve(), and DPMBase::speciesHandler.