mofem/html/brick__slice__test_8jou_source.html

reset


###### FILE DESCRIPTION ##############

#

# JOURNAL FILE FOR TEST BF420B of reference

# R P Birkett CSDMC P71 AEA-RS-4168 Part 1

# Hinkley Point B and Hunterston Geometry

# External loading scenario

# Created 14/5/2014 by Owen Booler AMEC

# Last Modified 14/5/2014 by Owen Booler AMEC

#

###### GEOMETRY##################

#

# creating the 45 degree sector Hunterston B fuel brick Geometry

create curve location 0 0 0 location 131.625 0 0

create curve location 131.625 0 0 location 192.18 0 0

create curve location 192.18 0 0 location 192.18 10.125 0

create curve location 192.18 10.125 0 location 229.9 10.125 0

create curve location 229.9 10.125 0 location 229.9 45.72995 0

curve 1  copy rotate 45 about z

create curve location at vertex 12   direction curve 6  length 60.555

create curve location at vertex 14   direction 1 -1 0 length 16.475

create curve location at vertex 16   direction 1 1 0 length 37.72

create curve location at vertex 18   direction 1 -1 0 length 29.255

create curve location at vertex 20   location at vertex 10

create curve arc center vertex 11  12  2

delete Curve 1

delete Curve 6

create surface curve 7 12 3 2 4 5 11 10 9 8

#

# Projecting 45 degree sector to make a 25mm thick octant

# change distance to vary the thickness units in mm

sweep surface 1  direction z  distance 25

#

# Copy octant to create full 360 degree slice

Volume 1  copy reflect y

Volume 1 2 copy rotate 90 about z

Volume 1 2 copy rotate 180 about z

Volume 1 2 copy rotate 270 about z

unite volume all

#

# Put fillets into keyway corners radius = 1.6

modify curve 17 99 27 53 38 188 203 234 219 128 143 174 159 68 114 83  blend radius 1.6

#

## Create crack

# Simulating crack from test BF420B from ref: R P Birkett CSDMC P71 AEA-RS-4168 Part 1

# Initial crack size will be 5mm long spanning full thickness of the brick slice.

create vertex on curve 277 279 midpoint

create curve vertex 239 240

create curve location at vertex 240 direction surface 105 length -5

create curve location at vertex 239 direction surface 105 length -5

create curve vertex 244 242

create surface curve 307 308 309 310

create cylinder height 25 radius 8

move Volume 10  location vertex 239  include_merged

move Volume 10  z 12.5 include_merged

chop volume 1 with volume 10

webcut volume 11 with plane surface 113

# Delete surface defining crack plane

delete body 9

#

#Creating Loading Keys (Interstitial)

create curve location 0 0 0 location 0 138 0

create curve location 0 138 0 location 138 138 0

create curve location 138 138 0 location 138 101 0

create curve location 0 0 0 location 37 0 0

create curve location at vertex 272   direction 1 1 0 length 55.6777849

create curve location at vertex 270   direction -1 -1 0 length 55.6777849

create curve location at vertex 276   direction 1 -1 0 length 34

create curve location at vertex 274   direction 1 -1 0 length 34

create curve vertex 278 280

Create surface curve 345 346 347 348 349 350 351 352 353

sweep surface 137 perpendicular distance 100.3

align Volume 14  surface 145  with surface 64

align Volume 14  surface 144  with surface 124

Volume 14  copy reflect y

Volume 14  copy reflect x

Volume 16  copy reflect y

#

# Creating Cross head plate for loading

brick x 30 y 520 z 120

align Volume 18 surface 184  with surface 159

move Volume 18  General location between location surface 159   location surface 170  except x z include_merged

#

# Creating extra vertices for loadpath

create vertex on curve 453  midpoint

create vertex on curve 457  midpoint

#

# Merge volumes

imprint volume all

merge volume all

#

###### MATERIAL PROPERTIES ########################

#

## Graphite Brick Slice

block 1 volume 11 12 13

block 1 name "MAT_ELASTIC_Graphite"

block 1 attribute count 10

#Youngs Modulus E=10.9GPa

#Applied Youngs Modulus = 10900*1e-4

block 1 attribute index 1 1.09

#Poissons ratio

block 1 attribute index 2 0.2

#

## Hardened steel keys

block 2 volume 18 14 17 15 16

block 2 name "MAT_ELASTIC_STEEL"

block 2 attribute count 10

#Youngs Modulus E=210GPa

#Applied Youngs Modulus = 210000*1e-4

block 2 attribute index 1 21

#Poissons ratio

block 2 attribute index 2 0.

#

####### BOUNDARY CONDITIONS ######################

#

Create Displacement name ' ' on surface 138 148 dof 1 dof 2 dof 3 fix 0

#

####### LOADING SCENARIO ########################

#

# BF420B Failed at a load of approximately 6200N

# Apply load either by force or pressure Hash (#) out the commands that

# are not appropriate

#

#Apply load as pressure P=F/A=6200/62400

#Create Pressure on surface 186 magnitude 1

#

#Apply load as a force Ff=6200N

#create force on surface 186 force value 6.2e-4 direction x

#

######## FRACTURE CONSTRAINTS ####################

#

# Create Block Set to plot load  vs dsplacement graphts

#block 3 vertex 347 350 354 351 346 355

#block 3 name 'LoadPath'

#

# Creating Nodesets and Sidesets for fracture constraints

#

## 101 - Vertices - Not crack

#nodeset 101 vertex all

#nodeset 101 vertex 257 258 259 260 remove

##

## 100 - Bounding Edges - Not crack

#sideset 100 curve all

#sideset 100 curve 329 330 331 332 343 335 341 336 313 315 remove

#sideset 100 curve 241 12 remove

##

## 102 - Bounding Surfaces - Not crack

#sideset 102 surface all

#sideset 102 surface 127 130 134 11 remove

##

## 10200 -

#sideset 10200  surface 11

##

## 200 - Crack Surface

#sideset 200 surface 127

##

## 201 - Crack front

#sideset 201 curve 329

##


####

brick x 20 y 20 z 120

move surface 212 location surface 186 include_merged

unite volume 18 19

#create force on surface 214 force value -6.2e-4 direction x

create force on surface 214 vector -0.00062 0 0 0 0 0  #added by zahur

####


###### MESH DEFINITION ######################

volume all scheme Tetmesh


#ADJUST FOR THE NEW VERSION OF MESH CUTTING

create sheet extended from surface 127

move curve 544 midpoint location curve 329 include_merged

surface 217 scale x 1 y 0.017 z 0.5 about vertex 257


surface 217 size auto factor 4

surface 217 scheme trimesh

mesh surface 217

sideset 400 surface 217

move surface 217 y 1 x 0.1  include_merged


unite volume 11 13


imprint volume 14 17 12 18 15 11 16

merge volume 14 17 12 18 15 11 16


block 5 vertex all in volume 14 12 16 18 17 15 11

block 5 vertex 258 257 remove

block 4 curve in volume 14 12 16 18 17 15 11

block 4 curve 336 341 329 343 335 315 313 278 255 remove


# Global refinement

volume all scheme tetmesh

volume 11 size auto factor 6

mesh volume 11

volume 14 17 12 18 15 16 size auto factor 6.5

mesh volume 12 to 18


###### SAVE MODEL ########################

#

# Adapt this to the file and directory you wish to save

#

save as "/Users/zahur/Desktop/Fracture/paper/brick_slice_test_paper.cub" overwrite

#

#

###### END OF FILE #########################