【可下载】用Deepseek做了一个随机生成锂离子电池二维多孔电极COMSOL模型

// ========== 参数定义部分 ==========
// 模型参数
double ELECTRODE_WIDTH = 50; // 电极宽度 (μm)
double ELECTRODE_HEIGHT = 85; // 电极高度 (μm)
double RIND = 0.5; // 边界保护层厚度 (μm)
double POROSITY = 0.3; // 目标孔隙率
// 颗粒参数
double MIN_RADIUS = 1; // 最小颗粒半径 (μm)
double MAX_RADIUS = 12.5; // 最大颗粒半径 (μm)
double number_particle = 150; // 最大颗粒数
// ========== 几何初始化部分 ==========
model.component("comp1").geom("geom1").lengthUnit("um");
// 设置几何单位为微米
// ========== 基底创建部分 ==========
// 创建电极基底(带边界保护)
model.component("comp1").geom("geom1").create("base", "Rectangle");
with(model.component("comp1").geom("geom1").feature("base")); //定义矩形尺寸
set("size", new double[]{
ELECTRODE_WIDTH,
ELECTRODE_HEIGHT});
endwith();
// ========== 颗粒生成部分 ==========
int successCount = 0;
double totalArea = (ELECTRODE_WIDTH)*(ELECTRODE_HEIGHT);
double targetSolidArea = totalArea*(1-POROSITY);
double currentSolidArea = 0;
while (successCount < number_particle && currentSolidArea < targetSolidArea) {
// 生成正态分布半径
double rand1 = Math.random();
double r = MIN_RADIUS+(MAX_RADIUS-MIN_RADIUS)*rand1;
r = Math.max(MIN_RADIUS, Math.min(MAX_RADIUS, r));
double area = Math.PI*r*r;

// 跳过超出剩余空间的颗粒
if (currentSolidArea+area > targetSolidArea) continue;

// 生成随机位置(考虑边界保护)
double rand2 = Math.random();
double rand3 = Math.random();
double x = RIND+r+(ELECTRODE_WIDTH-2*(RIND+r))*rand2;
double y = RIND+r+(ELECTRODE_HEIGHT-2*(RIND+r))*rand3;

// 碰撞检测
boolean collision = false;
for (int i = 0; i < successCount; i++) {
double[] existPos = model.component("comp1").geom("geom1").feature("part"+i).getDoubleArray("pos");
double existR = model.component("comp1").geom("geom1").feature("part"+i).getDouble("r");

double dx = x-existPos[0];
double dy = y-existPos[1];
double dist = Math.sqrt(dx*dx+dy*dy);

if (dist < (r+existR)) {
collision = true;
break;
}
}

// 创建颗粒
if (!collision) {
String ftName = "part"+successCount;
model.component("comp1").geom("geom1").create(ftName, "Circle");
model.component("comp1").geom("geom1").feature(ftName)
.set("r", r)
.set("pos", new double[]{x, y});

successCount++;
currentSolidArea += area;
}
}
// ========== 后处理部分 ==========
// 计算实际孔隙率
double actualPorosity = 1-(currentSolidArea/totalArea);
model.param().set("actual_porosity", actualPorosity);
// 生成几何
model.component("comp1").geom("geom1").run();