电力变压器铁芯柱截面优化设计

1、第一页答卷编号:论文题目：电力变压器铁心柱截面的优化设计姓名专业、班级有效联系电话参赛队员1XXXXX电气工程及其自动化2班XXXXXXXXX参赛队员2XXXXXXXXXXX电气工程及其自动化2班XXXXXXXXXXXXXX参赛队员3XXXXXXXXXX数学与应用数学XXXXXXXXXXXX扌旨导教师： XXXXXXXXX参赛学校：中国地质大学（北京）报名序号：XXXXXXXXXXXXX证书邮寄地址：北京市海淀区学院路29号中国地质大学（北京）16楼504宿舍（学校统一组织的请填写负责人）第二页答卷编号:阅卷专家1阅卷专家2阅卷专家3论文等级内容摘要以电力变压器铁心柱截面在圆形的线圈筒里面为例

2、，建立非线性整数规划数学模型，然后用 lingo编写程序求出级数和各级宽度和厚度的最优解，以提高铁心柱有效截面积。当铁心柱截面外接圆直径为650毫米时，级数为14.加入的油道应该位于第一.二级和第四.五级之间。关键词有效截面积非线性整数规划假设法lingo 编程百分误差1问题简述为了充分利用线圈内空间又便于生产管理，心式铁心柱截面常采用多级阶梯形结构，如 图1所示。截面在圆内上下轴对称，左右也轴对称。阶梯形的每级都是由许多同种宽度的硅 钢片迭起来的。由于制造工艺的要求，硅钢片的宽度一般取为5的倍数（单位：毫米）。因为在多级阶梯形和线圈之间需要加入一定的撑条来起到固定的作用，所以一般要求第一级的

3、厚铁心柱有效截面的面积，等于多级铁心柱的几何截面积（不包括油道）乘以叠片系数。设计时希望有效截面尽量大。但由于制造工艺限制，所以直径在某范围是，一般情况下铁心柱的级数可参 照表1选取。表1铁心柱截面级数的选择铁心柱直径mm级数80-1955-7200-2658-10270-39011400-74012-14760以上>15问题一：当铁心柱外接圆直径为650毫米时，如何确定铁心柱截面的级数、各级宽度和厚度, 才能使铁心柱的有效截面积最大。问题二：实际生产中线圈的内筒直径和铁心柱的外接圆直径不是精确地相等，而留有一定的间隙以便于安装和维修，设计的两个直径的取值范围称为各自的公差带。 因此可以

4、在设计铁心截面 时稍微增加铁心柱的外接圆的直径以使得铁心柱有更好的截面形状。 结合铁心柱截面的设计而设计 出二者的公差带。问题三：铜导线在电流流过时发热造成的功率损耗简称为铜损；铁心在磁力线通过时发热造 成的功率损耗简称为铁损。为了改善铁心内部的散热，铁心柱直径为380毫米以上时须设置冷却油 道。简单地说，就是在某些相邻阶梯形之间留下 6毫米厚的水平空隙（如图2所示），空隙里充满 油，变压器工作时油上下循环带走铁心里的热量。具体油道数可按表2选取。油道的位置应使其分 割的相邻两部分铁心柱截面积近似相等。分别针对问题一和问题二的情况，增加油道要求再设计，并指出油道的位置。油道图2带油道的铁心柱截

5、面表2冷却油道数的选择铁心柱直径mm半圆中6mm油道个数380-4100420-5001510-6902700-84032. 模型假设1 叠片系数通常与硅钢片厚度、表面的绝缘漆膜厚度、硅钢片的平整度以及压紧程度有关。在这里，视为定值，取0.97.2 有无油道对铁心柱截面级数无影响。3加入油道后，铁心主界面各级在圆形内的位置变化较小，用于初步确定油道大概位置。3. 符号约定由于电力变压器心式铁心柱截面是上下对称的，所以取上半部分考虑S铁心柱总截面积Si铁心柱第i级截面积i =1,2,.,n铁心柱的叠片系数（一般在0.97左右）Xi第i级宽度i =1,2,.,nyi第i级的厚度i=2，,n （i

6、=1表示第一级厚度的一半）尬 尾数，为5或10d外接圆直径r 外接圆半径4问题分析与模型建立:由于叠片系数一般情况下，取 0.97，为定值，所以要求有效面积最大，只要求总的截面面积 最大。根据题中的约束条件，分别在考虑有道和不考虑有道的情况下建立数学模型，用lingo编程求解。4.1模型一：不考虑有油道情况：由于制造工艺的要求，硅钢片的宽度一般取为 5或10的倍数（单位：毫米）。因为在多级阶梯 形和线圈之间需要加入一定的撑条来起到固定的作用，所以一般要求第一级的厚度最小为 26毫米, 硅钢片的宽度最小为20毫米。对直径4，以（xx2,，Xn）为其各级的宽度， （yy2,，yn ）为1i_1相应

7、的厚度。面积函数为Si =2Xiyi，x与y的关系利用勾股定理：y： = （r2 - （x/2）竽 八 i =1,2, n ,可得以下模型每一级的有效面积：S二diyii T则目标函数：max S=2'，Si经约束条件:各级宽度必须为 的倍数y1 -13第一级厚度必须大于26mm650-xn -20第n级宽度大于20mm各级宽度逐级递减Xr 丄 X2 _- xn此模型为非线性整数规划。4.1.1模型分析与求解求解:显然铁心柱的级数愈多，其截面愈接近于圆形，在一定的直径下铁心柱有效截面也愈大。但这样制造也工艺复杂，一般情况下铁心柱的级数可参照上面的表1选取。可知：i =12.13或14.

8、利用:lingo软件编写求解上述线性规划问题得到：当i =12时，截面面积最大值为：318927.2平方毫米当i =13时，截面面积最大值为：319923.9平方毫米当i =14时，截面面积最大值为：320753.1平方毫米又1 一般取0.97，则有效面积最大值为311133.999平方毫米。所以得到各级宽度和厚度如下表：表2.铁心柱截面设计最优解级数宽度（单位：毫米）厚度（单位：毫米）164056.78908262040.80702359533.24253457025.36635554024.69553651020.59389747520.35858843519.63976939516.61

9、3331035015.755201130014.452791224512.715421318011.26045141058.4416564.1.2公差带的确定:公差带是由标准公差和基本偏差两个基本要素确定的，标准公差确定公差带的大小；基本偏差确定公差带相对于零线的位置。1）标准公差：标准公差是由国家标准规定的，用于确定公差带大小的任一公差。公差等级确定尺寸的精确程度，国家标准把公差等组分为20个等级，分别用IT01、IT0、IT1IT18表示，称为标准公差，IT （InternationalToleranee ）表示标准公差。当基本尺寸一定时，公差等级愈高，标准公差值愈小，尺寸的精确度就愈高。

10、基本尺寸和公差等级相同的孔与轴，它们的标准公差相等。2）基本偏差；国家标排规定用来确定公差带相对于零线位置的上偏差或下偏差；一般为最靠 近零线的那个偏差为基水偏差。当公差带位于零线的上方时，基本偏差为下偏差；当公差带位于零 钱的下方时，基本偏差为上偏差。所以：参照附录三中的国际标准公差表，选取IT15等级进行计算。用.=3.2表示公差，则:对于铁心柱：d -： d ： d（650-3.2）毫米：d : （650+3.2）毫米。铁心柱的公差带为：（646.8，653.2 ）毫米。对于线圈的内筒直径：D =d 1662毫米.同样选取IT15等级进行计算。用=3.2表示线圈内筒直径。则根据附录四的表

11、格：Dmin 二d 12 -' =655.6毫米Dmax =d 12、=668.4毫米所以线圈内筒公差带为：（655.6,668.4）毫米。4.2模型二：考虑有油道的情况：根据题目意思，有油道情况下，既要满足所求各级宽度和厚度是最优解，又要满足油道分割相邻两部分面积近似。参见下表：表3冷却油道数的选择铁心柱直径mm半圆中6mm油道个数380-4100420-5001510-6902700-8403由此表，可知d=650毫米时，半圆中应选择两个油道，整个圆被分成五个部分.设S为油道平均分割各铁心柱截面各部分的的面积，且 Si ,S2 ,S3 ,S4 , S之间近似相等，即 允许与绝对分割

12、面积间存在误差。 先令此误差取较大的一个数 。 的取值不同，近似程度也就不 同，从而就会有不同的优化结果。由于是近似，所以可认为整个铁心柱被四个油道分割的五块面积3 ,S2 , S3 ,S4 , S5均近铁心柱在模型一种求得的面积S的1/5为64150.62平方毫米。下面以=10000平方毫米为例。此时以下求出各部分，百分误差最大为：（5 /S） 100%=15.883%.有：只要求出的优化结果中，S, ,S2 , S3 ,S4 , S5的值在（S/5） 一即可认为符合近似条件，并且满足截面总面积S最大即可。4.2.1运用假设法进行模型分析：不妨设半圆中，这两个油道分别在第m级与第m 1级之间

13、和第j级与第j 1级之间（其中m v j ）。下面求m和j即可确定油道位置。假设一：油道直接加入在模型一中的某两极之间，即加入油道后，模型一求出的优化结果仍然为此时的最优解！则由于：2S =72690.0224 平方毫米.S2 S3 S4 =59538.47725 平方毫米S5 &S14 =64493.06155 平方毫米所以：m=1, j=4.两个油道应该分别放在第一级和第二级之间，第四级和第五级之间。此时只要把第二级和第五级的厚度分别减去油道厚度后各级宽度和厚度即为放入油道后各级宽度和厚度!51 =20 =72690.0224 平方毫米52 -S4 =S2 S3 S4- 6X2=5

14、5818.47725 平方毫米FF53 = S5 =S5 S6. S14 -6x5 =61253.06155 平方毫米截面总面积为：S = S1S2S3S4S5 =306833.1平方毫米表4.假设一中铁心柱截面设计最优解级数宽度(单位：毫米)厚度(单位：毫米)164056.78908262034.80702359533.24253457025.36635554018.69553651020.59389747520.35858843519.63976939516.613331035015.755201130014.452791224512.715421318011.26045141058.44

15、1656假设二：加入油道后，模型一中求得的优化结果已经不是有油道的最优解了，那么要根据条件，重新建立模型求解：由假设一可基本确定油道一放在第一级和第二级之间，只不过此时得出的最优解的第一级和第二级的宽度和厚度不一定与假设一的结果相同。所以可以得到m的值为1.下面根据条件有：确定j和各级宽度和厚度。各级宽度和厚度根据勾股定理有如下关系：i =1 时：1 i 4yi =(r2 -(x/2)2)2ytt=1i =2，j 时：1 i 4yi =(r2 -(x/2)2)2yt -6t *i = j 1, n1 i Ayi =(-(x/2)2)2yt-12每一级的面积：Si =xi yii仝则目标函数：m

16、ax S= 2'-$i 4约束条件：s.t.各级宽度必须为的倍数yi -13第一级厚度必须大于26mm650_xn -20第n级宽度大于20mmS/5 -瓷 S : S/5确保Si ,S2 , S3符合近似条件各级宽度逐级递减S/5工 S2 : S/5S/5 -工 Ss :： S/5Xi 亠 X2 - xn4.2.2模型的求解：利用lingo软件编程求解，程序编完后，分别令j =2,3,41516.14，得到如下最优解:j =4时，S最大，此时，S =306943.4平方毫米F51 =2=72690.0224 平方毫米.FF52 二S4二S2S3S4-6x2=62680.66305 平

17、 方 毫FF53 二 S5 =S5 S6. S14 -6X5 -54446.02575 平方毫米则各级优化设计如下表：表5.假设二中铁心柱截面设计最优解级数宽度（单位：毫米）厚度（单位：毫米）164056.78908261542.41904358536.45609455527.50656552019.92923649018.54157745518.55486842015.94591938015.633601033514.835731128513.5820712230131.879951317.9.714067141007.443065综上可知，假设二中截面总面积：S = S S2 "

18、 S3 S4S5 =306943.4平方毫米大于假设一中截面总面积：S = s/ S2S3 s/ S5 =306833.1平方毫米，所以假设二更合理，油道分别位于第一二级和第四五级之间，有油道后的优化设计应采取假设二的优化结果。此时，公差带与模型一中一样，分别为：铁心柱的公差带为：（646.8，653.2 ）毫米。线圈内筒公差带为：（655.6, 668.4） 毫米。5. 模型评价及改进：1. 在考虑油道的情况下，采用了假设的方法，如若不采用此方法，也可以部分两种情况假设， 直接建立模型求解。2. 模型二中，假设二的前提是加入油道后，级数不变。如果级数变化，则应该重新建立模型求解。3. 可以在

19、分析有油道的情况下，由于是油道分割相邻两部分面积近似相等， 所以我们就近似的 取模型一中，总有效截面积的1/5作为基础数，把加入油道后各部分面积与之比较， 来确定百分误 差。其实可以定义一个较科学的相似度的概念，用来衡量二者的接近程度。6. 参考文献：1. 潘淑平王晓玲张秀兰电力变压器铁心柱截面的优化设计2006年3月2. 陈晔 何荣坚 郑可逵 电力变压器铁心柱截面的优化设计2004年12月3. 沈变力电压器设计手册铁心设计4. 国际标准公差表附录：一.无油道时，用Lingo软件编写求解最优的代码：（其中n分别为12.13.14时，运行求最优解。给出的是n=14的情况。）model :data

20、 :d,n=650,14;en ddatasets :n umber/1,2. n/:x,y,w;en dsets!目标函数；max=2* sur(number(j):x(j)*y(j);!约束条件；for (number(i):d>x(i);x14>20;y1>13;for (number(i):sum(number(i):y(i)<d/2);for (number(i):x(i)=5*w(i);for (number(i):dA2/4=x(i)A2/4+(sum(number(j)|j#le#i:y(j)A2);for (number(i)|i#le#n-1:x(i

21、)>=x(i+1)+5);for (number(i):gin (w(i);一.有油道时，用Lingo软件编写求解假设二中最优解的代码：(其中j分别为2.34514时，运行求最优解。给出的是j =4的情况。)model :data :d=650; n=14;j=4;p=64150.62;en ddatasets :S/1. n/:a,b,c;en dsets!目标函数；max =2* sum(S(l):a(l)*b(l);!约束条件；for (S(l):a(l)=5*c(l);for (S(l):(a(l)/2)A2+(if (l#le#j,6,12)+if (l#le#1,-6,0)+

22、sum(S(k)|k#le#l:b(k)A2=D*D/4);!满足近似条件；abs( sum(S(l)|l#le#j:a(l)*b(l)+p-sum(S(l):a(l)*b(l)<10000;abs(p+a(1)*b(1)-sum(S(l)|l#le#j:a(l)*b(l)<10000;abs (2*a1*b1-p)<10000;for (S(k)|k#le#(n-1):a(k+1)+5<=a(k);a(1)>20;b(1)>13;for (S(l):gin (c(l);.根据国际标准，以下为基本尺寸 0-1000mm 4-18级精度标准公差表基本尺寸公差值IT4lT5lT6lT7lT8lT9lT10lT11lT12lT13lT14lT15lT16lT17lT18大于到mm-331 4610142540600.100.140.250.400.601.01.43645812183048750.120.180.300.480.751.21.861046915223658900.150.220.360.580.901.52.2101858111827431 701100.180.270.430.701.101.82.