simple_shear.cc File Reference

#include "generic.h"
#include "solid.h"
#include "meshes/simple_cubic_mesh.template.h"

Classes
class	ElasticCubicMesh< ELEMENT >
	Simple cubic mesh upgraded to become a solid mesh. More...
class	SimpleShearProblem< ELEMENT >
	Boundary-driven elastic deformation of fish-shaped domain. More...

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
void	Global_Physical_Variables::body_force (const Vector< double > &xi, const double &t, Vector< double > &b)
	Body force vector: Vertically downwards with magnitude Gravity. More...
int	main ()
	Driver for simple elastic problem. More...

Function Documentation

main()

int main

(

)

Driver for simple elastic problem.

{
 //Initialise physical parameters
 Global_Physical_Variables::E  = 2.1; 
 Global_Physical_Variables::Nu = 0.4; 
 Global_Physical_Variables::C1 = 1.3; 
 
  for (unsigned i=0;i<2;i++)
  {
 
 // Define a strain energy function: Generalised Mooney Rivlin
 Global_Physical_Variables::Strain_energy_function_pt = 
  new GeneralisedMooneyRivlin(&Global_Physical_Variables::Nu,
                              &Global_Physical_Variables::C1,
                              &Global_Physical_Variables::E);
 
 // Define a constitutive law (based on strain energy function)
 Global_Physical_Variables::Constitutive_law_pt = 
  new IsotropicStrainEnergyFunctionConstitutiveLaw(
   Global_Physical_Variables::Strain_energy_function_pt);
 
 {
  //Set up the problem with pure displacement formulation
  SimpleShearProblem<QPVDElement<3,3> > problem(false);
  problem.run("RESLT");
 }
 
 //Discontinuous pressure
 {
  //Set up the problem with pure displacement formulation
  SimpleShearProblem<QPVDElementWithPressure<3> > problem(false);
  problem.run("RESLT_pres");
 }
 
 /*{
  //Set up the problem with pure displacement formulation
  SimpleShearProblem<QPVDElementWithPressure<3> > problem(true);
  problem.run("RESLT_pres_incomp");
  }*/
 
 
 {
  //Set up the problem with pure displacement formulation
  SimpleShearProblem<QPVDElementWithContinuousPressure<3> > problem(false);
  problem.run("RESLT_cont_pres");
 }
 
 /*{
  //Set up the problem with pure displacement formulation
  SimpleShearProblem<QPVDElementWithContinuousPressure<3> > problem(true);
  problem.run("RESLT_cont_pres_incomp");
  }*/
 
 
  }
 
 
}

References Global_Physical_Variables::C1, Global_Physical_Variables::Constitutive_law_pt, Global_Physical_Variables::E, i, Global_Physical_Variables::Nu, problem, and Global_Physical_Variables::Strain_energy_function_pt.