请教下FLAC2D中如何用interface模拟断层？感激涕零！

zhangjian321

小弟初学FLAC2D 想向大家请教下如何用interface模拟断层 最好有实例 在此谢过了先！！！！

sinlk

有吗？我也需要啊 小弟邮箱wanggangqiuzhi@163.com

feigle

用interface
命令 就行了

xiaoyehan

楼上太废话了

franco_chao

断层，是否可以直接建单剪实验模拟？？

我弄过一个单剪实验的模型。

franco_chao

单剪实验，速度加载，应力控制，一起分享。

;-------------- Contrainte impose ------------> OK
new

gen zone brick size 1 1 1 rat 1 1 1 &
p0 0 0.05 0 p1 1 0.05 0 p2 0 0.1 0 p3 0 0.05 0.05

group bas range y 0.05 0.1

int 1 fa ra y 0.05

gen zone brick size 1 1 1 rat 1 1 1 &
p0 -2 0 0 p1 -1 0 0 p2 -2 0.05 0 p3 -2 0 0.05

group haut range y 0 0.05

ini x add 2 ra gr haut

def propd
    Eyh=20e7
     Kh=Eyh/(3*(1-2*nuh))
    Gh=Eyh/(2*(1+nuh))
    Eyb=10e7
    nu=0.3
    Kb=Eyb/(3*(1-2*nub))
    Gb=Eyb/(2*(1+nub))
    knor=3850000000
    ktan=500000000
       sn0=-1e5
end
propd

model elas ra gr haut
prop bulk Kh shear Gh
model elas ra gr bas
prop bulk Kb shear Gb  

int 1 prop kn knor ks ktan tension 1e10

fix y ra y 0.1
apply sxx sn0 ra x 0
apply sxx sn0 ra x 1
apply szz sn0 ra z 0
apply szz sn0 ra z 0.05
apply syy sn0 ra y 0

solve rat 1e-10

sav ei3.sav

apply remove sxx
apply remove szz
apply remove syy
free x y z
fix x ra y 0
fix x ra y 0.1
fix x ra gr bas y 0.05
fix y
fix z
ini xdisp 0
ini ydisp 0
ini zdisp 0
ini xvel 0
ini yvel 0
ini zvel 0

    apply xvel dep ra gr haut x 0 y 0.05
    apply xvel dep ra gr haut x 1 y 0.05

wkh1102

这是3d的啊

franco_chao

这是3d的啊
wkh1102 发表于 2010-3-11 16:10

恩，是