这是单轴压缩试验的命令
n
gen zone cyl p0 0 0 0 p1 0.025 0 0 p2 0 0.1 0 p3 0 0 0.025 size 4 8 4
gen zone reflect norm 1,0,0
gen zone reflect norm 0,0,1
model mohr
pro den 2500 young 1.12e10 poisson 0.22 co 1.2e6 fric 50 ten 1.07e6
fix x y z range y -0.0001 0.0001
fix x y z range y 0.099 0.101
ini yvel 1e-5 range y -0.0001 0.0001
ini yvel -1e-5 range y 0.099 0.101
def ax_str
str = 0
sabcd=pi*0.025*0.025
pnt = gp_head
loop while pnt # null
if gp_ypos(pnt) <0.01 then _____________这行的含义是什么？ 0.01根据什么取得？
str = str + gp_yfunbal(pnt) _____________ 不平衡力的累加代表轴向力？
endif
pnt = gp_next(pnt)
endloop
ax_str = str / sabcd
end
要是想做三轴卸荷试验，怎么能把轴压增加的同时降低围压？
