TangentialSpringUnitTest.cpp File Reference

#include <iostream>
#include "Species/LinearViscoelasticSlidingFrictionSpecies.h"
#include "DPMBase.h"
#include "Particles/BaseParticle.h"
#include "Logger.h"

Go to the source code of this file.

Classes
class	TangentialSpringUnitTest
	In this file, the rolling behaviour of the tangential spring is tested. This is done by placing one normal partilce on top of a fixed partilce and letting graviry roll it over the other particle until it losses contact. More...

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

Variables
Logger< Log::ERROR >	testLogger ("Tangential Spring Unit Test")

Function Documentation

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

Definition at line 61 of file TangentialSpringUnitTest.cpp.

References BaseHandler< T >::begin(), BaseHandler< T >::copyAndAddObject(), FATAL, DPMBase::getTime(), logger, DPMBase::particleHandler, constants::pi, DPMBase::setGravity(), Files::setName(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), DPMBase::solve(), and DPMBase::speciesHandler.

{

    TangentialSpringUnitTest tangentialSpringUnitTestProblem;
    auto species = tangentialSpringUnitTestProblem.speciesHandler.copyAndAddObject(LinearViscoelasticSlidingFrictionSpecies());
    tangentialSpringUnitTestProblem.setName("TangentialSpringUnitTest");
    species->setDensity(6./constants::pi);
    species->setCollisionTimeAndRestitutionCoefficient(1e-3,.2,1./8);
    species->setSlidingStiffness(species->getStiffness()*2/7);
    species->setSlidingFrictionCoefficient(1e20);
    tangentialSpringUnitTestProblem.setTimeStep(2e-5);
    //tangentialSpringUnitTestProblem.setFileType(FileType::NO_FILE);
    tangentialSpringUnitTestProblem.setGravity(Vec3D(0,0,-1));
    tangentialSpringUnitTestProblem.setTimeMax(2.5);
    tangentialSpringUnitTestProblem.solve(argc,argv);
    
    std::vector<BaseParticle*>::iterator pIt = tangentialSpringUnitTestProblem.particleHandler.begin();
    
    if (!(*pIt)->getPosition().isEqualTo(Vec3D(0.5,0.5,0.25), 1e-7))  logger(FATAL, "First particles is in the wrong position. It is %",(*pIt)->getPosition());
    if (!(*pIt)->getVelocity().isEqualTo(Vec3D(0.0,0.0,0.0)  , 1e-7))  logger(FATAL, "First particle has the wrong velocity");
    ++pIt;

    if (!((*pIt)->getPosition().isEqualTo(Vec3D(0.961524052956078, 0.5, 0.450153103106219), 1e-7)))      logger(FATAL, "Second particle has the wrong position. It is % at time %",(*pIt)->getPosition(),tangentialSpringUnitTestProblem.getTime());
    if (!(*pIt)->getVelocity().isEqualTo(Vec3D(0.325869890236916, 0, -0.623251003973752)    , 1e-7))  logger(FATAL, "Second particle has the wrong velocity. It is %",(*pIt)->getVelocity());
    if (!(*pIt)->getAngularVelocity().isEqualTo(Vec3D(0.0,2.11896618998424,0.0)    , 1e-7))  logger(FATAL, "Second particles has the wrong angular velocity. It is %", (*pIt)->getAngularVelocity());

}

Variable Documentation

Logger<Log::ERROR> testLogger("Tangential Spring Unit Test")