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

    航空磁探中潛艇磁場建模方法分析

    2017-03-08 11:12:37周家新陳建勇單志超陳長康
    海軍航空大學學報 2017年1期
    關鍵詞:磁偶極子模擬法磁體

    周家新,陳建勇,單志超,陳長康

    (海軍航空工程學院a.研究生管理大隊;b.電子信息工程系,山東煙臺264001)

    航空磁探中潛艇磁場建模方法分析

    周家新a,陳建勇b,單志超b,陳長康a

    (海軍航空工程學院a.研究生管理大隊;b.電子信息工程系,山東煙臺264001)

    為提高航空磁探中潛艇磁場模型精度,分析現(xiàn)有潛艇磁場建模方法,對邊界積分法、有限元法、積分方程法、磁體模擬法等主要潛艇磁場建模方法的基本原理和推導過程進行分析。通過數(shù)值分析比較各方法的建模精度,最后歸納潛艇低空磁場建模方法的基本要點。

    航空磁探;潛艇磁場;邊界積分法;有限元法;積分方程法;磁體模擬法

    基于磁性探潛應用的潛艇磁場建模理論和技術研究始于一戰(zhàn)時期[1]。新倫敦海軍試驗站在1918年使用固定線圈型磁檢測器檢測磁性目標,并對潛艇的磁性進行研究。1944年,搭載Madcats航空磁探儀的63巡邏中隊在直布羅陀海峽準確追蹤到U-761潛艇[2]。20世紀90年代國外已開始使用FLUX3D、ANSYS等有限元分析軟件對潛艇進行磁場建模[3],并得到高模擬精度的數(shù)學模型。

    從PB報告中了解到美國早已對艦艇磁特性、艦艇低空磁場特征深入探索。日本最先采用磁偶極子近似推算潛艇的準靜態(tài)磁場低空磁異常信號[4]。德、英、法、俄等國通過磁性鋼板、薄壁艦艇殼層等[5-6]磁性模型,對艦艇磁場有限元建模分析方法進行了大量研究[7],用于艦艇感應磁場分析、渦流磁場分析和消磁繞組優(yōu)化設計等方面。

    我國有關艦艇磁場和航空磁性探潛的研究起始于20世紀50年代末,由于多種因素,發(fā)展并不順利。但經過多年研究,也取得了一些潛艇磁場建模的理論和技術成果[8-9]。潛艇空間磁場是航空磁探儀對潛探測和磁性武器對潛攻擊的探測源[10],其空間分布直接影響航空磁探儀對潛艇探測寬度和作用距離以及對磁性武器的磁隱身能力[11]。潛艇磁場建??煞譃槲锢斫:蛿?shù)學建模,與傳統(tǒng)的潛艇物理磁模型相比較,數(shù)學模型分析性價比更高,耗費更小,研制周期更短,效率更高。

    數(shù)學模型仿真方法的基本思路是構建潛艇磁場的延拓數(shù)學模型,然后根據(jù)部分測量數(shù)據(jù)作為延拓數(shù)學模型的輸入對其他空間的磁場進行換算,主要包括磁體模擬(又稱“等效源”)法[12-13]、積分方程法[14]、有限元(素)法[7]、邊界積分法[15]。

    圖1為常用的模型方法。

    1 邊界積分法

    邊界積分法在半空間區(qū)域基于格林函數(shù),通過標量磁位分布求拉普拉斯的邊值問題。根據(jù)麥克斯韋方程組可以得到潛艇靜磁場描述如下[16]:

    式(1)中:H為潛艇磁場的磁場強度;δ為電流密度;B為磁感應強度;μ0=4π×10-7H/m為真空磁導率。

    由潛艇附近的電流密度δ=0,可得

    式(3)中:Hx、Hy、Hz分別是磁場強度H在x、y、z軸向上的分量的幅值。

    圖2為潛艇磁場場域的示意圖。場源代表磁性目標潛艇,Γs為磁場數(shù)據(jù)觀測面,Q為觀測面以上的計算點,Ω為無源封閉區(qū)域,P為位于觀測面Γs上的流動積分點,?!逓闊o窮遠邊界曲面,r為計算點Q與流動點P之間的距離。場域為有限空間,觀測面為有限平面。

    在圖2所示的無源封閉區(qū)域中,代入標量磁位u,由靜磁場的麥克斯韋方程得到

    由??B=μ0??H=-μ0?2u=0,得到拉普拉斯方程

    根據(jù)格林函數(shù),解式(3)的拉普拉斯方程,則空間場點Q(x′,y′,z′)處的三分量磁場強度為:

    通過觀測面Γs上的實際測量數(shù)據(jù),再根據(jù)式(6)的潛艇磁場向上延拓公式便可得到潛艇磁場強度。

    邊界積分法對求解區(qū)域邊界進行單元化剖分,可對計算問題進行降維,加快模型求解速度;但其所需要的實際數(shù)據(jù)量大,觀測面必須滿足大平面要求。

    2 有限元法

    有限元法當前廣泛應用在工程實踐當中,在潛艇磁場建模中對潛艇幾何形狀進行元素化離散剖分,得到單元區(qū)域,進行偏微分數(shù)值求解。將潛艇磁場建模問題轉化成的泛函變分問題,離散求解該變分問題可以得到潛艇磁場的有限元數(shù)值解。

    將式(5)等效變換為潛艇高空磁場數(shù)學模型的泛函變分問題:

    式(7)中:Ω1為空間場域;μ為磁導率,對應不同的鐵磁性物質以及不同區(qū)域,取值不同。

    潛艇具有一定的磁矩,其磁場在無窮極限處近似衰減為零,即當|z|→∞時,得到無窮遠的邊界條件,邊界處的磁位為uB=-Hz0?z。

    將半空間場域Ω1離散為N個微元,各微元的體積為Vn(n=1,2,…,N),并形成m個節(jié)點,則磁位u可以使用單元上的節(jié)點磁位ui的三維差值函數(shù)Ni為:

    對式(7)進行離散化,并將式(8)代入,得到實際可計算的潛艇低空磁場數(shù)學模型的泛函變分為

    則變分問題可以變?yōu)?/p>

    當泛函F取極小值時,可以得到

    根據(jù)半空間無窮邊界面上的邊值條件,可得到各個剖分單元節(jié)點上的潛艇低空磁位ui。

    將磁位u代入式(4),即可得到潛艇磁場空間分布。

    反潛機一般飛行在一定高度上,場域空間極大,對場域空間剖分受到條件限制。

    當剖分的微元個數(shù)極大時,潛艇磁場模型求解必須有足夠的計算速度支持。

    3 積分方程法

    潛艇在空間場點Q(x′,y′,z′)處產生矢徑為r,磁場強度為H的磁場,且

    Hs為潛艇附近的空間電流因電磁感應在空間場點Q(x′,y′,z′)處產生的磁場,且

    式(14)中:j是海水中的電流密度;r′是海水中的電流元 jdv于空間場點處的矢徑。

    Hm是潛艇作為鐵磁物質磁化產生的磁場,且

    對式(13)進行離散化,得到:

    將式(16)寫成矩陣形式:

    空間場點Q(x′,y′,z′)處的磁場分量為:

    求解式(17)矩陣,得到磁化率χ′,聯(lián)立式(18)則可以得到潛艇磁場在空間中的分布。

    積分方程法依據(jù)嚴格的理論推導出來,其計算精度高,但需要對艦艇進行精確的剖分,計算量較大。

    4 磁體模擬法

    磁體模擬法又稱等效源法,用已知磁場分布的磁性物體來等效實際潛艇,模擬潛艇磁場。常用的磁性目標磁場建模磁體模擬模型有:旋轉橢球體模型、磁偶極子陣列模型、旋轉橢球體與磁偶極子陣列混合模型。裝載航空磁性探潛設備的巡邏機位于潛艇的遠場空間,潛艇低空磁場可使用如圖3所示的旋轉橢球體模型近似等效。圖3中,L=2a是潛艇長度,B=2b是潛艇寬度,半焦距。

    基于單個旋轉橢球體模型近似的潛艇磁場空間分布擬合精度不高,通過磁偶極子陣列模型擬合潛艇磁場,得到精度較高的模型。

    圖4所示為潛艇的磁偶極子陣列模型,通過一定數(shù)量的磁偶極子形成的矢量場疊加,擬合空間磁場分布。當磁偶極子個數(shù)足夠多時,預測磁場值將趨近真實值。

    針對航空磁探潛需求,為提高潛艇低空磁場擬合精度和求解速度,常使用圖5所示的旋轉橢球體與磁偶極子陣列混合模型。

    潛艇在空間場點Qj(x′j,y′j,z′j)處的磁場分量為Hx′j、Hy′j、Hz′j,可以得到:

    式(19)中:a、b、c根據(jù)磁體不同可查閱文獻[9]獲得;Mx′i、My′i、Mz′i分別為磁偶極子沿x軸、y軸、z軸方向的磁矩分量帶入觀測平面的實際測量磁場數(shù)據(jù)和坐標信息。

    當i=0時,式(19)表示旋轉橢球體模型;當i=1,2,…,n時,式(19)表示磁偶極子陣列模型;當i=0,1,2,…,n時,式(19)表示旋轉橢球體與磁偶極子陣列混合模型;j=1,2,…,n,表示觀測平面上的第j個實測磁場數(shù)據(jù)。

    式(19)的通用矩陣表示式為:

    基于磁體模擬思想的磁場建模方法能夠在低空遠場中擬合潛艇磁場,通過逐步回歸法、遺傳算法[17]、微粒群(PSO)算法[18]等方法對式(20)進行優(yōu)化求解,得到潛艇磁源參數(shù),再通過正演延拓便可得到低空磁場分布。磁體模擬法在計算量和實測磁場數(shù)據(jù)需求量等方面有突出優(yōu)勢。但使用這些優(yōu)化算法,通常計算時間長,可靠性也存在較大的問題。

    5 數(shù)值分析

    為驗證各方法建立的磁場模型精度,通過實測數(shù)據(jù)與預測值進行比較。使用一個雙層殼體結構的鋼管來模擬簡易的潛艇模型,進行數(shù)值分析。內鋼管的直徑為0.30 m,外鋼管直徑為0.46 m,長度為1 m,各層厚度均為0.002 m。模型指向磁北向,與x軸重合,其中x∈[-20 m,20 m]。以模型中心正上方50 m處為x軸的原點,在[-20 m,20 m]區(qū)間內每2 m進行一組測量。通過光泵測磁儀測量得到測量值與模型基數(shù)按得到的預測值進行比較,結果如圖6所示。

    定義最大相對誤差為:

    可得邊界積分法、有限元法、積分方程法、磁體模擬法各模型的最大相對誤差分別為15.3%、10.4%、7.6%、22.3%,其中磁體模擬法選用單個磁偶極子模型。

    6 結語

    航空磁異常探潛中,由于飛機距離潛艇的距離一般大于潛艇的線度,常應用簡單磁偶極子模型。但該模型沒有充分利用潛艇磁異常信號的全部信息,對潛艇磁場的有限擬合程度導致丟失部分更加精確的潛艇磁異常信號特征,誤差較大。本文通過對各種建模方法的比較發(fā)現(xiàn),積分方程法確定的磁場模型精度較高,邊界元法和有限元法的擬合精度受限于空間場域大小,磁體模擬法受限于計算量。潛艇磁場數(shù)學模型的選擇要求在低空遠場區(qū)域內的擬合精度高、實時性好、解算所需的磁場實測數(shù)據(jù)量少,并且簡單易行,計算量小。未來研究將基于積分方程法,進一步提高潛艇磁場模型的實時性等性能。

    [1]GEORGE IOANNIDIS.Identification of a ship or submarine from its magnetic signature[J].IEEE Transactions on Aerospace and Electronic Systems,1977,13(3):327-329.

    [2]JEFFREY MARK SCHWEIGER.Evaluation of geomagnetic activity in the mad frequency band[D].California:Naval Postgraduate School,1982.

    [3]XAVIER BRUNOTTE,GERARD MEUNIER,JEAN PAUL BONGIRAUD.Ship magnetizations modeling by the finite element method[J].IEEE Transactions on Magnetics,1993,29(2):1970-1975.

    [4]KREVET B.Nonlinear magnetic field calculation using dipole approximation[J].IEEE Transactions on Magnetics,1992,28(2):1060-1063.

    [5]OLIVIER CHADEBEC,JEAN LOUIS COULOMB,JEAN PAUL BONGIRAUD.Recent improvements for solving inverse magnetostatic problem applied to thin shells[J].IEEE Transactions on Magnetics,2002,38(2):1005-1008.

    [6]LEV FRUMKIS,BEN ZION KAPLAN.Spherical and spheroidal shells as models in magnetic detection[J]. IEEE Transactions on Magnetics,1999,35(5):4151-4158.

    [7]YAHIA EL OUAZZANI,F(xiàn)RANGOISE RIOUX DAMIDAU.Finite element modeling of unbounded problems:use of a geometrical transformation and comparison with the boundary integral method[J].IEEE Transactions on Magnetics,1996,32(3):1401-1404.

    [8]林春生,龔沈光.艦船物理場[M].北京:兵器工業(yè)出版社,2007:45-49. LIN CHUNSHENG,GONG SHENGUANG.Physical field of ship[M].Beijing:Weapon Industry Press,2007:45-49.(in Chinese)

    [9]周耀忠,張國友.艦船磁場分析計算[M].北京:國防工業(yè)出版社,2004:104-200. ZHOU YAOZHONG,ZHANG GUOYOU.Ship magnetic field analyze and compute[M].Beijing:National Defense Industry Press,2004:104-200.(in Chinese)

    [10]PALMSTROM Robert E,WANG GUNNAR.Detection ofmoving magneticdipolesby three-dimensional matched filter techniques[M].US:National Technical Information Service,1996:7-8.

    [11]翁行泰,曹梅芬.磁異探潛中潛艇的數(shù)學模型[J].上海交通大學學報,1995,29(3):27-32. WENG XINGTAI,CAO MEIFEN.Submarine math model of magnetic abnormity detection[J].Journal of Shanghai Jiaotong University,2003,29(3):27-32.(in Chinese)

    [12]ALASTAIR D McAULAY.Computerized model demon-strating magnetic submarine localization[J].IEEE Transactions on Aerospace and Electronic Systems,1977,13(3):246-254.

    [13]林春生.艦船磁場信號檢測與磁性目標定位[D].武漢:海軍工程大學,1996. LIN CHUNSHENG.Signal detection of ship magnetic field and magnetic target location[D].Wuhan:Naval Engineering Institute,1996.(in Chinese)

    [14]郭成豹,肖昌漢,劉大明.基于積分方程法和奇異值分解的磁性目標磁場延拓技術研究[J].物理學報,2008,57(7):4182-4188. GUO CHENGBAO,XIAO CHANGHAN,LIU DAMING.Research on the continuations of magnetic field of magnetic object based on integral equation method and singular value decomposition[J].Acta Physica Sinica,2008,57(7):4182-4188.(in Chinese)

    [15]VISHNEVSKY A,KRASNOV I,LAPOVOK A.Calculation of static magnetization for thin-walled constructions by boundary element method[J].IEEE Transactions on Magnetics,1993,29(3):2152-2155.

    [16]曲曉慧,楊日杰,單志超.潛艇磁場建模方法的分析與比較[J].艦船科學技術,2011,33(3):7-11. QU XIAOHUI,YANG RIJIE,SHAN ZHICHAO.Analysis and comparison on magnetic field modeling method of submarine[J].Ship Science and Technology,2011,33(3):7-11.(in Chinese)

    [17]劉勝道,劉大明,肖昌漢,等.基于遺傳算法的磁性目標磁模型[J].武漢理工大學學報,2008,32(6):1017-1020. LIU SHENGDAO,LIU DAMING,XIAO CHANGHAN,et al.Magnetic model of the magnetic object based on genetic algorithm[J].Journal of Wuhan University of Technology,2008,32(6):1017-1020.(in Chinese)

    [18]張朝陽,肖昌漢,徐杰.基于微粒群優(yōu)化算法的艦船磁模型分析[J].華中科技大學學報,2010,38(11):124-128. ZHANG ZHAOYANG,XIAO CHANGHAN,XU JIE. Analysis of magnetic model for ships using particle swarm optimization method[J].Journal of Huazhong University of Science&Technology,2010,38(11):124-128.(in Chinese)

    Analysis of Submarine Magnetic Field Modeling Method for Aero Magnetic Detection

    ZHOU Jiaxina,CHEN Jianyongb,SHAN Zhichaob,CHEN Changkanga
    (Naval Aeronautical and Astronautical University a.Graduate Students’Brigade; b.Department of Electronic and Information Engineering,Yantai Shandong 264001,China)

    To improve the accuracy of models of submarine magnetic fields in the aero magnetic detection,the existed mod?eling methods of submarine magnetic field was analyzed.Fundamental principles and derived process of main modeling methods of submarine magnetic field was derived,including boundary integral method,finite element method,integral equation method,magnet simulation method.Using numerical analysis,the accuracy of the existed modeling methods was calculated.Finally the basic points of modeling method of submarine magnetic field at low altitude were illustrated.

    aero magnetic detection;submarine magnetic field;boundary integral method;finite element method;integral equation method;magnet simulation method

    O441

    :A

    1673-1522(2017)01-0143-06

    10.7682/j.issn.1673-1522.2017.01.008

    2016-11-07;

    :2017-01-02

    部委預研基金資助項目(302020210)

    周家新(1992-),男,碩士生。

    猜你喜歡
    磁偶極子模擬法磁體
    DyF3熱擴滲技術提高NdFeB磁體矯頑力及其機理分析
    基于遞推更新卡爾曼濾波的磁偶極子目標跟蹤
    可控震源地震勘探中的數(shù)值模擬法應用
    磁偶極子跟蹤的漸進貝葉斯濾波方法
    自動化學報(2017年5期)2017-05-14 06:20:53
    蒙特卡洛模擬法計算電動汽車充電負荷
    含Ce燒結Nd-Fe-B磁體的腐蝕行為
    基于磁偶極子的磁場梯度張量縮并的試驗驗證及相關參數(shù)確定
    傳導冷卻高溫超導儲能磁體制作
    隨機模擬法求不規(guī)則圖形面積
    水平磁偶極子電磁場特征研究
    国产1区2区3区精品| 天堂动漫精品| 听说在线观看完整版免费高清| 男女视频在线观看网站免费 | 可以在线观看的亚洲视频| 国产成人啪精品午夜网站| 制服人妻中文乱码| 国产精品1区2区在线观看.| 久久欧美精品欧美久久欧美| 美女高潮到喷水免费观看| 十八禁人妻一区二区| 午夜成年电影在线免费观看| 岛国在线观看网站| 美女午夜性视频免费| 日韩欧美免费精品| 成熟少妇高潮喷水视频| 国产精品乱码一区二三区的特点| 国产精品香港三级国产av潘金莲| 麻豆av在线久日| 一卡2卡三卡四卡精品乱码亚洲| av电影中文网址| 好男人电影高清在线观看| 欧美午夜高清在线| 大型av网站在线播放| 国产成人欧美在线观看| 色播在线永久视频| 搡老岳熟女国产| 欧美激情高清一区二区三区| 色av中文字幕| 亚洲国产精品成人综合色| 非洲黑人性xxxx精品又粗又长| xxx96com| 一进一出抽搐gif免费好疼| 欧美黑人欧美精品刺激| 一本久久中文字幕| 日韩大尺度精品在线看网址| 丝袜在线中文字幕| 90打野战视频偷拍视频| 精品一区二区三区四区五区乱码| 久久 成人 亚洲| 每晚都被弄得嗷嗷叫到高潮| 淫秽高清视频在线观看| 可以在线观看的亚洲视频| 国产精品电影一区二区三区| 女人高潮潮喷娇喘18禁视频| 久久中文字幕人妻熟女| 美女扒开内裤让男人捅视频| 白带黄色成豆腐渣| 欧美激情 高清一区二区三区| 婷婷精品国产亚洲av| 此物有八面人人有两片| 性色av乱码一区二区三区2| 亚洲免费av在线视频| 国产亚洲精品综合一区在线观看 | 天堂动漫精品| 国产成人av教育| 精品高清国产在线一区| 亚洲国产毛片av蜜桃av| 欧美黑人巨大hd| 欧美中文综合在线视频| 日本黄色视频三级网站网址| 亚洲精品粉嫩美女一区| 欧美丝袜亚洲另类 | 欧美久久黑人一区二区| 一区二区三区高清视频在线| 国产av一区二区精品久久| 极品教师在线免费播放| 精品久久蜜臀av无| 亚洲专区国产一区二区| 日韩三级视频一区二区三区| 在线看三级毛片| 黑人操中国人逼视频| 女人高潮潮喷娇喘18禁视频| 男女视频在线观看网站免费 | 超碰成人久久| 黄色女人牲交| 午夜免费激情av| 91老司机精品| 97碰自拍视频| 丝袜美腿诱惑在线| 久久香蕉激情| 久久久久精品国产欧美久久久| 韩国av一区二区三区四区| 人人妻人人澡欧美一区二区| 久久亚洲真实| svipshipincom国产片| 手机成人av网站| 久热这里只有精品99| 老鸭窝网址在线观看| www.精华液| 亚洲欧美精品综合一区二区三区| 亚洲av成人av| 日本撒尿小便嘘嘘汇集6| 女人爽到高潮嗷嗷叫在线视频| 村上凉子中文字幕在线| 一个人观看的视频www高清免费观看 | 久久精品国产亚洲av香蕉五月| 夜夜爽天天搞| 日韩大尺度精品在线看网址| 在线视频色国产色| 老鸭窝网址在线观看| 最新美女视频免费是黄的| 亚洲精品在线美女| 天天添夜夜摸| 亚洲第一电影网av| 老司机靠b影院| 国产一区二区在线av高清观看| 国产成人av教育| 中文字幕av电影在线播放| 黄网站色视频无遮挡免费观看| 国产精品免费一区二区三区在线| 午夜精品在线福利| 韩国av一区二区三区四区| 亚洲国产精品sss在线观看| 满18在线观看网站| 国产在线精品亚洲第一网站| 波多野结衣高清无吗| 免费女性裸体啪啪无遮挡网站| 成人三级做爰电影| 国产精品98久久久久久宅男小说| 久久久久久久午夜电影| 亚洲国产中文字幕在线视频| 亚洲电影在线观看av| 99在线人妻在线中文字幕| 国产麻豆成人av免费视频| 女人爽到高潮嗷嗷叫在线视频| 国产一区在线观看成人免费| 精品久久久久久久末码| 国产激情偷乱视频一区二区| 久久亚洲精品不卡| 久久九九热精品免费| 日日爽夜夜爽网站| 在线观看免费视频日本深夜| 波多野结衣高清无吗| 男女之事视频高清在线观看| 国内精品久久久久精免费| 国产精品影院久久| 亚洲五月天丁香| 91老司机精品| 成人三级黄色视频| 久久午夜亚洲精品久久| 熟妇人妻久久中文字幕3abv| 黄片小视频在线播放| 精品熟女少妇八av免费久了| 亚洲午夜精品一区,二区,三区| 黄片小视频在线播放| 国产麻豆成人av免费视频| 欧美 亚洲 国产 日韩一| 国内揄拍国产精品人妻在线 | 精品福利观看| 18禁国产床啪视频网站| 久久精品91无色码中文字幕| 精品久久久久久成人av| 成人欧美大片| 国语自产精品视频在线第100页| 久久香蕉国产精品| 国产免费av片在线观看野外av| 一二三四社区在线视频社区8| 母亲3免费完整高清在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品 欧美亚洲| 黄色毛片三级朝国网站| 又紧又爽又黄一区二区| 亚洲真实伦在线观看| 亚洲成av片中文字幕在线观看| а√天堂www在线а√下载| 欧美成人一区二区免费高清观看 | 一区二区日韩欧美中文字幕| 国产麻豆成人av免费视频| 久久精品亚洲精品国产色婷小说| 亚洲精品美女久久久久99蜜臀| 99久久精品国产亚洲精品| 老熟妇乱子伦视频在线观看| 90打野战视频偷拍视频| 久久热在线av| 巨乳人妻的诱惑在线观看| 国产色视频综合| 搡老熟女国产l中国老女人| 一区福利在线观看| 麻豆成人午夜福利视频| 午夜激情av网站| 99热只有精品国产| 亚洲欧美精品综合久久99| 夜夜看夜夜爽夜夜摸| 久久久国产精品麻豆| 亚洲国产中文字幕在线视频| 欧美日韩黄片免| 久久中文看片网| 久久久国产成人免费| 国产成+人综合+亚洲专区| 色综合亚洲欧美另类图片| 久久香蕉激情| 侵犯人妻中文字幕一二三四区| 老汉色av国产亚洲站长工具| 国产久久久一区二区三区| 国产亚洲av高清不卡| 欧美黑人欧美精品刺激| 色av中文字幕| 韩国精品一区二区三区| 久久性视频一级片| 婷婷亚洲欧美| 成人亚洲精品一区在线观看| 国产高清videossex| 天天躁夜夜躁狠狠躁躁| 18禁裸乳无遮挡免费网站照片 | 欧美黑人欧美精品刺激| 精品久久久久久成人av| 日韩精品青青久久久久久| 人成视频在线观看免费观看| 黑人操中国人逼视频| 亚洲精品在线观看二区| 久久中文字幕一级| 99热这里只有精品一区 | 亚洲七黄色美女视频| 最新美女视频免费是黄的| 久久欧美精品欧美久久欧美| 99久久无色码亚洲精品果冻| 美国免费a级毛片| 黄色女人牲交| 可以免费在线观看a视频的电影网站| 一二三四在线观看免费中文在| 一二三四社区在线视频社区8| 久久草成人影院| 国产一级毛片七仙女欲春2 | 一进一出抽搐gif免费好疼| 欧美zozozo另类| 国产成人啪精品午夜网站| 亚洲中文字幕一区二区三区有码在线看 | 黑人巨大精品欧美一区二区mp4| a级毛片a级免费在线| 亚洲第一电影网av| 国产成人系列免费观看| 亚洲成人久久爱视频| 黄色成人免费大全| а√天堂www在线а√下载| 欧美一级a爱片免费观看看 | 高清毛片免费观看视频网站| 国产精品国产高清国产av| 国产精品九九99| 日本在线视频免费播放| 欧美在线黄色| 成人国语在线视频| 人人妻,人人澡人人爽秒播| 午夜免费激情av| 老司机深夜福利视频在线观看| 88av欧美| 中文字幕最新亚洲高清| 午夜a级毛片| 精品一区二区三区四区五区乱码| 成人特级黄色片久久久久久久| 91av网站免费观看| 亚洲av成人一区二区三| 欧美激情极品国产一区二区三区| 午夜久久久久精精品| 天天躁狠狠躁夜夜躁狠狠躁| 国产极品粉嫩免费观看在线| 精品一区二区三区视频在线观看免费| 精品久久久久久久末码| 757午夜福利合集在线观看| 久久久水蜜桃国产精品网| 亚洲国产中文字幕在线视频| 搡老岳熟女国产| 午夜精品久久久久久毛片777| 曰老女人黄片| 男女之事视频高清在线观看| 国产亚洲精品av在线| 18禁国产床啪视频网站| 国产精品久久电影中文字幕| 亚洲精品一区av在线观看| 国产色视频综合| 亚洲色图av天堂| 99国产综合亚洲精品| 一级毛片女人18水好多| 欧美黄色片欧美黄色片| 亚洲天堂国产精品一区在线| 两个人免费观看高清视频| 日本撒尿小便嘘嘘汇集6| 亚洲精品在线美女| 亚洲国产欧美一区二区综合| 757午夜福利合集在线观看| 国产精品爽爽va在线观看网站 | 巨乳人妻的诱惑在线观看| 国产成人一区二区三区免费视频网站| 久久久久久九九精品二区国产 | 精品久久久久久久久久免费视频| 欧美在线黄色| 欧美成人免费av一区二区三区| 午夜福利在线在线| 国产极品粉嫩免费观看在线| 欧美一区二区精品小视频在线| 99riav亚洲国产免费| 日韩精品中文字幕看吧| 欧美成人性av电影在线观看| 搞女人的毛片| 免费看美女性在线毛片视频| 欧美乱码精品一区二区三区| 90打野战视频偷拍视频| 一级毛片女人18水好多| 亚洲av成人av| 日韩高清综合在线| 白带黄色成豆腐渣| 亚洲一码二码三码区别大吗| 欧美性长视频在线观看| 女生性感内裤真人,穿戴方法视频| 国内久久婷婷六月综合欲色啪| 亚洲成人精品中文字幕电影| xxxwww97欧美| 国产av在哪里看| 亚洲人成电影免费在线| 老司机靠b影院| 午夜免费激情av| 亚洲欧美精品综合一区二区三区| 狠狠狠狠99中文字幕| 久久精品国产综合久久久| 亚洲精品美女久久av网站| www.熟女人妻精品国产| 一进一出好大好爽视频| 亚洲avbb在线观看| av在线播放免费不卡| 国产区一区二久久| 久久人人精品亚洲av| 免费在线观看黄色视频的| 欧美三级亚洲精品| 成年人黄色毛片网站| 国产视频一区二区在线看| 久久亚洲真实| 在线观看66精品国产| 成人午夜高清在线视频 | 日韩大码丰满熟妇| 成人特级黄色片久久久久久久| 久久天躁狠狠躁夜夜2o2o| 亚洲午夜精品一区,二区,三区| 免费女性裸体啪啪无遮挡网站| 亚洲国产日韩欧美精品在线观看 | 国内久久婷婷六月综合欲色啪| 久久国产乱子伦精品免费另类| bbb黄色大片| 每晚都被弄得嗷嗷叫到高潮| 亚洲人成电影免费在线| 亚洲 国产 在线| 精品熟女少妇八av免费久了| 亚洲全国av大片| 国产熟女午夜一区二区三区| 日本免费一区二区三区高清不卡| 精品熟女少妇八av免费久了| 自线自在国产av| 亚洲成人久久爱视频| 国产高清视频在线播放一区| 可以在线观看的亚洲视频| 亚洲精品中文字幕在线视频| 久9热在线精品视频| 国产午夜福利久久久久久| a在线观看视频网站| 久久国产精品影院| 嫁个100分男人电影在线观看| 少妇 在线观看| 欧美黑人巨大hd| 精品福利观看| a在线观看视频网站| 国产视频一区二区在线看| 视频在线观看一区二区三区| 日本黄色视频三级网站网址| 欧美色视频一区免费| 亚洲一区二区三区色噜噜| 成人手机av| 亚洲国产高清在线一区二区三 | 欧美日韩精品网址| 国产三级黄色录像| 最新美女视频免费是黄的| 国产私拍福利视频在线观看| 国产在线精品亚洲第一网站| 变态另类丝袜制服| 亚洲中文日韩欧美视频| 久久中文字幕一级| 757午夜福利合集在线观看| 最近最新中文字幕大全免费视频| 777久久人妻少妇嫩草av网站| 久久亚洲真实| 日韩中文字幕欧美一区二区| 1024手机看黄色片| 每晚都被弄得嗷嗷叫到高潮| 中文字幕精品亚洲无线码一区 | 成年人黄色毛片网站| 久久久久久免费高清国产稀缺| 国产亚洲av嫩草精品影院| 亚洲三区欧美一区| 男人操女人黄网站| 大型黄色视频在线免费观看| 97超级碰碰碰精品色视频在线观看| 99国产极品粉嫩在线观看| 免费在线观看成人毛片| 久久久国产欧美日韩av| 亚洲欧美精品综合久久99| 国产不卡一卡二| 午夜福利成人在线免费观看| 国产欧美日韩一区二区三| 亚洲天堂国产精品一区在线| 老汉色av国产亚洲站长工具| 色播在线永久视频| 日本a在线网址| 色婷婷久久久亚洲欧美| av超薄肉色丝袜交足视频| 在线观看66精品国产| 一本久久中文字幕| 一进一出抽搐gif免费好疼| 国产国语露脸激情在线看| 婷婷精品国产亚洲av在线| aaaaa片日本免费| 国产男靠女视频免费网站| 性欧美人与动物交配| 欧美在线黄色| 中文字幕最新亚洲高清| 啦啦啦 在线观看视频| 亚洲中文日韩欧美视频| 免费女性裸体啪啪无遮挡网站| 又黄又爽又免费观看的视频| 两个人视频免费观看高清| 韩国av一区二区三区四区| 久久久久久国产a免费观看| 亚洲午夜理论影院| 亚洲av第一区精品v没综合| 久久久久国产精品人妻aⅴ院| av超薄肉色丝袜交足视频| 中文亚洲av片在线观看爽| 精品欧美国产一区二区三| 大型黄色视频在线免费观看| www国产在线视频色| av在线播放免费不卡| 国产片内射在线| 深夜精品福利| 亚洲中文日韩欧美视频| а√天堂www在线а√下载| 国语自产精品视频在线第100页| www.www免费av| 国产成人欧美| 亚洲成a人片在线一区二区| 欧美在线黄色| 91国产中文字幕| 亚洲精品国产区一区二| 99热6这里只有精品| 久久人人精品亚洲av| 校园春色视频在线观看| 日本精品一区二区三区蜜桃| 精品卡一卡二卡四卡免费| 日本a在线网址| 在线观看66精品国产| tocl精华| 午夜影院日韩av| 亚洲国产中文字幕在线视频| 波多野结衣高清作品| 波多野结衣巨乳人妻| 制服丝袜大香蕉在线| 国产精品98久久久久久宅男小说| 1024手机看黄色片| 亚洲美女黄片视频| 在线播放国产精品三级| 久久香蕉国产精品| 亚洲va日本ⅴa欧美va伊人久久| 色精品久久人妻99蜜桃| 一区二区三区高清视频在线| 国产单亲对白刺激| 国产私拍福利视频在线观看| 国产激情欧美一区二区| 午夜亚洲福利在线播放| 国产亚洲精品久久久久久毛片| 国产真实乱freesex| 午夜影院日韩av| 久久伊人香网站| 国产精品一区二区三区四区久久 | 亚洲国产欧洲综合997久久, | 久久精品91无色码中文字幕| 午夜视频精品福利| 成人国语在线视频| 国产私拍福利视频在线观看| 看黄色毛片网站| 亚洲一码二码三码区别大吗| 国产99白浆流出| 亚洲人成伊人成综合网2020| 国内揄拍国产精品人妻在线 | 黄色女人牲交| 久久 成人 亚洲| 亚洲av第一区精品v没综合| 黑人巨大精品欧美一区二区mp4| 在线观看免费午夜福利视频| 成人三级黄色视频| 欧美久久黑人一区二区| 99国产综合亚洲精品| 国产精品98久久久久久宅男小说| 在线观看免费午夜福利视频| 99re在线观看精品视频| 精品国产乱子伦一区二区三区| 叶爱在线成人免费视频播放| 日本成人三级电影网站| 女性生殖器流出的白浆| 老熟妇乱子伦视频在线观看| 欧美人与性动交α欧美精品济南到| 2021天堂中文幕一二区在线观 | 88av欧美| 成人永久免费在线观看视频| 久久久久九九精品影院| 亚洲真实伦在线观看| 色尼玛亚洲综合影院| 国产一区二区激情短视频| 欧美国产精品va在线观看不卡| 91国产中文字幕| 亚洲真实伦在线观看| 久久伊人香网站| 女生性感内裤真人,穿戴方法视频| 日本熟妇午夜| 欧美成狂野欧美在线观看| 国产亚洲精品第一综合不卡| 一本精品99久久精品77| 国产av不卡久久| 日本熟妇午夜| 国产欧美日韩一区二区精品| 97人妻精品一区二区三区麻豆 | 两个人免费观看高清视频| 丁香六月欧美| 亚洲国产精品成人综合色| 少妇粗大呻吟视频| 久久久久国产一级毛片高清牌| 99热只有精品国产| 无限看片的www在线观看| 日韩高清综合在线| 男女午夜视频在线观看| 观看免费一级毛片| 丝袜美腿诱惑在线| 在线观看www视频免费| 在线av久久热| 国产精品二区激情视频| 欧美国产精品va在线观看不卡| 国产av又大| 亚洲国产毛片av蜜桃av| 人妻丰满熟妇av一区二区三区| 精品免费久久久久久久清纯| 麻豆久久精品国产亚洲av| 欧美一级毛片孕妇| 老汉色av国产亚洲站长工具| 成人国产综合亚洲| 日本一区二区免费在线视频| 黄色女人牲交| 在线观看舔阴道视频| 久久久久九九精品影院| 久久久国产欧美日韩av| 黄色a级毛片大全视频| 日本a在线网址| 12—13女人毛片做爰片一| 精品欧美一区二区三区在线| 久久精品影院6| 美女国产高潮福利片在线看| bbb黄色大片| 少妇被粗大的猛进出69影院| 欧美日本视频| 午夜福利视频1000在线观看| 精品久久久久久,| 亚洲欧洲精品一区二区精品久久久| 欧美黑人欧美精品刺激| 免费一级毛片在线播放高清视频| 国产精品自产拍在线观看55亚洲| 免费在线观看视频国产中文字幕亚洲| 精品久久久久久久末码| 亚洲一区二区三区色噜噜| 国产麻豆成人av免费视频| 精品久久久久久,| 精品久久久久久久久久久久久 | 老司机午夜十八禁免费视频| 亚洲美女黄片视频| 欧美一级a爱片免费观看看 | 国产精品 欧美亚洲| 精品久久久久久久毛片微露脸| 久久香蕉精品热| 狂野欧美激情性xxxx| 变态另类丝袜制服| 变态另类成人亚洲欧美熟女| 国产精品美女特级片免费视频播放器 | netflix在线观看网站| 日韩高清综合在线| 亚洲国产精品sss在线观看| 脱女人内裤的视频| 最近最新免费中文字幕在线| 国产亚洲精品第一综合不卡| 高潮久久久久久久久久久不卡| 亚洲aⅴ乱码一区二区在线播放 | 听说在线观看完整版免费高清| 久久精品91蜜桃| 亚洲av电影不卡..在线观看| videosex国产| 亚洲 欧美 日韩 在线 免费| 欧美一级毛片孕妇| 久久久国产精品麻豆| 夜夜看夜夜爽夜夜摸| 精品日产1卡2卡| 淫秽高清视频在线观看| 国产精品99久久99久久久不卡| 后天国语完整版免费观看| 国产亚洲精品第一综合不卡| 国产99久久九九免费精品| 久久精品91无色码中文字幕| 亚洲成人精品中文字幕电影| 国产激情欧美一区二区| 这个男人来自地球电影免费观看| 国产精品久久久av美女十八| 精品久久久久久久人妻蜜臀av| 嫁个100分男人电影在线观看| 国产又爽黄色视频| 人人妻人人看人人澡| 又大又爽又粗| 国产真人三级小视频在线观看| 看黄色毛片网站| 欧美成狂野欧美在线观看| 日本 av在线| 好看av亚洲va欧美ⅴa在| 一级作爱视频免费观看| 亚洲激情在线av| 亚洲最大成人中文|