COMSOL with MATLAB代码随机分布球/圆模型及代码。 包含二维三维，打包。 用于模拟多孔介质 二维COMSOL with MATLAB 接口代码 多孔介质生成 以及 互不相交小球生成程序 说明：本模型可以生成固定数目的互不相交的随机小球；也可以生成随机孔隙模型 一、若要生成固定数目的小球，则在修改小球个数count的同时，将n改为1 二、若要生成随机孔隙模型，则将count尽量调大，保证能生成足够多的小球 三维COMSOL with MATLAB代码：随机分布小球模型 功能： 1、本模型可以生成固定小球数量以及固定孔隙率的随机分布独立小球模型 2、小球半径服从正态分布，需要给定半径均值和标准差。 2、若要生成固定小球数量模型，则更改countsph，并将孔隙率n改为1 3、若要生成固定孔隙率模型，则更改孔隙率n，并将countsph改为一个极大值1e6.
搞多孔介质仿真的兄弟应该都懂，随机分布的小球模型简直是刚需。最近折腾了一套二维/三维通杀的MATLAB+COMSOL联动代码，能快速生成互不重叠的随机小球/圆模型，顺手分享几个实战技巧。
二维篇：画布上的小球战争
先看二维场景，核心逻辑就两条——要么固定数量的小球硬塞进去，要么按孔隙率让小球打到"挤不下为止"。上硬菜：
% 二维核心参数 R = 0.1; % 基准半径 count = 30; % 想塞多少个小球？ n = 1; % 模式开关：1=固定数量，0=按孔隙率 % 生成逻辑 for i = 1:count collision = 1; while collision x = rand()*(1-2*R) + R; % 边界留安全距离 y = rand()*(1-2*R) + R; % 碰撞检测 if all(sqrt((x - existing_x).^2 + (y - existing_y).^2) > 2*R) existing_x(i) = x; existing_y(i) = y; collision = 0; end end end这段代码的骚操作在于边界安全距离和暴力循环检测。特别提醒几个坑：
n=1时玩的是俄罗斯方块模式——必须塞满指定数量，这时候count别设太大否则死循环- 改成孔隙率模式时，把count调到500+，让小球疯狂生成直到塞不下
实战中建议加个尝试次数计数器，超过500次还没位置就自动break，避免程序卡死。
三维篇：球球大作战升级版
三维版本多了正态分布半径的玩法，代码更有内味儿：
% 三维参数暴击 mu = 0.08; % 平均半径 sigma = 0.02; % 标准差 porosity = 0.4; % 孔隙率 countsph = 1e6; % 数量模式开关 % 半径生成器 radii = abs(normrnd(mu, sigma, [countsph,1])); % 空间填充算法 maxAttempts = 1000; % 防卡死神器 for i = 1:countsph attempt = 0; while attempt < maxAttempts % 三维坐标生成 coord = rand(1,3).*(1-2*radii(i)) + radii(i); % 升级版碰撞检测 distances = sqrt(sum((coord - existing_coords).^2, 2)); if all(distances > (radii(i) + existing_radii)) % 记录新球参数 break; end attempt = attempt + 1; end end这里藏着几个骚操作：
abs(normrnd)防止出现负半径的智障情况- 三维碰撞检测用矩阵运算代替循环，速度直接起飞
maxAttempts是保命符——遇到实在塞不下的大球直接跳过
重点提醒：做孔隙率模式时，计算域体积要预先扣除孔隙部分，边塞球边计算当前总体积，超过阈值就停止。
避坑指南
- 二维边界处理：生成坐标时记得
(1-2*R)这个操作，否则靠近边界的球会被COMSOL判定为超出计算域 - 正态分布陷阱：当sigma设太大时，用
abs()强制非负可能导致实际均值偏移，建议加个均值补偿 - 性能玄学：三维代码里把existing_coords存成预设数组，比动态扩容快10倍不止
- 可视化彩蛋：在MATLAB里用scatter3画生成结果时，加上
alpha(0.5)半透明效果，重叠检测一目了然
这套代码实测在i7笔记本上，生成500个三维小球大概需要2分钟（没开并行）。想要更快的可以试试空间八叉树分割法，不过代码量得翻三倍...
最后扔个参数组合套餐给选择困难症：
- 新手套餐：二维+固定数量+均一半径
- 进阶套餐：三维+孔隙率模式+正态半径
- 魔鬼套餐：三维+动态调整sigma+实时孔隙率计算
模型生成后直接mphimport导回COMSOL，配上Darcy流模块，一套多孔介质仿真流水线直接拉满。下次遇到审稿人质疑模型随机性，反手就是20组不同随机种子的对比结果糊脸——就问你稳不稳？
