一种非平稳工况下IGBT模块寿命预测方法与流程

本发明涉及一种非平稳工况下igbt模块寿命预测方法，它基于器件电热模型和工作任务剖面计算出实际工况下功率模块的结温-时间曲线，提取出相应的结温波动信息；然后，基于功率模块的寿命预测模型，计算不同结温波动对器件造成的累积损伤；最后，结合各结温应力循环次数，对实际工况下igbt模块的mttf进行评估。此方法属于igbt模块寿命评价
技术领域：
。(二)
背景技术：
：由于绝缘栅双极型晶体管(igbt)模块具有电压应力高、电流密度大、导通压降低以及开关损耗低等优点，目前大多数功率变流器选择igbt模块作为功率器件。根据最近工业界的一份调查显示，以igbt模块为代表的功率器件失效占到了变流器失效的40％，因而针对实际工况下功率模块可靠性展开研究有助于提高功率变流器整体的使用可靠性。但由于igbt模块内部结构复杂，其可靠性问题涉及多学科的交叉，相关研究仍处于初级阶段，部分退化机理和模型尚不明确。通常，结温对igbt模块使用可靠性影响最为严重，而功率igbt模块内部芯片的结温往往无法直接获取。目前可通过以下四种方法测量得到：物理接触式测试法、光学非接触测试法、热阻抗模型预测法及热敏感电参数法。其中热敏感电参数法是通过测量器件随结温变化的电气特征参数反推出芯片结温，具有响应快、精度高、可在线测量等优点，但已有提出的不少热敏感电参数测量时需要引入外界激励源，影响功率模块的正常实际工作。此外，随着igbt模块的老化进程各热敏感电参数与结温的对应关系也会发生相应的改变，造成结温测量误差。因此，若能研究出和器件老化因素较小、且能实现非侵入式检测的热敏感电参数，从而建立出实际条件下的结温在线监测方法很有重要意义。随着igbt模块结构工艺的改进，已有的寿命预测模型在窄结温幅值下预测精度较差；同时，igbt模块常常应用于非平稳工况条件下，有必要对非平稳工况下的igbt模块进行寿命预测。为此，本方法以结温预测模型为基础，提出了一种考虑窄结温幅值的寿命预测方法，建立一种非平稳工况下igbt模块寿命预测方法。(三)技术实现要素：1、目的：本发明的目的是：提供一种非平稳工况下igbt模块寿命预测方法，该方法考虑了igbt模块窄结温幅值的情况。与传统的寿命预测评价方法相比，该仿真评价方法准确性更高，适用范围更广泛。2、技术方案：本发明一种非平稳工况下igbt模块寿命预测方法，它包括如下步骤：步骤一：igbt模块的退化机理调研典型键合封装的功率igbt模块的内部结构是由多层具有不同的热膨胀系数的材料组成。由于模块封装各层的材料不同，在温度循环的过程中材料之间热膨胀系数的不匹配会在不同层之间的界面上产生热机械应力，这种应力会造成igbt模块中键合线和焊料层的疲劳损伤。所以igbt的主要失效机理是键合线的裂纹和剥离，以及焊料层的空洞。步骤二：igbt模块电学模型建立igbt实际运行时，可在不借助额外的传感器或内置测量电路的前提下测量得到器件集电极-发射极两端电压vce和集电极电流ic。因此，建立igbt模块在导通情形下的饱和压降vcesat的数值模型，并针对该机理模型结合前述讨论的igbt模块退化机理，分析确定考虑退化特性的结温在线监测参数。该参数可用于在线监测模块内部igbt芯片正常工作时的结温。步骤三：基于simulink的igbt电学模型参数解耦在电学模型中，集电极电流ic、结温tj参数之间耦合关系复杂，若要考虑退化特性影响在线监测器件结温，还需对电学模型中的集电极电流ic、结温tj进行参数解耦。封装部分压降只与集电极电流ic成正比关系，只需对通态压降部分进行解耦。另一方面，由于难以通过纯数学变换方法对von进行参数解耦，因此在matlab/simulink中建立器件电学模型，考虑温度效应影响，仿真通态压降von同结温、集电极电流的关系，实现参数解耦。步骤四：考虑退化特性的结温在线监测搭建结温测试平台，为避免键合丝与铝金属层结构退化的影响，采用集电极开启电压vn作为热敏感参数用于在线监测功率模块内部igbt芯片的结温测量：vn＝kv·tj+bv集电极开启电压vn和结温tj呈负温度系数线性关系，能在功率模块实际工作电流较大的条件下对器件结温实现非侵入式的测量，能有效应用于复杂工况下结温的在线监测。步骤五：窄结温幅值寿命预测模型建立采用卡尔曼滤波的方法，结合建立的考虑退化特性的结温监测方法，得到功率igbt模块全开关周期的结温预测模型，从而获得功率器件工作时的结温-时间历程曲线，为后续功率循环试验的升温和降温时间提供依据，为实际工况下功率igbt模块的寿命预测提供参考。步骤六：功率循环试验平台搭建功率循环试验电路主要由功率器件大电流老化部分和小电流结温测量部分组成，根据卡尔曼滤波方法确定升温与降温时间，并确定功率循环试验终止条件。其中大电流老化部分主要是对器件通以大电流，使功率器件结温在短时间内快速上升，是功率循环试验平台电路的主要组成部分；小电流结温测量部分是在器件通以小电流10ma～100ma，可在大电流关断条件下测量器件饱和压降vcesat以间接获得器件结温。步骤七：窄结温幅值寿命预测模型建立根据多组不同应力条件的功率循环试验，基于试验数据的统计分析结果，利用basquin方程和arrhenius方程建立窄结温幅值情形下的功率igbt模块寿命预测模型，为后续实际非平稳工况下功率igbt模块剩余寿命的预测提供输入。步骤八：非平稳工况igbt模块寿命预测模型在寿命预测理论基本实现方法的基础上，使用雨流算法提取器件结温随机波动载荷信息。考虑功率igbt模块不同工况区间的分布特点及器件结温幅值对器件损伤影响，提出非平稳工况下功率igbt模块的寿命预测方法。(四)附图说明：图1是本发明的实施步骤流程示意图图2是simulink中建立的器件电学模型图3是结温tj保持不变，通态压降von与集电极电流ic的关系图4是集电极电流ic保持不变，通态压降von与结温tj的关系图5是结温测量试验平台电路原理图图6是结温测量试验平台图7是集电极开启电压vn与结温tj对应关系图8是卡尔曼滤波结温预测流程图图9是功率模块基于卡尔曼滤波结温预测方法整体思路图10是卡尔曼滤波的结温预测结果示例图11是功率循环试验平台电路示意图图12是搭建的功率循环试验平台示意图图13是窄结温幅值寿命预测模型拟合结果示意图图14是非平稳工况下功率igbt模块寿命评估流程图图15是基于雨流计数法的随机结温波动信息提取示意图(五)具体实施方式：下面将结合附图和非平稳工况下igbt结温预测方法，对本发明作进一步的详细说明。本发明一种典型cmos器件极限低温特性仿真方法，具体步骤如下：步骤一：igbt模块的退化机理调研igbt模块的退化机理可分为与芯片相关的退化机理(如时间相关介质电击穿、电迁移和铝金属化重构等)和与模块封装有关的退化机理。与芯片有关的退化通常和半导体物理相关，由硅材料内部缺陷引发并最终导致器件退化失效。这一失效虽不同于模块封装引起的退化失效，但和模块封装结构的变化可能有关。相比与芯片有关的退化机理，在众多应用领域科研人员更关注与模块封装有关的退化机理，因为模块封装的退化或失效会引起器件工作结温的升高，进一步导致器件安全工作区和性能参数的退化。由于模块封装各层的材料不同，在温度循环的过程中材料之间热膨胀系数的不匹配会在不同层之间的界面上产生热机械应力，这种应力会造成igbt模块中键合线和焊料层的疲劳损伤。当键合线的疲劳累积损伤到达一定程度时，就会导致铝键合线产生裂纹和剥离；而焊层的累积损伤也会造成硅芯片与直接覆铜基板、dcb与底板之间的焊料层产生空洞。这两种失效机理是功率igbt模块中最常见与模块封装相关的失效类型。步骤二：igbt模块电学模型建立igbt实际运行时，可在不借助额外的传感器或内置测量电路的前提下测量得到器件集电极-发射极两端电压vce和集电极电流ic。因此，建立igbt模块在导通情形下的饱和压降vcesat的数值模型，并针对该机理模型结合前述讨论的igbt模块退化机理，分析确定考虑退化特性的结温在线监测参数。该参数可用于在线监测模块内部igbt芯片正常工作时的结温。一般的，测量得到的器件饱和压降主要包括器件芯片上的通态压降von和集电极-发射极通道上封装部分的压降：vcesat＝von(tj,ic)+rcon·ic式中，von为器件芯片上的通态压降，是关于器件结温tj和集电极电流ic的一个复杂函数；rcon为铝金属层、键合线等电气互连结构形成的接触电阻。随着igbt模块键合线以及铝金属化层不断退化，rcon也会不断增大，因此在恒定ic下测量器件饱和压降表征器件结温时，键合线与铝金属化层的退化会对测量结果的准确度造成影响。步骤三：基于simulink的igbt电学模型参数解耦在电学模型中，集电极电流ic、结温tj参数之间耦合关系复杂，若要考虑退化特性影响在线监测器件结温，还需对电学模型中的集电极电流ic、结温tj进行参数解耦。封装部分压降只与集电极电流ic成正比关系，只需对通态压降部分进行解耦。另一方面，由于难以通过纯数学变换方法对von进行参数解耦，因此在matlab/simulink中建立器件电学模型，考虑温度效应影响，仿真通态压降von同结温、集电极电流的关系，实现参数解耦。选取典型键合封装的三菱cm150dy-12nf型(150a/600v)igbt模块为研究对象，电学模型仿真中所需的相关器件参数通过经验估计、器件手册外推得到，相应取值如表1所示：表1电学模型中相关参数取值建立的simulink器件电学模型如图2所示。在matlab中利用sim函数调用simulink进行数值仿真。仿真主要通过两个步骤实现参数解耦：①结温tj保持不变，仿真通态压降von与集电极电流ic的关系；②集电极电流ic保持不变，仿真通态压降von与结温tj的关系。根据igbt实际应用情况与本文所用igbt型号的器件手册规定，仿真中结温tj取值为300k到450k之间，集电极电流ic取值为10ma到200a之间，仿真结果如图3和图4所示。从图3可以看出，在上述四个不同的恒定结温tj取值下，当集电极电流ic取值较小时，通态压降von和集电极电流ic的非线性程度比较严重；而当ic超过一定值(5a)时，通态压降von和集电极电流ic呈现为较好的线性关系。因此，器件通态压降von可与结温tj和集电极电流ic作如下解耦：von(tj,ic)＝f(tj)·ic+vn(tj)式中，f(tj)和vn(tj)都是和ic无关、和tj相关的函数。当结温tj保持为恒定值时，通态压降von和集电极电流ic可近似认为呈线性关系。从图4可以看出，在六个不同的集电极电流ic取值下，通态压降von与结温tj均呈现为较好的线性关系。无论电流ic如何取值，f(tj)·ic+vn(tj)均为关于tj的线性函数，由此可知f(tj)和vn(tj)都是与集电极电流ic无关、关于结温tj的线性函数。即f(tj)＝kf·tj+bfvn(tj)＝kv·tj+bv式中，kf、bf、kv、bv均为常数，具体数值可通过实验得到。vcesat＝von(tj,ic)+rcon·ic＝f(tj)·ic+vn(tj)+rcon·ic＝(kf·tj+bf+rcon)·ic+kv·tj+bv步骤四：考虑退化特性的结温在线监测为避免键合线与铝金属层等结构退化，造成器件封装电阻rcon的增大对结温测量结果准确性的影响，以集电极开启电压vn作为热敏感参数用于在线监测功率模块内部igbt芯片的结温测量：vn＝kv·tj+bv集电极开启电压vn和结温tj呈负温度系数线性关系，且不受功率模块键合线与铝金属层疲劳退化影响，线性度好，能在功率模块实际工作电流较大的条件下对器件结温实现非侵入式的测量，能有效应用于复杂工况下结温的在线监测。结温测量试验平台电路原理图如图5所示。主电路由15v直流电源、恒流源和cm150dy-12nf型号的igbt模块构成，待测器件放于恒温箱中。采用高速数据采集卡，记录下待测器件开通瞬间的集射极压降变化情况，存储于pc机中。试验人员可从pc机中读取得到器件开通瞬间的饱和压降vcesat。搭建的结温测量试验平台如图6所示。进行结温测量试验时，将待测igbt模块置于恒温箱中，恒温20分钟后(使模块内部温度稳定在所设温度)施加大小为15v的正向偏压，通以不同集电极电流ic测量该温度下器件导通压降vcesat。参照本发明所用igbt器件手册规定及igbt实际应用情况，试验时电流取值为10ma、100ma、1a、10a、20a、30a、40a、50a、60a、70a、100a、120a、150a，饱和压降测量间隔为2分钟，以避免引入器件自热效应造成测量误差；恒温箱试验时温度设置范围为30℃～130℃，步长δt为10℃，测量igbt芯片在11个不同结温值下的饱和压降。试验测得的饱和压降如表2所示：表2饱和压降vcesat测量结果tj/℃3040506070809010011012013010ma0.3210.3030.2850.2640.2480.2310.2120.1940.1760.1600.144100ma0.4310.4130.3920.3720.3500.3320.3130.2910.2680.2490.2291a0.5690.5460.5230.5030.4790.4550.4330.4120.3900.3700.34910a0.8740.8460.8170.7900.7610.7330.7050.6750.6480.6210.59420a1.0020.9760.9490.9210.8950.8700.8440.8190.7910.7640.73830a1.0881.0641.0411.0180.9950.9710.9460.9210.9010.8790.85040a1.1651.1421.1191.0991.0751.0531.0311.0100.9900.9690.94650a1.2311.2111.1931.1711.1531.1321.1121.0931.0771.0601.04160a1.3011.2851.2691.2511.2331.2161.2011.1851.1671.1521.13570a1.3761.3631.3481.3351.3191.3051.2901.2771.2581.2421.227100a1.6481.6371.6281.6171.6081.5991.5921.5811.5731.5611.554120a1.7891.7851.7821.7781.7771.7741.7701.7681.7641.7601.758150a2.0382.0412.0422.0452.0482.0512.0532.0552.0562.0592.062从上述结果可以看出，当集电极电流ic超过一定范围时，器件的饱和压降vcesat和结温tj、集电极电流ic呈良好的线性关系。在电流较小(如10ma和100ma)且保持恒定时，vcesat和tj呈现非常好的线性关系，且灵敏度较高；但当器件工作不断增大时，vcesat和tj的灵敏度不断下降，尤其工作电流为120a～150a左右时，器件饱和压降vcesat与结温tj都近似成水平关系，灵敏度较低，在此种情形下测量器件饱和压降vcesat得到的器件结温值和实际结果可能存在较大的误差，这是大电流测量器件饱和压降方法的一个盲区。结合表2中饱和压降测量结果，依据线性回归关系得到器件在不同结温对应的集电极开启电压vn，如表3所示，并作出线性拟合结果如图7所示。表3不同结温的集电极开启电压vntj/℃304050607080vn/v0.82710.79590.76440.73390.70330.6730tj/℃90100110120130vn/v0.64200.61000.58200.55240.5210从图7可以看出，试验得到的集电极开启电压vn和结温tj呈良好的负温度系数线性关系。步骤五：基于卡尔曼滤波的结温预测模型为使用卡尔曼滤波对igbt模块进行全开关周期的结温预测，需获得模块热传导过程中的状态空间描述。热传导过程中的状态空间表达式可表达为如下形式：状态方程输出方程tj(t)＝cx(t)+du(t)其中x(t)代表rc网络每个网格的温度差，且u(t)＝[pdtc]t。状态空间表达式中的状态方程描述了rc网格温度差的变化趋势，而输出方程表示结温tj为fosterrc热网络各网格温度差与壳温之和。矩阵an×n为系统矩阵，bn×2为输入矩阵，c1×n为输出矩阵，d1×2为前馈矩阵。卡尔曼滤波通常用于线性离散情况下的系统状态预测，而上述建立的状态空间表达式是一个连续状态空间函数，因此需要利用向后欧拉法将此连续状态空间函数转换为离散的状态空间函数。取时间步长为ts，即假设k时刻对应时间为kts，并假设系统状态驱动和结温测量时均存在相互独立的白噪声，则得到的离散状态空间函数为：x[k]＝fx[k-1]+gu[k-1]+w[k-1]tj[k]＝cx[k]+du[k]+v[k]其中，c＝[111…d＝[01]x[k]指t＝kts时刻当前系统的状态变量，其余与此类似。考虑到实际复杂工况下功率模块焊料层退化造成rc热阻网络的热阻、热容发生变化，而卡尔曼滤波算法的更新特性使得其能有效适用于实际复杂工况下热流通道退化时功率模块的结温预测。基于卡尔曼滤波的结温预测方法可分为预测和校正两大过程。预测过程主要是利用时间更新方程，即建立的热传导状态空间模型计算当前结温的先验估计；校正过程负责反馈，利用观测更新方程在结温先验估计和热敏感电参数法测量得到的器件结温基础上得到了改进后的结温的后验估计。卡尔曼滤波结温预测的具体表达式如下：预测过程校正过程预测过程和校正过程的具体实施顺序如图8所示。可以看出，若当前时刻存在通过热敏感电参数得到的结温在线监测结果(即tj[k]存在)，则进行结温的预测与校正；反之，若当前时刻不存在结温的在线监测结果(即tj[k]不存在)，则只进行结温的预测。tj[k]不存在时可令c＝0，则此时系统状态估计和估计误差协方差p[k]的先验估计结果和后验估计结果保持一致。以cm150dy-12nf型igbt模块作为研究对象，开展结温预测方法验证。参考其器件手册，基于瞬态热阻抗曲线通过曲线拟合的方法得到模块fosterrc热网络中热阻、热容的值，当rc网格单元为3个时，模型已达到很好的拟合效果，对应的参数拟合结果如表4所示。表4fosterrc热网络参数拟合结果igbt模块的电学模型描述了器件饱和压降同集电极电流、结温的关系，由此可进一步根据器件结温、集电极电流得到当前功率模块的耗散功率，用于卡尔曼滤波结温预测的输入。根据2.3.2节饱和压降机理模型参数解耦结果，有器件饱和压降：vcesat(tj.ic)＝a1×ic×tj+a2×ic+a3×tj+a4其中，a1，a2，a3，a4均为常数。由此拟合结温测量试验得到的器件饱和压降数据，得到的拟合精度r2＝0.9987，且a1＝2.199×10-5,a2＝7.431×10-3a3＝-3.305×10-3,a4＝0.916在此基础上，结合结温在线监测的功率模块基于卡尔曼滤波的结温预测方法的整体思路如图9所示。由测量采集电路得到功率模块当前时刻的集电极电流与饱和压降，根据线性关系反推出集电极开启电压，获得当前时刻的结温在线监测结果，用于卡尔曼滤波的校正；卡尔曼滤波的预测结果即结温tj起到了反馈作用，由结温tj和集电极电流根据电学模型计算出当前器件功耗，结合器件壳温用于卡尔曼滤波的预测步骤。值得指出的是，器件功耗也可由测量采集电路采集的集电极电流和饱和压降计算得到，后者允许在功率模块退化情形下更为精确地进行器件结温预测。同时，假定在测量功率模块饱和压降和壳温时均存在相应的高斯噪声，造成结温在线监测结果误差。基于卡尔曼滤波预测功率模块结温时，均在器件集电极电流不少于10a的情形下得到，这也符合功率模块实际工作在大电流条件下的特点。图10为集电极电流为100a时，功率igbt模块工作周期为2s且占空比为0.5时，取时间步长为0.1ms，基于卡尔曼滤波的结温预测结果。可以看出，只有在功率igbt模块正向导通的情形下，才可进行功率模块结温的在线监测，从而对卡尔曼滤波的结温预测结果进行校正；当功率igbt模块反向阻断时，则可基于卡尔曼滤波只进行结温预测。步骤六：功率循环试验平台搭建功率循环试验平台的电路示意图如图11所示，该电路主要由功率器件大电流老化部分和小电流结温测量部分组成。其中大电流老化部分主要是对器件通以大电流，使功率器件结温在短时间内快速上升，是功率循环试验平台电路的主要组成部分；小电流结温测量部分是在器件通以小电流10ma～100ma，可在大电流关断条件下测量器件饱和压降vcesat以间接获得器件结温。大电流老化电路部分主要包括试验igbt模块(三菱cm150dy-12nf，150a/600v)、额定功率为4kw的安捷伦n5764a可程控直流功率电源、dh1718g-4型直流三路稳压稳流电源、栅极驱动电路、风冷散热系统、tp9000多路温度记录仪和数据采集卡及pc主机。其中，功率循环试验过程中由n5764a可程控直流功率电源向试验igbt模块通以恒定大电流导通，该功率电源最大可输出76a电流。试验中将两个功率电源并联以输出140a电流，并通过控制主电源的周期开关实现功率模块结温的交替变化。功率igbt模块栅极驱动电路主要实现控制电路与被驱动igbt栅极的电隔离。本文用dh1718g-4型直流三路稳压稳流电源在功率igbt模块栅极、发射极两端施加大小为15v的恒定电压，并选用富士公司的光耦驱动模块vla517-01rzip-13用于功率模块的栅极驱动。风冷散热系统利用功率igbt模块和散热器之间涂以导热硅脂进行传热，并且为了提高功率模块的降温速率以降低功率循环试验所需时间，在模块两侧分别放置轴流风机实现空气的平稳对流。采用拓普锐tp9000多通道温度记录仪，在功率模块基板底部放置热电偶监测采集功率循环试验过程中器件壳温变化情况。数据采集卡并与pc主机相连，在功率循环试验过程中定时监测igbt模块的饱和压降，一方面用于评估功率模块的老化程度，另一方面在功率模块大电流关断只通以小电流时获得功率模块结温数据，用于确定功率模块升温、降温时间。对于小电流结温测量电路，根据前面建立的饱和压降解耦模型可以看出，在电流ic较小的情形下，当功率模块键合线疲劳或芯片铝金属层重构造成接触电阻rcon增大时，器件饱和压降vcesat的退化影响较小，因而可近似认为小电流情形下饱和压降vcesat和结温tj的对应关系几乎不受模块老化影响发生变化。采用dh1718g-4型直流三路稳压稳流电源的另一路对功率igbt模块集电极、发射极两端通以恒定100ma小电流，从而可在大电流关断条件下监测器件饱和压降vcesat以间接获得器件结温。根据前面结温测量试验结果，使用matlab拟合函数线性拟合100ma时器件饱和压降vcesat和结温tj的关系，结果如下式所示。tj＝-491.6×vcesat+242.7最终，建立的功率循环试验平台如图12所示。功率模块的升温、降温时间利用卡尔曼滤波法进行粗略的估算，确定模块升降温时间ton和toff的值。设定功率循环试验的最大循环次数(考虑功率循环试验所需时间，这里设定为250000次)，并在arduinouno开发板中编程设定相应功率模块的升降温时间。在功率循环试验过程中功率模块不断加热冷却，监测器件的饱和压降vcesat和壳温，根据大电流条件下监测到的饱和压降判定器件是否发生达到失效阈值(饱和压降变化超过初始值15％)，若达到失效阈值则停止功率循环试验。同时，若器件功率循环次数已超过250000次而未发生失效，也停止功率循环试验。进行9组不同条件的窄结温幅值功率循环试验，有两只igbt管子进行同一试验条件进行试验。同时，为了和窄结温幅值功率循环试验结果进行对比，本文还进行了一组宽结温幅值功率循环试验。各组试验条件及升温、降温时间详细如表5所示。表5功率循环试验条件及升温、降温时间步骤七：窄结温幅值寿命预测模型建立该功率模块cm150dy-12nf在不同应力条件下失效前功率循环次数如表6所示，其中除器件j1外其余器件记录的功率循环次数为饱和压降vcesat相对初始值上升15％的功率循环次数。表6不同应力功率循环试验的失效前功率循环次数功率igbt模块在承受热机械应力冲击时，焊料层及键合线等受到的总应变δεtot可分为如下三类：与时间无关的弹性应变δεe、塑性应变δεp以及与时间相关的蠕性应变δεc，即δεtot＝δεe+δεp+δεc一般而言，当功率igbt模块经受功率循环试验的结温波动幅值较大时(低周疲劳)，材料遭受的应力超过其屈服强度，此时应变主要由塑形应变构成；反之，当功率igbt模块经受功率循环试验的结温波动幅值较小时(高周疲劳)，材料遭受的应力超过其屈服强度，此时应变主要由弹形应变构成。在通常情形下，由于功率igbt模块工作在高频条件下经历较为快速的温度循环，其蠕性应变基本可忽略不计。故而在功率igbt经历温度循环结温幅值较宽时，其寿命主要由塑性应变决定；在功率igbt经历温度循环结温幅值较窄时，其寿命主要由弹性应变决定。根据basquin方程，基于弹性应变的老化模型可将模块寿命和弹性应变用指数关系联系起来，其表达式为：式中，δεe代表器件在一个功率循环中的弹性应变范围；e代表模块的弹性模量；σf和b均为与材料相关的疲劳延性系数和疲劳延性指数。在弹性应变情形下，器件在一个功率循环中弹性应变范围δεe和结温幅值δtj成正比，上式可变为如下形式：nf＝a·(δtj)α参考arrhenius模型，考虑功率循环过程中结温平均值tm对功率循环寿命的影响，有：式中，nf为器件失效前功率循环次数；参数a1、α、a2均为常数，可通过表6的功率循环次数结果拟合得到。借助matlab，参数a1、α、a2的拟合结果分别为4521、-0.432和1758，拟合结果如图12所示。该模型适用于δtj＜80k的窄结温幅值情形，同时模型的相关参数仅限于当前型号的igbt模块。对于其它型号功率igbt模块，参数a1、α、a2可通过该型号器件的窄结温幅值功率循环试验结果拟合得到。步骤八：非平稳工况igbt模块寿命预测模型当功率igbt模块运用于风电变流器，或是机车牵引、电动汽车驱动变流器时，其往往工作于输入功率大范围波动造成结温随机波动的非平稳工况下。针对非平稳工况下igbt模块进行寿命预测时，需要计及结温幅值和工况类型的影响。同时，考虑到实际非平稳工况下难以准确确定功率变流器的工作任务剖面，本方法提出热敏感电参数法监测器件结温，并采用卡尔曼滤波方法获取器件的结温剖面，最终建立的非平稳工况下igbt模块寿命评估流程如图14所示。其评估步骤如下：1)分析功率igbt模块非平稳工况特点，划分各工况区间并确定各工况区间在实际工况下的概率分布。记工况区间序号为i(i＝1,2,…,n)，对应工况区间在非平稳工况所占比例为fi，显然有：2)针对各工况区间，根据热敏感电参数法在线监测器件结温，并基于卡尔曼滤波算法获取功率模块在该工况区间的结温波动剖面。3)采用雨流计数法提取各工况区间结温波动剖面相应的结温幅值、结温平均值和对应的应力循环次数。雨流计数法的步骤如下：a.将基于器件结温在线监测得到的结温-时间历程曲线沿着顺时针方向旋转90°，如图15所示，此时横坐标轴代表器件结温，纵坐标轴代表时间。b.每个峰值(或谷值)为雨流计数法计算的起始点，规定雨流从每个起始点内侧沿着斜坡往下流，进而根据雨点轨迹(图15虚线所示)提取结温波动信息：当雨流从某一峰值(或谷值)内侧开始流动，在流到峰值(或谷值)处竖直落下，直到流到比起始点更大的峰值(或更小的谷值)时停止流动；或者雨滴流动过程中遇到上层斜面留下的雨滴时便停止流动。当雨滴停止流动时，其轨迹就会形成一个完整的闭合曲线，即一个完整的结温波动循环。c.对提取到的某一特定结温波动循环，记其起点值为to，终点值为ts，对应的结温幅值δtj和结温平均值tm可由下式确定：δtj＝|to-ts|d.对应结温波动循环(δtj,tm)的应力循环次数n(δtj,tm)在原先统计结果上加1。e.重复步骤b～d步，直至遍历整个结温-时间历程曲线。4)基于线性疲劳累积损伤理论，计算在工况区间i下功率模块的mttfi：式中，n(δtj,tm)为器件在工况区间i时间ti内应力循环次数统计值；nf(δtj,tm)为根据器件寿命预测模型得到的在结温幅值为δtj、结温平均值为tm的寿命预测结果。在进行寿命预测时，应计及结温幅值对器件功率循环能力的影响，窄结温幅值情形(δtj≤80℃)下基于本方法前述建立的窄结温幅值寿命预测模型进行预测，较宽结温幅值情形(δtj＞80℃)下可基于lesit模型进行预测。5)根据非平稳工况下各工况区间的概率分布，计算实际非平稳工况下功率igbt模块的平均失效前时间mttf：非平稳工况下功率igbt模块的剩余使用寿命也可基于该寿命评估模型计算得到。当前第1页12