mofem/html/T3__horiz_8jou_source.html

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.0305

move Volume 7 location surface 55  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 64 63 52 51 49 50

block 2 name 'LoadPath'


#Create nodesets and sidesets for fracture constraints

#Boundary edges

sideset 100 curve all

sideset 100 curve 132 121 122 130 119 104 106 119 118 120 remove


#Boundary vertices

nodeset 101 vertex all

nodeset 101 vertex 65 66 56 57 55 58 remove


#Boundary surfaces

sideset 102 surface all

sideset 102 surface 73 74 81  7 17   remove


#Crack surface

sideset 200 surface 73


#Crack front

sideset 201 curve 119


#!Meshing

volume all scheme tetmesh

volume all size auto factor 7

volume 8 10 size auto factor 7

#!Mesh Crack Surface (those elemenst are streched)

surface 73  size auto factor 7


#!Mesh All

mesh surface 73

mesh volume all


save as "/Users/likask/MyBuild/mofem-bitbucket/mofem/meshes/T3_horiz.cub" overwrite