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

    球形燃料元件累積旋轉(zhuǎn)角度和角速度問題研究

    2016-04-19 05:48:29朱興望曲世祥中國科學(xué)院上海應(yīng)用物理研究所嘉定園區(qū)上海20800中國科學(xué)院大學(xué)北京00049
    核技術(shù) 2016年3期
    關(guān)鍵詞:角速度

    趙 穎 朱興望 曲世祥 陳 堃(中國科學(xué)院上海應(yīng)用物理研究所 嘉定園區(qū) 上海 20800)2(中國科學(xué)院大學(xué) 北京 00049)

    ?

    球形燃料元件累積旋轉(zhuǎn)角度和角速度問題研究

    趙 穎1,2朱興望1曲世祥1陳 堃1
    1(中國科學(xué)院上海應(yīng)用物理研究所 嘉定園區(qū)上海 201800)2(中國科學(xué)院大學(xué)北京 100049)

    摘要球床堆采用球形燃料元件,其在循環(huán)過程中存在滾動(dòng)磨損以及滑動(dòng)磨損并產(chǎn)生石墨粉塵。本文借助一種球形無線探測系統(tǒng)量化了角速度及累計(jì)旋轉(zhuǎn)角度,并量化了滾動(dòng)磨損對(duì)石墨粉塵產(chǎn)生的影響。實(shí)驗(yàn)結(jié)果表明,累計(jì)旋轉(zhuǎn)角度和角速度與初始徑向位置有關(guān),隨著初始徑向位置由里向外,呈V字型分布,與倉筒邊壁接觸的球累計(jì)旋轉(zhuǎn)角度和角速度遠(yuǎn)大于其他位置,邊壁滾動(dòng)效應(yīng)顯著。計(jì)算結(jié)果表明,滾動(dòng)磨損對(duì)石墨粉塵產(chǎn)生量的貢獻(xiàn)在11%-27%,且靠近倉筒邊壁處該百分比高于其他位置,邊壁效應(yīng)明顯。一般情況下,滾動(dòng)磨損產(chǎn)生石墨粉塵量要低于滑動(dòng)磨損產(chǎn)生的石墨粉塵量,以往在估算球床堆中的石墨粉塵量時(shí),只考慮滑動(dòng)磨損而忽略滾動(dòng)磨損,會(huì)導(dǎo)致石墨粉塵的估算偏高。

    關(guān)鍵詞球形燃料元件,累計(jì)旋轉(zhuǎn)角度,角速度,滾動(dòng)磨損,球形無線探測系統(tǒng)

    中國科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)項(xiàng)目(No.XDA02050100)資助

    第一作者:趙穎,女,1989年出生,2011年畢業(yè)于山東建筑大學(xué),現(xiàn)為博士研究生,研究領(lǐng)域?yàn)槎叹嚯x無線通信、微電子及球床堆堆芯顆粒流

    Supported by Strategic Priority Research Program of Chinese Academy of Sciences(No.XDA02050100)

    First author:ZHAO Ying,female,born in 1989,graduated from Shandong Architecture University in 2011,doctor student,focusing on short-range wireless communication technique,microelectronics,granular dynamics in pebble bed reactor cores

    Rotation angles and angular velocities study of pebble-shaped fuel element based on a detection system

    ZHAO Ying1,2ZHU Xingwang1QU Shixiang1CHEN Kun1
    1(Shanghai Institute of Applied Physics,Chinese Academy of Sciences,Jiading Campus,Shanghai 201800,China)2(University of Chinese Academy of Sciences,Beijing 100049,China)

    AbstractBackground:Pebble Bed Reactor(PBR),most often used in Generation IV high-temperature gas-cooled nuclear reactor,employs pebble-shaped fuel elements.As pebbles flow through the reactor,they will slide and rotate against each other and the wall,causing generation of graphite dust.Purpose:This study aims to investigate the rotation angles and angular velocities of pebble-shaped fuel elements(PSFE).Methods:The rotation angle and angular velocity were studied by using a simulated transparent cylinder to replace the reactor core.A pebble-shaped detection system is employed to test the angular velocities and accumulated rotation angles of PSFEs with wireless interface.The generation of graphite particles by rotational abrasion was also studied.Results:The experimental results show that the initial seeding position has influence on rotation angles and angular velocities.The distribution of rotation angles and angular velocities from the center to the wall can be described as an inverted V-type distribution.The calculation results show that wear mass by rotation accounts for 11%-27%of the total graphite dust.Conclusion:Sliding abrasion is stronger than rotation abrasion in general,graphite generation tends to be overestimated if we only take sliding abrasion into consideration and ignore rotation abrasion.

    Key wordsPebble-shaped fuel element,Rotation angles,Angular velocity,Rotation abrasion,Pebble-shaped detector system

    球床堆是一種采用球形燃料元件的反應(yīng)堆,冷卻劑一般采用氣體或熔鹽[1]。燃料球最外層的包覆材料為石墨,堆芯反射層一般也采用石墨材料[2-6]。球床堆運(yùn)行時(shí),燃料球由于在線裝卸料操作,在堆內(nèi)循環(huán)形成球流。球流穿過堆芯時(shí),燃料球相互之間、燃料球與石墨反射層和其它構(gòu)件之間發(fā)生摩擦、造成磨損并形成石墨粉塵。這些粉塵可能會(huì)積累在堆芯或熱交換系統(tǒng)中,影響反應(yīng)堆的安全和運(yùn)行。例如,AVR(Arbeitsgemeinschaft Versuchreaktor)在20年運(yùn)行期間產(chǎn)生的石墨粉塵50-60 kg[7]。

    針對(duì)石墨粉塵問題,各國學(xué)者開展了多年的模擬和實(shí)驗(yàn)研究。如Lee等[8]采用計(jì)算流體力學(xué)(Computational Fluid Dynamics,CFD)方法模擬燃料球的摩擦;Brendlé等[9]和Luo等[10]分別利用實(shí)驗(yàn)分析了堆芯溫度、壓強(qiáng)以及冷卻劑行為對(duì)石墨的磨損系數(shù)和摩擦系數(shù)的影響。研究表明,堆芯石墨粉塵的產(chǎn)生與燃料球在堆芯的磨損有關(guān)。燃料球在堆芯的磨損既包括滑動(dòng)磨損,又包括滾動(dòng)磨損?;瑒?dòng)磨損產(chǎn)生的石墨粉塵主要由燃料球的移動(dòng)距離和不同溫度對(duì)應(yīng)的磨損率決定。滾動(dòng)磨損則與燃料球在循環(huán)過程中的累計(jì)旋轉(zhuǎn)角度有關(guān)[11-12]。由于球流循環(huán)過程中的燃料元件的累計(jì)旋轉(zhuǎn)角度很難獲得,可供參考的研究數(shù)據(jù)較少,因此在分析石墨粉塵的產(chǎn)生時(shí),一般很難定量分析滾動(dòng)磨損的影響[13-14]。

    為研究燃料球運(yùn)動(dòng)過程中滾動(dòng)摩損對(duì)于石墨粉塵產(chǎn)生的影響,我們借助一種球形無線探測系統(tǒng),量化分析了燃料球的累積旋轉(zhuǎn)角度和角速度。為便于實(shí)驗(yàn),我們采用常溫下的模擬球來模擬實(shí)際堆芯中的球流。研究表明,經(jīng)過合理?;哪M球系統(tǒng)能夠較好反映實(shí)際堆芯中球流的運(yùn)動(dòng)學(xué)特性[15]。

    1 實(shí)驗(yàn)簡介

    1.1實(shí)驗(yàn)裝置

    實(shí)驗(yàn)研究對(duì)象選取的是麻省理工學(xué)院(Massachusetts Institute of Technology,MIT)設(shè)計(jì)的球床模塊高溫氣冷堆(Pebble-Bed Modular Reactor,PBMR)[16],這是一種有代表性的球床堆。實(shí)驗(yàn)的主要研究內(nèi)容是燃料球在卸料過程中的角速度和累計(jì)旋轉(zhuǎn)角度。通過?;治?,確定實(shí)驗(yàn)裝置的幾何尺寸依據(jù)PBMR的1/6 縮放[17]。PBMR采用的球形燃料元件直徑d1為60 mm,MIT設(shè)計(jì)的反應(yīng)堆半徑為1.75 m(即29d1),依據(jù)Rycroft等[17]的分析,經(jīng)1/6 縮放半徑四舍五入為5d1。PBMR出口半徑為5d1,但考慮到當(dāng)出口半徑小于2.5d1時(shí)就會(huì)出現(xiàn)卡球現(xiàn)象,故1/6縮放后出口半徑設(shè)定為2.5d1[18]。模擬球的直徑d為30 mm,實(shí)驗(yàn)裝置尺寸與PBMR幾何尺寸的對(duì)比如表1所示。實(shí)驗(yàn)裝置主容器為圓柱形倉筒結(jié)構(gòu),倉筒的直徑為10d。倉筒的底端連接一個(gè)坡度為30°的漏斗,漏斗的底端連接一個(gè)半徑為2.5d的圓柱形出口管,出口管的底端為一可閉合的擋板。

    表1 實(shí)驗(yàn)裝置尺寸與PBMR尺寸對(duì)比Table 1 Size of experiment apparatus vs.size of PBMR.

    實(shí)驗(yàn)倉筒等采用透明亞克力材質(zhì),高速攝像可記錄表層模擬球運(yùn)動(dòng)信息。模擬球采用聚丙烯材質(zhì),以便于模擬燃料球的摩擦系數(shù)[15]。實(shí)驗(yàn)共使用2 050個(gè)模擬球。

    1.2球形無線探測系統(tǒng)

    為實(shí)時(shí)記錄模擬球運(yùn)動(dòng)過程中的角速度和累積旋轉(zhuǎn)角度,我們采用了自主研發(fā)的球形無線探測系統(tǒng),如圖1所示。該系統(tǒng)主要包括探測球、接收端和主控機(jī)等。其中探測球的尺寸、密度及最外層材質(zhì)與模擬球完全一致。探測球內(nèi)置運(yùn)動(dòng)傳感器、溫度傳感器、無線收發(fā)裝置以及電池等。其中運(yùn)動(dòng)傳感器位于探測球球心,且通過配重方法,使球重心與幾何中心重合。探測球能夠?qū)崟r(shí)測量和無線傳輸實(shí)驗(yàn)過程中采集的加速度、角速度和溫度等信息,并給出實(shí)時(shí)姿態(tài)和累計(jì)旋轉(zhuǎn)角度。加速度和溫度的測量精度為相對(duì)誤差在3%以內(nèi),靜止時(shí)角速度的絕對(duì)誤差在±2°·s-1以內(nèi),運(yùn)動(dòng)狀態(tài)下角速度測量相對(duì)誤差在3%以內(nèi)。

    圖1 球形無線探測系統(tǒng)Fig.1 Pebble-shaped detection system.

    1.3實(shí)驗(yàn)過程

    實(shí)驗(yàn)時(shí)首先閉合倉筒底部擋板,令模擬球自倉筒頂部自由下落,形成高度約88 cm的不規(guī)則堆積球床,以模擬實(shí)際堆芯球床的堆積。在球床頂端放置探測球,探測球球心初始位置分為4個(gè)不同的徑向位置R,其半徑由內(nèi)向外依次為:0 cm、4 cm、8cm、13.5 cm,對(duì)應(yīng)的探測球個(gè)數(shù)依次為1、4、8、18顆。球心徑向位置位于13.5 cm處的探測球與倉筒內(nèi)壁直接接觸。實(shí)驗(yàn)時(shí)打開底部擋板,模擬球和探測球在重力作用下自由下落,直至全部流出倉筒,以模擬實(shí)際堆芯中燃料球的運(yùn)動(dòng)。在此過程中探測球以100 Hz的采樣率實(shí)時(shí)測量自身的加速度、角速度,并通過計(jì)算給出實(shí)時(shí)姿態(tài)和累計(jì)旋轉(zhuǎn)角度。

    圖2 圓柱形倉筒卸料實(shí)驗(yàn)Fig.2 Discharging experiment with cylindrical container.

    圖3 探測球初始位置Fig.3 Initial seeding position of detection pebbles.

    2 結(jié)果與分析

    由于探測球在倉筒中的運(yùn)動(dòng)具有隨機(jī)性,本文重復(fù)進(jìn)行了45次上述實(shí)驗(yàn)以獲取探測球物理量統(tǒng)計(jì)規(guī)律。

    2.1探測球滾動(dòng)的量化

    本文采用三維累積旋轉(zhuǎn)角度來量化探測球在下落過程中產(chǎn)生的旋轉(zhuǎn),其平均值的計(jì)算采用如下的方式進(jìn)行。此處以靠近筒壁(R=13.5 cm)的探測球?yàn)槔?,給出探測球累積旋轉(zhuǎn)角度平均值的計(jì)算過程。把符合該位置的全部實(shí)驗(yàn)數(shù)據(jù)(45×18=810組)進(jìn)行統(tǒng)計(jì),并用正態(tài)分布函數(shù)分別對(duì)三維累積旋轉(zhuǎn)角度分布進(jìn)行擬合(圖4),從而得到平均累積旋轉(zhuǎn)角度分別為:1 156.6°(X軸)、709.3°(Y軸)和1 026.2° (Z軸),其它擬合參數(shù)如表2所示。

    表2 擬合參數(shù)Table 2 Parameters of experiment apparatus.

    圖4 X軸(a)、Y軸(b)、Z軸(c)累計(jì)旋轉(zhuǎn)角度統(tǒng)計(jì)直方圖Fig.4 Statistical results of X-axis(a),Y-axis(b),Z-axis(c)rotation angle.

    同樣,位于R=0 cm、4 cm、8 cm的探測球也采用類似方法處理。由于計(jì)算過程相同,此處不再累述。不同位置R的旋轉(zhuǎn)累積角度如圖5所示。結(jié)合相同實(shí)驗(yàn)條件的滯留時(shí)間分布[18],給出不同徑向位置R的平均角速度如圖6所示。

    圖5 累計(jì)旋轉(zhuǎn)角度統(tǒng)計(jì)結(jié)果Fig.5 Statistical results of rotation angle.

    圖6 平均角速度Fig.6 Average angular velocity.

    由統(tǒng)計(jì)結(jié)果可看出,探測球在下落過程中產(chǎn)生的平均累積旋轉(zhuǎn)角度和角速度與初始位置相關(guān),整體分布呈V形,且隨初始位置半徑R而增大;靠近倉筒軸心處滾動(dòng)小而緩,而在與筒壁接觸的位置滾動(dòng)較為突出,約為中心位置探測球的兩倍。產(chǎn)生這一現(xiàn)象,主要是由于本次實(shí)驗(yàn)中R=13.5 cm處的探測球采用靠近筒壁的放置方式,實(shí)驗(yàn)開啟擋板后,探測球易貼近筒壁滾落,其角速度要比其它位置探測球大(圖7)。需要指出的是,目前在利用顆粒離散元方法(Discrete Element Method,DEM)等數(shù)值計(jì)算方法來模擬小球在球床中的運(yùn)動(dòng)規(guī)律時(shí),這一特殊的邊壁滾動(dòng)效應(yīng)通常并沒有被考慮到[18-20]。

    圖7 初始徑向位置為0 cm和13.5 cm處球的角速度Fig.7 Angular velocities for the initial radial location of 0 cm and 13.5 cm.

    2.2石墨粉塵的計(jì)算

    通過參考探測球滾動(dòng)的量化,可以進(jìn)一步建立在考慮滾動(dòng)摩擦情況下石墨粉塵的計(jì)算模型,并對(duì)滑動(dòng)與滾動(dòng)間進(jìn)行量化的比較。

    通常情況下,單個(gè)石墨因滾動(dòng)摩擦產(chǎn)生的粉塵可由式(1)計(jì)算:

    式中:N為Voronoi單元中與小球直接接觸的小球個(gè)數(shù);Prot為單位時(shí)間內(nèi)石墨球間因滾動(dòng)磨損產(chǎn)生的粉塵量(可通過實(shí)驗(yàn)測量),g·s-1;t為石墨球的平均滯留時(shí)間,s。

    而因滑動(dòng)摩擦產(chǎn)生的粉塵量通常由式(2)計(jì)算:

    式中:Pslid為單位距離上石墨粉塵因滑動(dòng)而產(chǎn)生的磨損量,g·m-1;Lslid為石墨球的平均滑動(dòng)距離,m。

    實(shí)驗(yàn)中受到重力作用而自由下落的探測球,其平均運(yùn)動(dòng)距離近似為流道長度H,探測球在下滑過程中發(fā)生“原地打轉(zhuǎn)”的概率較小,因此有:

    根據(jù)式(1)-(3)計(jì)算出下落過程中產(chǎn)生的石墨粉塵總量:

    探測球在下落期間的平均滾動(dòng)距離Lrot由本實(shí)驗(yàn)計(jì)算出:

    式中:θ為平均旋轉(zhuǎn)累積角度;r=1.5 cm為探測球半徑。

    利用理論模型并類比本實(shí)驗(yàn)裝置,在考慮滾動(dòng)和滑動(dòng)因素下對(duì)單個(gè)小球的粉塵產(chǎn)生量進(jìn)行分析。計(jì)算所需參數(shù)歸為表3。

    表3 計(jì)算使用的參數(shù)Table 3 Parameters in calculating.

    對(duì)表3中計(jì)算參數(shù)作簡要說明:

    1)滑動(dòng)磨損率Pslid。根據(jù)Troy等對(duì)石墨滑動(dòng)磨損率的實(shí)驗(yàn)研究[11],在一定范圍內(nèi),滑動(dòng)磨損率隨著角速度以及受力而增加,兩顆型號(hào)為GM-101(由Graphtek LLC制造)的石墨球接觸處的滑動(dòng)磨損率變化范圍在0.005-0.991 g·m-1之間。本實(shí)驗(yàn)結(jié)果表明,探測球的角速度及負(fù)載均小于Troy的實(shí)驗(yàn)條件,因此本文保守取Prot=0.005 g·m-1(實(shí)際上,此處Prot已被高估)。

    2)滾動(dòng)磨損率Prot。根據(jù)Troy等對(duì)石墨球滾動(dòng)磨損的實(shí)驗(yàn)研究[21],滾動(dòng)磨損率隨著角速度以及受力而增加,石墨球在450 r·min-1和1 500 r·min-1轉(zhuǎn)速下每個(gè)觸點(diǎn)的磨損率介于0.003-0.07 g·min-1之間。由于滑動(dòng)磨損率的研究采用的是Graphtek LLC 的 GM-101石墨球[11],為方便比較,本文采用GM-101且負(fù)載為15 kg左右的滾動(dòng)磨損率,對(duì)1500 r·min-1和450 r·min-1轉(zhuǎn)速下的滾動(dòng)磨損率求均值為0.00585-0.00365g·min-1,并外推求得30r·min-1時(shí)磨損率為0.00277g·min-1(4.617×10-5g·s-1)。

    3)N的取值。規(guī)則排布情況下,Voronoi單元中與球相互鄰近的球大于20個(gè),與之直接接觸的球的個(gè)數(shù)為12,如圖8所示。根據(jù)Yang等[22]的研究,在隨機(jī)排布的球床中,Voronoi單元中與球相互鄰近的球?yàn)?4個(gè)左右,故與之直接接觸的球的個(gè)數(shù)小于規(guī)則排布的球的個(gè)數(shù)12,此處保守取12個(gè),對(duì)于靠近倉筒邊壁的探測球,保守6個(gè)。排除重復(fù)計(jì)算的情況后,靠近倉筒邊壁的探測球N=3,其他位置N=6。

    圖8 規(guī)則排布時(shí)與中心球的接觸情況Fig.8 Contact conditions of spheres under rule configuration.

    經(jīng)過計(jì)算,小球的滾動(dòng)磨損量和滑動(dòng)磨損量如圖9所示。從圖9中可看出,石墨粉塵的產(chǎn)生量來自于滾動(dòng)磨損的貢獻(xiàn)在11%-27%,靠近倉筒中心位置處小,而靠近倉筒邊壁處貢獻(xiàn)較大。這意味著在計(jì)算球床堆中燃料球(石墨球)的石墨粉塵時(shí),若忽略滾動(dòng)磨損的貢獻(xiàn),則石墨粉塵的估算值偏大。另外,圖9也說明,邊壁處滾動(dòng)產(chǎn)生的石墨粉塵要比其它位置多,這就要求我們在分析時(shí)需要考慮燃料球(石墨球)的邊壁效應(yīng)。

    圖9 滾動(dòng)磨損質(zhì)量、滑動(dòng)磨損質(zhì)量以及滾動(dòng)磨損質(zhì)量與總量的百分比Fig.9 Rotation wear mass,sliding wear mass and the percentage of rotation wear mass and total wear mass.

    3 結(jié)語

    球床堆采用球形燃料元件,循環(huán)過程中存在球與球以及球與邊壁的滑動(dòng)摩擦和滾動(dòng)摩擦,從而導(dǎo)致石墨粉塵的產(chǎn)生。有關(guān)滑動(dòng)磨損的研究很多,且在估算石墨粉塵產(chǎn)生量時(shí),通常由滑動(dòng)磨損率和滑動(dòng)距離保守計(jì)算得到。本文利用自主開發(fā)的球形無線探測系統(tǒng),在模擬實(shí)驗(yàn)裝置中開展實(shí)驗(yàn),精確測量實(shí)時(shí)角速度,并給出累積旋轉(zhuǎn)角度(亦可求得累計(jì)滾動(dòng)距離),結(jié)合滯留時(shí)間給出平均角速度,從而量化了滾動(dòng)。通過類比文中實(shí)驗(yàn),首次量化了滾動(dòng)產(chǎn)生的石墨粉塵,計(jì)算結(jié)果表明,通常采用的忽略滾動(dòng)的做法會(huì)使得計(jì)算得到的石墨粉塵偏于高估。

    同時(shí),通過本次實(shí)驗(yàn),我們發(fā)現(xiàn)在采用數(shù)值計(jì)算等方法來模擬球床的運(yùn)動(dòng)規(guī)律時(shí),需要考慮球的邊壁效應(yīng)。利用本實(shí)驗(yàn)裝置采集的轉(zhuǎn)動(dòng)角度、角速度等信息可以進(jìn)一步幫助改善模型。

    參考文獻(xiàn)

    1 彭超,朱興望,張國慶,等.采用 SCALE 計(jì)算氟鹽冷卻高溫堆產(chǎn)氚量的一些問題[J].核技術(shù),2015,38(8):080601.DOI:10.11889/j.0253-3219.2015.hjs.38.080601 PENG Chao,ZHU Xingwang,ZHANG Guoqing,et al.Issues in the calculation of the tritium production of the fluoride-salt-cooled high-temperature reactors using SCALE[J].Nuclear Techniques,2015,38(8):080601.DOI:10.11889/j.0253-3219.2015.hjs.38.080601

    2 雒曉衛(wèi),于溯源,張振聲,等.HTR-10 產(chǎn)生石墨粉塵量的估算及其尺寸分布[J].核動(dòng)力工程,2005,26(2):203-208 LUO Xiaowei,YU Suyuan,ZHANG Zhensheng,et al.Estimation of graphite dust quantity and size distribution of graphite particle in HTR-10[J].Nuclear Power Engineering,2005,26(2):203-208

    3 Hong R,Huber S,Lee K,et al.Reactor safety and mechanical design for the annular pebble-bed advanced high temperature reactor[R].UC Berkeley:Department of Nuclear Engineering,Report UCBTH09-001,2009

    4 Bickel J E,Laufer M R,Li L,et al.Conceptual design,experiments,and analysis for the core of an FHR-16 test reactor[C].Proceedings of the 2010 International Congress on Advances in Nuclear Power Plants-ICAPP'10,2010

    5 Wu Z,Lin D,Zhong D.The design features of the HTR-10[J].Nuclear Engineering and Design,2002,218(1):25-32.DOI:10.1016/j.nucengdes.2009.02.023

    6 朱興望,王帥,彭超,等.固態(tài)釷基熔鹽堆中14C的產(chǎn)生及釋放探討[J].核技術(shù),2015,38(3):030603.DOI:10.11889/j.0253-3219.2015.hjs.38.030603 ZHU Xingwang,WANG Shuai,PENG Chao,et al.Production and release of14C in TMSR-SF1[J].Nuclear Techniques,2015,38(3):030603.DOI:10.11889/j.0253-3219.2015.hjs.38.030603

    7 Gottaut H,Krüger K.Results of experiments at the AVR reactor[J].Nuclear Engineering and Design,1990,121(2):143-153.DOI:10.1016/0029-5493(90)90099-J

    8 Lee J J,Park G C,Kim K Y,et al.Numerical treatment of pebble contact in the flow and heat transfer analysis of a pebble bed reactor core[J].Nuclear Engineering and Design,2007,237(22):2183-2196.DOI:10.1016/j.nucengdes.2007.03.046

    9 Brendlé M,Stempflé P.Triboreactions of graphite with moisture - a new model of triboreactor for integrating friction and wear[J].Wear,2003,254(9):818-826.DOI:10.1016/S0043-1648(03)00231-X

    10 Luo X,Yu S,Sheng X,et al.Temperature effect on IG-11 graphite wear performance[J].Nuclear Engineering and Design,2005,235(21):2261-2274.DOI:10.1016/j.nucengdes.2005.05.001

    11 Troy R S,Tompson R V,Ghosh T K,et al.Generation of graphite particles by sliding abrasion and their characterization[J].Nuclear Technology,2015,189(3):241-257.DOI:10.13182/NT14-25

    12 Troy R S,Tompson R V,Ghosh T K,et al.Particle production by rotational abrasion between graphite spheres[J].Nuclear Technology,2015,191(1):71-91.DOI:10.13182/NT14-109

    13 Rostamian M,Johnson G,Hiruta M,et al.Computational and experimental prediction of dust production in pebble bed reactors—Part I[J].Nuclear Engineering and Design,2013,263:500-508.DOI:10.1016/j.nucengdes.2013.04.019

    14 Hiruta M,Johnson G,Rostamian M,et al.Computational and experimental prediction of dust production in pebble bed reactors—Part II[J].Nuclear Engineering and Design,2013,263:509-514.DOI:10.1016/j.nucengdes.2013.04.032

    15 Griveau A.Modeling and transient analysis for the pebble bed advanced high temperature reactor(PB-AHTR)[R].MS Project Report,UCBTH-07-001,2007

    16 Kadak A C.MIT pebble bed reactor project[J].Nuclear Engineering and Technology,2007,39(2):95-102.DOI:10.5516/NET.2007.39.2.095

    17 Rycroft C H,Dehbi A,Lind T,et al.Granular flow in pebble-bed nuclear reactors:scaling,dust generation,and stress[J].Nuclear Engineering and Design,2013,265:69-84.DOI:10.1016/j.nucengdes.2013.07.010

    18 Zuriguel I,Pugnaloni L A,Garcimartin A,et al.Jamming during the discharge of grains from a silo described as a percolating transition[J].Physical Review E,2003,68(3):030301.DOI:10.1103/PhysRevE.68.030301

    19 Yu Y,Saxén H.Discrete element method simulation of properties of a 3D conical hopper with mono-sized spheres[J].Advanced Powder Technology,2011,22(3):324-331.DOI:10.1016/j.apt.2010.04.003

    20 Laufer M R.Granular dynamics in pebble bed reactor cores[D].Berkeley:University of California,2013

    21 Troy R S,Tompson R V,Ghosh T K,et al.Particle production by rotational abrasion between graphite spheres[J].Nuclear Technology,2015,191(1):71-91.DOI:dx.doi.org/10.13182/NT14-109

    22 Yang R Y,Zou R P,Yu A B.Voronoi tessellation of the packing of fine uniform spheres[J].Physical Review E,2002,65(4):041302.DOI:dx.doi.org/10.1103/PhysRevE.65.041302

    收稿日期:2015-11-26,修回日期:2016-01-27

    Corresponding author:CHEN Kun,E-mail:chenkun@sinap.ac.cn

    通信作者:陳堃,E-mail:chenkun@sinap.ac.cn

    DOI:10.11889/j.0253-3219.2016.hjs.39.030603

    中圖分類號(hào)TL37

    猜你喜歡
    角速度
    《液壓與氣動(dòng)》常用單位的規(guī)范
    MEMS偏航角速度傳感器在旋轉(zhuǎn)導(dǎo)彈穩(wěn)定回路中的應(yīng)用
    高中物理角速度矢量性問題的教學(xué)探究
    巧用一張紙分析角速度
    圓周運(yùn)動(dòng)角速度測量方法賞析
    半捷聯(lián)雷達(dá)導(dǎo)引頭視線角速度提取
    基于構(gòu)架點(diǎn)頭角速度的軌道垂向長波不平順在線檢測
    武術(shù)套路旋風(fēng)腳900°接馬步可行性試驗(yàn)研究
    對(duì)一道圓周問題的拓展
    一種無角速度信息的撓性航天器姿態(tài)控制方法
    日韩高清综合在线| 男人舔女人下体高潮全视频| 最近最新中文字幕大全电影3| 久久中文看片网| 日韩在线高清观看一区二区三区 | 黄片wwwwww| 色尼玛亚洲综合影院| 欧美一区二区亚洲| 一进一出抽搐gif免费好疼| 亚洲精品乱码久久久v下载方式| 亚洲av电影不卡..在线观看| or卡值多少钱| 18禁黄网站禁片免费观看直播| 久久6这里有精品| 日韩,欧美,国产一区二区三区 | 搡女人真爽免费视频火全软件 | 国产伦精品一区二区三区视频9| 可以在线观看毛片的网站| 国产亚洲av嫩草精品影院| 91av网一区二区| 精品人妻一区二区三区麻豆 | 亚洲性久久影院| 99精品久久久久人妻精品| 天堂动漫精品| 亚洲人成网站在线播放欧美日韩| 蜜桃久久精品国产亚洲av| 免费人成视频x8x8入口观看| av国产免费在线观看| 日本熟妇午夜| 亚洲av电影不卡..在线观看| 狠狠狠狠99中文字幕| 国产精品无大码| 亚洲精华国产精华液的使用体验 | 色av中文字幕| 观看美女的网站| 国产三级在线视频| 国产精品99久久久久久久久| 免费在线观看成人毛片| 露出奶头的视频| 无遮挡黄片免费观看| 少妇的逼水好多| 三级毛片av免费| 特级一级黄色大片| 能在线免费观看的黄片| 欧美日韩中文字幕国产精品一区二区三区| 久久久久久九九精品二区国产| 最近最新中文字幕大全电影3| 日韩高清综合在线| 亚州av有码| 国产女主播在线喷水免费视频网站 | av在线观看视频网站免费| 精品一区二区三区av网在线观看| 久久久久九九精品影院| 哪里可以看免费的av片| 亚洲成人免费电影在线观看| 精品久久久噜噜| 亚洲va在线va天堂va国产| 黄色女人牲交| 亚洲国产高清在线一区二区三| 国产色婷婷99| 亚洲人与动物交配视频| 欧美性感艳星| 一区福利在线观看| 欧美最新免费一区二区三区| 久久国产精品人妻蜜桃| 国产人妻一区二区三区在| 有码 亚洲区| 欧美+亚洲+日韩+国产| 一个人看视频在线观看www免费| 亚洲av不卡在线观看| 精品久久久久久久末码| 欧美黑人欧美精品刺激| 国产高清视频在线播放一区| 欧美一级a爱片免费观看看| 亚洲四区av| 国产av不卡久久| 2021天堂中文幕一二区在线观| 国内精品一区二区在线观看| 一本精品99久久精品77| 最后的刺客免费高清国语| 波多野结衣高清作品| 日韩一本色道免费dvd| 两性午夜刺激爽爽歪歪视频在线观看| 黄色丝袜av网址大全| 国产人妻一区二区三区在| a级毛片免费高清观看在线播放| 免费看a级黄色片| 成人高潮视频无遮挡免费网站| 老司机午夜福利在线观看视频| 国产中年淑女户外野战色| 超碰av人人做人人爽久久| av国产免费在线观看| 日韩中文字幕欧美一区二区| 午夜日韩欧美国产| 色尼玛亚洲综合影院| 在现免费观看毛片| 超碰av人人做人人爽久久| 我的女老师完整版在线观看| 久久久午夜欧美精品| 在线a可以看的网站| 免费一级毛片在线播放高清视频| 色视频www国产| 桃红色精品国产亚洲av| 两个人视频免费观看高清| 久久99热这里只有精品18| 久久精品91蜜桃| 天堂av国产一区二区熟女人妻| 热99re8久久精品国产| 免费人成视频x8x8入口观看| 久久精品国产亚洲av香蕉五月| 免费在线观看影片大全网站| 18+在线观看网站| 高清在线国产一区| 亚洲人与动物交配视频| 亚洲国产精品成人综合色| 久久精品综合一区二区三区| 美女免费视频网站| 身体一侧抽搐| 久久久久国产精品人妻aⅴ院| 联通29元200g的流量卡| 内地一区二区视频在线| 国产麻豆成人av免费视频| 久久久国产成人免费| 美女高潮喷水抽搐中文字幕| 成人av一区二区三区在线看| 国产毛片a区久久久久| 高清在线国产一区| 真人一进一出gif抽搐免费| 嫩草影院精品99| 黄色女人牲交| 日本与韩国留学比较| 免费观看在线日韩| 国产成人一区二区在线| 真实男女啪啪啪动态图| av女优亚洲男人天堂| 搡老妇女老女人老熟妇| 国产91精品成人一区二区三区| 三级国产精品欧美在线观看| 中文资源天堂在线| 九九爱精品视频在线观看| 久久久久性生活片| 最近最新免费中文字幕在线| 久久久精品欧美日韩精品| 97人妻精品一区二区三区麻豆| 欧美最新免费一区二区三区| 国产精品人妻久久久久久| 中文在线观看免费www的网站| xxxwww97欧美| 国内精品宾馆在线| 两个人视频免费观看高清| 99在线人妻在线中文字幕| 成人国产麻豆网| 婷婷精品国产亚洲av在线| 欧美三级亚洲精品| 久久精品国产亚洲av涩爱 | 亚洲av中文字字幕乱码综合| 女人被狂操c到高潮| 18禁裸乳无遮挡免费网站照片| 精品免费久久久久久久清纯| 国产黄a三级三级三级人| 亚洲综合色惰| 3wmmmm亚洲av在线观看| 91午夜精品亚洲一区二区三区 | 日日撸夜夜添| 精品久久久久久,| 精品一区二区三区视频在线观看免费| 午夜日韩欧美国产| 国产精品一区二区免费欧美| 国产欧美日韩一区二区精品| 校园人妻丝袜中文字幕| 国产欧美日韩精品亚洲av| 日韩欧美三级三区| 国产午夜精品久久久久久一区二区三区 | 色哟哟哟哟哟哟| 亚洲精品亚洲一区二区| av福利片在线观看| 亚洲,欧美,日韩| 真实男女啪啪啪动态图| 99热这里只有精品一区| 91精品国产九色| 啦啦啦啦在线视频资源| 国产真实乱freesex| a在线观看视频网站| 欧美激情在线99| 天美传媒精品一区二区| 欧美日韩中文字幕国产精品一区二区三区| 免费电影在线观看免费观看| 97人妻精品一区二区三区麻豆| 欧美又色又爽又黄视频| 国产精品久久久久久久电影| 小说图片视频综合网站| 国产精品98久久久久久宅男小说| 国产免费男女视频| 国产伦一二天堂av在线观看| 在线观看舔阴道视频| 日日夜夜操网爽| 亚洲人成网站高清观看| 高清在线国产一区| 日本三级黄在线观看| 欧美色欧美亚洲另类二区| 久久久久久久久大av| 婷婷亚洲欧美| 中文字幕精品亚洲无线码一区| 99久久九九国产精品国产免费| 午夜爱爱视频在线播放| 狂野欧美激情性xxxx在线观看| 亚洲精品国产成人久久av| 亚洲熟妇中文字幕五十中出| 国产男人的电影天堂91| 国产一区二区激情短视频| 午夜a级毛片| 男人舔女人下体高潮全视频| 他把我摸到了高潮在线观看| avwww免费| 久久国产乱子免费精品| www日本黄色视频网| 久久国产乱子免费精品| 婷婷色综合大香蕉| 又粗又爽又猛毛片免费看| 精品人妻偷拍中文字幕| 日本与韩国留学比较| 亚洲成人中文字幕在线播放| 此物有八面人人有两片| 亚洲av中文av极速乱 | 我要搜黄色片| 91久久精品国产一区二区成人| 欧美成人性av电影在线观看| 两个人视频免费观看高清| 国产黄色小视频在线观看| 少妇猛男粗大的猛烈进出视频 | 此物有八面人人有两片| 波多野结衣高清无吗| 午夜福利18| 热99在线观看视频| 3wmmmm亚洲av在线观看| 国产蜜桃级精品一区二区三区| 日韩中字成人| 国产一区二区三区视频了| 尾随美女入室| 最好的美女福利视频网| 亚洲经典国产精华液单| 精品国内亚洲2022精品成人| 99热这里只有精品一区| 精品一区二区三区视频在线| 欧美zozozo另类| 99久久中文字幕三级久久日本| 一级a爱片免费观看的视频| 老司机午夜福利在线观看视频| 美女高潮喷水抽搐中文字幕| 一a级毛片在线观看| 我的老师免费观看完整版| av在线亚洲专区| 观看免费一级毛片| 国产麻豆成人av免费视频| 亚洲国产精品sss在线观看| 在线观看66精品国产| 成人一区二区视频在线观看| 性色avwww在线观看| 国产精品久久久久久精品电影| 999久久久精品免费观看国产| 午夜精品一区二区三区免费看| 日本在线视频免费播放| 欧美黑人巨大hd| 九色成人免费人妻av| 午夜福利在线观看免费完整高清在 | 欧美高清成人免费视频www| 中文字幕精品亚洲无线码一区| 国产黄色小视频在线观看| 97超视频在线观看视频| 人人妻人人澡欧美一区二区| 日日摸夜夜添夜夜添av毛片 | 久久国内精品自在自线图片| 国产免费av片在线观看野外av| 色哟哟哟哟哟哟| 一个人看的www免费观看视频| 真人做人爱边吃奶动态| 亚洲不卡免费看| 18禁黄网站禁片午夜丰满| 久久久久久久精品吃奶| 床上黄色一级片| 男女做爰动态图高潮gif福利片| 黄色日韩在线| 亚洲欧美精品综合久久99| 男人的好看免费观看在线视频| 色综合站精品国产| 国产精品一区二区三区四区久久| 日韩欧美精品v在线| 一进一出抽搐gif免费好疼| 亚洲精华国产精华液的使用体验 | 一区二区三区高清视频在线| 日日夜夜操网爽| 色av中文字幕| 成年免费大片在线观看| 亚洲av免费高清在线观看| 一级毛片久久久久久久久女| 哪个播放器可以免费观看大片| 亚洲国产高清在线一区二区三| 国产一区二区三区综合在线观看 | 亚洲欧美精品自产自拍| 人妻 亚洲 视频| 伊人久久国产一区二区| 秋霞在线观看毛片| 精品一区二区三区视频在线| 大片免费播放器 马上看| 国产男人的电影天堂91| 国产伦精品一区二区三区视频9| 午夜福利在线在线| 日日摸夜夜添夜夜添av毛片| 国产日韩欧美在线精品| 免费黄色在线免费观看| 偷拍熟女少妇极品色| 又大又黄又爽视频免费| 国产爱豆传媒在线观看| 成人毛片a级毛片在线播放| 永久网站在线| 精品国产露脸久久av麻豆| 日韩国内少妇激情av| 91狼人影院| 国产在线男女| 99热这里只有是精品50| 中文在线观看免费www的网站| 91精品国产九色| 日韩 亚洲 欧美在线| 内地一区二区视频在线| 国产欧美日韩一区二区三区在线 | 久久精品国产亚洲网站| 亚洲av免费高清在线观看| 日韩免费高清中文字幕av| 毛片女人毛片| 国产老妇伦熟女老妇高清| 男人和女人高潮做爰伦理| 久久青草综合色| 国产亚洲欧美精品永久| 亚洲内射少妇av| 亚洲av成人精品一二三区| 免费播放大片免费观看视频在线观看| 亚洲激情五月婷婷啪啪| 我的老师免费观看完整版| 亚洲精品aⅴ在线观看| 边亲边吃奶的免费视频| 国产精品福利在线免费观看| 高清毛片免费看| 中文字幕亚洲精品专区| 18禁在线无遮挡免费观看视频| 国语对白做爰xxxⅹ性视频网站| 国产精品免费大片| 亚洲精品色激情综合| 色婷婷av一区二区三区视频| 亚洲av福利一区| 狂野欧美激情性xxxx在线观看| 成人特级av手机在线观看| 天堂8中文在线网| 国产成人a∨麻豆精品| 王馨瑶露胸无遮挡在线观看| 女人十人毛片免费观看3o分钟| 狂野欧美激情性xxxx在线观看| 日本猛色少妇xxxxx猛交久久| 久久久成人免费电影| 观看av在线不卡| 午夜福利在线在线| 熟女av电影| 欧美人与善性xxx| 亚洲人成网站在线观看播放| 国产69精品久久久久777片| 热re99久久精品国产66热6| 日本黄色片子视频| 欧美3d第一页| 久久久久久人妻| 久久久久久久大尺度免费视频| 国产久久久一区二区三区| 尤物成人国产欧美一区二区三区| 国产久久久一区二区三区| 国产在线一区二区三区精| 如何舔出高潮| 久久久久久人妻| 欧美成人午夜免费资源| 久久精品国产鲁丝片午夜精品| 亚洲国产精品国产精品| videossex国产| 美女xxoo啪啪120秒动态图| 视频中文字幕在线观看| 国产免费视频播放在线视频| 久久 成人 亚洲| 亚洲av男天堂| 国产欧美另类精品又又久久亚洲欧美| 七月丁香在线播放| 蜜臀久久99精品久久宅男| 国产精品偷伦视频观看了| 亚洲欧美成人综合另类久久久| 久久久久久伊人网av| 大码成人一级视频| 国产精品爽爽va在线观看网站| 在线观看三级黄色| 王馨瑶露胸无遮挡在线观看| 大香蕉97超碰在线| 欧美xxxx性猛交bbbb| 大陆偷拍与自拍| 麻豆精品久久久久久蜜桃| 成人无遮挡网站| 99久国产av精品国产电影| 国产人妻一区二区三区在| 少妇人妻精品综合一区二区| 久久久久国产精品人妻一区二区| 男人和女人高潮做爰伦理| 欧美日韩精品成人综合77777| 国产免费一级a男人的天堂| 草草在线视频免费看| 日本午夜av视频| 日本免费在线观看一区| 亚洲av国产av综合av卡| av在线观看视频网站免费| 又大又黄又爽视频免费| 大片电影免费在线观看免费| 久久久久人妻精品一区果冻| av一本久久久久| 2021少妇久久久久久久久久久| 日韩 亚洲 欧美在线| 国产亚洲精品久久久com| 91久久精品国产一区二区成人| 超碰av人人做人人爽久久| 超碰av人人做人人爽久久| 在线观看免费高清a一片| 亚洲婷婷狠狠爱综合网| 国产淫片久久久久久久久| 亚洲欧美清纯卡通| 五月天丁香电影| 国产精品一二三区在线看| 免费大片18禁| 久久国内精品自在自线图片| 男人爽女人下面视频在线观看| 99热全是精品| av线在线观看网站| 国产在线一区二区三区精| 国产日韩欧美在线精品| 日韩在线高清观看一区二区三区| 国产毛片在线视频| 高清午夜精品一区二区三区| 久久这里有精品视频免费| 亚洲欧美清纯卡通| 水蜜桃什么品种好| 超碰97精品在线观看| 久热久热在线精品观看| 成年女人在线观看亚洲视频| 91aial.com中文字幕在线观看| 美女国产视频在线观看| 99热国产这里只有精品6| 欧美日韩精品成人综合77777| 亚洲自偷自拍三级| 一级毛片电影观看| 一级毛片久久久久久久久女| 亚洲高清免费不卡视频| av不卡在线播放| 99久国产av精品国产电影| 欧美丝袜亚洲另类| 亚洲精品色激情综合| 国产综合精华液| 国产 精品1| 国产在线一区二区三区精| 国产欧美日韩精品一区二区| h日本视频在线播放| 亚洲精品久久久久久婷婷小说| 久久久久久久久大av| 日本爱情动作片www.在线观看| 简卡轻食公司| 99久久人妻综合| 不卡视频在线观看欧美| 亚洲中文av在线| 久久99精品国语久久久| 18禁裸乳无遮挡动漫免费视频| av在线播放精品| 免费黄色在线免费观看| 妹子高潮喷水视频| 欧美变态另类bdsm刘玥| 国产熟女欧美一区二区| 最黄视频免费看| 在线观看国产h片| 一区二区三区四区激情视频| 狂野欧美白嫩少妇大欣赏| 少妇精品久久久久久久| 尾随美女入室| 亚洲av欧美aⅴ国产| 日韩大片免费观看网站| 欧美zozozo另类| 一级黄片播放器| 国产一区有黄有色的免费视频| 免费看光身美女| 街头女战士在线观看网站| 麻豆国产97在线/欧美| 久久韩国三级中文字幕| 亚洲欧美日韩另类电影网站 | 黄色怎么调成土黄色| 色视频在线一区二区三区| 欧美xxxx性猛交bbbb| 国产深夜福利视频在线观看| 蜜桃在线观看..| 亚洲第一区二区三区不卡| 久久99热这里只有精品18| 精品一区二区免费观看| 国产高清国产精品国产三级 | 在线观看三级黄色| 亚洲成人手机| 狠狠精品人妻久久久久久综合| 亚洲图色成人| 亚洲欧美日韩东京热| 国产精品福利在线免费观看| 亚洲天堂av无毛| 国产亚洲最大av| 麻豆乱淫一区二区| av视频免费观看在线观看| 99热这里只有是精品50| 不卡视频在线观看欧美| 两个人的视频大全免费| 国产欧美另类精品又又久久亚洲欧美| 欧美成人一区二区免费高清观看| 亚洲国产最新在线播放| 欧美极品一区二区三区四区| 高清毛片免费看| 日韩电影二区| 欧美xxⅹ黑人| 久久久国产一区二区| 性色avwww在线观看| 少妇熟女欧美另类| 99久久精品一区二区三区| 精品人妻偷拍中文字幕| 激情五月婷婷亚洲| 乱系列少妇在线播放| 在线观看人妻少妇| 身体一侧抽搐| 亚洲精品乱久久久久久| 婷婷色综合www| 成人美女网站在线观看视频| av天堂中文字幕网| 一本色道久久久久久精品综合| 91在线精品国自产拍蜜月| 亚洲av综合色区一区| 丝瓜视频免费看黄片| 女人久久www免费人成看片| 在线观看一区二区三区| 国精品久久久久久国模美| 亚洲av在线观看美女高潮| 国产高清不卡午夜福利| 国产午夜精品久久久久久一区二区三区| 免费久久久久久久精品成人欧美视频 | 91狼人影院| 国产欧美亚洲国产| 又大又黄又爽视频免费| 在线观看免费高清a一片| 人妻少妇偷人精品九色| 亚洲av中文字字幕乱码综合| 大香蕉久久网| 少妇被粗大猛烈的视频| 亚洲精品中文字幕在线视频 | 亚洲人与动物交配视频| 久久这里有精品视频免费| 制服丝袜香蕉在线| 国产精品秋霞免费鲁丝片| 成人无遮挡网站| 欧美丝袜亚洲另类| 天堂俺去俺来也www色官网| av视频免费观看在线观看| 国产爱豆传媒在线观看| 各种免费的搞黄视频| 欧美xxxx黑人xx丫x性爽| 狠狠精品人妻久久久久久综合| 99精国产麻豆久久婷婷| 国产免费视频播放在线视频| 国产精品不卡视频一区二区| 日韩中文字幕视频在线看片 | 欧美日韩综合久久久久久| 亚洲,一卡二卡三卡| 亚洲性久久影院| 欧美另类一区| 成人综合一区亚洲| 久久久久久伊人网av| 国产久久久一区二区三区| av在线观看视频网站免费| 在线观看免费高清a一片| 国产高清三级在线| 国产av码专区亚洲av| 国产精品久久久久久久电影| 久久久久精品久久久久真实原创| 国模一区二区三区四区视频| 人妻系列 视频| 伊人久久国产一区二区| 久久精品国产亚洲av天美| 欧美日韩亚洲高清精品| av在线老鸭窝| 日韩在线高清观看一区二区三区| 97精品久久久久久久久久精品| 免费黄网站久久成人精品| 黄色配什么色好看| 亚洲综合色惰| 麻豆成人午夜福利视频| 人妻 亚洲 视频| 久久久久久久久久人人人人人人| 亚洲伊人久久精品综合| 亚洲国产色片| 久久久久久久久久成人| 精品人妻视频免费看| av卡一久久| 国产精品麻豆人妻色哟哟久久| 久久久亚洲精品成人影院| 成人18禁高潮啪啪吃奶动态图 | 一本色道久久久久久精品综合| 春色校园在线视频观看| 小蜜桃在线观看免费完整版高清| 久久99精品国语久久久| 舔av片在线| av线在线观看网站| 国产男女超爽视频在线观看| 又粗又硬又长又爽又黄的视频| 久久精品国产亚洲网站| 亚洲在久久综合| 久久久久人妻精品一区果冻| 欧美极品一区二区三区四区|