1 # Viscous layers construction
6 from salome.geom import geomBuilder
7 geompy = geomBuilder.New(salome.myStudy)
10 from salome.smesh import smeshBuilder
11 smesh = smeshBuilder.New(salome.myStudy)
13 X = geompy.MakeVectorDXDYDZ( 1,0,0 )
14 O = geompy.MakeVertex( 100,50,50 )
15 plane = geompy.MakePlane( O, X, 200 ) # plane YZ
17 box = geompy.MakeBoxDXDYDZ(200,100,100)
19 shape = geompy.MakeHalfPartition( box, plane )
21 faces = geompy.SubShapeAllSorted(shape, geompy.ShapeType["FACE"])
23 ignoreFaces = [ faces[0], faces[-1]]
25 geompy.addToStudy( shape, "shape" )
26 geompy.addToStudyInFather( shape, face1, "face1")
29 mesh = smesh.Mesh(shape, "CFD")
31 mesh.Segment().NumberOfSegments( 4 )
34 mesh.Quadrangle(face1)
36 algo3D = mesh.Tetrahedron()
41 layersHyp = algo3D.ViscousLayers(thickness,numberOfLayers,stretchFactor,ignoreFaces)
45 mesh.MakeGroup("Tetras",SMESH.VOLUME,SMESH.FT_ElemGeomType,"=",SMESH.Geom_TETRA)
46 mesh.MakeGroup("Pyras",SMESH.VOLUME,SMESH.FT_ElemGeomType,"=",SMESH.Geom_PYRAMID)
47 mesh.MakeGroup("Prims",SMESH.VOLUME,SMESH.FT_ElemGeomType,"=",SMESH.Geom_PENTA)