new set random wall id=1 nodes -5,7 -5,3 -2,0 wall id=2 nodes 2,0 5,3 5,7 wall id=3 nodes -6,0 -6,-7 6,-7 6,0 wall id=1 kn=1e16 ks=1e16 fric 0.1 wall id=2 kn=1e16 ks=1e16 fric 0.1 wall id=3 kn=1e16 ks=1e16 fric 0.2 ; set hist_rep 100 history wall yforce id 3 ; ; ball id 1 x 0. y 5. rad 0.4 ball id 2 x -0.1 y 4.5 rad 0.3 ball id 3 x 0.3 y 4.6 rad 0.3 clump add id 101 range id 1 3 ball id 4 x -3 y 3 rad 0.5 ball id 5 x -3.5 y 3 rad 0.4 clump add id 102 range id 4 5 ball id 6 x 2.9 y 3.5 rad 0.25 ball id 7 x 3.2 y 3.6 rad 0.3 ball id 8 x 3 y 4 rad 0.3 clump add id 103 range id 6 8 ; plot set background white plot add wall black ;plot add ball white shade on id on plot add clump yellow plot add cforce blue plot show ;pause ; property density 2000 kn 1e8 ks 1e8 property friction 0.25 ; set grav 0. -9.81 cycle 5000 ;pause ; ball id 11 x 0. y 5. rad 0.4 ball id 12 x -0.1 y 4.5 rad 0.3 ball id 13 x 0.3 y 4.6 rad 0.3 clump add id 201 range id 11 13 ball id 14 x -3 y 3 rad 0.5 ball id 15 x -3.5 y 3 rad 0.4 clump add id 202 range id 14 15 ball id 16 x 2.9 y 3.5 rad 0.25 ball id 17 x 3.2 y 3.6 rad 0.3 ball id 18 x 3 y 4 rad 0.3 clump add id 203 range id 16 18 property density 2000 kn 1e8 ks 1e8 property friction 0.25 cycle 5000 ;pause ; ball id 21 x 0. y 5. rad 0.4 ball id 22 x -0.1 y 4.5 rad 0.3 ball id 23 x 0.3 y 4.6 rad 0.3 clump add id 301 range id 21 23 ball id 24 x -3 y 3 rad 0.5 ball id 25 x -3.5 y 3 rad 0.4 clump add id 302 range id 24 25 ball id 26 x 2.9 y 3.5 rad 0.25 ball id 27 x 3.2 y 3.6 rad 0.3 ball id 28 x 3 y 4 rad 0.3 clump add id 303 range id 26 28 property density 2000 kn 1e8 ks 1e8 property friction 0.25 cycle 5000 ;pause ; ball id 31 x 0. y 5. rad 0.4 ball id 32 x -0.1 y 4.5 rad 0.3 ball id 33 x 0.3 y 4.6 rad 0.3 clump add id 401 range id 31 33 ball id 34 x -3 y 3 rad 0.5 ball id 35 x -3.5 y 3 rad 0.4 clump add id 402 range id 34 35 ball id 36 x 2.9 y 3.5 rad 0.25 ball id 37 x 3.2 y 3.6 rad 0.3 ball id 38 x 3 y 4 rad 0.3 clump add id 403 range id 36 38 property density 2000 kn 1e8 ks 1e8 property friction 0.25 cycle 5000 ;pause ; ball id 41 x 0. y 5. rad 0.4 ball id 42 x -0.1 y 4.5 rad 0.3 ball id 43 x 0.3 y 4.6 rad 0.3 clump add id 501 range id 41 43 ball id 44 x -3 y 3 rad 0.5 ball id 45 x -3.5 y 3 rad 0.4 clump add id 502 range id 44 45 ball id 46 x 2.9 y 3.5 rad 0.25 ball id 47 x 3.2 y 3.6 rad 0.3 ball id 48 x 3 y 4 rad 0.3 clump add id 503 range id 46 48 property density 2000 kn 1e8 ks 1e8 property friction 0.25 cycle 5000 ;pause ; ball id 51 x 0. y 5. rad 0.4 ball id 52 x -0.1 y 4.5 rad 0.3 ball id 53 x 0.3 y 4.6 rad 0.3 clump add id 601 range id 51 53 ball id 54 x -3 y 3 rad 0.5 ball id 55 x -3.5 y 3 rad 0.4 clump add id 602 range id 54 55 ball id 56 x 2.9 y 3.5 rad 0.25 ball id 57 x 3.2 y 3.6 rad 0.3 ball id 58 x 3 y 4 rad 0.3 clump add id 603 range id 56 58 property density 2000 kn 1e8 ks 1e8 property friction 0.25 cycle 5000 ;pause ; ball id 61 x 0. y 5. rad 0.4 ball id 62 x -0.1 y 4.5 rad 0.3 ball id 63 x 0.3 y 4.6 rad 0.3 clump add id 701 range id 61 63 ball id 64 x -3 y 3 rad 0.5 ball id 65 x -3.5 y 3 rad 0.4 clump add id 702 range id 64 65 ball id 66 x 2.9 y 3.5 rad 0.25 ball id 67 x 3.2 y 3.6 rad 0.3 ball id 68 x 3 y 4 rad 0.3 clump add id 703 range id 66 68 property density 2000 kn 1e8 ks 1e8 property friction 0.25 cycle 5000 ;pause ; ball id 71 x 0. y 5. rad 0.4 ball id 72 x -0.1 y 4.5 rad 0.3 ball id 73 x 0.3 y 4.6 rad 0.3 clump add id 801 range id 71 73 ball id 74 x -3 y 3 rad 0.5 ball id 75 x -3.5 y 3 rad 0.4 clump add id 802 range id 74 75 ball id 76 x 2.9 y 3.5 rad 0.25 ball id 77 x 3.2 y 3.6 rad 0.3 ball id 78 x 3 y 4 rad 0.3 clump add id 803 range id 76 78 property density 2000 kn 1e8 ks 1e8 property friction 0.25 cycle 5000 ;pause ; damp viscous normal 0.7 damp viscous shear 0.3 cycle 20000 , ;save pelda2_final.sav ;plot hist 1