計算結果如圖：

第一步：建立單軸壓縮模型

new

title 'Testing Bonded Particle Model'

def canshu 

    width = 0.05

    height = 0.1

    rmax= 1.66e-3

    rmin=  1e-3

    poro=  0.3

end

@canshu

; Set the domain extent

domain extent -0.03 0.03 -0.03 0.03 -0.06 0.06 condition destroy

cmat default model linear method deformability emod 1.0e9 kratio 1.0

cmat default property dp_nratio 0.5

; create walls extending past the edges of the sample

wall generate id 1 plane dip 0 dd 0   position 0 0 0.05

wall generate id 2 plane dip 0 dd 0   position 0 0 -0.05

wall generate id 3 plane dip 90 dd 90 position -0.025 0 0

wall generate id 4 plane dip 90 dd 90 position 0.025 0 0

wall generate id 5 plane dip 90 dd 0  position 0 -0.025 0

wall generate id 6 plane dip 90 dd 0  position 0 0.025 0

set random 10002

ball distribute porosity @poro ...

                resolution 1.0 ...

                numbins 1 ...

                bin 1 ...

                    radius [rmax*0.5] [rmin*0.5] ...

                    volumefraction 1 ...

                    group shazi ...  

                box -0.025 0.025 -0.025 0.025 -0.05 0.05

ball attribute density 2500 damp 0.7

; Calm the system

cycle 2000 calm 10

; Solve the system to a target limit (here the average force ratio)

; Use density scaling to quickly reach equilibrium

set timestep scale

solve aratio 1e-4

set timestep auto

calm

; Delete side walls

; Be careful to include the keyword 'walls' or else facets will be deleted

wall delete walls range id 3 6

; make cylinder

ball delete range cylinder end1 0 0 -0.05 end2 0 0 0.05 rad 0.025 not

wall generate id 7 cylinder axis 0 0 1  ...

                            base 0 0 -0.06 ...

                            height 0.12 ...

                            radius 0.025 ...

                            cap false false ...

                            onewall ...

                            resolution 0.1

cycle 2000 calm 10

set timestep scale

solve aratio 5e-4

set timestep auto

ball delete range cylinder end1 0 0 -0.05 end2 0 0 0.05 rad 0.025 not

save unbonded

第二步：導入xdf文件并生成clump，利用ball對其進行替換，生成cluster