refineable_periodic_load.cc File Reference

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

Classes
class	RefineablePeriodicLoadProblem< ELEMENT >
	Periodic loading problem. More...

Namespaces
	Global_Parameters
	Namespace for global parameters.

Functions
TimeHarmonicIsotropicElasticityTensor	Global_Parameters::E (Nu)
	The elasticity tensor for the solid. More...
void	Global_Parameters::exact_solution (const Vector< double > &x, Vector< double > &u)
	The exact solution for infinite depth case. More...
void	Global_Parameters::periodic_traction (const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &result)
	The traction function. More...
int	main (int argc, char *argv[])
	Driver code for PeriodicLoad linearly elastic problem. More...

Function Documentation

main()

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

Driver code for PeriodicLoad linearly elastic problem.

{
 // Number of elements in x-direction
 unsigned nx=2; 
 
 // Number of elements in y-direction (for (approximately) square elements)
 unsigned ny=unsigned(double(nx)*Global_Parameters::Ly/Global_Parameters::Lx);
 
 // Set up doc info
 DocInfo doc_info;
 
 // Set output directory
 doc_info.set_directory("RESLT");
 
 // Set up problem
 RefineablePeriodicLoadProblem<RefineableQLinearElasticityElement<2,3> >
  problem(nx,ny,Global_Parameters::Lx, Global_Parameters::Ly);
 
 // Run the simulation, allowing for up to 4 spatial adaptations
 unsigned max_adapt=4; 
 problem.newton_solve(max_adapt);
 
 // Output the solution
 problem.doc_solution(doc_info);
  
} // end_of_main

References Global_Parameters::Lx, Global_Parameters::Ly, Mesh_Parameters::nx, Mesh_Parameters::ny, problem, and oomph::DocInfo::set_directory().