请教高手：格栅拉拔试验中颗粒力与位移的关系

S160700105

  请教高手如何编程确定格栅拉拔试验中颗粒力与位移的关系，土工格栅的细观参数则是通过土工格栅拉伸
试验的颗粒流模拟来确定，模拟格栅长为100 mm，最后一个颗粒固定，对第1 个颗粒施加20 mm/min 的恒定速率，即对模拟格栅进行拉伸模拟，记录第1 个颗粒力与位移的关系，并与试验结果比较。
下面是我编的一个程序，但是得不到力与位移的关系曲线，还请高手指导。
  new
res Lb_model1.sav
;===========定义参数======================
def parameter_2
      width=500
      _mvel=-0.0020      ;加载速度      
    _ncheck=20
end
;======对土工格栅最左边的颗粒施加恒定的速率=========
def accel_vel

     command
       fix  x  y spin  range id=62
     end_command
     
     loop while 1#0

        if abs(x_dif)>_ncheck then
  
            exit
         end_if

        command
           prop xforce=_mvel range id=1
           cycle 100

       end_command
     

    end_loop

end
;============获得球地址===================
def find_bp1
  bp1=find_ball(1)
  cp=contact_head
   
end
;=======================================
parameter_2
find_bp1
x_dif
;============设置记录参量===================
history id=1 wexx
history id=2 s_for

plot create s_for-wexx
plot add hist 2 vs 1

accel_vel
save LbJg_model.sav
cy 11111111

S160700105

格栅是采用平行连接模型建立的。

ccb137731

楼主，程序不全呀，看着怪怪的。建议楼主把问题细化到某一点上。获得力与位移的关系曲线不难。

S160700105

这是前面生成格栅的程序，加上这部分就是全部了。
new
set echo on
set disk on
set random

macro zero 'ini xvel 0 yvel 0 spin 0'

;=========================
def parameter_1

   g_width=0.1              ;格栅宽度m
   g_dens=1000              ;格栅密度
   g_rbar=2e-3                  ;格栅颗粒半径

   gpb_ks=6e9              ;格栅切向粘结刚度Pa/m
   gpb_kn=6e9              ;格栅法向粘结刚度Pa/m
   g_kn=1e10                  ;格栅法向刚度
   g_ks=1e10                  ;格栅剪切刚度
   gpb_ns=1e10             ;格栅法向粘结强度Pa
   gpb_ss=1e10              ;格栅切向粘结强度Pa
   gpb_rad=1.0
   g_fric=0.5
   _mvel=-0.15e4      ;加载力   

end
;=======显示设置=============
plot create cantilever
plot set title text 'geshan_laba'
plot add ball id=on yellow
plot add pbond blue
plot add disp black
;=======生成格栅=============
def gen_geshan
   
    g_num = int(g_width/(2*g_rbar))   
    g_id=1
    g_x=0.0
    g_y=0.0
   
    loop n (1, g_num)
      command
        BALL id=g_id  x=g_x y=g_y rad=g_rbar
      end_command
        g_id =g_id + 1
        g_x = g_x + 2*g_rbar
    end_loop
   
   ii=out(string(g_num)+'gs_particles were created!')
   

end

liangpingf

history id=1 wexx
history id=2 s_for
你的wexx和s_for各表示什么意思？不是要测球颗粒1的力和位移的关系吗？

guoyedepaofan

楼主'我想请问您一个其它的问题'你的这个平行粘结强度和刚度参数是如何得出来的'请指教'谢谢

geotechnik

楼主，你的速度可能不对啊！
20 mm/min
换算成m/s是0.020/60=0.000333333333。
另外，拉拔试验中，一端受力，另一端是自由的吧？

lz12121212

学习:)

whiteboy

过来学习学习

yingkong2007

S160700105 发表于 2011-12-19 10:07
这是前面生成格栅的程序，加上这部分就是全部了。
new
set echo on

你做出来没  我感觉你想的过于简单了吧