periodic_load.cc File Reference

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

Classes
class	PeriodicLoadProblem< 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...

Variables
double	Global_Parameters::Amplitude = 1.0
	Amplitude of traction applied. More...
bool	Global_Parameters::Finite =false
double	Global_Parameters::Lx = 1.0
	Length of domain in x direction. More...
double	Global_Parameters::Ly = 2.0
	Length of domain in y direction. 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=5;
 
 // 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
 PeriodicLoadProblem<QLinearElasticityElement<2,3> > 
  problem(nx,ny,Global_Parameters::Lx, Global_Parameters::Ly);
 
 // Solve
 problem.newton_solve();
 
 // 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().