1 # This is an automation of the cylinder-box object, defined with the coordinates of its center, its radius, and the box's
3 # The pitch ratio is calculated automatically from the minimum of the box dimensions on x and y.
4 # This functions can take a groups input containing the group names of 4 sides in addition to the internal circular boundary
5 # in the following order : [South,North,West,East,Internal].
7 import sys, salome, geompy, smesh, SMESH, math, commands
8 CWD = commands.getoutput('pwd')
12 from MacObject import *
13 import Config, GenFunctions
15 def Cylinder (X0 , Y0 , D , DX , DY , LocalMeshing , **args) :
16 if args.__contains__('DLocal') : DLocal = float(args['DLocal'])
17 else : DLocal = float(min(DX,DY))
19 # K is the pitch ratio
20 K = float(D)/(DLocal-D)
21 print "A local pitch ratio of K =", K ," will be used. "
22 NumCuts = 2*GenFunctions.QuarCylParam(K)
23 InternalMeshing = int(math.ceil(math.pi*D/(4*NumCuts*LocalMeshing)))
24 if InternalMeshing == 0 : InternalMeshing = 1 # This sets a minimum meshing condition in order to avoid an error. The user is notified of the value considered for the local meshing
25 print "Possible Local meshing is :", math.pi*D/(4*NumCuts*InternalMeshing), "\nThis value is returned by this function for your convenience.\n"
26 if args.__contains__('groups') :
27 GroupNames = args['groups']
28 else : GroupNames = [None, None, None, None, None]
32 GN1 = [None,GroupNames[1],None,GroupNames[3],GroupNames[4]]
33 GN2 = [None,GroupNames[1],GroupNames[2],None,GroupNames[4]]
34 GN3 = [GroupNames[0],None,GroupNames[2],None,GroupNames[4]]
35 GN4 = [GroupNames[0],None,None,GroupNames[3],GroupNames[4]]
37 GN1 = [None,GroupNames[1],None,None,GroupNames[4]]
38 GN2 = [None,GroupNames[1],None,None,GroupNames[4]]
39 GN3 = [GroupNames[0],None,None,None,GroupNames[4]]
40 GN4 = [GroupNames[0],None,None,None,GroupNames[4]]
42 GN5 = [GroupNames[0],GroupNames[1],None,GroupNames[3]]
43 GN6 = [GroupNames[0],GroupNames[1],GroupNames[2],None]
46 GN1 = [None,None,None,GroupNames[3],GroupNames[4]]
47 GN2 = [None,None,GroupNames[2],None,GroupNames[4]]
48 GN3 = [None,None,GroupNames[2],None,GroupNames[4]]
49 GN4 = [None,None,None,GroupNames[3],GroupNames[4]]
50 GN7 = [GroupNames[0],None,GroupNames[2],GroupNames[3]]
51 GN8 = [None,GroupNames[1],GroupNames[2],GroupNames[3]]
53 GN1 = [None,None,None,None,GroupNames[4]]
54 GN2 = [None,None,None,None,GroupNames[4]]
55 GN3 = [None,None,None,None,GroupNames[4]]
56 GN4 = [None,None,None,None,GroupNames[4]]
58 GN5 = [None,None,None,GroupNames[3]]
59 GN6 = [None,None,GroupNames[2],None]
61 GN9 = [GroupNames[0],None,None,GroupNames[3]]
62 GN10 = [GroupNames[0],None,None,None]
63 GN11 = [GroupNames[0],None,GroupNames[2],None]
65 GN12 = [None,GroupNames[1],None,GroupNames[3]]
66 GN13 = [None,GroupNames[1],None,None]
67 GN14 = [None,GroupNames[1],GroupNames[2],None]
71 Obj.append(MacObject('QuartCyl',[(X0+DLocal/4.,Y0+DLocal/4.),(DLocal/2.,DLocal/2.)],[InternalMeshing,'NE',K], groups = GN1))
72 Obj.append(MacObject('QuartCyl',[(X0-DLocal/4.,Y0+DLocal/4.),(DLocal/2.,DLocal/2.)],['auto','NW',K], groups = GN2))
73 Obj.append(MacObject('QuartCyl',[(X0-DLocal/4.,Y0-DLocal/4.),(DLocal/2.,DLocal/2.)],['auto','SW',K], groups = GN3))
74 Obj.append(MacObject('QuartCyl',[(X0+DLocal/4.,Y0-DLocal/4.),(DLocal/2.,DLocal/2.)],['auto','SE',K], groups = GN4))
78 Obj.append(MacObject('CompBoxF',[(X0+DLocal/2.+dX/2.,Y0),(dX,DLocal)],['auto'], groups = GN5))
79 Obj.append(MacObject('CompBoxF',[(X0-DLocal/2.-dX/2.,Y0),(dX,DLocal)],['auto'], groups = GN6))
84 Obj.append(MacObject('CompBoxF',[(X0+DLocal/2.+dX/2.,Y0-DLocal/2.-dY/2.),(dX,dY)],['auto'], groups = GN9))
85 Obj.append(MacObject('CompBoxF',[(X0,Y0-DLocal/2.-dY/2.),(DLocal,dY)],['auto'], groups = GN10))
86 Obj.append(MacObject('CompBoxF',[(X0-DLocal/2.-dX/2.,Y0-DLocal/2.-dY/2.),(dX,dY)],['auto'], groups = GN11))
87 Obj.append(MacObject('CompBoxF',[(X0+DLocal/2.+dX/2.,Y0+DLocal/2.+dY/2.),(dX,dY)],['auto'], groups = GN12))
88 Obj.append(MacObject('CompBoxF',[(X0,Y0+DLocal/2.+dY/2.),(DLocal,dY)],['auto'], groups = GN13))
89 Obj.append(MacObject('CompBoxF',[(X0-DLocal/2.-dX/2.,Y0+DLocal/2.+dY/2.),(dX,dY)],['auto'], groups = GN14))
91 Obj.append(MacObject('CompBoxF',[(X0,Y0-DLocal/2.-dY/2.),(DLocal,dY)],['auto'], groups = GN7))
92 Obj.append(MacObject('CompBoxF',[(X0,Y0+DLocal/2.+dY/2.),(DLocal,dY)],['auto'], groups = GN8))