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

    電特大天線增益的高效求解

    2016-11-25 05:48:11孫旭敏潘小敏盛新慶
    北京理工大學(xué)學(xué)報 2016年7期
    關(guān)鍵詞:子層多極平面波

    孫旭敏,潘小敏,盛新慶

    (北京理工大學(xué) 信息與電子學(xué)院,電磁仿真中心,北京 100081)

    ?

    電特大天線增益的高效求解

    孫旭敏,潘小敏,盛新慶

    (北京理工大學(xué) 信息與電子學(xué)院,電磁仿真中心,北京 100081)

    為解決電特大天線輻射特性計算效率低的問題,提出并實現(xiàn)了一種基于多極子的聚集思想快速求解遠場信息的高效并行算法. 根據(jù)目標(biāo)在并行多層快速多極子中幾何樹的特征,提出高層遠場模式直接聚集,并且按行并行插值的混合MPI-OpenMP并行計算方案. 將該算法應(yīng)用于大射電天文望遠鏡FAST(five-hundred-meter aperture spherical telescope)的遠場方向圖和增益的計算,通過與等效電流直接積分法計算對比展示其計算精度與效率. 結(jié)果驗證該算法的遠場輻射特性計算精度,且同等遠場計算任務(wù)量下,平均單方向計算效率提高近萬倍.

    并行多層快速多極子;MPI并行;OpenMP并行;電特大;增益

    隨著計算機并行技術(shù)的發(fā)展,為了滿足大目標(biāo)尤其是電特大目標(biāo)特性的精確計算對高效算法的要求,基于集群系統(tǒng)的消息傳遞(message passage interface,MPI)模型的多進程并行技術(shù)以及多核處理器的共享內(nèi)存(OpenMP)多線程并行技術(shù)的應(yīng)用[1-5],并行MLFMA對目標(biāo)散射特性的計算能力已達到未知數(shù)為十億數(shù)量級.

    文獻[6]中利用并行MLFMA來分析500 m口徑射電天文望遠鏡(five-hundred-meter aperture spherical telescope,F(xiàn)AST)的輻射特性,此天線主瓣電平高波束窄,為了保證其計算精度,該文獻的做法是在最大增益附近處dθ=±10°內(nèi)計算更密的遠場信息,這樣需要較大遠場信息量的耗時計算,而且方向數(shù)量的多少直接影響其精度,如此人為的因素就需要考慮,在精度和效率之間需要找到一個平衡點,對于程序的使用者來說這無疑是比較繁瑣的.

    對于上述問題,文獻[7]中已實現(xiàn)其計算目標(biāo)散射的MPI并行算法的研究,該方法與電流直接積分求解相比,計算電大目標(biāo)效率提高很客觀. 為解決應(yīng)用目前的MPI-OpenMP并行架構(gòu)MLFMA計算電特大輻射目標(biāo)效率低的問題,采用該遠場計算方法,其實現(xiàn)過程中平面波聚集方案的選取與實現(xiàn)則直接影響算法的計算精度與效率. 多極子幾何樹結(jié)構(gòu)最高層(盒子尺寸最大的層)具有特殊性,計算輻射問題時需要重新考慮負載平衡的實現(xiàn),經(jīng)分析高層遠場模式不同的聚集方案所需的計算資源,包括按分層逐次/單層直接聚集以及改進型的按行/按列并行插值方案,選取最優(yōu)的MPI并行方案,并重點致力于對OpenMP并行的研究,最終實現(xiàn)高效匹配的MPI-OpenMP并行方案,將其應(yīng)用于求解電特大輻射問題,在保證所需內(nèi)存最小化的情況下,實現(xiàn)了對其輻射特性的精確高效計算.

    1 基于并行MLFMA高效求解輻射問題

    與并行MFLMA計算散射問題相比,應(yīng)用其計算輻射問題的主要區(qū)別就是饋源和計算指標(biāo)不同. 饋源部分選用簡化的電壓源模型,該模型采用δ函數(shù)源,饋源端口處的RWG邊作為激勵邊來初始化矩陣方程的右端項;在計算天線輻射特性時,以方向性系數(shù)近似作為天線軸向增益,其數(shù)值計算公式如下:

    (1)

    (2)

    此外,并行MLFMA中應(yīng)用多極子思想處理遠相互作用中,通過各基函數(shù)聚集而得的遠場模式公式為

    (3)

    對比上述兩式,可知計算天線方向性增益時所需的全向遠場模值可通過計算多極子中結(jié)合所有基函數(shù)貢獻的遠場模式,通過乘以常數(shù)因子可得,來提高計算效率.

    2 基于多極子快速求解增益的并行

    基于并行MLFMA的多極子思想計算輻射問題,實現(xiàn)重點就是應(yīng)用多極子的聚集思想,根據(jù)等效電流并行求解天線的遠場信息. 基于MPI-OpenMP并行架構(gòu)的遠相互作用處理的聚集、轉(zhuǎn)移操作高效所得的第2層遠場模式,已包含所有未知數(shù)的貢獻,因此需適當(dāng)選擇實現(xiàn)高層(第2層到第0層即目標(biāo)本身)遠場模式的聚集方案,并通過適當(dāng)選取零層截斷系數(shù)以保證高斯面積分方法所對應(yīng)的遠場模式方向數(shù)可滿足精度要求. 并對聚集過程所需計算資源進行分析驗證以確定最優(yōu)方案,其中插值系數(shù)矩陣的計算資源由父層和子層遠場模式數(shù)量決定,平移系數(shù)由父層與子層盒子空間關(guān)系及父層遠場模式?jīng)Q定,其實現(xiàn)中涉及的主要問題:高層聚集方式選取——分層逐次聚集與單層直接聚集;聚集過程的MPI并行方案——按盒子并行與按平面波并行;不同MPI并行方案下OpenMP并行實現(xiàn)的負載均衡問題.

    目標(biāo)總體遠場模式的計算,可以由第二層遠場模式通過分層逐次聚集與單層直接聚集兩種方式來實現(xiàn)(圖1中2D結(jié)構(gòu)樹). 而對于電大輻射問題高層的遠場模式數(shù)較多,單層聚集在構(gòu)造插值系數(shù)與平移系數(shù)及聚集過程中節(jié)省大量的計算資源,主要考察其與分層聚集的計算精度. 通過以sph-0242與sph-10(金屬球半徑電長度分別為2.42,10.00)的散射問題為例,比較兩種方案相對于電流直接積分法和解析所得的RCS結(jié)果,這兩種方案相對于解析解的相對誤差基本<5%,而誤差較大的方向均由于RCS接近于0 dB,并且直接聚集和電流積分的誤差相差不大,而對于其他數(shù)值算例以及電大目標(biāo)的散射計算也驗證單層直接聚集方案可以保證計算精度要求,故選取該方案,并在此基礎(chǔ)上討論并行方案的選取與實現(xiàn).

    2.1 基于多極子快速求解增益的MPI并行方案

    基于并行MLFMA采用的是高層平均分配平面波,低層平均分配盒子的混合MPI并行方案,由過渡層(iTrsLev)與第2層關(guān)系決定高層聚集過程的MPI并行方案. 當(dāng)iTrsLev<2,高層按盒子并行,除了進程間盒子相互依賴關(guān)系外,其插值和平移過程與原盒子并行類似,每層的全部遠場模式方向數(shù)同樣作用于所有進程;當(dāng)iTrsLev≥2,高層按平面波并行,針對大部分目標(biāo)尤其是電大目標(biāo),其高層遠場模式數(shù)量較多均為平面波并行,且這些層盒子尺寸大數(shù)量少,當(dāng)以盒子并行時隨著線程數(shù)增加其效率下降. 并行MLFMA的層間插值系數(shù)存儲方式,按盒子并行時以父層平面波作為負載,其關(guān)于子層平面波的插值系數(shù)共享,采用按所有父層平面波即按行存儲方式;按平面波并行則以該進程父層平面波作為負載,分別保存其關(guān)于所有對其有貢獻的所有子層平面波的插值系數(shù),采用按這些子層平面波即按列存儲方式,以此方案來節(jié)省計算資源. 而高層結(jié)構(gòu)樹的特殊性(父層即第0層只有一個盒子,子層即第2層所有盒子均是其子盒子)使得在采用相應(yīng)并行方案的同時,需要對其聚集過程中計算資源進行估計并改進算法,包括構(gòu)造相應(yīng)插值以及平移矩陣的內(nèi)存及時間,以及聚集過程中循環(huán)次序的調(diào)整,以實現(xiàn)線程并行的負載均衡及并行效率的提高.

    按盒子并行方案的插值系數(shù)依舊采用原始按行插值形式,其總平面波方向數(shù)較少,聚集后采用歸約操作(MPI_Reduce)對完整的遠場再進行數(shù)據(jù)后處理部分;按平面波并行方案是本文重點,由于高層結(jié)構(gòu)樹的特殊性以及此時父層遠場模式數(shù)較多,需要改進其并行聚集方案,其中對于插值系數(shù)矩陣、聚集過程以及數(shù)據(jù)通信部分的并行實現(xiàn)進行探討及驗證.

    其中分塊插值系數(shù)矩陣均為稀疏矩陣,改進型按行插值與按列插值的并行插值系數(shù)矩陣選取及存儲方案如圖2(a)所示,改進的按行方案先以父層所有平面波進行搜索,若對其插值有貢獻的子層平面波(數(shù)目≤16)有在本地的則進行保存,否則將不保存,則本進程所有被貢獻的父層平面波都將被保存,其中包含與之左右相鄰進程的部分平面波,以這些局部父層平面波為一維索引,得到對應(yīng)的子層平面波與插值系數(shù),相對于原按行方案,于父層平面波的局部性以節(jié)省內(nèi)存需求;圖2(b)所示,改進的按列方案為減少通信開銷,以本地子層平面波為貢獻項,即以本進程子層平面波為一維索引,得到其所有有貢獻的父層平面波與插值系數(shù). 相對于改進型的插值方案,平移過程均以父層局部平面波作為整體負載,按子層全部盒子進行循環(huán),平移之后按盒子貢獻疊加得到局部平面波. 從本質(zhì)來看,兩種插值方式的信息一致,子層都無需進行通信,父層均得到局部平面波信息,不過其相對應(yīng)的聚集過程中插值方案選取將直接影響并行計算資源,于線程并行中詳細討論.

    對于電特大輻射問題,遠場模式聚集方法的零層方向由截斷系數(shù)決定,該方向數(shù)很大,區(qū)別于散射問題對于遠場信息的歸約處理,為實現(xiàn)節(jié)省計算資源,需要對聚集得到的局部平面波進行通信(MPI_Sendrecv)得到按進程并行的平面波信息,即將局部平面波不屬于本進程的部分平面波發(fā)送給對應(yīng)進程,該進程對應(yīng)疊加接收到的平面波來進行處理. 最終將各進程的遠場模式通過數(shù)據(jù)收集(MPI_Gatherv)操作將主極化和交叉極化的場值均集中于主進程(0號進程),此時得到整個三維空間的全向遠場,通過天線增益計算公式即可得增益.

    2.2 基于多極子快速求解增益的OpenMP加速技術(shù)

    高層聚集應(yīng)用基于共享內(nèi)存的指導(dǎo)性調(diào)度線程并行OpenMP來加快計算,為了充分實現(xiàn)其最佳并行效率,需先確定計算密集區(qū),并保證達到負載均衡的任務(wù)量來實現(xiàn)線程同步. 高層聚集包括插值系數(shù)矩陣填充、方向向量填充、平面波聚集、遠場模式數(shù)據(jù)處理這四部分,其中平面波聚集最耗時,在數(shù)據(jù)后處理總時間中占60%多,且其遠場方向數(shù)足夠大可保證并行的任務(wù)量,因此對按行與按列兩種聚集方案的OpenMP并行實現(xiàn)均進行分析和討論.

    按行插值聚集中,單線程方案中以子層盒子作為循環(huán)對象,插值系數(shù)以父層平面波作為一維索引,插值過程以父層局部平面波作為整體進行. 線程并行化的實現(xiàn),首先需要確定負載任務(wù)的并行分配:由于父層為單盒子,子層盒子間為獨立計算,子層盒子數(shù)(3D最大為64)數(shù)目有限,而且根據(jù)插值系數(shù)的特征,父層平面波可作為獨立計算負載,此時需要調(diào)整循環(huán)次序,以父層局部平面波為一階循環(huán)來進行負載匹配,子層盒子為二階循環(huán),再應(yīng)用已獨立保存的平面波方向和盒子位置信息來進行平移操作. 因此應(yīng)用按行插值聚集方案通過循環(huán)次序的調(diào)整不僅可按最大任務(wù)量實現(xiàn)最佳并行,而且減少了平移的計算時間.

    按列插值聚集中,區(qū)別在于其插值系數(shù)的一維索引是子層平面波,故無法通過循環(huán)次序的變換將父層平面波負載作為線程匹配對象,因此考慮到線程并行共享內(nèi)存的模式與子層盒子的獨立性,只能仍以子層盒子數(shù)作為并行的循環(huán)對象,且需要對將寫入的父層局部平面波分配與線程數(shù)一致的內(nèi)存空間,獨立存儲,聚集完成后再以各線程結(jié)果求和來實現(xiàn)子層貢獻的疊加,該方案由于并行對象有限以及求和計算增加,并行效率明顯隨著線程數(shù)增加而下降,且對內(nèi)存的要求增大.

    由于父層平面波數(shù)目較大,為了節(jié)省內(nèi)存,采用了原遠相互作用中的內(nèi)存循環(huán)技術(shù),若結(jié)構(gòu)體中遠場模式除目前子層平面波后內(nèi)存足夠則無需另開內(nèi)存,而按列插值聚集并行過程中該內(nèi)存需求與線程數(shù)成正比,無法保證在提高并行計算效率的同時不提高內(nèi)存的消耗. 因此基于并行效率和內(nèi)存需求的多方面分析,最優(yōu)的并行方案為結(jié)合MPI-OpenMP并行技術(shù)的改進型按行插值聚集方案.

    3 算例的數(shù)值分析

    計算飛機散射問題來測試其遠場計算的線程并行效率,電尺寸為240λ,未知數(shù)為450萬,根據(jù)節(jié)點資源,線程數(shù)從2~12變化,其高層按平面波并行時線程并行效率如圖3,當(dāng)線程到達12時,并行效率仍高于83%,該并行策略較成功實現(xiàn)了并行加速.

    計算反射面天線FAST-300(反射面光滑,口徑為300 m)的輻射特性,與電流直接積分法計算結(jié)果比較來驗證本文并行算法的精度和效率. 該反射面天線于焦點以天文臺提供的階梯喇叭作為饋源,分析頻率為300 MHz,天線電尺寸為300λ,結(jié)構(gòu)樹層數(shù)為10層,平均剖分尺寸為0.13λ,總未知數(shù)為1 737萬,用64個進程6個線程并行計算,E-plane遠場方向圖結(jié)果如圖4,計算增益為58.19 dBi,與電流直接積分法和商業(yè)軟件FEKO(PO法)所得的58.28 dBi和58.22 dBi相比,其誤差均小于1‰,輻射特性結(jié)果吻合較好.

    其中兩種方法所需的計算資源如表1,快速的遠場聚集方案在遠場計算中應(yīng)用了遠場模式的內(nèi)存循環(huán)技術(shù),不過保存高層插值與平移系數(shù)需要額外的內(nèi)存,其中電流直接積分法為保證最大輻射方向的精度,仍計算兩次遠場信息,表中統(tǒng)計一次計算時間,方向選取(θ,φ)步進為1°,快速遠場計算中方向以高斯取點決定,其平均遠場計算時間提高約為1.6萬倍,總時間亦如預(yù)期提高,因此若計算結(jié)構(gòu)更復(fù)雜或精度要求更高的輻射問題時,同等計算任務(wù)量下該算法更具優(yōu)勢.

    表1 FAST-300天線計算資源

    4 結(jié) 論

    為實現(xiàn)應(yīng)用并行MLFMA高效準(zhǔn)確地計算電特大輻射問題,采用了快速的多極子遠場聚集求解方案. 基于并行結(jié)構(gòu)樹的高層特殊性,對于高層不同聚集方案的資源需求進行分析,并通過算例驗證,確定以改進型按行插值的單層直接聚集作為最優(yōu)方案,為進一步提高計算效率,于計算密集區(qū)通過循環(huán)次序調(diào)整結(jié)合線程并行,并應(yīng)用了內(nèi)存循環(huán)技術(shù)來節(jié)省額外內(nèi)存.

    通過計算理想金屬目標(biāo)的散射問題,驗證其遠場聚集的線程并行效率,在線程數(shù)為12時依然高達83%的并行效率;以電特大反射面天線FAST的計算,驗證了其對于多方向遠場需求計算的高效性,與電流直接積分法相比,輻射特性精度和效率均得到顯著提高,數(shù)據(jù)后處理部分將同等計算任務(wù)量條件下,效率提高近萬倍.

    [1] Song J M,Chew W C. Multilevel fast-multipole algorithm for solving combined field integral equations of electromagnetic scattering [J].Microwave and Optical Technology Letter,1995,10(1):14-19.

    [2] 潘小敏,盛新慶.一種多層快速多極子的高效并行方案[J].電子學(xué)報,2007,35(3):567-571.

    Pan Xiaomin,Sheng Xinqing. A highly efficient parallel approach of multi-level fast multipole algorithm[J]. Acta Electronica Sinica,2007,35(3):567-571. (in Chinese)

    [3] 潘小敏,盛新慶. 一種高性能并行多層快速多極子算法[J].電子學(xué)報,2010,38(3):580-584.

    Pan Xiaomin,Sheng Xinqing. A high-performance parallel multi-level fast multipole algorithm[J]. Acta Electronica Sinica,2010,38(3):580-584. (in Chinese)

    [4] 潘小敏,皮維超,盛新慶.基于共享內(nèi)存的高效OpenMP并行多層快速多極子算法[J].北京理工大學(xué)學(xué)報,2012,32(2):164-169.

    Pan Xiaomin,Pi Weichao,Sheng Xinqing. Efficient parallelization of multi-level fast multipole algorithm based on OpenMP[J]. Transactions of Beijing Institute of Technology,2012,32(2):164-169. (in Chinese)

    [5] Pan Xiaomin,Pi Weichao,Yang Minglin,et al. Solving problems with over one billion unknowns by the MLFMA[J]. IEEE Transactions on Antennas and Propagation,2012,60(5):2571-2574.

    [6] 段猛,李斌,盛新慶,等.基于并行多層快速多極子的大射電天文望遠鏡FAST的計算[C]∥2011年全國天線年會.南京:[s.n.],2011:666-669.

    Duan Meng,Li Bin,Sheng Xinqing,et al. The computation of large spherical radio telescope FAST antenna based on the parallel MLFMA[C]∥Proceedings of The 2011 National Conference on Antennas. Nanjing: [s.n.], 2011:666-669. (in Chinese)

    [7] Velamparambil S,Chew W C,Song J. 10 million unknowns: is it that big?[J]. IEEE Antennas and Propagation Magazine,2003,45(2):43-58.

    (責(zé)任編輯:劉芳)

    Efficient Solution for Extremely Electric Large Antenna Gain

    SUN Xu-min,PAN Xiao-min,SHENG Xin-qing

    (Center for Electromagnetic Simulation,School of Information and Electronics,Beijing Institute of Technology,Beijing 100081,China)

    An efficient parallel algorithm to accelerate the calculation of the far-fields was proposed and implemented for solving the electric large antenna radiation problem,based on the upward pass or the aggregation phase of the parallel multilevel fast multipole algorithm(MLFMA). According to the geometry tree structure of the parallel MLFMA,an efficient hybrid MPI-OpenMP parallelization algorithm was designed with the direct aggregation of coarser level far-field radiation pattern and parallel in-row interpolation. The computation of far-field radiation pattern and antenna gain of the extremely large parabolic reflector of FAST by this algorithm shows the accuracy and efficiency of the parallel program,compared with the computation by the equivalent current direct integration method. The numerical results demonstrate that the presented approach has good accuracy for the radiation characteristics,and the computational efficiency has improved about ten thousand times for single direction under same far-field calculation task compared to the equivalent current direct integration method.

    parallel multi-level fast multi-pole algorithm (parallel MLFMA);MPI parallel;OpenMP parallel;extremely electric large;gain

    2014-03-17

    國家“九七三”計劃項目(2012CB720702)

    孫旭敏(1989—),女,博士生,E-mail: sunxm1989@foxmail.com.

    盛新慶(1968—),男,教授,博士生導(dǎo)師,E-mail: xsheng@bit.edu.cn.

    O 441.4

    A

    1001-0645(2016)07-0718-05

    10.15918/j.tbit1001-0645.2016.07.011

    猜你喜歡
    子層多極平面波
    Landau-Lifshitz方程平面波解的全局光滑性
    5G OTA測量寬帶平面波模擬器的高效優(yōu)化方法與應(yīng)用
    復(fù)合材料厚層合板力學(xué)性能等效方法研究
    基于GPU并行運算的超聲平面波成像仿真
    電子制作(2016年11期)2016-11-07 08:43:45
    分離載荷近距協(xié)議的數(shù)據(jù)鏈路層仿真研究
    應(yīng)用于艦載天線電磁兼容性分析的改進多層快速多極子算法
    多極子陣列聲波測井在煤層氣勘探開發(fā)中的應(yīng)用
    中國煤層氣(2014年3期)2014-08-07 03:07:39
    以太網(wǎng)協(xié)議模型的演進分析
    積極為構(gòu)建“多點多極”戰(zhàn)略提供咨詢服務(wù)
    平面波展開法在聲子晶體中的應(yīng)用
    科技資訊(2012年19期)2012-04-29 08:23:50
    少妇熟女aⅴ在线视频| 中文字幕高清在线视频| 欧美日韩乱码在线| 欧美一区二区国产精品久久精品 | 天堂√8在线中文| 久久久精品欧美日韩精品| 久久久久久大精品| 搡老妇女老女人老熟妇| 亚洲五月婷婷丁香| 午夜精品在线福利| 国产97色在线日韩免费| 国产精品 欧美亚洲| 人妻丰满熟妇av一区二区三区| 一级毛片女人18水好多| 最近最新中文字幕大全电影3| 色尼玛亚洲综合影院| 怎么达到女性高潮| 国产熟女xx| 欧美日韩瑟瑟在线播放| 窝窝影院91人妻| 国产成人精品久久二区二区免费| 欧美绝顶高潮抽搐喷水| 一二三四在线观看免费中文在| 人妻夜夜爽99麻豆av| 国产亚洲精品第一综合不卡| 草草在线视频免费看| 五月玫瑰六月丁香| 99国产综合亚洲精品| 可以在线观看的亚洲视频| 亚洲欧美日韩高清在线视频| 色老头精品视频在线观看| 亚洲色图av天堂| 欧美日韩亚洲国产一区二区在线观看| 国产熟女午夜一区二区三区| 中出人妻视频一区二区| 日日摸夜夜添夜夜添小说| 亚洲专区国产一区二区| 18禁观看日本| 美女午夜性视频免费| 国产成年人精品一区二区| 亚洲男人天堂网一区| 19禁男女啪啪无遮挡网站| 国产精品精品国产色婷婷| 中国美女看黄片| 女人被狂操c到高潮| 亚洲欧美日韩高清在线视频| 成年版毛片免费区| 97碰自拍视频| 国产三级在线视频| 一区福利在线观看| 亚洲国产欧洲综合997久久,| 91av网站免费观看| 欧美成人一区二区免费高清观看 | 99久久无色码亚洲精品果冻| 成人av在线播放网站| 精品欧美一区二区三区在线| 99热这里只有是精品50| 91九色精品人成在线观看| 在线观看免费视频日本深夜| 嫩草影视91久久| 国产精品av视频在线免费观看| 哪里可以看免费的av片| 午夜福利视频1000在线观看| 91老司机精品| 久久久久久免费高清国产稀缺| 国产私拍福利视频在线观看| 欧美成狂野欧美在线观看| 国产在线精品亚洲第一网站| 身体一侧抽搐| 国产精品九九99| 午夜精品久久久久久毛片777| 人妻丰满熟妇av一区二区三区| 国产一级毛片七仙女欲春2| 深夜精品福利| 成人av在线播放网站| 国产精品野战在线观看| av天堂在线播放| 黄色视频,在线免费观看| 亚洲欧美激情综合另类| 欧美日韩国产亚洲二区| 啦啦啦免费观看视频1| 99国产精品一区二区蜜桃av| 高清在线国产一区| 精品不卡国产一区二区三区| 88av欧美| 国产麻豆成人av免费视频| 看黄色毛片网站| 成人国产一区最新在线观看| 极品教师在线免费播放| 婷婷丁香在线五月| 一二三四在线观看免费中文在| 亚洲激情在线av| a在线观看视频网站| 久久热在线av| 国产爱豆传媒在线观看 | 日韩国内少妇激情av| 神马国产精品三级电影在线观看 | 亚洲专区国产一区二区| 男女床上黄色一级片免费看| 麻豆国产av国片精品| 欧美最黄视频在线播放免费| 亚洲精品一区av在线观看| 我要搜黄色片| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品一卡2卡三卡4卡5卡| 两个人免费观看高清视频| 久久精品综合一区二区三区| 成人永久免费在线观看视频| 激情在线观看视频在线高清| 最好的美女福利视频网| 久久久国产欧美日韩av| 久久久久久久久久黄片| 亚洲无线在线观看| 中国美女看黄片| 精品久久蜜臀av无| 一级毛片精品| 女同久久另类99精品国产91| 亚洲真实伦在线观看| 久久精品影院6| 亚洲精品久久国产高清桃花| 成年免费大片在线观看| 精品国产乱子伦一区二区三区| 亚洲色图av天堂| av中文乱码字幕在线| www.999成人在线观看| a级毛片在线看网站| 美女 人体艺术 gogo| 老司机福利观看| 日本免费a在线| 麻豆av在线久日| 人成视频在线观看免费观看| 精品无人区乱码1区二区| 亚洲精品av麻豆狂野| 国产午夜福利久久久久久| 国产不卡一卡二| АⅤ资源中文在线天堂| 久久精品成人免费网站| 特大巨黑吊av在线直播| 日日爽夜夜爽网站| 老司机靠b影院| 亚洲精品美女久久av网站| tocl精华| 国产单亲对白刺激| 国产又黄又爽又无遮挡在线| 老熟妇乱子伦视频在线观看| 搡老岳熟女国产| 女生性感内裤真人,穿戴方法视频| 色综合站精品国产| 亚洲成av人片免费观看| 三级毛片av免费| 国产亚洲欧美98| 久久国产精品人妻蜜桃| 嫩草影视91久久| 中文字幕av在线有码专区| 色尼玛亚洲综合影院| 99国产精品99久久久久| 欧美一区二区精品小视频在线| 男女床上黄色一级片免费看| 1024手机看黄色片| 天堂影院成人在线观看| 国产亚洲av高清不卡| 99久久精品热视频| 好男人在线观看高清免费视频| 中文字幕最新亚洲高清| 国产亚洲av高清不卡| 黄片小视频在线播放| 欧美乱妇无乱码| 中文字幕久久专区| 国产精品久久久久久亚洲av鲁大| 一个人免费在线观看电影 | 一本综合久久免费| 亚洲av成人不卡在线观看播放网| 91成年电影在线观看| 桃色一区二区三区在线观看| 日韩精品中文字幕看吧| 1024香蕉在线观看| 日韩有码中文字幕| 黑人欧美特级aaaaaa片| 精品午夜福利视频在线观看一区| 久久中文字幕人妻熟女| 757午夜福利合集在线观看| 亚洲中文av在线| 国产三级在线视频| 在线观看舔阴道视频| 中文字幕久久专区| 久久久国产成人免费| 亚洲狠狠婷婷综合久久图片| 熟女电影av网| 亚洲成a人片在线一区二区| 久久久久亚洲av毛片大全| 女人爽到高潮嗷嗷叫在线视频| 少妇的丰满在线观看| 欧美日韩福利视频一区二区| 午夜精品久久久久久毛片777| 最近最新免费中文字幕在线| 亚洲精品久久国产高清桃花| 毛片女人毛片| av福利片在线| 在线看三级毛片| 亚洲avbb在线观看| 黄色a级毛片大全视频| 国产成人精品久久二区二区免费| 国产午夜精品论理片| 亚洲午夜精品一区,二区,三区| 国语自产精品视频在线第100页| 欧美精品亚洲一区二区| 亚洲va日本ⅴa欧美va伊人久久| 免费在线观看成人毛片| 欧美日韩中文字幕国产精品一区二区三区| 十八禁网站免费在线| 亚洲熟妇中文字幕五十中出| 亚洲自拍偷在线| 天堂av国产一区二区熟女人妻 | 欧美日韩亚洲国产一区二区在线观看| www.精华液| 色av中文字幕| 最新美女视频免费是黄的| 日日夜夜操网爽| 一本综合久久免费| 亚洲一区二区三区色噜噜| 午夜福利欧美成人| 欧美性猛交╳xxx乱大交人| 一进一出抽搐动态| 可以在线观看毛片的网站| 又紧又爽又黄一区二区| 无遮挡黄片免费观看| 女生性感内裤真人,穿戴方法视频| 在线观看免费午夜福利视频| 麻豆av在线久日| 久久精品国产99精品国产亚洲性色| 国产单亲对白刺激| 中文字幕最新亚洲高清| 香蕉久久夜色| 久久久久久九九精品二区国产 | 麻豆成人午夜福利视频| 亚洲第一电影网av| 国产伦一二天堂av在线观看| 老司机靠b影院| 丰满的人妻完整版| 岛国在线免费视频观看| 久久久久性生活片| 高潮久久久久久久久久久不卡| 在线a可以看的网站| 在线永久观看黄色视频| a级毛片在线看网站| 老司机午夜福利在线观看视频| 成人精品一区二区免费| 亚洲 欧美 日韩 在线 免费| 97碰自拍视频| 久久久久久免费高清国产稀缺| 欧美国产日韩亚洲一区| 国产亚洲av高清不卡| 国产精品久久久人人做人人爽| cao死你这个sao货| 女人被狂操c到高潮| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲欧洲精品一区二区精品久久久| 老司机在亚洲福利影院| 欧美午夜高清在线| 这个男人来自地球电影免费观看| 成年女人毛片免费观看观看9| 日日爽夜夜爽网站| 欧美性长视频在线观看| 亚洲精品在线观看二区| 亚洲中文av在线| tocl精华| www.熟女人妻精品国产| 日韩国内少妇激情av| 亚洲精品一区av在线观看| 欧美一区二区国产精品久久精品 | 91老司机精品| 午夜两性在线视频| 久久99热这里只有精品18| 欧美在线一区亚洲| 亚洲精品国产一区二区精华液| 色老头精品视频在线观看| 啦啦啦免费观看视频1| 黄色视频,在线免费观看| 久久久久精品国产欧美久久久| 日韩欧美国产一区二区入口| 欧美成人一区二区免费高清观看 | 午夜免费成人在线视频| 制服诱惑二区| 亚洲欧美精品综合一区二区三区| 亚洲精品中文字幕一二三四区| 精品国内亚洲2022精品成人| 一进一出好大好爽视频| 午夜福利成人在线免费观看| 女人高潮潮喷娇喘18禁视频| 久9热在线精品视频| 婷婷精品国产亚洲av| a级毛片a级免费在线| 国产aⅴ精品一区二区三区波| 色在线成人网| 国产麻豆成人av免费视频| 三级男女做爰猛烈吃奶摸视频| 亚洲性夜色夜夜综合| xxx96com| 久久精品国产清高在天天线| 人妻久久中文字幕网| 午夜福利在线在线| 久久香蕉激情| 亚洲国产高清在线一区二区三| 99国产精品一区二区蜜桃av| 在线观看舔阴道视频| 国产av不卡久久| 亚洲第一电影网av| 久久亚洲精品不卡| 亚洲第一欧美日韩一区二区三区| 神马国产精品三级电影在线观看 | 男女午夜视频在线观看| 三级毛片av免费| 国产一级毛片七仙女欲春2| 母亲3免费完整高清在线观看| 国产亚洲精品综合一区在线观看 | 香蕉久久夜色| 舔av片在线| 日日夜夜操网爽| 在线免费观看的www视频| 蜜桃久久精品国产亚洲av| 精品人妻1区二区| 香蕉国产在线看| 一级黄色大片毛片| 人成视频在线观看免费观看| 搞女人的毛片| 国产精品av久久久久免费| 欧美最黄视频在线播放免费| 88av欧美| 制服人妻中文乱码| 18禁国产床啪视频网站| 丰满人妻熟妇乱又伦精品不卡| 18禁黄网站禁片免费观看直播| 国产成人一区二区三区免费视频网站| 999久久久国产精品视频| 桃红色精品国产亚洲av| 麻豆成人av在线观看| 亚洲av第一区精品v没综合| 亚洲国产欧美网| 久久精品影院6| 国产精品一区二区免费欧美| 1024香蕉在线观看| 久久国产精品人妻蜜桃| 免费在线观看完整版高清| 法律面前人人平等表现在哪些方面| 欧美成人午夜精品| 黄片小视频在线播放| 波多野结衣高清作品| 一二三四社区在线视频社区8| 亚洲五月婷婷丁香| 怎么达到女性高潮| 日韩中文字幕欧美一区二区| 黄色成人免费大全| av天堂在线播放| 欧美日韩中文字幕国产精品一区二区三区| xxxwww97欧美| 曰老女人黄片| АⅤ资源中文在线天堂| 亚洲一区高清亚洲精品| 久久国产精品影院| 国产精品爽爽va在线观看网站| 亚洲av电影不卡..在线观看| 欧美日韩黄片免| 搡老妇女老女人老熟妇| 亚洲人成网站高清观看| 婷婷精品国产亚洲av在线| 免费看十八禁软件| 夜夜看夜夜爽夜夜摸| 国产一区二区三区在线臀色熟女| 亚洲avbb在线观看| 神马国产精品三级电影在线观看 | 看片在线看免费视频| 久久这里只有精品中国| 国产97色在线日韩免费| 午夜福利在线在线| 首页视频小说图片口味搜索| 天天躁狠狠躁夜夜躁狠狠躁| 精品欧美一区二区三区在线| 少妇粗大呻吟视频| a级毛片在线看网站| a级毛片在线看网站| 正在播放国产对白刺激| 波多野结衣高清无吗| 男女视频在线观看网站免费 | 一级黄色大片毛片| 听说在线观看完整版免费高清| 欧洲精品卡2卡3卡4卡5卡区| 悠悠久久av| 国产成人精品无人区| 国内精品久久久久精免费| 丁香欧美五月| 国产精华一区二区三区| 国产成人精品久久二区二区免费| 少妇被粗大的猛进出69影院| 亚洲欧美激情综合另类| 婷婷亚洲欧美| 美女高潮喷水抽搐中文字幕| 国内精品久久久久久久电影| 成人国产综合亚洲| 长腿黑丝高跟| 亚洲一区高清亚洲精品| 日韩中文字幕欧美一区二区| 国产亚洲精品久久久久5区| 国内精品久久久久精免费| 亚洲成人中文字幕在线播放| 国产精品美女特级片免费视频播放器 | 搡老妇女老女人老熟妇| 欧美丝袜亚洲另类 | a在线观看视频网站| 欧美一级a爱片免费观看看 | 亚洲人成电影免费在线| 亚洲午夜精品一区,二区,三区| 国产高清视频在线观看网站| 免费在线观看成人毛片| 国产精品一区二区三区四区免费观看 | 亚洲精品中文字幕在线视频| 两个人视频免费观看高清| 免费在线观看日本一区| 欧美黑人精品巨大| x7x7x7水蜜桃| 夜夜看夜夜爽夜夜摸| 国产高清有码在线观看视频 | 国产亚洲精品第一综合不卡| 亚洲成av人片免费观看| 手机成人av网站| 露出奶头的视频| 夜夜躁狠狠躁天天躁| 久久婷婷人人爽人人干人人爱| 在线看三级毛片| 欧美一区二区精品小视频在线| 国产亚洲av嫩草精品影院| 黄色成人免费大全| 午夜亚洲福利在线播放| 亚洲欧美日韩东京热| 久久久久久久精品吃奶| 夜夜看夜夜爽夜夜摸| 最近在线观看免费完整版| 国产精品综合久久久久久久免费| av片东京热男人的天堂| 国产一区在线观看成人免费| 久久午夜综合久久蜜桃| 亚洲中文av在线| 香蕉久久夜色| 亚洲五月天丁香| 国产视频一区二区在线看| 青草久久国产| 神马国产精品三级电影在线观看 | 宅男免费午夜| 国产99久久九九免费精品| 麻豆国产97在线/欧美 | 床上黄色一级片| 中文字幕av在线有码专区| 国产精品久久久久久久电影 | 天天添夜夜摸| 特大巨黑吊av在线直播| 伊人久久大香线蕉亚洲五| www.www免费av| 一级毛片高清免费大全| 床上黄色一级片| 欧美一级a爱片免费观看看 | 色综合亚洲欧美另类图片| 久久久久久大精品| 久久国产精品人妻蜜桃| 国产私拍福利视频在线观看| 亚洲欧美精品综合久久99| 国产男靠女视频免费网站| 国产成人av教育| ponron亚洲| 成年人黄色毛片网站| 熟女电影av网| 最近最新中文字幕大全免费视频| 黑人欧美特级aaaaaa片| 欧美黑人精品巨大| 中出人妻视频一区二区| 欧美不卡视频在线免费观看 | 久久久水蜜桃国产精品网| 久久久久久久午夜电影| 91大片在线观看| 日韩成人在线观看一区二区三区| 国产成人aa在线观看| 国产精品一区二区精品视频观看| 少妇被粗大的猛进出69影院| 天堂影院成人在线观看| 美女高潮喷水抽搐中文字幕| 最近最新免费中文字幕在线| 国产午夜精品论理片| 九色成人免费人妻av| 最新美女视频免费是黄的| 精品久久久久久,| 99精品久久久久人妻精品| 三级男女做爰猛烈吃奶摸视频| 51午夜福利影视在线观看| 一本综合久久免费| 国产精品av久久久久免费| 亚洲成av人片在线播放无| 欧美日韩中文字幕国产精品一区二区三区| 淫妇啪啪啪对白视频| 99热只有精品国产| 精品第一国产精品| 国产私拍福利视频在线观看| 高清在线国产一区| 一区福利在线观看| 亚洲,欧美精品.| 国产精品 欧美亚洲| 色噜噜av男人的天堂激情| 成年女人毛片免费观看观看9| 精品久久久久久,| 欧美中文日本在线观看视频| 国产三级在线视频| 性欧美人与动物交配| 2021天堂中文幕一二区在线观| 免费在线观看黄色视频的| 亚洲精品在线观看二区| 国内久久婷婷六月综合欲色啪| 亚洲精品国产精品久久久不卡| 成人永久免费在线观看视频| 黄片大片在线免费观看| 国产精品久久久久久亚洲av鲁大| svipshipincom国产片| 国产精品九九99| 岛国在线观看网站| 欧美人与性动交α欧美精品济南到| 天天躁狠狠躁夜夜躁狠狠躁| 两个人看的免费小视频| 天天添夜夜摸| 国产视频一区二区在线看| 久久这里只有精品19| 国产高清激情床上av| 啦啦啦观看免费观看视频高清| 极品教师在线免费播放| 一边摸一边做爽爽视频免费| xxxwww97欧美| 久久中文看片网| 91九色精品人成在线观看| 在线观看午夜福利视频| 日本五十路高清| 级片在线观看| 极品教师在线免费播放| 国产精品久久久久久久电影 | 日本免费a在线| 亚洲 欧美 日韩 在线 免费| 国产区一区二久久| 国产欧美日韩一区二区三| 亚洲av中文字字幕乱码综合| 日本 欧美在线| 国产爱豆传媒在线观看 | 国产高清视频在线观看网站| 特级一级黄色大片| 亚洲色图 男人天堂 中文字幕| 美女大奶头视频| 精品第一国产精品| 制服丝袜大香蕉在线| 在线观看www视频免费| 性欧美人与动物交配| 亚洲18禁久久av| 成人永久免费在线观看视频| 欧美成狂野欧美在线观看| 97人妻精品一区二区三区麻豆| 午夜精品一区二区三区免费看| 狠狠狠狠99中文字幕| 成人av一区二区三区在线看| 国产激情久久老熟女| 亚洲av日韩精品久久久久久密| 一二三四在线观看免费中文在| 一进一出抽搐动态| 国产精品久久久久久精品电影| 91在线观看av| 久久性视频一级片| 五月伊人婷婷丁香| 丁香六月欧美| 麻豆av在线久日| av天堂在线播放| 男女那种视频在线观看| 狠狠狠狠99中文字幕| 国产高清视频在线播放一区| 欧美极品一区二区三区四区| 欧美一区二区国产精品久久精品 | 男女之事视频高清在线观看| 亚洲成人久久性| 蜜桃久久精品国产亚洲av| 90打野战视频偷拍视频| 中文亚洲av片在线观看爽| 国产亚洲欧美98| 国产精品一区二区三区四区久久| 亚洲va日本ⅴa欧美va伊人久久| 别揉我奶头~嗯~啊~动态视频| 亚洲美女视频黄频| 亚洲午夜理论影院| 国产亚洲av嫩草精品影院| 巨乳人妻的诱惑在线观看| 黄色毛片三级朝国网站| 亚洲五月天丁香| 国产野战对白在线观看| 国产精品一区二区免费欧美| 岛国在线免费视频观看| 国产三级中文精品| 正在播放国产对白刺激| 搡老妇女老女人老熟妇| 999久久久国产精品视频| 亚洲欧洲精品一区二区精品久久久| 我要搜黄色片| 欧美一区二区精品小视频在线| 在线播放国产精品三级| 午夜亚洲福利在线播放| 日本免费一区二区三区高清不卡| 级片在线观看| 久久精品国产亚洲av高清一级| 999久久久国产精品视频| 一本精品99久久精品77| а√天堂www在线а√下载| 国产精品综合久久久久久久免费| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩乱码在线| 999久久久国产精品视频| 91大片在线观看| 欧美zozozo另类|