常方圓，方 峻，羅少敏

（南京理工大學 機械工程學院，江蘇 南京 210094）

1 引言

隨著三維數(shù)字化與信息技術(shù)的發(fā)展，國外各大先進制造企業(yè)已全面采用數(shù)字化設(shè)計制造和三維仿真技術(shù)的應(yīng)用。而我國兵器行業(yè)目前的制造裝配模式還處于數(shù)字化的起步階段，并沒有形成完整的數(shù)字化集成應(yīng)用[1]。對于傳統(tǒng)的、僅適于手工設(shè)計環(huán)境的尺寸公差分析已經(jīng)不能適應(yīng)現(xiàn)代產(chǎn)品設(shè)計、制造、仿真分析等要求。尤其是在CAD/CAM已高度普遍化的今天，它已成為制約三維數(shù)字化集成應(yīng)用技術(shù)發(fā)展的瓶頸。

公差分析（Tolerance Analysis）是指在公差制定后根據(jù)誤差的傳遞路線驗證設(shè)計公差是否滿足裝配功能和精度要求。而早期公差分析是以二維圖紙為基礎(chǔ)依賴于尺寸鏈和公差帶圖，通過傳統(tǒng)的尺寸鏈計算分析方法（如極值法、統(tǒng)計法）進行設(shè)計，同時常規(guī)的計算手段與計算方法工作量大，計算效率低并且校核困難。對此國內(nèi)外學者進行了大量的研究，文獻[2]首次提出通過數(shù)學方程式來描述零件的幾何特征，并以此來進行零件的尺寸和公差設(shè)計。文獻[3]在空間運動分析法的基礎(chǔ)上，提出了在幾何變動中利用Monte Carlo進行模擬的公差分析方法。文獻[4]對比分析了利用矢量環(huán)與雅克比矩陣進行三維公差分析求解的不同。隨著CAT技術(shù)在我國迅速推廣及應(yīng)用，國內(nèi)很多學者對此也進行了研究[5-7]。

重點以M16自動步槍為例進行裝配精度分析方法的研究，將特征要素離散化處理，利用矢量坐標矩陣描述物體特征要素在空間中的位置變動，物體的裝配偏差累積采用齊次特征矩陣變換來表示。通過在公差域范圍內(nèi)的隨機取值，模擬配合表面的誤差變動。結(jié)合小位移旋量理論，利用旋量參數(shù)與公差邊界的數(shù)學關(guān)系作為約束條件，建立符合自動武器結(jié)構(gòu)特點的三維公差分析模型，進行裝配精度分析，獲得了間隙的變動范圍。

2 理論分析

裝配精度分析也稱為公差累積分析，目的在于分析單個零件的變動范圍對裝配精度的影響，通過建立裝配公差模型，利用蒙特卡洛仿真分析零件的公差累積，從而得到裝配間隙變動范圍。

小位移旋量（Small Displacement Torsor，SDT）通過六個運動分量組成的矢量來表示剛體產(chǎn)生的微小位移。公差可以看作是幾何要素的變動約束，可以轉(zhuǎn)化為相關(guān)特征點在公差域內(nèi)的變動。而物體的點、線、面等基本幾何要素在公差域中的變動均可以通過旋量參數(shù)來表示。文獻[8]在1996年首次將SDT引入到公差領(lǐng)域，提出了基于SDT的公差數(shù)學表示方法。零件的幾何要素的變動可以用三個平移矢量D=［Δx，Δy，Δz］T和三個旋轉(zhuǎn)矢量R=［Δα，Δβ，Δγ］T來表示，則 SDT表示為τ=［R，D］＝（Δα，Δβ，Δγ；Δx，Δy，Δz）T。

幾何特征的變動可理解為名義要素上的點在公差帶約束條件下，經(jīng)過一定的空間運動后所形成的偏移。由此，特征點的空間運動可以利用機器人運動學中坐標變換理論，通過D-H矩陣的形式來表示

用上式的表述方法描述配合特征的變動形式，采用隨機數(shù)模擬實際特征要素的變動，通過矢量元素之間的約束方程，保證變動位于公差域內(nèi)。使每個裝配關(guān)系中配合面的變動符合實際裝配。

2.1 幾何特征公差域參數(shù)化建模

理想平面是通過幾何要素（點、線、面等）來描述的，在兩個理想平面間的任何相對位置均可以用小位移旋量的形式來表示。在此，考慮一個名義平面S0和一個擬合加工表面S1的變動關(guān)系。構(gòu)造一個平面S0，平面上任一點N（a，b，c）處的法線單位矢量為ZN，且為實體外法向矢量方向，如圖1所示。在S0附近設(shè)置一個間距為T的兩平行表面組成公差帶，使S0位于這兩個平行表面的約束范圍內(nèi)。擬合表面S1，只能在距離基準平面S0-TL～S0+TU間變動，這個空間區(qū)域是由S0的正負兩個偏移構(gòu)造而成。

圖1 有限個點擬合變動平面Fig.1 Finite Point Method for Fitting Plane

擬合表面S1偏移范圍如下：?N∈S0

任一平面能離散化為n個點Ni，這n個點在公差范圍內(nèi)變動可用上述方程表示。通過在名義平面S0上沿著法線矢量ZN方向構(gòu)造n個點，即n個點Ni的映射，可確定S0到Si的每一對應(yīng)點之間的距離。通過研究極限點在SDT的作用下產(chǎn)生的偏差必須限制在公差域，可以消除某些自由參數(shù)。

2.1.1 平面特征

在閉鎖機構(gòu)裝配中，配合特征多為圓形表面接觸或軸孔配合。故先建立圓形平面特征變動模型：在圓形平面圓周上任取一點 N，與 x 軸正方向夾角為 θ，則點坐標值為（rcosθ，rsinθ，0）。點N在公差域內(nèi)變動，通過齊次坐標變換如下：

故 Z 坐標軸方向的變化量為 δz=Δz-Δβ·rcosθ+Δα·rsinθ，由于δz的變動范圍必須限制在公差域［-TL，+TU］區(qū)間內(nèi)，則有：?θ∈［0，2π）

2.1.2 圓柱面特征

孔-軸配合公差可以看作是圓柱面的方位變動，即建立圓柱面的特征變動模型：設(shè)名義圓柱表面上任意一點M的坐標為（rcosθ，rsinθ，z）在旋量矩陣Δτ作用下M點轉(zhuǎn)換到新位置點M′：

圓柱度偏移的極限情況，如圖2所示。在XOZ平面內(nèi)圓柱度公差帶繞Y軸轉(zhuǎn)動的傾角β取得極大值。

圖2 圓柱度公差帶極限位置Fig.2 Limit Position of Cylindrical Tolerance Zone

圓柱度公差帶處于極限位置時，其邊界與XOZ平面的交線Lab的方程為：X=R+Δz+z·Δβmax

由于是微小位移的偏轉(zhuǎn)，故cosβmax≈1代入上式可得：

2.1.3 圓錐面特征

槍管與彈型樣柱之間的接觸面為圓錐面，其錐面的配合對裝配偏差累積有較大影響，為了滿足裝配要求需建立斜面圓跳動公差模型。斜面圓跳動的公差帶是在與基準同軸的任一半頂角的圓錐面上，距離為t的兩圓之間的區(qū)域，如圖3所示。

圖3 斜面圓跳動公差域Fig.3 Circular Run-Out Tolerance Region

設(shè)名義圓錐面上任意一點 P 的坐標為（tcosθ，tsinθ，t）在旋量矩陣Δτ作用下P點轉(zhuǎn)換到新位置點P′：

故 Z 坐標軸方向的變化量為 δz=Δz-tcosθ·Δβ+tsinθ·Δα，由于δz的變動范圍必須限制在距離為t的兩圓之間的區(qū)域內(nèi)，則有：

3 實例分析

3.1 閉鎖機構(gòu)結(jié)構(gòu)分析

閉鎖機構(gòu)是武器發(fā)射時關(guān)閉彈膛，承受彈殼底平面對槍機的軸向作用力并防止后逸的機構(gòu)，對自動武器工作的可靠性和射擊準確性都有重要的影響。M16自動步槍采用回轉(zhuǎn)式閉鎖機構(gòu)，閉鎖機構(gòu)由槍機、節(jié)套、機框和槍管的有關(guān)部件組成，如圖4、圖5所示。通過槍機上沿圓周均勻布置的閉鎖凸筍，作用在節(jié)套閉鎖槽內(nèi)的閉鎖支撐面上完成閉鎖。

圖4 閉鎖機構(gòu)結(jié)構(gòu)圖Fig.4 Structure Diagrams of Locking Mechanism

圖5 閉鎖機構(gòu)裝配爆炸圖Fig.5 Assembly Exploded Diagram of Locking Mechanism

表1 閉鎖機構(gòu)名義尺寸及其公差值Tab.1 Nominal Dimension and Tolerance Value of Locking Mechanism

表2 閉鎖機構(gòu)加工裝配要求Tab.2 Requirements for Processing and Assembling of Locking Mechanism

而在推彈進膛和閉鎖過程中，軸向間隙可能出現(xiàn)在不同部位。此時槍彈與彈膛定位面貼合，槍機與節(jié)套閉鎖支撐面貼合，在彈殼底平面與槍機彈底窩鏡面之間所形成的軸向間隙稱之為彈底間隙。為了保證武器在壽命期不發(fā)生斷殼故障，需要對產(chǎn)生間隙的各組成部分進行分析。根據(jù)加工圖紙，相關(guān)名義尺寸及其公差值，如表1所示。具體加工裝配要求，如表2所示。

3.2 分析計算

在公差分析中，幾何變動的誤差累積是按照運動學關(guān)系在裝配體中傳遞的，采用齊次坐標變換的方法來描述零件特征要素之間的幾何關(guān)系。首先以彈型樣柱底面的圓心為原點，以中心線為Z軸，建立固定坐標系。在每個接觸表面建立局部坐標系，幾何關(guān)系傳遞，如圖6所示。該裝配模型共4處誤差傳遞。從固定坐標系原點出發(fā)，經(jīng)過各個變動特征，再回到固定坐標點，形成一個封閉的回路。那么，裝配體中每個特征在固定N坐標系下的相對位置可以通過幾個齊次變換矩陣相乘得到：

將各個特征面的小位移旋量參數(shù)引入到旋量矩陣中，則有包含幾何變動信息的特征從局部坐標系向固定坐標系轉(zhuǎn)換，得到各參數(shù)的偏差累積模型：

式中：Ti—一個局部坐標系向另一個局部坐標系轉(zhuǎn)換的變換矩陣；Δτi—含有旋量參數(shù)的偏差傳遞矩陣；P—固定坐標系特征面上任一點的坐標矩陣；δP—沿三坐標軸的偏差累積。

圖6 閉鎖機構(gòu)裝配幾何特征傳遞關(guān)系圖Fig．6 Geometric Features of Locking Mechanism Transfer Relation Graph

蒙特卡洛模擬法是以概率統(tǒng)計為基礎(chǔ)，通過對隨機變量的統(tǒng)計實驗和隨機模擬來求解問題近似解。由于實際表面相對于名義表面的變動是無法預知的，具有隨機性，因此，公差分析中的旋量參數(shù)也是隨機量。而零件在實際加工的過程中通常會受眾多因素的影響，各種加工方法在隨機性因素的影響下所得加工尺寸的分散規(guī)律符合正態(tài)分布，其尺寸分散范圍為6σ。所以在一般情況下，加工方法的標準差σ與公差帶寬度T之間具有以下關(guān)系：

通過Matlab編程計算，在配合特征公差域內(nèi)產(chǎn)生隨機數(shù)，模擬零件特征要素的幾何變動，所取得的隨機數(shù)只有在滿足公差邊界條件與旋量參數(shù)約束方程的情況下才有效，否則重新取值直到滿足條件。進行10000次模擬，將每次分析結(jié)果存入分析矩陣，并進行統(tǒng)計分析。得出閉鎖間隙在基準坐標系中的空間誤差沿x，y，z軸坐標的誤差統(tǒng)計分布圖，如圖7所示。

由分析結(jié)果可知，閉鎖間隙沿三坐標軸的誤差變動具有一定的統(tǒng)計規(guī)律，近似服從正態(tài)分布，如圖8所示。且在x軸和y軸方向上變動較小，在有裝配功能要求的z軸方向上有較大的誤差變動，Δz變動范圍在（-0．06～0．16）mm 其中滿足裝配要求范圍的樣本占98．2%，與實際裝配情況相一致。由于沿x軸和y軸為實際配合表面的徑向方向，表現(xiàn)在同軸度的誤差變動。沿z軸為實際配合面的軸向方向，是全局誤差累積較大的地方，也是裝配要求重點保證的。

圖7 公差分析流程圖Fig．7 The Flow Chart of Tolerance Analysis

圖8 沿三坐標軸誤差統(tǒng)計分布圖Fig．8 Statistical Distribution of Deviation Along the Three Coordinate

3.3 軟件仿真分析

Vis VSA（Visualization Variation Simulation Analysis）軟件是基于三維模型，采用Monte Carlo模擬法進行公差分析，預測設(shè)計公差是否能夠滿足其關(guān)鍵的尺寸控制要求。分析步驟如下：（1）在UG中建立閉鎖機構(gòu)的三維模型，轉(zhuǎn)換成VSA軟件可識別的JT格式，導入分析軟件。定義裝配基準、功能特征、裝配關(guān)系和測量要求。（2）設(shè)定尺寸公差與形位公差等參數(shù)，將這些具有偏差的零件按照裝配順序組裝起來，在模擬多次裝配過程中，產(chǎn)品偏差的變化量被記錄下來，形成統(tǒng)計分布圖，如圖9所示。

圖9 公差分析模型Fig．9 Tolerance Analysis Model

本例采用10000個樣本點進行仿真模擬，得到閉鎖間隙沿Z軸方向的誤差統(tǒng)計分布圖，如圖10所示。當樣本點分布在±3σ范圍內(nèi)，我們認為結(jié)果是可信的。由分布圖可知間隙變動范圍為（-0．04～0．16）mm。

圖10 誤差統(tǒng)計分布圖Fig．10 Deviation Statistical Distribution

若只考慮線性尺寸公差可以得到極值法的計算結(jié)果：（-0．04～0．15）mm，將以上幾種方法進行分析對比，如表3所示。從結(jié)果可以看出，進行公差累積分析時，考慮形位公差和裝配間隙后計算出的變化范圍比極值法的范圍更大。這是由于形位誤差表現(xiàn)為配合表面形狀上的變化，在實際裝配中槍管和槍彈的配合是錐形面貼合，由于形狀誤差的存在，槍彈在彈膛中的實際接觸位置具有不確定性，從而對裝配精度造成影響。在線性尺寸鏈計算時考慮理想表面的接觸方式，忽略了幾何變動對接觸點的影響。而通過旋量矩陣法可以模擬特征面的實際變動，并同時分析沿三坐標軸的裝配誤差，相對于極值法來說真實反映了裝配的不確定性。同時，利用VisVSA公差仿真軟件分析可以快速、直觀、準確地對裝配精度進行預測。

表3 三種公差分析方法對比Tab.3 Comparisonof Three Kinds of Tolerance Analysis Method

4 結(jié)論

（1）基于三維模型的公差分析考慮了形位公差和裝配間隙等影響因素，提出了特征要素離散化處理的方法，解決了以往在建立約束條件時，因存在取點數(shù)量與位置不同而產(chǎn)生誤差的問題。（2）通過在公差范圍內(nèi)的隨機取值，模擬配合表面的誤差變動，在滿足公差邊界條件與旋量約束條件下進行裝配誤差累積分析，提出了基于三維模型的閉鎖機構(gòu)裝配精度預測方法。所得出的結(jié)果比線性尺寸鏈的分析結(jié)果較為嚴格，真實反映了加工的不確定性，更符合工程實際。

［1］羅少敏，王亞平，鄒衍．某型手槍擊發(fā)機構(gòu)裝配公差分析與信息化集成研究［J］．機械設(shè)計與制造，2014（28612）：132-136．（Luo Shao-min，Wang Ya-ping，Zou Yan．Research on assembly tolerance for firing mechanism of a pistoland solutions for information integration［J］．Machinery Design and Manufacture，2014（28612）：132-136．）

［2］HillyardR．C．Characterizing non-ideal shapes in termsof dimensionsan tolerances［J］．Computer Graphics，1983，12（3）：234-238．

［3］Turner，J．U．A．Feasibility space approach for automated tolerancing［J］．Journal of Engineering for Industry，1993，115（3）：325-341．

［4］MMarziale，WPolini．A review of two models for tolerance analysis of an assembly：vectorloopandmatrix［J］．InternationalJournalofAdvancedManufacturing Technology，2009，43（11）：1106-1123．

［5］吳昭同，楊將新．計算機輔助公差優(yōu)化設(shè)計［M］．杭州：浙江大學出版社，1999：3-9．（Wu Zhao-tong，Yang Jiang-xin．Computer Aided Tolerance Optimization Design［M］．Hangzhou：Zhejiang University Press，1999：3-9．）

［6］劉玉生，楊將新，吳昭同．基于數(shù)學定義的平面尺寸公差數(shù)學模型［J］．機械工程學報，2001，37（9）：12-17．（Liu Yu-sheng，Yang Jiang-xin，Wu Zhao-tong．Mathematical model of size tolerance for plane based on mathematical definition［J］．Chinese Journal of Mechanical Engineering，2001，37（9）：12-17．）

［7］胡潔，吳昭同，楊將新．基于旋量參數(shù)的三維公差累積的運動學模型［J］．中國機械工程，2003（2）：43-46．（Hu Jie，Wu Zhao-tong，Yang Jiang-xin．Kinematics model of 3D tolerance accumulation based on screw parameter［J］．China Mechanical Engineering，2003（2）：43-46．）

［8］Bourdet．A study of optimal-criteria identification based on the small displacement screw model［J］．Annals of the CIRP，1998（37）：503-506．