沥青混合料单轴压缩试验

zhaoqiqi

求助，请大神们看看下面的单轴压缩代码是否有问题？
restore C-3.sav
prop density 1000 friction 0.8
prop kn=2e8 n_bond=1000
prop ks=4e7 s_bond=800
prop pb_kn 2e8 pb_nstrength 1e5
prop pb_ks 4e7 pb_nstrength 2e4
prop pb_radius 0.0015
prop density 2700 friction 0.5 range group agg
prop kn 8e8 n_bond=10000 range group agg
prop ks 2e8 s_bond=6000 range group agg
;=====================================================================
def get_wall_ss
  command
    wall id=1 kn 1e+11 ks 2e+10
    wall id=3 kn 1e+11 ks 2e+10
  end_command
end
;=====================================================
get_wall_ss
delete wall 2,4
;=====================================================
def get_ss
  wadd = find_wall(3)
  weyy = - w_y(wadd)
  wsyy = w_yfob(wadd) / 0.15 / 1000000
end
;=================================================================
def p_ss
  while_stepping
  get_ss
  yingli = wsyy
  yingbian = weyy
end
;=================================================================
def c_t
    c_t = tdel
end
;=================================================================
def wall_loading
  command
    wall id 3 yvel -0.02;==============记载速度自己设置吧
  end_command
end
;=================================================================
def _ct
  _ct = time
end
;=================================================================
_ct
cyc 1
set dt 0.00001
p_ss
hist yingli
hist yingbian
hist dump 1,2
plot add hist 1 vs 2
wall_loading
hist write 1,2 file date.dat
cycle 80000
plot show

zhaoqiqi

谁有PFC3.10   二维的单轴压缩试验代码呀？