Tutorial1.cpp File Reference

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

Go to the source code of this file.

Classes
class	Tutorial1
	[T1:headers] More...

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

Function Documentation

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

[T1:class]

[T1:main]

[T1:problemSetup]

[T1:problemSetup]

[T1:speciesProp]

[T1:speciesProp]

[T1:output]

[T1:output]

[T1:visualOutput]

[T1:visualOutput]

[T1:solve]

[T1:solve]

Definition at line 61 of file Tutorial1.cpp.

References BaseHandler< T >::copyAndAddObject(), Files::dataFile, Files::eneFile, Files::fStatFile, File::getCounter(), 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
  Tutorial1 problem;

  problem.setName("Tutorial1");
  problem.setSystemDimensions(3);
  problem.setGravity(Vec3D(0.0,0.0,0.0));
  problem.setXMax(1.0);
  problem.setYMax(1.0);
  problem.setZMax(1.0);
  problem.setTimeMax(2.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);
  std::cout << problem.dataFile.getCounter() << std::endl;

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

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