1、直流母线在中压应用中的综合设计马丁 Buschendorf， 马可神户，罗德里戈阿尔瓦雷斯，斯特芬 Bernet电力电子集团 德累斯顿技术大学（TUD） 德国德累斯顿 martin.buschendorf TU-dresden.de摘要：在大功率中压变频器的应用中，低电感母线的设计是为了实现快速切换，较高效率，及过电压时的安全操作。本文将，采用 PEEC 仿真试验台的简化直流母线的分析计算，与最终设计的测量做以比较。它还提供了，在半导体故障及其合力的情况下短路的计算。.引言转换器在所有工作点的高效率和安全操作，是对作为可再生能源（4，5） ，工业加工，或牵引6的现代转换器（1，2，3）在应用中
2、的主要设计目标的。直流母线的漏电感直接影响 IGBT 的损耗7。漏电感在关断瞬间的影响尤为重要，关断瞬间的高 di/dt 结合大量漏电感产生一个电压峰值。这个峰值增大了 IGBT 的关断损耗，并可能因过电压故障（8，9）毁坏设备。用于保护设备的方法包括有源钳位10和减小电流斜率11。这两种方法通常会增加总的转换损耗和降低效率。因为这些原因，转换器低电感直流母线的设计就显得至关重要。实现低电感直流母线的一般设计规则如文献12, 13所示：最大限度地减少导体的厚度， 最大限度地提高铜的含量，从而降低电感和电阻 通过大量重叠使用多层母线。本文的重点是，对最多三个串联的 4.5 kV/1200AIGB
3、T 的中压试验台直流母线的综合设计。直流侧电容为 3.25mF 时，该系统的最大电压大约是 13.5kv 与的。为达到这些高能量和电压等级，有串并联连接电容器的电路是必要的。对最小的漏电感和最少量的铜，一个适当直流母线的几何形状是一个“U”形，其中，正负电极在图 3a 表示。该试验台的一个特点是漏电感可在 322 nH 的和595nH 之间连续变化。.发展现状文献显示了一些直流母线进行计算和模拟的方法。部分元等效电路（PEEC）法与有限元法（FEM）的比较所示于14。在该出版物中，两片之间的成果由 PEEC 法计算，并和通过有限元分析计算的结果进行比较。该 PEEC 计算是将棒细分成若干元素，
4、这样它可以看作是直的薄的导体并联。它表明，计算给出的结果一致，但 PEEC 法更容易和更快。中压变频器叠层母线的设计计算示于13和15。它表明，从一个简单的铜片开始，两个重叠的母线比单一母线的漏电感低得多。在其中，关于更复杂的几何形状和电流分布的多层母线被提及。据指出，母线设计在很大程度上依赖于，所研究的包括孔、穴的拓扑的特定结构，而且要准确预测漏电感，仿真模拟是十分必要的。大功率转换器直流母线的最优化记录在16。因此 Simulink 模型和 PEEC法被用来计算电流分布和考虑集肤效应及磁场的漏电感。叠层母线在大功率转换器应用的研究中被特别提到。文献17表示了如何通过使用 PEEC 优化一条
5、母线的方法，而且电路解算器仅和固定网格一起使用。它表明，电流平衡导致了漏电感的最小化。在简单结构中，它可由在母线设计中仅增加一条狭缝来完成。在更复杂的结构中，这种直观的方法可能会失败。.数学建模一 漏电感的计算如图 1 所示，对直流母线设计结构的数学建模可简化为一个单一带状导体。考虑到此带状导体的高度 h 和长度 d 比其宽度 W 和厚度大得多，且忽略顶部和底部边缘上的边界效应，这个问题可简化为二维问题。它可以通过，使用所提出的平而薄的差动导体环路，从而解析地解决，参见图 2。图 1 直流母线带状导体的简化模型图 2 直流母线的二维模型 根据毕奥 - 萨伐尔定律，则磁通密度 和磁通量 可被确定
6、，B其中， 是真空磁导率。求解此例的方程，则漏电感计算为，0机械设计提出的几何形状几何形状，如图 1 所示，最小宽度 W=25mm 和最大宽度 W=265mm 时，分析计算的漏电感所示于表。表 分析计算的漏电感这种方法的优点是使估算漏电感变得简单。最显著的缺点是，大幅度简化了机械设计且忽略了边界效应。电流密度，对更高频率的集肤效应和带状导体的高度也未代表。特别是对带状导体两个长边之间的短距离来说，有一个无边界效应的，几乎均匀的磁场的假设是比较好的。然而，由于磁耦合，在两个相对面板之间电流的均匀变得更重要。在铜板之间的长距离，假设和实际的偏差更加重要。该磁场和由分析计算给出的并不相同。因此，由于
7、带状导体顶部和底部的不均匀，导致计算漏电感的假性降低。二 短路电流和磁力的计算 直流短路时，直流母线的机械力与设计极为相关。直流母线的固有振动频率可用漏电感的最小计算值和直流连接电容计算。直流应用里瞬态短路电流 的最大振幅是18:pi根据18， 是母线和电容的电阻值， 是冲击系数。3.2R0.5计算两个平行导体的洛伦兹力 F 并考虑材料属性，可计算出直流母线的厚度是其中 是动态对静态弯曲应力比. 是主要导体应力系数， 是形状系数，V q它可以在19中找到， 是两个支撑点之间的最大距离， 是直流母线的高l h度， 是短路时最大变形为 0.2时，0.2 - 铜的明显拐点。短路的数0,2R学计算和它
8、的合力对简单方程有益。其缺点与它冲突，即考虑这几个方面影响的耦合。.采用 PEEC 的模拟为了更好的估计直流母线的实际漏电感，我们将之前提到的 PEEC 法几乎未简化的来模拟此漏电感。仿真工具 被应用于此。它包含20:InCaD集肤效应，邻近效应， 边界效应和 连接电阻。给母线搭配更精确的模型，则漏电感结果应更为准确。它还提供了一个全球性的解决方案，即允许对这个问题从总体上考虑漏电感，短路计算和磁力。图 3 直流母线的机械设计 a）CAD 草图 b） 采用 的3InCaDPEEC 网格 c）电流密度图 3 示出了直流母线机械设计的 CAD 草图（图 3a）和 PEEC 网格图（图 3b） ，包
9、括为了绝缘而必要的孔，为电容器和 IGBT 的前方端子串并联的板层。这个仿真模型提供了一个阻抗矩阵，它允许在该配置中，直流母线不同部分漏电感的计算。图 3c 给出了直流母线在 时电流分布的仿真结果。min25w对最终的仿真，电容器之间的距离可以从 变化到 。min25max265w表中几个距离下模拟漏电感的结果示于如下的分析计算中。表 采用 PEEC 模拟的漏电感最小和最大距离的计算结果是相比于漏电感的解析计算，-7的最小距离的的和 49的最大距离间的差异为我们所注意。这些差异是由于数学模型的简化造成的。尤其是短距离并行母线（ ）和长距离非均匀磁场在结果中引起了这些差异。正如已经提25wm到的
10、，由于母线顶部和的底部的边界效应，仿真中长距离的漏电感明显高于分析计算的漏电感。在图 3c 中可以看出，正如假定的分析计算，电流密度是不均匀的。在铜母线的前端，电流密度因导体宽度的减小而剧烈增加。这种影响不是能由分析计算表示。电容器端子电流密度的增加不能由在 0 段的简单带状导体表示。这两种影响导致漏电感的增加。. 实验验证测量是去核实分析计算，更准确的说是去核实模拟结果。测试电路和最终的机械结构示于图 4，直流电容器充电到 60v，并被两个 21.7kv/2.4 kA 的 IGBT短路。因此，直流电容器和电感的谐振电路振荡。其中电压和电流采用 10 GS / s 的示波器，600MHz 低电
11、压探针和 16 MHz 的 Rogowski 线圈测量，见图 5。考虑到阻尼 和短路电流频率 ，漏电感 可计算为：rfL图 4 a)漏电感的测量电路 b）最终机械设计漏电感在几个距离下的测量值于表和图 6 可观察到。该过程被重复数次，结果是几次测量的算术平均值。表 漏电感测量值图 5 , 时的电流电压测量值0DCUV3.25DCmF集电极 - 发射极电压， 集电极电流图 6 相比于电容器之间的距离的漏电感，分析计算，模拟， 测量 可以看出，在大距离时，测量和模拟漏电感值间的差别小于 10。很明显，PEEC 法非常适合估计该试验台的母线漏电感。对短距离，测量漏电感明显小与模拟结果。可以假定，母线
12、片之间距离短的接触电压降的影响比距离长的更相关。因此，阻尼常数增加且漏电感减小。这种情况的问题是如何在实际结构中有力的拧紧螺丝在下影响强。若不失败，实验设计的漏电感值总是比理想电气联系下的模拟值低一点。此外，PEEC 模拟使得 在各种情况下的估计值有L低于 10的错误率，相比与解析计算而言，这是一个巨大的进步。.结论文中给出了，确定一台中等电压试验台或一台变换器直流母线漏电感方法的概述。它展示了使用简化机械设计漏电感的解析计算。也计算出短路电流和机械强度。它将使用仿真程序 PEEC 分析的计算结果，和直流母线最终结构在各部分间不同距离的测量值做以比较。测量和解析计算之间的比较表明，因为机械设计
13、的必要简化及边界效应的忽略为，所以解析结果只是对大约 的10m粗略估计。对于较小的距离计算漏电感过高（ ） ，,amin(25)%L对较大的距离计算漏电感过低（ ） 。预定义为,amax(26)49L机械直流连接时，此模拟和它相比给出了一个相当不错的结果。对于小距离，模拟和测量漏电感值之间的差异降低到 ，且对于,ain1较大的距离，此值降低到 。,amax(0)5%.参考文献1 罗德里格斯,J.;Bernet,S.;吴斌;Pontt ,JO.;航路，S.;,“工业中压变频器多级电压源转换器的拓扑结构， ” 工业电子， IEEE 学报，第 54 卷，第 6期， pp.2930 -2945 ， 2
14、007 年 12 月 2 航路，S.;马林诺夫斯基，M.;戈帕库玛，K.;宝， J.; Franquelo，LG ;吴斌;罗德里格斯， J.;佩雷斯，MA;莱昂，姬; ， “多电平变换器的新进展和工业应用， ”工业电子，IEEE 学报， 第 57 卷，8 号，pp.2553 -2580 ，2010 年8 月 3 阿布 - 擦，H.;霍尔茨，J.;罗德里格斯， J.;葛宝明; ， “中压多电平换器艺术现状，挑战，以及在工业应用的要求， ”工业电子，IEEE 学报，第 57卷，8 号，pp.2581 -2596， 2010 年 8 月 4 卡拉斯科，JM ; Franquelo，LG ; Bial
15、asiewicz，JT ;高尔文，E.;Guisado，RCP，;普拉斯，Ma.AM ;莱昂，姬;莫雷诺 - 阿方索，N.;,“可再生能源并网：一项调查的电力电子系统， ” 工业电子，IEEE 学报，第 53 卷，第 4 期，pp.1002 -1016，2006 年 6 月 5 ahavi，JSA ; Kanagapriya，T.; Seyezhai，R.;， “可再生能源交错式升压转换器的设计与分析， ”计算机，电子和电气技术（ ICCEET ） ，2012 年国际会议，第 53 卷，第 4 期，p.447 - 451，2012 年 3 月 21-22 日 6 Bernet，S.; ， “大
16、功率转换器在工业和牵引应用的近期发展， ”电力电子技术， IEEE 学报，第 15 卷，第 6 期，pp.1102 -1117 ，2000 年 11 月 7 阿尔瓦雷斯，R.; Filsecker，F.; Bernet，S.;， “中压转换器新 4.5kv压力包 SPT + IGBT 的描述， ”能源转换会议博览会，2009 。ECCE2009 。IEEE学报，pp.3954 - 3962，2009 年 9 月 20-24 日 8 阿尔瓦雷斯，M. LOOI 和 H.Goh， “绝缘栅双极晶体管在过电压情况下的故障分析， ”国际可再生能源和电能质量会议，圣地亚哥德波斯特拉，2012 9 布劳恩
17、，D.; Pixler，D.;李梅，P.;， “ IGBT 模块的破裂分类和测试， ”工业应用研讨会，1997。第三十二届国际会计准则年度会议，国际会计准则97。IEEE 学报 1997 年会议记录，第 2 卷 pp.1259 -1266，1997 年 10 月 5-9 日 10 塞兹 J.; MERMET，M.;弗雷，D.; JEANNIN，PO ; Schanen ，JL ;Muszicki，P. “IGBT 串联的有源钳位电路的优化和整合， ”工业应用大会，2001 年。第三十六届国际会计准则届年会。I2001 年 EEE 会议记录，第 2 卷， pp.1046 -1051 2 卷， 2
18、001 年 9 月 30 日 11 Petterteig，A.;洛德，J.; Undeland ，T.M; ， “使用电容缓冲器硬开时IGBT的关断损耗， ” 1991年IEEE工业应用学会年会会议记录， 1991年，第2卷，pp.1501 -1507，1991年9月28- 10日 12 Schanen ，JL ;克拉韦尔，E.; Roudet，J.; ， “低电感母线连接的建模， ” IEEE学报，工业应用杂志，第2卷，第5号，pp.39 43，1996年9/10月 13 Guichon ，JM ;艾梅，J.; Schanen ，JL，马丁，C.; Roudet J.;Clavell，E.;
19、 Arpilliere，M.; Pasterczyk，R.;Le Floch，Y.; ， “母线设计：如何节约毫微亨？， ”工业应用会议，2006年国际第41届会计准则年度会议。2006年IEEE会议记录，第四卷，没，pp.1865 - 1869，2006年10月8-12日 14 皮埃特，N.;马雷夏尔，Y.;克拉韦尔，E.; ， “母线技术电动力的计算：部分元等效电路法（ PEEC ）和有限元法（FEM ）的比较 ， ”工业应用研讨会，1998。第三十三届国际会计年度会议。IEEE学报 ，第二卷，pp.921 -924，199810月12-15日15安藤，M.;和田，K.;隆夫 K.;金井，
20、T.;西泽，S.;大桥，H.; ， “中压逆变器母线结构的设计与分析， ”2011年第14届欧洲会议关于电力电子技术与应用的提议（EPE 2011）PP.1 - 10 ，2011年8月30日-2011年9月1日16科比 AT ; Vadlapati，KK ，斯塔基，J.P; Goldney，AP ;Kandilidis，SY ; Hinchley，DA ; ， “大功率层压母线的在高电流分布， ”2011年 - 第14届欧洲电力电子技术与应用（EPE 2011）会议，PP.1 - 10 ， 2011年8月30日- 2011年9月1日17 Pasterczyk ，RJ.;马丁，C.; Guich
21、on ,J.-M. ; Schanen ，J.-L. ; ， “关于分流和漏电感最小化平面母线的优化,”2005年欧洲会议电力电子与应用，pp9 18 DIN VDE 102-10: Kurzschlussstrme in GleichstromEigenbedarfsanlagen vonKraftwerken und Schaltanlagen, Teil1:Berechnung der Kurzschlussstrme, 1998.19 DIN VDE 103-10: Kurzschlussstrme in GleichstromEigenbedarfsanlagen von Kraft
22、werken und Schaltanlagen, Teil 2:Berechnung der Wirkung, 1998.20 Y. Le Floch, S. Cadeau-Belliard, J.-P. Gonnet and V. Mazauric, “Onthe InCa trail.,“Flux Magazine, pp. 10-11, 1 June 2006.Comprehensive Design of DC Busbars for Medium Voltage Applications Martin Buschendorf, Marco Kbe, Rodrigo Alvarez,
23、 Steffen Bernet Power Electronics Group Technical University of Dresden (TUD) Dresden, Germany martin.buschendorftu-dresden.de AbstractIn high power medium voltage converter applications a low inductance busbar is designed in order to achieve fast switching, high efficiency and safe operation with l
24、ow overvoltage. This paper compares an analytical calculation of a simplified DC busbar of a test bench with a PEEC simulation and the measurement of the final design. It also provides the short circuit calculation in case of semiconductor failure and its resulting forces. I. INTRODUCTION High effic
25、iency and safe operation of the converter at all operating points are the main design goals for modern converters ( 1, 2, 3) in applications as regenerative energy ( 4, 5), industrial processes, or traction 6. The stray inductance of the DC busbar directly influences the IGBT losses 7. The effect of
26、 the stray inductance is particularly important during the turn-off transient where a high combined with a large stray inductance results in a voltage peak. dti/This peak increases IGBT turn off losses and could even destroy the device due to overvoltage failure ( 8, 9). Methods used to protect the
27、device include active clamping 10 and reduction of the current slope 11. Both methods typically increase total converter losses and reduce efficiency. For these reasons the design of a low inductive DC busbar becomes essential for the converter. Some general design rules to achieve a low inductance
28、DC busbar are shown in 12, 13: minimize conductor thickness, maximize copper volume, which decreases inductance and resistance and use multilayer busbars with a high amount of overlap. This paper is focused on the comprehensive design of the DC busbar for a medium voltage test bench for the series c
29、onnection of up to three 4.5 kV/1200 A IGBTs. The maximum voltage of the system is about 13.5 kV with a DC-link capacitance of 3.25 mF . To reach these high energy and voltage levels a circuit of series and parallel connected capacitors is required. For minimal stray inductance and also minimal amou
30、nt of copper, a proper DC busbar geometry is a “U”-form, where the plus and minus poles are at the front- see Figure 3a. A special feature of the test bench is the continuously variable stray inductance between 322 nH and595 nH . II. STATE OF THE ART The literature shows some methods for calculation
31、 and simulation of DC busbars. The Partial Element Equivalent Circuit (PEEC) method compared to Finite Element Method (FEM) is shown in 14. In this publication the effort between two sheets is calculated with PEEC method and compared to the result calculated by a FEM-analysis. The PEEC calculation i
32、s made by subdividing the bars into several elements, which can be seen as straight thin conductors in parallel. It was shown that the computation gives coherent results, but the PEEC is easier and quicker. The design calculation of a laminated busbar for a medium voltage inverter is presented in 13
33、 and 15. It is shown that starting from a simple copper sheet two overlap busbars have a much lower stray inductance than a single one. For a more complex geometry and current distribution multilayer busbars are mentioned. It is stated that the busbar design strongly depends on the specific structur
34、e of the investigated topology including holes and apertures, and that simulation is necessary to accurately predict the stray inductance. The optimization of a DC busbar for high power converters is documented in 16. Therefore a simulink model and the PEEC method are used to calculate the current d
35、istribution and stray inductance with skin effect and magnetic fields in mind. Especially laminated busbars are mentioned in the investigated high power converter application. In 17 a methodology is shown for how a busbar can be optimized by using the PEEC and only the circuit solver is used with a
36、fixed meshing. It is shown that a current balancing leads to a minimization of the stray inductance. In simple structures this can be made by simply adding a slit into the busbar design. In more complex structures this intuitive approach can fail. III. MATHEMATICAL MODELING A. Calculation of the str
37、ay inductance For the mathematical modeling the DC busbar design structure was simplified to a single conductor ribbon as it is shown in Figure 1. Considering that the height and length of the hdribbon are much larger than its width and thickness and neglecting the boundary effects on the wtop and b
38、ottom edge, the problem can be simplified into two dimensions. This can be analytically solved by using a differential flat and thin conductor loop of the proposed ribbon, see Figure 2. Figure 1 Simplified DC busbar as a conductor ribbon Figure 2 Two dimensional model of the DC busbar Using the Biot
39、-Savart law the magnetic flux density and the magnetic flux can be determined Bwhere is the magnetic vacuum permeability. Solving the equation for this example the stray 0inductance can be calculated as For the proposed geometry of the mechanical design, as shown in Figure 1, the analytical calculat
40、ed stray inductance for a width distance of =25 mm and =265 mm are shown in minwmaxTABLE I . The advantage of this method is the simplicity to estimate the stray inductance. The most significant disadvantages are the substantial simplifications of the mechanical design and the neglect of the boundar
41、y effects. The current density, skin effect for higher frequencies and the height of the ribbon are not represented. Especially for short distances between the two long sides of the conductor ribbon the assumption to have a nearly homogenous magnetic field without boundary effects is comparatively g
42、ood. However, due to magnetic coupling, homogenous current in these two facing plates become more important. The deviation between assumption and reality is more important for long distances between the copper plates. The magnetic field is not as homogenous as given by the analytical calculation. He
43、nce the inhomogeneity on the top and the bottom of the ribbon leads to a falsely lower calculated stray inductance. B. Calculation of the short circuit current and the magnetic force The mechanical strength of the DC busbar in case of DC-short circuit is extremely relevant for the desighn.With the m
44、inimal calculated value of the stray inductance and the DC-link capacitance the natural frequency of the DC busbarcan be calculated. The maximum amplitude of the transient short circuit current for DC-piapplications is 18:where is the resistance of the busbar and the capacitors and is the shock 2.3R
45、 5.0coefficient,according 18 .Calculating the Lorentz force for two parallel conductors and considerig the material Fattributes the thickness of the DC busbar can be caoculated to Where is the ratio of dynamic to static bendig stress, is the factor representing the main Vconductor stress and the is
46、the form factor, which can be found in 19, is the maximal distance q lbetween two supporting points, is the height of the DC busbar and is the 0.2%apparent h2.0Ryielding point of the copper where the maximum is 0.2% in case of a short circuit. The mathematical calculation of the short circuit and it
47、s resulting force has the advantage of simple equations. In conflict to this the disadvantage is the uncoupled consideratio of the several effects. IV. SIMULATION WITH PEETo achieve a better estimation of the real stray inductance of the DC busbar, it was simulated with nearly no simplification by t
48、he already mentioned PEEC method .The simulation tool InCa3D is used for this. It contains20: skin effect, proximity effect, boundary effects and connection resistances. Paired with a more precise model of the busbar,the results for the stray inductance should be much more exact.It also provides a g
49、lobal solution methodology which allows a general view on the problem considering stray inductance,short circuit calculation and magnetic force. Figure 3 Mechanical design of the DC busbar a) as CAD sketch b) as PEEC grid in InCa3D c) current densityFigure 3 shows the mechanical DC busbar design as CAD sketch (Fig. 3a) and PEEC grid(Fig.3b) including the necessar
