【PFC6.0】泰森多边形区域划分及颗粒填充

model new

model domain extent -10 10 



wall generate box -5 5

ball distribute porosity 0.1 radius 0.8 1.5 box -5 5cmat default type ball-facet model linear method deform emod 100e9 kratio 1.5  cmat default type ball-ball model linear method deform emod 100e6 kratio 1.5   

ball attribute density 2e3 damp 0.7



model cycle 2000 calm 50ball delete range pos-x -5 5 notball delete range pos-y -5 5 notmodel solve model save "ball_dis"

rblock construct from-balls polydisperse true  

model save "rblock_dis"

model restore "rblock_dis"

ball deletedef GetGeo    geoname_count=1    loop foreach rb rblock.list        groupName=string.build("geo_%1",geoname_count)        idRblock=rblock.id(rb)        command            rblock export to-geometry @groupName range id @idRblock        endcommand        geoname_count =1    endloopend@GetGeorblock delete

def GenerateBallIn(geo_count_single)    groupName=string.build("geo_%1",geo_count_single)    command        ball distribute porosity 0.1 radius 0.03 0.06 group @groupName range geometry-space @groupName count 1         ball attribute density 2.7e3 damp 0.3    endcommandend

def genAllBall    loop gen_cur(1,geoname_count-1)        groupName=string.build("geo_%1",gen_cur)        command              wall import  from-geometry @groupName id 10000            model calm            ball fix velocity spin                   endcommand        GenerateBallIn(gen_cur)        command                       model cycle 2000 calm 20            ball delete range geometry-space @groupName count 1 not group @groupName        endcommand        command             model solve ratio-average 1e-5 or cycles 10000            model save @groupName            wall delete walls range id 10000        endcommand     endloopend@genAllBallmodel save "sample"