python控制修改comsol边界条件仿真方法

lqwh53542025-12-02 13:54

直接明了贴代码：

python 复制代码

from mph import Client
import mph
import numpy as np
import csv
import time


client = Client()

electrodes_1a = [144, 135, 126, 117, 108, 99, 90, 78]
electrodes_gnd = [149, 258, 364, 470, 576, 682, 788, 886]


model = client.load(r"E:\Wuxh\comsol\SR\CEP_basis\basis_model.mph")
OUTPUT_CSV = r"E:\Wuxh\comsol\SR\CEP_basis\R_matrix.csv"


Rmat = np.zeros((8, 8))
total = len(electrodes_1a) * len(electrodes_gnd)
count = 0
start_time = time.time()

for i, term_bd in enumerate(electrodes_1a):
    for j, gnd_bd in enumerate(electrodes_gnd):
        count += 1
        print(f"\n--- Running {count}/{total} : Terminal={term_bd}, Ground={gnd_bd} ---")

        model.java.component("comp1").physics("ec").feature("gnd1").selection().set(gnd_bd)
        model.java.component("comp1").physics("ec").feature("term1").selection().set(term_bd)

        print("gnd1 set:",model.java.component("comp1").physics("ec").feature("gnd1").selection())
        print("term1 set:",model.java.component("comp1").physics("ec").feature("term1").selection())

        model.java.sol("sol1").runAll()

        gev = model.java.result().numerical("gev1")

        tbl = model.java.result().table("tbl1")
        # tbl.clearTable()    # 重置
        gev.set("table", "tbl1")
        gev.setResult()
        print("计算结果：" ,tbl.getReal()[0][0])
        Rmat[i, j] = tbl.getReal()[0][0]

# 保存 CSV
with open(OUTPUT_CSV, "w", newline="") as f:
    writer = csv.writer(f)
    writer.writerow(["term\\gnd"] + electrodes_gnd)
    for idx, row in enumerate(Rmat):
        writer.writerow([electrodes_1a[idx]] + list(row))

elapsed = time.time() - start_time
print("\n==============================")
print("全部 64 组仿真完成！")
print(f"矩阵已保存至: {OUTPUT_CSV}")
print(f"总耗时: {elapsed:.1f} 秒")
print("==============================\n")

这是一个不断修改接地和接电流边界来不断跑仿真获得电阻的代码，需要准备的是已经能够跑通仿真的源文件，然后打开comsol server，运行即可，至于这里面的每一部分怎么设置，得看api，我是拿着ai给的建议一个一个试的，试了一下午终于试完。另外这里其实应该还要划分网格之类的操作，但是因为我原始模型都设置好了，然后改边界不影响网格分布，所以可以直接运行。