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

    DGTD用于RCS計算的初步研究

    2015-03-08 02:35:30謙魏李林茜葛德彪
    雷達學報 2015年3期
    關鍵詞:總場入射波四面體

    楊 謙魏 兵*李林茜葛德彪

    ①(西安電子科技大學物理與光電工程學院 西安 710071)

    ②(西安電子科技大學信息感知技術協(xié)同創(chuàng)新中心 西安 710071)

    DGTD用于RCS計算的初步研究

    楊 謙①②魏 兵*①②李林茜①②葛德彪①②

    ①(西安電子科技大學物理與光電工程學院 西安 710071)

    ②(西安電子科技大學信息感知技術協(xié)同創(chuàng)新中心 西安 710071)

    時域離散伽遼金法(Discontinuous Galerkin Time Domain, DGTD)同時具有時域有限元算法(Finite Element Time Domain, FETD)非結(jié)構(gòu)網(wǎng)格剖分和時域有限差分算法(Finite Difference Time Domain, FDTD)顯式迭代的優(yōu)點,是一種非常有前途的電磁計算方法,該文首先描述了基于矢量基函數(shù)的時域離散伽遼金法的基本原理。然后,給出了DGTD處理散射問題時平面波入射加入的具體實現(xiàn)方法。最后,給出了金屬球、介質(zhì)球和金屬彈頭寬帶散射的算例,算例結(jié)果的比較表明了該文算法的正確性和有效性。該文的研究,為復雜目標雷達散射截面RCS的準確預估打下了堅實的基礎。

    時域離散伽遼金方法;時域有限差分;有限元;雷達散射截面

    1 引言

    在目標特性的精確預估中,矩量法(Method of Moment, MoM)、有限元法(Finite Element Method, FEM)和時域有限差分方法(Finite Difference Time Domain, FDTD)均是應用廣泛的計算方法。MoM主要長于金屬目標散射問題的處理,F(xiàn)EM和FDTD方法則更適宜于復雜媒質(zhì)的計算。MoM和FEM需要求解大型線性方程組,應用于電大尺寸目標時求解效率低。FDTD采用顯式迭代的方式,易于并行計算,并且計算復雜度低,但其在對目標擬合上存在臺階誤差,即使采取亞網(wǎng)格或共形網(wǎng)格也很難精確擬合。近年來,離散伽遼金法(Discontinuous Galerkin method)算法受到了人們的重視[1]。該算法在流體力學及有限體積法

    (Finite volume method)中早有應用[2],近年來被引入有限元計算當中來形成了時域離散伽遼金算法(Discontinuous Galerkin Time Domain, DGTD)[3,4]。DGTD結(jié)合了FDTD顯式迭代和FEM網(wǎng)格精確擬合目標幾何性形狀的特點,是一種有廣闊應用前景的計算方法。近十余年,國際上的同行對該方法已經(jīng)進行了多方面的研究,但在國內(nèi)將基于矢量基函數(shù)的DGTD方法用于目標散射特性的分析還未見相關報道。本文首先簡要介紹了基于棱邊矢量基函數(shù)的3維DGTD方法的基本思想。然后,給出了在該方法實現(xiàn)中的平面波引入及總場邊界條件。最后,給出目標單站散射特性的算例表明本文算法的正確性。

    2 基于棱邊基函數(shù)的3 維DGTD

    DGTD可以看作由時域有限元方法(Finite Element Time Domain, FETD)改變邊界條件的處理方式得到。FETD是從支配方程和邊界條件出發(fā)將計算區(qū)域劃分為多個單元后導出矩陣方程并求解的方法,J. M. Jin, J. F. Lee等人在FETD方面做了很多工作[5–7]。一般來說,F(xiàn)ETD方程推導有兩種途徑,即變分法與Galerkin法。變分法尋找一個適當泛函對應于支配方程和邊界條件,Galerkin法通過加權余量來尋找支配方程和邊界條件的相應弱解形式。兩種方法都會得到一個大型矩陣方程。在FETD中的若干基本概念如基函數(shù)、單元積分、局域、全域、Galerkin加權等均可在DGTD中繼承使用,而兩種有限元算法的局域邊界處理方式有差異。

    下面從Maxwell方程組出發(fā)描述DGTD的基本思想,Maxwell旋度方程為:

    式中ε為介電系數(shù),μ為磁導系數(shù);σ為電導率,σm為導磁率;E和H分別為電場強度矢量與磁場強度矢量;J與Jm分別為電磁流。

    本文中采用四面體作為DGTD算法的基本離散單元(如圖1所示),采用Galerkin加權法[6],對式(1)在四面體單元內(nèi)積分,可得

    圖1 四面體單元示意圖Fig. 1 Schematic diagram of tetrahedron

    式中v為積分權函數(shù)。對式(2)中包含及的體積分項做代數(shù)變換

    在流體力學中為求解雙曲方程,邊界條件采用Numerical Flux[2],后又被引入了時域有限體積法(FVTD)[8],又被引入DGTD作核心思想[4]。在邊界上定義E*與H*為:

    式(4)中E*,H*稱為Numerical Flux,在兩個相鄰的四面體單元面兩側(cè)不相等但依賴于兩個相接面的切向場;為外法向單位矢量,E,H為待計算單元場值,E+,H+為與待計算單元的相鄰單元場值。

    由于E*,H*僅定義在邊界上,可得

    表1 Numerical Flux系數(shù)表達式Tab. 1 Coefficient of Numerical Flux

    經(jīng)整理后得矩陣形式的偏微分方程組,此方程組針對單個四面體,實際計算中在計算域內(nèi)循環(huán)所有四面體迭代求解。

    3 DGTD中的平面波引入及總場邊界條件

    用DGTD計算散射問題時可將計算區(qū)域劃分為總場區(qū)和散射場區(qū),如圖2所示。采用等效原理,類似FDTD中的總場邊界條件,在總場邊界上設置入射波的切向分量即可將入射波引入到總場區(qū)[10]。

    圖2 總場區(qū)與散射場區(qū)的劃分Fig. 2 The total field region and the scattered field region

    為便于說明總場邊界條件的加入,圖3中以三角形代替四面體,規(guī)定散射場區(qū)單元的場均屬于散射場,總場區(qū)單元的場均屬于總場。在總場區(qū)內(nèi)部及散射場區(qū)內(nèi)部,DGTD計算方式如前文所述,需要特殊處理的是邊界處相鄰的單元計算式,具體來講是在Numerical Flux表達式中增加入射波的計算部分。

    圖3 DGTD總場散射場邊界Fig. 3 TF/SF boundary for DGTD

    (1) 當四面體處于總場邊界的總場區(qū)一側(cè)時,應在Numerical Flux上加上入射波值,此時屬于總場,屬于散射場,Ei, Hi為入射波。

    (2) 當四面體處于總場邊界的散射場區(qū)一側(cè)時,應在Numerical Flux上減去入射波值,此時屬于散射場,屬于總場,Ei, Hi為入射波。

    采取上面的方法即可將平面波引入總場區(qū)域。

    4 算例

    為說明本文算法的正確性,用DGTD方法計算半徑為1 m的金屬球的后向RCS并與解析解結(jié)果相對比。計算區(qū)域被離散成為319308個四面體(計算域剖分截面圖如圖4所示),空間離散尺度為0.15 m,采用1階Mur吸收邊界。圖5是金屬球后向單站RCS的計算結(jié)果(圓圈),作為比較,圖中還給出了Mie級數(shù)計算結(jié)果(實線),由圖可見兩者的結(jié)果吻合較好。本算例DGTD消耗內(nèi)存約280 MB,計算耗時約20 min,計算平臺為Intel Core i5 3470;采用FETD計算所需內(nèi)存消耗為2 G,計算時間在6 h以上。顯然,本文方法有更高的計算效率。

    設介質(zhì)球參數(shù)為ε=4ε0,μ=μ0,σ=0.00015 S/m,其余參數(shù)同上。圖6是本文方法(圓圈)與

    Mie級數(shù)計算結(jié)果(實線)的比較,由圖可見,兩者相吻合。

    圖4 計算域剖分截面圖Fig. 4 Sectional view of the computational domain

    圖5 金屬球單站RCSFig. 5 The monostatic RCS of PEC sphere

    圖6 介質(zhì)球單站RCSFig. 6 The monostatic RCS of medium sphere

    為說明本文算法對于計算較為復雜目標的適用性,下面計算金屬彈頭的后向RCS。彈頭模型如圖7所示(底面半徑為1 m,高為2 m)。電磁波沿z軸負方向,電場沿x方向極化入射。圖8為彈頭的后向RCS,實線為DGTD算法的結(jié)果,斜十字為MoM算法掃頻結(jié)果。由圖可見,兩種算法的結(jié)果符合較好。

    圖7 彈頭模型Fig. 7 Model of PEC bullet

    圖8 金屬彈頭單站RCSFig. 8 The monostatic RCS of PEC bullet

    5 結(jié)論

    在散射問題的DGTD的分析中吸收邊界條件、時間離散方式、近遠場外推和四面體相鄰單元的快速判斷等方面都需要精細考慮,限于篇幅,本文沒有討論這些內(nèi)容。目前,由于吸收邊界采用1階Mur吸收邊界,本文算例中DGTD的計算精度有待提高。采用UPML吸收邊界、高階基函數(shù)及精度更高的時間離散方案將是我們下一步的工作??梢灶A期,作為一種新興算法,DGTD將在目標特性的精確預估,特別是含復雜媒質(zhì)的低RCS目標特性的精確預估方面扮演重要角色。

    [1]Ji X, Lu T, Cai W,et al.. Discontinuous Galerkin Time Domain (DGTD) methods for the study of 2-D waveguidecoupled microring resonators[J].Journal of Lightwave Technology, 2005, 23(11): 3864–3874.

    [2]Shu C W. A brief survey on discontinuous Galerkin methods in computational fluid dynamics[J].Advances in Mechanics, 2013, 43(6): 541–553.

    [3]Gedney S D, Kramer T, Luo C,et al.. The Discontinuous Galerkin Finite Element Time Domain method (DGFETD)[C]. IEEE International Symposium on Electromagnetic Compatibility, 2008: 1–4.

    [4]Alvarez J. A Discontinuous Galerkin Finite Element Method for the Time-Domain Solution of Maxwell Equations[D]. [Ph.D. dissertation], University of Granada (Spain), 2013: 31–39.

    [5]Lee Jin-Fa, Lee Robert, and Cangellaris Andreas. Timedomain finite-element methods[J].IEEE Transactions on Antennas and Propagation, 1997, 45(3): 430–442.

    [6]Jin Jian-ming. The Finite Element Method in Electromagnetic[M]. New York: John Wiley & Sons, 2002: 22–23.

    [7]Riley D J, Jin Jian-ming, Lou Z,et al.. Total-and scatteredfield decomposition technique for the finite-element timedomain method[J].IEEE Transactions on Antennas and Propagation, 2006, 54(1): 35–41.

    [8]Shankar V, Mohammadian A H, and Hall W F. A timedomain, finite-volume treatment for the Maxwell equations[J].Electromagnetics, 1990, 10(1/2): 127–145.

    [9]Bossavit A. Whitney forms: a class of finite elements for three-dimensional computations in electromagnetism[J].IEE Proceedings A(Physical Science,Measurement and Instrumentation,Management and Education,Reviews), 1988, 135(8): 493–500.

    [10]葛德彪, 魏兵. 電磁波時域計算方法(下冊)[M]. 西安電子科技大學出版社, 2014: 188–191. Ge D B and Wei B. Time Domain Computational Method for Electromagnetic Wave (Volume II)[M]. Xi’an: Xidian University Press, 2014: 188–191.

    楊 謙(1989–),男,博士研究生,主要研究方向為計算電磁學。

    E-mail: zijiangy@126.com

    魏 兵(1970–),男,教授,博士,主要研究方向為電磁理論、復雜系統(tǒng)中的場與波和計算電磁學等。

    E-mail: bwei@xidian.edu.cn

    李林茜(1985–),男,博士研究生,主要研究方向為計算電磁學。

    E-mail: 395106835@qq.com

    Preliminary Research on RCS Using DGTD

    Yang Qian①②Wei Bing①②Li Lin-qian①②Ge De-biao①②

    ①(School of Physics and Optoelectronic Engineering,XidianUniversity,Xi’an710071,China)

    ②(Collaborative Innovation Center of Information Sensing and Understanding,XidianUniversity,Xi’an710071,China)

    Discontinuous Galerkin Time Domain (DGTD) method appears to be very promising which combines the advantages of unstructured mesh in Finite Element Time Domain (FETD) and explicit scheme in Finite Difference Time Domain (FDTD). This paper first describes principle of DGTD base on vector basis function. Secondly, Specific method for incident plane wave is given for scattering problem. At last, the monostatic Radar Cross Section (RCS) of PEC sphere, medium sphere and the PEC bullet are computed by DGTD method. The numerical results illustrate the feasibility and correctness of the presented scheme. The study of this paper is a foundation for analyzing the RCS of complex target.

    Discontinuous Galerkin Time Domain (DGTD); Finite Difference Time Domain (FDTD); Finite Element Method (FEM); Radar Cross Section (RCS)

    O441.4

    A

    2095-283X(2015)03-0361-06

    10.12000/JR15052

    楊謙, 魏兵, 李林茜, 等. DGTD用于RCS計算的初步研究[J]. 雷達學報, 2015, 4(3): 361–366.

    10.12000/JR15052.

    Reference format:Yang Qian, Wei Bing, Li Lin-qian,et al.. Preliminary research on RCS using DGTD[J].Journal of Radars, 2015, 4(3): 361–366. DOI: 10.12000/JR15052.

    2015-05-05收到,2015-06-16改回

    國家自然科學基金項目(61231003, 61401344)資助課題

    *通信作者: 魏兵 bwei@xidian.edu.cn

    猜你喜歡
    總場入射波四面體
    四面體小把戲
    SHPB入射波相似律與整形技術的試驗與數(shù)值研究
    振動與沖擊(2022年6期)2022-03-27 12:18:26
    R3中四面體的幾個新Bonnesen型不等式
    R3中四面體的Bonnesen型等周不等式
    綜合施策打好棉花田管“組合拳”
    瞬態(tài)激勵狀態(tài)下樁身速度以及樁身內(nèi)力計算
    前向雷達目標回波成分與特性分析
    石總場早播棉花出苗顯行
    對機械波半波損失現(xiàn)象的物理解釋
    電子科技(2015年11期)2015-03-06 01:32:24
    后彎管式波力發(fā)電裝置氣室結(jié)構(gòu)的試驗研究*
    成人国产av品久久久| 99久久综合免费| 亚洲精品乱码久久久久久按摩| 妹子高潮喷水视频| 国产精品女同一区二区软件| 中文字幕av成人在线电影| 80岁老熟妇乱子伦牲交| 99久国产av精品国产电影| xxx大片免费视频| 国产av精品麻豆| 国产色爽女视频免费观看| 日本欧美国产在线视频| xxx大片免费视频| 日韩人妻高清精品专区| 亚洲av中文字字幕乱码综合| 中文乱码字字幕精品一区二区三区| 777米奇影视久久| 午夜福利高清视频| 免费av不卡在线播放| 欧美 日韩 精品 国产| 欧美zozozo另类| 国产伦精品一区二区三区视频9| 人体艺术视频欧美日本| 国产乱人偷精品视频| 精品少妇久久久久久888优播| 黄片wwwwww| 不卡视频在线观看欧美| 蜜桃久久精品国产亚洲av| 国产高清有码在线观看视频| 欧美精品一区二区免费开放| 亚洲久久久国产精品| 啦啦啦啦在线视频资源| 久久99热这里只频精品6学生| 婷婷色综合大香蕉| 色综合色国产| 国产免费又黄又爽又色| 国产亚洲精品久久久com| 色吧在线观看| 国产有黄有色有爽视频| 最近最新中文字幕大全电影3| 偷拍熟女少妇极品色| freevideosex欧美| 99热这里只有是精品在线观看| 亚洲国产欧美在线一区| 涩涩av久久男人的天堂| 国产黄片视频在线免费观看| 看十八女毛片水多多多| 三级国产精品片| 久久久午夜欧美精品| 在线播放无遮挡| 久久99精品国语久久久| 国产无遮挡羞羞视频在线观看| 国产精品女同一区二区软件| 日韩不卡一区二区三区视频在线| 精品少妇黑人巨大在线播放| 亚洲精品成人av观看孕妇| 高清av免费在线| 亚洲av中文av极速乱| 国产中年淑女户外野战色| 国产亚洲最大av| 国产乱人视频| 秋霞伦理黄片| 国产 精品1| 国产精品爽爽va在线观看网站| 男人舔奶头视频| 91精品一卡2卡3卡4卡| 国产人妻一区二区三区在| 男女边吃奶边做爰视频| 国产永久视频网站| 亚洲人成网站高清观看| 国产人妻一区二区三区在| 多毛熟女@视频| 美女国产视频在线观看| 亚洲综合色惰| 久久99热这里只有精品18| 人人妻人人爽人人添夜夜欢视频 | 99热这里只有是精品50| 好男人视频免费观看在线| 日韩三级伦理在线观看| 男女啪啪激烈高潮av片| 成人美女网站在线观看视频| 亚洲人与动物交配视频| 亚洲成人中文字幕在线播放| 欧美人与善性xxx| 纯流量卡能插随身wifi吗| 久久99热这里只频精品6学生| 久久精品久久久久久噜噜老黄| 大香蕉久久网| 免费黄频网站在线观看国产| 熟女人妻精品中文字幕| 国产精品精品国产色婷婷| 五月玫瑰六月丁香| www.色视频.com| 99热网站在线观看| 免费不卡的大黄色大毛片视频在线观看| 亚洲一级一片aⅴ在线观看| 十分钟在线观看高清视频www | 3wmmmm亚洲av在线观看| 国产又色又爽无遮挡免| 亚洲av在线观看美女高潮| 国产人妻一区二区三区在| 视频区图区小说| 国产亚洲5aaaaa淫片| 内射极品少妇av片p| 日韩 亚洲 欧美在线| 99九九线精品视频在线观看视频| 中文精品一卡2卡3卡4更新| 日日啪夜夜撸| 亚洲精品自拍成人| 人妻一区二区av| 春色校园在线视频观看| 亚洲av电影在线观看一区二区三区| 亚洲怡红院男人天堂| 午夜福利在线观看免费完整高清在| 小蜜桃在线观看免费完整版高清| 婷婷色av中文字幕| 亚洲婷婷狠狠爱综合网| 色婷婷av一区二区三区视频| 亚洲av国产av综合av卡| 免费大片黄手机在线观看| 蜜桃久久精品国产亚洲av| 国产老妇伦熟女老妇高清| 国产欧美另类精品又又久久亚洲欧美| 欧美精品一区二区大全| 亚洲久久久国产精品| 狠狠精品人妻久久久久久综合| 99久久综合免费| 国产黄片视频在线免费观看| 一级片'在线观看视频| 久久久久精品久久久久真实原创| 久久青草综合色| 天堂8中文在线网| 日韩大片免费观看网站| 久久久久久伊人网av| 欧美亚洲 丝袜 人妻 在线| 黑丝袜美女国产一区| 日韩,欧美,国产一区二区三区| 人人妻人人爽人人添夜夜欢视频 | 建设人人有责人人尽责人人享有的 | 最后的刺客免费高清国语| 欧美精品国产亚洲| 国产有黄有色有爽视频| 亚洲国产精品999| 国产亚洲av片在线观看秒播厂| 日本黄色片子视频| av国产精品久久久久影院| 久久精品久久久久久噜噜老黄| 成人毛片a级毛片在线播放| 观看免费一级毛片| 人妻系列 视频| 亚洲国产最新在线播放| 国产高清国产精品国产三级 | 99热网站在线观看| 亚洲精品乱码久久久久久按摩| 大陆偷拍与自拍| 欧美人与善性xxx| 国产成人精品一,二区| 亚洲av综合色区一区| 欧美精品亚洲一区二区| 大片免费播放器 马上看| 国产综合精华液| 久久 成人 亚洲| 久久人人爽av亚洲精品天堂 | 久久韩国三级中文字幕| 啦啦啦在线观看免费高清www| 亚洲第一区二区三区不卡| 18+在线观看网站| 免费黄频网站在线观看国产| 极品教师在线视频| 麻豆成人午夜福利视频| 卡戴珊不雅视频在线播放| 一级片'在线观看视频| 成人黄色视频免费在线看| av在线蜜桃| 高清视频免费观看一区二区| 久久久a久久爽久久v久久| 久久人人爽人人片av| 超碰av人人做人人爽久久| 国产精品一区二区性色av| 国产一区二区三区综合在线观看 | 成人毛片60女人毛片免费| 80岁老熟妇乱子伦牲交| 啦啦啦啦在线视频资源| 九九爱精品视频在线观看| 精品亚洲成国产av| 成年女人在线观看亚洲视频| 国产精品久久久久久精品电影小说 | 91久久精品国产一区二区成人| 久久久精品94久久精品| 成人国产av品久久久| 成人漫画全彩无遮挡| 欧美日韩综合久久久久久| 内地一区二区视频在线| 欧美精品国产亚洲| 国产又色又爽无遮挡免| 亚洲va在线va天堂va国产| 亚洲av.av天堂| 日韩av不卡免费在线播放| 国产在线视频一区二区| 99热这里只有是精品50| 成年免费大片在线观看| 午夜老司机福利剧场| 高清av免费在线| 亚洲欧洲国产日韩| 亚洲电影在线观看av| 人人妻人人爽人人添夜夜欢视频 | 精品人妻熟女av久视频| 久久97久久精品| 亚洲av日韩在线播放| 夫妻性生交免费视频一级片| 亚洲电影在线观看av| 亚洲激情五月婷婷啪啪| 亚洲精品视频女| 免费看日本二区| 国产黄色视频一区二区在线观看| 国产国拍精品亚洲av在线观看| 高清黄色对白视频在线免费看 | 亚洲av成人精品一二三区| 全区人妻精品视频| av黄色大香蕉| 少妇 在线观看| 美女内射精品一级片tv| 美女xxoo啪啪120秒动态图| 十分钟在线观看高清视频www | 精华霜和精华液先用哪个| 婷婷色综合大香蕉| 日本色播在线视频| 久久国产乱子免费精品| 精品人妻偷拍中文字幕| 1000部很黄的大片| av福利片在线观看| 欧美人与善性xxx| 亚洲人成网站在线观看播放| 大又大粗又爽又黄少妇毛片口| 精品亚洲成a人片在线观看 | 国产欧美日韩精品一区二区| 国产深夜福利视频在线观看| 一区二区三区乱码不卡18| 亚洲内射少妇av| 丰满人妻一区二区三区视频av| 欧美bdsm另类| 视频区图区小说| 看非洲黑人一级黄片| 一区在线观看完整版| 亚洲av国产av综合av卡| freevideosex欧美| 成人毛片60女人毛片免费| 寂寞人妻少妇视频99o| 老司机影院毛片| 免费av不卡在线播放| 岛国毛片在线播放| 亚洲色图av天堂| 五月开心婷婷网| 性高湖久久久久久久久免费观看| 国产精品.久久久| 极品教师在线视频| 久久国产亚洲av麻豆专区| 在线观看三级黄色| 亚州av有码| 亚洲怡红院男人天堂| 久久99热这里只有精品18| 久久久久精品久久久久真实原创| 日本黄大片高清| 99精国产麻豆久久婷婷| 色网站视频免费| 人人妻人人澡人人爽人人夜夜| 大香蕉久久网| 一个人看视频在线观看www免费| 亚洲人成网站在线观看播放| 国产久久久一区二区三区| 欧美高清性xxxxhd video| 天堂俺去俺来也www色官网| 国产精品国产三级国产专区5o| 久久久久性生活片| 国产精品久久久久久精品电影小说 | 久久精品国产亚洲av涩爱| 伦精品一区二区三区| 女性生殖器流出的白浆| 成人亚洲欧美一区二区av| 小蜜桃在线观看免费完整版高清| 免费av不卡在线播放| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品国产成人久久av| 91精品伊人久久大香线蕉| 人妻少妇偷人精品九色| 直男gayav资源| 国产一区二区三区综合在线观看 | 亚洲国产精品999| 网址你懂的国产日韩在线| 91午夜精品亚洲一区二区三区| 亚洲av男天堂| 久久6这里有精品| 亚洲精品第二区| 乱码一卡2卡4卡精品| 伦精品一区二区三区| 一级爰片在线观看| 国产探花极品一区二区| 国产精品人妻久久久久久| 波野结衣二区三区在线| av免费在线看不卡| 国产人妻一区二区三区在| 又爽又黄a免费视频| 亚洲最大成人中文| av.在线天堂| 亚洲av成人精品一二三区| 久久影院123| 最近的中文字幕免费完整| 又黄又爽又刺激的免费视频.| 不卡视频在线观看欧美| 1000部很黄的大片| 高清日韩中文字幕在线| 91在线精品国自产拍蜜月| 久久久久性生活片| 在线观看一区二区三区激情| 看十八女毛片水多多多| 日本av免费视频播放| 日韩强制内射视频| 热99国产精品久久久久久7| 在线亚洲精品国产二区图片欧美 | 又大又黄又爽视频免费| 91久久精品国产一区二区三区| 丰满少妇做爰视频| 免费少妇av软件| 成人亚洲精品一区在线观看 | 亚洲国产av新网站| 嫩草影院新地址| 色视频www国产| 在线观看一区二区三区激情| 岛国毛片在线播放| 亚洲美女黄色视频免费看| 国产高清国产精品国产三级 | 欧美区成人在线视频| 亚洲精品日韩av片在线观看| 激情五月婷婷亚洲| 毛片女人毛片| 亚洲国产精品国产精品| 黄色一级大片看看| 亚洲精品,欧美精品| 激情五月婷婷亚洲| 五月天丁香电影| 国产成人精品一,二区| 看十八女毛片水多多多| 各种免费的搞黄视频| 国产精品久久久久久精品电影小说 | 老女人水多毛片| 亚洲欧洲日产国产| 亚洲,一卡二卡三卡| 少妇猛男粗大的猛烈进出视频| 嫩草影院入口| 一级毛片我不卡| 在线 av 中文字幕| 少妇 在线观看| 一级黄片播放器| 99热这里只有是精品在线观看| 99热全是精品| 啦啦啦啦在线视频资源| 特大巨黑吊av在线直播| 亚洲不卡免费看| 久久精品国产鲁丝片午夜精品| 国产一级毛片在线| 国产无遮挡羞羞视频在线观看| 国产91av在线免费观看| 99热这里只有精品一区| 内地一区二区视频在线| 久久久久久久久久久丰满| 乱码一卡2卡4卡精品| 精品一区二区三卡| 熟女av电影| 岛国毛片在线播放| 国产av码专区亚洲av| 高清黄色对白视频在线免费看 | 日韩欧美 国产精品| 自拍欧美九色日韩亚洲蝌蚪91 | 国产精品久久久久久av不卡| 国产爽快片一区二区三区| 国产成人精品福利久久| 在线观看一区二区三区| 美女福利国产在线 | 国产成人一区二区在线| 久久精品国产a三级三级三级| 久久99热6这里只有精品| 国产成人一区二区在线| 亚洲成人中文字幕在线播放| 久久青草综合色| 女的被弄到高潮叫床怎么办| 岛国毛片在线播放| 久久久久精品久久久久真实原创| 男女国产视频网站| 最近的中文字幕免费完整| 久久婷婷青草| 两个人的视频大全免费| 国产精品欧美亚洲77777| 最近的中文字幕免费完整| 国产又色又爽无遮挡免| 我的老师免费观看完整版| 国产亚洲91精品色在线| 亚洲欧美中文字幕日韩二区| 在线观看免费视频网站a站| 久久久亚洲精品成人影院| 97在线视频观看| a级毛色黄片| 国产精品99久久久久久久久| 在线看a的网站| 亚洲久久久国产精品| 欧美精品一区二区免费开放| 免费观看性生交大片5| 亚洲美女搞黄在线观看| av福利片在线观看| 中国三级夫妇交换| 2021少妇久久久久久久久久久| 亚洲三级黄色毛片| 国产永久视频网站| 天堂8中文在线网| 极品教师在线视频| 亚洲av二区三区四区| 黑人高潮一二区| 麻豆成人午夜福利视频| 国产精品国产av在线观看| 免费人妻精品一区二区三区视频| 日本av手机在线免费观看| 亚洲美女搞黄在线观看| 网址你懂的国产日韩在线| 国产精品av视频在线免费观看| 我要看黄色一级片免费的| 日韩电影二区| 亚洲精品,欧美精品| 最近中文字幕高清免费大全6| 亚洲美女搞黄在线观看| 亚洲第一区二区三区不卡| 丰满乱子伦码专区| 婷婷色av中文字幕| 国产高清不卡午夜福利| 十分钟在线观看高清视频www | 免费观看在线日韩| 亚洲欧洲日产国产| 成年免费大片在线观看| 三级经典国产精品| 十分钟在线观看高清视频www | 日产精品乱码卡一卡2卡三| 亚洲内射少妇av| 亚洲精品成人av观看孕妇| 久久久久国产网址| 黑人高潮一二区| 韩国av在线不卡| 性高湖久久久久久久久免费观看| av免费在线看不卡| 久久精品夜色国产| 久久鲁丝午夜福利片| 亚洲综合色惰| av女优亚洲男人天堂| 国产又色又爽无遮挡免| 亚洲精品,欧美精品| 在线观看美女被高潮喷水网站| 老熟女久久久| 国内少妇人妻偷人精品xxx网站| 国产精品女同一区二区软件| 国产爽快片一区二区三区| 免费观看在线日韩| 亚洲av电影在线观看一区二区三区| 十分钟在线观看高清视频www | 欧美激情国产日韩精品一区| 一区二区av电影网| 春色校园在线视频观看| 久久毛片免费看一区二区三区| 色婷婷av一区二区三区视频| 亚洲av中文av极速乱| av国产精品久久久久影院| 精品亚洲成a人片在线观看 | 国产亚洲最大av| 国产欧美日韩一区二区三区在线 | 毛片一级片免费看久久久久| 又黄又爽又刺激的免费视频.| 亚洲真实伦在线观看| 亚洲国产最新在线播放| 亚洲丝袜综合中文字幕| 亚洲熟女精品中文字幕| 亚洲精品国产av成人精品| 久久久精品免费免费高清| 精华霜和精华液先用哪个| 久久99热这里只频精品6学生| 国产高清三级在线| 2022亚洲国产成人精品| 国产av国产精品国产| 在线观看人妻少妇| 91午夜精品亚洲一区二区三区| 亚洲色图综合在线观看| 国产淫语在线视频| 深夜a级毛片| 五月天丁香电影| 99re6热这里在线精品视频| 久久久久视频综合| 国产久久久一区二区三区| 亚洲熟女精品中文字幕| 97精品久久久久久久久久精品| 亚洲欧美日韩东京热| 久久久久久久久大av| 久久精品夜色国产| 亚洲成人av在线免费| 久久久久人妻精品一区果冻| av在线app专区| 深夜a级毛片| 精品人妻熟女av久视频| 亚洲经典国产精华液单| 男女国产视频网站| 亚洲精品国产av蜜桃| 国产爱豆传媒在线观看| 亚洲高清免费不卡视频| 国产又色又爽无遮挡免| 七月丁香在线播放| 99热全是精品| 国产成人精品久久久久久| 久久 成人 亚洲| 精品一品国产午夜福利视频| 国产精品av视频在线免费观看| 韩国av在线不卡| 男女免费视频国产| 毛片一级片免费看久久久久| 熟女av电影| 国产成人91sexporn| 国产精品伦人一区二区| 欧美日韩一区二区视频在线观看视频在线| 国产精品一区二区在线不卡| 乱码一卡2卡4卡精品| av免费观看日本| 免费黄色在线免费观看| 久久国产精品大桥未久av | 色吧在线观看| 香蕉精品网在线| 亚洲av国产av综合av卡| 亚洲欧美清纯卡通| 国产精品.久久久| 色视频在线一区二区三区| 国产精品国产av在线观看| 大片电影免费在线观看免费| 乱码一卡2卡4卡精品| 婷婷色综合大香蕉| 日本黄色片子视频| 中文字幕av成人在线电影| 男男h啪啪无遮挡| 又黄又爽又刺激的免费视频.| 日韩欧美 国产精品| 激情 狠狠 欧美| 国产一级毛片在线| 一个人看的www免费观看视频| 国产精品国产三级国产av玫瑰| 国产片特级美女逼逼视频| 午夜免费男女啪啪视频观看| 国产午夜精品久久久久久一区二区三区| 免费av中文字幕在线| 99热6这里只有精品| 亚洲在久久综合| 精品久久久久久电影网| 免费观看的影片在线观看| 久久久久久久亚洲中文字幕| 国产乱人视频| 99热国产这里只有精品6| 男男h啪啪无遮挡| 少妇精品久久久久久久| 男女边摸边吃奶| 欧美zozozo另类| 99热这里只有精品一区| 国产高清不卡午夜福利| 啦啦啦在线观看免费高清www| 最近的中文字幕免费完整| 国产毛片在线视频| 搡女人真爽免费视频火全软件| 亚洲综合精品二区| 亚洲丝袜综合中文字幕| 日本av免费视频播放| 精品久久久精品久久久| 性高湖久久久久久久久免费观看| 久久久久久久精品精品| 高清视频免费观看一区二区| av一本久久久久| 久久这里有精品视频免费| 亚州av有码| av黄色大香蕉| 亚洲伊人久久精品综合| 纯流量卡能插随身wifi吗| 日本av免费视频播放| 建设人人有责人人尽责人人享有的 | 成人午夜精彩视频在线观看| 啦啦啦在线观看免费高清www| 老师上课跳d突然被开到最大视频| 菩萨蛮人人尽说江南好唐韦庄| 极品教师在线视频| 久久国产乱子免费精品| 亚洲av综合色区一区| 日本午夜av视频| 美女视频免费永久观看网站| 在线观看免费日韩欧美大片 | 高清午夜精品一区二区三区| 97热精品久久久久久| 成人亚洲欧美一区二区av| 久久韩国三级中文字幕| 十八禁网站网址无遮挡 | a级毛片免费高清观看在线播放| 亚洲丝袜综合中文字幕| 国国产精品蜜臀av免费| 中文字幕久久专区| 多毛熟女@视频| 久久久精品免费免费高清| 97在线视频观看| 国产精品人妻久久久久久| 国产视频首页在线观看| 中国美白少妇内射xxxbb| 国产视频首页在线观看| 最新中文字幕久久久久| 最黄视频免费看| 极品教师在线视频| 99热国产这里只有精品6| 欧美高清性xxxxhd video| 亚洲国产精品专区欧美| 亚洲国产精品国产精品| 一本色道久久久久久精品综合|