New
清除所有颗粒
t random
/随机产生颗粒
t disk on ; treat balls as disks of unit thickness
title ’Idealized mine block-caving process’
为整个程序命名,在窗体上显示为“Job title 獾子油。。。”
wall id 1 ks 1e9 kn 1e9 fric 0.5 nodes -20,0 20,0 / 生成墙
wall id 墙编号 ;ks 切向刚度;kn 法向刚度;fric 摩擦系数;nodes x,y 产生墙的起始及结束坐标
gen x -20,20 y 0,3 id=20000,22550 rad 0.05,0.1 / 生成颗粒
x a1 ,b2 y b1,b2 ,产生颗粒范围;id 颗粒编号;rad r1, r2 产生颗粒半径范围;
ini rad mul 1.55 range id = 20000,22550
ini 初始条件设置;rad mul n 颗粒半径扩大倍数;id 颗粒编号;
plot create cave1
产生新的窗体
plot add ball lblue range (id 青春骄纵)
添加显示颗粒,range 确定范围,x y;id确定编号;if语句lblue 颗粒颜色;
(颗粒颜色:black,blue,green,cyan蓝绿色,red,magenta紫红色,orange字笑,brown,
lgray灰色,dgray,lblue,lgreen,lcyan,lred,lmagenta,lorange,yellow,white)
plot add clump
显示簇
plot add cf black
添加接触应力
wall id 6 xvel 5
设置墙体速度,x正方向数值为正,x负方向数值为负。
wall id 6 spin 5 x a1 y b1
墙体以(a1,b1)为中心,以5rad/s的角速度旋转,正值为顺时针旋转,负值为逆时针旋转
hist id 120 wall power id 6
3322监测墙体与球接触力的大小
hist id=1 diag mcf
监测平均接触力
hist id=2 diag muf
监测平均不平衡力
meas id 1 x -14 y 6 rad 6
设置测圆,measure;id,编号;x y ,测圆中心坐标;rad,测圆半径
hist id=3 measure porosity id=1
监测测圆孔隙度,id=3 梅州灵光寺监测编号;id=1测圆编号
prop ks 1e9 kn 1e9 dens 2000泛泛而谈
设置颗粒材料参数,ks kn,切向法向刚度;dens 材料密度单位kg/m3
t gravity 0,-10
设置重力加速度
cyc 6000
循环次数,或者step命令
save cave0.sav
save cave0.dat
保存文件
调用记事本编辑的程序
Restore abc.sav
Restore abc.dat
调用保存的数据文件
delete balls range id=20000,22550
删除颗粒 range 确定范围
hist id=101 ball yposition 0 4
监测单个颗粒y方向位移,此颗粒质心坐标距给定坐标(0,4)最近
t log on
print hist
t log off
记录文件及其数据存取的过程,并在与原文件相同的文件夹里产生pfc2d.log记事本文件
range name meiceng x 0 500 y 20 30
自定义名称代替范围
macro raft 'id=5000,5079'
property pb_kn 1e10 pb_ks 1e10 c_index 1 range raft
自定义任意字符,代替常用或重复使用的字符串
Solve
自动运行至合适的结果
delete wall 3
删除墙
property pb_rad 1.0 pb_s 1e20 pb_n 1e20 range raft
设置平行粘结半径,切向法向刚度
ini xv 0 range Raft
ini yforce -4e6.98 range Raft
x方向速度,y方向应力
t echo off
t x0=-19.75 y0=12.25 radius=0.25 /确定初始颗粒质心坐标及颗粒半径
t id_start=5000 n_col=80 n_row=1 /所产生颗粒的起始id膝盖英语怎么读编号,每排颗粒数及排数
chansheng
t echo on
产生规则排列的颗粒。
def chansheng
xc = x0
yc = y0
rc = radius
idc = id_start
r2 = 2.0 * radius
yinc = radius * sqrt(3.0)
loop row (1,n_row)
loop col (1,n_col)
command
ball id=idc x=xc y=yc rad=rc
end_command
idc = idc + 1
xc = xc + r2
end_loop
yc = yc + yinc
xc = x0 + radius * (row - (row/2) * 2)
end_loop
end
产生规则颗粒的fish程序
ball x 0 y 2 rad 1
产生单个颗粒,质心坐标,半径
ball x 0 y 2 rad 1
ball x 2 y 2 rad 1
ball x 4 y 2 rad 1
clump id=1 range x= -1 5 y=1,3
产生一系列单个颗粒,将这些颗粒连接成簇,x y 分别限定簇起作用的范围