chdharden 发表于 2017-10-10 12:16:01

抽水命令流有点问题

自己写了一个关于节点抽水的命令流,可以最后solve后发现没有任何变化,既不出现塑性区孔压等也不变,不知道是哪里出现问题了,希望大家帮忙解决一下,另外想知道“apply pwell -0.1 interior range x 16 17 y 0 3.5 z 9 10

ini pp 0 range id 1350
fix pp range id 1350“这两句命令的具体作用怎么描述,感谢了。命令流如下:

new
config fluid
gen zone brick p0 0 0 0 p1 40 0 0 p2 0 2 0 p3 0 0 20 p4 40 2 0 p5 0 2 20 &
p6 40 0 20 p7 40 2 20 &
size 40 2 20

fix y
fix z range z -0.1 0.1
fix x range x -0.1 0.1         
fix x range x 39.9 40.1

model moh
prop bulk 1e8 shear 5e9 coh 3e6 tens 1e6 fric 35 dila 20
set gravity 0 0 -10.0

ini density 2040 range z 18 20      
ini density 2000 range z 0 18

model fl_iso
water density 1000
water table ori 0 0 18 norm 0 0 1
prop perm 1.0e-11 poro 0.1         

ini fmod 1e9
set fluid off mech on
solve
sav 应力场.sav
pause
ini xdisp 0 ydisp 0 zdisp 0 xvel 0 yvel 0 zvel 0
;apply pwell -0.1 interior range x 16 17 y 0 3.5 z 9 10

ini pp 0 range id 1350
fix pp range id 1350   

set fluid on mech off

solve
sav 渗流场.sav

1012693062 发表于 2017-11-6 10:53:29

关于第一个fix pp 的命令,它是把那个节点处的孔隙水压力力设为0,再用fix固定,就是让他一直为零
后面关于孔压无变化,pwell的问题,我个人认为:pwell是关于定义节点流量的流入和流出问题,应该不能用于抽水吧,pwell只是描述的流入和流出 ,水没有消失,所以你的孔压没有变化,不知说的对不对,望指正!

j1381642 发表于 2018-11-27 12:24:32

如何把孔压分布导进tecplot绘图
页: [1]
查看完整版本: 抽水命令流有点问题