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

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-39011精品文档400-74012-14760以上>15问题一：当铁心柱外接圆直径为650毫米时,如何确定铁心柱截面的级数、各级宽度和厚度, 才能使铁心柱的有效截面积最大.问题二：实际生产中线圈的内筒直径和铁心柱的外接圆直径不是精确地相等,而留有一定的 间隙以便于安装和维修,设计的两个直径的取值范围称为各自的公差带. 因此可以

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

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

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

7、Xiyi,x与y的关系利用勾股定理： y： = r2 - x/2于 八 tvi =1,2, n ,可得以下模型每一级的有效面积：& = Xyii =n那么目标函数：max S=2 Si A约束条件：s.t.生二生各级宽度必须为的倍数上coy1 -13第一级厚度必须大于26mm650- xn -20第n级宽度大于20mmX1 一 X2 - Xn各级宽度逐级递减此模型为非线性整数规划.精品文档精品文档4.1.1模型分析与求解求解:显然铁心柱的级数愈多,其截面愈接近于圆形,在一定的直径下铁心柱有效截面也愈大.但这样制造也工艺复杂,一般情况下铁心柱的级数可参照上面的表1选取.可知：i =12.

8、13或14.利用: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

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

10、高.根本尺寸和公差等级相同的孔与轴,它们的标准公差相等.2根本偏差；国家标排规定用来确定公差带相对于零线位置的上偏差或下偏差；一般为最靠 近零线的那个偏差为基水偏差.当公差带位于零线的上方时,根本偏差为下偏差；当公差带位于零 钱的下方时,根本偏差为上偏差.所以：参照附录三中的国际标准公差表,选取IT15等级进行计算.用.=3.2表示公差,那么:对于铁心柱：(650-3.2)毫米：d :：: (650+3.2)毫米.铁心柱的公差带为：(646.8 , 653.2 )毫米. 对于线圈的内筒直径：D =d12 =662毫米.同样选取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为油道平均分割各铁心柱截面各局部的的面积,且 3占,S3 Q , S之间近似相等,即 允许与绝对分割面积间存在

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

13、 1级之间(其中m < j ).下面求m和j即可确定油道位置.精品文档假设一：油道直接参加在模型一中的某两极之间,即参加油道后,模型一求出的优化结果仍然为此时的最优解！那么由于：23 =72690.0224 平方毫米.S2S3 S4 =59538.47725 平方毫米Ss S6. 务=64493.06155平方毫米所以：m=1, j =4.两个油道应该分别放在第一级和第二级之间,第四级和第五级之间.此时只要把第二级和第五级的厚度分别减去油道厚度后各级宽度和厚度即为放入油道后各级 宽度和厚度！rS = 2Si =72690.0224平方毫米S2 = S4 =S2 S3 S4- 6x2 =5

14、5818.47725 平方毫米Ss 二 s5 =S5 S6 Sm -6x5 =61253.06155 平方毫米rrrrr截面总面积为：S = s S2 S3S4S5 =306833.1平方毫米表4.假设一中铁心柱截面设计最优解级数宽度单位：毫米厚度单位：毫米164056.78908262034.80702359533.24253457025.36635554018.69553651020.59389747520.35858843519.63976939516.613331035015.755201130014.452791224512.715421318011.26045141058.4416

15、56假设二:精品文档精品文档参加油道后,模型一中求得的优化结果已经不是有油道的最优解了,那么要根据条件,重新建立模型求解：由假设一可根本确定油道一放在第一级和第二级之间,只不过此时得出的最优解的第一级和第二级的宽度和厚度不一定与假设一的结果相同.所以可以得到m的值为1.下面根据条件有：确定j和各级宽度和厚度.各级宽度和厚度根据勾股定理有如下关系：i =1 时：1i J.yi =(r -(x/2) )2 -、yt12i =2, j 时：y = (r2 -(x/2)2)2yt -6t ii 二 j 1, ,nyi =(r -(x/2) )2yt-12每一级的面积：Si =Xi yii那么目标函数：

16、max S= 2 Sii 二约束条件：s.t.丸二彳各级宽度必须为的倍数曲coY1 -13第一级厚度必须大于26mm650-xn -20第n级宽度大于20mmS/5 -r:3:S/5S/5 -:：S2 :：S/5保证S1 ,S2 , S3符合近似条件S/5 -r:S3 :：S/5x-x2亠Xn各级宽度逐级递减4.2.2模型的求解:利用lingo软件编程求解,程序编完后,分别令j =2,3,4,5,6.14,得到如下最优解：j =4时,S最大,此时,S =306943.4平方毫米S =2=72690.0224平方毫米.FF52 二S4 二S2 S3 S4- 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 S3 S4 -

18、Ss =306943.4平方毫米大于假设一中截面总面积：S = S S2 S3 S4 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. 国际标准公差表附录:n分别为12.13.14时,运行求最无油道时,用Lingo软件编写求解最优的代码：其中 优解.给出的是n=14的情况.model :data :d

20、,n=650,14;en ddatasets :n umber/1,2. n/:x,y,w;en dsets!目标函数；max=2* sur(number(i):x(i)*y(i);!约束条件；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)&g

21、t;=x(i+1)+5);for (number(i):gin (w(i); 有油道时,用Lingo软件编写求解假设二中最优解的代码：(其中j分别为2.3.4.514时,运行求最优解.给出的是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)/2f2+(if (l#le#j,6,12)+if (l#le#1,-

22、6,0)+sum(S(k)|k#le#I: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(I)|I#le#j:a(I)*b(I)<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级精度标准公差表根本尺寸公差值IT4IT5IT6IT7IT8IT9IT10IT11IT12IT13IT14IT15IT16IT17IT18大于到mm-334610142540600.100.140.250.400.601.01.43645812183048750.120.180.300.480.751.21.861046915223658900.150.220.360.580.905 J2.210185811182743701100.180.270.430.701.101.82.71