Tutorial12_PrismWalls.cpp File Reference

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

Classes
class	Tutorial12
	[T12:headers] More...

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

Function Documentation

main()

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

[T12:class]

[T12:main]

[T12: constructor]

[T12:speciesProp]

[T12:speciesProp]

[T12: constructor]

[T12: Solve Problem]

[T12: Solve Problem]

{
    // Problem setup
    Tutorial12 problem; // instantiate an object of class Tutorial 12
    
    problem.setName("Tutorial12");
    problem.setGravity(Vec3D(0.0, 0.0, 0.0));
    problem.setXMax(0.5);
    problem.setYMax(0.5);
    problem.setZMax(0.1);
    problem.setTimeMax(5.0);
 
    // The normal spring stiffness and normal dissipation is computed and set as
    // For collision time tc=0.005 and restitution coefficient rc=1.0,
    LinearViscoelasticSpecies species;
    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.speciesHandler.copyAndAddObject(species);
 
    
    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.setParticlesWriteVTK(true);
    problem.wallHandler.setWriteVTK(FileType::ONE_FILE);
 
    problem.setTimeStep(.005 / 50.0); // (collision time)/50.0
    
    problem.solve(argc, argv);
    
    return 0;
}

References NO_FILE, ONE_FILE, problem, ParticleSpecies::setDensity(), LinearViscoelasticNormalSpecies::setDissipation(), and LinearViscoelasticNormalSpecies::setStiffness().