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

    孤立金屬尖端的電暈觸發(fā)閾值

    2014-07-06 10:54:38譚涌波朱俊儒李祥超梁忠武郭秀峰
    應(yīng)用氣象學(xué)報 2014年3期
    關(guān)鍵詞:電暈尖端電場

    譚涌波 朱俊儒 李祥超 梁忠武 郭秀峰

    1)(南京信息工程大學(xué)氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心,南京210044)2)(南京信息工程大學(xué)中國氣象局氣溶膠與云降水重點開放實驗室,南京210044)

    孤立金屬尖端的電暈觸發(fā)閾值

    譚涌波1)2)*朱俊儒1)2)李祥超1)2)梁忠武1)2)郭秀峰1)2)

    1)(南京信息工程大學(xué)氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心,南京210044)2)(南京信息工程大學(xué)中國氣象局氣溶膠與云降水重點開放實驗室,南京210044)

    為了得到金屬尖端在發(fā)生電暈放電時尖端處的電場強度,該文首先采用實驗室實驗得到不同高度、不同形狀、不同材質(zhì)的金屬尖端發(fā)生電暈放電時的環(huán)境電場閾值;再采用有限元法計算二維泊松方程,得到尖端處電暈觸發(fā)閾值,由此得出以下結(jié)論:環(huán)境電場閾值隨金屬尖端高度的增大基本呈線性減小趨勢,隨著尖端越來越尖,環(huán)境電場閾值呈先減小后增大的變化趨勢;高度、形狀對金屬物尖端處電暈觸發(fā)閾值無影響,尖端處電暈閾值為定值;給出尖端處電暈觸發(fā)閾值為158.75kV·m-1與空間分辨率的擬合公式,可為今后電暈放電數(shù)值模擬中判斷電暈放電的起始時刻提供參考。

    電暈放電;實驗室模擬;電暈觸發(fā)閾值;有限元;空間分辨率

    引 言

    高大建筑物在環(huán)境電場作用下,其尖端附近的電場會發(fā)生畸變[1-2];當(dāng)緊鄰建筑物尖端表面區(qū)域處的電場強度超過一定閾值時,便會發(fā)生電暈放電[34]。觀測發(fā)現(xiàn),電暈離子在電場力、風(fēng)場等作用下可以到達距離地面幾百米的高度處[5-7]。一般情況下,由于雷暴云云底的主電荷區(qū)為負(fù),所以電暈放電所釋放的離子通常帶正電。眾所周知,電暈放電對低層大氣電學(xué)的影響至關(guān)重要。首先,電暈離子聚集在尖端附近,會影響地表附近空間電荷層的分布,使尖端處電場強度減?。?-11],同時也改變了上行先導(dǎo)的觸發(fā)條件及上行先導(dǎo)向上傳輸條件。其次,電暈放電為對流層提供了豐富的離子源,這可能影響雷暴云起電過程及內(nèi)部電荷結(jié)構(gòu)[12-13]?;谝陨显?,電暈放電一直受到國內(nèi)外學(xué)者的關(guān)注[14-19]。到目前為止,對于電暈放電的研究主要集中于其形成的物理過程以及電參數(shù)相關(guān)特性方面。其中電暈放電的初始過程,尤其是放電閾值是此類研究中的重點和難點。

    目前關(guān)于電暈觸發(fā)閾值主要是通過野外觀測和實驗室實驗兩種手段獲得。野外觀測主要通過電暈電流測試儀上的電流與在其附近的同步電場儀上的電場,建立電暈電流與環(huán)境電場之間的擬合方程,并根據(jù)該方程求解電流為零時的環(huán)境電場值,即電暈放電的環(huán)境電場閾值[5,20-22]。而實驗室實驗是根據(jù)接地物內(nèi)產(chǎn)生的電暈電流達到人為規(guī)定的零電流時,在其附近的電場傳感器上顯示的電場值則被定為電暈觸發(fā)閾值[23-24]。然而,不論是野外觀測還是實驗室實驗,電暈觸發(fā)閾值的大小皆由電場儀所獲取的環(huán)境電場強度決定。由于電場儀測值大小會受其擺設(shè)位置、自身形狀以及金屬物自身的幾何尺寸等特征影響,目前已有的一些電暈初始閾值的計量結(jié)果其適用范圍都存在不確定性。因此,如何有效地給出電暈放電的觸發(fā)閾值仍是目前大氣電學(xué)研究中未解決的一個難題。

    基于此,本文通過實驗室實驗與理論計算相結(jié)合探討金屬物體尖端處觸發(fā)電暈放電的環(huán)境電場強度,擬采用尖端處電場大小作為判斷電暈放電的依據(jù),并分析其隨金屬物高度、形狀以及材質(zhì)的變化規(guī)律。

    1 實驗與計算方法

    首先,本文采用實驗室實驗得到不同類型的金屬尖端在發(fā)生電暈放電時上板施加的電壓值,即電暈電壓閾值;然后,將得到的電壓數(shù)據(jù)轉(zhuǎn)換成理論計算所需要的背景場;最后,在背景場下采用數(shù)值計算得到不同類型的尖端發(fā)生電暈放電時尖端處電場值。

    1.1 實驗方法

    實驗采用兩塊尺寸均為1.2m×1.2m×5mm金屬鐵板,在平板的上板施加一個負(fù)直流高壓,下板接地。實驗所采用的高壓設(shè)備是武漢大學(xué)恒安防雷接地技術(shù)中線研制的沖擊控制系統(tǒng)(ICGS),最高沖擊電壓可達到-60kV。在實驗中兩板間距離設(shè)為1m。尖端安置在下板正中間,其中尖端與下板之間絕緣。用于測量電暈電流的裝置由1個1MΩ電阻(在1μA的電暈電流通過1MΩ電阻時,電阻兩端的電壓為1V)與10MΩ電阻構(gòu)成反相運算放大電路。尖端放電電流通過1MΩ電阻后,將電流信號轉(zhuǎn)化成電壓信號,為了將更小的電暈電流信號被放大,將1MΩ電阻兩端的信號通過反相運算放大電路,將信號放大10倍。

    實驗采用的金屬尖端分為兩大類:①圓柱型,直徑D為2.0cm,高度Ha分別為25,30,39,50cm的鐵制材料尖端及高度Ha為30cm的鋁和銅材料尖端;②錐、柱型,尖端上部為圓錐形,下部分為圓柱形(直徑D=2.0cm),金屬尖端的總高度Ha均為30cm,材料均為鐵制,上端圓錐高度Hb分別為2.3,3.0,4.5,6.0cm (圖1)。

    圖1 實驗采用的不同高度的金屬尖端(a)及不同形狀的金屬尖端(b)Fig.1 Experiment with different heights of the lightning rod(a)and different shapes of lightning rod(b)

    1.2 計算方法

    由于金屬物尖端處的電場無法通過儀器測量得到,所以要得到尖端處的電場強度只能通過數(shù)值計算求得。本實驗測量的電壓數(shù)據(jù)為兩平行板之間的電位差,認(rèn)為兩板之間為均勻電場,在均勻電場背景下計算尖端處電暈觸發(fā)閾值,即本文分兩部分計算尖端處電暈電場閾值。

    1.2.1 均勻電場的計算

    由于實驗中上、下板邊緣部分的邊緣效應(yīng),使兩板之間的電場并不是均勻分布的;兩板之間的均勻電場要小于簡單計算得到的兩板間電壓U與距離d的比值,因此假設(shè)兩板間的均勻電場E0=αU/d,其中α為修正系數(shù),文章所選取的修正系數(shù)α=0.91[25]。

    1.2.2 電暈觸發(fā)閾值的計算

    在設(shè)定的均勻背景電場條件下求解模擬區(qū)域中的泊松方程,由于該區(qū)域內(nèi)無自由電荷存在,則泊松方程轉(zhuǎn)換為拉普拉斯方程。本文采用二維有限元法求解拉普拉斯方程最終計算出尖端處電場值[26]。本文理論計算設(shè)定的區(qū)域為1.2m×1.2m,將該區(qū)域劃分成若干個邊長為l互不重疊的四邊形有限單元,對每個四邊形單元進行線性插值得到四邊形單元及各頂點上的電位插值函數(shù),在線性四邊形單元基礎(chǔ)上采用里茲變分法導(dǎo)出電位方程組,采用迭代法計算出每個四邊形單元上的電位φ,最終根據(jù)E=-φ計算各節(jié)點上的電場[26-27]。對 于 邊 界 條件[28],下板、尖端與上板電位滿足第1類邊界條件(即狄里赫利邊界條件)為常數(shù),即上板電位為αU,下板與尖端的電位為零(下板與尖端是等電位)。模擬域的側(cè)邊界采用第3類邊界條件(即齊諾伊曼邊界條件),使用條件是物理量的法向?qū)?shù)在邊界值為零。

    2 實驗結(jié)果及計算結(jié)果

    2.1 實驗結(jié)果

    當(dāng)平行板上板施加的電壓致使金屬尖端發(fā)生尖端放電時,示波器會顯示出周期性的電暈放電脈沖信號。由于上板施加的電壓為負(fù),針上產(chǎn)生的電暈電流方向向上,因此反相運算放大電路的輸入電壓為負(fù),則最終的輸出電壓為正,即示波器上顯示的信號為正。對于電暈放電脈沖結(jié)束后伴隨著小的負(fù)脈沖出現(xiàn),本文認(rèn)為是由于電暈放電產(chǎn)生的正離子漂浮在金屬物體尖端處使尖端處電場極性發(fā)生反轉(zhuǎn)所致:假設(shè)電暈放電產(chǎn)生的離子為小離子且在針尖端附近僅受電場力作用,速度V=μE;其中,μ是小離子的遷移率,尖端處電場采用有限元法計算得到E=3.5×105V·m-1(分辨率f=0.2mm);最終得出電暈放電產(chǎn)生的小離子在正脈沖的持續(xù)時間內(nèi)所移動的最大距離d=μEt=2×10-4×3.5×105×4×10-5=2.8mm。由于運動距離為毫米量級,所以在電暈放電脈沖發(fā)生時,電暈放電產(chǎn)生的離子在尖端表層可以使尖端處的電場極性發(fā)生反轉(zhuǎn)。

    當(dāng)高壓設(shè)備施加在上板電壓上升到起先設(shè)定的電壓值后便會以0.1kV幅度下降,因此示波器上顯示出的電暈放電脈沖幅值會隨之減小,直至消失。由于存在干擾,無法根據(jù)示波器精確得到電暈脈沖起始觸發(fā)時刻。在本實驗中示波器上干擾信號的最大幅值為12mV,但12mV高頻脈沖信號所占比例甚少,所以實驗選取V=15mV作為判斷電暈放電的起始時刻。最終得出同一形狀、同一材質(zhì)、不同高度以及同一高度、同一材質(zhì)、不同形狀的尖端發(fā)生電暈放電時上板的電壓值,見表1;同一形狀、同一高度(Ha=30cm)、不同材質(zhì)的尖端發(fā)生電暈放電時上板的電壓值,見表2。

    表1 不同高度、不同形狀的鐵質(zhì)尖端發(fā)生電暈放電時上板電暈電壓值(單位:kV)Table 1 Corona voltage value of iron material rods with different heights and different shapes(unit:kV)

    表2 同一形狀、同一高度、不同材質(zhì)的尖端上板電暈電壓值(單位:kV)Table 2 Corona voltage value of lightning rod with the same shape,the same height and different materials(unit:kV)

    由實驗數(shù)據(jù)發(fā)現(xiàn),對于同為圓柱型鐵制材料的尖端在發(fā)生電暈放電時所需要的電暈電壓值隨高度的增長而減小;金屬尖端形狀變化對所需要觸發(fā)電暈放電的電壓影響較為強烈,隨著金屬物體上端越來越尖(即金屬物體尖端處圓錐高Hb與寬度D之比越大)所要觸發(fā)電暈放電的電壓值就越高,即表1。

    材料對電暈放電有微弱影響,鐵、鋁、銅電導(dǎo)率分別為0.1728,0.6321,1.0000,電導(dǎo)率相差比較明顯,而實驗中這3種材質(zhì)的尖端(相同高度、相同形狀)觸發(fā)電暈放電時所需要的電暈電壓只有輕微的減小(表2)。

    將實驗得到的鐵制尖端發(fā)生電暈放電時所需要電暈電壓值分別轉(zhuǎn)化成環(huán)境電場閾值,圖2a是鐵制尖端的環(huán)境電場閾值與高度之間的變化關(guān)系圖,圖2b為鐵制尖端的環(huán)境電場閾值與其形狀之間的變化規(guī)律。由圖2a可以看出,環(huán)境電場閾值隨高度增加基本上呈線性減小;由圖2b可以看出,實驗中所選尖端的尖端系數(shù)(Hb/D)分別為0,1.15,1.5,2.25,3,環(huán)境電場閾值分別為-21.13,-20.32,-21.28,-21.97,-22.95kV·m-1,電暈環(huán)境電場閾值隨尖端變得越來越尖出現(xiàn)先減小再增大的變化趨勢,這種變化與以往其他研究者得出的結(jié)論較為一致[29-31]。

    圖2 鐵制尖端觸發(fā)電暈放電時環(huán)境電場強度隨高度變化(a)以及Hb/D變化(b)Fig.2 The change of environmental field strength when iron lightning rod trigger corona discharge with heights(a)and withHb/D(b)

    2.2 計算結(jié)果

    在數(shù)值計算中,離散網(wǎng)格帶來的空間分辨率會給計算結(jié)果帶來一定的系統(tǒng)誤差,一般而言,分辨率越高,計算出來的結(jié)果就越接近連續(xù)空間中的真實值,系統(tǒng)誤差也就越??;分辨率越低,系統(tǒng)誤差就越大,計算就越不準(zhǔn)確。本文根據(jù)實驗得到的不同類型下的尖端在發(fā)生電暈放電時上板施加的電壓數(shù)據(jù),將其轉(zhuǎn)化成均勻環(huán)境電場時上板上電壓值,采用有限元理論在分辨率f=1cm情況下,計算出金屬尖端發(fā)生電暈放電時尖端處的電場閾值。表3為同一形狀、同一材質(zhì)、不同高度及同一高度、同一材質(zhì)、不同形狀下的尖端觸發(fā)電暈放電時尖端處的電場值。

    表3 不同高度、不同形狀的鐵質(zhì)尖端發(fā)生電暈放電時尖端處觸發(fā)閾值(單位:kV·m-1)Table 3 Corona field threshold of iron material rods with different heights and different shapes(unit:kV·m-1)

    由表3可以發(fā)現(xiàn),不論改變金屬尖端的高度,還是改變其形狀,最終通過理論計算得出的金屬尖端在發(fā)生電暈放電時尖端處電場分布都較為集中,即呈現(xiàn)鋸齒狀分布。所以本文認(rèn)為鐵質(zhì)材料的尖端不管是何種形狀、高度,在發(fā)生電暈放電時尖端處的電場是個定值,即當(dāng)鐵質(zhì)材料的金屬物體尖端處電場達到該值,尖端就會發(fā)生電暈放電。因此,對8種尖端發(fā)生電暈放電時尖端處電場閾值求平均值,即Ec=135.181kV·m-1(f=1cm)。

    2.3 結(jié)果分析

    表4 不同形狀、不同高度鐵制金屬物體尖端電暈電場閾值的相對誤差(單位:%)Table 4 The relative error of corona field threshold of iron metal objects with different shapes and different heights(unit:%)

    同時,采用有限元理論計算了材質(zhì)為鐵、鋁、銅尖端(Ha=30cm)發(fā)生電暈放電時尖端處電場值,見表5。若鋁、銅材質(zhì)的尖端選取了與鐵制材料一樣的尖端處電暈電場閾值Ec=135.181kV·m-1時,由此計算出每次實驗得到尖端處電暈電場閾值的相對誤差見表6。表6中尖端處電暈電場閾值的最大相對誤差為5.435%。由于最大誤差值小于10%,所以本文認(rèn)為對于銅、鋁材質(zhì)的尖端發(fā)生電暈放電時尖端處電暈電場值可以近似為鐵制材料的尖端發(fā)生電暈放電時尖端處電暈電場閾值,即Ec。

    表5 同一形狀、同一高度、不同材質(zhì)的尖端發(fā)生電暈放電時尖端處觸發(fā)值(單位:kV·m-1)Table 5 Corona field threshold of lightning rod with the same shape and the same height and different materials(unit:kV·m-1)

    表6 不同材質(zhì)的尖端電暈電場閾值的相對誤差(單位:%)Table 6 The relative error of corona field threshold of metal objects with different materials(unit:%)

    3 不同空間分辨率的電暈觸發(fā)閾值

    在采用數(shù)值計算方法計算金屬尖端在尖端處電暈電場閾值過程中,由于數(shù)值計算需要將連續(xù)的空間劃分為許多離散網(wǎng)格,在離散化的區(qū)域中計算得到每個單元上的電場值(其值表示為單元內(nèi)所有連續(xù)點上電場的平均值),這就會給計算結(jié)果帶來一定的誤差。此誤差值的大小僅與分辨率有關(guān),與尖端的尺寸大小無關(guān)。分辨率越低(網(wǎng)格越粗),計算結(jié)果的系統(tǒng)誤差值越大;分辨率越高(網(wǎng)格越細(xì)),計算結(jié)果的系統(tǒng)誤差越小。因此,在計算金屬物發(fā)生電暈電暈放電時尖端處電場閾值時,分辨率的選取對于計算尖端處電場強度的影響比較大。本文計算了在分辨率f分別為0.01,0.02,0.03,0.04,0.06,0.08,0.1,0.15,0.2,0.3,0.4m 情況下尖端處電暈電場閾值,計算結(jié)果分別為135.18,102.37,88.42,79,67.38,61.26,54.83,49.15,45.87,42.69,40.78kV·m-1(圖3)。由圖3可以看出,金屬體尖端處電暈電場閾值隨分辨率降低總體呈遞減趨勢,分辨率為0~0.1m時,隨分辨率降低,尖端電場閾值減小的比較明顯;當(dāng)分辨率f>0.1m時,隨分辨率降低尖端電場閾值減小的較為緩慢。

    圖3 不同空間分辨率的尖端處電暈觸發(fā)閾值Fig.3 The corona trigger threshold at the tip of metal cusps with different spatial resolutions

    計算不同分辨率尖端電暈電場閾值發(fā)現(xiàn),尖端處電暈電場閾值Ec和空間分辨率f之間呈負(fù)指數(shù)函數(shù)關(guān)系,則不同分辨率電暈電場閾值Ec(f)的擬合函數(shù)形式為

    式(1)中,常數(shù)a,b,c均大于0。圖3中顯示的黑線為尖端處電暈觸發(fā)閾值與分辨率f之間的擬合函數(shù)曲線,其擬合函數(shù)為

    擬合方程的決定系數(shù)R2(R是相關(guān)系數(shù),其值在0~1變化,若R2接近1時表明擬合效果好皆為1)為0.977,說明相關(guān)程度較高(達到0.05顯著性水平)。通過采用尖端處電暈觸發(fā)閾值作為判斷電暈放電的依據(jù),比用環(huán)境電場閾值作為判斷電暈放電的依據(jù)更可靠,主要是由于尖端處電場不受環(huán)境、自身形狀等不確定因素的影響。文中給出的尖端電暈觸發(fā)閾值與分辨率的擬合關(guān)系,可為今后電暈放電數(shù)值模擬研究中在判斷電暈放電的起始時刻提供參考。

    由于尖端頂部的實際電暈觸發(fā)閾值無法通過現(xiàn)有觀測手段獲得,而現(xiàn)有的網(wǎng)格化計算方法無法直接給出比擬實際連續(xù)空間的理論真值。因此本文通過擬合得到的尖端處電暈觸發(fā)閾值與網(wǎng)格間距之間的關(guān)系方程進行外推求極限:

    由式(3)可得,在連續(xù)空間中尖端發(fā)生電暈放電時尖端處電暈觸發(fā)閾值為158.75kV·m-1。

    4 結(jié)論與討論

    通過實驗得到同一材質(zhì)、同一形狀、不同高度,同一高度、同一形狀、不同材質(zhì),同一材質(zhì)、同一高度、不同形狀3種類型金屬尖端發(fā)生電暈時的電壓數(shù)據(jù),并采用有限元法計算金屬物發(fā)生電暈放電時尖端電場,可以得出以下結(jié)論:

    1)尖端尺寸對電暈環(huán)境電場閾值的影響隨尖端高度增高、尖端觸發(fā)電暈放電所需要的電暈環(huán)境電場閾值而減小,即尖端的高度越高,越容易觸發(fā)電暈放電;隨金屬物尖端形狀越來越尖,電暈環(huán)境電場閾值呈先減小再增大的變化趨勢。

    2)尖端形狀對尖端處電暈觸發(fā)閾值沒有影響,不管是改變尖端的高度還是改變其形狀,計算得到尖端觸發(fā)閾值為定值;對于銅、鋁、鐵材質(zhì)的金屬物在尖端發(fā)生電暈放電時尖端處的電場相差較小。因此,將這3種材質(zhì)的尖端在發(fā)生電暈放電時尖端處觸發(fā)閾值近似為158.75kV·m-1。

    3)空間分辨率變化影響最終計算尖端處電暈觸發(fā)閾值結(jié)果,最終給出尖端處電暈觸發(fā)閾值與分辨率f關(guān)系:Ec(f)=113.14e-f/0.0337+45.6。

    通過計算發(fā)現(xiàn),材質(zhì)選取對尖端處電暈電場閾值會產(chǎn)生影響,因此對于材料的選取將是以后實驗工作中所要探討的一部分。同時,實驗中濕度、氣壓、溫度等環(huán)境因素均為定值,因此改變實驗室的濕度、氣壓、溫度等環(huán)境因素,探討其對尖端處電暈觸發(fā)閾值的影響,也是下一步研究的重點。

    [1] 王道洪,郄秀書,郭昌明.雷電與人工引雷.上海:上海交通大學(xué)出版社,2000.

    [2] 吳亭,呂偉濤,劉曉陽,等.北京地區(qū)不同天氣條件下近地面大氣電場特征.應(yīng)用氣象學(xué)報,2009,20(4):394-401.

    [3] 張義軍,周秀驥.雷電研究的回顧和進展.應(yīng)用氣象學(xué)報,2006,17(6):829-834.

    [4] 張義軍,孟青,馬明,等.閃電探測技術(shù)發(fā)展和資料應(yīng)用.應(yīng)用氣象學(xué)報,2006,17(5):611-620.

    [5] Standler R B,Winn W P.Effects of coronae on electric fields beneath thunderstorms.Quart J Roy Meteor Soc,1979,105(443):285-302.

    [6] Chauzy S,Raizonville P.Space charge layers created by coronae at ground level below thunderclouds:Measurements and modeling.Journal of Geophysical Research:Oceans,1982,87(C4):3143-3148.

    [7] Chauzy S,Soula S.General interpretation of surface electric field variations between lightning flashes.Journal of Geophysical Research:Atmospheres,1987,92(D5):5676-5684.

    [8] Whipple F J W.Modern views on atmospheric electricity.Quarterly Journal of the Royal Meteorological Society,1938,64(275):199-222.

    [9] Toland R B,Vonnegut B.Measurement of maximum electric field intensities over water during thunderstorms.J Geophys Res,1977,82(3):438-440.

    [10] Whipple F J W,Scrase F J.Point discharge in the electric field of the earth.Geophysical Memoirs,H M S Meteorol Office,1936,68(7):1-20.

    [11] Chauzy S,Medale J C,Prieur S,et al.Multilevel measurement of the electric field underneath a thundercloud:1.A new system and the associated data processing.Journal of Geophysical Research:Atmospheres,1991,96(22):319-326.

    [12] Vonnegut B.Possible mechanism for the formation of thunderstorm electricity.Bull Amer Meteor Soc,1953,34:378-381.

    [13] Vonnegut B.Some facts and speculations concerning the origin and role of thunderstorm electricity.Meteorol Monogr,1963,5(27):224-241.

    [14] Chauzy S,Soula S,Despiau S.Ground coronae and lightning.Journal of Geophysical Research:Atmospheres,1989,94(D11):13115-13119.

    [15] 郄秀書.雷暴下地面自然尖端電暈放電離子時空演變的數(shù)值模擬.地球物理學(xué)報,1996,39(增刊):43-51.

    [16] Standler R B,Winn W P.Effects of coronae on electric fields beneath thunderstorms.Quart J Roy Meteor Soc,1979,105(443):285-302.

    [17] Soula S.Transfer of electrical space charge from corona between ground and thundercloud:Measurements and modeling.Journal of Geophysical Research:Atmospheres,1994,99(D5):10759-10765.

    [18] Aleksandrov N L,Bazelyan E M,Drabkin M M,et al.Corona discharge at the tip of a tall object in the electric field of a thundercloud.Plasma Physics Reports,2002,28(11):953-964.

    [19] Bazelyan E M,Raizer Y P,Aleksandrov N L,et al,Corona processes and lightning attachment:The effect of wind during thunderstorms.Atmospheric Research,2009,94(3):436-447.

    [20] 王才偉,陳茜,劉欣生,等.雷暴云下部正電荷中心產(chǎn)生的電場.高原氣象,1987,6(1):65-74.

    [21] 逯曦,張義軍,呂偉濤,等.近地面電暈電流組網(wǎng)觀測與數(shù)據(jù)分析.氣象科技,2010,38(6):746-751.

    [22] 趙中闊,郄秀書,張廣庶,等.雷暴云內(nèi)電場探測儀及初步實驗結(jié)果.高原氣象,2008,27(4):881-887.

    [23] Kasemir H W.Conference on Cloud Physics and Atmospheric Electricity of AMS.American Meteorological Society,Issaquah,Washington,1978.

    [24] D’Alessandro F,Berger G.Laboratory studies of corona emissions from air terminals.Journal of Physics D:Applied Physics,1999,32(21):2785-2790.

    [25] 張榴晨,徐松.有限元法在電磁計算中的應(yīng)用.北京:中國鐵道出版社,1996.

    [26] 史密斯(Smith I M).有限元方法編程(第三版).王崧,譯.北京:電子工業(yè)出版社,2003.

    [27] 金建銘,王建國,葛德彪.電磁場有限元方法.西安:西安電子科技大學(xué)出版社,1998.

    [28] 郭立新,李江挺,韓旭彪.計算物理學(xué).西安:西安電子科技大學(xué)出版社,2009.

    [29] Eriksson A J.The Lightning Ground Flash:An Engineering Study.University of KwaZulu-Natal,1979.

    [30] Moore C B.Improved configurations of lightning rods and air terminals.Journal of the Franklin Institute,1983,315(1):61-85.

    [31] 郭秀峰,譚涌波,郭鳳霞,等.建筑物尖端對大氣電場畸變影響的數(shù)值計算.應(yīng)用氣象學(xué)報,2012,12(1):136-151.

    [32] 言穆弘,郭昌明,劉欣生.雷暴對流起電機制理論分析.大氣科學(xué),1991,15(2):120-128.

    [33] 言穆弘,劉欣生.閃電先導(dǎo)靜電場波形理論分析.應(yīng)用氣象學(xué)報,1993,4(2):185-191.

    Onset Corona Field Based on Isolated Metal Tip

    Tan Yongbo1)2)Zhu Junru1)2)Li Xiangchao1)2)Liang Zhongwu1)2)Guo Xiufeng1)2)

    1)(Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters,Nanjing University of Information Science &Technology,Nanjing210044)
    2)(Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration,Nanjing University of Information Science &Technology,Nanjing210044)

    Corona discharges characterized by much lower current densities are often generated in a strong atmospheric electric field.What’s more,corona discharges are very important because the corona charge layers can change the distribution of the electric field near the ground and affect conditions of the initiation and development of an upward leader near different tall objects.Therefore,corona discharges have attracted the attention of scholars at home and abroad,especially in the onset corona field.However,how to effectively show corona emission threshold is still a difficult problem in the study of atmospheric electricity.

    In order to obtain the electric field strength at the tip of metal in the case of corona discharge,the onset corona voltages must be obtained by laboratory test for different heights,shapes and materials of metal tips which are arranged on the middle of lower plate of two parallel metal plates.Dimensions of two horizontal plates are 1.2mby 1.2mand the distance between them is 1m.Then environmental threshold electric fields are calculated for different heights,shapes and materials of metal tips.The two-dimensional Poisson equation is solved using finite element method and the corona emission thresholds at the tip of different types of metal tips are obtained.

    Effects on threshold electric fields by the height,material and shape are discussed,respectively.The following conclusions can be drawn from the experimental study.First,the environmental threshold fields are basically a linear decrease with heights of metal.The environmental threshold fields decrease first and then increase when metal tips are becoming sharp increasingly.Second,the height and shape of mental tips have no effects on corona emission threshold at the tip of metal tips.The corona threshold at the tip of mental tip is a constant.Finally,the corona triggering threshold is determined which isEc=158.75kV·m-1and the fitting function between corona emission threshold and spatial resolution is given.It will provide a reference for the judgment of starting moment of corona discharge for the future numerical simulation of corona discharge.

    corona discharge;laboratory simulation;corona field threshold;finite element;spatial resolution

    譚涌波,朱俊儒,李祥超,等.孤立金屬尖端的電暈觸發(fā)閾值.應(yīng)用氣象學(xué)報,2014,25(3):339-346.

    2013-05-28收到,2014-01-08收到再改稿。

    國家自然科學(xué)基金項目(41175003,41075003)

    *email:ybtan@ustc.edu

    猜你喜歡
    電暈尖端電場
    巧用對稱法 妙解電場題
    科學(xué)中國人(2018年8期)2018-07-23 02:26:56
    電場強度單個表達的比較
    電場中六個常見物理量的大小比較
    電暈籠內(nèi)導(dǎo)線交流電暈起始電壓判斷方法
    電暈籠內(nèi)導(dǎo)線三維電場強度分布淺析
    鏡頭看展
    500kV絕緣子串含零值絕緣子時的電暈放電分析
    電測與儀表(2014年6期)2014-04-04 11:59:46
    基于位移相關(guān)法的重復(fù)壓裂裂縫尖端應(yīng)力場研究
    斷塊油氣田(2014年5期)2014-03-11 15:33:49
    加速尖端機床國產(chǎn)化
    中文字幕久久专区| 女生性感内裤真人,穿戴方法视频| 亚洲成av人片免费观看| 乱人视频在线观看| 国产免费av片在线观看野外av| 国产一区二区在线观看日韩 | 欧美中文日本在线观看视频| 日本一二三区视频观看| 日韩av在线大香蕉| 久久久久久国产a免费观看| 久久中文看片网| 国产蜜桃级精品一区二区三区| 2021天堂中文幕一二区在线观| 97超级碰碰碰精品色视频在线观看| 亚洲欧美激情综合另类| 99在线视频只有这里精品首页| 亚洲av成人av| 精品国产三级普通话版| 久久久久亚洲av毛片大全| 亚洲 欧美 日韩 在线 免费| 国产一区二区三区视频了| 老熟妇仑乱视频hdxx| 午夜影院日韩av| 国产精品香港三级国产av潘金莲| 日韩成人在线观看一区二区三区| 午夜免费成人在线视频| 2021天堂中文幕一二区在线观| 热99在线观看视频| 婷婷六月久久综合丁香| 亚洲美女视频黄频| 欧美黄色片欧美黄色片| 亚洲国产欧美人成| 99国产综合亚洲精品| 国产精品免费一区二区三区在线| 免费高清视频大片| 国产黄a三级三级三级人| 久久久国产精品麻豆| 香蕉丝袜av| 欧美在线一区亚洲| www.www免费av| 天天躁日日操中文字幕| 久久久久九九精品影院| 99国产精品一区二区蜜桃av| 最近最新免费中文字幕在线| 黄片大片在线免费观看| 精品一区二区三区视频在线 | 精品一区二区三区人妻视频| 欧美一区二区精品小视频在线| 精品电影一区二区在线| 无遮挡黄片免费观看| 亚洲片人在线观看| 亚洲av成人av| 淫秽高清视频在线观看| 久久久久久大精品| 亚洲av一区综合| 最后的刺客免费高清国语| 男插女下体视频免费在线播放| 国产真实乱freesex| 日韩有码中文字幕| 精品午夜福利视频在线观看一区| 尤物成人国产欧美一区二区三区| 人妻丰满熟妇av一区二区三区| 亚洲国产精品久久男人天堂| 久久久久性生活片| 国产久久久一区二区三区| 99精品在免费线老司机午夜| 成人精品一区二区免费| 老司机深夜福利视频在线观看| 欧美日韩精品网址| 国内精品美女久久久久久| 日本黄大片高清| 欧美一区二区精品小视频在线| 成人国产一区最新在线观看| 亚洲自拍偷在线| 美女免费视频网站| 久久99热这里只有精品18| 亚洲国产欧洲综合997久久,| 九色成人免费人妻av| 天天添夜夜摸| 久久亚洲真实| 欧美乱妇无乱码| 性色avwww在线观看| 国产高清有码在线观看视频| 最近最新免费中文字幕在线| 69人妻影院| 久久天躁狠狠躁夜夜2o2o| 日韩高清综合在线| 久久精品91无色码中文字幕| 国产精华一区二区三区| 久久久久性生活片| 日本在线视频免费播放| 欧美一级a爱片免费观看看| 欧美乱色亚洲激情| 波野结衣二区三区在线 | 午夜免费成人在线视频| 亚洲美女视频黄频| www.999成人在线观看| 极品教师在线免费播放| 法律面前人人平等表现在哪些方面| 青草久久国产| 免费在线观看成人毛片| 精品99又大又爽又粗少妇毛片 | 美女黄网站色视频| 岛国视频午夜一区免费看| 亚洲18禁久久av| 少妇裸体淫交视频免费看高清| 免费人成视频x8x8入口观看| 国产亚洲精品久久久com| 一进一出抽搐动态| 成人亚洲精品av一区二区| 国产在线精品亚洲第一网站| 欧美+日韩+精品| 国产视频一区二区在线看| 成人国产一区最新在线观看| av国产免费在线观看| 精品乱码久久久久久99久播| 国产精品日韩av在线免费观看| 乱人视频在线观看| 色精品久久人妻99蜜桃| 国产淫片久久久久久久久 | 草草在线视频免费看| 有码 亚洲区| 极品教师在线免费播放| 中国美女看黄片| 99热6这里只有精品| 在线播放国产精品三级| 免费人成在线观看视频色| 久久精品亚洲精品国产色婷小说| 午夜福利在线观看免费完整高清在 | 国产精品亚洲一级av第二区| 国内少妇人妻偷人精品xxx网站| 精品久久久久久,| 99热只有精品国产| 免费观看人在逋| 久久九九热精品免费| 欧美3d第一页| 免费在线观看日本一区| 搞女人的毛片| 亚洲一区二区三区色噜噜| 欧美av亚洲av综合av国产av| 九色成人免费人妻av| 香蕉丝袜av| 黄片大片在线免费观看| 亚洲欧美精品综合久久99| 午夜精品一区二区三区免费看| 香蕉久久夜色| 女同久久另类99精品国产91| 欧美av亚洲av综合av国产av| 九色成人免费人妻av| 国内毛片毛片毛片毛片毛片| 亚洲av免费高清在线观看| 欧美日韩亚洲国产一区二区在线观看| 一级毛片女人18水好多| 激情在线观看视频在线高清| 一区二区三区高清视频在线| www日本黄色视频网| 国产成年人精品一区二区| 久久久精品欧美日韩精品| 桃红色精品国产亚洲av| 成年女人毛片免费观看观看9| 看免费av毛片| 欧美日韩综合久久久久久 | 成人国产一区最新在线观看| 精品国产三级普通话版| 嫩草影院精品99| www日本在线高清视频| 日本免费一区二区三区高清不卡| 国产伦精品一区二区三区视频9 | 日韩欧美精品免费久久 | 亚洲人与动物交配视频| 蜜桃久久精品国产亚洲av| 亚洲成av人片在线播放无| 在线观看66精品国产| 神马国产精品三级电影在线观看| 一区福利在线观看| 乱人视频在线观看| 最近在线观看免费完整版| 日本三级黄在线观看| 中文字幕久久专区| 国产老妇女一区| 美女免费视频网站| 久久久久久大精品| 日本一二三区视频观看| 桃色一区二区三区在线观看| 亚洲天堂国产精品一区在线| 亚洲国产精品合色在线| 精品人妻1区二区| 国产欧美日韩精品亚洲av| 一区二区三区免费毛片| 国产69精品久久久久777片| 国产激情欧美一区二区| 久久国产精品影院| 国产真实伦视频高清在线观看 | 欧美一级a爱片免费观看看| 国产毛片a区久久久久| 搡老熟女国产l中国老女人| av专区在线播放| 成年女人看的毛片在线观看| www日本黄色视频网| 欧美性猛交黑人性爽| a在线观看视频网站| 老司机午夜福利在线观看视频| 亚洲无线观看免费| 国产高潮美女av| 免费av毛片视频| 69人妻影院| 老司机午夜福利在线观看视频| 久久精品国产综合久久久| 身体一侧抽搐| 国产精品,欧美在线| 亚洲欧美日韩卡通动漫| 国产69精品久久久久777片| 美女被艹到高潮喷水动态| 在线观看66精品国产| 搡老岳熟女国产| 天天添夜夜摸| 美女高潮的动态| 在线播放无遮挡| 亚洲 国产 在线| 99精品欧美一区二区三区四区| 成人无遮挡网站| 老司机午夜福利在线观看视频| 国产亚洲精品av在线| 午夜福利在线观看免费完整高清在 | 欧美黑人巨大hd| 两人在一起打扑克的视频| 高潮久久久久久久久久久不卡| 亚洲中文字幕日韩| 久久久成人免费电影| 岛国在线免费视频观看| 在线看三级毛片| 99久久精品一区二区三区| 国内少妇人妻偷人精品xxx网站| 欧美黄色片欧美黄色片| 高清毛片免费观看视频网站| 国产精品永久免费网站| 国语自产精品视频在线第100页| 亚洲乱码一区二区免费版| 18禁在线播放成人免费| 国产精品亚洲av一区麻豆| 欧美高清成人免费视频www| 午夜免费成人在线视频| 噜噜噜噜噜久久久久久91| 亚洲成人久久爱视频| 脱女人内裤的视频| 欧美一区二区精品小视频在线| 欧美黑人欧美精品刺激| 久久精品91无色码中文字幕| 亚洲国产中文字幕在线视频| 无限看片的www在线观看| 欧美日本亚洲视频在线播放| 90打野战视频偷拍视频| 欧美日韩福利视频一区二区| 免费一级毛片在线播放高清视频| 每晚都被弄得嗷嗷叫到高潮| 国产精品免费一区二区三区在线| 日韩欧美精品免费久久 | 特大巨黑吊av在线直播| 国产探花极品一区二区| 久久欧美精品欧美久久欧美| 女生性感内裤真人,穿戴方法视频| 中国美女看黄片| 精品午夜福利视频在线观看一区| 亚洲在线观看片| 亚洲成人免费电影在线观看| 国产免费av片在线观看野外av| 久久久久久久精品吃奶| 久久伊人香网站| 国产黄片美女视频| 久久精品人妻少妇| 国产精品自产拍在线观看55亚洲| 两人在一起打扑克的视频| 亚洲一区二区三区不卡视频| 国内少妇人妻偷人精品xxx网站| 变态另类丝袜制服| 精品一区二区三区视频在线 | 日韩欧美一区二区三区在线观看| 最近最新免费中文字幕在线| 成人鲁丝片一二三区免费| 亚洲无线观看免费| 精品久久久久久,| 精品免费久久久久久久清纯| 中出人妻视频一区二区| 欧美+日韩+精品| 亚洲欧美日韩无卡精品| 99国产精品一区二区蜜桃av| 欧美成人免费av一区二区三区| 少妇高潮的动态图| 色噜噜av男人的天堂激情| 91久久精品国产一区二区成人 | 国产精品1区2区在线观看.| 男人舔女人下体高潮全视频| 偷拍熟女少妇极品色| 又爽又黄无遮挡网站| 熟妇人妻久久中文字幕3abv| 丁香欧美五月| 亚洲黑人精品在线| 欧洲精品卡2卡3卡4卡5卡区| 欧美激情久久久久久爽电影| 男人舔女人下体高潮全视频| 精品久久久久久久末码| 久久久久久久亚洲中文字幕 | 国产爱豆传媒在线观看| 亚洲无线在线观看| 精品欧美国产一区二区三| 真人做人爱边吃奶动态| 日本与韩国留学比较| 制服丝袜大香蕉在线| 国产精品久久久久久久久免 | 色综合欧美亚洲国产小说| 日本熟妇午夜| 久久久久久久久大av| 日本成人三级电影网站| 一区福利在线观看| 给我免费播放毛片高清在线观看| 中文字幕久久专区| 国产av在哪里看| а√天堂www在线а√下载| 亚洲国产中文字幕在线视频| 欧美三级亚洲精品| 亚洲精品色激情综合| 在线看三级毛片| 岛国在线观看网站| avwww免费| 在线观看午夜福利视频| 精品久久久久久久久久久久久| 国产99白浆流出| 99在线人妻在线中文字幕| 国产成人aa在线观看| 欧美高清成人免费视频www| 搡老岳熟女国产| 欧美成人一区二区免费高清观看| 一本精品99久久精品77| 国产高清视频在线播放一区| 精品久久久久久久毛片微露脸| 久久久国产成人精品二区| 三级男女做爰猛烈吃奶摸视频| 亚洲精品美女久久久久99蜜臀| 久99久视频精品免费| 韩国av一区二区三区四区| 久久久成人免费电影| 麻豆国产av国片精品| 亚洲午夜理论影院| 亚洲自拍偷在线| 十八禁人妻一区二区| 国产 一区 欧美 日韩| 在线观看av片永久免费下载| 天美传媒精品一区二区| 午夜久久久久精精品| 91字幕亚洲| 欧美zozozo另类| 一级毛片女人18水好多| 色在线成人网| 一本久久中文字幕| 亚洲成av人片在线播放无| 亚洲精品成人久久久久久| 久久久久九九精品影院| 成人精品一区二区免费| 男女下面进入的视频免费午夜| 美女cb高潮喷水在线观看| 美女被艹到高潮喷水动态| 搡老岳熟女国产| 久久精品91无色码中文字幕| 丁香六月欧美| 国产真人三级小视频在线观看| 狠狠狠狠99中文字幕| 中文字幕久久专区| 成人三级黄色视频| 欧美成狂野欧美在线观看| 丁香六月欧美| 久久精品国产自在天天线| 亚洲av免费高清在线观看| 亚洲一区二区三区不卡视频| 国产精品,欧美在线| 在线看三级毛片| 欧美日韩综合久久久久久 | 精品午夜福利视频在线观看一区| 日本黄色视频三级网站网址| 亚洲狠狠婷婷综合久久图片| 亚洲av电影不卡..在线观看| 精品国产超薄肉色丝袜足j| 99国产综合亚洲精品| 高清在线国产一区| 日韩人妻高清精品专区| 一夜夜www| 最近最新中文字幕大全电影3| 此物有八面人人有两片| 欧美成人一区二区免费高清观看| 怎么达到女性高潮| 美女大奶头视频| 欧美成人a在线观看| 99精品在免费线老司机午夜| 美女免费视频网站| 青草久久国产| 亚洲国产欧美网| 他把我摸到了高潮在线观看| 国产精品美女特级片免费视频播放器| 欧美绝顶高潮抽搐喷水| 波多野结衣高清作品| 女人被狂操c到高潮| 午夜福利在线观看免费完整高清在 | 1024手机看黄色片| 中文资源天堂在线| 18禁国产床啪视频网站| 午夜激情福利司机影院| 亚洲国产精品合色在线| 国产精品一区二区三区四区久久| 又爽又黄无遮挡网站| 欧美日韩黄片免| 最近最新中文字幕大全免费视频| 夜夜躁狠狠躁天天躁| 欧美日本亚洲视频在线播放| 亚洲精品久久国产高清桃花| 波多野结衣巨乳人妻| 久久精品国产自在天天线| 91在线精品国自产拍蜜月 | 国产私拍福利视频在线观看| 亚洲成a人片在线一区二区| 一级黄色大片毛片| 最近视频中文字幕2019在线8| 婷婷丁香在线五月| ponron亚洲| 午夜老司机福利剧场| 中文资源天堂在线| 国产精品免费一区二区三区在线| 国产97色在线日韩免费| 全区人妻精品视频| 日本a在线网址| 天堂av国产一区二区熟女人妻| av国产免费在线观看| 成人性生交大片免费视频hd| 国产精品1区2区在线观看.| 高清在线国产一区| 神马国产精品三级电影在线观看| 嫩草影院入口| 香蕉久久夜色| 变态另类丝袜制服| 99久久成人亚洲精品观看| 嫩草影视91久久| 99riav亚洲国产免费| 特级一级黄色大片| 老汉色av国产亚洲站长工具| 91麻豆av在线| 国产爱豆传媒在线观看| 久久天躁狠狠躁夜夜2o2o| 熟妇人妻久久中文字幕3abv| 在线观看午夜福利视频| 久久性视频一级片| 成人高潮视频无遮挡免费网站| 国产一区二区在线av高清观看| 法律面前人人平等表现在哪些方面| 最好的美女福利视频网| 禁无遮挡网站| 免费无遮挡裸体视频| 99国产综合亚洲精品| 亚洲欧美日韩高清在线视频| 亚洲18禁久久av| 国产亚洲欧美98| 一区二区三区高清视频在线| 国产精品久久视频播放| 日本熟妇午夜| 亚洲狠狠婷婷综合久久图片| 亚洲人成网站高清观看| 国产一区二区在线观看日韩 | 一级毛片女人18水好多| 丰满的人妻完整版| 精品午夜福利视频在线观看一区| 日韩欧美国产一区二区入口| 嫁个100分男人电影在线观看| www日本在线高清视频| 亚洲国产欧美人成| 国产高清三级在线| 婷婷六月久久综合丁香| 女人被狂操c到高潮| 伊人久久大香线蕉亚洲五| 久久九九热精品免费| 一个人看的www免费观看视频| 18禁在线播放成人免费| 国产探花在线观看一区二区| 日韩欧美一区二区三区在线观看| 欧美中文日本在线观看视频| 又粗又爽又猛毛片免费看| 国产亚洲av嫩草精品影院| 成年版毛片免费区| 日本免费a在线| 中文字幕av在线有码专区| www.色视频.com| 内地一区二区视频在线| 91av网一区二区| 国产麻豆成人av免费视频| av在线蜜桃| 一本久久中文字幕| 变态另类成人亚洲欧美熟女| 超碰av人人做人人爽久久 | 亚洲 国产 在线| 久9热在线精品视频| 亚洲中文字幕日韩| 最近最新免费中文字幕在线| 日本熟妇午夜| 欧美性感艳星| 少妇的逼水好多| 久久中文看片网| 欧洲精品卡2卡3卡4卡5卡区| 欧美bdsm另类| 色视频www国产| 欧美+亚洲+日韩+国产| 午夜精品久久久久久毛片777| 一区二区三区激情视频| 日韩av在线大香蕉| 一个人看视频在线观看www免费 | 人人妻人人看人人澡| 国产一区二区激情短视频| 亚洲av成人精品一区久久| 97人妻精品一区二区三区麻豆| 亚洲精品久久国产高清桃花| 观看免费一级毛片| 麻豆久久精品国产亚洲av| 亚洲av免费在线观看| 国模一区二区三区四区视频| 欧美国产日韩亚洲一区| 欧美区成人在线视频| 久久99热这里只有精品18| 日韩亚洲欧美综合| 亚洲午夜理论影院| 免费av毛片视频| 久久精品国产99精品国产亚洲性色| 99国产极品粉嫩在线观看| 男人的好看免费观看在线视频| 亚洲国产精品合色在线| 最近视频中文字幕2019在线8| 免费av观看视频| 极品教师在线免费播放| 国产亚洲精品综合一区在线观看| 国产一区二区激情短视频| 真人做人爱边吃奶动态| 少妇的逼好多水| 窝窝影院91人妻| 少妇的逼好多水| 国产激情欧美一区二区| 亚洲精品国产精品久久久不卡| 窝窝影院91人妻| 久久精品综合一区二区三区| 成人精品一区二区免费| 成人特级av手机在线观看| 少妇的逼水好多| 90打野战视频偷拍视频| 18+在线观看网站| 偷拍熟女少妇极品色| 99国产精品一区二区三区| 岛国在线免费视频观看| 国产乱人伦免费视频| 久久久久国内视频| 国产精品久久久久久久电影 | 在线a可以看的网站| 久久性视频一级片| 久久久久久九九精品二区国产| 午夜激情欧美在线| 人人妻人人澡欧美一区二区| 国产97色在线日韩免费| 午夜福利18| 精品日产1卡2卡| avwww免费| 中文在线观看免费www的网站| 在线免费观看的www视频| 黑人欧美特级aaaaaa片| 亚洲18禁久久av| 国产欧美日韩精品一区二区| 五月玫瑰六月丁香| 窝窝影院91人妻| 90打野战视频偷拍视频| 一本综合久久免费| 欧美成人性av电影在线观看| 亚洲欧美日韩东京热| 亚洲精华国产精华精| 亚洲中文字幕日韩| 最近在线观看免费完整版| 伊人久久大香线蕉亚洲五| e午夜精品久久久久久久| www日本黄色视频网| 精品99又大又爽又粗少妇毛片 | 免费一级毛片在线播放高清视频| 亚洲av电影在线进入| 成人永久免费在线观看视频| 最近最新中文字幕大全免费视频| 国产黄片美女视频| 国产成人影院久久av| 1000部很黄的大片| 老司机午夜十八禁免费视频| 欧美av亚洲av综合av国产av| 精品久久久久久成人av| 18禁国产床啪视频网站| 国产精品电影一区二区三区| or卡值多少钱| 成人高潮视频无遮挡免费网站| 久久久久久国产a免费观看| 在线国产一区二区在线| 国产高清有码在线观看视频| 久久国产精品影院| 精品欧美国产一区二区三| 国产精品影院久久| 欧美成人免费av一区二区三区| 国产欧美日韩一区二区三| 99精品在免费线老司机午夜| a级毛片a级免费在线| 欧美乱色亚洲激情| 亚洲精品影视一区二区三区av| 999久久久精品免费观看国产| 国产淫片久久久久久久久 | 内射极品少妇av片p| 天天躁日日操中文字幕| 亚洲中文字幕日韩| 51国产日韩欧美| 少妇人妻一区二区三区视频| 国产精品永久免费网站| 19禁男女啪啪无遮挡网站|