]> SALOME platform Git repositories - modules/smesh.git/commitdiff
Salome HOME
Prise en compte de la courbure de la face à mailler
authorGERALD NICOLAS <D68518@dsp0864451.atlas.edf.fr>
Fri, 16 Apr 2021 08:51:45 +0000 (10:51 +0200)
committerGERALD NICOLAS <D68518@dsp0864451.atlas.edf.fr>
Fri, 16 Apr 2021 08:51:45 +0000 (10:51 +0200)
src/Tools/blocFissure/CasTests/execution_Cas.py
src/Tools/blocFissure/gmu/genereMeshCalculZoneDefaut.py
src/Tools/blocFissure/gmu/insereFissureLongue_d.py
src/Tools/blocFissure/gmu/insereFissureLongue_e.py
src/Tools/blocFissure/gmu/mailleFacesFissure.py
src/Tools/blocFissure/gmu/mailleFacesPeau.py

index 31313568175f3f69a6c99d50d00eec701fc35d8d..983899999920b57277495b690fb298c5d003b7fb 100644 (file)
@@ -30,7 +30,7 @@ from blocFissure.gmu.casStandard import casStandard
 # -----------------------------------------------------------------------------------------------
 LOGFILE = os.path.join(tempfile.gettempdir(),"blocFissure.log")
 #initLog.setDebug(LOGFILE)     # debug = 10
-initLog.setVerbose(LOGFILE)    # info = 20
+#initLog.setVerbose(LOGFILE)   # info = 20
 #initLog.setRelease(LOGFILE)   # warning = 30
 #initLog.setPerfTests(LOGFILE) # critical = 50
 #initLog.setAlways(LOGFILE)    # critical = 50
index f245e9c8e216bf8dfaef9d54ac9c7c5df1cbcd2d..3ab7cf8149b9c362c4f9e12926b19f46e3dccf9e 100644 (file)
@@ -75,6 +75,9 @@ def genereMeshCalculZoneDefaut(facefiss, minSize, maxSize, \
     hypo2d.SetOptimize( 1 )
     hypo2d.SetFineness( 2 )
     hypo2d.SetMinSize( minSize )
+    hypo2d.SetChordalErrorEnabled (True)
+    hypo2d.SetChordalError( maxSize*0.25 )
+    hypo2d.SetUseSurfaceCurvature (True)
     hypo2d.SetQuadAllowed( 0 )
   putName(algo2d, "{}_2d_zoneFiss".format(mailleur), i_pref=nro_cas)
   putName(hypo2d, "hypo1d_zoneFiss", i_pref=nro_cas)
index 206b71c6b6ac8f85eef79c6d4f4d969ced91536b..cff54376abe18f9c5ac24156270415130c155545 100644 (file)
@@ -55,6 +55,9 @@ def insereFissureLongue_d (facePeau, edgePeauFiss, groupEdgesBordPeau, bordsLibr
     hypo2d.SetOptimize( 1 )
     hypo2d.SetFineness( 2 )
     hypo2d.SetMinSize( 2 )
+    hypo2d.SetChordalErrorEnabled (True)
+    hypo2d.SetChordalError( 250. )
+    hypo2d.SetUseSurfaceCurvature (True)
     hypo2d.SetQuadAllowed( 0 )
   putName(algo2d.GetSubMesh(), "facePeau", i_pref=nro_cas)
   putName(algo2d, "{}_2d_facePeau".format(mailleur), i_pref=nro_cas)
index 3a505845ebcf86252eaa9ce52e35d380d2fcf2b2..6cfd268ffad0ad3716358bc490ec85f9694d3bd3 100644 (file)
@@ -55,6 +55,9 @@ def insereFissureLongue_e (faceFiss, edgePeauFiss, groupEdgesPeauFiss, group_gen
     hypo2d.SetOptimize( 1 )
     hypo2d.SetFineness( 2 )
     hypo2d.SetMinSize( 2 )
+    hypo2d.SetChordalErrorEnabled (True)
+    hypo2d.SetChordalError( mesh_size*0.25 )
+    hypo2d.SetUseSurfaceCurvature (True)
     hypo2d.SetQuadAllowed( 0 )
   putName(algo2d.GetSubMesh(), "faceFiss", i_pref=nro_cas)
   putName(algo2d, "{}_2d_faceFiss".format(mailleur), i_pref=nro_cas)
index 0b63cfa9d27a62639f0cdd95acd058f9a2078cd3..ac01b8cd1c4c36032a9d7cc6726374e8b4251d7c 100644 (file)
@@ -37,7 +37,7 @@ def mailleFacesFissure(faceFissureExterne, \
   logging.info("Maillage avec %s pour le cas n°%s", mailleur, nro_cas)
 
   meshFaceFiss = smesh.Mesh(faceFissureExterne)
-  logging.info("Maillage avec %s", mailleur)
+  putName(meshFaceFiss.GetMesh(), "faceFiss", i_pref=nro_cas)
   if ( mailleur == "MeshGems"):
     algo2d = meshFaceFiss.Triangle(algo=smeshBuilder.MG_CADSurf)
     hypo2d = algo2d.Parameters()
@@ -54,22 +54,23 @@ def mailleFacesFissure(faceFissureExterne, \
     hypo2d.SetOptimize( 1 )
     hypo2d.SetFineness( 2 )
     hypo2d.SetMinSize( rayonPipe/float(nbsegRad) )
+    hypo2d.SetChordalErrorEnabled (True)
     hypo2d.SetChordalError( areteFaceFissure*0.25 )
+    hypo2d.SetUseSurfaceCurvature (True)
     hypo2d.SetQuadAllowed( 0 )
-  putName(algo2d.GetSubMesh(), "faceFiss", i_pref=nro_cas)
   putName(algo2d, "{}_2d_faceFiss".format(mailleur), i_pref=nro_cas)
   putName(hypo2d, "hypo2d_faceFiss", i_pref=nro_cas)
 
   logging.info("UseExisting1DElements depuis '%s'", edgesPipeFissureExterneC.GetName())
   algo1d = meshFaceFiss.UseExisting1DElements(geom=edgesPipeFissureExterneC)
-  hypo1d = algo1d.SourceEdges([ edgeFaceFissGroup ],0,0)
   putName(algo1d.GetSubMesh(), "edgeFissPeau", i_pref=nro_cas)
-  putName(algo1d, "algo1d_edgeFissPeau", i_pref=nro_cas)
-  putName(hypo1d, "hypo1d_edgeFissPeau", i_pref=nro_cas)
+  putName(algo1d, "edgeFissPeau_UseExisting1DElements", i_pref=nro_cas)
+  hypo1d = algo1d.SourceEdges([ edgeFaceFissGroup ],0,0)
+  putName(hypo1d, "edgeFissPeau_SourceEdges", i_pref=nro_cas)
 
-  grpFaceFissureExterne = meshFaceFiss.GroupOnGeom(faceFissureExterne, "fisOutPi", SMESH.FACE)
-  grpEdgesPeauFissureExterne = meshFaceFiss.GroupOnGeom(edgesPeauFissureExterneC,'edgesPeauFissureExterne',SMESH.EDGE)
   grpEdgesPipeFissureExterne = meshFaceFiss.GroupOnGeom(edgesPipeFissureExterneC,'edgesPipeFissureExterne',SMESH.EDGE)
+  grpEdgesPeauFissureExterne = meshFaceFiss.GroupOnGeom(edgesPeauFissureExterneC,'edgesPeauFissureExterne',SMESH.EDGE)
+  grpFaceFissureExterne = meshFaceFiss.GroupOnGeom(faceFissureExterne, "fisOutPi", SMESH.FACE)
 
   is_done = meshFaceFiss.Compute()
   text = "meshFaceFiss.Compute"
index e8289c22f9af556af4c827cea668230e30ab7341..d216abbdf79e39b9d4f83025c7e5b2ac5326e239 100644 (file)
@@ -79,7 +79,9 @@ def mailleFacesPeau(partitionsPeauFissFond, idFillingFromBout, facesDefaut, \
       hypo2d.SetOptimize( 1 )
       hypo2d.SetFineness( 2 )
       hypo2d.SetMinSize( rayonPipe/float(nbsegRad) )
+      hypo2d.SetChordalErrorEnabled (True)
       hypo2d.SetChordalError( dmoyen*0.25 )
+      hypo2d.SetUseSurfaceCurvature (True)
       hypo2d.SetQuadAllowed( 0 )
     putName(algo2d, "{}_2d_facePeau".format(mailleur), ifil, nro_cas)
     putName(hypo2d, "hypo2d_facePeau", ifil, nro_cas)