用ansys软件包模拟功率半导体器件用散热器稳态热阻.docVIP

用ANSYS软件包模拟功率半导体器件用散热器稳态热阻 1 引 言? ? 功率半导体器件是高发热的电器产品，功率半导体器件的结温高低直接关系到器件的工作状态和承受浪涌电流的能力[1]。功率半导体器件的结温一方面取决于器件损耗的大小，另一方面取决于这些损耗被带走的速度，即取决于散热冷却条件。功率半导体器件用的散热器的目的是控制半导体器件的结温Tj ，使其低于半导体器件的最大结温Tjmax ，从而提高半导体器件的稳定性和可靠性。因此，深入而准确地了解功率半导体器件用散热器的温度场分布，对于功率半导体器件的应用，无疑是一个十分重要的基本问题。温度场的仿真计算及分析，是世界各种优化设计其产品的一种重要手段，大到一艘潜艇，一架飞机，小到一台PC机机箱，一块冷却板均应用这一技术进行优化设计。在机车上，因此体积，重量受到严格限制，加上其电气发热密度大，因此通风冷却是机车关键设计技术之一。在我国机车通风冷却设计技术中，还没有或很少利用计算机进行仿真计算，优化，一般还是根据经验进行，设计较粗糙，不能满足生产实验的需要，必须对功率半导体器件用的散热器进行仿真计算，达到优化设计的目的。而ANSYS多物理场软件包恰恰可以对温度场，流场等方面内容进行仿真。因此，本文介绍用ANSYS多物理场软件包对功率半导体器件用散热器进行温度场仿真计算，且计算结果与测试结果吻合得很好。? ? 2 理论分析? ? 根据传热学模拟理论，不论是什么样的散热器模型都必须按照通一的原则加以分析和研究，即求解导热微分方程。我们知道，在传热学中，具有内热源和瞬态温度场的三维导热在直角坐标系中的微分方程式为[2]： t / τ=λ/ ρc （ 2t/ x2+ 2t/ y2+ 2t/ z2+θ/λ） ………… (1)? ? 式中：t：为物体的瞬态温度? ? τ：为时间? ? ρ：密度? ? c：为比热容? ? λ：为导热率? ? θ：为单位体积内热源的生成热? ? 求解上述导热问题，实质上归结为对导热微分方程式的求解。为了使微分方程得到唯一解，还必须给出边界条件和初始条件。? ? (Ⅰ) 规定了边界上的温度值，称为第一类边界条件，用公式描述为：? ? Tw = f1 ( x ，y ，z ，τ) ……… (2)? ? (Ⅱ) 规定了边界上的电流密度值，称为第二类边界条件，用公式描述为：? ? —λ( t/ n)w = f2 ( x ，y ，z ，τ ) ……… (3)? ? 式中：n为表面的法线方向。? ? (Ⅲ) 规定了边界上物体与周围流体间的换热系数a及周围流体的温度tf ，称为第三类边界条件。用公式描述为 ：? ? —λ( t/ n)w = a( tw—tf ) ………… (4)? ? 按照有限元法的理论[3]，将原始方程变换成积分形式，它是基于求一个量的最小值将导致真实解的概念。变分原理的等价定理指出，求解原始微分方程等价于求下式定义的泛函Ie (T)的最小值[4]。即使? ? Ie(Ti)/ Ti=∫∫∫ { ax· 2t/ x2 + ay· 2t/ y2 + az· 2t/ z2 +? ? △R? ? ? ? ( t/ τ- θ/ τ)}dxdydz +∫∫ （λ/ρc） t/ n·ds = 0 ………… (5)? ? △C ? ? 式中△R是单元所包含的子域，△C是在表面C上的面积，显然只有靠近边界C的单元才会出现这一项。(5)式对于所有网格点i成立，便构成一个联立方程系统，求解这个系统可得到整个问题的数值解，即温度场分布。? ? 3 ANSYS软件稳态热分析的基本过程? ? ANSYS软件包仿真热分析的基本过程包括三个步骤：? ? 第一步骤为前处理，其任务是建模。首先确定文件名、分析标题、单位，然后确定软件单元类型、实常数、材料热性能参数，最后创建几何模型并划分网格。? ? 第二步骤为加载计算。在这一过程中，首先确定文件分析类型，然后对该模型施加载荷，最后对该模型进行求解。? ? 第三步骤为后处理，其作用是查看结果。在这一过程中，ANSYS软件可用后处理软件查看由上一步计算好的数据。? ? 4 结果与验证? ? 本次我们仿真的对象为韶山7型（SS7D）电力机车主整流器功率半导体器件用散热器。由前面的分析可知，ANSYS软件包对热场的分析是基于有限元计算，若有限元分析计算量很大，轻则耗用的计算时间长，重则计算机无法正常工作。为此，在保证精度的情况下，对所模拟的散热器模型进行简化，是非常必要的。因此，我们结合韶山SS7D机车主整流器功率半导体器件用散热器的结构形状对称性，所选模拟模型为实际结构整体的1/4(如图1所示)。图2和图3分别是韶山SS7D用散热器在被施加1393W和1979W功率时，所得到的数值仿真温度场分