2018-05-25 22:04:48 +05:00
|
|
|
# Usage of Radial Prism 3D meshing algorithm
|
|
|
|
|
|
|
|
|
|
import salome
|
2021-08-12 14:38:10 +05:00
|
|
|
salome.salome_init_without_session()
|
2018-05-25 22:04:48 +05:00
|
|
|
from salome.geom import geomBuilder
|
2018-06-29 18:38:39 +05:00
|
|
|
geompy = geomBuilder.New()
|
2018-05-25 22:04:48 +05:00
|
|
|
import SMESH
|
|
|
|
|
from salome.smesh import smeshBuilder
|
2018-06-29 18:38:39 +05:00
|
|
|
smesh = smeshBuilder.New()
|
2018-05-25 22:04:48 +05:00
|
|
|
|
|
|
|
|
# Create geometry: hollow sphere
|
|
|
|
|
|
|
|
|
|
sphere_1 = geompy.MakeSphereR( 100 )
|
|
|
|
|
sphere_2 = geompy.MakeSphereR( 50 )
|
|
|
|
|
|
|
|
|
|
hollow_sphere = geompy.MakeCut( sphere_1, sphere_2, theName="hollow sphere")
|
|
|
|
|
|
|
|
|
|
faces = geompy.ExtractShapes( hollow_sphere, geompy.ShapeType["FACE"] )
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
# Create mesh
|
|
|
|
|
|
|
|
|
|
mesh = smesh.Mesh( hollow_sphere, "Mesh of hollow sphere" )
|
|
|
|
|
|
|
|
|
|
# assign Global Radial Prism algorithm
|
|
|
|
|
prism_algo = mesh.Prism()
|
|
|
|
|
|
|
|
|
|
# define projection between the inner and outer spheres
|
|
|
|
|
mesh.Triangle( smeshBuilder.NETGEN_1D2D, faces[0] ) # NETGEN on faces[0]
|
|
|
|
|
mesh.Projection1D2D( faces[1] ).SourceFace( faces[0] ) # projection faces[0] -> faces[1]
|
|
|
|
|
|
|
|
|
|
# define distribution of layers using Number of Segments hypothesis in logarithmic mode
|
|
|
|
|
prism_algo.NumberOfSegments( 4, 5. )
|
|
|
|
|
|
|
|
|
|
# compute the mesh
|
|
|
|
|
mesh.Compute()
|
