关于陈育民第九章导来导去法的问题讨论与分析

周玉5dc9b2cc

陈育民第九章导来导去法 命令流如下：
n
gen zone brick size 3 3 3
group 2 range x 1 2 y 1 2 z 1 2
group 1 range gr 2 not
save 1.sav
del ran group 2 not

interface 1 face range x 1 y 1 2 z 1 2
interface 1 face range x 2 y 1 2 z 1 2
interface 1 face range x 1 2 y 1 z 1 2
interface 1 face range x 1 2 y 1 2 z 1
interface 1 face range x 1 2 y 1 2 z 2
plot int yellow ske
save 2.sav

rest 1.sav
del ran group 2
expgrid 1.flac3d
rest 2.sav
impgrid 1.flac3d
plot int yellow ske

model ela
fix x y z ran z 0
ini den 2000
set grav 0 0 -10
interface 1 prop kn 20e6 ks 20e6 coh 10e3 fri 15
app nstr -200e3 ran x 0 1 y 1 2 z 3

solve运行后却出现了问题

周玉5dc9b2cc

今天在贴吧里看到同样的问题，看了一个师兄的命令流：
;导来导去法
;---------------------------------------------
n
gen zone brick size 3 3 3
group 2 range x 1 2 y 1 2 z 1 2
group 1 range gr 2 not
save 1.sav

del ran group 2 not
interface 1 face range x 1 y 1 2 z 1 2
interface 1 face range x 2 y 1 2 z 1 2
interface 1 face range x 1 2 y 1 z 1 2
interface 1 face range x 1 2 y 1 2 z 1
interface 1 face range x 1 2 y 1 2 z 2
plot int yellow ske
save 2.sav

rest 1.sav
del ran group 2
expgrid 1.flac3d
rest 2.sav
impgrid 1.flac3d
plot int yellow ske

model ela
fix x y z ran z 0
ini den 2000
set grav 0 0 -10
prop shear 30e8 bulk 20e8
interface 1 prop kn 20e6 ks 20e6 coh 10e3 fri 15
app nstr -200e3 ran x 0 1 y 1 2 z 3
hist gp zdisp 1 1 1
和教材中的唯一区别就是倒数第四行prop shear 30e8 bulk 20e8给模型进行了材料参数赋值，结果能运行出来
这是什么原因造成的？现在还不清楚，望各路大神
一起讨论共同交流
solve

周玉5dc9b2cc

怎么没有人回应呢

周玉5dc9b2cc

;导来导去法
;---------------------------------------------
n
gen zone brick size 3 3 3
group 2 range x 1 2 y 1 2 z 1 2
group 1 range gr 2 not
save 1.sav

del ran group 2 not
interface 1 face range x 1 y 1 2 z 1 2
interface 1 face range x 2 y 1 2 z 1 2
interface 1 face range x 1 2 y 1 z 1 2
interface 1 face range x 1 2 y 2 z 1 2
interface 1 face range x 1 2 y 1 2 z 1
interface 1 face range x 1 2 y 1 2 z 2
plot int yellow ske
save 2.sav

rest 1.sav
del ran group 2
expgrid 1.flac3d
rest 2.sav
impgrid 1.flac3d
plot int yellow ske

model ela
fix x y z ran z 0
ini den 2000
set grav 0 0 -10
prop shear 30e8 bulk 20e8
interface 1 prop kn 20e6 ks 20e6 coh 10e3 fri 15
app nstr -200e3 ran x 0 1 y 1 2 z 3
hist gp zdisp 1 1 1终于发现了问题所在，网上下载的命令流里接触面的建立少了一个面，红颜色标注的，大家要注意了哦
solve