机械优化设计案例.docVIP

機械優化設計案例1 1. 題目 對一對單級圓柱齒輪減速器，以體積最小為目標進行優化設計。 2.已知條件 已知數輸入功p=58kw，輸入轉速n1=1000r/min，齒數比u=5，齒輪的許用應力[]H=550Mpa，許用彎曲應力[]F=400Mpa。 建立優化模型 3.1問題分析及設計變數的確定 由已知條件得求在滿足零件剛度和強度條件下，使減速器體積最小的各項設計參數。由於齒輪和軸的尺寸（即殼體內的零件）是決定減速器體積的依據，故可按它們的體積之和最小的原則建立目標函數。 單機圓柱齒輪減速器的齒輪和軸的體積可近似的表示為： 式中符號意義由結構圖給出，其計算公式為 由上式知，齒數比給定之後，體積取決於b、z1 、m、l、dz1 和dz2 六個參數，則設計變數可取為 3.2目標函數為 3.3約束條件的建立 為避免發生根切，應有，得 2 )齒寬應滿足，和為齒寬係數的最大值和最小值，一般取=0.9，=1.4，得 動力傳遞的齒輪模數應大於2mm，得 為了限制大齒輪的直徑不至過大，小齒輪的直徑不能大於，得 齒輪軸直徑的範圍：得 軸的支撐距離按結構關係，應滿足條件：（可取=20），得 齒輪的接觸應力和彎曲應力應不大於許用值，得 齒輪軸的最大撓度不大於許用值，得 齒輪軸的彎曲應力不大於許用值，得 4.優化方法的選擇 由於該問題有6個設計變數，16個約束條件的優化設計問題，採用傳統的優化設計方法比較繁瑣，比較複雜，所以選用Matlab優化工具箱中的fmincon函數來求解此非線性優化問題，避免了較為繁重的計算過程。 5.數學模型的求解 5.1.1將已知及資料代入上式，該優化設計的數學優化模型表示為： Subject to: 5.1.2運用Matlab優化工具箱對數學模型進行程式求解 首先在Matlab優化工具箱中編寫目標函數的M檔 myfun.m,返回x處的函數值f： function f = myfun(x) f=0.785398*(4.75*x(1)*x(2)^2*x(3)^2+85*x(1)*x(2)*x(3)^2-85*x(1)*x(3)^2+0.92*x(1)*x(6)^2-x(1)*x(5)^2+0.8*x(1)*x(2)*x(3)*x(6)-1.6*x(1)*x(3)*x(6)+x(4)*x(5)^2+x(4)*x(6)^2+28*x(5)^2+32*x(6)^2) 由於約束條件中有非線性約束，故需要編寫一個描述非線性約束條件的M檔mycon.m： function[c,ceq]=myobj(x) c=[17-x(2);0.9-x(1)/(x(2)*x(3));x(1)/(x(2)*x(3))-1.4;2-x(3);x(2)*x(3)-300;100-x(5);x(5)-150;130-x(6);x(6)-200;x(1)+0.5*x(6)-x(4)-40;1486250/(x(2)*x(3)*sqrt(x(1)))-550; 7098/(x(1)*x(2)*x(3)^2*(0.169+0.006666*x(2)-0.0000854*x(2)^2))-400;7098/(x(1)*x(2)*x(3)^2*(0.2824+0.00177*x(2)-0.0000394*x(2)^2))-400;117.04*x(4)^4/(x(2)*x(3)*x(5)^4)-0.003*x(4);(1/(x(5)^3))*sqrt((2850000*x(4)/(x(2)*x(3)))^2+2.4*10^12)-5.5;(1/(x(6)^3))*sqrt((2850000*x(4)/(x(2)*x(3)))^2+6*10^13)-5.5]; ceq=[]; 最後在command window裡輸入： x0=[230;21;8;420;120;160];%給定初始值 [x,fval,exitflag,output]=fmincon(@myfun,x0,[],[],[],[],[],[],@myobj,output) %調用優化過程 5.1.3最優解以及結果分析 運行結果如下圖所示： 由圖可知，優化後的最終結果為 x=[123.3565 99.8517 1.7561 147.3157 150.4904 129.5096] f(x)=2.36e*107 由於齒輪模數應為標準值，齒數必須為整數，其它參數也要進行圓整,所以最優解不能直接採用，按設計規範，經標準化和圓整後： x=[124 100 2 148 150 130] f(x)=6.16 *107 結果對比分析 若按初始值減速器的體積V大約為6.32×107mm