elasticity.hpp

Go to the documentation of this file.

//==============================================================================
//==============================================================================
#ifndef ELASTICITY_HPP_
#define ELASTICITY_HPP_

#include <dune/common/fmatrix.hh>
#include <dune/common/fvector.hh>

namespace Opm {
namespace Elasticity {

  template<class GridType>
class Elasticity {
  public:
    static const int dim = GridType::dimension;

    typedef typename GridType::LeafGridView::ctype ctype;

    explicit Elasticity(const GridType& gv_) : gv(gv_) {}

      template<int funcdim>
    void getBVector(Dune::FieldVector<ctype,funcdim>& BVector,
                    const Dune::FieldVector<ctype,dim>& point);

      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