solid/three_d_cantilever/three_d_cantilever.cc File Reference

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

Classes
class	RefineableElasticQuarterTubeMesh< ELEMENT >
	Simple quarter tube mesh upgraded to become a solid mesh. More...
class	ElasticQuarterTubeMesh< ELEMENT >
	Simple quarter tube mesh upgraded to become a solid mesh. More...
class	CantileverProblem< ELEMENT >
	Problem class for the cantilever "beam" structure. 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...
int	main (int argc, char *argv[])
	Driver for 3D cantilever beam loaded by gravity. More...

Function Documentation

main()

int main	(	int argc	,
		char *	argv[]
	)

Driver for 3D cantilever beam loaded by gravity.

{
 
 // Run main demo code if no command line arguments are specified
 if (argc==1)
  {
   
   // Create incompressible Mooney Rivlin strain energy function
   Global_Physical_Variables::Strain_energy_function_pt = 
    new MooneyRivlin(&Global_Physical_Variables::C1,
                     &Global_Physical_Variables::C2);
   
   // 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 continous pressure/displacement
   CantileverProblem<RefineableQPVDElementWithContinuousPressure<3> > problem;
   
   // Doc solution
   problem.doc_solution();
   
   // Initial values for parameter values
   Global_Physical_Variables::Gravity=0.0;
   
   //Parameter incrementation
   unsigned nstep=10; 
   
   double g_increment=5.0e-4;   
   for(unsigned i=0;i<nstep;i++)
    {
     // Increment load
     Global_Physical_Variables::Gravity+=g_increment;
     
     // Solve the problem with Newton's method, allowing
     // up to max_adapt mesh adaptations after every solve.
     unsigned max_adapt=1;   
     problem.newton_solve(max_adapt);
     
     // Doc solution
     problem.doc_solution();    
    }
 
  } // end main demo code
 
 // Run self-tests
 else
  {
   
   // Additional test cases combining adaptive/non-adaptive
   // elements with displacement/displacement-pressure
   // discretisation and various constitutive equations
   
   
   // Initialise flag for FD evaluation
   bool use_fd=false; 
   
   // Number of cases per implementation
   unsigned ncase=5;
   
   // Is the material incomressible?
   bool incompress=true;
   
   // Loop over fd and analytical implementation
   for (unsigned i=0;i<2;i++) 
    {
     
     // Generalised Hookean constitutive equations
     //-------------------------------------------
     {
      Global_Physical_Variables::Constitutive_law_pt = 
       new GeneralisedHookean(&Global_Physical_Variables::Nu,
                              &Global_Physical_Variables::E);
      
      incompress=false;
      
#ifdef REFINE
      {
       //Set up the problem with pure displacement based elements
       CantileverProblem<RefineableQPVDElement<3,3> > problem; 
       problem.run_tests(0+i*ncase,incompress,use_fd);
      }
#else   
      {
       //Set up the problem with pure displacement based elements
       CantileverProblem<QPVDElement<3,3> > problem;
       problem.run_tests(0+i*ncase,incompress,use_fd);
      }
#endif
      
      
#ifdef REFINE
      {
       //Set up the problem with continous pressure/displacement
       CantileverProblem<RefineableQPVDElementWithContinuousPressure<3> > problem;
       problem.run_tests(1+i*ncase,incompress,use_fd);
      }
#else
      {
       //Set up the problem with continous pressure/displacement
       CantileverProblem<QPVDElementWithContinuousPressure<3> > problem;
       problem.run_tests(1+i*ncase,incompress,use_fd);
      }
#endif
      
      
#ifdef REFINE    
      {
       //Set up the problem with discontinous pressure/displacement
       CantileverProblem<RefineableQPVDElementWithPressure<3> > problem;
       problem.run_tests(2+i*ncase,incompress,use_fd);
     }
#else
      {
       //Set up the problem with discontinous pressure/displacement
       CantileverProblem<QPVDElementWithPressure<3> > problem;
       problem.run_tests(2+i*ncase,incompress,use_fd);
      }
#endif
      
      delete Global_Physical_Variables::Constitutive_law_pt;
      Global_Physical_Variables::Constitutive_law_pt=0;
     }
     
     
     
     // Incompressible Mooney Rivlin
     //-----------------------------
     {
      Global_Physical_Variables::Strain_energy_function_pt = 
       new MooneyRivlin(&Global_Physical_Variables::C1,
                        &Global_Physical_Variables::C2);
    
      // Define a constitutive law (based on strain energy function)
      Global_Physical_Variables::Constitutive_law_pt = 
       new IsotropicStrainEnergyFunctionConstitutiveLaw(
        Global_Physical_Variables::Strain_energy_function_pt);
    
      incompress=true;
 
 
#ifdef REFINE
      {
       //Set up the problem with continous pressure/displacement
       CantileverProblem<RefineableQPVDElementWithContinuousPressure<3> > problem;
       problem.run_tests(3+i*ncase,incompress,use_fd);
      }
#else
      {
       //Set up the problem with continous pressure/displacement
       CantileverProblem<QPVDElementWithContinuousPressure<3> > problem;     
       problem.run_tests(3+i*ncase,incompress,use_fd);
      }
#endif
 
 
    
#ifdef REFINE
      {
       //Set up the problem with discontinous pressure/displacement
       CantileverProblem<RefineableQPVDElementWithPressure<3> > problem;
       problem.run_tests(4+i*ncase,incompress,use_fd);
      }
#else
      {
       //Set up the problem with discontinous pressure/displacement
       CantileverProblem<QPVDElementWithPressure<3> > problem;
       problem.run_tests(4+i*ncase,incompress,use_fd);
      }
#endif
    
      delete  Global_Physical_Variables::Strain_energy_function_pt;
      Global_Physical_Variables::Strain_energy_function_pt=0;
    
      delete Global_Physical_Variables::Constitutive_law_pt;
      Global_Physical_Variables::Constitutive_law_pt=0;
     }
 
 
     use_fd=true;
     std::cout << "\n\n\n CR Total fill_in... : bla \n\n\n " << std::endl;
   
    }
  }
 
 
} //end of main

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