基于多项式拟合和寿命分布的锂离子电池动态可靠性评估方法技术

本发明专利技术公开了一种基于多项式拟合和寿命分布的锂离子电池动态可靠性评估方法，包括步骤一、获取工况折合因子关于虚拟失效阈值的函数关系；步骤二、得到产品的标准退化轨迹；步骤三、获取产品当前时刻的状态偏移量；步骤四、得到状态偏移度；步骤五、得到产品各时刻的动态实时可靠度。本发明专利技术充分利用了同类产品信息和待评估产品的实时监测信息，实现了对锂离子电池进行动态可靠性评估，评估效果显著，与现有的可靠性评估方法相比，本发明专利技术方法显著提高了通用性和精度。

【技术实现步骤摘要】

【技术保护点】
一种基于多项式拟合和寿命分布的锂离子电池动态可靠性评估方法，具体包括以下步骤：步骤一、获取工况折合因子，得到工况折合因子关于虚拟失效阈值的函数关系；当虚拟失效阈值的虚拟失效分布服从指数分布，具体为：若锂离子电池在标准工况下对应虚拟失效阈值Ln的虚拟失效分布服从指数分布，分布参数为μ1，虚拟失效阈值Ln即为当性能退化量达到某一值Ln时假设电池失效，标准工况即为工况一，μ1根据在工况一中运行得到的试验样本数据采用参数估计的方法估计得到，得到工况一下对应虚拟失效阈值Ln的虚拟可靠度函数为：R1(t1)=exp{-t1μ1}其中，t1表示工况一下性能退化量达到Ln所需的时间；锂离子电池在工况二下对应虚拟失效阈值Ln的虚拟失效分布也服从指数分布，分布参数为μ2，μ2根据在工况二中运行得到的试验样本数据采用参数估计的方法估计得到，则工况二下对应虚拟失效阈值Ln的虚拟可靠度函数为：R2(t2)=exp{-t2μ2}其中，t2表示工况二下性能退化量达到Ln所需的时间；由R1(t1)＝R2(t2)可得，在虚拟失效阈值Ln下工况二对应工况一的工况折合因子为：kn=t1t2=μ1μ2通过上述公式，获取工况折合因子与虚拟失效阈值的对应关系，即(Li,ki)，Li表示第i种虚拟失效阈值，ki表示Li的工况折合因子，i＝1,2,…，采用回归分析的方法建立工况二对应工况一的工况折合因子与虚拟失效阈值的函数关系k＝f(L)，其中L是虚拟失效阈值，k为对应的工况折合因子；锂离子电池在运行过程中经历了两个工况，设定锂离子电池当前的运行工况为工况一， 工况切换之前的运行工况为工况二，选定能够表征产品性能退化的性能退化特征量，并对其进行实时监测，得到产品的性能退化数据{(L1,t1),(L2,t2),…，(Lm,tm)}，ti表示产品第i个监测时刻，Li表示产品在ti时刻监测得到的性能退化特征量的值，将每一个Li都看做是其虚拟失效阈值，由于产品前后工况不一致，因此，根据运行时的工况不同，性能退化数据分为两组，工况二时的性能退化数据{(L1,t1),(L2,t2),…,(Ln,tn)}和工况一时的性能退化数据{(Ln+1，tn+1)，(Ln+2，tn+2)，…，(Lm，tm)}，将性能退化数据归一化到工况一中，具体分为以下三步：（1）将工况二中监测得到的性能退化特征量{L1,L2，…，Ln}，输入到k＝f(L)中，得到性能退化特征量{L1，L2,…，Ln}，对应的工况折合因子{k1,k2,…,kn}；（2）利用公式t工况一＝k×t工况二，其中，t工况一表示工况一对应时间，k表示工况折合因子，t工况二工况二对应时间，由工况折合因子{k1,k2,…,kn}和实际运行时间{t1，t2，…，tn}，可以得到折合后相当于在工况一中运行的时间{t′1，t′2，…，t′n}；（3）对于锂离子电池在工况一中的监测得到实际运行时刻，分为两部分，包含了工况切换前在工况二中运行的时间和工况切换后在工况一中运行的时间，因此对时间进行修正，修正公式为：ti′=ti-tn+tn′其中：ti为第i个实际监测时刻，tn工况切换的实际时间，t′n为tn经过工况折合因子修正之后的运行时间，t′i为修正后的监测时刻，i＝n+1,n+2,…，m；最终得到归一化之后的工况一下的性能退化数据{(L1,t′1),(L2,t′2),…，(Lm，t′m)}；步骤二、锂离子电池标准退化轨迹拟合，利用多个全寿命退化样本的性能退化量数据，采用多项式拟合的方法拟合得到产品的标准退化轨迹，其函数表达式为：y=f(t)=Σi=0naiti其中，y为性能退化量，t为性能退化量y对应的监测时刻，n为多项式的阶数，ai为i次项的系数；步骤三、状态偏移量获取，利用步骤二得到的标准退化轨迹和实时监测得到的性能退化数据，获取产品当前时刻的状态偏移量；设当前时刻监测得到的性能退化数据为(y0,t0)由当前监测时刻t0和步骤二得到的标准退化轨迹函数得到标准性能退化量的值y′0，值状态偏移量计算公式为：o=|y0-y0′|步骤四、状态偏移量归一化，选择指数归一化函数，根据同类产品在各时刻的可靠度和状态偏移量，确定归一化函数中的参数c的取值，根据得到的归一化函数，将状态偏移量归一化，得到状态偏移度；以一条样本数据为例，具体步骤如下：（1）获取样本状态偏移量将该实验样本各时刻的状态运行特征值{(y1,t1),(y2,t2),…,(yn,tn)}按照步骤三中的方法得到相应的状态偏...

【技术特征摘要】

【专利技术属性】
技术研发人员：栾家辉，刘正高，吴雷，曾辉，顾长鸿，单添敏，
申请(专利权)人：中国航天标准化研究所，
类型：发明
国别省市：