陳華，唐廣輝，陳志強(qiáng)，李志敏，金隼

（1.上海交通大學(xué)上海市復(fù)雜薄板結(jié)構(gòu)數(shù)字化制造重點(diǎn)實驗室，上海200240；2.上汽通用五菱汽車股份有限公司，廣西柳州545007）

基于雅可比旋量統(tǒng)計法的發(fā)動機(jī)三維公差分析

陳華1，唐廣輝2，陳志強(qiáng)2，李志敏1，金隼1

（1.上海交通大學(xué)上海市復(fù)雜薄板結(jié)構(gòu)數(shù)字化制造重點(diǎn)實驗室，上海200240；2.上汽通用五菱汽車股份有限公司，廣西柳州545007）

為了在發(fā)動機(jī)裝配體的公差分析中考慮傳統(tǒng)尺寸鏈無法處理的幾何公差信息，以發(fā)動機(jī)的做功機(jī)構(gòu)為分析對象，引入雅可比旋量模型，建立在尺寸和形位公差約束下的特征及其配合旋量的變動和約束方程。將其融入到蒙特卡洛算法中，實現(xiàn)發(fā)動機(jī)做功機(jī)構(gòu)的三維統(tǒng)計公差分析。與傳統(tǒng)分析方法相比，基于雅克比旋量模型的分析方法考慮了更多的幾何公差信息，分析結(jié)果也更加精確，該模型可以從尺寸鏈角度建立發(fā)動機(jī)零部件公差與整機(jī)性能的映射關(guān)系。

雅可比旋量；統(tǒng)計法；三維公差分析；做功機(jī)構(gòu)；發(fā)動機(jī)

對于發(fā)動機(jī)而言，燃油消耗率、功率、扭矩等性能指標(biāo)是汽車顧客最為關(guān)心的問題。研究表明，單缸燃燒室壓縮比的波動及其缸間的差異化是影響油耗的一個重要因素，零部件之間的運(yùn)動摩擦功耗對整機(jī)輸出功率有較大影響，這些都與零部件的公差設(shè)計和累積直接相關(guān)。零部件公差與整機(jī)性能的映射關(guān)系是發(fā)動機(jī)領(lǐng)域一個亟待解決的問題。

發(fā)動機(jī)的做功機(jī)構(gòu)，包括曲軸及軸瓦、連桿及連桿瓦、活塞銷和活塞、以及作為支撐件的氣缸體，是發(fā)動機(jī)最為重要的功能組件之一，也是壓縮比的重要構(gòu)件和摩擦功的主要制造者。以某型四缸汽油機(jī)為例，做功機(jī)構(gòu)的累積公差可導(dǎo)致壓縮比波動6%，僅次于缸蓋燃燒室和活塞燃燒室的容積公差導(dǎo)致的波動。同時，該機(jī)構(gòu)的摩擦功占整機(jī)摩擦功的70%以上［1］。

目前，該機(jī)構(gòu)的公差設(shè)計以經(jīng)驗為主，公差分析和分配基于傳統(tǒng)的平面尺寸鏈和極值法，分析對象僅限于零部件的尺寸公差。對于圖紙上大量的幾何公差信息，一直缺乏標(biāo)注依據(jù)和分析手段。

三維公差分析技術(shù)在空間領(lǐng)域描述特征在公差域的微小變動，以及這些變動矢量在裝配鏈中的傳遞關(guān)系，能有效的表達(dá)和傳遞除尺寸公差外的其他幾何公差。目前，基于三維幾何的公差分析模型主要有：以點(diǎn)集形式包含了滿足公差要求的所有公差域邊界及變動要素的T-Map模型［2］；用齊次變換矩陣表示實際特征在公差域中變動的矩陣模型［3］；將包含裝配偏差和幾何偏差的非線性隱式裝配關(guān)系一階線性展開的直接線性法（direct linearization method，DLM）模型［4，5］；視公差域為虛擬關(guān)節(jié)運(yùn)動空間，用雅可比矩陣計算各關(guān)節(jié)與末端之間的公差傳遞的雅可比模型［6］；基于工藝拓?fù)湎嚓P(guān)表面理論，適合公差表示的小位移旋量（small displacement torsor，SDT）模型［7］；以及結(jié)合了雅可比和小位移旋量模型的雅可比旋量模型［8］等。其中，雅可比旋量模型算法簡潔，工程語義豐富，特別適合于以圓柱和平面特征為主的復(fù)雜裝配體三維公差建模，除極值算法外，還能實現(xiàn)統(tǒng)計算法。

將三維公差分析技術(shù)與統(tǒng)計算法結(jié)合起來是當(dāng)前的研究熱點(diǎn)［9］。本文引入雅可比旋量模型，以發(fā)動機(jī)做功機(jī)構(gòu)為分析對象，建立在尺寸和形位公差約束下的特征及其配合旋量的變動和約束方程，將其融入到蒙特卡洛算法中，實現(xiàn)發(fā)動機(jī)做功機(jī)構(gòu)的三維統(tǒng)計公差分析，并詳細(xì)分析計算結(jié)果。

1 雅可比旋量統(tǒng)計算法

1.1 雅可比旋量模型

雅可比旋量模型結(jié)合了2個重要的概念，即適合公差傳遞的雅可比矩陣和適合公差表示的旋量模型。

在公差分析中，可將各公差域所蘊(yùn)含的變動視為虛擬關(guān)節(jié)的運(yùn)動，用雅可比矩陣將其傳遞到裝配尺寸鏈的封閉環(huán)。為方便分析計算，這里將公差分為兩類，即功能要求副（functional requirement，F(xiàn)R）和功能單元副（functional element，F(xiàn)E）。FR為設(shè)計目標(biāo)值，通常為封閉環(huán)；FE根據(jù)是否屬于同一個零件又可以分為內(nèi)部功能單元副（internal FE，IFE）和接觸功能單元副（contact FE，CFE）。雅可比矩陣用于公差分析的傳遞模型可表達(dá)為

式中：δS和δα分別代表FR的3個平動矢量和3個轉(zhuǎn)動矢量；δq[]FEi代表各FE中的小位移矢量；n為FE的數(shù)目；J[]FEi為各關(guān)節(jié)對應(yīng)的雅可比矩陣，其數(shù)學(xué)表達(dá)式為

旋量法是螺旋理論在公差領(lǐng)域的具體應(yīng)用，因其分析對象為小位移轉(zhuǎn)動或平動，故也稱小位移旋量。如圖1所示，公差域可視為實際特征S1對理想特征S0的變動空間，由3個沿坐標(biāo)軸平動的矢量和3個繞坐標(biāo)軸轉(zhuǎn)動的矢量組成。S1相對于S0的小位移旋量可表示為

圖1 面特征的旋量表示模型Fig.1 Screw model of a planar feature

對于S0上的任一點(diǎn)N，S1的變動空間，即公差域的旋量表示為

式中：ZN為N點(diǎn)的法向量，TSU為上公差，TSL為下公差。

小位移旋量模型用運(yùn)動旋量來表示公差變動域，有嚴(yán)格的數(shù)學(xué)定義及豐富的工程語義。但對于公差的傳遞和累積，卻不如雅可比矩陣方便直接［10］。

綜合雅可比矩陣和小位移旋量模型，就是雅可比旋量模型，其表達(dá)式為

從式（6）可以看出，基于雅可比旋量模型的公差分析方法具有簡潔的數(shù)學(xué)表達(dá)式、可量化位置和方向公差、豐富的工程語義等特點(diǎn)。然而，對于雅可比旋量模型，極值法的求解與公差標(biāo)準(zhǔn)相矛盾，因為同一FE的各矢量之間存在約束關(guān)系，而非彼此獨(dú)立。基于蒙特卡洛算法的統(tǒng)計法能在計算過程中考慮這種約束關(guān)系，因而可將統(tǒng)計法引入到雅可比旋量模型中，這就是雅可比旋量統(tǒng)計算法。

1.2 雅可比旋量統(tǒng)計算法

眾所周知，零件的加工質(zhì)量都呈現(xiàn)出一定的統(tǒng)計規(guī)律，其分布特征取決于工序能力和工藝技術(shù)。統(tǒng)計法也意味著一定的廢品率和不完全互換性，但相比于極值法，統(tǒng)計法更符合生產(chǎn)實際，且對組成環(huán)的公差分配較為寬松。

Ghie等［11］提出了雅可比旋量模型的統(tǒng)計算法，但在計算過程中忽略了屬于同一旋量的矢量之間的約束關(guān)系。本文提出了一種新的算法，能將旋量的變動和約束關(guān)系融入到計算流程中，使得計算結(jié)果更加合理。其流程如圖2所示。其中，i為實際仿真次數(shù)，N為設(shè)定仿真次數(shù)。

圖2 雅克比旋量模型統(tǒng)計法計算流程Fig.2 A statistical scheme for the Jacobian-Torsor model

該算法基于蒙特卡洛模擬方法，關(guān)鍵在于建立特征的變動和約束方程，并將其應(yīng)用到蒙特卡洛算法中。對于不同的特征和公差類型，以及特征之間的配合，都對應(yīng)不同的變動和約束方程［12-13］。

發(fā)動機(jī)做功機(jī)構(gòu)的各零部件以圓柱和平面特征為主，裝配約束主要是面匹配和軸孔間隙配合。在各種公差約束下，建立平面和圓柱特征，及其配合的變動和約束方程，然后代入到雅可比旋量統(tǒng)計算法流程中，就可對做功機(jī)構(gòu)進(jìn)行三維統(tǒng)計公差求解。

2 做功機(jī)構(gòu)三維統(tǒng)計公差分析

圖3所示為某型發(fā)動機(jī)第三缸活塞處于上止點(diǎn)時的功能要求示意圖，F(xiàn)R對發(fā)動機(jī)的壓縮比有直接的影響。

圖3 功能要求圖Fig.3 Functional requirement

圖4 為與FR相關(guān)的各零件設(shè)計公差圖。為方便計算，本文將曲軸瓦與曲軸孔視為一體，稱為缸體組件；將連桿大端孔與連桿瓦視為一體，稱為連桿組件。計算過程由基于Matlab編制的程序完成，其步驟如下：

1）建立雅可比旋量模型。做功機(jī)構(gòu)的裝配尺寸鏈傳遞見圖5所示，坐標(biāo)系固結(jié)于各零件公差域的中間位置，坐標(biāo)及方向如圖4所示，各坐標(biāo)方向均相同以避免投影換算。該裝配體包含兩條串聯(lián)環(huán)路：IFE1-IFE9-FR和CFE2-IFE3-CFE4-IFE5-CFE6-FE7-IFE8-FR，故需建立2個雅可比旋量模型。在本案例中忽略了曲軸和連桿的軸向止推約束，以及缸孔與活塞的圓周約束，以避免產(chǎn)生局部并聯(lián)而導(dǎo)致無法建立雅可比傳遞矩陣。

圖4 各零件的設(shè)計要求圖Fig.4 Detailed drawing of each part

圖5 功機(jī)構(gòu)的裝配鏈Fig.5 Connection graph of the working mechanism

由式（2）、（3）可得做功機(jī)構(gòu)的雅可比傳遞矩陣：

式中：dz1＝dz2＝186，dz3＝dz4＝146.5，dz5＝dz6＝dz7＝25.5，dy1＝dy2＝38，其余未說明向量為零。

2）建立旋量的變動范圍和約束。這里給出具有代表性的IFE3和IFE9為例說明。在IFE3中曲軸小頸對大頸有位置度TPO（Φ0.1）和平行度TPAR（Φ0.05））要求，其旋量中各矢量的變動范圍為

約束條件為

IFE9為缸體頂面對曲軸孔的位置TPO（Φ0.1）約束，其變動范圍為

約束條件為

通常，位置公差限制旋量的平動，而方向公差限制其轉(zhuǎn)動。

平面之間的匹配，一般情況下都默認(rèn)為完全貼合，故小位移旋量都為零。而圓柱的情況則不同，因為有間隙配合的情況，例如曲軸與曲軸孔的配合，其旋量的變動和約束方程與位置度相同，這里的“位置度”等于軸孔配合的最大間隙值。

3）計算各旋量中矢量的均值和標(biāo)準(zhǔn)差，按概率分布產(chǎn)生隨機(jī)數(shù)。其均值和標(biāo)準(zhǔn)差可由式（12）、（13）得到，其中K為標(biāo)準(zhǔn)化正態(tài)數(shù)，可依據(jù)預(yù)定廢品率查表得知；μT為變動范圍的中值：

這里假設(shè)各零部件的廢品率都為0.27%，且各矢量都為正態(tài)分布，故K＝3，隨機(jī)數(shù)由Matlab中的randn（）函數(shù)生成。

4）蒙特卡洛仿真計算。計算由程序?qū)崿F(xiàn)，需要注意的是仿真次數(shù)越多，結(jié)果越精確，但也越耗時，這里設(shè)置為10 000次。

5）建立FR的直方圖，并計算其均值和標(biāo)準(zhǔn)差。對于發(fā)動機(jī)做功機(jī)構(gòu)，只關(guān)注FR在活塞行程方向的公差累積，對應(yīng)于FR中的w。在程序中建立w的直方圖，并計算出均值和標(biāo)準(zhǔn)差。計算結(jié)果為：σFR1w＝0.031 1，σFR2w＝0.017 6。二者的均值都非常小，可置為0。

6）結(jié)果計算?；钊麎嚎s高至缸體上表面的公差累積為兩條公差傳遞環(huán)路之和。在統(tǒng)計法中，F(xiàn)R的標(biāo)準(zhǔn)差為FR1和FR2標(biāo)準(zhǔn)差的平方和的根。上一步計算出FR的標(biāo)準(zhǔn)差，從而得到其公差范圍為

即活塞處于上止點(diǎn)時，活塞壓縮高到缸體頂面在z軸方向的累積公差為±0.107 2。

3 結(jié)果分析

對于該做功機(jī)構(gòu)，僅考慮尺寸公差的傳統(tǒng)平面尺寸鏈統(tǒng)計算法結(jié)果為±0.128 9。對比可以看出，本文提出的基于雅可比旋量模型的統(tǒng)計算法在計算過程中考慮了平行度公差的約束，以及旋量中平動和轉(zhuǎn)動矢量間的約束，從而更加符合公差標(biāo)準(zhǔn)且精度更高，同時也為各組成零部件提供了更為寬松的公差設(shè)計空間。

對于做功機(jī)構(gòu)，不僅要知道零部件之間的公差累積，還需要知道各個公差對累積公差的貢獻(xiàn)度，為公差優(yōu)化提供依據(jù)。對于雅可比旋量的統(tǒng)計算法，其公差貢獻(xiàn)度的計算公式為

經(jīng)計算并歸一化處理，可得到各FE的貢獻(xiàn)度排列圖，如圖6所示。

圖6 公差貢獻(xiàn)度排列圖Fig.6 The pareto diagram of tolerance contribution

由圖6可以看出，曲軸孔的位置度對做功機(jī)構(gòu)的公差累積影響最為顯著，達(dá)到21.12%。這也體現(xiàn)了雅可比旋量模型的杠桿效應(yīng)。圖7所示為曲軸孔位置度（Φt）對缸體理論頂面（b）的影響示意圖，曲軸中心線在位置度范圍內(nèi)的轉(zhuǎn)動（a′）可導(dǎo)致頂面（b′）超出平動的變動范圍（±t／2），且距離d越大，杠桿效應(yīng)越明顯。故曲軸孔的位置度控制越嚴(yán)格越好。

圖7 曲軸孔位置度對頂面變動影響示意圖Fig.7 Graph for variational top plane influenced by the positional tolerance of hole of crankshaft

曲軸孔位置度、曲軸與曲軸孔的間隙配合、曲軸小徑對大頸的位置度、以及連桿小端孔對大端孔的位置度對FR的偏差貢獻(xiàn)度之和占總貢獻(xiàn)度的78.1%，而剩余的5項只占22.9%，因此，可嚴(yán)格控制前4項，而適當(dāng)放松后5項，做到質(zhì)量成本的優(yōu)化。

值得注意的是，以上流程都是假設(shè)零件的加工質(zhì)量和累積結(jié)果符合正態(tài)分布，但實際情況中可能還存在其他分布形式，這就需要對測量數(shù)據(jù)采樣分析。對于非正態(tài)分布的參數(shù)，尤其是計算結(jié)果，需分布檢驗程序。對分布的描述，除了均值和標(biāo)準(zhǔn)差，還應(yīng)該有偏度和峰度等參數(shù)［14］。另外，還可繼續(xù)分解式（15）以分別求出其中位置和方向公差的貢獻(xiàn)度，這里不做展開。

4 結(jié)束語

要理清諸如發(fā)動機(jī)等結(jié)構(gòu)復(fù)雜、制造成本相對較高的精密機(jī)構(gòu)零部件公差與性能映射關(guān)系，三維公差分析方法可從尺寸鏈角度提供思路和線索。雅可比旋量統(tǒng)計公差分析方法具有簡潔的數(shù)學(xué)表達(dá)式、可量化幾何公差信息、準(zhǔn)確的計算結(jié)果和豐富的工程語義等特點(diǎn)，尤其適合此類機(jī)構(gòu)的公差分析和優(yōu)化。此外，通過修正雅可比矩陣，實際工況下的熱、載荷等影響因素也可以考慮進(jìn)來［15］，在一定程度上可以減少物理實驗，降低開發(fā)成本。

［1］鐘明明.發(fā)動機(jī)制造精度與性能關(guān)系研究［D］.上海：上海交通大學(xué)，2012：6-26.ZHONG Mingming.Research on correlation between manufacturing precision and performance of engine［D］.Shanghai：Shanghai Jiao Tong University，2012：6-26.

［2］AMETA G，SERGE S，GIORDANO M.Comparison of spatial math models for tolerance analysis：tolerance-maps，deviation domain，and TTRS［J］.ASME Journal of Computing and Information Science in Engineering，2011，11：021004-1-8.

［3］WHITNEY D E，GILBERT O L，JASTRZEBSKI M.Representation of geometric variations using matrix transforms for statistical tolerance analysis in assemblies［J］.Research in Engineering Design，1994，6（4）：191-210.

［4］CHASE K W，MAGLEBY S P，GAO J.Tolerance analysis of 2-D and 3-D mechanical assemblies with small kinematic adjustments［C］／／Advanced Tolerancing Techniques.New York：John Wiley＆Sons，1997：103-137.

［5］CHASE K W，GAO J，MAGLEBY S P，et al.Including geometric feature variations in tolerance analysis of mechanical assembly［J］.IIE Transactions，1996，28（10）：795-807.

［6］LAPERRIèRE L，EIMAREGHY H A.Tolerance analysis and synthesis using Jacobian transforms［J］.CIRP Annals－Manufacturing Technology，2000，49（1）：359-362.

［7］TEISSANDIER D，COUéTARD Y，GéRARD A.A computer aided tolerancing model：proportioned assembly clearance volume［J］.Computer-Aided Design，1999，31：805-817.

［8］DESROCHERS A，GHIE W，LAPERRIèRE L.Application of a unified Jacobian-Torsor model for tolerance analysis［J］.ASME Journal of Computing and Information Science in Engineering，2003，3：1-13.

［9］HONG Y S，CHANG T C.A comprehensive review of tolerancing research［J］.International Journal of Production Research，2002，40（11）：2425-2459.

［10］MARZIALE M，POLINI W.A review of two models for tolerance analysis of an assembly：Jacobian and Torsor［J］.International Journal of Computer Integrated Manufacturing，2011，24（1）：74-86.

［11］GHIE W，LAPERRIèRE L，DESROCHERS A.Statistical tolerance analysis using the unified Jacobian-Torsor model［J］.International Journal of Production Research，2010，48（15）：4609-4630.

［12］ROY U，LI B，Representation and interpretation of geometric tolerances for polyhedral objects.II.Size，orientation and position tolerances［J］.Computer-Aided Design，1999，31（4）：273-285.

［13］蔡敏，楊將新，吳昭同.基于數(shù)學(xué)定義的圓柱要素形狀公差數(shù)學(xué)模型研究［J］.機(jī)械工程學(xué)報，2003，39（12）：86-90.CAI Ming.YANG Jiangxin，WU Zhaotong.Mathematical model of form tolerance for cylindrical feature based on mathematic definition［J］.Chinese Journal of Mechanical Engineering，2003，39（12）：86-90.

［14］NIGAM S D，TURNER J U.Review of statistical approaches to tolerance analysis［J］.Computer-Aided Design，1995，27（1）：5-15.

［15］張為民，陳燦，李鵬忠，等.基于雅可比旋量法的實際工況公差建模［J］.計算機(jī)集成制造系統(tǒng)，2011，17（1）：77-83.ZHANG Weimin，CHEN Can，LI Pengzhong，et al.Tolerance modeling in actual working condition based on Jacoian-Torsor theory［J］.Computer Integrated Manufacturing Systems，2011，17（1）：77-83.

Three-dimensional tolerance analysis of engine based on Jacobian-Torsor statistical model

CHEN Hua1，TANG Guanghui2，CHEN Zhiqiang2，LI Zhimin1，JIN Sun1
（1.Shanghai Key Laboratory of Digital Manufacture for Thin-walled Structures，Shanghai Jiao Tong University，Shanghai 200240，China；2.SAIC GM Wuling Automobile Co.，Ltd.，Liuzhou 545007，China）

The working mechanism of an engine is the target of analysis in order to consider geometric tolerances which cannot be handled by traditional tolerance chains in engine assemblies.The Jacobian-Torsor model is introduced，and variations and constraints of torsors for features and their joints under the restrictions of size and geometric tolerances are built.Based on Monte Carlo simulation，three-dimensional tolerance analysis for working mechanism with a statistical way is achieved.Compared with traditional tolerance analysis methods，the method based on the Jacobian-Torsor statistical model takes more geometric tolerances into consideration and obtains more accurate results.This model can build the relationship between tolerances of parts and engine performance from tolerance chains＇perspectives.

Jacobian-Torsor；statistical method；3D tolerance analysis；working mechanism；engine

10.3969／j.issn.1006-7043.201306058

http：／／www.cnki.net／kcms／doi／10.3969／j.issn.1006-7043.201306058.html

TP391

A

1006-7043（2014）11-1397-06

2013-06-18.網(wǎng)絡(luò)出版時間：2014-09-25.

國家自然科學(xué)基金資助項目（51121063，51175340）；“十二五”國家科技支撐計劃資助項目（2012BAF06B03）.

陳華（1981-），男，博士研究生；金隼（1973-），男，研究員，博士生導(dǎo)師.

金隼，E-mail：jinsun＠sjtu.edu.cn.