ForkMaster/smesh
1
0
Fork 0
mirror of https://git.salome-platform.org/gitpub/modules/smesh.git synced 2024-12-25 17:00:34 +05:00
Code Issues Packages Projects Releases Wiki Activity

IPAL54633: Bad viscous layers

Browse Source
This commit is contained in:
eap 2020-03-20 18:57:40 +03:00
parent 1f0895c06c
commit d487f882d9
1 changed files with 16 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
src/StdMeshers/StdMeshers_ViscousLayers.cxx
View File

@ -3304,12 +3304,22 @@ bool _ViscousBuilder::findShapesToSmooth( _SolidData& data )
if ( isC1 )
{
double maxEdgeLen = 3 * Min( eov._edges[0]->_maxLen, eov._hyp.GetTotalThickness() );
double eLen1 = SMESH_Algo::EdgeLength( TopoDS::Edge( dirOfEdges[i].first->_shape ));
double eLen2 = SMESH_Algo::EdgeLength( TopoDS::Edge( dirOfEdges[j].first->_shape ));
if ( eLen1 < maxEdgeLen ) eov._eosC1.push_back( dirOfEdges[i].first );
if ( eLen2 < maxEdgeLen ) eov._eosC1.push_back( dirOfEdges[j].first );
dirOfEdges[i].first = 0;
dirOfEdges[j].first = 0;
for ( int isJ = 0; isJ < 2; ++isJ ) // loop on [i,j]
{
size_t k = isJ ? j : i;
const TopoDS_Edge& e = TopoDS::Edge( dirOfEdges[k].first->_shape );
double eLen = SMESH_Algo::EdgeLength( e );
if ( eLen < maxEdgeLen )
{
TopoDS_Shape oppV = SMESH_MesherHelper::IthVertex( 0, e );
if ( oppV.IsSame( V ))
oppV = SMESH_MesherHelper::IthVertex( 1, e );
_EdgesOnShape* eovOpp = data.GetShapeEdges( oppV );
if ( dirOfEdges[k].second * eovOpp->_edges[0]->_normal < 0 )
eov._eosC1.push_back( dirOfEdges[k].first );
}
dirOfEdges[k].first = 0;
}
}
}
} // fill _eosC1 of VERTEXes