超级电容器建模

分类:分析 | 综合版 2026-04-06

超级电容器建模理论基础

🎓

电气双层电容器的离子传输和充放电特性。多孔质电极的PNP方程。能量存储装置。




支配方程




$$ C = \frac{\varepsilon A}{d} $$
$$ i = C\frac{dV}{dt} + \frac{V}{R_{leak}} $$



🧑🎓

等等,超级电容器建模...也适用于这种情况吗?


离散化方法

🧑🎓

这些方程在计算机上怎样实际求解呢?


🎓

使用有限元法(FEM)进行空间离散化。组装单元刚度矩阵,构造整体刚度方程。


🎓

转换为弱形式(变分形式),使用试验函数和形状函数进行Galerkin方法定义。单元类型的选择(低阶单元 vs. 高阶单元完全积分 vs. 降低积分)直接关系到解的精度与计算成本的权衡。




矩阵求解算法

🧑🎓

矩阵求解算法具体是什么意思呢?


🎓

使用直接法(LU分解Cholesky分解)或迭代法(CG法GMRES法)求解连立方程。对于大规模问题,有预处理的迭代法很有效。



求解法分类内存使用量适用规模
LU分解直接法O(n²)小~中规模
Cholesky分解直接法(对称正定值)O(n²)小~中规模
PCG法迭代法O(n)大规模
GMRES法迭代法O(n·m)大规模、非对称
AMG预处理预处理O(n)超大规模
🧑🎓

也就是说,在有限元法那部分偷工减料的话,后面会吃苦头对吧?记住了!


商用工具中的实现

🧑🎓

那超级电容器建模可以用什么软件呢?


工具名开发者/现在主要文件格式
COMSOL MultiphysicsCOMSOL AB.mph
Ansys Mechanical(旧ANSYS Structural)Ansys Inc..cdb, .rst, .db, .ans, .mac
Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
Ansys FluentAnsys Inc..cas, .dat, .msh, .jou

供应商系统和产品整合过程

🧑🎓

每个软件的发展历史都很有戏剧性吗?



COMSOL Multiphysics

🧑🎓

请给我讲讲"COMSOL Multiphysics"!


🎓

1986年在瑞典成立。作为MATLAB链接的FEMLAB开始,后改名为COMSOL。多物理是强项。

现在的所属公司:COMSOL AB



Ansys Mechanical(旧ANSYS Structural)

🧑🎓

请给我讲讲"Ansys Mechanical"!


🎓

1970年由Swanson Analysis Systems Inc. (SASI)开发。以APDL(Ansys参数设计语言)为基础。

现在的所属公司:Ansys Inc.




Simcenter STAR-CCM+

🧑🎓

接下来是Simcenter STAR吧。什么内容呢?


🎓

由CD-adapco开发。2016年被西门子收购,整合到Simcenter品牌。以多面体网格为特点。

现在的所属公司:Siemens Digital Industries Software


🧑🎓

啊,原来是这样啊!"年在瑞典成立"是这种机制啊。


文件格式和互操作性

🧑🎓

在不同软件间交换数据时有什么要注意的吗?


格式扩展名种类概述
STEP.stp/.step中立CAD符合ISO 10303的3D CAD数据交换格式。支持形状+PMI。
IGES.igs/.iges中立CAD初期CAD数据交换规范。曲面数据互操作性有问题。逐步迁移到STEP。
VTK.vtk/.vtu可视化Visualization Toolkit格式。ParaView等使用。
🎓

在不同求解器间转换模型时,要注意单元类型的对应关系、材料模型的兼容性、荷载和边界条件的表达差异。尤其是高阶单元或特殊单元(内聚单元、用户定义单元等)往往无法直接在求解器间转换。


🧑🎓

原来如此...文件格式看起来简单,其实深度很深啊。


实务注意事项

🧑🎓

教科书上没有的"现场智慧"有吗?


🎓

网格收敛性确认、边界条件的合理性验证、材料参数的敏感性分析非常重要。


🎓
  • 网格依赖性验证:至少用3个网格密度水平确认收敛性
  • 边界条件合理性:设置物理上有意义的约束条件
  • 结果验证:与理论解、实验数据、已知基准问题对比



  • 🎓

    嗯,很好的势头!实际动手操作才是最好的学习。有不懂的地方随时问我。


    Coffee Break 闲话

    为什么超级电容器"不是电池"——双层的物理学

    超级电容器(电气双层电容器、EDLC)的理论魅力在于其"储能方式"与电池完全不同。电池通过化学反应储能,但EDLC只是将离子物理吸附在电极表面形成的厚度仅1~2 nm的"电气双层"中。因此充放电极快,理论上化学反应引起的劣化不存在。活性炭电极的比表面积每克可达1,000~3,000 m²,这巨大的表面积是大容量的源头。建模时,双层静电容量以Gouy-Chapman-Stern理论为基础,但多孔质电极内的离子传输与双层的联成是理论的核心课题。

    超级电容器建模数值计算方法



    离散化表述



    🎓

    形状函数 $N_i$ 近似未知量:



    $$ u^h(\mathbf{x}) = \sum_{i=1}^{n} N_i(\mathbf{x}) \, u_i $$




    🎓

    这用数式表示就是这样。


    $$ K_e = \int_{\Omega_e} B^T \, D \, B \, d\Omega \approx \sum_{g=1}^{n_g} w_g \, B^T(\xi_g) \, D \, B(\xi_g) \, |J(\xi_g)| $$

    基本方程离散形式


    🎓

    用数式表示就是这样。


    $$ C = \frac{\varepsilon A}{d} $$
    $$ i = C\frac{dV}{dt} + \frac{V}{R_{leak}} $$

    🧑🎓

    嗯...只看公式不太明白呢...表示什么呢?


    🎓

    连续体的控制方程离散化后,得到如下代数方程组:



    $$ [K]\{u\} = \{F\} $$


    🎓

    这里 $[K]$ 是整体刚度矩阵(或等价的系统矩阵),$\{u\}$ 是未知节点变量向量,$\{F\}$ 是外力向量。


    🧑🎓

    啊,原来是这样!连续体的控制方程是这种机制啊。


    单元技术

    🧑🎓

    "单元技术"我听过,但可能没理解透......


    单元类型次数节点数(3D)精度计算成本
    四面体1阶线性4低(剪切锁定)
    四面体2阶二次10
    六面体1阶线性8
    六面体2阶二次20非常高
    楔形线性/二次6/15中~高

    积分方案

    🧑🎓

    积分方案具体是什么意思呢?


    🎓
    • 完全积分:精确积分所有项。刚度高估倾向(锁定
    • 降低积分:减少积分点数。计算效率提高,但有沙漏模式发生风险
    • 选择性降低积分(B-bar法):分离体积项和偏差项进行积分。规避锁定

    • 🧑🎓

      听到这里终于明白为什么单元类型这么重要了!


      收敛性和稳定性

      🧑🎓

      如果收敛不了,首先应该查什么?


      🎓
      • h-细化:细分网格(减小单元尺寸h)提高精度
      • p-细化:提高单元多项式次数提高精度
      • hp-细化:同时最优化h和p

      • 🎓

        收敛速度:二次单元以 $O(h^2)$ 阶减少误差(光滑解)


        🧑🎓

        原来如此...细分网格看起来简单,其实深度很深啊。


        求解器设置建议

        🧑🎓

        具体用什么算法求解超级电容器建模呢?


        参数推荐值备注
        迭代法收敛判定$10^{-6}$残差范数准则
        预处理方法ILU(0) or AMG依问题规模
        最大迭代次数1000不收敛时重新检查设置
        内存模式In-core尽可能采用

        单片法

        在1个连立方程系统中同时求解所有物理场。对强耦合稳定,但实现复杂,内存消耗大。

        分离法(分离迭代法

        各物理场独立求解,在界面处进行数据交换。实现容易,可利用现有求解器。适合弱耦合。

        界面数据转移

        最近邻法(最简单但精度低)、射影法(保守的)、RBF插值(对网格非一致性强)。保守性和精度的均衡重要。

        子迭代

        各耦合步内进行充分迭代,确保界面条件的一致性。残差准则应按各物理场的典型值进行缩放。

        Aitken松弛

        自动调整耦合迭代的松弛系数。防止过松弛导致的发散,加快收敛的自适应方法。

        稳定性条件

        注意added mass效应(流体-结构耦合时结构密度≈流体密度)。不稳定时采用Robin型界面条件或IQN-ILS法。

        超级电容器建模实务应用

        🎓

        讲解超级电容器建模的实务分析流程和注意事项。



        分析流程

        🧑🎓

        从最初的一步开始讲吧!应该从哪开始呢?


        🎓

        1. 预处理(Pre-processing)

        • CAD数据导入和形状简化
        • 材料特性定义
        • 网格生成(单元类型、尺寸确定)
        • 边界条件和荷载条件设置

        🎓

        2. 求解(Solving)

        • 求解器设置(求解法、收敛准则、输出控制)
        • 作业投入和计算实行
        • 收敛监控

        🎓

        3. 后处理(Post-processing)

        • 结果可视化(位移、应力、其他物理量)
        • 结果验证和合理性确认
        • 报告编制


        网格生成最佳实践

        🧑🎓

        网格的好坏怎样判断呢?



        单元品质指标

        🧑🎓

        请给我讲讲"单元品质指标"!


        指标理想值允许范围影响
        长宽比1.0< 5.0精度低下
        Jacobian比1.0> 0.3单元退化
        翘曲< 15°精度低下
        歪度< 45°收敛性恶化
        锥度比0< 0.5精度低下

        网格密度的确定

        🧑🎓

        网格密度的确定具体是什么意思呢?


        🎓
        • 应力集中部:配置最少3层以上单元
        • 应力梯度大的区域:单元尺寸缩小到周围的1/3~1/5
        • 荷载印加点附近:局部细分
        • 远处区域:粗网格确保计算效率


        • 边界条件设置指南

          🧑🎓

          听说边界条件搞错了全盘皆输......


          🎓
          • 注意过约束:刚体移动约束仅限6个自由度
          • 利用对称条件:减小计算规模
          • 荷载等价分配:集中荷载 vs. 分布荷载的选择

          • 🧑🎓

            啊,原来是这样!过约束注意仕組。


            各商用工具的实现步骤

            🧑🎓

            有很多软件对吧?分别讲讲各自特点吧!


            工具名开发者/现在主要文件格式
            COMSOL MultiphysicsCOMSOL AB.mph
            Ansys Mechanical(旧ANSYS Structural)Ansys Inc..cdb, .rst, .db, .ans, .mac
            Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
            Ansys FluentAnsys Inc..cas, .dat, .msh, .jou

            COMSOL Multiphysics

            🧑🎓

            请给我讲讲"COMSOL Multiphysics"!


            🎓

            1986年在瑞典成立。作为MATLAB链接的FEMLAB开始,后改名为COMSOL。多物理是强项。

            现在的所属公司:COMSOL AB



            Ansys Mechanical(旧ANSYS Structural)

            🧑🎓

            请给我讲讲"Ansys Mechanical"!


            🎓

            1970年由Swanson Analysis Systems Inc. (SASI)开发。以APDL(Ansys参数设计语言)为基础。

            现在的所属公司:Ansys Inc.


            🧑🎓

            老师讲的很好理解!工具名的模糊理解消除了。


            常见失败和对策

            🧑🎓

            初心者容易犯什么错误?想提前知道!


            现象原因对策
            计算不收敛网格品质不良、不恰当的边界条件网格改善、拘束条件见解
            应力异常大应力奇点、网格依赖避免奇点、局部网格细分
            位移不现实材料常数错误、单位系不一致确认输入数据
            计算时间过大不必要的细分、低效求解法网格最优化、并行计算

            质量保证清单

            🧑🎓

            教科书上没有的"现场智慧"有吗?


            🎓
            • 网格收敛性用3个水平以上确认了吗
            • 力的平衡(反力合计)验证了吗
            • 结果在物理合理范围内确认了吗
            • 与已知的理论解或基准问题对比了吗



            • 🎓

              嗯,很好的势头!实际动手操作才是最好的学习。有不懂的地方随时问我。


              Coffee Break 闲话

              电车回生制动中"吸收"超级电容器的真实力量

              超级电容器最有冲击力的实际应用是铁路回生制动能量的吸收。电车制动时,电动机发挥发电机作用,瞬间产生数百kW的电力。如果用普通锂离子电池吸收这个电力,充放电速率太高,寿命会大幅缩短。但超级电容器即使在10C~100C以上的高速率下,劣化也几乎不发生。东京都交通局在都营地下铁引入EDLC,有效利用约30~40%的回生电力,年度用电量削减约10%的成绩。这样的系统设计是由超级电容器建模仿真来支撑的。

              超级电容器建模软件比较

              🎓

              详述超级电容器建模对应的主要商用CAE工具的功能对比和各产品的历史背景。



              支持工具列表

              🧑🎓

              那超级电容器建模可以用什么软件呢?


              工具名开发者/现在主要文件格式
              COMSOL MultiphysicsCOMSOL AB.mph
              Ansys Mechanical(旧ANSYS Structural)Ansys Inc..cdb, .rst, .db, .ans, .mac
              Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
              Ansys FluentAnsys Inc..cas, .dat, .msh, .jou

              COMSOL Multiphysics

              🧑🎓

              请给我讲讲"COMSOL Multiphysics"!


              🎓

              1986年在瑞典成立。作为MATLAB链接的FEMLAB开始,后改名为COMSOL。多物理是强项。

              现在的所属公司:COMSOL AB



              Ansys Mechanical(旧ANSYS Structural)

              🧑🎓

              请给我讲讲"Ansys Mechanical"!


              🎓

              1970年由Swanson Analysis Systems Inc. (SASI)开发。以APDL(Ansys参数设计语言)为基础。

              现在的所属公司:Ansys Inc.




              Simcenter STAR-CCM+

              🧑🎓

              接下来是Simcenter STAR吧。什么内容呢?


              🎓

              由CD-adapco开发。2016年被西门子收购,整合到Simcenter品牌。以多面体网格为特点。

              现在的所属公司:Siemens Digital Industries Software



              Ansys Fluent

              🧑🎓

              接下来是Ansys Fluent吧。什么内容呢?


              🎓

              由Fluent Inc.开发。2006年被Ansys收购。基于非结构网格的通用CFD求解器。

              现在的所属公司:Ansys Inc.


              🧑🎓

              原来如此...年在瑞典成立看似简单,其实深度很深啊。


              功能对比矩阵

              🧑🎓

              预算和时间都有限,哪个最划算呢?


              功能COMSOLAnsys MechanicalStar-CCM+Fluent
              基本功能
              高级功能
              自动化/脚本
              并行计算
              GPU支持

              转换时的风险

              🧑🎓

              转换时的风险具体是什么意思呢?


              🎓
              • 单元类型的不兼容:求解器专用单元无法用中立格式表示
              • 材料模型的差异:同名但内部实现不同的情况
              • 边界条件的重新定义:多数情况需要手动重新设置
              • 结果数据的对比:输出变量的定义不同(节点值 vs. 单元值、积分点值)

              • 🧑🎓

                啊,原来是这样!不同工具间的模型转换是这种机制啊。


                许可证形式

                🧑🎓

                "许可证形式"我听过,但可能没理解透......


                工具许可证特点
                商用FEA节点锁定/浮动高费用但有官方支持
                OpenFOAMGPL免费但支持收费
                COMSOL节点锁定/浮动按模块购买
                Code_AsterGPLEDF开发的开源求解器

                选择指南

                🧑🎓

                最终选哪个,判断标准请告诉我?


                🎓

                在超级电容器建模工具选定中应考虑以下因素:


                🎓
                • 分析规模:数万~数亿DOF的可扩展性
                • 物理模型:必要的本构模型、单元类型支持状况
                • 工作流:与CAD的协同、自动化容易度
                • 成本:初期投资+年度维护+教育成本
                • 支持:技术支持的质量和响应速度



                • 🎓

                  嗯,很好的势头!实际动手操作才是最好的学习。有不懂的地方随时问我。


                  Coffee Break 闲话

                  电容器专用工具vs汎用多物理——实际选择

                  市场上专门针对超级电容器建模的工具几乎不存在,很多工程师都使用COMSOL Multiphysics或ANSYS Fluent这样的汎用多物理软件。将电池模块专门定制为EDLC是常见做法,COMSOL的"Batteries & Fuel Cells"模块也支持电气双层电容器。而在研究水平,电化学向代码的実装FEniCS使自作実装多。工具选择的关键是"实验参数同定工作流有多高效"。因为超级电容器的特性随材料变化很大,现场倾向选择能快速循环执行「实验→参数同定→模拟→再实验」的工具。

                  超级电容器建模先进研究

                  🎓

                  看看超级电容器建模中最新的研究动向和先进方法。



                  最新的数值方法

                  🧑🎓

                  接下来是最新数值方法吧。什么内容呢?



                  🧑🎓

                  嗯...只看公式不太明白呢...表示什么呢?


                  🎓
                  • 等几何分析(IGA):直接使用NURBS基函数,实现CAD-CAE无缝协同
                  • 粒子法(SPH, MPM):无网格方法追踪大变形、破坏
                  • 位场法(Phase-Field):通过隐式表示界面进行复杂界面追踪
                  • 机器学习辅助代理模型、物理信息神经网络(PINN)


                  • 高性能计算(HPC)支持


                    并行化方法概述适用求解器
                    MPI(领域分割)分布式内存型。大规模问题标准所有主要求解器
                    OpenMP共享内存型。节点内并行很多求解器
                    GPU(CUDA/OpenCL)GPGPU利用。特别在显式法有效LS-DYNA,Fluent等
                    混合MPI+OpenMP节点间+节点内并行大规模HPC环境

                    超级电容器建模故障排除



                    常见错误和对策

                    🧑🎓

                    老师也对超级电容器建模过夜调试吗?(笑)



                    1.收敛失败

                    🧑🎓

                    收敛失败具体是什么意思呢?


                    🎓

                    现象:求解器在指定迭代次数内无法收敛,异常终止


                    🎓

                    可能原因

                    • 网格品质不足(过度扭曲单元)
                    • 材料参数设置不恰当
                    • 不恰当的初始条件
                    • 非线性性太强(缺少荷载步)

                    🎓

                    对策

                    • 进行网格品质检查(长宽比、Jacobian)
                    • 确认材料参数的单位系
                    • 将荷载分割为多个步(增加子步数)
                    • 放宽收敛判定准则(但要注意精度)

                    🧑🎓

                    也就是说,在收敛失败那部分偷工减料,后面会吃苦头对吧?记住了!



                    2.非物理结果

                    🧑🎓

                    接下来是非物理结果吧。什么内容呢?


                    🎓

                    现象应力/位移/温度等出现物理上不现实的值


                    🎓

                    可能原因

                    • 边界条件设置误
                    • 单位系混用(SI单位与工程单位混用)
                    • 不恰当的单元类型选择
                    • 应力奇点存在

                    🎓

                    对策

                    • 确认反力合计(力的平衡)
                    • 检查单位系的一致性
                    • 重新考虑单元类型的恰当性
                    • 除去奇点或采用子建模

                    🧑🎓

                    前辈说"收敛失败务必做好"这话的意思现在明白了。




                    3.计算时间超过

                    🧑🎓

                    计算时间超过具体是什么意思呢?


                    🎓

                    现象:计算用时是预期时间的好几倍


                    🎓

                    对策

                    • 优化网格粗细分布
                    • 利用对称性(1/2、1/4模型)
                    • 最优化求解器设置(迭代法、预处理选择)
                    • 利用并行计算



                    4.内存不足

                    🧑🎓

                    请给我讲讲"内存不足"!


                    🎓

                    现象:Out of Memory错误


                    🧑🎓

                    前辈说"收敛失败务必做好"这话的意思现在明白了。


                    🎓

                    对策

                    • 采用核外求解法
                    • 减小网格规模
                    • 确认使用64位版求解器
                    • 增加内存分配

                    🧑🎓

                    哦~,收敛失败的话题超有意思!请多讲一些。


                    Nastran代表性错误

                    🧑🎓

                    代表性错误具体是什么意思呢?


                    🎓
                    • FATAL 2012:奇异刚度矩阵 → 重新检查约束条件
                    • USER WARNING 5291:单元品质不良 → 修正网格
                    • SYSTEM FATAL 3008:内存不足 → 调整MEM设置


                    • Abaqus代表性错误

                      🧑🎓

                      请给我讲讲"代表性错误"!


                      🎓
                      • Excessive distortion:单元过度变形 → 确认NLGEOM、改善网格
                      • Zero pivot:约束不足 → 追加边界条件
                      • Time increment too small