smesh/test/SMESH_Partition1_tetra.py

187 lines
5.8 KiB
Python
Raw Permalink Normal View History

2012-08-09 16:03:55 +06:00
# -*- coding: iso-8859-1 -*-
2024-01-20 19:16:41 +05:00
# Copyright (C) 2007-2024 CEA, EDF, OPEN CASCADE
2009-02-17 10:27:49 +05:00
#
2012-08-09 16:03:55 +06:00
# Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
# CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
#
2012-08-09 16:03:55 +06:00
# This library is free software; you can redistribute it and/or
# modify it under the terms of the GNU Lesser General Public
# License as published by the Free Software Foundation; either
2014-02-20 18:25:37 +06:00
# version 2.1 of the License, or (at your option) any later version.
#
2012-08-09 16:03:55 +06:00
# This library is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
# Lesser General Public License for more details.
#
2012-08-09 16:03:55 +06:00
# You should have received a copy of the GNU Lesser General Public
# License along with this library; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
#
2012-08-09 16:03:55 +06:00
# See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
2004-06-18 14:34:31 +06:00
#
2012-08-09 16:03:55 +06:00
2004-06-18 14:34:31 +06:00
# Tetrahedrization of the geometry generated by the Python script GEOM_Partition1.py
# Hypothesis and algorithms for the mesh generation are global
# -- Rayon de la bariere
2009-02-17 10:27:49 +05:00
#
2004-06-18 14:34:31 +06:00
import salome
salome.salome_init()
import GEOM
from salome.geom import geomBuilder
2017-06-13 15:01:10 +05:00
geompy = geomBuilder.New()
import SMESH, SALOMEDS
from salome.smesh import smeshBuilder
2017-06-13 15:01:10 +05:00
smesh = smeshBuilder.New()
2004-12-01 15:48:31 +05:00
from math import sqrt
2004-06-18 14:34:31 +06:00
#---------------------------------------------------------------
barier_height = 7.0
barier_radius = 5.6 / 2 # Rayon de la bariere
colis_radius = 1.0 / 2 # Rayon du colis
2017-06-13 15:01:10 +05:00
colis_step = 2.0 # Distance s<>parant deux colis
2004-06-18 14:34:31 +06:00
cc_width = 0.11 # Epaisseur du complement de colisage
# --
cc_radius = colis_radius + cc_width
colis_center = sqrt(2.0)*colis_step/2
# --
boolean_common = 1
boolean_cut = 2
boolean_fuse = 3
boolean_section = 4
# --
2004-12-01 15:48:31 +05:00
p0 = geompy.MakeVertex(0.,0.,0.)
vz = geompy.MakeVectorDXDYDZ(0.,0.,1.)
2004-06-18 14:34:31 +06:00
# --
2004-12-01 15:48:31 +05:00
barier = geompy.MakeCylinder(p0, vz, barier_radius, barier_height)
2004-06-18 14:34:31 +06:00
2004-12-01 15:48:31 +05:00
# --
2004-06-18 14:34:31 +06:00
2004-12-01 15:48:31 +05:00
colis = geompy.MakeCylinder(p0, vz, colis_radius, barier_height)
cc = geompy.MakeCylinder(p0, vz, cc_radius, barier_height)
2004-06-18 14:34:31 +06:00
2004-12-01 15:48:31 +05:00
colis_cc = geompy.MakeCompound([colis, cc])
colis_cc = geompy.MakeTranslation(colis_cc, colis_center, 0.0, 0.0)
2004-06-18 14:34:31 +06:00
2004-12-01 15:48:31 +05:00
colis_cc_multi = geompy.MultiRotate1D(colis_cc, vz, 4)
2004-06-18 14:34:31 +06:00
# --
Compound1 = geompy.MakeCompound([colis_cc_multi, barier])
SubShape_theShape = geompy.SubShapeAll(Compound1,geompy.ShapeType["SOLID"])
alveole = geompy.MakePartition(SubShape_theShape)
2004-06-18 14:34:31 +06:00
2017-03-20 17:27:30 +05:00
print("Analysis of the geometry to mesh (right after the Partition) :")
2004-06-18 14:34:31 +06:00
2004-12-01 15:48:31 +05:00
subShellList = geompy.SubShapeAll(alveole, geompy.ShapeType["SHELL"])
subFaceList = geompy.SubShapeAll(alveole, geompy.ShapeType["FACE"])
subEdgeList = geompy.SubShapeAll(alveole, geompy.ShapeType["EDGE"])
2004-06-18 14:34:31 +06:00
2017-03-20 17:27:30 +05:00
print("number of Shells in alveole : ", len(subShellList))
print("number of Faces in alveole : ", len(subFaceList))
print("number of Edges in alveole : ", len(subEdgeList))
2004-06-18 14:34:31 +06:00
2004-12-01 15:48:31 +05:00
subshapes = geompy.SubShapeAll(alveole, geompy.ShapeType["SHAPE"])
2004-06-18 14:34:31 +06:00
2012-08-09 16:03:55 +06:00
## there are 9 sub-shapes
2004-06-18 14:34:31 +06:00
2004-12-01 15:48:31 +05:00
comp1 = geompy.MakeCompound( [ subshapes[0], subshapes[1] ] )
comp2 = geompy.MakeCompound( [ subshapes[2], subshapes[3] ] )
comp3 = geompy.MakeCompound( [ subshapes[4], subshapes[5] ] )
comp4 = geompy.MakeCompound( [ subshapes[6], subshapes[7] ] )
2004-06-18 14:34:31 +06:00
2004-12-01 15:48:31 +05:00
compGOs = []
compGOs.append( comp1 )
compGOs.append( comp2 )
compGOs.append( comp3 )
compGOs.append( comp4 )
comp = geompy.MakeCompound( compGOs )
2004-06-18 14:34:31 +06:00
2004-12-01 15:48:31 +05:00
alveole = geompy.MakeCompound( [ comp, subshapes[8] ])
idalveole = geompy.addToStudy(alveole, "alveole")
2004-06-18 14:34:31 +06:00
2017-03-20 17:27:30 +05:00
print("Analysis of the geometry to mesh (right after the MakeCompound) :")
2004-06-18 14:34:31 +06:00
2004-12-01 15:48:31 +05:00
subShellList = geompy.SubShapeAll(alveole, geompy.ShapeType["SHELL"])
subFaceList = geompy.SubShapeAll(alveole, geompy.ShapeType["FACE"])
subEdgeList = geompy.SubShapeAll(alveole, geompy.ShapeType["EDGE"])
2004-06-18 14:34:31 +06:00
2017-03-20 17:27:30 +05:00
print("number of Shells in alveole : ", len(subShellList))
print("number of Faces in alveole : ", len(subFaceList))
print("number of Edges in alveole : ", len(subEdgeList))
2004-06-18 14:34:31 +06:00
2004-12-01 15:48:31 +05:00
status = geompy.CheckShape(alveole)
2017-03-20 17:27:30 +05:00
print(" check status ", status)
2004-06-18 14:34:31 +06:00
# ---- init a Mesh with the alveole
shape_mesh = salome.IDToObject( idalveole )
mesh = smesh.Mesh(shape_mesh, "MeshAlveole")
2017-03-20 17:27:30 +05:00
print("-------------------------- create Hypothesis (In this case global hypothesis are used)")
2004-06-18 14:34:31 +06:00
2017-03-20 17:27:30 +05:00
print("-------------------------- NumberOfSegments")
2004-06-18 14:34:31 +06:00
numberOfSegments = 10
regular1D = mesh.Segment()
hypNbSeg = regular1D.NumberOfSegments(numberOfSegments)
2017-03-20 17:27:30 +05:00
print(hypNbSeg.GetName())
print(hypNbSeg.GetId())
print(hypNbSeg.GetNumberOfSegments())
smesh.SetName(hypNbSeg, "NumberOfSegments_" + str(numberOfSegments))
2004-06-18 14:34:31 +06:00
2017-03-20 17:27:30 +05:00
print("-------------------------- MaxElementArea")
2004-06-18 14:34:31 +06:00
maxElementArea = 0.1
2022-04-15 20:49:22 +05:00
triangle = mesh.Triangle()
hypArea = triangle.MaxElementArea(maxElementArea)
2017-03-20 17:27:30 +05:00
print(hypArea.GetName())
print(hypArea.GetId())
print(hypArea.GetMaxElementArea())
smesh.SetName(hypArea, "MaxElementArea_" + str(maxElementArea))
2004-06-18 14:34:31 +06:00
2017-03-20 17:27:30 +05:00
print("-------------------------- MaxElementVolume")
2004-06-18 14:34:31 +06:00
maxElementVolume = 0.5
netgen3D = mesh.Tetrahedron(smeshBuilder.NETGEN)
hypVolume = netgen3D.MaxElementVolume(maxElementVolume)
2017-03-20 17:27:30 +05:00
print(hypVolume.GetName())
print(hypVolume.GetId())
print(hypVolume.GetMaxElementVolume())
smesh.SetName(hypVolume, "MaxElementVolume_" + str(maxElementVolume))
2004-06-18 14:34:31 +06:00
2017-03-20 17:27:30 +05:00
print("-------------------------- compute the mesh of alveole ")
ret = mesh.Compute()
if not ret:
raise Exception("Error when computing Mesh")
log = mesh.GetLog(0) # no erase trace
# for linelog in log:
# print(linelog)
print("Information about the Mesh_mechanic:")
print("Number of nodes : ", mesh.NbNodes())
print("Number of edges : ", mesh.NbEdges())
print("Number of faces : ", mesh.NbFaces())
print("Number of triangles : ", mesh.NbTriangles())
print("Number of volumes : ", mesh.NbVolumes())
print("Number of tetrahedrons: ", mesh.NbTetras())
2004-12-01 15:48:31 +05:00
2017-06-13 15:01:10 +05:00
salome.sg.updateObjBrowser()