isotropicelasticeqvar.h

Go to the documentation of this file.

#ifndef  _MY_IsotropicElasticEqVar_
#define _MY_IsotropicElasticEqVar_
#include <fe/fe_system.h>
#include <fe/fe_q.h>
#include "dealfewrapper.h"


class IsotropicElasticEqVar : public DealFEWrapper
{
 private:
  FE_Q<3> m_FE;
  FESystem<3> m_FES;
  QGauss<3> m_cQG;
  QGauss<2> m_fQG;
  VecIndex m_global_dofs_map;
  static const Index m_local_dof_face_cell_map[6][4];
  OrthoMesh &m_mesh;
  enum COMPONENTS{U0,U1,U2};
  double m_l1,m_l2; //Lame constants
 protected:

 public:
  unsigned n_max_dof_coupling(){return 81;} 
  Point3D get_global_dof_position(Index cDof){return getCell()->vertex(static_cast<VertexDirection3D>(getFE()->system_to_component_index(cDof).second));}

  const VecIndex& local_to_global_dofs_map(){return m_global_dofs_map;} 
  //Index local_face_to_local_cell_map(Index face,Index local_dof_face,unsigned component);  /*!<Get the local dof number relative to the cell given the local dof relative to the face*/
  void setCell(OrthoMesh::Cell_It &cell); 
  int nullProduct(unsigned i,unsigned j);


  //Functions to build the local system
  double localEqVar(unsigned i,unsigned j,unsigned qPoint);
  double localRHS(unsigned i,unsigned j);

  IsotropicElasticEqVar(OrthoMesh &mesh,  VecDouble &cK, double l1,double l2, Function3D &fD, Function3D &fN,unsigned debugLevel=0);

};

#endif

---