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

    利用并行八叉樹(shù)劃分技術(shù)開(kāi)展復(fù)雜構(gòu)型航天器放氣污染分析研究

    2017-01-11 07:23:21焦子龍姜利祥黃建國(guó)孫繼鵬朱云飛
    航天器環(huán)境工程 2016年6期
    關(guān)鍵詞:八叉樹(shù)構(gòu)型航天器

    焦子龍,姜利祥,黃建國(guó),孫繼鵬,朱云飛

    (北京衛(wèi)星環(huán)境工程研究所 可靠性與環(huán)境工程技術(shù)重點(diǎn)實(shí)驗(yàn)室,北京 100094)

    利用并行八叉樹(shù)劃分技術(shù)開(kāi)展復(fù)雜構(gòu)型航天器放氣污染分析研究

    焦子龍,姜利祥,黃建國(guó),孫繼鵬,朱云飛

    (北京衛(wèi)星環(huán)境工程研究所 可靠性與環(huán)境工程技術(shù)重點(diǎn)實(shí)驗(yàn)室,北京 100094)

    放氣污染對(duì)航天器造成的有害影響需要定量評(píng)估。放氣產(chǎn)物分子的輸運(yùn)過(guò)程可以采用 Monte-Carlo方法模擬,并用射線表示放氣分子輸運(yùn)軌跡。如此一來(lái),即轉(zhuǎn)換為射線追蹤問(wèn)題。由于航天器構(gòu)型復(fù)雜,射線追蹤的計(jì)算很耗時(shí),所以有必要研究其加速算法。文章將復(fù)雜構(gòu)型航天器表面用三角形非結(jié)構(gòu)網(wǎng)格表示,利用八叉樹(shù)劃分技術(shù)開(kāi)展分子運(yùn)動(dòng)軌跡與表面相交的加速模擬計(jì)算。同時(shí),在基于共享內(nèi)存多核計(jì)算平臺(tái),利用OpenMP API軟件實(shí)現(xiàn)了計(jì)算的并行化。該方法在對(duì)空間站的實(shí)例計(jì)算中顯著提高了運(yùn)算速度,并對(duì)空間站的污染情況進(jìn)行了初步評(píng)價(jià)分析。

    航天器;放氣污染;模擬計(jì)算;復(fù)雜構(gòu)型;八叉樹(shù);并行計(jì)算

    0 引言

    在航天器結(jié)構(gòu)中,廣泛應(yīng)用碳纖維/環(huán)氧復(fù)合材料、熱控涂層材料、隔熱材料、膠粘劑等[1-3]。航天器在軌運(yùn)行過(guò)程中,這些材料在真空、高低溫、粒子輻射等環(huán)境中將發(fā)生質(zhì)量損失,統(tǒng)稱為放氣[4]。這不僅影響材料本身的性能,還會(huì)對(duì)航天器造成污染。例如歐洲的Rosetta彗星探測(cè)器在利用四極質(zhì)譜計(jì)探測(cè)彗星周圍氣體成分時(shí),發(fā)現(xiàn)探測(cè)器材料放氣對(duì)探測(cè)結(jié)果產(chǎn)生了嚴(yán)重影響[5]。Cassini探測(cè)器在軌飛行中材料放氣使其窄視角相機(jī)受到嚴(yán)重污染,造成圖像出現(xiàn)光暈現(xiàn)象,不得不進(jìn)行了長(zhǎng)時(shí)間的加熱除氣[6]。因此,有必要對(duì)航天器用材料提出嚴(yán)格的放氣控制要求。

    除需對(duì)材料放氣進(jìn)行嚴(yán)格控制外,還需要定量評(píng)估放氣帶來(lái)的影響。國(guó)內(nèi)外對(duì)放氣污染評(píng)估分析方法進(jìn)行了廣泛研究。例如,針對(duì)國(guó)際空間站的污染分析與評(píng)估,ESA和 NASA分別開(kāi)發(fā)了COMOVA和NASAN污染分析軟件[7]。這2款軟件均采用Monte-Carlo方法,并用射線來(lái)表征分子的輸運(yùn)過(guò)程,包括分子從放氣源表面出射,在航天器周圍運(yùn)動(dòng)、碰撞,在敏感表面被吸附等過(guò)程。

    航天器構(gòu)型復(fù)雜,射線追蹤的計(jì)算很耗時(shí),因此有必要研究其加速算法。本文將計(jì)算機(jī)圖形學(xué)中常用到的射線追蹤加速方法即并行八叉樹(shù)劃分技術(shù)引入復(fù)雜構(gòu)型航天器的放氣污染分析中,并在OpenMP API軟件平臺(tái)上開(kāi)展了模擬分析,獲得了很好的加速效果。

    1 放氣污染計(jì)算

    1.1 蒙特卡羅方法

    粒子在軌道環(huán)境中傳輸分為直接流和背景散射返回流,對(duì)這些粒子流可以采用試驗(yàn)粒子蒙特卡羅方法(Test Particle Monte-Carlo, TPMC)進(jìn)行模擬。TPMC方法基于硬球碰撞(二體彈性碰撞)模型,采用統(tǒng)計(jì)方法對(duì)大量放氣分子的運(yùn)動(dòng)、碰撞和軌跡進(jìn)行模擬計(jì)算,以獲得相應(yīng)的散射分布和返回流比率,具體計(jì)算方法和步驟可參見(jiàn)文獻(xiàn)[8-10]。

    1.2 八叉樹(shù)劃分算法

    八叉樹(shù)(octa-tree)算法是一種常用的空間劃分技術(shù),具有數(shù)據(jù)結(jié)構(gòu)簡(jiǎn)單、集合運(yùn)算高效、節(jié)點(diǎn)遍歷運(yùn)算易于實(shí)現(xiàn)等優(yōu)點(diǎn),是射線追蹤普遍采用的加速方法,如圖1所示。

    圖1 八叉樹(shù)空間劃分示意圖Fig.1 Sketch of octa-tree subdivision

    其基本過(guò)程是:構(gòu)建包含所有空間對(duì)象的長(zhǎng)方體包圍盒作為根節(jié)點(diǎn);再對(duì)該包圍盒進(jìn)行細(xì)分,即沿3個(gè)坐標(biāo)軸方向?qū)⑵浞譃?個(gè)大小相等的子立方體;接著檢查每個(gè)子立方體,如果子立方體內(nèi)空間對(duì)象的個(gè)數(shù)不大于給定的閾值,則停止分解該子立方體,否則,需要對(duì)該子立方體繼續(xù)進(jìn)行8等分,直到內(nèi)部空間對(duì)象數(shù)低于給定閾值。

    采用八叉樹(shù)劃分后,當(dāng)需要求解射線與空間對(duì)象的交點(diǎn)時(shí),不需要判斷射線與每個(gè)空間對(duì)象是否相交,而只需對(duì)與射線相交的包圍盒中的空間對(duì)象進(jìn)行判斷(即先對(duì)包圍盒進(jìn)行簡(jiǎn)單的求交運(yùn)算,然后略過(guò)與射線不相交的包圍盒),從而達(dá)到減少求交次數(shù)和計(jì)算時(shí)間的目的。

    Agate等人[11]提出了一種自頂向下的八叉樹(shù)搜索算法,充分利用了八叉樹(shù)的空間結(jié)構(gòu)特性,與其他算法相比,計(jì)算速度更快。

    射線可用參數(shù)方程表示,即

    起點(diǎn)矢量為O=[Ox,Oy,Oz]T,方向矢量為D=[l,m,n]T。在判斷射線是否與坐標(biāo)軸平行的八叉樹(shù)子節(jié)點(diǎn)相交之前,可先求出射線與子節(jié)點(diǎn)邊界相交的t值:

    其中x0和x1分別為子節(jié)點(diǎn)的x坐標(biāo)值的下邊界和上邊界。

    類似地可以求出0yt、1yt、0zt、1zt。然后再求出tmin和tmax:

    當(dāng)沿八叉樹(shù)從根節(jié)點(diǎn)向下搜索時(shí),可以利用空間劃分的特性加速計(jì)算。八叉樹(shù)子節(jié)點(diǎn)劃分采用中間平面,因此其子節(jié)點(diǎn)的任一方向的t值恰好是對(duì)應(yīng)父節(jié)點(diǎn)t值的一半。然后,確定射線進(jìn)入節(jié)點(diǎn)的平面,再根據(jù)射線的參數(shù)確定后續(xù)的相交節(jié)點(diǎn)。如果該節(jié)點(diǎn)包含空間對(duì)象,則將其加入待求交點(diǎn)空間對(duì)象列表。遞歸調(diào)用搜索過(guò)程,直到射線穿出八叉樹(shù)為止。最后求出射線與前述空間對(duì)象列表每個(gè)元素相交的t值,則最小t值對(duì)應(yīng)的空間對(duì)象即為所求結(jié)果。

    1.3 并行處理技術(shù)

    隨著計(jì)算機(jī)技術(shù)的發(fā)展,為了解決單核心CPU處理能力難以繼續(xù)提高的瓶頸,多核心CPU技術(shù)得到了飛速發(fā)展和廣泛應(yīng)用。因此模擬計(jì)算中應(yīng)充分利用多核CPU的計(jì)算加速能力。在確定計(jì)算硬件條件后,接下來(lái)的主要工作是計(jì)算程序的并行化編程。

    多核CPU屬于共享存儲(chǔ)的對(duì)稱多處理器[12]。OpenMP是為共享存儲(chǔ)環(huán)境編寫并行程序而設(shè)計(jì)的一個(gè)應(yīng)用編程接口,支持Fortran、C/C++語(yǔ)言。OpenMP標(biāo)準(zhǔn)中包括一套編譯指導(dǎo)語(yǔ)句和一個(gè)支持函數(shù)庫(kù)。在并行執(zhí)行程序時(shí),它采用 Fork/Join方式。其基本流程是:程序開(kāi)始時(shí)只有一個(gè)主線程,即程序中的串行部分都由主線程執(zhí)行,而并行部分則通過(guò)派生(Fork)其他線程來(lái)執(zhí)行。待并行執(zhí)行程序結(jié)束后(Join)才能執(zhí)行后面非并行部分代碼,如圖2所示。

    分析TPMC算法計(jì)算過(guò)程,可以在網(wǎng)格循環(huán)部分和試驗(yàn)粒子循環(huán)部分進(jìn)行并行。若在試驗(yàn)粒子循環(huán)部分采用并行,則需要進(jìn)行多次的線程派生注銷操作,效率低。因此,需要在網(wǎng)格循環(huán)級(jí)別進(jìn)行并行,即在C/C++語(yǔ)言中具體實(shí)現(xiàn)方法是需要在循環(huán)前加入“#pragma omp parallel for”編譯指導(dǎo)語(yǔ)句。

    圖2 共享存儲(chǔ)并行程序的并行執(zhí)行流程Fig.2 Flowchart of parallel program on shared memory machine

    2 計(jì)算結(jié)果

    2.1 圓盤返回流

    對(duì)文獻(xiàn)[10]中介紹的圓盤返回流進(jìn)行計(jì)算,顯示了八叉樹(shù)劃分算法的加速效果。算例的模型如圖3所示。計(jì)算結(jié)果為返回流比率,即與背景分子碰撞后返回放氣表面的分子數(shù)與總放氣分子數(shù)的比。算例的輸入?yún)?shù)見(jiàn)表1。

    圖3 模型示意圖Fig.3 Model sketch

    表1 算例輸入?yún)?shù)Table 1 Input parameters for test case

    對(duì)粗網(wǎng)格(90個(gè)三角形網(wǎng)格)和細(xì)網(wǎng)格(1128個(gè)三角形網(wǎng)格)2種情況進(jìn)行了比較計(jì)算。網(wǎng)格劃分如圖4所示。

    圖4 兩種網(wǎng)格劃分示意圖Fig.4 Comparison between coarse meshes and fine meshes

    計(jì)算機(jī)硬件配置為Intel Core2 Duo E7500雙核CPU、內(nèi)存4 GB、操作系統(tǒng)Windows XP。文獻(xiàn)[10]中的計(jì)算結(jié)果為 1.367 7×10-6,本文的計(jì)算結(jié)果為1.326 2×10-6,二者偏差為3%,符合度較好。加速效果如表2所示。

    表2 計(jì)算時(shí)間對(duì)比Table 2 Comparison of calculation time

    可以看出,在網(wǎng)格數(shù)較多時(shí),八叉樹(shù)的加速效果極為顯著。并且,加速系數(shù)可用ncell/(klog8ncell)近似表示,其中k為八叉樹(shù)子節(jié)點(diǎn)中空間單元的最大個(gè)數(shù),本文中設(shè)置為 8,ncell表示總的空間單元個(gè)數(shù)。

    并行加速比定義為S=T串/T并,即串行執(zhí)行時(shí)間與并行執(zhí)行時(shí)間的比值,當(dāng)并行加速比等于計(jì)算內(nèi)核數(shù)時(shí),稱為線性加速。本文中加速比接近計(jì)算內(nèi)核數(shù)2,說(shuō)明算法具有較好的加速特性。當(dāng)計(jì)算任務(wù)復(fù)雜,需要較短時(shí)間獲得結(jié)果時(shí),可通過(guò)簡(jiǎn)單增加線程數(shù)目實(shí)現(xiàn)。

    2.2 空間站污染仿真分析

    對(duì)特定空間站構(gòu)型下的污染分布狀況進(jìn)行了初步分析。空間站構(gòu)型取自文獻(xiàn)[13],共劃分24 166個(gè)三角形網(wǎng)格,網(wǎng)格劃分情況如圖5所示??梢钥闯?,在表面曲率變化大的部位進(jìn)行了自動(dòng)加密,使得網(wǎng)格能夠較精確地貼合表面。

    圖5 空間站網(wǎng)格劃分Fig.5 Meshing of a space station

    已有試驗(yàn)證明太陽(yáng)電池板放氣較為嚴(yán)重[14],因此污染源選擇其中一塊面積較大的太陽(yáng)電池板。軌道大氣來(lái)流速度方向?yàn)?z,其他參數(shù)與2.1節(jié)中算例一致。該計(jì)算采用雙線程,耗時(shí) 387.822 s。計(jì)算結(jié)果如圖6所示,顯示的是從污染源釋放的分子直接運(yùn)動(dòng)或經(jīng)過(guò)與大氣分子碰撞后運(yùn)動(dòng)至其他表面的情況,即單位面積表面網(wǎng)格上的入射分子數(shù)。結(jié)果表明,污染分布與污染源及敏感表面之間的視角因子成正比。因此,與污染源距離越近,污染越嚴(yán)重。圖6(a)中正對(duì)污染源的位置所受污染最為嚴(yán)重。圖6(b)中,背對(duì)污染源的表面,未受到污染影響。從圖6中還可以發(fā)現(xiàn)由于表面遮擋造成的污染分布不均勻。污染源表面以及與其處于同平面的另一個(gè)太陽(yáng)電池板也有分子沉積,這是由于污染分子與軌道大氣分子碰撞后返回而形成的。

    圖6 污染分布計(jì)算結(jié)果Fig.6 The calculated contamination distribution

    3 結(jié)束語(yǔ)

    航天器的放氣污染需要定量評(píng)估其影響。采用Monte-Carlo方法進(jìn)行模擬,用射線表征放氣分子的輸運(yùn)過(guò)程,可將污染輸運(yùn)過(guò)程轉(zhuǎn)化為射線追蹤問(wèn)題。由于航天器構(gòu)型復(fù)雜,射線追蹤的計(jì)算很耗時(shí)。本文將復(fù)雜構(gòu)型航天器表面用三角形非結(jié)構(gòu)網(wǎng)格表示,用八叉樹(shù)劃分算法開(kāi)展分子運(yùn)動(dòng)軌跡與表面相交的加速計(jì)算。同時(shí),利用基于共享內(nèi)存多核計(jì)算平臺(tái)的OpenMP軟件實(shí)現(xiàn)了計(jì)算的并行化。計(jì)算實(shí)例證明該算法可大大地提高計(jì)算速度,且與串行計(jì)算結(jié)果一致。此外,上述算法還可應(yīng)用于輻射傳熱角系數(shù)計(jì)算、空間粒子輻射影響分析等。

    [1]陳烈民.航天器結(jié)構(gòu)與機(jī)構(gòu)[M].北京: 中國(guó)科學(xué)技術(shù)出版社, 2005: 63-70

    [2]侯增祺, 胡金剛.航天器熱控制技術(shù): 原理及其應(yīng)用[M].北京: 中國(guó)科學(xué)技術(shù)出版社, 2007: 129-140

    [3]陳樹(shù)海.衛(wèi)星結(jié)構(gòu)中的非金屬材料[J].上海航天, 2004, 21(3): 39-44 CHEN S H.Nonmetal materials in satellite structure[J].Aerospace Shanghai, 2004, 21(3): 39-44

    [4]MILLER S K R, BANKS B.Degradation of spacecraft materials in the space environment[J].MRS Bulletin, 2010, 35(1): 20-24

    [5]SCHLAPPI B, ALTWEGG K, BALSIGER H, et al.Characterization of the gaseous spacecraft environment of Rosetta by ROSINA[C]//3rdAIAA Atmospheric Space Environments Conference.Honolulu, HI, United States: American Institute for Aeronautics and Astronautics, 2011: 1-8

    [6]HAEMMERLE V R, GERHARD J H.Cassini camera contamination anomaly: experiences and lessons learned: AIAA 2006-5834[R]

    [7]ROUSSEL J F, SOARES C E, SCHMIDL W D.Development and validation of a new return flux model for International Space Station[C]//Optical System Contamination: Effects, Measurements, and Control VII, Proceedings of SPIE, 2002, 4774: 172-183

    [8]焦子龍, 龐賀偉.航天器分子污染返回流計(jì)算方法[J].宇航學(xué)報(bào), 2009, 30(2): 415-421 JIAO Z L, PANG H W.Simulation methods for return flux of molecular contamination[J].Journal of Astronautics, 2009, 30(2): 415-421

    [9]焦子龍.航天器分子污染計(jì)算方法綜述[J].航天器環(huán)境工程, 2010, 27(2): 181-189 JIAO Z L.Review on simulation method for spacecraft molecular contamination[J].Spacecraft Environment Engineering, 2010, 27(2): 181-189

    [10]FAN C, GEE C, FONG M C.Monte Carlo simulation of molecular flux on simple spacecraft surfaces due to selfand ambient-scatter of outgassing molecules: AIAA 1993-2867[R]

    [11]AGATE M, GRIMSDALE R L, LISTER P F.The HERO algorithm for ray tracing octa-trees[M]//GRIMSDALER L, STRASSER W.Advances in computer graphics hardware IV.New York: Springer-Verlag, 1991: 23-27

    [12]周偉明.多核計(jì)算與程序設(shè)計(jì)[M].武漢: 華中科技大學(xué)出版社, 2009: 5-10

    [13]翟睿瓊, 姜海富, 田東波, 等.空間站原子氧環(huán)境仿真研究[J].裝備環(huán)境工程, 2014, 11(3): 35-39 ZHAI R Q, JIANG H F, TIAN D B.et al.Simulation of atomic oxygen interaction with space station[J].Equipment Environmental Engineering, 2014, 11(3): 35-39

    [14]焦子龍, 龐賀偉, 易忠, 等.航天器真空熱試驗(yàn)污染物成分分析[J].航天器環(huán)境工程, 2010, 27(6): 711-714 JIAO Z L, PANG H W, YI Z.et al.Analysis of contaminant components of spacecraft thermal vacuum test[J].Spacecraft Environment Engineering, 2010, 27(6): 711-714

    (編輯:肖福根)

    Parallel and octa-tree-based acceleration of outgassing contamination simulation for spacecraft of complex geometries

    JIAO Zilong, JIANG Lixiang, HUANG Jianguo, SUN Jipeng, ZHU Yunfei
    (Science and Technology on Reliability and Environmental Engineering Laboratory, Beijing Institute of Spacecraft Environment Engineering, Beijing 100094, China)

    The organic materials used in spacecraft may outgas under the space environment, including, mainly, the vacuum, the high & low temperature, and the influence of the ionization irradiation.The outgassing product has harmful effects on the spacecraft itself.So a quantitative assessment of the contamination effects is necessary.The molecular transport process is modeled by the Monte-Carlo method, with the rays representing the molecule trajectories, to simulate the motion process of a large number of molecules.Thus the transport process is transformed into a ray tracing problem.Due to the complexity of the spacecraft configurations, the ray tracing computation might be very time consuming, therefore, it is extremely necessary to have an accelerating algorithm.In this paper, the complex geometry of the spacecraft surface is modeled with triangular unstructured grids, the octa-tree space subdivision technology is used to accelerate the molecular motion and surface intersection calculations.At the same time, based on the shared memory multi-core computing platform, the OpenMP API is used to implement parallel computing for the simulation.The test case shows a significant acceleration in the calculation, and the developed method is applied to the contamination analysis of a space station.

    spacecraft; outgassing contamination; simulation computation; complex geometries; octa-tree; parallel computing

    X502

    :A

    :1673-1379(2016)06-0601-05

    10.3969/j.issn.1673-1379.2016.06.005

    焦子龍(1979—),男,博士學(xué)位,高級(jí)工程師,從事航天器特殊空間環(huán)境效應(yīng)分析、試驗(yàn)及防護(hù)技術(shù)研究。E-mail: novanova2009@163.com。

    2016-04-01;

    :2016-11-10

    猜你喜歡
    八叉樹(shù)構(gòu)型航天器
    三維十字鏈表八叉樹(shù)的高效檢索實(shí)現(xiàn)
    2022 年第二季度航天器發(fā)射統(tǒng)計(jì)
    分子和離子立體構(gòu)型的判定
    2019 年第二季度航天器發(fā)射統(tǒng)計(jì)
    2018 年第三季度航天器發(fā)射統(tǒng)計(jì)
    2018年第二季度航天器發(fā)射統(tǒng)計(jì)
    航天器受迫繞飛構(gòu)型設(shè)計(jì)與控制
    遙感衛(wèi)星平臺(tái)與載荷一體化構(gòu)型
    兩個(gè)具stp三維拓?fù)錁?gòu)型的稀土配位聚合物{[Ln2(pda)3(H2O)2]·2H2O}n(Ln=Nd,La)
    散亂點(diǎn)云線性八叉樹(shù)結(jié)構(gòu)在GPU中的實(shí)現(xiàn)
    亚洲成av人片在线播放无| 在线观看日韩欧美| 欧美国产日韩亚洲一区| 曰老女人黄片| 午夜福利欧美成人| av在线天堂中文字幕| 久久久久国内视频| 国产精品永久免费网站| 白带黄色成豆腐渣| 美女午夜性视频免费| 欧美一区二区国产精品久久精品 | 久久久国产成人精品二区| 高潮久久久久久久久久久不卡| 不卡av一区二区三区| 国产一级毛片七仙女欲春2| 国产又色又爽无遮挡免费看| 嫩草影院精品99| a在线观看视频网站| 免费在线观看成人毛片| 黑人欧美特级aaaaaa片| 在线观看免费视频日本深夜| 中文字幕人成人乱码亚洲影| 亚洲性夜色夜夜综合| 18美女黄网站色大片免费观看| 亚洲免费av在线视频| 精品国内亚洲2022精品成人| 日韩中文字幕欧美一区二区| 国产在线观看jvid| 男女之事视频高清在线观看| 热99re8久久精品国产| 91大片在线观看| 亚洲最大成人中文| 97人妻精品一区二区三区麻豆| 99在线人妻在线中文字幕| 亚洲五月天丁香| 国产单亲对白刺激| 在线十欧美十亚洲十日本专区| 伊人久久大香线蕉亚洲五| 黄片小视频在线播放| 成人午夜高清在线视频| 久9热在线精品视频| 国产日本99.免费观看| 怎么达到女性高潮| 999久久久国产精品视频| 免费电影在线观看免费观看| 深夜精品福利| 成人三级黄色视频| 国产精品久久电影中文字幕| 国产片内射在线| 久久 成人 亚洲| 欧美不卡视频在线免费观看 | 久久婷婷人人爽人人干人人爱| 悠悠久久av| 女同久久另类99精品国产91| 国产激情偷乱视频一区二区| 亚洲中文日韩欧美视频| 久久久久免费精品人妻一区二区| 手机成人av网站| 亚洲 欧美 日韩 在线 免费| 久久精品人妻少妇| 97人妻精品一区二区三区麻豆| 国产欧美日韩一区二区三| 久久久水蜜桃国产精品网| 亚洲九九香蕉| 大型av网站在线播放| av视频在线观看入口| 精品久久久久久久毛片微露脸| 成人高潮视频无遮挡免费网站| 欧美乱码精品一区二区三区| netflix在线观看网站| 亚洲色图av天堂| 老汉色∧v一级毛片| 午夜福利18| 宅男免费午夜| 久久久久久国产a免费观看| 国产高清有码在线观看视频 | 91大片在线观看| 好看av亚洲va欧美ⅴa在| 全区人妻精品视频| 男女午夜视频在线观看| 99热只有精品国产| 国产精品爽爽va在线观看网站| 亚洲av电影在线进入| 变态另类丝袜制服| 黄色成人免费大全| 欧美三级亚洲精品| 熟妇人妻久久中文字幕3abv| 欧美成人一区二区免费高清观看 | 国产精品,欧美在线| 欧美性猛交黑人性爽| 欧美黑人巨大hd| 免费在线观看黄色视频的| 精品欧美一区二区三区在线| 日本一本二区三区精品| 国产亚洲精品综合一区在线观看 | 欧美zozozo另类| 禁无遮挡网站| www国产在线视频色| 久久久精品大字幕| 国产亚洲精品久久久久久毛片| 亚洲片人在线观看| 久久久久久久精品吃奶| 久久午夜综合久久蜜桃| 中文亚洲av片在线观看爽| 欧美黄色片欧美黄色片| 91国产中文字幕| 久久国产精品人妻蜜桃| 此物有八面人人有两片| 久久精品国产亚洲av高清一级| 午夜福利成人在线免费观看| 曰老女人黄片| 国产成年人精品一区二区| 久久午夜亚洲精品久久| 亚洲av片天天在线观看| 国产一区二区在线观看日韩 | 国产伦在线观看视频一区| 亚洲aⅴ乱码一区二区在线播放 | 中文字幕久久专区| 777久久人妻少妇嫩草av网站| 每晚都被弄得嗷嗷叫到高潮| 亚洲专区中文字幕在线| 国产精品影院久久| 黄片小视频在线播放| 欧美在线一区亚洲| 精品一区二区三区四区五区乱码| 黄色丝袜av网址大全| 一本综合久久免费| 精品国产超薄肉色丝袜足j| 国产在线观看jvid| 亚洲国产欧美人成| www.www免费av| 首页视频小说图片口味搜索| 天堂av国产一区二区熟女人妻 | 午夜免费观看网址| 日本a在线网址| 精品午夜福利视频在线观看一区| 成年版毛片免费区| 久久天堂一区二区三区四区| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人精品久久二区二区免费| 18禁黄网站禁片免费观看直播| 老鸭窝网址在线观看| 国产高清videossex| 成人特级黄色片久久久久久久| 波多野结衣高清作品| 国产精品美女特级片免费视频播放器 | 老司机午夜十八禁免费视频| 欧美 亚洲 国产 日韩一| 一a级毛片在线观看| 在线观看午夜福利视频| 天堂av国产一区二区熟女人妻 | 亚洲人与动物交配视频| 欧美精品啪啪一区二区三区| 丰满人妻一区二区三区视频av | 国产一区二区激情短视频| 很黄的视频免费| 成人高潮视频无遮挡免费网站| 国产精品一及| 又爽又黄无遮挡网站| 黄色成人免费大全| 听说在线观看完整版免费高清| 精品日产1卡2卡| 亚洲熟女毛片儿| 久久久国产精品麻豆| 欧美一级毛片孕妇| 99久久久亚洲精品蜜臀av| 中文字幕熟女人妻在线| 制服人妻中文乱码| 亚洲 国产 在线| 日韩av在线大香蕉| 精品高清国产在线一区| 啦啦啦免费观看视频1| 高清毛片免费观看视频网站| 韩国av一区二区三区四区| 久久精品夜夜夜夜夜久久蜜豆 | 精品不卡国产一区二区三区| 90打野战视频偷拍视频| 91大片在线观看| 久久久久亚洲av毛片大全| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品中文字幕一二三四区| 天堂影院成人在线观看| 久久人妻福利社区极品人妻图片| 日本三级黄在线观看| 999精品在线视频| 听说在线观看完整版免费高清| 熟女少妇亚洲综合色aaa.| 妹子高潮喷水视频| 国产成人精品久久二区二区91| 男女那种视频在线观看| 少妇人妻一区二区三区视频| 久久久久九九精品影院| 中文资源天堂在线| 久久久水蜜桃国产精品网| www日本黄色视频网| 午夜免费激情av| 看黄色毛片网站| 成人三级黄色视频| 午夜福利视频1000在线观看| 亚洲一区二区三区不卡视频| 国产蜜桃级精品一区二区三区| 麻豆国产av国片精品| 欧美乱色亚洲激情| 中文字幕人妻丝袜一区二区| 久99久视频精品免费| 天堂av国产一区二区熟女人妻 | 久9热在线精品视频| 久久久久性生活片| 淫妇啪啪啪对白视频| 超碰成人久久| 97碰自拍视频| 亚洲av美国av| cao死你这个sao货| 欧美中文日本在线观看视频| 国产亚洲精品综合一区在线观看 | 久9热在线精品视频| 亚洲在线自拍视频| 亚洲电影在线观看av| 国产免费av片在线观看野外av| 亚洲av熟女| 可以免费在线观看a视频的电影网站| 人人妻人人看人人澡| 亚洲 国产 在线| 亚洲av日韩精品久久久久久密| 亚洲一卡2卡3卡4卡5卡精品中文| 99精品在免费线老司机午夜| 黄色 视频免费看| 宅男免费午夜| 成人av在线播放网站| 欧美乱色亚洲激情| АⅤ资源中文在线天堂| 在线看三级毛片| 成人国产综合亚洲| 五月玫瑰六月丁香| 一卡2卡三卡四卡精品乱码亚洲| 欧美三级亚洲精品| 免费人成视频x8x8入口观看| 日日摸夜夜添夜夜添小说| 成人三级黄色视频| 久久性视频一级片| 欧美绝顶高潮抽搐喷水| 男女床上黄色一级片免费看| 精品日产1卡2卡| 我的老师免费观看完整版| 亚洲欧美一区二区三区黑人| 亚洲 欧美一区二区三区| 啦啦啦免费观看视频1| 国产av又大| 欧美最黄视频在线播放免费| 手机成人av网站| 欧美+亚洲+日韩+国产| 国产一区二区在线av高清观看| 两性夫妻黄色片| 精品第一国产精品| 妹子高潮喷水视频| www.999成人在线观看| 操出白浆在线播放| 亚洲午夜理论影院| 村上凉子中文字幕在线| 亚洲性夜色夜夜综合| 成人国语在线视频| 亚洲最大成人中文| 国产精品久久久人人做人人爽| 久久性视频一级片| 成人精品一区二区免费| 天天躁夜夜躁狠狠躁躁| 中文字幕熟女人妻在线| 99久久综合精品五月天人人| 18禁观看日本| 18美女黄网站色大片免费观看| 亚洲成av人片在线播放无| 1024香蕉在线观看| 99热这里只有精品一区 | 亚洲自拍偷在线| 久久九九热精品免费| 亚洲自偷自拍图片 自拍| 亚洲熟妇熟女久久| www.熟女人妻精品国产| 亚洲精品中文字幕在线视频| 免费高清视频大片| 黄频高清免费视频| 身体一侧抽搐| 一个人免费在线观看的高清视频| 两人在一起打扑克的视频| 99久久99久久久精品蜜桃| 亚洲成人中文字幕在线播放| 国产成人精品久久二区二区91| 伦理电影免费视频| 在线观看66精品国产| 美女午夜性视频免费| 国产精品精品国产色婷婷| 中出人妻视频一区二区| 国产伦在线观看视频一区| 男女下面进入的视频免费午夜| 国产精品久久久久久人妻精品电影| 成人国语在线视频| 99精品在免费线老司机午夜| av国产免费在线观看| 欧美成狂野欧美在线观看| 可以在线观看毛片的网站| 免费看a级黄色片| 日本免费a在线| 成年女人毛片免费观看观看9| 一区二区三区激情视频| 欧美日韩中文字幕国产精品一区二区三区| 成人三级做爰电影| 国产免费男女视频| 精品国产美女av久久久久小说| 俺也久久电影网| 亚洲男人的天堂狠狠| 两个人的视频大全免费| 久久久国产精品麻豆| 久久久久久亚洲精品国产蜜桃av| www日本在线高清视频| 日韩欧美一区二区三区在线观看| 俺也久久电影网| 亚洲av成人一区二区三| 亚洲,欧美精品.| 欧美色视频一区免费| 久久 成人 亚洲| netflix在线观看网站| 神马国产精品三级电影在线观看 | 老司机午夜福利在线观看视频| 国产精品免费视频内射| 黄色视频不卡| 欧美3d第一页| 成年版毛片免费区| 男人舔女人下体高潮全视频| 麻豆国产av国片精品| 久久精品国产综合久久久| 欧美在线黄色| 男插女下体视频免费在线播放| 免费一级毛片在线播放高清视频| 村上凉子中文字幕在线| 国产欧美日韩一区二区精品| 淫妇啪啪啪对白视频| 欧美一区二区精品小视频在线| 中文字幕熟女人妻在线| 亚洲免费av在线视频| 亚洲va日本ⅴa欧美va伊人久久| 露出奶头的视频| www.熟女人妻精品国产| 国产精品综合久久久久久久免费| 欧美成人性av电影在线观看| 精品国产乱子伦一区二区三区| 日本五十路高清| a级毛片在线看网站| 国产精品野战在线观看| 国产不卡一卡二| 狂野欧美白嫩少妇大欣赏| 国产精品国产高清国产av| 极品教师在线免费播放| 正在播放国产对白刺激| 精品福利观看| 亚洲av成人av| 国产成人欧美在线观看| 我的老师免费观看完整版| 国产精品久久久久久人妻精品电影| 在线观看舔阴道视频| 国产高清激情床上av| 首页视频小说图片口味搜索| 日本黄色视频三级网站网址| 少妇裸体淫交视频免费看高清 | 亚洲中文字幕日韩| 亚洲国产欧美人成| 亚洲欧美精品综合久久99| 黄色成人免费大全| 色综合站精品国产| 欧美国产日韩亚洲一区| 人人妻,人人澡人人爽秒播| 丝袜人妻中文字幕| 香蕉久久夜色| 亚洲精品色激情综合| 少妇的丰满在线观看| 亚洲精品在线美女| 久久精品国产亚洲av高清一级| xxx96com| 成人国产一区最新在线观看| 国产在线观看jvid| 18禁裸乳无遮挡免费网站照片| 中文在线观看免费www的网站 | 亚洲国产精品成人综合色| www.熟女人妻精品国产| 日本黄色视频三级网站网址| 免费观看人在逋| 又黄又粗又硬又大视频| 亚洲av电影在线进入| 亚洲av成人精品一区久久| 亚洲,欧美精品.| 校园春色视频在线观看| av天堂在线播放| 国产日本99.免费观看| 在线观看www视频免费| 床上黄色一级片| 久久人妻av系列| 日韩精品免费视频一区二区三区| 日韩欧美在线二视频| 欧美av亚洲av综合av国产av| 欧美人与性动交α欧美精品济南到| 日韩欧美 国产精品| 欧美成人一区二区免费高清观看 | netflix在线观看网站| 中文字幕熟女人妻在线| 中文字幕av在线有码专区| 欧美性猛交黑人性爽| 中文字幕久久专区| 日韩成人在线观看一区二区三区| 99久久综合精品五月天人人| 亚洲国产精品999在线| 亚洲人成电影免费在线| 久久人人精品亚洲av| 亚洲成人久久爱视频| 波多野结衣高清无吗| 亚洲欧美日韩高清专用| 麻豆成人av在线观看| 国产精华一区二区三区| 亚洲中文字幕日韩| 日韩av在线大香蕉| 日日爽夜夜爽网站| 国产精品香港三级国产av潘金莲| 国产精品久久久久久人妻精品电影| 国产激情欧美一区二区| 亚洲国产精品sss在线观看| 成人午夜高清在线视频| 国产1区2区3区精品| 亚洲美女视频黄频| 午夜精品在线福利| 99久久99久久久精品蜜桃| 国产精品av视频在线免费观看| 美女 人体艺术 gogo| 国产一区二区三区视频了| 90打野战视频偷拍视频| 亚洲国产精品久久男人天堂| 悠悠久久av| 亚洲七黄色美女视频| 日日爽夜夜爽网站| av在线播放免费不卡| 午夜福利视频1000在线观看| 欧美高清成人免费视频www| 欧美三级亚洲精品| 日日摸夜夜添夜夜添小说| 精品久久久久久,| 国产麻豆成人av免费视频| 白带黄色成豆腐渣| 亚洲男人天堂网一区| 成人三级做爰电影| 老熟妇仑乱视频hdxx| 这个男人来自地球电影免费观看| 在线观看舔阴道视频| 美女黄网站色视频| 国内精品久久久久久久电影| 动漫黄色视频在线观看| 亚洲七黄色美女视频| 日日摸夜夜添夜夜添小说| 亚洲美女视频黄频| 亚洲人成伊人成综合网2020| 久久九九热精品免费| 天天一区二区日本电影三级| 九九热线精品视视频播放| 天堂影院成人在线观看| 国产午夜福利久久久久久| cao死你这个sao货| 亚洲欧美日韩东京热| 少妇裸体淫交视频免费看高清 | 老汉色∧v一级毛片| bbb黄色大片| 男人舔女人的私密视频| 欧美中文日本在线观看视频| 搡老岳熟女国产| 一本精品99久久精品77| 精品乱码久久久久久99久播| 观看免费一级毛片| 久久久久久久午夜电影| 9191精品国产免费久久| 一个人免费在线观看的高清视频| 日韩欧美在线乱码| 极品教师在线免费播放| 国语自产精品视频在线第100页| av超薄肉色丝袜交足视频| 色综合婷婷激情| 亚洲av熟女| 成人国产综合亚洲| 一本一本综合久久| 三级国产精品欧美在线观看 | 精品日产1卡2卡| 搡老岳熟女国产| 91大片在线观看| 18禁国产床啪视频网站| 99国产精品99久久久久| 久久久精品欧美日韩精品| 中文资源天堂在线| 夜夜爽天天搞| 亚洲中文字幕日韩| 日本撒尿小便嘘嘘汇集6| av中文乱码字幕在线| 欧美国产日韩亚洲一区| 国产精品一及| 在线观看舔阴道视频| 少妇裸体淫交视频免费看高清 | 国产一区二区在线av高清观看| 国产99久久九九免费精品| 国产亚洲精品久久久久5区| 亚洲自偷自拍图片 自拍| 搡老妇女老女人老熟妇| 精品国产亚洲在线| 亚洲一卡2卡3卡4卡5卡精品中文| 丁香欧美五月| 国产精品 欧美亚洲| 国产av不卡久久| 黄色毛片三级朝国网站| 免费搜索国产男女视频| 在线观看免费视频日本深夜| 成年人黄色毛片网站| 欧美日本亚洲视频在线播放| 精品久久蜜臀av无| 五月玫瑰六月丁香| 丁香六月欧美| 人成视频在线观看免费观看| 中文亚洲av片在线观看爽| 一本一本综合久久| 999久久久国产精品视频| 欧美日韩一级在线毛片| 99re在线观看精品视频| 久久伊人香网站| 色综合站精品国产| 又黄又粗又硬又大视频| 给我免费播放毛片高清在线观看| 老司机靠b影院| 老司机福利观看| 国产欧美日韩一区二区三| 国产高清videossex| 国产一区二区三区在线臀色熟女| 看黄色毛片网站| 色精品久久人妻99蜜桃| 久久精品91蜜桃| 亚洲成a人片在线一区二区| 18禁黄网站禁片免费观看直播| 亚洲欧洲精品一区二区精品久久久| 亚洲精品美女久久久久99蜜臀| 可以在线观看毛片的网站| 人成视频在线观看免费观看| 国产亚洲av高清不卡| 国产区一区二久久| 国产真人三级小视频在线观看| 精品国内亚洲2022精品成人| 中文字幕熟女人妻在线| 日韩中文字幕欧美一区二区| 久9热在线精品视频| 亚洲国产看品久久| 久久午夜综合久久蜜桃| 岛国视频午夜一区免费看| 国产免费男女视频| 亚洲精华国产精华精| 97碰自拍视频| 99久久综合精品五月天人人| 亚洲精品粉嫩美女一区| 欧美乱码精品一区二区三区| 十八禁人妻一区二区| 悠悠久久av| 亚洲国产中文字幕在线视频| 18禁国产床啪视频网站| 老司机深夜福利视频在线观看| 亚洲黑人精品在线| 97碰自拍视频| 色综合欧美亚洲国产小说| 国产精品九九99| 国产午夜精品论理片| tocl精华| 日韩精品青青久久久久久| 午夜久久久久精精品| 男女床上黄色一级片免费看| 中文资源天堂在线| 色精品久久人妻99蜜桃| 少妇粗大呻吟视频| 欧美国产日韩亚洲一区| 老汉色av国产亚洲站长工具| 精品福利观看| 国产熟女午夜一区二区三区| 男人舔女人下体高潮全视频| 精品久久久久久久久久免费视频| 亚洲一码二码三码区别大吗| 久久久久久免费高清国产稀缺| 免费看日本二区| 在线永久观看黄色视频| 韩国av一区二区三区四区| 男女视频在线观看网站免费 | 久久香蕉激情| 村上凉子中文字幕在线| 欧美最黄视频在线播放免费| 叶爱在线成人免费视频播放| 久久久久久免费高清国产稀缺| 中文字幕精品亚洲无线码一区| 黄片小视频在线播放| 99久久99久久久精品蜜桃| 欧美乱妇无乱码| 神马国产精品三级电影在线观看 | 九色成人免费人妻av| 久久精品国产综合久久久| 黄片大片在线免费观看| 97碰自拍视频| 亚洲av电影不卡..在线观看| 五月伊人婷婷丁香| 看片在线看免费视频| 少妇裸体淫交视频免费看高清 | 一个人免费在线观看电影 | 妹子高潮喷水视频| 色av中文字幕| 国产欧美日韩一区二区三| 天天一区二区日本电影三级| 又大又爽又粗| 中文字幕人成人乱码亚洲影| 亚洲人成网站高清观看| 老鸭窝网址在线观看| 三级国产精品欧美在线观看 | 亚洲欧美日韩高清在线视频|