离散元后处理工具集：使用PFC数据绘制云图并导入MATLAB生成三维图形

离散元后处理，将PFC数据写出并导入到matlab中形成云图。 支持二维三维绘图。 内容包括：案例文件，fish代码和matlab代码。 pfc2d、3d_to_matlab.

最近在折腾离散元后处理，发现PFC自带的绘图功能实在不够用。刚好手头有个边坡破坏的案例，咱们试试用Matlab搞点炫酷的云图。先上效果图镇楼（此处假装有图）——从颗粒级应力场到速度场的动态可视化，Matlab玩起来是真的香。

先说二维情况。在PFC里跑模型的时候，记得在Fish脚本里加个数据记录器：

; 颗粒应力记录器 def save_stress fopen 'stress_data.txt' write ; 每次运行覆盖写入 loop foreach local gp ball.list ; 遍历所有颗粒 local pos = ball.pos(gp) local stress = ball.extra(gp,1) ; 假设应力存在extra 1 fprintf string.format("%.4f %.4f %.4f\n",pos.x,pos.y,stress) end_loop fclose end @save_stress ; 记得在需要时调用

这坨代码会把每个颗粒的坐标和应力值存成文本文件。注意文件路径别翻车，建议用绝对路径更稳妥。数据量大的时候记得分块保存，别一股脑全塞内存里。

转到Matlab这边，先整个暴力读取：

data = load('stress_data.txt'); x = data(:,1); y = data(:,2); stress = data(:,3); % 生成网格 xi = linspace(min(x),max(x),100); yi = linspace(min(y),max(y),100); [XI,YI] = meshgrid(xi,yi); % 插值上色 ZI = griddata(x,y,stress,XI,YI,'natural'); % 画图 figure('Color','white') pcolor(XI,YI,ZI) shading interp colorbar hold on contour(XI,YI,ZI,20,'k-')

这里有个坑——griddata插值遇到不规则区域会抽风。遇到这种情况建议上scatter直接画散点云图，加点透明度更带感：

scatter(x,y,30,stress,'filled','MarkerEdgeAlpha',0.3)

三维版本更刺激。PFC3D的Fish代码改改就能用：

def save_3d_stress fopen 'stress3d.txt' write loop foreach gp ball.list pos = ball.pos(gp) stress = ball.prop(gp,'stress_zz') fprintf string.format("%.4f %.4f %.4f %.4f\n",pos.x,pos.y,pos.z,stress) end_loop fclose end

Matlab这边处理三维数据得用点骚操作：

data3d = importdata('stress3d.txt'); x = data3d(:,1); y = data3d(:,2); z = data3d(:,3); s = data3d(:,4); % 三维插值 [XI,YI,ZI] = meshgrid(linspace(min(x),max(x),50),... linspace(min(y),max(y),50),... linspace(min(z),max(z),50)); VI = griddata(x,y,z,s,XI,YI,ZI); % 切片可视化 slice(XI,YI,ZI,VI,[],[],linspace(min(z),max(z),5)) shading interp colormap jet

遇到百万级颗粒数据时，建议上bin文件读写。Fish里用file.write写二进制，Matlab用fread读取，速度能快十倍。不过小心字节顺序问题，跨平台读写建议用同款操作系统。

最后甩个自用的PFC-Matlab直连工具（假装此处有GitHub链接）。这玩意儿能自动同步模型边界和颗粒数据，支持实时可视化更新。下次准备搞个GPU加速版本，有兴趣的可以留言区扣1（笑）。