1Particle4WallsRestartUnitTest.cpp File Reference

#include <Logger.h>
#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, 1 Particle and 4 Walls are loaded from files "1Particle4Walls.ini" and "1Particle4Walls.restart". More...

Function Documentation

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

In this file, 1 Particle and 4 Walls are loaded from files "1Particle4Walls.ini" and "1Particle4Walls.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 35 of file 1Particle4WallsRestartUnitTest.cpp.

References FATAL, BaseInteractable::getAngularVelocity(), BaseHandler< T >::getObject(), BaseInteractable::getVelocity(), Vec3D::isEqualTo(), logger, DPMBase::particleHandler, DPMBase::readRestartFile(), DPMBase::setAppend(), Files::setName(), DPMBase::solve(), DPMBase::writeRestartFile(), and helpers::writeToFile().

{
    helpers::writeToFile("1Particle4WallsRestartUnitTest.ini",
        "1 0 0 0 0 .96 1 .96\n"
        "0.48 0 0.48  0 0 0  0.5  0 0 0  0 0 0  0\n"
        );

    helpers::writeToFile("1Particle4WallsRestartUnitTest.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 5 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 yMin 0 yMax 1 zMin 0 zMax 1\n"
        "timeStep 1e-04 time 0 ntimeSteps 0 timeMax 0.01\n"
        "systemDimensions 3 particleDimensions 3 gravity 0 0 0\n"
        "Species 1\n"
        "LinearViscoelasticSlidingFrictionSpecies id 0 density 1.9098593 stiffness 200000 dissipation 0 slidingStiffness 57142.857 slidingDissipation 0 frictionCoefficient 0.5 frictionCoefficientStatic 0.5\n"
        "Walls 4\n"
        "InfiniteWall id 0 indSpecies 0 position 0    0 0    orientation 0 0 0 1 velocity 0 0 0 angularVelocity 0 0 0 0 force 0 0 0 torque 0 0 0 normal -1 0 0 factor 1\n"
        "InfiniteWall id 1 indSpecies 0 position 0.96 0 0    orientation 0 0 0 1 velocity 0 0 0 angularVelocity 0 0 0 0 force 0 0 0 torque 0 0 0 normal  1 0 0 factor 1\n"
        "InfiniteWall id 2 indSpecies 0 position 0    0 0    orientation 0 0 0 1 velocity 0 0 0 angularVelocity 0 0 0 0 force 0 0 0 torque 0 0 0 normal 0 0 -1 factor 1\n"
        "InfiniteWall id 3 indSpecies 0 position 0    0 0.96 orientation 0 0 0 1 velocity 0 0 0 angularVelocity 0 0 0 0 force 0 0 0 torque 0 0 0 normal 0 0  1 factor 1\n"
        "Boundaries 0\n"
        "Particles 1\n"
        "BaseParticle id 0 indSpecies 0 position 0.48 0 0.48 orientation 0 0 0 1 velocity 0 0 0 angularVelocity 0 60 0 0 force 0 0 0 torque 0 0 0 radius 0.5 invMass 1 invInertia 10\n"
        "Interactions 4\n"
        "LinearViscoelasticSlidingFrictionInteraction particleWallIds 0 0 timeStamp 0 force 0 0 0 torque 0 0 0 slidingSpring 0 0 0\n"
        "LinearViscoelasticSlidingFrictionInteraction particleWallIds 0 1 timeStamp 0 force 0 0 0 torque 0 0 0 slidingSpring 0 0 0\n"
        "LinearViscoelasticSlidingFrictionInteraction particleWallIds 0 2 timeStamp 0 force 0 0 0 torque 0 0 0 slidingSpring 0 0 0\n"
        "LinearViscoelasticSlidingFrictionInteraction particleWallIds 0 3 timeStamp 0 force 0 0 0 torque 0 0 0 slidingSpring 0 0 0\n"
        );
        
    DPMBase problem;
    problem.setName("1Particle4WallsRestartUnitTest");
    problem.readRestartFile();
    //problem.readDataFile("1Particle4WallsRestartUnitTest.ini");
    problem.setAppend(false);
    //problem.write(std::cout,true);
    problem.solve();
    //problem.write(std::cout,true);
    problem.writeRestartFile();
    
    BaseParticle* p = problem.particleHandler.getObject(0);
    if (!p->getAngularVelocity().isEqualTo(Vec3D(0.0,17.7251802668803,0.0), 1e-7))  logger(FATAL, "First particles has the wrong angular velocity it is %, while is should be %.",p->getAngularVelocity(),Vec3D(0.0,17.7251802668803,0.0));
    if (!p->getVelocity().isEqualTo(Vec3D(0.0,0.0,0.0)  , 1e-7))  logger(FATAL, "First particle has the wrong velocity, it is %, while it should be %",p->getVelocity(),0.0);

    //Now check the rotational energy in the system; it should loose some rotational energy in the first loop due to
    //sliding, then conserve energy as no tangential dissipation is used.
    //gnuplot:
    // p '1Particle4WallsRestartUnitTest_restart.ene' u 1:4 w l t 'rotational energy', '1Particle4WallsRestartUnitTest_restart.ene' u 1:($3+$4+$5-160) w l t 'total energy - normal elastic energy'

}