1 # -*- coding: utf-8 -*-
5 from geomsmesh import geompy
6 from geomsmesh import smesh
8 from putName import putName
10 def mailleAretesEtJonction(internalBoundary, aretesVivesCoupees, dmoyen):
12 edges de bord, faces défaut à respecter
16 aFilterManager = smesh.CreateFilterManager()
17 nbAdded, internalBoundary, _NoneGroup = internalBoundary.MakeBoundaryElements( SMESH.BND_1DFROM2D, '', '', 0, [ ])
19 unCritere = smesh.GetCriterion(SMESH.EDGE,SMESH.FT_FreeBorders,SMESH.FT_Undefined,0)
20 criteres.append(unCritere)
21 filtre = smesh.GetFilterFromCriteria(criteres)
22 bordsLibres = internalBoundary.MakeGroupByFilter( 'bords', filtre )
23 smesh.SetName(bordsLibres, 'bordsLibres')
25 # --- pour aider l'algo hexa-tetra à ne pas mettre de pyramides à l'exterieur des volumes repliés sur eux-mêmes
26 # on désigne les faces de peau en quadrangles par le groupe "skinFaces"
28 skinFaces = internalBoundary.CreateEmptyGroup( SMESH.FACE, 'skinFaces' )
29 nbAdd = skinFaces.AddFrom( internalBoundary.GetMesh() )
31 # --- maillage des éventuelles arêtes vives entre faces reconstruites
33 if len(aretesVivesCoupees) > 0:
34 aretesVivesC = geompy.MakeCompound(aretesVivesCoupees)
35 meshAretesVives = smesh.Mesh(aretesVivesC)
36 algo1d = meshAretesVives.Segment()
37 hypo1d = algo1d.LocalLength(dmoyen,[],1e-07)
38 putName(algo1d.GetSubMesh(), "aretesVives")
39 putName(algo1d, "algo1d_aretesVives")
40 putName(hypo1d, "hypo1d_aretesVives")
41 isDone = meshAretesVives.Compute()
42 logging.info("aretesVives fini")
43 grpAretesVives = meshAretesVives.CreateEmptyGroup( SMESH.EDGE, 'grpAretesVives' )
44 nbAdd = grpAretesVives.AddFrom( meshAretesVives.GetMesh() )
46 return (internalBoundary, bordsLibres, grpAretesVives)