Classes | Namespaces | Functions | Variables
cylinder.cc File Reference
#include <fenv.h>
#include "generic.h"
#include "axisym_linear_elasticity.h"
#include "meshes/rectangular_quadmesh.h"

Go to the source code of this file.

Classes

class  AxisymmetricLinearElasticityProblem< ELEMENT, TIMESTEPPER >
 Class to validate time harmonic linear elasticity (Fourier decomposed) More...
 

Namespaces

 Global_Parameters
 Namespace for global parameters.
 

Functions

void Global_Parameters::boundary_traction (const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &result)
 The traction function at r=Rmin: (t_r, t_z, t_theta) More...
 
void Global_Parameters::body_force (const double &time, const Vector< double > &x, Vector< double > &result)
 The body force function; returns vector of doubles in the order (b_r, b_z, b_theta) More...
 
void Global_Parameters::exact_solution_th (const Vector< double > &x, Vector< double > &u)
 Helper function - spatial components of the exact solution in a vector. This is necessary because we need to multiply this by different things to obtain the velocity and acceleration 0: u_r, 1: u_z, 2: u_theta. More...
 
double Global_Parameters::u_r (const double &time, const Vector< double > &x)
 Calculate the time dependent form of the r-component of displacement. More...
 
double Global_Parameters::u_z (const double &time, const Vector< double > &x)
 Calculate the time dependent form of the z-component of displacement. More...
 
double Global_Parameters::u_theta (const double &time, const Vector< double > &x)
 Calculate the time dependent form of the theta-component of displacement. More...
 
double Global_Parameters::d_u_r_dt (const double &time, const Vector< double > &x)
 Calculate the time dependent form of the r-component of velocity. More...
 
double Global_Parameters::d_u_z_dt (const double &time, const Vector< double > &x)
 Calculate the time dependent form of the z-component of velocity. More...
 
double Global_Parameters::d_u_theta_dt (const double &time, const Vector< double > &x)
 Calculate the time dependent form of the theta-component of velocity. More...
 
double Global_Parameters::d2_u_r_dt2 (const double &time, const Vector< double > &x)
 Calculate the time dependent form of the r-component of acceleration. More...
 
double Global_Parameters::d2_u_z_dt2 (const double &time, const Vector< double > &x)
 Calculate the time dependent form of the z-component of acceleration. More...
 
double Global_Parameters::d2_u_theta_dt2 (const double &time, const Vector< double > &x)
 Calculate the time dependent form of the theta-component of acceleration. More...
 
void Global_Parameters::exact_solution (const double &time, const Vector< double > &x, Vector< double > &u)
 The exact solution in a vector: 0: u_r, 1: u_z, 2: u_theta and their 1st and 2nd derivs. More...
 
int main (int argc, char *argv[])
 Driver code. More...
 

Variables

double Global_Parameters::Nu = 0.3
 Define Poisson's ratio Nu. More...
 
double Global_Parameters::E = 1.0
 Define the non-dimensional Young's modulus. More...
 
double Global_Parameters::Lambda = E*Nu/(1.0+Nu)/(1.0-2.0*Nu)
 Lame parameters. More...
 
double Global_Parameters::Mu = E/2.0/(1.0+Nu)
 
double Global_Parameters::Omega_sq = 0.5
 Square of the frequency of the time dependence. More...
 
unsigned Global_Parameters::Nr = 5
 Number of elements in r-direction. More...
 
unsigned Global_Parameters::Nz = 10
 Number of elements in z-direction. More...
 
double Global_Parameters::Lr = 1.0
 Length of domain in r direction. More...
 
double Global_Parameters::Lz = 2.0
 Length of domain in z-direction. More...
 
double Global_Parameters::Rmin = 0.1
 Set up min r coordinate. More...
 
double Global_Parameters::Zmin = 0.3
 Set up min z coordinate. More...
 
double Global_Parameters::Rmax = Rmin+Lr
 Set up max r coordinate. More...
 
double Global_Parameters::Zmax = Zmin+Lz
 Set up max z coordinate. More...
 

Function Documentation

◆ main()

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