SinterBed2SelfTest.cpp File Reference

#include "SphericalIndenter.h"
#include "Species/SinterReversibleAdhesiveSpecies.h"
#include "Species/SinterFrictionSpecies.h"
#include "Walls/InfiniteWall.h"
#include "Boundaries/PeriodicBoundary.h"
#include "Walls/AxisymmetricIntersectionOfWalls.h"
#include "Logger.h"

Classes
class	Indenter
	Single particle, indented slowly by spherical indenter. More...

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

Function Documentation

main()

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

{
    Mdouble timeMax = 200;
    Mdouble indenterDiameter = 127e-6;
    Mdouble indentationForce = 4e-3; //4mN
    Mdouble indentationVelocity = readFromFile<Mdouble>("in","indentationVelocity",2e-6);
 
    //Ered 2e8
    //F=4e-3 = 4/3 E c del
    //F=4e-3 = 4/3 0.2e9 1e-6 1e-6  * 15
    //k=267
    Indenter i(indenterDiameter, indentationVelocity, indentationForce);
    if (i.getSystemDimensions()==2) {
        i.setIndentationForce(i.getIndentationForce()/3.0);
    }
    Mdouble frictionCoefficient = readFromFile<Mdouble>("in","frictionCoefficient",10);
    i.species->setSlidingFrictionCoefficient(frictionCoefficient);
    i.species->setRollingFrictionCoefficient(.1*frictionCoefficient);
    i.setTimeMax(timeMax);
 
    writeToFile("SinterBed2.gnu",
                         "set xlabel 'displacement [um]'\n"
                          "set ylabel 'force [mN]'\n"
                          "p 'SinterBed2.ene' u (-$2*1e6):($3*1e3) w lp\n"
    );
 
    i.setSaveCount(50);
    i.fStatFile.setFileType(FileType::NO_FILE);
    i.restartFile.setFileType(FileType::ONE_FILE);
    //i.setParticlesWriteVTK(true);
    i.solve();
    
    logger(INFO, "Execute 'gnuplot SinterBed2.gnu' to view output");
 
    return 0;
}

References e(), i, INFO, logger, NO_FILE, ONE_FILE, and helpers::writeToFile().