src/StdMeshers/StdMeshers_ViscousLayers.hxx

   1 // Copyright (C) 2007-2016  CEA/DEN, EDF R&D, OPEN CASCADE
   2 //
   3 // This library is free software; you can redistribute it and/or
   4 // modify it under the terms of the GNU Lesser General Public
   5 // License as published by the Free Software Foundation; either
   6 // version 2.1 of the License, or (at your option) any later version.
   7 //
   8 // This library is distributed in the hope that it will be useful,
   9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 // Lesser General Public License for more details.
  12 //
  13 // You should have received a copy of the GNU Lesser General Public
  14 // License along with this library; if not, write to the Free Software
  15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 //
  17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 //
  19
  20 // File      : StdMeshers_ViscousLayers.hxx
  21 // Created   : Wed Dec  1 15:15:34 2010
  22 // Author    : Edward AGAPOV (eap)
  23
  24 #ifndef __StdMeshers_ViscousLayers_HXX__
  25 #define __StdMeshers_ViscousLayers_HXX__
  26
  27 #include "SMESH_StdMeshers.hxx"
  28
  29 #include "SMESH_Hypothesis.hxx"
  30 #include "SMESH_ProxyMesh.hxx"
  31 #include "SMESH_ComputeError.hxx"
  32
  33 #include <vector>
  34
  35 /*!
  36  * \brief Hypothesis defining parameters of viscous layers
  37  */
  38 class STDMESHERS_EXPORT StdMeshers_ViscousLayers : public SMESH_Hypothesis
  39 {
  40 public:
  41   StdMeshers_ViscousLayers(int hypId, int studyId, SMESH_Gen* gen);
  42
  43   // Set boundary shapes (faces in 3D, edges in 2D) either to exclude from
  44   // treatment or to make the Viscous Layers on
  45   void   SetBndShapes(const std::vector<int>& shapeIds, bool toIgnore);
  46   std::vector<int> GetBndShapes() const { return _shapeIds; }
  47   bool   IsToIgnoreShapes() const { return _isToIgnoreShapes; }
  48
  49   // Set total thickness of layers of prisms
  50   void   SetTotalThickness(double thickness);
  51   double GetTotalThickness() const { return _thickness; }
  52
  53   // Set number of layers of prisms
  54   void   SetNumberLayers(int nb);
  55   int    GetNumberLayers() const { return _nbLayers; }
  56
  57   // Set factor (>1.0) of growth of layer thickness towards inside of mesh
  58   void   SetStretchFactor(double factor);
  59   double GetStretchFactor() const { return _stretchFactor; }
  60
  61   // Method of computing node translation
  62   enum ExtrusionMethod {
  63     // node is translated along normal to a surface with possible further smoothing
  64     SURF_OFFSET_SMOOTH,
  65     // node is translated along the average normal of surrounding faces till
  66     // intersection with a neighbor face translated along its own normal
  67     // by the layers thickness
  68     FACE_OFFSET,
  69     // node is translated along the average normal of surrounding faces
  70     // by the layers thickness
  71     NODE_OFFSET
  72   };
  73   void   SetMethod( ExtrusionMethod how );
  74   ExtrusionMethod GetMethod() const { return _method; }
  75
  76   // Computes temporary 2D mesh to be used by 3D algorithm.
  77   // Return SMESH_ProxyMesh for each SOLID in theShape
  78   SMESH_ProxyMesh::Ptr Compute(SMESH_Mesh&         theMesh,
  79                                const TopoDS_Shape& theShape,
  80                                const bool          toMakeN2NMap=false) const;
  81
  82   // Checks compatibility of assigned StdMeshers_ViscousLayers hypotheses
  83   static SMESH_ComputeErrorPtr
  84     CheckHypothesis(SMESH_Mesh&                          aMesh,
  85                     const TopoDS_Shape&                  aShape,
  86                     SMESH_Hypothesis::Hypothesis_Status& aStatus);
  87
  88   // Checks if viscous layers should be constructed on a shape
  89   bool IsShapeWithLayers(int shapeIndex) const;
  90
  91   virtual std::ostream & SaveTo(std::ostream & save);
  92   virtual std::istream & LoadFrom(std::istream & load);
  93
  94   /*!
  95    * \brief Initialize my parameter values by the mesh built on the geometry
  96     * \param theMesh - the built mesh
  97     * \param theShape - the geometry of interest
  98     * \retval bool - true if parameter values have been successfully defined
  99    */
 100   virtual bool SetParametersByMesh(const SMESH_Mesh* theMesh, const TopoDS_Shape& theShape);
 101
 102   /*!
 103    * \brief Initialize my parameter values by default parameters.
 104    *  \retval bool - true if parameter values have been successfully defined
 105    */
 106   virtual bool SetParametersByDefaults(const TDefaults& dflts, const SMESH_Mesh* theMesh=0)
 107   { return false; }
 108
 109   static const char* GetHypType() { return "ViscousLayers"; }
 110
 111  private:
 112
 113   std::vector<int> _shapeIds;
 114   bool             _isToIgnoreShapes;
 115   int              _nbLayers;
 116   double           _thickness;
 117   double           _stretchFactor;
 118   ExtrusionMethod  _method;
 119 };
 120
 121 class SMESH_subMesh;
 122 namespace VISCOUS_3D
 123 {
 124   // sets a sub-mesh event listener to clear sub-meshes of sub-shapes of
 125   // the main shape when sub-mesh of the main shape is cleared,
 126   // for example to clear sub-meshes of FACEs when sub-mesh of a SOLID
 127   // is cleared
 128   void ToClearSubWithMain( SMESH_subMesh* sub, const TopoDS_Shape& main);
 129 }
 130
 131 #endif