请教邓肯张模型

guihonglong

小弟现在学习用flac计算粗粒土加筋问题，想用邓肯张模拟粗粒土。flac本身带了一个mduncan.fis，userguide里面和刘波的那本书里面也都给出了一个fish。我在计算中碰到了两个问题，请教一下各位高手。用自带的mduncan.fis时，
set cd name C:\Program Files\Itasca\flac500\gui\fishlib\ConstitutiveModels
call Mduncan.fis
set cd back
model m_duncan notnull
prop dens 2000.0 d_k=700.0 d_ku=1820.0 d_kb=280.0 d_n=0.37 d_m=0.19 d_rf=0.8 d_pa=105840.0 d_coh=30000.0 d_fric=33.0 notnull
ini sxx -10000 syy -10000 szz -10000 notnull
边界条件已经确定，计算时就出现divide by zero。一步都不能计算。
用userguide(P246,3-128)里面和刘波的书（p499)中给出的fish，编译时会出现ex_1没有定义的错误。我也没搞懂fish中ex_1应该是什么，不过看来确实似乎没有定义。不知道问题出在什么地方。
小弟flac仍是个新手，请各位大侠不吝赐教。或者有好心人能提供一个可用的duncan.fis也可以。先谢谢了。

guihonglong

没有人回答啊？把fish手册里用邓肯张的例子（P204）第一句confix axi去掉的话，也就是用平面应变的情况算，也会出现divided by zero这种情况。

guihonglong

没有理？？？help !版主帮忙。。。。。。。高手指教一下。。。

Stphenche

邓肯张模型不能出现拉应力 估计你的出现了拉应力 可能没办法用邓肯张了