用户定义接触模型(UDM)

sjzdh

谢谢提供分享

w511967305

多谢楼主分享！

wxz

谢谢分享

hhhm258

非常感谢

王军牛

      
视频文件7M，转成gif的33M了。更传不上去了，还是截图吧。呵呵

源代码如下：这是最初的版本，后面有过几次修改，不过这个基本上能反映这个问题啦。呵呵
new
set random

plot create tunnel
plot set back white
plot add wall black ball yellow

def initial_parameters
    w_stiff=1e8
    b_stiff=1e7
end
initial_parameters

wall id 10 type circle center 0 0 rad 4 kn=w_stiff ks=w_stiff
wall id 1 nodes (0,-2.5) (0,-4) kn=w_stiff ks=w_stiff
wall id 2 nodes (2.5,0) (4,0)  kn=w_stiff ks=w_stiff
wall id 3 nodes (0,2.5) (0,4)  kn=w_stiff ks=w_stiff
wall id 4 nodes (-2.5,0) (-4,0)  kn=w_stiff ks=w_stiff
wall id 5 nodes (1,0) (-1,0)  kn=w_stiff ks=w_stiff
wall id 6 nodes (0,-1) (0,1)  kn=w_stiff ks=w_stiff

plot show
gen x 0 2 y -2.5 0 rad 0.06 0.1 id 1,50

prop density=2650 kn=b_stiff ks=b_stiff
set gravity 0 -9.81
cyc 20000

wall prop spin 1.1 x 0 y 0

set plot avi size 640 480
movie avi_open file tunnel_movie.avi
movie step 10000 1 file tunnel_movie.avi
set dt 1e-4
cyc 100000
movie avi_close file tunnel_movie.avi

whx1002102308

谢谢楼主，学习学习