求助一个关于flac3d的问题，知道的帮忙解决一下 [复制链接]

1.txt (3 K) 下载次数:15 ，这是程序，我不知道是哪里出了错

学习学习，呵呵

嘿嘿，恩，互相学习

贴出来吧！
要本人帮忙，还要别人下载附件，不大合适。
请斟酌。

;=========================
;剪模为0.09e8
;=========================
new
config dynamic
set dynamic dt 1.0e-5
set dynamic time 0.0
set dynamic multi on
set dynamic damp rayleigh 0.1 10
gen zone brick p0 0 0 0 p1 5 0 0 p2 0 1 0 p3 0 0 1 size 25 5 5 group 1
gen zone brick p0 5 0 0 p1 5.2 0 0 p2 5 1 0 p3 5 0 1 size 10 5 5 group 2
gen zone brick p0 5.2 0 0 p1 15.2 0 0 p2 5.2 1 0 p3 5.2 0 1 size 50 5 5 group 1
gen zone brick p0 15.2 0 0 p1 15.4 0 0 p2 15.2 1 0 p3 15.2 0 1 size 10 5 5 group 2
gen zone brick p0 15.4 0 0 p1 20.4 0 0 p2 15.4 1 0 p3 15.4 0 1 size 25 5 5 group 1
;plot block group
;pau
mo elas
prop bulk 3e8 shear 1.8e8
prop bulk 0.15e8 shear 0.09e8 range group 2
ini dens 2500
ini dens 1500 range group 2
;plot block dens
;pau
fix x y z range x 0.1 -.1
;计算sin（wt）
def wave
    if dytime>1./2000.0
       wave=0.0
    else
       wave=-1.0e6*sin(2*pi*1000.0*dytime)
    endif
end
;施加动态应力
;app nquiet dquiet squiet range x 20.09 20.11
;app nquiet dquiet squiet range x -.1 .1
apply sxx 1.0 hist wave range x 20.39 20.41
set hist_rep 1
hist gp xdisp 0 0.5 1 ;1
hist gp xdisp 5 0.5 1 ;2
hist gp xdisp 5.2 0.5 1 ;3
hist gp xdisp 15.2 0.5 1 ;4
hist gp xdisp 15.4 0.5 1 ;5
hist gp xdisp 20.4 0.5 1 ;6
hist dytime ;7
plot create hhh
plot add hist 1 2 3 4 5 6 v 7
plot show
solve age 1
hist write 1 2 3 4 5 6 v 7 file 剪模0.09e8.txt
save 剪模0.09e8.sav
;=========================
;剪模为1.44e8
;=========================
new
config dynamic
set dynamic dt 1.0e-5
set dynamic time 0.0
;set dynamic multi on
set dynamic damp rayleigh 0.1 10
gen zone brick p0 0 0 0 p1 5 0 0 p2 0 1 0 p3 0 0 1 size 25 5 5 group 1
gen zone brick p0 5 0 0 p1 5.2 0 0 p2 5 1 0 p3 5 0 1 size 10 5 5 group 2
gen zone brick p0 5.2 0 0 p1 15.2 0 0 p2 5.2 1 0 p3 5.2 0 1 size 50 5 5 group 1
gen zone brick p0 15.2 0 0 p1 15.4 0 0 p2 15.2 1 0 p3 15.2 0 1 size 10 5 5 group 2
gen zone brick p0 15.4 0 0 p1 20.4 0 0 p2 15.4 1 0 p3 15.4 0 1 size 25 5 5 group 1
;plot block group
;pau
mo elas
prop bulk 3e8 shear 1.8e8
prop bulk 2.4e8 shear 1.44e8 range group 2
ini dens 2500
ini dens 1500 range group 2
;plot block dens
;pau
fix x y z range x 0.1 -.1
;计算sin（wt）
def wave
    if dytime>1./2000.0
       wave=0.0
    else
       wave=-1.0e6*sin(2*pi*1000.0*dytime)
    endif
end
;施加动态应力
;app nquiet dquiet squiet range x 20.09 20.11
;app nquiet dquiet squiet range x -.1 .1
apply sxx 1.0 hist wave range x 20.39 20.41
set hist_rep 1
hist gp xdisp 0 0.5 1 ;1
hist gp xdisp 5 0.5 1 ;2
hist gp xdisp 5.2 0.5 1 ;3
hist gp xdisp 15.2 0.5 1 ;4
hist gp xdisp 15.4 0.5 1 ;5
hist gp xdisp 20.4 0.5 1 ;6
hist dytime ;7
plot create hhh
plot add hist 1 2 3 4 5 6 v 7
plot show
solve age 1
hist write 1 2 3 4 5 6 v 7 file 剪模1.44e8.txt
save 剪模1.44e8.sav