亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于經(jīng)驗(yàn)?zāi)B(tài)與奇異值分解的振動源數(shù)估計(jì)方法

        2022-09-20 02:04:04劉維新葉超
        機(jī)床與液壓 2022年10期
        關(guān)鍵詞:振源轉(zhuǎn)輪減速器

        劉維新,葉超

        (中國工程物理研究院,機(jī)械制造工藝研究所,四川綿陽 621900)

        0 前言

        轉(zhuǎn)輪模塊是某大型光學(xué)實(shí)驗(yàn)裝置空間濾波器的核心部件,對裝置實(shí)現(xiàn)預(yù)定功能具有重要作用,其可靠性問題備受關(guān)注。轉(zhuǎn)輪模塊由電機(jī)、減速器、聯(lián)軸器、圓光柵、轉(zhuǎn)輪等部件組成,運(yùn)行在真空環(huán)境中,對可靠性的要求極高。

        本文作者在某大型光學(xué)實(shí)驗(yàn)裝置部件可靠性實(shí)驗(yàn)時發(fā)現(xiàn),由于機(jī)械結(jié)構(gòu)緊湊、低速重載、工況轉(zhuǎn)換頻繁、沖擊頻率高,轉(zhuǎn)輪模塊中的行星減速器是其主要故障源,是影響裝置可靠性的關(guān)鍵因素。

        目前,行星減速器故障診斷領(lǐng)域已有較多學(xué)者取得良好的研究成果,但大多針對均勻載荷和伺服電機(jī)驅(qū)動下的減速器,對非均勻載荷和步進(jìn)電機(jī)驅(qū)動下的行星減速器的研究較少。由于振動傳感器只能布置在行星減速器箱體外,行星減速箱內(nèi)一系列激勵源產(chǎn)生的振動響應(yīng)與干擾、噪聲等耦合在一起,振動傳感器測量得到的是統(tǒng)計(jì)意義上的耦合源信號,很難直接測量得到單一振動激勵源。非均勻載荷和步進(jìn)電機(jī)推力輸出特性造成的干擾與轉(zhuǎn)輪模塊行星減速器故障特征耦合在一起難以剝離,實(shí)現(xiàn)振動監(jiān)測信號內(nèi)激勵源(振源)與外激勵源(干擾)的盲源分離,是實(shí)現(xiàn)行星減速器故障檢測與辨識的重要前提。

        傳統(tǒng)盲源分離中,為計(jì)算簡單,一般假設(shè)源的數(shù)量等于觀測傳感器數(shù),則混合/解耦矩陣是方陣,容易估計(jì)。然而,在文中研究對象轉(zhuǎn)輪模塊行星減速器中,這種假設(shè)無法成立,因?yàn)槭馨惭b位置和空間限制,傳感器數(shù)量只能是有限的,而源的數(shù)量很容易超過傳感器數(shù)量。

        為實(shí)現(xiàn)轉(zhuǎn)輪模塊行星減速器振動監(jiān)測信號盲源分離,需準(zhǔn)確估計(jì)內(nèi)部激勵振源數(shù)。針對振動傳感器數(shù)小于本底振源數(shù)的源數(shù)估計(jì)問題,本文作者提出一種基于經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition,EMD)的虛擬通道擴(kuò)展方法。通過經(jīng)驗(yàn)?zāi)B(tài)分解得到的固有模態(tài)函數(shù)構(gòu)建振動信號觀測矩陣,以擴(kuò)充振動傳感器觀測通道數(shù)量。針對奇異值分解(Singular Value Decomposition,SVD)特征矩陣中噪聲和數(shù)據(jù)觀測誤差以及不確定性導(dǎo)致的源數(shù)估計(jì)不準(zhǔn)確問題,提出一種基于類內(nèi)散度與類間距離比值優(yōu)化的聚類分析方法。通過對奇異值分解后特征值矩陣中對角線特征值的聚類分析,獲得盲源數(shù)估計(jì)結(jié)果。通過仿真實(shí)驗(yàn)和模擬實(shí)驗(yàn),驗(yàn)證文中所提振源數(shù)估計(jì)方法的有效性。

        1 基于經(jīng)驗(yàn)?zāi)B(tài)分解的虛擬通道擴(kuò)展方法

        經(jīng)驗(yàn)?zāi)B(tài)分解能夠把復(fù)雜信號分解為一系列內(nèi)稟函數(shù)(Intrinsic Mode Function,IMF)之和,每個IMF分量代表原信號的一個固有信號模態(tài),且包含了不同的特征時間尺度,讓原信號特征以不同的分辨率表現(xiàn)出來。由于每個IMF所包含的頻率成分不僅與采樣頻率有關(guān),還隨著信號的變化而變化,所以EMD方法是一種自適應(yīng)的信號處理方法,適用于時變、非平穩(wěn)信號分解。

        本文作者分析轉(zhuǎn)輪模塊振動監(jiān)測系統(tǒng)的特點(diǎn)和振動機(jī)制,分別選取行星減速器2個振動傳感器監(jiān)測信號的前5個IMF,將原有的2個振動觀測通道擴(kuò)展為10個振動觀測通道。

        虛擬通道擴(kuò)展方法具體實(shí)現(xiàn)過程如下:

        (1)采集轉(zhuǎn)輪模塊行星減速器箱體表面振動信號();

        (2)辨識信號()上所有局部極大值和極小值點(diǎn),用三次樣條曲線分別連接所有局部極大值點(diǎn)和極小值點(diǎn),形成上下包絡(luò)線;

        (3)計(jì)算上下包絡(luò)線的均值,然后計(jì)算=()-,并按條件1、條件2判斷是否滿足IMF的條件,若滿足則為振動信號()的1個IMF分量;

        條件1:具有相同數(shù)量的極值點(diǎn)和過零點(diǎn),或最多相差一個;

        條件2:任意時刻極大值和極小值的包絡(luò)線關(guān)于時間軸局部對稱,平均值為0;

        (4)若不滿足IMF的條件,則將作為原始數(shù)據(jù),重復(fù)步驟(1)—(3)次,直到1滿足IMF的條件,記為=1,為振動信號()的1個IMF分量;

        (5)將IMF分量從振動信號()中分離出來,得到=()-,將作為原始數(shù)據(jù)重復(fù)步驟(1)—(4),重復(fù)循環(huán)次,得到個IMF分量,即得到=-;…;=-1-;

        (7)選取1路振動傳感器經(jīng)驗(yàn)?zāi)B(tài)分解后的前5個IMF,構(gòu)建該路傳感器新的5路虛擬通道觀測矩陣;

        (8)循環(huán)結(jié)束。

        在經(jīng)驗(yàn)?zāi)B(tài)分解過程中,原振動信號()相鄰2個極大值點(diǎn)或極小值點(diǎn)之間的時間跨度被稱為時間特征尺度,反映的是振動信號不同模態(tài)的局部振動特性。經(jīng)驗(yàn)?zāi)B(tài)分解過程就是按照特征時間尺度從小到大的順序依次分離各IMF,包含原振動信號從高到低不同頻率段的振動成分,也可以把經(jīng)驗(yàn)?zāi)B(tài)分解方法理解為是一組濾波器。

        2 基于奇異值分解的盲源數(shù)估計(jì)方法

        轉(zhuǎn)輪模塊行星減速器可視為一個多輸入多輸出系統(tǒng)(MIMO)。行星減速器中個振源(=1,…,)作為MIMO系統(tǒng)的輸入,經(jīng)過從源到觀測的頻率響應(yīng)函數(shù)(FRF)矩陣={}(=1,…,;=1,…,)的作用,在振動傳感器處混合形成系統(tǒng)輸出(=1,…,)。各源分量之間可能相關(guān)也可能不相關(guān)。

        由譜分析理論,個振源的自譜和互譜可構(gòu)成一個×的譜矩陣,譜矩陣的線性無關(guān)行(列)的數(shù)目,即矩陣的秩,即為振動源中的不相關(guān)源數(shù)。譜矩陣通常不可能直接觀測,它與傳感器觀測矩陣的關(guān)系為=。由線性代數(shù)相關(guān)定理可知,如果的秩rank()=,則rank()≥rank(),工程應(yīng)用中可近似認(rèn)為rank()=rank(),其中:分別對應(yīng)的頻率響應(yīng)矩陣。因此,工程中可以利用傳感器觀測信號的譜矩陣的秩估計(jì)不相關(guān)振源數(shù)。

        奇異值分解是一種對矩陣的正交化分解方法,是一種數(shù)值穩(wěn)定的矩陣求秩方法。一個×矩陣的SVD表達(dá)為=,其中:是一個由矩陣的正交特征向量組成的×對角矩陣;是一個由矩陣的奇異值組成的×對角矩陣;是一個由矩陣的正交特征向量組成的×矩陣。

        理論上,譜矩陣經(jīng)過奇異值分解后,非零的奇異值數(shù)目即為待求的秩。然而,對于實(shí)際采集的傳感器觀測數(shù)據(jù),譜矩陣的奇異值大多不會出現(xiàn)為零的情況,原因是工程應(yīng)用中采集的傳感器觀測信號,由于有限數(shù)據(jù)采樣精度及奇異值分解計(jì)算的誤差累積效應(yīng),觀測信號譜矩陣的獲得涉及對自譜和互譜密度的估計(jì),導(dǎo)致奇異值分解的非零誤差。一般來說,由數(shù)模轉(zhuǎn)換引入的誤差占比較大。

        通常,工程應(yīng)用中采集到的傳感器觀測信號奇異值分解后得到的奇異值呈現(xiàn)大小不等的兩端分布,其中較大的奇異值與待估計(jì)的振源有對應(yīng)關(guān)系。文獻(xiàn)[14]中建議采用閾值方法對不同奇異值進(jìn)行分類,即設(shè)置一固定閾值,大于該閾值的值認(rèn)為是有效的奇異值,對應(yīng)某一振源;小于該閾值的值認(rèn)為是無效的奇異值,對應(yīng)誤差或噪聲。然而,實(shí)際應(yīng)用中,閾值的設(shè)定必須綜合考慮信號處理與參數(shù)估計(jì)誤差對譜矩陣奇異值分解的影響,而影響譜估計(jì)精度的參數(shù)和因素較多且相互間并無明確的函數(shù)關(guān)系,導(dǎo)致閾值的合理設(shè)定較為困難。

        針對上述奇異值分解特征矩陣中噪聲和數(shù)據(jù)觀測誤差以及不確定性導(dǎo)致的源數(shù)估計(jì)不準(zhǔn)確問題,提出一種基于類內(nèi)散度與類間距離比值優(yōu)化的聚類分析方法。通過對奇異值分解后特征值矩陣中對角線特征值的聚類分析,獲得振動源數(shù)估計(jì)結(jié)果。所提方法的實(shí)施流程如圖1所示。

        圖1 所提振源數(shù)估計(jì)方法流程

        所提振源數(shù)估計(jì)方法的具體實(shí)現(xiàn)過程如下:

        (1)采集轉(zhuǎn)輪模塊減速箱箱體表面90°夾角布置的2個振動傳感器的振動信號()、();

        (2)采用第1節(jié)所述的基于經(jīng)驗(yàn)?zāi)B(tài)分解的虛擬通道擴(kuò)展方法,得到振動信號的固有模態(tài);采用每個單獨(dú)振動信號的前5個固有模態(tài)擴(kuò)展振動觀測信號矩陣;

        (3)計(jì)算振動觀測信號矩陣中各固有模態(tài)的自/互相關(guān)系數(shù),構(gòu)建自/互相關(guān)系數(shù)矩陣:

        其中:IMF,IMF,代表原始振動信號的第個固有模態(tài)函數(shù)與第個固有模態(tài)函數(shù)的互相關(guān)系數(shù),、取值范圍均為[1,2,3,4,5]。

        (4)根據(jù)所構(gòu)建的自/互相關(guān)系數(shù)矩陣,分別計(jì)算矩陣和矩陣的特征值和特征向量;計(jì)算得到的矩陣的特征向量為矩陣奇異值分解的左奇異向量,計(jì)算得到的矩陣的特征向量為矩陣奇異值分解的右奇異向量,特征值組成的對角線矩陣為矩陣的奇異值矩陣,即=,rank()=10;

        (5)以對角矩陣的10個對角線特征元素作為均值聚類的對象,基于類內(nèi)散度與類間距離的比值,即用DB指標(biāo)來評價聚類效果,確定最優(yōu)聚類數(shù);

        (6)最優(yōu)聚類數(shù)減1(噪聲干擾形成的奇異值矩陣特征元素聚集為一類)即為振動信號盲源數(shù)估計(jì)結(jié)果。

        3 仿真實(shí)驗(yàn)驗(yàn)證

        為驗(yàn)證文中所提振動信號盲源數(shù)估計(jì)方法的有效性,采用仿真實(shí)驗(yàn)驗(yàn)證盲源數(shù)估計(jì)結(jié)果是否正確。仿真信號的振源已知,便于評價算法的有效性。構(gòu)造仿真信號如下:

        =sin(2π)+sin(2π)+sin(2π)+0.1

        (1)

        式中:=250;=100;=50;為隨機(jī)干擾,范圍為[0,1)。由式(1)可看出該仿真信號具有3個振源,其時域波形如圖2所示。

        圖2 仿真信號時域波形

        原仿真信號觀測傳感器數(shù)量為1,由盲源分離原理可知,最多只能估計(jì)出一個振動源。為準(zhǔn)確估計(jì)振動源,采用文中第1節(jié)所述的基于經(jīng)驗(yàn)?zāi)B(tài)分解的虛擬通道擴(kuò)展方法,擴(kuò)展傳感器觀測矩陣,得到仿真信號固有模態(tài)分解結(jié)果如圖3所示。

        圖3中IMF8為殘余分量,IMF1~I(xiàn)MF7為仿真信號固有模態(tài)。采用這7個固有模態(tài)函數(shù)擴(kuò)展原觀測傳感器通道,可擴(kuò)展為7通道觀測信號。采用文中第2節(jié)所述的基于奇異值分解的源數(shù)估計(jì)方法,構(gòu)建自/互相關(guān)系數(shù)矩陣,然后對矩陣進(jìn)行奇異值分解,得到奇異值矩陣。對奇異值矩陣的7個對角線特征值進(jìn)行聚類分析,基于類內(nèi)散度與類間距離的比值來評價聚類效果,結(jié)果如圖4所示。

        圖3 仿真信號經(jīng)驗(yàn)?zāi)B(tài)分解結(jié)果

        由DB指標(biāo)的定義可知,DB指標(biāo)越小,聚類效果越好,但DB指標(biāo)不應(yīng)為0。因此,由圖4可判斷仿真信號觀測矩陣7個奇異值特征最優(yōu)分類數(shù)為4,分類結(jié)果如圖5所示。

        圖4 仿真信號DB指標(biāo)聚類效果 圖5 仿真奇異值矩陣特征值聚類結(jié)果

        由圖5可以看出:仿真信號奇異值矩陣對角線特征值聚類為4類;有3個聚類結(jié)果分別對應(yīng)不同的振動源,第4個聚類結(jié)果數(shù)值接近為0,對應(yīng)噪聲和數(shù)據(jù)觀測誤差等因素;文中方法振源辨識結(jié)果與仿真信號真實(shí)情況一致。結(jié)果表明:對于仿真信號,文中所提方法可有效實(shí)現(xiàn)振動源數(shù)盲估計(jì)。

        4 故障模擬實(shí)驗(yàn)對比驗(yàn)證

        為進(jìn)一步驗(yàn)證文中所提振源數(shù)估計(jì)方法的有效性,采用某大型科學(xué)裝置上的轉(zhuǎn)輪模塊原型樣機(jī)開展實(shí)驗(yàn)驗(yàn)證。轉(zhuǎn)輪模塊結(jié)構(gòu)與構(gòu)建的原型樣機(jī)狀態(tài)監(jiān)測系統(tǒng)分別如圖6、圖7所示。

        圖6 轉(zhuǎn)輪模塊結(jié)構(gòu)

        圖7 轉(zhuǎn)輪模塊狀態(tài)監(jiān)測系統(tǒng)

        該監(jiān)測系統(tǒng)由轉(zhuǎn)輪模塊、8個振動傳感器、1臺DASP-V11振動信號采集分析儀組成。8個振動傳感器分別布置在前支撐軸承、后支撐軸承、行星減速器、步進(jìn)電機(jī)上,其中,本文作者重點(diǎn)研究的行星減速器上呈90°夾角布置有2個振動傳感器。為獲取減速器故障狀態(tài)下的振動信號,通過在減速器太陽輪、行星輪表面進(jìn)行加工的方式模擬減速器齒輪故障,并獲取對應(yīng)工況下的減速器振動監(jiān)測信號。實(shí)驗(yàn)過程:轉(zhuǎn)輪旋轉(zhuǎn)360°后反向往復(fù)運(yùn)動,采樣頻率1 000 Hz。齒輪故障模擬效果如圖8所示。

        圖8 一級行星輪崩齒故障模擬效果

        以崩齒故障為例,驗(yàn)證文中方法的有效性。一級行星齒輪崩齒故障振動信號監(jiān)測結(jié)果如圖9所示。

        圖9 崩齒故障行星減速器振動監(jiān)測信號

        行星減速器具有2個振動觀測傳感器,根據(jù)盲源分離原理分析,最多只能估計(jì)出2個振動源。為準(zhǔn)確估計(jì)振動源,采用第1節(jié)所述的基于經(jīng)驗(yàn)?zāi)B(tài)分解的虛擬通道擴(kuò)展方法,擴(kuò)展傳感器觀測矩陣??紤]到算法的復(fù)雜度和運(yùn)算速度,選取每個振動傳感器監(jiān)測信號的前5個固有模態(tài)函數(shù)擴(kuò)展觀測信號,結(jié)果如圖10所示。

        圖10 崩齒故障振動信號經(jīng)驗(yàn)?zāi)B(tài)分解結(jié)果

        采用垂向和橫向振動信號的前5個固有模態(tài)IMF1~I(xiàn)MF5擴(kuò)展原觀測傳感器通道,可擴(kuò)展為10通道觀測信號。采用文中第2節(jié)所述的基于奇異值分解的源數(shù)估計(jì)方法,得到奇異值矩陣。對奇異值矩陣的10個對角線特征值進(jìn)行聚類分析,基于類內(nèi)散度與類間距離的比值評價聚類效果,結(jié)果如圖11所示。

        由圖11可判斷仿真信號觀測矩陣10個奇異值特征最優(yōu)分類數(shù)為5,分類結(jié)果如圖12所示。

        圖11 振動信號DB指標(biāo)聚類效果 圖12 振動信號奇異值矩陣特征值聚類結(jié)果

        由圖12可以看出:奇異值矩陣對角線特征值聚類為5類,其中5個聚類結(jié)果分別對應(yīng)不同的振動源,第5個聚類結(jié)果數(shù)值接近0,對應(yīng)噪聲和數(shù)據(jù)觀測誤差等因素,即文中方法估計(jì)轉(zhuǎn)輪模塊崩齒故障振動源有4個。而采用傳統(tǒng)的閾值方法處理奇異值矩陣對角線特征時,特征閾值如設(shè)為0.4,則小于0.4的對角線特征值都被歸類于噪聲和觀測誤差而不被認(rèn)為是振動源。傳統(tǒng)方法振源估計(jì)結(jié)果為1,顯著不符合轉(zhuǎn)輪模塊的物理特性,并且該特征閾值只能人為估計(jì),很難準(zhǔn)確設(shè)定。結(jié)果表明,文中所提振動信號盲源數(shù)估計(jì)方法效果優(yōu)于傳統(tǒng)閾值方法。

        5 結(jié)論

        本文作者針對轉(zhuǎn)輪模塊行星減速器振動源數(shù)盲估計(jì)問題,提出一種基于經(jīng)驗(yàn)?zāi)B(tài)分解與奇異值分解相結(jié)合的振動信號盲源數(shù)估計(jì)方法。針對振動傳感器數(shù)小于本底振源數(shù)的源數(shù)估計(jì)問題,通過經(jīng)驗(yàn)?zāi)B(tài)分解得到的固有模態(tài)函數(shù)構(gòu)建振動信號觀測矩陣,以擴(kuò)充振動傳感器觀測通道數(shù)量。針對奇異值分解特征矩陣中噪聲和數(shù)據(jù)觀測誤差以及不確定性導(dǎo)致的源數(shù)估計(jì)不準(zhǔn)確問題,通過對奇異值分解后特征值矩陣中對角線特征值的聚類分析,獲得盲源數(shù)估計(jì)結(jié)果。通過仿真實(shí)驗(yàn)和故障模擬實(shí)驗(yàn)驗(yàn)證所提方法的有效性。結(jié)果表明:利用傳統(tǒng)直接采用原始振動信號進(jìn)行源數(shù)估計(jì)的方法只能計(jì)算得到2個振動源;傳統(tǒng)的基于閾值的奇異值分解方法難以準(zhǔn)確估計(jì)振動源數(shù);所提的基于虛擬通道擴(kuò)展與奇異值分解相結(jié)合的方法可準(zhǔn)確估計(jì)振動源數(shù),效果優(yōu)于傳統(tǒng)方法。

        猜你喜歡
        振源轉(zhuǎn)輪減速器
        Effects of Prandtl number in two-dimensional turbulent convection?
        詞語大轉(zhuǎn)輪
        ——“AABC”和“無X無X”式詞語
        寺廟里有座大書架——神奇的轉(zhuǎn)輪藏
        駝峰第三制動位減速器夾停鉤車問題的改進(jìn)
        考慮振源相互作用的船舶甲板減振方法
        船海工程(2019年1期)2019-03-04 01:58:04
        低密度超音速減速器
        軍事文摘(2018年24期)2018-12-26 00:57:40
        我國第一臺分半鑄造的30萬千瓦水輪發(fā)電機(jī)轉(zhuǎn)輪實(shí)制成功
        蘭臺世界(2017年12期)2017-06-22 14:17:25
        一種小型化低噪聲微波本振源設(shè)計(jì)
        電子制作(2017年9期)2017-04-17 03:00:36
        基于ANSYS Workbench 的ATB260 減速器箱體模態(tài)分析
        振源傳導(dǎo)樁隔振系統(tǒng)作用機(jī)理與性能
        又湿又紧又大又爽a视频国产| 久久国产精品99精品国产987| 亚洲中文字幕无线乱码va| 亚洲性日韩一区二区三区| 欧美老妇交乱视频在线观看| 国产熟妇人妻精品一区二区动漫| 国产鲁鲁视频在线播放| 国产免费一区二区av| 日本精品久久不卡一区二区| 日本阿v片在线播放免费| 免费人成年小说在线观看| 中文字幕午夜AV福利片| 国产一区二区三区的区| 东京热无码av一区二区| 在教室伦流澡到高潮hnp视频 | 18成人片黄网站www| 99久久国内精品成人免费 | 在线天堂av一区二区| 国产精品永久免费| 欧美一片二片午夜福利在线快| 精品久久久亚洲中文字幕| 国产av一卡二卡日韩av| 色噜噜久久综合伊人一本| 亚洲中文久久精品无码ww16| 国产呦系列视频网站在线观看 | 国产做床爱无遮挡免费视频| 91国产熟女自拍视频| 日本真人边吃奶边做爽电影| 性生交大片免费看淑女出招 | 欧美a级毛欧美1级a大片免费播放| 1769国产精品短视频| 亚洲一区二区女优视频| 男人吃奶摸下挵进去啪啪软件| 女厕厕露p撒尿八个少妇| 亚洲欧洲日产国码无码| 国产精品 视频一区 二区三区 | 香蕉久久福利院| 亚洲不卡电影| 麻豆视频av在线观看| 成人精品天堂一区二区三区| 尤物yw无码网站进入|