平面应变条件下试件剪切带分布请教
模型高l0cm，宽5cm，厚0.5cm。将模型划分成单元网格0.125，0.25，0.5cm大小, 试件的前后面的y向位移被约束，以形成平面应变状态。光滑时上下端面侧向位移自由，摩擦很大时，上下端面侧向位移被约束，加载应变率为le-7m/s，侧压为le5Pa.光滑时
title
axial compression test
set logfile wumoca
set log on
gen zone brick size 5 1 10 p0 0,0,0 p1 0.05,0,0 p2 0,0.005,0 p3 0,0,0.1
model ss
prop bulk 1.19e10 shear 1.1e10
prop coh 2.72e5 fric 44 ten 2e5
prop ctab 1 ftab 2
table 1 0,2.72e5 1e-4,2e5 2e-4,1.5e5 3e-4,1.03e5 1,1.03e5
table 2 0,44 1e-4,42 2e-4,40 3e-4,38 1,38
fixy range y -0.0001 0.0001
fixy range y 0.0049 0.0051
apply sxx =-1e5range x -0.001 0.001
apply sxx =-1e5range x 0.049 0.051
ini zvel 1e-7 range z -.01 0.01
ini zvel -1e-7 range z 0.09 0.11
step 3000
hist gp zdisp 0,0,0
hist zone szz 0,0,0.05
plot his -2 vs 1
print zone strain
save wumaca.sav
摩擦很大时
title
axial compression test
set logfile moca
set log on
gen zone brick size 5 1 10 p0 0,0,0 p1 0.05,0,0 p2 0,0.005,0 p3 0,0,0.1
model ss
prop bulk 1.19e10 shear 1.1e10
prop coh 2.72e5 fric 44 ten 2e5
prop ctab 1 ftab 2
table 1 0,2.72e5 1e-4,2e5 2e-4,1.5e5 3e-4,1.03e5 1,1.03e5
table 2 0,44 1e-4,42 2e-4,40 3e-4,38 1,38
fix x y range z -0.01 0.01
fix x y range z 0.09 0.11
fixy range y -0.0001 0.0001
fixy range y 0.0049 0.0051
apply sxx =-1e5range x -0.001 0.001
apply sxx =-1e5range x 0.0049 0.0051
ini zvel 1e-7 range z -.01 0.01
ini zvel -1e-7 range z 0.09 0.11
step 3000
hist gp zdisp 0,0,0
hist zone szz 0,0,0.05
plot his -2 vs 1
print zone strain
save moca.sav
图是别人算的一个例子，上面命令流是我编的，我怎么也模拟不出如图示的剪切带，望指教。
[ 本帖最后由 wazoch 于 2006-5-30 18:34 编辑 ] 我看了一下，这种模拟应该还是很难做的，尤其是对这种试验的模拟，因为试样小，边界效应明显，然后跟你选择的本构模型关系密切，所以得不到这种结果也很正常。
如果别人也是用flac算的，那么请仔细阅读他的文章，看他的参数、模型、计算方法有什么特点，不可以想当然的就编出一个很好的结果。 楼主你的问题解决没？那个例子是王学滨教授的。我按照他的一篇文章，同样的参数，模拟不出他的效果。不知道是怎么回事
页:
[1]