考慮不完全維修的風機齒輪箱優化檢修策略(2)
作者:石家莊風機 日期:2015-7-30 瀏覽:1435
風機廠石家莊風機廠石家莊風機石家莊市風機廠石家莊風機維修石家莊風機銷售
風機齒輪箱的維修方式針對石家莊風機廠齒輪箱這一復雜多部件可修系統而言,在壽命期限內發生故障,不需要對整個設備進行更換,而僅更換或維修部分部件,來維持所預定的功能,故在其整個壽命期內將發生多次故障及維修活動。根據維修類型和維修程度的不同,維修可分為如下幾種 [15-16](1)完全維修。指設備修復如新,如齒輪箱整機更換就屬于完全維修。其維修因子:M =0。(2)最小維修。指設備修復如舊,如齒輪箱零部件加固、潤滑油更換等屬于最小維修。其維修因子 M =1。(3)不完全維修。指設備維修后,功能得以恢復,設備狀態介于完全維修和最小維修之間。如齒輪箱零部件更換或維修屬于不完全維修。其維修因子 0< M <1。對于石家莊風機廠齒輪箱常見的故障,如輪齒損壞、軸承磨損等故障類型而言,通常僅對其部分部件進行更換或維修而不進行整機更換。故雖然其功能得以恢復,但其健康狀態通常處于“修復如新”與“修復如舊”之間,即不完全維修(0<M<1)。因此,本文針對不完全維修建模研究維修活動與實際情況更加吻合。
在 PIM 中,每次故障事件是隨機的,故障事件通常被視為一個非齊次 Possion 過程,則系統發生失效的強度函數為T0(1)其中: t 是齒輪箱的運行時間; ( )vt為齒輪箱在 t 時刻的強度函數;0vt為基本強度函數Nt為設備在區間(0, t )內失效次數; ( ) Z t 是系統在時刻t的協變量;γ為協變量回歸參數。比例強度模型包括全參數型和半參數型,當失效規律已知(0 ( )v t 已知)時,為全參數型PIM,能同時估計出基本強度函數的參數和協變量的回歸系數,參數估計的準確率相對較高;當失效規律未知(0 ( )v t 未知)時,采用半參數型PIM。1994年Kumar等綜述了對PIM參數的各種估計方法,以基本強度函數服從威布爾分布的PIM的估計結果最為準確。 研究表明, 齒輪箱的故障時間服從威布爾分布,其強度函數為(2)其中: β 為形狀參數; η 為尺度參數;( ) Z t 為系統在 t 時刻的協變量; γ 為協變量的回歸參數。協變量可分為外部協變量和內部協變量。外部協變量包括石家莊風機廠齒輪箱的運行環境溫度、工作負荷等,而內部協變量包括反映齒輪箱故障征兆的振動數據、溫度數據等??紤]維修效應對強度函數的影響,既可以將維修效應作為衰減因子得到虛擬壽命過程,也可以將其作為協變量疊加到強度函數中,本文采用后者。
1.3 石家莊風機廠齒輪箱 PIM 的參數估計利用石家莊風機廠齒輪箱故障數據對PIM 的參數進行極大似然估計。首先驗證故障數據是否符合威布爾分布以確定模型選擇的正確性,然后構建似然函數,并在 Matlab 中應用 Newton-Raphson 迭代算法解非線性方程組,計算出估計值,具體步驟如下。第一步:樣本分布的檢驗。采用最常見的威布爾分布檢驗圖,檢驗圖反映了威布爾分布的基本特征,是用于檢驗樣本數據是否符合威布爾分布規律的最直觀方法。設威布爾分布的可靠度函數為( ) 1對上式兩邊分別求兩次對數(3)可知,對于石家莊風機廠齒輪箱故障時刻 t ,繪制 lnt 與 ln( ln( ( )))R t −之間的函數曲線。 若樣本服從威布爾分布,則該曲線應近似地擬合為一條斜率大于0的直線, 直線的斜率即為形狀參數的估計值, 如圖2所示,其故障數據服從威布爾分布。第二步:似然函數的構造。(4)式中: n 為樣本總數; q 為觀測到失效時間的樣本個數;iδ 為截尾狀態量,值為0表示失效為齒輪箱的概率密度函數i kR t β η γ 為可靠度函數。第三步:PIM模型參數估計。對似然函數兩邊求對數后,并對 β , η ,kγ 分別求一階偏導,且置0。
然后,采用Newton-Raphson法進行數值計算,求解出各參數值,并代入模型(2),得ˆ1ˆˆ ( ) ( ) exp( ( ))ˆˆk ktv t Z tββγη η−= ⋅∑式中,ˆβ 、 ˆ η 、 ˆ k γ 為 β 、 η 、kγ 的估計值。2 最優維修決策利用比例強度模型,綜合監測數據、歷史故障數據及維修活動,確定設備在運行時刻 t 的強度函數 ( )v t 。將強度函數與其維修閾值*v 進行比較以決定是否實施維修。 當 t 時刻的強度*( ) v t v > 時, 進行維修;當*( ) v t v ≤ 時,不進行維修。本文應用物理規劃方法確定最優維修閾值*v ,具體詳見下一節。石家莊風機廠齒輪箱實施維修的決策條件 (5)由于強度函數中指數處理復雜且不直觀,故對式(5)取對數,進行化簡整理(6)式(6)描述了石家莊風機廠齒輪箱運行時間 t 與狀態變量之間的關系。對于特定的維修閾值*v ,以運行時間 t 為橫坐標, 以*ln( ) ( 1)ln( ) 為縱坐標得一條對 數 曲 線 , 即 最 優 維 修 決 策 控 制 限 , 為縱坐標,得到置信上下限。將 t 時刻觀測到的狀態變量代入T( ) Z t γ ⋅∑中,當觀測點在置信下限以下時,保持設備繼續運截尾,為1表示行;當觀測點在置信上下限內,密切監視并隨時進行維修;當觀測點在置信上限以上時,立刻停機維修。