src/Tools/blocFissure/gmu/mailleAretesEtJonction.py

   1 # -*- coding: utf-8 -*-
   2
   3 import logging
   4
   5 from .geomsmesh import geompy
   6 from .geomsmesh import smesh
   7 import SMESH
   8
   9 from .putName import putName
  10
  11 def mailleAretesEtJonction(internalBoundary, aretesVivesCoupees, lgAretesVives):
  12   """
  13   edges de bord, faces défaut à respecter
  14   """
  15   logging.info('start')
  16
  17   aFilterManager = smesh.CreateFilterManager()
  18   nbAdded, internalBoundary, _NoneGroup = internalBoundary.MakeBoundaryElements( SMESH.BND_1DFROM2D, '', '', 0, [  ])
  19   criteres = []
  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')
  25
  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"
  28
  29   skinFaces = internalBoundary.CreateEmptyGroup( SMESH.FACE, 'skinFaces' )
  30   nbAdd = skinFaces.AddFrom( internalBoundary.GetMesh() )
  31
  32   # --- maillage des éventuelles arêtes vives entre faces reconstruites
  33
  34   grpAretesVives = None
  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() )
  47
  48   return (internalBoundary, bordsLibres, grpAretesVives)