楊繼坤，徐廷學(xué)，趙 亮，陳 卓

(1.海軍航空工程學(xué)院 研究生管理大隊，山東 煙臺264001;2.海軍航空工程學(xué)院 兵器科學(xué)與技術(shù)系，山東 煙臺264001;3.中國人民解放軍92076 部隊，北京102202)

0 引 言

對導(dǎo)彈武器系統(tǒng)可用度的評估，通常是根據(jù)裝備某個時期內(nèi)統(tǒng)計的各項時間因素值來進行計算，由于統(tǒng)計時間周期長，影響因素眾多，這種方法不能很好地反映裝備真實可用度水平［1］。目前，研究導(dǎo)彈武器系統(tǒng)等可修系統(tǒng)的主要數(shù)學(xué)工具是隨機過程理論。文獻［2］評定了出現(xiàn)保障延誤時間的可修系統(tǒng)穩(wěn)態(tài)可用度。文獻［3-4］利用馬爾可夫更新過程分析了偵察系統(tǒng)的可靠性和裝甲裝備的使用可用度。文獻［5］基于更新過程進行了導(dǎo)彈武器系統(tǒng)的可用度分析。

在對武器系統(tǒng)進行可用度分析時，由于經(jīng)濟上或者試驗環(huán)境難以模擬等原因，并不能確切知道影響可用度的時間分布和參數(shù)值，而只能根據(jù)已有的實驗和專家的經(jīng)驗確定這些分布類型和參數(shù)區(qū)間，將其描述為隨機模糊變量是更客觀、更合理的。目前，隨機模糊理論在可靠性領(lǐng)域有一定的應(yīng)用。文獻［6-7］基于模糊隨機理論建立了結(jié)構(gòu)可靠性分析模型。文獻［8］基于模糊理論進行了導(dǎo)彈系統(tǒng)貯存可靠性仿真方法研究。

本文首先基于馬爾可夫更新過程和拉普拉斯變換工具，建立考慮技術(shù)陣地設(shè)備部分戰(zhàn)損情況下和具有保障延誤的串聯(lián)可修導(dǎo)彈武器系統(tǒng)可用度分析模型。針對可用度分析中存在模糊因素的特點，將導(dǎo)彈視作模糊系統(tǒng)，引入隨機模糊的概念，采用三角模糊數(shù)量化模糊信息，估算可用度水平。

1 馬爾可夫更新過程

設(shè)隨機變量Zn取值在E={0，1，…，K}中，隨機變量Tn取值在［0，∞)中，且有0=T0≤T1≤T2≤…≤Tn，n=0，1，…，如果對所有j ∈E，t ≥0，有P{Zn+1=j，Tn+1-Tn≤t| Z0，Z1，…，Zn，T0，T1，…，Tn}=P{Zn+1=j，Tn+1-Tn≤t| Zn}，則隨機過程(Z，T)={Zn，Tn，n=0，1，…}稱為狀態(tài)空間E 上的馬爾可夫更新過程。經(jīng)簡化求得馬爾可夫更新方程為［3-4］

即:

其中，Ai(t)為系統(tǒng)可用度;gi(t)為定義在［0，∞)上的非負、在任何有限區(qū)間上的有界函數(shù)，在計算可用度時，通常這個函數(shù)是不同裝備服從任意分布的壽命、維修、保障延誤的時間。

在本文中，G(t)的Laplace 變換為

Laplace-Stieltjes 變換為

R(s)表示復(fù)變數(shù)s 的實部，Gk(t)表示G(t)的k 重卷積［9］，即，且G(0)(t)=1，t ≥0。

2 隨機模糊理論

設(shè)(Θ，P(Θ)，Cr)為可信性空間，其中Θ 為非空集合，P(Θ)為Θ 的冪集，Cr 為定義在P(Θ)上的可信性測度［10-12］。

定義1:假設(shè)ξ 為一從可信性空間(Θ，P(Θ)，Cr)到實數(shù)集R 上的函數(shù)，則ξ 為模糊變量。

定義2:令ξ 為可信性空間(Θ，P(Θ)，Cr)上的模糊變量，α ∈(0，1］，稱和分別為ξ 的α-悲觀值和α-樂觀值。

定義3:設(shè)ξ 為模糊變量，它的期望值E［ξ］是有限的，則。

定義4:若a=(al，am，an)，0 ＜al≤am≤al，稱a 為一個三角模糊數(shù)，如果其隸屬度函數(shù)可以表示為

定義5:設(shè)a 為三角模糊數(shù)，對?λ ∈(0，1)，a的截集aλ為有界閉區(qū)間［(am-al)λ+al，(am-an)λ+an］。

3 導(dǎo)彈武器系統(tǒng)可用度分析模型

在導(dǎo)彈武器系統(tǒng)的可用度研究中，大多數(shù)文獻考慮的系統(tǒng)都是假定部件故障后立即得到修理。但實際情況并非如此，而是常常有一段隨機的等待時間(稱這段時間為保障延誤時間)，例如申請等待備件等。而且在現(xiàn)實戰(zhàn)場環(huán)境下，敵方主要打擊的目標(biāo)是導(dǎo)彈技術(shù)測試陣地，對技術(shù)陣地的打擊會損毀技術(shù)陣地的測試和維修設(shè)備(包括測試中的導(dǎo)彈)，使其部分戰(zhàn)損，需要對其進行更換(或修理)后再工作，這都將影響到導(dǎo)彈技術(shù)測試、對接任務(wù)的進行，使技術(shù)陣地不能及時提供裝配好的導(dǎo)彈供作戰(zhàn)使用。因此，在系統(tǒng)可用度分析中，考慮技術(shù)陣地測試設(shè)備部分戰(zhàn)損和保障延誤的可修系統(tǒng)十分必要，而且更有重要的理論意義和實用價值［13-14］。

3.1 模型構(gòu)建步驟

1)明確系統(tǒng)組成。設(shè)系統(tǒng)為由n 個分系統(tǒng)串聯(lián)而成的可修系統(tǒng);

2)確定隨機模糊變量。時間分布未知的分系統(tǒng)和修理設(shè)備，要通過專家的知識與經(jīng)驗來描述其分布類型和參數(shù)情況，這些量即為模糊變量;

3)建立基于馬爾可夫更新過程的可用度模型，得出可用度的解析模型(本章下面兩小節(jié)主要對這方面進行研究)。

4)運用三角模糊數(shù)估算可用度水平。將隨機模糊變量用三角模糊數(shù)來表示，根據(jù)隨機模糊相關(guān)定理，估算其α-悲觀值和α-樂觀值，再根據(jù)定義3 求解其穩(wěn)態(tài)可用度。

3.2 基本條件和假設(shè)

1)系統(tǒng)的第i 個部件壽命Xi服從指數(shù)分布Fi(t)=1-e-λit，t ≥0，λi＞0;

2)系統(tǒng)發(fā)生故障，如果不能及時得到維修，視為不可用狀態(tài)，產(chǎn)生的保障延誤時間Wi服從一般分布Hi(t)，其平均保障延誤時間為

3)進行故障維修時，修理時間Yi服從一般的分布Gi(t)，其平均修理時間為

修復(fù)后系統(tǒng)立即轉(zhuǎn)為工作狀態(tài);

4)當(dāng)1 個部件故障等待修理或者正在修理時，其他部件停止工作，不再發(fā)生故障，此時系統(tǒng)處于故障狀態(tài)。部件修理完恢復(fù)如新，由指數(shù)分布的性質(zhì)，知道在同一時刻不會有多個部件故障;

5)技術(shù)陣地修理設(shè)備在執(zhí)行任務(wù)過程中，會因部分戰(zhàn)損而出現(xiàn)失效，其服務(wù)期U 服從指數(shù)分布U(t)=1-e-αt，t ≥0，α ＞0，修理設(shè)備失效后立即進行更換，其更換時間V 統(tǒng)一服從一般分布V(t)，且平均更換時間為，修理設(shè)備更換如新;

6)Xi，Yi，Wi，U，V 相互獨立。

3.3 構(gòu)建系統(tǒng)可用度模型

定義系統(tǒng)的廣義修理時間γ 為:在導(dǎo)彈技術(shù)陣地，測試和修理設(shè)備從開始測試和修理故障系統(tǒng)的時刻起，直到該部件技術(shù)準(zhǔn)備完成的時間長度，其中包括測試與修理設(shè)備可能的戰(zhàn)損而對測試和修理設(shè)備進行更換的時間。令，則

其LS 變換和平均廣義修理時間為

A(t)=P{X(t)=1 | X(0)=1}。

若令τk，δk，γk分別表示系統(tǒng)的第k 個工作壽命、第k 個保障延誤時間和第k 個廣義修理時間，不難證明，{τk+δk+γk，k=1，2，…}是個更新過程。其分布函數(shù)為

對式(6)2 端作LS(Laplace-Stieltjes)變換得

式(8)右端第1 項中，當(dāng)t ＜Xj時，時刻t 系 統(tǒng)必處于工作狀態(tài)，自然有X(t)=1;第2 項中，當(dāng)Xj≤t ≤Xj+Wj+γj，即系統(tǒng)處于修理狀態(tài)，不可能與X(t)=1 同時發(fā)生，故第2 項為0;第3 項如下式推導(dǎo)所示。

從而

4 算例分析

假設(shè)某型導(dǎo)彈武器系統(tǒng)由n 個分系統(tǒng)串聯(lián)而成，通過可靠性分析，建立系統(tǒng)可靠性框圖如圖1所示。

圖1 系統(tǒng)可靠性框圖Fig.1 Reliablity block diagram of system

根據(jù)已有的實驗和專家經(jīng)驗確定，其分系統(tǒng)和測試設(shè)備壽命服從指數(shù)分布，保障延誤時間和維修時間服從一般分布，應(yīng)用三角模糊數(shù)來評定其隨機模糊參數(shù)如表1所示。

表1 不確定參數(shù)的三角模糊表示Tab.1 The uncertain parameters of triangular fuzzy representation

根據(jù)三角模糊數(shù)的相關(guān)定理，可以得出如下不同參數(shù)的α-悲觀值和α-樂觀值(此處用x 表示α):

其中x ∈(0，1］。

根據(jù)定義2 可知，可用度的悲觀值和樂觀值可分別表示為

由定義3 可得出

由算例分析結(jié)果可知，系統(tǒng)的可用度為0.5811，如果不考慮測試和維修設(shè)備故障問題，則β=0，系統(tǒng)的可用度為0.618 0，與經(jīng)過前期MC 仿真得到的0.592 3 比較接近。結(jié)果表明，所建模型能夠很好地反映裝備實際組成和運行情況，該模型和算法適用于導(dǎo)彈武器系統(tǒng)可用度的分析。

5 結(jié) 語

本文在隨機模糊參數(shù)的背景下，研究了導(dǎo)彈武器系統(tǒng)可用度評估模型，通過分析可得到如下結(jié)論:

1)建立的馬爾可夫更新過程模型考慮到技術(shù)陣地設(shè)備部分戰(zhàn)損和保障延誤時間等因素，符合實際使用維修和作戰(zhàn)要求。

2)所建可用度模型中存在模糊因素，引入隨機模糊理論，采用三角模糊數(shù)量化模糊信息(分布類型和參數(shù))能合理、有效地解決這一問題。

3)為復(fù)雜的導(dǎo)彈武器系統(tǒng)可用度研究提供了一條技術(shù)途徑，對改進武器系統(tǒng)，制導(dǎo)武器系統(tǒng)的維護和管理有促進作用。

［1］黃奇.兩不同部件串聯(lián)可修系統(tǒng)的可用度的一個新的計算方法［J］.數(shù)學(xué)理論與應(yīng)用，2004，24(2):123-126.

［2］張幗奮，武洪萍.出現(xiàn)保障延誤時間的可修系統(tǒng)穩(wěn)態(tài)可用度的評定［J］.浙江大學(xué)學(xué)報，2006，33(3):268-271.

［3］沈吉鋒，張永志，等.基于馬爾科夫更新過程的偵查系統(tǒng)可靠性分析［J］.兵工自動化，2010，29(3):7-10.

［4］劉福勝，吳緯，單志偉.基于馬爾科夫更新過程的裝甲裝備使用可用度模型［J］.裝甲兵工程學(xué)院學(xué)報，2010，24(5):15-17.

［5］楊繼坤，徐廷學(xué).基于更新過程的導(dǎo)彈武器系統(tǒng)可用度分析［J］.兵工自動化，2012，31(2):14-16.

［6］田艷芳，師義民，李鳳.基于隨機模糊理論的結(jié)構(gòu)可靠性分析［J］.模糊系統(tǒng)與數(shù)學(xué)，2008，22(1):163-166.

［7］王磊，劉文.基于隨機模糊參數(shù)的結(jié)構(gòu)模糊可靠性分析模型［J］.北京航空航天大學(xué)學(xué)報，2005，31(4):412-415.

［8］陳海建，滕克難，李波.基于模糊理論的導(dǎo)彈系統(tǒng)貯存可靠性仿真方法研究［J］.彈箭與制導(dǎo)學(xué)報，2011，31(2):33-36.

［9］黃卓，郭波.更新理論推論過程及應(yīng)用［J］.數(shù)學(xué)的實踐與認識，2006，36(9):200-204.

［10］楊倫標(biāo)，高英儀.模糊數(shù)學(xué)［M］.廣州:華南理工大學(xué)出版社，1995.

［11］LIU Y B，QIAO Z，WANG G Y.Fuzzy random reliability of structures based on fuzzy random variables［J］.Fuzzy Sets and Systems，1997，86:345-355.

［12］郭書祥，呂震宇.概率模型含有模糊參數(shù)時的模糊失效概率計算方法［J］.機械強度，2003，25(5):527-529.

［13］MACHERE Y，KOEHN P，SPARROW D.Improving reliability and operational availability of military systems［A］.IEEE Aerospace Conference［C］，2005.3489-3957.

［14］劉鯖潔，陳桂明，劉小方.BP 神經(jīng)網(wǎng)絡(luò)在導(dǎo)彈武器系統(tǒng)使用可用度預(yù)測中的應(yīng)用［J］.戰(zhàn)術(shù)導(dǎo)彈技術(shù)，2011(1):26-28.