def createAnodePlate(self):
#create the plate
rect1 = self.makeRectanglarMesh([0, 0],
self.__widths[1] - self.__widths[0],
self.__heights[1])
rect2 = self.makeRectanglarMesh([0, self.__heights[1]],
self.__widths[1],
self.__heights[0] - self.__heights[1])
#Need some more code:
a1 = rect1.CreateEmptyGroup(SMESH.FACE, '1')
nbAdd = a1.AddFrom(rect1.GetMesh())
a1.SetColor(SALOMEDS.Color(0, 0.666667, 1))
a2 = rect2.CreateEmptyGroup(SMESH.FACE, '1')
nbAdd = a2.AddFrom(rect2.GetMesh())
a2.SetColor(SALOMEDS.Color(0, 0.666667, 1))
smesh.SetName(a1, str(1))
smesh.SetName(a2, str(1))
#create the void
rect3 = self.makeRectanglarMesh(
[self.__widths[1] - self.__widths[0], 0], self.__widths[0],
self.__heights[1])
a3 = rect3.CreateEmptyGroup(SMESH.FACE, '2')
nbAdd = a3.AddFrom(rect3.GetMesh())
a3.SetColor(SALOMEDS.Color(0, 0.666667, 1))
smesh.SetName(a3, str(2))
return [rect1, rect2, rect3]
def RebuildData(theStudy):
smesh.SetCurrentStudy(theStudy)
import StdMeshers
import NETGENPlugin
Mesh_1 = smesh.Mesh(Cut_1)
Regular_1D = Mesh_1.Segment()
Max_Size_1 = Regular_1D.MaxSize(20.2301)
MEFISTO_2D = Mesh_1.Triangle()
Tetrahedron_Netgen = Mesh_1.Tetrahedron(algo=smesh.NETGEN)
isDone = Mesh_1.Compute()
Mesh_1.Clear()
isDone = Mesh_1.Compute()
Mesh_2 = smesh.Mesh(Y)
Regular_1D_1 = Mesh_2.Segment()
ns_1 = Regular_1D_1.NumberOfSegments(50, 10, [])
isDone = Mesh_2.Compute()
Mesh_2.ExtrusionSweepObject1D(
Mesh_2, SMESH.DirStruct(SMESH.PointStruct(0, 0, 100)), 5)
Mesh_2.RotationSweepObject1D(Mesh_2, SMESH.AxisStruct(0, 0, 0, 0, 0, 10),
0.0872665, 6, 1e-05)
## set object names
smesh.SetName(Mesh_1.GetMesh(), 'Mesh_1')
smesh.SetName(Regular_1D.GetAlgorithm(), 'Regular_1D')
smesh.SetName(Max_Size_1, 'Max Size_1')
smesh.SetName(MEFISTO_2D.GetAlgorithm(), 'MEFISTO_2D')
smesh.SetName(Tetrahedron_Netgen.GetAlgorithm(), 'Tetrahedron (Netgen)')
smesh.SetName(Mesh_2.GetMesh(), 'Mesh_2')
smesh.SetName(ns_1, 'ns_1')
if salome.sg.hasDesktop():
salome.sg.updateObjBrowser(0)
### Store presentation parameters of displayed objects
import iparameters
ipar = iparameters.IParameters(
theStudy.GetModuleParameters("Interface Applicative", "SMESH", 1))
#Set up entries:
# set up entry SMESH_4 (Mesh_2) parameters
ipar.setParameter("SMESH_4", "VTKViewer_0_Visibility", "On")
ipar.setParameter("SMESH_4", "VTKViewer_0_Representation", "2")
ipar.setParameter("SMESH_4", "VTKViewer_0_IsShrunk", "0")
ipar.setParameter("SMESH_4", "VTKViewer_0_Entities", "e:1:f:1:v:0")
ipar.setParameter(
"SMESH_4", "VTKViewer_0_Colors",
"surface:0:0.666667:1:backsurface:0:0:1:edge:0:0.666667:1:node:1:0:0")
ipar.setParameter("SMESH_4", "VTKViewer_0_Sizes", "line:1:shrink:0.75")
ipar.setParameter("SMESH_4", "VTKViewer_0_PointMarker", "std:1:9")
ipar.setParameter("SMESH_4", "VTKViewer_0_Opacity", "1")
ipar.setParameter("SMESH_4", "VTKViewer_0_ClippingPlane", "Off")
pass
def createCathodePlate(self):
rightCornerX = sum(
self.__widths) - self.__widths[0] - self.__widths[-1]
#create the plate
print "right corner", rightCornerX
print self.__widths[-2], self.__widths[-1]
#print "op", [rightCornerX[0] ,0], -(self.__widths[1] -self.__widths[0]), self.__heights[0]
rect1 = self.makeRectanglarMesh(
[rightCornerX - (self.__widths[-2] - self.__widths[-1]), 0],
self.__widths[-2] - self.__widths[-1], self.__heights[2])
rect2 = self.makeRectanglarMesh(
[rightCornerX - self.__widths[-2], self.__heights[2]],
self.__widths[-2], self.__heights[0] - self.__heights[2])
a1 = rect1.CreateEmptyGroup(SMESH.FACE, '11')
nbAdd = a1.AddFrom(rect1.GetMesh())
a1.SetColor(SALOMEDS.Color(0, 0.666667, 1))
a2 = rect2.CreateEmptyGroup(SMESH.FACE, '11')
nbAdd = a2.AddFrom(rect2.GetMesh())
a2.SetColor(SALOMEDS.Color(0, 0.666667, 1))
smesh.SetName(a1, str(1))
smesh.SetName(a2, str(1))
#create the void
rect3 = self.makeRectanglarMesh([rightCornerX - self.__widths[-2], 0],
self.__widths[-1], self.__heights[2])
a3 = rect3.CreateEmptyGroup(SMESH.FACE, '10')
nbAdd = a3.AddFrom(rect3.GetMesh())
a3.SetColor(SALOMEDS.Color(0, 0.666667, 1))
smesh.SetName(a3, str(2))
return [rect1, rect2, rect3]
def generate(self):
if self.__widths[1] > 0:
for el in self.createAnodePlate():
self.subMeshes.append(el)
if self.__widths[-2] > 0:
for el in self.createCathodePlate():
self.subMeshes.append(el)
x = self.__widths[1]
materialIDCount = 3
for w in self.__widths[2:-2]:
rect = self.makeRectanglarMesh([x, 0], w, self.__heights[0])
a1 = rect.CreateEmptyGroup(SMESH.FACE, str(materialIDCount))
nbAdd = a1.AddFrom(rect.GetMesh())
a1.SetColor(SALOMEDS.Color(0, 0.666667, 1))
smesh.SetName(a1, str(materialIDCount))
self.subMeshes.append(rect)
x += w
materialIDCount += 1
self.mkCompound()
Arithmetic_1D_2.SetEndLength(1)
Arithmetic_1D_2.SetReversedEdges([])
Arithmetic_1D_2.SetObjectEntry("0:1:1:24")
Arithmetic_1D_1.SetStartLength(1)
Arithmetic_1D_1.SetEndLength(1)
Arithmetic_1D_1.SetReversedEdges([])
Arithmetic_1D_1.SetObjectEntry("0:1:1:24")
isDone = Maillage_1.Compute()
Sous_maillage_1 = Regular_1D.GetSubMesh()
Sous_maillage_2 = Regular_1D_1.GetSubMesh()
Sous_maillage_3 = Regular_1D_2.GetSubMesh()
Sous_maillage_4 = Regular_1D_3.GetSubMesh()
Sous_maillage_5 = Quadrangle_2D.GetSubMesh()
## set object names
smesh.SetName(Maillage_1.GetMesh(), 'Maillage_1')
smesh.SetName(NETGEN_2D.GetAlgorithm(), 'NETGEN_2D')
smesh.SetName(NETGEN_2D_Parameters, 'NETGEN 2D Parameters')
smesh.SetName(Regular_1D.GetAlgorithm(), 'Regular_1D')
smesh.SetName(Nb_Segments_1, 'Nb. Segments_1')
smesh.SetName(Arithmetic_1D_1, 'Arithmetic 1D_1')
smesh.SetName(Arithmetic_1D_2, 'Arithmetic 1D_2')
smesh.SetName(Nb_Segments_2, 'Nb. Segments_2')
smesh.SetName(Quadrangle_2D.GetAlgorithm(), 'Quadrangle_2D')
smesh.SetName(Quadrangle_Parameters_1, 'Quadrangle Parameters_1')
smesh.SetName(Sous_maillage_1, 'Sous-maillage_1')
smesh.SetName(Sous_maillage_2, 'Sous-maillage_2')
smesh.SetName(Sous_maillage_3, 'Sous-maillage_3')
smesh.SetName(Sous_maillage_4, 'Sous-maillage_4')
smesh.SetName(Sous_maillage_5, 'Sous-maillage_5')
Arr.SetColor(SALOMEDS.Color(0, 0.666667, 1))
AxeAv = Maillage_1.CreateEmptyGroup(SMESH.NODE, 'AxeAv')
nbAdd = AxeAv.Add([5])
AxeAv.SetColor(SALOMEDS.Color(0, 0.666667, 1))
AxeAr = Maillage_1.CreateEmptyGroup(SMESH.NODE, 'AxeAr')
nbAdd = AxeAr.Add([9, 8])
AxeAr.SetColor(SALOMEDS.Color(0, 0.666667, 1))
Pedale = Maillage_1.CreateEmptyGroup(SMESH.NODE, 'Pedale')
nbAdd = Pedale.Add([4])
Pedale.SetColor(SALOMEDS.Color(0, 0.666667, 1))
Potence = Maillage_1.CreateEmptyGroup(SMESH.NODE, 'Potence')
nbAdd = Potence.Add([1])
Potence.SetColor(SALOMEDS.Color(0, 0.666667, 1))
isDone = Maillage_1.Compute()
smesh.SetName(Nb_Segments_1, 'Nb. Segments_1')
smesh.SetName(Regular_1D, 'Regular_1D')
smesh.SetName(Maillage_1.GetMesh(), 'Maillage_1')
smesh.SetName(TOUT, 'TOUT')
smesh.SetName(Fourche_1, 'Fourche')
smesh.SetName(Direct, 'Direct')
smesh.SetName(Arr, 'Arr')
smesh.SetName(AxeAv, 'AxeAv')
smesh.SetName(AxeAr, 'AxeAr')
smesh.SetName(Pedale, 'Pedale')
smesh.SetName(Potence, 'Potence')
if salome.sg.hasDesktop():
salome.sg.updateObjBrowser(1)
iparameters.getSession().restoreVisualState(1)
### ### SMESH component ### import smesh, SMESH, SALOMEDS smesh.SetCurrentStudy(theStudy) import StdMeshers Maillage_1 = smesh.Mesh(Treillis) Regular_1D = Maillage_1.Segment() Nb_Segments_1 = Regular_1D.NumberOfSegments(1) Nb_Segments_1.SetDistrType( 0 ) isDone = Maillage_1.Compute() Treillis_1 = Maillage_1.GroupOnGeom(Treillis,'Treillis',SMESH.EDGE) Liaison = Maillage_1.CreateEmptyGroup( SMESH.NODE, 'Liaison' ) nbAdd = Liaison.Add( [ 1, 4 ] ) Efforts = Maillage_1.CreateEmptyGroup( SMESH.NODE, 'Efforts' ) nbAdd = Efforts.Add( [ 2, 3 ] ) isDone = Maillage_1.Compute() ## set object names smesh.SetName(Maillage_1.GetMesh(), 'Maillage_1') smesh.SetName(Regular_1D.GetAlgorithm(), 'Regular_1D') smesh.SetName(Nb_Segments_1, 'Nb. Segments_1') smesh.SetName(Treillis_1, 'Treillis') smesh.SetName(Liaison, 'Liaison') smesh.SetName(Efforts, 'Efforts') if salome.sg.hasDesktop(): salome.sg.updateObjBrowser(1)
Regular_1D = Maillage_1.Segment() Nb_Segments_1 = Regular_1D.NumberOfSegments(20) Nb_Segments_1.SetDistrType( 0 ) isDone = Maillage_1.Compute() POUTRE = Maillage_1.GroupOnGeom(Ligne_1,'Ligne_1',SMESH.EDGE) POUTRE_1 = Maillage_1.GroupOnGeom(Ligne_1,'Ligne_1',SMESH.NODE) POUTRE_1.SetName( 'POU§TRE' ) POUTRE_1.SetName( 'POUTRE' ) POUTRE.SetName( 'POUTRE' ) Encastre = Maillage_1.CreateEmptyGroup( SMESH.NODE, 'Encastre' ) nbAdd = Encastre.Add( [ 1 ] ) Effort = Maillage_1.CreateEmptyGroup( SMESH.NODE, 'Effort' ) nbAdd = Effort.Add( [ 2 ] ) isDone = Maillage_1.Compute() [ POUTRE, POUTRE_1, Encastre, Effort ] = Maillage_1.GetGroups() smesh.SetName(Maillage_1, 'Maillage_1') ## set object names smesh.SetName(Maillage_1.GetMesh(), 'Maillage_1') smesh.SetName(Regular_1D.GetAlgorithm(), 'Regular_1D') smesh.SetName(Nb_Segments_1, 'Nb. Segments_1') smesh.SetName(POUTRE, 'POUTRE') smesh.SetName(POUTRE_1, 'POUTRE') smesh.SetName(Encastre, 'Encastre') smesh.SetName(Effort, 'Effort') ### if salome.sg.hasDesktop(): salome.sg.updateObjBrowser(1)
smesh.SetCurrentStudy(theStudy)
import StdMeshers
Maillage_1 = smesh.Mesh(Assemblage_1)
Regular_1D = Maillage_1.Segment()
Nb_Segments_1 = Regular_1D.NumberOfSegments(1)
Nb_Segments_1.SetDistrType(0)
isDone = Maillage_1.Compute()
TOUT = Maillage_1.CreateEmptyGroup(SMESH.EDGE, 'TOUT')
nbAdd = TOUT.AddFrom(Maillage_1.GetMesh())
Group_1 = Maillage_1.CreateEmptyGroup(SMESH.NODE, 'Group_1')
nbAdd = Group_1.Add([6])
Group_1.SetColor(SALOMEDS.Color(1, 0.666667, 0))
Group_2 = Maillage_1.CreateEmptyGroup(SMESH.NODE, 'Group_2')
nbAdd = Group_2.Add([5])
Group_2.SetColor(SALOMEDS.Color(1, 0.666667, 0))
Group_2 = Maillage_1.CreateEmptyGroup(SMESH.NODE, 'Group_3')
nbAdd = Group_2.Add([1])
Group_2.SetColor(SALOMEDS.Color(1, 0.666667, 0))
## set object names
smesh.SetName(Maillage_1.GetMesh(), 'Maillage_1')
smesh.SetName(Regular_1D.GetAlgorithm(), 'Regular_1D')
smesh.SetName(Nb_Segments_1, 'Nb. Segments_1')
smesh.SetName(TOUT, 'TOUT')
smesh.SetName(Group_1, 'Group_1')
smesh.SetName(Group_2, 'Group_2')
if salome.sg.hasDesktop():
salome.sg.updateObjBrowser(1)
algo3D = tetraN.Tetrahedron(smesh.FULL_NETGEN)
# define hypotheses
n23_params = algo3D.Parameters()
#START MESH PARAM SECTION
n23_params.SetMaxSize(MaxMeshSize / 1000) # here specify in m
n23_params.SetMinSize(MinMeshSize / 1000)
n23_params.SetNbSegPerEdge(MeshSegPerEdge)
n23_params.SetGrowthRate(MeshGrowthRate)
n23_params.SetOptimize(MeshOptimize)
#END MESH PARAM SECTION
# compute the mesh
tetraN.Compute()
smesh.SetName(tetraN.GetMesh(), 'Mesh_1')
print("\n" * 100)
print "creating mesh groups"
# create mesh groups
# Sequence defines body and face index in elmer.
for aGroup in range(0, len(p1_solids)):
tetraN.GroupOnGeom(g[aGroup], 'body{0}'.format(aGroup + 1), SMESH.VOLUME)
for aSolid in range(0, len(p1_solids)):
p1_sufaces = geompy.ExtractShapes(p1_solids[aSolid],
geompy.ShapeType["FACE"], False)
gf = []
for aGroup in range(0, len(p1_sufaces)):
NETGEN_3D_Parameters_1.GroupOnGeom(gf[aGroup],
# set up entry GEOM_1 (Sommet_1) parameters
import smesh, SMESH, SALOMEDS
# Definition du maillage
smesh.SetCurrentStudy(theStudy)
import StdMeshers
Maillage_1 = smesh.Mesh(Anneau)
Regular_1D = Maillage_1.Segment()
Nb_Segments_1 = Regular_1D.NumberOfSegments(15, [], [])
Nb_Segments_1.SetDistrType(0)
isDone = Maillage_1.Compute()
TOUT = Maillage_1.GroupOnGeom(Anneau, 'Anneau', SMESH.EDGE)
Bas_1 = Maillage_1.GroupOnGeom(Bas, 'Bas', SMESH.NODE)
Haut_1 = Maillage_1.GroupOnGeom(Haut, 'Haut', SMESH.NODE)
TOUT.SetName('TOUT')
[TOUT, Bas_1, Haut_1] = Maillage_1.GetGroups()
smesh.SetName(Maillage_1, 'Maillage_1')
Maillage_1.ExportMED(currentPath, 0, SMESH.MED_V2_2, 1)
## set object names
smesh.SetName(Maillage_1.GetMesh(), 'Maillage_1')
smesh.SetName(Regular_1D.GetAlgorithm(), 'Regular_1D')
smesh.SetName(Nb_Segments_1, 'Nb. Segments_1')
smesh.SetName(TOUT, 'TOUT')
smesh.SetName(Bas_1, 'Bas')
smesh.SetName(Haut_1, 'Haut')
if salome.sg.hasDesktop():
salome.sg.updateObjBrowser(1)
iparameters.getSession().restoreVisualState(1)
Mesh_1 = smesh.Mesh(Box_1)
CompositeSegment_1D = Mesh_1.Segment(algo=smesh.COMPOSITE)
Local_Length_1 = CompositeSegment_1D.LocalLength(0.2)
Local_Length_1.SetPrecision( 1e-07 )
Quadrangle_2D = Mesh_1.Quadrangle()
Hexa_3D = smesh.CreateHypothesis('Hexa_3D')
status = Mesh_1.AddHypothesis(Hexa_3D)
forc_1 = Mesh_1.GroupOnGeom(forc,'forc',SMESH.FACE)
toutface_1 = Mesh_1.GroupOnGeom(toutface,'toutface',SMESH.FACE)
tout_1 = Mesh_1.GroupOnGeom(tout,'tout',SMESH.VOLUME)
isDone = Mesh_1.Compute()
Local_Length_1.SetLength( 0.25 )
Local_Length_1.SetPrecision( 1e-07 )
## set object names
smesh.SetName(Mesh_1.GetMesh(), 'Mesh_1')
smesh.SetName(CompositeSegment_1D.GetAlgorithm(), 'CompositeSegment_1D')
smesh.SetName(Local_Length_1, 'Local Length_1')
smesh.SetName(Quadrangle_2D.GetAlgorithm(), 'Quadrangle_2D')
smesh.SetName(Hexa_3D, 'Hexa_3D')
smesh.SetName(forc_1, 'forc')
smesh.SetName(toutface_1, 'toutface')
smesh.SetName(tout_1, 'tout')
### Store presentation parameters of displayed objects
import iparameters
ipar = iparameters.IParameters(theStudy.GetModuleParameters("Interface Applicative", "SMESH", 1))
if salome.sg.hasDesktop():
geompy.addToStudy(Partition_2, 'Partition_2')
geompy.addToStudy(Tuyau_2, 'Tuyau_2')
geompy.addToStudyInFather(Partition_1, Face_5, 'Face_5')
geompy.addToStudyInFather(Partition_1, Face_6, 'Face_6')
###
### SMESH component
###
import smesh, SMESH, SALOMEDS
smesh.SetCurrentStudy(theStudy)
import NETGENPlugin
import StdMeshers
Maillage_1 = smesh.Mesh(joint)
NETGEN_2D = Maillage_1.Triangle(algo=smesh.NETGEN_1D2D)
NETGEN_2D_Simple_Parameters = NETGEN_2D.Parameters(smesh.SIMPLE)
NETGEN_2D_Simple_Parameters.SetNumberOfSegments(15)
NETGEN_2D_Simple_Parameters.LengthFromEdges()
NETGEN_2D_Simple_Parameters.SetAllowQuadrangles(0)
isDone = Maillage_1.Compute()
## set object names
smesh.SetName(Maillage_1.GetMesh(), 'Maillage_1')
smesh.SetName(NETGEN_2D.GetAlgorithm(), 'NETGEN_2D')
smesh.SetName(NETGEN_2D_Simple_Parameters, 'NETGEN 2D Simple Parameters_1')
if salome.sg.hasDesktop():
salome.sg.updateObjBrowser(1)
import smesh, SMESH, SALOMEDS
smesh.SetCurrentStudy(theStudy)
import NETGENPlugin
Mesh_1 = smesh.Mesh(Plaque)
NETGEN_2D = Mesh_1.Triangle(algo=smesh.NETGEN_1D2D)
NETGEN_2D_Parameters = NETGEN_2D.Parameters()
NETGEN_2D_Parameters.SetMaxSize(1)
NETGEN_2D_Parameters.SetMinSize(0.5)
isDone = Mesh_1.Compute()
Encastre_1 = Mesh_1.GroupOnGeom(Encastre, 'Encastre', SMESH.EDGE)
Moment_1 = Mesh_1.GroupOnGeom(Moment, 'Moment', SMESH.EDGE)
Plaque_1 = Mesh_1.GroupOnGeom(Plaque, 'Plaque', SMESH.FACE)
## set object names
smesh.SetName(Mesh_1.GetMesh(), 'Mesh_1')
smesh.SetName(NETGEN_2D.GetAlgorithm(), 'NETGEN_2D')
smesh.SetName(NETGEN_2D_Parameters, 'NETGEN 2D Parameters')
smesh.SetName(Encastre_1, 'Encastre')
smesh.SetName(Moment_1, 'Moment')
smesh.SetName(Plaque_1, 'Plaque')
### Store presentation parameters of displayed objects
import iparameters
ipar = iparameters.IParameters(
theStudy.GetModuleParameters("Interface Applicative", "SMESH", 1))
if salome.sg.hasDesktop():
salome.sg.updateObjBrowser(1)
iparameters.getSession().restoreVisualState(1)
import NETGENPlugin
Maillage_1 = smesh.Mesh(Plaque)
NETGEN_2D = Maillage_1.Triangle(algo=smesh.NETGEN_1D2D)
NETGEN_2D_Parameters = NETGEN_2D.Parameters()
NETGEN_2D_Parameters.SetMaxSize(0.1)
NETGEN_2D_Parameters.SetSecondOrder(0)
NETGEN_2D_Parameters.SetOptimize(1)
NETGEN_2D_Parameters.SetFineness(2)
NETGEN_2D_Parameters.SetMinSize(0.002)
NETGEN_2D_Parameters.SetQuadAllowed(0)
isDone = Maillage_1.Compute()
BORD_1 = Maillage_1.GroupOnGeom(BORD, 'BORD', SMESH.EDGE)
BORD_2 = Maillage_1.GroupOnGeom(BORD, 'BORD', SMESH.NODE)
Plaque_1 = Maillage_1.GroupOnGeom(Plaque, 'Plaque', SMESH.FACE)
[BORD_1, BORD_2, Plaque_1] = Maillage_1.GetGroups()
smesh.SetName(Maillage_1, 'Maillage_1')
Maillage_1.ExportMED(
r'/home/bordreuil/Enseignement/elementsFiNimes/Examples/Plaques/E_Modes/mode.mmed',
0, SMESH.MED_V2_2, 1)
## set object names
smesh.SetName(Maillage_1.GetMesh(), 'Maillage_1')
smesh.SetName(NETGEN_2D.GetAlgorithm(), 'NETGEN_2D')
smesh.SetName(NETGEN_2D_Parameters, 'NETGEN 2D Parameters')
smesh.SetName(BORD_1, 'BORD')
smesh.SetName(BORD_2, 'BORD')
smesh.SetName(Plaque_1, 'Plaque')
if salome.sg.hasDesktop():
salome.sg.updateObjBrowser(1)
Compound_Mesh_1 = smesh.Concatenate([
Mesh_1.GetMesh(),
Mesh_2.GetMesh(),
Mesh_3.GetMesh(),
Mesh_4.GetMesh(),
Mesh_5.GetMesh(),
Mesh_6.GetMesh(),
Mesh_7.GetMesh()
], 1, 1, 1e-05)
[a1_2, a1_3, a2_2, a2_3, a3_2, a4_2, a5_1, a6_1, a7_1, a3_3,
a4_3] = Compound_Mesh_1.GetGroups()
# isDone = Compound_Mesh_1.RemoveElements( [ 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124 ] )
## set object names
smesh.SetName(Mesh_1.GetMesh(), 'Mesh_1')
smesh.SetName(Regular_1D.GetAlgorithm(), 'Regular_1D')
smesh.SetName(Local_Length_1, 'Local Length_1')
smesh.SetName(Quadrangle_2D.GetAlgorithm(), 'Quadrangle_2D')
smesh.SetName(Mesh_2.GetMesh(), 'Mesh_2')
smesh.SetName(Mesh_3.GetMesh(), 'Mesh_3')
smesh.SetName(Mesh_4.GetMesh(), 'Mesh_4')
smesh.SetName(Mesh_5.GetMesh(), 'Mesh_5')
smesh.SetName(Mesh_6.GetMesh(), 'Mesh_6')
smesh.SetName(Mesh_7.GetMesh(), 'Mesh_7')
smesh.SetName(a1, '1')
smesh.SetName(a2, '2')
smesh.SetName(a3, '3')
smesh.SetName(a4, '4')
smesh.SetName(a5, '5')
smesh.SetName(a6, '6')
