pressure_loaded_cylinder.cc File Reference

#include "generic.h"
#include "time_harmonic_fourier_decomposed_linear_elasticity.h"
#include "meshes/rectangular_quadmesh.h"
#include "meshes/triangle_mesh.h"

Classes
class	FourierDecomposedTimeHarmonicLinearElasticityProblem< ELEMENT >

Namespaces
	Global_Parameters
	Namespace for global parameters.

Functions
std::complex< double >	Global_Parameters::Nu (0.3, 0.05)
	Define Poisson's ratio Nu. More...
std::complex< double >	Global_Parameters::E (1.0, 0.01)
	Define the non-dimensional Young's modulus. More...
std::complex< double >	Global_Parameters::Omega_sq (10.0, 0.0)
const std::complex< double >	Global_Parameters::I (0.0, 1.0)
	Define the imaginary unit. More...
void	Global_Parameters::boundary_traction (const Vector< double > &x, const Vector< double > &n, Vector< std::complex< double > > &result)
	The traction function at r=rmin: (t_r, t_z, t_theta) More...
int	main (int argc, char *argv[])
	Driver code. More...

Function Documentation

main()

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

Driver code.

{
 // Number of elements in r-direction
 unsigned nr=10;
 
 // Number of elements in z-direction (for (approximately) square elements)
 unsigned nz=unsigned(double(nr)*Global_Parameters::Lz/Global_Parameters::Lr);
 
 // Set up doc info
 DocInfo doc_info;
 
 // Set output directory
 doc_info.set_directory("RESLT");
 
#ifdef ADAPTIVE
 
 // Set up problem
 FourierDecomposedTimeHarmonicLinearElasticityProblem
  <ProjectableTimeHarmonicFourierDecomposedLinearElasticityElement
   <TTimeHarmonicFourierDecomposedLinearElasticityElement<3> > >
  problem(nr,nz,Global_Parameters::rmin,Global_Parameters::rmax,
          Global_Parameters::zmin,Global_Parameters::zmax);
 
 // Solve
 unsigned max_adapt=3;
 problem.newton_solve(max_adapt);
 
#else
 
 // Set up problem
 FourierDecomposedTimeHarmonicLinearElasticityProblem
  <QTimeHarmonicFourierDecomposedLinearElasticityElement<3> > 
  problem(nr,nz,Global_Parameters::rmin,Global_Parameters::rmax,
          Global_Parameters::zmin,Global_Parameters::zmax);
 
 // Solve
 problem.newton_solve();
 
#endif
 
 // Output the solution
 problem.doc_solution(doc_info);
  
} // end_of_main

References Global_Parameters::Lr, Global_Parameters::Lz, Mesh_Parameters::nz, problem, Global_Parameters::rmax, Global_Parameters::rmin, oomph::DocInfo::set_directory(), Global_Parameters::zmax, and Global_Parameters::zmin.