Tutorial4.cpp File Reference

#include <Species/LinearViscoelasticSpecies.h>
#include <Mercury3D.h>
#include <Particles/BaseParticle.h>
#include <Walls/InfiniteWall.h>

Go to the source code of this file.

Classes
class	Tutorial4

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

Function Documentation

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

[T4:speciesProp]

[T4:speciesProp]

Definition at line 58 of file Tutorial4.cpp.

References BaseHandler< T >::copyAndAddObject(), Files::dataFile, Files::eneFile, Files::fStatFile, NO_FILE, ONE_FILE, Files::restartFile, File::setFileType(), DPMBase::setGravity(), Files::setName(), Files::setSaveCount(), DPMBase::setSystemDimensions(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), DPMBase::setXBallsAdditionalArguments(), DPMBase::setXMax(), DPMBase::setYMax(), DPMBase::setZMax(), DPMBase::solve(), and DPMBase::speciesHandler.

{

  // Problem setup
  Tutorial4 problem;

  problem.setName("Tutorial4");
  problem.setSystemDimensions(3);
  problem.setGravity(Vec3D(0.0,0.0,-9.81));
  problem.setXMax(1.0);
  problem.setYMax(1.0);
  problem.setZMax(2.0);
  problem.setTimeMax(5.0);

  // The normal spring stiffness and normal dissipation is computed and set as 
  // For collision time tc=0.005 and restitution coefficeint rc=0.88, 
  auto species = problem.speciesHandler.copyAndAddObject(LinearViscoelasticSpecies());
  species->setDensity(2500.0); // sets the species type-0 density
  species->setStiffness(259.0159);// sets the spring stiffness
  species->setDissipation(0.0334);// sets the dissipation

  problem.setSaveCount(10);
  problem.dataFile.setFileType(FileType::ONE_FILE);
  problem.restartFile.setFileType(FileType::ONE_FILE);
  problem.fStatFile.setFileType(FileType::NO_FILE);
  problem.eneFile.setFileType(FileType::NO_FILE);

  problem.setXBallsAdditionalArguments("-solidf -v0");  

  problem.setTimeStep(0.005/50.0); // (collision time)/50.0
  problem.solve(argc, argv);

  return 0;
}