# Rigid Voronoi slope model collapse **URL:** \ **Category:** UDEC\ **Created:** [May 5, 2022, 9:34am UTC](https://forum.itascainternational.com/t/rigid-voronoi-slope-model-collapse/499 "2022-05-05T09:34:36Z")\ **Posts on this page:** 13\ **Page:** 1
**Author:** ![Akash77skyz](https://avatars.discourse-cdn.com/v4/letter/a/8dc957/32.png) [@Akash77skyz](https://forum.itascainternational.com/u/Akash77skyz)\ **Post date:** [May 5, 2022, 9:34am UTC](https://forum.itascainternational.com/t/rigid-voronoi-slope-model-collapse/499/1 "2022-05-05T09:34:36Z")
I generated a voronoi model in UDEC 6 and made the blocks rigid. When I try to solve it, the block just collapse downwards as shown in picture. I have already applied zero velocity on the lower boundaries. I have also sent the code. Can anyone tell me how I stop this phenomenon. Thank you. round 0.2 edge 0.4 block 0,0 0,300 100,300 476,100 576,100 576,0 vor edge 10 change mat 1 range -29.3828,605.9828 -14.9911,320.9908 property mat 1 density 2.654E3 group joint ‘weakj:ID104’ range -31.7827,611.3825 -39.5897,334.19 joint model area jks 2E5 jkn 2E8 jfriction 25 range group ‘weakj:ID104’ boundary xvel 0 yvel 0 range -28.7829,628.1816 -11.3913,2.408 boundary xvel 0 yvel 0 range 575.9844,603.5829 -43.1895,144.0003 boundary xvel 0 yvel 0 range -22.1832,1.2155 -22.1907,327.5904 set gravity=0.0 -9.81 solve ![3456](https://canada1.discourse-cdn.com/flex031/uploads/itasca/original/1X/a9ef25994ab7fae25276e8f3651406da6f6a8d2d.png) ![123445](https://canada1.discourse-cdn.com/flex031/uploads/itasca/original/1X/d2c92e7b5a77ec09863b686ff139e4f3637870e8.png) ---
**Author:** ![jscic](https://avatars.discourse-cdn.com/v4/letter/j/958977/32.png) [@jscic](https://forum.itascainternational.com/u/jscic)\ **Post date:** [May 5, 2022, 8:26pm UTC](https://forum.itascainternational.com/t/rigid-voronoi-slope-model-collapse/499/2 "2022-05-05T20:26:46Z")
Use [In Situ] → [Fix] for the boundaries; The jkn is not stiff enough and causing overlapping. I added an order on it. ![displacement](https://canada1.discourse-cdn.com/flex031/uploads/itasca/original/1X/cf3549774273ff382d0b84b0a8a7a1fb8907ed97.jpeg) ![Hist Unbal](https://canada1.discourse-cdn.com/flex031/uploads/itasca/original/1X/930e46dd155fe9a965aec90459afb33de524a7d0.jpeg) round 0.2 edge 0.4 block 0,0 0,300 100,300 476,100 576,100 576,0 vor edge 10 change mat 1 range -29.3828,605.9828 -14.9911,320.9908 property mat 1 density 2.654E3 group joint weakj:ID104 range -31.7827,611.3825 -39.5897,334.19 joint model area jks 2E5 jkn 2E9 jfriction 25 range group weakj:ID104 ;boundary xvel 0 yvel 0 range -28.7829,628.1816 -11.3913,2.408 ;boundary xvel 0 yvel 0 range 575.9844,603.5829 -43.1895,144.0003 ;boundary xvel 0 yvel 0 range -22.1832,1.2155 -22.1907,327.5904 fix range -11.8351,585.3726 -9.9694,8.2065 fix range -11.8351,9.1104 -23.1859,308.801 fix range 571.4983,584.0657 -9.5713,108.771 set gravity=0.0 -9.81 solve ---
**Author:** ![Akash77skyz](https://avatars.discourse-cdn.com/v4/letter/a/8dc957/32.png) [@Akash77skyz](https://forum.itascainternational.com/u/Akash77skyz)\ **Post date:** [May 6, 2022, 3:34am UTC](https://forum.itascainternational.com/t/rigid-voronoi-slope-model-collapse/499/3 "2022-05-06T03:34:11Z")
@jscic Thanks a lot. And I had some other questions regarding this model. This time I use the same slope model but using Mohr-coulomb deformable blocks. There is sometimes the overlapping errors (when stiffness are low or near to the young’s modulus), so I increase the stiffness as the software says, but the more I increase the stiffness, the model becomes more stronger and reduces the block deformations I want to see. What I want to see is like the voronoi blocks slide off or roll off as rockfall or slide off the slope like this paper I send below (like toppling or sliding). Can you give me some suggestions on I can achieve this. Thank you. Paper: [s00603-019-01855-z.pdf](https://forum.itascainternational.com/uploads/short-url/h9676CgB0mkfsidL5VXsDOWpBFk.pdf) (4.6 MB) My code: round 0.02 edge 0.04 block 0,0 0,300 100,300 476,100 576,100 576,0 vor edge 10 gen edge 100.0 range region (-25.6195,-47.4439) (-25.6195,319.6272) (633.2542,319.6272) (633.2542,-47.4439) group zone ‘User:ID143’ range -20.3834,616.7823 -41.3897,325.7905 zone model mohr density 2.654E3 bulk 1.66667E8 shear 7.69231E7 friction 29 cohesion 8.4E4 range group ‘User:ID143’ group joint ‘weakj:ID104’ range -55.7815,604.183 -33.5901,343.1896 joint model area jks 2E5 jkn 2E9 jfriction 25 range group ‘weakj:ID104’ boundary xvel 0 yvel 0 range -20.1652,601.074 -15.8092,2.7352 boundary xvel 0 yvel 0 range 573.8027,637.072 -40.3533,153.818 boundary xvel 0 yvel 0 range -33.2554,2.7427 -7.0824,319.6272 set gravity=0.0 -9.81 solve ![00](https://canada1.discourse-cdn.com/flex031/uploads/itasca/original/1X/1c1df16b84eab5b8c4356039b42411f4032c7c21.png) ![55](https://canada1.discourse-cdn.com/flex031/uploads/itasca/original/1X/4d2f6cc48e96982b50f4db27207637fade571343.png) ---
**Author:** ![jscic](https://avatars.discourse-cdn.com/v4/letter/j/958977/32.png) [@jscic](https://forum.itascainternational.com/u/jscic)\ **Post date:** [May 9, 2022, 7:45pm UTC](https://forum.itascainternational.com/t/rigid-voronoi-slope-model-collapse/499/4 "2022-05-09T19:45:28Z")
One of the Apply Examples in the Itasca’s website should give you a great idea. Please check the link: [https://www.itascacg.com/software/support/examples/rockslide-runout-prediction](https://www.itascacg.com/software/support/examples/rockslide-runout-prediction) ---
**Author:** ![itsmoho\_sen](https://yyz1.discourse-cdn.com/flex031/user_avatar/forum.itascainternational.com/itsmoho_sen/32/421_2.png) [@itsmoho\_sen](https://forum.itascainternational.com/u/itsmoho_sen)\ **Post date:** [July 26, 2022, 8:54am UTC](https://forum.itascainternational.com/t/rigid-voronoi-slope-model-collapse/499/5 "2022-07-26T08:54:36Z")
I am trying to re-produce this code in UDEC 7, I made little changes but still does not work, Here is my code: ; model new round 0.02 edge 0.04 block create polygon 0,0 0,300 100,300 476,100 576,100 576,0 vor edge 10 gen edge 100.0 range region (-25.6195,-47.4439) (-25.6195,319.6272) (533.2542,319.6272) (533.2542,-47.4439) group zone ‘ID1’ range -20.3834,516.7823 -41.3897,325.7905 zone model mohr density 2.654E3 bulk 1.66667E8 shear 7.69231E7 friction 29 cohesion 8.4E4 range group ‘ID1’ group joint ‘weakj1’ range -55.7815,604.183 -33.5901,343.1896 joint model area jks 2E5 jkn 2E9 jfriction 25 range group ‘weakj1’ boundary xvel 0 yvel 0 range -20.1652,601.074 -15.8092,2.7352 boundary xvel 0 yvel 0 range 573.8027,637.072 -40.3533,153.818 boundary xvel 0 yvel 0 range -33.2554,2.7427 -7.0824,319.6272 set gravity=0.0 -9.81 Cycle 1 ---
**Author:** ![udecsupport](https://avatars.discourse-cdn.com/v4/letter/u/e9c0ed/32.png) [@udecsupport](https://forum.itascainternational.com/u/udecsupport)\ **Post date:** [July 28, 2022, 2:51pm UTC](https://forum.itascainternational.com/t/rigid-voronoi-slope-model-collapse/499/6 "2022-07-28T14:51:01Z")
The entire command set has changed in UDEC 7. The commands you have listed are UDEC 6 commands (not UDEC 7 commands). There is an automatic command conversion tool in UDEC 7. Load the file into a text pane and press edit-\>command conversion. ---
**Author:** ![itsmoho\_sen](https://yyz1.discourse-cdn.com/flex031/user_avatar/forum.itascainternational.com/itsmoho_sen/32/421_2.png) [@itsmoho\_sen](https://forum.itascainternational.com/u/itsmoho_sen)\ **Post date:** [July 28, 2022, 4:42pm UTC](https://forum.itascainternational.com/t/rigid-voronoi-slope-model-collapse/499/8 "2022-07-28T16:42:06Z")
this is the updates now but i still have a problem with this line (block zone generate edge 100.0 range region -25.6195 -47.4439 -25.6195 319.6272 533.2542 319.6272 533.2542 -47.4439) ; ; model new block tolerance corner-round-length 0.02 block tolerance minimum-edge-length 0.04 block create polygon 0 0 0 300 100 300 476 100 576 100 576 0 block cut voronoi edge-maximum 10 block zone generate edge 100.0 range region -25.6195 -47.4439 -25.6195 319.6272 533.2542 319.6272 533.2542 -47.4439 block zone group ‘ID1’ range position-x -20.3834 516.7823 position-y -41.3897 325.7905 block zone cmodel assign mohr-coulomb density 2.654E3 bulk 1.66667E8 shear 7.69231E7 friction 29 cohesion 8.4E4 range group ‘ID1’ block contact group ‘weakj1’ range position-x -55.7815 604.183 position-y -33.5901 343.1896 block contact cmodel assign area stiffness-shear 2E5 stiffness-normal 2E9 friction 25 range group ‘weakj1’ block gridpoint apply velocity-x 0 velocity-y 0 range position-x -20.1652 601.074 position-y -15.8092 2.7352 block gridpoint apply velocity-x 0 velocity-y 0 range position-x 573.8027 637.072 position-y -40.3533 153.818 block gridpoint apply velocity-x 0 velocity-y 0 range position-x -33.2554 2.7427 position-y -7.0824 319.6272 model gravity 0.0 -9.81 ; block cycle 1 ; return ; ---
**Author:** ![udecsupport](https://avatars.discourse-cdn.com/v4/letter/u/e9c0ed/32.png) [@udecsupport](https://forum.itascainternational.com/u/udecsupport)\ **Post date:** [July 29, 2022, 3:24pm UTC](https://forum.itascainternational.com/t/rigid-voronoi-slope-model-collapse/499/9 "2022-07-29T15:24:05Z")
1. The quote marks (‘ID1’ ) you are using for strings are not accepted by UDEC. For example use ‘ID1’ instead. 2. You need to add a seed command to the block cut voronoi command (I used seed 10010), Otherwise you will get a different model every time you run the data file. That makes it impossible to repeat or debug. 3. You need to reduce the edge length in the block zone gen edge command to 7.0. 4. You still have rigid blocks outside of the block zone gen range. You have not assigned properties to these blocks. The model will not cycle. ---
**Author:** ![udecsupport](https://avatars.discourse-cdn.com/v4/letter/u/e9c0ed/32.png) [@udecsupport](https://forum.itascainternational.com/u/udecsupport)\ **Post date:** [July 29, 2022, 3:26pm UTC](https://forum.itascainternational.com/t/rigid-voronoi-slope-model-collapse/499/10 "2022-07-29T15:26:06Z")
I see that the forum is changing the quote marks that are displayed in the listings. You need to use the single quote marks that are to the left of the enter key. ---
**Author:** ![itsmoho\_sen](https://yyz1.discourse-cdn.com/flex031/user_avatar/forum.itascainternational.com/itsmoho_sen/32/421_2.png) [@itsmoho\_sen](https://forum.itascainternational.com/u/itsmoho_sen)\ **Post date:** [August 2, 2022, 6:30am UTC](https://forum.itascainternational.com/t/rigid-voronoi-slope-model-collapse/499/12 "2022-08-02T06:30:55Z")
do you have your code to share with me please? I think I also do not employ this properly (block cut voronoi command (I used seed 10010),) ---
**Author:** ![udecsupport](https://avatars.discourse-cdn.com/v4/letter/u/e9c0ed/32.png) [@udecsupport](https://forum.itascainternational.com/u/udecsupport)\ **Post date:** [August 2, 2022, 8:10pm UTC](https://forum.itascainternational.com/t/rigid-voronoi-slope-model-collapse/499/13 "2022-08-02T20:10:36Z")
model new block tolerance corner-round-length 0.02 block tolerance minimum-edge-length 0.04 block create polygon 0 0 0 300 100 300 476 100 576 100 576 0 block cut voronoi edge-maximum 10 seed 10010 block zone generate edge 7.0 range region -25.6195 -47.4439 -25.6195 319.6272 533.2542 319.6272 533.2542 -47.4439 block zone group ‘ID1’ range position-x -20.3834 516.7823 position-y -41.3897 325.7905 block zone cmodel assign mohr-coulomb density 2.654E3 bulk 1.66667E8 shear 7.69231E7 friction 29 cohesion 8.4E4 range group ‘ID1’ block contact group ‘weakj1’ range position-x -55.7815 604.183 position-y -33.5901 343.1896 block contact cmodel assign area stiffness-shear 2E5 stiffness-normal 2E9 friction 25 range group ‘weakj1’ block gridpoint apply velocity-x 0 velocity-y 0 range position-x -20.1652 601.074 position-y -15.8092 2.7352 block gridpoint apply velocity-x 0 velocity-y 0 range position-x 573.8027 637.072 position-y -40.3533 153.818 block gridpoint apply velocity-x 0 velocity-y 0 range position-x -33.2554 2.7427 position-y -7.0824 319.6272 model gravity 0.0 -9.81 ; block cycle 1 ---
**Author:** ![itsmoho\_sen](https://yyz1.discourse-cdn.com/flex031/user_avatar/forum.itascainternational.com/itsmoho_sen/32/421_2.png) [@itsmoho\_sen](https://forum.itascainternational.com/u/itsmoho_sen)\ **Post date:** [August 4, 2022, 4:23pm UTC](https://forum.itascainternational.com/t/rigid-voronoi-slope-model-collapse/499/15 "2022-08-04T16:23:16Z")
Thank you very much I had some error in the text and now this is fixed and people can use it. ; model new ; block tolerance corner-round-length 0.02 block tolerance minimum-edge-length 0.04 ; block create polygon 0 0 0 300 100 300 476 100 576 100 576 0 ; block cut voronoi edge-maximum 10 seed 10010 ; block zone generate edge 7.0 range region -25.6195 -47.4439 -25.6195 319.6272 633.2542 319.6272 633.2542 -47.4439 block zone group ‘ID1’ range position-x -20.3834 616.7823 position-y -41.3897 325.7905 ; block zone cmodel assign mohr-coulomb density 2.654E3 bulk 1.66667E8 shear 7.69231E7 friction 29 cohesion 8.4E4 range group ‘ID1’ ; block contact group ‘weakj1’ range position-x -55.7815 604.183 position-y -33.5901 343.1896 ; block contact cmodel assign area stiffness-shear 2E5 stiffness-normal 2E9 friction 25 range group ‘weakj1’ ; block gridpoint apply velocity-x 0 velocity-y 0 range position-x -20.1652 601.074 position-y -15.8092 2.7352 block gridpoint apply velocity-x 0 velocity-y 0 range position-x 573.8027 637.072 position-y -40.3533 153.818 block gridpoint apply velocity-x 0 velocity-y 0 range position-x -33.2554 2.7427 position-y -7.0824 319.6272 ; model gravity 0.0 -9.81 ; ;model cycle 1 model solve Also, I am a 3DEC license owner, shouldn’t I will be able to run this model for a few cycles? does not even allow me to run it for one cycle ![image](https://canada1.discourse-cdn.com/flex031/uploads/itasca/original/1X/ddbb371c98616917dc4f34a0abf12be1b0d7cf0f.png) ---
**Author:** ![udecsupport](https://avatars.discourse-cdn.com/v4/letter/u/e9c0ed/32.png) [@udecsupport](https://forum.itascainternational.com/u/udecsupport)\ **Post date:** [August 4, 2022, 9:28pm UTC](https://forum.itascainternational.com/t/rigid-voronoi-slope-model-collapse/499/16 "2022-08-04T21:28:06Z")
Mohsen, A 3DEC license does not allow you to use UDEC and vis versa. They are not the same product. In demo mode you can run small models. But this is much larger than that. ![~WRD000.jpg](https://canada1.discourse-cdn.com/flex031/uploads/itasca/original/1X/52d5f1c08649212fe6768d2cacdcf8df711810d9.jpeg) ![image001.jpg](https://canada1.discourse-cdn.com/flex031/uploads/itasca/original/1X/2485d26919316b1d80aa60f5ca8e2cb5bc0275fe.jpeg)