掌子面支护力逐级卸载问题

zjfoeg

以下是我用逐级卸载的方法在掌子面施加支护力的命令流，得出来的结果是当lumda减小到后面几个的时候（支护力小的时候），掌子面的位移就不怎么变了，是不是卸载的命令出了问题啊？正常应该是在支护降到比较小的时候，位移越来越大才对啊.


def relax_load
p_gp=gp_near(0,20,0)
lumda=1.0
o_press=1.5e7
o_grad=grad_syx
loop n (1,10)
    t=float(n)*5.0e4
    lumda=1.0-float(n-1)*0.1
    app_load=o_press*lumda
command
    app nstress @app_load grad 0 0 @o_grad range kwm
set fl off
    set mech rat 6.1e-5
    solve
set fluid on
    solve age @t
end_command
   xtable(1,n)=lumda
   ytable(1,n)=gp_ydisp(p_gp)

   endloop
end