T3__horiz__assym_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

rotate Volume 7  angle 45  about X include_merged

move Volume 7 midpoint location curve 91  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 7


#!Mesh All

mesh surface 70

mesh volume all


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