facesea 发表于 2011-5-22 22:52:54

求助,关于fish

定义了一个求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
不知道怎么回事,求助
页: [1]
查看完整版本: 求助,关于fish