#include "Mercury3D.h"
#include "Walls/InfiniteWall.h"
#include "Species/LinearViscoelasticFrictionSpecies.h"
#include "Particles/ClumpParticle.h"
#include "../ClumpHeaders/ClumpInput.h"
#include "../ClumpHeaders/Mercury3DClump.h"
#include <stdlib.h>

Classes
class	multiParticleT1

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

Variables
Mdouble	f_min = -2

Mdouble	f_max = 2

int	SAVECOUNT = 400

Function Documentation

◆ main()

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

 {
     multiParticleT1 problem;
     auto species = problem.speciesHandler.copyAndAddObject(LinearViscoelasticFrictionSpecies());
     species->setDensity(1.0); // sets the species type-0 density
     species->setDissipation(30.0);
     species->setStiffness(1e6);
  
     species->setSlidingFrictionCoefficient(0.6);
     species->setSlidingStiffness(5e5);
     species->setRollingFrictionCoefficient(0.0);
     species->setRollingStiffness(5e5);
  
     const Mdouble collisionTime = species->getCollisionTime(problem.getClumpMass());
     problem.setClumpDamping(0.01);
     problem.setTimeStep(collisionTime / 50.0);
     problem.setSaveCount(SAVECOUNT);
     problem.setHGridMaxLevels(1);
     problem.setTimeMax(30.0);
     problem.removeOldFiles();
     problem.solve();
     return 0;
 }