• <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ī)則圖形面積
    水平磁偶極子電磁場特征研究
    综合色av麻豆| 亚洲av成人精品一二三区| 色吧在线观看| 综合色丁香网| 国产老妇伦熟女老妇高清| 久久精品久久久久久久性| 国产精品乱码一区二三区的特点| 国产一区二区在线av高清观看| 久久99热这里只有精品18| 色吧在线观看| 国产精品一区二区三区四区久久| 搡女人真爽免费视频火全软件| 大又大粗又爽又黄少妇毛片口| 三级国产精品欧美在线观看| 熟妇人妻久久中文字幕3abv| 伊人久久精品亚洲午夜| 久久久a久久爽久久v久久| 亚洲精品一区蜜桃| av线在线观看网站| 免费看日本二区| 成人亚洲精品av一区二区| 亚洲精品456在线播放app| 一级黄片播放器| 亚洲成人久久爱视频| 2021少妇久久久久久久久久久| 69人妻影院| av在线播放精品| av在线播放精品| 国语对白做爰xxxⅹ性视频网站| 99在线视频只有这里精品首页| 欧美极品一区二区三区四区| 中文字幕av在线有码专区| 国产精品久久久久久久电影| 色综合亚洲欧美另类图片| 久久综合国产亚洲精品| 国产乱来视频区| 国产日韩欧美在线精品| 国产av码专区亚洲av| 黄片wwwwww| 免费大片18禁| 久久99热6这里只有精品| 国产激情偷乱视频一区二区| 日韩欧美 国产精品| 2022亚洲国产成人精品| 黑人高潮一二区| 亚洲精品国产成人久久av| 69人妻影院| 99久国产av精品| 日日摸夜夜添夜夜添av毛片| 亚洲精品久久久久久婷婷小说 | 色播亚洲综合网| 午夜激情福利司机影院| 国产激情偷乱视频一区二区| 天天躁夜夜躁狠狠久久av| 亚州av有码| 日本一二三区视频观看| videos熟女内射| 女人久久www免费人成看片 | 成人综合一区亚洲| 天堂中文最新版在线下载 | 亚洲天堂国产精品一区在线| 人体艺术视频欧美日本| 丰满人妻一区二区三区视频av| 欧美成人一区二区免费高清观看| 日韩欧美国产在线观看| 女人久久www免费人成看片 | 波多野结衣巨乳人妻| 22中文网久久字幕| 淫秽高清视频在线观看| 麻豆乱淫一区二区| 赤兔流量卡办理| 天天躁夜夜躁狠狠久久av| 男女国产视频网站| 人妻少妇偷人精品九色| 18+在线观看网站| 中文字幕av成人在线电影| 床上黄色一级片| 国产老妇伦熟女老妇高清| 黄色欧美视频在线观看| 三级国产精品欧美在线观看| 欧美一区二区国产精品久久精品| 久久精品夜色国产| 高清视频免费观看一区二区 | av在线天堂中文字幕| 国产视频内射| 久久久久久久久久成人| 3wmmmm亚洲av在线观看| 欧美成人一区二区免费高清观看| 国产免费又黄又爽又色| 人妻夜夜爽99麻豆av| 中文精品一卡2卡3卡4更新| 午夜精品在线福利| 精品久久久久久久人妻蜜臀av| 久久热精品热| 偷拍熟女少妇极品色| 桃色一区二区三区在线观看| 亚洲最大成人中文| 国产不卡一卡二| 国产精品爽爽va在线观看网站| 午夜激情欧美在线| 亚洲在久久综合| 精华霜和精华液先用哪个| 日本爱情动作片www.在线观看| 久久精品国产99精品国产亚洲性色| 欧美人与善性xxx| 精品久久久久久久久亚洲| 久久6这里有精品| 欧美高清性xxxxhd video| 国产白丝娇喘喷水9色精品| 久久亚洲精品不卡| 大香蕉97超碰在线| 麻豆一二三区av精品| 天天躁夜夜躁狠狠久久av| 噜噜噜噜噜久久久久久91| 99热全是精品| 日本熟妇午夜| 亚洲av男天堂| 国产人妻一区二区三区在| 国产激情偷乱视频一区二区| 国产一区二区三区av在线| 久久精品久久精品一区二区三区| 别揉我奶头 嗯啊视频| 国产熟女欧美一区二区| 九九热线精品视视频播放| 热99re8久久精品国产| 尾随美女入室| 少妇丰满av| 天天躁日日操中文字幕| 99热精品在线国产| 亚洲欧美日韩无卡精品| 美女内射精品一级片tv| 老师上课跳d突然被开到最大视频| av.在线天堂| 亚洲国产欧美在线一区| 久久99蜜桃精品久久| 岛国在线免费视频观看| 老女人水多毛片| 大话2 男鬼变身卡| 99九九线精品视频在线观看视频| 日本三级黄在线观看| 日本黄色视频三级网站网址| 色吧在线观看| 久久精品国产亚洲网站| 亚洲最大成人中文| av又黄又爽大尺度在线免费看 | 亚洲成人av在线免费| 最近最新中文字幕免费大全7| 99久国产av精品国产电影| 亚洲丝袜综合中文字幕| 亚洲成人精品中文字幕电影| 国产高清有码在线观看视频| av在线老鸭窝| 国产亚洲5aaaaa淫片| 久久99热这里只频精品6学生 | 久久99精品国语久久久| 18禁在线无遮挡免费观看视频| 国产亚洲av片在线观看秒播厂 | 汤姆久久久久久久影院中文字幕 | 国产欧美另类精品又又久久亚洲欧美| 国内精品美女久久久久久| 高清午夜精品一区二区三区| 亚洲电影在线观看av| 日韩欧美精品v在线| av免费观看日本| 蜜桃久久精品国产亚洲av| 亚洲成人久久爱视频| 亚洲一区高清亚洲精品| 久久久久性生活片| 精品人妻熟女av久视频| 嘟嘟电影网在线观看| 欧美极品一区二区三区四区| 中文字幕久久专区| av福利片在线观看| 亚洲欧洲国产日韩| 嫩草影院精品99| 99久国产av精品国产电影| 欧美成人免费av一区二区三区| 国产成人精品一,二区| 搞女人的毛片| 国产女主播在线喷水免费视频网站 | 综合色av麻豆| 简卡轻食公司| 亚洲国产精品国产精品| 欧美激情国产日韩精品一区| 91午夜精品亚洲一区二区三区| 丰满人妻一区二区三区视频av| 91在线精品国自产拍蜜月| 韩国av在线不卡| 欧美潮喷喷水| 精品国产一区二区三区久久久樱花 | 国产精品人妻久久久久久| 嫩草影院新地址| 久久精品夜色国产| 看片在线看免费视频| 国产免费一级a男人的天堂| 亚洲一级一片aⅴ在线观看| av卡一久久| 超碰97精品在线观看| 波多野结衣巨乳人妻| 边亲边吃奶的免费视频| 亚洲一级一片aⅴ在线观看| 日韩成人av中文字幕在线观看| 日本爱情动作片www.在线观看| 69av精品久久久久久| 亚洲成人久久爱视频| 国产淫语在线视频| 国产一级毛片七仙女欲春2| 欧美日韩精品成人综合77777| 偷拍熟女少妇极品色| 建设人人有责人人尽责人人享有的 | av视频在线观看入口| a级毛色黄片| 亚洲精华国产精华液的使用体验| 午夜日本视频在线| 欧美精品国产亚洲| 国产视频内射| 久久久久久久久久久免费av| 26uuu在线亚洲综合色| 岛国在线免费视频观看| 午夜视频国产福利| 国产精品av视频在线免费观看| 熟女人妻精品中文字幕| 不卡视频在线观看欧美| 国产成年人精品一区二区| 最近的中文字幕免费完整| 欧美+日韩+精品| 日本wwww免费看| 国产精品日韩av在线免费观看| 99热这里只有精品一区| 身体一侧抽搐| av国产久精品久网站免费入址| 成人综合一区亚洲| 国产三级中文精品| 春色校园在线视频观看| 欧美高清成人免费视频www| 自拍偷自拍亚洲精品老妇| 男女视频在线观看网站免费| 成年免费大片在线观看| 精品少妇黑人巨大在线播放 | 97超视频在线观看视频| 韩国高清视频一区二区三区| 日本黄色视频三级网站网址| 男女国产视频网站| 日韩视频在线欧美| 欧美日韩精品成人综合77777| 性色avwww在线观看| 欧美+日韩+精品| 七月丁香在线播放| 91精品国产九色| 夫妻性生交免费视频一级片| 舔av片在线| 在线播放国产精品三级| 三级男女做爰猛烈吃奶摸视频| 免费看a级黄色片| 色综合站精品国产| 午夜爱爱视频在线播放| 99热这里只有精品一区| 大又大粗又爽又黄少妇毛片口| 亚洲精品自拍成人| 激情 狠狠 欧美| 亚洲人与动物交配视频| 中文字幕久久专区| 日韩av不卡免费在线播放| 69av精品久久久久久| 午夜福利成人在线免费观看| 黄色配什么色好看| 亚洲成人av在线免费| 久久国产乱子免费精品| 麻豆国产97在线/欧美| 最近手机中文字幕大全| 成人性生交大片免费视频hd| 国产成人aa在线观看| 成人亚洲精品av一区二区| 97超碰精品成人国产| 最近中文字幕2019免费版| av.在线天堂| 国产成人免费观看mmmm| 综合色丁香网| 国产男人的电影天堂91| 一边摸一边抽搐一进一小说| 久久久久久久国产电影| 伦理电影大哥的女人| 精品人妻偷拍中文字幕| 国产高清视频在线观看网站| 亚洲av一区综合| 欧美3d第一页| 成人高潮视频无遮挡免费网站| 日产精品乱码卡一卡2卡三| 一边亲一边摸免费视频| 乱系列少妇在线播放| 99热精品在线国产| 欧美又色又爽又黄视频| 亚洲国产精品成人综合色| 久久精品国产亚洲av天美| 亚洲欧美日韩东京热| 国内少妇人妻偷人精品xxx网站| 97在线视频观看| 国产精品一区www在线观看| av在线亚洲专区| 成人毛片60女人毛片免费| 一边亲一边摸免费视频| 一级黄片播放器| 小说图片视频综合网站| 免费观看a级毛片全部| 99九九线精品视频在线观看视频| 国产亚洲一区二区精品| 日韩国内少妇激情av| 久久久亚洲精品成人影院| 久久这里只有精品中国| 毛片一级片免费看久久久久| 国产精品美女特级片免费视频播放器| 国产激情偷乱视频一区二区| 欧美又色又爽又黄视频| 日韩成人av中文字幕在线观看| 一本一本综合久久| 中文字幕久久专区| av国产免费在线观看| 成人漫画全彩无遮挡| 美女高潮的动态| 熟妇人妻久久中文字幕3abv| 搡女人真爽免费视频火全软件| 国产亚洲av嫩草精品影院| 国产精品熟女久久久久浪| 欧美成人精品欧美一级黄| 尾随美女入室| 亚洲精品日韩在线中文字幕| 亚洲av福利一区| 99在线视频只有这里精品首页| 直男gayav资源| ponron亚洲| 99九九线精品视频在线观看视频| 欧美成人精品欧美一级黄| 欧美激情久久久久久爽电影| 日韩成人伦理影院| 久久99蜜桃精品久久| 丰满人妻一区二区三区视频av| 能在线免费看毛片的网站| 直男gayav资源| h日本视频在线播放| 男女下面进入的视频免费午夜| 观看美女的网站| 久久精品国产亚洲网站| 久久精品熟女亚洲av麻豆精品 | 日本欧美国产在线视频| 欧美一级a爱片免费观看看| 国内精品美女久久久久久| 国产真实伦视频高清在线观看| 自拍偷自拍亚洲精品老妇| 丰满乱子伦码专区| 丝袜美腿在线中文| 真实男女啪啪啪动态图| 又粗又爽又猛毛片免费看| 亚洲欧美日韩高清专用| 国产白丝娇喘喷水9色精品| 99在线人妻在线中文字幕| 国产乱人偷精品视频| 久久99蜜桃精品久久| 免费搜索国产男女视频| 国产伦一二天堂av在线观看| 免费观看性生交大片5| 午夜视频国产福利| 亚洲熟妇中文字幕五十中出| 久久这里有精品视频免费| 国产精品一及| 别揉我奶头 嗯啊视频| 国产精品av视频在线免费观看| 能在线免费看毛片的网站| 久久久久久久久久成人| 欧美精品国产亚洲| 青春草视频在线免费观看| 精品国产三级普通话版| 秋霞在线观看毛片| 国产av在哪里看| 国产在视频线精品| 欧美区成人在线视频| 亚洲va在线va天堂va国产| 国产伦精品一区二区三区四那| 国产精品国产三级专区第一集| 97热精品久久久久久| 国产伦理片在线播放av一区| 国产精品人妻久久久影院| 91精品国产九色| 搞女人的毛片| 97超视频在线观看视频| 国产精品永久免费网站| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲精品一区蜜桃| 中文欧美无线码| 免费看日本二区| 日韩一本色道免费dvd| 18禁在线播放成人免费| www.av在线官网国产| 有码 亚洲区| av线在线观看网站| 91aial.com中文字幕在线观看| 蜜桃亚洲精品一区二区三区| 欧美日本亚洲视频在线播放| 在线免费观看不下载黄p国产| 寂寞人妻少妇视频99o| 性插视频无遮挡在线免费观看| 精品少妇黑人巨大在线播放 | 黄色欧美视频在线观看| 国产精品一及| 国产亚洲av片在线观看秒播厂 | 亚洲真实伦在线观看| 国产在视频线精品| 国产精品1区2区在线观看.| 狠狠狠狠99中文字幕| 亚洲国产精品国产精品| 九九久久精品国产亚洲av麻豆| 久热久热在线精品观看| 又粗又爽又猛毛片免费看| 在现免费观看毛片| 亚洲中文字幕日韩| 欧美xxxx黑人xx丫x性爽| 国产伦一二天堂av在线观看| 男女视频在线观看网站免费| 少妇高潮的动态图| 亚洲欧美日韩卡通动漫| 亚洲成色77777| 欧美一区二区亚洲| 老女人水多毛片| 亚洲怡红院男人天堂| 午夜福利视频1000在线观看| 在线观看av片永久免费下载| 欧美xxxx性猛交bbbb| 蜜臀久久99精品久久宅男| 国产中年淑女户外野战色| 亚洲av中文av极速乱| 简卡轻食公司| 日本五十路高清| 非洲黑人性xxxx精品又粗又长| 91aial.com中文字幕在线观看| 美女高潮的动态| 国产一区二区在线观看日韩| 亚洲美女搞黄在线观看| 国产私拍福利视频在线观看| 久久精品国产亚洲网站| 国产精品福利在线免费观看| 国产精品一二三区在线看| 中文字幕av成人在线电影| 禁无遮挡网站| 1000部很黄的大片| 国产69精品久久久久777片| 男人舔女人下体高潮全视频| 真实男女啪啪啪动态图| 床上黄色一级片| 色5月婷婷丁香| 青青草视频在线视频观看| 精品少妇黑人巨大在线播放 | 国产伦理片在线播放av一区| 日本免费a在线| kizo精华| 亚洲国产精品成人综合色| 三级国产精品片| 看免费成人av毛片| 免费看a级黄色片| 免费看av在线观看网站| 日产精品乱码卡一卡2卡三| 老司机福利观看| 婷婷色麻豆天堂久久 | 麻豆久久精品国产亚洲av| 一级毛片我不卡| av在线蜜桃| 久久精品国产亚洲网站| 国产午夜精品久久久久久一区二区三区| 成人性生交大片免费视频hd| 22中文网久久字幕| 国产精品美女特级片免费视频播放器| 自拍偷自拍亚洲精品老妇| 男人和女人高潮做爰伦理| 国产av在哪里看| 人妻系列 视频| 你懂的网址亚洲精品在线观看 | 亚洲美女搞黄在线观看| av又黄又爽大尺度在线免费看 | 国产成人午夜福利电影在线观看| 日韩精品有码人妻一区| 久久久久久久亚洲中文字幕| 精品午夜福利在线看| 国产亚洲精品av在线| 久久久精品欧美日韩精品| 中文资源天堂在线| 亚洲精品影视一区二区三区av| 国产 一区 欧美 日韩| 欧美精品一区二区大全| 一区二区三区四区激情视频| 久久久久久久久中文| 久久久久久久久久黄片| 国产精品一二三区在线看| 亚洲精品aⅴ在线观看| 亚洲国产最新在线播放| 国产探花极品一区二区| 99久久中文字幕三级久久日本| 一卡2卡三卡四卡精品乱码亚洲| av播播在线观看一区| 成人漫画全彩无遮挡| 免费人成在线观看视频色| 国产精品日韩av在线免费观看| 中文在线观看免费www的网站| 日本黄大片高清| 狠狠狠狠99中文字幕| 欧美97在线视频| www日本黄色视频网| 一级毛片我不卡| 伦理电影大哥的女人| 狂野欧美激情性xxxx在线观看| 亚洲av不卡在线观看| 91精品一卡2卡3卡4卡| 婷婷色综合大香蕉| 精品人妻熟女av久视频| 看十八女毛片水多多多| 18禁在线播放成人免费| 成人美女网站在线观看视频| 亚洲成人av在线免费| 69av精品久久久久久| 免费av不卡在线播放| or卡值多少钱| av线在线观看网站| av在线观看视频网站免费| 毛片女人毛片| 少妇高潮的动态图| 精品午夜福利在线看| 99久国产av精品| 日韩欧美精品v在线| 一边摸一边抽搐一进一小说| 欧美不卡视频在线免费观看| 毛片一级片免费看久久久久| 国产精品99久久久久久久久| 亚洲av成人av| 乱系列少妇在线播放| 99久久九九国产精品国产免费| 国产精品人妻久久久影院| 天堂影院成人在线观看| 亚洲欧美成人综合另类久久久 | 亚洲国产精品合色在线| 午夜精品在线福利| 亚洲精品aⅴ在线观看| 国产亚洲91精品色在线| 赤兔流量卡办理| 高清日韩中文字幕在线| 国产91av在线免费观看| 深夜a级毛片| 在线播放国产精品三级| 国产高清三级在线| 3wmmmm亚洲av在线观看| 久久精品国产鲁丝片午夜精品| 国内揄拍国产精品人妻在线| 99久久成人亚洲精品观看| 亚洲国产欧洲综合997久久,| 国产午夜精品一二区理论片| 天美传媒精品一区二区| 国产黄片美女视频| 日本欧美国产在线视频| 又粗又爽又猛毛片免费看| 热99re8久久精品国产| 在线观看av片永久免费下载| 日韩高清综合在线| 国产久久久一区二区三区| 久久鲁丝午夜福利片| 日日撸夜夜添| 啦啦啦啦在线视频资源| 一区二区三区高清视频在线| 1024手机看黄色片| 老司机影院毛片| 亚洲精品国产成人久久av| 国产高清三级在线| 国产女主播在线喷水免费视频网站 | 午夜日本视频在线| 亚洲国产欧美在线一区| 国产精品久久久久久精品电影小说 | 精品久久久久久久人妻蜜臀av| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品永久免费网站| 久久国产乱子免费精品| 国产老妇伦熟女老妇高清| 国产午夜精品一二区理论片| 国产私拍福利视频在线观看| 亚洲国产最新在线播放| 水蜜桃什么品种好| 中文天堂在线官网| 美女内射精品一级片tv| 久久亚洲国产成人精品v| 国产成年人精品一区二区| 久久久精品大字幕| 国产精品国产三级国产av玫瑰| 亚洲人成网站高清观看| 国产成人精品婷婷| 真实男女啪啪啪动态图| 国内精品美女久久久久久| 欧美色视频一区免费| 精品人妻偷拍中文字幕| 一级毛片我不卡| 中文亚洲av片在线观看爽| 在现免费观看毛片| 蜜臀久久99精品久久宅男| 精品99又大又爽又粗少妇毛片| 日韩欧美三级三区| 亚洲国产成人一精品久久久| 精品久久久噜噜| 噜噜噜噜噜久久久久久91| 91aial.com中文字幕在线观看| 婷婷色综合大香蕉| 亚洲aⅴ乱码一区二区在线播放| 秋霞伦理黄片| 中文精品一卡2卡3卡4更新| 亚洲熟妇中文字幕五十中出| 国产大屁股一区二区在线视频| 欧美三级亚洲精品| 久久午夜福利片| 日韩国内少妇激情av| 日本av手机在线免费观看| 免费av观看视频| or卡值多少钱| 人人妻人人看人人澡|