materials.hh

Go to the documentation of this file.

//==============================================================================
//==============================================================================
#ifndef MATERIALS_HH_
#define MATERIALS_HH_
 
#include "material.hh"
 
namespace Opm::Elasticity {
 
class Isotropic : public Material
{
public:
  Isotropic(int ID, double Emod, double Poisson, double rho_= 0.0)
    : Material(ID,rho_)
  {
    E = Emod;
    nu = Poisson;
  }
 
  ~Isotropic() override {}
 
  int numPar() const override
  {
    return 2;
  }
 
  double getPar(int ipar = 1) const override
  {
    return ipar == 1 ? E : nu;
  }
 
  void setE(double E_)
  {
    E = E_;
  }
 
  double getE() const
  {
    return E;
  }
 
  bool getConstitutiveMatrix(Dune::FieldMatrix<double,6,6>& C,
                             bool invers = false) const override;
 
  bool getConstitutiveMatrix(Dune::FieldMatrix<double,3,3>& C,
                             bool invers = false) const override;
 
protected:
  std::ostream& write(std::ostream& os) const override;
 
private:
  double E;  
  double nu; 
};
 
class OrthotropicD : public Material
{
public:
  OrthotropicD(int ID, double Ex, double Ey, double Ez,
               double Gxy, double Gxz = double(-1), double Gyz = double(-1));
 
  ~OrthotropicD() override {}
 
  int numPar() const override
  {
    return 6;
  }
 
  double getPar(int ipar = 1) const override;
 
  bool getConstitutiveMatrix(Dune::FieldMatrix<double,6,6>& C,
                             bool invers = false) const override;
 
  bool getConstitutiveMatrix(Dune::FieldMatrix<double,3,3>& C,
                             bool invers = false) const override;
 
protected:
  std::ostream& write(std::ostream& os) const override;
 
private:
  double E[6]; 
};
 
class OrthotropicSym : public Material
{
public:
  OrthotropicSym(int ID, const Dune::DynamicVector<double>& Cu);
 
  ~OrthotropicSym() override {}
 
  int numPar() const override
  {
    return 21;
  }
 
  double getPar(int ipar = 1) const override;
 
  bool getConstitutiveMatrix(Dune::FieldMatrix<double,6,6>& C,
                             bool invers = false) const override;
 
  bool getConstitutiveMatrix(Dune::FieldMatrix<double,3,3>& C,
                             bool invers = false) const override;
 
protected:
  std::ostream& write(std::ostream& os) const override;
 
private:
  double Cupper[21]; 
};
 
} // namespace Opm::Elasticity
 
#endif