FreeFallSelfTest.cpp File Reference

#include "Mercury2D.h"
#include "Particles/BaseParticle.h"
#include "Walls/InfiniteWall.h"
#include <iostream>
#include <Species/LinearViscoelasticSpecies.h>

Go to the source code of this file.

Classes
class	FreeFallSelfTest
	This case does a single elastic particle falling on an infinite plane. The k is chosen so that the maximum overlap with the wall is around 2% of the partcles dimater; whereas, the time is taken to ensure 50 steps with a collision. More...

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

Function Documentation

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

Start off my solving the default problem

Definition at line 54 of file FreeFallSelfTest.cpp.

References SpeciesHandler::addObject(), Files::fStatFile, NO_FILE, File::setFileType(), Files::setName(), DPMBase::setParticleDimensions(), Files::setSaveCount(), DPMBase::setTimeStep(), DPMBase::setXMax(), DPMBase::setYMax(), DPMBase::solve(), and DPMBase::speciesHandler.

{
    std::cout<<"Single particle bouncing vertically on the bottom plate"<<std::endl;
    FreeFallSelfTest freeFallSelfTestProblem;
    auto species = new LinearViscoelasticSpecies;
    freeFallSelfTestProblem.speciesHandler.addObject(species);
    species->setDensity(2000.0);
    freeFallSelfTestProblem.setParticleDimensions(3);
    species->setStiffness(8000000.0);

    freeFallSelfTestProblem.setName("FreeFallSelfTest");
    freeFallSelfTestProblem.setSaveCount(500);
    freeFallSelfTestProblem.fStatFile.setFileType(FileType::NO_FILE);
    freeFallSelfTestProblem.setTimeStep(1e-6);
    freeFallSelfTestProblem.setYMax(0.1);
    freeFallSelfTestProblem.setXMax(0.01);
    freeFallSelfTestProblem.solve(argc, argv);
}