#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;
 }

Functions

Function Documentation

◆ main()