shock_disk.cc File Reference

#include <iostream>
#include <fstream>
#include <cmath>
#include "generic.h"
#include "solid.h"
#include "meshes/quarter_circle_sector_mesh.h"

Classes
class	ElasticRefineableQuarterCircleSectorMesh< ELEMENT >
class	DiskShockWaveProblem< ELEMENT, TIMESTEPPER >

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
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 (int argc, char *argv[])
	Driver for simple elastic problem. More...

Function Documentation

main()

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

Driver for simple elastic problem.

{
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 //Initialise physical parameters
 Global_Physical_Variables::E = 1.0; // ADJUST 
 Global_Physical_Variables::Nu = 0.3; // ADJUST
 
  // "Big G" Linear constitutive equations:
  Global_Physical_Variables::Constitutive_law_pt = 
   new GeneralisedHookean(&Global_Physical_Variables::Nu,
                          &Global_Physical_Variables::E);
 
  //Set up the problem:
  unsigned case_number=0;
 
 
 // Pure displacement formulation
  {
   cout << "Running case " << case_number 
        << ": Pure displacement formulation" << std::endl;
   DiskShockWaveProblem<RefineableQPVDElement<2,3>, Newmark<1> > problem;
   problem.run(case_number);
   case_number++;
  }
   
 // Pressure-displacement with Crouzeix Raviart-type pressure
  {
   cout << "Running case " << case_number 
        << ": Pressure/displacement with Crouzeix-Raviart pressure" << std::endl;
   DiskShockWaveProblem<RefineableQPVDElementWithPressure<2>, Newmark<1> >
    problem;
   problem.run(case_number);
   case_number++;
  }
 
 
  // Pressure-displacement with Taylor-Hood-type pressure
  {
   cout << "Running case " << case_number 
        << ": Pressure/displacement with Taylor-Hood pressure" << std::endl;
   DiskShockWaveProblem<RefineableQPVDElementWithContinuousPressure<2>, 
    Newmark<1> > problem;
   problem.run(case_number);
   case_number++;
  }
 
 
 // Clean up 
 delete Global_Physical_Variables::Constitutive_law_pt;
 Global_Physical_Variables::Constitutive_law_pt=0;
 
}

References Global_Physical_Variables::Constitutive_law_pt, Global_Physical_Variables::E, Global_Physical_Variables::Nu, problem, and Flag_definition::setup().