model new
model random 10000
block tol 1e-10

; make parallelepiped to fill with voronois

geometry import '0301block1.dxf'
;geometry import '0217.dxf'

block generate voronoi max-edge 0.026 

block export filename '0222-block'
model save '0222-block'

program return

model new
model random 10000
block tol 1e-6


geometry import '0301block2.dxf'

block generate voronoi max-edge 0.026

block export filename '0222-try1'

model save '0222-try1'

program return

model new
model random 10000
block tol 1e-6

; make parallelepiped to fill with voronois
geometry import '0301block3.dxf'
block generate voronoi max-edge 0.026

block export filename '0222-try2'

model save '0222-try2'

program return

model new
model random 10000
block tol 1e-6

; make parallelepiped to fill with voronois
geometry import '0301block4.dxf'
block generate voronoi max-edge 0.026

block export filename '0222-try3'

model save '0222-try3'

program return

model new
model random 10000
block tol 1e-6

; make parallelepiped to fill with voronois
geometry import '0301block5.dxf'
block generate voronoi max-edge 0.0255

block export filename '0222-try4'

model save '0222-try4'

program return

model new
model random 10000
block tol 1e-8

;geometry import '0301block6.dxf'
geometry import '0301-two.dxf'
block generate voronoi max-edge 0.0255

block export filename '0222-try5'

model save '0222-try5'

program return

model new
model random 10000
block tol 1e-6
model large-strain on 

block import filename '0222-block'
block import filename '0222-try1'
block import filename '0222-try2'
block import filename '0222-try3'
block import filename '0222-try4'
block import filename '0222-try5'
;block import filename '0222-try6'

block zone generate edge 0.1

block contact group 'weak' range joint-set 3 4 5 6 7 8  

block create prism &
face-1 0.0685948,0,0.249574  0.062,0,0.26  0.15,0,0.31  0.156,0,0.3 &
face-2 0.0685948,0.2,0.249574  0.062,0.2,0.26  0.15,0.2,0.31  0.156,0.2,0.3

block group 'metal' range id 2745

block zone cmodel elastic
block zone property density 734 young 5.335e8 poiss 0.3 range position-x 0 0.55

block property density 1e-3 range group 'metal'

block contact material-table default property stiffness-norm=7e8 stiffness-shear=2.8e8 friction=60 

block contact jmodel assign mohr
block contact prop stiffness-norm=2.89e9 stiffness-shear=7.8e8 friction=9.463  cohesion 0 range group 'weak'
;block contact material-table add jmodel mohr property &
;stiffness-norm=2.89e9 stiffness-shear=3.12e7 friction=9.463 friction-residual 9.463 range group 'weak'

block contact prop stiffness-norm=4.545e8 stiffness-shear=1.729e7 cohesion 4e5 friction=32 range joint-set 2
;block contact material-table add jmodel mohr property &
;stiffness-norm=4.545e9 stiffness-shear=1.729e8 cohesion 4e5 friction=32 cohesion-residual 4e4 friction-residual 32 range joint-set 2

block gridpoint apply velocity 0 0 0 range position-z 0 0.004
block gridpoint apply velocity 0 0 0 range position-x 0 0.004
block gridpoint apply velocity 0 0 0 range position-x 0.546 0.55

block gridpoint apply-remove velocity range group 'metal'

block gridpoint apply velocity-y 0 range position-x 0 0.55
block apply velocity-y 0 range group 'metal'

block tolerance 0

model gravity 0,0,-9.81

block gridpoint apply force-x 122.625  range id-list 55282 55281 55280 55279
block gridpoint apply force-z -212.393 range id-list 55282 55281 55280 55279


model solve mechanical ratio 1e-5 
;model solve convergence 1
;model solve cycle 600

model save '0222-zeroday'

program return