bf3D_linearElasticity.hh

Go to the documentation of this file.

 
#ifndef HP3D_LINEARELASTICITY_HH
#define HP3D_LINEARELASTICITY_HH
 
#include "basics/typedefs.hh"
#include "operator/bilinearForm.hh"
 
namespace concepts {
  // forward declarations
  template <class F>
  class Element;
 
  template <class F>
  class ElementMatrix;
}
 
namespace vectorial {
  // forward declaration
  template<class F, class G>
  class BilinearForm;
}
 
namespace hp3D {
  // forward declaration
  class Hexahedron;
  
  using concepts::Real;
  
  template <class F>
  class LinearElasticity : public concepts::BilinearForm<F>
  {
  public:
    LinearElasticity(const uint &i, const uint &j,
                     const Real &lambda, const Real &mu);
    
    virtual ~LinearElasticity();
    
    virtual void operator()(const concepts::Element<Real> &elmX,
                            const concepts::Element<Real> &elmY,
                            concepts::ElementMatrix<F> &em) const;
    
    virtual LinearElasticity* clone() const
    {
      return new LinearElasticity(i_, j_, lambda_, mu_);
    }
    
    virtual std::ostream& info(std::ostream &os) const;
    
    static void setup(vectorial::BilinearForm<F,Real> &bf,
                      const Real &lambda, const Real &mu);
    
  private:
    
    void operator()(const Hexahedron &elmX, const Hexahedron &elmY,
                    concepts::ElementMatrix<F> &em) const;
    
    
    uint i_;
    uint j_;
    Real lambda_;
    Real mu_;
    
    concepts::BilinearForm<F> *bfA_;
    concepts::BilinearForm<F> *bfB_;
    Real coefA_;
    Real coefB_;
    
  };
    
} // namespace hp3D
 
#endif // HP3D_LINEARELASTICITY_HH