首页/文章/ 详情

PFC休止角（堆积角）滚筒计算——稻谷

lobby

4年前浏览5176

这个算例农业的同学可以看过来了。

这里介绍一下PFC实现非圆形颗粒的方法——clump

用户可以引入图形轮廓文件，然后生成一个clump模板。

之后便可以像生成ball一样在算例里面生成clump了。

如图：

用户可以指定这个clump的参数，如果计算速度足够的话，这个轮廓可以显示的更加好。

下面给出变形的动图：

我这里没有完全计算完，这个已经算了三天三夜了，有条件的同学可以继续去算，下面附上代码：

new
domain extent -1 1
wall generate circle position 0 0 radius 1geometry import daogu.dxf
clump template create name daogu ...
    geometry daogu ...
    surfcalc bubblepack ratio 0.3 distance 150
clump distribute porosity 0.18             ...
                  resolution 1.0            ...
                  numbin    1               ...
                  diameter                  ...
                  bin 1 size 0.03 0.045 ...
                        volumefrac 1.0     ...
                        template daogu       ...
                            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
clump delete range y -10 -1
solve time 300