張春幸 劉菊花

摘 要：分析了影響銑刀桿設(shè)計的強度、剛度等諸多因素，建立了優(yōu)化選取銑刀桿結(jié)構(gòu)參數(shù)的數(shù)學(xué)模型，并運用MATLAB語言對該數(shù)學(xué)模型進行求解。

關(guān)鍵字：銑刀桿；MATLAB 優(yōu)化設(shè)計

在切削加工過程中，優(yōu)化選擇切削參數(shù)的問題已成為現(xiàn)代化機械制造業(yè)中極為重要的經(jīng)濟問題之一，而傳統(tǒng)的依靠經(jīng)驗數(shù)據(jù)和查閱手冊獲取的切削參數(shù)，其生產(chǎn)效率低且加工成本偏高。因此，以提高切削加工效率、降低加工成本、獲得高質(zhì)量產(chǎn)品為目的而進行的切削參數(shù)優(yōu)化的研究，對提高經(jīng)濟效益具有重要意義。銑削加工是機械加工中一種重要的、常見的方法之一。銑削加工的精度高，其加工的質(zhì)量、安全性及可靠性與諸多因素有關(guān)，其中，銑刀桿的設(shè)計質(zhì)量就是一個重要因素。工廠里的經(jīng)驗法，雖然基本上滿足設(shè)計要求，但卻不是最優(yōu)、最理想的設(shè)計。隨著現(xiàn)代計算機的普及和發(fā)展，以及在計算機條件下數(shù)值計算方法的成熟，為銑刀桿結(jié)構(gòu)參數(shù)優(yōu)化設(shè)計提供了條件。本文通過一個銑刀桿結(jié)構(gòu)參數(shù)優(yōu)化設(shè)計實例，來說明工程銑刀桿結(jié)構(gòu)優(yōu)化設(shè)計的基本方法，建立優(yōu)化選取銑刀桿結(jié)構(gòu)參數(shù)的數(shù)學(xué)模型，并運用當(dāng)今流行的科學(xué)數(shù)學(xué)計算語言MATLAB對該數(shù)學(xué)模型進行求解。

1．問題的提出

普通銑刀桿自由端有集中載荷，即切削阻力P= 1000N，扭矩M=100N*m。懸臂伸出不小于100mm，刀桿直徑不大于50mm，己知使用的材料的許用彎曲應(yīng)力120N/mm2，許用扭轉(zhuǎn)剪切應(yīng)力80N/mm2，許用撓度[f]=0.1mm，彈性模量E=200000N/mm2，要求設(shè)計的銑刀桿用料最省，試確定銑刀桿的結(jié)構(gòu)參數(shù)。

2．優(yōu)化模型的建立

2.1模型的基本假設(shè)及約定

銑刀桿種類繁多，但刀桿主要的設(shè)計參數(shù)應(yīng)為刀桿的桿長L(mm)和刀桿的直徑D(mm)，其他的輔助參數(shù)可由已得的L、D及使用的銑刀等條件查相關(guān)手冊得到。設(shè)計變量，寫成數(shù)學(xué)表達式：

2.2模型的建立

2.2.1確定目標(biāo)函數(shù)

一般的刀桿設(shè)計都是在滿足強度、剛度等的條件下，以求用料最省，即體積最小。所以，此處選刀桿體積最小為目標(biāo)函數(shù)：

2.2.2確定約束條件

（1）一般來講，為保證工件與銑刀的相對位置，刀桿的長度應(yīng)大于等于某個極限值 （mm），而刀桿直徑應(yīng)小于等于某個極限值 （mm）：

（2）刀桿在承受彎曲和扭轉(zhuǎn)時應(yīng)具有一定的強度，才能保證加工的安全性，彎曲應(yīng)力應(yīng)不大于材料許用應(yīng)力，剪切應(yīng)力不大于材料剪切應(yīng)力：

（3）銑刀桿的剛度影響被加工工件表面質(zhì)量等，所以，為了保證加工質(zhì)量，刀桿在承受載荷時發(fā)生的撓度f應(yīng)不大于許用撓度[f]：

其中 為不同的物理模型中的最大彎曲應(yīng)力與L，D的函數(shù)， 為不同的物理模型中的最大剪切應(yīng)力與L，D的函數(shù)， 為不同的物理模型中的最大撓度與L，D的函數(shù)。

2.2.3進行模型的抽象化處理

將上一步得到的目標(biāo)函數(shù)以及約束條件進行整理，令X 為二元矩陣，令D=x1，L=x2,可得優(yōu)化選取銑刀桿結(jié)構(gòu)參數(shù)的數(shù)學(xué)模型為：

注：對于不同的銑削形式，有不同的刀桿，也就對應(yīng)不同的力學(xué)模型，但總可以用以上模型描述。

3．基于MATLAB優(yōu)化工具箱進行模型求解

根據(jù)上面建立的數(shù)學(xué)模型對優(yōu)化設(shè)計實例進行處理可得：

首先建立目標(biāo)函數(shù)“fun.m”文件，以下是 fun 函數(shù)的語句：

%目標(biāo)函數(shù)

function f=fun(x)

f=(pi*x(2)*x(1)*x(1)/4)

注：因MATLAB中不能識別下標(biāo)值，故上述語句中的x（1）、x（2）即x1、x2分別代表刀桿直徑，刀桿長度。

接下來建立約束條件“mycon.m”文件，以下是“mycon.m”文件里包含的語句：

function[c,ceq]=mycon(x)

c(1)=[100000*x(2)-12*x(1)^3];

c(2)=[50000-8*x(1)^3];

c(3)=[34*x(2)^3-10*x(1)^4];

ceq=[];

注：上述語句中的c 矩陣為求解模型中的非線性約束。

最后建立求解文件“jie.m”（也可以不用建立求解文件，直接在MATLAB 的 command 窗口輸入以下語句即可），以下是 jie.m 中的語句：

x0=[50,100];

options=optimset('Largescale','on','TolCon',1e-8,'TolX',1e-8);

lb=[0,0];

A=[1,0;0,-1;-1,0];

b=[50;-100;-18.42];

Aeq=[];

beq=[];

[x,fval]=fmincon (@fun,x0,A,b,Aeq,beq,lb,[],@mycon,options)

注：因為刀桿直徑，刀桿長度不能為負值，因此 x1、x2（亦即語句中的 x（1）、x（2））的下限值向量 lb=［0，0，0］；矩陣A、b 為不等式線性約束 g1、g2的系數(shù)矩陣，矩陣 Aeq、beq 為等式線性約束的系數(shù)矩陣。通過以上的求解得x1= 42.9408mm，x2= 100.000mm，V= 1.4482e+005mm3。

表1 優(yōu)化值與經(jīng)驗值對比結(jié)果

比較項目 D L 目標(biāo)函數(shù) 性能約束

/mm /mm V(mm3)

經(jīng)驗值 45 120 1.9085e+005 11.6872 -74.5130 0.0433

優(yōu)化值 42.9408 100 1.4482e+005 0.0000 -73.2314 -0.0058

4．結(jié)論

比較優(yōu)化設(shè)計得到的結(jié)果，可以看出42.9408和100這兩個數(shù)值已經(jīng)能滿足設(shè)計要求，所以沒有必要再進行放大了。傳統(tǒng)設(shè)計法因為放大而導(dǎo)致了材料的浪費，而利用本文所述的優(yōu)化設(shè)計方法則可以很好地在經(jīng)濟性和安全性中求得一個最優(yōu)解。應(yīng)用 MATLAB優(yōu)化設(shè)計的方法，設(shè)計者只需要通過一次計算，便可得到可靠的設(shè)計結(jié)果，因此避免了傳統(tǒng)設(shè)計方法因設(shè)計者經(jīng)驗不足而造成的過度放大，提高了效率，縮短了設(shè)計周期。設(shè)計實例的正確性和合理性充分說明該模型的正確性，而本文運用MATLAB優(yōu)化工具的求解方法亦不失為一種既簡單又科學(xué)的方法。所以，本文建立的數(shù)學(xué)模型和解決方法是實用的。

5.參考文獻

[1] 羅佑新、郭惠聽、張龍庭等，機械零件的穩(wěn)健可靠性設(shè)計[J]，農(nóng)業(yè)機械學(xué)報，2002

[2]李萬祥，工程優(yōu)化設(shè)計與 MATLAB 實現(xiàn)，北京：清華大學(xué)出版社，2010

[3]吳道全、石琳，關(guān)于圓柱形球頭立銑刀結(jié)構(gòu)的研究[J]，貴州工學(xué)院學(xué)報，1996

[4] 張曉紅，系列化可轉(zhuǎn)位銑刀參數(shù)化設(shè)計系統(tǒng)的研究[J]，煤炭技術(shù)，2006