matlab(四連桿優(yōu)化設(shè)計)

1、機械優(yōu)化設(shè)計在matlab中的應用東南大學機械工程學院*一優(yōu)化設(shè)計目的：在生活和工作中，人們對于同一個問題往往會提出多個解決方案，并通過各方面的論證從中提取最佳方案。最優(yōu)化方法就是專門研究如何從多個方案中科學合理地提取出最佳方案的科學。由于優(yōu)化問題無所不在，目前最優(yōu)化方法的應用和研究已經(jīng)深入到了生產(chǎn)和科研的各個領(lǐng)域，如土木工程、機械工程、化學工程、運輸調(diào)度、生產(chǎn)控制、經(jīng)濟規(guī)劃、經(jīng)濟管理等，并取得了顯著的經(jīng)濟效益和社會效益。二優(yōu)化設(shè)計步驟：機械優(yōu)化設(shè)計的全過程一般可以分為如下幾個步驟：1）建立優(yōu)化設(shè)計的數(shù)學模型；2）選擇適當?shù)膬?yōu)化方法；3）編寫計算機程序；4）準備必要的初始數(shù)據(jù)并傷及計算；5）對

2、計算機求得的結(jié)果進行必要的分析。其中建立優(yōu)化設(shè)計數(shù)學模型是首要的和關(guān)鍵的一步，它是取得正確結(jié)果的前提。優(yōu)化方法的選取取決于數(shù)學模型的特點，例如優(yōu)化問題規(guī)模的大小，目標函數(shù)和約束函數(shù)的性態(tài)以及計算精度等。在比較各種可供選用的優(yōu)化方法時，需要考慮的一個重要因素是計算機執(zhí)行這些程序所花費的時間和費用，也即計算效率。建立數(shù)學模型的基本原則與步驟設(shè)計變量的確定；設(shè)計變量是指在優(yōu)化設(shè)計的過程中，不斷進行修改，調(diào)整，一直處于變化的參數(shù)稱為設(shè)計變量。設(shè)計變量的全體實際上是一組變量，可用一個列向量表示：x二，.，-.:.II。目標函數(shù)的建立；選擇目標函數(shù)是整個優(yōu)化設(shè)計過程中最重要的決策之一。當對某以設(shè)計性能有特

3、定的要求，而這個要求有很難滿足時，則針對這一性能進行優(yōu)化會得到滿意的效果。目標函數(shù)是設(shè)計變量的函數(shù)，是一項設(shè)計所追求的指標的數(shù)學反映，因此它能夠用來評價設(shè)計的優(yōu)劣。目標函數(shù)的一般表達式為：f（x）=，要根據(jù)實際的設(shè)計要求來設(shè)計目標函數(shù)。約束條件的確定。一個可行性設(shè)計必須滿足某些設(shè)計限制條件，這些限制條件稱為約束條件，簡稱約束。由若干個約束條件構(gòu)成目標函數(shù)的可行域，而可行域內(nèi)的所有設(shè)計點都是滿足設(shè)計要求的，一般情況下，其設(shè)計可行域可表示為在可行域中，任意設(shè)計點滿足全部約束條件，稱為可行解，但不是最優(yōu)解，而優(yōu)化設(shè)計就是要求出目標函數(shù)在可行域的最優(yōu)解。三實例分析（機械優(yōu)化設(shè)計P241頁例8-5）設(shè)計

4、一曲柄搖桿機構(gòu)如圖，要求：曲柄.從且已范圍內(nèi)變化。分析：1）設(shè)計變量的確定決定機構(gòu)尺寸的各桿長度，以及當搖桿按已知運動規(guī)律開始運行時，曲柄所載的位置角應列為設(shè)計變量，即：X=考慮到機構(gòu)的桿長按比例變化時，不會改變其運動規(guī)律，因此在計算時常取，而其他桿長則按比例取為的倍數(shù)。若取曲柄的初始位置角為極位角，則及相應搖桿位置角均為桿長的函數(shù)，幾何圖形關(guān)系如右圖，其關(guān)系式為：TOC o 1-5 h z=arcos（1）=arcos（2）將的長度代入上式（1）,（2）得到：=arcos=arcos因此，只有為獨立變量，設(shè)計變量減少，故最后的設(shè)計變量為：:；7=目標函數(shù)的建立目標函數(shù)可根據(jù)已知的運動規(guī)律與機

5、構(gòu)實際運動規(guī)律之間的偏差最小為指標來建立，即f(x)=min式中一期望輸出角，m輸入角等分數(shù)；實際輸出角，由下圖得:a)0W式中b)(3)將輸入角分成30等分，并用近似公式計算，可得目標函數(shù)的表達式：f(x)=由題意知，傳動角的變化范圍是，則上式中變量的最后形式可以寫成:arcos將帶入（4）（5）得:為當時的理想輸出角，其值在題目中已經(jīng)給出:3）約束條件的確定曲柄搖桿機構(gòu)應滿足曲柄存在條件，可得曲柄搖桿機構(gòu)的傳動角可得把約束條件簡化（=5x=g7(x)二云十云1.414陽16youhua4.15742.2909%最優(yōu)解fval=5.1899e-004%目標函數(shù)最優(yōu)點的值exitflag=%標

6、志值，5表示重要方向?qū)?shù)小于規(guī)定的容許范圍并且約束違背小于options.ToLConoutput=iterations:12%迭代次數(shù)funcCount:40%函數(shù)的評價次數(shù)lssteplength:1stepsize:7.6955e-005algorithm:medium-scale:SQP,Quasi-Newton,line-search%采用的中型算法firstorderopt:1.0832e-006constrviolation:-1.0852e-006%一階最優(yōu)性條件message:1x780char%跳出信息lambda=lower:2x1doubleupper:2x1doubl

7、eeqlin:0 x1doubleeqnonlin:0 x1doubleineqlin:5x1doubleineqnonlin:2x1doublegrad=1.0e-003*0.48880.4445%函數(shù)在最優(yōu)點處梯度信息hessian=%函數(shù)在最優(yōu)點處海塞矩陣0.00160.00750.00750.0468結(jié)果分析采用fmincon求解的最優(yōu)值:：.=4.1574;2.2909;采用算法：中型算法(mediun-scale)。這與課本給出的最優(yōu)解：=4.1286;2.3325,=0.0156相比，計算精度更高，最優(yōu)解的數(shù)值更精確，故計算準確度高。用matlab繪制輸入輸出曲線關(guān)系圖上圖中(單

8、位為“度”藍色的線代表曲柄搖桿機構(gòu)的實際輸出角與輸入角的關(guān)系，紅色的線代表理想輸出角與輸入角的關(guān)系?？梢钥闯觯簩嶋H輸出和理論輸出曲線之間存在線性誤差，其最大線性誤差為，誤差在允許的范圍之內(nèi)，故結(jié)果的可信度也較大，運用matlab優(yōu)化工具箱計算所得結(jié)果正確。小結(jié)通過結(jié)合實際問題的分析，計算，求解，更加深入地了解和掌握機械優(yōu)化設(shè)計的過程和步驟，比較重要的步驟是數(shù)學模型的建立，以及設(shè)計變量的選取，以及數(shù)學模型的尺度變換，根據(jù)機構(gòu)實際工作需要，建立目標函數(shù)的約束條件等等，當數(shù)學模型建好以后，剩下的工作可以再matlab里面完成，而matlab里面的優(yōu)化工具箱，給用戶提供了多種優(yōu)化函數(shù)，使用者只需要將數(shù)

9、學模型按要求編寫成子程序嵌入已有的優(yōu)化程序即可。在設(shè)計過程中也遇到一些困難，比如說在在用matlab計算時，計算機已知處于busy狀態(tài)，得不到函數(shù)的最優(yōu)解，最后反復的檢查，終于找的了其原因，是由于初始點選擇不恰當引起的，如果初始點選擇得好，可以節(jié)省計算時間和計算空間，故初始點的選取比較重要。附錄1.編寫目標函數(shù)M文件myfun.m：functionf=myfun(x)f=0;aO=acos(1+x(1)入2-x(2)%函數(shù)f賦初值入2+25)/(10*(1+x(1);%初始計算點曲柄和搖桿的角度bO=acos(1+x(1)入2-x(2)i=2;入225)/(10*x(2);while(i=31

10、)a(i)=a0+(pi/2)*(i/30);b(i)=b0+2*(a(i)a0)入2/(3*pi);r=sqrt(2610*cos(a(i);c(i)=acos(r2+x(2)A2x(1)人2)/(2*x(2)*r);d(i)=acos(r入2+24)/(10*r);ifa(i)=pie(i)=pic(i)d(i);elseifa(i)=2*pie(i)=pic(i)+d(i);end%設(shè)置迭代次數(shù)為30次%計算曲柄各分度的角度值%計算搖桿各分度的角度值%計算搖桿輸出的實際值enda(1)=a0;f=f+(b(i)e(i)入2)*(a(i)a(i1);i=i+1;%目標函數(shù)的計算end2.編寫非線性不等式約束M文件constrain.m：functioncceq=constrain(x)c=36x(1)八2x(2)八21.414*x(1)*x(2);x(1)入2+x(2)入2-1.414*x(1)*x(2)-16;ceq=;3.調(diào)用fmincon優(yōu)化函數(shù)，建立youhua.m文件：lb=1;1;%非線性不等式約束x0=4;2;A=-1,0;0,-1;-1,-1;1,-1;-1,1;b=-1;-1;-6;4;4;%設(shè)計變量的下界%迭代初始點%線性不等