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

    大氣壓環(huán)環(huán)電極結(jié)構(gòu)射流放電模型建立及仿真

    2016-04-07 00:35:51姚正秋南京工業(yè)大學(xué)電氣工程與控制科學(xué)學(xué)院南京210009
    電工技術(shù)學(xué)報(bào) 2016年4期

    方 志 錢 晨 姚正秋(南京工業(yè)大學(xué)電氣工程與控制科學(xué)學(xué)院 南京 210009)

    ?

    大氣壓環(huán)環(huán)電極結(jié)構(gòu)射流放電模型建立及仿真

    方 志 錢 晨 姚正秋
    (南京工業(yè)大學(xué)電氣工程與控制科學(xué)學(xué)院 南京 210009)

    摘要建立射流放電的仿真模型并對(duì)其進(jìn)行研究,該工作對(duì)于深入研究其放電特性,進(jìn)而在實(shí)際應(yīng)用中優(yōu)化射流放電等離子反應(yīng)器設(shè)計(jì)具有重要意義。建立了大氣壓下環(huán)環(huán)電極結(jié)構(gòu)射流放電系統(tǒng),采用電壓、電流波形和Lissajous圖測(cè)量及發(fā)光圖像拍攝等實(shí)驗(yàn)手段,研究了其放電特性。結(jié)合對(duì)實(shí)驗(yàn)結(jié)果和放電特性的分析,建立了一種環(huán)環(huán)電極結(jié)構(gòu)射流放電的等效電氣模型,將兩電極之間的放電等效為DBD放電模型,將噴出管外的射流等效為可變阻抗,同時(shí)還考慮了溢流效應(yīng)的影響,將其等效為可變阻抗,更真實(shí)準(zhǔn)確地反映了射流放電的實(shí)際情況?;诖穗姎饽P?,進(jìn)一步得到了放電等效電路,并利用Simulink建立了射流放電的動(dòng)態(tài)仿真模型。仿真得到的電壓、電流波形圖和Lissajous圖形與實(shí)驗(yàn)結(jié)果對(duì)比顯示,二者是吻合的,驗(yàn)證了提出的等效電氣模型的準(zhǔn)確性。在此基礎(chǔ)上,進(jìn)一步仿真研究了電源頻率和環(huán)環(huán)電極間距對(duì)電氣特性及放電參量的影響,結(jié)果表明,放電功率和傳輸電荷都隨電源頻率的增加非線性地增大,隨環(huán)環(huán)間距的增加非線性減小。還利用所建立的模型進(jìn)一步得到了實(shí)驗(yàn)過程中無法直接測(cè)量獲得的介質(zhì)電壓、氣隙電壓和放電電流等放電參量。

    關(guān)鍵詞:射流放電 環(huán)環(huán)電極結(jié)構(gòu) 等效電氣模型 動(dòng)態(tài)仿真模型 仿真研究

    國(guó)家自然科學(xué)基金(51377075)和江蘇省自然科學(xué)基金(BK20131412)資助項(xiàng)目。

    The Model and Simulation Studies for Ring-Ring Electrode Structure Jet Discharge at Atmospheric Pressure

    Fang Zhi Qian Chen Yao Zhengqiu
    (College of Electrical Engineering and Control Science Nanjing Tech University Nanjing 210009 China)

    Abstract The model establishment for jet discharge and its simulation analysis are of great interesting for studying the discharge characteristics as well as for optimizing the design of jet discharge plasma reactors in real applications. In this paper, the experimental system with the ring-ring electrode structure for atmospheric pressure plasma jet discharge is established, and the discharge characteristics are studied, by measuring voltage and current waveforms, Lissajous figures and lighting emission images. Based on the experimental results and discharge characteristics analysis, an equivalent electrical model for ring-ring electrode structure jet discharge is established, which can reflect the discharge. The discharge space between two electrodes is equivalent to a DBD discharge model, and the jet outside the tube is equivalent to variable impedance. The influence of spillover effect is also taken into account, and it is modeled as variable impedance. Based on the equivalent electrical model, an equivalent circuit diagram is deduced, and a dynamic simulation model is established in Simulink software. Voltage and current waveforms, and Lissajous figures can be well coherence withthe experimental results, which verify the proposed electrical model. The influences of the power supply frequency and ring-ring distance on the electrical characteristics and discharge parameters are also studied. Results show that the discharge power and transported charges both increase nonlinearly as the applied frequency increases, while both decrease as the ring-ring distance increases. The dynamic behaviors of discharge parameters that can not measured directly during the experiment, such as dielectric voltage, gas gap voltage and discharge current, are also obtained by the established model.

    Keywords:Jet discharge, ring-ring electrode structure, equivalent electrical model, dynamic simulation model, simulation studies

    0 引言

    近年來,He、Ar等純惰性氣體中產(chǎn)生的大氣壓低溫等離子體射流受到廣泛關(guān)注[1-5]。與其他形式的大氣壓低溫等離子體源相比,其能量密度適中,易于操作,成本低廉,反應(yīng)器結(jié)構(gòu)靈活,處理區(qū)域不受放電空間限制,可處理三維物體,在殺菌消毒、等離子體醫(yī)學(xué)及材料處理等領(lǐng)域具有廣泛的應(yīng)用前景[6-10]。因此,對(duì)大氣壓等離子體射流的放電特性進(jìn)行研究對(duì)促進(jìn)其應(yīng)用具有重要意義。

    目前,研究射流放電大多采用實(shí)驗(yàn)手段,由于受外部運(yùn)行參數(shù)、反應(yīng)器結(jié)構(gòu)及氣體條件等多種因素影響,實(shí)驗(yàn)工作量較大,且受到測(cè)量和診斷手段的限制,許多微觀的放電參量無法直接測(cè)量得到。而通過建立合適的等效電氣模型來模擬其放電情況,進(jìn)而通過仿真研究其放電特性及放電參量的變化規(guī)律,不僅可以減小實(shí)驗(yàn)工作量,還可以得到實(shí)驗(yàn)無法直接獲得的放電參量,從而對(duì)放電特性進(jìn)行更深入的研究,為優(yōu)化反應(yīng)器設(shè)計(jì)提供參考。環(huán)環(huán)電極結(jié)構(gòu)是一種產(chǎn)生射流的典型電極結(jié)構(gòu),國(guó)內(nèi)外已有一些關(guān)于其放電特性的研究[11-15]。江南等實(shí)驗(yàn)采用光電倍增管和數(shù)碼相機(jī)等研究了環(huán)環(huán)電極結(jié)構(gòu)射流的電壓、電流特性和發(fā)光特性的變化規(guī)律,并研究了放電參數(shù)對(duì)“溢流”現(xiàn)象的影響[15]。J. S. Oh等采用電流探頭測(cè)量了環(huán)環(huán)電極結(jié)構(gòu)射流的放電電流及傳輸電荷等放電參量[13]。S. Yonemori等采用激光誘導(dǎo)熒光技術(shù)對(duì)環(huán)環(huán)電極結(jié)構(gòu)射流中粒子進(jìn)行診斷[12]。關(guān)于環(huán)環(huán)電極結(jié)構(gòu)射流放電仿真研究國(guó)內(nèi)外報(bào)道較少,張冠軍教授采用結(jié)合氣體動(dòng)力學(xué)的仿真模型與實(shí)驗(yàn)相結(jié)合的辦法,分析了管徑及電壓對(duì)雙環(huán)電極射流長(zhǎng)度的影響[16]。研究者通常采用與射流放電類似的介質(zhì)阻擋放電(Dieletric Barrier Discharge, DBD)的電氣模型來分析射流放電的電氣特性,但射流放電不同于DBD放電,除兩電極間發(fā)生放電外,產(chǎn)生的等離子體還會(huì)噴出放電電極空間,使周圍的空氣電離,在一定條件下,射流放電還存在“電荷溢流”效應(yīng)[15],因此采用DBD電氣模型來模擬分析射流放電尚存在一定誤差,得到的結(jié)果往往不能真實(shí)地反映射流放電的電氣特性??偟膩碚f,目前已報(bào)道環(huán)環(huán)電極機(jī)構(gòu)射流放電的研究結(jié)果大多是實(shí)驗(yàn)研究,而通過建立更能反映放電真實(shí)情況的等效電氣模型對(duì)環(huán)環(huán)電極結(jié)構(gòu)射流放電電氣特性進(jìn)行仿真研究,國(guó)內(nèi)外研究報(bào)道不多。

    本文通過對(duì)大氣壓環(huán)環(huán)電極結(jié)構(gòu)Ar等離子體射流放電特性分析的基礎(chǔ)上,提出一種能反映放電實(shí)際情況的等效電氣模型,并進(jìn)一步得到等效電路,進(jìn)而基于此等效電路采用Simulink建模工具建立其動(dòng)態(tài)仿真模型,仿真得到其放電電氣特性,和實(shí)驗(yàn)結(jié)果比較驗(yàn)證其準(zhǔn)確性,并進(jìn)一步利用該模型對(duì)環(huán)環(huán)電極結(jié)構(gòu)射流放電的影響因素進(jìn)行仿真研究,以及研究了其放電參量的動(dòng)態(tài)變化規(guī)律。

    1 實(shí)驗(yàn)裝置及測(cè)量系統(tǒng)

    圖1給出了研究所采用的射流實(shí)驗(yàn)裝置及測(cè)量系統(tǒng)接線。電源采用幅值在0~20kV范圍內(nèi)、頻率在5~20kHz范圍內(nèi)可調(diào)的高頻高壓交流電源。射流反應(yīng)器部分由石英玻璃管、高壓電極和地電極組成。其中,石英玻璃管外徑4mm,內(nèi)徑2mm,長(zhǎng)180mm。高壓電極和地電極都為寬10mm的銅皮,均套在玻璃管外壁,地電極外端距離管口10mm,高低壓電極間距離為20mm。通過流量計(jì)使高純度(99.999%)Ar由進(jìn)氣口進(jìn)入玻璃管內(nèi),固定其流速為6L/min。外加電壓波形采用高壓探頭Tek P6015A來測(cè)量,放電電流波形通過在放電回路中串聯(lián)一個(gè)阻值為50Ω的無感電阻R獲得,放電空間傳輸?shù)碾姾赏ㄟ^在放電回路中串聯(lián)一個(gè)0.1μF的測(cè)量電容C獲得,放電Lissajous圖形通過把高壓探頭測(cè)得的反應(yīng)器上的電壓和C兩端的電壓分別加在示波器的XY軸得到。實(shí)驗(yàn)時(shí)測(cè)得的電壓、電流波形及Lissajous圖形由TDS—3054c數(shù)字示波器記錄。發(fā)光圖像用置于放電空間側(cè)面與放電氣隙平行的數(shù)碼相機(jī)Canon G6(曝光時(shí)間為1s)拍攝得到。

    圖1 射流實(shí)驗(yàn)裝置及測(cè)量系統(tǒng)示意圖Fig.1 Schematic diagram of jet experimental set-up and measurement system

    2 射流放電演變規(guī)律

    圖2和圖3分別給出了固定電源頻率為10kHz時(shí)不同外加電壓下測(cè)得的射流放電的電壓、電流波形和Lissajous圖形以及發(fā)光圖像。如圖2a所示,電壓幅值為11kV時(shí),電流波形表現(xiàn)為多脈沖形式,電壓每個(gè)正、負(fù)半周期內(nèi)各出現(xiàn)四個(gè)電流脈沖,其持續(xù)時(shí)間在μs數(shù)量級(jí),電流脈沖幅值為十幾到幾十mA,Lissajous圖形中,其正、負(fù)半周期的放電階段各有四個(gè)明顯階躍。如圖2b~圖2d所示,隨電壓幅值由12kV增大到14kV,電壓每半個(gè)周期內(nèi)電流脈沖的個(gè)數(shù)逐漸增多,在Lissajous圖形中,其正、負(fù)半周期的階躍個(gè)數(shù)也隨之增加。從圖3a中可以看出,射流放電發(fā)生后,放電貫穿高壓電極和地電極之間的氣隙,地電極外側(cè)等離子體在玻璃管內(nèi)向管口方向延伸至管外,在敞開的空氣中形成大氣壓等離子體射流;同時(shí)可以觀察到射流放電還存在“電荷溢流”效應(yīng)[16-18],即高壓電極下的電荷積累區(qū)向內(nèi)發(fā)展,產(chǎn)生等離子體電離通道,其向內(nèi)發(fā)展的長(zhǎng)度也隨外加電壓的增加而增加。從圖3中可以看出,射流的長(zhǎng)度隨外加電壓增加而增加,說明在空氣中產(chǎn)生的電離通道長(zhǎng)度增加。隨外加電壓由12kV增加到14kV,射流長(zhǎng)度由10mm增加到20mm;高壓電極內(nèi)側(cè)因溢流效應(yīng)產(chǎn)生的等離子體長(zhǎng)度由6mm增加到8mm。

    圖2 不同外加電壓下測(cè)得的射流放電的電壓、電流波形和Lissajous圖形Fig.2 Voltage and current waveforms and Lissajous figures of jet discharge measured at different applied voltages

    圖3 射流放電的發(fā)光圖像隨外加電壓變化情況Fig.3 Changing of lighting emission images of jet discharge with applied voltage

    3 仿真模型建立

    3.1 等效電氣模型

    射流放電是以兩電極之間放電擊穿和電離通道向兩電極外側(cè)擴(kuò)展為主要特征的氣體放電,其放電電離過程主要發(fā)生在如上放電通道中,因此射流放電的電氣仿真模型主要建立在對(duì)其放電電離通道的等效基礎(chǔ)上。結(jié)合上文的實(shí)驗(yàn)結(jié)果和分析,本文建立了如圖4所示的等效電氣模型,由于激勵(lì)電源為高壓交流高頻電源,因此放電未發(fā)生時(shí),放電通道未形成,這個(gè)射流電極結(jié)構(gòu)的各個(gè)部分等效為電容,放電前的等效電氣模型如圖4a所示,圖中Cd1和Cd2表示阻擋介質(zhì)(石英玻璃管)等效電容;Cg表示兩電極間氣隙等效電容;Cp1表示回流等離子體等效電容;Cp2表示等離子體射流等效電容;Cs、Rs和Ls回路分別表示反應(yīng)器雜散電容、引線電阻和寄生串聯(lián)電感,考慮這些因素更能真實(shí)地反映放電實(shí)際情況。外加電壓未達(dá)到擊穿電壓時(shí),如上各元件參數(shù)均為定值。當(dāng)外加電壓增加到一定值時(shí),高低壓電極間的Ar氣間隙擊穿,這時(shí)兩電極之間的放電表現(xiàn)為介質(zhì)阻擋放電(DBD),根據(jù)本文作者以前的研究成果,可以用一個(gè)電壓控制電流源(Voltage Controlled Current Source, VCCS)和一個(gè)阻容串聯(lián)回路(Cg和Rg)等效,其中VCCS反映放電電流變化規(guī)律,而Cg和Rg反映放電擊穿后兩電極間等離子通道阻抗變化[19-21],由于采用交流電源激勵(lì),外加電壓在每個(gè)半周會(huì)改變極性,因此雙環(huán)電極的兩個(gè)電極每隔半周極性交替變化,而氣流方向一直由管內(nèi)流向管口,因此會(huì)導(dǎo)致外加電壓正、負(fù)半周兩電極之間電離擊穿通道的阻抗有所不同,反映在放電電流上,正、負(fù)半周測(cè)量得到的電流幅值有所不同??紤]這一因素,在氣隙回路中并聯(lián)了Rf1和Rf2分別表示正、負(fù)半周期極性效應(yīng)引起的阻抗變化,如圖4b所示。而同時(shí)放電空間產(chǎn)生的粒子隨氣流移動(dòng)到低壓電極外的放電空間,使空氣電離形成電離通道,產(chǎn)生等離子射流,改變了其阻抗。此時(shí)地電極外放電空間阻抗不再是容性Cp1,相關(guān)研究表明,隨射流在空氣中擴(kuò)展,其在空氣中形成的導(dǎo)電通道的阻抗和電位分布不同[2,5],因此,射流部分的等離子體阻抗可用可變的Cp2(t)和Rp2(t)組成的阻容性串聯(lián)電路表示。同理,內(nèi)側(cè)因溢流效應(yīng)而向內(nèi)產(chǎn)生的電離通道可用Cp1(t)和Rp1(t)串聯(lián)等效,如圖4b所示。由圖3中實(shí)驗(yàn)結(jié)果和相關(guān)研究可知[2,5],隨著外加電壓的升高,電離增強(qiáng),射流和溢流效應(yīng)產(chǎn)生的等離子體長(zhǎng)度不斷增加,Cp1(t)和Cp2(t)不斷減小,Rp1(t)和Rp2(t)不斷增大。因此,本文考慮了射流通道和溢流放電通道的阻抗并采用可變阻容串聯(lián)電路來等效射流和溢流效應(yīng),更能真實(shí)地反映放電特性。

    圖4 等效電氣模型Fig.4 Equivalent electrical model

    3.2 等效電路

    圖5 等效電路Fig.5 Equivalent circuit diagram

    根據(jù)圖4中等效電氣模型及相關(guān)分析,考慮到電壓、電流關(guān)系,進(jìn)一步建立了放電不同階段的等效電路,如圖5所示。圖中各個(gè)元器件與等效電氣模型中的元器件相對(duì)應(yīng),考慮到實(shí)驗(yàn)中外電路雜散電阻和變壓器等效漏抗等的影響,還在外加電源旁串聯(lián)了電阻R和電感L。其中,Va(t)表示外加電壓;Vt(t)表示反應(yīng)器電壓;ia(t)表示回路總電流;it(t)表示放電回路總電流;iccs(t)表征DBD的電壓控制電流源;id1(t)和id2(t)分別表示兩個(gè)阻擋介質(zhì)上的位移電流;ip1(t)和ip2(t)分別表示等離子體射流和溢流等離子體中的位移電流;ig(t)表示氣隙中的位移電流;Vd1(t)和Vd2(t)分別表示兩個(gè)阻擋介質(zhì)上的電壓;Vg(t)表示氣隙上的電壓。

    圖5b中的射流放電發(fā)生后等效電路由于含有Cp1(t)、Cp2(t)、Rp1(t)和Rp2(t)等非線性元件和電壓控制電流源,分析起來十分復(fù)雜。考慮到雙環(huán)電極結(jié)構(gòu)射流放電,在外加電壓增加到擊穿電壓時(shí),放電先是在兩電極之間產(chǎn)生DBD放電,隨外加電壓繼續(xù)增加到一定值,射流放電和溢流效應(yīng)發(fā)生。因此將上述電路進(jìn)行一定簡(jiǎn)化,先不考慮射流和溢流,根據(jù)基爾霍夫定律,可以得到各電壓Vt(t)、Vd1(t)、Vd2(t)和Vg(t)之間的關(guān)系,以及它們與各電流it(t)、iccs(t)之間的關(guān)系[25]

    經(jīng)過計(jì)算和整理得到iccs(t)、Vt(t)和it(t)之間的關(guān)系為

    由式(5)可知,放電電流iccs(t)受Vt(t)控制,CCS是由遵循Cg非線性動(dòng)態(tài)變化規(guī)律的電壓信號(hào)激勵(lì),對(duì)任一iccs(t)均有一個(gè)it(t)與之對(duì)應(yīng)。得到上述關(guān)系后,再考慮Cp1(t)、Cp2(t)、Rp1(t)和Rp2(t)等接入,通過選取合適的值,來反應(yīng)射流放電和溢流效應(yīng)。

    3.3 動(dòng)態(tài)仿真模型

    根據(jù)圖5中的等效電路,利用Simulink建立了如圖6所示的射流放電的動(dòng)態(tài)仿真模型,來對(duì)圖5中放電的電氣特性進(jìn)行仿真。模型中使用五個(gè)開關(guān)控制模塊Subsystem1、Subsystem2、Subsystem3、Subsystem4和Subsystem5控制使能開關(guān)來實(shí)現(xiàn)放電的發(fā)生和熄滅。射流放電未發(fā)生時(shí),開關(guān)Sw1導(dǎo)通,開關(guān)Sw4和Sw5斷開,開關(guān)Sw6和Sw7導(dǎo)通,模型中放電回路為介質(zhì)等效電容Cd1和Cd2先分別與回流等離子體等效電容Cp1和等離子體射流等效電容Cp2并聯(lián)之后再串聯(lián),然后再與氣隙等效電容Cg串聯(lián),與圖5a一致。射流放電發(fā)生后,開關(guān)Sw4和Sw5導(dǎo)通,開關(guān)Sw6和Sw7斷開,在Cd1兩端并聯(lián)了Cp1(t)和Rp1(t),它們的值為線性變化的函數(shù),反映了隨外加電壓變化放電過程中回流等離子體等效阻抗的變化過程,由于它們的值不能通過測(cè)量或計(jì)算得到,本文通過實(shí)驗(yàn)和仿真結(jié)果比較得到它們的取值和初始值,確定方法為:在不同外加電壓下選擇不同值的Cp1(t)和Rp1(t)進(jìn)行仿真,將得到的仿真結(jié)果和實(shí)驗(yàn)結(jié)果進(jìn)行比較,調(diào)節(jié)Cp1(t)和Rp1(t)的值,直到仿真和實(shí)驗(yàn)結(jié)果一致,將不同外加電壓幅值和對(duì)應(yīng)的Cp1(t)和Rp1(t)數(shù)據(jù)進(jìn)行曲線擬合,得到其函數(shù)表達(dá)式分別為

    圖6 射流放電的動(dòng)態(tài)仿真模型Fig.6 Dynamic simulation model of jet discharge

    式中,Va為外加電壓幅值,其與圖2中的Va相對(duì)應(yīng);D1為衰減常數(shù),D1=-2×10-16;Cp1(t)的初始值為1×10-8F,其值通過D1=0時(shí)比較仿真和實(shí)驗(yàn)結(jié)果確定;A1為增益常數(shù),A1=1×10-2;Rp1(t)的初始值為50Ω,其值通過A1=0時(shí)比較仿真和實(shí)驗(yàn)結(jié)果確定。同時(shí),在Cd2兩端并聯(lián)了Cp2(t)和Rp2(t),其值也為線性變化的函數(shù),反映了放電過程中等離子體射流溢流現(xiàn)象導(dǎo)致的等效阻抗的變化過程,其取值和初始值的確定方法同Cp1(t)和Rp1(t),函數(shù)表達(dá)式分別為

    式中,Cp2(t)的初始值為1μF;D2為衰減常數(shù),D2= -3×10-16;Rp2(t)的初始值為5kΩ;A2為增益常數(shù),A2=300。

    圖7 模塊4和模塊5Fig.7 The module of Subsystem4 and Subsystem5

    模型中放電電流的變化用如圖7a所示的電壓控制電流源模塊Subsystem4來體現(xiàn),而考慮了放電正、負(fù)半周極性效應(yīng)的兩電極之間電離通道阻抗,用如圖7b等效回路模塊Subsystem5來反映。Subsystem4其內(nèi)部結(jié)構(gòu)展開如圖7a所示,該模塊符合式(5),其功能為根據(jù)反應(yīng)器電容值以及Vt(t)和it(t)得到電壓控制電流源激勵(lì)信號(hào)iccs(t)。Subsystem5其內(nèi)部結(jié)構(gòu)展開如圖7b所示,電壓正半周放電時(shí),Sw2導(dǎo)通,Sw3斷開,在Rg(t)和Cg(t)兩端并聯(lián)正半周期極性效應(yīng)等效電阻Rf1。電壓負(fù)半周放電時(shí),Sw2斷開,Sw3導(dǎo)通,在Rg(t)和Cg(t)兩端并聯(lián)負(fù)半周期極性效應(yīng)等效電阻Rf2,與圖5b一致,其中Rf1和Rf2的值分別為1kΩ和1.5kΩ。而模塊Subsystem1和Subsystem2為開關(guān)控制模塊,產(chǎn)生脈沖信號(hào)用于控制圖7b中開關(guān)Sw2和Sw3的開閉,進(jìn)而實(shí)現(xiàn)每個(gè)周期放電的開始與熄滅,來控制每個(gè)放電周期出現(xiàn)脈沖放電電流,通過調(diào)節(jié)產(chǎn)生脈沖信號(hào)的脈寬和脈沖個(gè)數(shù)使放電電流持續(xù)時(shí)間與實(shí)際放電電流一致。根據(jù)環(huán)環(huán)電極和環(huán)板電極的射流反應(yīng)器結(jié)構(gòu),兩電極間加上電壓產(chǎn)生電場(chǎng)后,玻璃管相當(dāng)于同軸電容器,而管內(nèi)氣體相當(dāng)于一個(gè)圓柱型電容器,等效電氣模型中的仿真參數(shù)Cd1、Cd2和Cg和可以由以下公式計(jì)算得到。

    式中,ε0為真空介電常數(shù);εd為介質(zhì)相對(duì)介電常數(shù);εg為氣體相對(duì)介電常數(shù);ld為電極寬度;lg為介質(zhì)層厚度;R和r分別為玻璃管的內(nèi)、外半徑。仿真時(shí),外電路雜散參數(shù)R和L分別取500Ω和1mH,反應(yīng)器雜散參數(shù)Cs、Rs和Ls分別取500Ω、0.1mH和100pF。其中,Cs在pF數(shù)量級(jí),對(duì)仿真結(jié)果影響較大,其值根據(jù)實(shí)驗(yàn)測(cè)量結(jié)果估算得到。其確定方法為:利用圖1中裝置在氣隙放電擊穿前測(cè)量得到電源電壓V(t)和相應(yīng)的回路電流i(t),由于此時(shí)放電未發(fā)生,i(t)僅為容性電流,利用關(guān)系式i(t)=CeqdV(t)/dt,可求得等效電容Ceq,Ceq為反應(yīng)器電容和雜散電容并聯(lián)的總電容,根據(jù)計(jì)算的反應(yīng)器電容,可估算出Cs的值約為94.6pF。而R、L、Rs和Ls等對(duì)仿真結(jié)果影響不大,為了真實(shí)反映引線電阻等雜散參數(shù)的影響,在仿真模型中也加以考慮。放電發(fā)生階段,前一個(gè)半周期放電時(shí)介質(zhì)表面積聚的電荷對(duì)下一個(gè)半周期的放電起促進(jìn)作用,導(dǎo)致起始放電電壓降低,放電提前發(fā)生。考慮到這種情況,模型中采用“Transport Delay”模塊對(duì)電壓進(jìn)行了延遲,延遲的時(shí)間函數(shù)用Td(t)表示,其函數(shù)表達(dá)式為式中,A3為時(shí)間增益常數(shù),A3=2×10-9;Td(t)的初始值為4.5×10-5,其值通過A3=0時(shí)比較仿真和實(shí)驗(yàn)結(jié)果確定。

    模型中還添加了多個(gè)測(cè)量模塊和顯示模塊,來得到和顯示放電電氣特性,為了能在所建立的模型中得到Lissajous圖形,放電回路中串聯(lián)一個(gè)電容C0=0.1μF。利用電壓表V1來測(cè)量電容C0兩端的的電壓,電壓表V2測(cè)量電源電壓Va(t),電流表I1測(cè)量回路總電流ia(t),電源電壓和整個(gè)放電回路上的總電流連接到顯示器“VI”兩端,得到電壓、電流波形。電源電壓和C0兩端電壓連接到顯示器“Lissajous”上,得到放電Lissajous圖形。

    4 仿真結(jié)果

    采用圖6所建立的仿真模型對(duì)環(huán)環(huán)電極結(jié)構(gòu)射流放電電氣特性進(jìn)行仿真,射流反應(yīng)器結(jié)構(gòu)及電源參數(shù)等實(shí)驗(yàn)條件與前面實(shí)驗(yàn)部分介紹相同。仿真中用到的參量Cg、Cd1、Cd2和Cd可根據(jù)反應(yīng)器結(jié)構(gòu)由計(jì)算公式得到[23],它們的值分別為3.4pF、200.04pF、200.04pF和100.02pF。圖8給出了不同外加電壓下仿真得到的射流放電的電壓、電流波形和Lissajous圖形。比較圖8和圖2中本文仿真和實(shí)驗(yàn)結(jié)果可知,兩者在電流脈沖個(gè)數(shù)、脈沖幅值及脈沖持續(xù)時(shí)間上均符合。比較本文的仿真結(jié)果與相關(guān)文獻(xiàn)中與本文類似實(shí)驗(yàn)條件下得到的實(shí)驗(yàn)測(cè)量結(jié)果,本文仿真得到的電壓、電流波形的特征和隨外加電壓的變化規(guī)律與相關(guān)文獻(xiàn)中結(jié)果是一致的[2,5,15]。說明了本文所建立的電氣模型和動(dòng)態(tài)仿真模型的正確性。

    圖8 不同外加電壓下仿真得到的射流放電的電壓、電流波形和Lissajous圖形Fig.8 Voltage and current waveforms and Lissajous figures of jet discharge obtained by simulations at different applied voltages

    放電功率P和傳輸電荷Q是表征放電的兩個(gè)重要參量,其中,P可由放電Lissajous圖獲得,而Q可由測(cè)量得到的電壓、電流波形計(jì)算得到[24-28]。它們的計(jì)算式分別為

    式中,f為外加電壓的頻率;A為L(zhǎng)issajous圖形的面積;i(t)為瞬時(shí)放電電流;T為一個(gè)放電周期。

    根據(jù)圖2和圖8中實(shí)驗(yàn)和仿真結(jié)果,由式(13)和式(14)計(jì)算得到P和Q,如圖9和圖10所示。當(dāng)電壓幅值由11kV增加到14kV時(shí),實(shí)驗(yàn)得到的P 從13.6W增加到23.5W,Q從507.6nC增加到958.8nC。仿真得到的P從12.5W增加到22.5W,Q 從496.6nC增加到942.8nC,P和Q都隨著外加電壓的增加而非線性地增大。由仿真結(jié)果計(jì)算得到的P與Q的變化趨勢(shì)和數(shù)值與由實(shí)驗(yàn)結(jié)果計(jì)算得到的基本一致。

    圖9 實(shí)驗(yàn)和仿真得到的放電功率Fig.9 Discharge power obtained by experiment and simulation

    圖10 實(shí)驗(yàn)和仿真得到的傳輸電荷Fig.10 Transported charge obtained by experiment and simulation

    5 放電影響因素的仿真

    外部運(yùn)行條件、反應(yīng)器結(jié)構(gòu)參數(shù)及氣體條件等都會(huì)影響射流放電特性,實(shí)驗(yàn)研究中各種影響因素的變化范圍大,實(shí)驗(yàn)研究工作量較大。本文通過仿真研究了電源頻率和電極間距對(duì)放電電氣特性的影響,從而為反應(yīng)器的進(jìn)一步優(yōu)化和設(shè)計(jì)提供參考。

    5.1 電源頻率的影響

    仿真條件為:外加電壓幅值保持為11kV,環(huán)環(huán)間距固定為20mm,電源頻率在9.5~11kHz內(nèi)變化,Cg、Cd1、Cd2和Cd的值與第4節(jié)給出的值相同。圖11給出了電源頻率分別為9.5kHz、10.5kHz和11kHz時(shí)仿真得到的的電壓、電流波形和Lissajous圖形,電源頻率為10kHz時(shí)的仿真結(jié)果如圖8a所示??梢钥闯?,當(dāng)電源頻率由9.5kHz增加到11kHz時(shí),電流脈沖最大幅值由9mA增加到11.5mA。利用仿真得到的Lissajous圖形,由式(13)和式(14)計(jì)算得到P和Q隨電源頻率的變化曲線如圖12所示。當(dāng)電源頻率由9.5kHz增加到11kHz時(shí),P和Q都隨著電源頻率的增加非線性增大,P從11.9W增加到14.1W,Q從481.7nC增加到538.9nC。

    圖11 不同電源頻率下仿真得到的電壓、電流波形和Lissajous圖形Fig.11 Voltage and current waveforms and Lissajous figures obtained by simulations at different frequency

    圖12 仿真得到的放電功率與傳輸電荷隨電源頻率變化曲線Fig.12 Variation of the discharge power and the transported charge with applied frequency obtained by simulations

    5.2 環(huán)環(huán)間距的影響

    仿真條件為:外加電壓幅值保持為11kV,電源頻率保持為10kHz,環(huán)環(huán)間距在10~30mm內(nèi)變化,仿真用到的參數(shù)Cd1、Cd2和Cd的值與第4節(jié)給出的值相同,10mm和30mm間距下Cg的經(jīng)計(jì)算分別為4.95pF和2.59pF。圖13給出了環(huán)環(huán)間距分別為10mm和30mm時(shí)仿真得到的電壓、電流波形和Lissajous圖形,環(huán)環(huán)間距為20mm時(shí)的仿真結(jié)果如圖8a所示??梢钥闯觯S著環(huán)環(huán)間距的增大,電流脈沖幅值有所下降,放電起始時(shí)刻延遲。利用仿真得到的Lissajous圖形,由式(13)和式(14)計(jì)算得到的P和Q隨環(huán)環(huán)間距的變化曲線如圖14所示。當(dāng)環(huán)環(huán)間距由10mm增加到30mm時(shí),P和Q都隨著環(huán)環(huán)間距的增加而非線性減小,P從13.2W減小到11.6W,Q從518.9nC減小到470.7nC。

    圖13 不同環(huán)環(huán)間距下仿真得到的電壓、電流波形和Lissajous圖形Fig.13 Voltage and current waveforms and Lissajous figures obtained by simulations at different ring-ring distance

    圖14 仿真條件下得到的放電功率與傳輸電荷隨環(huán)環(huán)間距變化曲線Fig.14 Variation of the discharge power and the transported charge with ring-ring distance obtained by simulations

    6 放電參量計(jì)算

    利用本文所建立的仿真模型,還可以計(jì)算得到實(shí)驗(yàn)無法直接測(cè)量的放電參量,如電源電壓Va、介質(zhì)電壓Vd、氣隙電壓Vg和放電電流iccs等,從而深入地對(duì)放電過程進(jìn)行分析。圖15給出了外加電壓為11kV時(shí)計(jì)算得到的放電參量。從圖15中可以看出放電未發(fā)生時(shí),Vg隨著外加電壓的增大而增大。當(dāng)Vg達(dá)到6kV時(shí),放電發(fā)生,放電電流脈沖幅值持續(xù)增加,由于介質(zhì)表面聚積電荷,形成一個(gè)反向電壓,從而導(dǎo)致Vg減小,出現(xiàn)第一次抖動(dòng),當(dāng)放電空間內(nèi)的電壓小于擊穿電壓時(shí),放電熄滅。此后,由于外加電壓上升,氣隙電壓繼續(xù)上升,放電重燃,此后再次熄滅,如此反復(fù),抖動(dòng)次數(shù)增加,11kV條件下正、負(fù)半周各有4次。Vg抖動(dòng)次數(shù)對(duì)應(yīng)電壓正、負(fù)半周期電流脈沖的個(gè)數(shù),同時(shí)與圖2中的電流脈沖個(gè)數(shù)一致,而阻擋介質(zhì)集聚電荷的作用是影響放電發(fā)生和熄滅的主要因素。

    圖15 仿真得到的放電參量Fig.15 Discharge parameters obtained by simulations

    7 結(jié)論

    1)環(huán)環(huán)電極結(jié)構(gòu)射流放電電流表現(xiàn)為多脈沖的形式,高低壓電極間的放電表現(xiàn)為DBD放電形式,除了噴出管外的射流,玻璃管內(nèi)高壓電極向內(nèi)側(cè)會(huì)出現(xiàn)一定的溢流放電。

    2)建立了一種能反映環(huán)環(huán)電極結(jié)構(gòu)射流放電過程的等效電氣模型,高低壓電極間的放電用基于VCCS電路等效,而射流和溢流部分用一個(gè)可變電阻和可變電容串聯(lián)組成的阻抗等效。基于此電氣模型的等效電路能準(zhǔn)確地反映大氣壓環(huán)環(huán)電極結(jié)構(gòu)射流放電的電氣特性。

    3)基于等效電路,利用Simulink建立的動(dòng)態(tài)仿真模型仿真得到的電壓、電流波形圖和Lissajous圖在電流脈沖個(gè)數(shù)、脈沖幅值以及脈沖持續(xù)時(shí)間上均與實(shí)驗(yàn)結(jié)果相符合。利用仿真模型可以進(jìn)一步研究電源頻率和環(huán)環(huán)間距對(duì)環(huán)環(huán)電極結(jié)構(gòu)射流放電電氣特性的影響情況,由仿真結(jié)果計(jì)算得到,放電功率和傳輸電荷都隨電源頻率的增加非線性增大,隨環(huán)環(huán)間距的增加非線性減小。故在實(shí)際應(yīng)用中,為了獲得較為強(qiáng)烈的射流放電,在進(jìn)行反應(yīng)器設(shè)計(jì)時(shí),應(yīng)采用較小的環(huán)環(huán)間距,同時(shí)在安全范圍內(nèi)盡可能提高電源頻率。

    4)利用此模型可以進(jìn)一步得到實(shí)際實(shí)驗(yàn)過程中無法直接測(cè)量獲得的放電參量,如Va、Vd、Vg和iccs等,從而對(duì)射流放電特性進(jìn)行更深入的研究,為優(yōu)化反應(yīng)器設(shè)計(jì)、提高放電效率提供參考。

    參考文獻(xiàn)

    [1] Deng X L, Nikiforov A Y, Vanraes P, et al. Direct current plasma jet at atmospheric pressure operating in nitrogen and air[J]. Current Applied Physics, 2013,113(2): 023305(1-9).

    [2] Shashurin A, Shneider M, Keidar M. Measurements of streamer head potential and conductivity of streamer column in cold nonequilibrium atmospheric plasmas[J]. Plasma Sources Science and Technology, 2012, 21(3): 1-6.

    [3] Wu S, Huang Q, Wang Z, et al. On the magnetic field signal radiated by an atmospheric pressure room temperature plasma jet[J]. Journal of Applied Physics, 2013, 113(4): 043305(1-6).

    [4] Xian Yubin, Wu Shuqun, Wang Zhan, et al. Discharge dynamics and modes of an atmospheric pressure non-equilibrium air plasma jet[J]. Plasma Processes and Polymers, 2013, 10(4): 372-378.

    [5] Brian L, Shih K, Jared W, et al. Dynamic electric potential redistribution and its influence on the development of a dielectric barrier plasma jet[J]. Plasma Sources science and Technology, 2012, 21(3): 1-11.

    [6] Zhang Qian, Liang Yongdong, Feng Hongqing, et al. A study of oxidative stress induced by non-thermal plasma-activated water for bacterial damage[J]. Applied Physics Letter, 2013, 102(20): 203701(1-4).

    [7] David B. The emerging role of reactive oxygen and nitrogen species in redox biology and some implications for plasma applications to medicine and biology[J]. Journal of Physics D: Applied Physics, 2012, 45(26): 1-42.

    [8] Kakiuchi H, Higashida K, Shibata T, et al. High-rate HMDSO-based coatings in open air using atmosphericpressure plasma jet[J]. Journal of Non-Crystalline Solids, 2012, 358(17): 2462-2465.

    [9] Lee E, Kwon J, Song D, et al. The effects of non-thermal atmospheric pressure plasma jet on cellular activity at SLA-treated titanium surfaces[J]. Current Applied Physics, 2013, 13(S1): 36-41.

    [10] Yan Xu, Xiong Zilan, Zou Fei, et al. Plasma-induced etch of HepG2 cancer cells intracellular effects of reactive species[J]. Plasma Processes and Polymers, 2013, 9(1): 59-66.

    [11] Cordula M, Saskia M, Bienvenida G, et al. Impact of homogeneous and filamentary discharge modes on the efficiency of dielectric barrier discharge ionization mass spectrometry[J]. Analytical and Bioanalytical Chemistry, 2013, 405(14): 4729-4735.

    [12] Yonemori S, Nakagawa Y, One R, et al. Measurement of OH density and air–helium mixture ratio in an atmospheric-pressure helium plasma jet[J]. Journal of Physics D: Applied Physics, 2012, 45(22): 1-8.

    [13] Oh J S, Walsh J L, Bradley J W. Plasma bullet current measurements in a free-stream helium capillary jet[J]. Plasma Sources Science and Technology, 2012, 21(3): 034020(1-6).

    [14] 方志, 劉源, 蔡玲玲. 大氣壓氬等離子體射流的放電特性[J]. 高電壓技術(shù), 2012, 38(7): 1613-1622.

    Fang Zhi, Liu Yuan, Cai Lingling. Discharge characteristics of atmosphere pressure plasma jet in Ar[J]. High Voltage Engineering, 2012, 38(7): 1613-1622.

    [15] 江南, 曹則賢. 一種大氣壓放電氦等離子體射流的實(shí)驗(yàn)研究[J]. 物理學(xué)報(bào), 2010, 59(5): 3324-3330.

    Jiang Nan, Cao Zexian. Experimental studies on an atmospheric pressure He plasma jet[J]. Asta Physics Sinica, 2010, 59(5): 3324-3330.

    [16] 張冠軍, 詹江楊, 邵先軍, 等. 大氣壓氬氣等離子體射流長(zhǎng)度的影響因素[J]. 高電壓技術(shù), 2011, 37(6): 1432-1438.

    Zhang Guanjun, Zhan Jiangyang, Shao Xianjun, et al. Influence factor analysis on jet length of atmospheric pressure argon plasma jets[J]. High Voltage Engineering, 2011, 37(6): 1432-1438.

    [17] Hao Z Y, Ji S C, Qiu A C. Study on the influence of dielectric barrier materials on the characteristics of atmospheric plasma jet in Ar[J]. IEEE Transactions on Plasma Science, 2012, 40(11): 2822-2830.

    [18] 陳田, 葉齊政, 譚丹, 等. 兩相體介質(zhì)阻擋放電中的三種放電形式[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2012, 32(22): 182-187.

    Chen Tian, Ye Qizheng, Tan Dan, et al. Dielectric barrier discharge in a two-phase mixture[J]. Proceedings of the CSEE, 2012, 32(22): 182-187.

    [19] 李清泉, 許光可, 房新振, 等. 沿面型介質(zhì)阻擋放電的數(shù)值仿真計(jì)算[J]. 高電壓技術(shù), 2012, 38(7): 1548-1555.

    Li Qingquan, Xu Guangke, Fang Xinzhen, et al. Numerical simulation of surface dielectric barrier discharge[J]. High Voltage Engineering, 2012, 38(7): 1548-1555.

    [20] 章程, 方志, 趙龍章, 等. 介質(zhì)阻擋放電特性的仿真與實(shí)驗(yàn)研究[J]. 真空科學(xué)與技術(shù)學(xué)報(bào), 2008, 28(1): 32-36.

    Zhang Cheng, Fang Zhi, Zhao Longzhang, et al. Simulated and experimental studies of discharge characteristics in dielectric barrier discharge[J]. Chinese Journal of Vacuum Science and Technology, 2008, 28(1): 32-36.

    [21] 章程, 方志, 胡建杭, 等. 不同條件下介質(zhì)阻擋放電的仿真與實(shí)驗(yàn)研究[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2008, 28(34): 33-39.

    Zhang Cheng, Fang Zhi, Hu Jianhang, et al. Simulation and experiment study on dielectric barrier discharge under different conditions[J]. Proceedings of the CSEE, 2008, 28(34): 33-39.

    [22] 潘俊, 方志. 多脈沖均勻介質(zhì)阻擋放電特性的仿真及實(shí)驗(yàn)研究[J]. 高電壓技術(shù), 2012, 38(5): 1132-1140.

    Pan Jun, Fang Zhi. Simulation and experimental studies on discharge characteristics of multiple pulse homogeneous dielectric barrier discharge[J]. High Voltage Engineering, 2012, 38(5): 1132-1140.

    [23] Fang Z, Ji S, Pan J, et al. Electrical model and experimental analysis of the atmospheric pressure homogeneous dielectric barrier discharge in He[J]. IEEE Transactions on Plasma Science, 2012, 40(3): 883-891.

    [24] 郝艷捧, 王曉蕾, 陽(yáng)林. 大氣壓氦氣介質(zhì)阻擋放電多脈沖輝光放電的形成條件[J]. 電工技術(shù)學(xué)報(bào), 2009, 24(9): 28-32.

    Hao Yanpeng, Wang Xiaolei, Yang Lin. Formation of dielectric barrier multi-pulse glow discharges in helium at atmospheric pressure[J]. Transactions of China Electrotechnical Society, 2009, 24(9): 28-32.

    [25] 周亦驍, 方志, 邵濤. 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.

    [26] 郝艷捧, 涂恩來, 陽(yáng)林, 等. 基于氣隙伏安特性研究大氣壓氦氣輝光放電的模式和機(jī)理[J]. 電工技術(shù)學(xué)報(bào), 2010, 25(7): 24-30.

    Hao Yanpeng, Tu Enlai, Yang Lin, et al. Mode and mechanism of multi-pulse glow discharges in helium at atmospheric pressure based on voltage-current characteristics[J]. Transactions of China Electrotechnical Society, 2010, 25(7): 24-30.

    [27] 章程, 邵濤, 龍凱華, 等. 大氣壓空氣中納秒脈沖介質(zhì)阻擋放電均勻性的研究[J]. 電工技術(shù)學(xué)報(bào), 2010, 25(1): 30-36.

    Zhang Cheng, Shao Tao, Long Kaihua, et al. Uniform of unipolar nanosecond pulse DBD in atmospheric air[J]. Transactions of China Electrotechnical Society, 2010, 25(1): 30-36.

    [28] 章程, 邵濤, 于洋, 等. 納秒脈沖介質(zhì)阻擋放電特性及其聚合物材料表面改性[J]. 電工技術(shù)學(xué)報(bào), 2010, 25(5): 31-37.

    Zhang Cheng, Shao Tao, Yu Yang, et al. Characteristics of unipolar nanosecond pulse DBD and its application on surface treatment of polymer films[J]. Transactions of China Electrotechnical Society, 2010, 25(5): 31-37.

    方 志 男,1975年生,教授,碩士生導(dǎo)師,從事高電壓絕緣及氣體放電等離子體的基本理論及其用于材料表面改性研究工作。E-mail: myfz@263.net(通信作者)

    錢 晨 男,1989年生,碩士研究生,研究方向?yàn)榈入x子體射流的實(shí)驗(yàn)與仿真研究。

    E-mail: money19891008@126.com

    作者簡(jiǎn)介

    收稿日期2014-01-14 改稿日期 2014-07-05

    中圖分類號(hào):TM213

    午夜视频国产福利| 国产乱人视频| 亚洲精品色激情综合| 中文字幕熟女人妻在线| 简卡轻食公司| 又黄又爽又刺激的免费视频.| 中文资源天堂在线| 日本五十路高清| 午夜福利在线观看免费完整高清在| 日本-黄色视频高清免费观看| 国产v大片淫在线免费观看| 男人的好看免费观看在线视频| 中文乱码字字幕精品一区二区三区 | 国产乱人偷精品视频| 亚洲天堂国产精品一区在线| 亚洲人成网站在线观看播放| 国产白丝娇喘喷水9色精品| 黄色日韩在线| 在线观看av片永久免费下载| 欧美日本亚洲视频在线播放| 午夜a级毛片| 七月丁香在线播放| 欧美性猛交黑人性爽| 激情 狠狠 欧美| 久久99热这里只频精品6学生 | 秋霞在线观看毛片| 青春草国产在线视频| 春色校园在线视频观看| 成人二区视频| 秋霞伦理黄片| 亚洲av二区三区四区| 欧美日韩一区二区视频在线观看视频在线 | av播播在线观看一区| 久久久久久久久久久免费av| 九九在线视频观看精品| 亚洲激情五月婷婷啪啪| 色噜噜av男人的天堂激情| 久久精品夜夜夜夜夜久久蜜豆| 日本wwww免费看| 亚州av有码| 国产高清国产精品国产三级 | 大话2 男鬼变身卡| 欧美bdsm另类| 色综合站精品国产| 久久久久九九精品影院| 美女高潮的动态| 内地一区二区视频在线| 国产午夜精品一二区理论片| 日本黄色视频三级网站网址| 1000部很黄的大片| 亚洲欧美成人综合另类久久久 | 老司机影院成人| 中文字幕精品亚洲无线码一区| 久热久热在线精品观看| 亚洲av中文字字幕乱码综合| av福利片在线观看| 亚洲乱码一区二区免费版| 99热网站在线观看| 欧美一区二区亚洲| 欧美色视频一区免费| 色视频www国产| 国语对白做爰xxxⅹ性视频网站| 国产精品一区www在线观看| 久久久国产成人精品二区| 国产在线一区二区三区精 | 小说图片视频综合网站| 国产亚洲5aaaaa淫片| 久久久久久伊人网av| 成年女人永久免费观看视频| 毛片一级片免费看久久久久| 久久久久性生活片| 久久久精品大字幕| 看十八女毛片水多多多| 日产精品乱码卡一卡2卡三| 久久这里只有精品中国| 成人漫画全彩无遮挡| 天天躁日日操中文字幕| 欧美丝袜亚洲另类| 男女那种视频在线观看| 一个人看视频在线观看www免费| 最近中文字幕2019免费版| 国产三级中文精品| 成人三级黄色视频| 一级二级三级毛片免费看| 亚洲自偷自拍三级| 一个人观看的视频www高清免费观看| 极品教师在线视频| 成人亚洲欧美一区二区av| 久久久成人免费电影| av卡一久久| 级片在线观看| 熟女电影av网| 九九热线精品视视频播放| 熟女人妻精品中文字幕| 国产伦在线观看视频一区| 国产精品无大码| 免费观看人在逋| 又爽又黄无遮挡网站| 伦精品一区二区三区| 欧美日韩在线观看h| 我要搜黄色片| 成人亚洲精品av一区二区| 色综合色国产| 日韩三级伦理在线观看| 国产黄色视频一区二区在线观看 | 色综合站精品国产| 久久精品国产99精品国产亚洲性色| 亚洲在久久综合| 亚洲精品成人久久久久久| 亚洲欧美日韩卡通动漫| 夜夜爽夜夜爽视频| 只有这里有精品99| 天天躁日日操中文字幕| 美女被艹到高潮喷水动态| 久久久亚洲精品成人影院| 青春草亚洲视频在线观看| 国产精品久久久久久久久免| 精品不卡国产一区二区三区| 国产精品三级大全| 免费看日本二区| 国产精品一区二区三区四区免费观看| 国产男人的电影天堂91| 成年免费大片在线观看| 国产高清不卡午夜福利| 男女边吃奶边做爰视频| 99热6这里只有精品| 久久久久久久久久成人| 免费不卡的大黄色大毛片视频在线观看 | 热99re8久久精品国产| 久久亚洲国产成人精品v| 日本与韩国留学比较| 日本av手机在线免费观看| 久久久久久久国产电影| 纵有疾风起免费观看全集完整版 | 我要看日韩黄色一级片| 国产久久久一区二区三区| 插逼视频在线观看| 日韩强制内射视频| 国产乱人偷精品视频| 国产精品一区二区性色av| 精品久久国产蜜桃| av女优亚洲男人天堂| 久久久久久久久中文| 亚洲成人精品中文字幕电影| 丝袜美腿在线中文| 国产亚洲av片在线观看秒播厂 | 亚洲欧美精品自产自拍| 久久久久久大精品| 国产乱人偷精品视频| 国产一级毛片七仙女欲春2| 国产欧美另类精品又又久久亚洲欧美| 99热这里只有是精品在线观看| 国产精品乱码一区二三区的特点| 亚洲av中文av极速乱| 亚洲国产欧美人成| 亚洲四区av| 国产免费又黄又爽又色| 日本色播在线视频| 久久久久九九精品影院| 赤兔流量卡办理| 成人三级黄色视频| 99在线视频只有这里精品首页| 中文在线观看免费www的网站| 欧美一区二区精品小视频在线| 久久久精品欧美日韩精品| 成年av动漫网址| av.在线天堂| 久久精品影院6| 亚洲国产最新在线播放| 久久久久性生活片| 国产精品福利在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 91在线精品国自产拍蜜月| 男人的好看免费观看在线视频| 久久精品国产99精品国产亚洲性色| 久久久a久久爽久久v久久| 老司机福利观看| a级一级毛片免费在线观看| 插逼视频在线观看| 午夜福利在线在线| 中文字幕av成人在线电影| 黄片wwwwww| 国产午夜精品论理片| 91精品一卡2卡3卡4卡| 久久精品国产亚洲av涩爱| 久久鲁丝午夜福利片| 免费av观看视频| 内射极品少妇av片p| 欧美日韩综合久久久久久| 永久网站在线| 国产三级中文精品| 亚洲成人中文字幕在线播放| 亚洲欧美成人综合另类久久久 | 欧美zozozo另类| 国产v大片淫在线免费观看| 久久久久免费精品人妻一区二区| 亚洲欧美精品专区久久| 久久韩国三级中文字幕| 色视频www国产| 免费一级毛片在线播放高清视频| 午夜福利网站1000一区二区三区| 亚洲欧洲日产国产| 国产高清有码在线观看视频| 美女被艹到高潮喷水动态| 午夜福利高清视频| 欧美精品国产亚洲| 精品久久久久久成人av| 三级男女做爰猛烈吃奶摸视频| 国产精品av视频在线免费观看| 欧美激情在线99| 女的被弄到高潮叫床怎么办| 国产精品1区2区在线观看.| 最新中文字幕久久久久| 国产精品美女特级片免费视频播放器| 99久久精品一区二区三区| 午夜福利高清视频| 有码 亚洲区| 亚洲国产精品成人久久小说| 波多野结衣巨乳人妻| 亚洲怡红院男人天堂| 亚洲va在线va天堂va国产| 成人美女网站在线观看视频| 夜夜看夜夜爽夜夜摸| 人体艺术视频欧美日本| 亚洲国产最新在线播放| 欧美成人精品欧美一级黄| 久久久久久久久中文| 老司机影院毛片| 最近最新中文字幕免费大全7| 视频中文字幕在线观看| 卡戴珊不雅视频在线播放| 中文字幕免费在线视频6| 青春草视频在线免费观看| h日本视频在线播放| 日日啪夜夜撸| 亚洲自拍偷在线| 精华霜和精华液先用哪个| 神马国产精品三级电影在线观看| 亚洲成人中文字幕在线播放| 高清av免费在线| 最近中文字幕高清免费大全6| 女的被弄到高潮叫床怎么办| 亚洲人成网站在线观看播放| av天堂中文字幕网| 亚洲精品日韩av片在线观看| 国产高清视频在线观看网站| 日本av手机在线免费观看| 99久久精品国产国产毛片| 特大巨黑吊av在线直播| 久久久精品大字幕| 亚洲成色77777| 一区二区三区高清视频在线| 日韩制服骚丝袜av| 一级二级三级毛片免费看| 欧美成人a在线观看| 校园人妻丝袜中文字幕| 乱系列少妇在线播放| 精品欧美国产一区二区三| 国产老妇伦熟女老妇高清| 国产综合懂色| 亚洲精品久久久久久婷婷小说 | 国产老妇女一区| 久久久久久大精品| 国模一区二区三区四区视频| 国产在线男女| 亚洲精品aⅴ在线观看| 久久久久久国产a免费观看| 午夜精品在线福利| 亚洲人成网站在线观看播放| 久久久亚洲精品成人影院| 亚洲va在线va天堂va国产| 国产精品电影一区二区三区| 中文欧美无线码| 嫩草影院入口| 亚洲丝袜综合中文字幕| 免费看a级黄色片| 亚洲av免费高清在线观看| 日韩一区二区视频免费看| 汤姆久久久久久久影院中文字幕 | 精华霜和精华液先用哪个| 一区二区三区四区激情视频| 国产高清三级在线| 亚洲国产精品专区欧美| 亚洲欧美中文字幕日韩二区| 国产成人一区二区在线| 亚洲av男天堂| 久99久视频精品免费| 欧美成人午夜免费资源| 夫妻性生交免费视频一级片| ponron亚洲| 日本猛色少妇xxxxx猛交久久| 久久精品夜色国产| 久久久久国产网址| av播播在线观看一区| 在线免费十八禁| 18+在线观看网站| 久久99精品国语久久久| 日韩av在线大香蕉| 精品国内亚洲2022精品成人| 91午夜精品亚洲一区二区三区| 免费观看人在逋| 成年av动漫网址| 国产午夜福利久久久久久| 国产黄a三级三级三级人| 国产精品国产三级国产专区5o | 六月丁香七月| 一级爰片在线观看| 日产精品乱码卡一卡2卡三| 99国产精品一区二区蜜桃av| 国产激情偷乱视频一区二区| 亚洲国产高清在线一区二区三| 中文精品一卡2卡3卡4更新| 亚洲成人久久爱视频| 美女被艹到高潮喷水动态| 色吧在线观看| 青春草国产在线视频| 欧美激情久久久久久爽电影| 在线免费观看的www视频| 国产老妇伦熟女老妇高清| 乱码一卡2卡4卡精品| 少妇丰满av| 亚洲av免费高清在线观看| 亚洲国产精品久久男人天堂| 日本熟妇午夜| 亚洲精品亚洲一区二区| 九九久久精品国产亚洲av麻豆| 日本猛色少妇xxxxx猛交久久| 综合色av麻豆| 有码 亚洲区| 国产一区二区在线av高清观看| 国产精品女同一区二区软件| 美女高潮的动态| 一级黄色大片毛片| 淫秽高清视频在线观看| 亚洲精品乱码久久久v下载方式| 色综合色国产| 色播亚洲综合网| 国产免费又黄又爽又色| 日韩欧美精品v在线| 一区二区三区乱码不卡18| 欧美变态另类bdsm刘玥| 免费av毛片视频| 麻豆成人午夜福利视频| 亚洲,欧美,日韩| 我要看日韩黄色一级片| 日韩欧美国产在线观看| 中文欧美无线码| 精品久久久久久久末码| 秋霞伦理黄片| 国产精品av视频在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 一区二区三区高清视频在线| 日韩成人av中文字幕在线观看| 在线播放国产精品三级| 久久久午夜欧美精品| 国产免费视频播放在线视频 | 国产单亲对白刺激| 亚洲中文字幕一区二区三区有码在线看| 国产精品人妻久久久影院| 免费搜索国产男女视频| a级毛色黄片| 亚洲av成人精品一区久久| 国产精品电影一区二区三区| 看黄色毛片网站| 国产精品综合久久久久久久免费| 亚洲精品乱久久久久久| 欧美激情在线99| 国产精品伦人一区二区| 日韩欧美在线乱码| 亚洲电影在线观看av| 婷婷色综合大香蕉| 免费大片18禁| 秋霞在线观看毛片| 国产在线男女| av播播在线观看一区| 一区二区三区免费毛片| 老女人水多毛片| 高清日韩中文字幕在线| 久久久国产成人精品二区| 熟女电影av网| 激情 狠狠 欧美| 亚洲人成网站在线播| 日本-黄色视频高清免费观看| 成人特级av手机在线观看| 亚洲精品456在线播放app| 在线观看av片永久免费下载| 嫩草影院入口| 亚洲欧美成人精品一区二区| 99久久中文字幕三级久久日本| 久久99热6这里只有精品| 日本wwww免费看| kizo精华| 少妇高潮的动态图| 日产精品乱码卡一卡2卡三| 欧美一区二区国产精品久久精品| 精品欧美国产一区二区三| 午夜爱爱视频在线播放| 禁无遮挡网站| 国产成人aa在线观看| 国产又黄又爽又无遮挡在线| 三级经典国产精品| 少妇人妻一区二区三区视频| av卡一久久| 人人妻人人看人人澡| 午夜a级毛片| 老司机影院成人| 伊人久久精品亚洲午夜| 国产国拍精品亚洲av在线观看| 亚洲国产成人一精品久久久| 欧美激情国产日韩精品一区| АⅤ资源中文在线天堂| 午夜激情欧美在线| 亚洲成色77777| 亚洲熟妇中文字幕五十中出| 熟女人妻精品中文字幕| 久久久欧美国产精品| 汤姆久久久久久久影院中文字幕 | 男女下面进入的视频免费午夜| 免费大片18禁| 国产精品美女特级片免费视频播放器| 91精品国产九色| 高清视频免费观看一区二区 | 免费电影在线观看免费观看| 最近中文字幕高清免费大全6| 国产精品一二三区在线看| 日韩在线高清观看一区二区三区| 日韩欧美精品v在线| 一区二区三区四区激情视频| 欧美丝袜亚洲另类| 小蜜桃在线观看免费完整版高清| 久久久亚洲精品成人影院| 国产精品日韩av在线免费观看| 欧美xxxx黑人xx丫x性爽| 成人漫画全彩无遮挡| a级毛色黄片| 亚洲电影在线观看av| .国产精品久久| 九九久久精品国产亚洲av麻豆| 亚洲综合精品二区| 国产欧美日韩精品一区二区| 黄片无遮挡物在线观看| 女人十人毛片免费观看3o分钟| 亚洲精品乱码久久久v下载方式| 欧美日韩一区二区视频在线观看视频在线 | 我要搜黄色片| av女优亚洲男人天堂| 成人二区视频| 伦理电影大哥的女人| 日本wwww免费看| 国语对白做爰xxxⅹ性视频网站| 国产精品一二三区在线看| 中文字幕av在线有码专区| 久久99蜜桃精品久久| 最近中文字幕高清免费大全6| av黄色大香蕉| 永久免费av网站大全| 精品久久久久久久久av| 一个人免费在线观看电影| 精品99又大又爽又粗少妇毛片| 国产黄片美女视频| 我的女老师完整版在线观看| 在现免费观看毛片| 青春草视频在线免费观看| 搡女人真爽免费视频火全软件| 精品久久久久久久末码| 熟妇人妻久久中文字幕3abv| 最近最新中文字幕免费大全7| 极品教师在线视频| 国产极品天堂在线| 国产v大片淫在线免费观看| 免费无遮挡裸体视频| 性色avwww在线观看| 岛国毛片在线播放| 日韩高清综合在线| 国产毛片a区久久久久| 午夜老司机福利剧场| 黄色一级大片看看| 亚洲最大成人av| av黄色大香蕉| 色吧在线观看| 老师上课跳d突然被开到最大视频| 自拍偷自拍亚洲精品老妇| 99久国产av精品国产电影| 99久久九九国产精品国产免费| or卡值多少钱| 亚洲在线观看片| 色综合站精品国产| 青春草视频在线免费观看| 成人国产麻豆网| 亚洲精品乱码久久久v下载方式| 日本午夜av视频| 韩国av在线不卡| 免费观看a级毛片全部| 久久99热6这里只有精品| 日韩av在线免费看完整版不卡| 99热全是精品| 午夜日本视频在线| 狂野欧美白嫩少妇大欣赏| 午夜日本视频在线| 国产真实乱freesex| 国产精品一区二区三区四区免费观看| 国产亚洲午夜精品一区二区久久 | 国产欧美另类精品又又久久亚洲欧美| 天堂影院成人在线观看| 亚洲成人中文字幕在线播放| 国产精品久久电影中文字幕| 日本色播在线视频| 美女高潮的动态| 亚洲欧洲国产日韩| 99热这里只有是精品在线观看| 国产爱豆传媒在线观看| 日韩欧美精品免费久久| 男的添女的下面高潮视频| 欧美色视频一区免费| 99视频精品全部免费 在线| 一本一本综合久久| 午夜福利在线观看免费完整高清在| 97热精品久久久久久| 国产精品一区二区在线观看99 | 国产在视频线在精品| 国内揄拍国产精品人妻在线| 国产69精品久久久久777片| 中文字幕av成人在线电影| 最近2019中文字幕mv第一页| 91午夜精品亚洲一区二区三区| 综合色丁香网| 精品人妻一区二区三区麻豆| 国产精品无大码| 婷婷色av中文字幕| 国产精品无大码| 啦啦啦啦在线视频资源| 亚洲av日韩在线播放| 99久久精品国产国产毛片| 日本一本二区三区精品| 免费不卡的大黄色大毛片视频在线观看 | 三级国产精品片| 国产三级在线视频| 欧美丝袜亚洲另类| 国产精品福利在线免费观看| 成年免费大片在线观看| 久久热精品热| 亚洲精品影视一区二区三区av| 麻豆成人av视频| 久久精品综合一区二区三区| 天堂av国产一区二区熟女人妻| a级毛片免费高清观看在线播放| 麻豆一二三区av精品| 国产av不卡久久| 啦啦啦韩国在线观看视频| 综合色丁香网| 永久网站在线| 蜜桃亚洲精品一区二区三区| 国产精品一区二区三区四区久久| 国产伦在线观看视频一区| 国产在视频线精品| 日本与韩国留学比较| 尤物成人国产欧美一区二区三区| 欧美一区二区国产精品久久精品| av天堂中文字幕网| 韩国av在线不卡| 热99在线观看视频| 日韩一区二区视频免费看| 嫩草影院入口| 欧美又色又爽又黄视频| 九九爱精品视频在线观看| 亚洲国产欧美在线一区| 日韩,欧美,国产一区二区三区 | 边亲边吃奶的免费视频| 免费看av在线观看网站| 丰满人妻一区二区三区视频av| 国产精品电影一区二区三区| 日日干狠狠操夜夜爽| 身体一侧抽搐| 日韩亚洲欧美综合| 你懂的网址亚洲精品在线观看 | 国产麻豆成人av免费视频| 亚洲精品乱码久久久久久按摩| 色尼玛亚洲综合影院| 黑人高潮一二区| 精品国内亚洲2022精品成人| 综合色av麻豆| 色5月婷婷丁香| 午夜免费激情av| 寂寞人妻少妇视频99o| 国产人妻一区二区三区在| av在线观看视频网站免费| 国产单亲对白刺激| 天堂√8在线中文| 午夜视频国产福利| av在线天堂中文字幕| 国产伦在线观看视频一区| 久久精品久久精品一区二区三区| 国产高清不卡午夜福利| 一区二区三区免费毛片| 日本黄色视频三级网站网址| 1024手机看黄色片| 汤姆久久久久久久影院中文字幕 | 亚洲综合色惰| 91精品一卡2卡3卡4卡| 成人午夜高清在线视频| 国产老妇女一区| 国内精品宾馆在线| 91av网一区二区| 亚洲美女搞黄在线观看| 久久婷婷人人爽人人干人人爱| 真实男女啪啪啪动态图| 床上黄色一级片| 欧美最新免费一区二区三区| 精品欧美国产一区二区三| 亚洲精品久久久久久婷婷小说 | 日韩在线高清观看一区二区三区| 亚洲精品国产成人久久av| 特级一级黄色大片| 日本av手机在线免费观看| 永久免费av网站大全|