【技术实现步骤摘要】
一种单面PCB结构稳态热分析方法
本专利技术属于PCB热分析的
，具体涉及一种单面PCB结构稳态热分析方法。
技术介绍
PCB是电子元器件的主要载体，已经成为了电子设备的重要组成部分，其可靠性关系到整个电子设备的可靠运行。而影响PCB及其上元器件可靠工作的主要因素之一就是温度。PCB结构的设计要保证其上元器件的发热能够有效的扩散并与外环境进行热交换，即须保证元器件的温升在安全的范围内，PCB也成为了元器件散热的主要媒介。因此，在考虑实际PCB结构与线路布局以及热边界条件的情况下，对PCB进行热分析可实现辅助优化PCB设计与分析电子系统可靠性的作用。目前市场上应用较为广泛的PCB热仿真软件为Mentor公司的FloTHERM与ANSYS公司的Icepak。在算法原理上前者主要基于与PCB局部覆铜率相关的近似等效热导率，而后者则基于有限元法。而同样基于有限元法的COMSOL与SigrityPowerDC也可在完成PCB及器件建模的基础上进行热仿真分析。而目前国内尚无针对PCB热分析的相关软件。Men...
【技术保护点】
1.一种单面PCB结构稳态热分析方法，其特征在于，包括：/nS1、构建单面PCB绝缘层稳态热平衡拉普拉斯方程，并基于PCB绝缘层的边界条件求解傅里叶级数解析解系数，得到绝缘层温度解析解矩阵方程；/nS2、基于元器件热阻参数计入元器件热效应，修正基板绝缘层边界条件，并将元器件的发热率作为已知热源条件计入解析解的矩阵方程表达式；/nS3、基于有限体积法近似离散金属层热平衡方程，归纳得到金属层热扩散集总矩阵方程，并在其中计入金属层与其表面发热器件间的热传导；/nS4、根据修正后的绝缘层温度解析解矩阵方程和金属层基于数值离散的热扩散矩阵方程，构建单面PCB结构的热分析耦合方程组；/...
【技术特征摘要】
1.一种单面PCB结构稳态热分析方法，其特征在于，包括：
S1、构建单面PCB绝缘层稳态热平衡拉普拉斯方程，并基于PCB绝缘层的边界条件求解傅里叶级数解析解系数，得到绝缘层温度解析解矩阵方程；
S2、基于元器件热阻参数计入元器件热效应，修正基板绝缘层边界条件，并将元器件的发热率作为已知热源条件计入解析解的矩阵方程表达式；
S3、基于有限体积法近似离散金属层热平衡方程，归纳得到金属层热扩散集总矩阵方程，并在其中计入金属层与其表面发热器件间的热传导；
S4、根据修正后的绝缘层温度解析解矩阵方程和金属层基于数值离散的热扩散矩阵方程，构建单面PCB结构的热分析耦合方程组；
S5、基于有限体积法近似离散金属层电流连续性方程，计算归纳得到描述金属层单元电势间线性关系的集总矩阵方程，并求出金属层内的电势分布，从而进一步求出金属层的电流密度分布与焦耳发热分布，并可将焦耳发热分布计入金属层热扩散矩阵方程，与温度分布间开展迭代计算，以求得计入线路焦耳发热前提下的PCB表面的温度分布。
S6、采用COMSOL建模运算结果验证单面PCB结构的热分析方法的计算精度。
2.根据权利要求1所述的单面PCB结构稳态热分析方法，其特征在于，所述步骤S1构建单面PCB绝缘层稳态热平衡拉普拉斯方程，并基于PCB绝缘层的边界条件求解傅里叶级数解析解系数，得到绝缘层温度解析解矩阵方程，包括：
构建单面PCB绝缘层稳态热平衡拉普拉斯方程，并得到PCB绝缘层四周与上、下表面的边界条件：
其中，T为相比于环境温度的温度变化量或在环境温度为0℃时的温度，qiu(x,y)为上表面传入的热流密度分布函数，hu与hd分别为绝缘层上、下表面热传递系数，ki为绝缘层热导率，Lx、Ly、Lin分别为绝缘层在笛卡尔坐标系下x、y方向的长度以及在z方向的厚度；
稳态热平衡方程解析解在傅里叶级数下的表达形式为：
其中，C1，C2，Cn,m，Cγn,m为傅里叶级数系数，βn、μm、γn,m为傅里叶级数中的特征值；
将傅里叶级数解代入T对应的下表面z＝Lin热边界条件，求得解中部分系数与特征值的关系，当hd不为零时，C2及Cγn,m对应的表达式为：
其中，Hu和Hd为解析解推导过程的中间参数；当hd为零时，C1,Cγn,m为：
将级数解代入上表面的热边界条件，得到：
将等式两边同乘cos(βnx)cos(μmy)，并在表面发热区域进行积分,根据三角函数的特性，不同级数相乘后的积分将为零，因此方程左边只有包含cos2(βnx)cos2(μmy)的级数项积分后不为零，而方程右边则对发热区域内的各个发热单元qiu,i进行积分后再求和，经运算后求得系数C1与Cn,m：
其中，dq为将表面热源网格离散化后的网格单元长度，积分域Sqiu,i即代表传入热流密度为qiu,i的表面离散单元区域，δn与δm为解析解推导过程的无量纲中间参数，qiu,1～qiu,N1为将绝缘层表面热源网格离散化后的单元热流密度；若某热源传入的热流密度较为均匀时，则不用对该热源所在区域传递的纵向热流密度进行网格离散；进一步的，可采用向量运算的方式来表达傅里叶级数的求和运算，对于上面的解析解，qiu,1～qiu,N1即构建了热流密度的列向量qiu，则单面结构绝缘层中任一点的解析解可表达为：
求解指定区域的温度分布即为将热阻行向量Ru(x,y,z)根据区域对应分布点的坐标扩展成一个热阻矩阵Ru，其中的每一行即对应求解某点温度的热阻行向量，再将Ru与对应热流密度向量q相乘即可得到指定区域的温度分布向量，即对应q在该指定区域产生的热效应；
由此，得到单面PCB绝缘层温度解析解矩阵方程：
Tu＝RMqM+RIqu
其中，金属层传入绝缘层的热流密度分布向量为qM；而qu为表面发热元器件与其覆盖区域间传递的热流密度分布向量，qu中包括qiu(x,y)经离散后的热流密度分布；RM为qM在绝缘层表面产生的热效应所对应的热阻矩阵；RI为qu作用于绝缘层所对应的热阻矩阵；由于qu中通过金属层表面传递的热流将在金属层对应的热扩散矩阵方程中被计入，则上述解析解矩阵方程中将仅计入表面发热元器件与其直接覆盖的绝缘层表面区域间的传递的热流，因此RI中将包括部分全零列向量。
3.根据权利要求1所述的单面PCB结构稳态热分析方法，其特征在于，所述步骤S2基于元器件热阻参数计入元器件热效应，修正基板绝缘层边界条件，并将元器件的发热率作为已知热源条件计入解析解的矩阵方程表达式，包括：
设表面发热器件覆盖区域的热流密度向量为qh,u：
qh,u＝qu+huTJ,u
其中，TJ,u为表面器件占据区域的温度分布向量；由于被元器件覆盖区域不存在直接与外环境间在hu条件下计算的热交换，但该部分热交换却已在前述步骤S1的解析解推导过程中被计入，因此须对元器件所占区域根据hu多计算的热交换进行补偿；并利用元器件提供的热阻参数推导相关的元器件与PCB表面之间的热传递，具体推导过程如下：
联立上述两式得：
qh,u＝CP,uqJ,u-CT,uTJ,u
其中，表面元器件表面的温度分布为Ttop，其表面与外界的热传递系数为htop，其表面存在的散热器对应的等效热传导也被计入htop；结心平均温度为Tj；表面与结心间的平均热阻为ψjt；结心与器件底部间的平均热阻为ψjc；阻焊层涂层的等效热阻为ψc,u；结心向元器件表面传递的热流密度为qjt；表面元器件对应的发热体积密度向量为PJ,u；由元器件高度参数组成的对角阵为CJ,u；表面元器件对应的发热面密度向量为qJ,u；CP,u为计入器件热阻参数后对应qJ,u的系数；CT,u为计入器件热阻参数后对应TJ,u的系数；
用qh,u替换表面温度分布解析解Tu中的qu：
其中，RICP,uqJ,u项为已知常数向量，而Tu中包括TJ,u向量的全部元素，因此可将矩阵系数通过添零升维的方式进行变换，以实现在矩阵方程中用Tu代替TJ,u，以达到统一向量的目的，等效变换等式如下：
CJTu＝RICT,uTJ,u
其中，CJ为实现TJ,u与Tu相关项等效变换的系数，进而得到：
Tu＝RMqM+RICP,uqJ,u-CJ,uTu。
4.根据权利要求1所述的单面PCB结构稳态热分析方法，其特征在于，所述步骤S3基于有限体积法近似离散金属层热平衡方程，归纳得到金属层热扩散集总矩阵方程，并在其中计入金属层与其表面发热器件间的热传导，包括：
采用结构化正方形网格将金属层进行离散，并基于二阶Taylor级数近似对金属层内热传导的热平衡方程进行离散，得到离散金属层热平衡方程：
其中，km为金属层热导率；qc为金属层的单位体积焦耳发热率分布；Tc为金属单元温度；dc为经结构化离散后的金属单元的单位边长；dm为金属层的厚度；qc,u为Tc对应单元其上...
【专利技术属性】
技术研发人员：张亚斌，郭嘉，董秀成，吴昌东，张彼德，陈琦，
申请(专利权)人：西华大学，四川文化艺术学院，
类型：发明
国别省市：四川;51
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。