CN110232206A - 基于ansys的脉冲功率器件的热变形仿真分析方法 - Google Patents

基于ansys的脉冲功率器件的热变形仿真分析方法 Download PDF

Info

Publication number
CN110232206A
CN110232206A CN201910351916.4A CN201910351916A CN110232206A CN 110232206 A CN110232206 A CN 110232206A CN 201910351916 A CN201910351916 A CN 201910351916A CN 110232206 A CN110232206 A CN 110232206A
Authority
CN
China
Prior art keywords
power device
pulsed power
simulation
result
hot
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910351916.4A
Other languages
English (en)
Other versions
CN110232206B (zh
Inventor
王晓锋
何小琦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201910351916.4A priority Critical patent/CN110232206B/zh
Publication of CN110232206A publication Critical patent/CN110232206A/zh
Application granted granted Critical
Publication of CN110232206B publication Critical patent/CN110232206B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

为了解决现有技术中热变形仿真方法并未考虑脉冲功率的变化对温度的影响而导致的仿真结果与实际结果有着较大的差距的问题,本发明提供一种基于ANSYS的脉冲功率器件的热变形仿真分析方法,首先对脉冲功率器件进行热仿真过程,得到温度随时间变化的脉冲功率器件热仿真结果并形成存储文件;然后,进行结构仿真过程,并将上述的脉冲功率器件热仿真结果文件作为载荷施加到结构仿真中,得到脉冲功率器件在各个时间点的热变形结果。本发明所述方法能在物理过程上更接近脉冲功率器件的实际温变过程与散热过程,获得更精确的脉冲功率器件的在不同时间点的温度场、温度随脉冲功率的变化过程和热变形结果。

Description

基于ANSYS的脉冲功率器件的热变形仿真分析方法
技术领域
本发明涉及脉冲功率器件的可靠性领域,特别是一种基于ANSYS的脉冲功率器件的热变形仿真分析方法。
背景技术
随着微电子封装技术的快速发展,电子器件不断向小型化和高度集成化的方向发展。脉冲功率器件的设计、封装和材料等相关技术的不断进步,致使其尺寸更小,高I/O数,电性能好并且易于装配到电路板上。脉冲功率器件在制造和使用中会受到温度、湿气、振动等条件因素的影响,易与造成器件可靠性的失效问题并影响器件的正常使用。有关调查数据表明,温度是造成脉冲功率器件失效的主要原因。由于器件结构中不同材料之间热膨胀系数差异,脉冲功率器件在温变过程中会产生热变形。热变形会影响到器件的共面度[2],引发芯片断裂、芯片分层和焊点失效等质量问题。因此,对于脉冲功率器件建模技术的研究就显得十分必要。
目前,常采用的脉冲功率器件的热变形仿真方法,一般是根据器件的物理结构,分析其平均脉冲功率稳态时的温度分布情况,以便寻找器件的热关键点,获得器件温度关键点的热变形结果,如图7~8所示。这种热变形仿真方法得到的温度分布,并未考虑到脉冲功率的变化对温度的影响,因而最后结果与实际结果有着较大的差距,不利于脉冲功率器件的热变形研究。
发明内容
为了解决现有技术中热变形仿真方法并未考虑脉冲功率的变化对温度的影响而导致的仿真结果与实际结果有着较大的差距的问题,本发明提供一种基于ANSYS的脉冲功率器件的热变形仿真分析方法。
本发明为了解决上述技术体所采用的技术方案是:基于ANSYS的脉冲功率器件的热变形仿真分析方法,其技术方案在于:首先对脉冲功率器件进行热仿真过程,得到温度随时间变化的脉冲功率器件热仿真结果并形成存储文件;然后,进行结构仿真过程,并将上述的脉冲功率器件热仿真结果文件作为载荷施加到结构仿真中,得到脉冲功率器件在各个时间点的热变形结果。
本发明的有益效果是:本发明所述方法包括脉冲功率器件的热仿真和结构仿真;所述脉冲功率器件的热仿真方法,通过选择单元类型、建立仿真模型、材料参数设定、网格划分、热载荷加载与瞬态求解控制设置,计算得到温度随时间变化的仿真结果;所述脉冲功率器件的结构仿真分析,通过清除热分析的载荷约束、热—结构转变、分析类型、材料参数设置、参考温度设置和热分析结果加载、约束条件施加,计算得到脉冲功率器件在各个时间点的热变形结果。本发明所述仿真方法能在物理过程上更接近脉冲功率器件的实际温变过程与散热过程,获得脉冲功率器件在不同时间点的温度场、温度随脉冲功率的变化过程和热变形结果。
附图说明
图1为本发明的流程图。
图2为脉冲功率器件的几何模型。
图3是t=15s时脉冲功率器件的温度分布云图。
图4是t=15s时节点8484的温度变化示意图。
图5是定义约束后的网格模型。
图6是t=15s时脉冲功率器件的热变形示意图。
图7是稳态时脉冲功率器件的温度结果。
图8是稳态时脉冲功率器件的热变形图。
具体实施方式
下面结合附图对本申请进行进一步的说明。
如图1所示,热仿真过程采用间接法进行热-结构耦合分析,包括以下步骤:
S1.选择脉冲功率器件热仿真的热单元类型;
S2.对脉冲功率器件的几何模型进行简化,建立几何模型;
S3.将S2步骤得到的几何模型导到入Ansys中,添加材料参数,进行网格划分并进行局部细化;
S4.进行热载荷的加载和瞬态求解控制设置;
S5.进行热仿真分析求解,得到温度随时间变化的脉冲功率器件热仿真结果;
S6.将S5步骤得到的脉冲功率器件热仿真结果进行保存。
其中,S4步骤中的热载荷的加载和瞬态求解控制设置过程包括以下步骤:
S401.将脉冲功率器件的脉冲功率转换成热生成数据;
S402.制定txt文档,将S401中的热生成数据加载到表格载荷中;
S403.根据脉冲功率器件的脉冲功率的脉宽和频率,设置相对应的时间步长,时间步长的时间 t1取脉冲信号中高电平时间点t2的一半,即:t1=t2/2;输出结果中,设置需要输出结果的频率n,输出结果时间点t3是时间步长的倍数,即:t3=n*t1;其中,n为1,2,3……的自然数。 n的大小由存储输出结果的存储硬盘大小决定,如果硬盘存储器空间足够,可以将n取值1,此时,存储输出结果的频率最高,存储的数据最多。在实际操作中,可以根据情况确定n的取值。
其中,结构仿真过程为静态分析,包括以下步骤:
S7.获取热仿真过程得到的温度随时间变化的脉冲功率器件热仿真结果文件,然后进行热-结构单元转换,将热单元转换成对应的结构单元,并删除热仿真过程中的所有加载和约束;
S8.定义结构仿真的类型、材料参数设置、参考温度设置,施加约束载荷,将S7步骤获取的温度随时间变化的脉冲功率器件热仿真结果文件作为载荷施加到结构仿真中,选择时间点T 的热分析结果作为热载荷加载到结果模型上;按照模型的结构,将底面中心的四个节点进行全约束;
S9.进行脉冲功率器件的结构仿真分析,得到脉冲功率器件在各个时间点的热变形结果。
具体实施方式I:本实施例结合某脉冲功率器件为仿真对象,提出基于ANSYS的热变形仿真方法。脉冲功率器件的电参数:脉冲功率为2W,脉冲周期是1ms,脉冲宽度是100μs。参照图1所示流程对脉冲功率器件的热变形仿真进行分析与建模。具体步骤如下:
步骤一,本次分析选用间接法进行热-结构耦合分析,选用Solid70热单元;
步骤二,对脉冲功率器件几何模型进行简化,去掉器件中对热变形影响较小的结构,如细的引线、引脚、焊点等细小结构;如图2,仿真模型结构为:基板(D)、粘接层(C)、芯片(B)、塑封层(A);
步骤三,添加模型的材料参数,包括芯片的弹性模量131Gpa、泊松比0.28、密度2300kg/m3、热膨胀系数2.8×10-6m/K,热导率131W/(m·K);粘接层的弹性模量7GPa,泊松比0.34,密度2500kg/m3,热膨胀系数1×10-4m/K,热导率0.2W/(m·K);塑封层的弹性模量20GPa,泊松比0.25,密度1900kg/m3,热膨胀系数4×10-5m/K,热导率0.25W/(m·K);基板的弹性模量Ex=Ey=7Gpa、Ez=7.5Gpa,Gxy=Gyz=2.4Gpa、Gxz=3.6Gpa,泊松比 Uxy=Uyz=0.42、Uxz=0.11;密度1700kg/m3,热膨胀系数Ax=Az=15×10-6m/K、Ay=25× 10-6m/K,热导率0.2W/(m·K);
步骤四,热载荷的施加,包括芯片上脉冲载荷的施加、初始温度为25℃、器件外表面对流系数设定为125w/m2;瞬态求解控制设置,时间总长(Time at end of loadsteps)是15s,时间步长(time step size)是0.00048s,结果输出的频率数(Write every Nthsubsteps) 是n=2;
步骤五,按照上述方式定义完成后首先进行热分析求解,得到脉冲功率器件的瞬态热仿真结果。t=15s时脉冲功率器件的温度场如图3所示:芯片对应的上表面温度最高(34.7℃),器件边角处的温度最低(28℃);图4为t=15s时节点8484的温度变化示意图:初始时温度增长较快,最后逐渐趋于稳定;
步骤六,热分析完成后,保存热仿真结果文件(*.rth);进行热-结构单元的转换,将热单元Solid70转换成Solid185结构单元;删除热分析中定义的所有加载和约束;定义结构分析为静态分析、参考温度为25℃,定义约束载荷,将底面中心的四个节点进行全约束,如图5所示;将某一时间点热分析结果文件作为载荷施加到结构分析中;
步骤七,进行脉冲功率器件的结构分析仿真,得到热变形的仿真结果。进行脉冲功率器件的结构分析仿真,得到热变形的仿真结果。t=15s时脉冲功率器件的热变形如图6所示:器件热变形呈现出两边向上翘起、中心下凹的变形,变形从中心约束点到边角处逐渐增大,变形值为1.4μm。
至此,本流程结束。
总之,本发明公开了一种脉冲功率器件的热变形仿真分析方法,包括脉冲功率器件的热仿真和结构仿真;通过选择单元类型、建立仿真模型、材料参数设定、网格划分、热载荷加载与瞬态求解控制设置,计算得到温度(℃)随时间(t)变化的仿真结果;然后通过清除热分析的载荷约束、热—结构转变、分析类型、材料参数设置、参考温度设置和热分析结果加载、约束条件施加,计算得到脉冲功率器件在各个时间点的热变形结果。该仿真方法能在物理过程上更接近脉冲功率器件的实际温变过程与散热过程,获得脉冲功率器件在不同时间点的温度场、温度随脉冲功率的变化过程和热变形结果,使脉冲功率器件的热变形仿真更加接近实际情况,得到更精确的测量结果。
以上所述仅为发明的较佳实施例而己,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.基于ANSYS的脉冲功率器件的热变形仿真分析方法,其特征在于:首先对脉冲功率器件进行热仿真过程,得到温度随时间变化的脉冲功率器件热仿真结果并形成存储文件;然后,进行结构仿真过程,并将上述的脉冲功率器件热仿真结果文件作为载荷施加到结构仿真中,得到脉冲功率器件在各个时间点的热变形结果。
2.根据权利要求1所述的基于ANSYS的脉冲功率器件的热变形仿真分析方法,其特征在于:热仿真过程采用间接法进行热-结构耦合分析,包括以下步骤:
S1.选择脉冲功率器件热仿真的热单元类型;
S2.对脉冲功率器件的几何模型进行简化,建立几何模型;
S3.将S2步骤得到的几何模型导到入Ansys中,添加材料参数,进行网格划分并进行局部细化;
S4.进行热载荷的加载和瞬态求解控制设置;
S5.进行热仿真分析求解,得到温度随时间变化的脉冲功率器件热仿真结果;
S6.将S5步骤得到的脉冲功率器件热仿真结果进行保存。
3.根据权利要求2所述的基于ANSYS的脉冲功率器件的热变形仿真分析方法,其特征在于:S4步骤中的热载荷的加载和瞬态求解控制设置过程包括以下步骤:
S401.将脉冲功率器件的脉冲功率转换成热生成数据;
S402.制定txt文档,将S401中的热生成数据加载到表格载荷中;
S403.根据脉冲功率器件的脉冲功率的脉宽和频率,设置相对应的时间步长,时间步长的时间t1取脉冲信号中高电平时间点t2的一半,即:t1=t2/2;输出结果中,设置需要输出结果的频率n,输出结果时间点t3是时间步长的倍数,即:t3=n*t1;其中,n为其中,n为1,2,3……的自然数。
4.根据权利要求1所述的基于ANSYS的脉冲功率器件的热变形仿真分析方法,其特征在于:结构仿真过程为静态分析,包括以下步骤:
S7.获取热仿真过程得到的温度随时间变化的脉冲功率器件热仿真结果文件,然后进行热-结构单元转换,将热单元转换成对应的结构单元,并删除热仿真过程中的所有加载和约束;
S8.定义结构仿真的类型、材料参数设置、参考温度设置,施加约束载荷,将S7步骤获取的温度随时间变化的脉冲功率器件热仿真结果文件作为载荷施加到结构仿真中,选择时间点T的热分析结果作为热载荷加载到结果模型上;按照模型的结构,将底面中心的四个节点进行全约束;
S9.进行脉冲功率器件的结构仿真分析,得到脉冲功率器件在各个时间点的热变形结果。
CN201910351916.4A 2019-04-28 2019-04-28 基于ansys的脉冲功率器件的热变形仿真分析方法 Active CN110232206B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910351916.4A CN110232206B (zh) 2019-04-28 2019-04-28 基于ansys的脉冲功率器件的热变形仿真分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910351916.4A CN110232206B (zh) 2019-04-28 2019-04-28 基于ansys的脉冲功率器件的热变形仿真分析方法

Publications (2)

Publication Number Publication Date
CN110232206A true CN110232206A (zh) 2019-09-13
CN110232206B CN110232206B (zh) 2023-07-25

Family

ID=67860356

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910351916.4A Active CN110232206B (zh) 2019-04-28 2019-04-28 基于ansys的脉冲功率器件的热变形仿真分析方法

Country Status (1)

Country Link
CN (1) CN110232206B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111104741A (zh) * 2019-12-17 2020-05-05 南方电网科学研究院有限责任公司 一种mmc子模块中金属膜电容电热耦合仿真方法
CN111428408A (zh) * 2020-03-24 2020-07-17 哈尔滨奥瑞德光电技术有限公司 一种基于ansys仿真分析的热弯机用加热板高温变形的计算方法及加热板表面结构
CN112257303A (zh) * 2020-09-14 2021-01-22 中国航空工业集团公司洛阳电光设备研究所 一种基于热仿真模型的温度稳定时间的测试方法
CN112329361A (zh) * 2020-11-18 2021-02-05 天津光电通信技术有限公司 一种基于热仿真分析空气过滤网影响系统散热功能的方法
CN114580222A (zh) * 2022-05-09 2022-06-03 山东大学 用于电子器件的高速瞬态热仿真方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102592019A (zh) * 2012-01-06 2012-07-18 北京航空航天大学 一种石英挠性加速度计焊接残余应力仿真方法
CN103559418A (zh) * 2013-11-20 2014-02-05 北京航空航天大学 一种基于故障物理的电子产品任务可靠度计算方法
CN105631122A (zh) * 2015-12-25 2016-06-01 鼎奇(天津)主轴科技有限公司 一种机床大件的热变形仿真分析与建模方法
CN106295088A (zh) * 2016-11-03 2017-01-04 上海卫星工程研究所 航天器空间机热耦合载荷高精度自适应传递的通用方法
CN106874566A (zh) * 2017-01-18 2017-06-20 西安电子科技大学 一种电子装备结构优化设计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102592019A (zh) * 2012-01-06 2012-07-18 北京航空航天大学 一种石英挠性加速度计焊接残余应力仿真方法
CN103559418A (zh) * 2013-11-20 2014-02-05 北京航空航天大学 一种基于故障物理的电子产品任务可靠度计算方法
CN105631122A (zh) * 2015-12-25 2016-06-01 鼎奇(天津)主轴科技有限公司 一种机床大件的热变形仿真分析与建模方法
CN106295088A (zh) * 2016-11-03 2017-01-04 上海卫星工程研究所 航天器空间机热耦合载荷高精度自适应传递的通用方法
CN106874566A (zh) * 2017-01-18 2017-06-20 西安电子科技大学 一种电子装备结构优化设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
程刚等: "基于Ansys的光学器件热变形仿真与分析", 《应用光学》 *
程刚等: "基于Ansys的光学器件热变形仿真与分析", 《应用光学》, no. 05, 15 September 2008 (2008-09-15) *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111104741A (zh) * 2019-12-17 2020-05-05 南方电网科学研究院有限责任公司 一种mmc子模块中金属膜电容电热耦合仿真方法
CN111104741B (zh) * 2019-12-17 2024-01-26 南方电网科学研究院有限责任公司 一种mmc子模块中金属膜电容电热耦合仿真方法
CN111428408A (zh) * 2020-03-24 2020-07-17 哈尔滨奥瑞德光电技术有限公司 一种基于ansys仿真分析的热弯机用加热板高温变形的计算方法及加热板表面结构
CN112257303A (zh) * 2020-09-14 2021-01-22 中国航空工业集团公司洛阳电光设备研究所 一种基于热仿真模型的温度稳定时间的测试方法
CN112257303B (zh) * 2020-09-14 2024-04-09 中国航空工业集团公司洛阳电光设备研究所 一种基于热仿真模型的温度稳定时间的测试方法
CN112329361A (zh) * 2020-11-18 2021-02-05 天津光电通信技术有限公司 一种基于热仿真分析空气过滤网影响系统散热功能的方法
CN114580222A (zh) * 2022-05-09 2022-06-03 山东大学 用于电子器件的高速瞬态热仿真方法及系统
CN114580222B (zh) * 2022-05-09 2022-07-22 山东大学 用于电子器件的高速瞬态热仿真方法及系统

Also Published As

Publication number Publication date
CN110232206B (zh) 2023-07-25

Similar Documents

Publication Publication Date Title
CN110232206A (zh) 基于ansys的脉冲功率器件的热变形仿真分析方法
Christiaens et al. A generic methodology for deriving compact dynamic thermal models, applied to the PSGA package
Evans et al. Automated fast extraction of compact thermal models for power electronic modules
Hubbard High-resolution modeling of the advance of the Younger Dryas ice sheet and its climate in Scotland
CN102592020B (zh) 一种航空参数处理设备可靠性薄弱环节确定方法
CN111950220B (zh) 一种电热耦合模型建立方法
CN103792476A (zh) 用于半导体器件的热阻测试方法
CN105913494A (zh) 多尺度裂缝精细地质建模及数值模拟方法和装置
CN108614941A (zh) 一种针对集成qfn芯片的板级封装设计优化方法
Menrath et al. Shape optimization of a pin fin heat sink
CN108536931B (zh) 基于贪婪算法的三维片上网络中温度传感器的部署方法
Li et al. End-to-end analysis of integration for thermocouple-based sensors into 3-D ICs
Ali Modeling and simulation of MEMS components: Challenges and possible solutions
Weber et al. Toolset for measuring thermal behavior of FPGA devices
Kaïd et al. Novel method for modelling the local electro-thermal behaviour of discrete power diode using Verilog-A in a standard CAD environment
CN118298897B (zh) 一种内存模组自动断电压力测试方法及系统
Chang et al. Package warpage and stress evaluation for a plastic electronic package
Lasance et al. The use of time-dependent temperature response curves for the generation of (dynamic) compact thermal models
Janicki et al. Dynamic thermal modelling of a power integrated circuit with the application of structure functions
CN113155313B (zh) 一种扇出型封装温度分布原位模拟结构及方法
Janicki et al. Practical study of temperature distribution in a thermal test integrated circuit
Rencz et al. Thermal and Electro-thermal System Simulation
Staszak Heat Management
Huang et al. Physically-based compact thermal modeling—achieving parametrization and boundary condition independence
Poppe et al. Thermal Metrics

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant