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