- 积分
- 11
- 注册时间
- 2005-5-23
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2009-6-26 15:08:41
|
显示全部楼层
来自 安徽合肥
此为计算时的命令流,望大家指点
gen zone radtun p0 0,0,0 p1 2,0,0 p2 0,10,0 p3 0,0,2 size 6 30 6 8 dim 1 1 1 1 ratio 1 1 1 1.1 fill
gen zone reflect normal 1 0 0 ori 0 0 0
gen zone reflect normal 0 0 1 ori 0 0 0
;group trench range y -0.1,0.1 x -0.2,0.2 z -0.2,0.2
mo mohr
pro bulk=41.67e9 shear=19.23e9 fric 28 coh 2e6 ten 2.5e6
ini dens 2600
fix x range x -2.1,-1.9
fix x range x 1.9 2.1
fix y range y 9.9 10.1
fix z range z -2.1 -1.9
fix z range z 1.9 2.1
;apply nstress=-3.6e6 range z 1.9 2.1
;apply nstress=-3.6e6 range z -2.1 -1.9
;apply nstress=-3.6e6 range x -2.1 -1.9
;apply nstress=-3.6e6 range x 1.9 2.1
;apply nstress=-18e6 range y 9.9 10.1
;apply nstress=-18e6 range x -2.1 -0.2 y -0.1 0.1
;apply nstress=-18e6 range x 0.2 2.1 y -0.1 0.1
;apply nstress=-18e6 range x -0.2,0.2 y -0.1 0.1 z 0.2 2.1
;apply nstress=-18e6 range x -0.2,0.2 y -0.1 0.1 z -2.1 -0.2
;step 7000
;sav origin0.2.sav
sel cable id=1 beg 0,0,0 end 0,4,0 nseg 16
sel cable id=1 beg 0,4,0 end 0,7,0 nseg 12
sel cable prop emod 2e11 ytension 350e3 xcarea 0.000254 &
gr_coh 0 gr_k 0 gr_per 0.085 range cid 1,16
sel cable prop emod 2e11 ytension 350e3 xcarea 0.000254 &
gr_coh 5e5 gr_k 2e8 gr_f 25 gr_per 0.085 range cid 17,28
sel liner range y=-.01, .01 x=-0.2,0.2 z=-0.2,0.2
sel liner PROP iso=(2e11, 0.15) thick=0.1 ; concrete
sel liner PROP cs_nk=8e8 cs_sk=8e8 cs_ncut=8e8 cs_scoh=8e8 cs_scohres=8e6 cs_sfric=68
sel delete link range id 1
sel link id=1001 1 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 1001
sel cable id=1 pretension 100e3 range cid 1,16
step 2000
hist unb
pl hist 1
sav 20.sav |
|