1 # -*- coding: utf-8 -*-
4 from .geomsmesh import smesh
5 from salome.smesh import smeshBuilder
7 # -----------------------------------------------------------------------------
8 # --- maillage face de fissure pour identification zone de defaut
10 def genereMeshCalculZoneDefaut(facefiss, minSize, maxSize):
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]
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]"""
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")
40 coordsNoeudsFissure = []
41 nodeIds = meshFissure.GetNodesId()
43 coords = meshFissure.GetNodeXYZ(id)
44 coordsNoeudsFissure.append(coords[0])
45 coordsNoeudsFissure.append(coords[1])
46 coordsNoeudsFissure.append(coords[2])
47 return coordsNoeudsFissure