5ParticlesRestartUnitTest.cpp File Reference

#include "DPMBase.h"
#include "Math/Helpers.h"

Go to the source code of this file.

Functions
int	main (int argc UNUSED, char *argv[] UNUSED)
	In this file, 5 Particles (4 fixed) are loaded from files "5Particles.ini" and "5Particles.restart". More...

Function Documentation

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

In this file, 5 Particles (4 fixed) are loaded from files "5Particles.ini" and "5Particles.restart".

The particles are aligned such that the single nonfixed particle rotates sinusoidally without moving. This is to test the behaviour of the tangential spring and the file loading routines.

Definition at line 34 of file 5ParticlesRestartUnitTest.cpp.

References DPMBase::readDataFile(), DPMBase::readRestartFile(), DPMBase::setAppend(), Files::setName(), DPMBase::solve(), and helpers::writeToFile().

{
    
    helpers::writeToFile("5Particles.ini",
        "5 0 0 0 0 1.6 1 1.6 \n"
        "0   0 0    0 0 0  0.5  0 0 0  0 0 0  0\n"
        "0   0 1.4  0 0 0  0.5  0 0 0  0 0 0  0\n"
        "1.4 0 0    0 0 0  0.5  0 0 0  0 0 0  0\n"
        "1.4 0 1.4  0 0 0  0.5  0 0 0  0 0 0  0\n"
        "0.7 0 0.7  0 0 0  0.5  0 0 0  0 30 0  0\n"
        );

    helpers::writeToFile("5Particles.restart",
        "restart_version 1.0 name 1Particle4WallsRestartUnitTest_restart\n"
        "dataFile name 1Particle4WallsRestartUnitTest_restart.data fileType ONE_FILE saveCount 5 counter 0 nextSavedTimeStep 0\n"
        "fStatFile name 1Particle4WallsRestartUnitTest_restart.fstat fileType ONE_FILE saveCount 5 counter 0 nextSavedTimeStep 0\n"
        "eneFile name 1Particle4WallsRestartUnitTest_restart.ene fileType ONE_FILE saveCount 1 counter 0 nextSavedTimeStep 0\n"
        "restartFile name 1Particle4WallsRestartUnitTest_restart.restart fileType ONE_FILE saveCount 5 counter 0 nextSavedTimeStep 0\n"
        "statFile name 1Particle4WallsRestartUnitTest_restart.stat fileType ONE_FILE saveCount 5 counter 0 nextSavedTimeStep 0\n"
        "xMin 0 xMax 1.6 yMin 0 yMax 1 zMin 0 zMax 1.6\n"
        "timeStep 1e-05 time 0 ntimeSteps 0 timeMax 0.01\n"
        "systemDimensions 3 particleDimensions 3 gravity 0 0 0\n"
        "Species 1\n"
        "LinearViscoelasticSlidingFrictionSpecies id 0 rho 1.9098593 k 200000 disp 0 slidingStiffness 57142.857 slidingDissipation 0 frictionCoefficient 0.5 frictionCoefficientStatic 0.5\n"
        "Walls 0\n"
        "WallsPeriodic 0\n"
        "Particles 5\n"
        "BP 0 0 0 0 0 0 0.5 0 0 0 0 0 0 0 0 0\n"
        "BP 0 0 0 0 0 0 0.5 0 0 0 0 0 0 0 0 0\n"
        "BP 0 0 0 0 0 0 0.5 0 0 0 0 0 0 0 0 0\n"
        "BP 0 0 0 0 0 0 0.5 0 0 0 0 0 0 0 0 0\n"
        "BP 0 0 0 0 0 0 0.5 0 0 0 0 0 0 0 0 0\n"
        "Interactions 0\n"
        );
    
    DPMBase problem;
    problem.setName("5Particles");
    problem.readRestartFile();
    problem.readDataFile("5Particles.ini");
    problem.setAppend(false);
    problem.solve();
    
    //Now check the rotational energy in the system
    //gnuplot: p '5Particles_restart.ene' u 1:4
}