3 ###### FILE DESCRIPTION ##############
5 # JOURNAL FILE FOR TEST BF420B of reference
6 # R P Birkett CSDMC P71 AEA-RS-4168 Part 1
7 # Hinkley Point B and Hunterston Geometry
8 # External loading scenario
9 # Created 14/5/2014 by Owen Booler AMEC
10 # Last Modified 14/5/2014 by Owen Booler AMEC
12 ###### GEOMETRY##################
14 # creating the 45 degree sector Hunterston B fuel brick Geometry
15 create curve location 0 0 0 location 131.625 0 0
16 create curve location 131.625 0 0 location 192.18 0 0
17 create curve location 192.18 0 0 location 192.18 10.125 0
18 create curve location 192.18 10.125 0 location 229.9 10.125 0
19 create curve location 229.9 10.125 0 location 229.9 45.72995 0
20 curve 1 copy rotate 45 about z
21 create curve location at vertex 12 direction curve 6 length 60.555
22 create curve location at vertex 14 direction 1 -1 0 length 16.475
23 create curve location at vertex 16 direction 1 1 0 length 37.72
24 create curve location at vertex 18 direction 1 -1 0 length 29.255
25 create curve location at vertex 20 location at vertex 10
26 create curve arc center vertex 11 12 2
29 create surface curve 7 12 3 2 4 5 11 10 9 8
31 # Projecting 45 degree sector to make a 25mm thick octant
32 # change distance to vary the thickness units in mm
33 sweep surface 1 direction z distance 25
35 # Copy octant to create full 360 degree slice
36 Volume 1 copy reflect y
37 Volume 1 2 copy rotate 90 about z
38 Volume 1 2 copy rotate 180 about z
39 Volume 1 2 copy rotate 270 about z
42 # Put fillets into keyway corners radius = 1.6
43 modify curve 17 99 27 53 38 188 203 234 219 128 143 174 159 68 114 83 blend radius 1.6
46 # Simulating crack from test BF420B from ref: R P Birkett CSDMC P71 AEA-RS-4168 Part 1
47 # Initial crack size will be 5mm long spanning full thickness of the brick slice.
48 create vertex on curve 277 279 midpoint
49 create curve vertex 239 240
50 create curve location at vertex 240 direction surface 105 length -5
51 create curve location at vertex 239 direction surface 105 length -5
52 create curve vertex 244 242
53 create surface curve 307 308 309 310
54 create cylinder height 25 radius 5
55 move Volume 10 location vertex 239 include_merged
56 move Volume 10 z 12.5 include_merged
57 chop volume 1 with volume 10
58 webcut volume 11 with plane surface 113
59 # Delete surface defining crack plane
62 #Creating Loading Keys (Interstitial)
63 create curve location 0 0 0 location 0 138 0
64 create curve location 0 138 0 location 138 138 0
65 create curve location 138 138 0 location 138 101 0
66 create curve location 0 0 0 location 37 0 0
67 create curve location at vertex 272 direction 1 1 0 length 55.6777849
68 create curve location at vertex 270 direction -1 -1 0 length 55.6777849
69 create curve location at vertex 276 direction 1 -1 0 length 34
70 create curve location at vertex 274 direction 1 -1 0 length 34
71 create curve vertex 278 280
72 Create surface curve 345 346 347 348 349 350 351 352 353
73 sweep surface 137 perpendicular distance 100.3
74 align Volume 14 surface 145 with surface 64
75 align Volume 14 surface 144 with surface 124
76 Volume 14 copy reflect y
77 Volume 14 copy reflect x
78 Volume 16 copy reflect y
80 # Creating Cross head plate for loading
81 brick x 30 y 520 z 120
82 align Volume 18 surface 184 with surface 159
83 move Volume 18 General location between location surface 159 location surface 170 except x z include_merged
85 # Creating extra vertices for loadpath
86 create vertex on curve 453 midpoint
87 create vertex on curve 457 midpoint
93 ###### MATERIAL PROPERTIES ########################
95 ## Graphite Brick Slice
96 block 1 volume 11 12 13
97 block 1 name "MAT_ELASTIC_Graphite"
98 block 1 attribute count 10
99 #Youngs Modulus E=10.9GPa
100 #Applied Youngs Modulus = 10900*1e-4
101 block 1 attribute index 1 1.09
103 block 1 attribute index 2 0.2
105 ## Hardened steel keys
106 block 2 volume 18 14 17 15 16
107 block 2 name "MAT_ELASTIC_STEEL"
108 block 2 attribute count 10
109 #Youngs Modulus E=210GPa
110 #Applied Youngs Modulus = 210000*1e-4
111 block 2 attribute index 1 21
113 block 2 attribute index 2 0.
115 ####### BOUNDARY CONDITIONS ######################
117 Create Displacement name ' ' on surface 138 148 dof 1 dof 2 dof 3 fix 0
119 ####### LOADING SCENARIO ########################
121 # BF420B Failed at a load of approximately 6200N
122 # Apply load either by force or pressure Hash (#) out the commands that
123 # are not appropriate
125 #Apply load as pressure P=F/A=6200/62400
126 #Create Pressure on surface 186 magnitude 1
128 #Apply load as a force Ff=6200N
129 #create force on surface 186 force value 6.2e-4 direction x
131 ######## FRACTURE CONSTRAINTS ####################
133 # Create Block Set to plot load vs dsplacement graphts
134 block 3 vertex 347 350 354 351 346 355
135 block 3 name 'LoadPath'
137 # Creating Nodesets and Sidesets for fracture constraints
139 # 101 - Vertices - Not crack
140 nodeset 101 vertex all
141 nodeset 101 vertex 257 258 259 260 remove
143 # 100 - Bounding Edges - Not crack
144 sideset 100 curve all
145 sideset 100 curve 329 330 331 332 343 335 341 336 313 315 remove
146 sideset 100 curve 241 12 remove
148 # 102 - Bounding Surfaces - Not crack
149 sideset 102 surface all
150 sideset 102 surface 127 130 134 11 remove
153 sideset 10200 surface 11
155 # 200 - Crack Surface
156 sideset 200 surface 127
159 sideset 201 curve 329
163 brick x 20 y 20 z 120
164 move surface 212 location surface 186 include_merged
166 create force on surface 214 force value 6.2e-4 direction x
169 ###### MESH DEFINITION ######################
170 volume all scheme Tetmesh
172 volume all size auto factor 7
173 volume 11 13 size auto factor 7
174 # Crack Mesh refinement
175 #curve 329 interval 15
176 #curve 343 interval 2
177 #curve 335 interval 2
178 #curve 330 interval 2
179 #curve 332 interval 2
180 #curve 341 interval 2
181 #curve 336 interval 2
184 ###### SAVE MODEL ########################
186 # Adapt this to the file and directory you wish to save
188 save as "/projects/symphonyb/MOFEM_Dev/model dev/HNB_Brick/HNB_BF420B/Test_BF420B.cub" overwrite
191 ###### END OF FILE #########################
