1 // Copyright (C) 2007-2024 CEA, EDF, OPEN CASCADE
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.
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.
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
17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
20 // File : StdMeshers_ViscousLayers.hxx
21 // Created : Wed Dec 1 15:15:34 2010
22 // Author : Edward AGAPOV (eap)
24 #ifndef __StdMeshers_ViscousLayers_HXX__
25 #define __StdMeshers_ViscousLayers_HXX__
27 #include "SMESH_StdMeshers.hxx"
29 #include "SMESH_Hypothesis.hxx"
30 #include "SMESH_ProxyMesh.hxx"
31 #include "SMESH_ComputeError.hxx"
38 * \brief Hypothesis defining parameters of viscous layers
40 class STDMESHERS_EXPORT StdMeshers_ViscousLayers : public SMESH_Hypothesis
43 StdMeshers_ViscousLayers(int hypId, SMESH_Gen* gen);
45 // Set boundary shapes (faces in 3D, edges in 2D) either to exclude from
46 // treatment or to make the Viscous Layers on
47 void SetBndShapes(const std::vector<int>& shapeIds, bool toIgnore);
48 std::vector<int> GetBndShapes() const { return _shapeIds; }
49 bool IsToIgnoreShapes() const { return _isToIgnoreShapes; }
51 // Set total thickness of layers of prisms
52 void SetTotalThickness(double thickness);
53 double GetTotalThickness() const { return _thickness; }
55 // Set number of layers of prisms
56 void SetNumberLayers(int nb);
57 int GetNumberLayers() const { return _nbLayers; }
59 // Set factor (>1.0) of growth of layer thickness towards inside of mesh
60 void SetStretchFactor(double factor);
61 double GetStretchFactor() const { return _stretchFactor; }
63 // Method of computing node translation
64 enum ExtrusionMethod {
65 // node is translated along normal to a surface with possible further smoothing
67 // node is translated along the average normal of surrounding faces till
68 // intersection with a neighbor face translated along its own normal
69 // by the layers thickness
71 // node is translated along the average normal of surrounding faces
72 // by the layers thickness
75 void SetMethod( ExtrusionMethod how );
76 ExtrusionMethod GetMethod() const { return _method; }
78 // name of a group to create
79 void SetGroupName(const std::string& name);
80 const std::string& GetGroupName() const { return _groupName; }
81 static SMDS_MeshGroup* CreateGroup( const std::string& theName,
83 SMDSAbs_ElementType theType);
85 // Computes temporary 2D mesh to be used by 3D algorithm.
86 // Return SMESH_ProxyMesh for each SOLID in theShape
87 SMESH_ProxyMesh::Ptr Compute(SMESH_Mesh& theMesh,
88 const TopoDS_Shape& theShape,
89 const bool toMakeN2NMap=false) const;
91 // Checks compatibility of assigned StdMeshers_ViscousLayers hypotheses
92 static SMESH_ComputeErrorPtr
93 CheckHypothesis(SMESH_Mesh& aMesh,
94 const TopoDS_Shape& aShape,
95 SMESH_Hypothesis::Hypothesis_Status& aStatus);
97 // Compute thickness of the 1st layer
98 static double Get1stLayerThickness( double thickTotal,
102 // Checks if viscous layers should be constructed on a shape
103 bool IsShapeWithLayers(int shapeIndex) const;
105 virtual std::ostream & SaveTo(std::ostream & save);
106 virtual std::istream & LoadFrom(std::istream & load);
109 * \brief Initialize my parameter values by the mesh built on the geometry
110 * \param theMesh - the built mesh
111 * \param theShape - the geometry of interest
112 * \retval bool - true if parameter values have been successfully defined
114 virtual bool SetParametersByMesh(const SMESH_Mesh* theMesh, const TopoDS_Shape& theShape);
117 * \brief Initialize my parameter values by default parameters.
118 * \retval bool - true if parameter values have been successfully defined
120 virtual bool SetParametersByDefaults(const TDefaults& /*dflts*/, const SMESH_Mesh* /*theMesh*/=0)
123 static const char* GetHypType() { return "ViscousLayers"; }
127 std::vector<int> _shapeIds;
128 bool _isToIgnoreShapes;
131 double _stretchFactor;
132 ExtrusionMethod _method;
133 std::string _groupName;
139 // sets a sub-mesh event listener to clear sub-meshes of sub-shapes of
140 // the main shape when sub-mesh of the main shape is cleared,
141 // for example to clear sub-meshes of FACEs when sub-mesh of a SOLID
143 void ToClearSubWithMain( SMESH_subMesh* sub, const TopoDS_Shape& main);