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

    基于Geant4的質(zhì)子輸運(yùn)計(jì)算的加速方法研究

    2015-09-16 02:24:56姚志明黑東煒宋顧周馬繼明韓長(zhǎng)材段寶軍岳志勤
    黑龍江科學(xué) 2015年8期
    關(guān)鍵詞:質(zhì)子個(gè)數(shù)粒子

    姚志明,黑東煒,宋顧周,馬繼明,韓長(zhǎng)材,段寶軍,岳志勤

    (西北核技術(shù)研究所強(qiáng)脈沖輻射環(huán)境模擬與效應(yīng)國(guó)家重點(diǎn)實(shí)驗(yàn)室,西安 710024)

    基于Geant4的質(zhì)子輸運(yùn)計(jì)算的加速方法研究

    姚志明,黑東煒,宋顧周,馬繼明,韓長(zhǎng)材,段寶軍,岳志勤

    (西北核技術(shù)研究所強(qiáng)脈沖輻射環(huán)境模擬與效應(yīng)國(guó)家重點(diǎn)實(shí)驗(yàn)室,西安 710024)

    質(zhì)子在醫(yī)學(xué)、材料科學(xué)、輻射成像等領(lǐng)域有著廣闊的應(yīng)用前景,質(zhì)子輸運(yùn)過(guò)程的數(shù)值模擬是質(zhì)子應(yīng)用理論分析的重要手段。本研究基于Geant4軟件編寫(xiě)了蒙卡模擬程序,研究了程序的加速方法。結(jié)果表明,次級(jí)粒子截?cái)嗪筒⑿杏?jì)算可以有效加快程序的運(yùn)行速度。程序易于擴(kuò)展,為質(zhì)子應(yīng)用的快速模擬奠定了基礎(chǔ)。

    質(zhì)子輸運(yùn);蒙特卡羅;加速方法;次級(jí)粒子截?cái)?;并行?jì)算

    質(zhì)子在醫(yī)學(xué)、材料科學(xué)、輻射成像等領(lǐng)域有著廣泛的應(yīng)用。目前,國(guó)內(nèi)有多臺(tái)質(zhì)子加速器已經(jīng)建成或正在建設(shè)當(dāng)中,它們有著不同的應(yīng)用前景。例如,中國(guó)原子能科學(xué)研究院的回旋加速器上開(kāi)展了很多國(guó)防核技術(shù)和核物理的研究工作[1];上海應(yīng)用物理研究所正在研制用于惡性腫瘤治療的質(zhì)子治療裝置[2];在建的中國(guó)散裂中子源的質(zhì)子加速器計(jì)劃用于空間輻射效應(yīng)模擬和質(zhì)子照相;等等。

    數(shù)值模擬計(jì)算可以為粒子的應(yīng)用提供理論指導(dǎo)和依據(jù)。數(shù)值模擬中的蒙特卡羅方法通過(guò)隨機(jī)抽樣,能夠比較真實(shí)地反映粒子的輸運(yùn)過(guò)程[3]。目前,模擬質(zhì)子輸運(yùn)的軟件包括SRIM、MCNP、MCNPX、Geant4等,其中Geant4允許用戶二次開(kāi)發(fā),容易根據(jù)實(shí)際需要對(duì)程序進(jìn)行改進(jìn)和擴(kuò)展。

    基于 Geant4的模擬質(zhì)子輸運(yùn)的研究比較豐富[4-6],加速方法的研究卻較少。文獻(xiàn)[7,8]中給出了通用的并行計(jì)算加快模擬速度的方法。然而,針對(duì)質(zhì)子輸運(yùn)過(guò)程的特點(diǎn)對(duì)Geant4程序進(jìn)行加速的方法未見(jiàn)報(bào)道。

    本研究以10 GeV的質(zhì)子在FTO (French Test Object)中的輸運(yùn)過(guò)程為例,基于Geant4建立蒙卡模擬程序[9,10],探討加快質(zhì)子輸運(yùn)模擬速度的方法。

    1 質(zhì)子與物質(zhì)的作用規(guī)律

    質(zhì)子與物質(zhì)的作用過(guò)程主要有電離能損、多次庫(kù)侖散射、與原子核的彈性碰撞和非彈性碰撞[11]。電離能損是質(zhì)子與很多核外電子碰撞,每次碰撞電子獲得一部分能量,質(zhì)子的能量隨之減小。多次庫(kù)侖散射是當(dāng)質(zhì)子與原子最近距離比原子半徑小時(shí),受到原子核電場(chǎng)的作用,在庫(kù)侖力作用下運(yùn)動(dòng)方向發(fā)生改變。與原子核碰撞時(shí)質(zhì)子受到核力的作用,按質(zhì)子能量是否發(fā)生改變分為彈性和非彈性碰撞。彈性碰撞按照原子處于基態(tài)和激發(fā)態(tài)可細(xì)分為彈性碰撞和準(zhǔn)彈性碰撞,質(zhì)子的運(yùn)動(dòng)方向發(fā)生較大偏轉(zhuǎn);非彈性碰撞是原子核俘獲質(zhì)子,發(fā)生核反應(yīng),核反應(yīng)過(guò)程中釋放次級(jí)粒子,導(dǎo)致入射質(zhì)子的個(gè)數(shù)損失,并可能產(chǎn)生能量較低的次級(jí)質(zhì)子。本文10 GeV質(zhì)子在FTO輸運(yùn)過(guò)程中以上四種作用都會(huì)發(fā)生,導(dǎo)致質(zhì)子的能量損失、個(gè)數(shù)衰減和方向偏轉(zhuǎn)。其中非彈性碰撞產(chǎn)生的中子、γ射線等次級(jí)粒子對(duì)結(jié)果貢獻(xiàn)小,模擬耗時(shí)多,截?cái)嗪罂梢蕴岣哂?jì)算速度。

    2 模擬程序的建立

    Geant4軟件采用模塊化設(shè)計(jì),由不同的模塊分別實(shí)現(xiàn)探測(cè)器幾何結(jié)構(gòu)定義、粒子與物質(zhì)相互作用過(guò)程描述、跟蹤并記錄粒子的狀態(tài)和改變、中止次級(jí)粒子的產(chǎn)生等[10]。圖1給出了軟件的一些常用的類、子函數(shù)以及它們的功能。

    模擬程序的幾何布局如圖2所示。Construct函數(shù)中定義物體的材料和幾何結(jié)構(gòu)。物體選為FTO[11],它是高能閃光照相中使用的一種通用測(cè)試對(duì)象,是由鎢、銅和塑料組成的多層同心球,內(nèi)半徑分別為1 cm、4.5 cm、6.5 cm,塑料的外半徑為22.5 cm。物體后1 mm處放置1 mm厚的足夠大長(zhǎng)方體,定義為敏感探測(cè)器。Geant4中的物理過(guò)程比較豐富,采用的物理模型也有多種選擇,為方便用戶的使用,Geant4開(kāi)發(fā)人員將物理模型打包,定義了多種physic list,其中包含了質(zhì)子與物質(zhì)作用的全部物理過(guò)程和產(chǎn)生的中子、伽瑪?shù)却渭?jí)粒子的作用過(guò)程。本文選取了physic list中的QGSP_BERT物理過(guò)程列表。GeneratePrimaries函數(shù)中設(shè)置發(fā)射粒子的參數(shù)信息。粒子類型為質(zhì)子,初始能量為 10 GeV,發(fā)射方向沿 z軸正方向。ProcessHits、Print、EndOfEvent函數(shù)中定義敏感探測(cè)器的功能,用于粒子信息的記錄。當(dāng)粒子打到敏感探測(cè)器上,在ProcessHits函數(shù)中將該粒子的出射位置、運(yùn)動(dòng)方向、類型、動(dòng)能和輸運(yùn)編號(hào)存儲(chǔ)在一個(gè)Hit中。在Print函數(shù)中定義Hit中各信息的輸出格式。當(dāng)發(fā)射的一個(gè)質(zhì)子輸運(yùn)完成后,EndOfEvent函數(shù)調(diào)用Print函數(shù)將各信息存儲(chǔ)在文本文件中,供分析使用。

    3 加速方法與結(jié)果分析

    3.1次級(jí)粒子截?cái)?/p>

    質(zhì)子輸運(yùn)過(guò)程中發(fā)生核反應(yīng)會(huì)產(chǎn)生大量中子、γ射線和低能量的質(zhì)子。除低能量的質(zhì)子外,還會(huì)產(chǎn)生一定數(shù)量的電子、μ介子、π介子等其他帶電粒子。這些次級(jí)粒子的輸運(yùn)過(guò)程模擬會(huì)耗費(fèi)大量時(shí)間。Geant4中的G4VUserStackingAction類可以對(duì)粒子的信息進(jìn)行判斷,實(shí)現(xiàn)次級(jí)粒子的截?cái)?,達(dá)到節(jié)省計(jì)算時(shí)間的目的。然而,次級(jí)粒子的截?cái)嗫赡軙?huì)影響質(zhì)子輸運(yùn)的模擬,本文選取了敏感探測(cè)器記錄到的質(zhì)子總個(gè)數(shù)和直穿質(zhì)子個(gè)數(shù)兩個(gè)參數(shù)來(lái)評(píng)估截?cái)鄬?duì)模擬準(zhǔn)確性的影響。

    Geant4軟件運(yùn)行在聯(lián)想ThinkCentre M8380t型號(hào)的計(jì)算機(jī)上,安裝了windows7 32位操作系統(tǒng),CPU頻率為3.4 GHz。

    首先,不截?cái)啻渭?jí)粒子,程序?qū)⒛M所有直穿質(zhì)子和次級(jí)粒子的輸運(yùn)過(guò)程。發(fā)射105個(gè)質(zhì)子,時(shí)間由G4Timer類記錄,模擬用時(shí)2 538 s。表1給出了敏感探測(cè)器記錄的粒子種類和個(gè)數(shù)。直穿質(zhì)子和次級(jí)質(zhì)子,可以由輸運(yùn)編號(hào)區(qū)分,輸運(yùn)編號(hào)為1的是直穿質(zhì)子,大于1的是次級(jí)質(zhì)子。

    表1 不截?cái)喈a(chǎn)生的粒子個(gè)數(shù)Tab.1 Particle generation without cut

    圖1 Geant4常用類和函數(shù)結(jié)構(gòu)圖Fig.1 Common class and function of Geant4

    圖2 模擬程序的幾何布局示意圖Fig.2 Geometry of the simulation program

    然后,在ClassifyNewTrack函數(shù)中分別實(shí)現(xiàn)γ射線、中子和其他帶電粒子的截?cái)?,并考慮了只模擬質(zhì)子和只模擬直穿質(zhì)子的情況。以截?cái)唳蒙渚€為例,ClassifyNewTrack函數(shù)中的程序代碼如下:

    表2給出了各情況下模擬所需的時(shí)間,通過(guò)與表1比較,給出了直穿質(zhì)子個(gè)數(shù)和總質(zhì)子個(gè)數(shù)減少的百分比,并計(jì)算了截?cái)嗲昂髸r(shí)間的比值,即加速比。

    由表2可以得出以下結(jié)論:第一,截?cái)啻渭?jí)粒子可以有效地減少計(jì)算時(shí)間,截?cái)嗖煌牧W訒r(shí)加速比均達(dá)到了2倍以上。第二,截?cái)嗖煌牧W雍髮?duì)直穿質(zhì)子的個(gè)數(shù)的影響均小于1%,可以認(rèn)為不會(huì)影響直穿質(zhì)子的模擬結(jié)果。第三,截?cái)唳蒙渚€對(duì)總質(zhì)子個(gè)數(shù)的影響小于1%,截?cái)嘀凶雍推渌麕щ娏W雍罂傎|(zhì)子個(gè)數(shù)減少了10%以上,分析原因是由于γ射線與物質(zhì)作用后不會(huì)產(chǎn)生質(zhì)子,截?cái)嗖粫?huì)影響質(zhì)子的模擬結(jié)果,而中子和一些帶電粒子與物質(zhì)作用會(huì)進(jìn)一步產(chǎn)生次級(jí)質(zhì)子,截?cái)噙@些粒子會(huì)使得出射的次級(jí)質(zhì)子有所減少,在需要考慮次級(jí)質(zhì)子輸運(yùn)過(guò)程的應(yīng)用中不能采用。

    表2 截?cái)嗪蟮挠?jì)算速度測(cè)試結(jié)果Tab.2 Results of calculation speed testing with cut

    3.2并行化加速方法

    有關(guān)Geant4程序并行化的研究已較為豐富[7,8]。本研究將并行化方法應(yīng)用于質(zhì)子輸運(yùn)模擬,采用文獻(xiàn)[8]中的方法對(duì)程序進(jìn)行了并行化改造。并行計(jì)算在銀河計(jì)算機(jī)上進(jìn)行,該計(jì)算機(jī)安裝了紅帽公司的Linux服務(wù)器操作系統(tǒng),CPU頻率為1.6 GHz。表3給出了并行化后的測(cè)試結(jié)果。需要說(shuō)明的是,當(dāng)模擬的質(zhì)子個(gè)數(shù)較多時(shí),如果將所有質(zhì)子的信息存儲(chǔ)在一個(gè)文件中,得到的數(shù)據(jù)文件將會(huì)很大,不便于后續(xù)數(shù)據(jù)的統(tǒng)計(jì)分析。在print函數(shù)中可以首先讀取CPU的編號(hào),用編號(hào)為輸出文件命名,將不同CPU的計(jì)算結(jié)果存儲(chǔ)在不同的文件當(dāng)中,解決數(shù)據(jù)文件過(guò)大的問(wèn)題。實(shí)現(xiàn)語(yǔ)句如下:

    由表3可以看出,加速比隨CPU個(gè)數(shù)的增加近似成線性增長(zhǎng),與文獻(xiàn)[8]中的結(jié)果一致。與CPU為3.4 GHz的聯(lián)想計(jì)算機(jī)的結(jié)果比較,在只模擬直穿質(zhì)子的條件下,聯(lián)想計(jì)算機(jī)用時(shí)65 s,銀河計(jì)算機(jī)用時(shí)295 s。銀河計(jì)算機(jī)單個(gè)CPU運(yùn)行的速率不及CPU頻率較高的聯(lián)想計(jì)算機(jī),然而,銀河計(jì)算機(jī)上最多可提供的CPU數(shù)達(dá)1 024個(gè),可以根據(jù)需要進(jìn)一步增加CPU的使用個(gè)數(shù)來(lái)提高質(zhì)子輸運(yùn)的模擬速度。

    表3 并行程序的計(jì)算速度測(cè)試結(jié)果Tab.3 Results of calculation speed testing of parallel program

    4 結(jié)論

    本研究基于Geant4軟件建立了模擬10 GeV能量的質(zhì)子穿透FTO過(guò)程的蒙卡程序。重點(diǎn)探討了程序的加速方法,在程序中分別實(shí)現(xiàn)了次級(jí)粒子截?cái)嗪筒⑿杏?jì)算。截?cái)唳?、中子和其他帶電粒子后,加速比均達(dá)到了2倍以上,截?cái)嘀凶雍鸵恍щ娏W訒?huì)使透射的次級(jí)質(zhì)子個(gè)數(shù)有所減少,當(dāng)不需要考慮次級(jí)質(zhì)子輸運(yùn)過(guò)程時(shí)可以采用。質(zhì)子輸運(yùn)模擬程序并行化后的加速比隨CPU使用個(gè)數(shù)的增加呈線性增長(zhǎng),與文獻(xiàn)結(jié)果一致。與單個(gè)CPU不采用粒子截?cái)嗟那闆r相比,只記錄初級(jí)質(zhì)子的16個(gè)CPU并性程序加速比可達(dá)610倍。

    [1]張?zhí)炀?,樊明武,李振?guó),等.強(qiáng)流回旋加速器技術(shù)的研究進(jìn)展[J].科學(xué)通報(bào),2009,54(15):2169-2176.

    [2]楊朝霞,李德明,張滿洲.基于Geant4模擬的質(zhì)子治療束配系統(tǒng)的束流光學(xué)設(shè)計(jì)[J].核技術(shù),2013,36(7):070201.

    [3]許淑艷.蒙特卡羅方法在實(shí)驗(yàn)核物理中的應(yīng)用[M].北京:原子能出版社,2006.

    [4]張慧,王振超,關(guān)世榮,等.質(zhì)子能量歧離的Geant4模擬[J].黑龍江科學(xué),2014,1(5):14-16.

    [5]路偉.Geant4的開(kāi)發(fā)及其在空間輻射效應(yīng)分析中的應(yīng)用[D].長(zhǎng)沙:國(guó)防科學(xué)技術(shù)大學(xué)研究生院,2007:13-15.

    [6]聶鵬煊,汪一夫,鄭濤,等.CsI(TI)對(duì)高能質(zhì)子能量響應(yīng)的蒙特卡羅研究[J].武漢科技學(xué)院學(xué)報(bào).2010,23(2):27-32.

    [7]王力,李煜,岳騫,等.基于Geant4的射線探測(cè)仿真并行化研究[J].核電子學(xué)與探測(cè)技術(shù),2011,11(11):1237-1240.

    [8]謝紅剛,朱金輝,黃流興.Geant4在高性能計(jì)算機(jī)上的并行[C]//第十屆全國(guó)蒙特卡羅方法及其應(yīng)用學(xué)術(shù)會(huì)論文:148-151.

    [9]張慧,梅雪松,關(guān)世榮,等.4He離子盧瑟福背散射的Geant4模擬[J].黑龍江科學(xué),2013,1(4):22-24.

    [10]Geant4 user's guide for application developers[EB10L].(292-11-30)[2015-03-06]http://geant4.web.cern.ch/geant4/UserDocumentation/ UserGuides/ForApplication Developer/html/index.html.

    [11]ZIOCK H J,ADAMS K J,ALRICK K R,et al.The Proton Radiography Concept[R].Los Alamos National Report,LA-UR-98-1368,1998:1-17.

    Research on Accelerating Method of Proton Transportation Simulation Program Based on Geant4

    YAOZhi-ming,HEI Dong-wei,SONGGu-zhou,MAJi-ming,HANChang-cai,DUANBao-jun,YUE Zhi-qin
    (State Key Laboratory of Intense Pulsed Radiation Simulation and Effect,Northwest Institute ofNuclear Technology,Xi'an 710024,China)

    Proton has bright future in the field of medical,material analysis and radiography.Numerical simulation provides a good way to analyze the process protons transmitting through materials.This article establishes a Monte Carlo program based on Geant4.The accelerating method is introduced.The result implies that secondary particles cut and program paralleled can shorten the time cost effectively.The program can be easily extended which provides a good foundation for fast simulation ofproton transportation.

    Proton transportation;Monte Carlo;Acceleration method;Secondaryparticle cut;Parallel program

    TL99

    A

    1674-8646(2015)07-0007-03

    2015-04-09

    超快物理過(guò)程射線診斷及應(yīng)用研究(51311060502)

    姚志明(1989-),男,吉林輝南人,碩士研究生,助理工程師,主要從事質(zhì)子掃描成像的原理探究。

    猜你喜歡
    質(zhì)子個(gè)數(shù)粒子
    怎樣數(shù)出小正方體的個(gè)數(shù)
    質(zhì)子束放療在腫瘤中的研究新進(jìn)展
    等腰三角形個(gè)數(shù)探索
    怎樣數(shù)出小木塊的個(gè)數(shù)
    基于粒子群優(yōu)化的橋式起重機(jī)模糊PID控制
    怎樣數(shù)出小正方體的個(gè)數(shù)
    基于粒子群優(yōu)化極點(diǎn)配置的空燃比輸出反饋控制
    淺談質(zhì)子守恒
    “質(zhì)子”號(hào)一箭發(fā)雙星
    太空探索(2014年6期)2014-07-10 13:06:11
    基于Matlab的α粒子的散射實(shí)驗(yàn)?zāi)M
    物理與工程(2014年4期)2014-02-27 11:23:08
    国产国语露脸激情在线看| 亚洲中文字幕日韩| 中文资源天堂在线| 啪啪无遮挡十八禁网站| 亚洲精品国产精品久久久不卡| 日本精品一区二区三区蜜桃| 亚洲av成人av| 91成年电影在线观看| a级毛片a级免费在线| or卡值多少钱| 91av网站免费观看| 一级片免费观看大全| 国产精品久久视频播放| 欧美成人一区二区免费高清观看 | 99热只有精品国产| 一边摸一边抽搐一进一小说| 丰满人妻熟妇乱又伦精品不卡| 熟妇人妻久久中文字幕3abv| 国产野战对白在线观看| 在线天堂中文资源库| www.www免费av| 色av中文字幕| 成人国产综合亚洲| 在线观看一区二区三区| 人成视频在线观看免费观看| 99热这里只有精品一区 | 99热6这里只有精品| 日韩大尺度精品在线看网址| 亚洲成人久久性| 香蕉av资源在线| 免费人成视频x8x8入口观看| 国产激情久久老熟女| 女人爽到高潮嗷嗷叫在线视频| 老熟妇仑乱视频hdxx| 男女床上黄色一级片免费看| av福利片在线| 亚洲av片天天在线观看| 男人舔奶头视频| 啦啦啦 在线观看视频| 国产亚洲av嫩草精品影院| 久久伊人香网站| 99国产精品一区二区三区| 午夜免费成人在线视频| 一级片免费观看大全| 欧美黄色淫秽网站| 日本熟妇午夜| av欧美777| 一夜夜www| 欧美色欧美亚洲另类二区| 丁香六月欧美| 欧美又色又爽又黄视频| 叶爱在线成人免费视频播放| 可以免费在线观看a视频的电影网站| 色播在线永久视频| 色在线成人网| 久久精品国产亚洲av高清一级| 两人在一起打扑克的视频| 久久 成人 亚洲| 黄色 视频免费看| 久久精品影院6| 成人手机av| 91大片在线观看| 午夜激情av网站| 久久狼人影院| 亚洲精品久久国产高清桃花| 精品一区二区三区av网在线观看| 欧美成人午夜精品| 亚洲熟妇中文字幕五十中出| 亚洲成人久久爱视频| 黄色毛片三级朝国网站| 女人高潮潮喷娇喘18禁视频| 日日干狠狠操夜夜爽| 老司机午夜福利在线观看视频| 9191精品国产免费久久| 97人妻精品一区二区三区麻豆 | 欧美最黄视频在线播放免费| 此物有八面人人有两片| 美女高潮到喷水免费观看| 在线观看日韩欧美| 成人欧美大片| 免费看a级黄色片| 听说在线观看完整版免费高清| 中文字幕精品亚洲无线码一区 | 日本精品一区二区三区蜜桃| 两性夫妻黄色片| 女警被强在线播放| 色综合婷婷激情| 777久久人妻少妇嫩草av网站| 亚洲av成人不卡在线观看播放网| 日韩视频一区二区在线观看| 色综合欧美亚洲国产小说| 757午夜福利合集在线观看| 亚洲国产精品sss在线观看| 国产亚洲精品久久久久5区| 亚洲免费av在线视频| 一本综合久久免费| 亚洲精品久久成人aⅴ小说| 麻豆av在线久日| 久久热在线av| 丝袜在线中文字幕| 欧美性猛交╳xxx乱大交人| 性欧美人与动物交配| 搞女人的毛片| 色综合亚洲欧美另类图片| a在线观看视频网站| 欧美大码av| 999精品在线视频| 哪里可以看免费的av片| 欧美丝袜亚洲另类 | 一级毛片女人18水好多| 自线自在国产av| 国产一区二区三区在线臀色熟女| 欧美中文日本在线观看视频| 国语自产精品视频在线第100页| 国产精品香港三级国产av潘金莲| 精品久久久久久久末码| 免费在线观看日本一区| 日本五十路高清| 亚洲九九香蕉| 日本在线视频免费播放| 欧美一区二区精品小视频在线| 嫩草影院精品99| 国产区一区二久久| 又黄又爽又免费观看的视频| 亚洲精品久久成人aⅴ小说| 免费高清在线观看日韩| 听说在线观看完整版免费高清| 亚洲五月天丁香| 99久久国产精品久久久| 国产亚洲欧美精品永久| 欧美一级a爱片免费观看看 | 夜夜躁狠狠躁天天躁| 9191精品国产免费久久| 麻豆成人av在线观看| 嫩草影视91久久| 啦啦啦免费观看视频1| 欧美午夜高清在线| 亚洲精品国产精品久久久不卡| 亚洲一码二码三码区别大吗| 最新美女视频免费是黄的| 午夜激情av网站| 1024香蕉在线观看| 黄色片一级片一级黄色片| 91在线观看av| 人人妻人人澡欧美一区二区| 69av精品久久久久久| 成人亚洲精品av一区二区| 一本大道久久a久久精品| 天天躁狠狠躁夜夜躁狠狠躁| www.www免费av| 制服人妻中文乱码| 欧美一区二区精品小视频在线| 亚洲av日韩精品久久久久久密| 人人妻人人澡人人看| 一级毛片高清免费大全| 精品国产一区二区三区四区第35| 一本一本综合久久| 午夜日韩欧美国产| 欧美黑人精品巨大| 久久九九热精品免费| 人人澡人人妻人| 午夜福利高清视频| 精品久久久久久久毛片微露脸| 精品一区二区三区四区五区乱码| 白带黄色成豆腐渣| 欧美午夜高清在线| bbb黄色大片| 9191精品国产免费久久| 一级a爱片免费观看的视频| 最近最新中文字幕大全免费视频| 成人18禁在线播放| 伊人久久大香线蕉亚洲五| 午夜福利免费观看在线| 国产97色在线日韩免费| 久久久久国产精品人妻aⅴ院| 宅男免费午夜| 国产熟女午夜一区二区三区| 久久精品国产综合久久久| 波多野结衣av一区二区av| 成人午夜高清在线视频 | 国产av一区二区精品久久| 成人手机av| 精品高清国产在线一区| 日韩精品免费视频一区二区三区| 久久久久国内视频| 真人做人爱边吃奶动态| 国产免费av片在线观看野外av| 男人操女人黄网站| 成年版毛片免费区| 精品久久久久久久毛片微露脸| 老司机午夜福利在线观看视频| 一边摸一边做爽爽视频免费| 精品国产乱码久久久久久男人| 中文字幕最新亚洲高清| 国内少妇人妻偷人精品xxx网站 | 成人午夜高清在线视频 | 国内毛片毛片毛片毛片毛片| 中文字幕最新亚洲高清| 国产真人三级小视频在线观看| 桃红色精品国产亚洲av| 色婷婷久久久亚洲欧美| 亚洲中文av在线| 99久久综合精品五月天人人| 国产高清视频在线播放一区| 人成视频在线观看免费观看| 免费在线观看日本一区| 日本 av在线| 午夜免费激情av| 色综合婷婷激情| 在线永久观看黄色视频| 国产成人欧美在线观看| 午夜两性在线视频| 免费搜索国产男女视频| 成人永久免费在线观看视频| 久久久久久国产a免费观看| 精品日产1卡2卡| 国内少妇人妻偷人精品xxx网站 | 叶爱在线成人免费视频播放| 欧美久久黑人一区二区| 国产在线观看jvid| 日韩一卡2卡3卡4卡2021年| 一级黄色大片毛片| 免费观看人在逋| 黑人欧美特级aaaaaa片| 精品久久久久久,| 久久久久九九精品影院| 一本大道久久a久久精品| 亚洲国产欧洲综合997久久, | 欧美激情 高清一区二区三区| 可以在线观看的亚洲视频| 精品福利观看| 一夜夜www| 日日爽夜夜爽网站| 国产成年人精品一区二区| 精品欧美一区二区三区在线| 日韩欧美 国产精品| 满18在线观看网站| 99精品在免费线老司机午夜| 国产又色又爽无遮挡免费看| 亚洲一区高清亚洲精品| 大型黄色视频在线免费观看| 国产精华一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 国产精品自产拍在线观看55亚洲| 人人澡人人妻人| www日本黄色视频网| 成年人黄色毛片网站| 成人精品一区二区免费| 色综合亚洲欧美另类图片| 久久婷婷人人爽人人干人人爱| 黄色视频不卡| 日本一本二区三区精品| 久久精品国产亚洲av香蕉五月| 在线看三级毛片| 欧美三级亚洲精品| 亚洲精华国产精华精| 美国免费a级毛片| 长腿黑丝高跟| 成人18禁在线播放| 一a级毛片在线观看| 真人一进一出gif抽搐免费| 日日夜夜操网爽| 亚洲片人在线观看| 首页视频小说图片口味搜索| 人人妻,人人澡人人爽秒播| 又大又爽又粗| 日韩精品中文字幕看吧| 日韩欧美免费精品| 91九色精品人成在线观看| 美女午夜性视频免费| 男人操女人黄网站| 熟女电影av网| 51午夜福利影视在线观看| 国产激情偷乱视频一区二区| 免费在线观看完整版高清| 久久精品国产99精品国产亚洲性色| 亚洲中文字幕一区二区三区有码在线看 | 成人精品一区二区免费| 国产成人系列免费观看| 香蕉丝袜av| 久久久久国产一级毛片高清牌| 99国产精品99久久久久| 国产成人系列免费观看| 最近最新免费中文字幕在线| 黄色a级毛片大全视频| 免费在线观看影片大全网站| 午夜福利免费观看在线| www.精华液| 久久久久久免费高清国产稀缺| 女生性感内裤真人,穿戴方法视频| 在线看三级毛片| 变态另类成人亚洲欧美熟女| 亚洲精品国产精品久久久不卡| 少妇的丰满在线观看| 欧美一区二区精品小视频在线| 欧美激情极品国产一区二区三区| 99热6这里只有精品| 国产精品亚洲美女久久久| 日韩 欧美 亚洲 中文字幕| 国产成人影院久久av| 我的亚洲天堂| 18禁观看日本| 日本免费a在线| 91大片在线观看| 日韩欧美国产一区二区入口| 亚洲av第一区精品v没综合| 亚洲第一电影网av| 在线观看免费视频日本深夜| 久久久水蜜桃国产精品网| 国产成人一区二区三区免费视频网站| 午夜日韩欧美国产| 国产免费av片在线观看野外av| 欧美成狂野欧美在线观看| 亚洲av中文字字幕乱码综合 | 给我免费播放毛片高清在线观看| 99国产精品99久久久久| 国产真人三级小视频在线观看| 久久天躁狠狠躁夜夜2o2o| av片东京热男人的天堂| 好看av亚洲va欧美ⅴa在| 波多野结衣高清无吗| 69av精品久久久久久| 天堂影院成人在线观看| 久久精品人妻少妇| 国内毛片毛片毛片毛片毛片| 校园春色视频在线观看| 免费在线观看影片大全网站| 国产黄色小视频在线观看| 国产激情久久老熟女| 欧美亚洲日本最大视频资源| 琪琪午夜伦伦电影理论片6080| 亚洲成国产人片在线观看| 在线观看66精品国产| 十八禁网站免费在线| 亚洲欧美激情综合另类| 搡老熟女国产l中国老女人| 国产一区在线观看成人免费| 一边摸一边抽搐一进一小说| 免费观看精品视频网站| 99国产极品粉嫩在线观看| 人成视频在线观看免费观看| 亚洲精品av麻豆狂野| 国产又色又爽无遮挡免费看| 国产欧美日韩一区二区精品| 又黄又粗又硬又大视频| 午夜福利一区二区在线看| 好男人在线观看高清免费视频 | 亚洲精品美女久久av网站| 成年女人毛片免费观看观看9| 欧美精品亚洲一区二区| 国产亚洲欧美在线一区二区| 国内精品久久久久久久电影| 女性被躁到高潮视频| 国产精品,欧美在线| 午夜老司机福利片| 国产黄片美女视频| 亚洲人成77777在线视频| 国产久久久一区二区三区| 97碰自拍视频| 免费搜索国产男女视频| 成人欧美大片| av欧美777| 中文字幕人妻熟女乱码| 久久久久免费精品人妻一区二区 | 中文字幕另类日韩欧美亚洲嫩草| 亚洲一区高清亚洲精品| 成人特级黄色片久久久久久久| 国内毛片毛片毛片毛片毛片| xxx96com| 97人妻精品一区二区三区麻豆 | 日日爽夜夜爽网站| 欧美在线黄色| tocl精华| 欧美精品啪啪一区二区三区| 变态另类丝袜制服| 精品久久久久久成人av| 国产爱豆传媒在线观看 | 国产极品粉嫩免费观看在线| 美女国产高潮福利片在线看| 亚洲欧美一区二区三区黑人| 国语自产精品视频在线第100页| 啪啪无遮挡十八禁网站| 国产亚洲av高清不卡| 午夜福利在线观看吧| av免费在线观看网站| 啪啪无遮挡十八禁网站| 亚洲熟女毛片儿| 18禁裸乳无遮挡免费网站照片 | 国产成人精品无人区| 亚洲午夜精品一区,二区,三区| 啪啪无遮挡十八禁网站| 亚洲三区欧美一区| 99精品久久久久人妻精品| 成人av一区二区三区在线看| 一进一出抽搐动态| 午夜福利18| 婷婷丁香在线五月| 日韩精品免费视频一区二区三区| 中文在线观看免费www的网站 | 国产亚洲av嫩草精品影院| 亚洲精品中文字幕在线视频| 久久久精品国产亚洲av高清涩受| 男女午夜视频在线观看| 亚洲avbb在线观看| 免费电影在线观看免费观看| 怎么达到女性高潮| 免费高清在线观看日韩| www日本在线高清视频| 欧美zozozo另类| 精品福利观看| 亚洲成av人片免费观看| 欧美+亚洲+日韩+国产| 97超级碰碰碰精品色视频在线观看| 色在线成人网| 久久伊人香网站| 国产精品香港三级国产av潘金莲| e午夜精品久久久久久久| 亚洲成人免费电影在线观看| 丝袜美腿诱惑在线| √禁漫天堂资源中文www| 嫩草影视91久久| 禁无遮挡网站| 搞女人的毛片| 久久香蕉国产精品| www.熟女人妻精品国产| 美女大奶头视频| 久久久精品国产亚洲av高清涩受| 精品国产超薄肉色丝袜足j| 日日干狠狠操夜夜爽| 久久精品国产综合久久久| 波多野结衣高清作品| 日本三级黄在线观看| 色在线成人网| 久久久精品国产亚洲av高清涩受| 国产一级毛片七仙女欲春2 | 国产精品日韩av在线免费观看| 99精品久久久久人妻精品| 免费无遮挡裸体视频| 午夜福利欧美成人| 色精品久久人妻99蜜桃| 国产精品亚洲美女久久久| 国产精品一区二区三区四区久久 | 欧美午夜高清在线| 亚洲国产欧美一区二区综合| 国产精品99久久99久久久不卡| 久久久久久九九精品二区国产 | 最近最新中文字幕大全免费视频| 成人18禁高潮啪啪吃奶动态图| 精品国产一区二区三区四区第35| 制服人妻中文乱码| av有码第一页| 免费在线观看视频国产中文字幕亚洲| 久久精品成人免费网站| 国内少妇人妻偷人精品xxx网站 | 真人一进一出gif抽搐免费| 一级a爱视频在线免费观看| 最近最新中文字幕大全免费视频| 欧美最黄视频在线播放免费| 99re在线观看精品视频| 非洲黑人性xxxx精品又粗又长| 色哟哟哟哟哟哟| 一级片免费观看大全| 国产免费男女视频| 成在线人永久免费视频| 精品久久久久久久末码| www日本在线高清视频| 欧美日韩中文字幕国产精品一区二区三区| 久久午夜综合久久蜜桃| 亚洲第一av免费看| aaaaa片日本免费| 每晚都被弄得嗷嗷叫到高潮| 国产亚洲av嫩草精品影院| 国产精品1区2区在线观看.| 91av网站免费观看| 亚洲国产日韩欧美精品在线观看 | 黄片大片在线免费观看| 国产精品1区2区在线观看.| av视频在线观看入口| 91九色精品人成在线观看| 女生性感内裤真人,穿戴方法视频| 国内精品久久久久精免费| 亚洲激情在线av| 法律面前人人平等表现在哪些方面| 在线视频色国产色| 老司机午夜福利在线观看视频| 欧美av亚洲av综合av国产av| 在线av久久热| 日韩精品青青久久久久久| 免费在线观看黄色视频的| 老熟妇仑乱视频hdxx| 99久久无色码亚洲精品果冻| 日本撒尿小便嘘嘘汇集6| 国产成人欧美| 久久久久九九精品影院| 夜夜看夜夜爽夜夜摸| 夜夜爽天天搞| 曰老女人黄片| av超薄肉色丝袜交足视频| 久久国产精品影院| 国内毛片毛片毛片毛片毛片| 精品高清国产在线一区| 日韩视频一区二区在线观看| 亚洲av成人一区二区三| 欧美激情高清一区二区三区| 久久精品亚洲精品国产色婷小说| 90打野战视频偷拍视频| 一区二区三区国产精品乱码| 午夜福利免费观看在线| 日日摸夜夜添夜夜添小说| 亚洲精品国产精品久久久不卡| 日日爽夜夜爽网站| or卡值多少钱| 欧美乱码精品一区二区三区| 免费高清在线观看日韩| 一本大道久久a久久精品| 欧美日韩乱码在线| 久久精品亚洲精品国产色婷小说| 亚洲欧美一区二区三区黑人| 欧美精品啪啪一区二区三区| 99久久无色码亚洲精品果冻| 激情在线观看视频在线高清| 亚洲七黄色美女视频| 国语自产精品视频在线第100页| 日韩av在线大香蕉| 国产精品乱码一区二三区的特点| 亚洲七黄色美女视频| 亚洲国产精品合色在线| 国产精品 国内视频| 一二三四在线观看免费中文在| 国产精品久久视频播放| 亚洲国产高清在线一区二区三 | 精品久久蜜臀av无| 欧美av亚洲av综合av国产av| 人妻久久中文字幕网| 国产又爽黄色视频| 亚洲国产中文字幕在线视频| 国产精品香港三级国产av潘金莲| www.熟女人妻精品国产| 91麻豆精品激情在线观看国产| 美女扒开内裤让男人捅视频| 巨乳人妻的诱惑在线观看| 成年女人毛片免费观看观看9| 国产在线观看jvid| 成人18禁在线播放| 黄片小视频在线播放| 男女那种视频在线观看| 深夜精品福利| 怎么达到女性高潮| 色综合站精品国产| 国产亚洲精品一区二区www| 精品免费久久久久久久清纯| 两性午夜刺激爽爽歪歪视频在线观看 | 美女大奶头视频| 91麻豆av在线| 丁香欧美五月| 欧洲精品卡2卡3卡4卡5卡区| 久久婷婷人人爽人人干人人爱| 欧美大码av| 一进一出好大好爽视频| 国产精品一区二区免费欧美| 19禁男女啪啪无遮挡网站| 热re99久久国产66热| 俄罗斯特黄特色一大片| 亚洲最大成人中文| 亚洲国产欧美一区二区综合| 亚洲人成网站在线播放欧美日韩| 禁无遮挡网站| 成人国产一区最新在线观看| 久久久精品国产亚洲av高清涩受| 国产成人系列免费观看| 色综合欧美亚洲国产小说| 一进一出好大好爽视频| 91老司机精品| 国产精品爽爽va在线观看网站 | 人人妻,人人澡人人爽秒播| 久久久国产成人精品二区| 美女免费视频网站| 十分钟在线观看高清视频www| 午夜福利欧美成人| 亚洲男人的天堂狠狠| 国内揄拍国产精品人妻在线 | 夜夜看夜夜爽夜夜摸| 在线观看日韩欧美| 亚洲精品在线观看二区| 人妻久久中文字幕网| 免费av毛片视频| 淫妇啪啪啪对白视频| 19禁男女啪啪无遮挡网站| 欧美日韩亚洲国产一区二区在线观看| 国产亚洲精品第一综合不卡| 国产av又大| 无限看片的www在线观看| 免费女性裸体啪啪无遮挡网站| 天天一区二区日本电影三级| 免费高清在线观看日韩| 香蕉久久夜色| 欧美日韩乱码在线| 搞女人的毛片| 久久久久久久精品吃奶| 在线观看66精品国产| 亚洲国产欧美一区二区综合| 老司机在亚洲福利影院| 亚洲国产精品成人综合色| 欧美乱码精品一区二区三区| 香蕉丝袜av| 欧美精品亚洲一区二区| 日韩高清综合在线| www日本在线高清视频| 亚洲av第一区精品v没综合| 女性被躁到高潮视频| 99国产精品一区二区蜜桃av| 亚洲欧美精品综合一区二区三区| 这个男人来自地球电影免费观看| 欧美一区二区精品小视频在线| 欧美色欧美亚洲另类二区|