Tutorial6.cpp File Reference

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

Go to the source code of this file.

Classes
class	Tutorial6
	[T6:headers] More...

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

Function Documentation

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

[T6:class]

[T6:speciesProp]

[T6:speciesProp]

Definition at line 69 of file Tutorial6.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
  Tutorial6 problem; // instantiate an object of class Tutorial 6

  problem.setName("Tutorial6");
  problem.setSystemDimensions(3);
  problem.setGravity(Vec3D(0.0,0.0,0.0));
  problem.setXMax(0.5);
  problem.setYMax(0.25);
  problem.setZMax(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(0.005/50.0);// (collision time)/50.0
  problem.solve(argc, argv);
  
  return 0;
}