请教下FLAC2D中如何用interface模拟断层?感激涕零!
小弟初学FLAC2D 想向大家请教下如何用interface模拟断层 最好有实例 在此谢过了先!!!! 有吗?我也需要啊 小弟邮箱wanggangqiuzhi@163.com 用interface命令 就行了 楼上太废话了 断层,是否可以直接建单剪实验模拟??
我弄过一个单剪实验的模型。 本帖最后由 franco_chao 于 2010-3-10 23:13 编辑
单剪实验,速度加载,应力控制,一起分享。
;-------------- 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 这是3d的啊 这是3d的啊
wkh1102 发表于 2010-3-11 16:10 http://forum.simwe.com/images/common/back.gif
恩,是
页:
[1]