- 相關推薦
電力變壓器鐵心柱截面的優化設計(一)
摘要
本文針對變壓器鐵心柱截面優化設計，建立數學模型并用matlab函數constr編程以及使用LINGO編程求解。
最優鐵心柱截面設計方案指的是在達到使用標準的情況下，讓鐵心柱橫截面
的有效截面面積最大，從而提高產品的使用性能，提高產品的競爭力。
問題一中，鐵心柱外接圓直徑為650Mm，根據截面級數的選擇表，可以確定其級數的范圍。根據制造工藝要求的第一級厚度最小26和最后一級寬度要求最小為20，各級之間寬度依次減少為約束，以有效截面積最大為目標，建立起含多個整數變量的非線性整數規劃模型——模型一，優化求解以確定最優的多級硅鋼片階梯形組成方案。
利用MATLAB軟件規劃求解得到最優的多級階梯形硅鋼片組成方案為：
鐵心柱有效截面積320240，鐵心利用率96.50%，級數為13，
各級的寬度和厚度為：
（640，114）（620，82）（595，66）（570，51）（535，57）（500，46）（460，44）（415，41）（370，34）（315，34）（260，27）（190，26）（110，19）
對于問題二和問題三，我們將第一問的模型進行拓展和延伸。
第二問要求線圈內筒和鐵心柱外接圓直徑的公差。以有效截面積利用率最大為目標，以變壓器設計結構設計要求為約束，我們再次建立了一個非線性整數規劃模型。
在考慮添加油道的問題中，需要確定油道位置，并且分割的各個面積要求基本相等。我們首先確定油道的位置，然后在分割出來的各個小塊里，分別去進行多級階梯形的最優化設計。各個分塊里實現了最優化，整個設計就實現了最優化。
該過程中，油道和各個分塊之間會產生間隙，針對此問題，我們將油道的位置向圓心靠近調整，然后新產生的間隙，用下一個分塊的第一級硅鋼片進行補充。
由于鐵心柱疊片系數是一個常數，在模型求解的過程中，可以忽略其影響。
在建模的過程中，我們運用了MATLAB和LINGO軟件進行求解。
關鍵詞:非線整數性規劃，鐵心利用率，分布填充
問題重述
電力變壓器的設計中很重要的一個環節就是鐵心柱的截面如何設計。我國變壓器制造業通常采用全國統一的標準鐵心設計圖紙。根據多年的生產經驗，各生產廠存在著對已有設計方案的疑問：能否改進及如何改進這些設計，才能在提高使用效果的同時降低變壓器的成本。
現在以心式鐵心柱為例試圖進行優化設計。
電力變壓器鐵心柱截面在圓形的線圈筒里面。為了充分利用線圈內空間又便于生產管理，心式鐵心柱截面常采用多級階梯形結構，如圖1所示。截面在圓內上下軸對稱，左右也軸對稱。階梯形的每級都是由許多同種寬度的硅鋼片迭起來的。由于制造工藝的要求，硅鋼片的寬度一般取為5的倍數（單位：毫米）。因為在多級階梯形和線圈之間需要加入一定的撐條來起到固定的作用，所以一般要求第一級的厚度最小為26毫米，硅鋼片的寬度最小為20毫米。
鐵心柱有效截面的面積，等于多級鐵心柱的幾何截面積（不包括油道）乘以疊片系數。而疊片系數通常與硅鋼片厚度、表面的絕緣漆膜厚度、硅鋼片的平整度以及壓緊程度有關。設計時希望有效截面盡量大，既節省材料又減少能量損耗。顯然鐵心柱的級數愈多，其截面愈接近于圓形，在一定的直徑下鐵心柱有效截面也愈大。但這樣制造也工藝復雜，一般情況下鐵心柱的級數可參照表1選取。
表1 鐵心柱截面級數的選擇
鐵心柱直徑mm 級數
80-195 5-7
200-265 8-10
270-390 11
400-740 12-14
760以上 >15
問題一：當鐵心柱外接圓直徑為650毫米時，如何確定鐵心柱截面的級數、各級寬度和厚度，才能使鐵心柱的有效截面積最大。
問題二：實際生產中線圈的內筒直徑和鐵心柱的外接圓直徑不是精確地相等，而留有一定的間隙以便于安裝和維修，設計的兩個直徑的取值范圍稱為各自的公差帶。因此可以在設計鐵心截面時稍微增加鐵心柱的外接圓的直徑以使得鐵心柱有更好的截面形狀。請結合鐵心柱截面的設計而設計出二者的公差帶。
問題三：銅導線在電流流過時發熱造成的功率損耗簡稱為銅損；鐵心在磁力線通過時發熱造成的功率損耗簡稱為鐵損。為了改善鐵心內部的散熱，鐵心柱直徑為380毫米以上時須設置冷卻油道。簡單地說，就是在某些相鄰階梯形之間留下6毫米厚的水平空隙（如圖2所示），空隙里充滿油，變壓器工作時油上下循環帶走鐵心里的熱量。具體油道數可按表2選取。油道的位置應使其分割的相鄰兩部分鐵心柱截面積近似相等。
分別針對問題一和問題二的情況，增加油道要求再給出設計，并指出油道的位置。
表2 冷卻油道數的選擇
鐵心柱直徑mm 半圓中6mm油道個數
380-410 0
420-500 1
510-690 2
700-840 3
2 條件假設
1.整個鐵心柱的硅鋼片出了長度，其它如厚度，表面絕緣漆膜厚度，平整度都相同；
2.硅鋼片之間是沒有形變的壓緊；
3.疊片系數是確定的已知數；
4.油道對稱分布；
5.不考慮工藝過程的影響。
3 符號說明
第i級疊片的的厚度；
第i級疊片的寬度；
疊片系數；
鐵心柱理論外接圓的直徑，也等于理論線圈內筒直徑；
第i根油道與直徑之間的；
油道分割出來的分塊面積；
多級階梯形前i級厚度之和；
線圈內筒的公差值；
4 問題分析
4. 問題背景
變壓器是一種應用電磁感應原理把電能從一個電路傳到另一個電路的電磁裝置。它在電路中起變壓，變流，變電阻的作用。它由三部分組成：鐵心，起導磁，助磁作用；而是初級線圈，接電源，起激磁作用；三是次級線圈接負載，利用不同次初級線圈匝數比，實現變壓，變流，變電阻的作用。
在變壓器的構成里，鐵心柱是很一個十分重要的組成部件。
因為鐵心柱的形狀，截面積，疊片的選擇，疊片的相關工藝過程都會影響將來變壓器的使用效果和壽命，以及使用成本。
我國變壓器制造業通常采用全國統一的標準鐵心設計圖紙，根據多年的生產經驗，在長期的生產研究過程中，各生產廠產生了對已有設計方案的疑問：能否改進及如何改進這些設計，才能在提高使用效益的同時降低變壓器成本。
4.1.
鐵心柱是安裝在線圈筒里面的，理論上，在線圈的直徑確定了的情況下，鐵心柱的有效面積越大，鐵心的電阻越大，使鐵心的鐵損最小，因而可以減少能量損耗，變壓器的使用性能會越好，使用壽命較長。為了充分利用空間和便于生產，鐵心柱截面長采用多級階梯形結構，用不同長度的硅鋼片，堆疊成不同厚度的級，并且選擇合適的級數去逼近與之配合的鐵心線圈圓，期望得到的有效面積最大，獲得更大的電阻。
截面優化設計是以保證到達設計標準為前提，盡可能改善和提高產品的使用效果，使產品競爭力提升。
4.2
公差是生產中允許工件尺寸和幾何形狀變動的范圍，用來限制誤差。工件的誤差在公差范圍內才為合格。規定公差是為了保證產品使用性能的前提下，給出盡可能大的公差。因為公差越小，精度越高，生產成本就越高。去較大的公差，有利于控制成本。
設計時要求留有一定的間隙邊緣和維修，因此鐵心柱和外圍線圈的配合應該是間隙配合。
我們根據《變壓器設計手冊》查詢得到的GB的線圈內筒和鐵心柱的直徑差
公差
直徑（mm） 70-200 205-255 260-500 500-650 650-800
公差 (mm) 5 8 10 12 14
理論數值的鐵心柱外接圓和線圈內筒直徑，它們在配合的時候，采用間隙配合。由于生產硅鋼片的時候，是以滿足鐵心柱的設計尺寸為參照的，用來裝配鐵心柱的時候，它們的寬度是一系列已經確定的數值，不可能因為包裹鐵心柱的線圈內筒的直徑的增大而增大。
即當線圈內筒的直徑變化的時候，多級階梯形的各級的寬度是不可能變化的。為了獲得更好的截面積，或者增大截面積的值，我們可以將多級階梯形的某一級的疊片數增加，從而增加該級厚度，去逼近實際尺寸的線圈內筒。
增大線圈內筒的直徑，再增大某級的厚度，獲得更大的有效截面積。我們認定有讓效截面積增加率最大的那個偏差，就是鐵心柱公差的最大值。
據此選擇合適的公差帶，來獲得最大的有效面積增加率。
4.3.1
由于變壓器工作時鐵心柱有鐵心損耗，外圍線圈有電阻，電流通過的時候會產生熱量。為了改善鐵心散熱，在鐵心柱的直徑達到380mm以上，就需要增加油道，讓油循環達到散熱的目的。
在考慮添加油道的時候，需要讓油道分割的各部分的面積近似相等。為了達到這個要求，我們以各個分割面積相等為約束，對稱布置油道，建立起方程組，用以先確定下油道的具體位置。
確定油道位置之后，再根據鐵心柱有效截面積最大的要求，可以分步用硅鋼片將各個油道之間的分塊堆疊。在各個分區里面，在滿足各項設計要求的情況下，求堆疊的最優方案，使得各個分區的截面積最大。最后將各個面積進行疊加，各個分區的硅鋼片進行參數匯總，就可得最終的最優方案。
4.3.2
對于變壓器鐵心柱的級寬，都可以通過已知的直徑算出相應的厚度，因此可以轉化為以各級的寬度為變量的數學模型。由于現有的各直徑鐵心柱的級數都是參照生產經驗，認為要達到一定的級數才能滿足一定的幾何截面積，而沒有科學的精確的求解過程。因此設想通過對模型的計算所得出的數據，只要在滿足一定的幾何截面積的前提下，如果能減少級數，便可減少生產成本。
5 模型建立求解
5.1.1
問題一的數學模型：
當鐵心柱的外接圓直徑為650mm時，根據要求可知級數應該在12到14之間。由變壓器鐵心設計的各項要求為約束條件，我們以鐵心柱有效截面積最大為目標函數，建立了一個非線性整數規劃的模型。
約束條件
硅鋼片的寬度為5mm的倍數：
（ k為正整數）
多級階梯形的第一級厚度最小為26mm:
硅鋼片的寬度最小為20mm，就是最后一級的寬度最小為20mm
最大寬度為直徑：
多級階梯形中，每級的寬度必須是遞減的：
（）
I級及i級之前的厚度之和與i級的寬度之間還必須滿足勾股定理：
該多級階梯形的幾何截面積為：
其中（n=12，13，14）
目標函數：
我們認為所有硅鋼片的疊片系數是相同的，即是確定的常數，
求有效截面積最大：
綜上可知：
為確定鐵心柱直徑為650mm時候的最佳設計方案，以有效截面積最大為目標，在滿足各項結構設計要求的約束下，我們建立了模型一：
目標函數：
5.1.2 模型一的求解：
將以上模型代入編寫的matlab程序中予以求解（程序見附件），得到如下最優化
表 I 直徑650mm方案對比
硅鋼片參數 1 2 3 4 5 6 7 8 9 10 11
12 13 14 面積 面積使用率(%)
寬度（mm） 640 620 590 560 520 480 435 385 330 270 200 115 319405 96.26
長度（mm） 114 82 78 57 60 48 45 41 36 31 27 21
寬度（mm） 640 620 595 570 535 500 460 415 370 315 260 190 110 320240 96.51
長度（mm） 114 82 66 51 57 46 44 41 34 34 27 26 19
寬度（mm） 640 625 600 575 545 515 480 440 400 350 300 245 180 105 320230 96.50
長度（mm） 114 65 71 53 51 42 42 40 34 35 29 25 23 17
對比可知，所要求的鐵心柱有效截面積最大的最優方案為14級，各級寬度和長度可由上表可得
5.2.1問題二模型
公差等級是確定尺寸精度的等級，國家標準是按照不同的基本尺寸分為20個等級。由于電力變壓器不像機械設備那樣需要承受大的負荷，比如說軸承那樣的高速旋轉，像發動機那樣的劇烈復雜的機械運動，承受如此復雜的力，振動，高溫高壓等復雜嚴苛的工作狀況，我們認為電力變壓器的公差等級為18，已經足夠滿足它的使用了。
根據變壓器設計即為合格產品的最大范圍區間，然后在這個范圍內根據最優化截面的思想，即求最大的有效截面增加率的所對應的es值，就是它的公差最大值。
考慮鐵心柱和線圈內筒的配合的時候，我們先以鐵心柱的尺寸完全符合設計尺寸，即沒有誤差。以鐵心柱為參考，考察線圈內筒直徑的增加，以及直徑增加引起的多級階梯形某一級厚度的增加，最終導致的鐵心柱有效截面積的增加。能夠讓有效截面積增加到達最大的那個值就是公差的最大值。
據此我們建立起了非線性整數規劃模型二：
目標函數，線圈內同直徑增加，鐵心柱各級疊片厚度變化之后有效截面積達到最大 ：
約束條件：
重新計算的公差要在GB范圍之內：
其中D為某個確定的直徑，
表 公差取值
直徑（mm） 70-200 205-255 260-500 500-650 650-800
公差 (mm) 5 8 10 12 14
硅鋼片的級數不會變化，各級硅鋼片的寬度和理論最優設計尺寸相等，且不會變化
其中，D為之間的任意值，是直徑D所對應的最大公差。
是對應直徑為D的鐵心柱最優方案中，多級階梯形各級的寬度。
我們根據鐵心柱有效截面積最大的設計要求，讓處于直徑為D的時候，選擇合適的多級階梯形的級數，各級的厚度和寬度的最優方案，此計算結果我們放在附錄里。
改變多級階梯形各級厚度之后，鐵心柱在增大直徑后的線圈內筒之內，即各級尺寸滿足勾股定理：
因此，建立數學模型如下
s.t
模型求解：
運用LINGO軟件，求得的結果如下：
表 公差
直徑（mm） 70-200 205-255 260-500 500-650 650-800
公差 (mm) 4.8 7.5 9.7 10.9 13.5
5.3.1
因為我們采用的方法是：先根據油道分割出來的各小塊的面積基本相等，確定出油道的具體位置，然后在滿足第一小問中關于制造工藝的約束的情況下，用硅鋼片填充各個分區塊，使各個區塊的鐵心面積最大，即以各個區塊鐵心柱面積最大為目標，然后面積之和最大為目標函數，求最優方案。
Setp.1
確定油道的位置和寬度，根據幾何關系可得如下方程組：
方程組
聯立求解求得油道的位置：
求得的第一根油道的尺寸為：
，
求得的第二根油道的尺寸為：
，
，
在獲得有油道的精確位置之后，然后用硅鋼片去填充每一個由油道分割出來的小區塊，并且讓填充的面積最大：
第一個區塊的填充面積為： ,并讓
第二個區塊的填充面積為：, 并讓
第三個區塊的填充面積為：,并讓
總的填充面積，即為鐵心柱幾何截面積的二分之一
目標函數：
Setp.2
第一根油道與水平的那根直徑之間的填充最優方案：
求解上述模型就可得到第一分塊的多級階梯截面積最大的最優方案。
Setp.3
求第一根油道和第二根油道之間，最優填充方案：
求解上述模型可以得到第一根油道和第二根油道之間區塊最佳填充方案。
Setp.4
求解第二根油道和圓頂部之間區塊的最優填充方案：
求解上述模型可求得第二根油道和圓頂之間的最優填充方案。
5.3.2 模型求解
利用MATLAB進行優化計算得到如下最優設計方案
1級 2級 3級 4級 5級 6級 7級 8級 9級 10級
第一分塊 645
40
第二分塊 620 600 570
82 55 62
第三分塊 635 630 625 615 610 600 595 585 580 570
25 21 19 32 14 26 12 22 10 19
第一模塊面積
51600
第二模塊面積
mj=59590
第三模塊面積
mj =60720
總面積（不包括油道）
292220
面積占用比
88.1%
7 模型評價與擴展
本文建立的模型屬于含多個整數變量的非線性整數規劃模型。對于問題中的約束，針對性地將其表示為程序容易讀取和識別的約束條件，在設立變量時，考慮了變量下標與已知數據以及與算法的匹配，使得程序有較大的可移植性和通用性。例如僅將第模型一中的的變量下標數字改變就能移植到第三問中進行優化求解。同時，對于題目中所建立的含多變量的非線性整數規劃模型，采用matlab 和LINGO軟件進行編程求解。
建立模型的時候，我們沒有充分考慮
模型的不足在于，由于LINGO軟件的限制，沒有效率較高的算法，使得軟件在求解程序時的速度受限。此外，可以采用MATLAB軟件進行遍歷搜索，結果與LINGO結果進行比較，從而判斷模型和LINGO算法的正確性。
模型推廣
本模型還可以推廣到鋼材，布料分割等要求原材料能夠充分利用，節約成本，改善產品使用性能等其他的生產活動中去，也可以應用的有關分配的問題中，具有較強的實用性，可操作性。
8 參考文獻
[1] 朱德通. 最優化模型與實驗. 上海：同濟大學出版社， 2003.
[2] 韓中庚 數學建模方法及其應用. 北京：高等教育出版社， 2005
[3] 鄧成梁 運籌學的原理和方法（第二版） 武漢：華中科技大學出版社，2006
[4] 王全保等 變壓器手冊-電子分冊 沈陽：遼寧科學技術出版社 1987
[5] 鄧英劍 楊冬生 公差配合與測量技術 北京：國防工業出版社 2007
[6]羅建軍 楊琦 MATLAB 西安：西安交通大學出版社 2002
[7]張立衛 最優化問題的擾動分析 北京：科學出版社 2007
9 附錄
第一問MATLAB程序
第二問MATLAB程序
第二問鐵心柱多級階梯形設計最優方案
第三問LINGO 程序
第一問程序
程序：
funf='f=-1*(x(1)*sqrt(325^2-x(1)^2)+x(2)*(sqrt(325^2-x(2)^2)-sqrt(325^2-x(1)^2))+x(3)*(sqrt(325^2-x(3)^2)-sqrt(325^2-x(2)^2))+x(4)*(sqrt(325^2-x(4)^2)-sqrt(325^2-x(3)^2))+x(5)*(sqrt(325^2-x(5)^2)-sqrt(325^2-x(4)^2))+x(6)*(sqrt(325^2-x(6)^2)-sqrt(325^2-x(5)^2))+x(7)*(sqrt(325^2-x(7)^2)-sqrt(325^2-x(6)^2))+x(8)*(sqrt(325^2-x(8)^2)-sqrt(325^2-x(7)^2))+x(9)*(sqrt(325^2-x(9)^2)-sqrt(325^2-x(8)^2))+x(10)*(sqrt(325^2-x(10)^2)-sqrt(325^2-x(9)^2))+x(11)*(sqrt(325^2-x(11)^2)-sqrt(325^2-x(10)^2))+x(12)*(sqrt(325^2-x(12)^2)-sqrt(325^2-x(11)^2)));' ;% 最大鐵芯柱柱截面面積函數
fung='g=[x(2)-x(1)+2.5;x(3)-x(2)+2.5;x(4)-x(3)+2.5;x(5)-x(4)+2.5;x(6)-x(5)+2.5;x(7)-x(6)+2.5;x(8)-x(7)+2.5;x(9)-x(8)+2.5;x(10)-x(9)+2.5;x(11)-x(10)+2.5;x(12)-x(11)+2.5];';%寬度逐級遞減的約束條件
fun=[funf fung];
x0=[ 37.5 35 32.5 30 27.5 25 22.5 20 17.5 15 12.5 10]; %初始值
options=[];
vlb=[ 37.5 35 32.5 30 27.5 25 22.5 20 17.5 15 12.5 10]; %取值下限
vub=[322.5 320 317.5 315 312.5 310 307.5 305 302.5 300 297.5 295]; %取值上限
[x,options]=constr(fun,x0,options,vlb,vub);
y=zeros(1,12);
x=x.*2;
x=(round(x./5)).*5; %以5為倍數的寬度
for i=1:12
yy=sum(y',1);
y(i)=sqrt(650^2-x(i)^2)-yy;
end
y=round(y);
x
save 1-12x.txt x -ASCII -TABS
y
save 1-12y.txt y -ASCII -TABS
s=x.*y;
mj=sum(s',1);
mj
save 1-12mj.txt mj -ASCII -TABS
zhan=mj/((325^2)*pi);
zhan
save 1-12zhan.txt zhan -ASCII -TABS
結果輸出：
> In C:\MATLAB6p5\toolbox\optim\constr.m at line 55
第二問程序
1）
MODEL:
DATA:
N=12;
d=650;
ENDDATA
SETS:
can/1..N/:x ,Y,T;
ENDSETS
@FOR(can:X^2+Y^2<0.25*d^2);
@FOR(can(I)|I#GT#1:X(I-1)>X(I););
@FOR(can(I)|I#GT#1:Y(I-1)<Y(I););
@FOR(can:T=X/2.5);
@FOR(can:@GIN(T));
Y(1)>13;
X(1)<325;
X(N)>20;
)| I #GT# 1:X(I)*Y(I-1));
END
表 鐵心柱多級階梯形最優方案
直徑D 1級寬 2級寬 3級寬 4級寬 5級寬 6級寬 7級寬 8級寬 9級寬 10級寬 11級寬 12級寬 13級寬 14級寬 Ｓ Ｓ 面積利用率
80 65 70 65 60 55 50 40 5027 4520 90%
110 105 100 95 85 75 60 40 9503 8796 93%
140 135 130 115 100 85 65 40 15394 14320 93%
170 165 150 130 115 95 70 40 22698 21144 93%
195 185 165 145 130 110 90 55 29865 27740 93%
200 195 158 170 155 135 115 95 70 40 31416 29936 95%
230 225 215 205 190 170 150 130 110 80 45 41548 39588 95%
265 260 250 235 195 185 150 125 90 55 55154 52564 95%
270 265 250 230 230 105 180 155 155 125 95 55 57256 54236 95%
330 325 315 300 285 285 260 230 200 165 125 75 85530 81456 95%
390 395 370 350 330 305 275 245 210 170 125 75 119459 114360 96%
400 395 385 370 350 350 335 300 275 245 215 195 150 110 65 125664 121104 96%
530 520 505 485 455 425 390 355 355 315 270 220 160 95 220618 211996 96%
740 650 630 605 580 540 505 505 470 450 385 340 280 205 120 430084 399728 93%
2）
model:
data
l=37.5,35,32.5,30,27.5,25,20;
d=80;
enddata
sets:
m/1..7/:w,l;
endsets
));
@for(m(i):4w(i)^2+4*@sum(m(i):d(i))*@sum(m(i):d(I))<=d^2+2*d*es+es^2);
es<=5;
end
第三問程序
1）
s1-a1*r*r-0.5*w1*h1=0;
s2-a2*r*r-0.5*w2*h2+s1=0;
s3-0.5*pi*r*r+s1+s2=0;
2*s1-s2=0;
s2-s3=0;
0.25*w1*w1+h1*h1-r*r=0;
0.25*w2*w2+h2*h2-r*r=0;
pi=3.1416;
r=325;
s1>0;
a1>0;
a2>0;
a1<0.5*pi;
a2<0.5*pi;
h1-w1* @TAN(a1)=0;
h2-w2* @TAN(a2)=0;
2）
程序：
1.求油道大致位置
eq1=sym('x2*y2*1/2+atan(y2/x2)*325^2=0.3*((325^2)*pi)/5');
eq2=sym('x1*y1*1/2+atan(y1/x1)*325^2=0.1*((325^2)*pi)/5');
eq3=sym('x1^2+y1^2=325^2');
eq4=sym('x2^2+y2^2=325^2');
[x1,y1,x2,y2]=solve(eq1,eq2,eq3,eq4)
2.最優化取數
funf='f=-1*(x(1)*sqrt(325^2-x(1)^2)+x(2)*(sqrt(325^2-x(2)^2)-sqrt(325^2-x(1)^2))+x(3)*(sqrt(325^2-x(3)^2)-sqrt(325^2-x(2)^2))+x(4)*(sqrt(325^2-x(4)^2)-sqrt(325^2-x(3)^2))+x(5)*(sqrt(325^2-x(5)^2)-sqrt(325^2-x(4)^2))+x(6)*(sqrt(325^2-x(6)^2)-sqrt(325^2-x(5)^2))+x(7)*(sqrt(325^2-x(7)^2)-sqrt(325^2-x(6)^2))+x(8)*(sqrt(325^2-x(8)^2)-sqrt(325^2-x(7)^2))+x(9)*(sqrt(325^2-x(9)^2)-sqrt(325^2-x(8)^2)));' ;% 最大鐵芯柱柱截面面積函數
fung='g=[x(2)-x(1)+2.5;x(3)-x(2)+2.5;x(4)-x(3)+2.5;x(5)-x(4)+2.5;x(6)-x(5)+2.5;x(7)-x(6)+2.5;x(8)-x(7)+2.5;x(9)-x(8)+2.59];';%寬度逐級遞減的約束條件
fun=[funf fung];
x0=[305 302.5 300 297.5 295 292.5 290 287.5 285]; %初始值
options=[];
vlb=[305 302.5 300 297.5 295 292.5 290 287.5 285]; %取值下限
vub=[320 317.5 315 312.5 310 307.5 305 302.5 300]; %取值上限
[x,options]=constr(fun,x0,options,vlb,vub);
y=zeros(1,9);
x=x.*2;
x=(round(x./5)).*5; %以5為倍數的寬度
for i=1:9
yy=sum(y',1);
y(i)=sqrt(650^2-x(i)^2)-yy;
end
y=round(y);
x
y
s=x.*（y-156）;
mj=sum(s',1);
mj
【電力變壓器鐵心柱截面的優化設計(一)】相關文章：
電力工程施工組織設計及其優化研究03-14
淺談電力變壓器故障原因及處理方法05-20
DOS用戶界面的設計03-20
電力工程施工組織設計及其優化研究理工論文11-17
建筑設計中板柱結構的適用高度03-07
提高建筑給排水設計節能優化03-14
談減速器設計的優化03-18
網站界面的功能介紹及UI設計06-04
