1 # Usage of Medial Axis Projection algorithm
3 # for meshing a ring face with quadrangles
6 salome.salome_init_without_session()
8 from salome.geom import geomBuilder
9 from salome.smesh import smeshBuilder
11 geom_builder = geomBuilder.New()
12 smesh_builder = smeshBuilder.New()
15 circleEdge1 = geom_builder.MakeCircleR( 3 )
16 circleEdge2 = geom_builder.MakeCircleR( 7 )
17 ring = geom_builder.MakeFaceWires( [ circleEdge1, circleEdge2 ], True, theName='Ring' )
18 circleLen1 = geom_builder.BasicProperties( circleEdge1 )[0]
19 circleLen2 = geom_builder.BasicProperties( circleEdge2 )[0]
23 mesh = smesh_builder.Mesh( ring )
26 algo1d = mesh.Segment()
27 algo1d.NumberOfSegments( circNbSeg ) # division of circle edges
29 algo2d = mesh.Quadrangle( smeshBuilder.QUAD_MA_PROJ )
30 algo2d.StartEndLength( circleLen2 / circNbSeg, circleLen1 / circNbSeg ) # radial division
32 if not mesh.Compute(): raise Exception("Error when computing Mesh")