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

    霍爾推力器磁場(chǎng)位形及其優(yōu)化的數(shù)值研究

    2011-10-23 12:13:46鄧立赟藍(lán)紅梅
    物理學(xué)報(bào) 2011年2期
    關(guān)鍵詞:位形磁導(dǎo)率推力器

    鄧立赟 藍(lán)紅梅 劉 悅

    (大連理工大學(xué)物理與光電工程學(xué)院,大連 116024)

    (2010年4月20日收到;2010年5月6日收到修改稿)

    霍爾推力器磁場(chǎng)位形及其優(yōu)化的數(shù)值研究

    鄧立赟 藍(lán)紅梅 劉 悅?

    (大連理工大學(xué)物理與光電工程學(xué)院,大連 116024)

    (2010年4月20日收到;2010年5月6日收到修改稿)

    基于麥克斯韋方程,在軸對(duì)稱假設(shè)下建立了霍爾推力器磁場(chǎng)的數(shù)學(xué)模型.用有限差分方法對(duì)模型進(jìn)行了離散.給出了數(shù)值求解模型的迭代法.通過對(duì)模型的數(shù)值求解,得到了相應(yīng)的數(shù)值結(jié)果.通過對(duì)所得數(shù)值結(jié)果的分析,研究了磁場(chǎng)線圈電流變化對(duì)霍爾推力器磁場(chǎng)位形的影響.通過調(diào)整磁場(chǎng)線圈電流的大小找到了理想磁場(chǎng)位形.研究表明,對(duì)于理想磁場(chǎng)位形,內(nèi)通道的磁鏡比在3—3.5之間,外通道的磁鏡比在0.4—0.9之間;增加磁場(chǎng)線圈的電流,出口的磁場(chǎng)強(qiáng)度隨著增加,但不能增加磁鏡比.通道內(nèi)部的磁場(chǎng)強(qiáng)度幾乎不隨著磁場(chǎng)線圈電流的變化而變化.

    霍爾推力器,磁場(chǎng)位形,磁場(chǎng)線圈電流,磁鏡比

    PACS:52.75.Di,41.20.Gz

    1.引 言

    進(jìn)入21世紀(jì)以來,隨著人類探索宇宙空間范圍的拓展,傳統(tǒng)的化學(xué)火箭已不能滿足日益增加的對(duì)于深空探測(cè)的要求.為了克服傳統(tǒng)推進(jìn)技術(shù)的缺點(diǎn),目前許多發(fā)達(dá)國(guó)家都在研究各種新型推進(jìn)技術(shù)以滿足未來太空探索的需要.新型推進(jìn)技術(shù)當(dāng)中,電推進(jìn)技術(shù)是最具應(yīng)用前景的,因此在目前研究中最為廣泛.霍爾推力器是最典型的電推進(jìn)器,由于其具有結(jié)構(gòu)緊湊、比沖大、效率高和壽命長(zhǎng)等優(yōu)點(diǎn),可以在未來的深空探測(cè)中應(yīng)用于衛(wèi)星姿態(tài)控制、航天器軌道修正、軌道轉(zhuǎn)移、動(dòng)力補(bǔ)償?shù)热蝿?wù).

    霍爾推力器主要是由同軸的圓柱形通道、陶瓷器壁、以及進(jìn)氣口處的陽極與噴射口處的陰極、用于產(chǎn)生磁場(chǎng)的線圈等組成.霍爾推力器的工作氣體主要選用氙氣,也有人選用氪氣.工作氣體通過進(jìn)氣口進(jìn)入通道,同時(shí)陰極產(chǎn)生的電子進(jìn)入通道對(duì)工作氣體進(jìn)行電離.陽極和陰極之間的電勢(shì)差產(chǎn)生電場(chǎng),使氣體充分電離并加速等離子體中的帶電粒子(電子和離子).磁場(chǎng)線圈產(chǎn)生的磁場(chǎng)通過霍爾效應(yīng)磁化等離子體中的電子,使電子留在通道內(nèi),而沒有磁化的離子噴出通道,從而產(chǎn)生推力.這就是霍爾推力器的工作原理.

    由于霍爾推力器在未來的深空探測(cè)任務(wù)中有很好的應(yīng)用前景,因此受到廣泛研究.人們分別對(duì)霍爾推力器中等離子體的震蕩[1—5]、電子的反常遷移 和 輸 運(yùn)[6,7]、等 離 子 體 與 器 壁 材 料 的 相 互 作用[8—12]、等離子體鞘層[13—16]、二次電子發(fā)射[17—19]和磁場(chǎng)位形的效應(yīng)[20—25]等進(jìn)行了深入研究.

    眾所周知,霍爾推力器是通過霍爾效應(yīng)來產(chǎn)生推力,而磁場(chǎng)是產(chǎn)生霍爾效應(yīng)的關(guān)鍵.同時(shí),磁場(chǎng)位形對(duì)于霍爾推力器放電特性、使用壽命和效率等許多方面產(chǎn)生很大影響,所以對(duì)霍爾推力器磁場(chǎng)位形的研究非常重要,既對(duì)霍爾推力器的機(jī)理研究有理論指導(dǎo)意義,也對(duì)霍爾推力器的設(shè)計(jì)和實(shí)驗(yàn)提供理論基礎(chǔ).2002年,Peterson等[20]研究了霍爾推力器的磁場(chǎng)結(jié)構(gòu).2003年,Garrigues等[21]用數(shù)值模擬的方法研究了霍爾推力器的磁場(chǎng)位形對(duì)壽命的影響.2005 年,Keidar和 Boyd[22]研究了霍爾推力器中的磁鏡效應(yīng).2009年,鄂鵬等[23]研究了霍爾推力器中的磁場(chǎng)位形,于達(dá)仁等[24]研究了霍爾推力器磁場(chǎng)位形對(duì)羽流的影響,耿少飛等[25]研究了霍爾推力器磁場(chǎng)對(duì)電子和離子的影響.

    本文基于麥克斯韋方程,在軸對(duì)稱假設(shè)下建立了霍爾推力器磁場(chǎng)的數(shù)學(xué)模型,模型中包含有鐵物質(zhì),使得其具有很強(qiáng)的非線性,很難求解.本文用有限差分法對(duì)模型進(jìn)行了離散,給出了數(shù)值求解模型的迭代法,并對(duì)模型進(jìn)行了數(shù)值求解,得到了相應(yīng)的數(shù)值結(jié)果.通過對(duì)所得數(shù)值結(jié)果的分析,研究了磁場(chǎng)線圈電流變化對(duì)霍爾推力器磁場(chǎng)位形的影響.

    2.霍爾推力器磁場(chǎng)的數(shù)學(xué)模型

    對(duì)于霍爾推力器的磁場(chǎng),由麥克斯韋方程有

    (1)式中磁場(chǎng)B可以用磁矢勢(shì)A表示

    根據(jù)霍爾推力器的結(jié)構(gòu),霍爾推力器的磁場(chǎng)可以認(rèn)為是軸對(duì)稱的.在軸對(duì)稱的假設(shè)下,在柱坐標(biāo)系 (r,θ,z)中,磁矢勢(shì)A只有θ分量,而且只是r和z的函數(shù),因此,磁矢勢(shì)A可以表示為

    電流密度J可以表示為

    將(5)式代入(4)式有

    由(7)式有

    將(3),(6)和(7)式代入(2)式有

    方程(10)的邊界條件為

    圖1 霍爾推力器剖面圖

    方程(8),(9)和(10)與上述邊界條件組成了霍爾推力器磁場(chǎng)的數(shù)學(xué)模型.

    在實(shí)際計(jì)算中,計(jì)算區(qū)域是一個(gè)包含推力器剖面的一個(gè)矩形區(qū)域,計(jì)算區(qū)域的邊界離推力器足夠遠(yuǎn),我們?nèi)【嗤屏ζ髯銐蜻h(yuǎn)處為無窮遠(yuǎn),并在這個(gè)邊界上采用上面的邊界條件.

    推力器區(qū)域如圖1所示.其中,有3個(gè)磁場(chǎng)線圈分別定義為1號(hào)線圈,2號(hào)線圈,3號(hào)線圈,相應(yīng)的電流分別為 NI1,NI2和 NI3,灰色區(qū)域則為鐵物質(zhì)區(qū)域.

    計(jì)算區(qū)域中有鐵物質(zhì),而鐵物質(zhì)的相對(duì)磁導(dǎo)率需要通過圖2所示的μ—B曲線得到.方程(10)中的磁導(dǎo)率μ在計(jì)算區(qū)域中是一個(gè)有空間分布的函數(shù).在沒有鐵物質(zhì)的區(qū)域里,我們?nèi)〈艑?dǎo)率為真空磁導(dǎo)率;在有鐵物質(zhì)的區(qū)域里,先通過μ—B曲線計(jì)算出相對(duì)磁導(dǎo)率,再乘以真空磁導(dǎo)率作為鐵物質(zhì)區(qū)域的磁導(dǎo)率.

    圖2 μ—B曲線

    在計(jì)算霍爾推力器磁場(chǎng)的時(shí)候,需要在上述邊界條件下求解方程(8),(9)和(10),這時(shí)需要用到磁導(dǎo)率μ;而利用μ—B曲線計(jì)算相對(duì)磁導(dǎo)率μ時(shí),需要知道磁場(chǎng)B.因此,這是一個(gè)具有很強(qiáng)非線性的偏微分方程邊值問題,需要數(shù)值求解,并且需要用迭代法.強(qiáng)非線性使得求解很難.

    3.數(shù)值方法

    首先,用有限差分法對(duì)模型進(jìn)行離散.方程(10)的差分格式如下:

    方程(8)的差分格式如下:

    方程(9)的差分格式如下:

    在離散的過程中,假設(shè)r方向和z方向的空間步長(zhǎng)相等,都是h.整理(11)式得

    其中

    方程(14)是個(gè)具有很強(qiáng)非線性的代數(shù)方程組.為了求解這個(gè)方程組,我們構(gòu)造如下的迭代法:

    7)輸出結(jié)果.

    4.數(shù)值結(jié)果與討論

    我們用上面的迭代法計(jì)算了3個(gè)磁場(chǎng)線圈在不同電流情況下產(chǎn)生的磁場(chǎng),通過改變3個(gè)磁場(chǎng)線圈電流的大小試圖找到理想磁場(chǎng)位形.

    理想磁場(chǎng)位形的判斷條件為:1)在通道的中間存在磁場(chǎng)零點(diǎn);2)磁場(chǎng)在通道出口處沿軸向增加;3)存在磁透鏡位形的磁場(chǎng).

    圖3表示當(dāng)3個(gè)磁場(chǎng)線圈電流分別為NI1=200 A,NI2=350 A和NI3=35 A時(shí)整個(gè)計(jì)算區(qū)域的磁場(chǎng)位形.

    圖3 整個(gè)計(jì)算區(qū)域磁通 NI1=200 A,NI2=350 A,NI3=35 A

    圖4 霍爾推力器磁通 NI1=200 A,NI2=350 A,NI3=35 A

    圖4表示當(dāng)3個(gè)磁場(chǎng)線圈電流分別為NI1=200 A,NI2=350 A和NI3=35 A時(shí)霍爾推力器的磁場(chǎng)位形.可以看到,在通道中間存在磁場(chǎng)零點(diǎn).磁場(chǎng)零點(diǎn)的存在使電子被束縛在磁阱里,能夠更好地增加電子對(duì)氣體的電離,提高電離效率.

    圖5表示當(dāng)3個(gè)磁場(chǎng)線圈電流分別為NI1=200 A,NI2=350 A和NI3=35 A時(shí)霍爾推力器通道內(nèi)的磁場(chǎng)位形.由圖5可以看出,在通道的出口處為磁透鏡位形的磁場(chǎng),磁透鏡磁場(chǎng)能夠?qū)﹄x子起到匯聚作用,減小羽流散射角,提高推進(jìn)器的效率.

    根據(jù)文獻(xiàn)[22],磁鏡比定義為 Bw/B0,其中 Bw表示器壁處的磁場(chǎng)強(qiáng)度,B0表示通道出口處中心磁場(chǎng)強(qiáng)度.此時(shí)通過計(jì)算,內(nèi)外通道磁鏡比分別為3.2730,0.6077,出口處中心磁場(chǎng)大小為 234×10-4T.

    圖5 霍爾推力器通道內(nèi)的磁通 NI1=200 A,NI2=350 A,NI3=35 A

    圖6 霍爾推力器通道中心的磁場(chǎng) NI1=200 A,NI2=350 A,NI3=35 A

    圖6為磁場(chǎng)通道中心從陽極區(qū)域沿軸向到出口處磁場(chǎng)強(qiáng)度的變化.可以看出磁場(chǎng)沿軸向從陽極區(qū)域到出口處不斷增加,即在此組磁場(chǎng)線圈電流下找到了理想磁場(chǎng)位形.

    經(jīng)過大量數(shù)值計(jì)算,我們?cè)?0組電流強(qiáng)度下找到了理想的磁場(chǎng)位形,表1給出了相應(yīng)的內(nèi)外通道磁鏡比和通道出口中心處的磁場(chǎng)強(qiáng)度.從表中可以看出,在理想磁場(chǎng)位形下,隨著磁場(chǎng)線圈電流強(qiáng)度的增加,通道出口處磁場(chǎng)強(qiáng)度也隨之增加,但是內(nèi)外通道磁鏡比并沒有隨著電流強(qiáng)度的增加而增加.

    圖7分別給出了5組參數(shù)下通道中央沿軸向的磁場(chǎng)分布.由圖7可以看出,在理想磁場(chǎng)位形下,磁場(chǎng)隨著電流強(qiáng)度的增加而增加,但磁場(chǎng)強(qiáng)度的增加主要在通道出口處附近,通道內(nèi)部磁場(chǎng)幾乎不會(huì)受到磁場(chǎng)線圈電流大小的影響.

    表1 理想磁場(chǎng)位形中的磁鏡比和出口磁場(chǎng)強(qiáng)度

    5.結(jié) 論

    了霍爾推力器磁場(chǎng)的數(shù)學(xué)模型.由于模型中包含有鐵物質(zhì),使得其具有很強(qiáng)的非線性,很難求解.本文采用有限差分法對(duì)模型進(jìn)行了離散,給出了數(shù)值求解模型的迭代法,并對(duì)模型進(jìn)行了數(shù)值求解,得到了相應(yīng)的數(shù)值結(jié)果.通過對(duì)所得數(shù)值結(jié)果的分析,研究了3個(gè)磁場(chǎng)線圈中電流的變化對(duì)霍爾推力器磁場(chǎng)位形的影響.通過調(diào)整磁場(chǎng)線圈電流的大小可以找到理相的磁場(chǎng)位形,在理想的磁場(chǎng)位形下內(nèi)通道的磁鏡比在3—3.5之間,外通道的磁鏡比在0.4—0.9之間.理想磁場(chǎng)位形下,增加磁場(chǎng)線圈的電流,通道出口處的磁場(chǎng)強(qiáng)度隨著增加,但不能增加磁鏡比.同時(shí)通道內(nèi)部的磁場(chǎng)強(qiáng)度幾乎不會(huì)隨著磁場(chǎng)線圈電流的變化而變化.

    感謝哈爾濱工業(yè)大學(xué)于達(dá)仁教授及其研究團(tuán)隊(duì)的老師和學(xué)生們給予的幫助和支持,感謝大連理工大學(xué)等離子體理論與模擬科技創(chuàng)新團(tuán)隊(duì)的老師和同學(xué)提出的寶貴意見.

    本文基于麥克斯韋方程,在軸對(duì)稱假設(shè)下建立

    [1]Boeuf J P,Garrigues L 1998Appl.Phys.Lett.84 3541

    [2]Guerrini G,Michaut C 1999Phys.Plasmas6 343

    [3]Choueiri E Y 2001Phys.Plasmas8 1411

    [4]Yu D R,Wei L Q,Zhao Z Y,Han K,Yan G Y 2008Phys.Plasmas15 043502

    [5]Yu D R,Wang C S,Wei L Q,Gao C,Yu G 2008Phys.Plasmas15 113503

    [6]Meezank N B,Hargus W A,Cappelli Jr M A 2001Phys.Rev.E 63 026401

    [7]Smirnov A,Raitses Y,F(xiàn)isch N J 2004Phys.Plasmas11 4922

    [8]Meezan N B,Cappelli M A 2002Phys.Rev.E 66 036401

    [9]Ahedo E,Gallardo J M,Martinez-Sanchez M 2003Phys.Plasmas10 3397

    [10]Gascon N,Dudeck M,Barral S 2003Phys.Plasmas10 4123

    [11]Barral S,Makowski K,Peradzynski Z,Gascon N,Dudeck M 2003Phys.Plasmas10 4137

    [12]Yu D R,Li Y Q,Song S H 2006J.Phys.D:Applied Physics39 2205

    [13]Raitses Y,Staack D,Smirnov A,F(xiàn)isch N J 2005Phys.Plasmas12 073507

    [14]Francesco T,Savino L,Mario C 2005Phys.Plasmas12 093506

    [15]Francesco T,Savino L,Mario C,Ralf S 2009Appl.Phys.Lett.94 251502

    [16]Yu D R,Zhang F K,Li H,Liu H 2009Acta Phys.Sin.58 1844(in Chinese)[于達(dá)仁、張鳳奎、李 鴻、劉 輝 2009物理學(xué)報(bào) 58 1184]

    [17]Sydorenkoa D,Smolyakov A,Kaganovich I,Raitses Y 2006Phys.Plasmas13 014501

    [18]Raitses Y,Smirnov A,Staack D,F(xiàn)isch N J 2006Phys.Plasmas13 014502

    [19]Garrigues L,Hagelaar G J M,Boniface C,Boeuf J P 2006J.Appl.Phys.100 123301

    [20]Peterson P Y,Gallimore A D,Haas J M 2002Phys.Plasmas9 4354

    [21]Garrigues L,Hagelaar G J M,Bareilles J,Boniface C,Boeuf J P 2003Phys.Plasmas10 4886

    [22]Keidar M,Boyd I D 2005Appl.Phys.Lett.87 121501

    [23]E P,Yu D R,Wu Z W,Han K 2009Acta Phys.Sin.58 2535(in Chinese)[鄂 鵬、于達(dá)仁、武志文、韓 軻 2009物理學(xué)報(bào) 58 2535]

    [24]Yu D R,Li J,Li H,Li Y,Jiang B H 2009Plasma Sci.Technol.11 714

    [25]Geng S F,Tang D L,Zhao J,Qiu X M 2009Acta Phys.Sin.58 5520(in Chinese)[耿少飛、唐德禮、趙 杰、邱孝明 2009物理學(xué)報(bào)58 5520]

    PACS:52.75.Di,41.20.Gz

    Numerical study on Hall thruster magnetic configuration and its optimization

    Deng Li-Yun Lan Hong-Mei Liu Yue?
    (School of Physics and Optoelectronic Technology,Dalian University of Technology,Dalian 116024,China)
    (Received 20 April 2010;revised manuscript received 6 May 2010)

    Based on Maxwell's equations,under axis-symmetric assumption,a mathematical model of the magnetic field of a Hall thruster is established.A finite difference method is used for discretizing the model.An iterative method for numerically solving the model is given.The numerical results are obtained.From analysis of the results,the effect of the current in the magnetic field coils on the magnetic configuration of a Hall thruster is investigated.Through adjusting the current,the ideal magnetic configurations are found.It was shown that,for ideal magnetic configuration,the magnetic mirror ratio in the inner channel is between 3 and 3.5,and the magnetic mirror ratio in the outer channel is between 0.4 and 0.9.With the increasing of the current in the magnetic field coils,the magnetic field at the exit of the channel increases,but the magnetic mirror ratio cannot be increased.The magnetic field in the inner channel hardly changes with change of the currents in the magnetic field coils.

    Hall thruster,magnetic configuration,current in magnetic field coil,magnetic mirror ratio

    ?通訊聯(lián)系人.E-mail:liuyue@dlut.edu.cn

    ?Corresponding author.E-mail:liuyue@dlut.edu.cn

    猜你喜歡
    位形磁導(dǎo)率推力器
    單組元推力器倒置安裝多余物控制技術(shù)驗(yàn)證
    寬頻高磁導(dǎo)率R10k軟磁材料的開發(fā)
    山東冶金(2022年2期)2022-08-08 01:50:52
    中間支撐剛度對(duì)雙跨梁屈曲穩(wěn)定性的影響
    基于FEMM的永磁電機(jī)動(dòng)態(tài)凍結(jié)磁導(dǎo)率并行仿真及程序
    基于旋量理論的四自由度抓取機(jī)械手奇異位形分析
    Advanced Frozen Permeability Technique and Applications in Developing High Performance Electrical Machines
    用于小行星探測(cè)的離子推力器技術(shù)研究
    基于可操作度的機(jī)器人最優(yōu)初始位形研究
    大眾科技(2015年11期)2015-11-24 01:57:16
    基于最優(yōu)初始位形的冗余度機(jī)器人可操作度優(yōu)化*
    離子推力器和霍爾推力器的異同
    太空探索(2015年6期)2015-07-12 12:48:42
    亚洲国产色片| 中文欧美无线码| videosex国产| 九草在线视频观看| 五月开心婷婷网| 欧美三级亚洲精品| 日本-黄色视频高清免费观看| 国产高清不卡午夜福利| 亚洲成人av在线免费| 你懂的网址亚洲精品在线观看| 免费观看在线日韩| 日韩av不卡免费在线播放| 国产精品国产三级专区第一集| 黄片无遮挡物在线观看| 国产精品秋霞免费鲁丝片| 亚洲内射少妇av| 最近中文字幕高清免费大全6| 日韩免费高清中文字幕av| 99九九在线精品视频| 人妻一区二区av| 国产一区二区三区综合在线观看 | 亚洲欧美中文字幕日韩二区| 啦啦啦啦在线视频资源| 久久精品国产a三级三级三级| 婷婷成人精品国产| 欧美日韩亚洲高清精品| 老熟女久久久| 亚洲精华国产精华液的使用体验| 亚洲av在线观看美女高潮| 99re6热这里在线精品视频| 精品酒店卫生间| 欧美三级亚洲精品| 免费高清在线观看视频在线观看| 一级,二级,三级黄色视频| 久久久久久久久久久久大奶| 久久久国产精品麻豆| 亚洲av二区三区四区| 亚洲国产精品一区二区三区在线| 制服丝袜香蕉在线| 国产精品国产三级专区第一集| 亚洲熟女精品中文字幕| 亚洲精品亚洲一区二区| 如日韩欧美国产精品一区二区三区 | 一本久久精品| 国产亚洲精品久久久com| av国产久精品久网站免费入址| 成人18禁高潮啪啪吃奶动态图 | 免费黄网站久久成人精品| 亚洲美女搞黄在线观看| 亚洲第一av免费看| 另类亚洲欧美激情| 成人免费观看视频高清| 夜夜看夜夜爽夜夜摸| 男的添女的下面高潮视频| 观看av在线不卡| 精品午夜福利在线看| 中国国产av一级| 啦啦啦视频在线资源免费观看| 久热久热在线精品观看| 老司机影院毛片| 亚洲国产av新网站| 一区二区三区免费毛片| 夜夜爽夜夜爽视频| 天天操日日干夜夜撸| 亚洲综合色网址| 老司机影院成人| 亚洲三级黄色毛片| 汤姆久久久久久久影院中文字幕| 久久午夜综合久久蜜桃| 18禁观看日本| 日韩熟女老妇一区二区性免费视频| 看非洲黑人一级黄片| 国产av国产精品国产| 国产欧美日韩一区二区三区在线 | 亚洲综合精品二区| h视频一区二区三区| 国产精品 国内视频| 亚洲精品国产色婷婷电影| 一级黄片播放器| 国产精品不卡视频一区二区| 最新中文字幕久久久久| 国语对白做爰xxxⅹ性视频网站| videos熟女内射| 你懂的网址亚洲精品在线观看| 国产精品人妻久久久影院| 色婷婷久久久亚洲欧美| 日韩成人伦理影院| 国产欧美日韩综合在线一区二区| 国产精品久久久久成人av| 欧美精品高潮呻吟av久久| 人妻夜夜爽99麻豆av| 国产亚洲午夜精品一区二区久久| 亚洲五月色婷婷综合| 18+在线观看网站| 黄片无遮挡物在线观看| 欧美日韩国产mv在线观看视频| 一个人免费看片子| 亚洲一区二区三区欧美精品| 日日摸夜夜添夜夜爱| 黑人巨大精品欧美一区二区蜜桃 | 色吧在线观看| 国产精品嫩草影院av在线观看| 男人添女人高潮全过程视频| 亚洲美女黄色视频免费看| 啦啦啦中文免费视频观看日本| 九九爱精品视频在线观看| 亚洲人成77777在线视频| av免费在线看不卡| 曰老女人黄片| 十八禁网站网址无遮挡| 亚州av有码| 精品人妻熟女av久视频| 日韩中文字幕视频在线看片| 中文字幕精品免费在线观看视频 | av网站免费在线观看视频| 国国产精品蜜臀av免费| 中文字幕亚洲精品专区| 国产成人精品无人区| a级毛片免费高清观看在线播放| 一级黄片播放器| 国产在线视频一区二区| 日本wwww免费看| 亚洲综合精品二区| 免费看不卡的av| 蜜臀久久99精品久久宅男| 国产片内射在线| 国产精品久久久久久精品古装| 性高湖久久久久久久久免费观看| xxx大片免费视频| 亚洲人与动物交配视频| 精品人妻熟女毛片av久久网站| 精品国产乱码久久久久久小说| 青春草国产在线视频| 国产高清三级在线| 一区二区三区精品91| 黄色毛片三级朝国网站| 美女中出高潮动态图| 蜜臀久久99精品久久宅男| av在线播放精品| 中文字幕人妻熟人妻熟丝袜美| 三级国产精品欧美在线观看| 亚洲精品久久午夜乱码| 国产一区有黄有色的免费视频| 亚洲av欧美aⅴ国产| 青春草亚洲视频在线观看| 中文字幕制服av| 九九在线视频观看精品| 日韩强制内射视频| 精品人妻一区二区三区麻豆| 纵有疾风起免费观看全集完整版| 一区二区三区四区激情视频| 久久精品熟女亚洲av麻豆精品| 亚洲精品,欧美精品| 卡戴珊不雅视频在线播放| 日本91视频免费播放| 久久免费观看电影| 国产免费一区二区三区四区乱码| 赤兔流量卡办理| 男人操女人黄网站| 在线 av 中文字幕| 亚洲国产日韩一区二区| 桃花免费在线播放| 建设人人有责人人尽责人人享有的| 五月开心婷婷网| 亚州av有码| 青春草国产在线视频| 国产免费现黄频在线看| 欧美另类一区| videos熟女内射| 亚洲激情五月婷婷啪啪| 精品少妇内射三级| 国产精品国产av在线观看| 五月伊人婷婷丁香| 久久精品夜色国产| 九九在线视频观看精品| 校园人妻丝袜中文字幕| 男女无遮挡免费网站观看| 国产黄频视频在线观看| 色婷婷久久久亚洲欧美| 成人免费观看视频高清| 中国国产av一级| av线在线观看网站| 欧美亚洲 丝袜 人妻 在线| 精品少妇内射三级| 午夜福利,免费看| 色婷婷av一区二区三区视频| 性色av一级| 99视频精品全部免费 在线| 欧美少妇被猛烈插入视频| 欧美 日韩 精品 国产| 人人妻人人澡人人爽人人夜夜| 国产国拍精品亚洲av在线观看| 国产亚洲一区二区精品| 黄色配什么色好看| 精品国产乱码久久久久久小说| 老司机亚洲免费影院| 91精品三级在线观看| 十分钟在线观看高清视频www| 免费黄色在线免费观看| 高清视频免费观看一区二区| 亚洲精品乱久久久久久| 欧美日韩亚洲高清精品| 天天操日日干夜夜撸| 成人漫画全彩无遮挡| 精品99又大又爽又粗少妇毛片| 国产免费一级a男人的天堂| 国产黄片视频在线免费观看| 在线观看免费视频网站a站| 3wmmmm亚洲av在线观看| 在线免费观看不下载黄p国产| 国产精品一区二区三区四区免费观看| 成人影院久久| 亚洲第一av免费看| www.色视频.com| 97精品久久久久久久久久精品| 成人免费观看视频高清| 国产精品三级大全| 极品少妇高潮喷水抽搐| 最黄视频免费看| 99九九在线精品视频| 在线播放无遮挡| 中国三级夫妇交换| 九九爱精品视频在线观看| 日韩一本色道免费dvd| 91精品伊人久久大香线蕉| 99久久精品一区二区三区| 日韩不卡一区二区三区视频在线| 亚洲精品第二区| 国产av精品麻豆| 9色porny在线观看| 久久婷婷青草| 精品国产一区二区久久| 妹子高潮喷水视频| 亚洲av欧美aⅴ国产| 男男h啪啪无遮挡| av视频免费观看在线观看| 亚洲av成人精品一二三区| 免费黄色在线免费观看| 男女免费视频国产| 人妻人人澡人人爽人人| 久久久久国产网址| 美女xxoo啪啪120秒动态图| 午夜影院在线不卡| 日韩制服骚丝袜av| 婷婷色av中文字幕| 天天操日日干夜夜撸| 午夜福利视频精品| 国产一区二区三区综合在线观看 | 国产片特级美女逼逼视频| 久久久久国产网址| 观看美女的网站| 日韩视频在线欧美| 日本vs欧美在线观看视频| 性色av一级| 亚洲精品aⅴ在线观看| 少妇猛男粗大的猛烈进出视频| 国产又色又爽无遮挡免| 九色成人免费人妻av| 黑人巨大精品欧美一区二区蜜桃 | 免费大片18禁| 久久久欧美国产精品| 一级毛片 在线播放| 香蕉精品网在线| 精品人妻偷拍中文字幕| 亚洲美女黄色视频免费看| 亚洲国产精品成人久久小说| 一级毛片黄色毛片免费观看视频| 久久午夜福利片| 伦精品一区二区三区| 爱豆传媒免费全集在线观看| 极品少妇高潮喷水抽搐| 精品国产一区二区久久| 99热这里只有精品一区| 亚洲精品久久久久久婷婷小说| 2021少妇久久久久久久久久久| 国产精品一区二区三区四区免费观看| 亚洲精品久久午夜乱码| 欧美人与善性xxx| 国产亚洲午夜精品一区二区久久| av在线app专区| 中文欧美无线码| 曰老女人黄片| 久久99一区二区三区| 美女内射精品一级片tv| 精品酒店卫生间| 在线观看国产h片| av在线老鸭窝| 国产精品国产三级国产av玫瑰| 韩国av在线不卡| 久久99一区二区三区| 成人无遮挡网站| 国产在线视频一区二区| 午夜免费鲁丝| 国产精品人妻久久久久久| 在线观看三级黄色| 亚洲第一av免费看| 亚洲综合色惰| 新久久久久国产一级毛片| 日韩伦理黄色片| 在线观看www视频免费| 国产永久视频网站| 亚洲国产精品成人久久小说| 纵有疾风起免费观看全集完整版| 免费日韩欧美在线观看| 寂寞人妻少妇视频99o| 亚洲国产最新在线播放| 国产精品一区二区在线观看99| 观看美女的网站| 在线播放无遮挡| 国国产精品蜜臀av免费| 黄色欧美视频在线观看| 久久久国产一区二区| 久久久a久久爽久久v久久| 草草在线视频免费看| 飞空精品影院首页| 亚洲美女黄色视频免费看| 少妇人妻精品综合一区二区| 日韩av不卡免费在线播放| 五月天丁香电影| 中文字幕人妻丝袜制服| 丰满少妇做爰视频| 欧美日韩综合久久久久久| 亚洲色图综合在线观看| 欧美激情国产日韩精品一区| 黄色怎么调成土黄色| 免费看不卡的av| 七月丁香在线播放| 极品少妇高潮喷水抽搐| 日韩伦理黄色片| 美女视频免费永久观看网站| 免费观看a级毛片全部| 久久精品国产亚洲网站| 国产精品偷伦视频观看了| 亚洲精品第二区| 日韩人妻高清精品专区| 久久久久国产精品人妻一区二区| av网站免费在线观看视频| 热99国产精品久久久久久7| 精品国产一区二区三区久久久樱花| 日日摸夜夜添夜夜爱| 国产在线视频一区二区| 啦啦啦视频在线资源免费观看| 国产探花极品一区二区| 亚洲色图综合在线观看| 成人免费观看视频高清| 一区二区三区免费毛片| 插逼视频在线观看| 夜夜骑夜夜射夜夜干| 国产av精品麻豆| 精品卡一卡二卡四卡免费| 九草在线视频观看| 热re99久久国产66热| 女人久久www免费人成看片| 亚洲精品久久午夜乱码| 国产成人精品福利久久| 精品卡一卡二卡四卡免费| 成人亚洲精品一区在线观看| 能在线免费看毛片的网站| 看十八女毛片水多多多| 丝袜美足系列| 国产av码专区亚洲av| 精品卡一卡二卡四卡免费| 亚洲av二区三区四区| 毛片一级片免费看久久久久| 波野结衣二区三区在线| 美女cb高潮喷水在线观看| 熟女电影av网| av不卡在线播放| av.在线天堂| 2021少妇久久久久久久久久久| 免费观看在线日韩| 久久久久精品性色| 亚洲精品自拍成人| 丰满饥渴人妻一区二区三| 婷婷色麻豆天堂久久| 麻豆乱淫一区二区| 亚洲精品日韩av片在线观看| 免费黄网站久久成人精品| 亚洲欧洲国产日韩| 中文字幕人妻熟人妻熟丝袜美| 国产成人免费观看mmmm| tube8黄色片| 国产精品人妻久久久影院| 国产免费又黄又爽又色| 多毛熟女@视频| 国产av一区二区精品久久| 伦理电影免费视频| 伦精品一区二区三区| 久久久久久久精品精品| 免费看光身美女| 久久97久久精品| 国产在线免费精品| 亚洲高清免费不卡视频| 99国产综合亚洲精品| 免费观看的影片在线观看| 永久网站在线| 欧美日韩av久久| 一级,二级,三级黄色视频| 99热6这里只有精品| 在线观看免费日韩欧美大片 | 一级毛片 在线播放| 又黄又爽又刺激的免费视频.| 菩萨蛮人人尽说江南好唐韦庄| 国产日韩欧美亚洲二区| 国产亚洲欧美精品永久| kizo精华| 日日爽夜夜爽网站| 人体艺术视频欧美日本| 中文乱码字字幕精品一区二区三区| 欧美xxⅹ黑人| 最近的中文字幕免费完整| 久久99精品国语久久久| 亚洲精品国产av蜜桃| 欧美日韩精品成人综合77777| 蜜桃国产av成人99| 国产一区二区在线观看av| 午夜福利视频在线观看免费| 国产精品免费大片| 人人妻人人爽人人添夜夜欢视频| 久久 成人 亚洲| 免费高清在线观看日韩| 成年人午夜在线观看视频| 欧美精品人与动牲交sv欧美| 麻豆精品久久久久久蜜桃| 如日韩欧美国产精品一区二区三区 | 麻豆成人av视频| 亚洲,一卡二卡三卡| 国产精品一区二区在线观看99| 春色校园在线视频观看| 国产永久视频网站| 欧美亚洲日本最大视频资源| 色婷婷av一区二区三区视频| 美女中出高潮动态图| 亚洲欧洲精品一区二区精品久久久 | 天天躁夜夜躁狠狠久久av| 一边摸一边做爽爽视频免费| av免费观看日本| 久久精品国产鲁丝片午夜精品| 亚州av有码| 国产av一区二区精品久久| 青青草视频在线视频观看| 女人久久www免费人成看片| 人妻人人澡人人爽人人| 日韩人妻高清精品专区| 国产精品久久久久久久久免| 极品人妻少妇av视频| 日韩视频在线欧美| 9色porny在线观看| 免费av不卡在线播放| 久久国内精品自在自线图片| 天天躁夜夜躁狠狠久久av| 免费看av在线观看网站| www.色视频.com| 大又大粗又爽又黄少妇毛片口| 男女高潮啪啪啪动态图| www.av在线官网国产| 国语对白做爰xxxⅹ性视频网站| 欧美精品一区二区大全| 国产69精品久久久久777片| 国产欧美日韩综合在线一区二区| 久久久久久久亚洲中文字幕| 99久久中文字幕三级久久日本| 久久亚洲国产成人精品v| 丰满乱子伦码专区| 国产探花极品一区二区| 人妻人人澡人人爽人人| 久久免费观看电影| 18禁动态无遮挡网站| 日韩大片免费观看网站| 热99久久久久精品小说推荐| 日日摸夜夜添夜夜爱| 免费观看在线日韩| 麻豆精品久久久久久蜜桃| 久久久久精品久久久久真实原创| 久久久久久久亚洲中文字幕| 在线观看国产h片| 女人久久www免费人成看片| 99久久精品国产国产毛片| 亚洲精品日韩在线中文字幕| 欧美日韩亚洲高清精品| 国产成人免费观看mmmm| 国产午夜精品一二区理论片| 黄色毛片三级朝国网站| 大香蕉久久网| 91在线精品国自产拍蜜月| 精品酒店卫生间| 2021少妇久久久久久久久久久| 在线天堂最新版资源| 美女福利国产在线| 精品熟女少妇av免费看| 久久久精品区二区三区| 久久久久久久久大av| 一区二区日韩欧美中文字幕 | 国内精品宾馆在线| 日本与韩国留学比较| 久久久久久久久久久丰满| 午夜激情久久久久久久| 不卡视频在线观看欧美| 免费观看的影片在线观看| 汤姆久久久久久久影院中文字幕| 亚洲图色成人| 欧美少妇被猛烈插入视频| 岛国毛片在线播放| 一本一本综合久久| 成人黄色视频免费在线看| 新久久久久国产一级毛片| 久久久精品94久久精品| 国产欧美亚洲国产| 中文乱码字字幕精品一区二区三区| 欧美精品一区二区免费开放| 一个人免费看片子| 天堂8中文在线网| 全区人妻精品视频| 国产精品一二三区在线看| 国产伦精品一区二区三区视频9| av专区在线播放| 亚洲欧美日韩卡通动漫| 纵有疾风起免费观看全集完整版| 欧美日韩亚洲高清精品| 久久人人爽人人爽人人片va| 久久久久国产网址| 男人爽女人下面视频在线观看| 国产午夜精品久久久久久一区二区三区| 久久精品国产a三级三级三级| 多毛熟女@视频| 国产精品.久久久| 国产 精品1| 日韩中字成人| h视频一区二区三区| 一二三四中文在线观看免费高清| 边亲边吃奶的免费视频| 国产爽快片一区二区三区| 一边亲一边摸免费视频| 看免费成人av毛片| 亚洲欧美中文字幕日韩二区| 国产国语露脸激情在线看| 午夜免费观看性视频| 美女福利国产在线| 亚洲精品视频女| 日日啪夜夜爽| 99九九线精品视频在线观看视频| 日本免费在线观看一区| 久久精品久久久久久久性| 亚洲精品乱久久久久久| 免费看光身美女| 男人爽女人下面视频在线观看| 国产一区二区在线观看日韩| 日韩一区二区视频免费看| 国产成人免费无遮挡视频| 国产精品 国内视频| 亚洲图色成人| 久久影院123| 日产精品乱码卡一卡2卡三| 色婷婷久久久亚洲欧美| 伊人久久精品亚洲午夜| 少妇猛男粗大的猛烈进出视频| 成年人免费黄色播放视频| 视频中文字幕在线观看| 在线播放无遮挡| 久久人人爽av亚洲精品天堂| 成年人免费黄色播放视频| 好男人视频免费观看在线| 国产精品久久久久久久电影| 热99国产精品久久久久久7| 母亲3免费完整高清在线观看 | 在线观看免费视频网站a站| 狂野欧美激情性xxxx在线观看| 亚洲,欧美,日韩| 高清不卡的av网站| 久久久久久久久久人人人人人人| 最近手机中文字幕大全| 边亲边吃奶的免费视频| 久久人人爽人人片av| 欧美人与性动交α欧美精品济南到 | 亚洲人成网站在线观看播放| 国产精品女同一区二区软件| 看十八女毛片水多多多| 欧美 亚洲 国产 日韩一| 母亲3免费完整高清在线观看 | 18禁在线播放成人免费| 国产成人精品久久久久久| 精品99又大又爽又粗少妇毛片| 亚洲精品国产色婷婷电影| 日本免费在线观看一区| 亚洲av不卡在线观看| 91久久精品国产一区二区三区| 国产有黄有色有爽视频| 91精品国产九色| 一级毛片电影观看| 日本免费在线观看一区| 午夜精品国产一区二区电影| 成人影院久久| 亚洲av日韩在线播放| 亚洲不卡免费看| 久久综合国产亚洲精品| 国产免费一级a男人的天堂| 老熟女久久久| 综合色丁香网| 久久久亚洲精品成人影院| 国产亚洲av片在线观看秒播厂| 免费观看性生交大片5| 亚洲精品视频女| 在线观看一区二区三区激情| 在线观看免费高清a一片| 亚洲精品,欧美精品| 日韩电影二区| 亚洲婷婷狠狠爱综合网| 午夜av观看不卡| 成年女人在线观看亚洲视频| xxxhd国产人妻xxx| 欧美日韩一区二区视频在线观看视频在线| 国产片特级美女逼逼视频| 国产成人午夜福利电影在线观看| 成人国产麻豆网| 亚洲欧美色中文字幕在线|