fangweino1 发表于 2010-1-25 09:14:56

各位大侠,请问我的预应力锚索中,预应力大小改变为什么对结果毫无影响啊?

各位大侠,请问我的预应力锚索中,预应力大小改变为什么对结果毫无影响啊?

jihezhilv 发表于 2010-1-25 22:49:39

请教楼主对预应力锚索的模拟是:压力分散型锚索还是拉力集中型锚索。不知楼主是否看过tg0215的帖子没有?http://forum.simwe.com/viewthread.php?tid=758239&extra=page%3D1%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D1
能否贴上你的命令流大家好交流,看看能不能帮助你。

fangweino1 发表于 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 0p3 37.5 0 5   size 4 3 4 ratio 1 1 1
gen zone brick p0 0 0 5p1 37.5 0 5p2 0 3 5p3 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-3emod = 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-3emod = 45e9 ytension = 60e4 gr_k = 2.75e7 gr_coh = 2.2e4 gr_fric = 20 gr_per = 0.314 range cid 13,16
sel delete linkrange id=1
sel link id=1 1 target zone
sel link attach xdir = rigid ydir = rigid zdir = rigid   xrdir = rigid yrdir = rigid zrdir = rigidrange 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-3emod = 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-3emod = 45e9 ytension = 60e4 gr_k = 2.75e7 gr_coh = 2.2e4 gr_fric = 20 gr_per = 0.314 range cid 27,30
sel delete linkrange id=19
sel link id=19 18 target zone
sel link attach xdir = rigid ydir = rigid zdir = rigid   xrdir = rigid yrdir = rigid zrdir = rigidrange 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-3emod = 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-3emod = 45e9 ytension = 60e4 gr_k = 2.75e7 gr_coh = 2.2e4 gr_fric = 20 gr_per = 0.314 range cid 39,41
sel delete linkrange id=35
sel link id=35 33 target zone
sel link attach xdir = rigid ydir = rigid zdir = rigid   xrdir = rigid yrdir = rigid zrdir = rigidrange 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

fangweino1 发表于 2010-1-31 09:30:42

哎~!
还是没人回答。。。。。:sleepy:

achilles918 发表于 2011-7-11 10:14:20

非关联流动法则求解命令是不也是 associated?
页: [1]
查看完整版本: 各位大侠,请问我的预应力锚索中,预应力大小改变为什么对结果毫无影响啊?