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

※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.120.102.105 ※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Civil/M.1479199286.A.B6C.html