求助：Ducan-Chang模型的压缩试验

gaodengdeng

最近在学习FISH语言自定义本构模型， C:\Program Files\Itasca\flac500\gui\fishlib\ConstitutiveModels
有个Ducan-Chang模型的压缩试验，有一行命令行看不太懂，特来请教下各位朋友，命令行如下：

; mduncan.dat
config axi
g 5 20
gen 0 0 0 1 .25 1 .25 0
ca mduncan.fis
; --- model and properties ---
mo m_duncan
prop den .00202
prop d_pa=1.0584 d_k=700 d_n=0.37 d_rf=0.8 d_ku=1820
prop d_kb=280 d_m=0.19
prop d_coh=0.31 d_fric=33
; --- initialisation ---
ini sxx -2 syy -2 szz -2
set d_ms3=2.
ini_duncan
; --- boundary conditions ---
fix y j=1
fix y j=21
ini yv -0.5e-6 j=21
ini yv  0.5e-6 j=1
apply pres 2 i=6
; --- fish functions ---
def s1_s3
   _area=pi*x(igp,jgp)^2
   sum=yforce(1,jgp)*x(2,jgp)*0.25
   loop i (2,igp)
      sum=sum+yforce(i,jgp)*x(i,jgp)
   end_loop
   sigmav=2*pi*sum/_area
   s1_s3=sigmav-2.0
   s1_s3_ult=7.4076
   trans_s1_s3=ev/(sigmav-2.0)         (什么含义？？？？？？？)
end
def ev
   ev=(ydisp(3,1)-ydisp(3,21))/(y(3,21)-y(3,1))
end
; --- histories ---
hist nstep 20
hist s1_s3
hist s1_s3_ult
hist trans_s1_s3
hist ev
hist unbal
hist syy i=3 j=10

gaodengdeng

请教各位大侠，有人用过内置的Ducan-Chang模型算过吗？
即使一个简单的沉降问题，计算时都会产生"Divided by zero"的错误，不知道怎么回事

wmy_32552316wmy

呵呵，大家都很热心的呃说

yushu1983

是会出现像二楼所说的问题，FLAC里面自带的duncan-chang模型是不能够直接算填筑的。

yushu1983

是会出现像二楼所说的问题，FLAC里面自带的duncan-chang模型是不能够直接算填筑的。

recally

因为应力必须初始化
一般都是假设0.1倍大气压