CN104318060A - 基于可视化算法的电抗器温升计算方法 - Google Patents
基于可视化算法的电抗器温升计算方法 Download PDFInfo
- Publication number
- CN104318060A CN104318060A CN201410494881.7A CN201410494881A CN104318060A CN 104318060 A CN104318060 A CN 104318060A CN 201410494881 A CN201410494881 A CN 201410494881A CN 104318060 A CN104318060 A CN 104318060A
- Authority
- CN
- China
- Prior art keywords
- reactor
- core
- temperature rise
- formula
- reactor temperature
- 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.)
- Pending
Links
Landscapes
- Soft Magnetic Materials (AREA)
Abstract
本发明公开了一种基于可视化算法的电抗器温升计算方法,本发明中将电抗器铁芯磁路长度l、电抗器铁芯截面积S和电抗器线圈匝数N作为三个变量,在满足电感量和体积优化的条件下由可视化图形可得到所求得的温升数值。本发明的有益效果是:本发明中将电抗器铁芯磁路长度l、电抗器铁芯截面积S和电抗器线圈匝数N作为三个变量,同时满足线圈匝数N和铁芯截面积S的限制条件下获得电抗器温升值,结合变化趋势图获得电抗器铁芯磁路长度l、电抗器铁芯截面积S和电抗器线圈匝数N,本发明通过可视化算法解决了单一化的问题,设计过程相对直观、简单、方便,且可方便观察函数的变化情况。
Description
技术领域
本发明涉及一种电抗器温升计算方法,具体涉及一种基于可视化算法的电抗器温升计算方法。
背景技术
随着电力电子技术的飞速发展,各种电力电子装置在电力系统、工业、交通、家庭等众多领域中的应用日益广泛,由此带来的谐波问题日益严重,所以解决谐波问题是当前最重要的问题。与其它电气元件不同,使用者很难采购到符合自己要求的电感,具体设计一个磁性元件在满足电气性能条件下,需要综合考虑其成本、体积、重量等困难。
常用的整流电路中采用三相桥式结构,虽然电容滤波电路受到许多设计者的要求,但是其输出的电压畸变率较大,会影响用电设备。
发明内容
为解决现有技术的不足,本发明的目的在于提供一种基于可视化算法的电抗器温升计算方法,直观明了,方法简单。
为了实现上述目标,本发明采用如下的技术方案:
基于可视化算法的电抗器温升计算方法,其特征是,包括以下步骤,
1)测量电抗器铁芯磁路长度l、电抗器铁芯截面积S和电抗器线圈匝数N;计算电抗器的磁芯损失公式PFe和铜芯损失公式PCu,其中,磁芯损失公式为: 式中me为磁芯的质量,k、m、n均为材料的系数;
2)计算电抗器的总损耗PΣ,PΣ=PFe+PCu;
3)将步骤2)的结果带入电抗器温升计算公式进行计算;
4)根据法拉第定律转化后得到的工作磁密公式 确定电抗器线圈匝数N和电抗器铁芯截面积S的限制条件NS≥LIPK/Bw,其中Bw为电抗器的磁密,IPK为电流峰峰值;
5)根据步骤3)中的电抗器温升计算公式可得到变化趋势图,将变化趋势图进行平面剖析,在满足步骤4)中的限制条件下找出解集。
前述的基于可视化算法的电抗器温升计算方法,其特征是,所述步骤3)中的电抗器温升公式为Tr=450ψ0.826(℃);其中功率耗散密度 At为表面积,其表达式为:At=KsAP 0.5;AP为面积积,Ks是与磁芯结构有关的常数。
前述的基于可视化算法的电抗器温升计算方法,其特征是,所述步骤5)中的变化趋势图根据步骤4)中的计算结果以及以电抗器铁芯磁路长度l、电抗器铁芯截面积S、电抗器线圈匝数N为变量绘出。
前述的基于可视化算法的电抗器温升计算方法,其特征是,所述变化趋势图的约束条件是:1)所述电抗器铁芯磁路长度l的范围约束条件以硅钢EI型材料最小和最大尺寸为依据,约束条件为0.001m~0.5m;2)所述电抗器线圈匝数N的约束条件为1~100;3)所述电抗器铁芯截面积S的约束条件为0.000001m2~0.005m2。
前述的基于可视化算法的电抗器温升计算方法,其特征是,所述步骤5)中获得电抗器温升平面图的过程为,在满足步骤4)中的限制条件下,依据选择的电抗器铁芯磁路长度l作为定量,在电抗器温升变化趋势图中找出相关满足域,即可得到变量参数值。本发明的有益效果是:本发明中将电抗器铁芯磁路长度l、电抗器铁芯截面积S和电抗器线圈匝数N作为三个变量,同时满足线圈匝数N 和铁芯截面积S的限制条件下获得电抗器温升值,结合变化趋势图获得电抗器铁芯磁路长度l、电抗器铁芯截面积S和电抗器线圈匝数N,本发明通过可视化算法解决了单一化的问题,设计过程相对直观、简单、方便,且可方便观察函数的变化情况。
附图说明
图1是电抗器温升变化趋势图。
图2是电抗器工作磁密数值的条件下,交集域中电抗器温升变化平面趋势图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
常用的整流电路中采用三相桥式结构,虽然电容滤波电路受到许多设计者的要求,但是其输出的电压畸变率较大,会影响用电设备。为了防止整流电路中受到大电流的冲击,在电容滤波电路的基础上串入直流铁芯电抗器,构成滤波电路,滤波电路结构包括整流桥、直流铁芯电抗器、滤波电容和负载。整流桥输出电压不仅含有直流分量,而且含有交流谐波,该谐波对于负载工作及其不利,如果滤波器设计的不好,则电压纹波系数较大,其输出波形不平直,故必须使得整流输出电压纹波系数达到较小。滤波电感在整流电路滤波器中所占的体积相当大,其质量也较重,制约着滤波器的体积小型化。
直流铁芯电抗器主要包括线圈、铁芯、气隙组成。电感是一个能量储存元器件,为了有效地存储和返回能量到电路中去,并要求体积最小,体积最小,成本最低的电感是设计追求的目标。体积最小意味着磁芯利用最好,损耗最小。针对此,设计出了一种基于可视化算法的电抗器温升计算方法,以便于对电抗 器温升进行分析,包括以下步骤,
1)测量电抗器铁芯磁路长度l、电抗器铁芯截面积S和电抗器线圈匝数N;计算电抗器的磁芯损失公式PFe和铜芯损失公式PCu,其中,磁芯损失公式为: 式中me为磁芯的质量,k、m、n均为材料的系数;
2)计算电抗器的总损耗PΣ,PΣ=PFe+PCu;
3)将步骤2)的结果带入电抗器温升计算公式进行计算;电抗器温升公式为Tr=450ψ0.826(℃);其中功率耗散密度 At为表面积,其表达式为:At=KsAP 0.5;AP为面积积,Ks是与磁芯结构有关的常数;
4)根据法拉第定律转化后得到的工作磁密公式 确定电抗器线圈匝数N和电抗器铁芯截面积S的限制条件NS≥LIPK/Bw,其中Bw为电抗器的磁密,IPK为电流峰峰值;
5)根据步骤3)中的电抗器温升计算公式可得到变化趋势图,将变化趋势图进行平面剖析,在满足步骤4)中的限制条件下找出解集。根据步骤4)中的电抗器温升计算公式以及以电抗器铁芯磁路长度l、电抗器铁芯截面积S、电抗器线圈匝数N为变量可得到变化趋势图。其中,变化趋势图的约束条件是:
(1)电抗器铁芯磁路长度l的范围约束条件以硅钢EI型材料最小和最大尺寸为依据,约束条件为0.001m~0.5m;
(2)电抗器线圈匝数N的约束条件为1~100;
(3)电抗器铁芯截面积S的约束条件为0.000001m2~0.005m2。
将变化趋势图进行平面剖析,在满足步骤5)中的限制条件下依据选择的电抗器铁芯磁路长度l作为定量,在电抗器温升变化趋势图中找出相关满足域,即可得到变量参数值,找出解集。
具体的我们带入数据进行实际计算:
根据步骤1),设计采用EI系列厚度为0.3556mm的硅钢片,查铁合金表可知:k=0.000557,m=1.68,n=1.86,BAC=0.25T;其中BAC为交流磁通密度,其频率为100Hz,硅钢的密度ρ为7.3g/cm3,则硅钢的质量为:
m=ρVL≈0.001S·l,VL为电抗器铁芯体积。
所以由上面所提供参数可以求得磁芯损失为:PFe=(0.001S·l)kf(m)BAC (n)=0.000097Sl导线裸面积: 式中I为流过电感器的电流(A);j为电流密度(A/mm2)。
由上式Aw可得Aw=0.0529cm2,查美国线规可知,应选择AWG=#9:裸面积Aw=0.0663cm2,电阻率为
计算绕组电阻:
式中MLT为平均匝长,此种情况下MLT≈1.2l,则铜损Pcu=I2RL=0.01815NlPFe=0.000097Sl,PCu=0.01815Nl。
带入步骤2)得到PΣ=Pcu+PFe=0.01815Nl+0.000097Sl。
电抗器温升公式为Tr=450ψ0.826(℃);其中功率耗散密度 At为表面积,其表达式为:At=KsAP 0.5;AP为面积积,Ks是与磁芯结构有关的常数。在本情况下,查硅钢型号数据表可知叠片磁芯系数Ks=41.3;AP=wa·s,wa=0.0208l2,由此可得:At=41.3(0.0208Sl2)0.5,温升Tr计算公式为:Tr=450ψ0.826(℃),将前述公式代入可得
本方法将电抗器铁芯磁路长度l、电抗器铁芯截面积S和电抗器线圈匝数N作为三个变量,同时满足线圈匝数N和铁芯截面积S的限制条件下获得电抗器 温升值,结合变化趋势图获得电抗器铁芯磁路长度l、电抗器铁芯截面积S和电抗器线圈匝数N的关系图。
结合如上的具体实例,并根据步骤5)采用MATLAB软件编程可得图1,所示为得到的电抗器温升变化趋势图,图中以电抗器铁芯磁路长度l、电抗器铁芯截面积S和电抗器线圈匝数N作为三个变量建立三维坐标,将Tr的数值以颜色的深浅表现出来,对应的数值在图2右侧体现。于是,得到了如图1所示的可视化图形。
图2是对图1进行平面剖析后,在满足步骤5)中的限制条件下依据选择的电抗器铁芯磁路长度l作为定量,在电抗器温升变化趋势图中找出相关满足域,找出的解集,解集用虚线框出。图1和图2的结果很清晰,分析起来直观简单。
本方法通过可视化算法解决了单一化的问题,设计过程相对直观、简单、方便,且可方便观察函数的变化情况。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (5)
1.基于可视化算法的电抗器温升计算方法,其特征是,包括以下步骤,
1)测量电抗器铁芯磁路长度l、电抗器铁芯截面积S和电抗器线圈匝数N;计算电抗器的磁芯损失公式PFe和铜芯损失公式PCu,其中,磁芯损失公式为: 式中me为磁芯的质量,k、m、n均为材料的系数;
2)计算电抗器的总损耗PΣ,PΣ=PFe+PCu;
3)将步骤2)的结果带入电抗器温升计算公式进行计算;
4)根据法拉第定律转化后得到的工作磁密公式 确定电抗器线圈匝数N和电抗器铁芯截面积S的限制条件NS≥LIPK/Bw,其中Bw为电抗器的磁密,IPK为电流峰峰值;
5)根据步骤3)中的电抗器温升计算公式可得到变化趋势图,将变化趋势图进行平面剖析,在满足步骤4)中的限制条件下找出解集。
2.根据权利要求1所述的基于可视化算法的电抗器温升计算方法,其特征是,所述步骤3)中的电抗器温升公式为Tr=450ψ0.826(℃);其中功率耗散密度 At为表面积,其表达式为:At=KsAP 0.5;AP为面积积,Ks是与磁芯结构有关的常数。
3.根据权利要求1所述的基于可视化算法的电抗器温升计算方法,其特征是,所述步骤5)中的变化趋势图根据步骤4)中的计算结果以及以电抗器铁芯磁路长度l、电抗器铁芯截面积S、电抗器线圈匝数N为变量绘出。
4.根据权利要求3所述的基于可视化算法的电抗器温升计算方法,其特征是,所述变化趋势图的约束条件是:1)所述电抗器铁芯磁路长度l的范围约束条件以硅钢EI型材料最小和最大尺寸为依据,约束条件为0.001m~0.5m;2)所述电抗器线圈匝数N的约束条件为1~100;3)所述电抗器铁芯截面积S的约束条 件为0.000001m2~0.005m2。
5.根据权利要求1所述的基于可视化算法的电抗器温升计算方法,其特征是,所述步骤5)中获得电抗器温升平面图的过程为,在满足步骤4)中的限制条件下,依据选择的电抗器铁芯磁路长度l作为定量,在电抗器温升变化趋势图中找出相关满足域,即可得到变量参数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410494881.7A CN104318060A (zh) | 2014-09-24 | 2014-09-24 | 基于可视化算法的电抗器温升计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410494881.7A CN104318060A (zh) | 2014-09-24 | 2014-09-24 | 基于可视化算法的电抗器温升计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104318060A true CN104318060A (zh) | 2015-01-28 |
Family
ID=52373291
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410494881.7A Pending CN104318060A (zh) | 2014-09-24 | 2014-09-24 | 基于可视化算法的电抗器温升计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104318060A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112800569A (zh) * | 2019-11-13 | 2021-05-14 | 天津经纬正能电气设备有限公司 | 干式空心电抗器的优化方法 |
CN116136557A (zh) * | 2023-04-20 | 2023-05-19 | 北京电科能创技术有限公司 | 柔性直流输电用桥臂电抗器温升测量方法、设备及介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11157946A (ja) * | 1997-12-01 | 1999-06-15 | Nippon Cement Co Ltd | セラミックス製造過程のシミュレーション方法、およびそれを用いたセラミックス部材の設計方法 |
CN103956259A (zh) * | 2014-04-18 | 2014-07-30 | 江苏南瑞帕威尔电气有限公司 | 一种基于可视化算法的电抗器设计方法 |
CN104036125A (zh) * | 2014-05-30 | 2014-09-10 | 河北省电力建设调整试验所 | 一种油浸式变压器内部温度场的精确计算方法 |
-
2014
- 2014-09-24 CN CN201410494881.7A patent/CN104318060A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11157946A (ja) * | 1997-12-01 | 1999-06-15 | Nippon Cement Co Ltd | セラミックス製造過程のシミュレーション方法、およびそれを用いたセラミックス部材の設計方法 |
CN103956259A (zh) * | 2014-04-18 | 2014-07-30 | 江苏南瑞帕威尔电气有限公司 | 一种基于可视化算法的电抗器设计方法 |
CN104036125A (zh) * | 2014-05-30 | 2014-09-10 | 河北省电力建设调整试验所 | 一种油浸式变压器内部温度场的精确计算方法 |
Non-Patent Citations (2)
Title |
---|
刘刚: "磁阀式可控电抗器的研究及计算机辅助设计", 《中国优秀硕士学位论文全文数据库工程科技II辑》 * |
刘旭光等: "直流铁心电感器优化设计", 《低压电器》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112800569A (zh) * | 2019-11-13 | 2021-05-14 | 天津经纬正能电气设备有限公司 | 干式空心电抗器的优化方法 |
CN116136557A (zh) * | 2023-04-20 | 2023-05-19 | 北京电科能创技术有限公司 | 柔性直流输电用桥臂电抗器温升测量方法、设备及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106250626B (zh) | 一种基于l-i曲线的特高压变压器空载直流偏磁快速计算方法 | |
CN103956259B (zh) | 一种基于可视化算法的电抗器设计方法 | |
Muhlethaler et al. | Core losses under the DC bias condition based on Steinmetz parameters | |
Kovačević et al. | 3-D electromagnetic modeling of EMI input filters | |
Njafi et al. | Evaluating and derating of three-phase distribution transformer under unbalanced voltage and unbalance load using finite element method | |
CN102592062A (zh) | 一种变压器直流偏磁动态漏电感计算方法 | |
Burkart et al. | Optimal inductor design for 3-phase voltage-source PWM converters considering different magnetic materials and a wide switching frequency range | |
Zimmanck et al. | Efficient calculation of winding-loss resistance matrices for magnetic components | |
CN111914413A (zh) | 一种对称/不对称矩形电压激励下磁芯高频损耗计算方法 | |
Kaymak et al. | Comparison of analytical methods for calculating the AC resistance and leakage inductance of medium-frequency transformers | |
Delgado et al. | Equivalent conductor layer for fast 3-D finite element simulations of inductive power transfer coils | |
CN104318060A (zh) | 基于可视化算法的电抗器温升计算方法 | |
Reinlein et al. | Optimization of analytical iron loss approaches for electrical machines | |
CN106546811A (zh) | 一种谐波电流下变压器负载损耗的检测方法及系统 | |
Su et al. | Estimation of the core losses of induction machines based on the corrected Epstein test method | |
He et al. | An improvement of core losses estimation model in power electronic transformer | |
Kefalas et al. | Reduction of cost and losses of transformers by using composite magnetic cores | |
Ranta et al. | Inclusion of hysteresis and eddy current losses in nonlinear time-domain inductance models | |
Sato et al. | Study on an accurate iron loss calculation method considering the non-uniformity of the magnetic flux density | |
CN103630746A (zh) | 一种测量和评估有气隙磁元件绕组交流电阻的方法 | |
Luo et al. | 2-D winding losses calculation for round conductor coil | |
Beyranvand et al. | Finite element study of ferroresonance in single-phase transformers considering magnetic hysteresis | |
Zhao et al. | Harmonic analysis of the DC biased Epstein frame-like core model by the harmonic balance finite element method | |
Susnjic et al. | Stray losses computation in power transformer | |
Choudhary | 3D modelling and finite element analysis of three phase transformer with integrated inductors in COMSOL Multiphysics |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150128 |