elasticity.hpp

Go to the documentation of this file.

//==============================================================================
//==============================================================================
#ifndef ELASTICITY_HPP_
#define ELASTICITY_HPP_
 
#include <dune/common/fmatrix.hh>
 
namespace Opm {
namespace Elasticity {
 
  template<class GridType>
class Elasticity {
  public:
    static const int dim = GridType::dimension;
 
    typedef typename GridType::LeafGridView::ctype ctype;
 
    Elasticity(const GridType& gv_) : gv(gv_) {}
 
      template<int components, int funcdim>
    void getBmatrix(Dune::FieldMatrix<ctype,components,funcdim>& B,
                    const Dune::FieldVector<ctype,dim>& point,
                    const Dune::FieldMatrix<ctype,dim,dim>& Jinv);
 
      template<int comp, int funcdim>
    void getStiffnessMatrix(Dune::FieldMatrix<ctype,funcdim,funcdim>& A,
                            const Dune::FieldMatrix<ctype,comp,funcdim>& B,
                            const Dune::FieldMatrix<ctype,comp,comp>& C,
                            ctype detJW);
 
      template<int comp, int funcdim>
    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);
  protected:
    const GridType& gv;
};
 
Dune::FieldVector<double,3> waveSpeeds(const Dune::FieldMatrix<double,6,6>& C, double phi,
                                       double theta, double density);
 
}
}
 
#include "elasticity_impl.hpp"
 
#endif