time_harmonic_fourier_decomposed_linear_elasticity/cylinder/cylinder.cc File Reference

#include "generic.h"
#include "time_harmonic_fourier_decomposed_linear_elasticity.h"
#include "meshes/rectangular_quadmesh.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, 5.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...
void	Global_Parameters::body_force (const Vector< double > &x, Vector< std::complex< double > > &b)
	Body force. More...
void	Global_Parameters::exact_solution (const Vector< double > &x, Vector< double > &u)
	The exact solution for infinite depth case. More...
int	main (int argc, char *argv[])
	Driver code. More...

Variables
std::complex< double >	Global_Parameters::lambda = E*Nu/(1.0+Nu)/(1.0-2.0*Nu)
std::complex< double >	Global_Parameters::mu = E/2.0/(1.0+Nu)
double	Global_Parameters::rmin = 0.1
double	Global_Parameters::zmin = 0.3
double	Global_Parameters::rmax = rmin+Lr
double	Global_Parameters::zmax = zmin+Lz

Function Documentation

main()

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

Driver code.

{
 // Number of elements in r-direction
 unsigned nr=5;
 
 // 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");
 
 // 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();
 
 // 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.