unstructured_solid/unstructured_three_d_solid.cc File Reference

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

Classes
class	ElasticTetMesh< ELEMENT >
	Triangle-based mesh upgraded to become a solid mesh. More...
class	UnstructuredSolidProblem< ELEMENT, MESH >
	Unstructured solid problem. More...

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
void	Global_Physical_Variables::gravity (const double &time, const Vector< double > &xi, Vector< double > &b)
	Non-dimensional gravity as body force. More...
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...
int	main ()
	Demonstrate how to solve an unstructured solid problem. More...

Function Documentation

main()

int main

(

)

Demonstrate how to solve an unstructured solid problem.

{
 
 // Label for output
 DocInfo doc_info;
 
 // Output directory
 doc_info.set_directory("RESLT");
 
 // Create generalised Hookean constitutive equations
 Global_Physical_Variables::Constitutive_law_pt = 
  new GeneralisedHookean(&Global_Physical_Variables::Nu);
 
 //Set up the problem
 UnstructuredSolidProblem<TPVDElement<3,3> > problem;
 
 //Output initial configuration
 problem.doc_solution(doc_info);
 doc_info.number()++;
 
 // Parameter study
 Global_Physical_Variables::Gravity=12.0e-3;
 Global_Physical_Variables::P=0.0;
 double pressure_increment=-8.0e-3;
 
 unsigned nstep=2; // 10;
 for (unsigned istep=0;istep<nstep;istep++)
  {
   // Solve the problem
   problem.newton_solve();
   
   //Output solution
   problem.doc_solution(doc_info);
   doc_info.number()++;
   
   // Bump up suction
   Global_Physical_Variables::P+=pressure_increment;
  }
 
}

References Global_Physical_Variables::Constitutive_law_pt, Global_Physical_Variables::Gravity, Global_Physical_Variables::Nu, oomph::DocInfo::number(), Global_Physical_Variables::P, problem, and oomph::DocInfo::set_directory().