- 相關(guān)推薦
電力變壓器鐵心柱截面的優(yōu)化設(shè)計(一)
摘要
本文針對變壓器鐵心柱截面優(yōu)化設(shè)計，建立數(shù)學(xué)模型并用matlab函數(shù)constr編程以及使用LINGO編程求解。
最優(yōu)鐵心柱截面設(shè)計方案指的是在達(dá)到使用標(biāo)準(zhǔn)的情況下，讓鐵心柱橫截面
的有效截面面積最大，從而提高產(chǎn)品的使用性能，提高產(chǎn)品的競爭力。
問題一中，鐵心柱外接圓直徑為650Mm，根據(jù)截面級數(shù)的選擇表，可以確定其級數(shù)的范圍。根據(jù)制造工藝要求的第一級厚度最小26和最后一級寬度要求最小為20，各級之間寬度依次減少為約束，以有效截面積最大為目標(biāo)，建立起含多個整數(shù)變量的非線性整數(shù)規(guī)劃模型——模型一，優(yōu)化求解以確定最優(yōu)的多級硅鋼片階梯形組成方案。
利用MATLAB軟件規(guī)劃求解得到最優(yōu)的多級階梯形硅鋼片組成方案為：
鐵心柱有效截面積320240，鐵心利用率96.50%，級數(shù)為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）
對于問題二和問題三，我們將第一問的模型進行拓展和延伸。
第二問要求線圈內(nèi)筒和鐵心柱外接圓直徑的公差。以有效截面積利用率最大為目標(biāo)，以變壓器設(shè)計結(jié)構(gòu)設(shè)計要求為約束，我們再次建立了一個非線性整數(shù)規(guī)劃模型。
在考慮添加油道的問題中，需要確定油道位置，并且分割的各個面積要求基本相等。我們首先確定油道的位置，然后在分割出來的各個小塊里，分別去進行多級階梯形的最優(yōu)化設(shè)計。各個分塊里實現(xiàn)了最優(yōu)化，整個設(shè)計就實現(xiàn)了最優(yōu)化。
該過程中，油道和各個分塊之間會產(chǎn)生間隙，針對此問題，我們將油道的位置向圓心靠近調(diào)整，然后新產(chǎn)生的間隙，用下一個分塊的第一級硅鋼片進行補充。
由于鐵心柱疊片系數(shù)是一個常數(shù)，在模型求解的過程中，可以忽略其影響。
在建模的過程中，我們運用了MATLAB和LINGO軟件進行求解。
關(guān)鍵詞:非線整數(shù)性規(guī)劃，鐵心利用率，分布填充
問題重述
電力變壓器的設(shè)計中很重要的一個環(huán)節(jié)就是鐵心柱的截面如何設(shè)計。我國變壓器制造業(yè)通常采用全國統(tǒng)一的標(biāo)準(zhǔn)鐵心設(shè)計圖紙。根據(jù)多年的生產(chǎn)經(jīng)驗，各生產(chǎn)廠存在著對已有設(shè)計方案的疑問：能否改進及如何改進這些設(shè)計，才能在提高使用效果的同時降低變壓器的成本。
現(xiàn)在以心式鐵心柱為例試圖進行優(yōu)化設(shè)計。
電力變壓器鐵心柱截面在圓形的線圈筒里面。為了充分利用線圈內(nèi)空間又便于生產(chǎn)管理，心式鐵心柱截面常采用多級階梯形結(jié)構(gòu)，如圖1所示。截面在圓內(nèi)上下軸對稱，左右也軸對稱。階梯形的每級都是由許多同種寬度的硅鋼片迭起來的。由于制造工藝的要求，硅鋼片的寬度一般取為5的倍數(shù)（單位：毫米）。因為在多級階梯形和線圈之間需要加入一定的撐條來起到固定的作用，所以一般要求第一級的厚度最小為26毫米，硅鋼片的寬度最小為20毫米。
鐵心柱有效截面的面積，等于多級鐵心柱的幾何截面積（不包括油道）乘以疊片系數(shù)。而疊片系數(shù)通常與硅鋼片厚度、表面的絕緣漆膜厚度、硅鋼片的平整度以及壓緊程度有關(guān)。設(shè)計時希望有效截面盡量大，既節(jié)省材料又減少能量損耗。顯然鐵心柱的級數(shù)愈多，其截面愈接近于圓形，在一定的直徑下鐵心柱有效截面也愈大。但這樣制造也工藝復(fù)雜，一般情況下鐵心柱的級數(shù)可參照表1選取。
表1 鐵心柱截面級數(shù)的選擇
鐵心柱直徑mm 級數(shù)
80-195 5-7
200-265 8-10
270-390 11
400-740 12-14
760以上 >15
問題一：當(dāng)鐵心柱外接圓直徑為650毫米時，如何確定鐵心柱截面的級數(shù)、各級寬度和厚度，才能使鐵心柱的有效截面積最大。
問題二：實際生產(chǎn)中線圈的內(nèi)筒直徑和鐵心柱的外接圓直徑不是精確地相等，而留有一定的間隙以便于安裝和維修，設(shè)計的兩個直徑的取值范圍稱為各自的公差帶。因此可以在設(shè)計鐵心截面時稍微增加鐵心柱的外接圓的直徑以使得鐵心柱有更好的截面形狀。請結(jié)合鐵心柱截面的設(shè)計而設(shè)計出二者的公差帶。
問題三：銅導(dǎo)線在電流流過時發(fā)熱造成的功率損耗簡稱為銅損；鐵心在磁力線通過時發(fā)熱造成的功率損耗簡稱為鐵損。為了改善鐵心內(nèi)部的散熱，鐵心柱直徑為380毫米以上時須設(shè)置冷卻油道。簡單地說，就是在某些相鄰階梯形之間留下6毫米厚的水平空隙（如圖2所示），空隙里充滿油，變壓器工作時油上下循環(huán)帶走鐵心里的熱量。具體油道數(shù)可按表2選取。油道的位置應(yīng)使其分割的相鄰兩部分鐵心柱截面積近似相等。
分別針對問題一和問題二的情況，增加油道要求再給出設(shè)計，并指出油道的位置。
表2 冷卻油道數(shù)的選擇
鐵心柱直徑mm 半圓中6mm油道個數(shù)
380-410 0
420-500 1
510-690 2
700-840 3
2 條件假設(shè)
1.整個鐵心柱的硅鋼片出了長度，其它如厚度，表面絕緣漆膜厚度，平整度都相同；
2.硅鋼片之間是沒有形變的壓緊；
3.疊片系數(shù)是確定的已知數(shù)；
4.油道對稱分布；
5.不考慮工藝過程的影響。
3 符號說明
第i級疊片的的厚度；
第i級疊片的寬度；
疊片系數(shù)；
鐵心柱理論外接圓的直徑，也等于理論線圈內(nèi)筒直徑；
第i根油道與直徑之間的；
油道分割出來的分塊面積；
多級階梯形前i級厚度之和；
線圈內(nèi)筒的公差值；
4 問題分析
4. 問題背景
變壓器是一種應(yīng)用電磁感應(yīng)原理把電能從一個電路傳到另一個電路的電磁裝置。它在電路中起變壓，變流，變電阻的作用。它由三部分組成：鐵心，起導(dǎo)磁，助磁作用；而是初級線圈，接電源，起激磁作用；三是次級線圈接負(fù)載，利用不同次初級線圈匝數(shù)比，實現(xiàn)變壓，變流，變電阻的作用。
在變壓器的構(gòu)成里，鐵心柱是很一個十分重要的組成部件。
因為鐵心柱的形狀，截面積，疊片的選擇，疊片的相關(guān)工藝過程都會影響將來變壓器的使用效果和壽命，以及使用成本。
我國變壓器制造業(yè)通常采用全國統(tǒng)一的標(biāo)準(zhǔn)鐵心設(shè)計圖紙，根據(jù)多年的生產(chǎn)經(jīng)驗，在長期的生產(chǎn)研究過程中，各生產(chǎn)廠產(chǎn)生了對已有設(shè)計方案的疑問：能否改進及如何改進這些設(shè)計，才能在提高使用效益的同時降低變壓器成本。
4.1.
鐵心柱是安裝在線圈筒里面的，理論上，在線圈的直徑確定了的情況下，鐵心柱的有效面積越大，鐵心的電阻越大，使鐵心的鐵損最小，因而可以減少能量損耗，變壓器的使用性能會越好，使用壽命較長。為了充分利用空間和便于生產(chǎn)，鐵心柱截面長采用多級階梯形結(jié)構(gòu)，用不同長度的硅鋼片，堆疊成不同厚度的級，并且選擇合適的級數(shù)去逼近與之配合的鐵心線圈圓，期望得到的有效面積最大，獲得更大的電阻。
截面優(yōu)化設(shè)計是以保證到達(dá)設(shè)計標(biāo)準(zhǔn)為前提，盡可能改善和提高產(chǎn)品的使用效果，使產(chǎn)品競爭力提升。
4.2
公差是生產(chǎn)中允許工件尺寸和幾何形狀變動的范圍，用來限制誤差。工件的誤差在公差范圍內(nèi)才為合格。規(guī)定公差是為了保證產(chǎn)品使用性能的前提下，給出盡可能大的公差。因為公差越小，精度越高，生產(chǎn)成本就越高。去較大的公差，有利于控制成本。
設(shè)計時要求留有一定的間隙邊緣和維修，因此鐵心柱和外圍線圈的配合應(yīng)該是間隙配合。
我們根據(jù)《變壓器設(shè)計手冊》查詢得到的GB的線圈內(nèi)筒和鐵心柱的直徑差
公差
直徑（mm） 70-200 205-255 260-500 500-650 650-800
公差 (mm) 5 8 10 12 14
理論數(shù)值的鐵心柱外接圓和線圈內(nèi)筒直徑，它們在配合的時候，采用間隙配合。由于生產(chǎn)硅鋼片的時候，是以滿足鐵心柱的設(shè)計尺寸為參照的，用來裝配鐵心柱的時候，它們的寬度是一系列已經(jīng)確定的數(shù)值，不可能因為包裹鐵心柱的線圈內(nèi)筒的直徑的增大而增大。
即當(dāng)線圈內(nèi)筒的直徑變化的時候，多級階梯形的各級的寬度是不可能變化的。為了獲得更好的截面積，或者增大截面積的值，我們可以將多級階梯形的某一級的疊片數(shù)增加，從而增加該級厚度，去逼近實際尺寸的線圈內(nèi)筒。
增大線圈內(nèi)筒的直徑，再增大某級的厚度，獲得更大的有效截面積。我們認(rèn)定有讓效截面積增加率最大的那個偏差，就是鐵心柱公差的最大值。
據(jù)此選擇合適的公差帶，來獲得最大的有效面積增加率。
4.3.1
由于變壓器工作時鐵心柱有鐵心損耗，外圍線圈有電阻，電流通過的時候會產(chǎn)生熱量。為了改善鐵心散熱，在鐵心柱的直徑達(dá)到380mm以上，就需要增加油道，讓油循環(huán)達(dá)到散熱的目的。
在考慮添加油道的時候，需要讓油道分割的各部分的面積近似相等。為了達(dá)到這個要求，我們以各個分割面積相等為約束，對稱布置油道，建立起方程組，用以先確定下油道的具體位置。
確定油道位置之后，再根據(jù)鐵心柱有效截面積最大的要求，可以分步用硅鋼片將各個油道之間的分塊堆疊。在各個分區(qū)里面，在滿足各項設(shè)計要求的情況下，求堆疊的最優(yōu)方案，使得各個分區(qū)的截面積最大。最后將各個面積進行疊加，各個分區(qū)的硅鋼片進行參數(shù)匯總，就可得最終的最優(yōu)方案。
4.3.2
對于變壓器鐵心柱的級寬，都可以通過已知的直徑算出相應(yīng)的厚度，因此可以轉(zhuǎn)化為以各級的寬度為變量的數(shù)學(xué)模型。由于現(xiàn)有的各直徑鐵心柱的級數(shù)都是參照生產(chǎn)經(jīng)驗，認(rèn)為要達(dá)到一定的級數(shù)才能滿足一定的幾何截面積，而沒有科學(xué)的精確的求解過程。因此設(shè)想通過對模型的計算所得出的數(shù)據(jù)，只要在滿足一定的幾何截面積的前提下，如果能減少級數(shù)，便可減少生產(chǎn)成本。
5 模型建立求解
5.1.1
問題一的數(shù)學(xué)模型：
當(dāng)鐵心柱的外接圓直徑為650mm時，根據(jù)要求可知級數(shù)應(yīng)該在12到14之間。由變壓器鐵心設(shè)計的各項要求為約束條件，我們以鐵心柱有效截面積最大為目標(biāo)函數(shù)，建立了一個非線性整數(shù)規(guī)劃的模型。
約束條件
硅鋼片的寬度為5mm的倍數(shù)：
（ k為正整數(shù)）
多級階梯形的第一級厚度最小為26mm:
硅鋼片的寬度最小為20mm，就是最后一級的寬度最小為20mm
最大寬度為直徑：
多級階梯形中，每級的寬度必須是遞減的：
（）
I級及i級之前的厚度之和與i級的寬度之間還必須滿足勾股定理：
該多級階梯形的幾何截面積為：
其中（n=12，13，14）
目標(biāo)函數(shù)：
我們認(rèn)為所有硅鋼片的疊片系數(shù)是相同的，即是確定的常數(shù)，
求有效截面積最大：
綜上可知：
為確定鐵心柱直徑為650mm時候的最佳設(shè)計方案，以有效截面積最大為目標(biāo)，在滿足各項結(jié)構(gòu)設(shè)計要求的約束下，我們建立了模型一：
目標(biāo)函數(shù)：
5.1.2 模型一的求解：
將以上模型代入編寫的matlab程序中予以求解（程序見附件），得到如下最優(yōu)化
表 I 直徑650mm方案對比
硅鋼片參數(shù) 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
對比可知，所要求的鐵心柱有效截面積最大的最優(yōu)方案為14級，各級寬度和長度可由上表可得
5.2.1問題二模型
公差等級是確定尺寸精度的等級，國家標(biāo)準(zhǔn)是按照不同的基本尺寸分為20個等級。由于電力變壓器不像機械設(shè)備那樣需要承受大的負(fù)荷，比如說軸承那樣的高速旋轉(zhuǎn)，像發(fā)動機那樣的劇烈復(fù)雜的機械運動，承受如此復(fù)雜的力，振動，高溫高壓等復(fù)雜嚴(yán)苛的工作狀況，我們認(rèn)為電力變壓器的公差等級為18，已經(jīng)足夠滿足它的使用了。
根據(jù)變壓器設(shè)計即為合格產(chǎn)品的最大范圍區(qū)間，然后在這個范圍內(nèi)根據(jù)最優(yōu)化截面的思想，即求最大的有效截面增加率的所對應(yīng)的es值，就是它的公差最大值。
考慮鐵心柱和線圈內(nèi)筒的配合的時候，我們先以鐵心柱的尺寸完全符合設(shè)計尺寸，即沒有誤差。以鐵心柱為參考，考察線圈內(nèi)筒直徑的增加，以及直徑增加引起的多級階梯形某一級厚度的增加，最終導(dǎo)致的鐵心柱有效截面積的增加。能夠讓有效截面積增加到達(dá)最大的那個值就是公差的最大值。
據(jù)此我們建立起了非線性整數(shù)規(guī)劃模型二：
目標(biāo)函數(shù)，線圈內(nèi)同直徑增加，鐵心柱各級疊片厚度變化之后有效截面積達(dá)到最大 ：
約束條件：
重新計算的公差要在GB范圍之內(nèi)：
其中D為某個確定的直徑，
表 公差取值
直徑（mm） 70-200 205-255 260-500 500-650 650-800
公差 (mm) 5 8 10 12 14
硅鋼片的級數(shù)不會變化，各級硅鋼片的寬度和理論最優(yōu)設(shè)計尺寸相等，且不會變化
其中，D為之間的任意值，是直徑D所對應(yīng)的最大公差。
是對應(yīng)直徑為D的鐵心柱最優(yōu)方案中，多級階梯形各級的寬度。
我們根據(jù)鐵心柱有效截面積最大的設(shè)計要求，讓處于直徑為D的時候，選擇合適的多級階梯形的級數(shù)，各級的厚度和寬度的最優(yōu)方案，此計算結(jié)果我們放在附錄里。
改變多級階梯形各級厚度之后，鐵心柱在增大直徑后的線圈內(nèi)筒之內(nèi)，即各級尺寸滿足勾股定理：
因此，建立數(shù)學(xué)模型如下
s.t
模型求解：
運用LINGO軟件，求得的結(jié)果如下：
表 公差
直徑（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
因為我們采用的方法是：先根據(jù)油道分割出來的各小塊的面積基本相等，確定出油道的具體位置，然后在滿足第一小問中關(guān)于制造工藝的約束的情況下，用硅鋼片填充各個分區(qū)塊，使各個區(qū)塊的鐵心面積最大，即以各個區(qū)塊鐵心柱面積最大為目標(biāo)，然后面積之和最大為目標(biāo)函數(shù)，求最優(yōu)方案。
Setp.1
確定油道的位置和寬度，根據(jù)幾何關(guān)系可得如下方程組：
方程組
聯(lián)立求解求得油道的位置：
求得的第一根油道的尺寸為：
，
求得的第二根油道的尺寸為：
，
，
在獲得有油道的精確位置之后，然后用硅鋼片去填充每一個由油道分割出來的小區(qū)塊，并且讓填充的面積最大：
第一個區(qū)塊的填充面積為： ,并讓
第二個區(qū)塊的填充面積為：, 并讓
第三個區(qū)塊的填充面積為：,并讓
總的填充面積，即為鐵心柱幾何截面積的二分之一
目標(biāo)函數(shù)：
Setp.2
第一根油道與水平的那根直徑之間的填充最優(yōu)方案：
求解上述模型就可得到第一分塊的多級階梯截面積最大的最優(yōu)方案。
Setp.3
求第一根油道和第二根油道之間，最優(yōu)填充方案：
求解上述模型可以得到第一根油道和第二根油道之間區(qū)塊最佳填充方案。
Setp.4
求解第二根油道和圓頂部之間區(qū)塊的最優(yōu)填充方案：
求解上述模型可求得第二根油道和圓頂之間的最優(yōu)填充方案。
5.3.2 模型求解
利用MATLAB進行優(yōu)化計算得到如下最優(yōu)設(shè)計方案
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 模型評價與擴展
本文建立的模型屬于含多個整數(shù)變量的非線性整數(shù)規(guī)劃模型。對于問題中的約束，針對性地將其表示為程序容易讀取和識別的約束條件，在設(shè)立變量時，考慮了變量下標(biāo)與已知數(shù)據(jù)以及與算法的匹配，使得程序有較大的可移植性和通用性。例如僅將第模型一中的的變量下標(biāo)數(shù)字改變就能移植到第三問中進行優(yōu)化求解。同時，對于題目中所建立的含多變量的非線性整數(shù)規(guī)劃模型，采用matlab 和LINGO軟件進行編程求解。
建立模型的時候，我們沒有充分考慮
模型的不足在于，由于LINGO軟件的限制，沒有效率較高的算法，使得軟件在求解程序時的速度受限。此外，可以采用MATLAB軟件進行遍歷搜索，結(jié)果與LINGO結(jié)果進行比較，從而判斷模型和LINGO算法的正確性。
模型推廣
本模型還可以推廣到鋼材，布料分割等要求原材料能夠充分利用，節(jié)約成本，改善產(chǎn)品使用性能等其他的生產(chǎn)活動中去，也可以應(yīng)用的有關(guān)分配的問題中，具有較強的實用性，可操作性。
8 參考文獻
[1] 朱德通. 最優(yōu)化模型與實驗. 上海：同濟大學(xué)出版社， 2003.
[2] 韓中庚 數(shù)學(xué)建模方法及其應(yīng)用. 北京：高等教育出版社， 2005
[3] 鄧成梁 運籌學(xué)的原理和方法（第二版） 武漢：華中科技大學(xué)出版社，2006
[4] 王全保等 變壓器手冊-電子分冊 沈陽：遼寧科學(xué)技術(shù)出版社 1987
[5] 鄧英劍 楊冬生 公差配合與測量技術(shù) 北京：國防工業(yè)出版社 2007
[6]羅建軍 楊琦 MATLAB 西安：西安交通大學(xué)出版社 2002
[7]張立衛(wèi) 最優(yōu)化問題的擾動分析 北京：科學(xué)出版社 2007
9 附錄
第一問MATLAB程序
第二問MATLAB程序
第二問鐵心柱多級階梯形設(shè)計最優(yōu)方案
第三問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)));' ;% 最大鐵芯柱柱截面面積函數(shù)
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為倍數(shù)的寬度
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
結(jié)果輸出：
> 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
表 鐵心柱多級階梯形最優(yōu)方案
直徑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.最優(yōu)化取數(shù)
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)));' ;% 最大鐵芯柱柱截面面積函數(shù)
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為倍數(shù)的寬度
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
【電力變壓器鐵心柱截面的優(yōu)化設(shè)計(一)】相關(guān)文章：
電力工程施工組織設(shè)計及其優(yōu)化研究03-14
淺談電力變壓器故障原因及處理方法05-20
DOS用戶界面的設(shè)計03-20
電力工程施工組織設(shè)計及其優(yōu)化研究理工論文11-17
建筑設(shè)計中板柱結(jié)構(gòu)的適用高度03-07
提高建筑給排水設(shè)計節(jié)能優(yōu)化03-14
談減速器設(shè)計的優(yōu)化03-18