mechanics_stress.hh

Go to the documentation of this file.

 
#ifndef MECHANICS_STRESS_HH
#define MECHANICS_STRESS_HH
 
#include "space/function.hh"
 
namespace mechanics {
 
  using concepts::Real;
  using concepts::Real2d;
  using concepts::Real3d;
 
  // ************************************************************** Stress3D **
 
  template<class F, class G = Real>
  class Stress3D : public concepts::ElementFunction<F,G>
  {
  public:
    Stress3D(const F lambda, const F mu);
 
    virtual ~Stress3D(){}
 
    virtual Stress3D<F,G>* clone() const
    {
      return new Stress3D(lambda_, mu_);
    }
 
    virtual uint n() const {
      return 9;
    }
 
    virtual void operator()(const concepts::Element<G>& elm,
                            const concepts::Array<F>& coeff, concepts::Array<F>& val,
                            const uint *i) const;
 
    virtual void operator()(const concepts::Element<G>& elm, const uint* j,
                            concepts::Array<F>& val, const uint* i) const;
 
    virtual void operator()(const concepts::Element<G>& elm,
                            const concepts::Array<F>& coeff, concepts::Array<F>& val,
                            const Real p, const Real t = 0.0) const;
 
    virtual void operator()(const concepts::Element<G>& elm,
                            const concepts::Array<F>& coeff, concepts::Array<F>& val,
                            const Real2d& p, const Real t = 0.0) const;
 
    virtual void operator()(const concepts::Element<G>& elm,
                            const concepts::Array<F>& coeff, concepts::Array<F>& val,
                            const Real3d& p, const Real t = 0.0) const;
 
  protected:
    std::ostream& info(std::ostream &os) const
    {
      return os << concepts::typeOf(*this);
    }
 
  private:
    const F lambda_;
    const F mu_;
 
    void process_(concepts::Array<F> &stressVal,
                  const concepts::Array<F> &strainVal) const;
 
  };
 
} // namespace mechanics
 
#endif // MECHANICS_STRESS_HH