Tutorial2.cpp File Reference

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

Go to the source code of this file.

Classes
class	Tutorial2

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

Function Documentation

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

[T2:speciesProp]

[T2:speciesProp]

Definition at line 53 of file Tutorial2.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
  Tutorial2 problem;

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

  // The normal spring stiffness and normal dissipation is computed and set as 
  // For collision time tc=0.005 and restitution coefficeint rc=1.0, 
  auto species = problem.speciesHandler.copyAndAddObject(LinearViscoelasticSpecies());
  species->setDensity(2500.0); // sets the species type-0 density
  species->setStiffness(258.5);// sets the spring stiffness
  species->setDissipation(0.0);// 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(.005/50.0);// (collision time)/50.0
  problem.solve(argc, argv);
  
  return 0;
}