Tutorial3_BouncingBallElastic.cpp File Reference

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

Classes
class	Tutorial3
	[T3:headers] More...

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

Function Documentation

main()

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

[T3:class]

[T3:main]

[T3:speciesProp]

[T3:speciesProp]

[T3:output]

[T3:output]

[T3:visualOutput]

[T3:visualOutput]

[T3:solve]

[T3:solve]

{
    
    // Problem setup
    Tutorial3 problem;
 
    // name the output files
    problem.setName("Tutorial3");
    // remove old output files with the same name
    problem.removeOldFiles();
    // set gravivational acceleration
    problem.setGravity(Vec3D(0.0, 0.0, -9.81));
    // set domain size (xMin = yMin = zMin = 0 by default)
    problem.setXMax(0.5);
    problem.setYMax(0.5);
    problem.setZMax(0.5);
    problem.setTimeMax(2.0);
 
    // Sets a linear spring-damper contact law.
    // The normal spring stiffness and normal dissipation is computed and set as
    // For collision time tc=0.005 and restitution coefficient rc=1.0,
    auto species = problem.speciesHandler.copyAndAddObject(LinearViscoelasticSpecies());
    species->setDensity(2500.0);
    double mass = species->getMassFromRadius(0.005);
    double collisionTime = 0.005;
    double restitution = 1.0;
    species->setCollisionTimeAndRestitutionCoefficient(collisionTime,restitution,mass);
    logger(INFO, "Stiffness %, dissipation %", species->getStiffness(), species->getDissipation());
 
    problem.setSaveCount(25);
 
    // Output vtk files for ParaView visualisation
    problem.wallHandler.setWriteVTK(FileType::ONE_FILE);
    problem.setParticlesWriteVTK(true);
 
    // Sets time step to 1/50th of the collision time
    problem.setTimeStep(0.005 / 50.0);
    // Start the solver (calls setupInitialConditions and initiates time loop)
    problem.solve(argc, argv);
    return 0;
}

References INFO, logger, ONE_FILE, and problem.