立即咨詢(xún)

基于敏感度分析的機床關(guān)鍵性幾何誤差源識別方法新算法

來(lái)源:互聯(lián)網(wǎng)   作者:機床世界    發(fā)表時(shí)間:2019-04-02    瀏覽量:7724



0 前言


隨著(zhù)航空航天、軍工、船舶、汽車(chē)等行業(yè)對精密零件加工的要求越來(lái)越高,機床的精度性能顯得更加重要,精密數控機床的設計制造也以高精度、高智能化、高效率為主要目標。精度設計是保證和提高機床精度的重要一環(huán),傳統的機床設計中,主要依靠經(jīng)驗的方法設計機床各部件的公差等級。由于各環(huán)節誤差對機床整體精度的影響程度不同,精度控制實(shí)現的難易程度也不同。因此,為節約機床的設計及制造成本,分析影響機床加工精度的各主要因素,建立機床的精度模型并進(jìn)行誤差敏感度分析顯得尤其重要。精度設計對精密數控機床的設計與制造具有重要的實(shí)際應用意義。


影響機床整體加工精度的各類(lèi)誤差主要有機床零部件的幾何誤差、熱誤差、載荷誤差、伺服誤差和插補誤差等,而幾何誤差對加工精度的影響又是精度設計的主要研究?jì)热?。機床的綜合幾何誤差體現在各零部件上,最終均將反映在被加工工件的加工誤差上,而通過(guò)準確建立機床的精度模型,可以反映出整體精度與各零部件的精度關(guān)系。國內外許多學(xué)者對數控機床空間誤差建模方法進(jìn)行了較廣泛而深入的研究,先后出現了幾何建模法、誤差矩陣法、二次關(guān)系模型法、機構學(xué)建模法、剛體運動(dòng)學(xué)法和多體系統理論法等。多體系統是一般機械系統最為全面的完整抽象、高度概括和有效描述,是分析和研究機械系統的最優(yōu)模型形式。目前國內外眾多學(xué)者已應用多體系統理論對機床進(jìn)行精度建模與分析。POTT等通過(guò)一種基于簡(jiǎn)化矢力的方法分析了并聯(lián)機構運動(dòng)靈敏度,并通過(guò)分析6自由度并聯(lián)機構Linapod,驗證了方法的正確性。王冰采用綜合誤差敏感度、絕對誤差敏感度和誤差方向敏感度這三個(gè)誤差評價(jià)指標,將各個(gè)機床誤差影響因素對機床終端誤差的影響進(jìn)行了分析。李小雷等對刻楦機進(jìn)行了誤差參數靈敏度分析。程剛等針對以3-RPS并聯(lián)機構可達空間內各誤差來(lái)源,提出了一種通過(guò)歸一化描述靈敏度的新方法,并建立基于數學(xué)統計意義的靈敏度系數數學(xué)模型。黃強等以滾齒機YK3610為對象,介紹基于多體系統理論和齊次坐標變換的機床誤差建模方法,并對機床敏感誤差辨識方法、步驟和關(guān)鍵點(diǎn)進(jìn)行闡述。但是,通過(guò)分析可以看出,目前針對機床誤差進(jìn)行敏感度分析的文獻還較少,對機床誤差溯源的研究深度不夠。


本文是在國家重大科技專(zhuān)項“精密立、臥式加工中心研究”和國家自然科學(xué)基金項目“基于模糊穩健設計的機床動(dòng)態(tài)精度分析與反演方法研究”的課題背景下,對精密加工中心進(jìn)行精度設計。為了有效地識別出對機床加工精度及空間誤差影響較大的幾何誤差參數,首先基于多體系統理論的方法,分析各部件幾何誤差對機床精度的影響,建立精密臥式加工中心的精度模型;然后提出了基于矩陣微分法的四軸機床誤差敏感度分析數學(xué)方法,根據模型計算精密臥式加工中心的誤差敏感度和敏感度系數,從而進(jìn)行誤差溯源。



1 基于多體系統理論的精度建模


根據精密臥式加工中心的結構與運動(dòng)特點(diǎn),將機床抽象成多體系統,用拓撲結構和低序體陣列描述機床各部分的關(guān)聯(lián)性;在多體系統中建立廣義坐標系,用齊次變換矩陣把機床各部件的誤差量之間的耦合關(guān)系表述出來(lái);推導出機床中兩相鄰體之間相對運動(dòng)的特征矩陣和運動(dòng)方程,從而建立加工中心的精度模型。


1.1 廣義坐標系的設置及機床的特征矩陣


精密臥式加工中心是一臺四軸數控機床,具有高剛度、高精度、高速度的特點(diǎn),根據設計要求,其技術(shù)指標如表1所示,其刀柄型號為BT50。


加工中心的結構簡(jiǎn)圖如圖1所示,它由床身,X、Y、Z三軸運動(dòng)部件,主軸箱(刀具)和繞B軸旋轉的工作臺組成。



根據多體系統理論的方法,將機床各組成部件分成相應的“體”。并按床身—X軸運動(dòng)部件—Y軸運動(dòng)部件—主軸箱(刀具)分支,和床身—Z軸運動(dòng)部件—工作臺—工件分支分別對其編號,圖2為機床的拓撲結構,表2是其低序體陣列,表2中Ln(j)是體Bj的n階低序體陣。



一般來(lái)說(shuō),機床的每一個(gè)單元部件在處于靜止或運動(dòng)狀態(tài)時(shí)都存在空間6個(gè)自由度方向的基本誤差,通過(guò)分析精密臥式加工中心的結構及其運動(dòng)關(guān)系可知,加上工件的裝夾誤差,其中影響機床加工精度的幾何誤差主要有35項,這部分幾何誤差在機床總誤差中占了較大的比重,表3列出了機床各部件的幾何誤差(部分)。



多體系統中各體之間的位置和運動(dòng)關(guān)系可以用相應的坐標系的位姿變換來(lái)表示,為了方便機床的精度建模,需要對坐標系進(jìn)行特殊設置處理。設置如下:①在床身(B0)和所有運動(dòng)部件(Bj)上均建立起與其固定聯(lián)接的右手笛卡兒三維坐標系,這些坐標系的集合成為廣義坐標系(又稱(chēng)參考坐標系),各體坐標系稱(chēng)為子坐標系(又稱(chēng)動(dòng)坐標系)。每個(gè)坐標系的3個(gè)正交基按右手定則分別為Xj、Yj、Zj軸;②廣義坐標系各元素Xj、Yj、Zj軸分別對應平行;③X軸部件、Y軸部件(主軸箱)、Z軸部件、工作臺的體運動(dòng)參考系與其對應相鄰低序體的體坐標系重合;④刀具坐標系原點(diǎn)與主軸端面中心重合;⑤工件上設定工件坐標系。


根據加工中心的結構和各部件之間的運動(dòng)關(guān)系,可建立各相鄰體間的變換特征矩陣??紤]到加工中心運動(dòng)部件的一些靜止誤差或者非運動(dòng)部件的運動(dòng)誤差相對很小,可令其誤差特征矩陣為單陣I4×4。表4為精密臥式加工中心各部件的特征矩陣T。



表4中帶下標p的量表示靜止狀態(tài)特征矩陣,帶下標s的量表示運動(dòng)狀態(tài)的特征矩陣。


式中,x、y、z分別為X軸部件、Y軸部件、Z軸部件的位移,α、β、γ分別為軸X、Y、Z的相位角。


1.2機床精度模型


設刀具成形點(diǎn)在刀具坐標系內的坐標為



工件上成形點(diǎn)在工件坐標系中的坐標為



理想運動(dòng)條件下有
式中



那么,機床在無(wú)誤差情況(理想情況)下,刀具成形點(diǎn)在工件坐標系內的理想成形函數為



在機械加工中,機床加工精度最終是由機床上刀具成形點(diǎn)與工件成形點(diǎn)之間的相對位移誤差決定的。在實(shí)際加工過(guò)程中,刀具成型點(diǎn)的實(shí)際位置不可避免地會(huì )偏離理想位置,從而產(chǎn)生空間位置誤差。實(shí)際成形點(diǎn)與理想刀具成形點(diǎn)的綜合空間位置誤差為



式中E——機床的空間誤差矢量


Ex,Ey,Ez——X軸、Y軸、Z軸方向上的誤差式(5)就是精密臥式加工中心的綜合精度模型,它由機床各部件的幾何誤差組成,根據此精度模型,基于敏感度分析可以識別對機床加工精度具有重要影響的關(guān)鍵性誤差源。


2關(guān)鍵誤差源識別


根據精密臥式加工中心的精度模型,利用矩陣微分法建立四軸機床誤差敏感度分析的數學(xué)模型。通過(guò)計算各個(gè)部件單元幾何誤差的敏感度,比較各個(gè)誤差元素對總的空間誤差的影響程度,最終識別出了影響機床加工精度的關(guān)鍵性誤差源,從而為合理經(jīng)濟地提高機床的精度提供重要的理論依據。


2.1機床空間誤差敏感度分析模型的建立


根據式(5),可以建立一般的四軸數控機床的誤差計算模型,該模型可以顯性地表示為



式中G——n個(gè)機床各零部件幾何誤差組成的誤差矢量
G=(△ei,△e2,...,△en)T
△ei——機床部件的幾何誤差,i=1,2,...
Pw——工件上成形點(diǎn)在工件坐標系中的坐標矢量
Pw=(pwxpwypwz1)T
U——機床各運動(dòng)軸的位置矢量
U=(x,y,z,B)T
Uw——工件位置坐標矢量
Uw=(xw,ywzw,1)T
Ut——刀具位置坐標矢量
Ut=(xt,ytt,zt1)T
根據式(5)的形式,可知F是變量G、Pw、U、Uw、Ut的連續可微函數,因此式(6)按一階泰勒級數展開(kāi),可以表達為



式中,ΔG、ΔPw、ΔU、ΔUw、ΔUt分別為G、Pw、U、Uw、Ut在其理想值處的微小波動(dòng)量,去掉高階項,可以得到



在實(shí)際的分析過(guò)程中,Pw、U、Uw、Ut可設為定值,即在固定的工件、刀具安裝位置下,模擬加工相同的工件軌跡,分析由機床各部件的幾何誤差變化時(shí)對總的空間誤差影響程度的大小。因此式(8)簡(jiǎn)化為



進(jìn)而可以得到機床加工精度的敏感度分析模型為



式中,S=δF/δG為雅可比矩陣,稱(chēng)為誤差敏感度矩陣,并具有如下的形式



至此,用矩陣微分的方法,得到了一般四軸數控機床的誤差敏感度矩陣的表達式。



2.2零部件幾何誤差參數的檢定


本論文以圖1所示的精密臥式加工中心為研究對象,以圖3為檢驗加工對象,利用QC20-W球桿儀和APIXD-6D雙頻激光干涉儀對各幾何誤差進(jìn)行測試檢定,如圖4所示。測得精密臥式加工中心的35項幾何誤差,考慮到篇幅有限,表5僅列出了部分幾何誤差值。



將表5所示的檢測幾何誤差以及工件的安裝位置(xwd,ywd,zwd,1)T、刀具的安裝位置(xtd,ytd,ztd,1)T、各導軌的運動(dòng)位移x、y、z代入到式(5),就可以得到工件上每個(gè)加工位置(pwx,pwy,pwz,1)T的誤差。根據精密加工中心的幾何精度檢驗標準GB/T20957.7—2007《精密加工中心檢驗條件第7部分:精加工試件精度檢驗》,將加工試件安裝在工作臺中心,在x-y平面上切削一個(gè)直徑d=218mm的圓,如圖3所示。



在機床的工作空間區域中,各坐標軸處于中間位置的工作區域使用最頻繁,因此主要對機床中間工作區域進(jìn)行誤差敏感度的分析。根據精密臥式加工中心的設計結構、技術(shù)參數以及多體系統中廣義坐標系的設置,可設定工件坐標系原點(diǎn)在工作臺坐標系中的位置坐標為
(xwdywdzwd1)T=(500-1901181)T
刀具坐標系原點(diǎn)在主軸坐標系中的坐標為
(xtdytdztd1)T=(00-1901)T
工件上加工點(diǎn)的坐標取值范圍是
pwx=[-109,109]pwy=[553.5,771.5]pwz=-150
各個(gè)運動(dòng)軸的取值范圍是
x=[391,609]y=[363.5,581.5]z=118
工作臺轉角b=0并且有
pwx=x-500pwy=y+190
其中,根據典型試件加工軌跡的特點(diǎn),可得x、y的關(guān)系為
(x-500)2+(y-472.5)2=1092
本文中,取機床加工試件時(shí)運動(dòng)到最高點(diǎn)(x,y,z)=(500,581.5,118)的位置對機床進(jìn)行誤差
敏感度分析(以上各參數值的單位為mm)。


2.3關(guān)鍵性誤差源參數的識別


機床的加工精度受到零部件幾何誤差耦合而成的空間誤差的影響,每一項誤差對綜合加工精度影響程度大小不一樣,究竟哪些空間幾何誤差參數對機床總加工精度的影響較大,主要影響著(zhù)機床的加工精度,為了找出各個(gè)誤差元素對總的空間誤差的影響程度大小,需要對機床的精度模型進(jìn)行誤差敏感度分析。針對本文分析的精密臥式加工中心,其精度模型包含了機床各個(gè)運動(dòng)部件產(chǎn)生的幾何誤差共35項。



誤差敏感度反映的是各個(gè)誤差元素產(chǎn)生微小變化時(shí)總空間誤差的變化程度,根據式(5)的精度模型和式(11)的誤差敏感度矩陣,可以定義機床的空間誤差E對各個(gè)幾何誤差Δei的偏導數的絕對值為精密臥式加工中心的誤差敏感度



對于的誤差敏感度S△γz,由式(5)、(11),可得



將上面的誤差參數值代入到式(13),并進(jìn)行整理,計算得



進(jìn)一步可計算得到△γz對應的誤差敏感度為



用同樣的計算方法,可以計算出精密臥式加工中心的精度模型中各誤差元素的敏感度。表6列出了加工中心各誤差元素對應的空間誤差E的表達式(部分),而從E的表達式中很容易得到相應的Si。


由表4可知,精度模型中的幾何誤差元素值按表3的數值取,并按試件的加工設定坐標變量后,總誤差E與各項誤差元素Δei的關(guān)系為
E=ai+ki△ei(14)
在分別對Em(m=x,y,z)進(jìn)行分析時(shí),式(14)中系數ki的絕對值越大,表明誤差元素Δei相對應的敏感度Si越大,其對空間誤差Em的變化影響也越大。為了更好地識別和分析關(guān)鍵性誤差源,將各誤差源參數的敏感度系數進(jìn)行歸一化處理,定義



式中,smi為Δei對應的誤差敏感度系數,所以敏感度系數之和為1。


對各幾何誤差對應的敏感度系數進(jìn)行分析可知,誤差敏感度系數越大,表明總誤差E對誤差元素Δei的變化敏感度越大,也即該幾何誤差源是影響機床總加工精度的一個(gè)關(guān)鍵因素,其對應的機床零部件對機床整機的加工精度影響越大。其中對Ex影響較大的關(guān)鍵性誤差源有6項,其對應的誤差表達式、敏感度和敏感度系數如表7所示。



表7中,Δγx、Δγz、ΔγB、Δγxy、ΔγxB這幾項誤差為機床本身的設計誤差,Δγw為工件安裝時(shí)產(chǎn)生的誤差。圖5a為此六項誤差元素對應的誤差敏感度系數比較圖。從圖5中可以看出這6項誤差的敏感度系數之和為0.89,其他誤差對應的敏感度系數僅為0.11。通過(guò)計算分析,可見(jiàn)機床的x軸運動(dòng)部件、z軸運動(dòng)部件、工作臺B軸轉動(dòng)部件繞z軸的轉角誤差,x、y軸的垂直度誤差,x、B軸的垂直度誤差,以及工件安裝時(shí)繞z軸的轉角誤差在加工零件時(shí)對總的空間誤差在x軸方向的誤差分量Ex影響較大。



同理,可以識別出對Ey、Ez產(chǎn)生重要影響的關(guān)鍵性誤差源參數。如圖5b所示,對Ey變化敏感度大的誤差元素有Δγz、ΔγB、ΔγxB、Δγw,這4項誤差的敏感度系數之和為0.82,其他誤差對應的敏感度系數僅為0.18??梢?jiàn)機床的z軸運動(dòng)部件、工作臺B軸轉動(dòng)部件繞z軸的轉角誤差,x、B軸的垂直度誤差,以及工件安裝時(shí)繞x軸產(chǎn)生的轉角誤差在機床加工過(guò)程中對總的空間誤差在y軸方向的誤差分量Ey影響較大。


如圖5c所示,對Ez變化敏感度大的誤差元素,有Δαx、Δαz、ΔαB、Δαw、Δαyz、ΔαzB、Δβz、ΔβB、Δβxz。這9項誤差的敏感度系數之和為0.9,其他誤差對應的敏感度系數僅為0.1??芍獧C床的x軸運動(dòng)部件繞x軸的轉角誤差,z軸運動(dòng)部件、工作臺B軸轉動(dòng)部件繞x和y軸的轉角誤差,y、z軸的垂直度誤差,x、z軸的垂直度誤差,z、B軸的垂直度誤差,以及工件安裝時(shí)繞x軸的轉角誤差在加工零件時(shí)對總的空間誤差在z軸方向的誤差分量Ez影響較大。


2.4靈敏度分析結果


通過(guò)對精密臥式加工中心進(jìn)行誤差靈敏度分析,可以得到以下幾點(diǎn)。


(1)機床的z軸運動(dòng)部件、工作臺B軸轉動(dòng)部件繞z軸的轉角誤差(Δγz:z軸導軌在垂直平面內平行度誤差的反映,ΔγB:工作臺沿x軸方向的水平度誤差的反映),以及x、B軸的垂直度誤差ΔγxB等這三項誤差對Ex、Ey的變化靈敏度都很大,因此,在機床的精度設計時(shí)應重點(diǎn)控制這幾個(gè)誤差的值。

(2)在機床的各個(gè)部件中,z軸導軌的直線(xiàn)度和平行度誤差、工作臺B軸轉動(dòng)部件產(chǎn)生的誤差對總的空間誤差E產(chǎn)生的影響較大。

(3)影響Ex的主要誤差有6項,影響Ey的主要誤差有4項,而影響Ez的主要誤差有9項。相比較而言,Ez受到的基本誤差元素的影響更復雜。


3結論


(1)精度建模與誤差敏感度分析是機床精度設計的核心內容之一。用多體系統理論的方法,綜合考慮機床各部件的幾何誤差,建立了精密臥式加工中心的精度模型。

(2)根據精度模型,提出了四軸數控機床的誤差敏感度分析方法,用矩陣微分法建立四軸數控機床的敏感度分析模型,并推導出精密臥式加工中心的敏感度分析模型。

(3)以誤差檢測試驗和加工典型標準試件為例對加工中心進(jìn)行誤差敏感度分析,得到z軸導軌在垂直平面內的平行度誤差、工作臺沿x軸方向的水平度誤差以及x、B軸的垂直度誤差ΔγxB等這幾項誤差因素對Ex、Ey的影響較關(guān)鍵,最終識別出了影響機床加工精度的關(guān)鍵性誤差,實(shí)現了誤差溯源,為精密數控機床的設計提供了重要的理論依據。