TangentialSpringUnitTest.cpp File Reference

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

Classes
class	TangentialSpringUnitTest

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

Function Documentation

main()

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

{
    TangentialSpringUnitTest tangentialSpringUnitTestProblem;
    tangentialSpringUnitTestProblem.setName("TangentialSpringUnitTest");
 
    //set the species properties, with a very large sliding friction coefficient
    auto species = tangentialSpringUnitTestProblem.speciesHandler.copyAndAddObject(
            LinearViscoelasticSlidingFrictionSpecies());
    species->setDensity(6. / constants::pi);
    species->setCollisionTimeAndRestitutionCoefficient(1e-3, .2, 1. / 8);
    species->setSlidingStiffness(species->getStiffness() * 2 / 7);
    species->setSlidingFrictionCoefficient(1e20);
 
    //set simulation-parameters: time, time step, gravity, output-file properties
    tangentialSpringUnitTestProblem.setTimeStep(2e-5);
    tangentialSpringUnitTestProblem.setFileType(FileType::NO_FILE);
    tangentialSpringUnitTestProblem.setGravity(Vec3D(0, 0, -1));
    tangentialSpringUnitTestProblem.setTimeMax(2.5);
 
    //solve the system
    tangentialSpringUnitTestProblem.solve(argc, argv);
 
    //Test if particles have the correct position and velocity
    BaseParticle* p = tangentialSpringUnitTestProblem.particleHandler.getObject(0);
    logger.assert_always(p->getPosition().isEqualTo(Vec3D(0.5, 0.5, 0.25), 1e-7),
                         "First particles is in the wrong position. It is %", p->getPosition());
    logger.assert_always(p->getVelocity().isEqualTo(Vec3D(0.0, 0.0, 0.0), 1e-7),
                         "First particle has the wrong velocity. It is %", p->getVelocity());
    
    p = tangentialSpringUnitTestProblem.particleHandler.getObject(1);
    logger.assert_always((p->getPosition().isEqualTo(Vec3D(0.961524052956078, 0.5, 0.450153103106219), 1e-7)),
                         "Second particle has the wrong position. It is % at time %", p->getPosition(),
                         tangentialSpringUnitTestProblem.getTime());
    logger.assert_always(p->getVelocity().isEqualTo(Vec3D(0.325869890236916, 0, -0.623251003973752), 1e-7),
                         "Second particle has the wrong velocity. It is %", p->getVelocity());
    logger.assert_always(p->getAngularVelocity().isEqualTo(Vec3D(0.0, 2.11896618998424, 0.0), 1e-7),
            "Second particles has the wrong angular velocity. It is %", p->getAngularVelocity());
}

References BaseHandler< T >::copyAndAddObject(), BaseInteractable::getAngularVelocity(), BaseHandler< T >::getObject(), BaseInteractable::getPosition(), DPMBase::getTime(), BaseInteractable::getVelocity(), Vec3D::isEqualTo(), logger, NO_FILE, DPMBase::particleHandler, constants::pi, DPMBase::setFileType(), DPMBase::setGravity(), DPMBase::setName(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), DPMBase::solve(), and DPMBase::speciesHandler.