# -*- coding: iso-8859-1 -*- # Copyright (C) 2007-2024 CEA, EDF, OPEN CASCADE # # Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN, # CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS # # This library is free software; you can redistribute it and/or # modify it under the terms of the GNU Lesser General Public # License as published by the Free Software Foundation; either # version 2.1 of the License, or (at your option) any later version. # # This library is distributed in the hope that it will be useful, # but WITHOUT ANY WARRANTY; without even the implied warranty of # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU # Lesser General Public License for more details. # # You should have received a copy of the GNU Lesser General Public # License along with this library; if not, write to the Free Software # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA # # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com # # ======================================= # import salome salome.salome_init() import GEOM from salome.geom import geomBuilder geompy = geomBuilder.New() import SMESH, SALOMEDS from salome.smesh import smeshBuilder smesh = smeshBuilder.New() # Geometrie # ========= # Construire avec des hexahedres une piece faite de 2 cylindres de diametre different # et dont les axes se coupent orthogonalement, puis mailler. # Valeurs # ------- cx = 0 cy = 0 cz = 0 g_rayon = 100.0 g_hauteur = 500 p_rayon = 50.0 p_hauteur = 500 g_trim = 1000 # Gros cylindre # ------------- cpd = [] g_base = geompy.MakeVertex(cx, cy, cz) g_dir = geompy.MakeVectorDXDYDZ(0, 0, 1) t_hauteur = p_rayon+10.0 g_cyl = geompy.MakeCylinder(g_base, g_dir, g_rayon, g_hauteur) g_coupe = geompy.MakeVectorDXDYDZ(1, 0, 0) g_tools = [] g_tools.append(geompy.MakePlane(geompy.MakeVertex(cx+t_hauteur, cy, cz), g_coupe, g_trim)) g_tools.append(geompy.MakePlane(geompy.MakeVertex(cx-t_hauteur, cy, cz), g_coupe, g_trim)) g_partie = geompy.MakePartition([g_cyl], g_tools, [], [], geompy.ShapeType["SOLID"]) g_bas, g_centre, g_haut = geompy.SubShapeAllSorted(g_partie, geompy.ShapeType["SOLID"]) # Partie basse du gros cylindre # ----------------------------- b_hauteur = 10 b_base = 20 b_boite = geompy.MakeBox(cx-t_hauteur, cy-b_base, cz, cx-t_hauteur-b_hauteur, cy+b_base, cz+g_hauteur) cpd.append(b_boite) b_cyl = geompy.MakeCut(g_bas, b_boite) b_tools = [] b_tools.append(geompy.MakePlane(geompy.MakeVertex(cx-t_hauteur-b_hauteur, cy+b_base, cz), geompy.MakeVectorDXDYDZ( 1, 1, 0), g_trim)) b_tools.append(geompy.MakePlane(geompy.MakeVertex(cx-t_hauteur-b_hauteur, cy-b_base, cz), geompy.MakeVectorDXDYDZ(-1, 1, 0), g_trim)) b_partie = geompy.MakePartition([b_cyl], b_tools, [], [], geompy.ShapeType["SOLID"]) b_element = geompy.SubShapeAll(b_partie, geompy.ShapeType["SOLID"]) cpd = cpd + b_element # Partie haute du gros cylindre # ----------------------------- h_plan = geompy.MakePlane(g_base, g_coupe, g_trim) cpd.append(geompy.MakeMirrorByPlane(b_boite, h_plan)) for h in b_element: h_symetrie = geompy.MakeMirrorByPlane(h, h_plan) cpd.append(h_symetrie) # Petit cylindre # -------------- z_arete = p_rayon/2 x_arete = z_arete*t_hauteur*2/g_hauteur px = cx-x_arete py = cy-1.5*g_rayon pz = cz+g_hauteur/2 p_base = geompy.MakeVertex(cx, py, pz) p_dir = geompy.MakeVectorDXDYDZ(0, 1, 0) p_cyl = geompy.MakeCylinder(p_base, p_dir, p_rayon, p_hauteur) p_boite = geompy.MakeBox(px, py, pz-z_arete, cx+x_arete, py+p_hauteur, pz+z_arete) # Partie interieure du petit cylindre # ----------------------------------- i_cyl = geompy.MakeCommon(p_cyl, g_cyl) i_tuyau = geompy.MakeCut(i_cyl, p_boite) i_boite = geompy.MakeCommon(p_boite, g_cyl) # Partie exterieure du petit cylindre # ----------------------------------- e_cyl0 = geompy.MakeCut(p_cyl, g_cyl) e_cyl = geompy.SubShapeAllSorted(e_cyl0, geompy.ShapeType["SOLID"]) e_tuyau = geompy.MakeCut(e_cyl[1], p_boite) e_boite0 = geompy.MakeCut(p_boite, g_cyl) e_boite = geompy.SubShapeAllSorted(e_boite0, geompy.ShapeType["SOLID"]) cpd.append(e_boite[1]) # Partie centrale du gros cylindre # -------------------------------- c_cyl = geompy.MakeCut(g_centre, p_cyl) # Partitionner # ------------ p_tools = [] p_tools.append(geompy.MakePlane(geompy.MakeVertex(px, py, pz-z_arete), geompy.MakeVectorDXDYDZ(-z_arete, 0, x_arete), g_trim)) p_tools.append(geompy.MakePlane(geompy.MakeVertex(px, py, pz+z_arete), geompy.MakeVectorDXDYDZ( z_arete, 0, x_arete), g_trim)) p_partie = geompy.MakePartition([e_tuyau], p_tools, [], [], geompy.ShapeType["SOLID"]) p_element = geompy.SubShapeAll(p_partie, geompy.ShapeType["SOLID"]) cpd = cpd + p_element q_partie = geompy.MakePartition([i_tuyau, c_cyl], p_tools, [], [], geompy.ShapeType["SOLID"]) q_element = geompy.SubShapeAll(q_partie, geompy.ShapeType["SOLID"]) q_element = q_element + [i_boite] q_tools = [] q_tools.append(geompy.MakePlane(geompy.MakeVertex(cx, cy-b_base, cz), geompy.MakeVectorDXDYDZ(0, 1, 0), g_trim)) q_tools.append(geompy.MakePlane(geompy.MakeVertex(cx, cy+b_base, cz), geompy.MakeVectorDXDYDZ(0, 1, 0), g_trim)) r_element = [] for e in q_element: r_partie = geompy.MakePartition([e], q_tools, [], [], geompy.ShapeType["SOLID"]) r_element = r_element + geompy.SubShapeAll(r_partie, geompy.ShapeType["SOLID"]) cpd = cpd + r_element # Compound # -------- piece = geompy.RemoveExtraEdges(geompy.MakeCompound(cpd), True) # Ajouter la piece dans l'etude # ----------------------------- piece_id = geompy.addToStudy(piece, "ex15_cyl2geometry") # Meshing # ======= # Create a hexahedral mesh # ------------------------ hexa = smesh.Mesh(piece, "ex15_cyl2geometry:hexa") algo = hexa.Segment() algo.NumberOfSegments(12) hexa.Quadrangle() hexa.Hexahedron() # Mesh calculus # ------------- isDone = hexa.Compute() if not isDone: raise Exception("Error when computing Mesh") # Update object browser # --------------------- salome.sg.updateObjBrowser()