1.一种基于降阶方法的IGBT瞬态温度快速计算方法,其特征在于:所述计算方法依次包括以下步骤:
S1、建立IGBT电热耦合强形式及其耦合机制;
S2、建立IGBT电热耦合强形式及其耦合机制的有限维空间有限元形式;
S3、对IGBT电热耦合强形式及其耦合机制的有限维空间有限元形式做离散化处理,将其由偏微分方程转化为矩阵方程;
S4、采集IGBT温度样本矩阵X2,再采用奇异值分解法对温度样本矩阵X2进行降阶处理,并根据截断误差选取特征值数量用于描述整个温度样本矩阵X2;
S5、将降阶处理后的温度样本矩阵X2代入由步骤S3得到的矩阵方程中,对矩阵方程进行整体降阶处理,最后求解经整体降阶处理后的矩阵方程,即可实现IGBT瞬态温度的快速计算。
2.根据权利要求1所述的一种基于降阶方法的IGBT瞬态温度快速计算方法,其特征在于:步骤S1中,所述IGBT电热耦合强形式及其耦合机制为:γ(T)=ρ0((1+α(T‑Tref))上式中, 表示在 这个有界区域内的向量, 为梯度, 为 处的电势,ρ为材料密度,Cp为定压比热,k为导热系数,t为时刻, 表示在 处的温度,Q为热量,γ为电导率, 分别表示 为x、y、z的偏导,γ(T)表示温度为T时的电导率,ρ0为参考密度,α为电阻温度系数,Tref为参考温度。
3.根据权利要求2所述的一种基于降阶方法的IGBT瞬态温度快速计算方法,其特征在于:所述步骤S2依次按照以下步骤进行:
S21、基于IGBT电热耦合强形式及其耦合机制求 推导出IGBT电热耦合强形式及其耦合机制的无限维空间有限元形式,所述无限维空间有限元形式为:(Tt,v2)+a(T,,v2)=(Q,v2)(Tt,v2)=∫ΩρCpT·v2dxdydz1
上式中,(Q,v2)表示Q和v2内积,H为希尔伯特空间,H (Ω)上标表示原函数及其1阶导数1
平方可积, 中 下标表示v1在H (Ω)边界为0, 中 下标表示v2在1
H (Ω)边界为0,表示电势,v1为电势测试函数,v2为温度测试函数,Q表示热量,k表示温度方程的系数, 表示温度变化量, 表示电势测试函数变化量, 表示温度测试函数变化量,Tt为温度T对时刻t求导;
S22、引入时间T1, 基于IGBT电热耦合强形式及其耦合机制的无限维空间有限元形式求 T1h∈U2h,推导出其有限维空间有限元形式,所述有限维空间有限元形式为:(Tt,v2h)+a(T,v2h)=(Q,v2h)上式中, 表示在U1h这个时空空间下的电势,v1h表示U2h这个时空空间内的电势测试函数,v2h表示U2h这个时空空间内的温度测试函数, 为有限元子空间,表示U2h为由0到T1在U1h空间的时空空间, T1h表示在U2h这个时空空间的温度,φj表示有限元节点值,Ns为有限元总节点数,Tj为j这个节点的温度, 表示 与t内积, 为节点的电势, 表示v1h在U2h的边界为0,表示v2h在U2h的边界为0。
4.根据权利要求3所述的一种基于降阶方法的IGBT瞬态温度快速计算方法,其特征在于:所述步骤S3依次按照以下步骤进行:
S31、选择测试函数V1h,V2h=φi(i=1,2,…,Ns),对由步骤S2得到的IGBT电热耦合方程及其耦合机制的有限维空间有限元形式进行数值积分,经数值积分得到以下方程:上式中,c1、c2分别为方程系数,Tj、分别代表j这个有限元节点的温度和电势, 为测试函数φi的梯度, 为测试函数φj的梯度,Tj(t)表示在t时刻的Tj;
S32、根据经数值积分得到的方程,定义电势刚度矩阵A1、温度刚度矩阵A2、关于温度的质量矩阵M、热源载荷向量 其中,所述电势刚度矩阵A1为:所述温度刚度矩阵A2为:
所述关于温度的质量矩阵M为:
所述热源载荷向量 为:
S33、先将电势刚度矩阵A1、温度刚度矩阵A2、关于温度的质量矩阵M、热源载荷向量 代入到标量电位的矩阵方程和温度的常微分方程中,然后采用向后欧拉法对温度的常微分方程进行离散化处理得到方程一,整理方程一得到方程二,将方程二左右同乘Δt得到矩阵方程,所述标量电位的矩阵方程为:A1X1=0
上式中,
所述温度的常微分方程为:
上式中,X2(t)表示t时刻的温度向量,X2′(t)表示为t时刻的温度向量的一阶导数,所述方程一具体为:上式中,Δt表示m到m+1时刻的时间差, 表示m+1时刻的温度向量, 表示m时刻的温度向量,pm‑1为终止时刻,A(tm+1)为在m+1时刻的温度刚度矩阵, 为在m+1时刻的热源载荷向量;
所述方程二为:
所述矩阵方程为:
5.根据权利要求4所述的一种基于降阶方法的IGBT瞬态温度快速计算方法,其特征在于:所述步骤S4依次按照以下步骤进行:
S41、采集IGBT温度样本矩阵X2,所述温度样本矩阵X2为:上式中, 表示第NS个节点第Nt时刻的温度;
S42、将温度样本矩阵X2按T=Uα形式表达,并采用奇异值分解法选取温度样本矩阵X2的前Nd个特征值代表降阶后的温度样本矩阵X2,U为上一节得到的特征向量形成的标准正交基,U具体为:上式中,uNsNd表示为Ns行Nd列的特征向量;
α为降阶后的温度向量,α具体为:
上式中, 为由Ns行降低为Nd行的温度向量;
所述降阶后的温度样本矩阵X2具体为:
6.根据权利要求5所述的一种基于降阶方法的IGBT瞬态温度快速计算方法,其特征在于:所述步骤S5依次按照以下步骤进行:
S51、将降阶后温度样本矩阵X2作为 代入由步骤S3得到的矩阵方程中得到方程三,所述方程三为:T T
S52、将方程三左乘U ,U为U的转置,使矩阵方程整体阶数下降到Nd,得到方程四,所述方程四为:S53、根据方程四求解降阶后温度向量α,再根据T=Uα计算得到瞬态温度。