- 积分
- 7
- 注册时间
- 2005-8-5
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2010-3-26 11:15:17
|
显示全部楼层
来自 湖北武汉
锚杆抗拔试验flac3d程序
gen zone brick p0 -3,0,-3 p1 3,0,-3 p2 -3,10,-3 p3 -3,0,3 size 12 20 12
mod elastic
def derive
s_mod = y_mod / (2.0 * (1.0 + p_ratio))
b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))
end
set y_mod = 4.66e9 p_ratio = 0.34
derive
pro bulk b_mod she s_mod
ini dens 2300
; =======================================================
; Create a single cable and set its associated properties
;
sel cable id=1 beg 0, 0, 0 end 0 ,8, 0 nseg 40
sel cable prop Xcarea=803e-6 Emod=25e9 YTension=0.249e6 &
GR_K=0.78e9 GR_Coh=0.52e5 GR_Fric=0.0 GR_Per=2.83e-1
; =======================================================
; Fix free end of rock block and apply velocity to cable end
;
def initval
pull_vel = -1.0e-6
end
initval
;
fix y range y=(-0.01 0.01)
sel node FIX x range y=(-0.01 0.01)
sel node INI yvel=pull_vel range y=(-0.01 0.01)
; Create FISH function for monitoring total axial force built up
; in the rock per unit cable length [force_len] and
; displacement of cable-end [disp_end].
;
def force_len
;
_sum = 0.0
gpp = gp_head
loop while gpp # null
if gp_ypos(gpp) <= 0.02 then
_sum = _sum + gp_yfunbal(gpp)
end_if
gpp = gp_next(gpp)
end_loop
force_len = _sum
disp_end = step * (-1.0)*pull_vel
end
; Set up histories for monitoring model behavior
;
history nstep=10
hist id=1 unbal
hist id=10 force_len
hist id=20 sel node ydisp id=1
plot create model_view
plot add sketch
plot add axes lgray
plot add sel geom black black scale 0.005
plot add his -10 vs -20
plot set color off
plot set center (1.874e-001,2.919e-001,1.020e-001)
plot set rotation (30.00, 0.00, 30.00)
plot set distance 2.1e+001
plot set angle 22.50
plot set magnification 5.12e-001
plot show
plot create bolt_forces_displ
plot add sketch
plot add axes black
plot add sel cable force
plot add sel cable grout slip
plot set center (1.874e-001,2.919e-001,1.020e-001)
plot set rotation (30.00, 0.00, 30.00)
plot set distance 1.4560e+000
plot set angle 22.50
plot set magnification 5.12e-001
; =======================================================
; Apply velocity to achieve total displacement of 1 cm
;
cycle 30000
sav kblny-a6-1.sav
return |
评分
-
1
查看全部评分
-
|