CoupledBeamDemo.cpp File Reference

#include "Species/LinearViscoelasticSlidingFrictionSpecies.h"
#include "Oomph/SCoupling/SCoupledSolidProblem.h"

Functions
int	main ()
	[CoupledBeamDemo] More...

Function Documentation

main()

int main

(

)

[CoupledBeamDemo]

{
    SCoupledSolidProblem<RefineableQDPVDElement<3, 2>> problem;
    //set name
    problem.setName("CoupledBeamUnitTest");
 
    //setup mercury
    problem.setDomain({0,0,0},{0.16,0.04,0.01});
    // add species
    auto species = problem.speciesHandler.copyAndAddObject(LinearViscoelasticSpecies());
    species->setDensity(1000);
    species->setStiffness(5.2);
    // add particle
    SphericalParticle p(species);
    p.setRadius(5e-3);
    p.setPosition({problem.getXCenter(),problem.getYCenter(),problem.getZMax()+p.getRadius()});
    p.setVelocity({0,0,-0.05});
    auto particle = problem.particleHandler.copyAndAddObject(p);
    // set time step
    double tc = species->getCollisionTime(species->getMassFromRadius(p.getRadius()));
    problem.setTimeStep(0.02*tc);
 
    // setup oomph
    problem.setElasticModulus(1e6);
    problem.setDensity(2500);
    problem.setSolidCubicMesh(7, 5, 5, 0, 0.16, 0, 0.04, 0, 0.01);
    problem.pinBoundaries({problem.Boundary::X_MIN, problem.Boundary::X_MAX});
    problem.setNewtonSolverTolerance(3e-8);
    problem.prepareForSolve();
    problem.coupleBoundary(problem.Boundary::Z_MAX);
    // sync time steps
    problem.setOomphTimeStep(problem.getTimeStep());
 
    // setup time
    problem.setTimeMax(100*problem.getTimeStep());
 
    // Solve the  problem and save mesh
    problem.solveSurfaceCoupling();
    problem.saveSolidMesh();
    return 0;
}

References e(), p, and problem.