- 积分
- 0
- 注册时间
- 2008-3-17
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2009-11-13 17:20:13
|
显示全部楼层
来自 重庆沙坪坝区
new
gen zone brick size 10 10 10 &
p0 0 0 0 p1 10 0 0 p3 0 0 10
model mohr
prop bulk 1e8 shear 0.3e8 fric 35 coh 1e10 ten 1e10
set grav 0 0 -9.81
ini dens 1000
fix x range x 0
fix x range x 10
fix y range y 0
fix y range y 10
fix z range z 0
hist unbal
set hist_rep 10
solve
save 初始平衡.sav
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
model null range x 0 2 y 0 2 z 0 2
solve
save 1次开挖.sav
restore 初始平衡.sav
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
model null range x 0 2 y 0 2 z 0 2
def place_cables
loop m(1,2)
caseof m
case 1
x_d1=0.5
x_d2=0.5
x_d3=0.5
y_d1=0.5
y_d2=0.5
y_d3=0.5
z_d1=1.5
z_d2=2.0
z_d3=3.5
case 2
x_d1=0.5
x_d2=0.5
x_d3=0.5
y_d1=1.5
y_d2=1.5
y_d3=1.5
z_d1=1.5
z_d2=2.0
z_d3=3.5
endcase
loop n(1,2)
x_d1=x_d1+1
x_d2=x_d2+1
x_d3=x_d3+1
y_d1=y_d1
y_d2=y_d2
y_d3=y_d3
z_d1=z_d1
z_d2=z_d2
z_d3=z_d3
command
sel cable beg x_d1 y_d1 z_d1 end x_d2 y_d2 z_d2 nseg 2
sel cable beg x_d2 y_d2 z_d2 end x_d3 y_d3 z_d3 nseg 4
sel cable prop emod 2e10 ytension 3.1e5 xcarea 0.0004906 &
gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,2
sel cable prop emod 2e10 ytension 3.1e5 xcarea 0.0004906 &
gr_coh 10e5 gr_k 2e7 gr_per 0.0785 range cid 3,6
sel delete link range id 1
sel link id=100 1 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 1
sel cable pretension 60e3 range cid 1,6
endcommand
end_loop
end_loop
end
place_cables |
|