有關PFC2D的FISH使用方法 - 土木

Table of Contents

我簡化一些東西之後剩這樣
new
set random

def aa
bp=ball_head
loop while bp#null
b_rad(bp)=F
endloop
end
aa

wall id=1 node 0,0 10,0
wall id=2 node 10,0 10,10
gen id=1,1 rad 0.2 0.2 x 0.2 0.2 y 0.2 0.2
wall id=1 kn=1e8 ks=1e8
wall id=2 kn=1e8 ks=1e8

property density 1000 kn 1e8 ks 1e8
set grav 0 -9.81
prop fric 0.5

plot show
plot add ball y
plot add wall b

property xforce=F range id 1,1
cyc 1000
print ball position

情況還是跟之前一樣,可以跑,但是就是一顆球站在那裏不動,照理說
我是期望他會受到本身自己粒徑大小的力,也就是我設的球半徑0.2,
問題在xforce=F(不動),xforce=0.2(就會動了),再請各位高手指教了。

--

All Comments

Gary avatarGary2016-11-18
把aa放到倒數第4行再執行