※ 引述《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，逻辑一定有出问题，
可以请高手帮我指正一下，或是解说一下方向。感恩。