blob: bfa5a3a0c4f9229e1bdc214fe65901d1234adafe (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
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
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
|
*DECK GeoHom
*----
* Name : GeoHom.c2m
* Type : DRAGON procedure
* Use : Geometry generation homogenization
* Author : R. Chambon
*
* Procedure called as:
*
*ALLGEO:= GeoHom :: <<Cote>> <<CoteL>> ;
*
* Input data :
* Cote : Lattice pitch (cm)
* CoteL : Lattice pitch + Water space (cm)
*
* Output data :
* ALLGEO : All Geometries for homogenizations
*
PARAMETER ALLGEO ::
EDIT 0
::: LINKED_LIST ALLGEO ;
;
LINKED_LIST GEOTMP ASSMB_HOM ;
*----
* Modules used in this procedure
*----
MODULE GEO: DELETE: BACKUP: END: ;
*----
* Input data recovery
*----
* ;
REAL Cote CoteL ;
:: >>Cote<< >>CoteL<< ;
STRING namgeohom namgeoht1 namgeoht2 namgeopin :=
"GEOHOM" "GEOHT1" "GEOHT2" "GEOPIN" ;
STRING namsphhom namsphht1 namsphht2 namsphpin :=
"SPHHOM" "SPHHT1" "SPHHT2" "SPHPIN" ;
*
* GEOMETRY for MULTICOMPO
*
*
*----
* Homogeneous assembly
*----
REAL Lassp Lassm ;
EVALUATE Lassp := Cote 15. * CoteL 2. * + 2. / ;
EVALUATE Lassm := Lassp -1. * ;
GEOTMP := GEO: :: CAR2D 1 1
EDIT 0
X- REFL X+ REFL
Y- REFL Y+ REFL
MIX 1
MESHX <<Lassm>> <<Lassp>>
MESHY <<Lassm>> <<Lassp>>
;
ALLGEO := BACKUP: GEOTMP ::
EDIT 0
STEP UP <<namgeohom>>
;
GEOTMP := DELETE: GEOTMP ;
*----
* Heterogeneous assembly
*----
REAL LHetp LHetm ;
EVALUATE LHetp := Cote 15. * 2. / ;
EVALUATE LHetm := LHetp -1. * ;
GEOTMP := GEO: :: CAR2D 3 3
EDIT 0
X- REFL X+ REFL
Y- REFL Y+ REFL
MIX 3 2 3
2 1 2
3 2 3
MESHX <<Lassm>> <<LHetm>> <<LHetp>> <<Lassp>>
MESHY <<Lassm>> <<LHetm>> <<LHetp>> <<Lassp>>
;
ALLGEO := BACKUP: ALLGEO GEOTMP ::
EDIT 0
STEP UP <<namgeoht1>>
;
GEOTMP := DELETE: GEOTMP ;
*
EVALUATE LHetp := Cote 13. * 2. / ;
EVALUATE LHetm := LHetp -1. * ;
GEOTMP := GEO: :: CAR2D 3 3
EDIT 0
X- REFL X+ REFL
Y- REFL Y+ REFL
MIX 3 2 3
2 1 2
3 2 3
MESHX <<Lassm>> <<LHetm>> <<LHetp>> <<Lassp>>
MESHY <<Lassm>> <<LHetm>> <<LHetp>> <<Lassp>>
;
ALLGEO := BACKUP: ALLGEO GEOTMP ::
EDIT 0
STEP UP <<namgeoht2>>
;
GEOTMP := DELETE: GEOTMP ;
*
*----
* Pin by pin assembly
*----
REAL Lp1p Lp1m Lp2p Lp2m Lp3p Lp3m Lp4p Lp4m
Lp5p Lp5m Lp6p Lp6m Lp7p Lp7m Lp8p Lp8m ;
EVALUATE Lp1p := Cote 1. * 2. / ;
EVALUATE Lp1m := Lp1p -1. * ;
EVALUATE Lp2p := Cote 3. * 2. / ;
EVALUATE Lp2m := Lp2p -1. * ;
EVALUATE Lp3p := Cote 5. * 2. / ;
EVALUATE Lp3m := Lp3p -1. * ;
EVALUATE Lp4p := Cote 7. * 2. / ;
EVALUATE Lp4m := Lp4p -1. * ;
EVALUATE Lp5p := Cote 9. * 2. / ;
EVALUATE Lp5m := Lp5p -1. * ;
EVALUATE Lp6p := Cote 11. * 2. / ;
EVALUATE Lp6m := Lp6p -1. * ;
EVALUATE Lp7p := Cote 13. * 2. / ;
EVALUATE Lp7m := Lp7p -1. * ;
EVALUATE Lp8p := Cote 15. * 2. / ;
EVALUATE Lp8m := Lp8p -1. * ;
GEOTMP := GEO: :: CAR2D 17 17
EDIT 0
X- REFL X+ REFL
Y- REFL Y+ REFL
MIX
45 44 42 39 35 30 24 17 9 17 24 30 35 39 42 44 45
44 43 41 38 34 29 23 16 8 16 23 29 34 38 41 43 44
42 41 40 37 33 28 22 15 7 15 22 28 33 37 40 41 42
39 38 37 36 32 27 21 14 6 14 21 27 32 36 37 38 39
35 34 33 32 31 26 20 13 5 13 20 26 31 32 33 34 35
30 29 28 27 26 25 19 12 4 12 19 25 26 27 28 29 30
24 23 22 21 20 19 18 11 3 11 18 19 20 21 22 23 24
17 16 15 14 13 12 11 10 2 10 11 12 13 14 15 16 17
9 8 7 6 5 4 3 2 1 2 3 4 5 6 7 8 9
17 16 15 14 13 12 11 10 2 10 11 12 13 14 15 16 17
24 23 22 21 20 19 18 11 3 11 18 19 20 21 22 23 24
30 29 28 27 26 25 19 12 4 12 19 25 26 27 28 29 30
35 34 33 32 31 26 20 13 5 13 20 26 31 32 33 34 35
39 38 37 36 32 27 21 14 6 14 21 27 32 36 37 38 39
42 41 40 37 33 28 22 15 7 15 22 28 33 37 40 41 42
44 43 41 38 34 29 23 16 8 16 23 29 34 38 41 43 44
45 44 42 39 35 30 24 17 9 17 24 30 35 39 42 44 45
MESHX <<Lassm>> <<Lp8m>> <<Lp7m>> <<Lp6m>> <<Lp5m>> <<Lp4m>>
<<Lp3m>> <<Lp2m>> <<Lp1m>> <<Lp1p>> <<Lp2p>> <<Lp3p>>
<<Lp4p>> <<Lp5p>> <<Lp6p>> <<Lp7p>> <<Lp8p>> <<Lassp>>
MESHY <<Lassm>> <<Lp8m>> <<Lp7m>> <<Lp6m>> <<Lp5m>> <<Lp4m>>
<<Lp3m>> <<Lp2m>> <<Lp1m>> <<Lp1p>> <<Lp2p>> <<Lp3p>>
<<Lp4p>> <<Lp5p>> <<Lp6p>> <<Lp7p>> <<Lp8p>> <<Lassp>>
;
ALLGEO := BACKUP: ALLGEO GEOTMP ::
EDIT 0
STEP UP <<namgeopin>>
;
GEOTMP := DELETE: GEOTMP ;
*
* GEOMETRY for SPH calculations
*
*
*
* Original geometry for flux calculation
* CELL TI C01 C03 T1 C08 C13 T2 C23 C31 L1
* C02 C04 C06 C09 C14 C18 C24 C32 L1
* C05 C07 C10 C15 C19 C25 C33 L1
* T1 C11 C16 T2 C26 C34 L1
* C12 C17 C20 C27 C35 L1
* T2 C21 C28 C36 L1
* C22 C29 C37 L1
* C30 C38 L1
* C39 L1
* L2
REAL Lp9p ;
EVALUATE Lp9p := Cote 17. * 2. / ;
*----
* Homogeneous assembly
*----
ASSMB_HOM := GEO: :: CAR2D 10 10
X- DIAG X+ REFL Y- SYME Y+ DIAG
MIX 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1
1 1 1 1 1 1
1 1 1 1 1
1 1 1 1
1 1 1
1 1
1
MESHX <<Lp1m>> <<Lp1p>> <<Lp2p>> <<Lp3p>> <<Lp4p>> <<Lp5p>>
<<Lp6p>> <<Lp7p>> <<Lp8p>> <<Lp9p>> <<Lassp>>
MESHY <<Lp1m>> <<Lp1p>> <<Lp2p>> <<Lp3p>> <<Lp4p>> <<Lp5p>>
<<Lp6p>> <<Lp7p>> <<Lp8p>> <<Lp9p>> <<Lassp>>
;
ALLGEO := BACKUP: ALLGEO ASSMB_HOM ::
EDIT 0
STEP UP <<namsphhom>>
;
ASSMB_HOM := DELETE: ASSMB_HOM ;
*----
* Heterogeneous assembly
*----
*(1 pin thick outer layer in DRAGON)
ASSMB_HOM := GEO: :: CAR2D 10 10
X- DIAG X+ REFL Y- SYME Y+ DIAG
MIX 1 1 1 1 1 1 1 1 2 2
1 1 1 1 1 1 1 2 2
1 1 1 1 1 1 2 2
1 1 1 1 1 2 2
1 1 1 1 2 2
1 1 1 2 2
1 1 2 2
1 2 2
3 3
3
MESHX <<Lp1m>> <<Lp1p>> <<Lp2p>> <<Lp3p>> <<Lp4p>> <<Lp5p>>
<<Lp6p>> <<Lp7p>> <<Lp8p>> <<Lp9p>> <<Lassp>>
MESHY <<Lp1m>> <<Lp1p>> <<Lp2p>> <<Lp3p>> <<Lp4p>> <<Lp5p>>
<<Lp6p>> <<Lp7p>> <<Lp8p>> <<Lp9p>> <<Lassp>>
;
ALLGEO := BACKUP: ALLGEO ASSMB_HOM ::
EDIT 0
STEP UP <<namsphht1>>
;
ASSMB_HOM := DELETE: ASSMB_HOM ;
*(2 pins thick outer layer in DRAGON)
ASSMB_HOM := GEO: :: CAR2D 10 10
X- DIAG X+ REFL Y- SYME Y+ DIAG
MIX 1 1 1 1 1 1 1 2 2 2
1 1 1 1 1 1 2 2 2
1 1 1 1 1 2 2 2
1 1 1 1 2 2 2
1 1 1 2 2 2
1 1 2 2 2
1 2 2 2
3 3 3
3 3
3
MESHX <<Lp1m>> <<Lp1p>> <<Lp2p>> <<Lp3p>> <<Lp4p>> <<Lp5p>>
<<Lp6p>> <<Lp7p>> <<Lp8p>> <<Lp9p>> <<Lassp>>
MESHY <<Lp1m>> <<Lp1p>> <<Lp2p>> <<Lp3p>> <<Lp4p>> <<Lp5p>>
<<Lp6p>> <<Lp7p>> <<Lp8p>> <<Lp9p>> <<Lassp>>
;
ALLGEO := BACKUP: ALLGEO ASSMB_HOM ::
EDIT 0
STEP UP <<namsphht2>>
;
ASSMB_HOM := DELETE: ASSMB_HOM ;
*----
* Pin by pin assembly
*----
ASSMB_HOM := GEO: :: CAR2D 10 10
X- DIAG X+ REFL Y- SYME Y+ DIAG
MIX 1 2 3 4 5 6 7 8 9 9
10 11 12 13 14 15 16 17 17
18 19 20 21 22 23 24 24
25 26 27 28 29 30 30
31 32 33 34 35 35
36 37 38 39 39
40 41 42 42
43 44 44
45 45
45
MESHX <<Lp1m>> <<Lp1p>> <<Lp2p>> <<Lp3p>> <<Lp4p>> <<Lp5p>>
<<Lp6p>> <<Lp7p>> <<Lp8p>> <<Lp9p>> <<Lassp>>
MESHY <<Lp1m>> <<Lp1p>> <<Lp2p>> <<Lp3p>> <<Lp4p>> <<Lp5p>>
<<Lp6p>> <<Lp7p>> <<Lp8p>> <<Lp9p>> <<Lassp>>
;
ALLGEO := BACKUP: ALLGEO ASSMB_HOM ::
EDIT 0
STEP UP <<namsphpin>>
;
ASSMB_HOM := DELETE: ASSMB_HOM ;
END: ;
QUIT .
|
