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

    弓網(wǎng)電弧等離子體光譜特性實(shí)驗(yàn)

    2017-01-21 06:28:16怡魏文賦雷棟高國強(qiáng)吳廣寧
    電工技術(shù)學(xué)報(bào) 2016年24期
    關(guān)鍵詞:弓網(wǎng)譜線滑板

    胡 怡魏文賦雷 棟高國強(qiáng)吳廣寧

    (1. 西南交通大學(xué)電氣工程學(xué)院 成都 610031 2. 中國鐵道科學(xué)研究院標(biāo)準(zhǔn)計(jì)量研究所 北京 100081)

    弓網(wǎng)電弧等離子體光譜特性實(shí)驗(yàn)

    胡 怡1魏文賦1雷 棟2高國強(qiáng)1吳廣寧1

    (1. 西南交通大學(xué)電氣工程學(xué)院 成都 610031 2. 中國鐵道科學(xué)研究院標(biāo)準(zhǔn)計(jì)量研究所 北京 100081)

    弓網(wǎng)電弧等離子體具有高溫度、高能量的特點(diǎn),對高速鐵路弓網(wǎng)電接觸的性能造成了威脅?;诠W(wǎng)電弧模擬實(shí)驗(yàn)平臺(tái),利用光譜診斷法對弓網(wǎng)電弧等離子體進(jìn)行研究。對輻射的主要特征譜線進(jìn)行標(biāo)識(shí)和歸屬,并計(jì)算弓網(wǎng)電弧等離子體的激發(fā)溫度、轉(zhuǎn)動(dòng)溫度、振動(dòng)溫度和電子密度。此外,還研究了電源輸出電流對激發(fā)溫度和電子密度的影響。實(shí)驗(yàn)結(jié)果表明,弓網(wǎng)電弧等離子體對接觸網(wǎng)銅導(dǎo)線和浸金屬碳滑板強(qiáng)烈的燒蝕作用會(huì)產(chǎn)生豐富的銅原子譜線、鐵原子譜線和CN B2∑+→X2∑+譜線。同時(shí),基于Boltzmann斜線法,發(fā)現(xiàn)了弓網(wǎng)電弧等離子體激發(fā)溫度隨電流的增加而增加。通過擬合電弧等離子體中CN B2∑+→X2∑+譜線,可知在30A電流條件下,弓網(wǎng)電弧等離子體的轉(zhuǎn)動(dòng)溫度和振動(dòng)溫度分別達(dá)到6 800K和9 000K。最后,討論了特征譜線的展寬機(jī)制,并利用Cu I 521.82nm譜線,探討了弓網(wǎng)電弧等離子體的電子密度隨電流的增加而上升的情況。

    弓網(wǎng)電弧等離子體 發(fā)射光譜 激發(fā)溫度 轉(zhuǎn)動(dòng)溫度 振動(dòng)溫度

    0 引言

    弓網(wǎng)關(guān)系作為高速鐵路關(guān)系安全運(yùn)行的三大基礎(chǔ)關(guān)系之一,保障著動(dòng)車組的穩(wěn)定受流[1,2]。隨著動(dòng)車組運(yùn)行速度的提升,受電弓與接觸網(wǎng)之間的振動(dòng)加劇,導(dǎo)致弓網(wǎng)離線頻繁,弓網(wǎng)電弧頻發(fā)。電弧等離子體有高能量、高溫度的特點(diǎn),不僅會(huì)使電力系統(tǒng)或設(shè)備承受高頻振蕩過電壓[3,4],還會(huì)嚴(yán)重?zé)g受電弓和接觸網(wǎng),影響列車的行車安全[5,6]。由于目前在高速列車的運(yùn)營現(xiàn)場,并不具備對弓網(wǎng)電弧等離子體進(jìn)行實(shí)驗(yàn)的條件,主要側(cè)重于外部電氣參數(shù)的測量,并不能直觀地反映出燃弧的強(qiáng)度。對弓網(wǎng)電弧等離子體的燃燒機(jī)理進(jìn)行研究,有利于對其進(jìn)行抑制,其中對溫度的測量就是上述研究的關(guān)鍵技術(shù)之一。

    弓網(wǎng)電弧等離子體處于高速列車運(yùn)行產(chǎn)生的強(qiáng)氣流場中,內(nèi)部包含大量的分子、原子和離子。光譜是粒子狀態(tài)發(fā)生躍遷時(shí),所釋放出波長不同的電磁輻射,按電磁輻射的本質(zhì)將光譜區(qū)分為:原子光譜、分子光譜、X射線能譜、γ 射線能譜等種類。目前國內(nèi)外相關(guān)研究人員針對弓網(wǎng)電弧的實(shí)驗(yàn)開展了大量的工作,但主要都集中于對其電氣特性的探索。S. Midya等[7,8]研究了高速鐵路情況下弓網(wǎng)電弧的現(xiàn)象,主要分析了速度和電流等因素對弓網(wǎng)電弧的影響和作用機(jī)制。L. Morin等[9]研究了不同負(fù)載對電弧能量的影響機(jī)制,結(jié)果表明感性負(fù)載回路中電弧能量明顯高于阻性負(fù)載回路中的電弧能量。吳廣寧等[5,10]基于自行研制的弓網(wǎng)電弧實(shí)驗(yàn)系統(tǒng)研究了電弧對受電弓滑板和接觸網(wǎng)導(dǎo)線的侵蝕現(xiàn)象,并對弓網(wǎng)電弧的抑制措施進(jìn)行了初步探索。

    然而對于弓網(wǎng)電弧等離子體特性的實(shí)驗(yàn)研究,目前還較少見諸報(bào)端。近年來,隨著光譜儀等光學(xué)設(shè)備的發(fā)展,使用光學(xué)的方法診斷等離子體參數(shù)變得更加簡便、有效、精確。目前針對大氣壓等離子體參數(shù)診斷方式常用Langmuir探針法和光學(xué)診斷法,但探針法會(huì)對等離子體本身造成干擾,影響診斷結(jié)果準(zhǔn)確性。斯紅等[11]利用發(fā)射光譜法,分析了TIG焊接電弧的計(jì)算譜線選擇規(guī)則和其溫度場分布。張喬根等[12]對染污絕緣子表面的局部電弧,基于光譜診斷手段,分析了其內(nèi)部溫度和電子密度隨泄漏電流的變化情況。對于等離子體輻射的分子光譜,彭志敏等[13]基于火焰發(fā)射光譜,利用光學(xué)多通道分析儀,對OH自由基和CH自由基的實(shí)驗(yàn)光譜和理論計(jì)算光譜進(jìn)行對比,分析了轉(zhuǎn)動(dòng)溫度和振動(dòng)溫度。屠昕等[14]基于N2+的B2∑+→X2∑+譜帶對電弧等離子體的振動(dòng)溫度和轉(zhuǎn)動(dòng)溫度進(jìn)行了測量,并分析了轉(zhuǎn)動(dòng)溫度和振動(dòng)溫度對譜帶輪廓的影響機(jī)制。周亦驍?shù)萚15]測量了大氣壓下Ar/O2和Ar/H2O等離子體射流的發(fā)射光譜,并分析了轉(zhuǎn)動(dòng)溫度和振動(dòng)溫度隨外加電壓幅值的變化情況。

    本文基于弓網(wǎng)電弧實(shí)驗(yàn)系統(tǒng),模擬列車運(yùn)行時(shí)的弓網(wǎng)振動(dòng)情況,產(chǎn)生電弧并利用光譜儀采集其發(fā)射光譜,選取合適的譜線,基于Boltzmann斜線法,計(jì)算了在不同電流下,弓網(wǎng)電弧等離子體的激發(fā)溫度。同時(shí),通過實(shí)驗(yàn)光譜與理論計(jì)算光譜進(jìn)行擬合,得到了等離子體的轉(zhuǎn)動(dòng)溫度和振動(dòng)溫度。討論了弓網(wǎng)電弧等離子體的特征譜線展寬機(jī)制,并利用譜線的展寬,分析了在不同電流下弓網(wǎng)電弧等離子體的電子密度。

    1 弓網(wǎng)電弧實(shí)驗(yàn)系統(tǒng)

    課題組研制了一套弓網(wǎng)電弧模擬系統(tǒng)[5],為保證模擬的準(zhǔn)確性,受電弓滑板和接觸線材料與實(shí)際弓網(wǎng)系統(tǒng)一致,分別為銅合金接觸網(wǎng)導(dǎo)線和浸金屬碳滑板,滑板的主要組分見表1。弓網(wǎng)電弧光譜實(shí)驗(yàn)系統(tǒng)裝置如圖1所示。

    表1 浸金屬碳滑板組分含量Tab.1 Component content of pantograph strips

    圖1中箭頭方向分別表示接觸線導(dǎo)線CW和受電弓滑板P的運(yùn)動(dòng)方向。變頻電機(jī)驅(qū)動(dòng)轉(zhuǎn)盤W帶動(dòng)接觸導(dǎo)線。受電弓滑板的運(yùn)動(dòng)受工控機(jī)IPC控制,由伺服電動(dòng)機(jī)進(jìn)行操作,可縱向模擬弓網(wǎng)系統(tǒng)的振動(dòng),整個(gè)弓網(wǎng)系統(tǒng)由直流電源S供電,其電壓可調(diào)范圍為0~80V,電流范圍為0~60A。同時(shí)利用示波器采集電弧電壓、電流波形,示波器型號為Tektronix MDO3024,帶寬200MHz。

    圖1 弓網(wǎng)電弧光譜實(shí)驗(yàn)裝置S—電源R—列車負(fù)載電阻L—列車負(fù)載電感 W—轉(zhuǎn)盤CW—接觸網(wǎng)導(dǎo)線 P—受電弓滑板 PC Plasma—弓網(wǎng)電弧等離子體SP—光譜儀 F—光纖 IPC—工控機(jī)Fig.1 The pantograph-catenary arc test device

    弓網(wǎng)電弧是一種典型氣體放電現(xiàn)象,其電壓主要由陰極電壓降、陽極電壓降和弧柱電壓共同構(gòu)成,其中陰極和陽極電壓降與電弧所處氣體和電極材料有關(guān),在空氣中分別取值約為8~11V和3~12V[16],而弧柱電壓與電弧長度有關(guān),其單位約為15V/mm[17]。所以當(dāng)弓網(wǎng)間隙為4mm時(shí),弓網(wǎng)電弧的電壓為80V左右。在牽引供電系統(tǒng)中,27.5kV的電壓主要消耗在車體、牽引變壓器等大阻抗器件,實(shí)際的小間隙弓網(wǎng)電弧電壓值也較小,與本文所采用的實(shí)驗(yàn)裝置仍具一定等效性。

    當(dāng)實(shí)驗(yàn)系統(tǒng)產(chǎn)生弓網(wǎng)電弧時(shí),通過光譜儀對弓網(wǎng)電弧等離子體的發(fā)射光譜進(jìn)行采集,將采集的光譜存儲(chǔ)至工控機(jī),以便處理。實(shí)驗(yàn)所采用的是PG-2000光譜儀,其具體參數(shù)見表2。

    表2 PG-2000光譜儀參數(shù)Tab.2 Parameters of PG-2000 spectroscope

    本實(shí)驗(yàn)?zāi)M電力機(jī)車在高速運(yùn)行時(shí),由于弓網(wǎng)系統(tǒng)振動(dòng),產(chǎn)生電?。河晒た貦C(jī)控制伺服電動(dòng)機(jī),驅(qū)動(dòng)受電弓滑板在接觸線下方振動(dòng),滑板的最高位置是與接觸線恰好接觸上。振動(dòng)頻率為1Hz,振動(dòng)間距固定設(shè)置為4mm,光譜儀的積分時(shí)間為100μs。

    在弓網(wǎng)電弧實(shí)驗(yàn)平臺(tái),進(jìn)行受電弓滑板和接觸線從接觸到分離的一次完整拉弧過程,采集到弓網(wǎng)電弧典型的電壓和電流波形如圖2所示。在拉弧過程中,在圖中Q點(diǎn)處,滑板和接觸線分離,電弧產(chǎn)生,隨著弓網(wǎng)間隙增大,電壓逐漸增大,電流始終保持恒定,在1s后弓網(wǎng)電弧熄滅,此時(shí)弓網(wǎng)間隙的電壓等于電源電壓,電流為零。

    圖2 弓網(wǎng)電弧電壓、電流波形Fig.2 The voltage and current waveforms of pantograph-catenary arc

    2 弓網(wǎng)電弧輻射譜

    2.1 弓網(wǎng)電弧等離子體光譜特性分析

    經(jīng)光譜儀校正曲線校正后,得到直流電流30A的情況下,弓網(wǎng)電弧等離子體波長范圍為350~500nm和500~800nm的光譜特征譜線分別如圖3a和圖3b所示。對特征譜線進(jìn)行比對和分析,其粒子類型和波長均如圖3所示。

    從圖3中可以看到弓網(wǎng)電弧等離子體輻射光譜存在以下信息:

    (1)弓網(wǎng)電弧等離子體輻射光譜主要由Cu原子譜線構(gòu)成,這是因弓網(wǎng)電弧等離子體對銅接觸線和受電弓碳滑板強(qiáng)烈的燒蝕作用導(dǎo)致。除Cu原子譜線外,還觀測到了其他金屬元素(如Fe)的特征譜線,這是由于受電弓碳滑板和銅鎳導(dǎo)線中均含有微量的鐵。除此之外,還有Cu+離子特征譜線,但譜線強(qiáng)度較低,其原因是Cu+離子的激發(fā)能遠(yuǎn)大于激發(fā)Cu原子的能量。

    圖3 弓網(wǎng)電弧等離子體輻射光譜Fig.3 The emission spectrum of pantographcatenary arc plasma

    (2)此外,還觀測到了CN自由基B2∑+→X2∑+電子帶系譜帶,其產(chǎn)生主要是受電弓碳滑板在弓網(wǎng)電弧等離子體燒蝕下會(huì)產(chǎn)生C2分子,與空氣中的氮?dú)獍l(fā)生反應(yīng)C2+N2=2CN。

    (3)在弓網(wǎng)電弧等離子體輻射光譜中還有O I 777.53nm的特征譜線,O原子的激發(fā)也是由于高能電子的碰撞,可能還會(huì)影響弓網(wǎng)電弧等離子體的燃燒強(qiáng)度。

    2.2 Boltzmann法測激發(fā)溫度

    弓網(wǎng)電弧等離子體是一種由中性粒子和帶電粒子共同組成,整體呈電中性的氣體。根據(jù)統(tǒng)計(jì)動(dòng)力學(xué)觀點(diǎn),溫度是一個(gè)統(tǒng)計(jì)的概念,在等離子體中,某種粒子各自由度下的平均能量對應(yīng)著不同的溫度,對于分子的平動(dòng)、振動(dòng)、轉(zhuǎn)動(dòng)自由度分別有平動(dòng)溫度、振動(dòng)溫度和轉(zhuǎn)動(dòng)溫度。電子溫度是指等離子體中自由電子的平動(dòng)溫度,而激發(fā)溫度對應(yīng)的是原子或離子(重粒子)內(nèi)部束縛態(tài)電子的激發(fā)自由度。只有當(dāng)系統(tǒng)達(dá)到完全熱力學(xué)平衡時(shí),即各自由度的溫度相同時(shí),系統(tǒng)才具有一個(gè)統(tǒng)一的熱力學(xué)溫度[18]。

    等離子體的激發(fā)溫度是影響輻射譜線強(qiáng)度的重要因素[19]。不考慮共振光電子激發(fā)過程,輻射譜線強(qiáng)度與等離子體溫度的關(guān)系為

    式中,I為譜線強(qiáng)度;T為等離子體溫度;n為等離子體的中性粒子密度;Aul為躍遷幾率;h為普朗克常數(shù),h=6.63×10-34J·s;v為譜線的頻率;gu為上能級統(tǒng)計(jì)權(quán)重;Z為原子配分函數(shù);Eu為高能態(tài)激發(fā)能;k為玻耳茲曼常數(shù),k=1.38×10-23J/K。

    對式(1)進(jìn)行變換后,左右兩端取對數(shù)得

    式中,λ為譜線波長。

    當(dāng)式(2)中的激發(fā)能單位為eV(1eV=11 600K)時(shí);玻耳茲曼常數(shù)k經(jīng)換算后為8.618×10-5eV/K;從式(2)可知與Eu是線性關(guān)系,對于同種原子的多條譜線,以和Eu分別作為橫、縱坐標(biāo)描出各點(diǎn),利用最小二乘法進(jìn)行擬合,擬合直線的斜率為-1/(kT)。

    由圖3可知,弓網(wǎng)電弧等離子體中含有大量的Cu原子譜線,通過NIST查詢到Cu原子特征譜線的參數(shù),見表3。

    表3 Cu原子特征譜線參數(shù)Tab.3 Parameters of copper atomic characteristic spectral lines

    由于在利用譜線強(qiáng)度計(jì)算電子溫度時(shí),需要使用躍遷幾率,但理論計(jì)算躍遷幾率值往往存在較大的不確定度,影響弓網(wǎng)電弧等離子體激發(fā)溫度計(jì)算的準(zhǔn)確性,為減小誤差,通常會(huì)選擇多譜線進(jìn)行擬合[20]。除此之外,對于計(jì)算溫度的譜線選擇而言,還需考慮其強(qiáng)度情況,強(qiáng)度越低的譜線也越容易受背景噪聲干擾。I. L. Babich等[21]針對Cu原子譜線的躍遷幾率的準(zhǔn)確性情況進(jìn)行了研究,發(fā)現(xiàn)Cu I 529.25nm譜線強(qiáng)度可能受N I發(fā)射譜線強(qiáng)度的干擾,會(huì)較為偏離Boltzmann圖的擬合結(jié)果,所以在弓網(wǎng)電弧等離子體激發(fā)溫度的計(jì)算中不采用。此外,Cu I 793.31nm在NIST中未能查詢到躍遷幾率,也同樣不能用于計(jì)算。

    綜上所述,針對弓網(wǎng)電弧等離子電子激發(fā)溫度的計(jì)算,選擇Cu I 465.11nm、Cu I 510.55nm、Cu I 515.32nm、Cu I 521.82nm、Cu I 570.02nm和Cu I 578.21nm六條譜線,其擬合結(jié)果如圖4所示。

    圖4 弓網(wǎng)電弧等離子體Boltzmann圖Fig.4 The Boltzmann plot of pantographcatenary arc plasma

    通過線性擬合,最后得到在電流為30A,如圖4的實(shí)驗(yàn)條件情況下,弓網(wǎng)電弧等離子體的Boltzmann圖的斜率為-1.508 7×10-4,計(jì)算得到激發(fā)溫度為4 163K。其不確定度可根據(jù)式(3)[22]得出。

    在弓網(wǎng)系統(tǒng)離線過程中,弓網(wǎng)電弧等離子體是弓網(wǎng)系統(tǒng)唯一的導(dǎo)電通道。接觸網(wǎng)、弓網(wǎng)系統(tǒng)和電力機(jī)車構(gòu)成了牽引供電回路,電弧中的電流會(huì)發(fā)生變化。由于弓網(wǎng)電弧燃弧不穩(wěn)定,為保證實(shí)驗(yàn)結(jié)果的準(zhǔn)確性,在同一電流下,共進(jìn)行了20次重復(fù)實(shí)驗(yàn)。利用弓網(wǎng)電弧實(shí)驗(yàn)平臺(tái),分析了在不同電流情況下,弓網(wǎng)電弧等離子體激發(fā)溫度實(shí)驗(yàn)結(jié)果分布情況如圖5所示。

    從圖5中可以看出,弓網(wǎng)電弧等離子體激發(fā)溫度隨電流的增大有上升趨勢,這是因?yàn)殡娀≈械碾娏髟龃?,注入弓網(wǎng)電弧等離子體的能量也隨之增大,導(dǎo)致激發(fā)溫度上升。弓網(wǎng)電弧等離子體燃弧隨機(jī)性較強(qiáng),隨著電流增大,激發(fā)溫度范圍也變大。當(dāng)電流值從30A增至60A,弓網(wǎng)電弧等離子體的激發(fā)溫度平均值由4 258K上升到4 763K。

    圖5 弓網(wǎng)電弧等離子體激發(fā)溫度與電流的關(guān)系Fig.5 Relationship of excitation temperature and current of pantograph-catenary arc plasma

    2.3 電弧等離子體的振動(dòng)溫度和轉(zhuǎn)動(dòng)溫度

    在不同電流情況下的弓網(wǎng)電弧等離子體光譜中,均檢測到了CN分子自由基B2∑+→X2∑+電子帶系,CN分子屬于雙原子分子,其分子結(jié)構(gòu)簡單,轉(zhuǎn)動(dòng)常數(shù)大,較易識(shí)別,特別是在高溫(6 000K以上)情況,相對其他分子自由基更為穩(wěn)定[23]。

    對于CN分子而言,由于內(nèi)部不同電子能級的躍遷產(chǎn)生譜帶群,弓網(wǎng)電弧等離子體中的CN的B2∑+→X2∑+躍遷形成了CN violet光譜群,而兩個(gè)電子能級B2∑+和X2∑+中的振動(dòng)能級間發(fā)生躍遷,形成譜帶(v′,v′′),而每一個(gè)譜帶間振動(dòng)能級的所有轉(zhuǎn)動(dòng)能級發(fā)生躍遷,形成譜線[24]。而譜線的波長極為靠近,受光譜儀分辨率的限制,只能呈現(xiàn)具有一定輪廓的譜線。對于弓網(wǎng)電弧等離子體而言,通過對其輻射光譜的輪廓進(jìn)行擬合,可以分析其振動(dòng)溫度和轉(zhuǎn)動(dòng)溫度。

    首先假定弓網(wǎng)電弧等離子體的轉(zhuǎn)動(dòng)溫度和振動(dòng)溫度,通過LIFBASE仿真軟件計(jì)算該轉(zhuǎn)動(dòng)和振動(dòng)溫度下CN violet分子的輻射光譜的譜線形狀,與實(shí)驗(yàn)測量的譜帶形狀進(jìn)行對比,通過Peak correlation和R-square兩個(gè)參數(shù)判斷擬合程度,最終確定CN分子的振動(dòng)溫度和轉(zhuǎn)動(dòng)溫度。

    對CN自由基B2∑+→X2∑+電子帶系譜帶在波長390~425nm處進(jìn)行擬合,該譜帶區(qū)域包含了CN violet帶系的Δν=0和Δν=-1兩個(gè)振動(dòng)能級變化帶,圖6為利用LIFBASE進(jìn)行譜帶擬合的CN自由基譜帶和實(shí)驗(yàn)測量的弓網(wǎng)電弧等離子體譜帶比對圖,經(jīng)LIFBASE擬合的結(jié)果,擬合參數(shù)最優(yōu)情況為Peak correlation=0.98和R-square=59,此時(shí)弓網(wǎng)電弧等離子體的振動(dòng)溫度和轉(zhuǎn)動(dòng)溫度分別為9 000K和6 800K。

    由圖6可以看出,CN B2∑+→X2∑+自由基仿真的光譜結(jié)果在兩個(gè)振動(dòng)帶都與實(shí)驗(yàn)的譜線符合較好。由于CN分子的轉(zhuǎn)動(dòng)能級間的能級差較小,其布居主要決定于重粒子的非彈性碰撞,轉(zhuǎn)動(dòng)能級的激發(fā)只需要較短時(shí)間,即所對應(yīng)的弛豫時(shí)間也較小。所以通常在大氣壓等離子體中,認(rèn)為轉(zhuǎn)動(dòng)溫度約等于氣體溫度[25,26]。

    對于弓網(wǎng)電弧等離子體而言,轉(zhuǎn)動(dòng)溫度和激發(fā)溫度不相等,原因可能是,在該實(shí)驗(yàn)條件下,弓網(wǎng)電弧等離子體偏離局部熱平衡狀態(tài),Boltzmann法計(jì)算得到的激發(fā)溫度僅能反映等離子體中原子內(nèi)部束縛態(tài)電子能級的布居情況,并不能反映等離子體中自由電子的溫度[27]。

    圖6 CN自由基B2∑+→X2∑+光譜擬合圖Fig.6 Fitting spectrum of CN violet(B2∑+→X2∑+)

    2.4 弓網(wǎng)電弧等離子體的電子密度

    本文利用了譜線展寬效應(yīng),計(jì)算弓網(wǎng)電弧等離子體中的電子密度。輻射原子和離子在周圍電子和離子擾動(dòng)下會(huì)引起譜線展寬,通過測量譜線展寬計(jì)算電子密度是一種簡單且可定量分析的手段,目前在電弧等離子體、激光誘導(dǎo)等離子體等都得到了大量應(yīng)用,并適用于當(dāng)局部熱平衡不占優(yōu)勢的等離子體[28]。

    弓網(wǎng)電弧等離子體特征譜線的線型主要是由幾種不同的加寬機(jī)制共同作用導(dǎo)致的,包括自然展寬、儀器展寬、多普勒展寬和碰撞展寬[29]。

    其中多普勒展寬來源于原子或離子的熱運(yùn)動(dòng),對于處于局部熱平衡情況的等離子體,認(rèn)為原子和離子的分布滿足Maxwell速率分布,該展寬機(jī)制作用下的譜線線型為Gauss線型。儀器展寬由光譜儀的狹縫寬度和散射情況決定,可通過汞燈等標(biāo)準(zhǔn)光源進(jìn)行標(biāo)定。而碰撞展寬主要是因電子或離子對等離子體中的原子或離子進(jìn)行碰撞,所引起的的展寬效應(yīng),稱為Stark展寬,該展寬機(jī)制作用下的譜線線型為Lorenz線型[30]。

    在Doppler展寬和Stark展寬的共同作用下,等離子體的譜線線型可近似由高斯線型和洛倫茲線型共同組成的Voigt函數(shù)擬合[31]。

    利用Origin 9.0中的Voigt函數(shù)對Cu I 521.82nm處的譜線展寬進(jìn)行擬合,擬合結(jié)果如圖7所示,擬合參數(shù)表征量R2=0.979,說明譜線形狀與Voigt線型擬合良好。

    圖7 Cu I 521.82nm擬合曲線Fig.7 The fitting curve of copper atomic line at 521.82nm

    譜線實(shí)際Stark半展寬為

    由于多普勒展寬反映的是原子和離子的熱運(yùn)動(dòng),所以與離子(重粒子)的溫度有關(guān),計(jì)算弓網(wǎng)電弧等離子體Cu I 521.82nm處的多普勒半寬,其計(jì)算公式為

    式中,Ti為重粒子溫度(即Boltzmann斜線法所對應(yīng)溫度);λ0為譜線中心波長;c為光速;m為原子質(zhì)量。經(jīng)過計(jì)算得到的多普勒展寬僅為0.002 4nm,并不是弓網(wǎng)電弧等離子體中譜線展寬的主要因素,可忽略其影響。

    輻射譜線的Stark半展寬和電子密度的關(guān)系為

    式中,第一項(xiàng)主要對應(yīng)的是電子的碰撞效應(yīng),第二項(xiàng)是考慮了準(zhǔn)靜態(tài)的離子電場作用。在實(shí)際計(jì)算中,通常都是忽略第二項(xiàng)的作用,則式(6)簡化為

    式中,w為電子對特征譜線的展寬系數(shù),與溫度的關(guān)聯(lián)性極小[32],利用文獻(xiàn)[33]中Cu I 521.82nm的展寬系數(shù),對于擬合曲線的FWHM扣除儀器展寬,最后計(jì)算得到弓網(wǎng)電弧等離子體電子密度為1.02× 1017cm-3左右。

    在不同電流下,弓網(wǎng)電弧等離子體的電子密度的變化情況如圖8所示。

    圖8 在不同電流下弓網(wǎng)電弧等離子體電子密度Fig.8 Electron density of pantograph-catenary arc plasma at different current

    從圖8中可以看出,隨著弓網(wǎng)電弧等離子體電流的增加,注入電弧的能量增加,內(nèi)部原子和分子粒子發(fā)生電離增多,導(dǎo)致電子密度增加,直流電流由30A變化至60A時(shí),電子密度的平均值由1.02×1017cm-3上升至1.52×1017cm-3。

    3 結(jié)論

    1)弓網(wǎng)電弧等離子體的光譜中含有豐富的特征譜線,主要包括銅原子、銅離子、氧原子、鐵原子和CN B2∑+→X2∑+自由基分子,譜線和譜帶并存,所輻射譜線表明弓網(wǎng)電弧等離子體對浸金屬碳滑板和接觸網(wǎng)銅鎳導(dǎo)線均存在著強(qiáng)烈的燒蝕作用。

    2)通過Boltzmann法分析了弓網(wǎng)電弧等離子體的激發(fā)溫度,并分析了電流對激發(fā)溫度的影響。當(dāng)電流由30A增至60A,激發(fā)溫度從4 258K上升到4 763K。利用弓網(wǎng)電弧等離子體輻射光譜中CN violet分子進(jìn)行擬合,最終得到在電流30A的情況下,弓網(wǎng)電弧等離子體的振動(dòng)溫度和轉(zhuǎn)動(dòng)溫度分別為9 000K和6 800K。轉(zhuǎn)動(dòng)溫度和激發(fā)溫度不等表明弓網(wǎng)電弧等離子體在該實(shí)驗(yàn)條件下偏離局部熱平衡狀態(tài)。

    3)弓網(wǎng)電弧等離子體中的譜線展寬的作用機(jī)制主要是Stark展寬,隨著弓網(wǎng)電弧等離子體的電流從30A變?yōu)?0A,電子密度由1.02×1017cm-3上升至1.52×1017cm-3。

    [1] 吳廣寧, 古圳, 高國強(qiáng), 等. 弓網(wǎng)電弧形態(tài)特性試驗(yàn)研究[J]. 高電壓技術(shù), 2015, 41(11): 3532-3537. Wu Guangning, Gu Zhen, Gao Guoqiang, et al. Experimental research on morphological characteristics of pantograph-catenary arc[J]. High Voltage Engineering, 2015, 41(11): 3532-3537.

    [2] 郭鳳儀, 王喜利, 王智勇, 等. 弓網(wǎng)電弧輻射電場噪聲實(shí)驗(yàn)研究[J]. 電工技術(shù)學(xué)報(bào), 2015, 30(14): 220-225. Guo Fengyi, Wang Xili, Wang Zhiyong, et al. Research on radiated electric field noise of pantograph arc[J]. Transactions of China Electrotechnical Society, 2015, 30(14): 220-225.

    [3] 楊柳林, 李德奎, 陳延明, 等. 基于RTDS的電弧接地故障自定義建模及仿真分析[J]. 電力系統(tǒng)保護(hù)與控制, 2016, 44(16): 138-142. Yang Liulin, Li Dekui, Chen Yanming, et al. User defined modeling and simulation of arc grounding fault based on RTDS[J]. Power System Protection and Control, 2016, 44(16): 138-142.

    [4] 陳博博, 屈衛(wèi)鋒, 楊宏宇, 等. 小電流接地系統(tǒng)單相接地綜合電弧模型與選線方法的研究[J]. 電力系統(tǒng)保護(hù)與控制, 2016, 44(16): 1-7. Chen Bobo, Qu Weifeng, Yang Hongyu, et al. Research on single phase grounding arc model and line selection for neutral ineffectively grounding system[J]. Power System Protection and Control, 2016, 44(16): 1-7.

    [5] 王萬崗, 吳廣寧, 高國強(qiáng), 等. 高速鐵路弓網(wǎng)電弧試驗(yàn)系統(tǒng)[J]. 鐵道學(xué)報(bào), 2012, 34(4): 22-27. Wang Wangang, Wu Guangning, Gao Guoqiang, et al. The pantograph-catenary arc test system for highspeed railways[J]. Journal of the China Railway Society, 2012, 34(4): 22-27.

    [6] 吳積欽, 錢清泉. 弓網(wǎng)系統(tǒng)電弧侵蝕接觸線時(shí)的熱分析[J]. 鐵道學(xué)報(bào), 2008, 30(3): 31-34. Wu Jiqin, Qian Qingquan. Thermal analysis of arc erosion of contact wire of the pantograph & catenary system[J]. Journal of the China Railway Society, 2008, 30(3): 31-34.

    [7] Surajit M, Dierk B, Thorsten S, et al. Pantograph arcing in electrified railways—mechanism and influence of various parameters—part II: with AC traction power supply[J]. IEEE Transactions on Power Delivery, 2009, 24(4): 1940-1950.

    [8] Surajit M, Dierk B, Anders L, et al. Understanding pantograph arcing in electrified railways-influence of various parameters[C]//IEEE International Symposium on Electromagnetic Compatibility, 2008: 1-6.

    [9] Morin L, Jemaa N B, Jeannot D, et al. Contacts materials performances under break arc in automotive applications[J]. IEEE Transactions on Components and Packaging Technologies, 2000, 23(2): 367-375.

    [10] 王萬崗. 高速鐵路弓網(wǎng)電弧動(dòng)態(tài)特性研究[D]. 成都: 西南交通大學(xué), 2013.

    [11] 斯紅, 華學(xué)明, 張旺, 等. 基于Boltzmann光譜法的焊接電弧溫度場測量計(jì)算[J]. 光譜學(xué)與光譜分析, 2012, 32(9): 2311-2313. Si Hong, Hua Xueming, Zhang Wang, et al. Welding arc temperature field measurements based on Boltzmann spectrometry[J]. Spectroscopy and Spectral Analysis, 2012, 32(9): 2311-2313.

    [12] 楊昊, 張喬根, 龐磊, 等. 染污絕緣表面交流電弧發(fā)射光譜及等離子體特性研究[J]. 中國電機(jī)工程學(xué)報(bào), 2015, 35(18): 4808-4816. Yang Hao, Zhang Qiaogen, Pang Lei, et al. Study on AC arc discharge characteristics over the polluted insulation surface using optical emission spectroscopy[J]. Proceedings of the CSEE, 2015, 35(18): 4808-4816.

    [13] 彭志敏, 丁艷軍, 楊乾鎖, 等. 基于OH自由基A2∑+→X2∏r電子帶系發(fā)射光譜的溫度測量技術(shù)[J]. 物理學(xué)報(bào), 2011, 60(5): 250-259. Peng Zhimin, Ding Yanjun, Yang Qiansuo, et al. Emission spectra of OH radical (A2∑+→X2∏r) and its application on high temperature gas[J]. Acta Physica Sinica, 2011, 60(5): 250-259.

    [14] 屠昕, 嚴(yán)建華, 馬增益, 等. 基于N2+(B2∑+u→X2∑+g)的電弧等離子體振動(dòng)溫度和轉(zhuǎn)動(dòng)溫度測量[J]. 光譜學(xué)與光譜分析, 2006, 26(12): 2161-2165. Tu Xin, Yan Jianhua, Ma Zengyi, et al. Measurement of rotational and vibrational temperatures in arc plasma based on the first negative system of N2+(B2∑+u→X2∑+g)[J]. Spectroscopy and Spectral Analysis, 2006, 26(12): 2161-2165.

    [15] 周亦驍, 方志, 邵濤. Ar/O2和Ar/H2O中大氣壓等離子體射流放電特性的比較[J]. 電工技術(shù)學(xué)報(bào), 2014, 29(11): 229-238. Zhou Yixiao, Fang Zhi, Shao Tao. Comparison of discharge characteristics of atmospheric pressure plasma jet in Ar/O2and Ar/H2O mixtures[J]. Transactions of China Electrotechnical Society, 2014, 29(11): 229-238.

    [16] 王其平. 電器電弧理論[M]. 北京: 機(jī)械工業(yè)出版社, 1982.

    [17] Kapetanovic M, Knol P, Van der Sluis L. Analysis of cooling power and network interaction of a switchingarc in an SF6high-voltage circuit-breaker[C]//The Colloquium of CIGRE SC 13, Sarajevo, Yugoslavia, 1989.

    [18] 榮命哲, 劉定新, 李美, 等. 非平衡態(tài)等離子體的仿真研究現(xiàn)狀與新進(jìn)展[J]. 電工技術(shù)學(xué)報(bào), 2014, 29(6): 271-282. Rong Mingzhe, Liu Dingxin, Li Mei, et al. Research status and new progress on the numerical simulation of non-equilibrium plasmas[J]. Transactions of China Electrotechnical Society, 2014, 29(6): 271-282.

    [19] 董華軍, 廖敏夫, 鄒積巖, 等. 基于CCD真空開關(guān)電弧等離子體參數(shù)診斷方法[J]. 電工技術(shù)學(xué)報(bào), 2007, 22(6): 65-68. Dong Huajun, Liao Minfu, Zou Jiyan, et al. Methods of diagnosing the plasma parameters in vacuum switching arcs based on CCD[J]. Transactions of China Electrotechnical Society, 2007, 22(6): 65-68.

    [20] 吳蓉, 李燕, 朱順官, 等. 等離子體電子溫度的發(fā)射光譜法診斷[J]. 光譜學(xué)與光譜分析, 2008, 28(4): 731-735. Wu Rong, Li Yan, Zhu Shunguan, et al. Emission spectroscopy diagnostics of plasma electron yemperature[J]. Spectroscopy and Spectral Analysis, 2008, 28(4): 731-735.

    [21] Babich I L, Boretskij V F, Veklich A N, et al. Spectroscopic data and Stark broadening of Cu I and Ag I spectral lines: selection and analysis[J]. Advances in Space Research, 2014, 54(7): 1254-1263.

    [22] Magesh Thiyagarajan, John Scharer. Experimental investigation of ultraviolet laser induced plasma density and temper- ature evolution in air[J]. Journalof Applied Physics, 2008, 104(1): 013303 (1-12).

    [23] Kurosawa K, Sugita S, Fujita K, et al. Rotationaltemperature measurements of chemically reacting CN using band-tail spectra[J]. Journal of Thermophysics & Heat Transfer, 2009, 23(3): 463-472.

    [24] 梁卓, 羅海云, 王新新, 等. 氣流對氮?dú)饨橘|(zhì)阻擋放電氣體溫度及放電模式的影響[J]. 物理學(xué)報(bào), 2010, 59(12): 8739-8746. Liang Zhuo, Luo Haiyun, Wang Xinxin, et al. Influences of gas flow on gas temperature and discharge mode in dielectric barrier discharge of nitrogen at atmospheric pressure[J]. Acta Physica Sinica, 2010, 59(12): 8739-8746.

    [25] Laux C O, Spence T G, Kruger C H, et al. Optical diagnostics of atmospheric pressure air plasmas[J]. Plasma Sources Science & Technology, 2003, 12(2): 125-138.

    [26] Staack D, Farouk B, Gutsol A, et al. Characterization of a DC atmospheric pressure normal glow discharge[J]. Plasma Sources Science & Technology, 2005, 14(4): 700-711.

    [27] 趙文華, 唐皇哉, 沈巖, 等. 譜線強(qiáng)度法所測得溫度的物理意義[J]. 光譜學(xué)與光譜分析, 2007, 27(11): 2145-2149. Zhao Wenhua, Tang Huangzai, Shen Yan, et al. Physical meaning of temperature measured by spectral line intensity method[J]. Spectroscopy and Spectral Analysis, 2007, 27(11): 2145-2149.

    [28] Iwasaki M, Inui H, Matsudaira Y, et al. Nonequilibrium atmospheric pressure plasma with ultrahigh electron density and high performance for glass surface cleaning[J]. Applied Physics Letters, 2008, 92(8): 081503.

    [29] 唐曉閂, 李春燕, 朱光來, 等. 激光誘導(dǎo)Al等離子體中電子密度和溫度的實(shí)驗(yàn)研究[J]. 中國激光, 2004, 31(6): 687-692. Tang Xiaoshuan, Li Chunyan, Zhu Guanglai, et al. Experimental investigation on the electron density and electron temperature of laser-induced Al plasmas[J]. Chinese Journal of Lasers, 2004, 31(6): 687-692.

    [30] Hofmann S, Van Gessel A F H, Verreycken T, et al. Power dissipation, gas temperatures and electron densities of cold atmospheric pressure helium and argon RF plasma jets[J]. Plasma Sources Science and Technology, 2011, 20(6): 065010(1-12).

    [31] Rautian S G, Sobel'man I I. The effect of collisions on the doppler broadening of spectral lines[J]. Soviet Physics Uspekhi, 1967, 9(5): 701-716.

    [32] 潘成剛, 華學(xué)明, 張旺, 等. 傅里葉變換法計(jì)算焊接電弧光譜Stark展寬研究[J]. 光譜學(xué)與光譜分析, 2012, 32(7): 1739-1743. Pan Chenggang, Hua Xueming, Zhang Wang, et al. Calculating the stark broadening of welding arc spectra by Fourier transform method[J]. Spectroscopy and Spectral Analysis, 2012, 32(7): 1739-1743.

    [33] Konjevi? N, Wiese W L. Experimental Stark widths and shifts for spectral lines of neutral and ionized atoms[J]. Journal of Physical and Chemical Reference Data, 1990, 19(6): 1307-1385.

    Experimental Investigation on Spectral Characteristics of Pantograph-Catenary Arc Plasma

    Hu Yi1Wei Wenfu1Lei Dong2Gao Guoqiang1Wu Guangning1
    (1. School of Electrical Engineering Southwest Jiaotong University Chengdu 610031 China 2. Standards & Metrology Research Institute China Academy of Railway Sciences Beijing 100081 China)

    Pantograph-catenary arc plasma, with high temperature and high energy density, can severely threaten the electrical contact performance of pantograph-catenary system in high speed railway. In this paper, pantograph-catenary arc experiments were conducted based on arc simulative platform in laboratory. The characteristic of pantograph-catenary arc plasma was studied by spectral diagnosis method, and then the primary characteristic lines were recognized and classified. The excitation temperature, rotational temperature, vibrational temperature and electron density were also calculated, respectively. Meanwhile, the influence of input current on excitation temperature and electron density was analyzed. The results show that abundant copper, iron and CN B2∑+→X2∑+characteristic lines would occur due to the intense erosion effects of pantograph-catenary arc plasma on copper wire and metal impregnated strips. In addition, it is found that the excitation temperature of pantograph-catenary arc plasma increases with the increase of current based on Boltzmann plot. The rotational temperature and vibrational temperature reach up to 6 800K and 9 000K respectively when the current is 30A, according to the theoretical calculation with the spectral line of CN B2∑+→X2∑+. Furthermore, the broadening mechanism of characteristic lines was discussed, and by Cu I 521.82nmcharacteristic line the increasing tendency of electron density with the input current was analyzed.

    Pantograph-catenary arc plasma, emission spectrum, excitation temperature, rotational temperature, vibration temperature

    TM89

    胡 怡 女,1993年生,碩士,研究方向?yàn)楦咚勹F路弓網(wǎng)電弧。

    E-mail: yhuswjtu@foxmail.com

    魏文賦 男,1987年生,博士,講師,研究方向?yàn)榉烹姷入x子體。

    E-mail: wfwei@home.swjtu.edu.cn(通信作者)

    國家自然科學(xué)基金資助項(xiàng)目(U1234202、51325704、51577158、51607147)。

    2016-05-30 改稿日期 2016-08-14

    猜你喜歡
    弓網(wǎng)譜線滑板
    基于HITRAN光譜數(shù)據(jù)庫的合并譜線測溫仿真研究
    高速列車弓網(wǎng)故障響應(yīng)研究
    鐵合金光譜譜線分離實(shí)驗(yàn)研究
    電子測試(2018年11期)2018-06-26 05:56:00
    快樂滑板
    湖上的滑板
    滑板改造記
    鍶原子光鐘鐘躍遷譜線探測中的程序控制
    去玩滑板
    高速鐵路弓網(wǎng)電弧仿真研究進(jìn)展
    電氣化鐵道(2016年4期)2016-04-16 05:59:44
    弓網(wǎng)系統(tǒng)接觸電阻特性的研究
    国产精品嫩草影院av在线观看 | 久久人人精品亚洲av| 久久精品综合一区二区三区| 最近在线观看免费完整版| 国产精品1区2区在线观看.| 乱系列少妇在线播放| av国产免费在线观看| 色综合婷婷激情| 99久久成人亚洲精品观看| 国产成人av教育| 我要搜黄色片| АⅤ资源中文在线天堂| 自拍偷自拍亚洲精品老妇| 91午夜精品亚洲一区二区三区 | 日韩人妻高清精品专区| 国产精品98久久久久久宅男小说| 欧美日本视频| 国产成人a区在线观看| 校园春色视频在线观看| 国产亚洲欧美98| 久久精品国产亚洲网站| 国产麻豆成人av免费视频| 最近在线观看免费完整版| 1024手机看黄色片| 国产精品爽爽va在线观看网站| 搡老妇女老女人老熟妇| 国产精品自产拍在线观看55亚洲| a级毛片a级免费在线| 十八禁网站免费在线| 精品一区二区免费观看| 午夜a级毛片| 精品乱码久久久久久99久播| 午夜a级毛片| 国产高清激情床上av| 国产男靠女视频免费网站| 特级一级黄色大片| 观看免费一级毛片| 高清在线国产一区| 如何舔出高潮| 男女视频在线观看网站免费| 精品久久久噜噜| 97超级碰碰碰精品色视频在线观看| 欧美日韩乱码在线| 99久久久亚洲精品蜜臀av| 久久99热6这里只有精品| 国产精品美女特级片免费视频播放器| 男插女下体视频免费在线播放| 18禁裸乳无遮挡免费网站照片| 美女cb高潮喷水在线观看| 1024手机看黄色片| aaaaa片日本免费| x7x7x7水蜜桃| 人妻少妇偷人精品九色| 久久久久国内视频| 69人妻影院| 午夜亚洲福利在线播放| 听说在线观看完整版免费高清| a级毛片a级免费在线| 亚洲性夜色夜夜综合| 狂野欧美激情性xxxx在线观看| 九九热线精品视视频播放| 床上黄色一级片| 国内少妇人妻偷人精品xxx网站| 又粗又爽又猛毛片免费看| 日本色播在线视频| 亚洲精品成人久久久久久| eeuss影院久久| 午夜福利欧美成人| 日本 av在线| 色在线成人网| 久久人妻av系列| 麻豆国产av国片精品| 高清毛片免费观看视频网站| 亚洲最大成人av| 哪里可以看免费的av片| 欧美日韩黄片免| 亚洲一区高清亚洲精品| 一个人观看的视频www高清免费观看| 国产主播在线观看一区二区| 欧美日韩国产亚洲二区| 免费人成在线观看视频色| 美女被艹到高潮喷水动态| 午夜视频国产福利| 内射极品少妇av片p| 中文资源天堂在线| 中出人妻视频一区二区| 狠狠狠狠99中文字幕| 国产精品,欧美在线| 99国产极品粉嫩在线观看| 在线a可以看的网站| av在线亚洲专区| av黄色大香蕉| 久久久久精品国产欧美久久久| 国产色婷婷99| 国产精品久久久久久精品电影| 他把我摸到了高潮在线观看| 亚洲狠狠婷婷综合久久图片| 一区二区三区免费毛片| 三级毛片av免费| 99久国产av精品| 大又大粗又爽又黄少妇毛片口| 国产白丝娇喘喷水9色精品| 国产高潮美女av| 最新中文字幕久久久久| 97超视频在线观看视频| 久久久午夜欧美精品| 亚洲欧美日韩高清在线视频| 美女高潮的动态| 久久亚洲真实| av国产免费在线观看| 好男人在线观看高清免费视频| 日本五十路高清| 麻豆国产97在线/欧美| 日韩av在线大香蕉| 免费人成视频x8x8入口观看| 美女xxoo啪啪120秒动态图| 精品乱码久久久久久99久播| 99国产极品粉嫩在线观看| 国产伦精品一区二区三区四那| 婷婷丁香在线五月| 亚洲性夜色夜夜综合| 岛国在线免费视频观看| а√天堂www在线а√下载| 亚洲精华国产精华精| 性色avwww在线观看| 给我免费播放毛片高清在线观看| 成人性生交大片免费视频hd| 国产私拍福利视频在线观看| 国产蜜桃级精品一区二区三区| 91麻豆精品激情在线观看国产| 丰满人妻一区二区三区视频av| 亚洲精品粉嫩美女一区| 婷婷六月久久综合丁香| 一卡2卡三卡四卡精品乱码亚洲| 久久天躁狠狠躁夜夜2o2o| 蜜桃亚洲精品一区二区三区| 不卡视频在线观看欧美| 久久热精品热| 色综合站精品国产| 给我免费播放毛片高清在线观看| 成人美女网站在线观看视频| 国产老妇女一区| 成年女人看的毛片在线观看| 亚洲自拍偷在线| 女生性感内裤真人,穿戴方法视频| 18禁裸乳无遮挡免费网站照片| 欧美成人免费av一区二区三区| 欧美日韩国产亚洲二区| 国产乱人视频| 五月玫瑰六月丁香| 日韩高清综合在线| av黄色大香蕉| 真人做人爱边吃奶动态| 不卡一级毛片| 九色成人免费人妻av| 精品久久久久久久末码| 欧美国产日韩亚洲一区| 国产高清不卡午夜福利| 美女黄网站色视频| 在现免费观看毛片| 看片在线看免费视频| 成人av一区二区三区在线看| 亚洲美女视频黄频| 一边摸一边抽搐一进一小说| 少妇人妻精品综合一区二区 | 国内精品久久久久精免费| 麻豆成人av在线观看| 亚洲美女视频黄频| 午夜福利视频1000在线观看| 日韩大尺度精品在线看网址| 亚洲电影在线观看av| 免费av毛片视频| 狠狠狠狠99中文字幕| 久久精品国产99精品国产亚洲性色| av专区在线播放| 国内少妇人妻偷人精品xxx网站| 久久精品影院6| 久久国内精品自在自线图片| 联通29元200g的流量卡| 永久网站在线| 欧美黑人欧美精品刺激| 国产精品一区www在线观看 | 日韩欧美三级三区| 亚洲最大成人手机在线| 亚洲四区av| 国产在线精品亚洲第一网站| 亚洲中文字幕日韩| 亚洲一级一片aⅴ在线观看| 国产成人影院久久av| 亚洲男人的天堂狠狠| 亚洲天堂国产精品一区在线| 国产伦一二天堂av在线观看| 久久热精品热| 国产黄a三级三级三级人| 国产一区二区三区在线臀色熟女| 99riav亚洲国产免费| 直男gayav资源| 成年版毛片免费区| 一区福利在线观看| 亚洲四区av| 免费在线观看日本一区| 色噜噜av男人的天堂激情| 可以在线观看的亚洲视频| 久久久成人免费电影| 亚洲国产色片| 亚洲专区国产一区二区| 精品免费久久久久久久清纯| 窝窝影院91人妻| 久久精品久久久久久噜噜老黄 | 看片在线看免费视频| 欧美bdsm另类| 一进一出抽搐动态| 看十八女毛片水多多多| 亚洲av一区综合| 十八禁国产超污无遮挡网站| 在线观看一区二区三区| 麻豆av噜噜一区二区三区| 欧美三级亚洲精品| 免费人成视频x8x8入口观看| 国产伦在线观看视频一区| 国产私拍福利视频在线观看| 校园春色视频在线观看| 人人妻人人澡欧美一区二区| 麻豆国产av国片精品| 免费看光身美女| 五月玫瑰六月丁香| 精品久久久久久久久久久久久| av在线观看视频网站免费| 亚洲欧美清纯卡通| 国语自产精品视频在线第100页| 色综合站精品国产| 啦啦啦观看免费观看视频高清| www日本黄色视频网| 最近在线观看免费完整版| 天天一区二区日本电影三级| 黄色视频,在线免费观看| 老女人水多毛片| 亚洲精品日韩av片在线观看| 午夜激情欧美在线| 永久网站在线| 99riav亚洲国产免费| 国产精品99久久久久久久久| 男女边吃奶边做爰视频| 十八禁国产超污无遮挡网站| 国产精品久久久久久精品电影| 夜夜夜夜夜久久久久| 夜夜看夜夜爽夜夜摸| 少妇裸体淫交视频免费看高清| 国产高清三级在线| 国产免费av片在线观看野外av| 在线免费观看不下载黄p国产 | 精品无人区乱码1区二区| 免费av不卡在线播放| 国产精品爽爽va在线观看网站| 久久人妻av系列| 欧美日韩精品成人综合77777| 日韩 亚洲 欧美在线| 免费无遮挡裸体视频| 久久久久久久午夜电影| 成人综合一区亚洲| 国产精品无大码| 国产精品一区www在线观看 | 可以在线观看的亚洲视频| 波野结衣二区三区在线| 天美传媒精品一区二区| 99热精品在线国产| 美女高潮喷水抽搐中文字幕| 精品日产1卡2卡| 久久欧美精品欧美久久欧美| 久9热在线精品视频| 老熟妇乱子伦视频在线观看| 69人妻影院| 在线国产一区二区在线| 99热只有精品国产| 日本三级黄在线观看| 波多野结衣高清无吗| 男插女下体视频免费在线播放| 黄色一级大片看看| 少妇裸体淫交视频免费看高清| 国产精品电影一区二区三区| 亚州av有码| 如何舔出高潮| av中文乱码字幕在线| 在线免费十八禁| h日本视频在线播放| 偷拍熟女少妇极品色| 啦啦啦韩国在线观看视频| 久久婷婷人人爽人人干人人爱| 少妇人妻精品综合一区二区 | 欧美极品一区二区三区四区| 亚洲精品一区av在线观看| 久久久色成人| 精品一区二区三区视频在线观看免费| 成人无遮挡网站| 韩国av一区二区三区四区| 国产精品一区二区三区四区久久| 国产精品永久免费网站| 亚洲自拍偷在线| 午夜激情福利司机影院| 亚洲第一区二区三区不卡| 日本 欧美在线| 成人永久免费在线观看视频| 中文字幕熟女人妻在线| 亚洲av中文字字幕乱码综合| 亚洲av成人av| 精品久久久久久久久久免费视频| 午夜a级毛片| 久久草成人影院| 此物有八面人人有两片| 可以在线观看的亚洲视频| 不卡一级毛片| 国产欧美日韩一区二区精品| 永久网站在线| 日韩亚洲欧美综合| 国产午夜福利久久久久久| 嫩草影院入口| 欧美日韩瑟瑟在线播放| or卡值多少钱| 亚洲中文日韩欧美视频| 丰满乱子伦码专区| 看黄色毛片网站| 久久久久久伊人网av| 干丝袜人妻中文字幕| 又粗又爽又猛毛片免费看| 午夜精品在线福利| 12—13女人毛片做爰片一| a在线观看视频网站| 白带黄色成豆腐渣| 亚洲精华国产精华精| 国产毛片a区久久久久| 欧美性感艳星| 两性午夜刺激爽爽歪歪视频在线观看| 国产男人的电影天堂91| 老女人水多毛片| 国产精品久久久久久精品电影| 欧美日韩乱码在线| 伦精品一区二区三区| 综合色av麻豆| 18+在线观看网站| 九色成人免费人妻av| 国产色婷婷99| 国产精品国产三级国产av玫瑰| 亚洲国产精品久久男人天堂| 国内揄拍国产精品人妻在线| 亚洲久久久久久中文字幕| 欧美bdsm另类| 国产精品永久免费网站| 国产精品,欧美在线| 国内久久婷婷六月综合欲色啪| 久久精品久久久久久噜噜老黄 | 亚洲熟妇中文字幕五十中出| 啦啦啦观看免费观看视频高清| 五月玫瑰六月丁香| 免费人成在线观看视频色| 免费高清视频大片| 免费观看精品视频网站| 色综合站精品国产| 搡老熟女国产l中国老女人| 噜噜噜噜噜久久久久久91| 日韩一区二区视频免费看| 亚洲av成人av| 亚洲中文日韩欧美视频| 免费看美女性在线毛片视频| 特级一级黄色大片| а√天堂www在线а√下载| 天天一区二区日本电影三级| 男女视频在线观看网站免费| 国产午夜精品论理片| 十八禁网站免费在线| 精品日产1卡2卡| 日本免费一区二区三区高清不卡| 久久久久久国产a免费观看| 我要搜黄色片| 亚洲av电影不卡..在线观看| 久久久久久九九精品二区国产| 亚洲av熟女| 欧美一级a爱片免费观看看| 自拍偷自拍亚洲精品老妇| 91在线精品国自产拍蜜月| 九色国产91popny在线| 日韩欧美国产一区二区入口| 亚洲五月天丁香| 看黄色毛片网站| 国产精品,欧美在线| 如何舔出高潮| 日韩人妻高清精品专区| 亚洲精品影视一区二区三区av| 在线国产一区二区在线| 亚州av有码| 中文字幕免费在线视频6| 大型黄色视频在线免费观看| 99riav亚洲国产免费| 伊人久久精品亚洲午夜| 精品久久国产蜜桃| 午夜久久久久精精品| av天堂在线播放| 99riav亚洲国产免费| 黄色一级大片看看| 欧美性感艳星| 日韩一区二区视频免费看| 日韩精品有码人妻一区| 中文字幕精品亚洲无线码一区| 欧美又色又爽又黄视频| 欧美3d第一页| 在线观看一区二区三区| 久久人妻av系列| 国产亚洲91精品色在线| 日韩在线高清观看一区二区三区 | 非洲黑人性xxxx精品又粗又长| 又粗又爽又猛毛片免费看| 欧美色欧美亚洲另类二区| 最好的美女福利视频网| 极品教师在线视频| 国产伦精品一区二区三区视频9| 久久久久久久久久久丰满 | 免费黄网站久久成人精品| 久久亚洲精品不卡| 国产成人a区在线观看| 亚洲av成人av| 在线观看舔阴道视频| 免费观看人在逋| 伊人久久精品亚洲午夜| 最后的刺客免费高清国语| 欧美日韩黄片免| bbb黄色大片| 色在线成人网| 中文字幕免费在线视频6| 婷婷六月久久综合丁香| 人妻丰满熟妇av一区二区三区| 精品久久久久久久久av| 精品久久久久久久人妻蜜臀av| 露出奶头的视频| 国产精品不卡视频一区二区| 制服丝袜大香蕉在线| 亚洲最大成人中文| 精品久久久久久久末码| 国产女主播在线喷水免费视频网站 | 变态另类成人亚洲欧美熟女| 亚洲av日韩精品久久久久久密| 国产午夜精品论理片| 日韩av在线大香蕉| 国产一区二区三区视频了| 成人性生交大片免费视频hd| 国产一区二区三区av在线 | 夜夜夜夜夜久久久久| 少妇被粗大猛烈的视频| 国产aⅴ精品一区二区三区波| 亚洲内射少妇av| 18+在线观看网站| 91在线精品国自产拍蜜月| 三级男女做爰猛烈吃奶摸视频| 午夜福利视频1000在线观看| 国产精品久久久久久亚洲av鲁大| 国产v大片淫在线免费观看| 他把我摸到了高潮在线观看| 啦啦啦观看免费观看视频高清| 狂野欧美激情性xxxx在线观看| 伦精品一区二区三区| 亚洲精品影视一区二区三区av| 不卡一级毛片| 最新在线观看一区二区三区| 热99在线观看视频| 亚洲在线观看片| 亚洲av日韩精品久久久久久密| 一个人看视频在线观看www免费| 亚洲精华国产精华精| av在线蜜桃| 色精品久久人妻99蜜桃| 国产免费一级a男人的天堂| 午夜激情欧美在线| 最新中文字幕久久久久| av在线老鸭窝| 俄罗斯特黄特色一大片| 久久久久国产精品人妻aⅴ院| 性色avwww在线观看| 亚洲专区中文字幕在线| 成熟少妇高潮喷水视频| 在线免费观看的www视频| 亚洲av中文字字幕乱码综合| 日本a在线网址| 在线天堂最新版资源| 美女 人体艺术 gogo| 亚洲精品成人久久久久久| 九色国产91popny在线| 淫秽高清视频在线观看| 国产黄片美女视频| 可以在线观看的亚洲视频| 搡老岳熟女国产| 日本一本二区三区精品| 狂野欧美激情性xxxx在线观看| 伦理电影大哥的女人| 最后的刺客免费高清国语| 国产乱人视频| 九九爱精品视频在线观看| 亚洲国产欧洲综合997久久,| 日韩在线高清观看一区二区三区 | 国产精品永久免费网站| 久久精品综合一区二区三区| 欧美激情久久久久久爽电影| 听说在线观看完整版免费高清| 婷婷精品国产亚洲av| 国产人妻一区二区三区在| 亚洲真实伦在线观看| 中文字幕免费在线视频6| 十八禁网站免费在线| 亚洲,欧美,日韩| 一本久久中文字幕| 2021天堂中文幕一二区在线观| 在线播放国产精品三级| 亚洲av日韩精品久久久久久密| 伦理电影大哥的女人| 亚洲国产色片| 少妇丰满av| 大又大粗又爽又黄少妇毛片口| 天美传媒精品一区二区| 国内精品美女久久久久久| 久久久国产成人精品二区| 美女黄网站色视频| 色5月婷婷丁香| 高清毛片免费观看视频网站| 午夜免费成人在线视频| 亚洲第一电影网av| 噜噜噜噜噜久久久久久91| 变态另类成人亚洲欧美熟女| 热99re8久久精品国产| 一a级毛片在线观看| 中文字幕久久专区| 日本爱情动作片www.在线观看 | 哪里可以看免费的av片| 国产精品国产高清国产av| 日韩亚洲欧美综合| 黄片wwwwww| 一本精品99久久精品77| 91久久精品国产一区二区三区| 色综合色国产| 亚洲欧美日韩高清专用| 成人三级黄色视频| 国产高清三级在线| 日本爱情动作片www.在线观看 | 国产精品三级大全| 日韩精品中文字幕看吧| 麻豆成人午夜福利视频| 国产一区二区三区av在线 | 久久久久精品国产欧美久久久| 日韩强制内射视频| 全区人妻精品视频| 97超级碰碰碰精品色视频在线观看| 久久精品国产亚洲av香蕉五月| 欧美潮喷喷水| 亚洲精品影视一区二区三区av| 亚洲av中文av极速乱 | 亚洲精品粉嫩美女一区| 久久久久国产精品人妻aⅴ院| 国产欧美日韩一区二区精品| 一a级毛片在线观看| 天堂av国产一区二区熟女人妻| 色播亚洲综合网| 久久久久久九九精品二区国产| 精品久久久久久成人av| a级毛片a级免费在线| 欧美最新免费一区二区三区| 热99re8久久精品国产| 他把我摸到了高潮在线观看| 色5月婷婷丁香| 久久久色成人| 日韩精品有码人妻一区| 国产精品人妻久久久久久| 夜夜爽天天搞| 成人av一区二区三区在线看| 啦啦啦啦在线视频资源| 国产av不卡久久| 亚洲经典国产精华液单| 美女 人体艺术 gogo| 丰满人妻一区二区三区视频av| 国产成人一区二区在线| 久久久国产成人免费| 婷婷六月久久综合丁香| 色5月婷婷丁香| 国产欧美日韩一区二区精品| 精品国产三级普通话版| 亚洲精华国产精华液的使用体验 | 久久亚洲精品不卡| 久久国产精品人妻蜜桃| 波多野结衣巨乳人妻| a级毛片免费高清观看在线播放| 欧美激情国产日韩精品一区| 一边摸一边抽搐一进一小说| 悠悠久久av| 亚洲成a人片在线一区二区| 欧美一区二区精品小视频在线| 亚洲欧美日韩高清专用| 成人国产一区最新在线观看| 欧美色欧美亚洲另类二区| 成人三级黄色视频| 直男gayav资源| 亚洲av成人av| 尾随美女入室| 69av精品久久久久久| 欧美三级亚洲精品| 国产成人av教育| 少妇被粗大猛烈的视频| 日本一本二区三区精品| 嫩草影院精品99| 国产精品一区二区三区四区免费观看 | 波多野结衣高清无吗| 欧美日韩亚洲国产一区二区在线观看| 日日撸夜夜添| 亚洲三级黄色毛片| 欧美性感艳星| 变态另类成人亚洲欧美熟女| 简卡轻食公司| 一级a爱片免费观看的视频| 色尼玛亚洲综合影院| 免费高清视频大片| 午夜日韩欧美国产|