prescribed_displ_lagr_mult2.cc File Reference

#include "generic.h"
#include "solid.h"
#include "constitutive.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	WarpedLine
	Warped line in 2D space. More...
class	PrescribedBoundaryDisplacementProblem< ELEMENT >

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main ()
	Driver code. More...

Variables
WarpedLine	Global_Physical_Variables::Boundary_geom_object (0.0)
	GeomObject specifying the shape of the boundary: Initially it's flat. More...

Function Documentation

main()

int main

(

)

Driver code.

{
 
 //Set up the problem
 PrescribedBoundaryDisplacementProblem<RefineableQPVDElement<2,3> > problem;
 
 // Doc initial domain shape
 problem.doc_solution();
 
 // Max. number of adaptations per solve
 unsigned max_adapt=1;
 
 //Parameter incrementation
 unsigned nstep=2; // 16;  
 for(unsigned i=0;i<nstep;i++)
  {
   // Increment imposed boundary displacement
   Global_Physical_Variables::Boundary_geom_object.ampl()+=0.0125;
 
   // Solve the problem with Newton's method, allowing
   // up to max_adapt mesh adaptations after every solve.
   problem.newton_solve(max_adapt);
   
   // Doc solution
   problem.doc_solution();
 
   // For maximum stability: Reset the current nodal positions to be
   // the "stress-free" ones -- this assignment means that the
   // parameter study no longer corresponds to a physical experiment
   // but is what we'd do if we wanted to use the solid solve
   // to update a fluid mesh in an FSI problem, say.
   problem.solid_mesh_pt()->set_lagrangian_nodal_coordinates();
   
  }
 
} //end of main

References WarpedLine::ampl(), Global_Physical_Variables::Boundary_geom_object, i, and problem.