clamped_or_pinned_shell.cc File Reference

#include "generic.h"
#include "shell.h"
#include "meshes/circular_shell_mesh.h"

Classes
class	ShellProblem< ELEMENT >

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
double &	Global_Physical_Variables::external_pressure ()
	Access function to value of external pressure. More...
void	Global_Physical_Variables::press_load (const Vector< double > &xi, const Vector< double > &x, const Vector< double > &N, Vector< double > &load)
	Load function: Perturbation pressure to force non-axisymmetric deformation. More...
int	main (int argc, char *argv[])
	Driver. More...

Function Documentation

main()

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

Driver.

{
 
 //Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 //Length of domain
 double L = 10.0;
 double L_phi=0.5*MathematicalConstants::Pi;
 
 // Choose problem type: 0: Clamped; 1: Pinned; 2: Pinned & symmetry
 unsigned problem_id=0;
 string dir_name;
 
 if (CommandLineArgs::Argc==2)
  {
   if (atoi(CommandLineArgs::Argv[1])==0)
    {
     std::cout << "Doing clamped case " << std::endl;
     problem_id=0;
     dir_name="RESLT_clamped";
    }
   else if (atoi(CommandLineArgs::Argv[1])==1)
    {
     std::cout << "Doing pinned case " << std::endl;
     problem_id=1;
     dir_name="RESLT_pinned";
    }
   else if (atoi(CommandLineArgs::Argv[1])==2)
    {
     std::cout << "Doing pinned periodic case " << std::endl;
     problem_id=2;
     dir_name="RESLT_pinned_periodic";
    }
   else
    {
     std::cout << "Wrong command line argument -- using default (0) instead" 
               << std::endl;
     problem_id=0;
     dir_name="RESLT_clamped";
    }
  }
 else
  {
   std::cout 
    << "No (or more than one) input argument(s) -- using default (0) instead" 
    << std::endl;  
     problem_id=0;
     dir_name="RESLT_clamped";  
  }
 
 
 // Run with pinned boundary conditions
 //------------------------------------
 if (problem_id==1)
 {
  //Set up the problem
  ShellProblem<StorableShapeSolidElement<DiagHermiteShellElement> >
   problem(5,3,L,L_phi,problem_id);
  
  //Solve the problem
  problem.solve(dir_name);
 }
 // Run with proper clamped boundary conditions
 //--------------------------------------------
 else if (problem_id==0)
  {
   {
    //Set up the problem
    ShellProblem<StorableShapeSolidElement<DiagHermiteShellElement> >
     problem(5,3,L,L_phi,problem_id);
   
    //Solve the problem
    problem.solve(dir_name);   
   }
//    {    
//     //Set up the problem
//     ShellProblem<StorableShapeSolidElement<DiagHermiteShellElement> >
//      problem(5,3,L,L_phi,problem_id);
    
//     //Solve the problem
//     problem.solve(dir_name);   
//    }
  }
 // Run with pinned periodic boundary conditions
 //---------------------------------------------
 else if (problem_id==2)
  {
   //Set up the problem
   ShellProblem<StorableShapeSolidElement<DiagHermiteShellElement> >
    problem(5,3,L,L_phi,problem_id);
   //Solve the problem
   problem.solve(dir_name);
   
  }
 else
  {
   std::cout  << "Wrong problem id " << problem_id << std::endl;
  }
 
 
 
}

References oomph::CommandLineArgs::Argc, oomph::CommandLineArgs::Argv, L, BiharmonicTestFunctions2::Pi, problem, and Flag_definition::setup().