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