- 积分
- 0
- 注册时间
- 2008-5-5
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2010-1-26 13:20:30
|
显示全部楼层
来自 日本
本帖最后由 fangweino1 于 2010-1-26 13:27 编辑
楼上兄弟,你好!~
tg0215的帖子我看过,
我的建模方式正是按照所推荐的第一种方式来做的。命令流如下,麻烦您有空再帮忙看看:
new
;;;;;;建立网格(倾角1:1.5,边坡前缘距坡脚7.5m,边坡后缘距坡顶10m,前、后缘相距40m,深5m,坡体左右截面间距10m,下卧基岩坡度为1:2.5)
gen zone brick p0 0 0 0 p1 37.5 0 0 p2 0 3 0 p3 0 0 5 size 20 3 4 ratio 1 1 1
gen zone brick p0 37.5 0 0 p1 45 0 0 p2 37.5 3 0 p3 37.5 0 5 size 4 3 4 ratio 1 1 1
gen zone brick p0 0 0 5 p1 37.5 0 5 p2 0 3 5 p3 0 0 20 p4 37.5 3 5 p5 0 3 20 p6 15 0 20 p7 15 3 20 size 20 3 12 ratio 1 1 1
;;;;plo block model id on
;;;分组
group hnt red range plane above ori 5 0 16 nor 2 0 5
group bedrock green range group hnt not
;;;;;;设置边界条件
fix x y z range z -0.1 0.1
fix x range x -0.1 0.1
fix x range x 44.9 45.1
fix y range y -0.1 0.1
fix y range y 2.9 3.1
;;;;;;;生成初始应力
m e
prop bulk 3e7 shear 1e7
prop density 1970
set grav 0 0 -9.81
solve
save mg1.sav
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
;;;;;;;;安全系数求解(采用关联流动法则:fi=dila)
m m
ini dens 1970 range group hnt
prop bulk 1.15E7 shear 5.31E6 coh 20000 fri 12.1 ten 1000 dila 0 range group hnt
ini dens 2000 range group bedrock
prop bulk 22.6E9 shear 11.1E9 coh 6.72E6 fri 42 ten 1.58E6 dila 10.0 range group bedrock
sel cable id=1 begin=(20.63,1.5,16.24) end=(10.11,1.5,13.03) nseg=12
sel cable id=1 begin=(10.11,1.5,13.03) end=(5.33,1.5,11.57) nseg=4
sel cable id=1 prop xcarea = 1.256e-3 emod = 45e9 ytension = 60e4 gr_k = 0 gr_coh = 0 gr_fric 0 gr_per = 0.314 range cid 1,12
sel cable id=1 prop xcarea = 1.256e-3 emod = 45e9 ytension = 60e4 gr_k = 2.75e7 gr_coh = 2.2e4 gr_fric = 20 gr_per = 0.314 range cid 13,16
sel delete link range id=1
sel link id=1 1 target zone
sel link attach xdir = rigid ydir = rigid zdir = rigid xrdir = rigid yrdir = rigid zrdir = rigid range id 1
sel cable id=1 pretension 1e5 range cid 1 12
sel cable id=2 begin=(26.25,1.5,12.5) end=(15.73,1.5,9.28) nseg=10
sel cable id=2 begin=(15.73,1.5,9.28) end=(10.95,1.5,7.82) nseg=4
sel cable id=2 prop xcarea = 1.256e-3 emod = 45e9 ytension = 60e4 gr_k = 0 gr_coh = 0 gr_fric 0 gr_per = 0.314 range cid 17,26
sel cable id=2 prop xcarea = 1.256e-3 emod = 45e9 ytension = 60e4 gr_k = 2.75e7 gr_coh = 2.2e4 gr_fric = 20 gr_per = 0.314 range cid 27,30
sel delete link range id=19
sel link id=19 18 target zone
sel link attach xdir = rigid ydir = rigid zdir = rigid xrdir = rigid yrdir = rigid zrdir = rigid range id 19
sel cable id=2 pretension 1e5 range cid 17 26
sel cable id=3 begin=(31.88,1.5,8.74) end=(21.36,1.5,5.53) nseg=8
sel cable id=3 begin=(21.36,1.5,5.53) end=(16.58,1.5,4.07) nseg=3
sel cable id=3 prop xcarea = 1.256e-3 emod = 45e9 ytension = 60e4 gr_k = 0 gr_coh = 0 gr_fric 0 gr_per = 0.314 range cid 31,38
sel cable id=3 prop xcarea = 1.256e-3 emod = 45e9 ytension = 60e4 gr_k = 2.75e7 gr_coh = 2.2e4 gr_fric = 20 gr_per = 0.314 range cid 39,41
sel delete link range id=35
sel link id=35 33 target zone
sel link attach xdir = rigid ydir = rigid zdir = rigid xrdir = rigid yrdir = rigid zrdir = rigid range id 35
sel cable id=3 pretension 1e5 range cid 31 38
plo sk sel geo id on cid on node on scale 0.005
solve fos file mg2.sav associated ;;;;;;Fs=1.37
plot fos con ssi outline on vel red ;;;;;;绘制安全系数,剪应变增量及速度矢量图
;;;;;;;;;;监测采样
hist id=6 gp disp 20.63,1.5,16.24
hist id=7 gp disp 26.25,1.5,12.5
hist id=8 gp disp 31.88,1.5,8.74
plo hist 6 7 8
solve
save mg3.sav |
|