Tutorial9_InclinedPlane.cpp File Reference

#include <Species/LinearViscoelasticFrictionSpecies.h>
#include <Mercury3D.h>
#include <Walls/InfiniteWall.h>

Classes
class	Tutorial9
	[T9:headers] More...

Functions
int	main (int argc, char *argv[])
	[T9:class] More...

Function Documentation

main()

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

[T9:class]

[T9:main]

{
    
    // Problem setup
    Tutorial9 problem;//instantiate an object of class Tutorial 9
    
    double angle = constants::pi / 180.0 * 20.0;
    
    problem.setName("Tutorial9");
    problem.setSystemDimensions(3);
    problem.setGravity(Vec3D(sin(angle), 0.0, -cos(angle)) * 9.81);
    problem.setXMax(0.3);
    problem.setYMax(0.3);
    problem.setZMax(0.05);
    problem.setTimeMax(0.5);
 
    // The normal spring stiffness and normal dissipation is computed and set as
    // For collision time tc=0.005 and restitution coefficeint rc=0.88,
 
    //Properties index 0
    LinearViscoelasticFrictionSpecies species0;
 
    species0.setDensity(2500.0);//sets the species type_0 density
    species0.setStiffness(259.018);//sets the spring stiffness.
    species0.setDissipation(0.0334);//sets the dissipation.
    species0.setSlidingStiffness(2.0 / 7.0 * species0.getStiffness());
    species0.setRollingStiffness(2.0 / 5.0 * species0.getStiffness());
    species0.setSlidingFrictionCoefficient(0.0);
    species0.setRollingFrictionCoefficient(0.0);
    auto ptrToSp0=problem.speciesHandler.copyAndAddObject(species0);
 
    //Properties index 1
    LinearViscoelasticFrictionSpecies species1;
 
    species1.setDensity(2500.0);//sets the species type-1 density
    species1.setStiffness(259.018);//sets the spring stiffness
    species1.setDissipation(0.0334);//sets the dissipation
    species1.setSlidingStiffness(2.0 / 7.0 * species1.getStiffness());
    species1.setRollingStiffness(2.0 / 5.0 * species1.getStiffness());
    species1.setSlidingFrictionCoefficient(0.5);
    species1.setRollingFrictionCoefficient(0.0);
    auto ptrToSp1=problem.speciesHandler.copyAndAddObject(species1);
 
    //Combination of properties index 0 and index 1
    auto species01 = problem.speciesHandler.getMixedObject(ptrToSp0,ptrToSp1);
 
    species01->setStiffness(259.018);//sets the spring stiffness
    species01->setDissipation(0.0334);//sets the dissipation
    species01->setSlidingStiffness(2.0 / 7.0 * species01->getStiffness());
    species01->setRollingStiffness(2.0 / 5.0 * species01->getStiffness());
    species01->setSlidingFrictionCoefficient(0.5);
    species01->setRollingFrictionCoefficient(0.0);
 
    //Properties index 2
    LinearViscoelasticFrictionSpecies species2;
 
    species2.setDensity(2500.0);//sets the species type-2 density
    species2.setStiffness(258.5);//sets the spring stiffness
    species2.setDissipation(0.0);//sets the dissipation
    species2.setSlidingStiffness(2.0 / 7.0 * species2.getStiffness());
    species2.setRollingStiffness(2.0 / 5.0 * species2.getStiffness());
    species2.setSlidingFrictionCoefficient(0.5);
    species2.setRollingFrictionCoefficient(0.5);
    auto ptrToSp2 = problem.speciesHandler.copyAndAddObject(species2);
 
    //Combination of properties index 0 and index 2
    auto species02 = problem.speciesHandler.getMixedObject(ptrToSp0, ptrToSp2);
 
    species02->setStiffness(259.018);//sets the stiffness
    species02->setDissipation(0.0334);//sets the dissipation
    species02->setSlidingStiffness(2.0 / 7.0 * species02->getStiffness());
    species02->setRollingStiffness(2.0 / 5.0 * species02->getStiffness());
    species02->setSlidingFrictionCoefficient(0.5);
    species02->setRollingFrictionCoefficient(0.5);
 
    //Output
    problem.setSaveCount(10);
    problem.dataFile.setFileType(FileType::ONE_FILE);
    problem.restartFile.setFileType(FileType::ONE_FILE);
    problem.fStatFile.setFileType(FileType::ONE_FILE);
    problem.eneFile.setFileType(FileType::NO_FILE);
    
    problem.setXBallsAdditionalArguments("-solidf -v0 -s 8");
    
    problem.setTimeStep(0.005 / 50.0);
    problem.solve(argc, argv);
    
    return 0;
}

References BaseHandler< T >::copyAndAddObject(), mathsFunc::cos(), DPMBase::dataFile, DPMBase::eneFile, DPMBase::fStatFile, SpeciesHandler::getMixedObject(), NO_FILE, ONE_FILE, constants::pi, DPMBase::restartFile, ParticleSpecies::setDensity(), File::setFileType(), DPMBase::setGravity(), DPMBase::setName(), DPMBase::setSaveCount(), DPMBase::setSystemDimensions(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), DPMBase::setXBallsAdditionalArguments(), DPMBase::setXMax(), DPMBase::setYMax(), DPMBase::setZMax(), mathsFunc::sin(), DPMBase::solve(), and DPMBase::speciesHandler.