PFC休止角（堆积角）滚筒计算方法——碎石

碎石的堆积角的关键点不是它容易碎，而是它表面的不规则轮廓。

做数值模拟首先要分析，保留你所需要的关键信息。

如果你做碎石的压缩力学之类的东西，那你主要考虑碎石的易碎性。

如果你做这种堆积角，那么碎石的表面不规则性则是你所要考虑的关键！

这里给出随时都额生成方式——clump

如图

这里用clump生成一个不规则的碎石模型，然后后面就跟ball一样，这里没有计算到最终状态，已经算了三天三夜了，伤不起。

代码如下：

new
domain extent -1 1
wall generate circle position 0 0 radius 1geometry import kuaiti.dxf
clump template create name gangti ...
    geometry kuaiti ...
    surfcalc bubblepack ratio 0.3 distance 150
clump distribute porosity 0.12             ...
                  resolution 1.0            ...
                  numbin    1               ...
                  diameter                  ... 
                  bin 1 size 0.03 0.045 ...
                        volumefrac  1.0     ...
                        template gangti       ...
                            range annulus center 0 0 radius 0 [1-0.06]...
                                        plane o 0 0 dip 0 below;clump distribute porosity 0.28 radius 0.006 0.009 range annulus center 0 0 radius 0 [1-0.009]...
                                        ;plane o 0 0 dip 0 belowclump attribute density 3e3 damp 0.7
cmat default model linear method deformability emod     100e6  kratio 1.74 property fric 0.5cycle 2000 calm 50
set gravity 9.8
solvesave sample
clump attribute displacement multiply 0
cmat add model linear method deformability emod   100e6  kratio 1.74 property fric 1.5 range contact type pebble-facetset mech age 0
[time_record=mech.age]
[baocunpinlv=2]
def savefile
   
    if mech.age-time_record > baocunpinlv then
        filename=string.build("jieguo%1",count)
        command
            save @filename
           
        endcommand
        time_record=mech.age
        count +=1
    endif
   
end
set fish callback -1.0 @savefilewall attribute centrotation  0 0 spin 50
solve time 70

这里将稻谷和此处的轮廓图附上，我这里是用adobe illustrate画的二维图，然后导出的dxf

链接：https://pan.baidu.com/s/1oRelDk669LZLtLBeKTSANA
提取码：43il
复制这段内容后打开百度网盘手机App，操作更方便哦