Tutorial8.cpp File Reference

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

Go to the source code of this file.

Classes
class	Tutorial8
	[T8:headers] More...

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

Function Documentation

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

[T8:class]

[T8:speciesProp]

[T8:speciesProp]

Definition at line 84 of file Tutorial8.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::solve(), and DPMBase::speciesHandler.

{

  // Problem setup
  Tutorial8 problem; // instantiate an object of class Tutorial 8
 
  problem.setName("Tutorial8");
  problem.setSystemDimensions(2);
  problem.setGravity(Vec3D(0.0,0.0,0.0));
  problem.setXMax(0.5);
  problem.setYMax(0.5);
  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=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 -s .85");   

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

  return 0;
}