lpfvsfl 发表于 2010-3-26 10:18:47

谁有锚杆拉拔试验的命令流，共享下啊

谁有锚杆拉拔试验的命令流，共享下啊，不甚感激

doctorluo 发表于 2010-3-26 11:15:17

锚杆抗拔试验flac3d程序
gen zone brick p0 -3,0,-3   p1 3,0,-3p2 -3,10,-3p3 -3,0,3size 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 bulkb_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 FIXx            range y=(-0.01 0.01)
sel node INIyvel=pull_velrange y=(-0.01 0.01)

; Create FISH function for monitoring total axial force built up
; in the rock per unit cable length and
; displacement of cable-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=20sel 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 magnification5.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 magnification5.12e-001
; =======================================================
; Apply velocity to achieve total displacement of 1 cm
;
cycle 30000
sav kblny-a6-1.sav
return

lpfvsfl 发表于 2010-3-26 13:23:47

谢谢您了，感激的无以言表 2# doctorluo

abblcn 发表于 2010-5-25 16:31:59

非常感谢，正在算拉拔试验。

struggle222 发表于 2010-9-15 16:52:51

很有用，学习了，谢谢。

liuzhiqiang06 发表于 2010-9-15 17:11:25

学习了，谢谢楼主

123456-123456 发表于 2010-9-20 17:20:19

很好的资料，谢谢楼主！

liuniankd 发表于 2012-7-7 20:46:58

多谢楼主啊！

猪猪侠88 发表于 2013-4-14 19:54:16

研究下 谢谢

drillcore 发表于 2014-4-24 22:28:38

感谢楼主的无私奉献~

zbbz110 发表于 2014-5-5 20:38:05

谢谢楼主

chenpf09 发表于 2014-5-7 21:22:18

果断收藏~~~~~~~~~~~~~

shimeng19890710 发表于 2014-7-24 00:37:00

很好的帖子 ，尤其谢谢二楼

sunlang 发表于 2016-1-3 21:42:37

这个命令有问题吧

无心人呜呜呜 发表于 2017-6-28 09:44:22

好东西！谢谢！

刘小波1 发表于 2018-5-16 19:26:30

很好的资料，十分感谢

ma5629075 发表于 2018-10-15 09:24:16

doctorluo 发表于 2010-3-26 11:15
锚杆抗拔试验flac3d程序
gen zone brick p0 -3,0,-3   p1 3,0,-3p2 -3,10,-3p3 -3,0,3size 12 20 12
...

您好，我想问你点问题可以么？
页: [1]
查看完整版本: 谁有锚杆拉拔试验的命令流，共享下啊