- 积分
- 0
- 注册时间
- 2009-1-19
- 仿真币
-
- 最后登录
- 1970-1-1
|
定义了一个求pore pressure 增量的函数
zpp=z_near(0.5,0.5,0.5)
if tstep=1
z_extra(zpp,1)=z_pp(zpp)
biotan1=5000;(1-a1)*z_pp(zpp)
biotan2=5000;(1-a2)*z_pp(zpp)
biotan3=5000;(1-a3)*z_pp(zpp)
else
biotan1=0;(1-a1)*(z_pp(zpp)-z_extra(zpp,1))
biotan2=0;(1-a2)*(z_pp(zpp)-z_extra(zpp,1))
biotan3=0;(1-a3)*(z_pp(zpp)-z_extra(zpp,1))
z_extra(zpp,1)=z_pp(zpp)
endif
在pore pressure 不变的情况下,我把增量设成第一步5000,其他设为0,可计算出来的结果每一部的biotan1,biotan2,biotan3都是0
如果换成后面的表达式,比如第一步计算出 biotan1,2,3是5000,后面的也是5000
不知道怎么回事,求助 |
|