config cppudm
g 10 10
gen 0 0 0 .6 .6 .6 .6 0
model load DuncanChang.dll
model duncan
prop cohesion 11.7e3 friction 34 fricDel 0.0 ratiofail 0.59 ke 140 ne 0.45 kb 150 mb 0.39 kur 844.8
prop den 1750
set grav=10
fix y j=1

def sigma3
  sigma3=-90e4
  sigma3_2=90e4
end
sigma3

app press sigma3_2 i 1     ;若此處用nstress則后面的應(yīng)力要用sigma3，即為負(fù)號(hào)。
app press sigma3_2 i 10    
app press sigma3_2 j 10    ;

ini szz sigma3
ini syy sigma3
ini sxx sigma3
sol
ini xd 0 yd 0
ini xv 0 yv 0
;hist id=2 gp xd 0 0 0.6
;hist id=3 gp yd 0 0 0.6
hist  syy i=1 j=1
hist  sxx i=1 j=1
;tab 1 name loads
;tab 2 name ev_ea
;第1次加卸載
def load1
  loop n(1,2000)
    zss_load=sigma3_2+ float(n)*1e3      ;若下面用nstress則此處應(yīng)為sigma3-float(n)*1e3
    z_dis = -1*ydisp(5,10) / 0.6         ;后面的也應(yīng)改為z_load = (sigma3-zss_load)
    z_load = (zss_load-sigma3_2)
    command
      app press zss_load j=10
      solve
      tab 1 z_dis z_load
    end_command
  end_loop
end
load1

def unload1
  loop m(1,1500)
    zss_load=sigma3_2+2000e3 - float(m)*1e3
    z_dis = -1*ydisp(5,10) / 0.6
    z_load = (zss_load-sigma3_2)
    command
      app press zss_load j=10
      solve
      tab 1 z_dis z_load
    endcommand
  endloop
end
unload1
;==============

;第2次加卸載
def load2
  loop n(1,2500)
    zss_load=sigma3_2+500e3 + float(n)*1e3
    z_dis = -1*ydisp(5,10) / 0.6
    z_load = (zss_load-sigma3_2)
    command
      app press zss_load j=10
      solve
      tab 1 z_dis z_load   
    end_command
  end_loop
end
load2

def unload2
  loop m(1,2500)
    zss_load=sigma3_2+3000e3 - float(m)*1e3
    z_dis = -1*ydisp(5,10) / 0.6
    z_load = (zss_load-sigma3_2)
    command
      app press zss_load j=10
      solve
      tab 1 z_dis z_load
    endcommand
  endloop
end
unload2
;==============

;第3次加卸載
def load3
  loop n(1,3500)
    zss_load=sigma3_2+500e3 + float(n)*1e3
    z_dis = -1*ydisp(5,10) / 0.6
    z_load = (zss_load-sigma3_2)
    command
      app press zss_load j=10
      solve
      tab 1 z_dis z_load   
    end_command
  end_loop
end
load3

def unload3
  loop m(1,3500)
    zss_load=sigma3_2+4000e3 - float(m)*1e3
    z_dis = -1*ydisp(5,10) / 0.6
    z_load = (zss_load-sigma3_2)
    command
      app press zss_load j=10
      solve
      tab 1 z_dis z_load
    endcommand
  endloop
end
unload3
;==============

;第4次加卸載
def load4
  loop n(1,4500)
    zss_load=sigma3_2+500e3 + float(n)*1e3
    z_dis = -1*ydisp(5,10) / 0.6
    z_load = (zss_load-sigma3_2)
    command
      app press zss_load j=10
      solve
      tab 1 z_dis z_load   
    end_command
  end_loop
end
load4

def unload4
  loop m(1,4500)
    zss_load=sigma3_2+5000e3 - float(m)*1e3
    z_dis = -1*ydisp(5,10) / 0.6
    z_load = (zss_load-sigma3_2)
    command
      app press zss_load j=10
      solve
      tab 1 z_dis z_load
    endcommand
  endloop
end
unload4

;==============

;第5次加卸載
def load5
  loop n(1,5500)
    zss_load=sigma3_2+500e3 + float(n)*1e3
    z_dis = -1*ydisp(5,10) / 0.6
    z_load = (zss_load-sigma3_2)
    command
      app press zss_load j=10
      solve
      tab 1 z_dis z_load   
    end_command
  end_loop
end
load5