mirror of
https://git.salome-platform.org/gitpub/modules/smesh.git
synced 2024-12-26 01:10:35 +05:00
Check result of Compute() in test and examples
This commit is contained in:
parent
d1392dc034
commit
1ade5b1884
@ -48,15 +48,19 @@ algo_2D.SetMinSize( 7. )
|
|||||||
|
|
||||||
# 2) Compute 2D mesh
|
# 2) Compute 2D mesh
|
||||||
isDone = Mesh_1.Compute()
|
isDone = Mesh_1.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# 3) Split quadrangles into triangles
|
# 3) Split quadrangles into triangles
|
||||||
isDone = Mesh_1.SplitQuadObject( Mesh_1, Diag13=True )
|
Mesh_1.SplitQuadObject( Mesh_1, Diag13=True )
|
||||||
|
|
||||||
# 4) Define a 3D mesher
|
# 4) Define a 3D mesher
|
||||||
Mesh_1.Tetrahedron()
|
Mesh_1.Tetrahedron()
|
||||||
|
|
||||||
# 5) Compute 3D mesh
|
# 5) Compute 3D mesh
|
||||||
Mesh_1.Compute()
|
isDone = Mesh_1.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
if salome.sg.hasDesktop():
|
if salome.sg.hasDesktop():
|
||||||
salome.sg.updateObjBrowser()
|
salome.sg.updateObjBrowser()
|
||||||
|
@ -25,7 +25,7 @@ coords = list(range(-100,100,10))
|
|||||||
cartHyp = cartAlgo.SetGrid( coords,coords,coords, 1000000)
|
cartHyp = cartAlgo.SetGrid( coords,coords,coords, 1000000)
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
print("nb hexahedra",mesh.NbHexas())
|
print("nb hexahedra",mesh.NbHexas())
|
||||||
print("nb tetrahedra",mesh.NbTetras())
|
print("nb tetrahedra",mesh.NbTetras())
|
||||||
print("nb polyhedra",mesh.NbPolyhedrons())
|
print("nb polyhedra",mesh.NbPolyhedrons())
|
||||||
@ -34,7 +34,7 @@ print()
|
|||||||
# define the grid by setting constant spacing
|
# define the grid by setting constant spacing
|
||||||
cartHyp = cartAlgo.SetGrid( "10","10","10", 1000000)
|
cartHyp = cartAlgo.SetGrid( "10","10","10", 1000000)
|
||||||
|
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
print("nb hexahedra",mesh.NbHexas())
|
print("nb hexahedra",mesh.NbHexas())
|
||||||
print("nb tetrahedra",mesh.NbTetras())
|
print("nb tetrahedra",mesh.NbTetras())
|
||||||
print("nb polyhedra",mesh.NbPolyhedrons())
|
print("nb polyhedra",mesh.NbPolyhedrons())
|
||||||
@ -44,7 +44,7 @@ print()
|
|||||||
# activate creation of faces
|
# activate creation of faces
|
||||||
cartHyp.SetToCreateFaces( True )
|
cartHyp.SetToCreateFaces( True )
|
||||||
|
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
print("nb hexahedra",mesh.NbHexas())
|
print("nb hexahedra",mesh.NbHexas())
|
||||||
print("nb tetrahedra",mesh.NbTetras())
|
print("nb tetrahedra",mesh.NbTetras())
|
||||||
print("nb polyhedra",mesh.NbPolyhedrons())
|
print("nb polyhedra",mesh.NbPolyhedrons())
|
||||||
@ -53,7 +53,7 @@ print()
|
|||||||
|
|
||||||
# enable consideration of shared faces
|
# enable consideration of shared faces
|
||||||
cartHyp.SetToConsiderInternalFaces( True )
|
cartHyp.SetToConsiderInternalFaces( True )
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
print("nb hexahedra",mesh.NbHexas())
|
print("nb hexahedra",mesh.NbHexas())
|
||||||
print("nb tetrahedra",mesh.NbTetras())
|
print("nb tetrahedra",mesh.NbTetras())
|
||||||
print("nb polyhedra",mesh.NbPolyhedrons())
|
print("nb polyhedra",mesh.NbPolyhedrons())
|
||||||
@ -64,7 +64,7 @@ print()
|
|||||||
spaceFuns = ["5","10+10*t"]
|
spaceFuns = ["5","10+10*t"]
|
||||||
cartAlgo.SetGrid( [spaceFuns, [0.5]], [spaceFuns, [0.5]], [spaceFuns, [0.25]], 10 )
|
cartAlgo.SetGrid( [spaceFuns, [0.5]], [spaceFuns, [0.5]], [spaceFuns, [0.25]], 10 )
|
||||||
|
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
print("nb hexahedra",mesh.NbHexas())
|
print("nb hexahedra",mesh.NbHexas())
|
||||||
print("nb tetrahedra",mesh.NbTetras())
|
print("nb tetrahedra",mesh.NbTetras())
|
||||||
print("nb polyhedra",mesh.NbPolyhedrons())
|
print("nb polyhedra",mesh.NbPolyhedrons())
|
||||||
@ -85,24 +85,24 @@ spc = "0.1" # spacing
|
|||||||
mesh = smesh_builder.Mesh( box, "custom axes")
|
mesh = smesh_builder.Mesh( box, "custom axes")
|
||||||
algo = mesh.BodyFitted()
|
algo = mesh.BodyFitted()
|
||||||
algo.SetGrid( spc, spc, spc, 10000 )
|
algo.SetGrid( spc, spc, spc, 10000 )
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
print("Default axes")
|
print("Default axes")
|
||||||
print(" nb hex:",mesh.NbHexas())
|
print(" nb hex:",mesh.NbHexas())
|
||||||
|
|
||||||
# set axes using edges of the box
|
# set axes using edges of the box
|
||||||
algo.SetAxesDirs( xDir, [-0.1,1,0], zDir )
|
algo.SetAxesDirs( xDir, [-0.1,1,0], zDir )
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
print("Manual axes")
|
print("Manual axes")
|
||||||
print(" nb hex:",mesh.NbHexas())
|
print(" nb hex:",mesh.NbHexas())
|
||||||
|
|
||||||
# set optimal orthogonal axes
|
# set optimal orthogonal axes
|
||||||
algo.SetOptimalAxesDirs( isOrthogonal=True )
|
algo.SetOptimalAxesDirs( isOrthogonal=True )
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
print("Optimal orthogonal axes")
|
print("Optimal orthogonal axes")
|
||||||
print(" nb hex:",mesh.NbHexas())
|
print(" nb hex:",mesh.NbHexas())
|
||||||
|
|
||||||
# set optimal non-orthogonal axes
|
# set optimal non-orthogonal axes
|
||||||
algo.SetOptimalAxesDirs( isOrthogonal=False )
|
algo.SetOptimalAxesDirs( isOrthogonal=False )
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
print("Optimal non-orthogonal axes")
|
print("Optimal non-orthogonal axes")
|
||||||
print(" nb hex:",mesh.NbHexas())
|
print(" nb hex:",mesh.NbHexas())
|
||||||
|
@ -45,15 +45,14 @@ Mesh_1 = smesh.Mesh(Sphere_1,'Mesh_1')
|
|||||||
status = Mesh_1.AddHypothesis( Sphere_1, NETGEN_3D_Parameters_1 )
|
status = Mesh_1.AddHypothesis( Sphere_1, NETGEN_3D_Parameters_1 )
|
||||||
NETGEN_1D_2D_3D = Mesh_1.Tetrahedron(algo=smeshBuilder.NETGEN_1D2D3D)
|
NETGEN_1D_2D_3D = Mesh_1.Tetrahedron(algo=smeshBuilder.NETGEN_1D2D3D)
|
||||||
isDone = Mesh_1.Compute()
|
isDone = Mesh_1.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Creating Dual mesh
|
# Creating Dual mesh
|
||||||
dual_Mesh_1 = smesh.CreateDualMesh( Mesh_1, 'dual_Mesh_1', True)
|
dual_Mesh_1 = smesh.CreateDualMesh( Mesh_1, 'dual_Mesh_1', True)
|
||||||
|
|
||||||
|
|
||||||
assert(dual_Mesh_1.NbPolyhedrons() > 0)
|
assert(dual_Mesh_1.NbPolyhedrons() > 0)
|
||||||
assert(dual_Mesh_1.NbTetras() == 0)
|
assert(dual_Mesh_1.NbTetras() == 0)
|
||||||
|
|
||||||
|
|
||||||
if salome.sg.hasDesktop():
|
if salome.sg.hasDesktop():
|
||||||
salome.sg.updateObjBrowser()
|
salome.sg.updateObjBrowser()
|
||||||
|
@ -27,8 +27,7 @@ algo3D.MaxElementVolume(900.)
|
|||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
ret = tetra.Compute()
|
ret = tetra.Compute()
|
||||||
if ret == 0:
|
if not ret:
|
||||||
print("problem when computing the mesh")
|
raise Exception("problem when computing the mesh")
|
||||||
else:
|
|
||||||
print("mesh computed")
|
print("mesh computed")
|
||||||
pass
|
|
||||||
|
@ -52,11 +52,11 @@ mesh.Hexahedron()
|
|||||||
surfaces = geom_builder.SubShapeAll(box, geom_builder.ShapeType["FACE"])
|
surfaces = geom_builder.SubShapeAll(box, geom_builder.ShapeType["FACE"])
|
||||||
|
|
||||||
# method 1: no sub-mesh is created
|
# method 1: no sub-mesh is created
|
||||||
mesh.Compute( surfaces[0] )
|
if not mesh.Compute( surfaces[0] ): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# method 2: a sub-mesh is created
|
# method 2: a sub-mesh is created
|
||||||
submesh = mesh.GetSubMesh( surfaces[2], "submesh 2" )
|
submesh = mesh.GetSubMesh( surfaces[2], "submesh 2" )
|
||||||
submesh.Compute()
|
if not submesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# compute the whole mesh
|
# compute the whole mesh
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -35,14 +35,17 @@ mesh.Triangle(geom=Face_3)
|
|||||||
# get existing sub-mesh priority order: F1 -> F2 -> F3
|
# get existing sub-mesh priority order: F1 -> F2 -> F3
|
||||||
[[SubMesh_F1, SubMesh_F3, SubMesh_F2]] = mesh.GetMeshOrder()
|
[[SubMesh_F1, SubMesh_F3, SubMesh_F2]] = mesh.GetMeshOrder()
|
||||||
isDone = mesh.Compute()
|
isDone = mesh.Compute()
|
||||||
|
if not isDone: raise Exception("Error when computing Mesh")
|
||||||
print("Nb elements at initial order of sub-meshes:", mesh.NbElements())
|
print("Nb elements at initial order of sub-meshes:", mesh.NbElements())
|
||||||
|
|
||||||
# set new sub-mesh order: F2 -> F1 -> F3
|
# set new sub-mesh order: F2 -> F1 -> F3
|
||||||
isDone = mesh.SetMeshOrder([[SubMesh_F2, SubMesh_F1, SubMesh_F3]])
|
isDone = mesh.SetMeshOrder([[SubMesh_F2, SubMesh_F1, SubMesh_F3]])
|
||||||
isDone = mesh.Compute()
|
isDone = mesh.Compute()
|
||||||
|
if not isDone: raise Exception("Error when computing Mesh")
|
||||||
print("Nb elements at new order of sub-meshes:", mesh.NbElements())
|
print("Nb elements at new order of sub-meshes:", mesh.NbElements())
|
||||||
|
|
||||||
# compute with other sub-mesh order: F3 -> F2 -> F1
|
# compute with other sub-mesh order: F3 -> F2 -> F1
|
||||||
isDone = mesh.SetMeshOrder([[SubMesh_F3, SubMesh_F2, SubMesh_F1]])
|
isDone = mesh.SetMeshOrder([[SubMesh_F3, SubMesh_F2, SubMesh_F1]])
|
||||||
isDone = mesh.Compute()
|
isDone = mesh.Compute()
|
||||||
|
if not isDone: raise Exception("Error when computing Mesh")
|
||||||
print("Nb elements at another order of sub-meshes:", mesh.NbElements())
|
print("Nb elements at another order of sub-meshes:", mesh.NbElements())
|
||||||
|
@ -41,19 +41,19 @@ hyp3 = algo_local.Arithmetic1D(1, 6)
|
|||||||
hyp4 = algo_local.Propagation()
|
hyp4 = algo_local.Propagation()
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
tria.Compute()
|
if not tria.Compute(): raise Exception("Error when computing Mesh")
|
||||||
PrintMeshInfo(tria)
|
PrintMeshInfo(tria)
|
||||||
|
|
||||||
# remove a local hypothesis
|
# remove a local hypothesis
|
||||||
tria.RemoveHypothesis(hyp4, edge)
|
tria.RemoveHypothesis(hyp4, edge)
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
tria.Compute()
|
if not tria.Compute(): raise Exception("Error when computing Mesh")
|
||||||
PrintMeshInfo(tria)
|
PrintMeshInfo(tria)
|
||||||
|
|
||||||
# change the value of the 2D hypothesis
|
# change the value of the 2D hypothesis
|
||||||
hyp2.SetMaxElementArea(2.)
|
hyp2.SetMaxElementArea(2.)
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
tria.Compute()
|
if not tria.Compute(): raise Exception("Error when computing Mesh")
|
||||||
PrintMeshInfo(tria)
|
PrintMeshInfo(tria)
|
||||||
|
@ -24,7 +24,7 @@ tetra.Triangle()
|
|||||||
tetra.Tetrahedron()
|
tetra.Tetrahedron()
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
tetra.Compute()
|
if not tetra.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# export the mesh in a MED file
|
# export the mesh in a MED file
|
||||||
medFile = tempfile.NamedTemporaryFile(suffix=".med").name
|
medFile = tempfile.NamedTemporaryFile(suffix=".med").name
|
||||||
|
@ -113,7 +113,7 @@ algo.Propagation()
|
|||||||
hexa.Quadrangle()
|
hexa.Quadrangle()
|
||||||
hexa.Hexahedron()
|
hexa.Hexahedron()
|
||||||
|
|
||||||
hexa.Compute()
|
if not hexa.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
hexa.Group(group_a)
|
hexa.Group(group_a)
|
||||||
hexa.Group(group_b)
|
hexa.Group(group_b)
|
||||||
|
@ -46,7 +46,7 @@ algo1D_1=Mesh_inf.Segment()
|
|||||||
algo1D_1.NumberOfSegments(10)
|
algo1D_1.NumberOfSegments(10)
|
||||||
algo2D_1=Mesh_inf.Quadrangle()
|
algo2D_1=Mesh_inf.Quadrangle()
|
||||||
algo3D_1=Mesh_inf.Hexahedron()
|
algo3D_1=Mesh_inf.Hexahedron()
|
||||||
Mesh_inf.Compute()
|
if not Mesh_inf.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# create a group on the top face
|
# create a group on the top face
|
||||||
Gsup1=Mesh_inf.Group(Fsup1, "Sup")
|
Gsup1=Mesh_inf.Group(Fsup1, "Sup")
|
||||||
@ -59,7 +59,7 @@ algo1D_2=Mesh_sup.Segment()
|
|||||||
algo1D_2.NumberOfSegments(5)
|
algo1D_2.NumberOfSegments(5)
|
||||||
algo2D_2=Mesh_sup.Quadrangle()
|
algo2D_2=Mesh_sup.Quadrangle()
|
||||||
algo3D_2=Mesh_sup.Hexahedron()
|
algo3D_2=Mesh_sup.Hexahedron()
|
||||||
Mesh_sup.Compute()
|
if not Mesh_sup.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# create a group on the top face
|
# create a group on the top face
|
||||||
Gsup2=Mesh_sup.Group(Fsup2, "Sup")
|
Gsup2=Mesh_sup.Group(Fsup2, "Sup")
|
||||||
|
@ -20,7 +20,7 @@ localAlgo = mesh.Triangle(face)
|
|||||||
mesh.Segment().NumberOfSegments( 3 )
|
mesh.Segment().NumberOfSegments( 3 )
|
||||||
mesh.Quadrangle()
|
mesh.Quadrangle()
|
||||||
mesh.Prism()
|
mesh.Prism()
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# objects to copy
|
# objects to copy
|
||||||
fGroup = mesh.GroupOnGeom( face, "2D on face")
|
fGroup = mesh.GroupOnGeom( face, "2D on face")
|
||||||
@ -56,4 +56,4 @@ smallBox = geom_builder.MakeScaleAlongAxes( box, None, 1, 0.5, 0.5 )
|
|||||||
cutBox = geom_builder.MakeCut( box, smallBox, theName="box - smallBox" )
|
cutBox = geom_builder.MakeCut( box, smallBox, theName="box - smallBox" )
|
||||||
|
|
||||||
ok, newMesh, groups, submehses, hyps, invIDs = smesh_builder.CopyMeshWithGeom( mesh, cutBox, "cutBox" )
|
ok, newMesh, groups, submehses, hyps, invIDs = smesh_builder.CopyMeshWithGeom( mesh, cutBox, "cutBox" )
|
||||||
newMesh.Compute()
|
if not newMesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -113,14 +113,18 @@ def run_test(nbox=2, boxsize=100):
|
|||||||
print("Starting sequential compute")
|
print("Starting sequential compute")
|
||||||
start = time.monotonic()
|
start = time.monotonic()
|
||||||
is_done = seq_mesh.Compute()
|
is_done = seq_mesh.Compute()
|
||||||
assert is_done
|
if not is_done:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
stop = time.monotonic()
|
stop = time.monotonic()
|
||||||
time_seq = stop-start
|
time_seq = stop-start
|
||||||
|
|
||||||
print("Starting parallel compute")
|
print("Starting parallel compute")
|
||||||
start = time.monotonic()
|
start = time.monotonic()
|
||||||
is_done = par_mesh.Compute()
|
is_done = par_mesh.Compute()
|
||||||
assert is_done
|
if not is_done:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
stop = time.monotonic()
|
stop = time.monotonic()
|
||||||
time_par = stop-start
|
time_par = stop-start
|
||||||
|
|
||||||
|
@ -31,4 +31,4 @@ deflection = 0.05
|
|||||||
mesh = smesh_builder.Mesh( shape )
|
mesh = smesh_builder.Mesh( shape )
|
||||||
mesh.Segment().Adaptive( minSize, maxSize, deflection )
|
mesh.Segment().Adaptive( minSize, maxSize, deflection )
|
||||||
mesh.Triangle().MaxElementArea( 300 )
|
mesh.Triangle().MaxElementArea( 300 )
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -40,4 +40,4 @@ hexa.Quadrangle()
|
|||||||
hexa.Hexahedron()
|
hexa.Hexahedron()
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
hexa.Compute()
|
if not hexa.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -43,4 +43,4 @@ algo_local = hexa.Segment(e_arc)
|
|||||||
algo_local.Deflection1D(1.0)
|
algo_local.Deflection1D(1.0)
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
hexa.Compute()
|
if not hexa.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -40,4 +40,4 @@ algo_local.StartEndLength(1, 6)
|
|||||||
algo_local.Propagation()
|
algo_local.Propagation()
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
hexa.Compute()
|
if not hexa.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -40,4 +40,4 @@ algo_local.LocalLength(2.)
|
|||||||
algo_local.Propagation()
|
algo_local.Propagation()
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
hexa.Compute()
|
if not hexa.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -40,4 +40,4 @@ algo = tria_mesh.Triangle()
|
|||||||
algo.MaxElementArea(100)
|
algo.MaxElementArea(100)
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
tria_mesh.Compute()
|
if not tria_mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -31,6 +31,6 @@ algo3D.MaxElementVolume(200.)
|
|||||||
# compute the mesh
|
# compute the mesh
|
||||||
ret = tetra.Compute()
|
ret = tetra.Compute()
|
||||||
if ret == 0:
|
if ret == 0:
|
||||||
print("problem when computing the mesh")
|
raise Exception("problem when computing the mesh")
|
||||||
else:
|
|
||||||
print("Computation succeeded")
|
print("Computation succeeded")
|
||||||
|
@ -33,4 +33,4 @@ algo2D = tria.Triangle()
|
|||||||
algo2D.LengthFromEdges()
|
algo2D.LengthFromEdges()
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
tria.Compute()
|
if not tria.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -39,7 +39,7 @@ algo_local.Arithmetic1D(1, 4)
|
|||||||
algo_local.Propagation()
|
algo_local.Propagation()
|
||||||
|
|
||||||
# compute the mesh which contains prisms
|
# compute the mesh which contains prisms
|
||||||
hexa.Compute()
|
if not hexa.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# create another mesh on the box
|
# create another mesh on the box
|
||||||
mesh = smesh_builder.Mesh(box, "Propagation of distribution of nodes")
|
mesh = smesh_builder.Mesh(box, "Propagation of distribution of nodes")
|
||||||
@ -60,4 +60,4 @@ algo_local.Arithmetic1D(1, 4)
|
|||||||
algo_local.PropagationOfDistribution()
|
algo_local.PropagationOfDistribution()
|
||||||
|
|
||||||
# compute the mesh which contains hexahedra only
|
# compute the mesh which contains hexahedra only
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -29,7 +29,7 @@ algo3D = hexa.Hexahedron()
|
|||||||
algo1D.Arithmetic1D(1, 4)
|
algo1D.Arithmetic1D(1, 4)
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
hexa.Compute()
|
if not hexa.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# 2. Create a tetrahedral mesh on the box
|
# 2. Create a tetrahedral mesh on the box
|
||||||
tetra = smesh_builder.Mesh(box, "Box : tetrahedrical mesh")
|
tetra = smesh_builder.Mesh(box, "Box : tetrahedrical mesh")
|
||||||
@ -48,4 +48,4 @@ algo1D.Arithmetic1D(1, 4)
|
|||||||
algo2D.LengthFromEdges()
|
algo2D.LengthFromEdges()
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
tetra.Compute()
|
if not tetra.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -47,7 +47,7 @@ src_mesh.Segment().NumberOfSegments(9,10)
|
|||||||
src_mesh.Quadrangle()
|
src_mesh.Quadrangle()
|
||||||
src_mesh.Hexahedron()
|
src_mesh.Hexahedron()
|
||||||
src_mesh.Triangle(f1) # triangular sub-mesh
|
src_mesh.Triangle(f1) # triangular sub-mesh
|
||||||
src_mesh.Compute()
|
if not src_mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Mesh the box using projection algorithms
|
# Mesh the box using projection algorithms
|
||||||
|
|
||||||
@ -75,7 +75,7 @@ proj2D.SourceFace( f2 )
|
|||||||
# 3D hypotheses to project prisms from the source to the target mesh
|
# 3D hypotheses to project prisms from the source to the target mesh
|
||||||
proj3D = tgt_mesh.Projection3D()
|
proj3D = tgt_mesh.Projection3D()
|
||||||
proj3D.SourceShape3D( box, src_mesh, v1F1, v1F2, v2F1, v2F2 )
|
proj3D.SourceShape3D( box, src_mesh, v1F1, v1F2, v2F1, v2F2 )
|
||||||
tgt_mesh.Compute()
|
if not tgt_mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Move the source mesh to visually compare the two meshes
|
# Move the source mesh to visually compare the two meshes
|
||||||
src_mesh.TranslateObject( src_mesh, smesh_builder.MakeDirStruct( 210, 0, 0 ), Copy=False)
|
src_mesh.TranslateObject( src_mesh, smesh_builder.MakeDirStruct( 210, 0, 0 ), Copy=False)
|
||||||
|
@ -29,9 +29,9 @@ geom_builder.addToStudyInFather( box, Face_2, 'Face_2' )
|
|||||||
src_mesh = smesh_builder.Mesh(Face_1, "Source mesh")
|
src_mesh = smesh_builder.Mesh(Face_1, "Source mesh")
|
||||||
src_mesh.Segment().NumberOfSegments(15)
|
src_mesh.Segment().NumberOfSegments(15)
|
||||||
src_mesh.Triangle()
|
src_mesh.Triangle()
|
||||||
src_mesh.Compute()
|
if not src_mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Mesh the target mesh using the algorithm Projection1D2D
|
# Mesh the target mesh using the algorithm Projection1D2D
|
||||||
tgt_mesh = smesh_builder.Mesh(Face_2, "Target mesh")
|
tgt_mesh = smesh_builder.Mesh(Face_2, "Target mesh")
|
||||||
tgt_mesh.Projection1D2D().SourceFace(Face_1,src_mesh)
|
tgt_mesh.Projection1D2D().SourceFace(Face_1,src_mesh)
|
||||||
tgt_mesh.Compute()
|
if not tgt_mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -37,4 +37,4 @@ Mesh_1.AddHypothesis(Fixed_points_1D_1)
|
|||||||
Regular_1D = Mesh_1.Segment()
|
Regular_1D = Mesh_1.Segment()
|
||||||
Quadrangle_2D = Mesh_1.Quadrangle()
|
Quadrangle_2D = Mesh_1.Quadrangle()
|
||||||
# Compute mesh
|
# Compute mesh
|
||||||
Mesh_1.Compute()
|
if not Mesh_1.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -24,7 +24,7 @@ radial_Quad_algo = mesh.Quadrangle(algo=smeshBuilder.RADIAL_QUAD)
|
|||||||
# The Radial Quadrange algorithm can work without any hypothesis
|
# The Radial Quadrange algorithm can work without any hypothesis
|
||||||
# In this case it uses "Default Nb of Segments" preferences parameter to discretize edges
|
# In this case it uses "Default Nb of Segments" preferences parameter to discretize edges
|
||||||
# So by default there will be 15 segments in both radial and circular directions
|
# So by default there will be 15 segments in both radial and circular directions
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# The Radial Quadrange uses global or local 1d hypotheses if it does
|
# The Radial Quadrange uses global or local 1d hypotheses if it does
|
||||||
# not have its own hypotheses.
|
# not have its own hypotheses.
|
||||||
@ -32,9 +32,9 @@ mesh.Compute()
|
|||||||
# So that there will be 5 radial layers and 10 circular segments
|
# So that there will be 5 radial layers and 10 circular segments
|
||||||
global_Nb_Segments = mesh.Segment().NumberOfSegments(5)
|
global_Nb_Segments = mesh.Segment().NumberOfSegments(5)
|
||||||
local_Nb_Segments = mesh.Segment(circle).NumberOfSegments(10)
|
local_Nb_Segments = mesh.Segment(circle).NumberOfSegments(10)
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Define own parameters of Radial Quadrange algorithm
|
# Define own parameters of Radial Quadrange algorithm
|
||||||
# The number of radial layers will be 4
|
# The number of radial layers will be 4
|
||||||
radial_Quad_algo.NumberOfLayers( 4 )
|
radial_Quad_algo.NumberOfLayers( 4 )
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -28,4 +28,4 @@ Nb_Segments_1 = Regular_1D.NumberOfSegments(10)
|
|||||||
Quadrangle_2D = Mesh_1.Quadrangle().TriangleVertex( triaVertex )
|
Quadrangle_2D = Mesh_1.Quadrangle().TriangleVertex( triaVertex )
|
||||||
|
|
||||||
# Compute the mesh
|
# Compute the mesh
|
||||||
Mesh_1.Compute()
|
if not Mesh_1.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -36,8 +36,8 @@ Regular_1D_1 = Mesh_1.Segment(geom=Edge_2)
|
|||||||
Nb_Segments_2 = Regular_1D_1.NumberOfSegments(10)
|
Nb_Segments_2 = Regular_1D_1.NumberOfSegments(10)
|
||||||
|
|
||||||
# Compute mesh (with Quadrangle Preference type)
|
# Compute mesh (with Quadrangle Preference type)
|
||||||
isDone = Mesh_1.Compute()
|
if not Mesh_1.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Change type to Reduced and compute again
|
# Change type to Reduced and compute again
|
||||||
Quadrangle_Parameters_1.SetQuadType( smeshBuilder.QUAD_REDUCED )
|
Quadrangle_Parameters_1.SetQuadType( smeshBuilder.QUAD_REDUCED )
|
||||||
isDone = Mesh_1.Compute()
|
if not Mesh_1.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -33,7 +33,7 @@ geom_builder.addToStudyInFather( boxes[1], midFace1, "middle Face")
|
|||||||
srcMesh = smesh_builder.Mesh(boxes[0], "source mesh") # box coloser to CS origin
|
srcMesh = smesh_builder.Mesh(boxes[0], "source mesh") # box coloser to CS origin
|
||||||
nSeg1 = srcMesh.Segment().NumberOfSegments(4)
|
nSeg1 = srcMesh.Segment().NumberOfSegments(4)
|
||||||
srcMesh.Quadrangle()
|
srcMesh.Quadrangle()
|
||||||
srcMesh.Compute()
|
if not srcMesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
srcFaceGroup = srcMesh.GroupOnGeom( midFace0, "src faces", SMESH.FACE )
|
srcFaceGroup = srcMesh.GroupOnGeom( midFace0, "src faces", SMESH.FACE )
|
||||||
|
|
||||||
# Import faces from midFace0 to the target mesh
|
# Import faces from midFace0 to the target mesh
|
||||||
@ -43,8 +43,8 @@ importAlgo = tgtMesh.UseExisting2DElements(midFace1)
|
|||||||
import2hyp = importAlgo.SourceFaces( [srcFaceGroup] )
|
import2hyp = importAlgo.SourceFaces( [srcFaceGroup] )
|
||||||
tgtMesh.Segment().NumberOfSegments(3)
|
tgtMesh.Segment().NumberOfSegments(3)
|
||||||
tgtMesh.Quadrangle()
|
tgtMesh.Quadrangle()
|
||||||
tgtMesh.Compute()
|
if not tgtMesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Import the whole source mesh with groups
|
# Import the whole source mesh with groups
|
||||||
import2hyp.SetCopySourceMesh(True,True)
|
import2hyp.SetCopySourceMesh(True,True)
|
||||||
tgtMesh.Compute()
|
if not tgtMesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -43,7 +43,7 @@ layersHyp = algo3D.ViscousLayers(thickness,numberOfLayers,stretchFactor,
|
|||||||
ignoreFaces, # optional
|
ignoreFaces, # optional
|
||||||
groupName = groupName) # optional
|
groupName = groupName) # optional
|
||||||
|
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# retrieve boundary prisms created by mesh.Compute()
|
# retrieve boundary prisms created by mesh.Compute()
|
||||||
boundaryGroup = mesh.GetGroupByName( layersHyp.GetGroupName() )[0]
|
boundaryGroup = mesh.GetGroupByName( layersHyp.GetGroupName() )[0]
|
||||||
@ -65,7 +65,7 @@ mesh.Segment().NumberOfSegments( 5 )
|
|||||||
vlHyp = mesh.Triangle().ViscousLayers2D( 2, 3, 1.5,
|
vlHyp = mesh.Triangle().ViscousLayers2D( 2, 3, 1.5,
|
||||||
edgeIds, isEdgesToIgnore=True, # optional
|
edgeIds, isEdgesToIgnore=True, # optional
|
||||||
groupName=groupName) # optional
|
groupName=groupName) # optional
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# retrieve boundary elements created by mesh.Compute()
|
# retrieve boundary elements created by mesh.Compute()
|
||||||
quadrangles = mesh.GetGroupByName( vlHyp.GetGroupName() )[0]
|
quadrangles = mesh.GetGroupByName( vlHyp.GetGroupName() )[0]
|
||||||
@ -74,4 +74,4 @@ print( "Nb boundary quadrangles", quadrangles.Size() )
|
|||||||
# viscous layers will be created on 3 edges, as we pass isEdgesToIgnore=False
|
# viscous layers will be created on 3 edges, as we pass isEdgesToIgnore=False
|
||||||
vlHyp.SetEdges( edgeIds, False )
|
vlHyp.SetEdges( edgeIds, False )
|
||||||
|
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -28,4 +28,4 @@ mesh.Hexahedron()
|
|||||||
# define refinement near vertices
|
# define refinement near vertices
|
||||||
algo1d.LengthNearVertex( 0.2 )
|
algo1d.LengthNearVertex( 0.2 )
|
||||||
|
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -15,7 +15,7 @@ face = geom_builder.MakeFaceHW(100, 100, 1, theName="quadrangle")
|
|||||||
mesh = smesh_builder.Mesh(face)
|
mesh = smesh_builder.Mesh(face)
|
||||||
mesh.Segment().NumberOfSegments(10)
|
mesh.Segment().NumberOfSegments(10)
|
||||||
mesh.Triangle().MaxElementArea(25)
|
mesh.Triangle().MaxElementArea(25)
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# get all free borders
|
# get all free borders
|
||||||
filter = smesh_builder.GetFilter(SMESH.EDGE, SMESH.FT_FreeBorders)
|
filter = smesh_builder.GetFilter(SMESH.EDGE, SMESH.FT_FreeBorders)
|
||||||
|
@ -16,7 +16,7 @@ geom_builder.addToStudy( face, "quadrangle" )
|
|||||||
mesh = smesh_builder.Mesh(face)
|
mesh = smesh_builder.Mesh(face)
|
||||||
mesh.Segment().NumberOfSegments(10)
|
mesh.Segment().NumberOfSegments(10)
|
||||||
mesh.Triangle().MaxElementArea(25)
|
mesh.Triangle().MaxElementArea(25)
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# get all faces with free edges
|
# get all faces with free edges
|
||||||
filter = smesh_builder.GetFilter(SMESH.FACE, SMESH.FT_FreeEdges)
|
filter = smesh_builder.GetFilter(SMESH.FACE, SMESH.FT_FreeEdges)
|
||||||
|
@ -16,7 +16,7 @@ mesh = smesh_builder.Mesh( box, "Box" )
|
|||||||
mesh.Segment().NumberOfSegments(10)
|
mesh.Segment().NumberOfSegments(10)
|
||||||
mesh.Quadrangle()
|
mesh.Quadrangle()
|
||||||
mesh.Hexahedron()
|
mesh.Hexahedron()
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
# copy all elements with translation and Merge nodes
|
# copy all elements with translation and Merge nodes
|
||||||
mesh.TranslateObject( mesh, smesh_builder.MakeDirStruct( 10,0,0), Copy=True )
|
mesh.TranslateObject( mesh, smesh_builder.MakeDirStruct( 10,0,0), Copy=True )
|
||||||
mesh.MergeNodes( mesh.FindCoincidentNodes(1e-7) )
|
mesh.MergeNodes( mesh.FindCoincidentNodes(1e-7) )
|
||||||
|
@ -16,7 +16,7 @@ mesh = smesh_builder.Mesh( box, "Box" )
|
|||||||
mesh.Segment().NumberOfSegments(10)
|
mesh.Segment().NumberOfSegments(10)
|
||||||
mesh.Quadrangle()
|
mesh.Quadrangle()
|
||||||
mesh.Hexahedron()
|
mesh.Hexahedron()
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
# copy all elements with translation
|
# copy all elements with translation
|
||||||
mesh.TranslateObject( mesh, [10,0,0], Copy=True )
|
mesh.TranslateObject( mesh, [10,0,0], Copy=True )
|
||||||
# create a filter to find nodes equal within tolerance of 1e-5
|
# create a filter to find nodes equal within tolerance of 1e-5
|
||||||
|
@ -16,7 +16,7 @@ mesh = smesh_builder.Mesh( box, "Box" )
|
|||||||
mesh.Segment().NumberOfSegments(10)
|
mesh.Segment().NumberOfSegments(10)
|
||||||
mesh.Quadrangle()
|
mesh.Quadrangle()
|
||||||
mesh.Hexahedron()
|
mesh.Hexahedron()
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
# copy all elements with translation and merge nodes
|
# copy all elements with translation and merge nodes
|
||||||
mesh.TranslateObject( mesh, [10,0,0], Copy=True )
|
mesh.TranslateObject( mesh, [10,0,0], Copy=True )
|
||||||
mesh.MergeNodes( mesh.FindCoincidentNodes( 1e-5 ))
|
mesh.MergeNodes( mesh.FindCoincidentNodes( 1e-5 ))
|
||||||
|
@ -34,7 +34,7 @@ v000 = geom_builder.MakeVertex( 100,100,0, theName='v000' ) # can use box sub-ve
|
|||||||
v001 = geom_builder.GetVertexNearPoint( box1, geom_builder.MakeVertex(100,100,100), theName='v001')
|
v001 = geom_builder.GetVertexNearPoint( box1, geom_builder.MakeVertex(100,100,100), theName='v001')
|
||||||
ijkAlgo.Renumber([ smeshBuilder.BlockCS( box1, v000, v001 ) ])
|
ijkAlgo.Renumber([ smeshBuilder.BlockCS( box1, v000, v001 ) ])
|
||||||
|
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Create filters with FT_ConnectedElements criterion by pointing a domain in different ways:
|
# Create filters with FT_ConnectedElements criterion by pointing a domain in different ways:
|
||||||
|
|
||||||
|
@ -44,7 +44,7 @@ Mesh_1 = smesh_builder.Mesh(Partition_1)
|
|||||||
Mesh_1.Segment().NumberOfSegments(15)
|
Mesh_1.Segment().NumberOfSegments(15)
|
||||||
Mesh_1.Triangle().LengthFromEdges()
|
Mesh_1.Triangle().LengthFromEdges()
|
||||||
Mesh_1.Tetrahedron()
|
Mesh_1.Tetrahedron()
|
||||||
Mesh_1.Compute()
|
if not Mesh_1.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# relevant groups of volumes and faces
|
# relevant groups of volumes and faces
|
||||||
|
|
||||||
|
@ -33,7 +33,7 @@ quadra.Quadrangle()
|
|||||||
algo1D.NumberOfSegments(7)
|
algo1D.NumberOfSegments(7)
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
quadra.Compute()
|
if not quadra.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# create SMESH group on the face with name "SMESHGroup1"
|
# create SMESH group on the face with name "SMESHGroup1"
|
||||||
aSmeshGroup1 = quadra.GroupOnGeom(face, "SMESHGroup1")
|
aSmeshGroup1 = quadra.GroupOnGeom(face, "SMESHGroup1")
|
||||||
|
@ -16,7 +16,7 @@ box = geom_builder.MakeBoxDXDYDZ(10,10,10)
|
|||||||
mesh = smesh_builder.Mesh(box,"Quad mesh")
|
mesh = smesh_builder.Mesh(box,"Quad mesh")
|
||||||
hyp1D = mesh.Segment().StartEndLength( 1, 4 )
|
hyp1D = mesh.Segment().StartEndLength( 1, 4 )
|
||||||
mesh.Quadrangle()
|
mesh.Quadrangle()
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# create a group on filter selecting faces of medium size
|
# create a group on filter selecting faces of medium size
|
||||||
critaria = [ \
|
critaria = [ \
|
||||||
@ -30,7 +30,7 @@ print("Group on filter contains %s elements" % filtGroup.Size())
|
|||||||
# group on filter is updated if the mesh is modified
|
# group on filter is updated if the mesh is modified
|
||||||
hyp1D.SetStartLength( 2.5 )
|
hyp1D.SetStartLength( 2.5 )
|
||||||
hyp1D.SetEndLength( 2.5 )
|
hyp1D.SetEndLength( 2.5 )
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
print("After mesh change, group on filter contains %s elements" % filtGroup.Size())
|
print("After mesh change, group on filter contains %s elements" % filtGroup.Size())
|
||||||
|
|
||||||
# set a new filter defining the group
|
# set a new filter defining the group
|
||||||
|
@ -27,8 +27,8 @@ mesh2 = smesh_builder.Mesh(box2, 'box2')
|
|||||||
mesh2.Segment().NumberOfSegments(5)
|
mesh2.Segment().NumberOfSegments(5)
|
||||||
mesh2.Triangle().MaxElementArea(20)
|
mesh2.Triangle().MaxElementArea(20)
|
||||||
|
|
||||||
mesh1.Compute()
|
if not mesh1.Compute(): raise Exception("Error when computing Mesh")
|
||||||
mesh2.Compute()
|
if not mesh2.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# compute min distance from grp1 to the origin (not available yet)
|
# compute min distance from grp1 to the origin (not available yet)
|
||||||
smesh_builder.MinDistance(mesh1)
|
smesh_builder.MinDistance(mesh1)
|
||||||
|
@ -27,8 +27,8 @@ mesh2 = smesh_builder.Mesh(box2, 'box2')
|
|||||||
mesh2.Segment().NumberOfSegments(5)
|
mesh2.Segment().NumberOfSegments(5)
|
||||||
mesh2.Triangle().MaxElementArea(20)
|
mesh2.Triangle().MaxElementArea(20)
|
||||||
|
|
||||||
mesh1.Compute()
|
if not mesh1.Compute(): raise Exception("Error when computing Mesh")
|
||||||
mesh2.Compute()
|
if not mesh2.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# compute bounding box for mesh1
|
# compute bounding box for mesh1
|
||||||
mesh1.BoundingBox()
|
mesh1.BoundingBox()
|
||||||
|
@ -19,7 +19,7 @@ submesh_1d = mesh.Segment().NumberOfSegments(5)
|
|||||||
submesh_2d = mesh.Quadrangle()
|
submesh_2d = mesh.Quadrangle()
|
||||||
submesh_3d = mesh.Hexahedron()
|
submesh_3d = mesh.Hexahedron()
|
||||||
submesh_2d_face = mesh.Triangle(face)
|
submesh_2d_face = mesh.Triangle(face)
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# create a group
|
# create a group
|
||||||
|
|
||||||
|
@ -83,4 +83,4 @@ mesh.Quadrangle(sub_face4)
|
|||||||
|
|
||||||
# finally compute whole mesh
|
# finally compute whole mesh
|
||||||
|
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -73,7 +73,7 @@ algo3d = tetra.Tetrahedron()
|
|||||||
algo3d.MaxElementVolume(100)
|
algo3d.MaxElementVolume(100)
|
||||||
|
|
||||||
# Compute the mesh
|
# Compute the mesh
|
||||||
tetra.Compute()
|
if not tetra.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Create a mesh group of all triangles generated on geom faces present in faces_group
|
# Create a mesh group of all triangles generated on geom faces present in faces_group
|
||||||
group = tetra.Group(faces_group)
|
group = tetra.Group(faces_group)
|
||||||
|
@ -16,7 +16,7 @@ mesh = smesh_builder.Mesh(box, 'box')
|
|||||||
mesh.Segment().NumberOfSegments(10)
|
mesh.Segment().NumberOfSegments(10)
|
||||||
mesh.Triangle().MaxElementArea(5)
|
mesh.Triangle().MaxElementArea(5)
|
||||||
|
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
print("After Compute(): %s nodes, %s faces" % ( mesh.NbNodes(), mesh.NbFaces()))
|
print("After Compute(): %s nodes, %s faces" % ( mesh.NbNodes(), mesh.NbFaces()))
|
||||||
|
|
||||||
# remove nodes #246 and #255
|
# remove nodes #246 and #255
|
||||||
|
@ -14,7 +14,7 @@ box = geom_builder.MakeBoxDXDYDZ(200, 200, 200)
|
|||||||
mesh = smesh_builder.Mesh( box )
|
mesh = smesh_builder.Mesh( box )
|
||||||
mesh.Segment().AutomaticLength(0.1)
|
mesh.Segment().AutomaticLength(0.1)
|
||||||
mesh.Quadrangle()
|
mesh.Quadrangle()
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# find node at (0,0,0) which is located on a geom vertex
|
# find node at (0,0,0) which is located on a geom vertex
|
||||||
node000 = None
|
node000 = None
|
||||||
|
@ -46,7 +46,7 @@ def Mesh1D(shape1d, nbSeg, name):
|
|||||||
algo = mesh1d_tool.Segment()
|
algo = mesh1d_tool.Segment()
|
||||||
hyp = algo.NumberOfSegments(nbSeg)
|
hyp = algo.NumberOfSegments(nbSeg)
|
||||||
isDone = mesh1d_tool.Compute()
|
isDone = mesh1d_tool.Compute()
|
||||||
if not isDone: print('Mesh ', name, ': computation failed')
|
if not isDone: raise Exception('Mesh %s: computation failed' % name)
|
||||||
return mesh1d_tool
|
return mesh1d_tool
|
||||||
|
|
||||||
# Create a mesh with six nodes, seven edges and two quadrangle faces
|
# Create a mesh with six nodes, seven edges and two quadrangle faces
|
||||||
|
@ -28,7 +28,7 @@ algo1D.NumberOfSegments(3)
|
|||||||
Mesh_1.Quadrangle()
|
Mesh_1.Quadrangle()
|
||||||
|
|
||||||
isDone = Mesh_1.Compute()
|
isDone = Mesh_1.Compute()
|
||||||
if not isDone: print('Mesh Mesh_1 : computation failed')
|
if not isDone: raise Exception('Mesh Mesh_1 : computation failed')
|
||||||
|
|
||||||
# build a triangle mesh on Face_2
|
# build a triangle mesh on Face_2
|
||||||
Mesh_2 = smesh_builder.Mesh(Face_2)
|
Mesh_2 = smesh_builder.Mesh(Face_2)
|
||||||
@ -39,7 +39,7 @@ algo2D = Mesh_2.Triangle()
|
|||||||
algo2D.MaxElementArea(240)
|
algo2D.MaxElementArea(240)
|
||||||
|
|
||||||
isDone = Mesh_2.Compute()
|
isDone = Mesh_2.Compute()
|
||||||
if not isDone: print('Mesh Mesh_2 : computation failed')
|
if not isDone: raise Exception('Mesh Mesh_2 : computation failed')
|
||||||
|
|
||||||
# create a 2d pattern
|
# create a 2d pattern
|
||||||
pattern = smesh_builder.GetPattern()
|
pattern = smesh_builder.GetPattern()
|
||||||
@ -60,7 +60,7 @@ Mesh_3.Segment().NumberOfSegments(1)
|
|||||||
Mesh_3.Quadrangle()
|
Mesh_3.Quadrangle()
|
||||||
Mesh_3.Hexahedron()
|
Mesh_3.Hexahedron()
|
||||||
isDone = Mesh_3.Compute()
|
isDone = Mesh_3.Compute()
|
||||||
if not isDone: print('Mesh Mesh_3 : computation failed')
|
if not isDone: raise Exception('Mesh Mesh_3 : computation failed')
|
||||||
|
|
||||||
# create a 3d pattern (hexahedrons)
|
# create a 3d pattern (hexahedrons)
|
||||||
pattern_hexa = smesh_builder.GetPattern()
|
pattern_hexa = smesh_builder.GetPattern()
|
||||||
@ -103,7 +103,7 @@ Mesh_4.Segment().NumberOfSegments(1)
|
|||||||
Mesh_4.Quadrangle()
|
Mesh_4.Quadrangle()
|
||||||
Mesh_4.Hexahedron()
|
Mesh_4.Hexahedron()
|
||||||
isDone = Mesh_4.Compute()
|
isDone = Mesh_4.Compute()
|
||||||
if not isDone: print('Mesh Mesh_4 : computation failed')
|
if not isDone: raise Exception('Mesh Mesh_4 : computation failed')
|
||||||
|
|
||||||
# create another 3d pattern (pyramids)
|
# create another 3d pattern (pyramids)
|
||||||
pattern_pyra = smesh_builder.GetPattern()
|
pattern_pyra = smesh_builder.GetPattern()
|
||||||
|
@ -23,7 +23,7 @@ Mesh.Tetrahedron()
|
|||||||
|
|
||||||
# compute mesh
|
# compute mesh
|
||||||
|
|
||||||
Mesh.Compute()
|
if not Mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# convert to quadratic
|
# convert to quadratic
|
||||||
# theForce3d = 1; this results in the medium node lying at the
|
# theForce3d = 1; this results in the medium node lying at the
|
||||||
|
@ -37,7 +37,7 @@ algo3D = tetra.Tetrahedron()
|
|||||||
algo3D.MaxElementVolume("MaxElementVolume")
|
algo3D.MaxElementVolume("MaxElementVolume")
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
ret = tetra.Compute()
|
if not tetra.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# translate the mesh
|
# translate the mesh
|
||||||
point = SMESH.PointStruct("Offset", 0., 0.)
|
point = SMESH.PointStruct("Offset", 0., 0.)
|
||||||
|
@ -72,4 +72,4 @@ mesh.Segment(bigQuad).LocalLength( 3 )
|
|||||||
mesh.Triangle(bigQuad)
|
mesh.Triangle(bigQuad)
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -29,4 +29,4 @@ algo1d.NumberOfSegments( circNbSeg ) # division of circle edges
|
|||||||
algo2d = mesh.Quadrangle( smeshBuilder.QUAD_MA_PROJ )
|
algo2d = mesh.Quadrangle( smeshBuilder.QUAD_MA_PROJ )
|
||||||
algo2d.StartEndLength( circleLen2 / circNbSeg, circleLen1 / circNbSeg ) # radial division
|
algo2d.StartEndLength( circleLen2 / circNbSeg, circleLen1 / circNbSeg ) # radial division
|
||||||
|
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -20,7 +20,7 @@ mesh = smesh_builder.Mesh( Fuse, "Deflection_2D")
|
|||||||
algo = mesh.Segment()
|
algo = mesh.Segment()
|
||||||
algo.LocalLength(35)
|
algo.LocalLength(35)
|
||||||
algo = mesh.Triangle()
|
algo = mesh.Triangle()
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# get min and max deflection
|
# get min and max deflection
|
||||||
minMax = mesh.GetMinMax( SMESH.FT_Deflection2D )
|
minMax = mesh.GetMinMax( SMESH.FT_Deflection2D )
|
||||||
|
@ -23,7 +23,7 @@ algo = mesh.Segment()
|
|||||||
algo.NumberOfSegments(5)
|
algo.NumberOfSegments(5)
|
||||||
algo = mesh.Triangle()
|
algo = mesh.Triangle()
|
||||||
algo.MaxElementArea(20.)
|
algo.MaxElementArea(20.)
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# criterion : free borders
|
# criterion : free borders
|
||||||
aFilter = smesh_builder.GetFilter(SMESH.EDGE, SMESH.FT_FreeBorders)
|
aFilter = smesh_builder.GetFilter(SMESH.EDGE, SMESH.FT_FreeBorders)
|
||||||
|
@ -23,7 +23,7 @@ algo = mesh.Segment()
|
|||||||
algo.NumberOfSegments(5)
|
algo.NumberOfSegments(5)
|
||||||
algo = mesh.Triangle()
|
algo = mesh.Triangle()
|
||||||
algo.MaxElementArea(20.)
|
algo.MaxElementArea(20.)
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Criterion : Borders at multi-connection
|
# Criterion : Borders at multi-connection
|
||||||
nb_conn = 2
|
nb_conn = 2
|
||||||
|
@ -23,7 +23,7 @@ algo = mesh.Segment()
|
|||||||
algo.NumberOfSegments(5)
|
algo.NumberOfSegments(5)
|
||||||
algo = mesh.Triangle()
|
algo = mesh.Triangle()
|
||||||
algo.MaxElementArea(20.)
|
algo.MaxElementArea(20.)
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Criterion : Length > 3.
|
# Criterion : Length > 3.
|
||||||
length_margin = 3.
|
length_margin = 3.
|
||||||
|
@ -18,7 +18,7 @@ idbox = geom_builder.addToStudy(box, "box")
|
|||||||
mesh = smesh_builder.Mesh(box, "Mesh_free_nodes")
|
mesh = smesh_builder.Mesh(box, "Mesh_free_nodes")
|
||||||
mesh.Segment().NumberOfSegments(10)
|
mesh.Segment().NumberOfSegments(10)
|
||||||
mesh.Triangle().MaxElementArea(150.)
|
mesh.Triangle().MaxElementArea(150.)
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Remove some elements to obtain free nodes
|
# Remove some elements to obtain free nodes
|
||||||
# Criterion : AREA < 80.
|
# Criterion : AREA < 80.
|
||||||
|
@ -42,7 +42,7 @@ Mesh_1 = smesh_builder.Mesh(Partition_1)
|
|||||||
Mesh_1.Segment().MaxSize(34.641)
|
Mesh_1.Segment().MaxSize(34.641)
|
||||||
Mesh_1.Triangle()
|
Mesh_1.Triangle()
|
||||||
Mesh_1.Tetrahedron()
|
Mesh_1.Tetrahedron()
|
||||||
Mesh_1.Compute()
|
if not Mesh_1.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# create a group of free faces
|
# create a group of free faces
|
||||||
aFilter = smesh_builder.GetFilter(SMESH.FACE, SMESH.FT_FreeFaces )
|
aFilter = smesh_builder.GetFilter(SMESH.FACE, SMESH.FT_FreeFaces )
|
||||||
|
@ -16,7 +16,7 @@ geom_builder.addToStudy( box, "box" )
|
|||||||
mesh = smesh_builder.Mesh(box)
|
mesh = smesh_builder.Mesh(box)
|
||||||
mesh.Segment().NumberOfSegments(3)
|
mesh.Segment().NumberOfSegments(3)
|
||||||
mesh.Quadrangle()
|
mesh.Quadrangle()
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# remove 2 faces
|
# remove 2 faces
|
||||||
allFaces = mesh.GetElementsByType(SMESH.FACE)
|
allFaces = mesh.GetElementsByType(SMESH.FACE)
|
||||||
|
@ -23,7 +23,7 @@ algo = mesh.Segment()
|
|||||||
algo.NumberOfSegments(5)
|
algo.NumberOfSegments(5)
|
||||||
algo = mesh.Triangle()
|
algo = mesh.Triangle()
|
||||||
algo.MaxElementArea(20.)
|
algo.MaxElementArea(20.)
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Criterion : Length 2D > 5.7
|
# Criterion : Length 2D > 5.7
|
||||||
length_margin = 5.7
|
length_margin = 5.7
|
||||||
|
@ -23,7 +23,7 @@ algo = mesh.Segment()
|
|||||||
algo.NumberOfSegments(5)
|
algo.NumberOfSegments(5)
|
||||||
algo = mesh.Triangle()
|
algo = mesh.Triangle()
|
||||||
algo.MaxElementArea(20.)
|
algo.MaxElementArea(20.)
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Criterion : MULTI-CONNECTION 2D = 3
|
# Criterion : MULTI-CONNECTION 2D = 3
|
||||||
nb_conn = 3
|
nb_conn = 3
|
||||||
|
@ -33,4 +33,4 @@ mesh.Projection1D2D( faces[1] ).SourceFace( faces[0] ) # projection faces[0] ->
|
|||||||
prism_algo.NumberOfSegments( 4, 5. )
|
prism_algo.NumberOfSegments( 4, 5. )
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -24,7 +24,7 @@ mesh = smesh_builder.Mesh( shape )
|
|||||||
mesh.Segment().NumberOfSegments(1)
|
mesh.Segment().NumberOfSegments(1)
|
||||||
mesh.Quadrangle()
|
mesh.Quadrangle()
|
||||||
mesh.Triangle( ff[1] )
|
mesh.Triangle( ff[1] )
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# make group of quadrangles and extrude them into a hexahedron
|
# make group of quadrangles and extrude them into a hexahedron
|
||||||
quadGroup = mesh.Group( ff[0], "quads")
|
quadGroup = mesh.Group( ff[0], "quads")
|
||||||
|
@ -40,7 +40,7 @@ smeshhomard.SetLogInFile(True)
|
|||||||
smeshhomard.SetLogFile(log_file_1)
|
smeshhomard.SetLogFile(log_file_1)
|
||||||
smeshhomard.SetRemoveLogOnSuccess(False)
|
smeshhomard.SetRemoveLogOnSuccess(False)
|
||||||
smeshhomard.SetVerboseLevel(3)
|
smeshhomard.SetVerboseLevel(3)
|
||||||
smeshhomard.Compute()
|
if not smeshhomard.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
if osp.isfile(output_med_1):
|
if osp.isfile(output_med_1):
|
||||||
os.remove(output_med_1)
|
os.remove(output_med_1)
|
||||||
@ -79,7 +79,7 @@ smeshhomard.SetLogInFile(True)
|
|||||||
smeshhomard.SetLogFile(log_file_2)
|
smeshhomard.SetLogFile(log_file_2)
|
||||||
smeshhomard.SetRemoveLogOnSuccess(True)
|
smeshhomard.SetRemoveLogOnSuccess(True)
|
||||||
smeshhomard.SetVerboseLevel(0)
|
smeshhomard.SetVerboseLevel(0)
|
||||||
smeshhomard.Compute()
|
if not smeshhomard.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
if osp.isfile(output_med_2):
|
if osp.isfile(output_med_2):
|
||||||
print("Test Uniform refinement Case 2: Error: output med file has not been removed")
|
print("Test Uniform refinement Case 2: Error: output med file has not been removed")
|
||||||
|
@ -42,7 +42,7 @@ smeshhomard.SetLogInFile(True)
|
|||||||
smeshhomard.SetLogFile(log_file_1)
|
smeshhomard.SetLogFile(log_file_1)
|
||||||
smeshhomard.SetRemoveLogOnSuccess(False)
|
smeshhomard.SetRemoveLogOnSuccess(False)
|
||||||
smeshhomard.SetVerboseLevel(3)
|
smeshhomard.SetVerboseLevel(3)
|
||||||
smeshhomard.Compute()
|
if not smeshhomard.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
if osp.isfile(output_med_1):
|
if osp.isfile(output_med_1):
|
||||||
os.remove(output_med_1)
|
os.remove(output_med_1)
|
||||||
@ -85,7 +85,7 @@ smeshhomard.SetLogInFile(True)
|
|||||||
smeshhomard.SetLogFile(log_file_2)
|
smeshhomard.SetLogFile(log_file_2)
|
||||||
smeshhomard.SetRemoveLogOnSuccess(True)
|
smeshhomard.SetRemoveLogOnSuccess(True)
|
||||||
smeshhomard.SetVerboseLevel(0)
|
smeshhomard.SetVerboseLevel(0)
|
||||||
smeshhomard.Compute()
|
if not smeshhomard.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
if osp.isfile(output_med_2):
|
if osp.isfile(output_med_2):
|
||||||
print("Test Uniform refinement Case 2: Error: output med file has not been removed")
|
print("Test Uniform refinement Case 2: Error: output med file has not been removed")
|
||||||
|
@ -18,7 +18,7 @@ Regular_1D = Mesh1.Segment()
|
|||||||
Nb_Segments_1 = Regular_1D.NumberOfSegments(3)
|
Nb_Segments_1 = Regular_1D.NumberOfSegments(3)
|
||||||
Nb_Segments_1.SetDistrType( 0 )
|
Nb_Segments_1.SetDistrType( 0 )
|
||||||
Quadrangle_2D = Mesh1.Quadrangle()
|
Quadrangle_2D = Mesh1.Quadrangle()
|
||||||
isDone = Mesh1.Compute()
|
if not Mesh1.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
#Perform scale operation for the whole mesh and creation of a new mesh
|
#Perform scale operation for the whole mesh and creation of a new mesh
|
||||||
newMesh = Mesh1.ScaleMakeMesh(Mesh1,SMESH.PointStruct(100,100,200),[0.5,0.3,0.7],True,"ScaledMesh")
|
newMesh = Mesh1.ScaleMakeMesh(Mesh1,SMESH.PointStruct(100,100,200),[0.5,0.3,0.7],True,"ScaledMesh")
|
||||||
|
@ -40,7 +40,7 @@ algo1D.NumberOfSegments(6)
|
|||||||
algo2D = trias.Triangle()
|
algo2D = trias.Triangle()
|
||||||
algo2D.LengthFromEdges()
|
algo2D.LengthFromEdges()
|
||||||
|
|
||||||
trias.Compute()
|
if not trias.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# create a group of all triangles currently present in the mesh
|
# create a group of all triangles currently present in the mesh
|
||||||
faceTriGroup = trias.Group( face1, "face triangles" )
|
faceTriGroup = trias.Group( face1, "face triangles" )
|
||||||
@ -49,7 +49,7 @@ faceTriGroup = trias.Group( face1, "face triangles" )
|
|||||||
circlemesh = smesh_builder.Mesh(circle, "Path mesh")
|
circlemesh = smesh_builder.Mesh(circle, "Path mesh")
|
||||||
algo = circlemesh.Segment()
|
algo = circlemesh.Segment()
|
||||||
algo.NumberOfSegments(10)
|
algo.NumberOfSegments(10)
|
||||||
circlemesh.Compute()
|
if not circlemesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# extrusion of the mesh
|
# extrusion of the mesh
|
||||||
trias.ExtrusionAlongPath([], circlemesh, circle, 1, MakeGroups=True )
|
trias.ExtrusionAlongPath([], circlemesh, circle, 1, MakeGroups=True )
|
||||||
|
@ -35,7 +35,7 @@ algo_local = mesh.Segment(edge1)
|
|||||||
algo_local.Arithmetic1D(1, 4)
|
algo_local.Arithmetic1D(1, 4)
|
||||||
algo_local.Propagation()
|
algo_local.Propagation()
|
||||||
|
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# sew border to side
|
# sew border to side
|
||||||
# FirstNodeIDOnFreeBorder, SecondNodeIDOnFreeBorder, LastNodeIDOnFreeBorder,
|
# FirstNodeIDOnFreeBorder, SecondNodeIDOnFreeBorder, LastNodeIDOnFreeBorder,
|
||||||
|
@ -35,7 +35,7 @@ algo_local = mesh.Segment(edge1)
|
|||||||
algo_local.Arithmetic1D(1, 4)
|
algo_local.Arithmetic1D(1, 4)
|
||||||
algo_local.Propagation()
|
algo_local.Propagation()
|
||||||
|
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# sew conform free borders
|
# sew conform free borders
|
||||||
# FirstNodeID1, SecondNodeID1, LastNodeID1, FirstNodeID2, SecondNodeID2
|
# FirstNodeID1, SecondNodeID1, LastNodeID1, FirstNodeID2, SecondNodeID2
|
||||||
|
@ -24,7 +24,7 @@ mesh = smesh_builder.Mesh(shape, "Two faces : quadrangle mesh")
|
|||||||
mesh.Segment().Arithmetic1D( 0.1, 0.4 )
|
mesh.Segment().Arithmetic1D( 0.1, 0.4 )
|
||||||
mesh.Segment(q1).NumberOfSegments( 5 )
|
mesh.Segment(q1).NumberOfSegments( 5 )
|
||||||
mesh.Quadrangle()
|
mesh.Quadrangle()
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# sew free borders
|
# sew free borders
|
||||||
|
|
||||||
|
@ -30,7 +30,7 @@ algo_local = mesh.Segment(EdgesList[8])
|
|||||||
algo_local.NumberOfSegments(4)
|
algo_local.NumberOfSegments(4)
|
||||||
algo_local.Propagation()
|
algo_local.Propagation()
|
||||||
|
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# sew side elements
|
# sew side elements
|
||||||
|
|
||||||
|
@ -21,7 +21,7 @@ mesh.Quadrangle()
|
|||||||
mesh.Hexahedron()
|
mesh.Hexahedron()
|
||||||
|
|
||||||
# Compute mesh
|
# Compute mesh
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Duplicate nodes only
|
# Duplicate nodes only
|
||||||
|
|
||||||
|
@ -25,7 +25,7 @@ mesh = smesh_builder.Mesh( shape, "test_Reorient2D")
|
|||||||
mesh.AutomaticHexahedralization(0.5)
|
mesh.AutomaticHexahedralization(0.5)
|
||||||
localAlgo = mesh.Segment(faces[0])
|
localAlgo = mesh.Segment(faces[0])
|
||||||
localAlgo.NumberOfSegments( 11 )
|
localAlgo.NumberOfSegments( 11 )
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
group = mesh.Group( faces[1] )
|
group = mesh.Group( faces[1] )
|
||||||
|
|
||||||
vec = geom_builder.MakeVectorDXDYDZ( 1, 1, 1 )
|
vec = geom_builder.MakeVectorDXDYDZ( 1, 1, 1 )
|
||||||
|
@ -101,7 +101,7 @@ geom_builder.addToStudy( f2, "f2" )
|
|||||||
# compute 1D mesh
|
# compute 1D mesh
|
||||||
mesh = smesh_builder.Mesh( box )
|
mesh = smesh_builder.Mesh( box )
|
||||||
mesh.Segment().NumberOfSegments( 5 )
|
mesh.Segment().NumberOfSegments( 5 )
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# compute 2D mesh
|
# compute 2D mesh
|
||||||
mesh.Quadrangle()
|
mesh.Quadrangle()
|
||||||
@ -109,8 +109,8 @@ mesh.UseExistingFaces(f1) # UseExistingFaces() allows using my2DMeshing();
|
|||||||
mesh.UseExistingFaces(f2) # assign UseExistingFaces() BEFORE calling my2DMeshing()!
|
mesh.UseExistingFaces(f2) # assign UseExistingFaces() BEFORE calling my2DMeshing()!
|
||||||
my2DMeshing(f1, mesh)
|
my2DMeshing(f1, mesh)
|
||||||
my2DMeshing(f2, mesh)
|
my2DMeshing(f2, mesh)
|
||||||
assert mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# compute 3D mesh
|
# compute 3D mesh
|
||||||
mesh.Prism()
|
mesh.Prism()
|
||||||
assert mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
@ -32,11 +32,11 @@ Regular_1D_1_1 = tetra.Segment(geom=Face_1)
|
|||||||
Nb_Segments_1 = Regular_1D_1_1.NumberOfSegments(5)
|
Nb_Segments_1 = Regular_1D_1_1.NumberOfSegments(5)
|
||||||
Nb_Segments_1.SetDistrType( 0 )
|
Nb_Segments_1.SetDistrType( 0 )
|
||||||
Quadrangle_2D = tetra.Quadrangle(geom=Face_1)
|
Quadrangle_2D = tetra.Quadrangle(geom=Face_1)
|
||||||
isDone = tetra.Compute()
|
if not tetra.Compute(): raise Exception("Error when computing Mesh")
|
||||||
submesh = Regular_1D_1_1.GetSubMesh()
|
submesh = Regular_1D_1_1.GetSubMesh()
|
||||||
|
|
||||||
# compute the mesh
|
# compute the mesh
|
||||||
tetra.Compute()
|
if not tetra.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Creation of group
|
# Creation of group
|
||||||
group = tetra.CreateEmptyGroup( SMESH.FACE, 'Group' )
|
group = tetra.CreateEmptyGroup( SMESH.FACE, 'Group' )
|
||||||
|
@ -16,7 +16,7 @@ box = geom_builder.MakeBoxDXDYDZ(100,100,100)
|
|||||||
# Create a mesh
|
# Create a mesh
|
||||||
mesh = smesh_builder.Mesh(box,"Mesh")
|
mesh = smesh_builder.Mesh(box,"Mesh")
|
||||||
mesh.AutomaticHexahedralization()
|
mesh.AutomaticHexahedralization()
|
||||||
mesh.Compute()
|
if not mesh.Compute(): raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Create a point
|
# Create a point
|
||||||
x,y,z = 0, 0, 1
|
x,y,z = 0, 0, 1
|
||||||
|
@ -203,6 +203,8 @@ def TEST_createBoxMesh():
|
|||||||
Hexa_3D = smesh.CreateHypothesis('Hexa_3D')
|
Hexa_3D = smesh.CreateHypothesis('Hexa_3D')
|
||||||
status = boxmesh.AddHypothesis(Hexa_3D)
|
status = boxmesh.AddHypothesis(Hexa_3D)
|
||||||
isDone = boxmesh.Compute()
|
isDone = boxmesh.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
smesh.SetName(boxmesh.GetMesh(), 'boxmesh')
|
smesh.SetName(boxmesh.GetMesh(), 'boxmesh')
|
||||||
if salome.sg.hasDesktop():
|
if salome.sg.hasDesktop():
|
||||||
|
@ -100,8 +100,11 @@ smesh.SetName(hypArea200, "Max. Element Area")
|
|||||||
|
|
||||||
print("---------------------Compute the mesh")
|
print("---------------------Compute the mesh")
|
||||||
|
|
||||||
ret = mesh.Compute()
|
isDone = mesh.Compute()
|
||||||
print(ret)
|
print(isDone)
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
|
|
||||||
salome.sg.updateObjBrowser()
|
salome.sg.updateObjBrowser()
|
||||||
|
|
||||||
|
@ -76,8 +76,13 @@ hypNbSeg2 = algoReg2.NumberOfSegments(34)
|
|||||||
smesh.SetName(hypNbSeg2, "NumberOfSegments 2")
|
smesh.SetName(hypNbSeg2, "NumberOfSegments 2")
|
||||||
|
|
||||||
# compute meshes
|
# compute meshes
|
||||||
mesh1.Compute()
|
isDone = mesh1.Compute()
|
||||||
mesh2.Compute()
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
|
isDone = mesh2.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
|
|
||||||
# ---- udate object browser
|
# ---- udate object browser
|
||||||
|
@ -86,8 +86,13 @@ smesh.SetName(hypNbSeg2, "NumberOfSegments_" + str(numberOfSegments2))
|
|||||||
|
|
||||||
|
|
||||||
# compute meshes
|
# compute meshes
|
||||||
mesh1.Compute()
|
isDone = mesh1.Compute()
|
||||||
mesh2.Compute()
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
|
isDone = mesh2.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# ---- update object browser
|
# ---- update object browser
|
||||||
salome.sg.updateObjBrowser()
|
salome.sg.updateObjBrowser()
|
||||||
|
@ -29,7 +29,8 @@ def CheckBelongToGeomFilterOld(theMeshGen, theMesh, theShape, theSubShape, theEl
|
|||||||
aName = str(theSubShape)
|
aName = str(theSubShape)
|
||||||
geompy.addToStudyInFather(theShape,theSubShape,aName)
|
geompy.addToStudyInFather(theShape,theSubShape,aName)
|
||||||
|
|
||||||
theMeshGen.Compute(theMesh,theShape)
|
if not theMeshGen.Compute(theMesh,theShape):
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
aFilterMgr = theMeshGen.CreateFilterManager()
|
aFilterMgr = theMeshGen.CreateFilterManager()
|
||||||
aFilter = aFilterMgr.CreateFilter()
|
aFilter = aFilterMgr.CreateFilter()
|
||||||
@ -48,7 +49,9 @@ def CheckBelongToGeomFilter(theMesh, theShape, theSubShape, theElemType):
|
|||||||
aName = str(theSubShape)
|
aName = str(theSubShape)
|
||||||
geompy.addToStudyInFather(theShape,theSubShape,aName)
|
geompy.addToStudyInFather(theShape,theSubShape,aName)
|
||||||
|
|
||||||
theMesh.Compute()
|
if not theMesh.Compute():
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
aFilter = smesh.GetFilter(theElemType, SMESH.FT_BelongToGeom, theSubShape)
|
aFilter = smesh.GetFilter(theElemType, SMESH.FT_BelongToGeom, theSubShape)
|
||||||
return aFilter.GetElementsId(theMesh.GetMesh())
|
return aFilter.GetElementsId(theMesh.GetMesh())
|
||||||
|
|
||||||
|
@ -73,7 +73,9 @@ algo1D_1=Mesh_inf.Segment()
|
|||||||
algo1D_1.NumberOfSegments(10)
|
algo1D_1.NumberOfSegments(10)
|
||||||
algo2D_1=Mesh_inf.Quadrangle()
|
algo2D_1=Mesh_inf.Quadrangle()
|
||||||
algo3D_1=Mesh_inf.Hexahedron()
|
algo3D_1=Mesh_inf.Hexahedron()
|
||||||
Mesh_inf.Compute()
|
isDone = Mesh_inf.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# create a group on the top face
|
# create a group on the top face
|
||||||
Gsup1=Mesh_inf.Group(Fsup1, "Sup")
|
Gsup1=Mesh_inf.Group(Fsup1, "Sup")
|
||||||
@ -86,7 +88,9 @@ algo1D_2=Mesh_sup.Segment()
|
|||||||
algo1D_2.NumberOfSegments(5)
|
algo1D_2.NumberOfSegments(5)
|
||||||
algo2D_2=Mesh_sup.Quadrangle()
|
algo2D_2=Mesh_sup.Quadrangle()
|
||||||
algo3D_2=Mesh_sup.Hexahedron()
|
algo3D_2=Mesh_sup.Hexahedron()
|
||||||
Mesh_sup.Compute()
|
isDone = Mesh_sup.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# create a group on the top face
|
# create a group on the top face
|
||||||
Gsup2=Mesh_sup.Group(Fsup2, "Sup")
|
Gsup2=Mesh_sup.Group(Fsup2, "Sup")
|
||||||
|
@ -28,7 +28,9 @@ from SMESH_test1 import *
|
|||||||
|
|
||||||
# Compute the mesh created in SMESH_test1
|
# Compute the mesh created in SMESH_test1
|
||||||
|
|
||||||
mesh.Compute()
|
isDone = mesh.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Create geometry groups on plane:
|
# Create geometry groups on plane:
|
||||||
aGeomGroup1 = geompy.CreateGroup(face , geompy.ShapeType["FACE"])
|
aGeomGroup1 = geompy.CreateGroup(face , geompy.ShapeType["FACE"])
|
||||||
|
@ -34,7 +34,9 @@ from SMESH_test1 import *
|
|||||||
|
|
||||||
# Compute the mesh created in SMESH_test1
|
# Compute the mesh created in SMESH_test1
|
||||||
|
|
||||||
mesh.Compute()
|
isDone = mesh.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Create geometry groups on plane:
|
# Create geometry groups on plane:
|
||||||
aGeomGroup1 = geompy.CreateGroup(face , geompy.ShapeType["FACE"])
|
aGeomGroup1 = geompy.CreateGroup(face , geompy.ShapeType["FACE"])
|
||||||
@ -70,7 +72,10 @@ print("aGroupOnShell ids :", aGroupOnShell.GetListOfID())
|
|||||||
print(" ")
|
print(" ")
|
||||||
|
|
||||||
print("Re-compute mesh, contents of aGroupOnShell changes again:")
|
print("Re-compute mesh, contents of aGroupOnShell changes again:")
|
||||||
mesh.Compute()
|
isDone = mesh.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
print("aGroupOnShell size =", aGroupOnShell.Size())
|
print("aGroupOnShell size =", aGroupOnShell.Size())
|
||||||
print("aGroupOnShell ids :", aGroupOnShell.GetListOfID())
|
print("aGroupOnShell ids :", aGroupOnShell.GetListOfID())
|
||||||
|
|
||||||
|
@ -39,7 +39,9 @@ def BuildGroupLyingOn(theMesh, theElemType, theName, theShape):
|
|||||||
#Example
|
#Example
|
||||||
from SMESH_test1 import *
|
from SMESH_test1 import *
|
||||||
|
|
||||||
mesh.Compute()
|
isDone = mesh.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# First way
|
# First way
|
||||||
BuildGroupLyingOn(mesh.GetMesh(), SMESH.FACE, "Group of faces lying on edge #1", edge )
|
BuildGroupLyingOn(mesh.GetMesh(), SMESH.FACE, "Group of faces lying on edge #1", edge )
|
||||||
|
@ -147,7 +147,9 @@ smesh.SetName(hVolume, "MaxElementVolume_"+str(theMaxElementVolume))
|
|||||||
|
|
||||||
|
|
||||||
print("-------------------------- compute the mesh of the mechanic piece")
|
print("-------------------------- compute the mesh of the mechanic piece")
|
||||||
mesh.Compute()
|
isDone = mesh.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
print("Information about the Nut:")
|
print("Information about the Nut:")
|
||||||
print("Number of nodes : ", mesh.NbNodes())
|
print("Number of nodes : ", mesh.NbNodes())
|
||||||
|
@ -169,8 +169,9 @@ smesh.SetName(hypVolume, "MaxElementVolume_" + str(maxElementVolume))
|
|||||||
|
|
||||||
print("-------------------------- compute the mesh of alveole ")
|
print("-------------------------- compute the mesh of alveole ")
|
||||||
ret = mesh.Compute()
|
ret = mesh.Compute()
|
||||||
|
if not ret:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
if ret != 0:
|
|
||||||
log = mesh.GetLog(0) # no erase trace
|
log = mesh.GetLog(0) # no erase trace
|
||||||
# for linelog in log:
|
# for linelog in log:
|
||||||
# print(linelog)
|
# print(linelog)
|
||||||
@ -181,7 +182,5 @@ if ret != 0:
|
|||||||
print("Number of triangles : ", mesh.NbTriangles())
|
print("Number of triangles : ", mesh.NbTriangles())
|
||||||
print("Number of volumes : ", mesh.NbVolumes())
|
print("Number of volumes : ", mesh.NbVolumes())
|
||||||
print("Number of tetrahedrons: ", mesh.NbTetras())
|
print("Number of tetrahedrons: ", mesh.NbTetras())
|
||||||
else:
|
|
||||||
print("problem when computing the mesh")
|
|
||||||
|
|
||||||
salome.sg.updateObjBrowser()
|
salome.sg.updateObjBrowser()
|
||||||
|
@ -116,6 +116,8 @@ algo = my_hexa.Segment()
|
|||||||
algo.NumberOfSegments(NbSeg)
|
algo.NumberOfSegments(NbSeg)
|
||||||
my_hexa.Quadrangle()
|
my_hexa.Quadrangle()
|
||||||
my_hexa.Hexahedron()
|
my_hexa.Hexahedron()
|
||||||
my_hexa.Compute()
|
isDone = my_hexa.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
salome.sg.updateObjBrowser()
|
salome.sg.updateObjBrowser()
|
||||||
|
@ -68,6 +68,8 @@ alg3D.SetName("algo3D")
|
|||||||
|
|
||||||
# compute mesh
|
# compute mesh
|
||||||
|
|
||||||
box_mesh.Compute()
|
isDone = box_mesh.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
sg.updateObjBrowser()
|
sg.updateObjBrowser()
|
||||||
|
@ -124,7 +124,9 @@ smesh.SetName(hypVolume, "MaxElementVolume_" + str(maxElementVolume))
|
|||||||
print("-------------------------- compute shell")
|
print("-------------------------- compute shell")
|
||||||
ret = mesh.Compute()
|
ret = mesh.Compute()
|
||||||
print(ret)
|
print(ret)
|
||||||
if ret != 0:
|
if not ret:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
log = mesh.GetLog(0) # no erase trace
|
log = mesh.GetLog(0) # no erase trace
|
||||||
# for linelog in log:
|
# for linelog in log:
|
||||||
# print(linelog)
|
# print(linelog)
|
||||||
@ -135,7 +137,5 @@ if ret != 0:
|
|||||||
print("Number of triangles : ", mesh.NbTriangles())
|
print("Number of triangles : ", mesh.NbTriangles())
|
||||||
print("Number of volumes : ", mesh.NbVolumes())
|
print("Number of volumes : ", mesh.NbVolumes())
|
||||||
print("Number of tetrahedrons: ", mesh.NbTetras())
|
print("Number of tetrahedrons: ", mesh.NbTetras())
|
||||||
else:
|
|
||||||
print("probleme when computing the mesh")
|
|
||||||
|
|
||||||
salome.sg.updateObjBrowser()
|
salome.sg.updateObjBrowser()
|
||||||
|
@ -134,7 +134,9 @@ smesh.SetName(hypVolume, "MaxElementVolume_" + str(maxElementVolume))
|
|||||||
print("-------------------------- compute shell")
|
print("-------------------------- compute shell")
|
||||||
ret = mesh.Compute()
|
ret = mesh.Compute()
|
||||||
print(ret)
|
print(ret)
|
||||||
if ret != 0:
|
if not ret:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
log = mesh.GetLog(0) # no erase trace
|
log = mesh.GetLog(0) # no erase trace
|
||||||
# for linelog in log:
|
# for linelog in log:
|
||||||
# print(linelog)
|
# print(linelog)
|
||||||
@ -145,7 +147,5 @@ if ret != 0:
|
|||||||
print("Number of triangles : ", mesh.NbTriangles())
|
print("Number of triangles : ", mesh.NbTriangles())
|
||||||
print("Number of volumes : ", mesh.NbVolumes())
|
print("Number of volumes : ", mesh.NbVolumes())
|
||||||
print("Number of tetrahedrons: ", mesh.NbTetras())
|
print("Number of tetrahedrons: ", mesh.NbTetras())
|
||||||
else:
|
|
||||||
print("probleme when computing the mesh")
|
|
||||||
|
|
||||||
salome.sg.updateObjBrowser()
|
salome.sg.updateObjBrowser()
|
||||||
|
@ -93,7 +93,9 @@ smesh.SetName(hypVolume, "MaxElementVolume_" + str(maxElementVolume))
|
|||||||
print("-------------------------- compute the mesh of the box")
|
print("-------------------------- compute the mesh of the box")
|
||||||
ret = mesh.Compute()
|
ret = mesh.Compute()
|
||||||
print(ret)
|
print(ret)
|
||||||
if ret != 0:
|
if not ret:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
log = mesh.GetLog(0) # no erase trace
|
log = mesh.GetLog(0) # no erase trace
|
||||||
# for linelog in log:
|
# for linelog in log:
|
||||||
# print(linelog)
|
# print(linelog)
|
||||||
@ -104,7 +106,5 @@ if ret != 0:
|
|||||||
print("Number of triangles : ", mesh.NbTriangles())
|
print("Number of triangles : ", mesh.NbTriangles())
|
||||||
print("Number of volumes : ", mesh.NbVolumes())
|
print("Number of volumes : ", mesh.NbVolumes())
|
||||||
print("Number of tetrahedrons: ", mesh.NbTetras())
|
print("Number of tetrahedrons: ", mesh.NbTetras())
|
||||||
else:
|
|
||||||
print("probleme when computing the mesh")
|
|
||||||
|
|
||||||
salome.sg.updateObjBrowser()
|
salome.sg.updateObjBrowser()
|
||||||
|
@ -73,6 +73,8 @@ top_1 = Mesh_1.GroupOnGeom(top,'top',SMESH.FACE)
|
|||||||
middle_1 = Mesh_1.GroupOnGeom(middle,'middle',SMESH.FACE)
|
middle_1 = Mesh_1.GroupOnGeom(middle,'middle',SMESH.FACE)
|
||||||
bottom_1 = Mesh_1.GroupOnGeom(bottom,'bottom',SMESH.FACE)
|
bottom_1 = Mesh_1.GroupOnGeom(bottom,'bottom',SMESH.FACE)
|
||||||
isDone = Mesh_1.Compute()
|
isDone = Mesh_1.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
[ top_1, middle_1, bottom_1 ] = Mesh_1.GetGroups()
|
[ top_1, middle_1, bottom_1 ] = Mesh_1.GetGroups()
|
||||||
|
|
||||||
|
@ -83,6 +83,8 @@ algo3d = Mesh_1.Tetrahedron(algo=smeshBuilder.NETGEN_3D)
|
|||||||
algo3d.MaxElementVolume(0.0002)
|
algo3d.MaxElementVolume(0.0002)
|
||||||
|
|
||||||
isDone = Mesh_1.Compute()
|
isDone = Mesh_1.Compute()
|
||||||
|
if not isDone:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
# Create groups
|
# Create groups
|
||||||
d_geom_groups = {}
|
d_geom_groups = {}
|
||||||
|
@ -181,9 +181,10 @@ salome.sg.updateObjBrowser()
|
|||||||
print("-------------------------- compute the mesh of the volume")
|
print("-------------------------- compute the mesh of the volume")
|
||||||
|
|
||||||
ret = mesh.Compute()
|
ret = mesh.Compute()
|
||||||
|
|
||||||
print(ret)
|
print(ret)
|
||||||
if ret != 0:
|
if not ret:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
##log=mesh.GetLog(0) # no erase trace
|
##log=mesh.GetLog(0) # no erase trace
|
||||||
##for linelog in log:
|
##for linelog in log:
|
||||||
## print linelog
|
## print linelog
|
||||||
@ -194,7 +195,5 @@ if ret != 0:
|
|||||||
print("Number of triangles : ", mesh.NbTriangles())
|
print("Number of triangles : ", mesh.NbTriangles())
|
||||||
print("Number of volumes : ", mesh.NbVolumes())
|
print("Number of volumes : ", mesh.NbVolumes())
|
||||||
print("Number of tetrahedrons: ", mesh.NbTetras())
|
print("Number of tetrahedrons: ", mesh.NbTetras())
|
||||||
else:
|
|
||||||
print("problem when Computing the mesh")
|
|
||||||
|
|
||||||
salome.sg.updateObjBrowser()
|
salome.sg.updateObjBrowser()
|
||||||
|
@ -85,7 +85,9 @@ hexa3D.SetName("Hexa_3D")
|
|||||||
print("-------------------------- compute compshell")
|
print("-------------------------- compute compshell")
|
||||||
ret = mesh.Compute()
|
ret = mesh.Compute()
|
||||||
print(ret)
|
print(ret)
|
||||||
if ret != 0:
|
if not ret:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
log = mesh.GetLog(0) # no erase trace
|
log = mesh.GetLog(0) # no erase trace
|
||||||
# for linelog in log:
|
# for linelog in log:
|
||||||
# print(linelog)
|
# print(linelog)
|
||||||
@ -96,7 +98,5 @@ if ret != 0:
|
|||||||
print("Number of quadrangles : ", mesh.NbQuadrangles())
|
print("Number of quadrangles : ", mesh.NbQuadrangles())
|
||||||
print("Number of volumes : ", mesh.NbVolumes())
|
print("Number of volumes : ", mesh.NbVolumes())
|
||||||
print("Number of hexahedrons : ", mesh.NbHexas())
|
print("Number of hexahedrons : ", mesh.NbHexas())
|
||||||
else:
|
|
||||||
print("problem when Computing the mesh")
|
|
||||||
|
|
||||||
salome.sg.updateObjBrowser()
|
salome.sg.updateObjBrowser()
|
||||||
|
@ -64,7 +64,9 @@ netgen.SetFineness( smeshBuilder.Fine )
|
|||||||
print("-------------------------- compute mesh")
|
print("-------------------------- compute mesh")
|
||||||
ret = mesh.Compute()
|
ret = mesh.Compute()
|
||||||
print(ret)
|
print(ret)
|
||||||
if ret != 0:
|
if not ret:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
print("Information about the MeshcompShel:")
|
print("Information about the MeshcompShel:")
|
||||||
print("Number of nodes : ", mesh.GetMesh().NbNodes())
|
print("Number of nodes : ", mesh.GetMesh().NbNodes())
|
||||||
print("Number of edges : ", mesh.GetMesh().NbEdges())
|
print("Number of edges : ", mesh.GetMesh().NbEdges())
|
||||||
@ -73,7 +75,4 @@ if ret != 0:
|
|||||||
print("Number of volumes : ", mesh.GetMesh().NbVolumes())
|
print("Number of volumes : ", mesh.GetMesh().NbVolumes())
|
||||||
print("Number of tetrahedrons : ", mesh.GetMesh().NbTetras())
|
print("Number of tetrahedrons : ", mesh.GetMesh().NbTetras())
|
||||||
|
|
||||||
else:
|
|
||||||
print("problem when computing the mesh")
|
|
||||||
|
|
||||||
salome.sg.updateObjBrowser()
|
salome.sg.updateObjBrowser()
|
||||||
|
@ -96,7 +96,9 @@ smesh.SetName(hypVolume, "MaxElementVolume_" + str(maxElementVolume))
|
|||||||
print("-------------------------- compute compshell")
|
print("-------------------------- compute compshell")
|
||||||
ret = mesh.Compute(mesh)
|
ret = mesh.Compute(mesh)
|
||||||
print(ret)
|
print(ret)
|
||||||
if ret != 0:
|
if not ret:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
log = mesh.GetLog(0) # no erase trace
|
log = mesh.GetLog(0) # no erase trace
|
||||||
# for linelog in log:
|
# for linelog in log:
|
||||||
# print(linelog)
|
# print(linelog)
|
||||||
@ -108,7 +110,4 @@ if ret != 0:
|
|||||||
print("Number of volumes : ", mesh.NbVolumes())
|
print("Number of volumes : ", mesh.NbVolumes())
|
||||||
print("Number of tetrahedrons : ", mesh.NbTetras())
|
print("Number of tetrahedrons : ", mesh.NbTetras())
|
||||||
|
|
||||||
else:
|
|
||||||
print("problem when computing the mesh")
|
|
||||||
|
|
||||||
salome.sg.updateObjBrowser()
|
salome.sg.updateObjBrowser()
|
||||||
|
@ -94,7 +94,9 @@ smesh.SetName(hypLengthFromEdge,"LengthFromEdge")
|
|||||||
print("-------------------------- compute the skin flight")
|
print("-------------------------- compute the skin flight")
|
||||||
ret = mesh.Compute()
|
ret = mesh.Compute()
|
||||||
print(ret)
|
print(ret)
|
||||||
if ret != 0:
|
if not ret:
|
||||||
|
raise Exception("Error when computing Mesh")
|
||||||
|
|
||||||
log = mesh.GetLog(0) # no erase trace
|
log = mesh.GetLog(0) # no erase trace
|
||||||
# for linelog in log:
|
# for linelog in log:
|
||||||
# print(linelog)
|
# print(linelog)
|
||||||
@ -104,7 +106,5 @@ if ret != 0:
|
|||||||
print("Number of faces : ", mesh.NbFaces())
|
print("Number of faces : ", mesh.NbFaces())
|
||||||
print("Number of triangles : ", mesh.NbTriangles())
|
print("Number of triangles : ", mesh.NbTriangles())
|
||||||
print("Number of volumes : ", mesh.NbVolumes())
|
print("Number of volumes : ", mesh.NbVolumes())
|
||||||
else:
|
|
||||||
print("probleme when computing the mesh")
|
|
||||||
|
|
||||||
salome.sg.updateObjBrowser()
|
salome.sg.updateObjBrowser()
|
||||||
|
Some files were not shown because too many files have changed in this diff Show More
Loading…
Reference in New Issue
Block a user