• <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)
    男插女下体视频免费在线播放| 国产69精品久久久久777片| 久99久视频精品免费| 久久精品国产亚洲网站| 午夜精品国产一区二区电影 | 精品久久国产蜜桃| 边亲边吃奶的免费视频| 国产毛片a区久久久久| 国产亚洲精品av在线| 日韩一区二区视频免费看| 内地一区二区视频在线| 特大巨黑吊av在线直播| 久久久久网色| videossex国产| 少妇裸体淫交视频免费看高清| 中国国产av一级| 69人妻影院| 久久久久久久久久黄片| 边亲边吃奶的免费视频| 国产爱豆传媒在线观看| 精品久久久久久成人av| 精品人妻一区二区三区麻豆| 夜夜爽夜夜爽视频| 亚洲av二区三区四区| 久久精品人妻少妇| av又黄又爽大尺度在线免费看 | 国产色爽女视频免费观看| 亚洲成人av在线免费| 久久久久精品久久久久真实原创| 卡戴珊不雅视频在线播放| 国产单亲对白刺激| 国产亚洲精品av在线| 久久久午夜欧美精品| 国产精品日韩av在线免费观看| 国产在视频线精品| 99久久人妻综合| 少妇高潮的动态图| 精品一区二区免费观看| 亚洲综合精品二区| 青春草视频在线免费观看| 国产午夜福利久久久久久| av免费在线看不卡| 久久久久九九精品影院| 国产免费又黄又爽又色| 欧美又色又爽又黄视频| 少妇裸体淫交视频免费看高清| 久久精品久久久久久久性| 欧美变态另类bdsm刘玥| 亚洲乱码一区二区免费版| 国产亚洲最大av| 日韩 亚洲 欧美在线| 国产高清三级在线| 欧美三级亚洲精品| 久久久精品94久久精品| 男人狂女人下面高潮的视频| 国产伦一二天堂av在线观看| 草草在线视频免费看| 国产免费视频播放在线视频 | 日韩国内少妇激情av| 中文欧美无线码| av在线天堂中文字幕| 99热全是精品| 亚洲最大成人中文| 国产精品久久久久久精品电影| 久久韩国三级中文字幕| 丝袜喷水一区| 国内少妇人妻偷人精品xxx网站| 欧美一区二区亚洲| 国产91av在线免费观看| 在线天堂最新版资源| 丰满人妻一区二区三区视频av| 一区二区三区乱码不卡18| 国产大屁股一区二区在线视频| 人人妻人人澡人人爽人人夜夜 | 少妇的逼好多水| 亚洲欧美日韩高清专用| 啦啦啦观看免费观看视频高清| 麻豆久久精品国产亚洲av| 欧美精品国产亚洲| 简卡轻食公司| 内射极品少妇av片p| 欧美日韩精品成人综合77777| 久久6这里有精品| 又爽又黄无遮挡网站| 免费看光身美女| 婷婷色综合大香蕉| 国产精品久久视频播放| 国产伦精品一区二区三区四那| 久久欧美精品欧美久久欧美| 丝袜美腿在线中文| 国产成人精品一,二区| 免费观看在线日韩| 嫩草影院精品99| 日本黄色片子视频| 黄色配什么色好看| 日本与韩国留学比较| 日本av手机在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 黄色一级大片看看| 亚洲精品乱久久久久久| 久久久久性生活片| 内射极品少妇av片p| 少妇丰满av| 亚洲婷婷狠狠爱综合网| 日韩一区二区视频免费看| 亚洲精品久久久久久婷婷小说 | 国产一区有黄有色的免费视频 | 干丝袜人妻中文字幕| 51国产日韩欧美| 精品久久久久久成人av| 亚洲av成人精品一区久久| 亚洲国产日韩欧美精品在线观看| 一个人看的www免费观看视频| 三级国产精品欧美在线观看| 黄片无遮挡物在线观看| 丰满人妻一区二区三区视频av| 成人漫画全彩无遮挡| 69人妻影院| 国产老妇伦熟女老妇高清| 欧美极品一区二区三区四区| 国产成人精品一,二区| 在线播放无遮挡| 熟妇人妻久久中文字幕3abv| 22中文网久久字幕| 久久久成人免费电影| 变态另类丝袜制服| 国产黄色小视频在线观看| 1024手机看黄色片| 国产69精品久久久久777片| 十八禁国产超污无遮挡网站| 人人妻人人看人人澡| 老司机影院毛片| 久久欧美精品欧美久久欧美| 国产一区亚洲一区在线观看| 一个人免费在线观看电影| 国产精品熟女久久久久浪| 免费看日本二区| 少妇的逼水好多| 免费观看在线日韩| 日韩欧美三级三区| 亚洲欧美清纯卡通| 亚洲人成网站在线播| 最近中文字幕2019免费版| 九色成人免费人妻av| 国产视频内射| 亚洲中文字幕日韩| 观看美女的网站| 内地一区二区视频在线| 亚洲精品自拍成人| 一边摸一边抽搐一进一小说| 中文资源天堂在线| av免费在线看不卡| 欧美3d第一页| 在线观看66精品国产| 一边亲一边摸免费视频| 亚洲av不卡在线观看| 99视频精品全部免费 在线| 国产国拍精品亚洲av在线观看| 麻豆乱淫一区二区| av天堂中文字幕网| 国产视频首页在线观看| 老师上课跳d突然被开到最大视频| 欧美三级亚洲精品| 亚洲av福利一区| av在线播放精品| 精品久久久久久久久久久久久| 午夜免费男女啪啪视频观看| 小说图片视频综合网站| 国产精品国产三级国产av玫瑰| 丰满少妇做爰视频| 国产免费又黄又爽又色| 天堂中文最新版在线下载 | 插阴视频在线观看视频| 国产av一区在线观看免费| 国产高清不卡午夜福利| 91av网一区二区| 国产精品精品国产色婷婷| 狂野欧美白嫩少妇大欣赏| 老师上课跳d突然被开到最大视频| 成年版毛片免费区| 亚洲真实伦在线观看| a级毛片免费高清观看在线播放| 禁无遮挡网站| 中文字幕制服av| 久久久久久久久中文| 精品一区二区三区人妻视频| 亚洲av成人精品一二三区| 我要搜黄色片| 欧美xxxx性猛交bbbb| 国产成人精品一,二区| 亚洲成色77777| 国产精品野战在线观看| 在线天堂最新版资源| 成人国产麻豆网| 色播亚洲综合网| 婷婷色av中文字幕| 欧美人与善性xxx| 97超视频在线观看视频| 成人鲁丝片一二三区免费| 精品国产露脸久久av麻豆 | 69人妻影院| 亚洲在线自拍视频| 精品人妻视频免费看| 一本久久精品| 日韩,欧美,国产一区二区三区 | 91狼人影院| 九九久久精品国产亚洲av麻豆| 亚洲国产精品成人久久小说| 99视频精品全部免费 在线| 中国国产av一级| 国产在视频线在精品| 亚洲精品亚洲一区二区| 国产亚洲av嫩草精品影院| 中国美白少妇内射xxxbb| 七月丁香在线播放| 三级国产精品片| 欧美区成人在线视频| 国产精品爽爽va在线观看网站| 在线播放国产精品三级| 亚洲精品,欧美精品| 97热精品久久久久久| kizo精华| 国产亚洲午夜精品一区二区久久 | 久久精品久久久久久久性| 亚洲国产高清在线一区二区三| 免费观看人在逋| 亚洲精品一区蜜桃| 全区人妻精品视频| 久久婷婷人人爽人人干人人爱| 国模一区二区三区四区视频| 国产亚洲精品av在线| 97超视频在线观看视频| 亚洲精品乱码久久久v下载方式| 日韩精品青青久久久久久| 日韩欧美精品免费久久| 欧美不卡视频在线免费观看| 97在线视频观看| 国产在线一区二区三区精 | 国产片特级美女逼逼视频| 国产精品一区二区三区四区免费观看| 免费看美女性在线毛片视频| 精品久久久久久久人妻蜜臀av| 色综合亚洲欧美另类图片| 免费在线观看成人毛片| 天天躁夜夜躁狠狠久久av| 18禁在线无遮挡免费观看视频| 国产亚洲av片在线观看秒播厂 | 亚洲精品aⅴ在线观看| 自拍偷自拍亚洲精品老妇| 亚洲av电影在线观看一区二区三区 | 成人鲁丝片一二三区免费| 国产探花极品一区二区| 97人妻精品一区二区三区麻豆| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产成人一精品久久久| 国产精品1区2区在线观看.| 美女黄网站色视频| 国产老妇伦熟女老妇高清| 亚洲人成网站在线观看播放| 中文欧美无线码| 91精品伊人久久大香线蕉| 啦啦啦观看免费观看视频高清| 亚洲乱码一区二区免费版| 少妇猛男粗大的猛烈进出视频 | 青春草亚洲视频在线观看| 特大巨黑吊av在线直播| 99视频精品全部免费 在线| 亚洲国产欧美人成| 成人亚洲精品av一区二区| 在线天堂最新版资源| 麻豆国产97在线/欧美| 青春草国产在线视频| kizo精华| av免费观看日本| av国产久精品久网站免费入址| 欧美一级a爱片免费观看看| 久久精品人妻少妇| 日日干狠狠操夜夜爽| 亚洲久久久久久中文字幕| 国产精品久久久久久久电影| 麻豆国产97在线/欧美| 青春草国产在线视频| 看十八女毛片水多多多| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩,欧美,国产一区二区三区 | 精品人妻一区二区三区麻豆| av国产免费在线观看| 国产中年淑女户外野战色| 久久精品91蜜桃| 91在线精品国自产拍蜜月| 国产精品人妻久久久久久| av又黄又爽大尺度在线免费看 | 日韩欧美精品v在线| 国产精品精品国产色婷婷| 亚洲国产精品国产精品| 夜夜看夜夜爽夜夜摸| 国产精品三级大全| 如何舔出高潮| 亚洲电影在线观看av| 欧美性猛交黑人性爽| 久久人人爽人人片av| 不卡视频在线观看欧美| 久久婷婷人人爽人人干人人爱| 午夜精品在线福利| 国产精品人妻久久久久久| 日韩欧美精品v在线| 日韩欧美精品免费久久| 成人一区二区视频在线观看| 91aial.com中文字幕在线观看| videossex国产| 秋霞伦理黄片| 91aial.com中文字幕在线观看| 女人十人毛片免费观看3o分钟| 欧美高清性xxxxhd video| 九九在线视频观看精品| 国产亚洲91精品色在线| 欧美一级a爱片免费观看看| 亚洲av男天堂| 最近最新中文字幕大全电影3| 久久婷婷人人爽人人干人人爱| 老司机福利观看| 国产精品无大码| 三级毛片av免费| av在线蜜桃| 97热精品久久久久久| 少妇高潮的动态图| 国产免费视频播放在线视频 | 亚洲精品一区蜜桃| 精品久久久久久电影网 | 高清av免费在线| 亚洲在线自拍视频| 精品午夜福利在线看| 中文字幕人妻熟人妻熟丝袜美| eeuss影院久久| 美女高潮的动态| 午夜福利在线在线| 久久精品久久久久久久性| 午夜福利在线在线| 日本黄色片子视频| 国产一区二区三区av在线| 青春草国产在线视频| 一区二区三区高清视频在线| 99久国产av精品国产电影| 最近中文字幕2019免费版| 亚洲五月天丁香| 99久久九九国产精品国产免费| 免费观看性生交大片5| 国产精品99久久久久久久久| 一个人免费在线观看电影| 欧美日韩国产亚洲二区| 免费看美女性在线毛片视频| 国产探花极品一区二区| 日日摸夜夜添夜夜爱| 欧美激情久久久久久爽电影| 男女视频在线观看网站免费| .国产精品久久| 色播亚洲综合网| 国产毛片a区久久久久| av视频在线观看入口| 亚洲国产精品合色在线| 嫩草影院入口| 亚洲中文字幕日韩| 欧美色视频一区免费| 国产免费福利视频在线观看| 午夜爱爱视频在线播放| 日韩av在线免费看完整版不卡| 麻豆成人午夜福利视频| 久久久亚洲精品成人影院| 国产精品美女特级片免费视频播放器| 国产乱人视频| 国产乱来视频区| 日本与韩国留学比较| 欧美丝袜亚洲另类| 久久国产乱子免费精品| 内地一区二区视频在线| 高清日韩中文字幕在线| 中文字幕久久专区| 色综合站精品国产| 国产伦理片在线播放av一区| 观看美女的网站| 亚洲欧美日韩东京热| 搡老妇女老女人老熟妇| 国产精品.久久久| 一个人观看的视频www高清免费观看| 日韩,欧美,国产一区二区三区 | 永久网站在线| 简卡轻食公司| 亚洲真实伦在线观看| 国产69精品久久久久777片| 日韩欧美精品v在线| 嘟嘟电影网在线观看| 少妇人妻一区二区三区视频| 18禁动态无遮挡网站| 国产一区二区在线观看日韩| 亚洲国产精品国产精品| 亚洲国产欧美在线一区| 赤兔流量卡办理| 一级毛片电影观看 | 国内少妇人妻偷人精品xxx网站| 99久久无色码亚洲精品果冻| 99久久成人亚洲精品观看| 国产午夜福利久久久久久| 国产精品熟女久久久久浪| 国产精品一区二区性色av| 成人性生交大片免费视频hd| 久久久午夜欧美精品| 美女高潮的动态| 精华霜和精华液先用哪个| 亚洲欧美精品综合久久99| 热99在线观看视频| 亚洲av成人精品一区久久| 高清av免费在线| 赤兔流量卡办理| 一级二级三级毛片免费看| 久久精品国产亚洲av涩爱| 日日干狠狠操夜夜爽| 国语自产精品视频在线第100页| 亚洲五月天丁香| 高清午夜精品一区二区三区| 国产精品久久久久久精品电影小说 | 亚洲五月天丁香| 国产熟女欧美一区二区| 能在线免费观看的黄片| 日本与韩国留学比较| 如何舔出高潮| 我的女老师完整版在线观看| 国产淫语在线视频| 我要看日韩黄色一级片| 欧美不卡视频在线免费观看| 久久国内精品自在自线图片| 91午夜精品亚洲一区二区三区| 久久精品国产亚洲网站| 男女下面进入的视频免费午夜| 亚洲国产精品久久男人天堂| 国产欧美另类精品又又久久亚洲欧美| 秋霞伦理黄片| 一夜夜www| 亚洲经典国产精华液单| 老司机影院成人| 亚洲精品日韩在线中文字幕| 中文字幕制服av| 七月丁香在线播放| 成人午夜精彩视频在线观看| 久久99热这里只频精品6学生 | 69人妻影院| 夫妻性生交免费视频一级片| 亚洲国产精品成人久久小说| 只有这里有精品99| 精品99又大又爽又粗少妇毛片| 七月丁香在线播放| 久久这里有精品视频免费| 欧美一区二区亚洲| 欧美最新免费一区二区三区| 男人和女人高潮做爰伦理| 久久久久久九九精品二区国产| 美女高潮的动态| 五月伊人婷婷丁香| 最近最新中文字幕大全电影3| 精品久久久久久久久av| 久久精品国产亚洲av涩爱| 亚洲五月天丁香| 国内精品美女久久久久久| 亚洲欧美一区二区三区国产| av国产免费在线观看| 午夜久久久久精精品| 国产黄a三级三级三级人| 久久亚洲国产成人精品v| 亚洲aⅴ乱码一区二区在线播放| 中文亚洲av片在线观看爽| 国产精品久久久久久精品电影小说 | 久久久国产成人免费| 日本免费一区二区三区高清不卡| 免费看av在线观看网站| 久久99热这里只有精品18| 国产成人精品一,二区| 狂野欧美激情性xxxx在线观看| 九草在线视频观看| 老司机影院成人| 自拍偷自拍亚洲精品老妇| 欧美高清性xxxxhd video| 精品午夜福利在线看| 国产亚洲精品av在线| 免费搜索国产男女视频| 身体一侧抽搐| 麻豆精品久久久久久蜜桃| 少妇丰满av| 十八禁国产超污无遮挡网站| 欧美性猛交╳xxx乱大交人| 国产精品国产高清国产av| 国产亚洲精品久久久com| 蜜臀久久99精品久久宅男| 亚洲不卡免费看| 久久人人爽人人爽人人片va| 久久精品国产亚洲av涩爱| 一本久久精品| 国产亚洲5aaaaa淫片| av线在线观看网站| 国产亚洲午夜精品一区二区久久 | 欧美日本视频| 午夜福利在线观看吧| 精品人妻偷拍中文字幕| 99久国产av精品| 国产成年人精品一区二区| 国产伦理片在线播放av一区| 一区二区三区乱码不卡18| 亚洲av中文av极速乱| 精品99又大又爽又粗少妇毛片| 国产探花极品一区二区| 亚洲va在线va天堂va国产| 边亲边吃奶的免费视频| 91精品国产九色| 亚洲一级一片aⅴ在线观看| 在线播放无遮挡| av国产久精品久网站免费入址| 久久久亚洲精品成人影院| 中文字幕免费在线视频6| 变态另类丝袜制服| 精品酒店卫生间| 麻豆成人午夜福利视频| 欧美bdsm另类| av专区在线播放| 亚洲aⅴ乱码一区二区在线播放| 国产免费一级a男人的天堂| 国产一级毛片在线| 亚洲国产欧洲综合997久久,| 亚洲欧美日韩东京热| 久久久色成人| www日本黄色视频网| 欧美一级a爱片免费观看看| 校园人妻丝袜中文字幕| 午夜福利在线观看免费完整高清在| 韩国高清视频一区二区三区| 七月丁香在线播放| 亚洲性久久影院| 一级毛片我不卡| 日产精品乱码卡一卡2卡三| 级片在线观看| 欧美一区二区亚洲| 日本-黄色视频高清免费观看| 干丝袜人妻中文字幕| av又黄又爽大尺度在线免费看 | 亚洲精品国产成人久久av| 波野结衣二区三区在线| 亚洲人与动物交配视频| 色综合亚洲欧美另类图片| 97超碰精品成人国产| 18禁裸乳无遮挡免费网站照片| 只有这里有精品99| 99热这里只有是精品50| 久久精品久久精品一区二区三区| 一区二区三区高清视频在线| 内地一区二区视频在线| 国产亚洲最大av| 美女脱内裤让男人舔精品视频| 一级黄片播放器| 亚洲伊人久久精品综合 | 精品久久久久久久久av| 久久久久久久久久久免费av| 秋霞在线观看毛片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 秋霞在线观看毛片| 美女黄网站色视频| 国产成人freesex在线| 男人和女人高潮做爰伦理| 国产一区二区在线观看日韩| 麻豆一二三区av精品| videossex国产| 久久99蜜桃精品久久| 国产免费又黄又爽又色| 99久久无色码亚洲精品果冻| 婷婷色麻豆天堂久久 | 青春草亚洲视频在线观看| 亚洲国产高清在线一区二区三| 国产一级毛片在线| 亚洲欧洲日产国产| 亚洲经典国产精华液单| 99热这里只有是精品在线观看| 国产精品人妻久久久影院| 天天一区二区日本电影三级| 国产精品美女特级片免费视频播放器| 免费av毛片视频| 欧美人与善性xxx| 国产精品久久久久久精品电影小说 | 桃色一区二区三区在线观看| 国产高清视频在线观看网站| 伦理电影大哥的女人| 狂野欧美白嫩少妇大欣赏| 色吧在线观看| 亚洲国产精品合色在线| 非洲黑人性xxxx精品又粗又长| 亚洲一级一片aⅴ在线观看| 国内精品一区二区在线观看| 在线免费观看不下载黄p国产| 亚洲一级一片aⅴ在线观看| 人妻制服诱惑在线中文字幕| 午夜激情欧美在线| 久久久久久久久久黄片| 精品酒店卫生间| 国产乱来视频区| 久久人人爽人人爽人人片va| 夫妻性生交免费视频一级片| 又爽又黄无遮挡网站| 久久久久久国产a免费观看| 精品国产露脸久久av麻豆 | av国产久精品久网站免费入址| 国产午夜福利久久久久久| 午夜激情欧美在线| 免费观看a级毛片全部| 国产又黄又爽又无遮挡在线| 国产男人的电影天堂91| 可以在线观看毛片的网站| 亚洲av一区综合| 色网站视频免费| 午夜激情福利司机影院| 热99在线观看视频|