内容发布更新时间 : 2024/11/9 23:44:24星期一 下面是文章的全部内容请认真阅读。

在英文输入的状态下，按m是放大视图，按shift 和m是缩小视图。 按x/y/z分别是视图绕x/y/z轴旋转。 set hist rep=5

hist ball zvel（z轴方向的速度） 3 1 0（球心的坐标） hist diagnostic muf

set dt dscale（设置密度）

set dt=value set dt auto

PLOT ball id on（the ID numbers will be largeenough to see） poros = 0.40; desired final porosity

rat = 1.5; ratio of largest to smallest radii mult_a = 0 ; storage for l.h. radius multiplier setup，什么意思，是设置参数的意思吗 前面必须有def吗 def make_block

tot_vol = (x2 - x1) * (y2 - y1) * (z2 - z1) get_poros; Now measure actual porosity

plot set title text ’Basic collection of particles in a box’设置标题 density-scaling选项（简写为dscale） Cycle 3000 设定求解步数

当写入这句程序后，计算就开始了，计算的步数，时间步，总时间，不平衡力平均值和最大不平衡力这些信息均显示在屏幕上，并且每隔5个计算步更新一次（这是set hist_rep=5的作用），图像动画每隔20个计算步更新一次（默认），这个更新时间可用set pinterval 来改变。另外，阻尼已经在运动方程中被默认设定，用于吸收振动能量，更快地达到平衡状态。 10查看结果

求解完毕后，用户可输入命令plot his 1来查看历史记录中第一个变量的数值曲线。 wall id 6 spin 5 x a1 y b1墙体以（a1,b1）为中心，以5rad/s的角速度旋转，正值为顺时针旋转，负值为逆时针旋转

hist id 120 wall power id 6 监测墙体与球接触力的大小 hist id=2 diag muf 监测平均不平衡力

hist id=101 ball yposition 0 4监测单个颗粒y方向位移，此颗粒质心坐标距给定坐标（0，4）最近

set log on print hist set 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 自定义任意字符，代替常用或重复使用的字符串

property pb_rad 1.0 pb_s 1e20 pb_n 1e20 range raft 设置平行粘结半径，切向法向刚度 Solve自动运行至合适的结果

ini xv 0 range Raft ini yforce -4e6.98 range Raft set x0=-19.75 y0=12.25 radius=0.25 /确定初始颗粒质心坐标及颗粒半径

set id_start=5000 n_col=80 n_row=1 /所产生颗粒的起始id编号，每排颗粒数及排数 plot hist 3 显示第三个变量的监测结果