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