TractionField.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _TRACTION_FIELD_H
#define _TRACTION_FIELD_H
 
#include "Function.h"
#include "Tensor.h"
 
class TensorFunc;
class STensorFunc;
 
 
class TractionField : public TractionFunc
{
  const STensorFunc* sigma; 
  const TensorFunc* sigmaN; 
 
public:
  explicit TractionField(const STensorFunc& field);
  explicit TractionField(const TensorFunc& field);
  virtual ~TractionField() {}
 
  virtual bool isZero() const;
 
  virtual Vec3 deriv(const Vec3& x, const Vec3& n) const;
 
protected:
  virtual Vec3 evaluate(const Vec3& x, const Vec3& n) const;
};
 
 
class PressureField : public TractionFunc
{
  const RealFunc* pressure; 
  char            pdir;     
  const VecFunc*  pdfn;     
 
public:
  explicit PressureField(Real p, int dir = 0);
  explicit PressureField(const RealFunc* p, int dir = 0)
    : pressure(p), pdir(dir), pdfn(nullptr) {}
  PressureField(const RealFunc* p, const VecFunc* dir)
    : pressure(p), pdir(-1), pdfn(dir) {}
  virtual ~PressureField() { delete pressure; delete pdfn; }
 
  virtual bool isNormalPressure() const { return !pdfn && (pdir<1 || pdir>3); }
  virtual bool isZero() const { return pressure ? pressure->isZero() : true; }
 
  virtual Vec3 deriv(const Vec3& x, const Vec3& n) const;
 
protected:
  virtual Vec3 evaluate(const Vec3& x, const Vec3& n) const;
};
 
 
class ForceDirField : public TractionFunc
{
  typedef utl::Function<Real,Vec3> Vec3Func; 
 
  Vec3   X0;  
  Tensor Tlg; 
 
  const ScalarFunc* force; 
  const ScalarFunc* angle; 
  const Vec3Func*   fdir;  
  const RealFunc*   shape; 
 
  bool dirVec;  
  char rotAxis; 
 
public:
  ForceDirField(const ScalarFunc* f, const Vec3Func* d, const RealFunc* s,
                const Vec3& Xaxis = Vec3(Real(1),Real(0),Real(0)),
                const Vec3& Zaxis = Vec3(Real(0),Real(0),Real(1)),
                const Vec3& Xorig = Vec3(), bool v = false)
    : X0(Xorig), Tlg(Xaxis,Zaxis,false,true), force(f), angle(nullptr), fdir(d),
      shape(s), dirVec(v), rotAxis(0) {}
 
  ForceDirField(const ScalarFunc* f, const ScalarFunc* a, char x,
                const RealFunc* s,
                const Vec3& Xaxis = Vec3(Real(1),Real(0),Real(0)),
                const Vec3& Zaxis = Vec3(Real(0),Real(0),Real(1)),
                const Vec3& Xorig = Vec3())
    : X0(Xorig), Tlg(Xaxis,Zaxis,false,true), force(f), angle(a), fdir(nullptr),
      shape(s), dirVec(false), rotAxis(x) {}
 
  virtual ~ForceDirField();
 
  virtual bool isZero() const { return force ? force->isZero() : true; }
 
protected:
  virtual Vec3 evaluate(const Vec3& x, const Vec3&) const;
};
 
#endif