矩形三维随机裂隙网络。 使用COMSOL with Matlab接口编程。 可以直接导入COMSOL中，无需CAD，无需提取数据，方便快捷可以直接计算。 裂隙由matlab编程生成，能够生成两组不同产状的裂隙。 裂隙长度的分布律可以为确定的裂隙长度，也可以为在一定范围内随机均匀分布的长度。 注释十分详细，有运行的示范视频，可以直接改数据生成需要的三维裂隙网格。 三维随机裂隙网络模型均为自己编程，保证能够运行。 可以生成多组不同产状的裂隙。 （任意组数都行）

COMSOL-Matlab 协同下的三维随机裂隙网络全自动建模方案

矩形三维随机裂隙网络。 使用COMSOL with Matlab接口编程。 可以直接导入COMSOL中，无需CAD，无需提取数据，方便快捷可以直接计算。 裂隙由matlab编程生成，能够生成两组不同产状的裂隙。 裂隙长度的分布律可以为确定的裂隙长度，也可以为在一定范围内随机均匀分布的长度。 注释十分详细，有运行的示范视频，可以直接改数据生成需要的三维裂隙网格。 三维随机裂隙网络模型均为自己编程，保证能够运行。 可以生成多组不同产状的裂隙。 （任意组数都行）

——从脚本驱动到几何 partition 的零 CAD 工作流

一、背景与需求

在岩土、水文及石油工程中，裂隙网络的统计特征直接决定渗流、力学与传热仿真的可信度。传统做法先在外部 CAD 平台生成裂隙片，再逐片导入 COMSOL 做布尔分割，过程繁琐、易丢参、难复现。本文介绍的脚本化方案完全基于 COMSOL with Matlab 接口，无需任何外部文件，即可在内存中一次性完成“统计参数→几何实体→可计算网格”的闭环，实现：

1. 两组不同产状裂隙的批量随机投放；
2. 裂隙尺度支持“定长”“均匀分布”“正态分布”三种律；
3. 裂隙形态可在“圆盘”与“矩形”之间一键切换；
4. 所有裂隙自动归入同一 CumulativeSelection，随后一次 Partition 即生成可网格化的复合体；
5. 全程参数化，可嵌入优化或反演循环。

二、整体架构

脚本采用“模型句柄 + 几何工厂 + 选择集 + 一次性布尔”四级架构：

• 模型句柄：ModelUtil.create('Model') 在 Matlab 工作区开出 COMSOL 进程，返回 model 对象，后续所有操作均为内存级 API 调用，无磁盘 IO。
• 几何工厂：以 WorkPlane 为裂隙载体，通过 transdispl 与 transeulerang 实现三维空间刚体变换；圆盘/矩形在各自 WorkPlane 内原生创建，保证法向与尺寸精度。
• 选择集：所有裂隙统一 contributeto 到名为 csel1 的 CumulativeSelection，后续 Partition 仅需引用该命名选择集，无需关心裂隙数量。
• 一次性布尔：Partition 特征将 blk1 与 csel1 做布尔分割，生成共享拓扑的复合实体，直接具备网格划分条件。

三、核心功能拆解

3.1 参数化岩体容器

仅需修改 c,k,g 三变量即可定义长方体范围；脚本自动把容器命名为 blk1，作为后续 Partition 的“被切”对象。

3.2 分组裂隙统计

• count 控制组数；num 数组按序存放每组裂隙数量，支持任意正整数。
• 走向/倾向均值与方差独立可配，例如 zx_mu=[90 0] 表示第一组平均走向 90°、第二组走向 0°，方差为 0 时完全平行。
• 圆盘半径或矩形迹长通过 flag 切换分布律，其余参数自动失效，减少误配。

3.3 空间随机投放

中心点坐标在 [0,c]×[0,k]×[0,g] 内均匀抽样；欧拉角采用 normrnd 生成，符合 Fisher 或正态假设。所有随机种子继承 Matlab 当前状态，便于重复试验。

3.4 形态切换

• 圆盘脚本：在 WorkPlane 内创建 Circle，半径变量 r 按分布律生成。
• 矩形脚本：创建 Rectangle，长/宽可独立抽样，支持非等比矩形。

通过替换脚本文件即可在两种形态间切换，无需改动主流程。

3.5 自动选择与布尔

裂隙生成后立即 contributeto 到 csel1，Partition 阶段仅两行代码：

feature('par1').selection('input').set({'blk1'});
feature('par1').selection('tool').named('csel1');

无论裂隙数量是 10 还是 10 000，脚本行数不变；COMSOL 内部采用并行布尔算法，显著优于逐片循环布尔。

3.6 可视化与持久化

runPre('fin') 触发几何最终化；mphgeom 在 Matlab 图形窗弹出三维视图；mphsave 将模型保存为 .mph 文件，用户可继续在 GUI 中添加物理场、材料及研究步骤。

四、运行环境

• COMSOL Multiphysics ≥ 5.5，带 Matlab LiveLink；
• Matlab ≥ R2019b（支持 COMSOL 6.x 的 Java API）；
• 内存建议：每条裂隙约占用 1–2 KB 几何数据，百万级裂隙需 8 GB 以上内存。

五、使用示例

1. 在 Matlab 命令行切换至脚本目录；
2. 调用 model = randomfracture3Dround; 或 randomfracture3Drectangle;
3. 脚本结束后自动弹出几何视图，并在当前目录生成 randomfracture3D.mph；
4. 打开 COMSOL GUI，导入该文件，即可基于“裂隙网络”继续定义渗流或固体力学仿真。

六、扩展方向

• 裂隙开度：在 WorkPlane 内再建一圆/矩形，与原始面做 Difference，即可生成带真实厚度的裂隙通道；
• 非贯通裂隙：引入截断正态或 Beta 分布，控制半径上限，避免裂隙穿透整个岩块；
• 非均匀各向异性场：将 zxmu/qxmu 改为三维随机场数组，实现 P32 密度空间变异；
• 优化耦合：把脚本封装为函数句柄，目标输出设为等效渗透系数，即可与遗传算法或贝叶斯反演结合，实现“参数-响应”自动校准。

七、小结

该脚本化建模方案将“统计参数→几何→网格”的完整链路压缩到数十行 Matlab 代码，彻底摆脱 CAD 中间格式，显著降低人为干预。其“选择集 + 一次布尔”的设计模式，可线性扩展到数万条裂隙，为后续多物理场仿真提供高保真、可复现、可版本控制的几何输入。