PFC3D 5.0代码 可破碎的cluster,可模拟碎石、矿渣混凝土材料
计算结果如图:
第一步:建立单轴压缩模型
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
该付费内容为:PFC3D 5.0代码 可破碎的cluster,可模拟碎石、矿渣混凝土材料,cluster替换clump 单轴压缩试验模拟
1人购买