PFC2D-FISH
;fname: fishb5.DAT
new
wall id 5 kn 1e8 nodes -5,0 5,0
ball id 1 x 0 y 1.1 rad 1
prop dens 1000 kn 1e8
set grav 0 -10
set dt max=2e-3
def down_force
down_force = b_yfob(ball_head)
end
hist down_force
cyc 500
plot his 1
**缺少 set dt max=2e-3,则
;fname: fishb6.DAT
new
def box_geometry
n_stiff = 2e8
s_stiff = 1e8
xx = 5.0
yy = 3.0
end
box_geometry
macro wall1nodes 'nodes ( 0, 0) (xx, 0)'
macro wall2nodes 'nodes (xx, 0) (xx, yy)'
macro wall3nodes 'nodes xx yy 0 yy'
macro wall4nodes 'nodes 0 yy 0 0'
wall id=1 ks=s_stiff kn=n_stiff wall1nodes
wall id=2 ks=s_stiff kn=n_stiff wall2nodes
wall id=3 ks=s_stiff kn=n_stiff wall3nodes
wall id=4 ks=s_stiff kn=n_stiff wall4nodes
plot wall red id on
print wall prop
print macro
**将墙的属性设于文件开头可减少数据修
改和犯错。
评论0
最新资源