Domino.cpp File Reference

#include "Mercury3D.h"
#include "Walls/InfiniteWall.h"
#include "Species/LinearViscoelasticFrictionSpecies.h"
#include "Particles/ClumpParticle.h"
#include "../ClumpHeaders/ClumpInput.h"
#include "../ClumpHeaders/Mercury3DClump.h"
#include <stdlib.h>
#include "Boundaries/PeriodicBoundary.h"

Classes
struct	dominoes
class	multiParticleT1

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

Variables
dominoes	D

Function Documentation

main()

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

{
 
    int NumParams = 1;
 
    // Get parameters passed through the command line
    /*
    std::vector <Mdouble> params(0);
    for (int i = 0; i<NumParams; i++) {
        std::string a;
        Mdouble param = stod(helpers::readFromCommandLine(argc, argv, "-p"+std::to_string(i), a));
        params.push_back(param);
        std::cout<<params[i]<<std::endl;
    }
*/
    multiParticleT1 problem;
 
    // D.S_dom = 0.2 * (2 * D.k_peb * D.S_peb) + 0.8 * (2 * D.k_peb * D.S_peb) * (params[0]/100);
 
    // Domino species
    auto species0 = problem.speciesHandler.copyAndAddObject(LinearViscoelasticFrictionSpecies());
 
    // Wall species
    auto species1 = problem.speciesHandler.copyAndAddObject(LinearViscoelasticFrictionSpecies());
 
    // Mixed
    auto species01 = problem.speciesHandler.getMixedObject(species0, species1);
 
 
 
    species0->setDensity(1.0); // sets the species type-0 density
    //species->setConstantRestitution(0);
    species0->setSlidingFrictionCoefficient(0.0);
    species0->setSlidingStiffness(5e5);
    species0->setRollingFrictionCoefficient(0.0);
    species0->setRollingStiffness(5e5);
    species0->setDissipation(1.0);
    species0->setStiffness(1e6);
    const Mdouble collisionTime = species0->getCollisionTime(D.mass);
 
    species01->setSlidingFrictionCoefficient(0.2);
    species01->setSlidingStiffness(5e5);
    species01->setRollingFrictionCoefficient(0.1);
    species01->setDissipation(1.0);
    species01->setRollingStiffness(5e5);
    species01->setStiffness(1e6);
 
 
 
    problem.setClumpDamping(5);
    problem.setTimeStep(collisionTime / 50.0);
 
    // Quick demonstration
    problem.setSaveCount(500);
    problem.setTimeMax(10);
 
 
    // For time averaging featuring equipartition - takes a while
    // problem.setSaveCount(50000);
    // problem.setTimeMax(100000);
 
    problem.removeOldFiles();
    problem.solve();
 
    if ((D.started)&&(D.finished)) D.velocity = D.S_peb * (D.N_fin - D.N_ini)  / (D.T2-D.T1);
 
    std::cout<<"DOMINO VELOCITY: "<<D.velocity<<std::endl;
    std::cout<<"N_ini: "<<D.N_ini<<std::endl;
    std::cout<<"N_fin: "<<D.N_fin<<std::endl;
    std::cout<<"T1: "<<D.T1<<std::endl;
    std::cout<<"T2: "<<D.T2<<std::endl;
 
    // Add velocity to log
    std::ofstream funct("velocity.txt", std::ios_base::app | std::ios_base::out);
    funct << D.velocity <<"\n";
    funct.close();
 
 
 
    return 0;
}

References BaseHandler< T >::copyAndAddObject(), D, dominoes::finished, SpeciesHandler::getMixedObject(), dominoes::mass, dominoes::N_fin, dominoes::N_ini, DPMBase::removeOldFiles(), dominoes::S_peb, multiParticleT1::setClumpDamping(), DPMBase::setSaveCount(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), DPMBase::solve(), DPMBase::speciesHandler, dominoes::started, dominoes::T1, dominoes::T2, and dominoes::velocity.

Variable Documentation

D

dominoes D

Referenced by multiParticleT1::actionsAfterTimeStep(), DomainHandler::addObject(), IntersectionOfWalls::getDistanceAndNormal(), main(), multiParticleT1::multiParticleT1(), AxisymmetricHopper::setupInitialConditions(), multiParticleT1::setupInitialConditions(), test1(), and test2().