• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于容錯(cuò)齒磁通的五相永磁電機(jī)匝間短路故障診斷

    2023-07-26 05:29:40尹佐生王天乙宋志翌
    電工技術(shù)學(xué)報(bào) 2023年14期
    關(guān)鍵詞:繞線匝間匝數(shù)

    黃 建 尹佐生 王天乙 宋志翌 隋 義

    基于容錯(cuò)齒磁通的五相永磁電機(jī)匝間短路故障診斷

    黃 建1尹佐生1王天乙1宋志翌1隋 義2

    (1. 北京自動(dòng)化控制設(shè)備研究所 北京 100074 2. 哈爾濱工業(yè)大學(xué)電氣工程及自動(dòng)化學(xué)院 哈爾濱 150001)

    匝間短路故障是一種常見(jiàn)的電機(jī)繞組故障,具有危害大、易擴(kuò)散、診斷難的特點(diǎn),是多相容錯(cuò)永磁同步電機(jī)研究過(guò)程中關(guān)注的重點(diǎn)。該文以一臺(tái)高故障隔離能力五相永磁同步電機(jī)為基礎(chǔ),對(duì)基于齒磁通的故障診斷方法進(jìn)行了研究。建立了五相永磁同步電機(jī)定子磁路模型,給出了匝間短路故障對(duì)繞線齒及容錯(cuò)齒磁通的影響規(guī)律。研究了不同繞制方式下探測(cè)線圈感應(yīng)電動(dòng)勢(shì)隨電機(jī)運(yùn)行狀態(tài)的變化規(guī)律,對(duì)各繞制方案的故障識(shí)別能力、故障靈敏度、復(fù)雜程度等進(jìn)行了綜合評(píng)價(jià),提出基于容錯(cuò)齒磁通的匝間短路故障診斷方法,給出基于故障標(biāo)志的故障識(shí)別和定位方法。研究表明,當(dāng)電機(jī)發(fā)生匝間短路故障時(shí),故障相值明顯降低,剩余健康相值接近于1,根據(jù)這一特征可以對(duì)故障進(jìn)行準(zhǔn)確識(shí)別和定位。該方法不依賴傅里葉變換、不受負(fù)載和電機(jī)轉(zhuǎn)速影響,具有準(zhǔn)確度高、魯棒性強(qiáng)的優(yōu)點(diǎn)。該文進(jìn)一步分析了槽口寬度、容錯(cuò)齒長(zhǎng)度對(duì)故障敏感度和電機(jī)性能的影響,給出兼顧高低匝數(shù)短路故障診斷靈敏度和電機(jī)性能的尺寸選取方法。

    五相永磁同步電機(jī) 匝間短路 探測(cè)線圈 故障診斷

    0 引言

    相比傳統(tǒng)三相永磁同步電機(jī),多相容錯(cuò)永磁同步電機(jī)具備相冗余特性,故障下通過(guò)剩余相電流重構(gòu)可以保證電機(jī)繼續(xù)運(yùn)行,極大地提高了系統(tǒng)的可靠性和容錯(cuò)性,在航空航天、電動(dòng)汽車等對(duì)可靠性要求較高的場(chǎng)合具有良好的應(yīng)用前景[1-4]。

    繞組開(kāi)路和短路故障是電機(jī)的易發(fā)故障,對(duì)于多相容錯(cuò)永磁電機(jī)而言,想要做到開(kāi)路、短路等故障情況下的快速響應(yīng)和容錯(cuò)運(yùn)行,電機(jī)系統(tǒng)本身需要具備能夠識(shí)別繞組故障的快速準(zhǔn)確故障診斷功能[5-7]。電機(jī)繞組故障會(huì)導(dǎo)致轉(zhuǎn)矩、氣隙磁通及電樞電壓、電流、阻抗等物理量的畸變或不平衡,進(jìn)而引起額外的振動(dòng)噪聲,現(xiàn)有方法一般通過(guò)對(duì)上述物理量的監(jiān)測(cè)實(shí)現(xiàn)故障診斷。其中,繞組開(kāi)路故障的診斷相對(duì)簡(jiǎn)單,當(dāng)電機(jī)某一相繞組出現(xiàn)開(kāi)路故障時(shí),該相電流持續(xù)為零,并且諧波平面電流分量會(huì)出現(xiàn)畸變[8-10],因此可以采取監(jiān)測(cè)各相電流過(guò)零點(diǎn),或檢測(cè)諧波平面電流分量的方法對(duì)開(kāi)路故障進(jìn)行診斷,當(dāng)前研究的重點(diǎn)在于縮短故障診斷周期,提高響應(yīng)速度[11]。

    相比開(kāi)路故障,電機(jī)繞組短路故障危害更大、診斷更為困難,尤其是匝間短路故障,當(dāng)短路匝數(shù)較少時(shí),故障特征不明顯,診斷難度大[12-13]。當(dāng)電機(jī)繞組出現(xiàn)短路故障時(shí),一方面繞組阻抗特性會(huì)發(fā)生變化,影響電壓激勵(lì)下的電流響應(yīng);另一方面會(huì)導(dǎo)致各相電流或電壓不均衡,使幅值、有效值和諧波成分等出現(xiàn)差異,進(jìn)而導(dǎo)致零序電流電壓和旋轉(zhuǎn)坐標(biāo)系下的電流電壓分量出現(xiàn)畸變。文獻(xiàn)[14-15]提出了基于高頻電壓注入的電機(jī)繞組故障診斷方法,在逆變器端施加額外的高頻載波或調(diào)制波電壓,通過(guò)測(cè)量繞組中高頻電流分量對(duì)繞組狀態(tài)進(jìn)行監(jiān)測(cè),該方法不受基波電流頻率的影響,對(duì)電機(jī)運(yùn)行狀態(tài)不敏感,適用于繞組匝間短路故障診斷。文獻(xiàn)[16-17]建立了匝間短路故障數(shù)學(xué)模型,研究指出,匝間短路狀態(tài)下定子電流和dq平面電流分量中會(huì)出現(xiàn)一定的諧波成分,通過(guò)采取快速傅里葉變換(Fast Fourier Transformation, FFT)提取特定次諧波可以實(shí)現(xiàn)對(duì)匝間短路故障的診斷,但這類基于FFT的故障特征提取方法易受電機(jī)工況影響,多工況下的匝間短路故障診斷存在難度。文獻(xiàn)[18-19]指出,匝間短路故障會(huì)破壞電機(jī)各相電壓的對(duì)稱性,使零序電壓發(fā)生畸變,通過(guò)搭建額外的三相阻抗平衡網(wǎng)絡(luò)檢測(cè)中性點(diǎn)電壓可以實(shí)現(xiàn)對(duì)匝間短路故障的診斷,該方法適用于星形聯(lián)結(jié)電機(jī),且需要附加額外的電壓檢測(cè)電路。在電機(jī)定子齒部繞制探測(cè)線圈也是一種常用的匝間短路故障診斷方法,具有準(zhǔn)確度高、響應(yīng)快等優(yōu)勢(shì)[20]。文獻(xiàn)[21-22]中電機(jī)每個(gè)定子繞線齒上均繞制一個(gè)探測(cè)線圈,各線圈相互串聯(lián),研究表明,電機(jī)未發(fā)生故障時(shí),探測(cè)線圈電壓只包含3次諧波,當(dāng)出現(xiàn)短路故障時(shí),電壓中包含明顯的1次分量,通過(guò)對(duì)該1次分量的監(jiān)測(cè)可以實(shí)現(xiàn)對(duì)匝間短路故障的診斷,但該故障診斷方法同樣需要采用FFT對(duì)故障特征進(jìn)行處理,易受電機(jī)工況影響,并且對(duì)故障位置的判斷比較困難。文獻(xiàn)[23]采用的探測(cè)線圈為各齒獨(dú)立布置,并未對(duì)線圈進(jìn)行串聯(lián),該方法需要更多的檢測(cè)通道,優(yōu)點(diǎn)是對(duì)故障的定位更容易。文獻(xiàn)[21-23]均選用探測(cè)線圈感應(yīng)電動(dòng)勢(shì)低頻分量作為故障特征,對(duì)低匝數(shù)匝間短路故障敏感度較弱,針對(duì)這一問(wèn)題,文獻(xiàn)[24]提出采用逆變器側(cè)PWM供電引起的探測(cè)線圈高頻電壓分量作為故障判據(jù),該方法對(duì)低匝數(shù)匝間短路故障更為敏感,但同樣需要采用FFT對(duì)故障特征進(jìn)行處理,計(jì)算較為復(fù)雜,實(shí)時(shí)性難以保證。

    本文對(duì)基于探測(cè)線圈的五相永磁同步電機(jī)匝間短路故障診斷方法進(jìn)行研究。分析了正常運(yùn)行及匝間短路狀態(tài)下的電機(jī)定子齒磁通分布規(guī)律,為故障診斷提供理論依據(jù)。研究了不同繞制方式下探測(cè)線圈感應(yīng)電動(dòng)勢(shì)隨電機(jī)運(yùn)行狀態(tài)的變化規(guī)律,針對(duì)傳統(tǒng)方法對(duì)低匝數(shù)故障敏感度低、FFT實(shí)時(shí)性差的問(wèn)題,提出了基于容錯(cuò)齒磁通的匝間短路故障診斷方法,給出了故障識(shí)別和定位方法。分析了槽口寬度等關(guān)鍵尺寸對(duì)故障敏感度和電機(jī)性能的影響規(guī)律,給出了兼顧高低匝數(shù)短路故障敏感度和電機(jī)性能的尺寸選取方法。

    1 電機(jī)拓?fù)浣Y(jié)構(gòu)簡(jiǎn)介

    本文所研究的五相永磁同步電機(jī)結(jié)構(gòu)如圖1所示。同一相繞組的兩個(gè)線圈空間相鄰排布,磁通相互耦合,從根源上降低相間磁耦合。同一相的兩個(gè)線圈及不同相的線圈間通過(guò)容錯(cuò)齒實(shí)現(xiàn)物理和熱隔離,不同相線圈間容錯(cuò)齒為槽漏磁通提供路徑,避免一相繞組產(chǎn)生的槽漏磁通與相鄰相繞組耦合,進(jìn)一步阻絕槽漏磁引起的相間磁耦合。此外,還可以通過(guò)調(diào)整容錯(cuò)齒和繞線齒寬度配比提高電機(jī)基波繞組因數(shù)、降低諧波繞組因數(shù),提高電機(jī)功率密度的同時(shí)實(shí)現(xiàn)對(duì)定子磁動(dòng)勢(shì)諧波的抑制,降低轉(zhuǎn)子損耗。該電機(jī)主要結(jié)構(gòu)參數(shù)及關(guān)鍵指標(biāo)見(jiàn)表1。

    表1 五相永磁同步電機(jī)主要結(jié)構(gòu)參數(shù)及關(guān)鍵指標(biāo)

    Tab.1 Main structual parameters of the five-phase permanent-magnet synchronous machine

    2 正常運(yùn)行及匝間短路狀態(tài)下的定子齒磁通分布

    2.1 正常運(yùn)行狀態(tài)

    首先建立正常運(yùn)行狀態(tài)下該電機(jī)主磁路的等效磁路模型,如圖2a所示。定義同相線圈間容錯(cuò)齒為容錯(cuò)齒Ⅰ,不同相線圈間容錯(cuò)齒為容錯(cuò)齒Ⅱ。為了簡(jiǎn)化分析,此處作如下假設(shè):

    (1)鐵心磁導(dǎo)率無(wú)窮大,忽略鐵心磁壓降。

    (2)忽略定子開(kāi)槽影響。

    (3)永磁體磁導(dǎo)率等于空氣磁導(dǎo)率。

    (4)磁通沿徑向方向穿過(guò)氣隙。

    圖2 主磁通等效磁路和簡(jiǎn)化磁路模型

    在圖2a中,A1、A2分別為A相A1、A2齒上線圈產(chǎn)生的磁動(dòng)勢(shì);d0為繞線齒齒頂?shù)睫D(zhuǎn)子鐵心之間的磁阻;d1為容錯(cuò)齒Ⅰ齒頂?shù)睫D(zhuǎn)子鐵心之間的磁阻;d2為容錯(cuò)齒Ⅱ齒頂?shù)睫D(zhuǎn)子鐵心之間的磁阻。A1、A2、d0、d1和d2分別為

    式中,A1A2分別為A1、A2齒上線圈串聯(lián)匝數(shù);0為相電流幅值;ef0、ef1、ef2分別為繞線齒、容錯(cuò)齒Ⅰ、容錯(cuò)齒Ⅱ齒頂?shù)睫D(zhuǎn)子鐵心外表面間徑向距離;t0、t1、t2分別為繞線齒、容錯(cuò)齒Ⅰ、容錯(cuò)齒Ⅱ兩側(cè)槽口中心線間弧長(zhǎng);ef為電機(jī)定轉(zhuǎn)子鐵心軸向長(zhǎng)度;0為相對(duì)磁導(dǎo)率。

    進(jìn)一步對(duì)圖2a中的等效磁路模型進(jìn)行簡(jiǎn)化,將A相繞線齒和容錯(cuò)齒Ⅰ以外的并聯(lián)磁路磁阻等效為0,簡(jiǎn)化后的磁路模型如圖2b所示。在圖2b中,0為A相繞線齒和容錯(cuò)齒Ⅰ以外磁路對(duì)應(yīng)的磁阻并聯(lián)值,即0=d2//d2//d0//d0//d1…,計(jì)算得到0表達(dá)式如式(2)所示;A為定子鐵心軛部到轉(zhuǎn)子鐵心之間的磁動(dòng)勢(shì)。根據(jù)簡(jiǎn)化磁路模型,由節(jié)點(diǎn)電流定律得到式(3)所示方程。

    根據(jù)式(3)可以得到A為

    進(jìn)一步得到A1、A2齒及其他繞線齒和容錯(cuò)齒上的主磁路磁通分別為

    在正常運(yùn)行狀態(tài)下,A1、A2齒上線圈電流相等、匝數(shù)相同,即A1A2、A1+A2=0。此時(shí),定子鐵心軛部到轉(zhuǎn)子鐵心之間的磁動(dòng)勢(shì)A=0,說(shuō)明在正常運(yùn)行狀態(tài)下,除A1、A2齒外,A相線圈在其他所有繞線齒和容錯(cuò)齒上產(chǎn)生的磁通為零,即各相繞組產(chǎn)生的磁場(chǎng)僅分布在該相繞線齒上。

    除了圖2所示的穿過(guò)電機(jī)氣隙的主磁路磁通外,定子齒上還存在經(jīng)由槽口和槽內(nèi)區(qū)域閉合的槽漏磁通。以A相為例,槽漏磁通路徑及等效磁路如圖3所示。由圖3中等效磁路可以得到A相繞線齒、容錯(cuò)齒Ⅰ和容錯(cuò)齒Ⅱ上徑向磁通分別為

    圖3 槽漏磁通等效磁路

    圖4 正常運(yùn)行狀態(tài)下磁通分布

    2.2 匝間短路狀態(tài)

    有限元仿真得到A1≠A2狀態(tài)下僅A相繞組通電時(shí)的磁通分布,如圖5所示??梢钥闯觯珹相繞組產(chǎn)生的磁通不止分布在A相線圈所在繞線齒和容錯(cuò)齒Ⅱ上,A相容錯(cuò)齒Ⅰ上也有徑向磁通分布,此外,其他繞線齒和容錯(cuò)齒上也存在少量徑向磁通。

    圖5 匝間短路狀態(tài)下磁通分布

    匝間短路故障導(dǎo)致繞線齒、容錯(cuò)齒上磁通出現(xiàn)不同程度變化。從式(5)和式(6)中可以看出,A2齒線圈匝數(shù)輕微變化對(duì)A1齒磁通的影響較小,但是會(huì)引起A2齒磁通的變化。A1、A2齒上線圈匝數(shù)變化均會(huì)導(dǎo)致A相容錯(cuò)齒Ⅰ上出現(xiàn)定子側(cè)磁通,并且該磁通由主磁路磁通和槽漏磁通兩部分構(gòu)成,該磁通除了與A1和A2差值有關(guān)外,還受s1影響。

    3 基于容錯(cuò)齒磁通的匝間短路故障診斷方法

    3.1 探測(cè)線圈方案選取

    匝間短路狀態(tài)下故障相繞線齒和容錯(cuò)齒的磁通均會(huì)發(fā)生變化,可以通過(guò)探測(cè)線圈檢測(cè)齒磁通變化實(shí)現(xiàn)對(duì)匝間短路故障的診斷。由于容錯(cuò)齒Ⅱ中的磁通由兩側(cè)不同相線圈的磁通疊加而成,當(dāng)探測(cè)線圈繞制在容錯(cuò)齒Ⅱ上時(shí),難以對(duì)故障發(fā)生位置作出判斷,因此,本文對(duì)繞線齒繞制探測(cè)線圈和容錯(cuò)齒Ⅰ繞制探測(cè)線圈方案進(jìn)行研究,探測(cè)線圈繞制方式示意圖如圖6所示,定義單個(gè)齒上探測(cè)線圈匝數(shù)為T。在繞線齒繞制探測(cè)線圈方案Ⅰ中,五相電機(jī)的10個(gè)定子繞線齒上各繞制一個(gè)T=5匝的探測(cè)線圈;在繞線齒繞制探測(cè)線圈方案Ⅱ中,各相1、2號(hào)繞線齒上的探測(cè)線圈同向串聯(lián);在繞線齒繞制探測(cè)線圈方案Ⅲ中,各相1、2號(hào)繞線齒上的探測(cè)線圈反向串聯(lián);在容錯(cuò)齒繞制探測(cè)線圈方案中,電機(jī)的5個(gè)容錯(cuò)齒Ⅰ上各繞制一個(gè)T=5匝的探測(cè)線圈。

    圖6 探測(cè)線圈繞制方式示意圖

    有限元仿真得到轉(zhuǎn)速為1 000 r/min、相電流幅值為30 A時(shí),電機(jī)正常運(yùn)行狀態(tài)下上述探測(cè)線圈中感應(yīng)電動(dòng)勢(shì)波形如圖7所示。

    為了通過(guò)有限元仿真模擬電機(jī)繞組匝間短路故障,建模時(shí)將電機(jī)A相繞組分為兩部分,剩余健康繞組和短路繞組,如圖8所示。在有限元軟件中將短路繞組設(shè)置為“電壓激勵(lì)”,根據(jù)短路匝數(shù)給定短路繞組電阻值,可以模擬繞組匝間短路情況。

    圖7 正常運(yùn)行狀態(tài)下探測(cè)線圈感應(yīng)電動(dòng)勢(shì)波形

    圖8 匝間短路有限元仿真設(shè)置示意圖

    當(dāng)A相A2齒上1匝線圈發(fā)生短路時(shí),有限元仿真得到轉(zhuǎn)速為1 000 r/min、相電流幅值為30 A時(shí)探測(cè)線圈感應(yīng)電動(dòng)勢(shì)波形,如圖9所示。需要注意的是,匝間短路狀態(tài)下,除了定子齒磁通引起的定子分量,探測(cè)線圈感應(yīng)電動(dòng)勢(shì)中也包含轉(zhuǎn)子磁場(chǎng)引起的轉(zhuǎn)子分量,轉(zhuǎn)子分量?jī)H與轉(zhuǎn)子狀態(tài)有關(guān),不能表征定子繞組狀態(tài),本文不對(duì)其進(jìn)行單獨(dú)研究。

    在正常運(yùn)行狀態(tài)下,不同相探測(cè)線圈感應(yīng)電動(dòng)勢(shì)幅值均相等,當(dāng)電機(jī)A相A2齒上線圈的1匝發(fā)生短路故障時(shí),繞線齒繞制探測(cè)線圈方案Ⅰ中的TA2+-TA2-線圈感應(yīng)電動(dòng)勢(shì)幅值變小,TA1+-TA1-~TE1+-TE1-線圈感應(yīng)電動(dòng)勢(shì)幅值無(wú)明顯差異;繞線齒繞制探測(cè)線圈方案Ⅱ中的TA1+-TA2+~TE1+-TE2+線圈感應(yīng)電動(dòng)勢(shì)幅值無(wú)明顯差異;繞線齒繞制探測(cè)線圈方案Ⅲ中的TA1+-TA2-線圈感應(yīng)電動(dòng)勢(shì)幅值變?。蝗蒎e(cuò)齒繞制探測(cè)線圈方案中的TA+-TA-線圈感應(yīng)電動(dòng)勢(shì)幅值變小。可以看出,繞線齒繞制探測(cè)線圈方案Ⅰ僅能檢測(cè)出發(fā)生在同一繞線齒上的線圈匝間短路故障,需要10組探測(cè)線圈才能實(shí)現(xiàn)對(duì)匝間短路故障的診斷,并且故障特征不明顯;繞線齒繞制探測(cè)線圈方案Ⅱ無(wú)法對(duì)匝間短路故障進(jìn)行檢測(cè)。繞線齒繞制探測(cè)線圈方案Ⅲ和容錯(cuò)齒繞制探測(cè)線圈方案均對(duì)匝間短路故障比較敏感,但繞線齒繞制探測(cè)線圈方案Ⅲ中線圈個(gè)數(shù)較多、方案復(fù)雜,對(duì)槽滿率的影響也更大。

    進(jìn)一步仿真得到轉(zhuǎn)速為1 000 r/min、相電流幅值為30 A時(shí),A相繞組開(kāi)路故障下探測(cè)線圈感應(yīng)電動(dòng)勢(shì)波形,如圖10所示。繞線齒繞制探測(cè)線圈方案Ⅰ、Ⅱ、Ⅲ中探測(cè)線圈感應(yīng)電動(dòng)勢(shì)明顯不均衡;而容錯(cuò)齒繞制探測(cè)線圈方案中各線圈感應(yīng)電動(dòng)勢(shì)幅值相等,可以有效區(qū)分開(kāi)路故障和短路故障,并且對(duì)匝間短路故障的診斷無(wú)需采用FFT,診斷耗時(shí)短、操作簡(jiǎn)單。因此,本文選取容錯(cuò)齒繞制探測(cè)線圈方案對(duì)該電機(jī)的匝間短路故障進(jìn)行診斷。

    3.2 故障識(shí)別和故障位置判斷

    為了定量表征不同狀態(tài)下各探測(cè)線圈感應(yīng)電動(dòng)勢(shì)幅值的變化規(guī)律,定義TA+-TA-~TE+-TE-線圈感應(yīng)電動(dòng)勢(shì)幅值分別為TA、TB、TC、TD、TE。從第2節(jié)中的分析可以看出,未發(fā)生匝間短路故障時(shí),正常相繞組在對(duì)應(yīng)探測(cè)線圈中產(chǎn)生的磁通基本為零,探測(cè)線圈中感應(yīng)電動(dòng)勢(shì)由轉(zhuǎn)子磁場(chǎng)產(chǎn)生,只與電機(jī)轉(zhuǎn)速有關(guān),與轉(zhuǎn)速成正比,可以根據(jù)電機(jī)轉(zhuǎn)速得到TA~TE的正常參考值,記作T0,表達(dá)式為

    定義故障標(biāo)志為T+-T-探測(cè)線圈感應(yīng)電動(dòng)勢(shì)幅值測(cè)量值Ti與參考值T0之比,表達(dá)式為

    當(dāng)電機(jī)未發(fā)生匝間短路故障時(shí),TA+-TA-~TE+-TE-線圈感應(yīng)電動(dòng)勢(shì)幅值相等,均等于參考值,各相均等于1。當(dāng)電機(jī)第相繞組出現(xiàn)匝間短路故障時(shí),T+-T-探測(cè)線圈感應(yīng)電動(dòng)勢(shì)幅值明顯降低,導(dǎo)致<1,其余健康相探測(cè)線圈仍為1,根據(jù)這一特點(diǎn)可以實(shí)現(xiàn)對(duì)故障的識(shí)別和定位。

    有限元仿真得到轉(zhuǎn)速為1 000 r/min、相電流幅值為30 A情況下,各探測(cè)線圈值隨A相短路匝數(shù)的變化趨勢(shì),結(jié)果如圖11所示。從圖11中可以看出,A相出現(xiàn)匝間短路故障后,A<1,隨著短路匝數(shù)的增加,A值呈先減小后增大的趨勢(shì),B~E值基本不變,接近1。由于故障相值不隨短路匝數(shù)的增加單調(diào)變化,所以很難根據(jù)的值判斷短路匝數(shù)。

    圖11 不同短路匝數(shù)情況下ai變化趨勢(shì)

    以A相短路匝數(shù)為1的情況為例,有限元仿真得到轉(zhuǎn)速為1 000 r/min時(shí)各探測(cè)線圈值隨剩余相電流幅值的變化趨勢(shì),結(jié)果如圖12所示。從圖12中可以看出,隨著剩余相電流幅值的增大,A值逐漸減小,B~E值基本不變,接近1。

    圖12 不同相電流幅值情況下ai變化趨勢(shì)

    以A相短路匝數(shù)為1的情況為例,有限元仿真得到剩余相電流幅值為30 A時(shí)各探測(cè)線圈值隨電機(jī)轉(zhuǎn)速的變化趨勢(shì),結(jié)果如圖13所示。從圖13中可以看出,隨著電機(jī)轉(zhuǎn)速的增大,A值呈先減小后增大的趨勢(shì),B~E值基本不變,接近1。當(dāng)電機(jī)轉(zhuǎn)速較低時(shí),A值與其他各相值差異不明顯,說(shuō)明低轉(zhuǎn)速下本文所提出的故障診斷方法靈敏度會(huì)有所降低,但大部分轉(zhuǎn)速區(qū)間內(nèi)故障相值與健康相存在明顯區(qū)分。

    圖13 不同轉(zhuǎn)速下ai變化趨勢(shì)

    從上述分析可以看出,當(dāng)電機(jī)發(fā)生匝間短路故障時(shí),故障相值會(huì)明顯降低,剩余健康相值基本不變,接近于1。剩余健康相值基本不隨短路匝數(shù)、剩余相電流幅值和電機(jī)轉(zhuǎn)速變化,極大地簡(jiǎn)化了診斷流程,提高了故障診斷魯棒性和準(zhǔn)確度。而故障相值隨短路匝數(shù)非單調(diào)變化,還受負(fù)載和轉(zhuǎn)速影響,很難基于值對(duì)故障程度做出判斷。

    有限元仿真得到轉(zhuǎn)速為1 000 r/min、相電流幅值為30 A情況下,A~E相繞組分別發(fā)生1匝匝間短路故障時(shí),探測(cè)線圈雷達(dá)圖,不同相短路時(shí)的探測(cè)線圈感應(yīng)電動(dòng)勢(shì)幅值如圖14所示。當(dāng)A~E某一相繞組發(fā)生匝間短路故障,該相容錯(cuò)齒上的探測(cè)線圈值降低,其他容錯(cuò)齒上的值基本不變,根據(jù)這一特點(diǎn)可以實(shí)現(xiàn)對(duì)故障位置的判斷。

    圖14 不同相短路時(shí)的探測(cè)線圈感應(yīng)電動(dòng)勢(shì)幅值

    綜上所述,本文所提出的基于容錯(cuò)齒磁通的故障診斷方法可以實(shí)現(xiàn)對(duì)匝間短路故障的識(shí)別和定位。需要注意的是,受電機(jī)加工誤差、轉(zhuǎn)速計(jì)算精度等因素的影響,實(shí)際系統(tǒng)中,非匝間短路相的值可能在1附近輕微波動(dòng),為了防止誤判,實(shí)際操作過(guò)程中,需要為設(shè)置一定的閾值。本文結(jié)合上述分析給出了基于容錯(cuò)齒磁通的故障診斷方法流程,如圖15所示,具體流程總結(jié)如下:

    圖15 故障診斷流程

    (1)采集TA+-TA-~TE+-TE-探測(cè)線圈端電壓。探測(cè)線圈外接電壓采集電路,采集電機(jī)運(yùn)行過(guò)程中的各線圈端電壓。

    (2)濾波得到TA、TB、TC、TD、TE。對(duì)采集到的端電壓進(jìn)行濾波,消除高頻干擾成分,計(jì)算得到TA、TB、TC、TD、TE。

    (3)根據(jù)當(dāng)前轉(zhuǎn)速計(jì)算得到T0。根據(jù)位置傳感器反饋計(jì)算當(dāng)前轉(zhuǎn)速,根據(jù)式(7)計(jì)算得到當(dāng)前轉(zhuǎn)速下的T0值。

    (4)計(jì)算得到A、B、C、D、E。以步驟(2)中得到的TA~TE值和步驟(3)中得到的T0值為輸入,根據(jù)式(8)計(jì)算得到A~E的值。

    (5)判斷A~E的值是否低于閾值。若值均在閾值以上,判定未發(fā)生匝間短路故障;若出現(xiàn)低于閾值的相,判定該相發(fā)生匝間短路故障。

    此外,本文所采用的探測(cè)線圈繞制在電機(jī)容錯(cuò)齒上,會(huì)占用一部分槽面積。相同狀態(tài)下,探測(cè)線圈線匝數(shù)越多,線圈感應(yīng)電動(dòng)勢(shì)幅值越高,占用的槽面積也越大,探測(cè)線圈匝數(shù)選取應(yīng)綜合考慮電機(jī)工況、后端電壓采集電路接口設(shè)計(jì)和電機(jī)功率密度。下面以本文所研究的電機(jī)為例,對(duì)探測(cè)線圈匝數(shù)選取及引入探測(cè)線圈對(duì)電機(jī)性能的影響進(jìn)行分析。

    對(duì)于圖1所示電機(jī),最高轉(zhuǎn)速為10 000 r/min,單匝探測(cè)線圈感應(yīng)電動(dòng)勢(shì)系數(shù)為ET=8.15× 10-5V·s/rad,可以得到單匝探測(cè)線圈感應(yīng)電動(dòng)勢(shì)最大值為5.12 V,當(dāng)后端測(cè)量電路輸入電壓范圍為0~30 V時(shí),探測(cè)線圈匝數(shù)選取為5匝。后端測(cè)量電路輸入電流通常在20 mA以內(nèi),按照電流密度5 A/mm2標(biāo)準(zhǔn)選取,探測(cè)線圈最小線徑為0.06 mm,考慮線圈強(qiáng)度,可選擇0.12 mm線徑漆包線,外徑為0.14 mm。按照電機(jī)槽滿率估算方法,5匝線圈占用面積約為0.1 mm2,考慮探測(cè)線圈外層包覆絕緣套管等,實(shí)際線圈占用面積在0.5 mm2以內(nèi)。該電機(jī)定子單個(gè)槽槽面積為130 mm2,探測(cè)線圈占槽面積比例不到0.4%,探測(cè)線圈對(duì)槽滿率和電機(jī)性能的影響基本可以忽略。

    3.3 定子結(jié)構(gòu)參數(shù)對(duì)故障診斷敏感性的影響分析

    從2.2節(jié)中的分析可以看出,匝間短路狀態(tài)下,s1影響流經(jīng)容錯(cuò)齒Ⅰ的槽漏磁通大小,d1影響流經(jīng)容錯(cuò)齒Ⅰ的主磁通大小,這兩個(gè)參數(shù)對(duì)探測(cè)線圈感應(yīng)電動(dòng)勢(shì)幅值影響較大,進(jìn)而對(duì)故障診斷產(chǎn)生影響。而s1和d1主要由槽口寬度01和容錯(cuò)齒長(zhǎng)度t決定,01和t的定義如圖16所示。01和t還會(huì)對(duì)電機(jī)定子磁導(dǎo)分布、漏磁等產(chǎn)生影響,進(jìn)而影響電機(jī)短路電流、轉(zhuǎn)矩等性能。下面采用有限元仿真研究01和t變化對(duì)、短路電流、平均轉(zhuǎn)矩和轉(zhuǎn)矩波動(dòng)的影響規(guī)律。

    圖16 槽口寬度b01和容錯(cuò)齒Ⅰ長(zhǎng)度Lt示意圖

    仿真得到轉(zhuǎn)速為1 000 r/min、相電流幅值為30 A情況下,、短路電流、平均轉(zhuǎn)矩和轉(zhuǎn)矩波動(dòng)隨槽口寬度01的變化,結(jié)果如圖17所示。隨著01的增大,低匝數(shù)匝間短路情況下的值逐漸增大,高匝數(shù)匝間短路情況下的值逐漸減小,選取適中的01值有助于兼顧高低匝數(shù)匝間短路故障診斷靈敏度。此外,隨著01的增大,短路電流幅值逐漸增大,平均轉(zhuǎn)矩先增大后減小,轉(zhuǎn)矩波動(dòng)逐漸降低。為了更好地兼顧高低匝數(shù)短路故障診斷靈敏度、電機(jī)出力能力和短路電流抑制能力,本文選取01= 1.5 mm。

    圖17 不同b01下ai、短路電流和轉(zhuǎn)矩變化趨勢(shì)

    仿真得到轉(zhuǎn)速為1 000 r/min、相電流幅值為30 A情況下,、短路電流、平均轉(zhuǎn)矩和轉(zhuǎn)矩波動(dòng)隨容錯(cuò)齒長(zhǎng)度t的變化,結(jié)果如圖18所示。隨著t的減小,低匝數(shù)匝間短路情況下的值變化不大,高匝數(shù)匝間短路情況下的值逐漸增大,會(huì)出現(xiàn)故障相值與非故障相接近的情況,無(wú)法識(shí)別故障。此外,隨著t的減小,短路電流幅值和平均轉(zhuǎn)矩略微增大,轉(zhuǎn)矩波動(dòng)變化不明顯。綜上所述,為了保證故障識(shí)別能力,本文選取t=19.5 mm。

    圖18 不同Lt下ai、短路電流和轉(zhuǎn)矩變化趨勢(shì)

    4 結(jié)論

    本文建立了五相永磁同步電機(jī)定子磁路模型,給出了匝間短路狀態(tài)下電機(jī)定子齒磁通的變化規(guī)律,研究表明,匝間短路狀態(tài)下,電機(jī)繞線齒和容錯(cuò)齒中的定子磁通均會(huì)發(fā)生變化,其中,繞線齒Ⅰ中的磁通“由無(wú)到有”,這一特征可以用于故障診斷。對(duì)基于探測(cè)線圈的電機(jī)故障診斷方法進(jìn)行了研究,給出了不同繞制方式下探測(cè)線圈感應(yīng)電動(dòng)勢(shì)的變化規(guī)律,提出了基于容錯(cuò)齒磁通的匝間短路故障診斷方法,給出了基于故障標(biāo)志的故障識(shí)別和定位方法。研究表明,當(dāng)電機(jī)發(fā)生匝間短路故障時(shí),故障相值明顯降低,剩余健康相值接近于1,剩余健康相值基本不隨短路匝數(shù)、剩余相電流幅值和電機(jī)轉(zhuǎn)速變化。該方法不依賴傅里葉變換、不受負(fù)載和電機(jī)轉(zhuǎn)速影響,具有準(zhǔn)確度高、魯棒性強(qiáng)的優(yōu)點(diǎn)。給出了槽口寬度、容錯(cuò)齒長(zhǎng)度對(duì)故障敏感度和電機(jī)性能的影響規(guī)律,結(jié)果表明,為了兼顧高低匝數(shù)匝間短路故障診斷靈敏度和電機(jī)性能,槽口寬度01取值不宜過(guò)大或過(guò)小,且不宜采用短齒設(shè)計(jì)。

    [1] 陶濤, 趙文祥, 程明, 等. 多相電機(jī)容錯(cuò)控制及其關(guān)鍵技術(shù)綜述[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2019, 39(2): 316-326, 629.

    Tao Tao, Zhao Wenxiang, Cheng Ming, et al. Review on fault-tolerant control of multi-phase machines and their key technologies[J]. Proceedings of the CSEE, 2019, 39(2): 316-326, 629.

    [2] Zhao Wenxiang, Xu Liang, Liu Guohai. Overview of permanent-magnet fault-tolerant machines: topology and design[J]. CES Transactions on Electrical Machines and Systems, 2018, 2(1): 51-64.

    [3] 李修東, 鄭曉欽, 王海峰, 等. 最大轉(zhuǎn)矩范圍內(nèi)九相永磁同步電機(jī)缺相容錯(cuò)運(yùn)行銅耗優(yōu)化策略[J]. 電工技術(shù)學(xué)報(bào), 2022, 37(17): 4355-4363.

    Li Xiudong, Zheng Xiaoqin, Wang Haifeng, et al. Copper loss optimization strategy for nine-phase permanent magnet synchronous motors fault-tolerant operation in maximum torque range[J]. Transactions of China Electrotechnical Society, 2022, 37(17): 4355-4363.

    [4] 陳前, 夏雨航, 趙文祥, 等. 采用無(wú)差拍電流跟蹤的五相梯形反電動(dòng)勢(shì)永磁電機(jī)開(kāi)路容錯(cuò)控制[J]. 電工技術(shù)學(xué)報(bào), 2022, 37(2): 368-379.

    Chen Qian, Xia Yuhang, Zhao Wenxiang, et al. Open- circuit fault-tolerant control for five-phase permanent magnet motors with trapezoidal back-EMF by deadbeat current tracking[J]. Transactions of China Electrotechnical Society, 2022, 37(2): 368-379.

    [5] 隋義. 純電動(dòng)汽車用五相容錯(cuò)永磁同步電機(jī)的關(guān)鍵技術(shù)研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2015.

    [6] 尹佐生. 航天伺服用四相容錯(cuò)永磁電機(jī)研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2018.

    [7] 張業(yè)成, 劉國(guó)海, 陳前. 基于電流波動(dòng)特征的永磁同步電機(jī)匝間短路與局部退磁故障分類診斷研究[J]. 電工技術(shù)學(xué)報(bào), 2022, 37(7): 1634-1643, 1653.

    Zhang Yecheng, Liu Guohai, Chen Qian. Discrimi- nation of interturn short-circuit and local demag- netization in permanent magnet synchronous motor based on current fluctuation characteristics[J]. Transactions of China Electrotechnical Society, 2022, 37(7): 1634-1643, 1653.

    [8] Hang Jun, Zhang Jianzhong, Cheng Ming, et al. Detection and discrimination of open-phase fault in permanent magnet synchronous motor drive system[J]. IEEE Transactions on Power Electronics, 2016, 31(7): 4697-4709.

    [9] Estima J O, Marques Cardoso A J. A new algorithm for real-time multiple open-circuit fault diagnosis in voltage-fed PWM motor drives by the reference current errors[J]. IEEE Transactions on Industrial Electronics, 2013, 60(8): 3496-3505.

    [10] 張瀟丹, 段斌, 劉昌杰, 等. 基于本體追蹤的永磁同步電機(jī)驅(qū)動(dòng)系統(tǒng)開(kāi)路故障診斷方法[J]. 電工技術(shù)學(xué)報(bào), 2022, 37(13): 3341-3349, 3360.

    Zhang Xiaodan, Duan Bin, Liu Changjie, et al. Open circuit fault diagnosis method for permanent magnet synchronous machine drive system based on ontology tracking[J]. Transactions of China Electrotechnical Society, 2022, 37(13): 3341-3349, 3360.

    [11] Duran M J, Gonzalez-Prieto I, Rios-Garcia N, et al. A simple, fast, and robust open-phase fault detection technique for six-phase induction motor drives[J]. IEEE Transactions on Power Electronics, 2018, 33(1): 547-557.

    [12] 謝穎, 胡圣明, 陳鵬, 等. 永磁同步電機(jī)匝間短路故障溫度場(chǎng)分析[J]. 電工技術(shù)學(xué)報(bào), 2022, 37(2): 322-331.

    Xie Ying, Hu Shengming, Chen Peng, et al. Thermal field analysis on inter-turn short circuit fault of permanent magnet synchronous motor[J]. Transa- ctions of China Electrotechnical Society, 2022, 37(2): 322-331.

    [13] 杜博超. 電動(dòng)汽車永磁同步電機(jī)驅(qū)動(dòng)系統(tǒng)故障診斷與容錯(cuò)技術(shù)研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2016.

    [14] Briz F, Degner M W, Zamarron A, et al. Online stator winding fault diagnosis in inverter-fed AC machines using high-frequency signal injection[J]. IEEE Transa- ctions on Industry Applications, 2003, 39(4): 1109- 1117.

    [15] Hu Rongguang, Wang Jiabin, Mills A R, et al. High-frequency voltage injection based stator inter- turn fault detection in permanent magnet machines[J]. IEEE Transactions on Power Electronics, 2021, 36(1): 785-794.

    [16] Ebrahimi B M, Faiz J. Feature extraction for short- circuit fault detection in permanent-magnet syn- chronous motors using stator-current monitoring[J]. IEEE Transactions on Power Electronics, 2010, 25(10): 2673-2682.

    [17] Kim K H. Simple online fault detecting scheme for short-circuited turn in a PMSM through current harmonic monitoring[J]. IEEE Transactions on Industrial Electronics, 2011, 58(6): 2565-2568.

    [18] Urresty J C, Riba J R, Romeral L. Diagnosis of inter- turn faults in PMSMs operating under nonstationary conditions by applying order tracking filtering[J]. IEEE Transactions on Power Electronics, 2013, 28(1): 507-515.

    [19] Hang Jun, Zhang Jianzhong, Cheng Ming, et al. Online interturn fault diagnosis of permanent magnet synchronous machine using zero-sequence com- ponents[J]. IEEE Transactions on Power Electronics, 2015, 30(12): 6731-6741.

    [20] 孫宇光, 杜威, 桂林, 等. 用于多相無(wú)刷勵(lì)磁機(jī)開(kāi)路與短路故障檢測(cè)的磁極探測(cè)線圈設(shè)計(jì)[J]. 電工技術(shù)學(xué)報(bào), 2022, 37(14): 3542-3554.

    Sun Yuguang, Du Wei, Gui Lin, et al. Design of pole detection coils for open-circuit and short-circuit faults in multiphase brushless exciter[J]. Transactions of China Electrotechnical Society, 2022, 37(14): 3542- 3554.

    [21] Chai Jie, Wang Jiabin, Atallah K, et al. Performance comparison and winding fault detection of duplex 2-phase and 3-phase fault-tolerant permanent magnet brushless machines[C]//2007 IEEE Industry Appli- cations Annual Meeting, New Orleans, LA, USA, 2007: 566-572.

    [22] Kim K T, Lee S T, Hur J. Diagnosis technique using a detection coil in BLDC motors with interturn faults[J]. IEEE Transactions on Magnetics, 2014, 50(2): 885- 888.

    [23] Lee S T, Kim K T, Hur J. Diagnosis technique for stator winding inter-turn fault in BLDC motor using detection coil[C]//2015 9th International Conference on Power Electronics and ECCE Asia (ICPE-ECCE Asia), Seoul, Korea (South), 2015: 2925-2931.

    [24] 曾沖. 基于定子齒磁通的永磁同步電機(jī)故障診斷研究[D]. 重慶: 重慶大學(xué), 2019.

    Diagnosis of Inter-Turn Short-Circuit Fault for Five-Phase Permanent Magnet Machine Based on Fault-Tolerant Tooth Flux

    11112

    (1. Beijing Institute of Automatic Control Equipment Beijing 100074 China 2. School of Electrical Engineering and Automation Harbin Institute of Technology Harbin 150001 China)

    Inter-turn short-circuit fault is a common machine winding fault with great harm and easy- diffusion characteristics. The fast inter-turn short-circuit fault diagnosis is necessary to avoid further fatal failure, which is a research focus for the multi-phase fault-tolerant permanent magnet synchronous machine (PMSM). The existing short-circuit fault diagnosis methods based on high-frequency signal injection, current harmonic component extraction, zero-sequence voltage monitoring, and detection coil are mostly bound with fast Fourier transform (FFT). The dependence on FFT brings problems of increased computation, longer diagnosis time, and sensitivity to operating conditions. It is hard to achieve the balance of diagnosis speed, accuracy, and simplicity. Therefore, an inter-turn short-circuit fault diagnosis method based on fault-tolerant tooth flux is proposed for a five-phase PMSM. The proposed method is independent of FFT, load, and machine speed, which has the advantages of high accuracy, fast speed, simple process, and robustness.

    Firstly, the stator magnetic circuit model of the five-phase PMSM is established, and the influence of inter-turn short-circuit fault on the flux of winding and fault-tolerant teeth is given. Secondly, variations of the induced electromotive force (EMF) of different arranged detection coils with the operating state of the machine are investigated. The fault identification ability, sensitivity and complexity of each winding scheme are comprehensively evaluated. Thirdly, the fault identification and locating method by fault signis presented for the fault diagnosis based on the fault-tolerant tooth flux. Finally, the influence law of slot-opening width and fault-tolerant tooth length on fault sensitivity and machine performances is further analyzed, and the size selection method is given. There is flux in the fault-tolerant teeth between in-phase coils only when the inter-turn short-circuit fault occurs in the phase. Compared with other arranging schemes, the fault detection coils arranged on the fault-tolerant tooth can distinguish the normal condition, open-circuit fault, and short-circuit fault explicitly, which is adopted for fault diagnosis. The decreased induced EMF amplitude in a specific detection coil indicates an inter-turn short-circuit fault. The fault signis introduced as the ratio of detection coil-induced EMF amplitude to the reference value to characterize the decreasing detection coil induced EMF in the fault phase. When the inter-turn short-circuit fault occurs to the machine,αvalue of the fault phase is significantly reduced,values of the remaining healthy phases are close to 1. According to this feature, faults can be accurately identified and located. With the increase of slot-opening width01, the sensitivities to small and large turn short-circuit faults are lower and higher, respectively. The fault diagnosis is more complex with the decrease of fault-tolerant tooth lengtht. With additional consideration of machine performances, values of01andtare chosen as 1.5 mm and 19.5 mm, respectively.

    The following conclusions can be drawn from the analytical and simulation analysis: (1) The flux in the fault-tolerant teeth between in-phase coils appears only when the inter-turn short-circuit fault occurs, which is a key characteristic for fault diagnosis. (2) An inter-turn short-circuit fault diagnosis method based on fault-tolerant tooth flux is proposed, with independent detection coils in each fault-tolerant tooth. (3) When the inter-turn short-circuit fault occurs, the fault signvalue of the fault phase is significantly reduced, andvalues of the remaining healthy phases are close to 1, which can indicate fault appearance and fault position. (4) The proposed method is independent of FFT, load, and machine speed, with the advantages of high accuracy, fast speed, simple process, and robustness. (5) To balance the sensitivities of small and large turn short-circuit fault diagnosis and machine performances, slot-opening width01should not be too large or too small, and a short tooth design is not recommended.

    Five-phase permanent magnet synchronous machine, inter-turn short-circuit, detection coil, fault-diagnosis

    國(guó)家自然科學(xué)基金(52175038)和國(guó)家自然科學(xué)基金青年基金(51907009)資助項(xiàng)目。

    2022-07-19

    2022-09-08

    10.19595/j.cnki.1000-6753.tces.221382

    TM351

    黃 建 男,1982年生,博士,研究員,研究方向?yàn)槎嘞嗳蒎e(cuò)永磁電機(jī)及其控制、電動(dòng)伺服控制技術(shù)。E-mail: huangjian421@126.com

    尹佐生 男,1994年生,博士,研究方向?yàn)槎嘞嗳蒎e(cuò)永磁電機(jī)及其控制、電動(dòng)伺服控制技術(shù)。E-mail: yinzuosheng0102@163.com(通信作者)

    (編輯 崔文靜)

    猜你喜歡
    繞線匝間匝數(shù)
    轉(zhuǎn)鼓過(guò)濾機(jī)繞線故障分析及防范措施
    一種無(wú)刷雙饋電機(jī)匝數(shù)配合的半解析優(yōu)化方法*
    高壓電機(jī)定子線圈匝間膠化工藝研究
    抽水蓄能發(fā)電電動(dòng)機(jī)轉(zhuǎn)子匝間短路的重復(fù)脈沖法診斷研究
    交流繞線電機(jī)電刷牌號(hào)及數(shù)量選擇
    不同側(cè)出線的變壓器線圈匝數(shù)計(jì)算
    次級(jí)線圈繞線錐度對(duì)LVDT靜態(tài)特性的影響
    馬達(dá)繞線參數(shù)對(duì)其功能的影響研究
    電子制作(2017年20期)2017-04-26 06:57:36
    基于EEMD-HHT的雙饋感應(yīng)發(fā)電機(jī)定子匝間短路故障研究
    大型水輪發(fā)電機(jī)中性點(diǎn)CT匝間短路故障分析處理
    欧美中文综合在线视频| 黄色片一级片一级黄色片| 久久久久国产一级毛片高清牌| 国产激情欧美一区二区| 老熟妇仑乱视频hdxx| 国产私拍福利视频在线观看| 69av精品久久久久久| 麻豆成人午夜福利视频| 成人一区二区视频在线观看| 久久婷婷成人综合色麻豆| 国产免费av片在线观看野外av| 日本a在线网址| 99国产精品99久久久久| 波多野结衣高清无吗| 国产精品久久久久久精品电影 | 伦理电影免费视频| 精品不卡国产一区二区三区| 精品国产美女av久久久久小说| 搡老妇女老女人老熟妇| 亚洲成国产人片在线观看| 99久久99久久久精品蜜桃| 每晚都被弄得嗷嗷叫到高潮| 99在线人妻在线中文字幕| a在线观看视频网站| 法律面前人人平等表现在哪些方面| 又黄又粗又硬又大视频| 巨乳人妻的诱惑在线观看| АⅤ资源中文在线天堂| 制服诱惑二区| 亚洲精品av麻豆狂野| 精品国产超薄肉色丝袜足j| 日韩欧美一区二区三区在线观看| 99riav亚洲国产免费| 99热这里只有精品一区 | 亚洲一区高清亚洲精品| 日本熟妇午夜| 日本 av在线| 搡老妇女老女人老熟妇| √禁漫天堂资源中文www| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品国产综合久久久| 一区二区三区高清视频在线| 91成年电影在线观看| 美女 人体艺术 gogo| 99精品在免费线老司机午夜| 精品人妻1区二区| 久久久国产精品麻豆| 久99久视频精品免费| 看免费av毛片| 精品人妻1区二区| a级毛片a级免费在线| 女警被强在线播放| 久久久久免费精品人妻一区二区 | 麻豆国产av国片精品| 男女那种视频在线观看| 午夜视频精品福利| 久久久久久久午夜电影| av片东京热男人的天堂| 又黄又爽又免费观看的视频| 亚洲精品久久成人aⅴ小说| 狠狠狠狠99中文字幕| 久久精品国产99精品国产亚洲性色| 国产精品亚洲美女久久久| 欧美日本亚洲视频在线播放| 18禁美女被吸乳视频| 国产精品亚洲av一区麻豆| 国产精品亚洲美女久久久| 日韩成人在线观看一区二区三区| 亚洲国产欧美一区二区综合| 人人澡人人妻人| 亚洲精品在线美女| 国产精品日韩av在线免费观看| 欧美日韩福利视频一区二区| 亚洲一区二区三区不卡视频| 久久亚洲精品不卡| 亚洲精品久久国产高清桃花| 免费在线观看日本一区| 人人妻,人人澡人人爽秒播| 亚洲狠狠婷婷综合久久图片| 校园春色视频在线观看| 夜夜看夜夜爽夜夜摸| 亚洲av美国av| 可以在线观看的亚洲视频| 国产激情偷乱视频一区二区| 91成人精品电影| 91av网站免费观看| 午夜成年电影在线免费观看| 国产一级毛片七仙女欲春2 | 欧美在线一区亚洲| 午夜成年电影在线免费观看| 悠悠久久av| 亚洲精品粉嫩美女一区| 在线av久久热| 精品少妇一区二区三区视频日本电影| 久久精品影院6| 免费在线观看黄色视频的| 精品久久久久久久人妻蜜臀av| 久久欧美精品欧美久久欧美| 亚洲一码二码三码区别大吗| svipshipincom国产片| 侵犯人妻中文字幕一二三四区| 久久精品国产清高在天天线| 黄色片一级片一级黄色片| 在线观看舔阴道视频| 国产在线观看jvid| a级毛片在线看网站| 在线播放国产精品三级| 丝袜在线中文字幕| 少妇裸体淫交视频免费看高清 | 九色国产91popny在线| 一级黄色大片毛片| 国产一区二区激情短视频| 国产亚洲av嫩草精品影院| 99热这里只有精品一区 | 伊人久久大香线蕉亚洲五| 99在线人妻在线中文字幕| 99riav亚洲国产免费| 久久久久久久久免费视频了| 一区二区三区激情视频| 久久婷婷人人爽人人干人人爱| 禁无遮挡网站| 国产99久久九九免费精品| 国产精品影院久久| 国产精品日韩av在线免费观看| 欧美日韩一级在线毛片| 亚洲午夜理论影院| 丁香欧美五月| 久久精品成人免费网站| 黄色a级毛片大全视频| 老司机深夜福利视频在线观看| 婷婷精品国产亚洲av在线| 满18在线观看网站| 国产激情久久老熟女| 桃红色精品国产亚洲av| 97超级碰碰碰精品色视频在线观看| 99国产精品99久久久久| 91成人精品电影| 日韩av在线大香蕉| 狂野欧美激情性xxxx| 国产黄片美女视频| 波多野结衣高清作品| 久久中文字幕一级| svipshipincom国产片| 亚洲欧美精品综合一区二区三区| 久久精品国产清高在天天线| 免费在线观看亚洲国产| 国产精品99久久99久久久不卡| 日韩欧美三级三区| 国产成人欧美在线观看| 精品少妇一区二区三区视频日本电影| 国产精品自产拍在线观看55亚洲| 老司机午夜福利在线观看视频| 免费无遮挡裸体视频| 在线观看一区二区三区| 人成视频在线观看免费观看| 免费在线观看日本一区| 久久人人精品亚洲av| 啪啪无遮挡十八禁网站| 国产成人av激情在线播放| 亚洲av成人不卡在线观看播放网| 欧美性猛交╳xxx乱大交人| av片东京热男人的天堂| 国产精品一区二区精品视频观看| 日本精品一区二区三区蜜桃| 国产真实乱freesex| 久久久久免费精品人妻一区二区 | 国产野战对白在线观看| 国产高清有码在线观看视频 | 熟女少妇亚洲综合色aaa.| 两性午夜刺激爽爽歪歪视频在线观看 | 国产av在哪里看| 亚洲性夜色夜夜综合| 国产成人一区二区三区免费视频网站| 少妇 在线观看| 精品国产美女av久久久久小说| 久久香蕉激情| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩 欧美 亚洲 中文字幕| 极品教师在线免费播放| 国产亚洲欧美精品永久| 久久精品91蜜桃| 国产欧美日韩精品亚洲av| 国产色视频综合| 9191精品国产免费久久| 在线观看一区二区三区| 别揉我奶头~嗯~啊~动态视频| 中文字幕av电影在线播放| 18禁观看日本| 啪啪无遮挡十八禁网站| 91老司机精品| 一a级毛片在线观看| 视频在线观看一区二区三区| 亚洲 欧美一区二区三区| 国产在线观看jvid| 亚洲三区欧美一区| 黄频高清免费视频| 欧美黑人精品巨大| 午夜福利在线观看吧| 999精品在线视频| 亚洲国产中文字幕在线视频| 亚洲精品中文字幕一二三四区| 亚洲午夜理论影院| 久久精品夜夜夜夜夜久久蜜豆 | 午夜a级毛片| 亚洲国产精品sss在线观看| 国产精品免费视频内射| 久久热在线av| 这个男人来自地球电影免费观看| 黄片小视频在线播放| 国产不卡一卡二| 十八禁人妻一区二区| 99热这里只有精品一区 | 国产在线观看jvid| 又大又爽又粗| www.熟女人妻精品国产| 色婷婷久久久亚洲欧美| 99国产极品粉嫩在线观看| 国产又色又爽无遮挡免费看| 久久九九热精品免费| √禁漫天堂资源中文www| 超碰成人久久| 亚洲中文av在线| 中文字幕人妻熟女乱码| 久久国产亚洲av麻豆专区| 国产精品一区二区精品视频观看| 久久精品亚洲精品国产色婷小说| 最近最新免费中文字幕在线| 国产精品精品国产色婷婷| 国产av一区二区精品久久| 色婷婷久久久亚洲欧美| 在线av久久热| 人妻久久中文字幕网| 国产精品国产高清国产av| 在线免费观看的www视频| 给我免费播放毛片高清在线观看| 欧美最黄视频在线播放免费| 久久久久亚洲av毛片大全| 亚洲狠狠婷婷综合久久图片| 在线看三级毛片| 亚洲av五月六月丁香网| 好男人电影高清在线观看| 久久久国产精品麻豆| 热99re8久久精品国产| 精品高清国产在线一区| 丁香欧美五月| 中亚洲国语对白在线视频| 女警被强在线播放| 免费在线观看视频国产中文字幕亚洲| 无遮挡黄片免费观看| 午夜亚洲福利在线播放| 法律面前人人平等表现在哪些方面| 国产黄片美女视频| 国产欧美日韩精品亚洲av| 最近在线观看免费完整版| 国产1区2区3区精品| 大型黄色视频在线免费观看| 露出奶头的视频| 非洲黑人性xxxx精品又粗又长| 精品国产国语对白av| 99国产极品粉嫩在线观看| 欧美日韩亚洲综合一区二区三区_| 亚洲国产欧美一区二区综合| 国产单亲对白刺激| 看片在线看免费视频| 精品电影一区二区在线| 国产黄色小视频在线观看| 日韩一卡2卡3卡4卡2021年| 欧美国产精品va在线观看不卡| 亚洲自拍偷在线| 满18在线观看网站| 精品久久久久久久久久久久久 | 日韩欧美国产在线观看| 麻豆久久精品国产亚洲av| 欧美激情极品国产一区二区三区| 欧美日韩一级在线毛片| 国产成人精品久久二区二区免费| 国内少妇人妻偷人精品xxx网站 | 欧美激情久久久久久爽电影| 国产av一区二区精品久久| 一边摸一边抽搐一进一小说| 亚洲欧洲精品一区二区精品久久久| 岛国在线观看网站| 欧美激情 高清一区二区三区| 国产色视频综合| 怎么达到女性高潮| 精品人妻1区二区| 男女之事视频高清在线观看| 午夜激情福利司机影院| 麻豆久久精品国产亚洲av| 国产一区二区三区在线臀色熟女| 婷婷丁香在线五月| 一边摸一边做爽爽视频免费| 亚洲无线在线观看| 午夜久久久在线观看| 中文亚洲av片在线观看爽| 亚洲国产欧洲综合997久久, | 美女国产高潮福利片在线看| 嫩草影院精品99| 色综合欧美亚洲国产小说| 国语自产精品视频在线第100页| 国产精品久久久久久精品电影 | 欧美+亚洲+日韩+国产| 欧美性猛交黑人性爽| 国产单亲对白刺激| 欧美三级亚洲精品| 亚洲九九香蕉| 亚洲男人的天堂狠狠| 日韩精品青青久久久久久| 午夜亚洲福利在线播放| 精品一区二区三区av网在线观看| 村上凉子中文字幕在线| 国内精品久久久久久久电影| 亚洲va日本ⅴa欧美va伊人久久| 美女高潮到喷水免费观看| x7x7x7水蜜桃| 99久久久亚洲精品蜜臀av| 极品教师在线免费播放| svipshipincom国产片| 午夜影院日韩av| 亚洲国产欧美一区二区综合| 日日摸夜夜添夜夜添小说| 亚洲成人国产一区在线观看| 99久久精品国产亚洲精品| 老熟妇乱子伦视频在线观看| a级毛片a级免费在线| 国产激情久久老熟女| 亚洲自拍偷在线| 大香蕉久久成人网| 手机成人av网站| 啦啦啦 在线观看视频| 亚洲成a人片在线一区二区| 精品电影一区二区在线| 亚洲精品久久国产高清桃花| 亚洲av五月六月丁香网| or卡值多少钱| 午夜成年电影在线免费观看| 欧美黑人精品巨大| 久久香蕉激情| 久久国产精品影院| 亚洲性夜色夜夜综合| 99国产精品99久久久久| 国产精品乱码一区二三区的特点| 日韩视频一区二区在线观看| 人人澡人人妻人| 自线自在国产av| 欧美色视频一区免费| 欧美性长视频在线观看| 叶爱在线成人免费视频播放| 91老司机精品| 亚洲熟妇中文字幕五十中出| 色av中文字幕| 妹子高潮喷水视频| 一本一本综合久久| 亚洲午夜精品一区,二区,三区| 日本a在线网址| 黄色a级毛片大全视频| 国产爱豆传媒在线观看 | 麻豆av在线久日| 成人欧美大片| АⅤ资源中文在线天堂| 欧美另类亚洲清纯唯美| 男女下面进入的视频免费午夜 | 国产精品 欧美亚洲| 亚洲国产欧美日韩在线播放| 国产精品九九99| 丁香六月欧美| 少妇被粗大的猛进出69影院| e午夜精品久久久久久久| 美国免费a级毛片| 亚洲狠狠婷婷综合久久图片| 国产成人av激情在线播放| 久久人人精品亚洲av| 亚洲五月婷婷丁香| 国产免费男女视频| 国产高清有码在线观看视频 | 成人三级做爰电影| 亚洲中文字幕一区二区三区有码在线看 | 欧美成人性av电影在线观看| 亚洲av中文字字幕乱码综合 | 亚洲欧美激情综合另类| 午夜福利高清视频| 国产麻豆成人av免费视频| 好男人电影高清在线观看| 日本一区二区免费在线视频| 久久欧美精品欧美久久欧美| 99久久精品国产亚洲精品| 国产av在哪里看| 亚洲,欧美精品.| 欧美中文日本在线观看视频| 在线观看免费日韩欧美大片| 国产极品粉嫩免费观看在线| 在线观看免费视频日本深夜| 露出奶头的视频| 高清在线国产一区| 老汉色∧v一级毛片| av有码第一页| 老司机午夜福利在线观看视频| 午夜激情av网站| 亚洲狠狠婷婷综合久久图片| 伦理电影免费视频| 免费看a级黄色片| 琪琪午夜伦伦电影理论片6080| 亚洲人成伊人成综合网2020| 精品国产亚洲在线| 久久天堂一区二区三区四区| 欧美在线一区亚洲| 国产精品野战在线观看| 久久国产精品男人的天堂亚洲| 日日干狠狠操夜夜爽| 美女扒开内裤让男人捅视频| 欧美国产精品va在线观看不卡| 1024香蕉在线观看| 国产97色在线日韩免费| 中文在线观看免费www的网站 | 欧美人与性动交α欧美精品济南到| 我的亚洲天堂| 亚洲精品在线美女| 91字幕亚洲| 日韩高清综合在线| 最近最新中文字幕大全电影3 | 一进一出抽搐动态| 久久久久免费精品人妻一区二区 | 久久人妻av系列| 亚洲va日本ⅴa欧美va伊人久久| 久久婷婷人人爽人人干人人爱| 麻豆国产av国片精品| 日日爽夜夜爽网站| 亚洲 国产 在线| 国产精品久久久久久人妻精品电影| 日本一本二区三区精品| 男人舔奶头视频| 成年人黄色毛片网站| 无限看片的www在线观看| 亚洲成人国产一区在线观看| 狂野欧美激情性xxxx| 亚洲熟女毛片儿| 老司机午夜福利在线观看视频| 午夜福利高清视频| 国产色视频综合| 国产精品乱码一区二三区的特点| 他把我摸到了高潮在线观看| 亚洲片人在线观看| 青草久久国产| 香蕉久久夜色| 亚洲成av片中文字幕在线观看| 精品久久蜜臀av无| 成人国产综合亚洲| 欧美性猛交╳xxx乱大交人| 国产在线精品亚洲第一网站| 国产主播在线观看一区二区| 国产成人一区二区三区免费视频网站| 一级a爱视频在线免费观看| 国产黄色小视频在线观看| 可以在线观看毛片的网站| 亚洲精品色激情综合| 国产亚洲欧美98| 妹子高潮喷水视频| 男女视频在线观看网站免费 | av视频在线观看入口| 日韩欧美在线二视频| 老司机在亚洲福利影院| 国产精品自产拍在线观看55亚洲| 18禁黄网站禁片午夜丰满| 夜夜躁狠狠躁天天躁| 午夜福利在线在线| tocl精华| 搞女人的毛片| 搡老岳熟女国产| 夜夜夜夜夜久久久久| 日本成人三级电影网站| 嫁个100分男人电影在线观看| 久久婷婷成人综合色麻豆| 国产熟女午夜一区二区三区| 人人妻人人澡欧美一区二区| 国产97色在线日韩免费| 国产又色又爽无遮挡免费看| 国产成人啪精品午夜网站| 中文字幕人妻丝袜一区二区| 精品国产乱子伦一区二区三区| 成人亚洲精品一区在线观看| 成人永久免费在线观看视频| 欧美黑人巨大hd| 少妇 在线观看| 在线观看免费午夜福利视频| 久久久久亚洲av毛片大全| 欧美日韩亚洲国产一区二区在线观看| 一区二区日韩欧美中文字幕| 免费观看精品视频网站| 国产精品,欧美在线| 人人澡人人妻人| 久久久久久久久久黄片| 精品卡一卡二卡四卡免费| 国产一区二区激情短视频| 欧美国产精品va在线观看不卡| 国产91精品成人一区二区三区| 国内久久婷婷六月综合欲色啪| 波多野结衣av一区二区av| 99热6这里只有精品| 12—13女人毛片做爰片一| 亚洲免费av在线视频| 无遮挡黄片免费观看| 国产人伦9x9x在线观看| 啦啦啦 在线观看视频| 色播在线永久视频| 欧美黑人巨大hd| 亚洲成人免费电影在线观看| 国产视频内射| 搡老岳熟女国产| 黄片播放在线免费| av在线播放免费不卡| 人人澡人人妻人| 88av欧美| 国产激情欧美一区二区| 国产伦人伦偷精品视频| 手机成人av网站| 别揉我奶头~嗯~啊~动态视频| ponron亚洲| 免费看日本二区| 久久狼人影院| 国产精品野战在线观看| 欧美性猛交╳xxx乱大交人| 在线观看66精品国产| 18美女黄网站色大片免费观看| 国产爱豆传媒在线观看 | 美女 人体艺术 gogo| 欧美在线黄色| 高清在线国产一区| 日韩视频一区二区在线观看| 一级黄色大片毛片| 成年人黄色毛片网站| 天天躁狠狠躁夜夜躁狠狠躁| 熟女少妇亚洲综合色aaa.| 一二三四社区在线视频社区8| 久久久久久久久中文| 最近最新免费中文字幕在线| 日本黄色视频三级网站网址| 亚洲成国产人片在线观看| 亚洲国产精品999在线| 成人18禁高潮啪啪吃奶动态图| 久久久水蜜桃国产精品网| 色播在线永久视频| 这个男人来自地球电影免费观看| 一进一出抽搐gif免费好疼| 精品欧美一区二区三区在线| 国产蜜桃级精品一区二区三区| 黄色毛片三级朝国网站| 亚洲九九香蕉| 久久久久久免费高清国产稀缺| 一级作爱视频免费观看| 欧美黑人精品巨大| 国内精品久久久久精免费| 亚洲一码二码三码区别大吗| 每晚都被弄得嗷嗷叫到高潮| 两人在一起打扑克的视频| 91麻豆精品激情在线观看国产| 99精品在免费线老司机午夜| 老鸭窝网址在线观看| 国产一区二区三区视频了| 特大巨黑吊av在线直播 | 精品第一国产精品| 欧美黑人欧美精品刺激| 免费在线观看视频国产中文字幕亚洲| 免费看十八禁软件| 伊人久久大香线蕉亚洲五| 一卡2卡三卡四卡精品乱码亚洲| 免费在线观看完整版高清| 国产免费男女视频| 日韩高清综合在线| 欧美av亚洲av综合av国产av| 精品一区二区三区四区五区乱码| 久久中文看片网| 50天的宝宝边吃奶边哭怎么回事| 性色av乱码一区二区三区2| 欧美亚洲日本最大视频资源| 日韩一卡2卡3卡4卡2021年| 免费在线观看完整版高清| 在线国产一区二区在线| 啪啪无遮挡十八禁网站| 久久精品91蜜桃| 午夜激情福利司机影院| 神马国产精品三级电影在线观看 | 欧美大码av| 国产野战对白在线观看| 99久久国产精品久久久| 欧美av亚洲av综合av国产av| 国产成人欧美在线观看| 国产三级在线视频| 精品日产1卡2卡| 亚洲成人精品中文字幕电影| 午夜激情福利司机影院| 波多野结衣av一区二区av| 欧美丝袜亚洲另类 | 国产久久久一区二区三区| 国产爱豆传媒在线观看 | 91成年电影在线观看| 欧美国产日韩亚洲一区| 美女午夜性视频免费| 黄色丝袜av网址大全| 日本 av在线| 国产aⅴ精品一区二区三区波| 日韩精品免费视频一区二区三区| 亚洲成国产人片在线观看| 国产男靠女视频免费网站| x7x7x7水蜜桃| 国产不卡一卡二| 丰满人妻熟妇乱又伦精品不卡| 国产成年人精品一区二区| 久久久久久久久免费视频了| 欧美黑人欧美精品刺激| 国产高清videossex| 午夜a级毛片| 国产99久久九九免费精品| 久久香蕉国产精品| 欧美国产精品va在线观看不卡| 人成视频在线观看免费观看|