NurbsWallUnitTest.cpp File Reference

#include "Nurbs/NurbsSurface.h"
#include "Walls/NurbsWall.h"
#include "Math/Helpers.h"
#include <sstream>
#include <DPMBase.h>
#include <Interactions/NormalForceInteractions/LinearViscoelasticInteraction.h>
#include <Species/LinearViscoelasticSpecies.h>

Functions
int	main ()

Function Documentation

main()

int main

(

)

{
    //define quarter circle as nurbs surface
    std::vector<double> knotsU = {0,0,0,1,1,1};
    std::vector<double> knotsV = {0,0,1,1};
    std::vector<std::vector<Vec3D>> controlPoints = {{{1,0,0},{1,0,2}},{{1,1,0},{1,1,2}},{{0,1,0},{0,1,2}}} ;
    std::vector<std::vector<Mdouble>> weights = {{1,1},{1,1},{2,2}};
    NurbsSurface nurbsSurface(knotsU,knotsV,controlPoints,weights);
 
    //create a basic dpm class
    DPMBase dpm;
    dpm.setName("NurbsSurfaceUnitTest");
    auto s = dpm.speciesHandler.copyAndAddObject(LinearViscoelasticSpecies());
 
    //define nurbs wall and add to dpm class
    NurbsWall nurbsWall;
    nurbsWall.setSpecies(s);
    nurbsWall.set(nurbsSurface);
    auto w = dpm.wallHandler.copyAndAddObject(nurbsWall);
 
    SphericalParticle particle;
    particle.setSpecies(s);
    particle.setRadius(.99);
    for (double x=-.5; x<=4.5; x+=.125) {
        for (double z = -.5; z <= 2.5; z += .5) {
            particle.setPosition({x, 2-x, z});
            dpm.particleHandler.copyAndAddObject(particle);
        }
    }
 
    //test if contact can be found
    Vec3D normal;
    double distance;
    std::stringstream ss;
    for (auto p : dpm.particleHandler) {
        if (w->getDistanceAndNormal(*p, distance, normal)) {
            p->setPosition(distance * normal + p->getPosition());
            ss << p->getPosition() << '\t' << distance << '\t' << distance * normal + p->getPosition() << '\n';
        } else {
            p->setPosition({0,0,0});
        }
    }
    helpers::writeToFile("NurbsWallUnitTest.txt",ss.str());
 
    //write vtk file
    dpm.setName("NurbsWallUnitTest");
    dpm.wallHandler.setWriteVTK(FileType::ONE_FILE);
    dpm.setParticlesWriteVTK(true);
    dpm.forceWriteOutputFiles();
 
    return 0;
}

References BaseHandler< T >::copyAndAddObject(), DPMBase::forceWriteOutputFiles(), WallFunction::normal(), ONE_FILE, p, DPMBase::particleHandler, s, NurbsWall::set(), DPMBase::setName(), DPMBase::setParticlesWriteVTK(), BaseInteractable::setPosition(), BaseParticle::setRadius(), BaseParticle::setSpecies(), BaseWall::setSpecies(), WallHandler::setWriteVTK(), DPMBase::speciesHandler, w, DPMBase::wallHandler, calibrate::weights, helpers::writeToFile(), and plotDoE::x.