12#ifndef ELASTICITY_HPP_
13#define ELASTICITY_HPP_
15#include <dune/common/fmatrix.hh>
21 template<
class Gr
idType>
25 static const int dim = GridType::dimension;
28 typedef typename GridType::LeafGridView::ctype
ctype;
38 template<
int components,
int funcdim>
39 void getBmatrix(Dune::FieldMatrix<ctype,components,funcdim>& B,
40 const Dune::FieldVector<ctype,dim>& point,
41 const Dune::FieldMatrix<ctype,dim,dim>& Jinv);
48 template<
int comp,
int funcdim>
50 const Dune::FieldMatrix<ctype,comp,funcdim>& B,
51 const Dune::FieldMatrix<ctype,comp,comp>& C,
60 template<
int comp,
int funcdim>
62 const Dune::FieldVector<ctype,funcdim>& v,
63 const Dune::FieldVector<ctype,comp>& eps0,
64 const Dune::FieldMatrix<ctype,comp,funcdim>& B,
65 const Dune::FieldMatrix<ctype,comp,comp>& C);
76Dune::FieldVector<double,3>
waveSpeeds(
const Dune::FieldMatrix<double,6,6>& C,
double phi,
77 double theta,
double density);
Elasticity helper class.
Definition elasticity.hpp:22
Elasticity(const GridType &gv_)
Default constructor.
Definition elasticity.hpp:32
void getStressVector(Dune::FieldVector< ctype, comp > &sigma, const Dune::FieldVector< ctype, funcdim > &v, const Dune::FieldVector< ctype, comp > &eps0, const Dune::FieldMatrix< ctype, comp, funcdim > &B, const Dune::FieldMatrix< ctype, comp, comp > &C)
Return the stress vector in a quadrature point.
Definition elasticity_impl.hpp:92
void getStiffnessMatrix(Dune::FieldMatrix< ctype, funcdim, funcdim > &A, const Dune::FieldMatrix< ctype, comp, funcdim > &B, const Dune::FieldMatrix< ctype, comp, comp > &C, ctype detJW)
Return the stiffness matrix contributions in a quadrature point.
Definition elasticity_impl.hpp:76
GridType::LeafGridView::ctype ctype
A basic number.
Definition elasticity.hpp:28
const GridType & gv
Const reference to our grid.
Definition elasticity.hpp:68
static const int dim
The dimension of our grid.
Definition elasticity.hpp:25
void getBmatrix(Dune::FieldMatrix< ctype, components, funcdim > &B, const Dune::FieldVector< ctype, dim > &point, const Dune::FieldMatrix< ctype, dim, dim > &Jinv)
Returns the B matrix in a quadrature point.
Definition elasticity_impl.hpp:22
Dune::FieldVector< double, 3 > waveSpeeds(const Dune::FieldMatrix< double, 6, 6 > &C, double phi, double theta, double density)
Compute the elastic wave velocities.
Definition elasticity_impl.hpp:101
Elasticity helper class - template implementations.
Inverting small matrices.
Definition ImplicitAssembly.hpp:43