麻 凱,高繼東,謝書港

(1.天津大學(xué)機械工程學(xué)院，天津 300072； 2.中國汽車技術(shù)研究中心，天津 300300)

2016057

懸架運動學(xué)特性的區(qū)間可靠性控制*

麻 凱1,2,高繼東2,謝書港2

(1.天津大學(xué)機械工程學(xué)院，天津 300072； 2.中國汽車技術(shù)研究中心，天津 300300)

本文中結(jié)合區(qū)間控制理論、概率論和可靠性理論，提出了計及硬點坐標(biāo)不確定性的懸架運動學(xué)特性區(qū)間可靠性控制方法。該方法將與懸架運動學(xué)特性可靠性相關(guān)的硬點坐標(biāo)參數(shù)區(qū)間作為控制對象，將懸架運動學(xué)特性區(qū)間可靠性作為控制目標(biāo)，進(jìn)行可靠性控制。最后，以前束角曲線和外傾角曲線為區(qū)間可靠性控制目標(biāo)，靈敏度較高的硬點坐標(biāo)區(qū)間作為區(qū)間控制參數(shù)，進(jìn)行仿真。結(jié)果進(jìn)一步驗證了該方法的可行性。

懸架運動學(xué)特性；靈敏度分析；可靠性控制；區(qū)間控制；概率論

前言

近年來，隨著產(chǎn)品質(zhì)量標(biāo)準(zhǔn)的不斷提高，對結(jié)構(gòu)可靠性的控制成為研究的熱點。工程中所有的機械結(jié)構(gòu)都受到各種因素帶來的不確定性影響，而這些不確定性可能來自于邊界條件、結(jié)構(gòu)尺寸和材料屬性等因素，這些不確定因素會對結(jié)構(gòu)的應(yīng)力、應(yīng)變、疲勞壽命和運動特性產(chǎn)生影響。這些影響對可靠性評估十分重要，不可忽略。在可靠性的控制中，結(jié)構(gòu)的可靠性按照傳統(tǒng)的定義方法規(guī)定為，在給定的可預(yù)期的使用環(huán)境下，結(jié)構(gòu)完成預(yù)定工作的能力[1-2]。隨著研究的深入，可靠性又被細(xì)化定義為概率可靠性和非概率可靠性，它們分別基于概率模型和非概率模型對結(jié)構(gòu)可靠性進(jìn)行研究。

基于概率論的可靠性研究已比較成熟。概率論即隨機理論[3-5]，指在規(guī)定的使用環(huán)境下，機構(gòu)完成指定工作的概率，在工程領(lǐng)域內(nèi)它已經(jīng)被廣泛使用。而隨著人們對不確定性認(rèn)識的深入，隨機方法被發(fā)現(xiàn)已經(jīng)不能準(zhǔn)確地描述所有的不確定性問題。因此,僅用概率方法研究可靠性問題是不夠的。于是在該模型基礎(chǔ)上，發(fā)展了不同的研究可靠性的方法，例如：灰色理論[6-7]、模糊理論[8-10]、凸模型理論[11-15]和區(qū)間理論[14-15]。并出現(xiàn)了使用以上兩種及兩種以上方法組合描述不確定性的可靠性分析方法[16-19]。這些方法彌補了用傳統(tǒng)隨機方法分析可靠性的不足。在各種描述不確定性的方法中，由于區(qū)間方法對數(shù)據(jù)樣本信息量要求較低，用區(qū)間方法評價結(jié)構(gòu)的可靠性得到了廣泛應(yīng)用。

文獻(xiàn)[20]中提出通過控制懸架結(jié)構(gòu)參數(shù)容差保證其運動學(xué)特性曲線的一致性的方法,提出利用區(qū)間方法建立懸架運動學(xué)特性曲線一致性的不確定性結(jié)構(gòu)參數(shù)模型。它把懸架結(jié)構(gòu)中參數(shù)容差看作不確定參數(shù),即區(qū)間參數(shù),通過區(qū)間控制方法控制其外特性的一致性。

本文中將區(qū)間控制方法與概率方法相結(jié)合，設(shè)計了一種新的可靠性控制方法，該方法可以控制目標(biāo)的分布區(qū)間和分布規(guī)律。

本文中將所設(shè)計的方法應(yīng)用在懸架運動學(xué)特性曲線可靠性控制中，提出了懸架運動學(xué)特性區(qū)間可靠性控制方法。使用該方法，設(shè)定合理的結(jié)構(gòu)硬點坐標(biāo)參數(shù)的中值和概率分布規(guī)律，達(dá)到控制懸架運動學(xué)特性可靠性的目的。

最后，通過某款車的雙橫臂懸架運動學(xué)特性曲線可靠性的控制，進(jìn)一步證明了這種方法的可行性。

1 實數(shù)區(qū)間及其運算

(1)

(2)

(3)

為區(qū)間XI的中值和半徑，即不確定量。

[XC,XC]+[-ΔX,ΔX]=XC+ΔXI=

XC+ΔX[-1,1]=XC+ΔXeΔ

(4)

式中：ΔXI=[-ΔX,ΔX]；eΔ=[-1,1]。這種表示方法稱為區(qū)間的中心區(qū)間表示法。

2 區(qū)間控制方法

考慮3個函數(shù)f1(X),f2(X)和f3(X)，都含有4個自變量X=(x1,x2,x3,x4)，略去高階項，1階泰勒展開式為

(5)

可知,df1/dx1≈a1,df2/dx1≈a2,df3/dx1≈a3,當(dāng)x1增量為Δx1=(df1/dx1)-1=1/a1,式(5)中的3個函數(shù)的增量Δf1,Δf2和Δf3可分別表示為

(6)

(7)

(8)

i=1,2,3,4；k=1,2,3

(9)

如果給定Δf1(X),Δf2(X),Δf3(X)，求滿足條件的ΔX=(Δx1,Δx2,Δx3,Δx4)，則可利用式(10)先求出λ1,λ2,λ3。

(10)

再將求出的λ1,λ2,λ3代入式(11)，求出Δx1,Δx2,Δx3,Δx4。

(11)

通過上面提出的區(qū)間半徑控制方法來實現(xiàn)區(qū)間控制，具體步驟如下。

步驟1：設(shè)一個極小值Δx，即初始的區(qū)間半徑，差分計算函數(shù)靈敏度。

(12)

步驟2：計算

(13)

步驟3：用如下方程組求λ1,λ2,…,λm。

(14)

步驟4：計算

(15)

步驟5：計算

Δfk(X)=fk(X+ΔX)-fk(X)

(16)

如果‖Δfk(X)-Δfkmax(X)‖f≤ε，停止計算，(Δfkmax(X)為控制目標(biāo)半徑)，這時式(15)得到的Δx就是滿足Δfkmax(X)要求的參數(shù)區(qū)間半徑。否則，令xi=xi+Δxi，計算

(17)

并回到第2步繼續(xù)迭代計算。這里，ε是一個給定的閾值。

3 正態(tài)分布

若隨機變量X服從一個位置參數(shù)為μ、尺度參數(shù)為σ的概率分布，且其概率密度函數(shù)為

(18)

則這個隨機變量就稱為正態(tài)隨機變量，正態(tài)隨機變量服從的分布就稱為正態(tài)分布，記作X～N(μ,σ2)。

正態(tài)分布是有兩個參數(shù)μ和σ2的連續(xù)型隨機變量的分布，參數(shù)μ是服從正態(tài)分布的隨機變量的均值，參數(shù)σ2是此隨機變量的方差，所以正態(tài)分布記作N(μ,σ2)。服從正態(tài)分布的隨機變量取離μ越近的值的概率越大，而取離μ越遠(yuǎn)的值的概率越?。沪以叫?，分布越集中在μ附近，σ越大，分布越分散。

σ描述正態(tài)分布數(shù)據(jù)分布的離散程度，σ越大，數(shù)據(jù)分布越分散；σ越小，數(shù)據(jù)分布越集中。也稱為是正態(tài)分布的形狀參數(shù)，σ越大，曲線越扁平；反之，σ越小，曲線越瘦高。

4 懸架運動學(xué)特性曲線區(qū)間可靠性控 制方法

在生產(chǎn)過程中，加工誤差是不可避免的。要保證懸架運動學(xué)特性曲線的可靠性，必須合理控制加工精度，即保證一定比例的懸架運動學(xué)特性分布在給定的懸架運動學(xué)特性的不確定區(qū)間中。因為工程中硬點坐標(biāo)和懸架運動學(xué)特性都滿足正態(tài)分布X～N(μ,σ2)，所以σ可看作σ分布在[-σ,σ]區(qū)間，即σ也可視為區(qū)間參數(shù)。

因此，可靠性控制主要控制懸架運動學(xué)特性曲線分布范圍fI=[-f,f]和分布率σI=[-σ,σ]，合并向量FI=[fI，σI]。懸架運動學(xué)特性曲線一般是一條參數(shù)隨輪跳變化的曲線?？刂七@些曲線較理想的方法是控制曲線的兩端。所以，本文中將對懸架運動學(xué)特性曲線的兩個端點的分布區(qū)間和分布率控制來實現(xiàn)對懸架運動學(xué)特性可靠性的控制。主要步驟如下。

5 算例

以雙橫臂獨立懸架模型為例，懸架模型和硬點名稱見圖1?？刂魄?，前束角隨輪跳而變化的曲線如圖2所示，中間曲線為確定性參數(shù)下的懸架模型的前束角隨輪跳變化曲線，另兩條線為由結(jié)構(gòu)參數(shù)容差引起的前束角隨輪跳變化曲線的上下邊界線，即曲線區(qū)間上界和下界。由于結(jié)構(gòu)參數(shù)與前束角變化曲線具有單調(diào)性，且輪跳為零時，所有前束角的值相同，所以前束角邊界曲線如果在輪跳正向為上邊界時，在輪跳負(fù)向為下邊界，另一條邊界曲線相反。同理，外傾角隨輪跳而變化的曲線如圖3所示。

圖1 雙橫臂獨立懸架模型及控制點

圖2 控制前雙橫臂獨立懸架車輪前束角曲線區(qū)間

圖3 控制前雙橫臂獨立懸架車輪外傾角曲線區(qū)間

由結(jié)構(gòu)的對稱性，選取Z方向7對靈敏度最大的硬點坐標(biāo)參數(shù)容差為要控制的區(qū)間參數(shù)的半徑，見表1。初始參數(shù)下的前束曲線區(qū)間和尺度參數(shù)見表2，初始參數(shù)下的外傾曲線區(qū)間和尺度參數(shù)見表3。

表1 控制前初始硬點坐標(biāo)區(qū)間半徑和尺度參數(shù)

表2 控制前的前束角區(qū)間半徑和尺度參數(shù)

表3 控制前的外傾角區(qū)間半徑

設(shè)計制造結(jié)構(gòu)參數(shù)區(qū)間和各自尺度參數(shù)，控制懸架外特性一致性要求如下：

(1) 車輪前束角θtoe在輪跳-50mm～50mm變化區(qū)間為±0.05°，尺度參數(shù)為0.1；

(2) 車輪外傾角θcamber在輪跳-50mm～50mm變化區(qū)間為±0.05°，尺度參數(shù)為0.1。

通過對前束角和外傾角的靈敏度分析可知，同一區(qū)間參數(shù)下前束角和外傾角的區(qū)間半徑都隨輪跳增加而增大，所以控制前束角和外傾角曲線一致性時，可以通過控制正、負(fù)向最大輪跳點的前束角和外傾角來實現(xiàn)。

區(qū)間可靠性控制后結(jié)構(gòu)參數(shù)區(qū)間半徑和尺度參數(shù)見表4，前束角曲線區(qū)間見圖4，外傾角曲線區(qū)間見圖5，輪跳極值點前束角區(qū)間半徑和尺度參數(shù)見表5，輪跳極值點外傾角區(qū)間半徑和尺度參數(shù)見表6。

表4 控制后硬點坐標(biāo)區(qū)間和尺度參數(shù)

圖4 控制后雙橫臂獨立懸架車輪前束角曲線區(qū)間

圖5 控制后雙橫臂獨立懸架車輪外傾角曲線區(qū)間

輪跳/mm前束角中值/(°)前束角區(qū)間半徑/(°)尺度參數(shù)σ-502.27440.0550.10450-1.41210.0470.0107

表6 控制后的外傾角區(qū)間半徑和尺度參數(shù)

從結(jié)果可以看出，控制后前束曲線和外傾曲線都在目標(biāo)區(qū)間范圍內(nèi)，其分布規(guī)律也滿足控制要求。因此，所有控制目標(biāo)均滿足區(qū)間可靠性控制要求。

6 結(jié)論

本文中將區(qū)間控制理論與概率論相結(jié)合，建立應(yīng)用于工程的區(qū)間可靠性方法，并應(yīng)用該方法合理設(shè)計汽車懸架硬點的坐標(biāo)區(qū)間和均方差，從而保證了懸架運動學(xué)特性的區(qū)間和分布規(guī)律，實現(xiàn)了可靠性控制。算例中通過對雙橫臂獨立前懸架前束角和外傾角的區(qū)間可靠性控制，進(jìn)一步說明了這種區(qū)間可靠性控制方法是可行的。

[1] 董聰.現(xiàn)代結(jié)構(gòu)系統(tǒng)可靠性理論及其應(yīng)用[M].北京:科學(xué)出版社,2001.

[2] ELISHAKOFF I. Discussion on: A Non-probabilistic Concept of Reliability[J]. Structural Safety,1995,17(3):195-199.

[3] 呂震宙,馮蘊雯.結(jié)構(gòu)可靠性問題研究的若干進(jìn)展[J].力學(xué)進(jìn)展,2000,30(1):21-28.

[4] ZHAO Y G, LU Z J. Fourth-moment Standardization for Structural Reliabilityassessment[J]. J Struct Eng,2007,133(7):916-924.

[5] ZHAO Y G, ONO T. A General Procedure for First/second-order Reliabilitymethod[J]. Structural Safety,1999,21(2):95-112.

[6] 肖芳淳.復(fù)合材料結(jié)構(gòu)灰色可靠性分析與設(shè)計[J].強度與環(huán)境,2000(1):60-65.

[7] 段齊駿,張福祥.復(fù)雜系統(tǒng)可靠性研究中的灰色問題[J].系統(tǒng)工程理論與實踐,1999,119(12):96-98.

[8] HELTON J C, OBERKAMPF W L. Alternative Representations of Epistemicuncertainty[J]. Reliability Engineering & System Safety,2004,85(1-3):1-10.

[9] KLIR G J, SMITH R M. On Measuring Uncertainty and Uncertainty-basedinformation: Recent Developments[J]. Annal Math Art Int,2001,32(1-4):5-33.

[10] BERND M. Discussion on: Structural Reliability Analysis Through Fuzzynumber Approach, with Application to Stability[J]. Computers & Structures,2004,82(2-3):325-327.

[11] BEN-HAIM Y. A Non-probabilistic Concept of Reliability[J]. Structural Safety,1994,14(4):227-245.

[12] BEN-HAIM Y. A Non-probabilistic Measure of Reliability of Linear Systemsbased on Expansion of Convex Models[J]. Structural Safety,1995,17(2):91-109.

[13] AU F, CHENG Y S, LG THAM E A. Robust Design of Structures Usingconvex Models[J]. Computers & Structures,2003,81(28-29):2611-2619.

[14] 郭書祥,呂震宙,馮元生.基于區(qū)間分析的結(jié)構(gòu)非概率可靠性模型[J].計算力學(xué)學(xué)報,2001,18(1):56-60.

[15] 王登剛,李杰.結(jié)構(gòu)可靠性分析的區(qū)間方法[C].第12屆全國反應(yīng)堆結(jié)構(gòu)力學(xué)會議論文專輯,2002:244-247.

[16] 董玉革.隨機應(yīng)力模糊強度時機械零件模糊可靠性設(shè)計方法的研究[J].機械設(shè)計,1999,16(3):11-13.

[17] 黃洪鐘,孫占全,郭東明,等.隨機應(yīng)力模糊強度時模糊可靠性的計算理論[J].機械強度,2001,23(3):305-307.

[18] 麻凱,陳塑寰.結(jié)構(gòu)區(qū)間參數(shù)的二階近似動態(tài)優(yōu)化方法[J].東北大學(xué)學(xué)報(科學(xué)版),2008(S2):400-405.

[19] JIANG Chao, LIU Lixin, TANG Yifei. A Structural Reliability Analysis Method for Hybrid Model with Interval Variables[J]. Automotive Engineering,2012(34):727-732.

[20] 麻凱,管欣,逄淑一,等.懸架運動學(xué)特性一致性區(qū)間控制方法[J].吉林大學(xué)學(xué)報(工學(xué)版),2011,4(41):910-914.

Interval Reliability Control for Suspension Kinematic Characteristics

Ma Kai1,2, Gao Jidong2& Xie Shugang2

1.SchoolofMechanicalEngineering,TianjinUniversity,Tianjin300072; 2.ChinaAutomotiveTechnology&ResearchCenter,Tianjin300300

By combining the theories of interval control, probability and reliability, an interval reliability control method for suspension kinematic characteristics is proposed with consideration of the uncertainty of hard point coordinates in this paper. In the method, reliability control is performed with the interval of hard point coordinate parameters related to the reliability of suspension kinematic characteristics as control object, and the interval reliability of suspension kinematic characteristics as control objective. Finally a specific simulation is conducted with toe-in angle and camber angle curves as the objectives of interval reliability control and the interval of hard point coordinates with higher sensibility as interval control parameters. The results further verify the feasibility of the method proposed.

suspension kinematic characteristics; sensitivity analysis; reliability control; interval control; probability theory

*吉林省科技廳青年科研基金(20130522180JH)資助。

原稿收到日期為2015年7月15日。