model new model large-strain on model random 12001 model domain extent -0.01 0.36 -0.01 0.05 -0.01 0.2 model domain condition stop stop stop contact cmat proximity 0.05 contact cmat add 1 model linear property kn 1e8 ks 1e8 ... range group 'wall1' matches 2 contact cmat modify 1 model linearpbond ... property kn 5e7 fric 0.25 dp_nratio 0.5 ... pb_kn 1e8 pb_ks 1e8 ... pb_ten 1e12 pb_coh 1e12 pb_fa 30 ... range group 'wall1' matches 2 contact cmat apply contact method bond ;is it required? contact cmat add 2 model rrlinear ... property kn 1e9 dp_nratio 0.2 ... range group 'poly1' matches 2 contact cmat add 3 model rrlinear ... property kn 1e9 dp_nratio 0.2 ... range group 'extra1' matches 2 contact cmat default type ball-ball model rrlinear ... property dp_nratio 0.2 contact cmat default type ball-facet model linear ... property kn 1e6 dp_nratio 0.2 geometry polygon create set 'poly' by-positions (0 0 0) (0 0 0.12) (0.12 0 0.12)... (0.12 0 0.045) (0.125 0 0.045) (0.125 0 0.18) (0.36 0 0.18) (0.36 0 0)... extrude (0 0.05 0) group 'poly1' geometry generate group 'extra1' box 0.12 0.125 0 0.05 0.045 0.12 set 'extra' ;front wall wall generate group 's1' id 1 polygon (0 0 0) (0 0 0.12) (0.12 0 0.12)... (0.12 0 0.045) (0.125 0 0.045) (0.125 0 0.18) (0.36 0 0.18) (0.36 0 0) triangulate ;back wall wall generate group 'ss' id 2 polygon (0 0.05 0) (0 0.05 0.12) (0.12 0.05 0.12)... (0.12 0.05 0.045) (0.125 0.05 0.045) (0.125 0.05 0.18) (0.36 0.05 0.18) (0.36 0.05 0) triangulate wall generate group 's3' id 3 polygon (0 0 0) (0 0.05 0) (0.36 0.05 0) (0.36 0 0) ; bottom wall wall generate group 's4' id 4 polygon (0 0 0) (0 0.05 0) (0 0.05 0.12) (0 0 0.12) ;lateral left wall wall generate group 's5' id 5 polygon (0.36 0 0) (0.36 0.05 0) (0.36 0.05 0.18) (0.36 0 0.18) ;lateral right wall wall generate group 's6'id 6 polygon (0 0 0.12) (0 0.05 0.12) (0.12 0.05 0.12) (0.12 0 0.12) ;top wall wall generate group 's7' id 7 polygon (0.125 0 0.18) (0.125 0.05 0.18) (0.36 0.05 0.18) (0.36 0 0.18) ;top wall wall generate group 's9' id 9 polygon (0.12 0 0.045) (0.12 0 0.12) (0.125 0 0.12) (0.125 0 0.045) ;front small rectangle wall generate group 's10' id 10 polygon (0.12 0.05 0.045) (0.12 0.05 0.12) (0.125 0.05 0.12) (0.125 0.05 0.045) ;back small rectangle wall generate group 's11' id 11 polygon (0.12 0 0.12) (0.12 0.05 0.12) (0.125 0.05 0.12) (0.125 0 0.12) model save '05_08_soildomain' ball distribute porosity 0.4 bin 1 radius 0.0015 0.0025 group 'poly1'... range position-z 0 0.12 geometry-space 'poly' count 1 ball distribute porosity 0.4 bin 1 radius 0.0015 0.0025 group 'extra1'... range position-y 0 0.05 geometry-space 'extra' count 1 ball attribute density 1600 damp 0.7 range geometry-space 'poly' count 1 ball attribute density 1600 range geometry-space 'extra' count 1 model clean model cycle 1000 calm 100 model save 'initial' model restore 'initial' measure create id 20 position 0.18 0.002 0.09 radius 0.025 measure history name '01' porosity id 20 ball delete range geometry-space 'poly' count 1 not wall delete range set id 9 wall delete range set id 10 wall delete range set id 11 geometry delete set 'extra' geometry generate group 'wall1' box 0.12 0.125 0 0.05 0.045 0.18 set 'wall' wall generate group 'w1' box 0.12 0.125 0 0.05 0.045 0.18 model save '08_05_wall' ball distribute porosity 0.4 bin 1 radius 0.00125 group 'wall1'... range geometry-space 'wall' count 1 ball distribute porosity 0.4 bin 1 radius 0.0015 0.0025 group 'poly1'... range position-z 0.12 0.14 geometry-space 'poly' count 1 ball distribute porosity 0.4 bin 1 radius 0.0015 0.0025 group 'poly1'... range position-z 0.14 0.16 geometry-space 'poly' count 1 ball distribute porosity 0.4 bin 1 radius 0.0015 0.0025 group 'poly1'... range position-z 0.16 0.18 geometry-space 'poly' count 1 ball attribute density 1600 damp 0.7 range geometry-space 'poly' count 1 ball attribute density 2400 range geometry-space 'wall' count 1 model clean measure create id 21 position 0.126 0.0025 0.17 radius 0.0035 measure delete range id 21 measure create id 21 position 0.129 0.0025 0.175 radius 0.004 measure create id 22 position 0.129 0.0025 0.165 radius 0.004 measure create id 23 position 0.129 0.0025 0.155 radius 0.004 measure create id 24 position 0.129 0.0025 0.145 radius 0.004 measure create id 25 position 0.129 0.0025 0.135 radius 0.004 measure create id 26 position 0.129 0.0025 0.125 radius 0.004 measure create id 27 position 0.129 0.025 0.175 radius 0.005 measure create id 28 position 0.129 0.025 0.165 radius 0.005 measure create id 29 position 0.129 0.025 0.155 radius 0.005 measure create id 30 position 0.129 0.025 0.145 radius 0.005 measure create id 31 position 0.129 0.025 0.135 radius 0.005 measure create id 32 position 0.129 0.025 0.125 radius 0.005