3brick x 0.27 y 0.1 z 0.1
4webcut volume 1 with plane xplane offset 0.1225 noimprint merge
5webcut volume 2 with plane xplane offset -0.1225 noimprint merge
6webcut volume 2 with plane xplane offset 0 noimprint merge
7webcut volume 2 with plane xplane offset 0.0025 noimprint merge
8webcut volume 4 with plane xplane offset -0.0025 noimprint merge
10sweep surface 51 perpendicular distance 0.0025
14create Cylinder height 0.3 radius 0.033
15rotate Volume 7 angle 45 about X include_merged
16#move Volume 7 midpoint location curve 91 include_merged
17move volume 7 y -0.113669 include_merged
18chop volume 6 with volume 7
19webcut volume 8 with plane xplane offset 0 noimprint nomerge
27block 1 name "MAT_ELASTIC_Graphite"
28block 1 attribute count 10
29#Elastic modulus E = 10.9GPa
30#Poisson's ratio v = 0.2
31#block 1 attribute index 1 1.09e+10
32block 1 attribute index 1 109
33block 1 attribute index 2 0.2
35#Apply displacement bcs
36create Displacement on curve 31 dof 1 dof 2 dof 3 fix 0.0
37create displacement on curve 14 dof 2 dof 3 fix 0.0
40#Apply pressure with magnitude = F/(pad_w x pad_l)
41#F = 18kN, pad_w = 0.005m, pad_l = 0.100m
43create pressure on surface 51 magnitude 0.36000000
45#Create Block Set with Vertices to ploat load-disp. path
46block 2 vertex 52 51 49 50
47block 2 name 'LoadPath'
49#Create nodesets and sidesets for fracture constraints
52sideset 100 curve 115 113 121 119 112 111 110 remove
56nodeset 101 vertex 60 56 55 59 remove
59sideset 102 surface all
60sideset 102 surface 76 70 72 7 17 remove
69volume all scheme tetmesh
70volume all size auto factor 7
71volume 8 10 size auto factor 8
72#!Mesh Crack Surface (those elemenst are streched)
73#surface 70 size auto factor 6
79#save as "/Users/likask/MyBuild/mofem-bitbucket/mofem/meshes/T3_corner_notch.cub" overwrite