inclined_plane.cpp File Reference

#include "DPMBase.h"
#include "Particles/BaseParticle.h"
#include "Walls/InfiniteWall.h"
#include <iostream>
#include <Species/LinearViscoelasticSlidingFrictionSpecies.h>

Go to the source code of this file.

Classes
class	incl_plane

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 50 of file inclined_plane.cpp.

References BaseHandler< T >::copyAndAddObject(), BaseHandler< T >::getObject(), BaseInteractable::getVelocity(), DPMBase::particleHandler, constants::pi, DPMBase::setGravity(), Files::setName(), DPMBase::setParticleDimensions(), Files::setSaveCount(), DPMBase::setSystemDimensions(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), DPMBase::setXMax(), DPMBase::setYMax(), DPMBase::setZMax(), DPMBase::solve(), DPMBase::speciesHandler, and DPMBase::write().

{
    std::cout<<"Single particle rolling on a bottom plate"<<std::endl;
    incl_plane problem;
    
    auto species = problem.speciesHandler.copyAndAddObject(LinearViscoelasticSlidingFrictionSpecies());
    problem.setName("inclined_plane");
    problem.setSaveCount(5000);
    problem.setSystemDimensions(2);
    problem.setParticleDimensions(3);
    Mdouble theta = 15.0*constants::pi/180.0;
    problem.setGravity(9.8*Vec3D(std::sin(theta),-std::cos(theta),0.0));
    species->setDensity(2000);
    species->setStiffness(10000);
    species->setDissipation(0);
    species->setSlidingStiffness(10000);
    species->setSlidingDissipation(0);
    species->setSlidingFrictionCoefficient(0.5);
    //species->setRollingStiffness(10000);
    //species->setRollingDissipation(0);
    //species->setRollingFrictionCoefficient(0.5);
    //tc=7e-5
    problem.setTimeStep(2e-6);
    problem.setTimeMax(1.0);
    problem.setXMax(0.1);
    problem.setYMax(0.1);
    problem.setZMax(0.1);
    problem.solve(argc, argv);
    problem.write(std::cout, true);
    std::cout << problem.particleHandler.getObject(0)->getVelocity() << std::endl;

}