1 # CEA/LGLS 2004-2005, Francis KLOSS (OCC)
2 # =======================================
11 # A small cube centered and put on a great cube build by points, edges, faces and solids
16 greatPoint111 = MakeVertex( 0, 0, 0)
17 greatPoint211 = MakeVertex(10, 0, 0)
18 greatPoint311 = MakeVertex(20, 0, 0)
19 greatPoint411 = MakeVertex(30, 0, 0)
21 greatPoint121 = MakeVertex( 0, 10, 0)
22 greatPoint221 = MakeVertex(10, 10, 0)
23 greatPoint321 = MakeVertex(20, 10, 0)
24 greatPoint421 = MakeVertex(30, 10, 0)
26 greatPoint112 = MakeVertex( 0, 0, 10)
27 greatPoint212 = MakeVertex(10, 0, 10)
28 greatPoint312 = MakeVertex(20, 0, 10)
29 greatPoint412 = MakeVertex(30, 0, 10)
31 greatPoint122 = MakeVertex( 0, 10, 10)
32 greatPoint222 = MakeVertex(10, 10, 10)
33 greatPoint322 = MakeVertex(20, 10, 10)
34 greatPoint422 = MakeVertex(30, 10, 10)
36 greatPoint113 = MakeVertex( 0, 0, 20)
37 greatPoint213 = MakeVertex(10, 0, 20)
38 greatPoint313 = MakeVertex(20, 0, 20)
39 greatPoint413 = MakeVertex(30, 0, 20)
41 greatPoint123 = MakeVertex( 0, 10, 20)
42 greatPoint223 = MakeVertex(10, 10, 20)
43 greatPoint323 = MakeVertex(20, 10, 20)
44 greatPoint423 = MakeVertex(30, 10, 20)
46 greatPoint114 = MakeVertex( 0, 0, 30)
47 greatPoint214 = MakeVertex(10, 0, 30)
48 greatPoint314 = MakeVertex(20, 0, 30)
49 greatPoint414 = MakeVertex(30, 0, 30)
51 greatPoint124 = MakeVertex( 0, 10, 30)
52 greatPoint224 = MakeVertex(10, 10, 30)
53 greatPoint324 = MakeVertex(20, 10, 30)
54 greatPoint424 = MakeVertex(30, 10, 30)
57 smallPoint111 = greatPoint222
58 smallPoint211 = greatPoint322
59 smallPoint121 = MakeVertex(10, 20, 10)
60 smallPoint221 = MakeVertex(20, 20, 10)
62 smallPoint112 = greatPoint223
63 smallPoint212 = greatPoint323
64 smallPoint122 = MakeVertex(10, 20, 20)
65 smallPoint222 = MakeVertex(20, 20, 20)
70 smallEdgeX11 = MakeEdge(smallPoint111, smallPoint211)
71 smallEdgeX21 = MakeEdge(smallPoint121, smallPoint221)
72 smallEdgeX12 = MakeEdge(smallPoint112, smallPoint212)
73 smallEdgeX22 = MakeEdge(smallPoint122, smallPoint222)
75 smallEdgeY11 = MakeEdge(smallPoint111, smallPoint121)
76 smallEdgeY21 = MakeEdge(smallPoint211, smallPoint221)
77 smallEdgeY12 = MakeEdge(smallPoint112, smallPoint122)
78 smallEdgeY22 = MakeEdge(smallPoint212, smallPoint222)
80 smallEdgeZ11 = MakeEdge(smallPoint111, smallPoint112)
81 smallEdgeZ21 = MakeEdge(smallPoint211, smallPoint212)
82 smallEdgeZ12 = MakeEdge(smallPoint121, smallPoint122)
83 smallEdgeZ22 = MakeEdge(smallPoint221, smallPoint222)
86 greatEdgeX111 = MakeEdge(greatPoint111, greatPoint211)
87 greatEdgeX211 = MakeEdge(greatPoint211, greatPoint311)
88 greatEdgeX311 = MakeEdge(greatPoint311, greatPoint411)
89 greatEdgeX121 = MakeEdge(greatPoint121, greatPoint221)
90 greatEdgeX221 = MakeEdge(greatPoint221, greatPoint321)
91 greatEdgeX321 = MakeEdge(greatPoint321, greatPoint421)
93 greatEdgeX112 = MakeEdge(greatPoint112, greatPoint212)
94 greatEdgeX212 = MakeEdge(greatPoint212, greatPoint312)
95 greatEdgeX312 = MakeEdge(greatPoint312, greatPoint412)
96 greatEdgeX122 = MakeEdge(greatPoint122, greatPoint222)
97 greatEdgeX222 = smallEdgeX11
98 greatEdgeX322 = MakeEdge(greatPoint322, greatPoint422)
100 greatEdgeX113 = MakeEdge(greatPoint113, greatPoint213)
101 greatEdgeX213 = MakeEdge(greatPoint213, greatPoint313)
102 greatEdgeX313 = MakeEdge(greatPoint313, greatPoint413)
103 greatEdgeX123 = MakeEdge(greatPoint123, greatPoint223)
104 greatEdgeX223 = smallEdgeX12
105 greatEdgeX323 = MakeEdge(greatPoint323, greatPoint423)
107 greatEdgeX114 = MakeEdge(greatPoint114, greatPoint214)
108 greatEdgeX214 = MakeEdge(greatPoint214, greatPoint314)
109 greatEdgeX314 = MakeEdge(greatPoint314, greatPoint414)
110 greatEdgeX124 = MakeEdge(greatPoint124, greatPoint224)
111 greatEdgeX224 = MakeEdge(greatPoint224, greatPoint324)
112 greatEdgeX324 = MakeEdge(greatPoint324, greatPoint424)
114 greatEdgeY11 = MakeEdge(greatPoint111, greatPoint121)
115 greatEdgeY21 = MakeEdge(greatPoint211, greatPoint221)
116 greatEdgeY31 = MakeEdge(greatPoint311, greatPoint321)
117 greatEdgeY41 = MakeEdge(greatPoint411, greatPoint421)
119 greatEdgeY12 = MakeEdge(greatPoint112, greatPoint122)
120 greatEdgeY22 = MakeEdge(greatPoint212, greatPoint222)
121 greatEdgeY32 = MakeEdge(greatPoint312, greatPoint322)
122 greatEdgeY42 = MakeEdge(greatPoint412, greatPoint422)
124 greatEdgeY13 = MakeEdge(greatPoint113, greatPoint123)
125 greatEdgeY23 = MakeEdge(greatPoint213, greatPoint223)
126 greatEdgeY33 = MakeEdge(greatPoint313, greatPoint323)
127 greatEdgeY43 = MakeEdge(greatPoint413, greatPoint423)
129 greatEdgeY14 = MakeEdge(greatPoint114, greatPoint124)
130 greatEdgeY24 = MakeEdge(greatPoint214, greatPoint224)
131 greatEdgeY34 = MakeEdge(greatPoint314, greatPoint324)
132 greatEdgeY44 = MakeEdge(greatPoint414, greatPoint424)
134 greatEdgeZ111 = MakeEdge(greatPoint111, greatPoint112)
135 greatEdgeZ211 = MakeEdge(greatPoint211, greatPoint212)
136 greatEdgeZ311 = MakeEdge(greatPoint311, greatPoint312)
137 greatEdgeZ411 = MakeEdge(greatPoint411, greatPoint412)
139 greatEdgeZ121 = MakeEdge(greatPoint121, greatPoint122)
140 greatEdgeZ221 = MakeEdge(greatPoint221, greatPoint222)
141 greatEdgeZ321 = MakeEdge(greatPoint321, greatPoint322)
142 greatEdgeZ421 = MakeEdge(greatPoint421, greatPoint422)
144 greatEdgeZ112 = MakeEdge(greatPoint112, greatPoint113)
145 greatEdgeZ212 = MakeEdge(greatPoint212, greatPoint213)
146 greatEdgeZ312 = MakeEdge(greatPoint312, greatPoint313)
147 greatEdgeZ412 = MakeEdge(greatPoint412, greatPoint413)
149 greatEdgeZ122 = MakeEdge(greatPoint122, greatPoint123)
150 greatEdgeZ222 = smallEdgeZ11
151 greatEdgeZ322 = smallEdgeZ21
152 greatEdgeZ422 = MakeEdge(greatPoint422, greatPoint423)
154 greatEdgeZ113 = MakeEdge(greatPoint113, greatPoint114)
155 greatEdgeZ213 = MakeEdge(greatPoint213, greatPoint214)
156 greatEdgeZ313 = MakeEdge(greatPoint313, greatPoint314)
157 greatEdgeZ413 = MakeEdge(greatPoint413, greatPoint414)
159 greatEdgeZ123 = MakeEdge(greatPoint123, greatPoint124)
160 greatEdgeZ223 = MakeEdge(greatPoint223, greatPoint224)
161 greatEdgeZ323 = MakeEdge(greatPoint323, greatPoint324)
162 greatEdgeZ423 = MakeEdge(greatPoint423, greatPoint424)
167 smallFaceX1 = MakeQuad(smallEdgeY11, smallEdgeZ11, smallEdgeY12, smallEdgeZ12)
168 smallFaceX2 = MakeQuad(smallEdgeY21, smallEdgeZ21, smallEdgeY22, smallEdgeZ22)
169 smallFaceY1 = MakeQuad(smallEdgeX11, smallEdgeZ11, smallEdgeX12, smallEdgeZ21)
170 smallFaceY2 = MakeQuad(smallEdgeX21, smallEdgeZ12, smallEdgeX22, smallEdgeZ22)
171 smallFaceZ1 = MakeQuad(smallEdgeX11, smallEdgeY11, smallEdgeX21, smallEdgeY21)
172 smallFaceZ2 = MakeQuad(smallEdgeX12, smallEdgeY12, smallEdgeX22, smallEdgeY22)
175 greatFaceX11 = MakeQuad(greatEdgeY11, greatEdgeZ111, greatEdgeY12, greatEdgeZ121)
176 greatFaceX21 = MakeQuad(greatEdgeY21, greatEdgeZ211, greatEdgeY22, greatEdgeZ221)
177 greatFaceX31 = MakeQuad(greatEdgeY31, greatEdgeZ311, greatEdgeY32, greatEdgeZ321)
178 greatFaceX41 = MakeQuad(greatEdgeY41, greatEdgeZ411, greatEdgeY42, greatEdgeZ421)
180 greatFaceX12 = MakeQuad(greatEdgeY12, greatEdgeZ112, greatEdgeY13, greatEdgeZ122)
181 greatFaceX22 = MakeQuad(greatEdgeY22, greatEdgeZ212, greatEdgeY23, greatEdgeZ222)
182 greatFaceX32 = MakeQuad(greatEdgeY32, greatEdgeZ312, greatEdgeY33, greatEdgeZ322)
183 greatFaceX42 = MakeQuad(greatEdgeY42, greatEdgeZ412, greatEdgeY43, greatEdgeZ422)
185 greatFaceX13 = MakeQuad(greatEdgeY13, greatEdgeZ113, greatEdgeY14, greatEdgeZ123)
186 greatFaceX23 = MakeQuad(greatEdgeY23, greatEdgeZ213, greatEdgeY24, greatEdgeZ223)
187 greatFaceX33 = MakeQuad(greatEdgeY33, greatEdgeZ313, greatEdgeY34, greatEdgeZ323)
188 greatFaceX43 = MakeQuad(greatEdgeY43, greatEdgeZ413, greatEdgeY44, greatEdgeZ423)
190 greatFaceY111 = MakeQuad(greatEdgeX111, greatEdgeZ111, greatEdgeX112, greatEdgeZ211)
191 greatFaceY211 = MakeQuad(greatEdgeX211, greatEdgeZ211, greatEdgeX212, greatEdgeZ311)
192 greatFaceY311 = MakeQuad(greatEdgeX311, greatEdgeZ311, greatEdgeX312, greatEdgeZ411)
193 greatFaceY121 = MakeQuad(greatEdgeX121, greatEdgeZ121, greatEdgeX122, greatEdgeZ221)
194 greatFaceY221 = MakeQuad(greatEdgeX221, greatEdgeZ221, greatEdgeX222, greatEdgeZ321)
195 greatFaceY321 = MakeQuad(greatEdgeX321, greatEdgeZ321, greatEdgeX322, greatEdgeZ421)
197 greatFaceY112 = MakeQuad(greatEdgeX112, greatEdgeZ112, greatEdgeX113, greatEdgeZ212)
198 greatFaceY212 = MakeQuad(greatEdgeX212, greatEdgeZ212, greatEdgeX213, greatEdgeZ312)
199 greatFaceY312 = MakeQuad(greatEdgeX312, greatEdgeZ312, greatEdgeX313, greatEdgeZ412)
200 greatFaceY122 = MakeQuad(greatEdgeX122, greatEdgeZ122, greatEdgeX123, greatEdgeZ222)
201 greatFaceY222 = smallFaceY1
202 greatFaceY322 = MakeQuad(greatEdgeX322, greatEdgeZ322, greatEdgeX323, greatEdgeZ422)
204 greatFaceY113 = MakeQuad(greatEdgeX113, greatEdgeZ113, greatEdgeX114, greatEdgeZ213)
205 greatFaceY213 = MakeQuad(greatEdgeX213, greatEdgeZ213, greatEdgeX214, greatEdgeZ313)
206 greatFaceY313 = MakeQuad(greatEdgeX313, greatEdgeZ313, greatEdgeX314, greatEdgeZ413)
207 greatFaceY123 = MakeQuad(greatEdgeX123, greatEdgeZ123, greatEdgeX124, greatEdgeZ223)
208 greatFaceY223 = MakeQuad(greatEdgeX223, greatEdgeZ223, greatEdgeX224, greatEdgeZ323)
209 greatFaceY323 = MakeQuad(greatEdgeX323, greatEdgeZ323, greatEdgeX324, greatEdgeZ423)
211 greatFaceZ11 = MakeQuad(greatEdgeX111, greatEdgeY11, greatEdgeX121, greatEdgeY21)
212 greatFaceZ21 = MakeQuad(greatEdgeX211, greatEdgeY21, greatEdgeX221, greatEdgeY31)
213 greatFaceZ31 = MakeQuad(greatEdgeX311, greatEdgeY31, greatEdgeX321, greatEdgeY41)
215 greatFaceZ12 = MakeQuad(greatEdgeX112, greatEdgeY12, greatEdgeX122, greatEdgeY22)
216 greatFaceZ22 = MakeQuad(greatEdgeX212, greatEdgeY22, greatEdgeX222, greatEdgeY32)
217 greatFaceZ32 = MakeQuad(greatEdgeX312, greatEdgeY32, greatEdgeX322, greatEdgeY42)
219 greatFaceZ13 = MakeQuad(greatEdgeX113, greatEdgeY13, greatEdgeX123, greatEdgeY23)
220 greatFaceZ23 = MakeQuad(greatEdgeX213, greatEdgeY23, greatEdgeX223, greatEdgeY33)
221 greatFaceZ33 = MakeQuad(greatEdgeX313, greatEdgeY33, greatEdgeX323, greatEdgeY43)
223 greatFaceZ14 = MakeQuad(greatEdgeX114, greatEdgeY14, greatEdgeX124, greatEdgeY24)
224 greatFaceZ24 = MakeQuad(greatEdgeX214, greatEdgeY24, greatEdgeX224, greatEdgeY34)
225 greatFaceZ34 = MakeQuad(greatEdgeX314, greatEdgeY34, greatEdgeX324, greatEdgeY44)
230 smallBlock = MakeHexa(smallFaceX1, smallFaceX2, smallFaceY1, smallFaceY2, smallFaceZ1, smallFaceZ2)
232 greatBlock11 = MakeHexa(greatFaceX11, greatFaceX21, greatFaceY111, greatFaceY121, greatFaceZ11, greatFaceZ12)
233 greatBlock21 = MakeHexa(greatFaceX21, greatFaceX31, greatFaceY211, greatFaceY221, greatFaceZ21, greatFaceZ22)
234 greatBlock31 = MakeHexa(greatFaceX31, greatFaceX41, greatFaceY311, greatFaceY321, greatFaceZ31, greatFaceZ32)
236 greatBlock12 = MakeHexa(greatFaceX12, greatFaceX22, greatFaceY112, greatFaceY122, greatFaceZ12, greatFaceZ13)
237 greatBlock22 = MakeHexa(greatFaceX22, greatFaceX32, greatFaceY212, greatFaceY222, greatFaceZ22, greatFaceZ23)
238 greatBlock32 = MakeHexa(greatFaceX32, greatFaceX42, greatFaceY312, greatFaceY322, greatFaceZ32, greatFaceZ33)
240 greatBlock13 = MakeHexa(greatFaceX13, greatFaceX23, greatFaceY113, greatFaceY123, greatFaceZ13, greatFaceZ14)
241 greatBlock23 = MakeHexa(greatFaceX23, greatFaceX33, greatFaceY213, greatFaceY223, greatFaceZ23, greatFaceZ24)
242 greatBlock33 = MakeHexa(greatFaceX33, greatFaceX43, greatFaceY313, greatFaceY323, greatFaceZ33, greatFaceZ34)
248 c_l.append(smallBlock)
249 c_l.append(greatBlock11)
250 c_l.append(greatBlock21)
251 c_l.append(greatBlock31)
252 c_l.append(greatBlock12)
253 c_l.append(greatBlock22)
254 c_l.append(greatBlock32)
255 c_l.append(greatBlock13)
256 c_l.append(greatBlock23)
257 c_l.append(greatBlock33)
259 piece = MakeCompound(c_l)
264 piece_id = addToStudy(piece, "ex01_cube2build")
269 # Create hexahedrical mesh on piece
270 # ---------------------------------
272 hexa = smesh.Mesh(piece, "ex01_cube2build:hexa")
274 algo = hexa.Segment()
275 algo.NumberOfSegments(4)
281 # Create local hypothesis
282 # -----------------------
284 algo = hexa.Segment(greatEdgeX111)
286 algo.Arithmetic1D(1, 4)