mirror of
https://git.salome-platform.org/gitpub/modules/smesh.git
synced 2024-12-27 09:50:34 +05:00
112 lines
3.0 KiB
Python
112 lines
3.0 KiB
Python
# Usage of Body Fitting algorithm
|
|
|
|
|
|
import salome
|
|
salome.salome_init_without_session()
|
|
import GEOM
|
|
from salome.geom import geomBuilder
|
|
geompy = geomBuilder.New()
|
|
|
|
import SMESH, SALOMEDS
|
|
from salome.smesh import smeshBuilder
|
|
smesh = smeshBuilder.New()
|
|
|
|
|
|
# create a sphere
|
|
sphere = geompy.MakeSphereR( 50 )
|
|
|
|
# cut the sphere by a box
|
|
box = geompy.MakeBoxDXDYDZ( 100, 100, 100 )
|
|
partition = geompy.MakePartition([ sphere ], [ box ], theName="partition")
|
|
|
|
# create a mesh and assign a "Body Fitting" algo
|
|
mesh = smesh.Mesh( partition )
|
|
cartAlgo = mesh.BodyFitted()
|
|
|
|
# define a cartesian grid using Coordinates
|
|
coords = list(range(-100,100,10))
|
|
cartHyp = cartAlgo.SetGrid( coords,coords,coords, 1000000)
|
|
|
|
# compute the mesh
|
|
mesh.Compute()
|
|
print("nb hexahedra",mesh.NbHexas())
|
|
print("nb tetrahedra",mesh.NbTetras())
|
|
print("nb polyhedra",mesh.NbPolyhedrons())
|
|
print()
|
|
|
|
# define the grid by setting constant spacing
|
|
cartHyp = cartAlgo.SetGrid( "10","10","10", 1000000)
|
|
|
|
mesh.Compute()
|
|
print("nb hexahedra",mesh.NbHexas())
|
|
print("nb tetrahedra",mesh.NbTetras())
|
|
print("nb polyhedra",mesh.NbPolyhedrons())
|
|
print("nb faces",mesh.NbFaces())
|
|
print()
|
|
|
|
# activate creation of faces
|
|
cartHyp.SetToCreateFaces( True )
|
|
|
|
mesh.Compute()
|
|
print("nb hexahedra",mesh.NbHexas())
|
|
print("nb tetrahedra",mesh.NbTetras())
|
|
print("nb polyhedra",mesh.NbPolyhedrons())
|
|
print("nb faces",mesh.NbFaces())
|
|
print()
|
|
|
|
# enable consideration of shared faces
|
|
cartHyp.SetToConsiderInternalFaces( True )
|
|
mesh.Compute()
|
|
print("nb hexahedra",mesh.NbHexas())
|
|
print("nb tetrahedra",mesh.NbTetras())
|
|
print("nb polyhedra",mesh.NbPolyhedrons())
|
|
print("nb faces",mesh.NbFaces())
|
|
print()
|
|
|
|
# define the grid by setting different spacing in 2 sub-ranges of geometry
|
|
spaceFuns = ["5","10+10*t"]
|
|
cartAlgo.SetGrid( [spaceFuns, [0.5]], [spaceFuns, [0.5]], [spaceFuns, [0.25]], 10 )
|
|
|
|
mesh.Compute()
|
|
print("nb hexahedra",mesh.NbHexas())
|
|
print("nb tetrahedra",mesh.NbTetras())
|
|
print("nb polyhedra",mesh.NbPolyhedrons())
|
|
print()
|
|
|
|
# Example of customization of dirtections of the grid axes
|
|
|
|
# make a box with non-orthogonal edges
|
|
xDir = geompy.MakeVectorDXDYDZ( 1.0, 0.1, 0.0, "xDir" )
|
|
yDir = geompy.MakeVectorDXDYDZ(-0.1, 1.0, 0.0, "yDir" )
|
|
zDir = geompy.MakeVectorDXDYDZ( 0.2, 0.3, 1.0, "zDir" )
|
|
face = geompy.MakePrismVecH( xDir, yDir, 1.0 )
|
|
box = geompy.MakePrismVecH( face, zDir, 1.0, theName="box" )
|
|
|
|
spc = "0.1" # spacing
|
|
|
|
# default axes
|
|
mesh = smesh.Mesh( box, "custom axes")
|
|
algo = mesh.BodyFitted()
|
|
algo.SetGrid( spc, spc, spc, 10000 )
|
|
mesh.Compute()
|
|
print("Default axes")
|
|
print(" nb hex:",mesh.NbHexas())
|
|
|
|
# set axes using edges of the box
|
|
algo.SetAxesDirs( xDir, [-0.1,1,0], zDir )
|
|
mesh.Compute()
|
|
print("Manual axes")
|
|
print(" nb hex:",mesh.NbHexas())
|
|
|
|
# set optimal orthogonal axes
|
|
algo.SetOptimalAxesDirs( isOrthogonal=True )
|
|
mesh.Compute()
|
|
print("Optimal orthogonal axes")
|
|
print(" nb hex:",mesh.NbHexas())
|
|
|
|
# set optimal non-orthogonal axes
|
|
algo.SetOptimalAxesDirs( isOrthogonal=False )
|
|
mesh.Compute()
|
|
print("Optimal non-orthogonal axes")
|
|
print(" nb hex:",mesh.NbHexas())
|