馮志杰，肖慧婷，楊永鋒

1.航空工業(yè)航宇救生裝備有限公司航空防護(hù)救生技術(shù)航空科技重點(diǎn)實驗室，湖北 襄陽 441000

2.西北工業(yè)大學(xué)，陜西 西安 710072

機(jī)構(gòu)在使用過程中會受到各種因素的影響，其中振動環(huán)境會影響機(jī)構(gòu)關(guān)鍵部件的使用壽命，甚至影響整個機(jī)構(gòu)的安全性及可靠性。而載荷譜作為機(jī)構(gòu)關(guān)鍵部件疲勞壽命及環(huán)境可靠性研究的重要組成部分，能否對其進(jìn)行正確的數(shù)據(jù)處理與編制，對于保障整個機(jī)構(gòu)的安全性及可靠性具有十分重要的意義。

李雷[1]等建立了某型彈射座椅的有限元計算模型，提出了飛機(jī)彈射座椅振動響應(yīng)的有限元計算方法?；谠摲椒ǎT志杰[2]等對某型飛機(jī)彈射座椅椅載設(shè)備進(jìn)行了隨機(jī)振動仿真研究，獲取了椅載設(shè)備兩個不同測點(diǎn)處的振動載荷譜，并對其進(jìn)行數(shù)據(jù)歸納。當(dāng)選取少量測點(diǎn)對振動數(shù)據(jù)進(jìn)行歸納時，出現(xiàn)偶然誤差的概率較大，可能導(dǎo)致得到的載荷譜過考核或欠考核，因此需對部件多個測點(diǎn)響應(yīng)進(jìn)行歸納分析。豐志強(qiáng)[3]等針對飛機(jī)機(jī)載設(shè)備在振動環(huán)境可靠性試驗剖面的編制需要，依據(jù)標(biāo)準(zhǔn)對機(jī)載設(shè)備振動試驗數(shù)據(jù)進(jìn)行了歸納。田永衛(wèi)[4]等對某型飛機(jī)多個艙位和主要結(jié)構(gòu)部位進(jìn)行了振動量值的飛行實測，并對多個不同測點(diǎn)的試驗數(shù)據(jù)采用參數(shù)歸納法進(jìn)行了數(shù)據(jù)處理。依據(jù)GJB 150A的準(zhǔn)則，參數(shù)歸納法僅適用于總體近似服從正態(tài)分布的數(shù)據(jù)集[5]，但由于測點(diǎn)位置的選取可能會影響響應(yīng)數(shù)據(jù)集概率分布情況，因此當(dāng)測得的數(shù)據(jù)集總體不完全服從正態(tài)分布時，參數(shù)歸納法便不再適用，應(yīng)采用非參數(shù)歸納法處理不同測點(diǎn)數(shù)據(jù)。

本文對于機(jī)構(gòu)部件載荷譜獲取方法進(jìn)行了研究，首先對機(jī)構(gòu)進(jìn)行隨機(jī)振動分析，然后提取機(jī)構(gòu)關(guān)鍵部件上不同測點(diǎn)的響應(yīng)數(shù)據(jù)，再通過非參數(shù)上限統(tǒng)計歸納方法處理響應(yīng)數(shù)據(jù)?；谟行院蜔o偏性分析，評價了不同歸納方法的載荷譜歸納結(jié)果，最終得出較優(yōu)的機(jī)構(gòu)關(guān)鍵部件載荷譜歸納結(jié)果，為后續(xù)關(guān)于部件的可靠性及壽命的研究提供幫助。

1 機(jī)構(gòu)部件動力學(xué)響應(yīng)獲取方法研究

若直接將機(jī)構(gòu)安裝位置載荷譜作為關(guān)鍵部件安裝位置處的載荷譜對部件進(jìn)行動力學(xué)研究[6-7]，則不能考慮從激勵輸入點(diǎn)到關(guān)鍵部件之間的傳力路徑，因此提出更符合實際的機(jī)構(gòu)關(guān)鍵部件載荷譜獲取方法[8-12]。

將機(jī)構(gòu)安裝位置處的載荷譜作為機(jī)構(gòu)載荷譜，對機(jī)構(gòu)進(jìn)行隨機(jī)振動仿真分析，提取關(guān)鍵部件安裝位置附近不同測點(diǎn)處的響應(yīng)，對不同測點(diǎn)處的響應(yīng)進(jìn)行歸納包絡(luò)，將歸納后的曲線作為關(guān)鍵部件的載荷譜，以探究關(guān)鍵部件在振動環(huán)境下的動力學(xué)響應(yīng)及壽命等。機(jī)構(gòu)關(guān)鍵部件載荷譜獲取過程如圖1所示。

圖1 機(jī)構(gòu)關(guān)鍵部件載荷譜獲取流程Fig.1 Load spectrum acquisition process of key components of mechanism

隨機(jī)振動的計算過程包括兩個階段[13]：（1）計算系統(tǒng)響應(yīng)特性，在時域內(nèi)用脈沖響應(yīng)函數(shù)h(t)表示，在頻域內(nèi)用復(fù)頻響函數(shù)H(ω)描述，其物理意義為系統(tǒng)響應(yīng)與激勵之比；（2）通過系統(tǒng)響應(yīng)特性和隨機(jī)激勵計算功率譜密度矩陣Sxx(ω)。

通過隨機(jī)振動分析得到關(guān)鍵部件不同測點(diǎn)的響應(yīng)后，對結(jié)果進(jìn)行歸納。以往在對機(jī)構(gòu)關(guān)鍵部件響應(yīng)數(shù)據(jù)進(jìn)行處理時，直接用仿真或?qū)崪y的某一個點(diǎn)或某兩個點(diǎn)的包絡(luò)作為關(guān)鍵部件的載荷譜[1-2]。取點(diǎn)數(shù)目較少時，得到的部件載荷譜不夠準(zhǔn)確，因此需對關(guān)鍵部件多個測點(diǎn)響應(yīng)進(jìn)行歸納研究。

2 振動數(shù)據(jù)歸納方法研究

GJB 150A 的18A[5]中規(guī)定了5 種數(shù)據(jù)集上限的統(tǒng)計分析方法，包含兩種參數(shù)上限統(tǒng)計估計方法和三種非參數(shù)上限統(tǒng)計估計方法。兩種參數(shù)上限統(tǒng)計估計方法分別是正態(tài)單邊容差上限法（NTL)和正態(tài)預(yù)測上限法（NPL)；三種非參數(shù)上限統(tǒng)計估計方法分別是包絡(luò)上限法（ENV)、無驗前分布容差上限法（DFL)和經(jīng)驗容差上限法（ETL)。參數(shù)上限統(tǒng)計法適用于滿足正態(tài)分布或者經(jīng)轉(zhuǎn)化后滿足正態(tài)分布的數(shù)據(jù)集，而非參數(shù)上限統(tǒng)計法適用于不滿足正態(tài)分布或者轉(zhuǎn)化后不滿足正態(tài)分布的數(shù)據(jù)集。在對部件不同測點(diǎn)響應(yīng)進(jìn)行數(shù)據(jù)歸納時，測點(diǎn)位置的選取會影響數(shù)據(jù)集分布，因此采用非參數(shù)上限統(tǒng)計方法歸納部件不同測點(diǎn)數(shù)據(jù)集。

2.1 ENV法

ENV 法的計算方法為選取數(shù)據(jù)集中最大估計值作為最大上限，以譜曲線一個頻段上的譜線作為研究對象，所有樣本在該段譜線上的容差上限為

式中：N表示數(shù)據(jù)集的樣本容量，xi表示第i個樣本的譜線值。

該方法計算簡單，忽略了樣本的分布特性，無法給出超過該值的概率。當(dāng)估計集有異常的情況下，包絡(luò)上限法可能過于保守，對于譜線帶寬也很敏感。

2.2 DFL法

DFL法計算方法同ENV法一致，取數(shù)據(jù)集的最大值作為上限估計值，但不同于ENV 法的是，DFL 法引入了置信度和覆蓋率。

式中：β表示覆蓋率，γ表示置信度。用xβ代表數(shù)據(jù)集樣本的真實上限，則數(shù)據(jù)集樣本小于xβ的概率為β，即

若要使樣本中的最大值xmax小于xβ，則數(shù)據(jù)集中每一個樣本都應(yīng)小于xβ。每個樣本之間相互獨(dú)立，故有

xmax≥xβ的概率為

工程中通常使用的置信度至少為50%，覆蓋率為95%[11]。若要使估計上限置信度γ≥50%且β= 95%，由式（5)計算得N≥14；若要保持置信度大小不變，取覆蓋率較大時(β＞0.95)，計算所需的樣本數(shù)量較大。

2.3 ETL法

經(jīng)驗容差上限一般表示為ETL(β)。假設(shè)共有N個測量點(diǎn)，每個測量點(diǎn)的輸出響應(yīng)譜曲線有M段譜線，則總數(shù)據(jù)集中共包含NM個樣本

式中：j= 1,2,…,M，i= 1,2,…,N，mj表示第j段譜線的樣本均值。在M段譜線上構(gòu)造的正則化估計集為

將正則化估計集{u}中所有數(shù)據(jù)按從升序排列，數(shù)據(jù)集中第k個元素以u(k)來表示。構(gòu)造{u}的經(jīng)驗分布函數(shù)FN

取覆蓋率為β時，有

采用經(jīng)驗分布函數(shù)FN代替實際分布函數(shù)F，取覆蓋率為β，由式（9)計算得k值，用u(k)作為容差上限uβ的估計值。每段譜線上的經(jīng)驗容差上限為

取置信度為50%，此時xβj表示在第j段譜線上，EBL(β)以50%的置信度超過所有值的100β%。若選擇大于xβj的值，置信度會增加。當(dāng)樣本容量N達(dá)到一定要求時，不同頻率分辨率帶寬上的正則化估計集應(yīng)服從同一分布。

ETL法是根據(jù)整個頻帶分布計算各頻率段譜線上的容差上限，因此可以消除一些偶然因素帶來的誤差，可以得到相對穩(wěn)定的結(jié)果，但其樣本量通常應(yīng)大于10[10]。

3 上限估計方法優(yōu)劣性評價

不同方法歸納得到的載荷譜必然存在優(yōu)良之分，故需要評價哪個估計方法相對較好。評價標(biāo)準(zhǔn)主要有無偏性評價和有效性評價。（1）無偏性評價：對于真實容差上限θ，其估計量?滿足E=θ，則稱估計量?為真實容差上限θ的無偏估計量。（2）有效性評價：設(shè)總體X～F(x;θ)，對于真實容差上限θ的兩個無偏估計量?，，如果，則比有效。

本文通過對不同歸納方法得到的載荷譜結(jié)果的有效性和無偏性評價，對比分析歸納結(jié)果的優(yōu)良性。取總樣本容量為N，將每段譜線樣本按樣本觀測值升序排列，取第95%N個樣本作為總體覆蓋率為95%的理論上限?；跓o偏性和有效性，評價各方法得到的容差上限估計量的優(yōu)劣性。當(dāng)各容差上限估計量均值越接近理論上限，方差越小，則說明該方法估計得到的載荷譜越接近真實工況。

4 實例分析

以某型彈射座椅為例，通過有限元分析獲得仿真結(jié)果，將椅載關(guān)鍵設(shè)備安裝位置不同測點(diǎn)處的響應(yīng)結(jié)果進(jìn)行載荷譜歸納。能準(zhǔn)確獲得人椅系統(tǒng)彈射座艙的各種參數(shù)是提高救生性能的關(guān)鍵[12]，但在緊急狀況下彈射座椅幫助飛行員順利出艙的首要前提是保證彈射座椅關(guān)鍵設(shè)備在飛機(jī)正常飛行過程中的安全性和可靠性。采用程序控制器作為彈射座椅關(guān)鍵部件，討論不同數(shù)據(jù)歸納方法得到的載荷譜歸納結(jié)果的優(yōu)良性。

4.1 模型建立及響應(yīng)計算

以保證主結(jié)構(gòu)完整性、傳力路徑不變?yōu)樵瓌t，對某型飛機(jī)彈射座椅模型零部件進(jìn)行簡化。該型飛機(jī)彈射座椅主要材料為鋁合金，其各項材料參數(shù)見表1。

表1 材料參數(shù)Table 1 Material parameters

座椅激勵譜是通過加速度傳感器實測得到，再根據(jù)GJB 150[5]中相關(guān)規(guī)定對試驗數(shù)據(jù)歸納處理得到。以某型彈射座椅振動耐久垂向載荷譜為例進(jìn)行計算分析，載荷曲線如圖2所示。

圖2 飛機(jī)座艙振動耐久垂向加速度載荷譜Fig.2 Aircraft cockpit vibration endurance vertical acceleration load spectrum

完成有限元模型處理后，先對模型進(jìn)行頻率響應(yīng)分析獲取頻響函數(shù)，再對其進(jìn)行隨機(jī)響應(yīng)分析，獲取測點(diǎn)處振動響應(yīng)。在程序控制器安裝位置附近各取900個節(jié)點(diǎn)作為響應(yīng)輸出測點(diǎn)，座椅簡化模型及程序控制器安裝位置如圖3所示。

圖3 彈射座椅簡化模型及程控器安裝位置Fig.3 Ejection seat simplified model and installation position of program controller

4.2 正態(tài)性檢驗

當(dāng)譜密度曲線樣本的每個頻率分辨率帶寬數(shù)據(jù)集服從正態(tài)分布時，加速度均方根值（RMS）也近似服從正態(tài)分布[10]。各測試點(diǎn)的加速度RMS值可由式（12）計算

式中：Gj和Δfj分別為第j個譜線上的譜密度和頻率帶寬。

加速度RMS值的分布特性可以反映信號的采樣點(diǎn)數(shù)、平均值和方差，因此可以通過對不同測點(diǎn)RMS值的分布研究譜密度曲線樣本是否服從正態(tài)分布。分別繪制程控器安裝位置不同測點(diǎn)處振動響應(yīng)的加速度均方根值Q-Q 圖和直方圖，如圖4和圖5所示。

圖4 程控器各測加速度點(diǎn)RMS值Q-Q圖Fig.4 Quantile-Quantile diagram of RMS value of each accel eration point measured by the program controller

圖5 程控器各測點(diǎn)加速度RMS值直方圖Fig.5 Histogram of acceleration RMS value of each measuring point of the program controlle

通過程控器安裝位置附近不同測點(diǎn)處的加速度RMS值Q-Q圖和直方圖，可判斷該譜密度曲線集不完全滿足正態(tài)分布，因此需對其進(jìn)行非參數(shù)上限統(tǒng)計歸納。

4.3 歸納結(jié)果

對于ENV、DFL 法，計算方法較為簡單，按照第2 節(jié)所述連接各段譜線上最大值，最終得到載荷歸納譜；ETL方法需將各個測試點(diǎn)的所有頻段上的譜線數(shù)據(jù)提取出來作為估計集，依據(jù)式（6）將所有樣本值進(jìn)行正則化，得到正則化估計集，對該正則化估計集取覆蓋率為β= 95%，利用式（9）～式（11）計算得到每段譜線寬上的經(jīng)驗容差上限值，終得到95%覆蓋率下的載荷歸納譜。采用ENV、DFL和ETL方法繪制得到的歸納曲線如圖6～圖8所示。

圖6 ENV法歸納譜曲線Fig.6 The spectrum curve induced by ENV method

圖7 DFL法歸納譜曲線Fig.7 The spectrum curve induced by DFL method

圖8 ETL法歸納譜曲線Fig.8 The spectrum curve induced by ETL method

圖6和圖7中的載荷歸納譜一樣，原因是ENV法和DFL法計算容差上限估計值的方法是一樣的。兩種方法的區(qū)別在于DFL法引入了參數(shù)置信度和覆蓋率，可以通過樣本容量和其中一個參數(shù)計算另一個參數(shù)。在本文后面討論中將ENV 法和DFL 法歸并討論。對比ENV/DFL 法和ETL 法的歸納結(jié)果（見圖9）可看出，由于ETL法考慮了整個頻率范圍內(nèi)的分布，消除了一部分由取點(diǎn)位置隨機(jī)性等因素產(chǎn)生的譜線值波動，導(dǎo)致其在峰值點(diǎn)的上限估計值均小于ENV/DFL法峰值點(diǎn)的上限估計值。為了驗證計算的準(zhǔn)確性，對比參考文獻(xiàn)[2]中的試驗結(jié)果與ETL 法歸納結(jié)果（見圖10）可看出，試驗結(jié)果與歸納結(jié)果幅值及趨勢基本一致。計算得歸納譜加速度RMS=4.8844g，試驗測得加速度RMS=5.0893g，計算得誤差為4.02%。

圖9 ENV/DFL法和ETL法歸納結(jié)果對比Fig.9 Comparison of inductive results between ENV/DFL and ETL

圖10 ETL法歸納結(jié)果與試驗結(jié)果對比Fig.10 Comparison between ETL inductive results and experimental results

4.4 歸納結(jié)果優(yōu)劣性評價

對三種方法的歸納結(jié)果進(jìn)行無偏性和有效性評價。取總體樣本容量N為900，按升序排列各譜線值，取第855 個數(shù)據(jù)（95%覆蓋率）作為理論上限值。樣本容量分別設(shè)定為5、13、30、45、60、75 和90，每個容量下重復(fù)抽樣10 次，計算10次抽樣的載荷譜歸納結(jié)果的均值和方差。載荷譜線的加速度RMS值能反映信號的平均能量，峰值點(diǎn)的波動大小能反映載荷對機(jī)構(gòu)部件的破壞性強(qiáng)弱，因此取載荷譜線的加速度RMS 值和峰值f= 62.78Hz 處的歸納結(jié)果作為討論對象進(jìn)行優(yōu)良性評價。首先分別計算各樣本容量10 次抽樣歸納結(jié)果的均值和方差，然后將計算所得均值除以理論上限進(jìn)行歸一化處理。計算得到的結(jié)果見表2～表5，對比不同歸納方法得到的峰值點(diǎn)和加速度RMS 值的歸一化均值及方差可以發(fā)現(xiàn)：（1）無論是對于峰值點(diǎn)還是加速度RMS值，ENV/DFL 法得到的歸納結(jié)果隨著樣本容量的增加，歸一化平均值不斷增加，且方差波動幅度較大，無法趨于穩(wěn)定；（2）ETL 法由于考慮了各測點(diǎn)每段譜線整體分布，其歸納結(jié)果消除了一部分波動。所以隨樣本量增加，其歸納結(jié)果的歸一化均值變化幅度較小，方差逐漸趨于穩(wěn)定。

表2 不同抽樣次數(shù)下各歸納方法在f=62.78Hz 處的歸一化上限平均值Table 2 The normalized upper bound mean of each inductive method at f=62.78Hz under different sampling times

表3 不同抽樣次數(shù)下各歸納方法在f=62.78Hz處的上限值方差（(g2/Hz)2?10?6）Table 3 The upper bound variance of each inductive method at f=62.78Hz under different sampling times

表4 不同抽樣次數(shù)下各歸納方法加速度RMS值的歸一化上限平均值Table 4 The normalized upper bound mean of each inductive method of acceleration RMS under different sampling times

表5 不同抽樣次數(shù)下各歸納方法加速度RMS值的上限值方差（(g2/Hz)2?10?4）Table 5 The upper bound variance of each inductive method of acceleration RMS under different sampling times

綜上所述，ETL 方法得到的載荷譜歸納結(jié)果的無偏性和有效性均優(yōu)于ENV/DFL方法得到的歸納結(jié)果。

5 結(jié)束語

本文建立了一種機(jī)構(gòu)關(guān)鍵部件載荷譜獲取方法。首先對機(jī)構(gòu)進(jìn)行隨機(jī)振動分析，獲取關(guān)鍵部件測點(diǎn)響應(yīng)，然后對響應(yīng)數(shù)據(jù)集進(jìn)行歸納處理，最終得出部件的載荷譜。通過對比不同歸納方法的理論分析以及實例結(jié)果，得出以下結(jié)論：當(dāng)樣本容量N≤13時采用ENV/DFL方法得到的估計結(jié)果波動較大，故采用ENV/DFL 法時樣本量不宜過小；ETL方法相對ENV/DFL方法得到的估計結(jié)果更穩(wěn)定，但當(dāng)樣本量較大時ETL法對于載荷譜波動較小的頻段上限估計值較大，而波動較大的峰值附近其上限估計值較小，因此采用ETL法時樣本量不宜過大。