T3_corner_notch.jou

Go to the documentation of this file.

reset
 
brick x 0.27 y 0.1 z 0.1
webcut volume 1  with plane xplane offset 0.1225 noimprint merge  
webcut volume 2  with plane xplane offset -0.1225 noimprint merge  
webcut volume 2  with plane xplane offset 0 noimprint merge  
webcut volume 2  with plane xplane offset 0.0025 noimprint merge  
webcut volume 4  with plane xplane offset -0.0025 noimprint merge  
unite volume 4 5
sweep surface 51 perpendicular distance 0.0025
unite volume 6 4 2  
 
##Create notch
create Cylinder height 0.3 radius 0.033
rotate Volume 7  angle 45  about X include_merged  
#move Volume 7 midpoint location curve 91  include_merged
move volume 7 y -0.113669 include_merged
chop volume 6  with volume 7   
webcut volume 8 with plane xplane offset 0 noimprint nomerge  
  
#Merge Volumes  
imprint volume all   
merge volume all   
  
##Material Properties   
block 1 volume all    
block 1 name "MAT_ELASTIC_Graphite"     
block 1 attribute count 10     
#Elastic modulus E = 10.9GPa     
#Poisson's ratio v = 0.2     
#block 1 attribute index 1 1.09e+10  
block 1 attribute index 1 109  
block 1 attribute index 2 0.2     
  
#Apply displacement bcs  
create Displacement  on curve 31 dof 1 dof 2 dof 3 fix 0.0  
create displacement on curve 14 dof 2 dof 3 fix 0.0  
  
#Loading    
#Apply pressure with magnitude = F/(pad_w x pad_l)     
#F = 18kN, pad_w = 0.005m, pad_l = 0.100m     
#P = 36MN/m2 * 10e-8  
create pressure on surface 51 magnitude 0.36000000     
  
#Create Block Set with Vertices to ploat load-disp. path  
block 2 vertex 52 51 49 50   
block 2 name 'LoadPath'  
  
#Create nodesets and sidesets for fracture constraints    
#Boundary edges    
sideset 100 curve all   
sideset 100 curve 115 113 121 119 112 111 110  remove 
 
#Boundary vertices    
nodeset 101 vertex all  
nodeset 101 vertex 60 56 55 59 remove    
  
#Boundary surfaces    
sideset 102 surface all    
sideset 102 surface 76 70 72 7 17 remove  
  
#Crack surface 
sideset 200 surface 70 
 
#Crack front 
sideset 201 curve 110 
 
#!Meshing  
volume all scheme tetmesh     
volume all size auto factor 7  
volume 8 10 size auto factor 8  
#!Mesh Crack Surface (those elemenst are streched) 
#surface 70  size auto factor 6 
 
#!Mesh All 
#mesh surface 70 
mesh volume all  
  
#save as "/Users/likask/MyBuild/mofem-bitbucket/mofem/meshes/T3_corner_notch.cub" overwrite