solid/airy_cantilever/airy_cantilever.cc File Reference

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

Classes
class	oomph::MySolidElement< ELEMENT >
	Wrapper class for solid element to modify the output. More...
class	oomph::FaceGeometry< MySolidElement< ELEMENT > >
	FaceGeometry of wrapped element is the same as the underlying element. More...
class	CantileverProblem< ELEMENT >
	Problem class for the cantilever "beam" structure. More...

Namespaces
	oomph
	DRAIG: Change all instances of (SPATIAL_DIM) to (DIM-1).
	Global_Physical_Variables
	Global variables.

Functions
void	Global_Physical_Variables::constant_pressure (const Vector< double > &xi, const Vector< double > &x, const Vector< double > &n, Vector< double > &traction)
	Constant pressure load. More...
void	Global_Physical_Variables::gravity (const double &time, const Vector< double > &xi, Vector< double > &b)
	Non-dimensional gravity as body force. More...
int	main ()

Function Documentation

main()

int main

(

)

Driver for cantilever beam loaded by surface traction and/or gravity

{
 
 // Create generalised Hookean constitutive equations
 Global_Physical_Variables::Constitutive_law_pt = 
  new GeneralisedHookean(&Global_Physical_Variables::Nu,
                         &Global_Physical_Variables::E);
 
 //Set up the problem
 CantileverProblem<MySolidElement<RefineableQPVDElement<2,3> > > problem;
 
 
 // Uncomment these as an exercise
 
 // CantileverProblem<MySolidElement<
 //  RefineableQPVDElementWithContinuousPressure<2> > > problem;
 
 // CantileverProblem<MySolidElement<
 //  RefineableQPVDElementWithPressure<2> > > problem;
 
 // Initial values for parameter values
 Global_Physical_Variables::P=0.0; 
 Global_Physical_Variables::Gravity=0.0;
 
 // Max. number of adaptations per solve
 unsigned max_adapt=3;
 
 //Parameter incrementation
 unsigned nstep=5; 
 double p_increment=1.0e-5;   
 for(unsigned i=0;i<nstep;i++)
  {
   // Increment pressure load
   Global_Physical_Variables::P+=p_increment;
 
   // 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();
 
  }
 
} //end of main

References Global_Physical_Variables::Constitutive_law_pt, Global_Physical_Variables::E, Global_Physical_Variables::Gravity, i, Global_Physical_Variables::Nu, Global_Physical_Variables::P, and problem.