手把手玩转COMSOL+MATLAB生成多孔介质模型

cvcNYgAdnd

708人浏览 · 2026-03-16 08:30:00

cvcNYgAdnd · 2026-03-16 08:30:00 发布

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

这里藏着几个骚操作：