數控機床服役態全行程尺寸關系模型研究*

馬軍旭，金　濤，胡　敏，趙萬華

(西安交通大學 機械制造系統工程國家重點實驗室，西安　710054)

數控機床服役態全行程尺寸關系模型研究*

馬軍旭，金濤，胡敏，趙萬華

(西安交通大學 機械制造系統工程國家重點實驗室，西安710054)

摘要：重力、熱等造成機床在裝配時的精度發生變化，因此在機床精度設計時應考慮服役態因素的影響。文章提出一種考慮服役態影響的尺寸關系模型用于機床精度設計。首先對基礎件進行離散表征，利用尺寸鏈理論，建立了服役態全行程尺寸關系模型。其次給出了服役態下尺寸關系模型中各環的等效尺寸計算方法。最后以某型號機床Y軸尺寸為例，進行了服役態全行程尺寸關系模型實驗驗證。結果顯示：在全行程范圍內，封閉環尺寸實測值與利用服役態尺寸關系模型得到的計算值最大相差2.2μm，封閉環尺寸實測值與幾何態封閉環尺寸計算值最大相差3.8μm，服役態尺寸關系模型更接近于實際情況。

關鍵詞：服役態；全行程；離散表征；尺寸關系模型

0引言

數控機床幾何精度占數控機床精度的70%以上。通常通過設計制造(主動)和補償(被動)兩種措施保證機床幾何精度。在幾何精度設計時，使用尺寸鏈進行公差分析與設計。

從1978年Hillyard R.C.博士[1]首次提出計算機輔助確定零件的幾何形狀、尺寸和公差的概念開始，公差模型經歷了屬性模型[2]、漂移模型[3]、參數化模型[4]、運動學模型[5]、自由度模型[6]等時期。利用公差模型進行公差分析，提出了基于齊次轉換矩陣方法[7]、矢量環及直接線性化方法[8]、基于雅克比矩陣的虛擬接點法[9]、約束自由度模型[10]、矢量公差法[11]、Gap-Space方法[12]和擴展空間技術[13]等方法或技術。利用公差模型進行公差設計，建立了一系列的公差-成本模型[14-15]、公差優化算法[16-18]、公差綜合模型[19-20]等，出現了一系列的公差設計理論，如：矢量公差設計[21]、并行公差設計[22]、魯棒公差設計[23]和動態公差控制[24]等。

上述文章主要針對公差模型、公差分析和公差設計，利用多體系統理論、小位移旋量理論、齊次坐標變換方法、模擬退火算法、遺傳算法和模糊神經網絡等建立了不同的數學模型，研究了公差的建模方法，公差大小是否滿足幾何功能要求、是否滿足可裝配性、可制造性以及制造成本等CAT問題。

圖1　服役態與幾何態對比示意圖

如圖2所示，圖2a表示固定件是剛體，移動件在不同位置時，固定件尺寸、移動件尺寸和封閉環尺寸關系，此時A1=A3。圖2b表示固定件是柔體，移動件在不同位置時，三者間的尺寸關系，此時A1≠A3。

圖2　在全行程不同位置對比剛體與柔體變形示意圖

因此，數控機床在精度設計時，需要考慮全行程不同位置時服役態的影響，上述模型或方法不能很好地滿足這種情況。

本文的目的是建立用于公差設計的服役態尺寸關系模型。利用離散點尺寸描述全行程內固定件的尺寸變化，根據尺寸鏈理論建立全行程服役態尺寸關系模型。計算在全行程不同位置處，服役態因素造成的封閉環尺寸變化。

2服役態尺寸變化模型

2.1幾何要素的離散與表征

由于傳統的尺寸鏈中，直線和平面簡化成理想情況所造成的不足，建立服役態尺寸關系模型時，需要對零部件中與尺寸相關的幾何要素(線、面)重新進行表征。如圖3所示，把固定件沿X方向離散成若干份。離散原則如下：

(1)測量基準視為理想要素，沒有誤差。

(2)運動方向等距離離散。

在離散后的物體上放置一個坐標系XOY，稱為零件坐標系，任一位置的尺寸利用特征點坐標表示。點P的坐標為Pi(xi,yi)，測試基準與X軸的距離為y0，則在全行程內不同位置處，固定件尺寸為

Ly(xi)=yi+y0i=1…n

圖3　尺寸與坐標關系示意圖

2.2服役態尺寸關系模型

如圖4所示的數控機床直線運動軸部件，通常由固定件(床身等)、導軌滑塊系統(直線導軌或靜壓導軌燈)和移動部件(工作臺等)組成，在進給方向(Y方向)上由螺母絲杠或者直線電機驅動。

圖4　機床直線軸結構示意圖

設在Y方向全行程內分成n段，則工作臺在全行程有(n+1)個位置。在導軌1和導軌2上，分別有2(n+1)個位置與滑塊相對應。取2(n+1)個點的坐標代替床身的Z向尺寸。如圖5所示。

圖5　床身尺寸離散示意圖

根據尺寸鏈理論，在全行程任一位置，建立工作臺與地面(理想基準)的服役態尺寸關系模型，如圖6所示。則服役態全行程尺寸關系模型為：

AΣi=A1i+A2ii=1…n

(1)

圖6　服役態尺寸關系模型

2.3封閉環尺寸的計算

由服役態尺寸關系模型可知，組成環為固定件尺寸A1i和移動件尺寸A2i；封閉環為AΣi。

在全行程的任一位置，床身尺寸(固定件尺寸)由滑塊尺寸決定。

在Y軸不同位置處，固定件尺寸A1i的等效尺寸為：

(2)

A滑塊是滑塊位置對應的床身尺寸。由加工后尺寸和服役態變形共同決定，在全行程內隨滑塊位置變化。

本文僅以俯仰造成的尺寸變化為例等效移動件尺寸。在Y軸不同位置處，移動件尺寸A2i的等效尺寸為：

(3)

A2為移動件的理想尺寸；L為滑塊在Y向的間距，如圖7所示。

圖7　尺寸A2i等效示意圖

由公式(1)、(2)、(3)即可得到封閉環尺寸。

3實驗方法

本實驗的目的是驗證在全行程不同位置時，服役態因素(此時只考慮重力)造成的封閉環尺寸變化，以及本文所提出的方法有效性。

圖8呈現的是本文的實驗裝置。包括兩個部分：某型號機床和激光跟蹤儀。激光跟蹤儀Laser Tracker Leica TN90精度為3μm。

圖8　實驗測試圖

本文中，實驗分兩步進行。第一步，拆下工作臺，分別在兩個導軌上推動滑塊，建立固定坐標系。在Y方向全行程分成5段，在離散的不同位置處，利用激光跟蹤儀采集滑塊在固定坐標系下的坐標值。根據離散數量，兩根導軌共采集24個點。第二步，安裝工作臺。激光跟蹤儀保持不動，保證測試在同一個坐標系進行。在離散的6個位置處，以工作臺上一點代替工作臺表面，分別采集其在固定坐標系的坐標值。不論導軌滑塊坐標值測試點還是工作臺上測試點，每個測試點重復采集五次坐標值，各點的坐標值由五次平均得到。

利用游標卡尺(精度為0.01mm)分別在四個滑塊位置測試工作臺尺寸，每個位置測試五次，取平均值作為工作臺尺寸(移動件尺寸)。

由尺寸和坐標的關系可知，固定件尺寸A1i和封閉環尺寸AΣi由坐標測試值和坐標軸到尺寸測試面的距離決定。本文采用坐標測試值代替A1i和AΣi，相當于測試基準由固定件底面變為坐標系坐標軸，不影響計算結果的準確性。

測試時，環境溫度變化在24.2±0.2℃以內。

4結果和討論

表1所示為工作臺拆下后(幾何態)，在Y軸行程不同位置處，由公式(2)計算得到的床身尺寸等效坐標值；以及工作臺重力造成床身不同位置處變形后(服役態)，床身尺寸等效坐標值。由表可知，工作臺安裝到床身后，在工作臺重力作用下，Y軸不同位置處的床身尺寸不同。床身尺寸變化的最小值為2.5μm，最大值為2.9μm。

游標卡尺測試得到工作臺平均尺寸為95.81mm。由公式(3)計算得到工作臺尺寸等效坐標值，如表2所示。

由公式(1)計算分別得到幾何態和服役態下封閉環尺寸如表3所示。

由激光跟蹤儀測試得到的裝配后工作臺尺寸(封閉環尺寸)坐標值如表3所示。

圖9為封閉環尺寸實測值、服役態計算值和幾何態計算值的對比圖。

圖9　封閉環尺寸實測值、服役態計算值及幾何態計算值對比

Y軸位置(mm)0100200300400500床身坐標(mm)幾何態-17.8951-17.8938-17.8938-17.8915-17.8941-17.8982服役態-17.8892-17.8900-17.8931-17.8906-17.8925-17.8937床身尺寸變化量(μm)5.93.80.70.91.64.5

表2　Y軸不同位置處工作臺尺寸等效后的坐標值

表3　服役態、幾何態計算值與實測值對比表

在Y軸全行程不同位置處，封閉環尺寸是變化的。只考慮幾何狀態時，封閉環尺寸為77.9185mm，考慮服役態時，全行程內，封閉環尺寸最大差值為6.9μm，實測值最大差值為4.5μm。只考慮幾何狀態時，床身尺寸是由不同位置的測試值平均得到的一個數值。由于床身尺寸采用離散表征，服役態造成在全行程不同位置時，床身尺寸的變化是不同的。由此也會造成移動件的姿態變化，馮更新[25]等的研究證明了磨床中確實存在這種姿態變化。

在服役態和幾何態兩種情況下，封閉環尺寸計算值是不相同的。幾何態計算值和服役態計算值最大相差6.0μm，平均值相差3.1μm。這是由于重力作用造成的尺寸變化。文中只考慮了重力的影響，機床在使用過程中，其他因素(熱、裝配應力等)也會造成服役態尺寸變化。

測試值與服役態計算值最大相差2.2μm，各點偏差的平均值1.03μm；與幾何態計算值最大相差3.8μm，各點偏差的平均值2.78μm。這說明考慮服役態影響的計算結果比幾何態的計算結果更符合實際情況。

5結論

通過理論分析和實驗測試得到如下結論：

(1)本文建立了服役態全行程尺寸關系模型。利用此模型計算出的服役態全行程尺寸變化量，是進行尺寸設計和精度設計的依據，后續可以依此為基礎，研究公差分配和誤差補償問題。

(2)本文利用所建立的模型計算了工作臺在重力作用下造成裝配后封閉環尺寸的變化量。結果表明考慮服役態影響的尺寸計算比純幾何態的尺寸計算更符合實際情況。

[參考文獻]

[1] Hillyard R C. Dimensions and Tolerances in Shape Design [D].London: University of Cambridge, 1978.

[2] Tsai J C, Wang W W. Development of a Computer-Aided Tolerancing System in CAD Environment[M].Global Consistency of Tolerances. Springer Netherlands, 1999: 47-54.

[3] Requicha A A G. Toward a theory of geometric tolerancing[J]. The International Journal of Robotics Research, 1983, 2(4): 45-60.

[4] Krishnan K K, Eyada O K, Ong J B. Modeling of manufacturing processes characteristics for automated tolerance analysis[J]. INTERNATIONAL JOURNAL OF INDUSTRIAL ENGINEERING-APPLICATIONS AND PRACTICE, 1997, 4(3): 187-196.

[5] Desrochers A, Rivière A. A matrix approach to the representation of tolerance zones and clearances[J]. The International Journal of Advanced Manufacturing Technology, 1997, 13(9): 630-636.

[6] Kandikjan T,Shah J J, Davidson J K. A mechanism for validating dimensioning and tolerancing and tolerancing schemes in CAD systems[J]. Computer-Aided Design,2001, 33(10):721-737.

[7] Whitney D E,Gibert O L. Representation of geometric variations using matrix transforms for statistical tolerance analysis in assemblies[J].IEEE International Conference on Robotics and Automation.1993,93(2):314-321.

[8] Chase K W, Magleby S P, Glancy C G. A comprehensive system for computer-aided tolerance analysis of 2-D and 3-D mechanical assemblies[M]//Geometric Design Tolerancing: Theories, Standards and Applications. Springer US, 1998: 294-307.

[9] Laperriere L, Lafond P.Tolerance analysis and synthesis using virtual joints[J].In Proceedings of 6th CIRP Seminars on Computer Aided Tolerancing.Enschede,Nether lands,1999:452-458.

[10] 胡潔，熊光楞，吳昭同.基于變動幾何約束網格的公差設計研究.機械工程學報.2003,39(5):20-26.

[11] Weber C,Britten W,Thome O. Conversion of geometrical tolerances into vectorial tolerance representations-A major step towards computer aided tolerancing[C].In International Design Conference, Dubrovnik, 1998.

[12] Zou Z, Morse E P. A gap-based approach to capture fitting conditions for mechanical assembl[J]y.Computer-Aided Design,2004, 36(8):691-700.

[13] Lee S,Yi C. Tolerance analysis for multi-chain assemblies with sequence and functionality constraints[J].In Proceedings of IEEE1997,International Conference on robotics and automation.Albuqyerque,New Mexico,1997.

[14] Bjorke O.Computer-Aided Tolerancing[M].2nded. NewYork:ASME Press,1989.

[15] Dong Z, Hu W, Xue D. New production cost-tolerance models for tolerance synthesis[J].Journal of Engineering for Industry. 1994, 116(3):199-206.

[16] Parkinson D B. Assessment and optimization of dimensional tolerances. Computer-Aided Design. 1985,17(4):191-199.

[17] Lee W J, Johnson G E. Optimal tolerance allotment using a genetic algorithm and truncated Monte-Carlo Simulation[J]. Computer-Aided Design. 1993,25:601-611.

[18] Chase K W,Greenwood W H,Loosti B G,et al. Least cost tolerance allocation for mechanical assemblies with automated process selection[J].Manufacturing Review. 1990.3:49-59.

[19] Pramanik N. A generic framework for design evolution:From function to form mapping to synthesis of tolerances[D]. Syracuse: Syracuse University,2003.

[20] Dantan J Y,Mathieu L,Ballu A,et al.Tolerance synthesis:quantifier notion and virtual boundary[J]. Computer-Aided Design.2005,37(2):231-240.

[21] Humienny Z,Thome O,Weber C.Vectorial dimensioning and tolerancing in computer aided tolerancing[C]. Poland:Proceedings of International Conference Mechatronics,2000, 47-450.

[22] Peng Heping,Jiang Xiangqian,Liu Xiaojun.Concurrent optimal allocation of design and process tolerance for mechanical assemblies with interelated dimension chains[J].International Journal of Journal of Production Research,2008,46(24):6963-6979.

[23] Jeang A.Robust tolerance design by response surface methodology[J].International Journal of Advanced Manufacturing Technology,1999,15:399-403.

[24] McGarvey R G,Lehtihet E A,Castillo E D,et al.On the frequency and location of set point adjustments in sequential tolerance control[J].Internatioanal Journal of Production Research,2001,39(12):2659-2674.

[25] 馮更新.淺析機床動態下的裝配尺寸鏈問題[J].磨床與磨削，1994(2):48-50.

(編輯趙蓉)

Research on Dimension Relationship Model Under Service Conditions in Whole Stroke of CNC Machine Tool

MA Jun-xu, JIN Tao, HU Min, ZHAO Wan-hua

(State Key Laboratory for Manufacturing System Engineering，Xi’an Jiaotong University，Xi’an 710054，China)

Abstract：The precision of CNC machine tools changes in the assembly process due to gravity and thermal effect. So, the effect of the service state on the precision should be considered during the design of CNC machine tools. This paper presents a dimension relationship model under service conditions, which is used for the precision design of CNC machine tools. Firstly, the geometrical feature of fixed part is discretized as a set of coordinate points. Then a dimension relationship model under service conditions in the whole stroke of operation is constructed based on the dimension chain theory. Secondly, the calculation method of the dimensions in the model is developed. Finally, Y-axis dimension of a CNC machine tool is tested to verify the dimension relationship model under service conditions. The results show that at any position in the whole stroke of operation, the biggest different is 2.2μm between the tested values and the calculated values; and the biggest different is 3.8μm between the tested values and the calculated values.

Key words：service conditions; whole stroke; discrete; dimension relationship model

文章編號：1001-2265(2016)06-0092-04

DOI:10.13462/j.cnki.mmtamt.2016.06.024

收稿日期：2015-06-16

*基金項目：國家自然科學基金資助項目(51235009)；國家“高檔數控機床與基礎制造裝備”科技重大專項(2011ZX04016-101，2012ZX04005011)

作者簡介：馬軍旭(1981—)，男，河南周口人，西安交通大學博士研究生，研究方向為數控機床精度保持性，(E-mail)mx810606@sohu.com；通訊作者：趙萬華(1965—)，男，遼寧本溪人，西安交通大學教授，博士研究生導師，研究方向為高速、超高速以及精密數控機床集成理論和機電耦合，(E-mail)whzhao@mail.xjtu.edu.cn。

中圖分類號：TH162;TG506

文獻標識碼：A