PlasticForceUnitTest.cpp File Reference

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

Go to the source code of this file.

Classes
class	PlasticForceUnitTest
	This code tests our plastic force model, as published in Luding 2008. More...

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

Function Documentation

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

Definition at line 82 of file PlasticForceUnitTest.cpp.

References BaseHandler< T >::begin(), FATAL, PlasticForceUnitTest::get_chi(), INFO, logger, DPMBase::particleHandler, PlasticForceUnitTest::set_chi(), Files::setName(), DPMBase::setParticleDimensions(), DPMBase::setSystemDimensions(), DPMBase::solve(), PlasticForceUnitTest::species, and DPMBase::writeRestartFile().

{   
    PlasticForceUnitTest PlasticForceUnitTestProblem;
    double k1=100.0;
    PlasticForceUnitTestProblem.species->setPlasticParameters(k1, 5.0*k1, k1, 0.05);
    PlasticForceUnitTestProblem.species->setDissipation(0);
    PlasticForceUnitTestProblem.setParticleDimensions(3);
    PlasticForceUnitTestProblem.setSystemDimensions(3);

    logger(INFO,"Testing particle particles collision for elastic plastic forces. \n"
        "This will be done for serveral values of scaled relative velocity chi");
    
    //PlasticForceUnitTestProblem.restartFile.getFstream().precision(20);
    
    //Set up constant data that will be used
    const std::vector<double> chi = {0.34, 0.69, 1.1, 1.37};
    const std::vector<Vec3D> leftFinalVecloity = {
        Vec3D(-0.032721738352012,0.0,0.0),
        Vec3D(-0.0138683231953154,0.0,0.0),
        Vec3D(-0.0204655358555405,0.0,0.0),
        Vec3D(-0.163049415300304,0.0,0.0)};
    const std::vector<Vec3D> leftFinalPosition = {
        Vec3D(0.995546292935715,1.0,1.0),
        Vec3D(1.00695193269955,1.0,1.0),
        Vec3D(1.00840467123501,1.0,1.0),
        Vec3D(0.969386085767181,1.0,1.0)};

    //Loop over all test cases
    for (int i=0; i<4; i++)
    {
        logger(INFO, "Running for chi=%",chi[i]);
        PlasticForceUnitTestProblem.set_chi(chi[i]);
        std::stringstream ss("");
        ss << "PlasticForceUnitTest" << PlasticForceUnitTestProblem.get_chi();
        PlasticForceUnitTestProblem.setName(ss.str().c_str());
        PlasticForceUnitTestProblem.solve();
        PlasticForceUnitTestProblem.writeRestartFile();

        //Now check the particles are in the right place for each of the 4 cases
        auto pIt = PlasticForceUnitTestProblem.particleHandler.begin();
        if (!(*pIt)->getPosition().isEqualTo(leftFinalPosition[i], 1e-10))
                logger(FATAL,"Left particle is in the wrong position. It is at % and should be %",(*pIt)->getPosition(),leftFinalPosition[i]);
        if (!(*pIt)->getVelocity().isEqualTo(leftFinalVecloity[i]  , 1e-10))
                logger(FATAL,"Left particle has the wrong velocity. It is at % and should be %",(*pIt)->getVelocity(),leftFinalVecloity[i]);
    }

    //to display use
    //gnuplot> plot 'plastic/collinear1.1.fstat' u 7:9 w lp
}