陳一凡 張根保,3 冉 琰 李宇龍 庾 輝
1.重慶大學(xué)機(jī)械工程學(xué)院,重慶,400044 2.重慶大學(xué)機(jī)械傳動(dòng)國(guó)家重點(diǎn)實(shí)驗(yàn)室,重慶,400044 3.重慶文理學(xué)院機(jī)械工程學(xué)院,重慶,402160
可靠性分配是產(chǎn)品早期設(shè)計(jì)和開(kāi)發(fā)階段的關(guān)鍵步驟之一,其目的是將產(chǎn)品整機(jī)的可靠性指標(biāo)分配給每個(gè)基本單元,以確保組成整機(jī)的各基本單元的可靠性滿足可靠性設(shè)計(jì)要求[1-3]。目前,傳統(tǒng)的可靠性分配方法主要側(cè)重于兩個(gè)方面:一是基于模糊理論的可靠性分配[4-5];二是基于優(yōu)化算法的可靠性優(yōu)化設(shè)計(jì)[6-8]。雖然這些可靠性分配技術(shù)得到了廣泛的應(yīng)用,但是,前者方法主要依賴專家評(píng)分,主觀性較強(qiáng),而后者方法忽略了機(jī)械系統(tǒng)的結(jié)構(gòu)特性,此外,它們有一個(gè)共性問(wèn)題,即它們都沿襲了電子產(chǎn)品的可靠性分配模式,因此,為了提高機(jī)械系統(tǒng)的可靠性,研究一種適用于機(jī)械系統(tǒng)的可靠性分配方法迫在眉睫。
近幾年來(lái),國(guó)內(nèi)外學(xué)者對(duì)機(jī)械系統(tǒng)的可靠性分配工作進(jìn)行了相關(guān)研究。張強(qiáng)等[9]考慮機(jī)械系統(tǒng)分配中的多層影響因素,建立了兩層因素的可靠性分配模型。王昊等[10]提出一種考慮子系統(tǒng)故障相關(guān)性的多因素可靠性分配方法。張義民等[11]通過(guò)引入模糊數(shù)學(xué)理論,提出了一種基于用戶反饋的數(shù)控機(jī)床故障信息的可靠性模糊分配方法。張根保等[12]引入“任務(wù)”概念,提出了一種基于“任務(wù)”的數(shù)控機(jī)床模糊可靠性分配方法。YU等[13]提出了一種考慮失效影響和可靠性成本的綜合模糊分配方法,隨后,他們?cè)谀繕?biāo)可行性(feasibility-of-objectives,F(xiàn)OO)方法的基礎(chǔ)上,提出了一種考慮依賴因子的分配方法[14]。CHENG等[15]將模糊分配與最大熵有序加權(quán)平均(ME-OWA)方法相結(jié)合并將其應(yīng)用于數(shù)控機(jī)床的子系統(tǒng),提高了分配的靈活性。ZHANG等[16]提出了一種綜合考慮運(yùn)動(dòng)機(jī)構(gòu)的組成和三類失效模式的可靠性分配方法,并成功將其應(yīng)用于飛機(jī)齒輪門(mén)鎖機(jī)構(gòu)。LIU等[6]為了找到全局最優(yōu)的可靠性指標(biāo),提出了一種基于多學(xué)科設(shè)計(jì)優(yōu)化的可靠性分配建模求解方法。
然而,現(xiàn)階段機(jī)械系統(tǒng)的可靠性分配研究還存在以下關(guān)鍵問(wèn)題需要進(jìn)一步解決:①現(xiàn)有的可靠性分配仍然按照“系統(tǒng)-部件-零件”的路徑分解,而實(shí)際的機(jī)械系統(tǒng)為由運(yùn)動(dòng)單元組成的傳動(dòng)鏈,傳統(tǒng)的自頂向下的分解會(huì)導(dǎo)致分配任務(wù)越來(lái)越繁瑣,分配效率較低[14,17];②運(yùn)動(dòng)單元的可靠性決定了整機(jī)系統(tǒng)的可靠性,在不確定性條件下,各運(yùn)動(dòng)單元對(duì)系統(tǒng)可靠性的影響存在一定的差異,這一點(diǎn)在傳統(tǒng)的可靠性分配中并未體現(xiàn)。
針對(duì)以上問(wèn)題,本文提出了一種基于元?jiǎng)幼骷捌潇`敏度的機(jī)械傳動(dòng)系統(tǒng)可靠性分配方法。從運(yùn)動(dòng)和傳動(dòng)的角度,將機(jī)械系統(tǒng)分解為最基本的運(yùn)動(dòng)單元(即元?jiǎng)幼?,并將其作為可靠性分配的最小粒度,然后,考慮到處于傳動(dòng)系統(tǒng)結(jié)構(gòu)不同位置的元?jiǎng)幼鲗?duì)整個(gè)系統(tǒng)可靠性的影響存在很大的差異,本文引入靈敏度指標(biāo)對(duì)其進(jìn)行不確定性量化。該方法能夠?qū)㈦y以定量評(píng)估的問(wèn)題進(jìn)行量化,以靈敏度指標(biāo)作為可靠性分配的有力依據(jù),從而為可靠性分配工作提供新的思路。
針對(duì)傳統(tǒng)的機(jī)械傳動(dòng)系統(tǒng)可靠性分配的不足,本文提出一種基于元?jiǎng)幼骷捌潇`敏度的機(jī)械傳動(dòng)系統(tǒng)可靠性分配方法,其分配框架如圖1所示。首先,采用“功能-運(yùn)動(dòng)-動(dòng)作(function-motion-action,F(xiàn)MA)”的分解方法[17-19]對(duì)機(jī)械系統(tǒng)進(jìn)行分解得到元?jiǎng)幼骷?,并將其作為分配的最小粒?然后,利用拉丁超立方體抽樣(Latin hypercube sampling,LHS)方法計(jì)算低階逼近(low-rank approximations,LRA)元模型的系數(shù),并通過(guò)后處理LRA的系數(shù)求解各元?jiǎng)幼髟趥鲃?dòng)鏈中的靈敏度;最后,基于獲得的元?jiǎng)幼黛`敏度計(jì)算結(jié)果識(shí)別出系統(tǒng)的關(guān)鍵元?jiǎng)幼?,并將傳?dòng)系統(tǒng)的整機(jī)可靠性指標(biāo)分配到各元?jiǎng)幼?。雖然基于LRA的Sobol’方法在其他領(lǐng)域已有相關(guān)研究[20],但是將其應(yīng)用于可靠性分配領(lǐng)域還未見(jiàn)相關(guān)報(bào)道。
圖1 可靠性分配框架Fig.1 The reliability allocation framework
注:不帶箭頭的實(shí)線連接表示組成關(guān)系;帶有箭頭的實(shí)線連接表示映射關(guān)系;帶有箭頭的虛線連接表示傳動(dòng)關(guān)系圖2 FMA分解示意圖Fig.2 Diagram of FMA decomposition
功能分解是將產(chǎn)品功能分解為相應(yīng)的子系統(tǒng)或部件的運(yùn)動(dòng)。傳統(tǒng)的功能分解方法[21-22]的最后節(jié)點(diǎn)處是針對(duì)零件而言的,如果分解對(duì)象為復(fù)雜的大型機(jī)電產(chǎn)品,則不利于產(chǎn)品的性能分析。FMA方法是從運(yùn)動(dòng)的角度將產(chǎn)品功能分解到最基本的運(yùn)動(dòng)單元,即元?jiǎng)幼?FMA分解過(guò)程如圖2所示。從圖2中可以看出,元?jiǎng)幼鲗訉?shí)際上是由一條傳動(dòng)鏈組成的,只要確保各元?jiǎng)幼髂軌蛘_\(yùn)行,那么就能實(shí)現(xiàn)產(chǎn)品所需的功能,這與可靠性分配的目的是一致的。因此,功能分解的目的是明確各元?jiǎng)幼鞯脑敿?xì)可靠性輸入值,以進(jìn)一步合理設(shè)計(jì)各元?jiǎng)幼鳌?/p>
由于數(shù)控機(jī)床是一個(gè)串聯(lián)結(jié)構(gòu),各功能的正常輸出取決于所有元?jiǎng)幼鞯恼9ぷ?,因此,可以將它們視為一個(gè)串聯(lián)系統(tǒng)。根據(jù)可靠性理論,元?jiǎng)幼鞯目煽啃耘c系統(tǒng)的可靠性關(guān)系表示為
(1)
式中,Ri為元?jiǎng)幼鰽i的可靠度;RS為傳動(dòng)系統(tǒng)的可靠度;n為元?jiǎng)幼鲾?shù)目。
隨著時(shí)間的推移,系統(tǒng)可靠性呈單調(diào)遞減趨勢(shì),這可以理解為系統(tǒng)中某個(gè)或多個(gè)元?jiǎng)幼鞒霈F(xiàn)退化,而在這個(gè)退化過(guò)程中存在不確定性因素,如應(yīng)力、磨損和員工的操作方式等,這些不確定性因素造成元?jiǎng)幼鞯男阅芴幱谝粋€(gè)區(qū)間范圍,因此,每個(gè)元?jiǎng)幼骺梢钥闯墒窍到y(tǒng)的一個(gè)隨機(jī)變量,則式(1)可以等效為
Y=f(X1,X2,…,Xn)
(2)
將式(1)轉(zhuǎn)換成式(2)的目的是為了后續(xù)分析每個(gè)隨機(jī)變量Xi(i=1,2,…,n)對(duì)函數(shù)Y的靈敏度,即元?jiǎng)幼鲗?duì)系統(tǒng)可靠性的影響程度。
1.2.1全局靈敏度分析
全局靈敏度分析(global sensitivity analysis,GSA)定義為響應(yīng)模型中輸入變量的不確定性對(duì)模型輸出響應(yīng)的不確定性的貢獻(xiàn)程度。Sobol’方法[23]是工程實(shí)際應(yīng)用中最廣泛的GSA方法之一,它可以分析單一輸入變量作用或多個(gè)輸入變量之間相互作用對(duì)系統(tǒng)輸出的影響。
為不失一般性,設(shè)函數(shù)f(X)在單位超立方體Ω(n)中,且輸入為n維,則Ω(n)可以表示為
Ω(n)={X|0≤Xi≤1;i=1,2,…,n}
(3)
Sobol’法的關(guān)鍵步驟是將函數(shù)f(X)進(jìn)行分解,即
(4)
其中,f0為常數(shù),fi是Xi的函數(shù),fij是Xi和Xj的函數(shù),其余各項(xiàng)對(duì)其所包含的每一個(gè)因素的積分為0,即
(5)
1≤k≤s
由于函數(shù)分解中的所有項(xiàng)之間都為正交,且式(4)的分解存在唯一性,那么所有展開(kāi)式的和都可以通過(guò)積分以遞歸的方式計(jì)算得到,即
(6)
式中,fi表示Xi單獨(dú)變化對(duì)Y的影響(稱為Xi的主效應(yīng));fij表示除了i、j各自變化的影響之外,Xi和Xj同時(shí)變化對(duì)Y的影響;E(·)為期望算子。
對(duì)于f2(X)可積函數(shù),可以對(duì)函數(shù)分解進(jìn)行平方和積分,進(jìn)而得到f(X)的總方差,即
(7)
式(7)等式左邊為Y的方差,右邊為方差項(xiàng)。根據(jù)Xi的集合分解式,方差表達(dá)式的分解表示為
(8)
(9)
式中,X~i為去除Xi之外的所有變量的集合,即X~i={X1,…,Xi-1,Xi+1,…,Xn}。
因此,一階Sobol’指標(biāo)(first-order Sobol’indices,F(xiàn)SI)和總Sobol’指標(biāo)(total Sobol’indices,TSI)分別表示為
(10)
(11)
FSI反映了單一輸入變量Xi的不確定性對(duì)輸出變量的影響,Si越大,說(shuō)明該輸入變量對(duì)系統(tǒng)輸出的影響越顯著。
TSI用于評(píng)估單一輸入變量及其與其他輸入變量的相互作用對(duì)系統(tǒng)輸出的共同影響。如果輸入變量間相互獨(dú)立,則FSI值總是小于TSI值;若兩者相差較大,則說(shuō)明輸入變量間的相互作用對(duì)系統(tǒng)輸出的影響更顯著。
通常,計(jì)算Sobol’法指標(biāo)的常用方法有兩種,即基于蒙特卡羅(Monte Carlo,MC)方法和基于元模型方法。相對(duì)來(lái)說(shuō),元模型(即替代模型)比MC方法計(jì)算效率更高,因?yàn)樗怯媒y(tǒng)計(jì)上的等價(jià)模型來(lái)代替復(fù)雜模型的。常用的元模型有多項(xiàng)式混沌展開(kāi)(polynomial chaos expansion,PCE)和LRA,鑒于LRA比PCE的計(jì)算結(jié)果誤差更小[24],本文提出一種基于LRA的靈敏度后處理方法,該方法具有計(jì)算效率高、耗時(shí)短等優(yōu)勢(shì)。
1.2.2基于LRA的全局靈敏度指標(biāo)計(jì)算方法
LRA是元模型中處理高維模型的一種有效方法,其目的是用少量的秩一張量的和來(lái)逼近模型的輸出響應(yīng)值,并利用多元基的張量積形式減少未知系數(shù)的數(shù)目。因此,Sobol’靈敏度指標(biāo)可以通過(guò)后處理LRA的系數(shù)來(lái)計(jì)算得到[20],多項(xiàng)式基的LRA元模型可以表示為
(12)
因此,基于LRA的Sobol’方法的計(jì)算分別取決于以下FSI和TSI的表達(dá)式:
(13)
(14)
由于式(13)、式(14)都涉及均值和方差,因此,將模型Y代入LRA元模型中,然后基于LRA系數(shù)計(jì)算出這些值,則模型輸出的均值和方差分別近似為
(15)
(16)
Si表達(dá)式中的E[E(Y|Xi)2]可以近似為
(17)
STi表達(dá)式中的E[E(Y|X~i)2]可以近似表示為
(18)
對(duì)于更高階指標(biāo)的表達(dá)式也可以根據(jù)LRA系數(shù)獲得。通過(guò)利用多項(xiàng)式基的張量積形式可以將未知系數(shù)的數(shù)量減少一個(gè)數(shù)量級(jí),進(jìn)而使用相對(duì)較少的秩數(shù)以足夠的精度逼近響應(yīng)值。因此,基于LRA的Sobol’方法能夠快速計(jì)算出靈敏度指標(biāo)。
為了使整個(gè)傳動(dòng)系統(tǒng)可靠性保持相對(duì)穩(wěn)定,應(yīng)該降低元?jiǎng)幼鲗?duì)系統(tǒng)可靠性的靈敏度。以靈敏度指標(biāo)作為分配權(quán)重,將具有較高靈敏度的元?jiǎng)幼鞣峙湟愿叩目煽啃?。因此,元?jiǎng)幼髋c傳動(dòng)系統(tǒng)可靠性的映射關(guān)系表示為
(19)
式中,S′i為元?jiǎng)幼鞯南鄬?duì)權(quán)重,記為ωi。
然后,將ωi取反(即1-ωi),且對(duì)其進(jìn)行歸一化處理,從而得到元?jiǎng)幼鞯淖罱K分配權(quán)重ω′為
(20)
最后,各元?jiǎng)幼鞣峙涞玫降目煽慷葹?/p>
(21)
數(shù)控機(jī)床作為復(fù)雜機(jī)電產(chǎn)品的經(jīng)典類型之一,在制造業(yè)加工領(lǐng)域扮演著重要的角色,例如汽車、航空工業(yè)和船舶工業(yè)等。特別是早期的產(chǎn)品設(shè)計(jì),將整機(jī)或系統(tǒng)的可靠性目標(biāo)值合理地分配給各運(yùn)動(dòng)單元是非常關(guān)鍵的,合理的可靠性分配能提高產(chǎn)品的質(zhì)量,并進(jìn)一步提高產(chǎn)品的生產(chǎn)效率。鑒于數(shù)控轉(zhuǎn)臺(tái)是加工中心的一個(gè)關(guān)鍵功能部件[17,25],本文將其作為實(shí)際應(yīng)用案例來(lái)驗(yàn)證提出方法的可行性和有效性,數(shù)控轉(zhuǎn)臺(tái)的結(jié)構(gòu)示意圖請(qǐng)掃描本文首頁(yè)二維碼查看。根據(jù)產(chǎn)品可靠性要求,數(shù)控轉(zhuǎn)臺(tái)的分度回轉(zhuǎn)運(yùn)動(dòng)可靠度RS=0.89。
根據(jù)數(shù)控轉(zhuǎn)臺(tái)運(yùn)動(dòng)特性,利用FMA分解方法將數(shù)控轉(zhuǎn)臺(tái)的分度運(yùn)動(dòng)分解為6個(gè)元?jiǎng)幼鱗26],共組成4條傳動(dòng)鏈,如圖3所示。
圖3 數(shù)控轉(zhuǎn)臺(tái)FMA分解Fig.3 FMA decomposition of CNC turntable
運(yùn)用Sobol’法計(jì)算全局靈敏度時(shí),需要先給定參數(shù)的變化區(qū)間及概率分布函數(shù)。數(shù)控機(jī)床常見(jiàn)的可靠性指標(biāo)主要有可靠度、失效率和平均故障間隔時(shí)間(mean time between failure,MTBF),它們之間往往是同一概念的不同表達(dá)形式。為了便于分析,這里使用MTBF作為統(tǒng)計(jì)的指標(biāo)。分析現(xiàn)有少量的數(shù)控轉(zhuǎn)臺(tái)元?jiǎng)幼鞯墓收蠑?shù)據(jù)發(fā)現(xiàn),大多數(shù)元?jiǎng)幼鞯姆植碱愋头木鶆蚍植?。由文獻(xiàn)[27]可知,在不確定性因素條件下,如果隨機(jī)變量的分布不是特別明確,使用均勻分布計(jì)算得到的靈敏度結(jié)果不會(huì)發(fā)生顯著性影響。此外,均勻分布能夠避免工程實(shí)際工程中數(shù)據(jù)量少、無(wú)法確定分布的問(wèn)題,只需要確定各元?jiǎng)幼鱉TBF的最小值和最大值便可進(jìn)行靈敏度分析。因此,本文選擇均勻分布作為元?jiǎng)幼鞯募僭O(shè)分布類型。具體取值情況如表1所示。
表1 數(shù)控轉(zhuǎn)臺(tái)的元?jiǎng)幼鞣植己腿≈登闆rTab.1 Distribution and value of the meta-actions of the CNC turntable
本文利用UQLab軟件[28]建立了基于LRA的Sobol’法的元?jiǎng)幼黛`敏度模型,其概率輸入模型由6個(gè)獨(dú)立的隨機(jī)變量組成。利用LHS方法對(duì)元模型進(jìn)行隨機(jī)抽樣,并基于200個(gè)樣本的實(shí)驗(yàn)設(shè)計(jì)創(chuàng)建數(shù)控轉(zhuǎn)臺(tái)可靠性模型的LRA元模型。通過(guò)LRA系數(shù)快速求解元?jiǎng)幼鞯腇SI值和TSI值,并計(jì)算出總靈敏度貢獻(xiàn)率(TSR)。各元?jiǎng)幼鲗?duì)系統(tǒng)可靠性的靈敏度結(jié)果如表2所示。
表2 元?jiǎng)幼鲗?duì)數(shù)控轉(zhuǎn)臺(tái)可靠性靈敏度結(jié)果Tab.2 Reliability sensitivity results of meta-actions to the CNC turntable
由表2可知,靈敏度大小排序依次為SA4、SA1、SA5、SA2、SA3、SA6,其中A4對(duì)數(shù)控轉(zhuǎn)臺(tái)的可靠性影響起決定性作用,即是數(shù)控轉(zhuǎn)臺(tái)可靠性的薄弱環(huán)節(jié)。其原因可解釋為:由圖3可知,數(shù)控轉(zhuǎn)臺(tái)由4條傳動(dòng)鏈組成,其中M342是最長(zhǎng)的傳動(dòng)鏈,即由三級(jí)元?jiǎng)幼鹘M成,從而增大了不確定性因素的影響,這正是在復(fù)雜機(jī)電產(chǎn)品設(shè)計(jì)中縮短傳動(dòng)鏈的主要原因之一。另外,A4位于傳動(dòng)鏈M342的末端,則上游元?jiǎng)幼鞯男阅芊€(wěn)定性輸出會(huì)對(duì)下游造成某種程度的影響[18],而A4的性能穩(wěn)定性又直接決定M342是否能夠完成指定的運(yùn)動(dòng)。因此,為了提高數(shù)控轉(zhuǎn)臺(tái)的可靠性,將A4分配更高的可靠性是有必要的。A1、A4、A5的靈敏度指標(biāo)為正,且TSR明顯比較大,這表明它們的影響很大程度上會(huì)降低數(shù)控轉(zhuǎn)臺(tái)的可靠性,在設(shè)計(jì)環(huán)節(jié)必須采取適當(dāng)措施提高其可靠性。另外,可以觀察到各個(gè)元?jiǎng)幼鞯腡SI值與FSI值的差值小于0.1,這說(shuō)明元?jiǎng)幼髦g存在較弱的相互作用。為了簡(jiǎn)化計(jì)算過(guò)程,這里不考慮相互作用。
根據(jù)式(19)~式(21)計(jì)算出各元?jiǎng)幼鞯目煽慷?,如?所示。
表3 元?jiǎng)幼骺煽啃杂成浣Y(jié)果Tab.3 Result of reliability mapping of meta-actions
此外,以MC方法模擬計(jì)算結(jié)果作為基準(zhǔn)檢驗(yàn)LRA元模型計(jì)算結(jié)果的準(zhǔn)確性和收斂性。圖4給出了3個(gè)元?jiǎng)幼鞯腇SI值和TSI值隨MC方法采樣次數(shù)增加的收斂過(guò)程,可以發(fā)現(xiàn)采樣規(guī)模N=100 000時(shí)MC方法的計(jì)算結(jié)果具有較高的準(zhǔn)確性。
(a)FSI收斂過(guò)程
在上述實(shí)例基礎(chǔ)上,本文將提出的方法與MC方法和PCE方法作對(duì)比分析,仿真結(jié)果如圖5所示。由仿真結(jié)果可知,基于LRA的Sobol’法的靈敏度結(jié)果與其他兩種方法的結(jié)果基本是一致的,這說(shuō)明抽樣規(guī)模取200時(shí)滿足收斂條件。由圖5可知,基于MC的Sobol’方法的樣本規(guī)模為100 000,而本文提出的方法在較少的抽樣規(guī)模(N=200)條件下保證了計(jì)算結(jié)果的精度,且計(jì)算效率更高。
(a)FSI計(jì)算結(jié)果
由上述分析結(jié)果可知,A4、A1和A5是數(shù)控轉(zhuǎn)臺(tái)的關(guān)鍵元?jiǎng)幼?,根?jù)“木桶效應(yīng)”,為了提高數(shù)控轉(zhuǎn)臺(tái)的可靠性,應(yīng)盡可能使靈敏度較大的元?jiǎng)幼鞯目煽啃圆话l(fā)生顯著變化。因此,提高它們的MTBF,即A4的平均故障間隔時(shí)間tMTBFA4=[658,1363]、tMTBFA1=[837,1821]和tMTBFA5=[889,1682],然后再次對(duì)它們的可靠性靈敏度進(jìn)行分析,其結(jié)果如圖6所示。
圖6 優(yōu)化前后元?jiǎng)幼鞯撵`敏度結(jié)果Fig.6 Sensitivity results of meta-actions before and after optimization
由圖6可以看出,優(yōu)化后A4的靈敏度明顯比優(yōu)化前下降了,A1和A5優(yōu)化前后的靈敏度差異不大。此外,為了檢驗(yàn)優(yōu)化效果,以靈敏度標(biāo)準(zhǔn)差作為波動(dòng)程度的評(píng)價(jià)標(biāo)準(zhǔn),通過(guò)計(jì)算得到優(yōu)化前后的標(biāo)準(zhǔn)差分別為0.160和0.096,這就說(shuō)明優(yōu)化后的整體靈敏度對(duì)系統(tǒng)的影響要小于優(yōu)化前整體靈敏度對(duì)系統(tǒng)的影響。
本文利用Sobol’方法分析了元?jiǎng)幼鲗?duì)整個(gè)傳動(dòng)系統(tǒng)可靠性的影響程度,利用“功能-運(yùn)動(dòng)-動(dòng)作”(FMA)分解得到的元?jiǎng)幼髯鳛闄C(jī)械傳動(dòng)系統(tǒng)可靠性分配的最小粒度,這有利于減小可靠性分配的工作量,提高分配效率。同時(shí),以元?jiǎng)幼鲗?duì)系統(tǒng)的輸出影響作為分配的依據(jù),避免了影響較小的單元分配較高的可靠性設(shè)計(jì)值,從而使得系統(tǒng)可靠性達(dá)到最優(yōu)?;诘碗A逼近(LRA)的Sobol’方法在較少樣本條件下能夠準(zhǔn)確得到元?jiǎng)幼鲗?duì)系統(tǒng)可靠性的不確定性輸出響應(yīng)。在此基礎(chǔ)上,為提高數(shù)控轉(zhuǎn)臺(tái)的運(yùn)行穩(wěn)定性,對(duì)關(guān)鍵元?jiǎng)幼鞯目煽啃赃M(jìn)行了優(yōu)化,從而為可靠性設(shè)計(jì)工作提供了指導(dǎo)。
本文從元?jiǎng)幼鞯慕嵌冉⒘丝煽啃苑峙淠P?,在今后的研究中,為了獲取組成機(jī)械系統(tǒng)所有零件的可靠性設(shè)計(jì)值,組成元?jiǎng)幼鞯钠渌慵?如支撐件、緊固件和傳動(dòng)件等)的可靠性設(shè)計(jì)值也有待進(jìn)一步研究。