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

    蒙特卡羅自動建模中的空腔簡化算法研究

    2015-03-20 08:21:14俞盛朋程夢云龍鵬程胡麗琴
    原子能科學技術(shù) 2015年1期
    關(guān)鍵詞:蒙特卡羅空腔復雜度

    俞盛朋,程夢云,龍鵬程,胡麗琴

    (中國科學院 核能安全技術(shù)研究所,中國科學院 中子輸運理論與輻射安全重點實驗室,安徽 合肥 230031)

    蒙特卡羅粒子輸運計算方法通過在幾何空間中追蹤粒子與材料之間發(fā)生的物理反應來近乎真實地模擬粒子的運動規(guī)律。為能精確高效進行蒙特卡羅粒子輸運計算,需準確簡潔描述蒙特卡羅計算的幾何模型。傳統(tǒng)的手工蒙特卡羅建模方法費時費力,并可能引入幾何錯誤,限制了蒙特卡羅程序的使用。為解決蒙特卡羅的建模困難,國內(nèi)外已研發(fā)了一系列基于CAD的蒙特卡羅自動建模工具,如德國KIT研發(fā)的McCad[1]、日本JAEA研發(fā)的GEOMIT[2]和FDS團隊研發(fā)的SuperMC/MCAM[3-4]。其中SuperMC/MCAM(multi-physics coupling analysis modeling program)是FDS團隊研發(fā)的超級蒙特卡羅核計算仿真軟件系統(tǒng)SuperMC[5-6]的建模模塊,是一款基于現(xiàn)代軟件工程理念[7]設計開發(fā)的多物理耦合分析自動建模軟件系統(tǒng),其目的是解決核分析中的多物理耦合建模挑戰(zhàn)。目前SuperMC/MCAM 已廣泛應用于國內(nèi)外核相關(guān)領(lǐng)域中[8-17]。

    出于計算效率和正確性的考慮,在典型蒙特卡羅程序,如MCNP[18]中,所有的空間包括實體間空隙(空腔)均被描述為柵元。由于一般的工程CAD 模型僅建立代表實體的幾何,對于空腔幾何,需要由建模軟件自動生成,否則需要大量的人力手工創(chuàng)建代表空腔的CAD 模型[3],并且在手工建立空腔模型的過程中仍有導入錯誤的風險[19]。由于CAD 模型的轉(zhuǎn)換表達式相對固定,自動建模軟件中自動空腔生成方法的優(yōu)劣決定了整體模型轉(zhuǎn)換效率以及最終蒙特卡羅輸運計算的效率。

    目前在自動建模軟件中已實現(xiàn)了基本的自動空腔填充方法,如在SuperMC/MCAM 中已實現(xiàn)了基于CAD實體的自動空腔填充算法[3-4,19],該方法能自動使用子空間表達式減去與其相交的實體表達式描述子空間內(nèi)的空腔,但其對空腔的描述復雜度受限于所相交實體的描述復雜度,在復雜模型中無法生成高效簡潔的空腔描述,限制了其在復雜模型上的應用。為進一步優(yōu)化空腔生成的效率,本文發(fā)展一種基于凸實體的自動空腔生成算法并基于SuperMC/MCAM 最新版本SuperMC/MCAM5.2進行實現(xiàn),基于國際熱核聚變實驗堆ITER 中子學基準模型,對新的空腔生成算法進行測試。

    1 基于凸實體的空腔生成算法

    基于凸實體的空腔簡化算法是使用較實體幾何體更小單位的幾何體進行空腔相交計算,通過描述更少的相交幾何體部分使空腔描述更簡潔,以提高空腔生成的速度及輸運計算中粒子徑跡跟蹤的效率。該算法的基本流程如圖1所示,對其描述如下:

    1)創(chuàng)建包裹所有實體的長方體a,形成初始的空腔長方體集合A={a},假設所有的實體s組成實體集合S。

    2)將所有實體切分為凸實體c,形成整體凸實體集合,若凸實體由原實體外殼面組成,則凸實體為正凸實體,記其集合為Cout,否則凸實體屬于內(nèi)凸實體集合Cin。凸實體可描述實體,如對 于s∈S,s=(∪c∩s=candc∈CoutDesc(c))-(∪c∩s=sandc∈CinDesc(c)),其中,Desc(c)指實 體c的蒙特卡羅幾何描述。

    圖1 基于凸實體的空腔生成算法流程Fig.1 Flow chart of void generation algorithm based on convex solids

    3)執(zhí)行如下步驟,直到A=?。

    (1)從空腔長方體集合中取出空腔長方體a∈A,檢查與a 相交的兩類凸實體,形成兩類凸實體集合:與a相交的外殼集合Cint,out={c|c∈Coutand c∩a=?}以及與a 相交的內(nèi)殼集合Cint,in={c|c∈Cinand c∩a≠?},a 中的空腔部分設為空腔va,其表達式為:Desc(va)=(Desc(a)- (∪c∈CoutDesc(c)))∪(Desc(a)∩(∪c∈Cint,inDesc(c)))。

    (2)判斷a 內(nèi)空腔的描述Desc(va),若達到預定的要求(空腔表達式整體描述長度小于閾值),則生成該空腔的表達式,并將a 從A 中刪除并轉(zhuǎn)到步驟1,否則轉(zhuǎn)到步驟3。

    (3)取a最長一邊作為切分邊,以通過該邊中點且垂直于該邊的平面切分空腔長方體形成空腔長方體anew1和anew2,加入A,轉(zhuǎn)到步驟1。

    相比SuperMC/MCAM 中原空腔算法,基于凸實體的空腔生成算法(新算法)提升了模型轉(zhuǎn)化的效率,主要原因如下:1)新算法只需判斷凸實體與空腔長方體是否相交,而原算法需判斷空腔長方體與整個實體相交與否,相交的判斷和參與判斷的實體復雜度成正比,所以相比原算法,新算法減少了空腔生成過程的幾何計算量,提升了空腔生成效率;2)新算法采用空腔長方體減凸實體描述方式,該描述方式較原描述方式減少了實體中不與空腔長方體相交的部分的描述,故新算法的空腔描述復雜度降低,因此采用較少的分割即可達到與原算法同樣的空腔復雜度要求,這也提升了空腔生成效率。同時,新空腔生成算法相比原算法提升了蒙特卡羅計算效率,主要原因如下:1)新空腔算法相比原空腔算法,建立的空腔的極限復雜度更低,降低了粒子與空腔柵元單次計算的復雜度;2)新空腔算法相比原空腔算法,在相同復雜度下,減少了粒子所需計算的空腔數(shù)量,減少了計算次數(shù)。

    2 測試驗證

    本文在SuperMC/MCAM 中實現(xiàn)新的空腔算法。為證實新空腔算法的正確性以及效率提升,采用ITER Alite模型,通過MCNP建模與MCNP計算兩方面對比新算法與基本算法,進行正確性和效率的測試。圖2為ITER Alite模型,該模型實體部分包含4 194個實體柵元,32 683 個 面,是ITER 于2009 年 發(fā) 布 的 用 于ITER 中子學分析的基準模型。該模型包括了托卡馬克本體和部分生物屏蔽層,考慮到模型的計算效率和托卡馬克裝置的對稱性,該模型僅一個40°的扇形區(qū)域,通過反射面模擬全反應堆裝置計算。

    圖2 SuperMC/MCAM 中的ITER Alite三維視圖Fig.2 3Dview of ITER Alite model in SuperMC/MCAM

    2.1 正確性測試

    原空腔生成算法的正確性已得到證明且已應 用于ITER 建 模 分 析 中[10,17,19],并 已 成 功 應用于ITER Alite的自動建模中[20]。為驗證新算法的正確性,本文基于ITER Alite的CAD模型,分別使用基于凸實體的空腔生成算法和原空腔生成算法生成MCNP模型,在兩個模型上分別進行MCNP計算。在計算模型中,于坐標(600,0,0)處設置各向同性能量為14 MeV的點源,進行輸運計算并統(tǒng)計模型各柵元中的粒子徑跡長度。使用原空腔建模算法和新空腔建模算法的計算結(jié)果的相對偏差如圖3所示。

    在兩種空腔建模算法得到的模型的計算結(jié)果中,柵元徑跡長度的計算結(jié)果的最大相對偏差為2.80×10-4,在允許范圍內(nèi)。實驗結(jié)果證實了新空腔算法建模的正確性。

    圖3 兩種空腔算法建模后計算結(jié)果的相對偏差Fig.3 Relative deviation of calculation result of models generated by two void generation algorithms

    2.2 效率測試

    為對比新、舊空腔建模算法的建模效率,分別使用不同的復雜度參數(shù)生成一系列模型。其中空腔復雜度指描述空腔使用的面的數(shù)量,此參數(shù)限制了生成的空腔的最大復雜度。兩種空腔生成算法的空腔生成時間隨空腔復雜度變化的趨勢如圖4所示。

    圖4 兩種算法空腔建模時間與空腔復雜度的關(guān)系Fig.4 Relationship between void modeling time and void complexity of two algorithms

    空腔復雜度為500以下時,新空腔生成算法的建模時間均小于原空腔算法,在二者均能達到的空腔復雜度區(qū)域,新空腔生成算法的效率是原算法的1.1~25倍;原空腔建模算法在平均空腔復雜度為251時已達極限,而新空腔算法的建模極限是空腔復雜度為49,這證實了新空腔算法對空腔復雜度的可擴展性。

    為驗證新、舊空腔算法所建立模型的計算效率,對比了同樣空腔復雜度情況下新、舊算法生成的空腔數(shù)量,從模型整體復雜度的角度來側(cè)面對比兩種算法生成模型的計算效率。圖5為新、舊算法的空腔復雜度與空腔數(shù)量的關(guān)系。

    圖5 兩種算法空腔數(shù)量與空腔復雜度的關(guān)系Fig.5 Relationship between void number and void complexity of two algorithms

    在同樣的空腔復雜度下,新算法的空腔數(shù)量均低于舊算法。在二者有交集的空腔復雜度范圍內(nèi),新空腔生成算法的空腔數(shù)量是原空腔算法的1/3~1/60。

    為正面驗證新空腔建模算法得到的模型在計算性能上的提升,對不同復雜度的模型的計算速度進行測試。測試均使用了真空材料,點源為設置在(600,0,0)處能量為14 MeV 的中子源,模擬1千萬粒子歷史的測試結(jié)果如圖6所示。

    圖6 兩種算法空腔復雜度與計算時間的關(guān)系Fig.6 Relationship between void complexity andcalculation time of two algorithms

    從圖6可看出,在同樣的空腔復雜度下,由于新空腔算法可生成更少的空腔數(shù)量,其計算效率優(yōu)于原空腔算法;而由于新空腔算法可在很好地控制空腔數(shù)量的前提下得到更低的空腔復雜度,所以其全局最優(yōu)計算效率要高于原空腔算法。試驗得到,新空腔算法生成的模型中計算速度最快的時間為26.9min,而原空腔算法生成的模型最快計算時間為64.91 min,新空腔算法提高了1.4倍計算效率。

    3 總結(jié)與展望

    本文提出了一種新的基于凸實體的自動空腔生成算法,將參與空腔描述的單位從柵元縮小為組成柵元的凸實體,這使得空腔自動建模過程可在控制空腔數(shù)量的前提下得到較優(yōu)的空腔復雜度。該方法實現(xiàn)后已集成進SuperMC/MCAM 中。通過測試對比發(fā)現(xiàn),對應不同的空腔復雜度,新算法的建模效率提升最多達25倍,使生成的空腔復雜度簡化至1/60。對比兩種算法各自生成的最高效的計算模型,新算法較原算法提升了1.4 倍計算效率。

    感謝FDS團隊其他成員提供的幫助和支持。

    [1] TSIGE-TAMIRAT H,F(xiàn)ISCHER U,SERIKOV A,et al.Automatic generation and validation of an ITER neutronics model from CAD data[J].Fusion Engineering and Design,2007,82:1 956-1 959.

    [2] SATO S,IIDA H,NISHITANI T.Development of CAD/MCNP interface program prototype for fusion reactor nuclear analysis[J].Fusion Engineering and Design,2006,81:2 767-2 772.

    [3] WU Y.CAD-based interface programs for fusion neutron transport simulation[J].Fusion Engineering and Design,2009,84(7-11):1 987-1 992.

    [4] 吳宜燦,李瑩,盧磊,等.蒙特卡羅粒子輸運計算自動建模程序系統(tǒng)的研究與發(fā)展[J].核科學與工程,2006,26(1):20-27.WU Yican,LI Ying,LU Lei,et al.Research and development of the automatic modeling system for Monte Carlo particle transport simulation[J].Chinese Journal of Nuclear Science and Engineering,2006,26(1):20-27(in Chinese).

    [5] WU Y C,SONG J,ZHENG H Q,et al.CADbased Monte Carlo program for integrated simulation of nuclear system SuperMC[J].Annals of Nuclear Energy,doi:10.1016/j.anucene.2014.08.058.

    [6] SONG Jing,SUN Guangyao,CHEN Zhenping,et al.Benchmarking of CAD-based SuperMC with ITER benchmark model[J].Fusion Engineering and Design,2014,89(11):2 499-2 503.

    [7] 吳宜燦,胡麗琴,龍鵬程,等.先進核能軟件發(fā)展與核信息學實踐[M]∥中國科研信息化藍皮書.北京:科學出版社,2013:232-244.

    [8] HU H,WU Y,CHEN M,et al.Benchmarking of SNAM with the ITER 3D model[J].Fusion Engineering and Design,2007,82:2 867-2 871.

    [9] BOTH J P,MAZZOLO A,PENELIAU Y,et al.User manual for version 4.3of the TRIPOLI-4 Monte-Carlo method particle transport computer code[M].France:CEA,2003.

    [10]LI Y,LU L,DING A,et al.Benchmarking of MCAM 4.0with the ITER 3Dmodel[J].Fusion Engineering and Design,2007,82(15-24):2 861-2 866.

    [11]吳宜燦,李靜驚,李瑩,等.大型集成多功能中子學計算與分析系統(tǒng)VisualBUS 的研究與發(fā)展[J].核科學與工程,2007,27(4):365-373.WU Yican,LI Jingjing,LI Ying,et al.An integrated multi-functional neutronics calculation and analysis code system:VisualBUS[J].Chinese Journal of Nuclear Science and Engineering,2007,27(4):365-373(in Chinese).

    [12]WU Y,F(xiàn)DS Team.Conceptual design activities of FDS series fusion power plants in China[J].Fusion Engineering and Design,2006,81(23-24):2 713-2 718.

    [13]WU Y,QIAN J,YU J.The fusion-driven hybrid system and its material selection[J].Journal of Nuclear Materials,2002,307-311:1 629-1 636.

    [14]WU Y,F(xiàn)DS Team.Design analysis of the China Dual-Functional Lithium Lead(DFLL)test blanket module in ITER[J].Fusion Engineering and Design,2007,82:1 893-1 903.

    [15]CHEN Y,WU Y.Conceptual study on high performance blanket in a spherical Tokamak fusiondriven transmitter[J].Fusion Engineering and Design,2000,49-50:507-512.

    [16]吳宜燦,柏云清,宋勇,等.中國鉛基研究反應堆概念設計研究[J].核科學與工程,2014,34(2):56-63.WU Yican,BAI Yunqing,SONG Yong,et al.Conceptual design of China Lead-based Research Reactor CLEAR-Ⅰ[J].Nuclear Science and Engineering,2014,34(2):56-63(in Chinese).

    [17]王國忠,黨同強,熊健,等.MCAM4.8 在ITER建筑大廳中子學建模中的應用[J].核科學與工程,2011,31(4):351-355.WANG Guozhong,DANG Tongqiang,XIONG Jian,et al.Application of MCAM4.8in creating neutronics model for ITER building[J].Chinese Journal of Nuclear Science and Engineering,2011,31(4):351-355(in Chinese).

    [18]BRIESMEISTER J F.MCNP:A general Monte Carlo N-particle transport code,version 4C,LA-13709-M[M].US:LANL,2000.

    [19]李瑩.蒙特卡羅粒子輸運計算自動建模系統(tǒng)預處理關(guān)鍵技術(shù)研究[D].合肥:中國科學院合肥物質(zhì)科學研究院,2009.

    [20]LU L,LEE Y K,ZHANG J,et al.Development of Monte Carlo automatic modeling functions of MCAM for TRIPOLI-ITER application[J].Nuclear Instruments and Methods in Physics Research Section A,2009,605(3):384-387.

    猜你喜歡
    蒙特卡羅空腔復雜度
    基于邊光滑有限元法的二維復合彈性空腔聲振特性分析
    利用蒙特卡羅方法求解二重積分
    智富時代(2019年6期)2019-07-24 10:33:16
    一種低復雜度的慣性/GNSS矢量深組合方法
    求圖上廣探樹的時間復雜度
    空腔參數(shù)對重力壩穩(wěn)定的影響分析
    前置污水去油池
    前置污水去油池
    某雷達導51 頭中心控制軟件圈復雜度分析與改進
    出口技術(shù)復雜度研究回顧與評述
    探討蒙特卡羅方法在解微分方程邊值問題中的應用
    男女啪啪激烈高潮av片| 午夜福利在线观看免费完整高清在 | h日本视频在线播放| 2022亚洲国产成人精品| 九色成人免费人妻av| 国模一区二区三区四区视频| 嘟嘟电影网在线观看| 国产在线男女| ponron亚洲| 国产精品国产三级国产av玫瑰| 免费搜索国产男女视频| 别揉我奶头 嗯啊视频| 少妇熟女aⅴ在线视频| 12—13女人毛片做爰片一| 国产精品av视频在线免费观看| 噜噜噜噜噜久久久久久91| 能在线免费观看的黄片| 亚洲精华国产精华液的使用体验 | 深爱激情五月婷婷| 国产精品av视频在线免费观看| 又黄又爽又刺激的免费视频.| 狂野欧美白嫩少妇大欣赏| 蜜桃亚洲精品一区二区三区| 乱码一卡2卡4卡精品| kizo精华| 亚洲成人av在线免费| 日本免费一区二区三区高清不卡| 最好的美女福利视频网| 高清日韩中文字幕在线| 免费不卡的大黄色大毛片视频在线观看 | 色吧在线观看| 亚洲一级一片aⅴ在线观看| 国产高清三级在线| 麻豆精品久久久久久蜜桃| 欧美成人免费av一区二区三区| 中文字幕免费在线视频6| 丝袜喷水一区| 国产成人freesex在线| 国产成人精品久久久久久| 亚洲国产精品sss在线观看| 久久精品国产自在天天线| 国产一区二区三区在线臀色熟女| 亚洲高清免费不卡视频| 卡戴珊不雅视频在线播放| 一级毛片我不卡| 男人和女人高潮做爰伦理| 亚洲人与动物交配视频| 国产麻豆成人av免费视频| 免费在线观看成人毛片| 中出人妻视频一区二区| 国产日韩欧美在线精品| 99久久精品一区二区三区| 亚洲成av人片在线播放无| 亚洲在线观看片| 国产探花在线观看一区二区| 大又大粗又爽又黄少妇毛片口| 99久久精品一区二区三区| 久久久精品欧美日韩精品| 国产视频首页在线观看| 免费大片18禁| 亚洲av一区综合| 国产黄片视频在线免费观看| 日本爱情动作片www.在线观看| 18禁在线播放成人免费| 日日干狠狠操夜夜爽| 91麻豆精品激情在线观看国产| 国产一区亚洲一区在线观看| 国产精品一及| АⅤ资源中文在线天堂| 99久久无色码亚洲精品果冻| 夜夜夜夜夜久久久久| 日韩视频在线欧美| 亚洲,欧美,日韩| 亚洲人成网站在线播放欧美日韩| 欧美日韩在线观看h| a级毛片免费高清观看在线播放| 男人和女人高潮做爰伦理| 真实男女啪啪啪动态图| 成人高潮视频无遮挡免费网站| 日本与韩国留学比较| 五月玫瑰六月丁香| 91久久精品国产一区二区成人| 婷婷色av中文字幕| av在线亚洲专区| 欧美+亚洲+日韩+国产| 麻豆国产97在线/欧美| 日本欧美国产在线视频| 丝袜喷水一区| 亚洲中文字幕一区二区三区有码在线看| 美女黄网站色视频| 午夜老司机福利剧场| 直男gayav资源| 老熟妇乱子伦视频在线观看| 看非洲黑人一级黄片| 国产高清有码在线观看视频| 国产单亲对白刺激| 熟妇人妻久久中文字幕3abv| 亚洲乱码一区二区免费版| 听说在线观看完整版免费高清| 三级毛片av免费| 国产探花在线观看一区二区| 成人毛片60女人毛片免费| 亚洲国产欧洲综合997久久,| 国内精品宾馆在线| 久久久久久久亚洲中文字幕| 日本黄色视频三级网站网址| 91麻豆精品激情在线观看国产| 久久久久九九精品影院| 国产 一区 欧美 日韩| 熟女电影av网| 国产高清视频在线观看网站| 亚洲四区av| 久久久a久久爽久久v久久| 女同久久另类99精品国产91| 99热网站在线观看| 国产成人影院久久av| 久久久a久久爽久久v久久| 亚洲中文字幕日韩| 深夜精品福利| 国产淫片久久久久久久久| 成人一区二区视频在线观看| 91久久精品电影网| 成年女人永久免费观看视频| 乱系列少妇在线播放| 国产av麻豆久久久久久久| 午夜福利高清视频| 国产精品日韩av在线免费观看| 能在线免费观看的黄片| 人人妻人人看人人澡| 自拍偷自拍亚洲精品老妇| 深爱激情五月婷婷| 男人狂女人下面高潮的视频| 国内精品一区二区在线观看| 色噜噜av男人的天堂激情| 嫩草影院入口| 级片在线观看| 国产一级毛片七仙女欲春2| 国产成人精品婷婷| 欧美xxxx性猛交bbbb| 国产亚洲精品av在线| 久久99精品国语久久久| 国产精品人妻久久久影院| 18禁裸乳无遮挡免费网站照片| 99热这里只有精品一区| 欧美人与善性xxx| 日日啪夜夜撸| 激情 狠狠 欧美| 日韩亚洲欧美综合| 免费黄网站久久成人精品| 午夜免费男女啪啪视频观看| 天堂中文最新版在线下载 | 免费无遮挡裸体视频| 自拍偷自拍亚洲精品老妇| 久久久精品94久久精品| 一本一本综合久久| 国产精品蜜桃在线观看 | 国产成人精品久久久久久| 我的女老师完整版在线观看| 国产黄色小视频在线观看| 51国产日韩欧美| 啦啦啦观看免费观看视频高清| 高清午夜精品一区二区三区 | 亚洲国产色片| 极品教师在线视频| 日韩视频在线欧美| 亚洲成人av在线免费| 黄片无遮挡物在线观看| 18禁黄网站禁片免费观看直播| 久久精品国产亚洲av香蕉五月| 久久久久久久午夜电影| videossex国产| 亚洲av电影不卡..在线观看| 在线播放无遮挡| 国产三级中文精品| 国产国拍精品亚洲av在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲最大成人中文| 免费观看精品视频网站| 国产在线男女| 久久精品国产亚洲av天美| 亚洲av不卡在线观看| 91精品国产九色| 国产精品久久电影中文字幕| 精品久久久久久成人av| 精品日产1卡2卡| 国产久久久一区二区三区| 久久精品国产鲁丝片午夜精品| 国产乱人偷精品视频| 日本一本二区三区精品| 国产在视频线在精品| 成人综合一区亚洲| 日韩一区二区视频免费看| 赤兔流量卡办理| 国产精品乱码一区二三区的特点| 免费看光身美女| 国产精品.久久久| 欧美区成人在线视频| 亚洲成人久久性| 亚洲精品乱码久久久久久按摩| 亚洲欧美精品自产自拍| 亚洲av.av天堂| 六月丁香七月| 直男gayav资源| 婷婷精品国产亚洲av| 日韩精品有码人妻一区| 欧美一区二区亚洲| 嘟嘟电影网在线观看| 中国美白少妇内射xxxbb| 男女啪啪激烈高潮av片| 亚洲自偷自拍三级| 最近手机中文字幕大全| 男人狂女人下面高潮的视频| 亚洲性久久影院| 99国产极品粉嫩在线观看| 黄色视频,在线免费观看| 51国产日韩欧美| 精品人妻一区二区三区麻豆| 看片在线看免费视频| 国产熟女欧美一区二区| 亚洲av.av天堂| 少妇猛男粗大的猛烈进出视频 | 一级毛片电影观看 | 边亲边吃奶的免费视频| 欧美bdsm另类| 国产黄色小视频在线观看| 51国产日韩欧美| 天堂av国产一区二区熟女人妻| 22中文网久久字幕| 国产男人的电影天堂91| 毛片一级片免费看久久久久| 99久国产av精品国产电影| 国内久久婷婷六月综合欲色啪| 久久久久免费精品人妻一区二区| 国产精品综合久久久久久久免费| 国产私拍福利视频在线观看| 国模一区二区三区四区视频| 亚洲天堂国产精品一区在线| 黄色日韩在线| 国产不卡一卡二| 97在线视频观看| 老司机影院成人| 啦啦啦观看免费观看视频高清| 91久久精品电影网| 长腿黑丝高跟| 性插视频无遮挡在线免费观看| 18禁在线无遮挡免费观看视频| 91在线精品国自产拍蜜月| 亚洲av一区综合| 在线观看美女被高潮喷水网站| 偷拍熟女少妇极品色| 青春草亚洲视频在线观看| 一夜夜www| 一级毛片久久久久久久久女| 色播亚洲综合网| 午夜老司机福利剧场| 亚洲欧美日韩高清专用| 欧美色欧美亚洲另类二区| 欧美色欧美亚洲另类二区| 欧美成人免费av一区二区三区| 一本精品99久久精品77| 大又大粗又爽又黄少妇毛片口| 中文字幕av在线有码专区| 久久久久久久久中文| 欧美激情国产日韩精品一区| 2021天堂中文幕一二区在线观| 91久久精品国产一区二区三区| 最近的中文字幕免费完整| 欧美人与善性xxx| 99久久精品热视频| 美女cb高潮喷水在线观看| 日韩一区二区视频免费看| 日本与韩国留学比较| 99久久无色码亚洲精品果冻| 国产亚洲av片在线观看秒播厂 | 亚洲精品色激情综合| 久久国内精品自在自线图片| 亚洲在线观看片| 国产熟女欧美一区二区| 国产 一区精品| 激情 狠狠 欧美| 最近2019中文字幕mv第一页| 日韩欧美在线乱码| 波多野结衣高清无吗| av在线天堂中文字幕| 日本撒尿小便嘘嘘汇集6| 中文字幕久久专区| 日本成人三级电影网站| 99国产极品粉嫩在线观看| 搡女人真爽免费视频火全软件| 淫秽高清视频在线观看| 在线免费观看不下载黄p国产| 91精品一卡2卡3卡4卡| 国产免费一级a男人的天堂| 亚洲精品色激情综合| 国产亚洲av嫩草精品影院| 国产久久久一区二区三区| 免费一级毛片在线播放高清视频| 国产极品精品免费视频能看的| 精品一区二区免费观看| 女人十人毛片免费观看3o分钟| 51国产日韩欧美| 国产日韩欧美在线精品| 成年版毛片免费区| 欧美日韩在线观看h| av在线天堂中文字幕| 在线免费观看不下载黄p国产| a级一级毛片免费在线观看| 人人妻人人澡欧美一区二区| av在线天堂中文字幕| 麻豆乱淫一区二区| 欧美最黄视频在线播放免费| 全区人妻精品视频| 看非洲黑人一级黄片| 日本免费a在线| 成年av动漫网址| 99视频精品全部免费 在线| 国产精品,欧美在线| 精品免费久久久久久久清纯| 少妇熟女欧美另类| 99久久精品一区二区三区| 国产精品久久久久久久电影| 精华霜和精华液先用哪个| 欧美一区二区亚洲| 色视频www国产| 伊人久久精品亚洲午夜| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲欧美日韩东京热| 久久婷婷人人爽人人干人人爱| 亚洲最大成人手机在线| 欧美变态另类bdsm刘玥| 国产精品无大码| 日日摸夜夜添夜夜爱| 在线播放国产精品三级| 欧美日韩一区二区视频在线观看视频在线 | 乱系列少妇在线播放| 久久久精品大字幕| 亚洲av中文av极速乱| 久久6这里有精品| 国产老妇伦熟女老妇高清| 亚洲av成人精品一区久久| 天天躁夜夜躁狠狠久久av| 美女xxoo啪啪120秒动态图| 免费看日本二区| 小说图片视频综合网站| 热99在线观看视频| 久久久色成人| kizo精华| 国产精品久久久久久精品电影| 午夜激情欧美在线| 国产黄片美女视频| 久久人人爽人人爽人人片va| 男人狂女人下面高潮的视频| 亚洲国产高清在线一区二区三| 欧美性感艳星| 午夜精品国产一区二区电影 | 亚洲国产日韩欧美精品在线观看| 国产高清激情床上av| 成人av在线播放网站| 国产精品国产高清国产av| 亚洲精品粉嫩美女一区| 久久久久久久久久久丰满| 天天躁夜夜躁狠狠久久av| 日日摸夜夜添夜夜爱| av专区在线播放| 久久精品国产亚洲av香蕉五月| 国产中年淑女户外野战色| 国产老妇女一区| 18禁在线播放成人免费| 十八禁国产超污无遮挡网站| 亚洲欧美日韩无卡精品| 久久久色成人| 久久久久久久久久久免费av| 国产一区二区三区在线臀色熟女| 国产v大片淫在线免费观看| 女的被弄到高潮叫床怎么办| 在线播放无遮挡| 国产精品免费一区二区三区在线| 高清毛片免费观看视频网站| 精品人妻一区二区三区麻豆| 亚洲真实伦在线观看| 一本久久中文字幕| 天堂网av新在线| 国产精品久久久久久亚洲av鲁大| 午夜老司机福利剧场| 赤兔流量卡办理| 国产精品一及| 桃色一区二区三区在线观看| 国产成人午夜福利电影在线观看| 午夜视频国产福利| 五月伊人婷婷丁香| 最近手机中文字幕大全| 午夜视频国产福利| 久久精品久久久久久噜噜老黄 | 日韩中字成人| 青青草视频在线视频观看| 蜜桃亚洲精品一区二区三区| 欧美日韩综合久久久久久| АⅤ资源中文在线天堂| 99热全是精品| 看免费成人av毛片| 免费不卡的大黄色大毛片视频在线观看 | 国产亚洲精品久久久com| 午夜精品一区二区三区免费看| 蜜臀久久99精品久久宅男| 99久久精品国产国产毛片| 国产精品99久久久久久久久| 国国产精品蜜臀av免费| 中文字幕免费在线视频6| 亚洲天堂国产精品一区在线| 夜夜看夜夜爽夜夜摸| 日韩强制内射视频| 久久久久久久久久久丰满| 在线观看av片永久免费下载| 岛国在线免费视频观看| 久久久a久久爽久久v久久| 国产av在哪里看| 高清毛片免费观看视频网站| 国产免费男女视频| 成人特级av手机在线观看| 亚洲第一电影网av| 亚洲第一区二区三区不卡| 国产精品女同一区二区软件| 又黄又爽又刺激的免费视频.| av免费观看日本| 一进一出抽搐动态| 插逼视频在线观看| 亚洲国产精品成人综合色| 成人综合一区亚洲| 国产不卡一卡二| 国产真实乱freesex| 少妇熟女aⅴ在线视频| 美女国产视频在线观看| 精品一区二区三区视频在线| 特级一级黄色大片| 99国产极品粉嫩在线观看| 国产人妻一区二区三区在| 九草在线视频观看| 淫秽高清视频在线观看| 国产精品麻豆人妻色哟哟久久 | 免费观看的影片在线观看| 看十八女毛片水多多多| 校园人妻丝袜中文字幕| 久久这里有精品视频免费| 国产精品一及| 女的被弄到高潮叫床怎么办| 国产亚洲5aaaaa淫片| 国产黄色视频一区二区在线观看 | 全区人妻精品视频| 岛国在线免费视频观看| 国产精品电影一区二区三区| 日本五十路高清| 欧美xxxx性猛交bbbb| 欧美精品一区二区大全| 97热精品久久久久久| 伦精品一区二区三区| 日韩视频在线欧美| 女同久久另类99精品国产91| 一区二区三区四区激情视频 | 有码 亚洲区| 久久人妻av系列| 菩萨蛮人人尽说江南好唐韦庄 | 六月丁香七月| 国产探花极品一区二区| 成人特级av手机在线观看| 亚洲在线自拍视频| 精品一区二区免费观看| 久久精品国产99精品国产亚洲性色| 日本黄大片高清| 97热精品久久久久久| 日本熟妇午夜| 色哟哟·www| 丝袜喷水一区| av免费在线看不卡| 亚洲国产精品久久男人天堂| 午夜福利高清视频| 好男人在线观看高清免费视频| www.av在线官网国产| 亚洲av一区综合| 国产伦精品一区二区三区四那| av免费在线看不卡| 99热这里只有精品一区| 一区二区三区高清视频在线| 美女cb高潮喷水在线观看| 夫妻性生交免费视频一级片| 女人十人毛片免费观看3o分钟| 色5月婷婷丁香| 男女边吃奶边做爰视频| 十八禁国产超污无遮挡网站| 久久久久久久久久黄片| 中出人妻视频一区二区| 亚洲欧美日韩无卡精品| 日韩一区二区三区影片| 91久久精品电影网| 麻豆国产av国片精品| 亚洲欧美精品自产自拍| 亚洲va在线va天堂va国产| 三级男女做爰猛烈吃奶摸视频| 超碰av人人做人人爽久久| 亚洲成av人片在线播放无| 丰满乱子伦码专区| 午夜精品一区二区三区免费看| 精品久久久久久久久亚洲| 免费在线观看成人毛片| 精品人妻偷拍中文字幕| 国产精品久久久久久精品电影小说 | 国产精品久久久久久久久免| 久久99精品国语久久久| 午夜福利视频1000在线观看| 夜夜看夜夜爽夜夜摸| 日韩av在线大香蕉| av福利片在线观看| 性插视频无遮挡在线免费观看| 国产精品乱码一区二三区的特点| 一级毛片电影观看 | 国产人妻一区二区三区在| 夜夜看夜夜爽夜夜摸| 免费人成视频x8x8入口观看| 亚洲国产精品合色在线| а√天堂www在线а√下载| 国产伦一二天堂av在线观看| 最近视频中文字幕2019在线8| 国产欧美日韩精品一区二区| 久久精品影院6| 黄色欧美视频在线观看| 国产蜜桃级精品一区二区三区| 亚洲av电影不卡..在线观看| 亚洲精品久久久久久婷婷小说 | 毛片女人毛片| 国内精品美女久久久久久| 国产人妻一区二区三区在| 波多野结衣高清无吗| 性欧美人与动物交配| 精品人妻视频免费看| 极品教师在线视频| 日韩中字成人| 成人国产麻豆网| 一卡2卡三卡四卡精品乱码亚洲| 少妇被粗大猛烈的视频| 亚洲欧美精品综合久久99| av在线观看视频网站免费| 久久久久性生活片| 中文欧美无线码| 国产v大片淫在线免费观看| 男女做爰动态图高潮gif福利片| 国产精品永久免费网站| 色综合亚洲欧美另类图片| 国产av一区在线观看免费| 色综合亚洲欧美另类图片| 一本精品99久久精品77| 黑人高潮一二区| 亚洲在线自拍视频| 国产av不卡久久| 看片在线看免费视频| 听说在线观看完整版免费高清| 国产伦理片在线播放av一区 | 女的被弄到高潮叫床怎么办| www.色视频.com| 啦啦啦观看免费观看视频高清| 丰满乱子伦码专区| 亚洲欧美日韩东京热| 日韩高清综合在线| 插逼视频在线观看| 国内精品宾馆在线| 成年版毛片免费区| 99久久精品国产国产毛片| 一区二区三区高清视频在线| 三级毛片av免费| 深夜a级毛片| 国产精品久久久久久久久免| 亚洲av男天堂| 国产色爽女视频免费观看| 99热网站在线观看| 国产av一区在线观看免费| 日韩欧美三级三区| 又黄又爽又刺激的免费视频.| 亚洲国产欧美在线一区| 午夜精品一区二区三区免费看| 亚洲精品色激情综合| 亚洲三级黄色毛片| 精品国产三级普通话版| 久久精品国产亚洲av天美| 国产精品国产三级国产av玫瑰| 联通29元200g的流量卡| 久久久久久大精品| 久久久久久久久久成人| 欧美zozozo另类| 女同久久另类99精品国产91| 亚洲最大成人手机在线| 国内少妇人妻偷人精品xxx网站| 免费在线观看成人毛片| 久久国产乱子免费精品| 内地一区二区视频在线| 色吧在线观看| 一级毛片aaaaaa免费看小| 噜噜噜噜噜久久久久久91| 搞女人的毛片| 中国美女看黄片| 成人毛片a级毛片在线播放| 日韩一区二区三区影片| 草草在线视频免费看| 日韩大尺度精品在线看网址| 青青草视频在线视频观看| 波野结衣二区三区在线| 久久久色成人| 国产一区二区三区av在线 | 久久国产乱子免费精品| 久久精品国产鲁丝片午夜精品| 亚洲国产欧美在线一区| 久久亚洲国产成人精品v| 中文字幕熟女人妻在线| 伊人久久精品亚洲午夜| videossex国产| 永久网站在线| 非洲黑人性xxxx精品又粗又长| 3wmmmm亚洲av在线观看|