materials.hh

Go to the documentation of this file.

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