src/Tools/blocFissure/gmu/genereMeshCalculZoneDefaut.py

   1 # -*- coding: utf-8 -*-
   2
   3 import logging
   4 from geomsmesh import smesh
   5 from salome.smesh import smeshBuilder
   6
   7 # -----------------------------------------------------------------------------
   8 # --- maillage face de fissure pour identification zone de defaut
   9
  10 def genereMeshCalculZoneDefaut(facefiss, minSize, maxSize):
  11   """
  12   -Permet de générer un maillage sur l'objet géométrique 'facefiss' via
  13    l'algorithme NETGEN_1D2D :
  14       -SetMaxSize     = dimension max d'un élément (maxSize)
  15       -SetSecondOrder = élément quadratique (Y=1, N=0)
  16       -SetOptimize    = élément régulier (Y=1, N=0)
  17       -SetFineness    = finesse du maillage
  18        [very_coarse, coarse, moderate, fine, very_fine, custom]
  19        [0,           1,      2,        3,    4,         5     ]
  20       -SetMinSize     = dimension min d'un élément (minSize)
  21       -SetQuadAllowed = permission quadrangle dans maillage triangle
  22   -On récupère les coordonnées de chaque noeud de la fissure qu'on stocke
  23    dans une liste sous la forme : [X0, Y0, Z0, ..., Xn, Yn, Zn]"""
  24
  25   logging.info('start')
  26
  27   meshFissure = smesh.Mesh(facefiss)
  28   algo2d = meshFissure.Triangle(algo=smeshBuilder.NETGEN_1D2D)
  29   hypo2d = algo2d.Parameters()
  30   hypo2d.SetMaxSize( maxSize )
  31   hypo2d.SetSecondOrder( 0 )
  32   hypo2d.SetOptimize( 1 )
  33   hypo2d.SetFineness( 2 )
  34   hypo2d.SetMinSize( minSize )
  35   hypo2d.SetQuadAllowed( 0 )
  36   isDone = meshFissure.Compute()
  37   smesh.SetName(algo2d, "algo2d_zoneFiss")
  38   smesh.SetName(hypo2d, "hypo1d_zoneFiss")
  39
  40   coordsNoeudsFissure = []
  41   nodeIds = meshFissure.GetNodesId()
  42   for id in nodeIds:
  43     coords = meshFissure.GetNodeXYZ(id)
  44     coordsNoeudsFissure.append(coords[0])
  45     coordsNoeudsFissure.append(coords[1])
  46     coordsNoeudsFissure.append(coords[2])
  47   return coordsNoeudsFissure