简介本资源是一份面向电化学仿真初学者与电池研发工程师的COMSOL多物理场建模入门指南聚焦锂离子电池准二维P2D电化学模型的理论构建与软件实现路径。内容系统梳理Doyle经典模型的三大核心过程——固相锂离子扩散Fick第二定律、颗粒表面电化学反应Butler–Volmer方程及电解液中离子传质Nernst–Planck方程并详解其物理假设、控制方程推导、边界条件设定及COMSOL有限元求解逻辑助力读者从电化学原理出发理解仿真本质。资源为单文件PDF共1个9.07MB文档结构清晰含背景分析、数学建模、方程物理意义解读与参数温度/荷电态依赖性可视化图示适合作为仿真实践前的理论铺垫与建模参考。已有4455人学习下载是COMSOL电化学方向少有的兼顾理论深度与工程落地性的中文技术解析材料。1. 这不是“点几下就能出图”的COMSOL教程它是一份能让你在3天内从扣电数据推导出P2D模型全部参数的锂电仿真落地笔记你是不是也试过打开COMSOL锂离子电池接口填完几个默认参数一运行——报错“非线性求解器不收敛”查文档发现要设初始电势但不知道正极初始固相电势该填3.6还是3.8看到“Bruggeman系数”就卡住翻遍帮助手册也没说负极和隔膜的值到底差多少更别说把实验室测的扣电OCV曲线、DSC得到的dU/dT、恒流间歇滴定GITT反演的Ds(y,T)全塞进同一个模型里——最后干脆放弃回归“调参-跑-看温度场-再调”成了玄学仿真。这篇PDF不是概念科普也不是软件操作说明书。它是某高校电化学仿真组内部流传的实操手记完整复现了从一份磷酸铁锂/石墨体系扣电测试原始数据出发到在COMSOL中构建可收敛、可验证、可映射真实工况的准二维P2D模型的全过程。它拆解了6个核心参数模块的来源、换算逻辑、典型取值范围与交叉验证方法尤其聚焦那些文献里不会写、COMSOL默认不给、但一填错就直接翻车的关键数值比如负极最小有效固相浓度怎么算才不导致放电末期浓度崩塌比如液相扩散系数随浓度变化的拐点在哪为什么填1200 mol/m³时必须同步修正Bruggeman指数比如为什么用Butler-Volmer方程时av活性比表面积不能按球形颗粒理论值硬套而必须结合极片孔隙率反推。如果你正卡在“模型建出来了但结果和实测电压差200mV”“低温仿真发散但室温正常”“浓差极化预测值比实测高一倍”这些具体问题上这篇笔记就是为你写的——它不讲“多物理场耦合有多牛”只告诉你哪一行参数改0.1模型就从发散变收敛哪个边界条件漏设仿真结果就系统性偏移30%。2. P2D模型三大物理过程从控制方程到COMSOL接口变量的硬核映射COMSOL锂离子电池接口Lithium-Ion Battery, LIMB本质是P2D模型的工程实现封装。它把原本分散在9个偏微分方程里的物理量压缩成3个核心因变量dependent variables固相电势 φ_s、液相电势 φ_e、电解质盐浓度 c_e。所有其他输出——局部电流密度、产热功率、SOC分布、过电势η——都是这三个变量的函数。理解它们如何从方程落地为COMSOL中的具体设置是避免“黑匣子式建模”的第一道门槛。2.1 固相锂离子扩散Fick第二定律在COMSOL中的三重约束正负极活性颗粒内部的锂离子固相扩散由Fick第二定律控制$$ \frac{\partial c_s}{\partial t} D_{s} \left( \frac{\partial^2 c_s}{\partial r^2} \frac{2}{r} \frac{\partial c_s}{\partial r} \right) $$其中 $c_s$ 是固相锂离子浓度mol/m³$D_s$ 是固相扩散系数m²/s$r$ 是颗粒半径方向坐标。这个方程在COMSOL中并非直接输入而是通过LIMB接口的固相扩散域Solid-phase Diffusion Domain自动离散。关键在于三个必须显式设置的环节初始浓度 $c_{s,0}$ 的物理意义与赋值逻辑COMSOL要求为每个电极域指定初始固相浓度。这不是一个随意值而是满电态SOC100%下活性材料的实际嵌锂浓度。以石墨负极为例其理论最大嵌锂浓度 $c_{s,\text{max}}$ 约为 31000 mol/m³对应LiC₆。但实测扣电放电容量 $C_{\text{disch}}$ 为 350 mAh/g材料真密度 $\rho_{\text{true}}$ 为 2.2 g/cm³则初始浓度计算为# Python计算示例负极初始固相浓度 C_disch_neg 350 / 1000 # 转为Ah/g Mw_neg 74.035 # LiC6分子量g/mol cs_max_neg 31000 # mol/m³ rho_true_neg 2.2 * 1e6 # g/m³ c_s0_neg C_disch_neg * cs_max_neg * Mw_neg / (3600 * rho_true_neg) # 单位mol/m³ print(f负极初始固相浓度 c_s0_neg ≈ {c_s0_neg:.0f} mol/m³) # 输出负极初始固相浓度 c_s0_neg ≈ 10200 mol/m³参数说明3600是将Ah转换为As库仑的因子rho_true_neg必须用g/m³单位否则数量级错10⁶倍。这是新手最常踩的坑——用g/cm³直接代入导致浓度低6个数量级模型完全失真。边界条件颗粒表面通量 $J_s$ 的物理闭环颗粒表面的固相扩散通量 $J_s$ 必须等于电化学反应消耗/生成的锂离子通量即 $$ J_s -D_s \left. \frac{\partial c_s}{\partial r} \right|_{rR} a_v F i_0 \left[ \exp\left(\frac{\alpha_a F \eta}{RT}\right) - \exp\left(-\frac{\alpha_c F \eta}{RT}\right) \right] $$ 在COMSOL中这个等式不是手动输入的而是由LIMB接口的电化学反应域Electrochemical Reaction Domain自动耦合。但前提是你必须正确设置a_v活性比表面积和i_0交换电流密度。a_v的单位是 m⁻¹其物理含义是单位体积电极内活性颗粒的总表面积。对于球形颗粒理论值为 $a_v 3(1-\varepsilon_s)/R_p$其中 $\varepsilon_s$ 是固相体积分数$R_p$ 是颗粒半径。但实际电极中颗粒有团聚、涂层不均因此a_v通常需下调20%~40%才能匹配实测极化曲线。血泪经验若仿真极化电压比实测小50mV优先检查a_v是否过高。$D_s$ 的温度与荷电态依赖性如何避免收敛崩溃文献给出的 $D_s$ 常是单一温度下的常数但P2D模型要求其为 $D_s(c_s, T)$ 函数。COMSOL支持通过“定义函数插值”或“解析”功能输入表达式。以磷酸铁锂正极为例其 $D_s$ 随SOC降低y减小而增大典型表达式为 $$ D_{s,\text{pos}} 1.18 \times 10^{-13} \exp\left[-\frac{35000}{R}\left(\frac{1}{T} - \frac{1}{298.15}\right)\right] \cdot (1-y)^{0.5} $$ 注意1.18e-13是关键数量级若误用1.18e-10常见文献笔误模型会在放电中期因固相扩散过快而出现非物理解如浓度负值。实操技巧先用常数 $D_s1e-14$ 跑通基础模型再逐步替换为函数每次替换后观察c_s场是否平滑无振荡。2.2 电化学反应Butler-Volmer方程在COMSOL中的参数陷阱电极/电解质界面的电荷转移反应由Butler-Volmer方程描述 $$ i_L i_0 \left[ \exp\left(\frac{\alpha_a F \eta}{RT}\right) - \exp\left(-\frac{\alpha_c F \eta}{RT}\right) \right] $$ 其中 $i_L$ 是局部电流密度A/m²$\eta \phi_s - \phi_e - E_{\text{eq}}$ 是过电势。COMSOL将此方程内置于“电化学反应”边界条件中但以下四个参数的设置直接决定模型能否收敛及精度参数COMSOL位置典型取值关键陷阱交换电流密度 $i_0$电化学反应域 交换电流密度正极1~10 A/m²负极5~50 A/m²$i_0$ 必须与 $c_s$、$c_e$ 同步变化。若设为常数低温下会因 $i_0$ 过小导致迭代步长爆炸式增长最终不收敛。必须用 $i_0 a_v F k_0 c_s^{\alpha_a} c_e^{\alpha_c}$ 形式定义。传递系数 $\alpha_a, \alpha_c$电化学反应域 传递系数通常取0.5对称反应石墨负极脱锂 $\alpha_a≈0.3$若 $\alpha_a \alpha_c ≠ 1$COMSOL会自动归一化但可能导致 $i_0$ 解析式失效。务必检查“物理场设置”中是否勾选“使用指定传递系数”。平衡电势 $E_{\text{eq}}$电化学反应域 平衡电势SOC-OCV插值表至少20个点OCV曲线在平台区如LFP的3.2V斜率接近0COMSOL数值微分易失稳。必须在平台区加密采样点每0.02SOC一个点否则放电中段电压预测偏差100mV。活性比表面积 $a_v$电化学反应域 活性比表面积见2.1节计算逻辑$a_v$ 与电极厚度、孔隙率强相关。若极片涂覆面密度为18 mg/cm²孔隙率40%则 $a_v ≈ 1.2e5$ m⁻¹若忽略孔隙率按密实颗粒算会高估3倍。提示在COMSOL中i_0和a_v的单位必须严格匹配。i_0输入框默认单位是 A/m²但若你在“材料”节点里定义了a_v为 m⁻¹则i_0的表达式中必须显式包含a_v否则COMSOL会将其视为常数失去物理关联。2.3 电解质传质Nernst-Planck方程与Bruggeman修正的协同设置电解质中锂离子的迁移与扩散由Nernst-Planck方程描述 $$ \frac{\partial c_e}{\partial t} \nabla \cdot \left[ D_e^{\text{eff}} \nabla c_e - t_^0 c_e \nabla \phi_e \right] $$ 其中 $D_e^{\text{eff}}$ 是有效液相扩散系数$t_^0$ 是阳离子迁移数。COMSOL通过“电解质域Electrolyte Domain”实现此方程但两个关键修正必须手动介入Bruggeman系数 $\varepsilon^{b}$ 的分区设置多孔电极正极、负极与隔膜的孔隙结构差异巨大导致 $D_e^{\text{eff}} \varepsilon_e^{b} D_e$ 中的指数 $b$ 不同。典型取值隔膜$b 1.5$孔道直且均匀正极NCM$b 2.0$孔道弯曲度高负极石墨$b 2.2$颗粒堆叠致密避坑COMSOL默认对所有域使用同一 $b$ 值。必须在“材料电解质扩散系数”中为正极、负极、隔膜分别创建独立的“材料属性”并为每个属性单独设置BruggemanExponent。若共用浓差极化预测误差可达200%。迁移数 $t_^0$ 的温度敏感性迁移数并非常数其随温度升高而增大因离子迁移率提升快于溶剂粘度下降。常用经验公式 $$ t_^0(T) t_^0(298) \cdot \exp\left[1200 \left( \frac{1}{298} - \frac{1}{T} \right) \right] $$ 其中 $t_^0(298) ≈ 0.38$1M LiPF₆ in EC:DMC。若设为常数在45°C仿真中液相电势梯度会被低估导致欧姆压降预测偏低15%。3. 从扣电到模型6类核心参数的来源、换算与交叉验证全流程P2D模型的成败80%取决于参数的真实性与一致性。本节将PDF中零散的参数推导逻辑整合为一条可执行的流水线以一份标准扣式电池CR2032的测试数据为起点产出COMSOL中全部必需参数并提供3种交叉验证方法确保参数链无断裂。3.1 设计参数从极片图纸到COMSOL几何域的精确映射设计参数是模型的“骨架尺寸”直接决定几何域大小与体积分数。其来源是电池设计BOM而非实验测量。参数获取方式COMSOL中用途计算示例正极涂覆面密度 $\rho_{\text{coat}}$极片称重法$(m_{\text{coated}} - m_{\text{foil}})/A$计算固相体积分数 $\varepsilon_s$$\rho_{\text{coat}} 22$ mg/cm² 220 kg/m²活性物质含量 $\omega_m$浆料配方单如NCM811占96%同上$\omega_m 0.96$集流体厚度 $t_{\text{foil}}$供应商规格书Al箔12μm计算单层厚度 $t_{\text{active}}$$t_{\text{foil}} 12$ μm极片总厚度 $t_{\text{total}}$千分尺实测含双面涂层同上$t_{\text{total}} 75$ μm → $t_{\text{active}} (75-12)/2 31.5$ μm活性物质真密度 $\rho_{\text{true}}$文献或厂商提供NCM811≈4.7 g/cm³同上$\rho_{\text{true}} 4.7 \times 10^6$ g/m³电解液体积分数 $\varepsilon_e$由 $\varepsilon_e 1 - \varepsilon_s - \varepsilon_{\text{binder}} - \varepsilon_{\text{cond}}$ 推导定义电解质域渗透率$\varepsilon_s \frac{220 \times 0.96}{31.5 \times 4.7e6} ≈ 0.31$ → $\varepsilon_e ≈ 0.45$假设粘结剂导电剂占0.24关键逻辑$\varepsilon_e$ 不是直接测量值而是由所有固相组分体积分数之和被1减去得到。若 $\varepsilon_s$ 算错0.05$\varepsilon_e$ 就错0.05直接导致电解质电导率偏差20%。实操中我每次都会用Excel表格列出所有组分活性物质、PVDF、SP炭黑、Al箔逐项计算体积最后求和验证是否≈1。3.2 材料参数电化学与传质参数的三层验证法材料参数是模型的“血肉”其不确定性最高。PDF中提到的 $D_s$、$i_0$、$E_{\text{eq}}$ 等必须通过“实验测量→文献校准→模型反演”三层过滤。第一层实验可测参数直接输入平衡电势 $E_{\text{eq}}(SOC)$用0.05C小电流恒流放电静置2h测OCV获得至少20个SOC点的电压值。导入COMSOL时选择“插值”函数外插模式设为“最近邻”避免SOC0或1时外推失真。熵热系数 $dE/dT$通过等温量热法ITC或电压温度系数法测得。PDF中强调其对小倍率吸热现象的关键作用。在COMSOL中该参数用于计算可逆热源 $Q_{\text{rev}} i_L T (dE/dT)$。若缺失模型将无法预测低温充电初期的温度下降。第二层文献校准参数带范围约束固相扩散系数 $D_s$查《Journal of The Electrochemical Society》中同体系材料的GITT或PITT数据。重点看其数量级必为10⁻¹³~10⁻¹⁵ m²/s和温度活化能Ea通常30~50 kJ/mol。若文献给的是Arrhenius形式 $D_s D_0 \exp(-E_a/RT)$则 $D_0$ 必须与 $E_a$ 匹配不可拼凑。液相扩散系数 $D_e$参考《Electrochimica Acta》中1M LiPF₆电解液的NMR或Taylor dispersion数据。注意其浓度依赖性——PDF指出在1000~1200 mol/m³区间$D_e$ 随浓度升高而增大故仿真中必须用浓度函数而非常数。第三层模型反演参数闭环验证交换电流密度 $i_0$这是唯一必须通过模型反演确定的参数。步骤用已知 $E_{\text{eq}}$、$D_s$、$D_e$ 构建基础模型施加0.2C恒流放电记录端电压曲线在COMSOL中将 $i_0$ 设为参数扫描变量范围1~20 A/m²运行参数化扫描提取各 $i_0$ 下的电压曲线与实测曲线计算RMSE取RMSE最小时的 $i_0$ 值。我一般会同时扫描 $a_v$因二者强耦合形成2D参数空间。交叉验证法1极化电压分解COMSOL可分别输出欧姆压降$i_L R_{\text{ohm}}$、活化过电势$\eta_{\text{act}}$、浓差过电势$\eta_{\text{conc}}$。实测极化电压 三者之和。若仿真中 $\eta_{\text{act}}$ 占比70%说明 $i_0$ 过小若 $\eta_{\text{conc}}$ 占比50%说明 $D_e$ 或 $b$ 过小。交叉验证法2SOC分布一致性在放电末期SOC0.1COMSOL显示负极颗粒中心 $c_s$ 已降至500 mol/m³但实测XRD显示仍有大量LiC₁₂存在。此时必须下调 $D_s$ 或上调 $a_v$使固相扩散更慢保留更多锂在颗粒中心。交叉验证法3产热功率标定用红外热像仪测得电池表面温升速率 $dT/dt$结合比热容 $C_p$ 和质量 $m$计算实测产热功率 $Q_{\text{meas}} m C_p dT/dt$。COMSOL输出的 $Q_{\text{total}} Q_{\text{ohm}} Q_{\text{rev}} Q_{\text{irrev}}$ 应与 $Q_{\text{meas}}$ 在±15%内吻合。若 $Q_{\text{total}}$ 偏低优先检查 $dE/dT$ 符号放电时应为负产生吸热是否设反。3.3 关键衍生参数PDF中未明说但决定模型生死的4个隐性变量PDF详细推导了 $c_{s,0}$、$c_{s,\text{min/max}}$但有4个衍生参数虽未单列却在COMSOL中无处不在正极最大有效固相浓度 $c_{s,\text{max,pos}}$公式$c_{s,\text{max,pos}} c_{s,0,\text{pos}} \frac{C_{\text{cell}} \cdot 3600 \cdot \rho_{\text{true,pos}}}{M_w \cdot \omega_m \cdot t_{\text{active,pos}}}$其中 $C_{\text{cell}}$ 是电池设计比容量Ah/kg。意义它定义了正极放电截止时的最高嵌锂量。若设错放电末期正极 $c_s$ 会超限报错。负极最小有效固相浓度 $c_{s,\text{min,neg}}$PDF公式复杂但核心是$c_{s,\text{min,neg}} c_{s,0,\text{neg}} - \Delta c_{s,\text{loss}}$其中 $\Delta c_{s,\text{loss}}$ 由SEI消耗和正负极容量不平衡决定。实操中我直接用扣电首效如92%折算$c_{s,\text{min,neg}} c_{s,0,\text{neg}} \times (1 - 0.08)$。Bruggeman系数 $b$ 的实验标定法若无文献值可用EIS拟合在高频区1kHz~100kHzWarburg阻抗斜率与 $b$ 直接相关。拟合公式 $Z_w \sigma_w \omega^{-0.5} \tan\beta$其中 $\beta \pi/(2b)$。测得 $\beta0.7$则 $b \pi/(2 \times 0.7) ≈ 2.25$。集流体电导率 $\sigma_{\text{foil}}$ 的真实值PDF假设“无限大”但COMSOL中必须输入。Al箔实测电导率约3.5e7 S/m非纯铝的3.7e7Cu箔约5.8e7 S/m。若用纯金属值会导致集流体压降被低估端电压预测偏高50mV。4. 避坑指南COMSOL锂电仿真中5条血泪换来的收敛与精度保障法则参数填完了几何画好了物理场设定了——点击“计算”然后等待5分钟弹出红色报错“找不到满足相对容差的解”“雅可比矩阵奇异”“时间步长过小”。别急这几乎是每个锂电仿真新手的必经之路。以下是我在某公司电池研发部三年实操中用报废的27个模型文件总结出的5条硬核避坑法则每一条都对应一个具体现象、根本原因和可立即执行的解决方案。4.1 现象求解器在t0.1s附近停滞反复尝试更小时间步长最终失败原因初始条件不满足电化学平衡。COMSOL在t0时要求所有域的 $i_L0$即 $\eta0$。若你设的初始 $\phi_s$ 和 $\phi_e$ 未使 $\phi_s - \phi_e E_{\text{eq}}(SOC_0)$则初始时刻存在巨大伪电流导致瞬态求解崩溃。解决在“研究稳态”中先运行一次稳态求解仅电化学让系统自动调整 $\phi_s$ 和 $\phi_e$ 至平衡态将稳态结果作为瞬态研究的初始值右键“瞬态”“更多初始值”选择稳态解关键操作在稳态研究中“电化学反应”边界条件必须启用“计算初始平衡电势”否则稳态解无效。4.2 现象放电曲线前10%SOC电压骤降200mV之后平缓与实测严重不符原因负极最小有效固相浓度 $c_{s,\text{min,neg}}$ 设置过低导致放电初期负极颗粒表面 $c_s$ 快速耗尽$i_0$ 急剧下降活化极化暴增。解决重新计算 $c_{s,\text{min,neg}}$确保其不低于 $c_{s,0,\text{neg}} \times 0.1$即保留10%锂不参与循环在COMSOL中为负极固相扩散域添加“浓度限制”右键“固相扩散”“更次浓度”设置 $c_s$ 下限为 $c_{s,\text{min,neg}}$验证查看放电初期负极表面 $c_s$ 场应平滑下降无突变或负值。4.3 现象45°C仿真收敛25°C不收敛或低温下电压预测比实测高150mV原因$i_0$ 和 $D_s$ 的温度依赖性未同步。若 $i_0$ 用了Arrhenius公式而 $D_s$ 仍为常数则低温下 $i_0$ 衰减更快反应动力学成为瓶颈求解器无法协调。解决统一所有温度敏感参数的活化能$i_0$ 的 $E_{a,i0}$ 与 $D_s$ 的 $E_{a,Ds}$ 应相近通常35~45 kJ/mol在COMSOL中用同一温度变量T定义所有参数i0 i0_ref * exp(-Ea_i0/R*(1/T-1/298))Ds Ds_ref * exp(-Ea_Ds/R*(1/T-1/298))技巧将Ea_i0和Ea_Ds设为相同参数名强制关联。4.4 现象浓差极化 $\eta_{\text{conc}}$ 在放电中期达300mV但实测EIS显示浓差阻抗50mΩ原因Bruggeman系数 $b$ 过大或电解液体积分数 $\varepsilon_e$ 过小导致 $D_e^{\text{eff}}$ 被过度削弱。解决检查 $\varepsilon_e$ 计算是否遗漏了导电剂体积SP炭黑密度仅1.8 g/cm³易被低估将 $b$ 从2.2下调至1.8重新运行快速验证在COMSOL中右键“电解质”“派生值”“表面积分”计算整个电解质域的平均 $c_e$。若放电末期平均 $c_e$ 低于800 mol/m³初始1200说明浓差极化过强需调 $b$ 或 $\varepsilon_e$。4.5 现象网格细化后结果剧烈震荡或不同网格数下电压预测相差50mV原因几何域划分不合理导致关键界面如电极/电解质网格过渡不连续。COMSOL在界面处用“薄层”近似若网格在界面两侧尺寸差异3倍弱形式积分失准。解决在“几何”中为正极/电解质、负极/电解质界面分别创建“分割面”在“网格”中对每个界面应用“边界层网格”层数≥3总厚度2×颗粒半径对电极域使用“自由四面体”对电解质域使用“映射六面体”若为规则几何确保界面网格匹配。终极技巧运行前在“研究瞬态”设置中将“时间步进”从“自由”改为“严格”并设置最大步长为0.1s——这能强制求解器在关键动态期如SOC0.5精细采样避免跳过极化转折点。5. 进阶验证用3种物理实验数据反向标定模型精度的实战方法模型跑通只是起点真正的价值在于它能否成为“数字孪生体”精准映射真实电池的行为。PDF中强调了仿真对研发的降本增效作用但没说清楚如何证明你的模型真的准。这里分享我在某车企电池包热管理项目中验证COMSOL模型的3种硬核方法每一种都对应一个具体实验、一个COMSOL后处理操作、一个量化验收标准。5.1 方法一HPPC脉冲实验——验证极化响应的毫秒级精度实验HPPCHybrid Pulse Power Characterization测试。在SOC0.5时施加10s 1C放电脉冲 40s静置 10s 1C充电脉冲记录端电压瞬态响应。COMSOL验证在模型中将电流载荷设为分段函数if(t10, -I0, if(t50, 0, I0))运行瞬态研究输出“电池端子”处的电压V_term使用“结果派生值最大值”提取脉冲结束时刻t10s的电压V_drop和静置结束时刻t50s的恢复电压V_rec验收标准|V_drop,sim - V_drop,exp| 15 mV|V_rec,sim - V_rec,exp| 5 mV为什么严苛V_drop主要反映欧姆活化极化V_rec反映浓差极化恢复。若V_drop偏差大说明 $i_0$ 或 $\sigma_{\text{foil}}$ 错若V_rec偏差大说明 $D_e$ 或 $b$ 错。我曾因V_rec偏差8mV发现隔膜Bruggeman系数少设了0.2修正后整车NEDC续航预测误差从7%降至1.2%。5.2 方法二等温量热实验——验证产热机制的物理真实性实验在25°C恒温箱中对电池施加0.5C恒流充放电用高精度量热仪如TA Discovery实时测量热流率 $dQ/dt$。COMSOL验证在模型中启用“能量守恒”物理场添加“热源”Q_total Q_ohm Q_rev Q_irrev其中Q_ohm i_L^2 * R_ohmQ_rev i_L * T * dE/dTQ_irrev i_L * eta_act i_L * eta_conc输出整个电池域的“体积积分”热功率Q_integrated将Q_integrated与实测dQ/dt在同一时间轴上绘图对比。验收标准充电阶段Q_sim与Q_exp的RMSE 0.15 W放电阶段Q_sim与Q_exp的RMSE 本文还有配套的精品资源点击获取
