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

    UN-U3Si2復(fù)合燃料的多尺度蠕變行為研究

    2022-12-22 14:36:10楊青峰謝哲驍陳平張靜高士鑫丁國琛周毅尹春雨丁淑蓉何梁孫
    核技術(shù) 2022年12期
    關(guān)鍵詞:機(jī)制實(shí)驗(yàn)模型

    楊青峰謝哲驍陳 平張 靜高士鑫丁國琛周 毅尹春雨丁淑蓉何 梁孫 丹

    1(中國核動(dòng)力研究設(shè)計(jì)院核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室 成都 610213)

    2(復(fù)旦大學(xué)航空航天系力學(xué)與工程仿真研究所 上海 200433)

    自福島發(fā)生核事故以來,開發(fā)具有更高安全性的燃料一直是核領(lǐng)域關(guān)注的焦點(diǎn)。耐事故燃料應(yīng)能提升反應(yīng)堆在嚴(yán)重事故工況下的安全性,同時(shí)也能保持甚至提升正常運(yùn)行狀態(tài)下的性能[1?2]。與傳統(tǒng)的UO2燃料相比,UN、U3Si2等燃料不僅具有更高的鈾密度和更高的熱導(dǎo)率,且其熱導(dǎo)率隨著溫度的升高而增加[3?4]。UN、U3Si2燃料的高導(dǎo)熱性能不僅保障了燃料棒在正常運(yùn)行時(shí)具有較低的熱儲(chǔ)能,而且在事故工況下可以較快地導(dǎo)出衰變熱;高鈾密度則提高了中子經(jīng)濟(jì)性,可全面提高反應(yīng)堆的運(yùn)行利潤[5?8]。因此,UN、U3Si2高鈾密度燃料及其復(fù)合燃料是極具研發(fā)潛力的耐事故燃料。

    UN-U3Si2復(fù)合燃料芯塊在服役過程中會(huì)產(chǎn)生復(fù)雜的輻照-熱-力耦合行為,其宏觀蠕變性能直接影響其安全性。燃料的蠕變會(huì)引起應(yīng)力松弛效應(yīng),影響復(fù)合燃料芯塊內(nèi)部的力學(xué)場及芯塊與包殼之間的力學(xué)相互作用[9?10]。在不同的應(yīng)力和輻照條件下,UN-U3Si2復(fù)合燃料的宏觀蠕變應(yīng)變及其內(nèi)部主導(dǎo)的蠕變機(jī)制將會(huì)產(chǎn)生變化,與UN和U3Si2多晶燃料的蠕變性能密切相關(guān)。因此,對(duì)復(fù)合燃料的多尺度蠕變行為進(jìn)行研究非常重要,可以為復(fù)合燃料芯塊的優(yōu)化設(shè)計(jì)及性能分析奠定基礎(chǔ)。

    已有的燃料蠕變性能研究主要針對(duì)UO2燃料展開,關(guān)于UN、U3Si2高鈾密度燃料蠕變的理論及實(shí)驗(yàn)研究較少,關(guān)于UN-U3Si2復(fù)合燃料的蠕變模型研究尚未見報(bào)道。Zeisser等[11]對(duì)UN燃料進(jìn)行了輻照條件下的單軸壓縮蠕變?cè)囼?yàn),獲得了溫度520~1 050℃、應(yīng)力5~63 MPa范圍內(nèi)的實(shí)驗(yàn)結(jié)果,但未進(jìn)一步建立蠕變率模型。Hayes等[12]通過對(duì)UN高溫壓縮蠕變實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合,得到了UN蠕變率模型,此模型已被應(yīng)用于愛達(dá)荷國家工程和環(huán)境實(shí)驗(yàn)室開展的UN燃料堆內(nèi)熱-力學(xué)行為的計(jì)算模擬。此外,Brucklacher等[13]也針對(duì)UN燃料開展了蠕變實(shí)驗(yàn)研究,得到了不同條件下的實(shí)驗(yàn)數(shù)據(jù),但未進(jìn)一步發(fā)展蠕變率模型。Mercado[14]對(duì)U3Si2進(jìn)行了高溫下的單軸壓縮蠕變?cè)囼?yàn),獲得了溫度范圍為850~1 050℃,應(yīng)力范圍為30~80 MPa下的實(shí)驗(yàn)數(shù)據(jù),并擬合得到了與溫度、應(yīng)力相關(guān)的蠕變率模型;Freeman[15]也利用Mercado得到的實(shí)驗(yàn)數(shù)據(jù),采用另外的函數(shù)形式得到了新的蠕變率模型;Yingling等[16]進(jìn)一步考慮了晶粒尺寸的影響,得到了與溫度、應(yīng)力和晶粒尺寸關(guān)聯(lián)的蠕變率模型,并引入U(xiǎn)3Si2-SiC燃料堆內(nèi)熱-力學(xué)行為的數(shù)值模擬研究。Metzger[17]基于已有的蠕變變形主導(dǎo)機(jī)制的分析,提出了U3Si2的蠕變率模型,尚未得到實(shí)驗(yàn)的驗(yàn)證。目前的這些蠕變率模型是經(jīng)驗(yàn)公式,沒有綜合考慮輻照蠕變及熱蠕變的貢獻(xiàn),也沒有全面考慮不同應(yīng)力及溫度條件下的蠕變機(jī)制進(jìn)行建模,因此需要結(jié)合已有報(bào)道的所有實(shí)驗(yàn)結(jié)果建立新的蠕變率模型。

    利用文獻(xiàn)中已有的蠕變實(shí)驗(yàn)數(shù)據(jù),基于不同條件下的主控蠕變機(jī)制建立了UN和U3Si2的蠕變率模型,并基于均勻化理論獲得UN-U3Si2復(fù)合燃料的宏觀等效蠕變應(yīng)變與各組分貢獻(xiàn)的關(guān)聯(lián)模型。建立了組分隨機(jī)分布的有限元模型,計(jì)算模擬了UN-U3Si2復(fù)合燃料輻照條件下的拉伸蠕變實(shí)驗(yàn),獲得了各組分對(duì)宏觀蠕變應(yīng)變的定量貢獻(xiàn)及主控的蠕變機(jī)制,并考察了燃料輻照腫脹對(duì)復(fù)合燃料多尺度蠕變行為的影響。

    1 UN和U3Si2燃料的輻照熱-力學(xué)性能模型

    UN-U3Si2復(fù)合燃料由UN和U3Si2組成,本節(jié)給出其彈性性能和輻照腫脹模型。

    1.1 彈性常數(shù)

    UN的彈性模量(EUN)和泊松比(νUN)分別為[12]:

    式中:ρ是理論密度比,%,本文取100;T是溫度,K。

    U3Si2的彈性模量(EU3Si2)和泊松比(νU3Si2)分別為120 GPa和0.177[14]。

    1.2 輻照腫脹模型

    UN和U3Si2的輻照體積腫脹采用如下經(jīng)驗(yàn)?zāi)P兔枋觯?8?19]:

    式中:Bu表示燃耗,單位為FIMA(Fissions Per Initial Heavy Metal Atoms)。

    2 復(fù)合燃料的多尺度蠕變模型

    2.1 多晶材料相的蠕變率模型

    由于目前缺乏完善的UN和U3Si2蠕變率模型,本研究結(jié)合文獻(xiàn)報(bào)道的實(shí)驗(yàn)結(jié)果,綜合考慮輻照蠕變和熱蠕變的貢獻(xiàn),并基于不同應(yīng)力及輻照工況條件下的蠕變機(jī)制分析,建立兩種組分材料的新蠕變率模型。

    UN和U3Si2均為多晶材料,多晶材料的蠕變變形主要有以下幾種典型機(jī)制[20]:空位擴(kuò)散機(jī)制[21]、位錯(cuò)攀移輔助滑移機(jī)制[22]、位錯(cuò)湮滅機(jī)制[23]。這些機(jī)制與晶粒內(nèi)部各類缺陷的運(yùn)動(dòng)以及晶粒之間的相互作用有關(guān),在不同條件下的主導(dǎo)蠕變機(jī)制不同,也因材料而異。

    對(duì)于UN材料,根據(jù)實(shí)驗(yàn)結(jié)果[11?12]獲知,在中低溫下空位擴(kuò)散蠕變機(jī)制主控;在低溫下主要發(fā)生輻照蠕變,在中等溫度下輻照蠕變和熱擴(kuò)散蠕變均有貢獻(xiàn);在高溫下以位錯(cuò)蠕變機(jī)制為主,不受裂變率的影響。通過蠕變率與應(yīng)力之間的關(guān)系可以判斷,高溫蠕變主要由位錯(cuò)湮滅機(jī)制所導(dǎo)致??瘴粩U(kuò)散機(jī)制下,蠕變率與應(yīng)力成正比關(guān)系,與晶粒尺寸的平方成反比關(guān)系;位錯(cuò)湮滅機(jī)制下,蠕變率與應(yīng)力的4.5次方成正比關(guān)系,與晶粒尺寸無關(guān)?;谶@兩種蠕變機(jī)制的典型特征,結(jié)合實(shí)驗(yàn)結(jié)果建立如下等效蠕變率模型:

    式中:ε?cr是宏觀蠕變 率,s?1;f?是裂變率,fissions·(cm3·s)?1;R是理想氣體常數(shù),取8.314 J·(mol·K)?1;σ是應(yīng)力,MPa;d是晶粒尺寸,μm;A1=6.592 6×10?23,A2=1.613 8×10?2,QA=114 007 J·mol?1,B1=2.285 6×10?3,QB=329 040 J·mol?1。

    對(duì)于U3Si2材料,根據(jù)實(shí)驗(yàn)結(jié)果[16]分析得出,其蠕變變形由空位擴(kuò)散機(jī)制和位錯(cuò)攀移輔助滑移機(jī)制共同作用,且在低溫下主要來自于輻照蠕變。對(duì)于這種燃料,其空位熱擴(kuò)散蠕變和位錯(cuò)攀移輔助滑移蠕變兩部分均受裂變率和溫度的影響[24]。目前缺乏U3Si2燃料的輻照蠕變實(shí)驗(yàn)數(shù)據(jù),根據(jù)文獻(xiàn)[17]得出,特定裂變率下的輻照蠕變率可由轉(zhuǎn)變溫度下的熱蠕變推導(dǎo)得到??瘴粩U(kuò)散機(jī)制下,蠕變率與應(yīng)力成正比關(guān)系,與晶粒尺寸的平方成反比關(guān)系;位錯(cuò)攀移輔助滑移機(jī)制下,蠕變率與應(yīng)力的3次方成正比關(guān)系,與晶粒尺寸無關(guān)。所以,結(jié)合實(shí)驗(yàn)結(jié)果得到了如下的等效蠕變率模型:

    式中:A1=6.524 8×10?23,A2=4.792,QA=172 647 J·mol?1,B1=3.361 6×10?29,B2=3.578 6×10?6,QB=171 453 J·mol?1。

    式(5)和式(6)中的各系數(shù)均通過文獻(xiàn)中的蠕變實(shí)驗(yàn)數(shù)據(jù)[11?12,16]及反推的輻照蠕變數(shù)據(jù)擬合得到,為解決多參數(shù)擬合的問題,針對(duì)低溫和中高溫區(qū)采用了分批次擬合的方法。在低溫下,溫度相關(guān)項(xiàng)的量級(jí)遠(yuǎn)小于輻照相關(guān)項(xiàng),所以先擬合出與溫度無關(guān)的輻照主導(dǎo)項(xiàng)的系數(shù),再利用中高溫下的實(shí)驗(yàn)數(shù)據(jù),將總?cè)渥兟蕼p去輻照相關(guān)項(xiàng)的貢獻(xiàn)后,得到熱蠕變率,進(jìn)一步擬合得到溫度相關(guān)項(xiàng)的系數(shù)。

    圖1(a)、(b)分別展示了采用式(5)和式(6)計(jì)算得到的UN和U3Si2蠕變率,并同時(shí)給出了文獻(xiàn)上的實(shí)驗(yàn)結(jié)果[11?12,16]。其中的實(shí)驗(yàn)結(jié)果包含了輻照條件下的數(shù)據(jù)及堆外蠕變實(shí)驗(yàn)結(jié)果,也具有不同晶粒尺寸的核燃料試樣的蠕變實(shí)驗(yàn)結(jié)果。可以看出,模型計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果吻合較好,證明了所建立模型的合理有效性。

    圖1 UN蠕變(a)和U3Si2蠕變(b)模型計(jì)算值與實(shí)驗(yàn)值的對(duì)比情況(彩色見網(wǎng)絡(luò)版)Fig.1 Comparison of the calculated and experimental values of the UN creep model(a)and the U3Si2 creep model(b)(color online)

    2.2 復(fù)合燃料的宏觀等效蠕變模型

    在宏觀外載作用下,UN-U3Si2復(fù)合燃料試樣的宏觀蠕變變形來自于UN和U3Si2的蠕變貢獻(xiàn)??梢曰诰鶆蚧碚摚?5?27],通過選取代表性體元模擬單軸拉伸蠕變實(shí)驗(yàn),通過對(duì)細(xì)觀非均勻分布的蠕變應(yīng)變進(jìn)行體積平均,計(jì)算得到復(fù)合燃料的宏觀等效蠕變應(yīng)變,并同時(shí)獲得UN和U3Si2兩相對(duì)復(fù)合燃料宏觀蠕變的定量貢獻(xiàn)及主控蠕變機(jī)制。

    對(duì)于計(jì)算模擬中的每個(gè)增量步,以增量步起始時(shí)刻的構(gòu)形作為參考構(gòu)形,各應(yīng)變?cè)隽繚M足:

    式中:Δεij為復(fù)合燃料芯塊的總應(yīng)變?cè)隽浚?Δεij)f和(Δεij)m分別為彌散相U3Si2和基體相UN對(duì)各自體積所平均得到的總應(yīng)變?cè)隽?;fν和fm分別為彌散相和基體相的現(xiàn)時(shí)體積分?jǐn)?shù)。

    由于蠕變應(yīng)變?yōu)槠睆埩?,而輻照腫脹應(yīng)變?yōu)榍蛐螐埩浚瑒t復(fù)合燃料的宏觀體積變化主要來自各組分材料輻照腫脹的貢獻(xiàn)。由于彈性應(yīng)變的貢獻(xiàn)可以忽略,可得復(fù)合燃料拉伸方向(x方向)的宏觀等效蠕變應(yīng)變?cè)隽繚M足:

    3 復(fù)合燃料單軸拉伸蠕變實(shí)驗(yàn)?zāi)M的有限元模型

    基于§1中的UN和U3Si2燃料彈性模型和輻照腫脹模型,以及§2.1中建立的UN和U3Si2多晶燃料蠕變率模型,建立了旋轉(zhuǎn)坐標(biāo)系下的三維大變形增量本構(gòu)關(guān)系,推導(dǎo)了應(yīng)力更新算法和一致剛度模量,并編寫了UMAT子程序,在ABAQUS有限元計(jì)算平臺(tái)實(shí)現(xiàn)了UN-U3Si2復(fù)合燃料的單軸拉伸蠕變實(shí)驗(yàn)的有限元模擬。

    目前,UN-U3Si2復(fù)合燃料芯塊主要是U3Si2和UN粉末共同燒結(jié)制備而成[28],芯塊的金相結(jié)構(gòu)剖面如圖2所示[28]。UN-U3Si2復(fù)合燃料可以視作以U3Si2為彌散相、以UN為基體相的彌散型核燃料,在結(jié)構(gòu)上類似顆粒復(fù)合燃料,但U3Si2顆粒的形狀不規(guī)則、大小不均勻,在UN基體中也呈現(xiàn)非均勻分布。

    圖2 UN-U3Si2復(fù)合燃料芯塊組織Fig.2 Microstructure of UN-U3Si2 composite fuel pellets

    根據(jù)復(fù)合燃料的微觀結(jié)構(gòu)特征和兩相分布情況,暫不考慮孔隙的存在,本文針對(duì)立方體形代表性體元發(fā)展了如下的隨機(jī)建模方法:

    1)將整個(gè)立方體體元?jiǎng)澐志鶆虻牧骟w網(wǎng)格;

    2)在計(jì)算前,借助外部程序進(jìn)行隨機(jī)建模,獲得兩相材料的空間分布情況,并輸出文本文檔;

    3)在計(jì)算開始時(shí)刻,通過UMAT讀取上述文本文檔和單元積分點(diǎn)坐標(biāo),根據(jù)積分點(diǎn)坐標(biāo)判斷材料屬性,并記錄材料狀態(tài)變量;

    4)在每一個(gè)增量步開始時(shí),讀取最初記錄的材料狀態(tài)變量,并調(diào)用對(duì)應(yīng)的材料子程序。

    圖3展示了通過外部程序進(jìn)行隨機(jī)分布建模的流程示意圖。20%彌散相含量的復(fù)合燃料隨機(jī)分布有限元模型如圖4所示。

    圖3 彌散相隨機(jī)分布建模的流程示意圖Fig.3 Flowchart diagram of random distribution modeling

    圖4 隨機(jī)分布模型示意圖Fig.4 Schematic diagram of random distribution model

    根據(jù)細(xì)觀力學(xué)的思想,復(fù)合燃料在宏觀上由代表性體元周期性排列得到。在宏觀外載作用下,代表性體元的變形亦周期性分布[27]。對(duì)于內(nèi)部彌散相隨機(jī)分布的代表性體元,由于其沒有幾何上的對(duì)稱性,所以需要施加位移約束使其滿足邊界條件的周期性。周期性位移邊界條件可表示為[29]:

    其中:如圖5所示,N0和Nj(j=1,2,3)是4個(gè)主節(jié)點(diǎn)和是相對(duì)面上的對(duì)應(yīng)點(diǎn),正負(fù)號(hào)代表相應(yīng)點(diǎn)所在面的外法線方向沿坐標(biāo)軸的正向或負(fù)向。uj+、uj?、Uj、U0分別是對(duì)應(yīng)點(diǎn)和主節(jié)點(diǎn)的位移。將主節(jié)點(diǎn)N0設(shè)置為固定,約束主節(jié)點(diǎn)Nj在其他兩個(gè)方向上的位移,用于約束剛體平動(dòng)和轉(zhuǎn)動(dòng)。以x方向?yàn)槔旆较?,體元在法線為x方向的兩個(gè)面上分別施加大小為P的均布拉力。需要指出,由于需要同時(shí)滿足位移周期性條件,雖然在周期性邊界上施加了均勻拉力,周期性邊界上的實(shí)際應(yīng)力非均勻分布且滿足應(yīng)力的周期性。

    圖5 代表性體元的周期性邊界條件及施加載荷示意圖Fig.5 Schematic diagram of periodic boundary conditions and applied loads for cube-shaped voxels

    4 結(jié)果與討論

    4.1 輻照腫脹對(duì)復(fù)合燃料多尺度蠕變行為的影響

    本文模擬復(fù)合燃料輻照條件下的單軸拉伸蠕變實(shí)驗(yàn),考慮了兩種材料的輻照腫脹貢獻(xiàn),以更好地模擬復(fù)合燃料的輻照效應(yīng)。雖然在計(jì)算宏觀等效蠕變應(yīng)變時(shí),已經(jīng)去除了輻照腫脹所引起的體積應(yīng)變貢獻(xiàn),但是各組分材料的輻照腫脹會(huì)影響復(fù)合燃料細(xì)觀的應(yīng)力狀態(tài),進(jìn)而影響多尺度的蠕變變形。為了研究輻照腫脹對(duì)復(fù)合燃料多尺度蠕變行為的影響,本文也對(duì)不考慮組分材料輻照腫脹的算例進(jìn)行了計(jì)算模擬,獲得了在均布拉力50 MPa、溫度500 K、裂變率1020fissions·m?3·s?1條件下的結(jié)果。

    圖6給出了考慮與不考慮輻照腫脹兩種情況下的宏觀等效蠕變應(yīng)變隨時(shí)間的演化結(jié)果??梢钥吹?,考慮輻照腫脹情況下的結(jié)果略大于不考慮輻照腫脹下的結(jié)果,其相對(duì)偏差隨著燃耗的增長而逐漸增大。在裂變密度達(dá)到4.32×1027fissions·m?3時(shí),考慮輻照腫脹計(jì)算得到的等效蠕變應(yīng)變?yōu)?.037 78,不考慮輻照腫脹計(jì)算得到的等效蠕變應(yīng)變?yōu)?.035 772,偏差為5.2%,總體上差別不大。

    圖6 考慮與不考慮輻照腫脹情況下宏觀等效蠕變應(yīng)變的對(duì)比Fig.6 Comparison of macroscopic equivalent creep strain with and without irradiation swelling

    圖7(a)和(b)分別給出了考慮與不考慮輻照腫脹計(jì)算得到的細(xì)觀von Mises等效應(yīng)力分布圖,為復(fù)合燃料在裂變密度達(dá)到4.32×1027fissions·m?3時(shí)的結(jié)果??梢钥吹?,對(duì)于考慮輻照腫脹的情況,由于彌散相和基體相均會(huì)產(chǎn)生腫脹,且體積腫脹不一致,導(dǎo)致兩相的界面處發(fā)生強(qiáng)烈的力學(xué)相互作用,最大的von Mises應(yīng)力近300 MPa,大約為施加的外部均布拉力的6倍;最小的von Mises應(yīng)力約為30 MPa,低于外部施加的均布拉力。對(duì)于不考慮輻照腫脹的情況,兩相的von Mises應(yīng)力基本呈均勻分布,大小接近于施加的外部均布拉力。

    圖7 裂變密度達(dá)到4.32×1027 fissions·m?3時(shí)的von Mises應(yīng)力云圖(剖面圖)(a)考慮輻照腫脹,(b)不考慮輻照腫脹Fig.7 Von Mises stress contourplots(section view)at the fission density of 4.32×1027 fissions·m?3(a)With radiation swelling,(b)Without radiation swelling

    圖8分別針對(duì)考慮與不考慮輻照腫脹兩種情況,給出了裂變密度達(dá)到4.32×1027fissions·m?3時(shí)的細(xì)觀等效蠕變分布云圖。可以看到,對(duì)于考慮輻照腫脹的情況,非均勻分布的von Mises等效應(yīng)力導(dǎo)致了等效蠕變的非均勻分布,等效蠕變應(yīng)變的最大值約為最小值的3.2倍。對(duì)于不考慮輻照的情況,等效蠕變應(yīng)變雖然也非均勻,但最大和最小值相差不大。

    圖8 裂變密度達(dá)到4.32×1027 fissions·m?3時(shí)的等效蠕變?cè)茍D(剖面圖)(a)考慮輻照腫脹,(b)不考慮輻照腫脹Fig.8 Equivalent creep contour plots(section view)at the fission density of 4.32×1027 fissions·m?3(a)With radiation swelling,(b)Without radiation swelling

    結(jié)合以上結(jié)果可知,由于兩相材料的輻照腫脹不同,兩相材料之間需要變形協(xié)調(diào),導(dǎo)致兩相之間存在較為強(qiáng)烈的力學(xué)相互作用。但由于離開界面處的基體應(yīng)力及顆粒相的應(yīng)力被削弱了,所以兩相材料的輻照腫脹對(duì)復(fù)合燃料宏觀等效蠕變的影響可以忽略。

    4.2 復(fù)合燃料宏觀蠕變應(yīng)變的細(xì)觀主導(dǎo)機(jī)制分析

    通過均勻化方法,可以得到彌散相和基體相對(duì)總?cè)渥儜?yīng)變的貢獻(xiàn),也可以得到不同機(jī)制蠕變的占比情況,分別如圖9、10所示??梢姡w相對(duì)總?cè)渥兊呢暙I(xiàn)較大。裂變密度達(dá)到4.32×1027fissions·m?3時(shí),基體蠕變?cè)诳側(cè)渥冎械恼急燃s為79%,這主要由于基體相在復(fù)合燃料的體積分?jǐn)?shù)較大。由圖10可以看出,輻照擴(kuò)散蠕變占主導(dǎo),其占比隨時(shí)間逐漸減??;輻照位錯(cuò)蠕變也產(chǎn)生了可觀的貢獻(xiàn),其占比隨時(shí)間逐漸增大,這主要是由于組分材料的輻照腫脹使得復(fù)合燃料內(nèi)部產(chǎn)生了越來越強(qiáng)烈的力學(xué)相互作用;熱擴(kuò)散蠕變和熱位錯(cuò)蠕變基本沒有貢獻(xiàn)。本研究所考慮的溫度較低,熱蠕變基本可以忽略,所以在本文的條件下輻照蠕變機(jī)制占主導(dǎo)。

    圖9 彌散相和基體相對(duì)總?cè)渥儜?yīng)變的貢獻(xiàn)Fig.9 Contribution of dispersed phase and matrix to total creep strain

    圖10 各蠕變貢獻(xiàn)的占比Fig.10 Proportions of different creep contribution

    5 結(jié)語

    本文基于已有文獻(xiàn)中的蠕變實(shí)驗(yàn)結(jié)果,綜合考慮了輻照蠕變及熱蠕變的貢獻(xiàn),通過分析各種條件下的主導(dǎo)蠕變機(jī)制建立了UN和U3Si2多晶燃料的蠕變率模型,證明了模型的合理有效性;同時(shí)利用均勻化理論建立了UN-U3Si2復(fù)合燃料宏觀等效蠕變和兩相貢獻(xiàn)的關(guān)聯(lián)模型;并發(fā)展了隨機(jī)建模方法,實(shí)現(xiàn)了UN-U3Si2復(fù)合燃料輻照條件下的單軸拉伸蠕變實(shí)驗(yàn)的有限元模擬,獲得兩相材料在復(fù)合燃料宏觀蠕變中的定量貢獻(xiàn),并分析了主控的蠕變機(jī)制。同時(shí),考察了組分材料的輻照腫脹對(duì)復(fù)合燃料多尺度蠕變行為的影響,得到的主要結(jié)論如下:

    1)考慮和不考慮輻照腫脹的計(jì)算結(jié)果表明:由于組分材料的輻照腫脹行為的差異性,導(dǎo)致復(fù)合燃料的顆粒和基體相產(chǎn)生了強(qiáng)烈的力學(xué)相互作用,在裂變密度達(dá)到4.32×1027fissions·m?3時(shí)的最大von Mises應(yīng)力約為外部均布拉力的6倍;由于同時(shí)存在細(xì)觀應(yīng)力被削弱的區(qū)域,組分材料不同的輻照腫脹對(duì)宏觀等效蠕變應(yīng)變的影響很小,可以忽略。

    2)在均布拉力50 MPa、溫度500 K、裂變率1020fissions·m?3·s?1的條件下,基體相對(duì)總?cè)渥兊呢暙I(xiàn)較大,在裂變密度達(dá)到4.32×1027fissions·m?3時(shí),拉伸方向的基體蠕變?cè)诳側(cè)渥冎械恼急燃s79%;在此參數(shù)范圍內(nèi),輻照擴(kuò)散機(jī)制蠕變占主導(dǎo),其占比隨燃耗的增長逐漸減?。欢椪瘴诲e(cuò)蠕變的占比隨燃耗的增長逐漸增大,歸功于組分輻照腫脹對(duì)細(xì)觀應(yīng)力場的影響;熱擴(kuò)散機(jī)制蠕變和熱位錯(cuò)機(jī)制蠕變基本沒有貢獻(xiàn)。

    作者貢獻(xiàn)聲明楊青峰:論文整體設(shè)計(jì),起草文章;謝哲驍:起草文章,模型計(jì)算;陳平:對(duì)文章作批評(píng)性審閱,研究經(jīng)費(fèi)支持;張靜:模型計(jì)算指導(dǎo);高士鑫:對(duì)文章作批評(píng)性審閱;丁國琛:模型計(jì)算指導(dǎo);周毅:論文整體設(shè)計(jì);尹春雨:論文整體設(shè)計(jì);丁淑蓉:研究指導(dǎo),論文修改;何梁:數(shù)據(jù)處理;孫丹:數(shù)據(jù)分析。

    猜你喜歡
    機(jī)制實(shí)驗(yàn)模型
    一半模型
    記一次有趣的實(shí)驗(yàn)
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    做個(gè)怪怪長實(shí)驗(yàn)
    自制力是一種很好的篩選機(jī)制
    文苑(2018年21期)2018-11-09 01:23:06
    3D打印中的模型分割與打包
    NO與NO2相互轉(zhuǎn)化實(shí)驗(yàn)的改進(jìn)
    實(shí)踐十號(hào)上的19項(xiàng)實(shí)驗(yàn)
    太空探索(2016年5期)2016-07-12 15:17:55
    破除舊機(jī)制要分步推進(jìn)
    男人狂女人下面高潮的视频| 日韩在线高清观看一区二区三区| 国产综合懂色| 尤物成人国产欧美一区二区三区| 国产精品麻豆人妻色哟哟久久 | 少妇熟女欧美另类| 在线免费观看不下载黄p国产| 99热这里只有是精品在线观看| 精品无人区乱码1区二区| 极品教师在线视频| 特大巨黑吊av在线直播| 国产av一区在线观看免费| 成人综合一区亚洲| av专区在线播放| 亚洲丝袜综合中文字幕| 亚洲精华国产精华液的使用体验| 日韩,欧美,国产一区二区三区 | 能在线免费看毛片的网站| 日韩亚洲欧美综合| 非洲黑人性xxxx精品又粗又长| 嘟嘟电影网在线观看| 精品99又大又爽又粗少妇毛片| 欧美丝袜亚洲另类| 99热这里只有精品一区| 又黄又爽又刺激的免费视频.| 国产精品野战在线观看| 欧美+日韩+精品| 亚洲经典国产精华液单| 国产免费一级a男人的天堂| 日日啪夜夜撸| 成人三级黄色视频| 成人国产麻豆网| 亚洲熟妇中文字幕五十中出| 国产精品女同一区二区软件| 亚洲精品影视一区二区三区av| 国产亚洲精品久久久com| 国产精品美女特级片免费视频播放器| 男人狂女人下面高潮的视频| 久久人妻av系列| 国产精品人妻久久久久久| 欧美激情久久久久久爽电影| 成年女人看的毛片在线观看| 午夜福利网站1000一区二区三区| 国产三级在线视频| 九九久久精品国产亚洲av麻豆| 天堂√8在线中文| 一个人免费在线观看电影| 性色avwww在线观看| 国产麻豆成人av免费视频| 国产精品久久电影中文字幕| 精品久久久噜噜| 尾随美女入室| av天堂中文字幕网| 中文字幕av成人在线电影| 久久久久精品久久久久真实原创| 亚洲av福利一区| 我要搜黄色片| 久久久久九九精品影院| 少妇熟女欧美另类| 国产综合懂色| 久久久久久久久久成人| 精品国产露脸久久av麻豆 | 亚洲经典国产精华液单| 91久久精品电影网| 日本午夜av视频| 亚洲熟妇中文字幕五十中出| 69人妻影院| 男人和女人高潮做爰伦理| 免费看av在线观看网站| 老女人水多毛片| 久久久久免费精品人妻一区二区| 在线观看66精品国产| 亚洲精品亚洲一区二区| 日韩精品有码人妻一区| 国产成人福利小说| 狠狠狠狠99中文字幕| 亚洲综合精品二区| 三级国产精品片| 久久精品国产鲁丝片午夜精品| 丝袜喷水一区| 免费av观看视频| 亚洲国产欧洲综合997久久,| 自拍偷自拍亚洲精品老妇| 日日撸夜夜添| 久久综合国产亚洲精品| 国产精华一区二区三区| 三级国产精品片| 国产成人a区在线观看| 久久久久久久久久成人| 国产免费一级a男人的天堂| 在线观看一区二区三区| 国产精品日韩av在线免费观看| 少妇猛男粗大的猛烈进出视频 | av国产免费在线观看| 国产成人精品婷婷| 亚洲av熟女| 国产精品美女特级片免费视频播放器| 国产精品永久免费网站| 在线播放无遮挡| 你懂的网址亚洲精品在线观看 | 韩国高清视频一区二区三区| 国产高清三级在线| 不卡视频在线观看欧美| 永久免费av网站大全| 国产一区二区在线观看日韩| 成年女人永久免费观看视频| 亚洲高清免费不卡视频| 最近的中文字幕免费完整| 国产精品久久电影中文字幕| 日韩在线高清观看一区二区三区| 天堂中文最新版在线下载 | 麻豆av噜噜一区二区三区| 男女那种视频在线观看| 亚洲欧美日韩东京热| 中文天堂在线官网| 精品久久久久久电影网 | 国产精品99久久久久久久久| 身体一侧抽搐| 久久亚洲精品不卡| 国产精品女同一区二区软件| 国产人妻一区二区三区在| 99久久人妻综合| 欧美日韩国产亚洲二区| 欧美精品一区二区大全| 精品一区二区免费观看| 久久综合国产亚洲精品| 熟女人妻精品中文字幕| 日韩欧美国产在线观看| 内射极品少妇av片p| ponron亚洲| 美女黄网站色视频| 一区二区三区四区激情视频| 一级黄片播放器| 一级爰片在线观看| 成年版毛片免费区| 免费看a级黄色片| 色哟哟·www| 插阴视频在线观看视频| 婷婷色av中文字幕| 日韩一区二区视频免费看| 你懂的网址亚洲精品在线观看 | 中文字幕免费在线视频6| 日韩欧美在线乱码| 黄色配什么色好看| 国产乱人视频| 五月玫瑰六月丁香| 能在线免费看毛片的网站| 欧美极品一区二区三区四区| 97超碰精品成人国产| 两个人的视频大全免费| 国产精品av视频在线免费观看| 欧美性感艳星| 午夜爱爱视频在线播放| 国产大屁股一区二区在线视频| 天天一区二区日本电影三级| 欧美又色又爽又黄视频| 久久人人爽人人爽人人片va| 能在线免费看毛片的网站| 免费黄网站久久成人精品| 久久久久久久久久成人| 汤姆久久久久久久影院中文字幕 | 免费播放大片免费观看视频在线观看 | 超碰av人人做人人爽久久| 亚洲国产精品成人久久小说| 99久久人妻综合| 国产精品一区二区三区四区免费观看| 自拍偷自拍亚洲精品老妇| 久久久精品94久久精品| 国产老妇伦熟女老妇高清| 国产色爽女视频免费观看| 在线观看66精品国产| 黄片wwwwww| 少妇裸体淫交视频免费看高清| 在线免费观看不下载黄p国产| 欧美潮喷喷水| 亚洲欧美精品自产自拍| 国产精品久久久久久精品电影| 成人午夜高清在线视频| 91精品国产九色| 国产精品99久久久久久久久| 2021天堂中文幕一二区在线观| 国产一区二区在线av高清观看| 18禁裸乳无遮挡免费网站照片| 欧美成人午夜免费资源| 人妻制服诱惑在线中文字幕| 亚洲国产精品sss在线观看| 成人漫画全彩无遮挡| 人妻少妇偷人精品九色| 国产色婷婷99| 日本-黄色视频高清免费观看| 精品一区二区免费观看| 国产免费视频播放在线视频 | 夜夜爽夜夜爽视频| 中文亚洲av片在线观看爽| 热99在线观看视频| 国产精品久久久久久精品电影| 国产精品久久久久久av不卡| 国产69精品久久久久777片| 人人妻人人看人人澡| 22中文网久久字幕| 亚洲国产欧美人成| 成人av在线播放网站| 夜夜爽夜夜爽视频| 国内精品美女久久久久久| 日韩三级伦理在线观看| 免费观看的影片在线观看| 亚洲天堂国产精品一区在线| 美女脱内裤让男人舔精品视频| 99久久精品热视频| 国产精品久久视频播放| 别揉我奶头 嗯啊视频| 白带黄色成豆腐渣| 久久精品久久久久久噜噜老黄 | 国产精品.久久久| 热99在线观看视频| 免费观看性生交大片5| 色视频www国产| 尾随美女入室| 亚洲精品日韩av片在线观看| 午夜亚洲福利在线播放| 国产单亲对白刺激| 日本熟妇午夜| 五月伊人婷婷丁香| 久久国内精品自在自线图片| 国产黄片视频在线免费观看| 色播亚洲综合网| 欧美一区二区精品小视频在线| 久久久成人免费电影| 国产精品一及| 九草在线视频观看| 九九在线视频观看精品| 亚洲av成人av| 中文字幕制服av| 成年女人看的毛片在线观看| 国产一区二区在线av高清观看| 亚洲欧美中文字幕日韩二区| 国产成人一区二区在线| 99热网站在线观看| 简卡轻食公司| 黄色一级大片看看| 床上黄色一级片| 中文字幕制服av| 永久免费av网站大全| 青青草视频在线视频观看| 特大巨黑吊av在线直播| 男人舔奶头视频| 18禁在线无遮挡免费观看视频| av专区在线播放| 日韩欧美精品v在线| 99久久精品国产国产毛片| 中国美白少妇内射xxxbb| 色吧在线观看| 亚洲国产成人一精品久久久| 天天躁夜夜躁狠狠久久av| 亚洲av熟女| 欧美潮喷喷水| 美女国产视频在线观看| 久热久热在线精品观看| 亚洲精品亚洲一区二区| 插阴视频在线观看视频| 日本熟妇午夜| 激情 狠狠 欧美| 午夜久久久久精精品| 又粗又爽又猛毛片免费看| 久久99热这里只频精品6学生 | 精品一区二区免费观看| 国产单亲对白刺激| 日本免费一区二区三区高清不卡| 一区二区三区免费毛片| 国产精品99久久久久久久久| 一级毛片电影观看 | 老司机影院毛片| 欧美性猛交╳xxx乱大交人| 日韩欧美国产在线观看| 色噜噜av男人的天堂激情| 亚洲成人久久爱视频| 国语自产精品视频在线第100页| 久久久久久久午夜电影| 国产在线一区二区三区精 | 日本欧美国产在线视频| 亚洲性久久影院| 我的女老师完整版在线观看| 欧美另类亚洲清纯唯美| 亚洲自拍偷在线| 亚洲欧美日韩高清专用| 免费av不卡在线播放| 听说在线观看完整版免费高清| 久久久a久久爽久久v久久| 91av网一区二区| 欧美xxxx性猛交bbbb| 亚洲成色77777| 91久久精品国产一区二区三区| av黄色大香蕉| 最新中文字幕久久久久| 成人二区视频| 国产av在哪里看| 亚洲四区av| 黄色一级大片看看| 狠狠狠狠99中文字幕| 日韩一区二区视频免费看| 亚洲伊人久久精品综合 | 久久热精品热| 成人av在线播放网站| 亚洲成色77777| 搡老妇女老女人老熟妇| 欧美日本视频| 男女国产视频网站| 久久精品熟女亚洲av麻豆精品 | 色哟哟·www| 级片在线观看| 热99在线观看视频| 国产精品伦人一区二区| 日产精品乱码卡一卡2卡三| 国内精品宾馆在线| 精品国产露脸久久av麻豆 | 大香蕉97超碰在线| 国产老妇女一区| 亚洲aⅴ乱码一区二区在线播放| 日韩欧美在线乱码| 国产精品三级大全| 免费观看精品视频网站| 国产亚洲91精品色在线| 观看免费一级毛片| 日日撸夜夜添| 亚洲国产精品国产精品| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕精品亚洲无线码一区| 国产午夜福利久久久久久| 少妇猛男粗大的猛烈进出视频 | 欧美变态另类bdsm刘玥| 九色成人免费人妻av| 老师上课跳d突然被开到最大视频| 91av网一区二区| 韩国高清视频一区二区三区| 日日啪夜夜撸| 免费黄色在线免费观看| 日本免费一区二区三区高清不卡| 乱系列少妇在线播放| av天堂中文字幕网| 亚洲精品影视一区二区三区av| 高清视频免费观看一区二区 | 在线免费观看不下载黄p国产| 精品熟女少妇av免费看| 久久久久九九精品影院| 九九在线视频观看精品| 亚洲精品国产av成人精品| av天堂中文字幕网| 性色avwww在线观看| 国产精品一区二区三区四区免费观看| 久久99热这里只频精品6学生 | 亚洲av免费在线观看| 亚洲国产精品久久男人天堂| 免费播放大片免费观看视频在线观看 | 久久久久九九精品影院| 日本免费在线观看一区| 春色校园在线视频观看| 看片在线看免费视频| av又黄又爽大尺度在线免费看 | 免费观看a级毛片全部| 精品无人区乱码1区二区| 亚洲欧美成人综合另类久久久 | av又黄又爽大尺度在线免费看 | 成人国产麻豆网| 亚洲国产成人一精品久久久| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕亚洲精品专区| 国产一区二区在线观看日韩| 亚洲av电影不卡..在线观看| 3wmmmm亚洲av在线观看| 麻豆精品久久久久久蜜桃| 夫妻性生交免费视频一级片| 亚洲伊人久久精品综合 | 精品欧美国产一区二区三| 国产精品伦人一区二区| 免费观看精品视频网站| 最近中文字幕2019免费版| 91av网一区二区| 搡老妇女老女人老熟妇| 成年av动漫网址| 国产精品.久久久| 精品久久久久久电影网 | 九草在线视频观看| 国产成人免费观看mmmm| 欧美色视频一区免费| 天堂av国产一区二区熟女人妻| 亚洲真实伦在线观看| 一区二区三区四区激情视频| av线在线观看网站| 国产精品久久久久久久电影| 五月玫瑰六月丁香| 好男人在线观看高清免费视频| 九草在线视频观看| 欧美xxxx黑人xx丫x性爽| 国语自产精品视频在线第100页| 毛片一级片免费看久久久久| 夫妻性生交免费视频一级片| 99久久精品国产国产毛片| 午夜福利在线观看免费完整高清在| 国产在视频线精品| 国产午夜福利久久久久久| 久久久色成人| 一区二区三区高清视频在线| 中文字幕精品亚洲无线码一区| 国产精品国产三级专区第一集| 99热6这里只有精品| 午夜日本视频在线| 日本黄色片子视频| 哪个播放器可以免费观看大片| 男人狂女人下面高潮的视频| 国产一区二区在线观看日韩| 高清在线视频一区二区三区 | 中文字幕亚洲精品专区| 亚洲伊人久久精品综合 | 久久精品夜夜夜夜夜久久蜜豆| 高清毛片免费看| 久久久久久久久大av| av黄色大香蕉| 99在线人妻在线中文字幕| 国产淫语在线视频| 综合色av麻豆| 欧美激情国产日韩精品一区| 欧美zozozo另类| 免费av不卡在线播放| 亚洲欧美清纯卡通| 成人毛片60女人毛片免费| 亚洲国产精品成人久久小说| 亚洲精品自拍成人| 女人被狂操c到高潮| 成人三级黄色视频| 九九爱精品视频在线观看| av免费在线看不卡| 黄色欧美视频在线观看| 久久久久久久午夜电影| 搞女人的毛片| 毛片一级片免费看久久久久| 午夜福利在线观看免费完整高清在| av女优亚洲男人天堂| 在线观看av片永久免费下载| 桃色一区二区三区在线观看| 国产一区二区在线观看日韩| 亚洲国产精品久久男人天堂| 国产伦精品一区二区三区视频9| 久久精品国产鲁丝片午夜精品| 国产私拍福利视频在线观看| 真实男女啪啪啪动态图| 亚洲熟妇中文字幕五十中出| 亚洲精华国产精华液的使用体验| 日本av手机在线免费观看| 欧美激情久久久久久爽电影| 国模一区二区三区四区视频| 一区二区三区高清视频在线| 亚洲一区高清亚洲精品| 亚洲欧洲国产日韩| 国产高清三级在线| 一级毛片电影观看 | 日韩强制内射视频| 少妇熟女aⅴ在线视频| 最后的刺客免费高清国语| 大话2 男鬼变身卡| 韩国高清视频一区二区三区| 菩萨蛮人人尽说江南好唐韦庄 | 九九久久精品国产亚洲av麻豆| 一级爰片在线观看| 一级二级三级毛片免费看| 亚洲国产精品合色在线| 一本久久精品| 婷婷色麻豆天堂久久 | 亚洲精品乱久久久久久| 欧美变态另类bdsm刘玥| 丝袜美腿在线中文| 天堂中文最新版在线下载 | 青青草视频在线视频观看| h日本视频在线播放| 最近手机中文字幕大全| 日日啪夜夜撸| 直男gayav资源| 精品午夜福利在线看| 人人妻人人看人人澡| 校园人妻丝袜中文字幕| 国产老妇女一区| 国产精品美女特级片免费视频播放器| 婷婷色综合大香蕉| 欧美一级a爱片免费观看看| 91久久精品国产一区二区三区| 看免费成人av毛片| 男人的好看免费观看在线视频| 中文天堂在线官网| 亚洲国产成人一精品久久久| 国产高清视频在线观看网站| 日韩成人伦理影院| 亚洲精品影视一区二区三区av| 永久免费av网站大全| 又爽又黄无遮挡网站| 久久久成人免费电影| 国产伦精品一区二区三区四那| 一级毛片久久久久久久久女| 老司机影院毛片| 国产美女午夜福利| 亚洲欧美成人综合另类久久久 | 久久99热这里只频精品6学生 | 精品久久久久久久末码| 三级国产精品欧美在线观看| 韩国高清视频一区二区三区| 日日撸夜夜添| 中文资源天堂在线| av国产久精品久网站免费入址| 纵有疾风起免费观看全集完整版 | 国产高清三级在线| 一级毛片久久久久久久久女| 国产精品久久久久久久电影| 韩国高清视频一区二区三区| 色吧在线观看| 亚洲性久久影院| 亚洲综合精品二区| 欧美日本亚洲视频在线播放| 亚洲国产精品国产精品| 国产精品爽爽va在线观看网站| 国语自产精品视频在线第100页| 简卡轻食公司| 国内揄拍国产精品人妻在线| 联通29元200g的流量卡| 欧美zozozo另类| 日韩欧美 国产精品| 一级av片app| 亚洲自偷自拍三级| 两性午夜刺激爽爽歪歪视频在线观看| 国产精华一区二区三区| 日韩欧美精品v在线| 99久久人妻综合| 淫秽高清视频在线观看| 国产毛片a区久久久久| 午夜精品在线福利| 国产一区二区三区av在线| 99热这里只有精品一区| 亚洲内射少妇av| 国产精品一区二区三区四区免费观看| 简卡轻食公司| 婷婷六月久久综合丁香| 国产又色又爽无遮挡免| 亚洲五月天丁香| 波多野结衣高清无吗| 亚洲国产精品久久男人天堂| 国产视频首页在线观看| 九九在线视频观看精品| 麻豆精品久久久久久蜜桃| 最近最新中文字幕免费大全7| 欧美日韩国产亚洲二区| 午夜免费激情av| 欧美一区二区精品小视频在线| 亚洲人成网站在线播| 变态另类丝袜制服| 亚洲经典国产精华液单| 99热6这里只有精品| 高清日韩中文字幕在线| 国产精品不卡视频一区二区| 国产精华一区二区三区| 亚洲精品成人久久久久久| 日韩视频在线欧美| 色播亚洲综合网| 国产三级中文精品| 国产av一区在线观看免费| av在线天堂中文字幕| 在线天堂最新版资源| 日本黄色片子视频| 亚洲熟妇中文字幕五十中出| www.av在线官网国产| av女优亚洲男人天堂| 国产精品av视频在线免费观看| 成人午夜精彩视频在线观看| 能在线免费看毛片的网站| 91aial.com中文字幕在线观看| 国产精品麻豆人妻色哟哟久久 | 精品久久久久久成人av| 狂野欧美激情性xxxx在线观看| 在线a可以看的网站| 韩国av在线不卡| 高清日韩中文字幕在线| 伊人久久精品亚洲午夜| 精品一区二区三区视频在线| 搡女人真爽免费视频火全软件| 麻豆一二三区av精品| 亚洲国产最新在线播放| 久久精品人妻少妇| 级片在线观看| 久久99热这里只频精品6学生 | 亚洲最大成人手机在线| 日本午夜av视频| 美女被艹到高潮喷水动态| 国产成人a∨麻豆精品| 午夜日本视频在线| 亚洲av成人精品一区久久| 最近中文字幕高清免费大全6| 亚洲精品aⅴ在线观看| 亚洲久久久久久中文字幕| 91久久精品国产一区二区三区| 亚洲中文字幕日韩| 午夜福利成人在线免费观看| 97在线视频观看| 日本黄色视频三级网站网址| 男人狂女人下面高潮的视频| 久久久久久国产a免费观看| 夫妻性生交免费视频一级片| 国产单亲对白刺激| 久久久久久久午夜电影| 亚洲国产高清在线一区二区三| 亚洲av电影不卡..在线观看| 激情 狠狠 欧美| 1000部很黄的大片| 级片在线观看| 亚洲av熟女| 国产亚洲5aaaaa淫片| 国产成人精品一,二区| 亚洲五月天丁香| 亚洲精品日韩av片在线观看| 小蜜桃在线观看免费完整版高清|