Re: [請益] 有關PFC2D的FISH使用方法

看板Civil (土木)作者albert8211 (大GG)時間8年前 (2016/11/15 16:41)推噓1(1推 0噓 1→)

留言2則, 2人參與討論串3/4 (看更多)

※ 引述《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://www.ptt.cc/bbs/Civil/M.1479199286.A.B6C.html