round 0.0001
edge 0.0002
block 0,-0.005 0,0.105 0.05,0.105 0.05,-0.005
crack 0,0 0.05,0
crack 0,0.1 0.05,0.1
vor edge 0.0025 round 1.0E-4 range -2.093E-2,5.414E-2 3.085E-2,8.913E-2
 jdelete
gen edge 0.01 range -2.384E-2,6.108E-2 -1.758E-3,0.1025
gen edge 0.1 range -7.331E-3,6.011E-2 0.1016,0.1081
gen edge 0.1 range -1.219E-2,5.83E-2 -4.672E-3,-2.036E-3
group zone 'sandblockmat:ID94' range -1.857E-2,5.914E-2 -1.481E-3,0.1008
group zone 'Steel:ID97' range -6.638E-3,6.094E-2 -6.753E-3,-1.342E-3
group zone 'Steel:ID97' range -1.566E-2,6.08E-2 0.1012,0.1055
zone model mohr density 1.68E3 bulk 8.33333E7 shear 3.84615E7 friction 34 cohesion 2.2E4 range group 'sandblockmat:ID94'
zone model elastic density 7.85E3 bulk 1.51515E11 shear 7.8125E10 range group 'Steel:ID97'
group joint 'Vorjoint:ID116' range -7.768E-3,5.9E-2 4.89E-5,9.993E-2
group joint 'PLATECONTACT:ID124' range -2.739E-3,5.156E-2 9.993E-2,0.1013
group joint 'PLATECONTACT:ID124' range -4.29E-3,5.301E-2 -1.342E-3,4.89E-5
joint model area jks 2E7 jkn 2E9 jfriction 10 jcohesion 4E1 range group 'Vorjoint:ID116'
joint model area jks 1E7 jkn 1E9 jfriction 0 jcohesion 0 range group 'PLATECONTACT:ID124'
boundary yvelocity -1E-3 range -6.135E-3,6.033E-2 0.1038,0.1074
boundary yvelocity 0 range -1.169E-2,5.95E-2 -6.901E-3,-3.848E-3
boundary xvelocity 0 range -2.389E-3,5.839E-2 -6.901E-3,-2.044E-3

def sigmav
sum = 0.0
z_count = 0
i_b = block_head
loop while i_b # 0 
i_z = b_zone (i_b)
loop while i_z # 0
sum = sum - z_syy (i_z)
i_z = z_next (i_z)
z_count = z_count + 1
end_loop
i_b = b_next (i_b)
end_loop
sigmav = sum/z_count
end

def ev 
ev = (gp_ydis (i_gb) - gp_ydis (i_gt))/ (gp_y (i_gt) - gp_y (i_gb))
end 

def set_lim 
i_gt = gp_near (2.5,9.5)
i_gb = gp_near (2.5,0.5)
end set_lim

hist sigmav
label hist 1
axial stress
hist ev
label hist 2
axial strain

solve step 100000