Matlab根据flac、pfc或其他软件导出的坐标及应力、位移数据再现云图 案例包括导出在flac6.0中导出位移的fish代码（也可以自己先准备软件导出的坐标数据及对应点的位移或应力数据，可根据需要自行修改为自己需要的云图数据）、matlab中读取三维坐标及位移数据然后云图呈现的代码

一、功能概述

本套工具由FLAC脚本（main.f3dat）与Matlab脚本（main.m）组成，核心功能是实现从数值模拟结果到可视化位移云图的完整流程。通过FLAC软件完成三维模型的力学计算后，自动提取模型各关键点的坐标信息与Z方向位移数据，生成标准化数据文件；再借助Matlab脚本读取该数据文件，通过两种可视化方案（散点云图、多边形插值云图）实现位移分布的三维可视化呈现，清晰直观地展示模型在重力作用下的位移响应特征。

二、适用环境与前置要求

（一）软件环境

1. 数值模拟端：支持运行FLAC3D的相关版本（需兼容fish语言脚本执行）；
2. 可视化端：Matlab 21a及以上版本（低版本需注意文件编码兼容性，具体解决方案见下文注意事项）；
3. 运行路径：FLAC脚本运行需在纯英文路径下执行，Matlab脚本运行不受路径中英文影响。

（二）前置条件

1. 已安装上述所需软件，且熟悉FLAC基本建模流程与Matlab基础操作；
2. 确保FLAC模型计算完成后，模型关键点（grid point）的位置与位移数据已正常生成；
3. 数据文件（displacement.dat）在生成后需与Matlab脚本置于同一目录下，或在脚本中明确指定文件路径。

三、核心流程说明

（一）数据采集阶段（FLAC脚本执行）

1. 模型初始化与计算：FLAC脚本首先完成三维圆柱体模型的创建，通过反射操作扩展模型范围，定义模型边界条件与重力场（重力加速度沿Z轴负方向，大小9.8m/s²）；
2. 材料参数配置：为模型赋予摩尔-库仑本构模型及对应的物理力学参数（密度、体积模量、剪切模量、内聚力、摩擦角等），确保力学计算的合理性；
3. 计算执行：设置计算步数并运行，完成模型在重力作用下的力学平衡计算；
4. 数据提取与存储：通过fish语言自定义函数，遍历模型所有关键点，采集各点的X、Y、Z坐标及Z方向位移数据，按固定格式写入displacement.dat文件，完成数据导出。

（二）可视化呈现阶段（Matlab脚本执行）

1. 数据读取与预处理：Matlab脚本读取displacement.dat文件中的结构化数据，分别提取坐标列（X、Y、Z）与位移列（disp），并计算位移数据的最大值与最小值，为后续颜色映射提供范围依据；
2. 可视化方案一：三维散点云图

• 采用scatter3函数绘制三维散点图，每个散点代表模型的一个关键点；
• 以位移数据作为散点颜色映射依据，通过colormap（jet色系）直观区分不同位移大小；
• 配置3D视图、坐标轴比例、颜色条、坐标轴标签与标题，优化图形显示效果，字体与字号统一规范。

1. 可视化方案二：多边形插值云图

• 通过alphaShape函数构建三维点云的多边形包围结构，基于alphaTriangulation进行三角化处理；
• 利用patch函数绘制多边形面，将位移数据作为面颜色插值依据，实现位移分布的平滑过渡显示；
• 保持与散点云图一致的视图配置、颜色映射与标注风格，提供更贴近实体模型的位移分布可视化效果。

四、输出结果说明

1. 中间数据文件：displacement.dat，文本格式，存储模型关键点的X、Y、Z坐标及对应Z方向位移数据，每行代表一个关键点的完整信息，数据间以空格分隔，可直接用文本编辑器打开查看；
2. Matlab可视化图形：两个独立的三维图形窗口，分别为散点云图与多边形插值云图；

• 图形包含坐标轴（X、Y、Z）、颜色条（标注位移数值范围）、标题（Disp picture）；
• 颜色从蓝色（最小位移）到红色（最大位移）渐变，清晰反映位移分布的空间差异。

五、注意事项与常见问题解决

（一）文件编码兼容性

低版本Matlab打开高版本（如21a）编写的.m文件时可能出现乱码，解决方案：

1. 用记事本或其他文本编辑器打开乱码的.m文件；
2. 复制文本编辑器中显示正常的代码；
3. 粘贴到Matlab的.m文件中替换原有乱码内容，保存后即可正常运行。

（二）路径要求

FLAC脚本（main.f3dat）必须在纯英文路径下运行，否则可能导致数据文件导出失败；Matlab脚本无此限制，但需确保能正确读取displacement.dat文件（建议置于同一目录）。

（三）参数调整建议

1. 若散点图显示过于密集或稀疏，可调整scatter3函数中的“30”（散点大小）参数；
2. 若多边形云图的平滑度不符合需求，可修改alphaShape函数中的“0.6”（形状参数），数值越大面越平滑，数值越小越贴近原始点云；
3. 颜色映射方案可通过修改colormap函数（如将jet改为parula、viridis等）调整，适应不同可视化需求。

（四）数据有效性检查

若运行Matlab脚本后无图形显示或报错，需检查：

1. displacement.dat文件是否成功生成，且数据格式正确（无缺失行、无多余字符）；
2. 数据文件中的坐标与位移数据是否合理（无异常超大或超小值）；
3. Matlab脚本中数据读取的列索引是否与.dat文件的列顺序一致。

六、使用建议

1. 首次使用时，建议先运行FLAC脚本确认displacement.dat文件正常生成，再执行Matlab脚本；
2. 若需可视化其他方向（X或Y方向）的位移，可修改FLAC脚本中的数据提取逻辑，对应调整Matlab脚本中的位移列索引；
3. 可视化图形可通过Matlab的图形导出功能保存为高分辨率格式（如.png、.eps），适用于学术报告或技术文档展示；
4. 对于大规模模型数据，可适当优化FLAC的关键点采样密度或Matlab的图形渲染参数，平衡可视化效果与运行效率。