求助~！UDEC进行渗流分析中出错

jia_yi_

我现正在用UDEC 4.0做节理岩质边坡的渗流分析，遇上了自己查找不到原因的问题，但又找不到人问，非常着急。计算采用不可压缩的瞬态流(set flow transient)，首先让模型在自重作用下达到平衡，再进行渗流分析。在进行渗流计算的时候，就一直提示“warning contact overlap with zero force”，最后提示“contact overlap too great”，模型损坏得面目全非。我查了手册中的Error messages，没能找到关于第一条提示的解释；对照着手册写了段简单的命令流反复试算也没能找出问题究竟出在哪儿。下面是试算模型的命令流，请大家看看，万分感谢~！

new
config tflow
round 0.2


block -117.04 -300.09 -117.04 -0.09 0.00 0.00 17.94 6.60 119.21 76.16 &
      172.07 116.49 294.92 230.01 357.07 280.22 420.21 347.01 438.38 360.62 &
      463.55 368.01 539.95 411.20 574.85 440.93 622.83 479.01 653.70 505.02 &
      669.97 516.03 742.64 565.21 770.78 588.16 804.07 624.72 847.15 676.12 &
      883.38 696.94 895.79 700.89 953.77 704.93 1012.15 709.00 1063.85 714.18 1120.63 721.39 &
      1146.00 726.06 1170.96 740.78 1213.16 748..99 1342.96 759.49 1342.96 -300.09


jset (-75,0) (100,0) (0,0) (22.8,0) id=1


gen edge 30


chan cons=3 jcons=2
prop mat=1 dens=2700 shear=1.23e10 bulk=1.79e10 coh=3.2e6 fric=50 tens=2e6


prop jmat=1 jkn=3e10 jks=1.2e10 jfric=28 jcoh=0.15e6 jten=0
prop jmat=1 jperm=1e8 azero=5e-4 ares=1e-4


grav 0,-9.81
fluid dens 1000


table 1 0.00 0.00 79.20 5.80 221.53 16.93 426.02 36.26 655.31 58.95 &
        839.76 79.78 1032.30 102.74 1342.96 141.15
insitu ywtable table 1
insitu stress -12.72e6 0 0 ygrad 16.75e3 0 0


bound xvel 0 yvel 0 range x -118 -116
bound xvel 0 yvel 0 range x 1342 1343
bound xvel 0 yvel 0 range y -301 -299


hist unbal
hist xdisp 900.10 701.19
hist ydisp 900.10 701.19


damp auto
solve
save 1.sav
pause
reset hist time


set caprat 100
set flow transient
set dtflow 10 voltol 1e-3 maxmech 500


bound pp 0.8e6 pygrad -0.01e6 range -120,208 -299,80
bound pp 1e6 pygrad -0.01e6 range 1340,1345 -301,100
bound stress 0,0,-0.8e6 ygrad 0,0,0.01e6 range -120,208 -1,80


step ftime 1200
save 2.sav

周泽

请问你还在研究这一块的内容么。可以相互交流一下么？