1 # Usage of Medial Axis Projection algorithm
3 # for meshing a ring face with quadrangles
6 salome.salome_init_without_session()
7 from salome.geom import geomBuilder
8 geompy = geomBuilder.New()
9 from salome.smesh import smeshBuilder
10 smesh = smeshBuilder.New()
13 circleEdge1 = geompy.MakeCircleR( 3 )
14 circleEdge2 = geompy.MakeCircleR( 7 )
15 ring = geompy.MakeFaceWires( [ circleEdge1, circleEdge2 ], True, theName='Ring' )
16 circleLen1 = geompy.BasicProperties( circleEdge1 )[0]
17 circleLen2 = geompy.BasicProperties( circleEdge2 )[0]
21 mesh = smesh.Mesh( ring )
24 algo1d = mesh.Segment()
25 algo1d.NumberOfSegments( circNbSeg ) # division of circle edges
27 algo2d = mesh.Quadrangle( smeshBuilder.QUAD_MA_PROJ )
28 algo2d.StartEndLength( circleLen2 / circNbSeg, circleLen1 / circNbSeg ) # radial division