本研究题为“Electro-Thermo-Mechanical Analysis and Modeling of High-Power Intergrated Gate Commutated Thyristors”，由 Xin Yan, Zhanqing Yu, Lu Qu, Zhizheng Gan, Yulong Huang, 和 Jian Feng 共同完成。所有作者均来自清华大学电力系统及大型发电设备安全控制和仿真国家重点实验室。该研究作为一篇完整的原创研究论文，发表于 IEEE Transactions on Power Electronics 期刊，第39卷第6期，出版于2024年6月。
学术背景与研究动机 本研究属于电力电子器件物理特性分析与建模领域，具体聚焦于大功率集成门极换流晶闸管（IGCT, Intergrated Gate Commutated Thyristor）。IGCT作为一种全硅片（whole-wafer）器件，因其具有强关断能力、高浪涌电流能力和高耐压能力等优点，在直流断路器和逆变器等大功率场合得到越来越广泛的应用。然而，其全硅片结构也带来了一个关键挑战：器件内部的电流分布和结温分布无法直接测量。这种“黑箱”状态使得设计者无法准确获知芯片的工作状况，增加了器件失效的风险。为了确保可靠性，工程实践中通常采用增加器件裕量（如串并联）的方法，但这显著增加了系统的体积和成本。因此，精确理解IGCT芯片内部的物理状态，特别是电流、温度和机械压力之间的耦合关系，并建立能够预测其分布的模型，对于优化器件设计、提高利用率和可靠性、降低成本具有至关重要的意义。现有研究主要集中于绝缘栅双极型晶体管（IGBT），而IGCT与IGBT在结构和封装上存在根本差异，导致IGBT的研究方法和结论无法直接应用于IGCT。此前针对类似结构门极可关断晶闸管（GTO）的研究，也多关注动态过程或采用破坏性、精度有限的测量方法，缺乏对稳态导通状态下芯片内部多物理场耦合的综合研究。因此，本研究旨在填补这一空白，提出一种创新的方法来研究和建模IGCT芯片的内部物理特性。
详细研究流程与方法 本研究包含一系列相互关联的步骤，从芯片特性分析到多物理场耦合模型的建立与验证。
第一步：提出并验证分环建模方法的必要性与有效性。 研究首先分析了IGCT芯片的结构特点。与由多个独立小芯片组成的压接式IGBT不同，IGCT芯片是单一的、完整的圆形晶圆。其上的元胞（cell）沿圆周均匀分布，但在径向（由内到外）可能因制造工艺和位置不同而存在差异。基于此观察，研究提出了“分环建模”的核心思想：将整个芯片沿径向划分为若干个同心环（文中分为10环），并假设同一环上的元胞特性相同，而不同环的元胞特性可以不同。为了验证这一假设，研究团队对芯片上不同位置的单个元胞进行了电流-电压（I-V）特性测试。结果明确显示，不同环上元胞的I-V曲线存在显著差异，外围环的元胞等效电阻高于内部环，证实了径向非均匀性的存在。同时，研究考虑了封装压力分布的影响。通过仿真分析发现，即使在理想压接下，芯片径向的压力分布也是不均匀的，而阴极钼片的高度差异会进一步加剧这种压力分布的不均。这种压力不均会导致接触热阻的差异，进而影响结温和电流分布。这些实验和仿真结果共同证明了传统将芯片视为均一整体的模型是不精确的，而基于圆周均匀、径向分环的建模方法（Split-ring modeling）是必要且有效的。为了便于对每个环进行独立研究，研究团队创新性地采用了激光切割技术，将阴极钼片分割成环，使得可以对每个环的元胞特性进行单独测试。
第二步：多物理场耦合关系分析与数学建模。 在确立分环建模思路后，研究深入分析了电、热、力三个物理场之间的相互耦合关系，并建立了相应的数学表征。 1. 热电耦合：温度显著影响芯片的电学特性（如通态压降），而电学特性产生的损耗又决定了热源，从而影响温度分布。研究从半导体物理的PIN二极管模型出发，推导了通态压降与温度的理论关系。研究表明，描述PN结电压和基区压降的参数与温度近似呈二次关系。基于此，研究改进了IGCT数据手册中常用的ABCD经验模型（Vt = A + B*It + C*ln(It+1) + D*√It），将其中参数C和D表征为温度T的二次函数。通过TCAD软件仿真不同温度下的元胞I-V特性，并拟合参数C和D，验证了其与温度之间的二次关系。这样，就通过可测量的模型参数将温度对电特性的影响进行了量化。 2. 力热耦合：压力的影响主要通过接触热阻体现。接触热阻是IGCT总热阻的重要组成部分（可达50%），它受到接触压力和表面粗糙度的强烈影响。研究介绍了基于接触力学和传热学的理论模型，如Bahrami和Yovanovich的公式，阐述了接触热导与压力、材料微硬度、表面形貌参数之间的关系。仿真结果表明，随着压力增大或粗糙度降低，接触热阻显著减小。这意味着压力分布的不均会直接导致散热路径热阻的差异，从而影响结温分布。反之，温升引起的热应力在IGCT封装中通常由碟形弹簧吸收，因此温度对压力分布的影响可忽略。 3. 电力耦合：压力通过影响接触电阻来影响电学特性，主要体现于阴极铝层与钼片之间的界面。研究通过COMSOL仿真分析了该接触电阻与压力的关系，发现其服从经验公式。然而，仿真也表明，该接触电阻的量级（纳欧到微欧）远小于芯片体电阻（毫欧级），因此在实际建模中其影响可以忽略。同样，电场产生的应力远小于外部机械应力，故电特性对压力分布的影响也不予考虑。
第三步：多物理场耦合仿真模型的构建。 基于上述分析，研究在MATLAB/Simulink软件中构建了一个场-路耦合的多物理场仿真模型。 1. 电学特性建模：利用激光切割阴极钼片后单独测试获得的不同环元胞在不同温度下的I-V特性曲线，使用改进的ABCD模型（参数C、D为温度函数）进行拟合，得到了每个环独立的、包含温度效应的电学模型参数（论文表I）。 2. 热网络建模：基于Cauer模型，为整个IGCT封装建立了集总参数热阻-热容网络。网络考虑了几何结构和材料属性（论文表II）。特别关键的是，将第二步中分析的、与压力相关的接触热阻模型集成到热网络中。研究使用白光干涉仪实测了芯片阴极梳条等关键表面的粗糙度（论文图12，表III），作为接触热阻计算的输入。最终，为分环模型中的每个环建立了对应的热网络，并考虑了环间的横向热耦合（论文图13）。通过与三维有限元仿真结果对比（论文图14），验证了该热网络模型的准确性。 3. 耦合模型集成：在Simulink中搭建了完整的耦合仿真模型（论文图15）。模型的外部输入为每个环的压力、环境温度和总电流。压力决定了各环的接触热阻，从而影响热网络。热网络计算出的各环节点温度，反馈给各环的电学ABCD模型，决定其I-V特性。各环的电学特性决定了在总电流约束下的电流分配（即各环电流）。各环的电流和压降又决定了该环的功率损耗，作为热源输入回热网络。如此，电、热、力三个场实现了闭环耦合求解。
第四步：仿真分析与实验验证。 利用所建模型，研究对两种典型工况进行了仿真，并通过实验验证了模型的准确性。 1. 浪涌电流工况：仿真了环境温度110°C下10kA浪涌电流的情况（论文图16）。结果显示，大部分电流流经第8和第9环，但电流密度和最高结温却出现在第5环。这表明失效风险最高的区域不一定是电流绝对值最大的区域。在浪涌条件下，芯片内部电流分布不均，最高结温差可达20°C。 2. 长期导通工况：模拟了类似固态断路器的稳态工作条件。仿真表明，稳态下的电流分布与浪涌时不同，第9环电流显著高于其他环。由于热时间常数的作用，结温分布更为均匀，最大温差仅约6°C。 3. 实验验证：通过浪涌测试，将实测的器件整体I-V特性曲线、不同温度下的峰值通态压降与模型仿真结果进行对比（论文图17）。两者吻合度非常高，最大电压偏差小于0.2V。此外，长期导通下的通态压降仿真与实测结果误差也在10%以内（论文图18）。这些对比充分验证了所提出的多物理场耦合模型的准确性和有效性。 4. 失效定位验证：研究进行了关断失效测试（论文图20）。仿真预测的风险最高区域（第5环）与实际发生失效的位置完全一致，而电流绝对值较大的第9环并未失效。这强有力地证明了模型在预测器件薄弱点和失效位置方面的实用价值。
研究结论与意义 本研究成功提出并验证了一种用于分析大功率IGCT芯片内部物理状态的分环建模方法及多物理场耦合模型。主要结论包括：1）IGCT芯片在径向存在固有的电、热、力学特性不均匀性，分环建模是精确表征其行为的有效手段。2）电、热、力三场之间存在紧密耦合，其中热电耦合和力热耦合是关键，本研究通过改进的ABCD模型和接触热阻模型量化了这些耦合关系。3）所建立的MATLAB/Simulink场-路耦合模型能够高效、准确地仿真不同工况下芯片内部的电流分布和结温分布。4）仿真揭示，在浪涌条件下芯片内部电流和温度分布极不均匀，失效风险点（第5环）并非电流最大点（第8、9环）；而在稳态条件下，分布则相对均匀。
本研究的科学价值在于首次对IGCT芯片内部的多物理场耦合机制进行了系统性、量化的研究，建立了一套完整的建模与分析框架，填补了该领域的研究空白。其应用价值尤为突出：该模型使设计者能够“透视”器件内部状态，精确评估芯片在不同工作条件下的安全裕度和失效风险，从而指导优化器件设计（如优化元胞布局、压力分布）、提高现有器件的利用率和可靠性，并减少为保安全而过度设计的成本。它为IGCT在直流断路器、大功率变流器等关键设备中的更可靠、更经济应用提供了重要的理论工具和设计依据。
研究亮点 1. 方法创新性：提出了针对全硅片IGCT芯片的“分环建模”思想，并通过激光切割阴极钼片的实验方法实现了对单个环特性的独立研究，构思巧妙且切实可行。 2. 模型系统性：首次建立了综合考虑电、热、机械压力耦合效应的IGCT芯片多物理场模型，耦合关系分析全面、数学表征清晰。 3. 实用高效性：在MATLAB/Simulink中实现的场-路耦合模型，相比纯有限元仿真大幅提高了计算速度，更适用于工程设计和工况分析，实现了精度与效率的平衡。 4. 验证充分性：研究不仅通过特性曲线对比验证了模型精度，更通过失效定位测试，直观且有力地证明了模型预测实际失效风险的能力，使研究成果具有很高的可信度和实用价值。 5. 发现启发性：研究揭示了“最大电流环”与“最高风险环”可能分离的现象，强调了仅关注总电流或局部最大电流的局限性，指出结温与电流密度相结合的综合评估才是关键，这对电力电子器件的可靠性设计具有普遍启发意义。