• <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
    国产成人a∨麻豆精品| 日韩一本色道免费dvd| 久99久视频精品免费| 日韩人妻高清精品专区| 三级男女做爰猛烈吃奶摸视频| 九九久久精品国产亚洲av麻豆| 国产黄片美女视频| 国产大屁股一区二区在线视频| 免费观看性生交大片5| 亚洲成人中文字幕在线播放| av在线蜜桃| 伦理电影大哥的女人| 午夜精品一区二区三区免费看| 最近2019中文字幕mv第一页| 国产精品日韩av在线免费观看| 99热网站在线观看| 九九久久精品国产亚洲av麻豆| 国产伦精品一区二区三区视频9| 日韩欧美精品v在线| 亚洲乱码一区二区免费版| 亚洲欧美日韩无卡精品| 中文字幕av成人在线电影| 久久久久国产网址| 一级毛片黄色毛片免费观看视频| 中国美白少妇内射xxxbb| 成人美女网站在线观看视频| 亚洲成人一二三区av| 国内精品宾馆在线| 亚州av有码| 国产成人精品久久久久久| 波多野结衣巨乳人妻| 欧美成人精品欧美一级黄| 国模一区二区三区四区视频| 久久人人爽人人爽人人片va| 久久国产乱子免费精品| 熟女电影av网| 男女那种视频在线观看| 91精品国产九色| 干丝袜人妻中文字幕| 亚洲天堂国产精品一区在线| 亚洲不卡免费看| 欧美成人精品欧美一级黄| 婷婷色av中文字幕| 国产精品国产三级国产专区5o| 亚洲四区av| 麻豆成人午夜福利视频| 亚洲欧美中文字幕日韩二区| 麻豆成人午夜福利视频| 亚洲精品国产av成人精品| 国产精品久久视频播放| 国产成人a∨麻豆精品| 久久6这里有精品| 日韩av在线大香蕉| 免费看美女性在线毛片视频| 免费黄色在线免费观看| 18禁在线无遮挡免费观看视频| 亚洲国产欧美在线一区| 国产精品av视频在线免费观看| 日韩欧美 国产精品| 免费看光身美女| 国产亚洲精品av在线| 欧美日韩亚洲高清精品| 狂野欧美白嫩少妇大欣赏| 国产精品精品国产色婷婷| 亚洲国产精品sss在线观看| 国产成人一区二区在线| 国产69精品久久久久777片| 我要看日韩黄色一级片| 欧美日韩精品成人综合77777| 熟女电影av网| 直男gayav资源| 欧美成人a在线观看| 国产一区亚洲一区在线观看| 少妇丰满av| 亚洲av男天堂| 国产一区二区亚洲精品在线观看| 午夜福利在线观看吧| 亚洲成人久久爱视频| 欧美xxxx黑人xx丫x性爽| 最近手机中文字幕大全| 91在线精品国自产拍蜜月| 国内揄拍国产精品人妻在线| 欧美三级亚洲精品| 99re6热这里在线精品视频| 尤物成人国产欧美一区二区三区| 亚洲欧美中文字幕日韩二区| 免费看美女性在线毛片视频| 国产 亚洲一区二区三区 | 欧美3d第一页| 美女内射精品一级片tv| 2018国产大陆天天弄谢| 午夜日本视频在线| 人体艺术视频欧美日本| .国产精品久久| 午夜日本视频在线| 久久人人爽人人爽人人片va| 在现免费观看毛片| 欧美一级a爱片免费观看看| 亚洲色图av天堂| 久久久午夜欧美精品| 亚洲国产精品专区欧美| 一个人免费在线观看电影| 国产精品一区www在线观看| 日本三级黄在线观看| 国产伦在线观看视频一区| av.在线天堂| av免费观看日本| 亚洲欧洲日产国产| 午夜福利视频1000在线观看| 女的被弄到高潮叫床怎么办| 搞女人的毛片| 国产免费一级a男人的天堂| 精品久久国产蜜桃| 亚洲欧美中文字幕日韩二区| 国产一区有黄有色的免费视频 | 少妇人妻精品综合一区二区| 亚洲三级黄色毛片| 插逼视频在线观看| 国产精品av视频在线免费观看| 日本猛色少妇xxxxx猛交久久| 天天躁日日操中文字幕| 欧美3d第一页| 真实男女啪啪啪动态图| 精品国产一区二区三区久久久樱花 | 免费看av在线观看网站| 国产精品麻豆人妻色哟哟久久 | 亚洲精品久久久久久婷婷小说| 免费不卡的大黄色大毛片视频在线观看 | 免费观看性生交大片5| 亚洲av免费高清在线观看| 成年人午夜在线观看视频 | 丝袜喷水一区| 亚洲最大成人中文| 免费观看在线日韩| 精品一区二区免费观看| 日本与韩国留学比较| 天堂√8在线中文| 69av精品久久久久久| 能在线免费观看的黄片| 国产日韩欧美在线精品| 成人亚洲精品一区在线观看 | 久久久久久国产a免费观看| 国语对白做爰xxxⅹ性视频网站| 91精品国产九色| 在线播放无遮挡| 国产免费一级a男人的天堂| 色吧在线观看| av卡一久久| 免费看日本二区| 中文字幕av在线有码专区| 简卡轻食公司| 极品少妇高潮喷水抽搐| 国产精品1区2区在线观看.| 精品人妻偷拍中文字幕| 肉色欧美久久久久久久蜜桃 | 免费黄频网站在线观看国产| 狠狠精品人妻久久久久久综合| 熟妇人妻不卡中文字幕| 亚洲人与动物交配视频| 国产老妇伦熟女老妇高清| 秋霞在线观看毛片| 一级毛片aaaaaa免费看小| 国产熟女欧美一区二区| 亚洲伊人久久精品综合| av国产久精品久网站免费入址| 亚洲av不卡在线观看| 国产免费又黄又爽又色| 欧美精品一区二区大全| 午夜激情久久久久久久| 简卡轻食公司| 国产淫语在线视频| 91久久精品电影网| 高清av免费在线| 亚洲国产av新网站| 国产精品久久视频播放| 国产亚洲av片在线观看秒播厂 | 又黄又爽又刺激的免费视频.| 国产探花极品一区二区| 黄片无遮挡物在线观看| 国产精品精品国产色婷婷| 亚洲精品影视一区二区三区av| 两个人视频免费观看高清| 亚洲欧美成人精品一区二区| 插阴视频在线观看视频| 国产亚洲精品久久久com| 国产精品国产三级国产av玫瑰| 久久久成人免费电影| 最新中文字幕久久久久| 亚洲自拍偷在线| 亚洲精品,欧美精品| 久久久午夜欧美精品| 男女啪啪激烈高潮av片| 亚洲精品aⅴ在线观看| 欧美日韩一区二区视频在线观看视频在线 | 色播亚洲综合网| 狠狠精品人妻久久久久久综合| 免费大片黄手机在线观看| 国产黄频视频在线观看| 亚洲av电影在线观看一区二区三区 | 中文精品一卡2卡3卡4更新| 亚洲成人一二三区av| 精品熟女少妇av免费看| 亚洲欧美日韩卡通动漫| 可以在线观看毛片的网站| 哪个播放器可以免费观看大片| 午夜免费激情av| 18+在线观看网站| 成人国产麻豆网| 亚洲国产色片| 免费无遮挡裸体视频| 久久久国产一区二区| 国产精品伦人一区二区| 91午夜精品亚洲一区二区三区| 毛片一级片免费看久久久久| 亚洲av中文字字幕乱码综合| 色吧在线观看| 午夜精品在线福利| 国产精品久久久久久久久免| 精品久久久久久久人妻蜜臀av| 亚洲精品日韩在线中文字幕| 91在线精品国自产拍蜜月| 精品久久久久久成人av| 国产老妇女一区| 日韩欧美三级三区| 亚洲精品第二区| 国产又色又爽无遮挡免| 久久久午夜欧美精品| 深爱激情五月婷婷| 97精品久久久久久久久久精品| 18+在线观看网站| 免费看a级黄色片| 久久综合国产亚洲精品| 伊人久久国产一区二区| 三级国产精品欧美在线观看| 97在线视频观看| 午夜激情福利司机影院| 亚洲人成网站在线观看播放| 国产精品国产三级专区第一集| 综合色丁香网| 亚洲人与动物交配视频| av在线亚洲专区| 久久久午夜欧美精品| 国产成人精品婷婷| 欧美zozozo另类| 久久6这里有精品| 国产黄色小视频在线观看| 国产一区亚洲一区在线观看| 久久这里只有精品中国| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 纵有疾风起免费观看全集完整版 | 亚洲av免费高清在线观看| 国产黄色小视频在线观看| 舔av片在线| 欧美性猛交╳xxx乱大交人| 一级毛片我不卡| 国精品久久久久久国模美| 日日摸夜夜添夜夜爱| 久久精品久久精品一区二区三区| 18禁在线无遮挡免费观看视频| 日韩av在线大香蕉| 欧美成人a在线观看| 国国产精品蜜臀av免费| 有码 亚洲区| 我要看日韩黄色一级片| 欧美一区二区亚洲| 少妇裸体淫交视频免费看高清| av在线播放精品| 一本久久精品| 中文乱码字字幕精品一区二区三区 | 国产黄色免费在线视频| 纵有疾风起免费观看全集完整版 | 精品少妇黑人巨大在线播放| 欧美日本视频| 国产探花在线观看一区二区| 国产精品麻豆人妻色哟哟久久 | 好男人在线观看高清免费视频| 国产在视频线精品| 亚洲高清免费不卡视频| 丝袜美腿在线中文| 毛片女人毛片| 一级黄片播放器| 少妇熟女欧美另类| 韩国高清视频一区二区三区| 精品久久久久久久末码| 国产精品久久久久久av不卡| 久久精品国产亚洲网站| 欧美高清性xxxxhd video| 欧美变态另类bdsm刘玥| 日韩欧美国产在线观看| 国产午夜精品一二区理论片| 纵有疾风起免费观看全集完整版 | 日韩av不卡免费在线播放| 欧美性猛交╳xxx乱大交人| 日日干狠狠操夜夜爽| 一级毛片我不卡| 国产麻豆成人av免费视频| 免费看av在线观看网站| 美女cb高潮喷水在线观看| 国产精品久久久久久久久免| 欧美潮喷喷水| 国产熟女欧美一区二区| 欧美xxxx黑人xx丫x性爽| 熟女人妻精品中文字幕| 亚洲av成人精品一区久久| 久久久久国产网址| 精品熟女少妇av免费看| 美女主播在线视频| 美女xxoo啪啪120秒动态图| av女优亚洲男人天堂| 日韩欧美国产在线观看| 九草在线视频观看| 自拍偷自拍亚洲精品老妇| 18禁在线播放成人免费| 免费黄网站久久成人精品| 亚洲精华国产精华液的使用体验| 日本熟妇午夜| 亚洲精品视频女| 欧美精品国产亚洲| 大片免费播放器 马上看| 久久精品久久久久久久性| 亚洲丝袜综合中文字幕| h日本视频在线播放| 日韩一本色道免费dvd| 免费电影在线观看免费观看| 免费大片18禁| 国产成人aa在线观看| 十八禁国产超污无遮挡网站| h日本视频在线播放| 国产三级在线视频| 国产成人福利小说| 日韩av在线大香蕉| 秋霞在线观看毛片| 午夜久久久久精精品| 91精品一卡2卡3卡4卡| 精品久久久久久久人妻蜜臀av| 日韩av在线免费看完整版不卡| 国产精品爽爽va在线观看网站| a级毛色黄片| 91狼人影院| 日韩精品有码人妻一区| 一区二区三区免费毛片| 亚洲精品一区蜜桃| 国产男女超爽视频在线观看| 色综合亚洲欧美另类图片| 91在线精品国自产拍蜜月| 国产综合懂色| 高清午夜精品一区二区三区| 国产精品一区二区三区四区久久| 免费在线观看成人毛片| 看免费成人av毛片| 国产成人免费观看mmmm| 国产黄色视频一区二区在线观看| 国产黄片美女视频| 国产黄片视频在线免费观看| 国产成人精品福利久久| 一二三四中文在线观看免费高清| 精品久久久久久久人妻蜜臀av| 神马国产精品三级电影在线观看| 26uuu在线亚洲综合色| 亚洲自拍偷在线| 午夜福利视频精品| 亚洲不卡免费看| 精品国产一区二区三区久久久樱花 | 黄色配什么色好看| 九九久久精品国产亚洲av麻豆| 日本爱情动作片www.在线观看| 丝袜喷水一区| 久久久午夜欧美精品| 天堂√8在线中文| 99热6这里只有精品| 成人一区二区视频在线观看| 日本午夜av视频| 成年人午夜在线观看视频 | 亚洲成色77777| 亚洲精品乱久久久久久| 日日摸夜夜添夜夜爱| 亚洲人成网站在线播| 永久网站在线| 亚洲精品色激情综合| 免费在线观看成人毛片| 寂寞人妻少妇视频99o| 日本av手机在线免费观看| 亚洲人成网站在线播| 亚洲国产高清在线一区二区三| 久久综合国产亚洲精品| 日韩欧美国产在线观看| .国产精品久久| 草草在线视频免费看| 国产精品久久久久久精品电影小说 | 少妇猛男粗大的猛烈进出视频 | 国产白丝娇喘喷水9色精品| 视频中文字幕在线观看| 少妇的逼水好多| 美女黄网站色视频| 欧美日韩亚洲高清精品| 亚洲人成网站在线观看播放| av免费在线看不卡| 国产在线一区二区三区精| 3wmmmm亚洲av在线观看| 久久精品久久久久久久性| 成人毛片60女人毛片免费| 日韩三级伦理在线观看| 51国产日韩欧美| 国产精品1区2区在线观看.| 久久99精品国语久久久| 日本欧美国产在线视频| 网址你懂的国产日韩在线| 亚洲内射少妇av| 大又大粗又爽又黄少妇毛片口| 亚洲18禁久久av| 又爽又黄a免费视频| 人人妻人人澡欧美一区二区| 男人爽女人下面视频在线观看| 国产白丝娇喘喷水9色精品| 免费不卡的大黄色大毛片视频在线观看 | 日本爱情动作片www.在线观看| 国产精品一二三区在线看| 97超碰精品成人国产| 国产精品av视频在线免费观看| 国产精品一区www在线观看| 久久久久国产网址| 国产伦精品一区二区三区视频9| 国产精品美女特级片免费视频播放器| 日日摸夜夜添夜夜添av毛片| 精品一区二区免费观看| 禁无遮挡网站| 尾随美女入室| 国产一区亚洲一区在线观看| 国产精品久久久久久精品电影小说 | 欧美三级亚洲精品| 男女视频在线观看网站免费| 国产男女超爽视频在线观看| 亚洲综合精品二区| 欧美区成人在线视频| 日韩亚洲欧美综合| 亚洲av中文av极速乱| 建设人人有责人人尽责人人享有的 | 亚洲不卡免费看| 最近最新中文字幕免费大全7| 国产一区有黄有色的免费视频 | 国产精品麻豆人妻色哟哟久久 | av.在线天堂| 美女内射精品一级片tv| 国产熟女欧美一区二区| 人妻少妇偷人精品九色| 少妇人妻一区二区三区视频| 亚洲精品乱码久久久v下载方式| 免费av毛片视频| 免费观看精品视频网站| 免费大片18禁| 精品不卡国产一区二区三区| 男人爽女人下面视频在线观看| 日产精品乱码卡一卡2卡三| 如何舔出高潮| 国产精品一及| 亚洲av成人精品一区久久| 国产亚洲午夜精品一区二区久久 | 91aial.com中文字幕在线观看| 女的被弄到高潮叫床怎么办| 免费在线观看成人毛片| 成人毛片a级毛片在线播放| 插阴视频在线观看视频| 少妇高潮的动态图| 日本av手机在线免费观看| 亚洲国产日韩欧美精品在线观看| 夫妻午夜视频| 国产精品人妻久久久久久| 久久97久久精品| 久久久久久久久久人人人人人人| 精品一区二区三区视频在线| 亚洲欧美日韩无卡精品| 久久久久久久久久久丰满| 国产精品国产三级国产av玫瑰| 80岁老熟妇乱子伦牲交| 日本爱情动作片www.在线观看| 成人亚洲欧美一区二区av| 久久久久网色| 美女大奶头视频| 亚洲美女搞黄在线观看| 天堂俺去俺来也www色官网 | 日本黄色片子视频| 免费av毛片视频| 国产伦精品一区二区三区四那| 一级毛片久久久久久久久女| 插逼视频在线观看| 国产麻豆成人av免费视频| 婷婷色综合大香蕉| 国产免费又黄又爽又色| 在线观看美女被高潮喷水网站| 三级经典国产精品| 久99久视频精品免费| 亚洲乱码一区二区免费版| 亚洲美女搞黄在线观看| 男女边吃奶边做爰视频| 国产毛片a区久久久久| 亚洲av免费高清在线观看| 免费看a级黄色片| 高清毛片免费看| 国产成人福利小说| 国产精品一区二区三区四区免费观看| 亚洲精品日本国产第一区| 黄片无遮挡物在线观看| 极品少妇高潮喷水抽搐| 97在线视频观看| 少妇熟女aⅴ在线视频| 国语对白做爰xxxⅹ性视频网站| 亚洲激情五月婷婷啪啪| 丝瓜视频免费看黄片| 一级a做视频免费观看| 熟妇人妻久久中文字幕3abv| 18禁裸乳无遮挡免费网站照片| 午夜日本视频在线| 国产在线一区二区三区精| 久久久精品94久久精品| 人人妻人人澡人人爽人人夜夜 | 国产精品一区二区在线观看99 | 午夜福利在线观看免费完整高清在| 免费黄色在线免费观看| 最新中文字幕久久久久| 精品久久久噜噜| 床上黄色一级片| 美女被艹到高潮喷水动态| 九色成人免费人妻av| 大又大粗又爽又黄少妇毛片口| 日本猛色少妇xxxxx猛交久久| 免费人成在线观看视频色| 亚洲欧洲日产国产| 欧美变态另类bdsm刘玥| 综合色av麻豆| 中国美白少妇内射xxxbb| 乱人视频在线观看| 国产伦一二天堂av在线观看| 久久精品久久精品一区二区三区| 欧美不卡视频在线免费观看| 欧美成人一区二区免费高清观看| 夜夜看夜夜爽夜夜摸| 在现免费观看毛片| 边亲边吃奶的免费视频| 日韩欧美精品免费久久| 永久网站在线| 国产美女午夜福利| 久久99蜜桃精品久久| 国产69精品久久久久777片| 91aial.com中文字幕在线观看| 欧美成人精品欧美一级黄| 成人国产麻豆网| 亚洲天堂国产精品一区在线| 国产亚洲91精品色在线| 国产午夜精品一二区理论片| 听说在线观看完整版免费高清| 国产探花极品一区二区| 亚洲精品影视一区二区三区av| 白带黄色成豆腐渣| 成人无遮挡网站| 肉色欧美久久久久久久蜜桃 | 美女高潮的动态| 美女主播在线视频| 高清在线视频一区二区三区| 国产精品一区二区三区四区久久| 精品久久久久久久人妻蜜臀av| 99九九线精品视频在线观看视频| 岛国毛片在线播放| 亚洲va在线va天堂va国产| 午夜福利在线观看免费完整高清在| 国产人妻一区二区三区在| 男女边摸边吃奶| 韩国高清视频一区二区三区| 男女下面进入的视频免费午夜| 中文在线观看免费www的网站| 国产精品精品国产色婷婷| 亚洲电影在线观看av| 人妻系列 视频| 美女国产视频在线观看| 国产伦一二天堂av在线观看| 老司机影院成人| 最新中文字幕久久久久| 国产成人午夜福利电影在线观看| 欧美成人a在线观看| 亚洲激情五月婷婷啪啪| 乱码一卡2卡4卡精品| 欧美性猛交╳xxx乱大交人| 精品熟女少妇av免费看| 一级av片app| 又粗又硬又长又爽又黄的视频| 久久鲁丝午夜福利片| 水蜜桃什么品种好| 日日撸夜夜添| 精品少妇黑人巨大在线播放| 在线 av 中文字幕| 婷婷六月久久综合丁香| 国产又色又爽无遮挡免| 成人无遮挡网站| 精品人妻视频免费看| 啦啦啦韩国在线观看视频| 在线观看免费高清a一片| 国产极品天堂在线| 十八禁国产超污无遮挡网站| 欧美一区二区亚洲| 寂寞人妻少妇视频99o| 亚洲美女搞黄在线观看| 99九九线精品视频在线观看视频| 国产精品一区二区性色av| 亚洲美女搞黄在线观看| 18禁在线播放成人免费| 纵有疾风起免费观看全集完整版 | 97人妻精品一区二区三区麻豆| 亚洲三级黄色毛片| 在线观看免费高清a一片| 久久99精品国语久久久| 国产免费又黄又爽又色| 美女内射精品一级片tv| 国产一区二区三区av在线| 免费人成在线观看视频色| 亚洲国产最新在线播放| 日韩人妻高清精品专区|