2021-02-11 13:36:00 +05:00
|
|
|
import salome
|
|
|
|
|
salome.salome_init()
|
|
|
|
|
|
|
|
|
|
import GEOM
|
|
|
|
|
from salome.geom import geomBuilder
|
|
|
|
|
geompy = geomBuilder.New()
|
|
|
|
|
|
|
|
|
|
import SMESH, SALOMEDS
|
|
|
|
|
from salome.smesh import smeshBuilder
|
|
|
|
|
smesh = smeshBuilder.New()
|
|
|
|
|
|
|
|
|
|
import math
|
|
|
|
|
|
|
|
|
|
axes = [
|
|
|
|
|
geompy.MakeVectorDXDYDZ(1, 0, 0),
|
|
|
|
|
geompy.MakeVectorDXDYDZ(0, 1, 0),
|
|
|
|
|
geompy.MakeVectorDXDYDZ(0, 0, 1)
|
|
|
|
|
]
|
|
|
|
|
|
|
|
|
|
vtx = [
|
|
|
|
|
geompy.MakeVertex(1 / math.sqrt(2), -1 / math.sqrt(2), 0.5),
|
|
|
|
|
geompy.MakeVertex(-1 / math.sqrt(2), 1 / math.sqrt(2), -0.5),
|
|
|
|
|
geompy.MakeVertex(-1, -1, -0.5),
|
|
|
|
|
geompy.MakeVertex(-0.5, -0.5, 0)
|
|
|
|
|
]
|
|
|
|
|
|
|
|
|
|
box = geompy.MakeBoxTwoPnt(vtx[0], vtx[1])
|
|
|
|
|
box = geompy.MakeRotation(box, axes[2], 45 * math.pi / 180.0)
|
|
|
|
|
|
|
|
|
|
#alpha = [ x * 0.01 for x in range(1, 28 + 1) ]
|
|
|
|
|
alpha=[0.15]
|
|
|
|
|
|
|
|
|
|
#for coef in alpha:
|
|
|
|
|
spheres = geompy.MakeSpherePntR(vtx[2], math.sqrt(3) / 4 / (1 - alpha[0]))
|
|
|
|
|
spheres = geompy.MakeMultiTranslation2D(spheres, axes[0], 1, 3, axes[1], 1, 3)
|
|
|
|
|
spheres = geompy.MakeMultiTranslation1D(spheres, axes[2], 1, 2)
|
|
|
|
|
|
|
|
|
|
spheres2 = geompy.MakeSpherePntR(vtx[3], math.sqrt(3) / 4 / (1 - alpha[0]))
|
|
|
|
|
spheres2 = geompy.MakeMultiTranslation2D(spheres2, axes[0], 1, 3, axes[1], 1, 3)
|
|
|
|
|
|
2021-02-11 15:48:17 +05:00
|
|
|
PoreBC = geompy.MakeCutList(box, [spheres, spheres2], False)
|
2021-02-11 13:36:00 +05:00
|
|
|
|
2021-02-11 15:48:17 +05:00
|
|
|
geompy.addToStudy(PoreBC, 'PoreBC')
|
2021-02-11 13:36:00 +05:00
|
|
|
|
|
|
|
|
box2 = geompy.MakeBoxTwoPnt(vtx[0], geompy.MakeVertex(0, 0, 0))
|
|
|
|
|
box2 = geompy.MakeRotation(box2, axes[2], 45 * math.pi / 180.0)
|
2021-02-11 15:48:17 +05:00
|
|
|
Segment1_8 = geompy.MakeCommonList([PoreBC, box2], True)
|
2021-02-11 13:36:00 +05:00
|
|
|
|
|
|
|
|
geompy.addToStudy(Segment1_8, 'Segment1_8')
|
|
|
|
|
|
2021-02-11 15:48:17 +05:00
|
|
|
mesh = smesh.Mesh(PoreBC)
|
|
|
|
|
netgen = mesh.Tetrahedron(algo=smeshBuilder.NETGEN_1D2D3D)
|
|
|
|
|
|
|
|
|
|
param = netgen.Parameters()
|
|
|
|
|
param.SetSecondOrder( 0 )
|
|
|
|
|
param.SetOptimize( 1 )
|
|
|
|
|
param.SetChordalError( -1 )
|
|
|
|
|
param.SetChordalErrorEnabled( 0 )
|
|
|
|
|
param.SetUseSurfaceCurvature( 1 )
|
|
|
|
|
param.SetFuseEdges( 1 )
|
|
|
|
|
param.SetCheckChartBoundary( 0 )
|
|
|
|
|
param.SetMinSize( 0.01 )
|
|
|
|
|
param.SetMaxSize( 0.02 )
|
|
|
|
|
param.SetFineness( 5 )
|
|
|
|
|
param.SetGrowthRate( 0.1 )
|
|
|
|
|
param.SetNbSegPerEdge( 5 )
|
|
|
|
|
param.SetNbSegPerRadius( 10 )
|
|
|
|
|
param.SetQuadAllowed( 1 )
|
|
|
|
|
|
|
|
|
|
#vlayer = netgen.ViscousLayers(0.05, 3, 1.5, [15, 29, 54], 1, smeshBuilder.SURF_OFFSET_SMOOTH)
|
|
|
|
|
|
|
|
|
|
isDone = mesh.Compute()
|
|
|
|
|
|
|
|
|
|
if not isDone:
|
|
|
|
|
print("Mesh is not computed")
|
|
|
|
|
|
|
|
|
|
|
2021-02-11 13:36:00 +05:00
|
|
|
if salome.sg.hasDesktop():
|
|
|
|
|
salome.sg.updateObjBrowser()
|
|
|
|
|
|
