app-cwinkelm/equivalentStress.hh

Go to the documentation of this file.

#ifndef equivalentStress_hh
#define equivalentStress_hh

#include "basics/typedefs.hh"
#include "function/vector.hh"
#include "elasticityTensor.hh"
#include "vectorial.hh"
#include "linDG2D/space.hh"
#include "operator/compositions.hh"

using concepts::Real;

// ********************************************** EquivalentStressHypothesis **

template<int dim>
class EquivalentStressHypothesis {
 public:
  virtual Real operator() (concepts::Mapping<Real, dim>& sigma) const = 0;
  virtual concepts::Point<Real, dim>
  principalStresses(concepts::Mapping<Real, dim>& sigma) const;
};

// **************************************** MaximumPrincipalStressHypothesis **

template<int dim>
class MaximumPrincipalStressHypothesis : 
  public EquivalentStressHypothesis<dim> {
 public:
  virtual Real operator() (concepts::Mapping<Real, dim>& sigma) const;
};

// ******************************************** MaximumShearStressHypothesis **

template<int dim>
class MaximumShearStressHypothesis : 
  public EquivalentStressHypothesis<dim> {
 public:
  virtual Real operator() (concepts::Mapping<Real, dim>& sigma) const;
};

// ************************************** MaximumShearStrainEnergyHypothesis **

template<int dim>
class MaximumShearStrainEnergyHypothesis : 
  public EquivalentStressHypothesis<dim> {
 public:
  virtual Real operator() (concepts::Mapping<Real, dim>& sigma) const;
};

// ****************************************************** EquivalentStress2d **

class EquivalentStress2d : public concepts::Operator<Real> {
 public:
  EquivalentStress2d(const EquivalentStressHypothesis<2>& esh,
                     const ElasticityTensor<2>& et,
                     const vectorial::Space<Real>& spc)
    : concepts::Operator<Real>(spc.get(0)->dim(), spc.dim())
    , esh_(esh), et_(et), spc_(spc) {
      conceptsAssert(spc.vdim() == 2, concepts::Assertion());
    }
  void operator()
    (const concepts::Function<Real>& fncY, concepts::Function<Real>& fncX);
  const concepts::Space<Real>& spaceX() const { return *spc_.get(0); }
  const concepts::Space<Real>& spaceY() const { return spc_; }
protected:
  virtual std::ostream& info (std::ostream &os) const;
private:
  const EquivalentStressHypothesis<2>& esh_;
  const ElasticityTensor<2>& et_;
  const vectorial::Space<Real>& spc_;
};

#endif // equivalentStress_hh