model new

model domain extent -10 10geometry set "box"geometry generate box -5 5

def par    origin_rad_min=0.2    origin_rad_max=origin_rad_min*1.5    scale=3    block_radius_max=origin_rad_max*scale    block_radius_min=origin_rad_min*scaleend@parrblock construct from-geometry "box"  minimum-edge  @origin_rad_min ...        maximum-edge @origin_rad_max group "origin" slot "1"

def inbox(pos,pointer)    groupName=rblock.group(pointer)    if groupName#"1=origin" then        inbox=false        exit    endif    vl=vector(x_pos_random-sousuo_rad_x,y_pos_random-sousuo_rad_y)    vh=vector(x_pos_random+sousuo_rad_x,y_pos_random+sousuo_rad_y)    rbpos=rblock.pos(pointer)    if comp.x(rbpos)>comp.x(vl) & comp.x(rbpos)<comp.x(vh) &comp.y(rbpos)>comp.y(vl) & comp.y(rbpos)<comp.y(vh) then                inbox=true        exit    endif    inbox=falseend
def conbine_rblock    stop_flag=true    name_count=1    loop while stop_flag        x_pos_random=math.random.uniform*10-5        y_pos_random=math.random.uniform*10-5        sousuo_rad_x=math.random.uniform*(block_radius_max-block_radius_min)+block_radius_min        sousuo_rad_y=math.random.uniform*(block_radius_max-block_radius_min)+block_radius_min        ishaveRb=false        loop foreach rb rblock.list            vl=vector(x_pos_random-sousuo_rad_x,y_pos_random-sousuo_rad_y)            vh=vector(x_pos_random+sousuo_rad_x,y_pos_random+sousuo_rad_y)            if inbox(vl,rb) then                ishaveRb=true                exit loop            endif        endloop                group_name=string.build("zu_%1",name_count)        name_count+=1        if ishaveRb=true then                    command                                rblock merge group "new" slot "1" group @group_name slot "10"  ...                        range fish @inbox                model clean            endcommand        endif        orgin_count=0        loop foreach rb rblock.groupmap("origin",1)            orgin_count+=1        endloop        if orgin_count==0 then            stop_flag=false        endif    endloopend@conbine_rblock

def findMinRb    VLimit=origin_rad_min*origin_rad_min*10    loop foreach rb rblock.groupmap("new",1)        if rblock.vol(rb)<VLimit*4 then            rblock.group(rb,1)="minRB"        endif    endloopend@findMinRb

def findNearestRb(pos)    minPos=1e100    minId=0    loop foreach rbfind rblock.groupmap("new",1)        vecLenghth=rblock.pos(rbfind)-pos        if math.mag(vecLenghth)<minPos then            minPos=math.mag(vecLenghth)            minId=rblock.id(rbfind)        endif    endloop    findNearestRb=minIdend

def solveMinRb        loop foreach rb  rblock.groupmap("minRB",1)        rbNearest=rblock.find(findNearestRb(rblock.pos(rb)))        rblock.group(rb,20)="dai"        rblock.group(rbNearest,20)="dai"        group_name=string.build("zu_%1",name_count)        name_count+=1        command             rblock merge group "new" slot "1" group @group_name ...                slot "10" range group "20=dai"        endcommand        loop foreach rbdai rblock.groupmap("dai",20)            rblock.group(rbdai,20)="none"        endloop    endloopend@solveMinRb

model restore "rblock_slip"rblock export to-geometry "allgeo" cmat default type ball-facet model linear method deform emod 100e6 kratio 1.5 property fric 0.5cmat default type ball-ball model linear method deform emod 100e6 kratio 1.5  property fric 0.5def 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 ballDeleteGeo(geo_count_single)    if geo_count_single=1 then        exit    endif    thisGroup=string.build("geo_%1",geo_count_single)    loop n(1,geo_count_single-1)         groupName=string.build("geo_%1",n)        command            ball delete range geometry-space @groupName count 1 group @thisGroup        endcommand    endloopenddef 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    endcommand    ballDeleteGeo(geo_count_single)end

wall import  from-geometry "box" id 100 def genWalls(geo_count_single)    loop n(1,geo_count_single)         groupName=string.build("geo_%1",n)        command            wall import  from-geometry @groupName id [n+10000]         endcommand    endloopenddef deleWalls(geo_count_single)    loop n(1,geo_count_single)        command            wall delete walls range id [n+10000]         endcommand    endloopend

def genAllBall    loop gen_cur(1,geoname_count-1)                genWalls(gen_cur)        command                       model calm            ball fix velocity spin                   endcommand        GenerateBallIn(gen_cur)        groupName=string.build("geo_%1",gen_cur)        command                       model cycle 2000 calm 20            ball delete range geometry-space @groupName count 1 not group @groupName        endcommand        ballDeleteGeo(gen_cur)        command             model solve ratio-average 1e-5 or cycles 10000            model save @groupName        endcommand        deleWalls(gen_cur)    endloopend@genAllBallball free velocity spinmodel cycle 1model solvemodel save "sample"