compressed_square.cc File Reference

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

Classes
class	oomph::MySolidElement< ELEMENT >
	Wrapper class for solid element to modify the output. More...
class	CompressedSquareProblem< ELEMENT >
	Problem class. More...

Namespaces
	oomph
	DRAIG: Change all instances of (SPATIAL_DIM) to (DIM-1).
	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 ()
	Driver for compressed square. More...

Function Documentation

main()

int main

(

)

Driver for compressed square.

{
 
 //Flag to indicate if we want the material to be incompressible 
 bool incompress=false;
 
 // Label for different cases
 int case_number=-1;
 
 // Generalised Hookean constitutive equations
 //===========================================
 {
  // Nu values
  Vector<double> nu_value(3);
  nu_value[0]=0.45;
  nu_value[1]=0.499999;
  nu_value[2]=0.5;
  
  // Loop over nu values
  for (unsigned i=0;i<3;i++)
   {
    // Set Poisson ratio
    Global_Physical_Variables::Nu=nu_value[i];
  
    std::cout << "===========================================\n";
    std::cout << "Doing Nu=" <<  Global_Physical_Variables::Nu 
              << std::endl;
    std::cout << "===========================================\n";
 
    // Create constitutive equation
    Global_Physical_Variables::Constitutive_law_pt = 
     new GeneralisedHookean(&Global_Physical_Variables::Nu);
    
    // Displacement-based formulation
    //--------------------------------
    // (not for nu=1/2, obviously)
    if (i!=2) 
     {
      case_number++;
      std::cout 
       << "Doing Generalised Hookean with displacement formulation: Case "
       << case_number << std::endl;
      
      //Set up the problem with pure displacement based elements
      incompress=false;
      CompressedSquareProblem<MySolidElement<QPVDElement<2,3> > > 
       problem(incompress);
      
      // Run it
      problem.run_it(case_number,incompress);
     }
    
 
    // Generalised Hookean with continuous pressure, compressible
    //-----------------------------------------------------------
    {
     case_number++;
     std::cout 
      << "Doing Generalised Hookean with displacement/cont pressure "
      << "formulation: Case " << case_number << std::endl;     
     
     //Set up the problem with continous pressure/displacement
     incompress=false;
     CompressedSquareProblem<MySolidElement<
      QPVDElementWithContinuousPressure<2> > > 
      problem(incompress); 
     
     // Run it
     problem.run_it(case_number,incompress);
    }
     
    
    // Generalised Hookean with discontinuous pressure, compressible
    //--------------------------------------------------------------
    {
     case_number++;
     std::cout 
      << "Doing Generalised Hookean with displacement/discont pressure "
      << "formulation: Case " << case_number << std::endl; 
 
     //Set up the problem with discontinous pressure/displacement
     CompressedSquareProblem<MySolidElement<
      QPVDElementWithPressure<2> > > problem(incompress); 
     
     // Run it
     problem.run_it(case_number,incompress);
    }
 
 
 
    // Generalised Hookean with continous pressure/displacement; 
    //----------------------------------------------------------
    // incompressible
    //---------------
    {
     case_number++;
     std::cout 
      << "Doing Generalised Hookean with displacement/cont pressure, "
      << "incompressible formulation: Case " << case_number << std::endl;
     
     //Set up the problem with continous pressure/displacement
     incompress=true;
     CompressedSquareProblem<MySolidElement<
      QPVDElementWithContinuousPressure<2> > > 
      problem(incompress); 
     
     // Run it
     problem.run_it(case_number,incompress);
    }
 
 
    // Generalised Hookean with discontinous pressure/displacement; 
    //-------------------------------------------------------------
    // incompressible
    //---------------
    {
     case_number++;
     std::cout 
      << "Doing Generalised Hookean with displacement/discont pressure, "
      << "incompressible formulation: Case " << case_number << std::endl; 
     
     //Set up the problem with discontinous pressure/displacement
     incompress=true;
     CompressedSquareProblem<MySolidElement<
      QPVDElementWithPressure<2> > > problem(incompress); 
     
     // Run it
     problem.run_it(case_number,incompress);
    }
 
    // Clean up
    delete Global_Physical_Variables::Constitutive_law_pt;
    Global_Physical_Variables::Constitutive_law_pt=0;
    
   } 
 
 } // end generalised Hookean
 
 
 // Incompressible Mooney Rivlin
 //=============================
 {
 
  // Create 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);
    
 
  // Mooney-Rivlin with continous pressure/displacement; 
  //----------------------------------------------------
  // incompressible
  //---------------
  {
   case_number++;
   std::cout 
    << "Doing Mooney Rivlin with cont pressure formulation: Case "
    << case_number << std::endl; 
   
   //Set up the problem with continous pressure/displacement
   incompress=true;
   CompressedSquareProblem<MySolidElement<
    QPVDElementWithContinuousPressure<2> > > 
    problem(incompress); 
     
   // Run it
   problem.run_it(case_number,incompress);
  }
    
 
  // Mooney-Rivlin with discontinous pressure/displacement; 
  //-------------------------------------------------------
  // incompressible
  //---------------
  {
   case_number++;
   std::cout 
    << "Doing Mooney Rivlin with discont pressure formulation: Case "
    << case_number << std::endl; 
   
   //Set up the problem with discontinous pressure/displacement
   incompress=true;
   CompressedSquareProblem<MySolidElement<
    QPVDElementWithPressure<2> > > 
    problem(incompress); 
     
   // Run it
   problem.run_it(case_number,incompress);
  }
 
 
  // Clean up
  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;
  
 }
   
} //end of main

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