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

22526: SMESH 2864 - Projection and Extrusion

Browse Source 
Make attempts to fix the projected mesh using smoothing in case of non-homothetic
  source and target faces
This commit is contained in:
eap 2014-11-14 20:23:20 +03:00
parent 380c397cfd
commit bff3771eab
1 changed files with 39 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

40
src/StdMeshers/StdMeshers_Prism_3D.cxx
View File

@ -49,6 +49,7 @@
#include <Bnd_B3d.hxx> #include <Bnd_B3d.hxx>
#include <Geom2dAdaptor_Curve.hxx> #include <Geom2dAdaptor_Curve.hxx>
#include <Geom2d_Line.hxx> #include <Geom2d_Line.hxx>
#include <GeomLib_IsPlanarSurface.hxx>
#include <Geom_Curve.hxx> #include <Geom_Curve.hxx>
#include <TopExp.hxx> #include <TopExp.hxx>
#include <TopExp_Explorer.hxx> #include <TopExp_Explorer.hxx>
@ -2214,11 +2215,11 @@ bool StdMeshers_Prism_3D::projectBottomToTop( const gp_Trsf & bottom
case 3: { case 3: {
newFace = myHelper->AddFace(nodes[0], nodes[1], nodes[2]); newFace = myHelper->AddFace(nodes[0], nodes[1], nodes[2]);
break; break;
} }
case 4: { case 4: {
newFace = myHelper->AddFace( nodes[0], nodes[1], nodes[2], nodes[3] ); newFace = myHelper->AddFace( nodes[0], nodes[1], nodes[2], nodes[3] );
break; break;
} }
default: default:
newFace = meshDS->AddPolygonalFace( nodes ); newFace = meshDS->AddPolygonalFace( nodes );
} }
@ -2228,6 +2229,41 @@ bool StdMeshers_Prism_3D::projectBottomToTop( const gp_Trsf & bottom
myHelper->SetElementsOnShape( oldSetElemsOnShape ); myHelper->SetElementsOnShape( oldSetElemsOnShape );
// Check the projected mesh
if ( thePrism.myNbEdgesInWires.size() > 1 && // there are holes
topHelper.IsDistorted2D( topSM, /*checkUV=*/false ))
{
SMESH_MeshEditor editor( topHelper.GetMesh() );
// smooth in 2D or 3D?
TopLoc_Location loc;
Handle(Geom_Surface) surface = BRep_Tool::Surface( topFace, loc );
bool isPlanar = GeomLib_IsPlanarSurface( surface ).IsPlanar();
bool isFixed = false;
set<const SMDS_MeshNode*> fixedNodes;
for ( int iAttemp = 0; !isFixed && iAttemp < 10; ++iAttemp )
{
TIDSortedElemSet faces;
for ( faceIt = topSMDS->GetElements(); faceIt->more(); )
faces.insert( faces.end(), faceIt->next() );
SMESH_MeshEditor::SmoothMethod algo =
iAttemp ? SMESH_MeshEditor::CENTROIDAL : SMESH_MeshEditor::LAPLACIAN;
// smoothing
editor.Smooth( faces, fixedNodes, algo, /*nbIterations=*/ 10,
/*theTgtAspectRatio=*/1.0, /*the2D=*/!isPlanar);
isFixed = !topHelper.IsDistorted2D( topSM, /*checkUV=*/true );
}
if ( !isFixed )
return toSM( error( TCom("Projection from face #") << botSM->GetId()
<< " to face #" << topSM->GetId()
<< " failed: inverted elements created"));
}
return true; return true;
} }