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

    改性釓/硼/聚乙烯納米復合材料的制備及對中子和伽馬射線的屏蔽性能

    2022-07-18 02:38:04趙盛霍志鵬鐘國強張宏胡立群
    高等學?;瘜W學報 2022年6期
    關鍵詞:伽馬射線模擬計算中子

    趙盛, 霍志鵬, 鐘國強, 張宏, 胡立群

    改性釓/硼/聚乙烯納米復合材料的制備及對中子和伽馬射線的屏蔽性能

    趙盛1,2, 霍志鵬1, 鐘國強1, 張宏1, 胡立群1

    (1. 中國科學院合肥物質科學研究院等離子體物理研究所, 合肥 230026; 2. 中國科學技術大學, 合肥 230026)

    采用偶聯(lián)劑-甲基丙烯酰氧基丙基三甲氧基硅烷(KH570)對納米氧化釓(nanoGd2O3)表面進行包覆改性, 通過熱壓成型法制備了一種新型的改性納米氧化釓/碳化硼/高密度聚乙烯(M-nanoGd2O3/B4C/HDPE)復合材料. 傅里葉變換紅外光譜(FTIR)、 掃描電子顯微鏡(SEM)和X射線能譜(EDS)分析結果表明, nanoGd2O3成功被偶聯(lián)劑改性, 且改性nanoGd2O3在聚乙烯基體內的界面相容性和分散性顯著提高. 熱重分析(TGA)、 差示掃描量熱分析(DSC)和力學拉伸實驗表明, 改性nanoGd2O3的引入增強了復合材料的熱穩(wěn)定性, 提高了復合材料的拉伸強度、 楊氏模量和斷裂伸長率. 對復合材料的中子和伽馬射線屏蔽性能進行了實驗測試和蒙特卡羅模擬計算, 研究了nanoGd2O3改性、 材料形狀和材料厚度對屏蔽性能的影響. 結果表明, 界面相容性和分散性優(yōu)良的nanoGd2O3能夠有效提高中子及伽馬射線屏蔽率. 方形M-nanoGd2O3/B4C/HDPE材料在厚度為11.7 cm時中子屏蔽率達到90%, 在厚度為13.5 cm時伽馬射線屏蔽率達到70%.

    氧化釓; 偶聯(lián)劑; 表面改性; 聚乙烯; 復合材料; 輻射防護

    核能作為清潔、 低碳、 高效的基荷能源, 在降低化石能源依賴、 改善霧霾天氣及應對全球變暖中起到正面作用. 然而核反應堆中大量的核輻射粒子(如粒子,粒子, 中子和伽馬射線等)會對人體的生殖系統(tǒng)、 免疫系統(tǒng)和神經系統(tǒng)及自然環(huán)境等造成直接損害, 例如福島核電站爆炸造成的核污染給當?shù)鼐用窈铜h(huán)境帶來了嚴重災難. 中子和伽馬射線穿透能力很強, 必須采用專門的防護措施. 聚乙烯(PE)因具備優(yōu)異的化學穩(wěn)定性、 良好的中子慢化效果、 質輕和體積小等優(yōu)點而被廣泛應用于核輻射防護領域[1~3], 而將核輻射屏蔽無機功能填料引入聚乙烯基體可以提升其屏蔽性能和力學性能[4,5].10B的熱中子吸收截面很高, 天然B中10B的含量約為20%, B及含B化合物如碳化硼(B4C)及氮化硼(BN)等通常作為熱中子吸收劑摻雜到PE基體中[6~8]. 高原子序數(shù)元素如鉛(Pb)易和中低能量伽馬光子產生光電效應和康普頓散射, 具備優(yōu)異的伽馬射線吸收能力, 可作為伽馬射線吸收劑[9~11]. 鉛硼聚乙烯(Pb/B4C/PE)是中國核動力院研發(fā)的一種中子及伽馬射線復合屏蔽材料[12], 由于具備綜合屏蔽性能優(yōu)良、 質輕和能簡化屏蔽結構等優(yōu)點而被廣泛用于核輻射混合場的屏蔽. 但鉛具有生物毒性且對環(huán)境不友好, 這在一定程度上限制了含鉛材料的應用. 此外, 無機材料與PE基體間界面相容性差, 導致其力學性能下降和無機顆粒團聚, 進而影響復合材料的屏蔽性能. 研究表明, 利用偶聯(lián)劑對填料表面改性是改善無機-有機界面相容性的有效途徑[13~15].

    稀土元素釓(Gd)作為高Z元素具備一定的伽馬射線屏蔽能力, 其同位素155Gd和157Gd還具備超高的熱中子吸收截面[16], 由于兼具熱中子吸收和伽馬射線屏蔽的能力而被廣泛應用, 通常以氧化釓(Gd2O3)形式存在自然界中. Gd2O3具有熱穩(wěn)定性優(yōu)良、 毒性弱且無污染性等特點, 是一種良好的代鉛材料. ?rim等[17]通過向PE中摻雜納米Gd2O3成功制備了無鉛復合屏蔽材料, 發(fā)現(xiàn)隨著填料濃度的增加, 中子和伽馬射線的透射率降低,納米粒子的加入還提高了復合材料的拉伸強度和模量等力學性能. Baykara等[18]在聚酰亞胺(PI)內部引入了納米Gd2O3和納米BN, 結果表明, 納米復合材料的中子透過率和伽馬衰減系數(shù)隨著填料含量和材料厚度的增加而降低. 此外, 與純PI相比, 復合材料的抗拉強度和熱性能均有所提升. ?rim等和Baykara等均指出納米顆粒的引入使得復合材料的抗拉強度優(yōu)于原基體. 然而, 納米Gd2O3作為無機物與PE基體間的界面相容性較差、 比表面積大和表面能高等特性使得其易在PE內部團聚, 這影響了復合材料的屏蔽性能和力學性能. 因此, 有必要對納米Gd2O3粉末進行表面修飾來改善其在PE內部的分散性.

    本文采用偶聯(lián)劑-甲基丙烯酰氧基丙基三甲氧基硅烷(KH570)對nanoGd2O3進行表面改性, 并采用熱壓成型法制備了一種M-nanoGd2O3/B4C/HDPE無機-有機納米復合材料, 研究了nanoGd2O3表面改性對復合材料微觀結構、 熱學性能、 力學性能及中子和伽馬射線屏蔽性能的影響, 并結合蒙特卡羅模擬計算和輻射屏蔽實驗研究了復合材料的輻射屏蔽機制.

    1 實驗部分

    1.1 試劑與儀器

    碳化硼(B4C), 粒徑1~10 μm, 分析純, 阿拉丁試劑有限公司; 納米氧化釓(nanoGd2O3, 粒徑200~500 nm)、 KH570和無水乙醇(CH3CH2OH), 分析純, 麥克林生化科技有限公司; 乙酸溶液(CH3COOH), 1 mol/L, 廣州臻萃質檢技術服務有限公司; 高密度聚乙烯(HDPE), 500目,n=40000~300000,=0.96 g/cm3, 深圳華創(chuàng)塑膠原料有限公司.

    SFM-11型Ⅴ形混料機, 合肥科晶材料技術有限公司; R32212017型熱壓機, 武漢啟恩科技發(fā)展有限責任公司; Nicolet 8700型傅里葉變換紅外光譜儀(FTIR), 美國尼高利儀器公司, 分辨率為4 cm-1, 掃描范圍400~4000 cm-1; G500型肖特基場發(fā)射掃描電子顯微鏡(SEM), 美國FEI電子光學公司; 牛津 Aztec系列能量色散X射線能譜儀(EDS), 英國牛津儀器公司, 分辨率(Mn)127 eV, 分析元素范圍: B—U; TA Discovery TGA型熱重分析儀(TGA), 美國TA公司, 氮氣氛圍, 升溫速率10 ℃/min, 測試溫度范圍: 20~800 ℃; TA Q2000 DSC型差示掃描量熱儀(DSC), 美國TA公司, 氮氣氛圍, 升溫速率10 ℃/min, 測試溫度范圍: 20~250 ℃; Instron 9657型多功能電子萬能材料力學性能實驗機, 美國英斯特朗公司, 樣品規(guī)格為GBT 1040.2-2006 1BA型, 拉伸速率: 1 mm/min; Cf-252中子源、 混凝土屏蔽體和SmartREM型便攜式高靈敏中子測量儀, 美國福祿克公司, 其中Cf-252中子源的活度為6×107Bq, 平均能量為2.45 MeV, 中子測量儀的靈敏度為0.74 cps·h/μSv, 測量范圍: 10 nSv/h~100 mSv/h; Cs-137伽馬源、 鉛室和FHZ 672 E-10型多功能輻射測量儀, 美國賽默飛世爾公司, 其中Cs-137伽馬源的活度為3.7×104Bq, 光子能量為0.662 MeV, 探測器靈敏度2000 h/(μSv·s), 測量范圍1 nSv~100 μSv/h.

    1.2 改性納米氧化釓的制備

    將nanoGd2O3放入電熱恒溫鼓風干燥箱中, 于110 ℃干燥3 h; 將無水乙醇和超純水的混合液(體積比為8∶1)倒入燒杯中, 加入適量nanoGd2O3, 超聲分散2 h后, 將硅烷偶聯(lián)劑KH570倒入燒杯中 [m(KH570)∶m(nanoGd2O3)=1∶10]; 邊攪拌邊逐滴加入1 mol/L的乙酸溶液, 利用pH計測量混合溶液的pH值, 當pH值為3~4時, 停止攪拌和滴加. 將上述混合液轉移至三口燒瓶中并加入攪拌子, 利用恒溫磁力攪拌器在80 ℃條件下恒溫攪拌10 h; 在5000 r/min條件下離心5 min, 用乙醇溶解沉淀物, 重復以上醇洗操作, 直至試管上清液透明; 隨后加入超純水進行水洗、 離心操作, 直至試管上清液澄清; 將沉淀物置于于60 ℃ ?0.08 MPa真空度下真空干燥12 h, 得到表面改性的nanoGd2O3(M-nanoGd2O3). 改性過程如Scheme 1所示. 上述實驗流程的相關參數(shù)是經過多次實驗優(yōu)化后確定的優(yōu)選結果.

    Scheme 1synthetic route for the modified nanoGd2O3

    1.3 納米復合屏蔽材料的制備

    將M-nanoGd2O3(或nanoGd2O3), B4C和HDPE按照質量比1∶2∶7加入SFM-11型V形混料機的料筒中, 并加入研磨球, 以120 r/min轉速混料1 h; 隨后將混合粉末轉移至R32212017型熱壓機模具中, 并墊上聚四氟乙烯脫模紙進行熱壓: 在2 MPa的壓力下先預壓5 min; 然后加壓至30 MPa, 溫度設定為120 ℃, 保溫保壓20 min; 隨后加壓至40 MPa, 溫度設定為180 ℃, 保溫保壓30 min; 上述過程完成后, 進行冷卻、 卸壓、 脫模, 即可得到M-nanoGd2O3/B4C/HDPE(或nanoGd2O3/B4C/HDPE)納米復合屏蔽材料. 上述熱壓工藝流程的相關參數(shù)均為經過多次熱壓步驟后確定的優(yōu)選結果.

    1.4 計算方法

    屏蔽率(, %)是表征材料對中子或伽馬射線吸收程度的物理量, 其計算公式為[19,20]:

    式中:(μSv/h)為放置材料時中子或伽馬射線探測器的劑量讀數(shù);0(μSv/h)為未放置材料時中子或伽馬射線探測器的劑量讀數(shù).

    線性衰減系數(shù)(cm-1)表示單位厚度的屏蔽材料對射線的吸收程度, 其計算公式為[20]:

    式中:(cm)表示屏蔽材料的厚度.

    質量衰減系數(shù)m(cm2/g)表示單位質量的屏蔽材料對射線的吸收程度, 其計算公式為[20]:

    式中:(g/cm3)為屏蔽材料的密度.

    半值層厚度HVL(cm)表示使射線的強度減少到初始輻射值一半時所需要的屏蔽材料的厚度, 其計算公式[21]為:

    伽馬實驗中的數(shù)據(jù)擬合曲線公式為[22]:

    式中:(cm)為材料的厚度;為積累因子, 在理想窄束條件下,γ=1, 寬束條件下,>1.

    中子實驗中的數(shù)據(jù)擬合曲線公式為[23]:

    式中:n為寬束情況下中子的積累因子,(cm-1)為復合材料的宏觀總截面, 表示入射中子在單位厚度介質中發(fā)生相互作用的幾率.

    1.5 蒙特卡羅模擬計算

    應用中國科學院核能安全技術研究所FDS團隊開發(fā)的SuperMC(Super Monte Carlo Simulation Program for Nuclear and Radiation Process)軟件對中子和伽馬射線屏蔽實驗裝置進行幾何和材料建模, 并進行了中子和伽馬光子的輸運模擬計算和Cf-252中子源在出射孔平面處的能譜計算.

    2 結果與討論

    2.1 材料微結構表征

    圖1(A)給出KH570, nanoGd2O3粉末和M-nanoGd2O3粉末的紅外光譜. nanoGd2O3的和M-nanoGd2O3的譜圖均在550 cm-1附近處出現(xiàn)了Gd—O的伸縮振動峰[24], 這是氧化釓粉末的固有特征峰. KH570分子結構中存在甲基、 亞甲基、 碳碳雙鍵、 碳氧雙鍵、 硅氧單鍵及硅碳單鍵等特征官能團和化學鍵, 在其紅外譜圖上均能找到與上述官能團和化學鍵對應的特征峰. 在M-nanoGd2O3的紅外譜圖中除了Gd—O特征峰, 還在2962和2927 cm-1處出現(xiàn)了雙峰[25], 這是偶聯(lián)劑中R基團的—CH3峰和—CH2—峰, 1702 cm-1處為R基團中的C=O伸縮振動峰, 1176 cm-1處為Si—O鍵的伸縮振動峰[26], 878 cm-1處為Si—C伸縮振動峰[27], 1421 cm-1處為R基團中的—CH2—彎曲振動峰, 1557 cm-1處為R基團中的C=C伸縮振動峰[28]. 以上新出現(xiàn)的峰表明偶聯(lián)劑分子與Gd2O3表面發(fā)生了相互作用, Gd2O3表面的—OH與水解后的偶聯(lián)劑分子的—OH發(fā)生了脫水縮合, 實現(xiàn)了KH570對nanoGd2O3的偶聯(lián)接枝. 圖1(B)為nanoGd2O3和M-nanoGd2O3在0~600 ℃范圍內的熱重曲線. 氧化釓的熱分解溫度極高, 在2300 ℃以上, 但圖中顯示nanoGd2O3和M-nanoGd2O3在升溫過程中均發(fā)生了質量損失. 其中nanoGd2O3在升溫結束時的質量為初始的99.1%, 這是粉末試樣中含有未除凈的水造成的; M-nanoGd2O3在升溫結束時的質量為初始的83.3%, 這是由于粉末試樣中含有未除凈的水和包覆在nanoGd2O3表面的KH570分子在升溫過程中逐步降解造成的. 根據(jù)圖1(B)可得到偶聯(lián)劑的接枝率為15.8%.

    Fig.1 FTIR spectraof KH570(a), nanoGd2O3(b) and M?nanoGd2O3(c)(A) and TG curves of nanoGd2O3(a) and M?nanoGd2O3(b)(B)

    進一步對nanoGd2O3和M-nanoGd2O3進行EDS分析(圖2). 圖2(A)顯示觀測區(qū)域只含有元素Gd和O, 均來自于Gd2O3; 圖2(B)顯示觀測區(qū)域包含元素Gd, O, C和Si, 其中Gd來自于Gd2O3, C和Si來自于KH570, O元素來自Gd2O3和KH570. nanoGd2O3的表面改性有助于解決填料與HDPE基體的界面相容性問題.

    Fig.2 EDS spectra of nanoGd2O3(A) and M?nanoGd2O3(B)

    圖3(A)為nanoGd2O3的SEM照片. 可見nanoGd2O3呈碎屑狀, 由于納米粒子比表面積大和表面能高的特性, 部分nanoGd2O3發(fā)生了團聚. 圖3(B)為M-nanoGd2O3的SEM照片. 可以看出, M-nanoGd2O3的團聚體表面有一層有機物質附著. 圖3(C)為nanoGd2O3/B4C/HDPE(質量比1∶2∶7)的斷面SEM照片, 在紅圈內部區(qū)域出現(xiàn)了大量顆粒團聚現(xiàn)象. 圖3(D)為相同質量比的M-nanoGd2O3/B4C/HDPE斷面SEM照片, 基體內部較為明顯的Gd2O3和B4C顆粒已在圖中圈出, 與圖3(C)相比, 顆粒的團聚堆積現(xiàn)象得到一定程度的改善, 說明偶聯(lián)劑改性能夠提高填料在基體內部的分散性.

    Fig.3 SEM images of nanoGd2O3 and nanoGd2O3/B4C/HDPE fracture surface

    (A) NanoGd2O3; (B) M-nanoGd2O3; (C) nanoGd2O3/B4C/HDPE; (D) M-nanoGd2O3/B4C/HDPE.

    2.2 材料熱穩(wěn)定性

    圖4(A)給出復合材料及純HDPE的熱重分析結果. 以初始分解溫度5%(℃, 材料發(fā)生5%降解的溫度)作為評估熱穩(wěn)定性的特征溫度. 結果顯示, HDPE, nanoGd2O3/B4C/HDPE和M-nanoGd2O3/B4C/HDPE的5%分別為450, 458和461 ℃. 與HDPE相比, nanoGd2O3/B4C/HDPE和M-nanoGd2O3/B4C/HDPE的TGA曲線整體偏向高溫側, 初始分解溫度也更高, 這是由于在高溫下納米復合材料中的nanoGd2O3粒子和B4C粒子在一定程度上阻礙了HDPE鏈的鏈段運動, 延遲了聚合物鏈段在高溫下的熱效應. 3種材料的質量損失主要集中在400~510 ℃, 這是由于HDPE分子鏈被熱降解, 而納米復合材料中剩余的殘留物為熱穩(wěn)定性優(yōu)良的nanoGd2O3和B4C. 圖4(B)給出復合材料及純HDPE的DSC測試結果. 吸熱峰的前沿最大斜率處的切線與基線延長線的交點對應的溫度為初始熔融溫度onset(℃),p(℃)為峰值溫度. 與HDPE的p(p1, 129 ℃)相比, nanoGd2O3/B4C/HDPE和M-nanoGd2O3/B4C/HDPE的p得到提高; 與nanoGd2O3/B4C/HDPE的p(p2, 132 ℃)相比, M-nanoGd2O3/B4C/HDPE的p(p3, 138 ℃)得到進一步 提高, 表明偶聯(lián)劑的改性有利于提高復合材料的熱穩(wěn)定性. 這是由于nanoGd2O3改性后極性變弱, 增強了與非極性HDPE的親和性, 導致與HDPE界面的結合力增強, 進一步限制了HDPE鏈的相對運動.

    Fig.4 TGA(A) and DSC(B) curves of three kinds of materials

    2.3 材料的力學性能

    圖5為3種材料的應力-應變曲線, 表1列出了相關力學數(shù)據(jù).復合材料的模量和拉伸強度相較HDPE均有所提升, 而斷裂強度均有所降低, 說明向HDPE基體內摻雜適量nanoGd2O3和B4C有利于提高其拉伸強度和剛性. 這是因為適量的填料能夠阻礙聚乙烯鏈的鏈段運動, 同時對基體內部的裂紋具有釘扎作用[29], 有效阻止了裂紋的擴展. 此外, 填料與HDPE界面之間的應力傳導機制有助于分擔部分復合材料的應力載荷. 3種材料的楊氏模量分別為511.39, 1281.78和1357.90 MPa; 拉伸強度分別為13.23, 15.37和17.40 MPa; 斷裂伸長率分別為13.5%, 8.2%和8.8%. 其中M-nanoGd2O3/B4C/HDPE的拉伸強度相較nanoGd2O3/B4C/HDPE提升了13.2%, 楊氏模量和斷裂伸長率分別提升了5.9%和7.3%. 這是由于填料的改性能夠提升其在基體內部的分散性, 改善了顆粒團聚現(xiàn)象, 減輕了材料內部的應力集中效應, 從而起到了增強增韌的效果; 另一方面, 無機填料的改性能夠提升有機-無機界面的結合強度[30,31], 使得應力可以更有效地通過界面?zhèn)鬟f, 緩解了因無機填料加入聚合物基體中導致的相分離問題, 從而提高了改性復合材料的拉伸強度和斷裂伸長率.

    Fig.5 Stress?strain curves of three kinds of materials at the rate of 1 mm/min

    Table 1 Results of tensile stress test for three kinds of materials

    2.4 輻射屏蔽性能

    圖6為中子和伽馬射線屏蔽實驗裝置示意圖, 本文中所有關于中子和伽馬射線的屏蔽率或透過率的測試結果均為5次測試結果的平均值. 圖7為所制備的兩種形狀規(guī)格的M-nanoG2O3/B4C/HDPE納米復合材料的照片.

    Fig.6 Schematic diagram of shielding test devices for neutron(A) and gamma?ray(B)

    圖8(A)是Cf-252中子源的標準能譜[32], 平均能量為2.54 MeV. 圖8(B)是Cf-252中子源經混凝土屏蔽筒準直后在出射孔平面處的模擬計算能譜, 平均能量為1.17 MeV. 可以看出, 準直后在出射孔平面處的能譜分布相比標準能譜平均能量有所降低, 譜圖整體向低能段偏移, 慢中子(<1 KeV)份額增加, 快中子(0.1 MeV<<10 MeV)和中能中子 (1 KeV<<100 KeV)總體份額有所降低. 本研究中進行屏蔽性能測試的所有復合材料均緊貼于出射孔平面外側, 因此, 材料所受照射的中子能量分布更接近于圖8(B).

    Fig.8 Standard energy spectrum of Cf?252 source(A) and energy spectrum after collimation of Cf?252 source(B)

    Fig.7 Photographs of the M?nanoG2O3 /B4C /HDPE materials with different shapes

    圖9為探測器探測到的中子能量分布譜圖(SuperMC模擬計算結果), 其中圖9(A)為未加載屏蔽材料時從準直孔處出射并照射在探測器表面上的中子能譜, 圖9(B)和圖9(C)分別為透過4.5 cm和15 cm厚度方形復合材料的出射中子照射在探測器表面的中子能譜. 根據(jù)各能區(qū)中子份額統(tǒng)計數(shù)據(jù)可知圖 9(A)慢中子份額最低, 為7.8%, 說明從準直孔方向的出射中子在無材料慢化吸收的情況下大部分為快中子和中能中子. 在加載屏蔽材料之后探測器測得的中子能譜發(fā)生了明顯衰減, 測得的中子數(shù)在全部能區(qū)均大幅減少, 說明制備的復合屏蔽材料有優(yōu)良的中子屏蔽性能. 此外, 中子份額統(tǒng)計數(shù)據(jù)顯示圖9(C)的慢中子份額(25.48%)要高于圖9(B)的(10.03%). 這是由于更厚的材料使入射中子與原子核的發(fā)生的碰撞次數(shù)更多(例如中子與H和C核的彈性碰撞, 與Gd核的非彈性碰撞), 從材料中的原子核間隙透過的幾率更低, 中子被B和Gd吸收的幾率更大, 這就使得穿過15 cm復合材料的平均中子能量和中子數(shù)量低于穿過4.5 cm的復合材料, 在譜圖上表現(xiàn)為低能區(qū)中子份額上升.

    Fig.9 Neutron energy spectrum on the surface of the detector

    (A) the neutron spectrum without shielding material; (B) the neutron spectrum under the composite with the thickness of 4.5 cm; (C) the neutron spectrum under the composite with the thickness of 15 cm.

    圖10為方形(15 cm×15 cm)的M-nanoGd2O3/B4C/HDPE, nanoGd2O3/B4C/HDPE和HDPE在厚度為1.5 cm時的中子和伽馬射線屏蔽結果. 純HDPE缺少中子吸收劑, 且由輕核元素C和H組成, 因此中子和伽馬射線屏蔽率均低于納米復合材料. 表面改性后的M-nanoGd2O3/B4C/HDPE的中子和伽馬射線屏蔽性能均明顯優(yōu)于未改性的nanoGd2O3/B4C/HDPE, 這是由于M-nanoGd2O3相較nanoGd2O3在HDPE基體內部的界面相容性和分散性得到提升, 從而提升了填料與中子和伽馬射線的作用幾率, 進而導致改性復合材料的綜合屏蔽性能優(yōu)于未改性復合材料.

    Fig.10 Neutron and gamma?ray shielding rates of three kinds of materials

    Scheme 2為中子及伽馬射線與復合材料相互作用的機理圖. 在中子屏蔽過程中, 主要有兩種情況: 一種情況是入射中子與Gd發(fā)生非彈性散射, 使中子能量迅速降低到1~2 MeV[33], 并產生次級伽馬射線. 然后減速的中子與H, C和O核發(fā)生幾次彈性碰撞, 直至變成熱中子. 根據(jù)能量守恒和動能守恒定律, 在彈性碰撞過程中, 中子與H核的碰撞是降低中子能量最有效的途徑. 另一種情況是, 入射中子只與H核發(fā)生多次彈性碰撞, 使其能量迅速降低到熱中子的水平. 熱中子吸收也存在兩種可能, 一種是被B吸收, 另一種是被Gd吸收. 重核Gd通過輻射俘獲反應吸收熱中子[34], 輕核B通過(n,)核反應吸收熱中子, 并產生帶電粒子[35].

    Scheme 2Mechanism diagram of neutron andphoton interacting with the composite material

    圖11是圓形材料和方形材料的中子屏蔽性能的實測數(shù)據(jù)和SuperMC模擬計算結果. 采用式(6)對實驗數(shù)據(jù)和SuperMC模擬計算數(shù)據(jù)進行擬合, 相關數(shù)據(jù)見表2. 在實驗中獲得50%屏蔽率下圓形材料厚度需要達到7 cm, 而方形材料僅需3 cm. 材料面積的增大能夠顯著增加中子屏蔽性能, SuperMC模擬計算結果也得到類似規(guī)律. 這是由于較大的屏蔽面積能夠有效地遮擋出射口前端各個方向的散射中子, 使大部分出射中子從屏蔽材料內部穿過, 減少了材料與屏蔽筒縫隙處散射的中子份額, 從而降低了環(huán)境的中子場強, 使得讀數(shù)結果更加接近材料的真實屏蔽情況. 此外, 無論是方形材料還是圓形材料, 當材料較薄時, 厚度的增加能使中子屏蔽率顯著升高, 但隨著材料增加到一定厚度時, 屏蔽率的增加幅度減小. 例如方形材料厚度從1.5 cm增至3 cm時, 屏蔽率提升了47.0%, 而厚度從9 cm增至10.5 cm時屏蔽率僅提升了6.7%, 這一規(guī)律同樣在SuperMC模擬計算結果中得到驗證. 這是由于隨著材料厚度的增加, 大部分出射中子已被材料屏蔽, 環(huán)境周圍散射中子的份額越來越大, 此時對探測器造成讀數(shù)波動的已經由出射中子變成為散射中子, 在這種情況下可以考慮增加屏蔽材料的屏蔽面積來降低中子的散射程度.

    Table 2 Neutron shielding parameters of real measured data and SuperMC simulated results*

    *is the macroscopic cross section of composites, HVL is the half value layer of composites,2is the correlation coefficient of the fitting curve.

    Fig.11 Real measured data and simulated results of neutron shielding tests

    宏觀總截面是一個中子同單位體積內的原子核發(fā)生核反應的平均幾率大小的量度. 本文制備的方形M-nanoGd2O3/B4C/HDPE在厚度為7.5 cm和11.7 cm時分別達到了80%和90%的中子屏蔽率, 其的實測和模擬計算結果分別為0.332和0.395 cm-1, 是一種高性能的中子屏蔽復合材料.

    圖12是圓形材料和方形材料伽馬射線屏蔽性能的實測數(shù)據(jù)和SuperMC模擬計算結果. 采用式(5)對實驗數(shù)據(jù)和SuperMC模擬計算數(shù)據(jù)進行擬合, 相關數(shù)據(jù)見表3. 在實驗中達到50%伽馬射線屏蔽率需要加載厚度為9.5 cm的圓形材料和7.8 cm的方形材料, SuperMC的模擬計算結果分別為6.5和6.0 cm. 在達到相同屏蔽率時, 伽馬射線實驗中兩種規(guī)格材料厚度差距相較中子實驗小, 這是由于伽馬光子相較中子來說散射程度和透射能力較小, 且本實驗中鉛室的出射孔僅為1 cm, 在很大程度上降低了伽馬射線的散射效果. 因此, 實測數(shù)據(jù)的誤差棒相較中子實驗也更短, 材料面積的增大對光子屏蔽的提升效果低于中子, 這在SuperMC模擬計算結果中也得到證實. 表3為利用式(3)和式(4)計算出的m和HVL等相關屏蔽性能參數(shù), 其中HVL的計算值與實測值很接近, 數(shù)值相差不超過 0.2. 本文制備的方形M-nanoGd2O3/B4C/HDPE材料在厚度為13.5 cm時能達到70%的伽馬射線屏 蔽率.

    Table 3 Gamma-ray shielding parameters of the real measured data and SuperMC simulated results*

    *is the linear attenuation coefficient of composites,mis the mass attenuation coefficient of composites, HVL is the half value layer of composites,2is the correlation coefficient of the fitting curve.

    Fig.12 Real measured data and simulated results of gamma?ray shielding tests

    綜合表2和表3結果可發(fā)現(xiàn)實測結果和模擬計算結果有一定出入, 這是因為SuperMC模型中屏蔽材料的元素分布是均一的, 而在實際情況中材料內部的各種元素不可能達到完全均勻分布. 此外, 軟件建模無法完全還原真實實驗環(huán)境, 這也在一定程度上增加了模擬計算與實測值之間的差距.

    3 結 論

    以HDPE為基體, 改性nanoGd2O3和B4C為增強填料, 通過熱壓成型法制備了一種中子及伽馬射線復合屏蔽材料. 應用偶聯(lián)劑對nanoGd2O3進行了表面包覆改性, 改善了nanoGd2O3在基體內部的界面相容性和分散性. 通過對實驗結果和計算結果的分析, 得出以下結論: (1) 改性nanoGd2O3的引入有利于提高復合材料的熱分解溫度和熔融溫度, 這是由于改性nanoGd2O3與HDPE界面結合力增強, 進一步限制了HDPE鏈的相對運動. (2) 向HDPE基體內摻雜適量nanoGd2O3和B4C提高了材料的拉伸強度, 因為填料能夠阻礙聚乙烯鏈的鏈段運動, 同時對基體內部的裂紋具有釘扎作用. 其中nanoGd2O3的表面改性使得復合材料的拉伸強度進一步提升, 這是由于改性填料界面相容性和分散性提高, 從而減輕了材料內部的應力集中效應. (3) 輻射屏蔽實驗結果證明nanoGd2O3的表面改性有效提升了復合材料的中子和伽馬屏蔽率, 這是由于填料分散性的提高增大了與中子和光子的作用幾率. 此外, 復合材料面積的增大能夠有效提升材料對中子和伽馬射線的屏蔽率. 值得注意的是, 由于伽馬光子的散射程度較中子散射小, 因此材料面積的增大對伽馬射線屏蔽率的提升幅度不如中子屏蔽率的提升顯著. (4) 采用蒙特卡羅方法對中子和伽馬射線屏蔽實驗的模擬計算結果與實測結果的規(guī)律一致. 本文制備的方形M-nanoGd2O3/B4C/HDPE復合材料在11.7 cm厚度時能達到90%的中子屏蔽率, 在13.5 cm厚度時能夠達到70%伽馬射線屏蔽率, 有望成為一種具有良好應用前景的輕質無鉛的中子及伽馬射線復合屏蔽材料.

    [1] Shang Y., Yang G., Su F., Shen C.,, 2020,,147—153

    [2] Alsayed Z, Badawi M. S., Awad R., Thabet A. A.,, 2020,(8), 853—867

    [3] Nambiar S., Yeow J., El?Khatib A. M.,, 2012,(11), 5717—5726

    [4] Soltani Z., Beigzadeh A., Ziaie F., Asadi E.,, 2016,, 182—187

    [5] More C. V., Alsayed Z., Badawi M. S., Thabet A. A., Pawar P. P.,, 2021,(3), 2057—2090

    [6] Cherkashina N. I., Pavlenko V. I., Noskov A. V.,, 2019,, 111—117

    [7] Darwish A. A., Hassan M. H., Mandour M., Maarouf A. A.,, 2018,, 142—147

    [8] Abdelaziz M. M., Gwaily S. E., Makarious A. S., Abdo A. E. S.,, 1995,(2), 235—240

    [9] Akkurt I., Basyigit C., Kilincarslan S.,, 2005,(1), 1—11

    [10] Rammah Y. S., Olarinoye I. O., El?Agawany F. I., El?Adawy A.,, 2021,(2), 2547—2556

    [11] Samia B., Hotc D.,, 2019,(17), 23561—23571

    [12] Lv J. X., Chen J. Y.,1994,, 370—374(呂繼新, 陳建延. 核動力工程, 1994,, 370—374)

    [13] Gu Y. B., Hu Y. H., Du L., Deng S. F., Zhou Y., Ru L. A.,2017,(9), 1670—1677(顧淵博, 扈艷紅, 杜磊, 鄧詩峰, 周燕, 汝凌傲. 高等學校化學學報, 2017,(9), 1670—1677)

    [14] Fereiduni E., Ghasemi A., Elbestawi M., Jadhav S. D.,, 2021,, 1299—1307

    [15] Chen H. S., Wang W. X., Li Y. L., Zhang P., Nie H. H., Wu Q. C.,, 2015,, 23—29

    [16] Zhang P., Jia C., Li J., Wang W. X.,, 2020,, 1280—1293

    [17] ?rim S. G., Wis A. A., Keskin M. A., Baykara O., Ozkoc G., Avci A., Dogru M., Karakoc M.,, 2018,, 434—443

    [18] Baykara O., ?rim S. G., Wis A. A., Keskin M. A., Ozkoc G., Avc A., Dogru M.,, 2020,(11), 2466—2479

    [19] Donzella A., Bonomi G., Giroletti E., Zenoni A.,, 2012,(4), 414—420

    [20] Biswas R., Sahadath H., Mollah A. S., Huq M. F.,, 2016,(1), 26—34

    [21] Al?Buriahi M. S., Sayyed M. I., Al?Hadeethi Y.,2020,(9), 13622—13629

    [22] Eke C., Boztosun I.,2014,(1), 103—108

    [23] Blaauw M.,1996,(3), 431—435

    [24] Guay?Bégin A. A., Chevallier P., Faucher L., Turgeon S., Fortin M. A.,, 2012,(1), 774—782

    [25] Gong L. X., Zhou Z. F., Wang S. M., Wang B.,2013,(3), 1212—1217

    [26] Ni W. S., Wu S. P., Ren Q.,2013,, 272—277

    [27] Zhu Q., Zhou H. L., Song Y. X., Chang Z. D., Li W. J.,2017,(2), 208—215

    [28] Wang H. Y, Yao X. D, Sui G. X, Yin L. M., Wang L. H.,, 2015,(2), 164—170

    [29] Ochiai S., Uehara T., Osamura K.,1986,(3), 1020—1026

    [30] Gao J. M., Yang Y., Lei T., Wang J., Liu J., Zhang J. M.,, 2020,(1), 36—42

    [31] Gu J. W., Zhang Q. Y., Dang J., Zhang J. P., Yang Z. Y.,, 2010,(5), 1030—1034

    [32] Hajek M., Suarez R. C.,2016,(4), 265—268

    [33] Clark J. K., Yannello V., Samarakoon A. M., Ross C., Uible M. C., Garlea V. O., Shatruk M.,2021,(29), 16183—16190

    [34] Tokuuye K., Tokita N., Akine Y., Nakayama H., Sakurai Y., Kobayashi T., Kanda K.,2000,(2), 81—83

    [35] Mastromarco M., Manna A., Aberle O.,2019,(453), 45—59

    Preparation of Modified Gadolinium/Boron/Polyethylene Nanocomposite and Its Radiation Shielding Performance for Neutron and-ray

    ZHAOSheng1,2, HUOZhipeng1, ZHONGGuoqiang1, ZHANGHong1, HULiqun1

    (1,,,230031,;2,230026,)

    The surface of nano-sized gadolinium oxide(nanoGd2O3) was modified by the coupling agent, 3-(Trimethoxysilyl)propyl methacrylate(KH570), and a new type of nanocomposite(modified-gadolinia/boron carbide/polyethylene, M-nanoGd2O3/B4C/HDPE) was prepared by the hot-pressing method. The results of Fourier transform infrared spectroscopy(FTIR), scanning electron microscopy(SEM) and energy dispersive X-ray energy spectroscopy(EDS) show that the nanoGd2O3wassuccessfully modified bycoupling reagent and the surface modification of nanoGd2O3significantly improved the interfacial compatibility and dispersion in the polyethylene matrix. Thermogravimetric analysis(TGA), differential scanning calorimetry(DSC) analysis and mechanical tensile tests showed that the modification of nanoGd2O3improved the thermal stability, tensile strength, elongation at break and Young’s modulus of the composites. The neutron andray shielding performance of the M-nanoGd2O3/B4C/HDPE was studied by real tests and Monte Carlo simulation. The results showed that the modified nanoGd2O3with excellent interfacial compatibility and dispersion could enhance the neutron andray shielding rate effectively. The square shape M-nanoGd2O3/B4C/HDPE material could achieve a neutron shielding rate of 90% at thickness of 11.7 cm and aray shielding rate of 70% at thickness of 13.5 cm.

    Gadolinium oxide; Coupling agent; Surface modification; Polyethylene; Composite; Radiation protection

    O614; TL344

    A

    10.7503/cjcu20220039

    2022-01-15

    2022-03-20.

    霍志鵬, 男, 博士, 副研究員, 主要從事輻射防護與環(huán)境保護及輻射屏蔽材料研究. E-mail: zhipeng.huo@ipp.ac.cn

    國家自然科學基金(批準號: 21777165)、 中國科學院合肥大科學中心“高端用戶培育基金”(批準號: 2020HSC-UE012)、 聚變堆主機關鍵系統(tǒng)綜合研究設施(CRAFT)項目(批準號: 2018-000052-73-01-001228)和能源研究院(安徽省能源實驗室)項目(批準號: 21KZL401, 21KHH105, 21KZS205)資助.

    Supported by the National Natural Science Foundation of China(No.21777165), the Users with Excellence Program of Hefei Science Center of Chinese Academy of Sciences(No.2020HSC-UE012), the Comprehensive Research Facility for Fusion Technology Program of China(No. 2018-000052-73-01-001228) and the Institute of Energy(Hefei Comprehensive National Science Center), China(Nos.21KZL401, 21KHH105, 21KZS205).

    (Ed.: W, K, M)

    猜你喜歡
    伽馬射線模擬計算中子
    我國觀測到迄今最亮伽馬射線暴
    科學之友(2022年11期)2022-11-03 10:45:18
    R1234ze PVTx熱物性模擬計算
    能源工程(2022年1期)2022-03-29 01:06:26
    3D打印抗中子輻照鋼研究取得新進展
    7個有關伽馬射線的驚人事實
    科學之謎(2017年4期)2017-05-13 00:25:46
    擠出發(fā)泡片材褶皺分析及模擬計算
    中國塑料(2016年9期)2016-06-13 03:18:54
    基于PLC控制的中子束窗更換維護系統(tǒng)開發(fā)與研究
    DORT 程序進行RPV 中子注量率計算的可靠性驗證
    實際發(fā)射工況下底排藥柱結構完整性的模擬計算
    火炸藥學報(2014年3期)2014-03-20 13:17:43
    中子深度定量分析的相對分析法
    計算物理(2014年2期)2014-03-11 17:01:27
    伽馬射線波段的地球與天空
    飛碟探索(2014年1期)2014-03-04 09:08:53
    亚洲精华国产精华精| 麻豆久久精品国产亚洲av| 男人舔奶头视频| 久久这里只有精品中国| 91av网站免费观看| 99久久国产精品久久久| 欧美最黄视频在线播放免费| 中文亚洲av片在线观看爽| 色噜噜av男人的天堂激情| 亚洲无线在线观看| 热99re8久久精品国产| 精品一区二区三区av网在线观看| 老司机在亚洲福利影院| 亚洲成人精品中文字幕电影| 久久久久久大精品| 无人区码免费观看不卡| 真人做人爱边吃奶动态| 老汉色∧v一级毛片| 免费av毛片视频| 国产成人精品久久二区二区91| 一进一出抽搐动态| 看黄色毛片网站| 一边摸一边做爽爽视频免费| 午夜免费观看网址| 18禁裸乳无遮挡免费网站照片| 少妇的丰满在线观看| 国产激情欧美一区二区| 777久久人妻少妇嫩草av网站| 美女午夜性视频免费| 一二三四社区在线视频社区8| 又大又爽又粗| 国产激情久久老熟女| 一夜夜www| 欧美日韩中文字幕国产精品一区二区三区| 99精品欧美一区二区三区四区| 国产主播在线观看一区二区| 国产伦在线观看视频一区| 国产免费男女视频| 日日夜夜操网爽| 脱女人内裤的视频| 国产激情偷乱视频一区二区| aaaaa片日本免费| 欧美在线黄色| 国产亚洲av高清不卡| 午夜精品一区二区三区免费看| av福利片在线| 久久久久国内视频| 我要搜黄色片| 一二三四社区在线视频社区8| 精品电影一区二区在线| 日本成人三级电影网站| 国产精品久久电影中文字幕| 亚洲一区高清亚洲精品| 一进一出抽搐动态| 制服丝袜大香蕉在线| 一a级毛片在线观看| 久久这里只有精品中国| 免费搜索国产男女视频| 一卡2卡三卡四卡精品乱码亚洲| 韩国av一区二区三区四区| 999久久久国产精品视频| 变态另类丝袜制服| 国产精华一区二区三区| 精品一区二区三区av网在线观看| 老熟妇仑乱视频hdxx| 老汉色av国产亚洲站长工具| 麻豆成人午夜福利视频| 啦啦啦观看免费观看视频高清| 亚洲人成网站在线播放欧美日韩| 99re在线观看精品视频| 国产高清视频在线观看网站| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人欧美在线观看| 国产伦一二天堂av在线观看| 亚洲成人中文字幕在线播放| 男女视频在线观看网站免费 | 不卡一级毛片| 变态另类丝袜制服| 亚洲人成伊人成综合网2020| 亚洲一区二区三区不卡视频| 午夜福利18| 精品国产乱码久久久久久男人| 国产亚洲欧美98| 中文资源天堂在线| 国产又色又爽无遮挡免费看| 色精品久久人妻99蜜桃| 国内久久婷婷六月综合欲色啪| 欧美精品亚洲一区二区| 久久中文看片网| 色精品久久人妻99蜜桃| 亚洲人成网站高清观看| 欧美一区二区精品小视频在线| 日韩三级视频一区二区三区| 桃色一区二区三区在线观看| 国产免费男女视频| 他把我摸到了高潮在线观看| 一本一本综合久久| 欧美高清成人免费视频www| 国产三级中文精品| www.999成人在线观看| 欧美日韩一级在线毛片| 亚洲av中文字字幕乱码综合| 伊人久久大香线蕉亚洲五| 日本精品一区二区三区蜜桃| 国产麻豆成人av免费视频| 成人欧美大片| 欧美日韩黄片免| 黄片小视频在线播放| 桃红色精品国产亚洲av| 脱女人内裤的视频| 无遮挡黄片免费观看| 久久中文字幕一级| 免费看日本二区| 又黄又爽又免费观看的视频| 国内少妇人妻偷人精品xxx网站 | av片东京热男人的天堂| 无限看片的www在线观看| 欧美av亚洲av综合av国产av| 成人午夜高清在线视频| 伦理电影免费视频| 每晚都被弄得嗷嗷叫到高潮| 在线观看美女被高潮喷水网站 | www.自偷自拍.com| 亚洲欧美一区二区三区黑人| 高清在线国产一区| 999久久久国产精品视频| 久久中文字幕一级| 正在播放国产对白刺激| 18禁观看日本| 一夜夜www| www.熟女人妻精品国产| 亚洲 国产 在线| 我要搜黄色片| 观看免费一级毛片| 老熟妇仑乱视频hdxx| 叶爱在线成人免费视频播放| 日本在线视频免费播放| 18禁观看日本| 色播亚洲综合网| 日本一本二区三区精品| 亚洲精品在线观看二区| 精品免费久久久久久久清纯| 久久久久精品国产欧美久久久| 在线观看www视频免费| 九色成人免费人妻av| 三级毛片av免费| 韩国av一区二区三区四区| 午夜a级毛片| 香蕉国产在线看| 精品一区二区三区av网在线观看| 成人三级黄色视频| 青草久久国产| 国产成+人综合+亚洲专区| 国产乱人伦免费视频| 久久精品人妻少妇| 色播亚洲综合网| 国内精品久久久久精免费| av中文乱码字幕在线| 免费看美女性在线毛片视频| 国产成人av激情在线播放| 精品电影一区二区在线| 哪里可以看免费的av片| 欧美成人免费av一区二区三区| av福利片在线观看| 露出奶头的视频| 欧美一级毛片孕妇| 两个人的视频大全免费| 午夜亚洲福利在线播放| 亚洲av成人精品一区久久| 精品人妻1区二区| 久久精品国产亚洲av高清一级| www.精华液| 国产视频一区二区在线看| 欧美黄色片欧美黄色片| 亚洲欧美激情综合另类| 黄色成人免费大全| 久久香蕉国产精品| 色综合欧美亚洲国产小说| 在线a可以看的网站| 国产黄a三级三级三级人| 中文字幕精品亚洲无线码一区| 国产成+人综合+亚洲专区| 欧美午夜高清在线| 国产69精品久久久久777片 | 香蕉av资源在线| 18禁裸乳无遮挡免费网站照片| 人人妻人人看人人澡| www日本黄色视频网| 久久精品成人免费网站| 亚洲成av人片在线播放无| 国产av一区二区精品久久| 91国产中文字幕| 波多野结衣巨乳人妻| 国产一区在线观看成人免费| 亚洲一卡2卡3卡4卡5卡精品中文| 日本成人三级电影网站| а√天堂www在线а√下载| 亚洲av成人精品一区久久| 亚洲美女黄片视频| 成年人黄色毛片网站| 一本大道久久a久久精品| 村上凉子中文字幕在线| 欧美日韩福利视频一区二区| 日日摸夜夜添夜夜添小说| 国产亚洲精品第一综合不卡| 淫秽高清视频在线观看| 日本在线视频免费播放| 成人三级黄色视频| 中文字幕精品亚洲无线码一区| 久热爱精品视频在线9| 久久香蕉国产精品| 日本一本二区三区精品| 啦啦啦免费观看视频1| 夜夜爽天天搞| 99热这里只有是精品50| 色噜噜av男人的天堂激情| 午夜精品久久久久久毛片777| 看黄色毛片网站| 麻豆久久精品国产亚洲av| 日韩 欧美 亚洲 中文字幕| 国产精品影院久久| 日韩 欧美 亚洲 中文字幕| 欧美日韩乱码在线| 国产69精品久久久久777片 | 成人高潮视频无遮挡免费网站| 亚洲av成人不卡在线观看播放网| 欧美黑人欧美精品刺激| tocl精华| 精品久久蜜臀av无| 欧美久久黑人一区二区| 正在播放国产对白刺激| 小说图片视频综合网站| 天天躁夜夜躁狠狠躁躁| 国产欧美日韩精品亚洲av| 亚洲色图 男人天堂 中文字幕| 少妇的丰满在线观看| 99国产精品一区二区三区| 麻豆成人av在线观看| 久久精品人妻少妇| 黄色视频不卡| 中文字幕精品亚洲无线码一区| 免费观看精品视频网站| 999久久久国产精品视频| 国产主播在线观看一区二区| 国产成人影院久久av| 国产激情偷乱视频一区二区| 亚洲激情在线av| 俺也久久电影网| 免费在线观看视频国产中文字幕亚洲| 99在线人妻在线中文字幕| 午夜成年电影在线免费观看| 国产av在哪里看| 久久午夜亚洲精品久久| 白带黄色成豆腐渣| 叶爱在线成人免费视频播放| 看黄色毛片网站| 99热只有精品国产| 欧美黄色淫秽网站| 婷婷精品国产亚洲av在线| 成人手机av| 亚洲av五月六月丁香网| 国产伦人伦偷精品视频| 国产午夜精品论理片| 99精品欧美一区二区三区四区| 亚洲avbb在线观看| 欧美在线黄色| 中出人妻视频一区二区| 女人高潮潮喷娇喘18禁视频| 午夜老司机福利片| 亚洲中文字幕一区二区三区有码在线看 | 91字幕亚洲| 午夜福利在线观看吧| 国产熟女xx| 国产精品一区二区三区四区久久| 高清在线国产一区| 亚洲国产精品999在线| www.www免费av| 久久久精品大字幕| 不卡一级毛片| 88av欧美| 亚洲精品久久国产高清桃花| 久久久久精品国产欧美久久久| 又紧又爽又黄一区二区| 欧美黄色片欧美黄色片| 亚洲五月婷婷丁香| 99精品在免费线老司机午夜| 日韩欧美 国产精品| 亚洲成人国产一区在线观看| 成人手机av| 成人一区二区视频在线观看| 18禁观看日本| 51午夜福利影视在线观看| 午夜久久久久精精品| 国产一级毛片七仙女欲春2| 亚洲中文av在线| 亚洲第一欧美日韩一区二区三区| 国产精品精品国产色婷婷| 很黄的视频免费| 777久久人妻少妇嫩草av网站| 亚洲精品粉嫩美女一区| 久久草成人影院| 日韩国内少妇激情av| 精华霜和精华液先用哪个| 欧美成人一区二区免费高清观看 | 免费观看精品视频网站| 嫩草影视91久久| 免费看美女性在线毛片视频| 麻豆成人午夜福利视频| 1024香蕉在线观看| xxxwww97欧美| 女同久久另类99精品国产91| 九色成人免费人妻av| 99久久综合精品五月天人人| 久久精品91蜜桃| 日本免费a在线| 国产精品 欧美亚洲| 国产男靠女视频免费网站| 成人手机av| 白带黄色成豆腐渣| 欧美日韩中文字幕国产精品一区二区三区| 亚洲精品国产一区二区精华液| 两个人的视频大全免费| 久久久精品欧美日韩精品| 美女扒开内裤让男人捅视频| 国产一级毛片七仙女欲春2| 久久中文字幕一级| 国产精品久久电影中文字幕| www.自偷自拍.com| 午夜福利视频1000在线观看| 搞女人的毛片| 久久久国产精品麻豆| e午夜精品久久久久久久| 又黄又爽又免费观看的视频| 亚洲国产欧美人成| 欧美日韩亚洲国产一区二区在线观看| 三级男女做爰猛烈吃奶摸视频| 亚洲熟妇中文字幕五十中出| 欧美中文日本在线观看视频| 国产一区二区三区视频了| 亚洲一码二码三码区别大吗| 天堂av国产一区二区熟女人妻 | 一二三四社区在线视频社区8| 女警被强在线播放| 特级一级黄色大片| 国产主播在线观看一区二区| 俺也久久电影网| 九色国产91popny在线| 女同久久另类99精品国产91| 久久中文看片网| 久久精品成人免费网站| 亚洲天堂国产精品一区在线| 又粗又爽又猛毛片免费看| 国产精品免费一区二区三区在线| 黄片大片在线免费观看| 欧美另类亚洲清纯唯美| 免费观看精品视频网站| 久久人妻av系列| 亚洲精品在线观看二区| 国产精品自产拍在线观看55亚洲| 国产69精品久久久久777片 | 在线观看免费日韩欧美大片| 欧美中文日本在线观看视频| 亚洲中文字幕日韩| 亚洲精品久久成人aⅴ小说| 制服诱惑二区| 精品高清国产在线一区| 巨乳人妻的诱惑在线观看| 久久精品成人免费网站| 欧美一级a爱片免费观看看 | xxxwww97欧美| 亚洲熟女毛片儿| 免费电影在线观看免费观看| 女警被强在线播放| 麻豆成人av在线观看| 亚洲精品粉嫩美女一区| 国产69精品久久久久777片 | 99久久久亚洲精品蜜臀av| 一二三四社区在线视频社区8| 欧美色视频一区免费| 亚洲一区二区三区不卡视频| 国产精品av视频在线免费观看| 99热6这里只有精品| 中文字幕久久专区| 男女下面进入的视频免费午夜| 在线看三级毛片| 999久久久国产精品视频| 俺也久久电影网| 国产熟女xx| 国产精品九九99| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看66精品国产| 亚洲国产欧美人成| 国产久久久一区二区三区| 麻豆久久精品国产亚洲av| a级毛片在线看网站| 看黄色毛片网站| 女生性感内裤真人,穿戴方法视频| 麻豆国产av国片精品| 久久精品国产综合久久久| 天堂动漫精品| 午夜久久久久精精品| 两个人免费观看高清视频| 夜夜夜夜夜久久久久| 精品国产亚洲在线| 国产成人精品久久二区二区91| 一区二区三区国产精品乱码| 成人三级黄色视频| 久久久久国内视频| 国产在线精品亚洲第一网站| 日本五十路高清| 国产亚洲精品久久久久5区| 国产日本99.免费观看| 国产1区2区3区精品| 国产不卡一卡二| 久9热在线精品视频| 亚洲av中文字字幕乱码综合| 亚洲男人天堂网一区| 欧美一级毛片孕妇| 精品人妻1区二区| 午夜精品久久久久久毛片777| 麻豆成人av在线观看| 岛国在线免费视频观看| 99国产精品一区二区三区| 久久中文字幕一级| 天堂av国产一区二区熟女人妻 | 亚洲国产精品成人综合色| 欧美乱色亚洲激情| 亚洲av美国av| 亚洲国产精品久久男人天堂| 国产视频一区二区在线看| 黄色视频,在线免费观看| 99热这里只有是精品50| 99re在线观看精品视频| 亚洲人与动物交配视频| 久久久久久免费高清国产稀缺| 亚洲九九香蕉| 亚洲成av人片在线播放无| 免费av毛片视频| 老司机午夜福利在线观看视频| 亚洲熟女毛片儿| 亚洲精品一卡2卡三卡4卡5卡| 一夜夜www| 级片在线观看| 国产成年人精品一区二区| 久久午夜综合久久蜜桃| 日本精品一区二区三区蜜桃| 精品国内亚洲2022精品成人| 亚洲成人国产一区在线观看| 国产又黄又爽又无遮挡在线| 十八禁网站免费在线| 精品国产乱子伦一区二区三区| 国产又色又爽无遮挡免费看| 日本在线视频免费播放| 18禁观看日本| 日本 欧美在线| 国产伦在线观看视频一区| 国产不卡一卡二| 亚洲人成伊人成综合网2020| 免费在线观看视频国产中文字幕亚洲| 午夜免费观看网址| 久久久国产精品麻豆| 欧美成人一区二区免费高清观看 | 黄色成人免费大全| 51午夜福利影视在线观看| 99久久99久久久精品蜜桃| 最近最新中文字幕大全免费视频| 在线观看www视频免费| 美女 人体艺术 gogo| 日韩欧美精品v在线| 精品欧美一区二区三区在线| 草草在线视频免费看| xxx96com| 天堂√8在线中文| 亚洲午夜理论影院| 亚洲五月天丁香| 一区二区三区激情视频| 亚洲欧美激情综合另类| 黄色女人牲交| 看片在线看免费视频| 精品久久久久久久毛片微露脸| 国产精品国产高清国产av| 又大又爽又粗| 国产精品影院久久| 成人18禁在线播放| 少妇被粗大的猛进出69影院| 国产精品一区二区免费欧美| 久久中文字幕一级| 国产亚洲av高清不卡| 91国产中文字幕| 亚洲国产高清在线一区二区三| 在线观看免费日韩欧美大片| av中文乱码字幕在线| e午夜精品久久久久久久| 99国产极品粉嫩在线观看| 亚洲人成伊人成综合网2020| 99久久精品热视频| 人妻丰满熟妇av一区二区三区| 久久精品影院6| 一区二区三区国产精品乱码| 男女视频在线观看网站免费 | 热99re8久久精品国产| 美女午夜性视频免费| 中文字幕最新亚洲高清| 欧美日韩中文字幕国产精品一区二区三区| 手机成人av网站| 久久久国产成人免费| 国产在线观看jvid| 午夜精品一区二区三区免费看| 亚洲人成77777在线视频| 亚洲va日本ⅴa欧美va伊人久久| а√天堂www在线а√下载| 国产午夜福利久久久久久| 亚洲免费av在线视频| 两个人免费观看高清视频| 久久 成人 亚洲| 亚洲狠狠婷婷综合久久图片| 日韩精品中文字幕看吧| 麻豆成人午夜福利视频| 亚洲国产看品久久| 亚洲av成人精品一区久久| 国内毛片毛片毛片毛片毛片| 日本黄色视频三级网站网址| 欧美黄色淫秽网站| 欧美黑人欧美精品刺激| 老司机深夜福利视频在线观看| 久久性视频一级片| 国产亚洲精品第一综合不卡| 黑人欧美特级aaaaaa片| 亚洲国产欧美人成| 国产精品1区2区在线观看.| 久久午夜亚洲精品久久| 一级a爱片免费观看的视频| 正在播放国产对白刺激| 国产三级黄色录像| 亚洲va日本ⅴa欧美va伊人久久| 香蕉国产在线看| 国产精品久久久av美女十八| 中文字幕久久专区| 99国产精品一区二区蜜桃av| 国产一级毛片七仙女欲春2| 成年女人毛片免费观看观看9| 少妇被粗大的猛进出69影院| 亚洲国产日韩欧美精品在线观看 | 黄色视频,在线免费观看| 女同久久另类99精品国产91| 国产亚洲欧美98| 岛国视频午夜一区免费看| 熟妇人妻久久中文字幕3abv| 一本一本综合久久| 无遮挡黄片免费观看| xxxwww97欧美| 亚洲精品色激情综合| 韩国av一区二区三区四区| 男插女下体视频免费在线播放| 亚洲自拍偷在线| 精品第一国产精品| 听说在线观看完整版免费高清| 久久国产精品影院| 无人区码免费观看不卡| 91麻豆精品激情在线观看国产| 亚洲一区二区三区不卡视频| 亚洲色图 男人天堂 中文字幕| 国产一区二区在线观看日韩 | 小说图片视频综合网站| av片东京热男人的天堂| 宅男免费午夜| 国产伦一二天堂av在线观看| 欧美色视频一区免费| 一区福利在线观看| 99精品在免费线老司机午夜| av超薄肉色丝袜交足视频| 婷婷精品国产亚洲av在线| xxx96com| 久久久久久久精品吃奶| 男人舔奶头视频| 久久国产精品人妻蜜桃| 两个人视频免费观看高清| 三级国产精品欧美在线观看 | 变态另类成人亚洲欧美熟女| 亚洲成a人片在线一区二区| 国产欧美日韩一区二区三| 欧美日本视频| 久久性视频一级片| av在线播放免费不卡| 国产精品野战在线观看| 这个男人来自地球电影免费观看| 欧美av亚洲av综合av国产av| 成人欧美大片| 久久久久国内视频| 搡老岳熟女国产| 99国产精品99久久久久| 黄片大片在线免费观看| 日韩 欧美 亚洲 中文字幕| 人人妻人人看人人澡| 毛片女人毛片| 老熟妇仑乱视频hdxx| 亚洲第一欧美日韩一区二区三区| 久久精品91蜜桃| 黄色女人牲交| 欧美另类亚洲清纯唯美| 亚洲18禁久久av| 中文字幕久久专区| 久久久久久久午夜电影| 欧美日韩亚洲国产一区二区在线观看| 免费在线观看完整版高清| 免费在线观看日本一区| 免费观看人在逋| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产精品成人综合色| 日本黄色视频三级网站网址| 国产97色在线日韩免费| 久久精品aⅴ一区二区三区四区| 成人欧美大片| 亚洲九九香蕉| 国产精品久久久久久久电影 |