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

※ 引述《Jate (桔子)》之銘言：
: ※ 引述《albert8211 (大GG)》之銘言：
: : 小弟第一次發文，懇請各位高手幫忙。
: : 最近再用PFC2D 4.0版本，要用此軟體解決土壤顆粒受滲流力作用的模擬。
: : 因為滲流力設定時需對應各顆粒的粒徑。
: : 例如：1號顆粒受到1號顆粒粒徑的滲流力(F=i*rw*V)依此類推......。
: : 現在想到的是用指令b_rad(bp)，但是就是一直弄不出來，可能是指令輸入邏輯問題，
: : 或是其他問題，希望參考各位高手的建議或是想法。感激不盡。
: 沒用過PFC的滲流模組, 但我記得設定好流場他好像會自己算?
: def aa
: bp = ball_head
: loop while bp#null
: F = i*rw*V <-純量, 想辦法把它拆成x跟y方向分量Fx跟Fy
: b_xfap(bp) = Fx <-滲流力x方向分量
: b_yfap(bp) = Fy <-滲流力y方向分量
: bp = b_next(bp)
: endloop
: end
: 我猜你要的應該是類似這種的東西
: 很久沒用了, 有錯的話不要怪我嘿

我現在先以最簡單的一個顆粒做測試，力的方向只採用x方向，預期是顆粒受到自己
粒徑大小的力。

new
set random

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

wall id=1 node 0,0 10,0
wall id=2 node 10,0 10,10
wall id=3 node 10,10 0,10
wall id=4 node 0,10 0,0
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
wall id=3 kn=1e8 ks=1e8
wall id=4 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
plot add axes r
cyc 1000
print ball position
property xforce=F range id 1,1
以上是我整套語法，可是跑出來是可以，但是就是沒有xforce，邏輯一定有出問題，
可以請高手幫我指正一下，或是解說一下方向。感恩。

--