[python] ffff
Viewer
*** This page was generated with the meta tag "noindex, nofollow". This happened because you selected this option before saving or the system detected it as spam. This means that this page will never get into the search engines and the search bot will not crawl it. There is nothing to worry about, you can still share it with anyone.
- #!/usr/bin/env python
- ###
- ### This file is generated automatically by SALOME v9.9.0 with dump python functionality
- ###
- import sys
- import salome
- import math
- salome.salome_init()
- import salome_notebook
- notebook = salome_notebook.NoteBook()
- sys.path.insert(0, r'/home/mika/OpenFOAM/mika-9/run/nozzle_1/geometry/')
- # входные параметры
- P_input = 200000
- T_input = 1000
- P_output = 8000
- G = 2
- alpha_ = 14
- betta_ = 22
- R = 287
- k = 1.4
- # рассчет входного сечения
- betta_vh = 0.999
- # давление
- p_vh = betta_vh*P_input
- # температура
- T_vh = T_input*betta_vh**((k-1)/k)
- # удельный объем
- mu_vh = R*T_vh/p_vh
- # плотность
- ro_vh = 1/mu_vh
- # скорость потока
- c_vh = math.sqrt(2*k/(k-1)*R*(T_input-T_vh))
- # местная скорость звука
- a_vh = math.sqrt(k*R*T_vh)
- # площадь входного поперечного сечения
- F_vh = G/ro_vh/c_vh
- # диаметр входного сечения сопла
- d_vh = round(math.sqrt(4/3.14*F_vh)*1e3)
- betta_kr = (2/(k+1))**(k/(k-1))
- p_kr = betta_kr*P_input
- T_kr = T_input*betta_kr**((k-1)/k)
- mu_kr = R*T_kr/p_kr
- ro_kr = 1/mu_kr
- c_kr = math.sqrt(2*k/(k-1)*R*(T_input-T_kr))
- a_kr = math.sqrt(k*R*T_kr)
- F_kr = G/ro_kr/c_kr
- d_kr = round(math.sqrt(4/3.14*F_kr)*1e3)
- # расчет выходного сечения
- betta_out = P_output/P_input
- # давление
- p_out = betta_out*P_input
- # Temnepatypa
- T_out = T_input*betta_out**((k-1)/k)
- # удельный объем
- mu_out = R*T_out/p_out
- # плотность
- ro_out = 1/mu_out
- # скорость потока
- c_out = math.sqrt(2*k/(k-1)*R*(T_input-T_out))
- # местная скорость звука
- a_out = math.sqrt(k*R*T_out)
- # площадь выходного поперечного сечения
- F_out = G/ro_out/c_out
- # диаметр выходного сечения сопла
- d_out = round(math.sqrt(4/3.14*F_out)*1e3)
- # расчетная длина сопла до критического сечения и после
- l_left = (d_vh-d_kr)/(2*math.tan(alpha_/2*math.pi/180))
- l_right = (d_out-d_kr)/(2*math.tan(betta_/2*math.pi/180))
- d_inlet_ = d_vh
- d_kr_ = d_kr
- d_outlet_ = d_out
- r_ = d_kr_
- OS_ = r_+d_kr_/2
- l1_ = (l_left+l_right)*1.6
- l2_ = d_outlet_/5
- h_ = d_outlet_*4
- arc_left = math.pi * alpha_ * r_ / 180
- arc_right = math.pi * betta_ * r_ / 180
- # atmosphere height above nozzle
- k_height = 177*6/885
- # atmosphere width
- k_width = (1189.09+152.41+2146.4)/2146.4
- # left top atmosphere
- k_left_top = (2146.4+30.482)/30.482
- # nozzle width
- k_nozzle = (1189.09+152.41)/1189.09
- # arc left
- k_arc_left = (14*3.14*100)/180/628
- # arc right
- k_arc_right = (28*3.14*100)/180/628
- NbSg1 = 15
- NbSg2 = 1
- NbSg3 = round(l_left*k_nozzle/13.4)
- NbSg4 = round(arc_left * k_arc_left*2)
- NbSg5 = round(arc_right * k_arc_right)
- NbSg6 = round(l_right*k_nozzle/5.7)
- NbSg7 = round(l1_*k_width/7)
- NbSg8 = round(l2_*k_left_top/360)
- NbSg9 = round(h_*k_height/15)
- ###
- ### SHAPER component
- ###
- from salome.shaper import model
- model.begin()
- partSet = model.moduleDocument()
- ### Create Part
- Part_1 = model.addPart(partSet)
- Part_1_doc = Part_1.document()
- model.addParameter(Part_1_doc, "d_inlet", "d_inlet_")
- model.addParameter(Part_1_doc, "d_kr", "d_kr_")
- model.addParameter(Part_1_doc, "d_outlet", "d_outlet_")
- model.addParameter(Part_1_doc, "alpha", "alpha_")
- model.addParameter(Part_1_doc, "betta", "betta_")
- model.addParameter(Part_1_doc, "r", "r_")
- model.addParameter(Part_1_doc, "OS", "OS_")
- model.addParameter(Part_1_doc, "l1", "l1_")
- model.addParameter(Part_1_doc, "l2", "l2_")
- model.addParameter(Part_1_doc, "h", "h_")
- model.addParameter(Part_1_doc, "l_left", "l_left")
- model.addParameter(Part_1_doc, "l_right", "l_right")
- ### Create Sketch
- Sketch_1 = model.addSketch(Part_1_doc, model.defaultPlane("XOY"))
- ### Create SketchLine
- SketchLine_1 = Sketch_1.addLine(0, 0, -12.18693434051474, 0)
- ### Create SketchProjection
- SketchProjection_1 = Sketch_1.addProjection(model.selection("VERTEX", "PartSet/Origin"), False)
- SketchPoint_1 = SketchProjection_1.createdFeature()
- Sketch_1.setCoincident(SketchLine_1.startPoint(), SketchPoint_1.result())
- Sketch_1.setHorizontal(SketchLine_1.result())
- ### Create SketchLine
- SketchLine_2 = Sketch_1.addLine(-12.18693434051474, 0, -1195.19084049934, 0)
- Sketch_1.setCoincident(SketchLine_1.endPoint(), SketchLine_2.startPoint())
- Sketch_1.setHorizontal(SketchLine_2.result())
- ### Create SketchLine
- SketchLine_3 = Sketch_1.addLine(0, 49.99999999999998, 0, 0)
- ### Create SketchProjection
- SketchProjection_2 = Sketch_1.addProjection(model.selection("EDGE", "PartSet/OY"), False)
- SketchLine_4 = SketchProjection_2.createdFeature()
- Sketch_1.setCoincident(SketchLine_3.startPoint(), SketchLine_4.result())
- Sketch_1.setCoincident(SketchLine_1.startPoint(), SketchLine_3.endPoint())
- ### Create SketchLine
- SketchLine_5 = Sketch_1.addLine(24.19218955996678, 0, 0, 0)
- ### Create SketchProjection
- SketchProjection_3 = Sketch_1.addProjection(model.selection("EDGE", "PartSet/OX"), False)
- SketchLine_6 = SketchProjection_3.createdFeature()
- Sketch_1.setCoincident(SketchLine_5.startPoint(), SketchLine_6.result())
- Sketch_1.setCoincident(SketchLine_1.startPoint(), SketchLine_5.endPoint())
- ### Create SketchLine
- SketchLine_7 = Sketch_1.addLine(164.6881315646525, 0, 24.19218955996678, 0)
- Sketch_1.setCoincident(SketchLine_5.startPoint(), SketchLine_7.endPoint())
- Sketch_1.setHorizontal(SketchLine_7.result())
- ### Create SketchLine
- SketchLine_8 = Sketch_1.addLine(-1195.19084049934, 196, -1195.19084049934, 0)
- Sketch_1.setCoincident(SketchLine_2.endPoint(), SketchLine_8.endPoint())
- Sketch_1.setVertical(SketchLine_8.result())
- ### Create SketchLine
- SketchLine_9 = Sketch_1.addLine(-1195.19084049934, 196, -12.18693434051474, 50.74538483586778)
- Sketch_1.setCoincident(SketchLine_8.startPoint(), SketchLine_9.startPoint())
- ### Create SketchLine
- SketchLine_10 = Sketch_1.addLine(-12.18693434051474, 50.74538483586778, -12.18693434051474, 0)
- Sketch_1.setCoincident(SketchLine_9.endPoint(), SketchLine_10.startPoint())
- Sketch_1.setCoincident(SketchLine_1.endPoint(), SketchLine_10.endPoint())
- Sketch_1.setVertical(SketchLine_10.result())
- ### Create SketchLine
- SketchLine_11 = Sketch_1.addLine(24.19218955996678, 52.97042737240035, 24.19218955996678, 0)
- Sketch_1.setCoincident(SketchLine_5.startPoint(), SketchLine_11.endPoint())
- Sketch_1.setVertical(SketchLine_11.result())
- ### Create SketchLine
- SketchLine_12 = Sketch_1.addLine(164.6881315646525, 88, 164.6881315646525, 0)
- Sketch_1.setCoincident(SketchLine_7.startPoint(), SketchLine_12.endPoint())
- Sketch_1.setVertical(SketchLine_12.result())
- ### Create SketchLine
- SketchLine_13 = Sketch_1.addLine(164.6881315646525, 88, 24.19218955996678, 52.97042737240035)
- Sketch_1.setCoincident(SketchLine_12.startPoint(), SketchLine_13.startPoint())
- Sketch_1.setCoincident(SketchLine_11.startPoint(), SketchLine_13.endPoint())
- ### Create SketchConstraintAngle
- Sketch_1.setAngle(SketchLine_9.result(), SketchLine_2.result(), "alpha/2", type = "Direct")
- ### Create SketchLine
- SketchLine_14 = Sketch_1.addLine(0, 49.99999999999998, 0, 150)
- SketchLine_14.setAuxiliary(True)
- Sketch_1.setCoincident(SketchLine_3.startPoint(), SketchLine_14.startPoint())
- Sketch_1.setCoincident(SketchLine_14.endPoint(), SketchLine_4.result())
- Sketch_1.setVerticalDistance(SketchLine_14.endPoint(), SketchLine_1.startPoint(), "OS")
- Sketch_1.setVerticalDistance(SketchLine_14.endPoint(), SketchLine_3.startPoint(), "r")
- ### Create SketchConstraintAngle
- Sketch_1.setAngle(SketchLine_13.result(), SketchLine_7.result(), "betta/2", type = "Direct")
- Sketch_1.setVerticalDistance(SketchLine_13.startPoint(), SketchLine_7.startPoint(), "d_outlet/2")
- Sketch_1.setVerticalDistance(SketchLine_8.startPoint(), SketchLine_8.endPoint(), "d_inlet/2")
- ### Create SketchArc
- SketchArc_1 = Sketch_1.addArc(-1.102898468723179e-33, 150, 0, 49.99999999999998, 24.19218955996678, 52.97042737240035, False)
- Sketch_1.setCoincident(SketchLine_3.startPoint(), SketchArc_1.startPoint())
- Sketch_1.setCoincident(SketchLine_11.startPoint(), SketchArc_1.endPoint())
- ### Create SketchArc
- SketchArc_2 = Sketch_1.addArc(-1.102898468723179e-33, 150, 0, 49.99999999999998, -12.18693434051474, 50.74538483586778, True)
- Sketch_1.setCoincident(SketchLine_3.startPoint(), SketchArc_2.startPoint())
- Sketch_1.setCoincident(SketchLine_9.endPoint(), SketchArc_2.endPoint())
- Sketch_1.setTangent(SketchArc_2.results()[1], SketchLine_9.result())
- Sketch_1.setTangent(SketchArc_1.results()[1], SketchLine_13.result())
- Sketch_1.setCoincident(SketchArc_2.center(), SketchArc_1.center())
- Sketch_1.setRadius(SketchArc_1.results()[1], "r")
- Sketch_1.setCoincident(SketchArc_1.center(), SketchLine_14.result())
- ### Create SketchLine
- SketchLine_15 = Sketch_1.addLine(164.6881315646525, 88, 164.6881315627929, 792)
- Sketch_1.setCoincident(SketchLine_12.startPoint(), SketchLine_15.startPoint())
- Sketch_1.setVertical(SketchLine_15.result())
- ### Create SketchLine
- SketchLine_16 = Sketch_1.addLine(164.6881315627929, 792, 129.4881315646525, 792)
- Sketch_1.setCoincident(SketchLine_15.endPoint(), SketchLine_16.startPoint())
- ### Create SketchLine
- SketchLine_17 = Sketch_1.addLine(129.4881315646525, 792, 129.4881315646525, 88)
- Sketch_1.setCoincident(SketchLine_16.endPoint(), SketchLine_17.startPoint())
- Sketch_1.setVertical(SketchLine_17.result())
- ### Create SketchLine
- SketchLine_18 = Sketch_1.addLine(129.4881315646525, 88, 164.6881315646525, 88)
- Sketch_1.setCoincident(SketchLine_17.endPoint(), SketchLine_18.startPoint())
- Sketch_1.setCoincident(SketchLine_12.startPoint(), SketchLine_18.endPoint())
- Sketch_1.setHorizontal(SketchLine_18.result())
- ### Create SketchLine
- SketchLine_19 = Sketch_1.addLine(164.6881315646525, 0, 2311.058131564647, 0)
- Sketch_1.setCoincident(SketchLine_7.startPoint(), SketchLine_19.startPoint())
- Sketch_1.setHorizontal(SketchLine_19.result())
- ### Create SketchLine
- SketchLine_20 = Sketch_1.addLine(2311.058131564647, 0, 2311.058131564647, 792)
- Sketch_1.setCoincident(SketchLine_19.endPoint(), SketchLine_20.startPoint())
- Sketch_1.setVertical(SketchLine_20.result())
- ### Create SketchLine
- SketchLine_21 = Sketch_1.addLine(2311.058131564647, 792, 164.6881315627929, 792)
- Sketch_1.setCoincident(SketchLine_20.endPoint(), SketchLine_21.startPoint())
- Sketch_1.setCoincident(SketchLine_15.endPoint(), SketchLine_21.endPoint())
- Sketch_1.setHorizontal(SketchLine_21.result())
- Sketch_1.setHorizontal(SketchLine_16.result())
- ### Create SketchLine
- SketchLine_22 = Sketch_1.addLine(164.6881315646525, 88, 2311.058131564647, 88)
- Sketch_1.setCoincident(SketchLine_12.startPoint(), SketchLine_22.startPoint())
- Sketch_1.setCoincident(SketchLine_22.endPoint(), SketchLine_20.result())
- Sketch_1.setHorizontal(SketchLine_22.result())
- Sketch_1.setVerticalDistance(SketchLine_17.startPoint(), SketchLine_18.startPoint(), "h")
- Sketch_1.setHorizontalDistance(SketchLine_15.endPoint(), SketchLine_21.startPoint(), "l1")
- Sketch_1.setHorizontalDistance(SketchLine_15.endPoint(), SketchLine_17.startPoint(), "l2")
- model.do()
- Sketch_1.changeFacesOrder([[SketchLine_3.result(), SketchLine_5.result(), SketchLine_11.result(), SketchArc_1.results()[1]],
- [SketchLine_3.result(), SketchArc_2.results()[1], SketchLine_10.result(), SketchLine_1.result()],
- [SketchLine_12.result(), SketchLine_19.result(), SketchLine_20.result(), SketchLine_22.result()],
- [SketchLine_10.result(), SketchLine_9.result(), SketchLine_8.result(), SketchLine_2.result()],
- [SketchLine_12.result(), SketchLine_13.result(), SketchLine_11.result(), SketchLine_7.result()],
- [SketchLine_15.result(), SketchLine_16.result(), SketchLine_17.result(), SketchLine_18.result()],
- [SketchLine_22.result(), SketchLine_20.result(), SketchLine_21.result(), SketchLine_15.result()]
- ])
- model.do()
- ### Create Face
- Face_1_objects = [model.selection("EDGE", "Sketch_1/SketchLine_9"),
- model.selection("EDGE", "Sketch_1/SketchLine_8"),
- model.selection("EDGE", "Sketch_1/SketchLine_2"),
- model.selection("EDGE", "Sketch_1/SketchLine_10")]
- Face_1 = model.addFace(Part_1_doc, Face_1_objects)
- ### Create Face
- Face_2_objects = [model.selection("EDGE", "Sketch_1/SketchArc_2_2"),
- model.selection("EDGE", "Sketch_1/SketchLine_3"),
- model.selection("EDGE", "Sketch_1/SketchLine_1"),
- model.selection("EDGE", "Sketch_1/SketchLine_10")]
- Face_2 = model.addFace(Part_1_doc, Face_2_objects)
- ### Create Face
- Face_3 = model.addFace(Part_1_doc, [model.selection("FACE", "Sketch_1/Face-SketchLine_3f-SketchLine_5r-SketchLine_11r-SketchArc_1_2r")])
- ### Create Face
- Face_4 = model.addFace(Part_1_doc, [model.selection("FACE", "Sketch_1/Face-SketchLine_12r-SketchLine_13f-SketchLine_11f-SketchLine_7r")])
- ### Create Face
- Face_5 = model.addFace(Part_1_doc, [model.selection("FACE", "Sketch_1/Face-SketchLine_22f-SketchLine_20f-SketchLine_21f-SketchLine_15r")])
- ### Create Face
- Face_6 = model.addFace(Part_1_doc, [model.selection("FACE", "Sketch_1/Face-SketchLine_15f-SketchLine_16f-SketchLine_17f-SketchLine_18f")])
- ### Create Face
- Face_7 = model.addFace(Part_1_doc, [model.selection("FACE", "Sketch_1/Face-SketchLine_12f-SketchLine_19f-SketchLine_20f-SketchLine_22r")])
- ### Create Shell
- Shell_1_objects = [model.selection("FACE", "Face_1_1"),
- model.selection("FACE", "Sketch_1/Face-SketchLine_3r-SketchArc_2_2r-SketchLine_10f-SketchLine_1r"),
- model.selection("FACE", "Face_3_1"),
- model.selection("FACE", "Face_4_1"),
- model.selection("FACE", "Sketch_1/Face-SketchLine_12f-SketchLine_19f-SketchLine_20f-SketchLine_22r"),
- model.selection("FACE", "Face_5_1"),
- model.selection("FACE", "Face_6_1")]
- Shell_1 = model.addShell(Part_1_doc, Shell_1_objects)
- ### Create Revolution
- Revolution_1 = model.addRevolution(Part_1_doc, [model.selection("SHELL", "Shell_1_1")], model.selection("EDGE", "PartSet/OX"), 2.5, 2.5, "Faces")
- ### Create Group
- Group_1_objects = [model.selection("EDGE", "[Revolution_1_1_1/Generated_Face&Sketch_1/SketchLine_10][Revolution_1_1_1/Generated_Face&Sketch_1/SketchLine_9]"),
- model.selection("EDGE", "[Revolution_1_1_2/Generated_Face&Sketch_1/SketchLine_3][Revolution_1_1_2/Generated_Face&Sketch_1/SketchArc_2_2]"),
- model.selection("EDGE", "[Revolution_1_1_1/Generated_Face&Sketch_1/SketchLine_8][Revolution_1_1_1/Generated_Face&Sketch_1/SketchLine_9]"),
- model.selection("EDGE", "[Revolution_1_1_3/Generated_Face&Sketch_1/SketchLine_11][Revolution_1_1_3/Generated_Face&Sketch_1/SketchArc_1_2]"),
- model.selection("EDGE", "[Revolution_1_1_4/Generated_Face&Sketch_1/SketchLine_12][Revolution_1_1_4/Generated_Face&Sketch_1/SketchLine_13]"),
- model.selection("EDGE", "[Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_17][Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_18]"),
- model.selection("EDGE", "[Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_16][Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_17]"),
- model.selection("EDGE", "[Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_15][Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_16]"),
- model.selection("EDGE", "[Revolution_1_1_6/Generated_Face&Sketch_1/SketchLine_20][Revolution_1_1_6/Generated_Face&Sketch_1/SketchLine_21]"),
- model.selection("EDGE", "[Revolution_1_1_5/Generated_Face&Sketch_1/SketchLine_20][Revolution_1_1_5/Generated_Face&Sketch_1/SketchLine_22]")]
- Group_1 = model.addGroup(Part_1_doc, "Edges", Group_1_objects)
- ### Create Group
- Group_2_objects = [model.selection("FACE", "Revolution_1_1_1/Generated_Face&Sketch_1/SketchLine_8"),
- model.selection("FACE", "Revolution_1_1_2/Generated_Face&Sketch_1/SketchLine_10"),
- model.selection("FACE", "Revolution_1_1_3/Generated_Face&Sketch_1/SketchLine_3"),
- model.selection("FACE", "Revolution_1_1_4/Generated_Face&Sketch_1/SketchLine_11"),
- model.selection("FACE", "Revolution_1_1_5/Generated_Face&Sketch_1/SketchLine_12"),
- model.selection("FACE", "Revolution_1_1_5/Generated_Face&Sketch_1/SketchLine_20")]
- Group_2 = model.addGroup(Part_1_doc, "Faces", Group_2_objects)
- ### Create Group
- Group_3_objects = [model.selection("EDGE", "[Revolution_1_1_1/Generated_Face&Sketch_1/SketchLine_9][Revolution_1_1_1/From_Face]"),
- model.selection("EDGE", "[Revolution_1_1_1/Generated_Face&Sketch_1/SketchLine_9][Revolution_1_1_1/To_Face]"),
- model.selection("EDGE", "[Revolution_1_1_1/From_Face][Revolution_1_1_1/To_Face]")]
- Group_3 = model.addGroup(Part_1_doc, "Edges", Group_3_objects)
- ### Create Group
- Group_4_objects = [model.selection("EDGE", "[Revolution_1_1_2/Generated_Face&Sketch_1/SketchArc_2_2][Revolution_1_1_2/From_Face]"),
- model.selection("EDGE", "[Revolution_1_1_2/From_Face][Revolution_1_1_2/To_Face]"),
- model.selection("EDGE", "[Revolution_1_1_2/Generated_Face&Sketch_1/SketchArc_2_2][Revolution_1_1_2/To_Face]")]
- Group_4 = model.addGroup(Part_1_doc, "Edges", Group_4_objects)
- ### Create Group
- Group_5_objects = [model.selection("EDGE", "[Revolution_1_1_3/Generated_Face&Sketch_1/SketchArc_1_2][Revolution_1_1_3/From_Face]"),
- model.selection("EDGE", "[Revolution_1_1_3/Generated_Face&Sketch_1/SketchArc_1_2][Revolution_1_1_3/To_Face]"),
- model.selection("EDGE", "[Revolution_1_1_3/From_Face][Revolution_1_1_3/To_Face]")]
- Group_5 = model.addGroup(Part_1_doc, "Edges", Group_5_objects)
- ### Create Group
- Group_6_objects = [model.selection("EDGE", "[Revolution_1_1_4/Generated_Face&Sketch_1/SketchLine_13][Revolution_1_1_4/From_Face]"),
- model.selection("EDGE", "[Revolution_1_1_4/Generated_Face&Sketch_1/SketchLine_13][Revolution_1_1_4/To_Face]"),
- model.selection("EDGE", "[Revolution_1_1_4/From_Face][Revolution_1_1_4/To_Face]")]
- Group_6 = model.addGroup(Part_1_doc, "Edges", Group_6_objects)
- ### Create Group
- Group_7_objects = [model.selection("EDGE", "[Revolution_1_1_6/Generated_Face&Sketch_1/SketchLine_21][Revolution_1_1_6/From_Face]"),
- model.selection("EDGE", "[Revolution_1_1_5/Generated_Face&Sketch_1/SketchLine_22][Revolution_1_1_5/From_Face]"),
- model.selection("EDGE", "Face_7_1/Modified_Edge&Sketch_1/SketchLine_22"),
- model.selection("EDGE", "[Revolution_1_1_5/From_Face][Revolution_1_1_5/To_Face]"),
- model.selection("EDGE", "[Revolution_1_1_6/Generated_Face&Sketch_1/SketchLine_21][Revolution_1_1_6/To_Face]")]
- Group_7 = model.addGroup(Part_1_doc, "Edges", Group_7_objects)
- ### Create Group
- Group_8_objects = [model.selection("EDGE", "[Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_16][Revolution_1_1_7/From_Face]"),
- model.selection("EDGE", "[Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_16][Revolution_1_1_7/To_Face]"),
- model.selection("EDGE", "[Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_18][Revolution_1_1_7/From_Face]"),
- model.selection("EDGE", "[Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_18][Revolution_1_1_7/To_Face]")]
- Group_8 = model.addGroup(Part_1_doc, "Edges", Group_8_objects)
- ### Create Group
- Group_9_objects = [model.selection("EDGE", "[Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_17][Revolution_1_1_7/From_Face]"),
- model.selection("EDGE", "[Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_17][Revolution_1_1_7/To_Face]"),
- model.selection("EDGE", "[Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_15][Revolution_1_1_7/From_Face]"),
- model.selection("EDGE", "[Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_15][Revolution_1_1_7/To_Face]"),
- model.selection("EDGE", "[Revolution_1_1_6/Generated_Face&Sketch_1/SketchLine_20][Revolution_1_1_6/From_Face]"),
- model.selection("EDGE", "[Revolution_1_1_6/Generated_Face&Sketch_1/SketchLine_20][Revolution_1_1_6/To_Face]")]
- Group_9 = model.addGroup(Part_1_doc, "Edges", Group_9_objects)
- ### Create Group
- Group_10_objects = [model.selection("FACE", "Revolution_1_1_1/To_Face"),
- model.selection("FACE", "Revolution_1_1_3/To_Face"),
- model.selection("FACE", "Revolution_1_1_2/To_Face"),
- model.selection("FACE", "Revolution_1_1_4/To_Face"),
- model.selection("FACE", "Revolution_1_1_5/To_Face"),
- model.selection("FACE", "Revolution_1_1_6/To_Face"),
- model.selection("FACE", "Revolution_1_1_7/To_Face")]
- Group_10 = model.addGroup(Part_1_doc, "Faces", Group_10_objects)
- Group_10.setName("front")
- Group_10.result().setName("front")
- ### Create Group
- Group_11_objects = [model.selection("FACE", "Revolution_1_1_1/From_Face"),
- model.selection("FACE", "Revolution_1_1_2/From_Face"),
- model.selection("FACE", "Revolution_1_1_3/From_Face"),
- model.selection("FACE", "Revolution_1_1_4/From_Face"),
- model.selection("FACE", "Revolution_1_1_7/From_Face"),
- model.selection("FACE", "Revolution_1_1_6/From_Face")]
- Group_11 = model.addGroup(Part_1_doc, "Faces", Group_11_objects)
- Group_11.setName("back")
- Group_11.result().setName("back")
- ### Create Group
- Group_12_objects = [model.selection("FACE", "Revolution_1_1_1/Generated_Face&Sketch_1/SketchLine_9"),
- model.selection("FACE", "Revolution_1_1_2/Generated_Face&Sketch_1/SketchArc_2_2"),
- model.selection("FACE", "Revolution_1_1_3/Generated_Face&Sketch_1/SketchArc_1_2"),
- model.selection("FACE", "Revolution_1_1_4/Generated_Face&Sketch_1/SketchLine_13"),
- model.selection("FACE", "Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_18")]
- Group_12 = model.addGroup(Part_1_doc, "Faces", Group_12_objects)
- Group_12.setName("wall")
- Group_12.result().setName("wall")
- ### Create Group
- Group_13 = model.addGroup(Part_1_doc, "Faces", [model.selection("FACE", "Revolution_1_1_1/Generated_Face&Sketch_1/SketchLine_8")])
- Group_13.setName("inlet")
- Group_13.result().setName("inlet")
- ### Create Group
- Group_14 = model.addGroup(Part_1_doc, "Faces", [model.selection("FACE", "Revolution_1_1_7/Generated_Face&Sketch_1/SketchLine_17")])
- Group_14.setName("outletLeft")
- Group_14.result().setName("outletLeft")
- ### Create Group
- Group_15 = model.addGroup(Part_1_doc, "Faces", [model.selection("FACE", "Revolution_1_1_6/Generated_Face&Sketch_1/SketchLine_20"), model.selection("FACE", "Revolution_1_1_5/Generated_Face&Sketch_1/SketchLine_20")])
- Group_15.setName("outletRight")
- Group_15.result().setName("outletRight")
- model.end()
- ###
- ### SHAPERSTUDY component
- ###
- model.publishToShaperStudy()
- import SHAPERSTUDY
- Face_2_1, = SHAPERSTUDY.shape(model.featureStringId(Face_2))
- Face_7_1, Group_7_1, = SHAPERSTUDY.shape(model.featureStringId(Face_7))
- Revolution_1_1, Group_1_1, Group_2_1, Group_3_1, Group_4_1, Group_5_1, Group_6_1, Group_7_2, Group_8_1, Group_9_1, front, back, wall, inlet, outletLeft, outletRight, = SHAPERSTUDY.shape(model.featureStringId(Revolution_1))
- ###
- ### GEOM component
- ###
- import GEOM
- from salome.geom import geomBuilder
- import math
- import SALOMEDS
- geompy = geomBuilder.New()
- O = geompy.MakeVertex(0, 0, 0)
- OX = geompy.MakeVectorDXDYDZ(1, 0, 0)
- OY = geompy.MakeVectorDXDYDZ(0, 1, 0)
- OZ = geompy.MakeVectorDXDYDZ(0, 0, 1)
- geompy.addToStudy( O, 'O' )
- geompy.addToStudy( OX, 'OX' )
- geompy.addToStudy( OY, 'OY' )
- geompy.addToStudy( OZ, 'OZ' )
- ###
- ### SMESH component
- ###
- import SMESH, SALOMEDS
- from salome.smesh import smeshBuilder
- smesh = smeshBuilder.New()
- #smesh.SetEnablePublish( False ) # Set to False to avoid publish in study if not needed or in some particular situations:
- # multiples meshes built in parallel, complex and numerous mesh edition (performance)
- Mesh_1 = smesh.Mesh(Revolution_1_1,'Mesh_1')
- Regular_1D = Mesh_1.Segment()
- Number_of_Segments_1 = Regular_1D.NumberOfSegments(NbSg1, None, [])
- Prism_3D = Mesh_1.Prism()
- front_1 = Mesh_1.GroupOnGeom(front,'front',SMESH.FACE)
- back_1 = Mesh_1.GroupOnGeom(back,'back',SMESH.FACE)
- wall_1 = Mesh_1.GroupOnGeom(wall,'wall',SMESH.FACE)
- inlet_1 = Mesh_1.GroupOnGeom(inlet,'inlet',SMESH.FACE)
- outletLeft_1 = Mesh_1.GroupOnGeom(outletLeft,'outletLeft',SMESH.FACE)
- outletRight_1 = Mesh_1.GroupOnGeom(outletRight,'outletRight',SMESH.FACE)
- Regular_1D_1 = Mesh_1.Segment(geom=Group_1_1)
- Sub_mesh_1 = Regular_1D_1.GetSubMesh()
- Number_of_Segments_2 = Regular_1D_1.NumberOfSegments(NbSg2, None, [])
- Regular_1D_2 = Mesh_1.Segment(geom=Group_2_1)
- Sub_mesh_2 = Regular_1D_2.GetSubMesh()
- status = Mesh_1.AddHypothesis(Number_of_Segments_1,Group_2_1)
- RadialQuadrangle_1D2D = Mesh_1.Quadrangle(algo=smeshBuilder.RADIAL_QUAD,geom=Group_2_1)
- isDone = Mesh_1.SetMeshOrder( [ [ Sub_mesh_1, Sub_mesh_2 ] ])
- [ front_1, back_1, wall_1, inlet_1, outletLeft_1, outletRight_1 ] = Mesh_1.GetGroups()
- Regular_1D_3 = Mesh_1.Segment(geom=Group_3_1)
- Number_of_Segments_3 = Regular_1D_3.NumberOfSegments(NbSg3,None,[])
- Number_of_Segments_3.SetConversionMode( 0 )
- Number_of_Segments_3.SetExpressionFunction( '0.5*t' )
- [ front_1, back_1, wall_1, inlet_1, outletLeft_1, outletRight_1 ] = Mesh_1.GetGroups()
- Regular_1D_4 = Mesh_1.Segment(geom=Group_4_1)
- Number_of_Segments_4 = Regular_1D_4.NumberOfSegments(NbSg4, None, [])
- Regular_1D_5 = Mesh_1.Segment(geom=Group_5_1)
- Number_of_Segments_5 = Regular_1D_5.NumberOfSegments(NbSg5, None, [])
- Regular_1D_6 = Mesh_1.Segment(geom=Group_6_1)
- Number_of_Segments_6 = Regular_1D_6.NumberOfSegments(NbSg6, None, [])
- [ front_1, back_1, wall_1, inlet_1, outletLeft_1, outletRight_1 ] = Mesh_1.GetGroups()
- Regular_1D_7 = Mesh_1.Segment(geom=Group_7_2)
- Number_of_Segments_7 = Regular_1D_7.NumberOfSegments(NbSg7, None, [])
- [ front_1, back_1, wall_1, inlet_1, outletLeft_1, outletRight_1 ] = Mesh_1.GetGroups()
- Regular_1D_8 = Mesh_1.Segment(geom=Group_8_1)
- Number_of_Segments_8 = Regular_1D_8.NumberOfSegments(NbSg8, None, [])
- [ front_1, back_1, wall_1, inlet_1, outletLeft_1, outletRight_1 ] = Mesh_1.GetGroups()
- Regular_1D_9 = Mesh_1.Segment(geom=Group_9_1)
- Number_of_Segments_9 = Regular_1D_9.NumberOfSegments(NbSg9,None,[])
- Number_of_Segments_9.SetConversionMode( 0 )
- Number_of_Segments_9.SetExpressionFunction( '-t*t' )
- [ front_1, back_1, wall_1, inlet_1, outletLeft_1, outletRight_1 ] = Mesh_1.GetGroups()
- isDone = Mesh_1.Compute()
- [ front_1, back_1, wall_1, inlet_1, outletLeft_1, outletRight_1 ] = Mesh_1.GetGroups()
- Sub_mesh_3 = Regular_1D_3.GetSubMesh()
- Sub_mesh_4 = Regular_1D_4.GetSubMesh()
- Sub_mesh_5 = Regular_1D_5.GetSubMesh()
- Sub_mesh_6 = Regular_1D_6.GetSubMesh()
- Sub_mesh_7 = Regular_1D_7.GetSubMesh()
- Sub_mesh_8 = Regular_1D_8.GetSubMesh()
- Sub_mesh_9 = Regular_1D_9.GetSubMesh()
- ## Set names of Mesh objects
- smesh.SetName(front_1, 'front')
- smesh.SetName(back_1, 'back')
- #smesh.SetName(RadialQuadrangle_1D2D.GetAlgorithm(), 'RadialQuadrangle_1D2D')
- #smesh.SetName(Prism_3D.GetAlgorithm(), 'Prism_3D')
- #smesh.SetName(Regular_1D.GetAlgorithm(), 'Regular_1D')
- #smesh.SetName(Mesh_1.GetMesh(), 'Mesh_1')
- #smesh.SetName(Number_of_Segments_4, 'Number of Segments_4')
- #smesh.SetName(Number_of_Segments_5, 'Number of Segments_5')
- #smesh.SetName(Number_of_Segments_6, 'Number of Segments_6')
- #smesh.SetName(Number_of_Segments_7, 'Number of Segments_7')
- #smesh.SetName(Number_of_Segments_1, 'Number of Segments_1')
- #smesh.SetName(Number_of_Segments_2, 'Number of Segments_2')
- #smesh.SetName(Number_of_Segments_3, 'Number of Segments_3')
- #smesh.SetName(Number_of_Segments_8, 'Number of Segments_8')
- #smesh.SetName(Number_of_Segments_9, 'Number of Segments_9')
- #smesh.SetName(Sub_mesh_9, 'Sub-mesh_9')
- #smesh.SetName(Sub_mesh_8, 'Sub-mesh_8')
- #smesh.SetName(Sub_mesh_5, 'Sub-mesh_5')
- #smesh.SetName(Sub_mesh_4, 'Sub-mesh_4')
- #smesh.SetName(Sub_mesh_7, 'Sub-mesh_7')
- #smesh.SetName(Sub_mesh_6, 'Sub-mesh_6')
- #smesh.SetName(Sub_mesh_1, 'Sub-mesh_1')
- smesh.SetName(outletLeft_1, 'outletLeft')
- smesh.SetName(outletRight_1, 'outletRight')
- #smesh.SetName(Sub_mesh_3, 'Sub-mesh_3')
- smesh.SetName(wall_1, 'wall')
- #smesh.SetName(Sub_mesh_2, 'Sub-mesh_2')
- smesh.SetName(inlet_1, 'inlet')
- try:
- Mesh_1.ExportUNV(r'/home/mika/OpenFOAM/mika-9/run/nozzle_1/geometry/Mesh.unv')
- pass
- except:
- print('ExportUNV() failed. Invalid file name?')
- if salome.sg.hasDesktop():
- salome.sg.updateObjBrowser()
Editor
You can edit this paste and save as new:
File Description
- ffff
- Paste Code
- 02 Dec-2022
- 28.76 Kb
You can Share it:
Latest Code Pastes