• <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
    一本色道久久久久久精品综合| av在线老鸭窝| 亚洲精品乱久久久久久| 亚洲精品aⅴ在线观看| 极品少妇高潮喷水抽搐| 久久精品久久精品一区二区三区| 一级a做视频免费观看| 亚洲一区二区三区欧美精品| 一级毛片aaaaaa免费看小| 少妇 在线观看| 国产精品三级大全| 女人精品久久久久毛片| 哪个播放器可以免费观看大片| 久久女婷五月综合色啪小说| 亚洲av免费高清在线观看| 男女国产视频网站| 午夜激情久久久久久久| 丝袜美足系列| 亚洲美女黄色视频免费看| 久久免费观看电影| 久久精品人人爽人人爽视色| 国产乱人偷精品视频| 婷婷色综合www| 美女xxoo啪啪120秒动态图| 午夜福利网站1000一区二区三区| 国产探花极品一区二区| 国产色婷婷99| 中文字幕av电影在线播放| 黑人高潮一二区| 街头女战士在线观看网站| 两个人的视频大全免费| 国产成人精品无人区| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美成人精品一区二区| 中文字幕制服av| 秋霞在线观看毛片| 国产黄色视频一区二区在线观看| 婷婷色综合www| 少妇人妻 视频| 久久精品人人爽人人爽视色| 精品一区二区三卡| 中文字幕制服av| 婷婷色av中文字幕| 最近手机中文字幕大全| 汤姆久久久久久久影院中文字幕| 国产精品人妻久久久久久| 久久久a久久爽久久v久久| 王馨瑶露胸无遮挡在线观看| 在线观看国产h片| 18禁在线播放成人免费| 国产男人的电影天堂91| 国产国语露脸激情在线看| 视频区图区小说| 男女边摸边吃奶| 免费观看av网站的网址| 成人亚洲精品一区在线观看| 国产成人freesex在线| 久久这里有精品视频免费| 在线免费观看不下载黄p国产| 9色porny在线观看| 午夜久久久在线观看| 午夜激情久久久久久久| 欧美精品一区二区大全| a级片在线免费高清观看视频| 丁香六月天网| 超碰97精品在线观看| 熟女人妻精品中文字幕| 高清欧美精品videossex| 久久99热这里只频精品6学生| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品日韩在线中文字幕| 亚洲人成77777在线视频| 亚洲成人av在线免费| 免费黄频网站在线观看国产| xxxhd国产人妻xxx| 国产成人精品婷婷| 99热国产这里只有精品6| 国产精品国产三级专区第一集| 久久狼人影院| 午夜精品国产一区二区电影| 午夜精品国产一区二区电影| 一区二区三区四区激情视频| 女人久久www免费人成看片| av黄色大香蕉| av天堂久久9| 97超视频在线观看视频| 日本av免费视频播放| 寂寞人妻少妇视频99o| 国产熟女欧美一区二区| .国产精品久久| 人体艺术视频欧美日本| 岛国毛片在线播放| 在线观看免费高清a一片| 18禁裸乳无遮挡动漫免费视频| 日韩中字成人| 国产成人91sexporn| 亚洲一级一片aⅴ在线观看| 22中文网久久字幕| 国产又色又爽无遮挡免| 国产成人精品一,二区| 欧美3d第一页| 国产免费视频播放在线视频| 能在线免费看毛片的网站| 亚洲四区av| 人妻一区二区av| 韩国av在线不卡| 天堂俺去俺来也www色官网| 大香蕉97超碰在线| 99九九线精品视频在线观看视频| 亚洲欧洲精品一区二区精品久久久 | 啦啦啦视频在线资源免费观看| 大码成人一级视频| 亚洲av中文av极速乱| 美女脱内裤让男人舔精品视频| 性高湖久久久久久久久免费观看| 欧美人与善性xxx| 久久久欧美国产精品| 免费人成在线观看视频色| 亚洲欧美日韩另类电影网站| 少妇猛男粗大的猛烈进出视频| 在线观看免费高清a一片| 狂野欧美激情性xxxx在线观看| 精品久久蜜臀av无| 亚洲国产av影院在线观看| 欧美bdsm另类| 街头女战士在线观看网站| 99国产综合亚洲精品| 久久久久国产网址| av视频免费观看在线观看| 国产色婷婷99| av在线老鸭窝| 精品人妻一区二区三区麻豆| 纯流量卡能插随身wifi吗| 国产69精品久久久久777片| 久久精品夜色国产| 在现免费观看毛片| 午夜老司机福利剧场| 最后的刺客免费高清国语| 亚洲欧洲日产国产| 国产欧美另类精品又又久久亚洲欧美| 大码成人一级视频| 97超视频在线观看视频| av免费观看日本| 26uuu在线亚洲综合色| 九九久久精品国产亚洲av麻豆| 99久久中文字幕三级久久日本| 国产成人精品婷婷| 一级a做视频免费观看| 桃花免费在线播放| 亚洲综合色网址| 国产日韩一区二区三区精品不卡 | 国产男女内射视频| 国产精品国产三级专区第一集| 日韩av在线免费看完整版不卡| 制服丝袜香蕉在线| 亚洲av国产av综合av卡| 国产一区二区三区av在线| 少妇丰满av| 男女啪啪激烈高潮av片| 日韩中字成人| 日本wwww免费看| 国产一区亚洲一区在线观看| xxx大片免费视频| 2018国产大陆天天弄谢| 大片电影免费在线观看免费| 午夜免费鲁丝| 一区二区三区乱码不卡18| 2021少妇久久久久久久久久久| 99久久精品国产国产毛片| 国产精品女同一区二区软件| 丰满迷人的少妇在线观看| 国产亚洲欧美精品永久| 激情五月婷婷亚洲| 久久午夜福利片| 亚洲欧美中文字幕日韩二区| 亚洲欧美日韩另类电影网站| av女优亚洲男人天堂| 国产老妇伦熟女老妇高清| av在线播放精品| 夜夜骑夜夜射夜夜干| 色婷婷久久久亚洲欧美| 中文字幕精品免费在线观看视频 | 又粗又硬又长又爽又黄的视频| 欧美国产精品一级二级三级| 18禁在线播放成人免费| 高清不卡的av网站| 婷婷色综合www| 久久国内精品自在自线图片| 简卡轻食公司| 青春草国产在线视频| 国产精品久久久久成人av| 亚洲四区av| 丝袜脚勾引网站| 亚洲五月色婷婷综合| 天堂俺去俺来也www色官网| 久久综合国产亚洲精品| 99久久中文字幕三级久久日本| 亚洲精品乱久久久久久| 国国产精品蜜臀av免费| 9色porny在线观看| 高清不卡的av网站| 丰满少妇做爰视频| 亚洲av免费高清在线观看| 狂野欧美白嫩少妇大欣赏| 欧美丝袜亚洲另类| 成年av动漫网址| 中文字幕久久专区| 一级爰片在线观看| 久久精品熟女亚洲av麻豆精品| 少妇被粗大的猛进出69影院 | 午夜福利视频精品| 国语对白做爰xxxⅹ性视频网站| 欧美精品亚洲一区二区| 国产爽快片一区二区三区| 精品人妻在线不人妻| 精品国产露脸久久av麻豆| 久久99蜜桃精品久久| 99热国产这里只有精品6| 黑人高潮一二区| 国精品久久久久久国模美| 欧美精品亚洲一区二区| 亚洲精品一区蜜桃| 中文字幕人妻熟人妻熟丝袜美| 在线免费观看不下载黄p国产| 欧美成人精品欧美一级黄| 亚洲国产精品国产精品| 久久av网站| 久久免费观看电影| av免费在线看不卡| 国产成人免费观看mmmm| 久久人人爽人人爽人人片va| 国产黄片视频在线免费观看| av国产久精品久网站免费入址| 免费观看的影片在线观看| 蜜桃国产av成人99| 久热这里只有精品99| 夜夜看夜夜爽夜夜摸| 街头女战士在线观看网站| 亚洲国产av影院在线观看| 色视频在线一区二区三区| 亚洲欧美清纯卡通| 欧美 日韩 精品 国产| 亚洲人成网站在线播| 亚洲,欧美,日韩| 久久久久久久久久久久大奶| 亚洲精品色激情综合| 久久 成人 亚洲| 一本一本综合久久| 中文乱码字字幕精品一区二区三区| 久久精品久久久久久久性| 日韩电影二区| 全区人妻精品视频| 中文欧美无线码| 国产成人免费观看mmmm| 久久午夜综合久久蜜桃| 最近手机中文字幕大全| 婷婷色av中文字幕| 日韩精品免费视频一区二区三区 | 国产免费一区二区三区四区乱码| 九色亚洲精品在线播放| 亚洲第一av免费看| 免费看不卡的av| 午夜福利网站1000一区二区三区| 国产精品成人在线| 成人影院久久| 大片免费播放器 马上看| 51国产日韩欧美| 欧美成人午夜免费资源| 精品视频人人做人人爽| 岛国毛片在线播放| 男的添女的下面高潮视频| 少妇被粗大猛烈的视频| 久久久国产一区二区| 国产成人a∨麻豆精品| 国产老妇伦熟女老妇高清| 日本欧美国产在线视频| 蜜桃国产av成人99| 成人午夜精彩视频在线观看| 夫妻性生交免费视频一级片| 亚洲人成网站在线观看播放| .国产精品久久| 亚洲精品美女久久av网站| 一级,二级,三级黄色视频| 色哟哟·www| 老司机影院成人| 黑人巨大精品欧美一区二区蜜桃 | 国产日韩欧美亚洲二区| 久久女婷五月综合色啪小说| 美女主播在线视频| 在线观看国产h片| 国产老妇伦熟女老妇高清| 婷婷成人精品国产| 免费少妇av软件| 全区人妻精品视频| 亚洲国产av影院在线观看| 在线观看免费高清a一片| 婷婷色av中文字幕| 在线播放无遮挡| av.在线天堂| 黄片无遮挡物在线观看| 国产精品三级大全| 一本一本综合久久| 欧美成人午夜免费资源| 九草在线视频观看| 一本—道久久a久久精品蜜桃钙片| 日日摸夜夜添夜夜爱| 亚洲av不卡在线观看| 久久鲁丝午夜福利片| 日韩制服骚丝袜av| 最后的刺客免费高清国语| 蜜桃在线观看..| 亚洲综合色惰| 国产乱来视频区| 高清视频免费观看一区二区| 蜜臀久久99精品久久宅男| 亚洲av电影在线观看一区二区三区| 久久精品国产亚洲网站| 美女中出高潮动态图| 国产免费福利视频在线观看| 大片电影免费在线观看免费| 亚洲国产精品成人久久小说| 亚洲精品久久成人aⅴ小说 | av免费观看日本| 欧美性感艳星| 狂野欧美白嫩少妇大欣赏| 一个人免费看片子| 成人亚洲欧美一区二区av| 菩萨蛮人人尽说江南好唐韦庄| 大片免费播放器 马上看| 一个人免费看片子| 人体艺术视频欧美日本| 亚洲欧美成人综合另类久久久| 精品少妇久久久久久888优播| 亚洲精品第二区| 亚洲丝袜综合中文字幕| 在线看a的网站| 大码成人一级视频| 国产亚洲精品久久久com| 老女人水多毛片| 免费高清在线观看日韩| 日本wwww免费看| 午夜福利视频精品| 伊人久久国产一区二区| 久久女婷五月综合色啪小说| 免费播放大片免费观看视频在线观看| 99九九线精品视频在线观看视频| 亚洲国产精品成人久久小说| 亚洲欧美色中文字幕在线| 国产男女超爽视频在线观看| 观看av在线不卡| 国产综合精华液| 人妻系列 视频| 少妇被粗大的猛进出69影院 | 亚洲国产成人一精品久久久| 伦精品一区二区三区| 国产极品天堂在线| 国产亚洲午夜精品一区二区久久| 免费观看的影片在线观看| 国产在线一区二区三区精| 九草在线视频观看| 午夜av观看不卡| 成人18禁高潮啪啪吃奶动态图 | 成年av动漫网址| 久久久国产精品麻豆| 最近最新中文字幕免费大全7| 五月天丁香电影| kizo精华| 这个男人来自地球电影免费观看 | www.色视频.com| 国产亚洲精品久久久com| 国产av一区二区精品久久| 亚洲美女黄色视频免费看| av在线播放精品| 亚洲人成网站在线播| 久久久久久久久久人人人人人人| 亚洲av福利一区| 在线免费观看不下载黄p国产| 欧美日本中文国产一区发布| 亚洲怡红院男人天堂| 午夜91福利影院| 新久久久久国产一级毛片| 免费高清在线观看日韩| 国产精品免费大片| 91精品国产国语对白视频| 欧美3d第一页| 国产乱人偷精品视频| 国产在线视频一区二区| 欧美日韩成人在线一区二区| 青春草亚洲视频在线观看| 国产精品一二三区在线看| 免费高清在线观看日韩| 久久人人爽人人爽人人片va| 建设人人有责人人尽责人人享有的| 亚洲性久久影院| 日本色播在线视频| 日日爽夜夜爽网站| 极品少妇高潮喷水抽搐| 满18在线观看网站| 制服人妻中文乱码| 亚洲精品国产av蜜桃| 黑丝袜美女国产一区| 精品久久国产蜜桃| 亚洲国产成人一精品久久久| 国产探花极品一区二区| 丝袜喷水一区| 永久免费av网站大全| av一本久久久久| 久久久久久久国产电影| 寂寞人妻少妇视频99o| 制服人妻中文乱码| 免费黄网站久久成人精品| 高清午夜精品一区二区三区| 一级a做视频免费观看| 国产69精品久久久久777片| 五月开心婷婷网| 亚洲国产av影院在线观看| 国产精品嫩草影院av在线观看| 18+在线观看网站| 26uuu在线亚洲综合色| 亚洲精品日本国产第一区| 人妻一区二区av| 五月伊人婷婷丁香| 亚洲综合精品二区| 日韩电影二区| 欧美日韩在线观看h| 18禁在线无遮挡免费观看视频| 国产成人精品一,二区| 蜜臀久久99精品久久宅男| 黄色毛片三级朝国网站| 老女人水多毛片| 高清欧美精品videossex| 精品亚洲乱码少妇综合久久| 日本免费在线观看一区| 国国产精品蜜臀av免费| 蜜桃久久精品国产亚洲av| 在线看a的网站| 99九九线精品视频在线观看视频| 欧美日韩亚洲高清精品| 国产熟女欧美一区二区| 一级a做视频免费观看| 韩国高清视频一区二区三区| 女人久久www免费人成看片| 久久久久视频综合| 久久久久国产精品人妻一区二区| 欧美老熟妇乱子伦牲交| 99热国产这里只有精品6| 国产精品久久久久成人av| 国产精品国产av在线观看| 日韩中字成人| 一本大道久久a久久精品| 热99久久久久精品小说推荐| 成人国产麻豆网| 人妻少妇偷人精品九色| www.av在线官网国产| 亚洲综合色网址| 乱人伦中国视频| 满18在线观看网站| 亚洲精品av麻豆狂野| 国产日韩欧美在线精品| 蜜臀久久99精品久久宅男| 嫩草影院入口| 赤兔流量卡办理| 国产欧美日韩综合在线一区二区| 高清欧美精品videossex| 伊人久久精品亚洲午夜| 日本猛色少妇xxxxx猛交久久| 乱码一卡2卡4卡精品| 国产午夜精品一二区理论片| av福利片在线| 欧美精品亚洲一区二区| 免费大片黄手机在线观看| 九色成人免费人妻av| 亚洲不卡免费看| 久久久亚洲精品成人影院| 亚洲高清免费不卡视频| a级片在线免费高清观看视频| 日日撸夜夜添| 国产精品久久久久久久久免| 国产熟女欧美一区二区| 日本爱情动作片www.在线观看| 亚洲三级黄色毛片| 亚洲怡红院男人天堂| 欧美精品亚洲一区二区| videosex国产| 亚洲av电影在线观看一区二区三区| tube8黄色片| 69精品国产乱码久久久| 中文乱码字字幕精品一区二区三区| 亚洲国产日韩一区二区| 夫妻午夜视频| 91成人精品电影| av免费观看日本| 久久久久国产精品人妻一区二区| 精品熟女少妇av免费看| 日日摸夜夜添夜夜爱| 亚洲人成网站在线观看播放| 色吧在线观看| 热re99久久精品国产66热6| 肉色欧美久久久久久久蜜桃| 午夜福利视频精品| 尾随美女入室| 少妇的逼好多水| 一个人免费看片子| 日韩不卡一区二区三区视频在线| 99热全是精品| 插逼视频在线观看| av国产精品久久久久影院| 两个人免费观看高清视频| 七月丁香在线播放| 亚洲国产欧美日韩在线播放| 日韩人妻高清精品专区| 婷婷色综合www| 欧美激情极品国产一区二区三区 | 91aial.com中文字幕在线观看| 水蜜桃什么品种好| 交换朋友夫妻互换小说| av国产久精品久网站免费入址| 亚洲久久久国产精品| 国产精品欧美亚洲77777| av在线老鸭窝| 日韩伦理黄色片| 欧美丝袜亚洲另类| 男女边摸边吃奶| 国产片特级美女逼逼视频| 99精国产麻豆久久婷婷| 校园人妻丝袜中文字幕| 18+在线观看网站| 日韩中文字幕视频在线看片| 婷婷色av中文字幕| 国产永久视频网站| 性色avwww在线观看| 少妇被粗大的猛进出69影院 | 亚洲,一卡二卡三卡| 精品99又大又爽又粗少妇毛片| 人人妻人人澡人人爽人人夜夜| 亚洲欧美成人综合另类久久久| 王馨瑶露胸无遮挡在线观看| 老女人水多毛片| 交换朋友夫妻互换小说| 搡老乐熟女国产| 美女视频免费永久观看网站| 国产白丝娇喘喷水9色精品| www.色视频.com| 99久国产av精品国产电影| 在线观看免费视频网站a站| 亚洲精品乱码久久久久久按摩| 男女高潮啪啪啪动态图| 国产亚洲一区二区精品| 国产欧美亚洲国产| 亚洲国产av新网站| 亚洲欧美清纯卡通| 日韩精品免费视频一区二区三区 | 亚洲色图综合在线观看| 黑人巨大精品欧美一区二区蜜桃 | 嘟嘟电影网在线观看| 成人免费观看视频高清| 免费不卡的大黄色大毛片视频在线观看| 欧美精品亚洲一区二区| 日本-黄色视频高清免费观看| 丝袜美足系列| 高清视频免费观看一区二区| 欧美少妇被猛烈插入视频| 成年女人在线观看亚洲视频| 自拍欧美九色日韩亚洲蝌蚪91| 一级爰片在线观看| 一本色道久久久久久精品综合| 中文字幕人妻熟人妻熟丝袜美| 国产精品秋霞免费鲁丝片| 欧美成人午夜免费资源| 如日韩欧美国产精品一区二区三区 | 啦啦啦在线观看免费高清www| 国产成人精品一,二区| av福利片在线| 中文字幕亚洲精品专区| 3wmmmm亚洲av在线观看| 国产欧美日韩一区二区三区在线 | 亚洲欧美日韩另类电影网站| 精品人妻偷拍中文字幕| 国产精品一区www在线观看| 亚洲av日韩在线播放| 大香蕉97超碰在线| av.在线天堂| 日韩亚洲欧美综合| 在线看a的网站| 亚洲欧美中文字幕日韩二区| av网站免费在线观看视频| 国产精品99久久99久久久不卡 | 人人妻人人澡人人爽人人夜夜| 啦啦啦视频在线资源免费观看| 欧美激情国产日韩精品一区| 天堂8中文在线网| 18+在线观看网站| 久久午夜综合久久蜜桃| 亚洲精品一二三| 日韩强制内射视频| 在线播放无遮挡| 久久人人爽av亚洲精品天堂| 日韩在线高清观看一区二区三区| 男女边吃奶边做爰视频| 精品人妻在线不人妻| 激情五月婷婷亚洲| 91精品一卡2卡3卡4卡| 免费大片18禁| 欧美激情国产日韩精品一区| 久久久久久久精品精品| 伊人久久精品亚洲午夜| 久久97久久精品| 亚洲经典国产精华液单| 精品久久久精品久久久| 午夜福利在线观看免费完整高清在| 一本一本综合久久| 人成视频在线观看免费观看| 亚洲人与动物交配视频| av在线观看视频网站免费| 日本黄色片子视频|