【導(dǎo)讀】數(shù)處理，只考慮鐵心柱的幾何截面積，建立了簡(jiǎn)單的單目標(biāo)規(guī)劃模型。心柱截面的級(jí)數(shù)應(yīng)為14，各級(jí)的寬度和厚度見正文。外接圓直徑的公差帶。序列法將其轉(zhuǎn)化為只考慮截面形狀的單目標(biāo)規(guī)劃模型，用lingo軟件編程實(shí)現(xiàn)，考慮國(guó)標(biāo)中關(guān)于公差的要求，得出鐵心柱外接圓直徑的公差帶為650—652mm，線圈的內(nèi)筒直徑的公差帶為—。指出油道的位置。我們分別在前兩問中增加“被油道分割的相鄰兩部分鐵心柱截。方法，求得的結(jié)果完全滿足題目的要求，并有一定的實(shí)際應(yīng)用價(jià)值。壓器制造業(yè)通常采用全國(guó)統(tǒng)一的標(biāo)準(zhǔn)鐵心設(shè)計(jì)圖紙。上下軸對(duì)稱，左右也軸對(duì)稱。階梯形的每級(jí)都是由許多同種寬度的硅鋼片迭起來。求第一級(jí)的厚度最小為26毫米，硅鋼片的寬度最小為20毫米。而疊片系數(shù)通常與硅鋼片厚度、表面的絕緣漆膜厚度、硅鋼片的平整。程度均相同，即鐵心柱的疊片系數(shù)處處相等，為一常數(shù)；據(jù)我們的結(jié)論，為生產(chǎn)廠家提供一種更好的設(shè)計(jì)，提高變壓器的性價(jià)比。

　　

【正文】 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 模型 2的 0d? 時(shí)程序 MODEL: Title 2。 sets: !y 為厚度 ,x 為寬度的一半 ,臨時(shí)變量 k 用于保證寬度為 5 的倍數(shù) 。 27 Piece/1..14/:y, x, k。 d1/1..13/:r1。 endsets !目標(biāo)函數(shù) 。 min=@abs(dis)。 !d 為直徑增加的量 ,。 dis=(@sum(d1(i):r1(i)) + r2)/27 650/2。 r2= ((650+d)/2)*14。 !s 為幾何面積 。 s=(@sum(Piece(i): x(i)*y(i))*x(1)*y(1))*4。 r1(1) =(x(2)^2 + (@sum(Piece(i)|iLE1 :y(i))*y(1))^2 )^。 r1(2) =(x(3)^2 + (@sum(Piece(i)|iLE2 :y(i))*y(1))^2 )^。 r1(3) =(x(4)^2 + (@sum(Piece(i)|iLE3 :y(i))*y(1))^2 )^。 r1(4) =(x(5)^2 + (@sum(Piece(i)|iLE4 :y(i))*y(1))^2 )^。 r1(5) =(x(6)^2 + (@sum(Piece(i)|iLE5 :y(i))*y(1))^2 )^。 r1(6) =(x(7)^2 + (@sum(Piece(i)|iLE6 :y(i))*y(1))^2 )^。 r1(7) =(x(8)^2 + (@sum(Piece(i)|iLE7 :y(i))*y(1))^2 )^。 r1(8) =(x(9)^2 + (@sum(Piece(i)|iLE8 :y(i))*y(1))^2 )^。 r1(9) =(x(10)^2 + (@sum(Piece(i)|iLE9 :y(i))*y(1))^2 )^。 r1(10)=(x(11)^2 + (@sum(Piece(i)|iLE10 :y(i))*y(1))^2 )^。 r1(11)=(x(12)^2 + (@sum(Piece(i)|iLE11 :y(i))*y(1))^2 )^。 r1(12)=(x(13)^2 + (@sum(Piece(i)|iLE12 :y(i))*y(1))^2 )^。 r1(13)=(x(14)^2 + (@sum(Piece(i)|iLE13 :y(i))*y(1))^2 )^。 !約束第一級(jí)厚度和長(zhǎng)度 。 x(1)=10 。 y(1)=26。 !各級(jí)厚度和寬度的一半符合勾股定理 。 x(1)^2+(*y(1))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 2 :(y(i)))*y(1))^2 + (x(2))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 3 :(y(i)))*y(1))^2 + (x(3))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 4 :(y(i)))*y(1))^2 + (x(4))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 5 :(y(i)))*y(1))^2 + (x(5))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 6 :(y(i)))*y(1))^2 + (x(6))^2 = ((650+d)/2)^2。 28 (@sum(Piece(i)|i LE 7 :(y(i)))*y(1))^2 + (x(7))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 8 :(y(i)))*y(1))^2 + (x(8))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 9 :(y(i)))*y(1))^2 + (x(9))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 10:(y(i)))*y(1))^2 + (x(10))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 11:(y(i)))*y(1))^2 + (x(11))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 12:(y(i)))*y(1))^2 + (x(12))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 13:(y(i)))*y(1))^2 + (x(13))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 14:(y(i)))*y(1))^2 + (x(14))^2 = ((650+d)/2)^2。 !總幾何面積小于外接圓面積大于問題一所求解出的面積 。 s=*((650+d)/2)^2。 s=。 !d 為直徑增加的量 。 d=0。 !隨著級(jí)數(shù)的增加硅鋼片寬度依次減小 。 @for(Piece(i)|iGT1:x(i1)x(i))。 !約束硅鋼片寬度必須為 5 的倍數(shù) 。 @for(Piece(i):x(i)*2=5*k(i))。 @for(Piece(i):@Gin(k(i)))。 @free(dis)。 END 模型 2的最終模型的程序 MODEL: Title 2。 sets: !y 為厚度 ,x 為寬度的一半 ,臨時(shí)變量 k 用于保證寬度為 5的倍數(shù) 。 Piece/1..14/:y, x, k。 d1/1..13/:r1。 endsets !目標(biāo)函數(shù) 。 min=dis。 !d 為直徑增加的量 ,。 dis=@abs((@sum(d1(i):r1(i)) + r2)/27 650/2)。 r2= ((650+d)/2)*14。 !s 為幾何面積 。 29 s=(@sum(Piece(i): x(i)*y(i))*x(1)*y(1))*4。 r1(1) =(x(2)^2 + (@sum(Piece(i)|iLE1 :y(i))*y(1))^2 )^。 r1(2) =(x(3)^2 + (@sum(Piece(i)|iLE2 :y(i))*y(1))^2 )^。 r1(3) =(x(4)^2 + (@sum(Piece(i)|iLE3 :y(i))*y(1))^2 )^。 r1(4) =(x(5)^2 + (@sum(Piece(i)|iLE4 :y(i))*y(1))^2 )^。 r1(5) =(x(6)^2 + (@sum(Piece(i)|iLE5 :y(i))*y(1))^2 )^。 r1(6) =(x(7)^2 + (@sum(Piece(i)|iLE6 :y(i))*y(1))^2 )^。 r1(7) =(x(8)^2 + (@sum(Piece(i)|iLE7 :y(i))*y(1))^2 )^。 r1(8) =(x(9)^2 + (@sum(Piece(i)|iLE8 :y(i))*y(1))^2 )^。 r1(9) =(x(10)^2 + (@sum(Piece(i)|iLE9 :y(i))*y(1))^2 )^。 r1(10)=(x(11)^2 + (@sum(Piece(i)|iLE10 :y(i))*y(1))^2 )^。 r1(11)=(x(12)^2 + (@sum(Piece(i)|iLE11 :y(i))*y(1))^2 )^。 r1(12)=(x(13)^2 + (@sum(Piece(i)|iLE12 :y(i))*y(1))^2 )^。 r1(13)=(x(14)^2 + (@sum(Piece(i)|iLE13 :y(i))*y(1))^2 )^。 !約束第一級(jí)厚度和長(zhǎng)度 。 x(1)=10 。 y(1)=26。 !各級(jí)厚度和寬度的一半符合勾股定理 。 x(1)^2+(*y(1))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 2 :(y(i)))*y(1))^2 + (x(2))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 3 :(y(i)))*y(1))^2 + (x(3))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 4 :(y(i)))*y(1))^2 + (x(4))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 5 :(y(i)))*y(1))^2 + (x(5))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 6 :(y(i)))*y(1))^2 + (x(6))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 7 :(y(i)))*y(1))^2 + (x(7))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 8 :(y(i)))*y(1))^2 + (x(8))^2 = ((650+d)/2)^2。 (@sum(Piece(i)|i LE 9 :(y(i)))*y(1))^2 + (x(9))^2 = ((