v0.14.0
Loading...
Searching...
No Matches
T3_corner_notch.jou
Go to the documentation of this file.
1reset
2
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
9unite volume 4 5
10sweep surface 51 perpendicular distance 0.0025
11unite volume 6 4 2
12
13##Create notch
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
20
21#Merge Volumes
22imprint volume all
23merge volume all
24
25##Material Properties
26block 1 volume all
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
34
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
38
39#Loading
40#Apply pressure with magnitude = F/(pad_w x pad_l)
41#F = 18kN, pad_w = 0.005m, pad_l = 0.100m
42#P = 36MN/m2 * 10e-8
43create pressure on surface 51 magnitude 0.36000000
44
45#Create Block Set with Vertices to ploat load-disp. path
46block 2 vertex 52 51 49 50
47block 2 name 'LoadPath'
48
49#Create nodesets and sidesets for fracture constraints
50#Boundary edges
51sideset 100 curve all
52sideset 100 curve 115 113 121 119 112 111 110 remove
53
54#Boundary vertices
55nodeset 101 vertex all
56nodeset 101 vertex 60 56 55 59 remove
57
58#Boundary surfaces
59sideset 102 surface all
60sideset 102 surface 76 70 72 7 17 remove
61
62#Crack surface
63sideset 200 surface 70
64
65#Crack front
66sideset 201 curve 110
67
68#!Meshing
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
74
75#!Mesh All
76#mesh surface 70
77mesh volume all
78
79#save as "/Users/likask/MyBuild/mofem-bitbucket/mofem/meshes/T3_corner_notch.cub" overwrite
80