10 from salome.geom import geomBuilder
12 geompy = geomBuilder.New()
24 O = geompy.MakeVertex(0, 0, 0)
25 OX = geompy.MakeVectorDXDYDZ(1, 0, 0)
26 OY = geompy.MakeVectorDXDYDZ(0, 1, 0)
27 OZ = geompy.MakeVectorDXDYDZ(0, 0, 1)
28 geompy.addToStudy( O, 'O' )
29 geompy.addToStudy( OX, 'OX' )
30 geompy.addToStudy( OY, 'OY' )
31 geompy.addToStudy( OZ, 'OZ' )
33 Cylinder_1 = geompy.MakeCylinderRH(r1, h1)
34 Cylinder_2 = geompy.MakeCylinderRH(r2, h2)
35 Rotation_1 = geompy.MakeRotation(Cylinder_2, OY, -90*math.pi/180.0)
36 Translation_1 = geompy.MakeTranslation(Rotation_1, 0, 0, length_piquage)
38 tpipe = geompy.MakeFuseList([Cylinder_1, Translation_1], True, True)
39 geompy.addToStudy( tpipe, 'tpipe' )
41 Inlet_z = geompy.GetFaceNearPoint(tpipe, O)
42 geompy.addToStudyInFather( tpipe, Inlet_z, 'Inlet_z' )
44 p_inlet_x = geompy.MakeVertex(-h2, 0, length_piquage)
45 Inlet_x = geompy.GetFaceNearPoint(tpipe, p_inlet_x)
46 geompy.addToStudyInFather( tpipe, Inlet_x, 'Inlet_x' )
48 p_outlet = geompy.MakeVertex(0, 0, h1)
49 Outlet = geompy.GetFaceNearPoint(tpipe, p_outlet)
50 geompy.addToStudyInFather( tpipe, Outlet, 'Outlet' )
52 Wall = geompy.CreateGroup(tpipe, geompy.ShapeType["FACE"])
53 faces = geompy.SubShapeAll(tpipe, geompy.ShapeType["FACE"])
54 geompy.UnionList(Wall, faces)
55 geompy.DifferenceList(Wall, [Inlet_x, Inlet_z, Outlet])
56 geompy.addToStudyInFather( tpipe, Wall, 'Wall' )
58 p_corner = geompy.MakeVertex(-r2, 0, length_piquage+r2)
59 corner = geompy.GetVertexNearPoint(tpipe, p_corner)
60 geompy.addToStudyInFather( tpipe, corner, 'corner' )
62 geom_groups = [Inlet_x, Inlet_z, Outlet, Wall]
64 volumeGEOM = geompy.BasicProperties(tpipe)[2]
70 import SMESH, SALOMEDS
71 from salome.smesh import smeshBuilder
73 from salome.StdMeshers import StdMeshersBuilder
75 smesh = smeshBuilder.New()
77 # Coarse mesh with default hypothesis
78 Mesh_1 = smesh.Mesh(tpipe, "Mesh_coarse")
80 Mesh_1.Triangle(algo=smeshBuilder.NETGEN_1D2D)
82 Mesh_1.Tetrahedron(algo=smeshBuilder.NETGEN_3D)
83 isDone = Mesh_1.Compute()
87 for geom_group in geom_groups:
88 gr = Mesh_1.Group(geom_group)
89 gr_name = gr.GetName()
90 d_mesh_groups[gr_name] = gr
92 # Check tetra mesh volume
93 mesh_volume = smesh.GetVolume(Mesh_1)
94 shape_volume = geompy.BasicProperties(tpipe)[2]
95 assert abs(mesh_volume-shape_volume)/shape_volume < 0.05
97 dual_Mesh_raw_1 = smesh.CreateDualMesh(Mesh_1, 'dual_Mesh_raw_1', False)
100 assert dual_Mesh_raw_1.NbPolyhedrons() > 0
101 assert dual_Mesh_raw_1.NbPolyhedrons() == dual_Mesh_raw_1.NbVolumes()
103 # Check dual mesh volume
104 dual_raw_volume = dual_Mesh_raw_1.GetVolume()
105 assert abs(dual_raw_volume-shape_volume)/shape_volume < 0.11
108 dual_Mesh_raw_groups = dual_Mesh_raw_1.GetGroups()
109 dual_Mesh_raw_group_names = dual_Mesh_raw_1.GetGroupNames()
111 assert len(dual_Mesh_raw_groups) == 4
113 for gr_dual, gr_name in zip(dual_Mesh_raw_groups, dual_Mesh_raw_group_names):
114 gr_name = gr_name.strip()
115 gr_tri = d_mesh_groups[gr_name]
116 area_gr_tri = smesh.GetArea(gr_tri)
117 area_gr_dual = smesh.GetArea(gr_dual)
119 print("Area tri: ", area_gr_tri)
120 print("Area dual:", area_gr_dual)
121 assert abs(area_gr_tri-area_gr_dual)/area_gr_tri < 1e-3
123 #dual_Mesh_1 = smesh.CreateDualMesh(Mesh_1, 'dual_Mesh_1', True)
126 #dual_volume = dual_Mesh_1.GetVolume()
127 #dual_raw_volume = dual_Mesh_raw_1.GetVolume()
128 #print("dual_volume: ", dual_volume)
129 #print("dual_raw_volume: ", dual_raw_volume)
131 #assert (dual_volume >= dual_raw_volume)
133 if salome.sg.hasDesktop():
134 salome.sg.updateObjBrowser()