-# -*- coding: latin-1 -*-
+# -*- coding: utf8 -*-
# Copyright (C) 2009-2016 CEA/DEN, EDF R&D
#
# This library is free software; you can redistribute it and/or
Face_bas = geompy.GetFaceNearPoint(Bielle_geom, Pt_B)
-# Param�tre pour:
-# -r�cup�ration des points.
+# Paramètre pour:
+# -récupération des points.
# -et associations.
u_1 = 0.75#0.8
u_2 = 0.5
Propa.setLaw( Law )# appliquer la loi de discretisation sur tout le modele et generer le maillage
#====================================
-# G�n�rer des maillages
+# Générer des maillages
#====================================
-print " --- MAILLAGE HEXAHEDRIQUE --- "
+print(" --- MAILLAGE HEXAHEDRIQUE --- ")
mesh_hexas = hexablock.mesh(doc, "Bielle:hexas")
-print "Nombre d hexaedres:" , mesh_hexas.NbHexas()
-print "Nombre de quadrangles:", mesh_hexas.NbQuadrangles()
-print "Nombre de segments:" , mesh_hexas.NbEdges()
-print "Nombre de noeuds:" , mesh_hexas.NbNodes()
+print("Nombre d hexaedres:" , mesh_hexas.NbHexas())
+print("Nombre de quadrangles:", mesh_hexas.NbQuadrangles())
+print("Nombre de segments:" , mesh_hexas.NbEdges())
+print("Nombre de noeuds:" , mesh_hexas.NbNodes())
-print " --- MAILLAGE QUADRANGULAIRE --- "
+print(" --- MAILLAGE QUADRANGULAIRE --- ")
mesh_quads = hexablock.mesh(doc, "Bielle:quads", 2)
-print "Nombre d hexaedres:" , mesh_quads.NbHexas()
-print "Nombre de quadrangles:", mesh_quads.NbQuadrangles()
-print "Nombre de segments:" , mesh_quads.NbEdges()
-print "Nombre de noeuds:" , mesh_quads.NbNodes()
+print("Nombre d hexaedres:" , mesh_quads.NbHexas())
+print("Nombre de quadrangles:", mesh_quads.NbQuadrangles())
+print("Nombre de segments:" , mesh_quads.NbEdges())
+print("Nombre de noeuds:" , mesh_quads.NbNodes())
-print " --- MAILLAGE 1D --- "
+print(" --- MAILLAGE 1D --- ")
mesh_edges = hexablock.mesh(doc, "Bielle:edges", 1)
-print "Nombre d hexaedres:" , mesh_edges.NbHexas()
-print "Nombre de quadrangles:", mesh_edges.NbQuadrangles()
-print "Nombre de segments:" , mesh_edges.NbEdges()
-print "Nombre de noeuds:" , mesh_edges.NbNodes()
+print("Nombre d hexaedres:" , mesh_edges.NbHexas())
+print("Nombre de quadrangles:", mesh_edges.NbQuadrangles())
+print("Nombre de segments:" , mesh_edges.NbEdges())
+print("Nombre de noeuds:" , mesh_edges.NbNodes())
-print " --- MAILLAGE des sommets --- "
+print(" --- MAILLAGE des sommets --- ")
mesh_nodes = hexablock.mesh(doc, "Bielle:nodes", 0)
-print "Nombre d hexaedres:" , mesh_nodes.NbHexas()
-print "Nombre de quadrangles:", mesh_nodes.NbQuadrangles()
-print "Nombre de segments:" , mesh_nodes.NbEdges()
-print "Nombre de noeuds:" , mesh_nodes.NbNodes()
-xaedres:" , mesh_edges.NbHexas()
-print "Nombre de quadrangles:", mesh_edges.NbQuadrangles()
-print "Nombre de segments:"
+print("Nombre d hexaedres:" , mesh_nodes.NbHexas())
+print("Nombre de quadrangles:", mesh_nodes.NbQuadrangles())
+print("Nombre de segments:" , mesh_nodes.NbEdges())
+print("Nombre de noeuds:" , mesh_nodes.NbNodes())
# dans cette première version, on considère une grille cylindrique avec nr = 1 et nl = 1
if (nr !=1) or (nl != 1):
- print "nr ou nl invalide(s) !!"
+ print("nr ou nl invalide(s) !!")
return None
vertices = []
li_z = [nl, 0] # indice nl => face du haut, indice 0 => face du bas
li_x = [nr, 0] # indice nr => externe, indice 0 => interne
- li_y = range(na) # indices de 0 à (na-1)
+ li_y = list(range(na)) # indices de 0 à (na-1)
for k in li_z: # suivant z
for i in li_x: # suivant x
for j in li_y: # suivant y
geom_vert_grd.extend(geom_vert_grd_bas)
-print "nombre vertex grand cylindre geom = ", len(geom_vert_grd)
+print("nombre vertex grand cylindre geom = ", len(geom_vert_grd))
# ajout des vertex de la géométrie à l'étude :
for num, vertex in enumerate(geom_vert_grd):
geom_vert_pte_bas = [geompy.MakeVertexWithRef(vertex, 0.0, 0.0, -h) for vertex in geom_vert_pte]
geom_vert_pte.extend(geom_vert_pte_bas)
-print "nombre vertex petit cylindre geom = ", len(geom_vert_pte)
+print("nombre vertex petit cylindre geom = ", len(geom_vert_pte))
# ajout des vertex de la géométrie à l'étude :
for num, vertex in enumerate(geom_vert_pte):
# association :
set_association_from_cylindrical(mod_pte, geom_vert_pte)
-
-le meme ordre :
-# sommets externes puis internes) :
-geom_vert_pte_bas = [geompy.MakeVertexWithRef(vertex, 0.0, 0.0, -h) for vertex in geom_vert_pte]
-geom_vert_pte.extend(geom_vert_pte_bas)
-
-print "nombre vertex petit cylindre geom = ", len(geom_vert_pte)
-
-# ajout des vertex
\ No newline at end of file