有大神帮忙翻译下这些指令的意思吗 锚杆的
Simulation of pull-test for grouted reinforcement
;confinement 4 MPa
; File: Pull03.dat
;
; ==================================================================
new
title
Pull Test for Grouted Cable Anchor - Confinement 4 MPa, fric=30
; =======================================================
; Create a single rock block and set its material properties.
;
gen zone brick size 4 7 4
mod el
ini x mul .1
ini y mul .1
ini z mul .1
pro bulk 5e9 she 3e9
;model mohr
;property density 2320 bulk 1.98e9 shear 1.19e9 cohesion 7e4 friction 24
;
apply nstress -4.0e6 range x=-0.01 0.01
apply nstress -4.0e6 range x=0.39 0.41
apply nstress -4.0e6 range z=-0.01 0.01
apply nstress -4.0e6 range z=0.39 0.41
cy 100000
save pull031.sav
;
; =======================================================
; Create a single cable and set its associated properties
;
sel cable id=1 begin=(0.2, 0.0, 0.2) end=(0.2, 0.5, 0.2) nseg=10
sel cable prop Xcarea=362e-6 Emod=197.2e9 YTension=0.464e6 &
GR_K=2.24e7 GR_Coh=3.5e5 GR_Fric=60.0 GR_Per=15.7e-2
; =======================================================
; 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.05 then
_sum = _sum + gp_yfunbal(gpp)
end_if
gpp = gp_next(gpp)
end_loop
force_len = _sum / 0.5
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
hist id=100 sel cableSEL force cid=1
hist id=200 sel cableSEL force cid=3
hist id=300 sel cableSEL force cid=5
hist id=400 sel cableSEL force cid=7
hist id=500 sel cableSEL force cid=9
hist id=600 sel cableSEL force cid=10
plo show
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 xmax 5e-2 ymax 3.8e5
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 1.4560e+000
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 2.0 cm
;
cycle 45000
; =======================================================
sav pull032.sav
; =======================================================
return