锂离子电池固态扩散阻抗的Voigt建模与参数辨识.pdf
《锂离子电池固态扩散阻抗的Voigt建模与参数辨识.pdf》由会员分享,可在线阅读,更多相关《锂离子电池固态扩散阻抗的Voigt建模与参数辨识.pdf(7页珍藏版)》请在咨信网上搜索。
1、2023.7Vol.47No.7研 究 与 设 计收稿日期:2023-01-25作者简介:刘昆朋(1997),男,山东省人,硕士,主要研究方向为锂离子电池。通信作者:方小红锂离子电池固态扩散阻抗的Voigt建模与参数辨识刘昆朋1,2,朱广焱3,张鹏博3,方小红1,2(1.中国科学院 上海高等研究院,上海 201210;2.中国科学院大学,北京 100049;3.上海派能能源科技股份有限公司,上海 201210)摘要:固态扩散阻抗表征锂离子电池电极性能,为测试和改进电极材料提供了便利。将可兼顾分析电极参数影响和模拟过电压且具有物理意义的Voigt模型从表示正极平面、球形颗粒阻抗推广至圆柱形颗粒阻
2、抗模型用于石墨等负极材料,将颗粒粒径等电极参数与模型参数相关联,建立模型过电压表达式,分析电极参数对电压的影响,并引入具有全局搜索方式的花朵授粉算法将表达式拟合电压测量值,得到模型参数的最优解。进一步地将该模型用于商用磷酸铁锂电池的测试分析,预测电极在各类工况下的固态扩散过电压,为锂离子电池电极材料的开发和性能提升提供技术支撑。关键词:Voigt模型;固态扩散阻抗;锂离子电池中图分类号:TM 912文献标识码:A文章编号:1002-087 X(2023)07-0878-07DOI:10.3969/j.issn.1002-087X.2023.07.012Voigt modeling and pa
3、rameter identification of solid-state diffusionimpedance in lithium-ion batteriesLIU Kunpeng1,2,ZHU Guangyan3,ZHANG Pengbo3,FANG Xiaohong1,2(1.Shanghai Advanced Research Institute,Chinese Academy of Sciences,Shanghai 201210,China;2.University of Chinese Academy ofSciences,Beijing 100049,China;3.Pylo
4、n Technologies,Co.,Ltd.,Shanghai 201210,China)Abstract:Solid-state diffusion impedance characterizes the performance of lithium-ion battery electrodes,whichprovides convenience for testing and improving materials.The Voigt model,which could both analyze the influence ofelectrode parameters and simul
5、ate overvoltage and had physical significance,was extended from representing thepositive planar and spherical particle impedance to the cylindrical particle impedance model for anode materials suchas graphite.By associating electrode parameters such as particle size with model parameters,model overv
6、oltageexpression was established to analyze the influence of electrode parameters on voltage,and a flower pollinationalgorithm with global search method was introduced to fit the expression to voltage measurement values,and theoptimal solution of model parameters was obtained.Further,the model was a
7、pplied to the test and analysis ofcommercial lithium iron phosphate batteries to predict the solid-state diffusion overvoltage of electrodes under variousworking conditions,which provided technical support for the development and performance improvement of lithiumion battery electrode materials.Key
8、words:Voigt model;solid-state diffusion impedance;lithium-ion battery锂离子电池(lithium-ion battery,LIB)具有能量密度大,循环寿命长等诸多优点,在消费电子、能源存储等领域得到广泛应用1-2。LIB的性能与电极中锂离子固态扩散过程密切相关2,研究该过程对提升 LIB 性能有重要作用。固态扩散阻抗是常用研究对象3,从阻抗频率特性中可得出固态扩散过程的动力学信息,能够用于分析颗粒粒径、锂离子固态扩散系数等电极参数对 LIB性能的影响3-4,测量该阻抗对新型电极材料的开发与应用具有指导意义。阻抗测量方法包括频域
9、测量方法和时域测量方法(timedomain measurement,TDM)5-6。频域测量低频的固态扩散阻抗耗时较长,使LIB荷电状态(state of charge,SOC)改变降低测量准确性5。TDM对LIB施加激励信号和测量响应信号并分别做时-频域变换6,耗时较短更适合于测量固态扩散阻抗。常用于固态扩散阻抗的激励信号为脉冲电流信号7-9,其弛豫电压响应信号中包含正负极固态扩散过程和电解液中扩散过程的过电压,但无法直接从总阻抗中分离各过程阻抗以分别分析动力学等方面的信息9,还需要借助其他分析手段。目 前 常 用 分 析 手 段 是 构 建 等 效 电 路 模 型(equivalentc
10、ircuit model,ECM)拟合低频阻抗,由模型各元件参数分析各过程对 LIB 性能的影响4,10。常用表示固态扩散阻抗的元件包括RC元件和Warburg元件等,但RC元件不具备物理意义,无法从总过电压中识别固态扩散过程过电压和分析电极材料对该电压的影响10-11;具有物理意义的Warburg元件阻抗表达式较复杂,难以对其做频域-时域变换并与电流表达式相乘8782023.7Vol.47No.7研 究 与 设 计计算过电压,无法量化电极材料参数改变引起的过电压变化3,12。故上述模型在分析固态扩散过程对 LIB 性能的影响上存在局限性,难以为改进电极材料提供借鉴。Voigt模型9解决了其他
11、元件无法兼顾分析电极参数影响和模拟过电压的问题,建立其阻抗表达式需要对菲克扩散方程作拉普拉斯变换3,9,因此表达式中各模型参数均与扩散距离(颗粒粒径)、扩散系数等电极参数相关联,且阻抗表达式相对简单易于做频域-时域变换9,模拟不同参数电极材料的固态扩散过电压。目前 Voigt模型已应用于镍钴锰酸锂电池9,被用于分析 SOC、温度对固态扩散过电压的影响,以预测电池在各类工况下性能,根据需求有针对性优化材料参数;也用于分析电池老化过程及预测老化电池电压。目前已建立平面、球形颗粒阻抗模型用于具有该形状颗粒的镍钴锰酸锂等多数种类 LIB的正极,因均满足正极过电压远大于石墨负极这一条件而忽略负极。但应用
12、于磷酸铁锂电池时13,因其正负极过电压较接近,在引入平面颗粒阻抗模型用于磷酸铁锂正极的同时尚需新建适用于石墨负极的圆柱形颗粒阻抗模型4。辨识 Voigt模型参数需要基于模型建立弛豫电压表达式并拟合测量值9,常用拟合算法为信赖域反射(trust-region-reflective,TRR)算法,搜索方式为局部搜索,拟合非线性电压数据无法得到最优解14。花朵授粉算法(flower pollinationalgorithm,FPA)可全局搜索15-16,相较于信赖域反射算法有较大的搜索空间增强找到最优解的能力15,更适合于辨识Voigt模型参数。本文将基于 Voigt模型建立正负极固态扩散过电压表达
13、式,应用 FPA 从弛豫电压中辨识模型参数,将其应用于已建立模型不适用的商用磷酸铁锂电池,拟首次准确预测该电池在不同工况下固态扩散过电压并分析电极参数对该过电压的影响。1 研究方法1.1 基于Voigt模型建立正负极固态扩散过电压表达式Voigt模型表示单个电极颗粒的固态扩散阻抗,如图 1(a)所示,扩散电容 Clim表示锂离子在颗粒不可渗透边界处的扩散阻抗,表达式与电极颗粒形状有关3,一定数目 RC元件表示基于有限体积法将颗粒离散化后各体积段的扩散阻抗9,其系数 akbk(k=1,2,.n)与体积段数目、电极颗粒形状有关9。模型参数为极化电阻 RD和时间常数 D,与电极材料各参数相关联,如式
14、(1)、(2)所示3:RD=-crADZF(1)D=r2D(2)式中:-c为颗粒表面电位随锂离子浓度变化的偏导数;r为颗粒半径;A为颗粒活性表面积;D为锂离子扩散系数;Z为电极反应时单个分子转移电子数;F为法拉第常数。为建立正负极固态扩散过电压表达式,基于有限体积法将圆柱形颗粒离散化,基于菲克扩散方程计算表示各体积段阻抗的 RC 元件的系数,建立表示该颗粒阻抗的 Voigt 模型,以用于与平面、球形颗粒 Voigt模型分别表示 LIB负极、正极颗粒固态扩散阻抗。电极固态扩散阻抗是大量不同粒径颗粒阻抗之和17-18,因 颗 粒 间 存 在 补 偿 固 态 扩 散 阻 抗 差 异 的 均 匀 化(
15、homogenization)扩散过程,可通过串联Voigt模型和表示均匀化过程阻抗的 RC元件来表示9,如图1(b)所示。在参数识别时可通过homo数值更接近同电极的D9这一假设区分homo。对 ECM 频域阻抗表达式作拉普拉斯逆变换得到时域表达 式,Voigt 模 型 表 达式 ZVoigt(t)为k=1nakRDbkDe-tbkD(t),其中(t)为阶跃函数;RC 元件表达式 Zhomo(t)为Rhomohomoe-thomo(t),其中 Rhomo、homo分别是均匀化过程的极化电阻和时间常数。将表达式与各类电流表达式 I(t)相乘建立固态扩散过电压表达式 u(t),用于分析电极材料各
16、参数对LIB性能的影响,如式(3)所示:u(t)=I(t)ZVoigt(t)+Zhomo(t)(3)1.2 应用FPA从弛豫电压中识别模型参数应用 FPA将弛豫电压表达式拟合测量值辨识模型参数,由于弛豫电压中包含正负极固态扩散过程和电解液中扩散过程的过电压,需分别建立其过电压表达式。将脉冲电流表达式Ip(t-t0)-t-(t0+Tp)代入式(3)建立固态扩散过电压表达式,如式(4)、(5)所示:uVoigt(t)=k=1nakRDIp(1-e-t-t0bkD)-1-e-t-(t0+Tp)bkD (4)uhomo(t)=RhomoIp(1-e-t-t0homo)-(1-e-t-(t0+Tp)ho
17、mo)(5)式中:Ip为脉冲电流大小;t0为施加脉冲电流的时刻;Tp为脉冲电流持续时间。对电解液中扩散过程 RC 元件阻抗9做拉普拉斯逆变换并与脉冲电流表达式相乘建立过电压表达式,如式(6)所示:uDiff,e(t)=RDiff,eIp(1-e-t-t0Diff,e)-1-e-t-(t0+Tp)Diff,e(6)式中:RDiff,e、Diff,e分别为电解液中扩散过程的极化电阻和时间常数。上述过电压与LIB开路电压uOCV之和为弛豫电压U(t),如式(7)所示:U(t)=uVoigt(t)+uhomo(t)+uDiff,e(t)+uOCV(7)将式(7)拟合弛豫电压辨识模型参数,目标函数设为i
18、=1N|ui-ui2,其中 N 为电压采样点的数目;ui为弛豫电压测量值;ui为模拟值。图1(a)单个颗粒固态扩散阻抗的Voigt模型和(b)电极固态扩散阻抗的ECM8792023.7Vol.47No.7研 究 与 设 计应用FPA找到使目标函数最小的参数值,原理15为:设有n个产生花粉的植物,授粉方式为异花授粉(全局搜索)和自花授粉(局部搜索)。异花授粉公式为:xt+1i=xti+L()(g*-xti)(8)式中:xti为第t次迭代时第i个花粉;g*为已找到的最优解;为控制步长的尺度因子;L()为花粉在不同种类植物间的传输因子,表达式为()sin(2)1h1+,其中 ()为伽马函数;h为步长
19、。自花授粉的公式为:xt+1i=xti+k(xtj-xtk)(9)式中:xtj、xtk为第 t次迭代时不同植物的花粉;k为0,1间均匀分布的随机数。设定总迭代次数,在每次迭代时根据切换概率 P选择异花授粉或自花授粉来搜索最优解13。对比具有全局搜索方式的 FPA 与具有局部搜索方式的TRR14,19对模拟数据的拟合结果,拟合步骤为:设定式(7)中参数生成模拟数据;施加平均值为0.3 mV的随机噪声模拟采样设备误差;设定两组初值,相对误差范围分别为 0.007 11%265%与0.084 6%2 740%;分别使用两种算法拟合。1.3 研究商用磷酸铁锂电池固态扩散阻抗选择标称容量为 6.5 Ah
20、的32700型磷酸铁锂电池为实验电池,正极磷酸铁锂具有平面颗粒,负极石墨具有圆柱形颗粒,根据颗粒形状选择相应 Voigt 模型系数代入式(7)建立弛豫电压表达式。使用CT-4008-5V6A-S1型电池性能测试系统测量弛豫电压,流程为:在待测温度(15、20、25、30和35)下将电池以 0.2 C 倍率放电至待测 SOC(10%、20%、40%、60%、80%、95%),弛豫 24 h;以 0.5 C 倍率脉冲放电 20 s;再次弛豫直至电压趋于稳定并采样。应用 FPA 将弛豫电压表达式拟合测量值得到模型参数,分析参数随电池SOC、温度的变化。2 结果分析和讨论2.1 建立表示圆柱形颗粒固态
21、扩散阻抗的Voigt模型建立圆柱形颗粒的Voigt模型以用于石墨负极材料,运用有限体积法将颗粒离散化为一定数目的体积段,如图 2(a)所示。图中 r0为颗粒半径;h为颗粒高度;Jr为各体积段间锂离子通量。各体积段9的平均半径 rk为 r0(1-2 k-12 n)(k=1,2,.n),其中 n 为体积段数目;平均体积 Vk为 hr20(1-k-1n)2-(1-kn)2;各体积段电容为VkZF-c。假设锂离子在具有不可渗透边界的颗粒中进行3,9径向扩散,菲克扩散定律可表示为:c(r,t)r=-J(r,t)A(r)2 rhD(10)两个体积段间锂离子流量J(r,t)A(r)恒定9,将式(10)积分后
22、得到锂离子浓度表达式:c(r,t)=-J(r,t)A(r)2 Dhln(r)+C(11)各体积段间电阻为锂离子浓度与流量的比值,如式(12)所示:cn+1(t)-cn(t)J(r,t)A(r)=ln(rn-1rn)2 Dh(12)建立表示圆柱形颗粒固态扩散阻抗的链梯模型,如图2(b)所示。在状态空间中,以各电容两端电压UCk(k=1,2,.n)及流经各电阻的电流 IRk(k=1,2,.n)为状态变量,通过变换状态空间方程将链梯模型9降阶为Voigt模型,如图2(c)所示。模型系数akbk(k=1,2,.n)通过变换状态空间方程得到,与n相关9,如表1 所示,为确保准确性n应不小于1009。2.
23、2 对比FPA和TRR拟合结果对比全局搜索的 FPA 与局部搜索的 TRR 对模拟数据的拟合结果,如表 2 所示。初值 1(相对误差范围 0.007 11%265%)下,FPA 拟合结果中各参数相对误差均为 0%,与参考值相同;TRR 结果中除 R-D、Diff,e和 RDiff,e外各参数相对误差均超过 10%。初值 2(相对误差范围 0.084 6%2 740%)下,FPA拟合结果相对误差保持为 0%,不受初值变化影响且与参考值相同;TRR结果中除Diff,e和RDiff,e外各参数相对误差均超过50%,且-homo的相对误差达到 1 600%,与参考值相差较大。这表明FPA的结果更优,相
24、较于TRR更适合于辨识本文所建图2(a)离散化为三个体积段的圆柱形颗粒;表示圆柱形颗粒固态扩散阻抗的(b)链梯模型和(c)Voigt模型表 1 表示圆柱形颗粒固态扩散阻抗的 Voigt 模型系数(n=100)k ak/10-4 bk/10-4 k ak/10-4 bk/10-4 k ak/10-4 bk/10-4 k ak/10-4 bk/10-4 1 1 360.000 681.000 26 2.890 1.480 51 0.719 0.400 76 0.296 0.188 2 406.000 203.000 27 2.680 1.380 52 0.690 0.385 77 0.288 0.
25、184 3 193.000 96.6.000 28 2.490 1.280 53 0.662 0.371 78 0.279 0.180 4 113.000 56.300 29 2.320 1.200 54 0.636 0.358 79 0.271 0.175 5 73.700 36.900 30 2.160 1.120 55 0.611 0.346 80 0.263 0.171 6 51.900 26.000 31 2.020 1.050 56 0.588 0.334 81 0.255 0.168 7 38.600 19.300 32 1.900 0.988 57 0.566 0.323 82
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 锂离子电池 固态 扩散 阻抗 Voigt 建模 参数 辨识
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。