- 积分
- 1
- 注册时间
- 2002-6-16
- 仿真币
-
- 最后登录
- 1970-1-1
|
模型高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
fix y range y -0.0001 0.0001
fix y range y 0.0049 0.0051
apply sxx =-1e5 range x -0.001 0.001
apply sxx =-1e5 range 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
fix y range y -0.0001 0.0001
fix y range y 0.0049 0.0051
apply sxx =-1e5 range x -0.001 0.001
apply sxx =-1e5 range 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 编辑 ] |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|