# Attention! The scripts for Adding nodes and Elements have been written using the old approach basing on direct usage of SMESH idl interface.
# For the moment smesh package doesn't provide methods to add nodes and elements.
# In the next SALOME version the scripts will be updated to use only the commands from smesh package.
import SMESH
import SMESH_mechanic
mesh = SMESH_mechanic.mesh
aMeshEditor = mesh.GetMeshEditor()
# add node
new_id = aMeshEditor.AddNode(50, 10, 0)
print ""
if new_id == 0: print "KO node addition."
else: print "New Node has been added with ID ", new_id
import SMESH
import SMESH_mechanic
mesh = SMESH_mechanic.mesh
aMeshEditor = mesh.GetMeshEditor()
print ""
# add node
n1 = aMeshEditor.AddNode(50, 10, 0)
if n1 == 0: print "KO node addition."
# add edge
e1 = aMeshEditor.AddEdge([n1, 38])
if e1 == 0: print "KO edge addition."
else: print "New Edge has been added with ID ", e1
import SMESH
import SMESH_mechanic
mesh = SMESH_mechanic.mesh
aMeshEditor = mesh.GetMeshEditor()
print ""
# add node
n1 = aMeshEditor.AddNode(50, 10, 0)
if n1 == 0: print "KO node addition."
# add triangle
t1 = aMeshEditor.AddFace([n1, 38, 39])
if t1 == 0: print "KO triangle addition."
else: print "New Triangle has been added with ID ", t1
import SMESH
import SMESH_mechanic
mesh = SMESH_mechanic.mesh
aMeshEditor = mesh.GetMeshEditor()
print ""
# add node
n1 = aMeshEditor.AddNode(50, 10, 0)
if n1 == 0: print "KO node addition."
n2 = aMeshEditor.AddNode(40, 20, 0)
if n2 == 0: print "KO node addition."
# add quadrangle
q1 = aMeshEditor.AddFace([n2, n1, 38, 39])
if q1 == 0: print "KO quadrangle addition."
else: print "New Quadrangle has been added with ID ", q1
import SMESH
import SMESH_mechanic
mesh = SMESH_mechanic.mesh
aMeshEditor = mesh.GetMeshEditor()
print ""
# add node
n1 = aMeshEditor.AddNode(50, 10, 0)
if n1 == 0: print "KO node addition."
# add tetrahedron
t1 = aMeshEditor.AddVolume([n1, 38, 39, 246])
if t1 == 0: print "KO tetrahedron addition."
else: print "New Tetrahedron has been added with ID ", t1
import SMESH
import SMESH_mechanic
mesh = SMESH_mechanic.mesh
aMeshEditor = mesh.GetMeshEditor()
print ""
# add nodes
nId1 = aMeshEditor.AddNode(50, 10, 0)
nId2 = aMeshEditor.AddNode(47, 12, 0)
nId3 = aMeshEditor.AddNode(50, 10, 10)
nId4 = aMeshEditor.AddNode(47, 12, 10)
if nId1 == 0 or nId2 == 0 or nId3 == 0 or nId4 == 0: print "KO node addition."
# add hexahedron
vId = aMeshEditor.AddVolume([nId2, nId1, 38, 39, nId4, nId3, 245, 246])
if vId == 0: print "KO Hexahedron addition."
else: print "New Hexahedron has been added with ID ", vId
import math
import salome
import smesh
# create an empty mesh structure
gen = smesh.smesh
mesh = gen.CreateEmptyMesh()
MeshEditor = mesh.GetMeshEditor()
# a method to build a polygonal mesh element with <nb_vert> angles:
def MakePolygon (a_mesh, x0, y0, z0, radius, nb_vert):
al = 2.0 * math.pi / nb_vert
node_ids = []
# Create nodes for a polygon
for ii in range(nb_vert):
nid = MeshEditor.AddNode(x0 + radius * math.cos(ii*al),
y0 + radius * math.sin(ii*al),
z0)
node_ids.append(nid)
pass
# Create a polygon
return MeshEditor.AddPolygonalFace(node_ids)
# Create three polygons
f1 = MakePolygon(mesh, 0, 0, 0, 30, 13)
f2 = MakePolygon(mesh, 0, 0, 10, 21, 9)
f3 = MakePolygon(mesh, 0, 0, 20, 13, 6)
salome.sg.updateObjBrowser(1)
import salome
import math
import smesh
# create an empty mesh structure
gen = smesh.smesh
mesh = gen.CreateEmptyMesh()
MeshEditor = mesh.GetMeshEditor()
# Create nodes for 12-hedron with pentagonal faces
al = 2 * math.pi / 5.0
cosal = math.cos(al)
aa = 13
rr = aa / (2.0 * math.sin(al/2.0))
dr = 2.0 * rr * cosal
r1 = rr + dr
dh = rr * math.sqrt(2.0 * (1.0 - cosal * (1.0 + 2.0 * cosal)))
hh = 2.0 * dh - dr * (rr*(cosal - 1) + (rr + dr)*(math.cos(al/2) - 1)) / dh
dd = [] # top
cc = [] # below top
bb = [] # above bottom
aa = [] # bottom
for i in range(5):
cos_bot = math.cos(i*al)
sin_bot = math.sin(i*al)
cos_top = math.cos(i*al + al/2.0)
sin_top = math.sin(i*al + al/2.0)
nd = MeshEditor.AddNode(rr * cos_top, rr * sin_top, hh ) # top
nc = MeshEditor.AddNode(r1 * cos_top, r1 * sin_top, hh - dh) # below top
nb = MeshEditor.AddNode(r1 * cos_bot, r1 * sin_bot, dh) # above bottom
na = MeshEditor.AddNode(rr * cos_bot, rr * sin_bot, 0) # bottom
dd.append(nd) # top
cc.append(nc) # below top
bb.append(nb) # above bottom
aa.append(na) # bottom
pass
# Create a polyhedral volume (12-hedron with pentagonal faces)
MeshEditor.AddPolyhedralVolume([dd[0], dd[1], dd[2], dd[3], dd[4], # top
dd[0], cc[0], bb[1], cc[1], dd[1], # -
dd[1], cc[1], bb[2], cc[2], dd[2], # -
dd[2], cc[2], bb[3], cc[3], dd[3], # - below top
dd[3], cc[3], bb[4], cc[4], dd[4], # -
dd[4], cc[4], bb[0], cc[0], dd[0], # -
aa[4], bb[4], cc[4], bb[0], aa[0], # .
aa[3], bb[3], cc[3], bb[4], aa[4], # .
aa[2], bb[2], cc[2], bb[3], aa[3], # . above bottom
aa[1], bb[1], cc[1], bb[2], aa[2], # .
aa[0], bb[0], cc[0], bb[1], aa[1], # .
aa[0], aa[1], aa[2], aa[3], aa[4]], # bottom
[5,5,5,5,5,5,5,5,5,5,5,5])
salome.sg.updateObjBrowser(1)
# Attention! This script has been written using the old approach basing on direct usage of SMESH idl interface.
# For the moment smesh package doesn't provide methods to remove nodes and elements.
# In the next SALOME version the scripts will be updated to use only the commands from smesh package.
import SMESH_mechanic
mesh = SMESH_mechanic.mesh
aMeshEditor = mesh.GetMeshEditor()
# remove nodes #246 and #255
res = aMeshEditor.RemoveNodes([246, 255])
if res == 1: print "Nodes removing is OK!"
else: print "KO nodes removing."
import SMESH_mechanic
mesh = SMESH_mechanic.mesh
anEditor = mesh.GetMeshEditor()
# remove three elements: #850, #859 and #814
res = anEditor.RemoveElements([850, 859, 814])
if res == 1: print "Elements removing is OK!"
else: print "KO Elements removing."
# Attention! This script has been written using the old approach basing on direct usage of SMESH idl interface.
# For the moment smesh package doesn't provide methods to renumber nodes and elements.
# In the next SALOME version the scripts will be updated to use only the commands from smesh package.
import SMESH_mechanic
mesh = SMESH_mechanic.mesh
anEditor = mesh.GetMeshEditor()
anEditor.RenumberNodes()
# Attention! This script has been written using the old approach basing on direct usage of SMESH idl interface.
# For the moment smesh package doesn't provide methods to move nodes.
# In the next SALOME version the scripts will be updated to use only the commands from smesh package.
import SMESH
import SMESH_mechanic
mesh = SMESH_mechanic.mesh
aMeshEditor = mesh.GetMeshEditor()
# move node #38
aMeshEditor.MoveNode(38, 20., 10., 0.)
# Attention! This script has been written using the old approach basing on direct usage of SMESH idl interface.
# For the moment smesh package doesn't provide methods to produce a diagonal inversion.
# In the next SALOME version the scripts will be updated to use only the commands from smesh package.
import salome
import smesh
# create an empty mesh structure
gen = smesh.smesh
mesh = gen.CreateEmptyMesh()
aMeshEditor = mesh.GetMeshEditor()
# create the following mesh:
# .----.----.----.
# | /| /| /|
# | / | / | / |
# | / | / | / |
# |/ |/ |/ |
# .----.----.----.
bb = [0, 0, 0, 0]
tt = [0, 0, 0, 0]
ff = [0, 0, 0, 0, 0, 0]
bb[0] = aMeshEditor.AddNode( 0., 0., 0.)
bb[1] = aMeshEditor.AddNode(10., 0., 0.)
bb[2] = aMeshEditor.AddNode(20., 0., 0.)
bb[3] = aMeshEditor.AddNode(30., 0., 0.)
tt[0] = aMeshEditor.AddNode( 0., 15., 0.)
tt[1] = aMeshEditor.AddNode(10., 15., 0.)
tt[2] = aMeshEditor.AddNode(20., 15., 0.)
tt[3] = aMeshEditor.AddNode(30., 15., 0.)
ff[0] = aMeshEditor.AddFace([bb[0], bb[1], tt[1]])
ff[1] = aMeshEditor.AddFace([bb[0], tt[1], tt[0]])
ff[2] = aMeshEditor.AddFace([bb[1], bb[2], tt[2]])
ff[3] = aMeshEditor.AddFace([bb[1], tt[2], tt[1]])
ff[4] = aMeshEditor.AddFace([bb[2], bb[3], tt[3]])
ff[5] = aMeshEditor.AddFace([bb[2], tt[3], tt[2]])
# inverse the diagonal bb[1] - tt[2]
print "\nDiagonal inversion ... ",
res = aMeshEditor.InverseDiag(bb[1], tt[2])
if not res: print "failed!"
else: print "done."
salome.sg.updateObjBrowser(1)
# Attention! This script has been written using the old approach basing on direct usage of SMESH idl interface.
# For the moment smesh package doesn't provide methods to unite two triangles.
# In the next SALOME version the scripts will be updated to use only the commands from smesh package.
import salome
import smesh
# create an empty mesh structure
gen = smesh.smesh
mesh = gen.CreateEmptyMesh()
aMeshEditor = mesh.GetMeshEditor()
# create the following mesh:
# .----.----.----.
# | /| /| /|
# | / | / | / |
# | / | / | / |
# |/ |/ |/ |
# .----.----.----.
bb = [0, 0, 0, 0]
tt = [0, 0, 0, 0]
ff = [0, 0, 0, 0, 0, 0]
bb[0] = aMeshEditor.AddNode( 0., 0., 0.)
bb[1] = aMeshEditor.AddNode(10., 0., 0.)
bb[2] = aMeshEditor.AddNode(20., 0., 0.)
bb[3] = aMeshEditor.AddNode(30., 0., 0.)
tt[0] = aMeshEditor.AddNode( 0., 15., 0.)
tt[1] = aMeshEditor.AddNode(10., 15., 0.)
tt[2] = aMeshEditor.AddNode(20., 15., 0.)
tt[3] = aMeshEditor.AddNode(30., 15., 0.)
ff[0] = aMeshEditor.AddFace([bb[0], bb[1], tt[1]])
ff[1] = aMeshEditor.AddFace([bb[0], tt[1], tt[0]])
ff[2] = aMeshEditor.AddFace([bb[1], bb[2], tt[2]])
ff[3] = aMeshEditor.AddFace([bb[1], tt[2], tt[1]])
ff[4] = aMeshEditor.AddFace([bb[2], bb[3], tt[3]])
ff[5] = aMeshEditor.AddFace([bb[2], tt[3], tt[2]])
# delete the diagonal bb[1] - tt[2]
print "\nUnite two triangles ... ",
res = aMeshEditor.DeleteDiag(bb[1], tt[2])
if not res: print "failed!"
else: print "done."
salome.sg.updateObjBrowser(1)
# Attention! This script has been written using the old approach basing on direct usage of SMESH idl interface.
# For the moment smesh package doesn't provide methods to unite a set of triangles.
# In the next SALOME version the scripts will be updated to use only the commands from smesh package.
import salome
import smesh
# create an empty mesh structure
gen = smesh.smesh
mesh = gen.CreateEmptyMesh()
aMeshEditor = mesh.GetMeshEditor()
# create the following mesh:
# .----.----.----.
# | /| /| /|
# | / | / | / |
# | / | / | / |
# |/ |/ |/ |
# .----.----.----.
bb = [0, 0, 0, 0]
tt = [0, 0, 0, 0]
ff = [0, 0, 0, 0, 0, 0]
bb[0] = aMeshEditor.AddNode( 0., 0., 0.)
bb[1] = aMeshEditor.AddNode(10., 0., 0.)
bb[2] = aMeshEditor.AddNode(20., 0., 0.)
bb[3] = aMeshEditor.AddNode(30., 0., 0.)
tt[0] = aMeshEditor.AddNode( 0., 15., 0.)
tt[1] = aMeshEditor.AddNode(10., 15., 0.)
tt[2] = aMeshEditor.AddNode(20., 15., 0.)
tt[3] = aMeshEditor.AddNode(30., 15., 0.)
ff[0] = aMeshEditor.AddFace([bb[0], bb[1], tt[1]])
ff[1] = aMeshEditor.AddFace([bb[0], tt[1], tt[0]])
ff[2] = aMeshEditor.AddFace([bb[1], bb[2], tt[2]])
ff[3] = aMeshEditor.AddFace([bb[1], tt[2], tt[1]])
ff[4] = aMeshEditor.AddFace([bb[2], bb[3], tt[3]])
ff[5] = aMeshEditor.AddFace([bb[2], tt[3], tt[2]])
# unite a set of triangles
aFilterMgr = smesh.smesh.CreateFilterManager()
aFunctor = aFilterMgr.CreateMinimumAngle()
print "\nUnite a set of triangles ... ",
res = aMeshEditor.TriToQuad([ff[2], ff[3], ff[4], ff[5]], aFunctor, 60.)
if not res: print "failed!"
else: print "done."
salome.sg.updateObjBrowser(1)
# Attention! This script has been written using the old approach basing on direct usage of SMESH idl interface.
# For the moment smesh package doesn't provide methods to change orientation.
# In the next SALOME version the scripts will be updated to use only the commands from smesh package.
import salome
import smesh
# SMESH module
gen = smesh.smesh
mesh = gen.CreateEmptyMesh()
MeshEditor = mesh.GetMeshEditor()
# build five quadrangles:
dx = 10
dy = 20
n1 = MeshEditor.AddNode(0.0 * dx, 0, 0)
n2 = MeshEditor.AddNode(1.0 * dx, 0, 0)
n3 = MeshEditor.AddNode(2.0 * dx, 0, 0)
n4 = MeshEditor.AddNode(3.0 * dx, 0, 0)
n5 = MeshEditor.AddNode(4.0 * dx, 0, 0)
n6 = MeshEditor.AddNode(5.0 * dx, 0, 0)
n7 = MeshEditor.AddNode(0.0 * dx, dy, 0)
n8 = MeshEditor.AddNode(1.0 * dx, dy, 0)
n9 = MeshEditor.AddNode(2.0 * dx, dy, 0)
n10 = MeshEditor.AddNode(3.0 * dx, dy, 0)
n11 = MeshEditor.AddNode(4.0 * dx, dy, 0)
n12 = MeshEditor.AddNode(5.0 * dx, dy, 0)
f1 = MeshEditor.AddFace([n1, n2, n8 , n7 ])
f2 = MeshEditor.AddFace([n2, n3, n9 , n8 ])
f3 = MeshEditor.AddFace([n3, n4, n10, n9 ])
f4 = MeshEditor.AddFace([n4, n5, n11, n10])
f5 = MeshEditor.AddFace([n5, n6, n12, n11])
# Change the orientation of the second and the fourth faces.
MeshEditor.Reorient([2, 4])
salome.sg.updateObjBrowser(1)
# Attention! This script has been written using the old approach basing on direct usage of SMESH idl interface.
# For the moment smesh package doesn't provide methods to cut quadrangles.
# In the next SALOME version the scripts will be updated to use only the commands from smesh package.
import SMESH
import SMESH_mechanic
smesh = SMESH_mechanic.smesh
mesh = SMESH_mechanic.mesh
# cut two quadrangles: 405 and 406
aFilterMgr = smesh.CreateFilterManager()
aFunctor = aFilterMgr.CreateMinimumAngle()
aMeshEditor = mesh.GetMeshEditor()
aMeshEditor.QuadToTri([405, 406], aFunctor)
# Attention! This script has been written using the old approach basing on direct usage of SMESH idl interface.
# For the moment smesh package doesn't provide methods to perform smoothing.
# In the next SALOME version the scripts will be updated to use only the commands from smesh package.
import salome
import geompy
import SMESH
import SMESH_mechanic
mesh = SMESH_mechanic.mesh
aMeshEditor = mesh.GetMeshEditor()
# select the top face
faces = geompy.SubShapeAllSorted(SMESH_mechanic.shape_mesh, geompy.ShapeType["FACE"])
face = faces[3]
geompy.addToStudyInFather(SMESH_mechanic.shape_mesh, face, "face planar with hole")
# create a group of faces to be smoothed
GroupSmooth = mesh.CreateGroupFromGEOM(SMESH.FACE, "Group of faces (smooth)", face)
# perform smoothing
# boolean SmoothObject(Object, IDsOfFixedNodes, MaxNbOfIterations, MaxAspectRatio, Method)
res = aMeshEditor.SmoothObject(GroupSmooth, [], 20, 2., SMESH.SMESH_MeshEditor.CENTROIDAL_SMOOTH)
print "\nSmoothing ... ",
if not res: print "failed!"
else: print "done."
salome.sg.updateObjBrowser(1)
# Attention! This script has been written using the old approach basing on direct usage of SMESH idl interface.
# For the moment smesh package doesn't provide methods to perform extrusion.
# In the next SALOME version the scripts will be updated to use only the commands from smesh package.
import salome
import geompy
import SMESH
import SMESH_mechanic
mesh = SMESH_mechanic.mesh
aMeshEditor = mesh.GetMeshEditor()
# select the top face
faces = geompy.SubShapeAllSorted(SMESH_mechanic.shape_mesh, geompy.ShapeType["FACE"])
face = faces[7]
geompy.addToStudyInFather(SMESH_mechanic.shape_mesh, face, "face circular top")
# create a vector for extrusion
point = SMESH.PointStruct(0., 0., 5.)
vector = SMESH.DirStruct(point)
# create a group to be extruded
GroupTri = mesh.CreateGroupFromGEOM(SMESH.FACE, "Group of faces (extrusion)", face)
# perform extrusion of the group
aMeshEditor.ExtrusionSweepObject(GroupTri, vector, 5)
salome.sg.updateObjBrowser(1)
# Attention! This script has been written using the old approach basing on direct usage of SMESH idl interface.
# For the moment smesh package doesn't provide methods to perform extrusion along a path.
# In the next SALOME version the scripts will be updated to use only the commands from smesh package.
import math
import salome
# Geometry
import geompy
# 1. Create points
points = [[0, 0], [50, 30], [50, 110], [0, 150], [-80, 150], [-130, 70], [-130, -20]]
iv = 1
vertices = []
for point in points:
vert = geompy.MakeVertex(point[0], point[1], 0)
geompy.addToStudy(vert, "Vertex_" + `iv`)
vertices.append(vert)
iv += 1
pass
# 2. Create edges and wires
Edge_straight = geompy.MakeEdge(vertices[0], vertices[4])
Edge_bezierrr = geompy.MakeBezier(vertices)
Wire_polyline = geompy.MakePolyline(vertices)
Edge_Circle = geompy.MakeCircleThreePnt(vertices[0], vertices[1], vertices[2])
geompy.addToStudy(Edge_straight, "Edge_straight")
geompy.addToStudy(Edge_bezierrr, "Edge_bezierrr")
geompy.addToStudy(Wire_polyline, "Wire_polyline")
geompy.addToStudy(Edge_Circle , "Edge_Circle")
# 3. Explode wire on edges, as they will be used for mesh extrusion
Wire_polyline_edges = geompy.SubShapeAll(Wire_polyline, geompy.ShapeType["EDGE"])
for ii in range(len(Wire_polyline_edges)):
geompy.addToStudyInFather(Wire_polyline, Wire_polyline_edges[ii], "Edge_" + `ii + 1`)
pass
# Mesh
import smesh
import SMESH
gen = smesh.smesh
smeshgui = salome.ImportComponentGUI("SMESH")
smeshgui.Init(salome.myStudyId)
# 1D algorithm and three 1D hypotheses
Wire_discretisation = gen.CreateHypothesis('Regular_1D', 'libStdMeshersEngine.so')
Nb_Segments_3 = gen.CreateHypothesis('NumberOfSegments', 'libStdMeshersEngine.so')
Nb_Segments_7 = gen.CreateHypothesis('NumberOfSegments', 'libStdMeshersEngine.so')
Nb_Segments_8 = gen.CreateHypothesis('NumberOfSegments', 'libStdMeshersEngine.so')
Nb_Segments_3.SetNumberOfSegments(3)
Nb_Segments_7.SetNumberOfSegments(7)
Nb_Segments_8.SetNumberOfSegments(8)
# Mesh given shape with given 1d hypothesis
def Mesh1D(shape1d, hyp1d, name):
mesh1d_tool = smesh.Mesh(shape1d)
mesh1d = mesh1d_tool.GetMesh()
status = mesh1d.AddHypothesis(shape1d, hyp1d)
status = mesh1d.AddHypothesis(shape1d, Wire_discretisation)
isDone = mesh1d_tool.Compute()
if not isDone: print 'Mesh ', name, ': computation failed'
return mesh1d
# Create a mesh with six nodes, seven edges and two quadrangle faces
def MakeQuadMesh2(mesh_name):
quad_1 = gen.CreateEmptyMesh()
smeshgui.SetName(salome.ObjectToID(quad_1), mesh_name)
editor_1 = quad_1.GetMeshEditor()
# six nodes
n1 = editor_1.AddNode(0, 20, 10)
n2 = editor_1.AddNode(0, 40, 10)
n3 = editor_1.AddNode(0, 40, 30)
n4 = editor_1.AddNode(0, 20, 30)
n5 = editor_1.AddNode(0, 0, 30)
n6 = editor_1.AddNode(0, 0, 10)
# seven edges
editor_1.AddEdge([n1, n2]) # 1
editor_1.AddEdge([n2, n3]) # 2
editor_1.AddEdge([n3, n4]) # 3
editor_1.AddEdge([n4, n1]) # 4
editor_1.AddEdge([n4, n5]) # 5
editor_1.AddEdge([n5, n6]) # 6
editor_1.AddEdge([n6, n1]) # 7
# two quadrangle faces
editor_1.AddFace([n1, n2, n3, n4]) # 8
editor_1.AddFace([n1, n4, n5, n6]) # 9
return [quad_1, editor_1, [1,2,3,4,5,6,7], [8,9]]
# Path meshes
Edge_straight_mesh = Mesh1D(Edge_straight, Nb_Segments_7, "Edge_straight")
Edge_bezierrr_mesh = Mesh1D(Edge_bezierrr, Nb_Segments_7, "Edge_bezierrr")
Wire_polyline_mesh = Mesh1D(Wire_polyline, Nb_Segments_3, "Wire_polyline")
Edge_Circle_mesh = Mesh1D(Edge_Circle , Nb_Segments_8, "Edge_Circle")
# Initial meshes (to be extruded)
[quad_1, editor_1, ee_1, ff_1] = MakeQuadMesh2("quad_1")
[quad_2, editor_2, ee_2, ff_2] = MakeQuadMesh2("quad_2")
[quad_3, editor_3, ee_3, ff_3] = MakeQuadMesh2("quad_3")
[quad_4, editor_4, ee_4, ff_4] = MakeQuadMesh2("quad_4")
[quad_5, editor_5, ee_5, ff_5] = MakeQuadMesh2("quad_5")
[quad_6, editor_6, ee_6, ff_6] = MakeQuadMesh2("quad_6")
[quad_7, editor_7, ee_7, ff_7] = MakeQuadMesh2("quad_7")
# ExtrusionAlongPath
# IDsOfElements, PathMesh, PathShape, NodeStart,
# HasAngles, Angles, HasRefPoint, RefPoint
refPoint = SMESH.PointStruct(0, 0, 0)
a10 = 10.0*math.pi/180.0
a45 = 45.0*math.pi/180.0
# 1. Extrusion of two mesh edges along a straight path
error = editor_1.ExtrusionAlongPath([1,2], Edge_straight_mesh, Edge_straight, 1,
0, [], 0, refPoint)
# 2. Extrusion of one mesh edge along a curved path
error = editor_2.ExtrusionAlongPath([2], Edge_bezierrr_mesh, Edge_bezierrr, 1,
0, [], 0, refPoint)
# 3. Extrusion of one mesh edge along a curved path with usage of angles
error = editor_3.ExtrusionAlongPath([2], Edge_bezierrr_mesh, Edge_bezierrr, 1,
1, [a45, a45, a45, 0, -a45, -a45, -a45], 0, refPoint)
# 4. Extrusion of one mesh edge along the path, which is a part of a meshed wire
error = editor_4.ExtrusionAlongPath([4], Wire_polyline_mesh, Wire_polyline_edges[0], 1,
1, [a10, a10, a10], 0, refPoint)
# 5. Extrusion of two mesh faces along the path, which is a part of a meshed wire
error = editor_5.ExtrusionAlongPath(ff_5 , Wire_polyline_mesh, Wire_polyline_edges[2], 4,
0, [], 0, refPoint)
# 6. Extrusion of two mesh faces along a closed path
error = editor_6.ExtrusionAlongPath(ff_6 , Edge_Circle_mesh, Edge_Circle, 1,
0, [], 0, refPoint)
# 7. Extrusion of two mesh faces along a closed path with usage of angles
error = editor_7.ExtrusionAlongPath(ff_7, Edge_Circle_mesh, Edge_Circle, 1,
1, [a45, -a45, a45, -a45, a45, -a45, a45, -a45], 0, refPoint
salome.sg.updateObjBrowser(1)
# Attention! This script has been written using the old approach basing on direct usage of SMESH idl interface.
# For the moment smesh package doesn't provide methods to perform revolution.
# In the next SALOME version the scripts will be updated to use only the commands from smesh package.
import math
import SMESH
import SMESH_mechanic
mesh = SMESH_mechanic.mesh
aMeshEditor = mesh.GetMeshEditor()
# create a group of faces to be revolved
FacesRotate = [492, 493, 502, 503]
GroupRotate = mesh.CreateGroup(SMESH.FACE,"Group of faces (rotate)")
GroupRotate.Add(FacesRotate)
# define revolution angle and axis
angle45 = 45 * math.pi / 180
axisXYZ = SMESH.AxisStruct(-38.3128, -73.3658, -23.321, -13.3402, -13.3265, 6.66632)
# perform revolution of an object
aMeshEditor.RotationSweepObject(GroupRotate, axisXYZ, angle45, 4, 1e-5)
# Attention! This script has been written using the old approach basing on direct usage of SMESH idl interface.
# For the moment smesh package doesn't provide methods to perform pattern mapping.
# In the next SALOME version the scripts will be updated to use only the commands from smesh package.
import geompy
import smesh
# define the geometry
Box_1 = geompy.MakeBoxDXDYDZ(200., 200., 200.)
geompy.addToStudy(Box_1, "Box_1")
faces = geompy.SubShapeAll(Box_1, geompy.ShapeType["FACE"])
Face_1 = faces[0]
Face_2 = faces[1]
geompy.addToStudyInFather(Box_1, Face_1, "Face_1")
geompy.addToStudyInFather(Box_1, Face_2, "Face_2")
# build a quadrangle mesh 3x3 on Face_1
Mesh_1 = smesh.Mesh(Face_1)
algo1D = Mesh_1.Segment()
algo1D.NumberOfSegments(3)
Mesh_1.Quadrangle()
isDone = Mesh_1.Compute()
if not isDone: print 'Mesh Mesh_1 : computation failed'
# build a triangle mesh on Face_2
Mesh_2 = smesh.Mesh(Face_2)
algo1D = Mesh_2.Segment()
algo1D.NumberOfSegments(1)
algo2D = Mesh_2.Triangle()
algo2D.MaxElementArea(240)
isDone = Mesh_2.Compute()
if not isDone: print 'Mesh Mesh_2 : computation failed'
# create a pattern
pattern = smesh.smesh.GetPattern()
isDone = pattern.LoadFromFace(Mesh_2.GetMesh(), Face_2, 0)
if (isDone != 1): print 'LoadFromFace :', pattern.GetErrorCode()
# apply the pattern to a face of the first mesh
pattern.ApplyToMeshFaces(Mesh_1.GetMesh(), [17], 0, 0)
isDone = pattern.MakeMesh(Mesh_1.GetMesh(), 0, 0)
if (isDone != 1): print 'MakeMesh :', pattern.GetErrorCode()