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

    基于DSC的預(yù)浸料固化工藝評價和優(yōu)化

    2017-01-17 03:25:08陶家祥
    高科技纖維與應(yīng)用 2016年1期
    關(guān)鍵詞:環(huán)氧樹脂樹脂動力學(xué)

    陶家祥,李 煒

    (東華大學(xué) 紡織學(xué)院 紡織面料技術(shù)教育部重點實驗室,上海 201620)

    基于DSC的預(yù)浸料固化工藝評價和優(yōu)化

    陶家祥,李 煒

    (東華大學(xué) 紡織學(xué)院 紡織面料技術(shù)教育部重點實驗室,上海 201620)

    采用差示掃描量熱法(DSC)測量了預(yù)浸料固化反應(yīng)熱過程。并對預(yù)浸料固化的影響因素、固化參數(shù)的確定,以及通過DSC確定預(yù)浸料的固化參數(shù)對預(yù)浸料固化工藝進(jìn)行評價和通過預(yù)浸料固化動力學(xué)對預(yù)浸料的固化工藝進(jìn)行優(yōu)化。

    預(yù)浸料;差示掃描量熱法;固化工藝;評價;優(yōu)化

    0 引言

    預(yù)浸料固化工藝對復(fù)合材料性能至關(guān)重要,其工藝主要包括固化溫度、固化時間和固化壓力3 個重要參數(shù)[1]。預(yù)浸料固化工藝可以根據(jù)固化度、玻璃化轉(zhuǎn)變溫度(Tg)等指標(biāo)來進(jìn)行評定,幫助預(yù)測復(fù)合材料力學(xué)性能,減少測試工作量。同時,利用固化動力學(xué)模型對固化工藝中的溫度和時間進(jìn)行轉(zhuǎn)換,可以很好地反映樹脂固化過程,根據(jù)實際需求選擇合適的固化溫度和固化時間。通過DSC對預(yù)浸料進(jìn)行分析測試,了解樹脂化學(xué)反應(yīng)特點,測定固化度、Tg等指標(biāo)對預(yù)浸料固化工藝進(jìn)行評定,并通過固化動力學(xué)模型來制定和優(yōu)化固化工藝的方法。

    1 預(yù)浸料固化工藝參數(shù)

    差示掃描量熱法(Differential Scanning Calorimetry,簡稱DSC),是指在程序溫度控制下,測量試樣與參比物之間建立零溫度差所需單位時間內(nèi)能量差(或功率差)隨溫度變化的一種技術(shù)??梢詼y定一定升溫速率條件下樹脂發(fā)生固化反應(yīng)的溫度及固化時所放出的熱量,是定量測定樹脂固化反應(yīng)程度的有效方法。根據(jù)DSC測定預(yù)浸料固化溫度和固化時間。

    1.1 影響DSC測試結(jié)果的因素

    DSC作為一種動態(tài)的量熱技術(shù),有很多因素會影響試驗結(jié)果,一般主要受升溫速度、試樣大小和試樣均勻程度的影響[2]。毛如增[3]等人在測定環(huán)氧樹脂固化反應(yīng)溫度和反應(yīng)熱研究中研究了升溫速率和樣品量對DSC測試的環(huán)氧樹脂的固化反應(yīng)溫度和反應(yīng)熱測量結(jié)果的影響。

    1.1.1 升溫速度

    DSC測試中,升溫速度愈快,反應(yīng)愈快,固化反應(yīng)不完全,造成了固化放熱滯后現(xiàn)象。如圖1所示,升溫速度快的DSC曲線的放熱峰向高溫方向移動,峰形變得尖而陡,則測定的預(yù)浸料反應(yīng)熱偏小。應(yīng)盡量采用較慢的升溫速度,以確保固化反應(yīng)完全,但升溫速度越慢,需要的時間越長,一般選定5或10 K/min[2~3]。

    表1 樣品量對環(huán)氧樹脂固化反應(yīng)熱的影響[3]Table 1 Effect of sample volume on the epoxy resin curing reaction heat[3]

    圖1 升溫速率對DSC曲線的影響Fig.1 Effect of heating rate on the DSC curve

    圖2 DSC測定特征溫度Fig.2 The characteristic temperature determined by DSC

    圖3 不同升溫速度下的特征溫度Fig.3 The characteristic temperature at different heating rate

    1.1.2 樣品

    測試樣品量和樣品中樹脂質(zhì)量分?jǐn)?shù)的均勻性會直接影響DSC數(shù)據(jù)的準(zhǔn)確性。樣品量越多,樣品內(nèi)部傳熱越慢,形成溫度梯度越大,峰形擴(kuò)張,峰頂溫度移向高溫[3]。如表1所示,通過樣品量取5 mg與10 mg的DSC對比,10 mg的方差小,實驗精度高,所以一般取樣為10 mg[3]。同時,為了樣品有良好的熱接觸和樹脂質(zhì)量分?jǐn)?shù),則一般采用疊加式取樣。此外,在制備樣品中,樣品受污染也會使得有多余的放熱峰等出現(xiàn),在取樣時要帶一次性手套,防止污染樣品。

    1.2 固化溫度

    預(yù)浸料若采用直接升溫至固化溫度進(jìn)行固化的固化方式,則反應(yīng)活性太大,容易集聚放熱,造成復(fù)合材料成型品質(zhì)不良,所以,預(yù)浸料一般采用階梯式固化工藝。

    預(yù)浸料固化工藝是根據(jù)樹脂性能來確定的。如圖2所示,通過DSC對預(yù)浸料進(jìn)行測試,根據(jù)放熱峰可以得到3 個特征溫度:起始溫度(Ti)、峰值溫度(Tp)和峰末溫度(Tf)。圖3是試樣在5、10和15 ℃/min升溫速度下的DSC曲線,隨著升溫速度的提高,放熱峰向右移動,3 個特征溫度均會增大,這樣就得到了在不同升溫速度(β)下的特征溫度(T),通過T-β圖外推法,可以得到β=0時3 個特征溫度(Ti、Tp和Tf),3 個溫度分別對應(yīng)為近似凝膠溫度Tgl、固化溫度Tcure和后處理溫度Ttreat,這3 個溫度即可作為固化工藝階梯溫

    度[4~15]。

    外推法得到的溫度可以作為預(yù)浸料固化工藝溫度。如果僅對純樹脂進(jìn)行DSC研究,則在制定預(yù)浸料固化工藝時必須考慮增強(qiáng)纖維對實際工藝產(chǎn)生的影響[6,11~12]。

    1.3 固化時間

    預(yù)浸料固化時間決定了預(yù)浸料成型時的固化程度,固化時間越長,預(yù)浸料固化程度越高。固化時間根據(jù)恒溫時完全固化所需時間來確定,即在確定固化溫度后,通過恒溫DSC測試預(yù)浸料在近似凝膠溫度Tgl、固化溫度Tcure、后處理溫度Ttreat等不同溫度間的放熱峰距離來確定固化時間,如圖4是根據(jù)DSC恒溫曲線獲得固化時間[14~15]。

    一般來說,固化時間過長,不僅會提高加工成本,而且也增加了成型產(chǎn)品中的熱應(yīng)力,反而會降低產(chǎn)品性能。因此,需要選擇合適的固化時間。在預(yù)固化溫度期間,一般達(dá)到40%~50%的固化度[4],固化溫度期間達(dá)到60%~80%的固化度[5]。

    2 固化工藝評價

    通常復(fù)合材料力學(xué)性能隨固化度增加而提高,Tg隨固化度提高而上升,并且當(dāng)固化度達(dá)到90%時,其復(fù)合材料的熱力學(xué)性能隨著固化度的增大而變化緩慢,因此,可以通過DSC測定固化度和Tg對預(yù)浸料固化工藝進(jìn)行評定[14~25]。

    固化度α是熱固性聚合物材料一個很重要的指標(biāo),用DSC可以很方便的根據(jù)其反應(yīng)熱進(jìn)行測量,計算公式如式⑴

    其中,ΔH0為樹脂體系進(jìn)行完全固化時所放出的總熱量(J/g),ΔHR為固化后剩余反應(yīng)熱(J/g)。

    圖5為通過DSC測定Tg。Tg沒有固定值,會隨著預(yù)浸料試樣重量、升溫速率和材料室溫放置時間等條件而變化。在說明Tg時需要說明測試方法和條件。測定時,一般選取15 mg預(yù)浸料樣品,升溫速度為10 ℃/min。通常,Tg值越高,預(yù)浸料具有較高的交聯(lián)密度和分子網(wǎng)絡(luò)結(jié)構(gòu),說明該預(yù)浸料固化越充分,耐熱性能越好[22~25]。

    圖4 由DSC曲線獲得固化時間Fig.4 The curing time obtained by DSC curve

    圖5 DSC測定TgFig.5 The glass transition temperature determined by DSC

    3 固化工藝優(yōu)化

    預(yù)浸料固化反應(yīng)動力學(xué)可以用數(shù)學(xué)形式得出在某一時刻或具體溫度條件下,己轉(zhuǎn)化物質(zhì)量與轉(zhuǎn)化速率的關(guān)系,利用DSC曲線進(jìn)行動力學(xué)分析,獲得固化反應(yīng)模型,優(yōu)化預(yù)浸料固化工藝。

    利用DSC對固化反應(yīng)過程進(jìn)行研究,一般分為非等溫和等溫兩種方法[27~33]。等溫方法是指預(yù)浸料在不同固化溫度下做等溫DSC曲線,一般用于反應(yīng)復(fù)雜時判定預(yù)浸料反應(yīng)類型;非等溫方法是在不同升溫速度下獲得的預(yù)浸料DSC曲線。非等溫法較之等溫法有許多優(yōu)點,如所需試樣量少;能在反應(yīng)開始到結(jié)束的整個溫度區(qū)間內(nèi)連續(xù)計算動力學(xué)參數(shù);節(jié)省時間并可同時獲得多種信息;可以克服等溫 DSC 法升溫過程帶來的誤差等[33]。

    根據(jù)預(yù)浸料中樹脂固化機(jī)制不同,有不同的固化動力學(xué)模型??煞譃閚級動力學(xué)模型、自催化模型和Kamal模型[34]。這3 種模型都是采用一些經(jīng)驗性的模型方程作為固化反應(yīng)動力學(xué)的基礎(chǔ)方程,通過數(shù)學(xué)擬合求取模型方程中的各個參數(shù)值。通常預(yù)浸料固化都是其中樹脂經(jīng)過軟化、凝膠、固化等階段,其反應(yīng)速度呈現(xiàn)出先增大再減小的趨勢,而n級動力學(xué)模型是假設(shè)一開始反應(yīng)速度最大,固化速度逐漸減少,符合反應(yīng)后期的趨勢,符合固化度范圍大致在0.7~0.98;自催化模型則強(qiáng)調(diào)反應(yīng)速度有個過渡期,是先增大后減少的過程,反應(yīng)速度最大發(fā)生在一段時間后,較符合反應(yīng)中期即固化度在0.2~0.7范圍內(nèi),一般最大速度發(fā)生在固化度0.3~0.4期間內(nèi);而Kamal模型即復(fù)合模型是將n級動力學(xué)模型和自催化模型結(jié)合[29~40]。

    有時,模型并不能準(zhǔn)確表達(dá)出預(yù)浸料中樹脂的固化機(jī)制,特別是在預(yù)浸料中增強(qiáng)材料以及其表面的形狀和漿料基團(tuán)會對樹脂的反應(yīng)產(chǎn)生傳熱、阻位等影響,則需引入新的參數(shù)或模型對固化動力學(xué)模型加以修正,準(zhǔn)確表達(dá)預(yù)浸料在固化時固化度、溫度和時間的關(guān)系,從而對預(yù)浸料固化工藝做出最正確的優(yōu)化。

    借用固化反應(yīng)模型可以預(yù)測轉(zhuǎn)化率與溫度和時間的關(guān)系,一定時間內(nèi)到達(dá)特定轉(zhuǎn)化溫度所需時間,在不同速度下的DSC曲線模擬,從而對整個固化反應(yīng)過程進(jìn)行預(yù)估,而后根據(jù)所需要的預(yù)浸料固化后產(chǎn)品性能對預(yù)浸料整個固化工藝進(jìn)行修改與優(yōu)化,達(dá)到預(yù)期產(chǎn)品性能目的[40~41]。

    4 結(jié)論

    通過DSC測定,可以對預(yù)浸料固化特征進(jìn)行深入的研究,獲得預(yù)浸料固化工藝中的溫度和時間參數(shù)。通過DSC中預(yù)浸料固化放熱積分曲線獲得預(yù)浸料的固化度和Tg,可以對預(yù)浸料的固化工藝進(jìn)行評價。目前,對于預(yù)浸料固化動力學(xué)研究基本是沿用樹脂的動力學(xué)模型,雖然纖維在復(fù)合材料當(dāng)中對固化反應(yīng)機(jī)理沒有影響,但是纖維的存在使得預(yù)浸料固化變成非均體,會對固化時的樹脂流動和擴(kuò)散受到影響,所以研究預(yù)浸料固化時還需要引進(jìn)一些參數(shù),將纖維的作用體現(xiàn)到實際的固化中,從而更好地得到預(yù)浸料固化模型方程,控制和優(yōu)化溫度及時間等固化工藝參數(shù)。

    [1]張鳳翻. 復(fù)合材料用預(yù)浸料[J]. 高科技纖維與應(yīng)用, 1999, 24(5): 29-32.

    [2]高家武, 周福珍, 劉士昕, 等. 高分子材料熱分析曲線集[M]. 北京: 科學(xué)出版社, 1990: 19-20.

    [3]毛如增, 冀克儉, 張銀生, 等. DSC法測定環(huán)氧樹脂固化反應(yīng)溫度和反應(yīng)熱[J]. 工程塑料應(yīng)用, 2002, 30(11): 36-39.

    [4]謝旻, 顧軼卓, 李敏, 等. 碳纖維/雙馬樹脂預(yù)浸料固化過程動態(tài)力學(xué)性能[J]. 復(fù)合材料學(xué)報, 2010, (4): 52-58.

    [5]Makoto Okamoto, Shuntaro Kuriyama. Method and apparatus for molding thermosetting composite material: US, 20050140045[P]. 2005-06-30.

    [6]宋偉, 齊俊偉, 肖軍, 等. 碳纖維/環(huán)氧預(yù)浸料的熱固化工藝研究[J]. 玻璃鋼/復(fù)合材料, 2011, (1): 48-51.

    [7]紀(jì)雙英, 張元晉, 包建文, 等. 真空袋成型用環(huán)氧樹脂固化體系的研究[C]//中國航空學(xué)會. 第17屆全國復(fù)合材料學(xué)術(shù)會議(復(fù)合材料及其原材料分論壇)論文集. 2012.

    [8]歐秋仁, 宮元勛, 趙宏杰, 等. 改性氰酸酯及其復(fù)合材料的性能研究[C]//中國航空學(xué)會. 第17屆全國復(fù)合材料學(xué)術(shù)會議(智能與功能復(fù)合材料分論壇)論文集. 2012.

    [9]熊艷麗, 凌輝, 蔣文革, 等. 結(jié)構(gòu)復(fù)合材料用耐高溫環(huán)氧樹脂體系[J]. 宇航材料工藝, 2011, (2): 64-67.

    [10]烏云其其格, 廖子龍. 中溫固化樹脂/碳布復(fù)合材料工藝性能研究[J]. 航空材料學(xué)報, 2006, (5): 65-69.

    [11]羅立, 唐慶如. 玻璃纖維對環(huán)氧基體預(yù)浸料固化特征的影響研究[J]. 玻璃鋼/復(fù)合材料, 2014, (4): 13-17.

    [12]邱軍, 魏良明. 芳綸纖維/環(huán)氧樹脂復(fù)合材料固化工藝的研究[J]. 沈陽化工學(xué)院學(xué)報, 2000, (3): 201-204.

    [13]高相南, 孫志杰, 顧軼卓, 等. 碳纖維對環(huán)氧樹脂固化反應(yīng)過程的影響[J]. 玻璃鋼/復(fù)合材料, 2012, (5): 47-53.

    [14]M. L. Costa, E. C. Botelho, M. C. Rezende. Monitoring of cure kinetic prepreg and cure cycle modeling[J]. Journal of Materials Science, 2006.

    [15]Liangfeng Sun, SuSeng Pang, Arthur M. Sterling, Ioan I. Negulescu,Michael A. Stubblefield. Thermal analysis of curing process of epoxy prepreg[J]. J. Appl. Polym. Sci., 2001.

    [16]賈玉瑛, 王榮國, 梁成利. 固化度對碳纖維復(fù)合材料熱機(jī)械性能及耐熱性的影響[J]. 纖維復(fù)合材料, 1986, (5): 15-21.

    [17]李艷國, 丁和斌, 李志東. 固化度對界面粘結(jié)強(qiáng)度及耐熱沖擊性能的影響研究[J]. 印制電路信息, 2009, (S1): 99-104.

    [18]張紀(jì)奎, 酈正能, 關(guān)志東, 等. 固化度與固化收縮對非對稱復(fù)合材料層合板固化變形的影響[J]. 復(fù)合材料學(xué)報, 2007, (2): 120-124.

    [19]龔文化, 許亞洪. 中溫環(huán)氧樹脂體系9368固化工藝與固化度關(guān)系研究[C]//中國力學(xué)學(xué)會. 第十五屆全國復(fù)合材料學(xué)術(shù)會議論文集(上冊). 2008: 4.

    [20]王曉潔, 謝群煒, 張煒, 等. 環(huán)氧樹脂基體固化研究[J].玻璃鋼/復(fù)合材料, 2001, (2): 10-12.

    [21]D. Kim, T. Centea, S.R. Nutt. Out-Time Effects on Cure Kinetics and Viscosity for an Out-Of-Autoclave (OOA) Prepreg: Modelling and Monitoring[J]. Composites Science and Technology, 2014.

    [22]付偉. 預(yù)浸料玻璃化轉(zhuǎn)變溫度測試方法研究[J]. 機(jī)械管理開發(fā), 2010, (6): 32-33.

    [23]Laraib A. Khan, Ayesha Kausar, Syed Tajammul Hussain, et al. Cure characterization of Cycom 977-2A carbon/epoxy composites for quickstep processing[J]. Polym Eng Sci, 2014: 544.

    [24]S. M. Sabzevari, S. Alavi-Soltani, B. Minaie. New Method for Estimating the Extent of Curing of Thermosetting Prepregs[J]. Wiley Online Library, 2011: 4-25.

    [25]Vyazovkin S. Thermal analysis[J]. Anal Chem, 2008, 80(12): 430.

    [26]Kissinger D E. Reaction kinetic in differential thermal analysis[J]. Analytical Chemistry, 1957, 29(11): 1 702-1 706.

    [27]Loris Giorgini, Laura Mazzocchetti, Tiziana Benelli, et al. Kinetics and modeling of curing behavior for two different prepregs based on the same epoxy precursor: A case study for the industrial design of thick composites[J]. Polym. Compos., 2013.

    [28]Mehdy Vafayan, Mohammad Hosain Beheshty, Mir Hamid Reza Ghoreishy, Hossein Abedini. Advanced integral isoconversional analysis for evaluating and predicting the kinetic parameters of the curing reaction of epoxy prepreg[J]. Thermochimica Acta, 2013.

    [29]楊光, 龍國榮. NY9200G樹脂體系非等溫DSC法固化動力學(xué)研究[J]. 洪都科技, 2006, (3): 37-41.

    [30]祁國成, 張博明, 葉金蕊, 等. T800級碳纖維/環(huán)氧預(yù)浸料固化過程研究[C]//中國科學(xué)技術(shù)協(xié)會、河北省人民政府. 第十四屆中國科協(xié)年會第11分會場: 低成本、高性能復(fù)合材料發(fā)展論壇論文集. 2012: 6.

    [31]王洪恩, 楊洋, 陳璐圓, 等. 一種高溫碳纖/環(huán)氧預(yù)浸料的固化反應(yīng)動力學(xué)[J]. 航空制造技術(shù), 2013, (Z2): 111-113.

    [32]柯紅軍, 喬英杰, 陳俊英, 等. 等溫DSC法研究液態(tài)橡膠改性環(huán)氧樹脂固化動力學(xué)[J]. 玻璃鋼/復(fù)合材料, 2012, (3): 43-47.

    [33]Mehran Hayaty, Mohammad Hosain Beheshty, Masoud Esfandeh. Cure kinetics of a glass/epoxy prepreg by dynamic differential scanning calorimetry[J]. J. Appl. Polym. Sci., 2010.

    [34]楊立平, 李鵬, 薛忠民, 等. 聚砜納米纖維增韌環(huán)氧樹脂體系的固化動力學(xué)研究[J]. 玻璃鋼/復(fù)合材料, 2007, (5): 7-10.

    [35]Shin D D, Hahn H T. A consistent kinetic model for AS4/3502 graphite/epoxy[J]. Composites: Part A: applied science and manufacturing, 2000, 31(4): 991-999.

    [36]Fraga F, Soto V H, Rodriguez-Nunez E, et al. Cure kinetic of the epoxy network diglycidyl ether of bisphenol A (BADGE n=0)/Amantidine[J]. J Therm Anal Calorim, 2007, 87(1): 97.

    [37]Khanna U, Chanda M. A kinetic scheme for anhydride curing of diglycidyl ester with tertiary amine as catalyst[J]. J Appl Polym Sci, 1993, 50(9): 1 635.

    [38]劉明昌. 碳纖維增強(qiáng)環(huán)氧樹脂預(yù)浸料及其復(fù)合材料的制備與性能研究[D]. 華東理工大學(xué), 2012.

    [39]王洪恩, 楊洋, 陳璐圓, 等. 一種高溫碳纖/環(huán)氧預(yù)浸料的固化反應(yīng)動力學(xué)[J]. 航空制造技術(shù), 2013, (Z2):111-113.

    [40]Liangfeng Sun, SuSeng Pang, Arthur M. Sterling, Ioan I. Negulescu, Michael A. Stubblefield. Dynamic modeling of curing process of epoxy prepreg[J]. J. Appl. Polym. Sci., 2002,. [41]周正剛, 張寶艷, 劉天舒, 等. LT-01低溫固化復(fù)合材料的反應(yīng)動力學(xué)與工藝優(yōu)化研究[C]//中國宇航學(xué)會.復(fù)合材料—基礎(chǔ)、創(chuàng)新、高效: 第十四屆全國復(fù)合材料學(xué)術(shù)會議論文集(下). 2006: 5.

    Evaluation and optimization of the prepreg's curing process via DSC analysis

    TAO Jia-xiang, LI Wei
    ( College of Textiles, Donghua University, Key Laboratory of Textile Science &Technology, Ministry of Education, Shanghai 201620 China )

    The curing stage of prepreg, as the intermediate matrix composite, directly influence the properties of the finished composite materials when the properties of the composite are determined by the curing process of composite. Differential Scanning Calorimetry(DSC), as a kind of thermal analysis instrument, can measure the heat process of the curing well. This paper mainly introduces the influencing factors when DSC is in use, the determination of curing parameter, the evaluation of the curing process and how to determine, evaluate and optimize the curing process by the curing kinetics theory of prepreg.

    prepreg; differential scanning calorimetry; curing process; evaluation; optimization

    TB332

    A

    1007-9815(2016)01-0066-06

    定稿日期:2016-02-20

    上海市經(jīng)信委產(chǎn)學(xué)研合作計劃(滬CXY-2013-11)

    陶家祥(1991-),男,江蘇淮安人,在讀研究生,主要研究方向為碳纖維預(yù)浸料固化工藝,(電子信箱)jimjiaxiang @163.com。

    猜你喜歡
    環(huán)氧樹脂樹脂動力學(xué)
    《空氣動力學(xué)學(xué)報》征稿簡則
    2018年1—4月我國PVC樹脂產(chǎn)量
    聚氯乙烯(2018年5期)2018-02-18 03:30:28
    2018年3月PVC樹脂進(jìn)出口數(shù)據(jù)
    聚氯乙烯(2018年5期)2018-02-18 03:30:28
    粉末涂料用環(huán)氧樹脂的合成研究
    上海建材(2017年2期)2017-07-21 14:02:10
    磁性離子交換樹脂的制備及其對Cr3+的吸附
    碳納米管陣列/環(huán)氧樹脂的導(dǎo)熱導(dǎo)電性能
    可膨脹石墨對環(huán)氧樹脂的阻燃改性
    高導(dǎo)熱填充型環(huán)氧樹脂復(fù)合材料研究進(jìn)展
    中國塑料(2015年7期)2015-10-14 01:02:35
    基于隨機(jī)-動力學(xué)模型的非均勻推移質(zhì)擴(kuò)散
    TNAE的合成和熱分解動力學(xué)
    久久天堂一区二区三区四区| 90打野战视频偷拍视频| 久久精品91蜜桃| 亚洲最大成人中文| 国产激情欧美一区二区| 在线观看日韩欧美| 亚洲国产欧美一区二区综合| 国产三级黄色录像| 18禁国产床啪视频网站| 欧美日本亚洲视频在线播放| 国产精品免费视频内射| 搡老熟女国产l中国老女人| 亚洲在线自拍视频| 欧美久久黑人一区二区| 国产亚洲欧美98| 在线观看免费日韩欧美大片| 国产精品久久久久久人妻精品电影| 国产亚洲欧美在线一区二区| av在线播放免费不卡| 国产三级在线视频| 老汉色∧v一级毛片| av免费在线观看网站| 欧美日本视频| 高清在线国产一区| 国内久久婷婷六月综合欲色啪| 国产久久久一区二区三区| 欧美国产日韩亚洲一区| 国产激情久久老熟女| 桃色一区二区三区在线观看| 国产成人av教育| 真人一进一出gif抽搐免费| 久久九九热精品免费| 国产精品 国内视频| 国产亚洲精品久久久久5区| 一个人免费在线观看的高清视频| 亚洲精品色激情综合| 免费在线观看完整版高清| 久久狼人影院| 国产精品久久久久久亚洲av鲁大| 波多野结衣高清无吗| 国产精品久久久久久亚洲av鲁大| 色哟哟哟哟哟哟| 精品一区二区三区视频在线观看免费| 亚洲中文字幕日韩| 亚洲三区欧美一区| ponron亚洲| 日日爽夜夜爽网站| 欧洲精品卡2卡3卡4卡5卡区| 老司机午夜福利在线观看视频| 99在线视频只有这里精品首页| 青草久久国产| 99在线视频只有这里精品首页| 精品欧美国产一区二区三| 日日干狠狠操夜夜爽| 亚洲七黄色美女视频| 亚洲,欧美精品.| 国产私拍福利视频在线观看| 久久中文看片网| 美女高潮喷水抽搐中文字幕| 国产91精品成人一区二区三区| 日韩 欧美 亚洲 中文字幕| 男人舔奶头视频| 精品国产亚洲在线| 精品久久久久久久久久久久久 | 成人三级黄色视频| 久久香蕉精品热| 黑人欧美特级aaaaaa片| 国产片内射在线| 亚洲欧洲精品一区二区精品久久久| 久久这里只有精品19| 国产国语露脸激情在线看| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色片一级片一级黄色片| 最新在线观看一区二区三区| 十分钟在线观看高清视频www| 一级a爱片免费观看的视频| 欧美成人午夜精品| 亚洲中文字幕日韩| 亚洲 欧美 日韩 在线 免费| 精品少妇一区二区三区视频日本电影| 免费无遮挡裸体视频| 久久精品亚洲精品国产色婷小说| 99久久99久久久精品蜜桃| 美国免费a级毛片| 国产精品乱码一区二三区的特点| 黑人操中国人逼视频| 国产精品亚洲美女久久久| 757午夜福利合集在线观看| 国产精品一区二区精品视频观看| 亚洲 欧美 日韩 在线 免费| 悠悠久久av| 国产成人欧美| 伊人久久大香线蕉亚洲五| 嫩草影视91久久| 村上凉子中文字幕在线| 久久精品人妻少妇| 自线自在国产av| 999精品在线视频| 国产三级在线视频| 欧美黑人精品巨大| 深夜精品福利| 黄色视频不卡| 国产一区二区三区在线臀色熟女| 日本精品一区二区三区蜜桃| 日韩欧美三级三区| 在线天堂中文资源库| 欧美日韩乱码在线| 美女高潮到喷水免费观看| 亚洲av电影在线进入| 叶爱在线成人免费视频播放| 大型av网站在线播放| 久久久久国内视频| 亚洲人成伊人成综合网2020| 夜夜夜夜夜久久久久| 特大巨黑吊av在线直播 | 男女午夜视频在线观看| 久久99热这里只有精品18| 两个人免费观看高清视频| 亚洲男人天堂网一区| 19禁男女啪啪无遮挡网站| 亚洲av电影在线进入| 国产高清有码在线观看视频 | av免费在线观看网站| 欧美日韩黄片免| 亚洲国产精品合色在线| 免费无遮挡裸体视频| 亚洲三区欧美一区| 老汉色∧v一级毛片| 亚洲av成人不卡在线观看播放网| 一级a爱片免费观看的视频| 男人舔女人下体高潮全视频| 免费在线观看完整版高清| 亚洲中文av在线| 免费女性裸体啪啪无遮挡网站| 又大又爽又粗| 可以在线观看的亚洲视频| 在线av久久热| 级片在线观看| 亚洲 欧美一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 成年版毛片免费区| 老汉色∧v一级毛片| 男女做爰动态图高潮gif福利片| 免费av毛片视频| 亚洲男人天堂网一区| 国产欧美日韩精品亚洲av| 欧美性猛交╳xxx乱大交人| 又大又爽又粗| 免费在线观看亚洲国产| 国产久久久一区二区三区| 在线国产一区二区在线| 国语自产精品视频在线第100页| 免费在线观看影片大全网站| 亚洲成人精品中文字幕电影| av欧美777| 一二三四社区在线视频社区8| 亚洲精品美女久久av网站| 久久精品夜夜夜夜夜久久蜜豆 | 国产精品乱码一区二三区的特点| 女同久久另类99精品国产91| 神马国产精品三级电影在线观看 | 国产爱豆传媒在线观看 | 国产精品香港三级国产av潘金莲| 十八禁网站免费在线| 日韩欧美免费精品| 午夜免费激情av| 99久久综合精品五月天人人| 国产精品一区二区免费欧美| 搞女人的毛片| 精品国产超薄肉色丝袜足j| 欧美在线黄色| www.999成人在线观看| 搡老妇女老女人老熟妇| 嫁个100分男人电影在线观看| bbb黄色大片| 久久九九热精品免费| 欧美成人性av电影在线观看| 麻豆av在线久日| 国产国语露脸激情在线看| 在线播放国产精品三级| 窝窝影院91人妻| 婷婷精品国产亚洲av| 亚洲av成人不卡在线观看播放网| 免费av毛片视频| 久久中文字幕人妻熟女| 一边摸一边抽搐一进一小说| 免费在线观看成人毛片| 久久久久精品国产欧美久久久| 女人爽到高潮嗷嗷叫在线视频| 国产一级毛片七仙女欲春2 | 免费在线观看成人毛片| 国产片内射在线| 91在线观看av| 亚洲成人久久性| 两个人免费观看高清视频| 午夜福利免费观看在线| 夜夜爽天天搞| e午夜精品久久久久久久| 久久午夜综合久久蜜桃| 人人妻人人看人人澡| 成人特级黄色片久久久久久久| 无遮挡黄片免费观看| 久久中文看片网| 91av网站免费观看| 美女 人体艺术 gogo| 国产av一区在线观看免费| 精品高清国产在线一区| 熟女电影av网| 美女大奶头视频| 欧美黄色片欧美黄色片| 51午夜福利影视在线观看| 别揉我奶头~嗯~啊~动态视频| 高清在线国产一区| 亚洲av中文字字幕乱码综合 | 黑丝袜美女国产一区| 亚洲国产欧美一区二区综合| 十分钟在线观看高清视频www| 观看免费一级毛片| 欧美日韩瑟瑟在线播放| 日韩欧美国产一区二区入口| 老汉色av国产亚洲站长工具| 黄片大片在线免费观看| 老司机福利观看| 熟妇人妻久久中文字幕3abv| 国产99白浆流出| 欧美黑人精品巨大| 啪啪无遮挡十八禁网站| 十分钟在线观看高清视频www| 一级毛片女人18水好多| 黑人巨大精品欧美一区二区mp4| 日本五十路高清| 国产三级黄色录像| 国产成人影院久久av| 欧美成人免费av一区二区三区| 中文字幕高清在线视频| 亚洲三区欧美一区| 国产真人三级小视频在线观看| 国产一区二区三区在线臀色熟女| 桃红色精品国产亚洲av| 婷婷丁香在线五月| netflix在线观看网站| 身体一侧抽搐| 亚洲男人天堂网一区| 久久精品成人免费网站| av片东京热男人的天堂| 亚洲av中文字字幕乱码综合 | 国产精品 国内视频| 久久久国产成人精品二区| 亚洲av片天天在线观看| 丝袜在线中文字幕| 久久中文看片网| 这个男人来自地球电影免费观看| 欧美丝袜亚洲另类 | 亚洲avbb在线观看| 国产激情欧美一区二区| 欧美日韩精品网址| 亚洲av成人不卡在线观看播放网| 久久久久久久久久黄片| 国产精品久久视频播放| 美女扒开内裤让男人捅视频| 欧美乱码精品一区二区三区| 日韩欧美国产一区二区入口| 久久久久久大精品| 99久久国产精品久久久| 国产激情偷乱视频一区二区| 不卡一级毛片| 1024香蕉在线观看| 国产精品久久视频播放| 国产成人啪精品午夜网站| 两性午夜刺激爽爽歪歪视频在线观看 | 69av精品久久久久久| 99久久精品国产亚洲精品| 9191精品国产免费久久| 一区二区三区高清视频在线| www.999成人在线观看| 国产欧美日韩一区二区精品| 黄色a级毛片大全视频| 动漫黄色视频在线观看| 午夜免费观看网址| 一边摸一边抽搐一进一小说| 在线播放国产精品三级| 777久久人妻少妇嫩草av网站| 男人舔女人下体高潮全视频| 两个人视频免费观看高清| 欧美中文日本在线观看视频| 亚洲午夜理论影院| 亚洲美女黄片视频| 免费女性裸体啪啪无遮挡网站| 老司机午夜福利在线观看视频| 中文字幕精品免费在线观看视频| 一个人免费在线观看的高清视频| 免费看十八禁软件| 日韩欧美国产在线观看| 亚洲专区国产一区二区| 日本 av在线| 一级片免费观看大全| 老司机福利观看| 女同久久另类99精品国产91| 丰满的人妻完整版| 俺也久久电影网| 国产精品一区二区免费欧美| 免费在线观看影片大全网站| 变态另类丝袜制服| 无人区码免费观看不卡| 窝窝影院91人妻| 最好的美女福利视频网| 久久人妻福利社区极品人妻图片| 欧美黄色淫秽网站| 好看av亚洲va欧美ⅴa在| 欧美+亚洲+日韩+国产| 露出奶头的视频| 两人在一起打扑克的视频| 国产精品乱码一区二三区的特点| 欧美在线黄色| 视频在线观看一区二区三区| 国产在线精品亚洲第一网站| 亚洲精品中文字幕在线视频| 色综合站精品国产| 中文字幕精品亚洲无线码一区 | 一进一出抽搐gif免费好疼| 天堂√8在线中文| 精品无人区乱码1区二区| 好看av亚洲va欧美ⅴa在| 伊人久久大香线蕉亚洲五| 欧美日韩亚洲综合一区二区三区_| 可以免费在线观看a视频的电影网站| 无限看片的www在线观看| 一a级毛片在线观看| 一级片免费观看大全| 九色国产91popny在线| 18禁观看日本| 免费看美女性在线毛片视频| 首页视频小说图片口味搜索| 日韩国内少妇激情av| 久久亚洲精品不卡| 午夜福利视频1000在线观看| 我的亚洲天堂| www.自偷自拍.com| 国产亚洲精品久久久久5区| 日韩欧美在线二视频| 侵犯人妻中文字幕一二三四区| 天堂动漫精品| 免费在线观看成人毛片| 欧美日韩瑟瑟在线播放| 久久婷婷人人爽人人干人人爱| 禁无遮挡网站| 亚洲av中文字字幕乱码综合 | 99精品久久久久人妻精品| 色播在线永久视频| 无遮挡黄片免费观看| 日本 欧美在线| 国产激情偷乱视频一区二区| 亚洲成人精品中文字幕电影| 免费人成视频x8x8入口观看| 亚洲国产中文字幕在线视频| 亚洲熟女毛片儿| 色综合站精品国产| 免费看日本二区| 又黄又爽又免费观看的视频| 18禁国产床啪视频网站| bbb黄色大片| 两性午夜刺激爽爽歪歪视频在线观看 | 啦啦啦观看免费观看视频高清| 久久久久久人人人人人| 亚洲精品国产精品久久久不卡| 久久这里只有精品19| 欧美激情久久久久久爽电影| 免费电影在线观看免费观看| 日韩有码中文字幕| 在线十欧美十亚洲十日本专区| 欧美午夜高清在线| 在线观看免费午夜福利视频| 热re99久久国产66热| 丁香欧美五月| 国产色视频综合| 在线国产一区二区在线| 久久国产精品男人的天堂亚洲| 日韩高清综合在线| 国产激情欧美一区二区| 美女免费视频网站| 国产色视频综合| 人妻丰满熟妇av一区二区三区| 悠悠久久av| 亚洲av熟女| 久久久久久久久免费视频了| av中文乱码字幕在线| 搡老妇女老女人老熟妇| 欧美国产日韩亚洲一区| 国产黄片美女视频| 女性生殖器流出的白浆| 特大巨黑吊av在线直播 | 欧美成人性av电影在线观看| 国产在线精品亚洲第一网站| 男人舔奶头视频| 日本a在线网址| 99精品在免费线老司机午夜| 搡老岳熟女国产| 大香蕉久久成人网| 色综合婷婷激情| 日韩 欧美 亚洲 中文字幕| 十八禁人妻一区二区| 亚洲中文日韩欧美视频| 一a级毛片在线观看| 青草久久国产| 亚洲国产日韩欧美精品在线观看 | 亚洲第一av免费看| 性色av乱码一区二区三区2| 日韩欧美免费精品| 国产高清激情床上av| 免费一级毛片在线播放高清视频| 91麻豆av在线| 桃红色精品国产亚洲av| 亚洲av美国av| 国内毛片毛片毛片毛片毛片| 青草久久国产| 国产熟女xx| 国产欧美日韩一区二区三| 热99re8久久精品国产| av电影中文网址| 亚洲成av片中文字幕在线观看| 亚洲第一欧美日韩一区二区三区| 中文字幕最新亚洲高清| 伊人久久大香线蕉亚洲五| 99久久无色码亚洲精品果冻| 国产99久久九九免费精品| 午夜福利一区二区在线看| 三级毛片av免费| 在线观看66精品国产| 成人免费观看视频高清| 亚洲欧洲精品一区二区精品久久久| 精品人妻1区二区| 一区二区三区精品91| 女人爽到高潮嗷嗷叫在线视频| 国产成年人精品一区二区| 欧美一级毛片孕妇| 午夜视频精品福利| 免费搜索国产男女视频| 在线国产一区二区在线| 久久久久久久久久黄片| 日韩欧美在线二视频| 最好的美女福利视频网| 国产精华一区二区三区| 熟妇人妻久久中文字幕3abv| 十八禁网站免费在线| 久久精品国产99精品国产亚洲性色| 亚洲国产高清在线一区二区三 | 嫁个100分男人电影在线观看| 国产免费男女视频| 黄色视频不卡| 色综合欧美亚洲国产小说| 18禁观看日本| 18禁裸乳无遮挡免费网站照片 | 老熟妇乱子伦视频在线观看| 女性生殖器流出的白浆| 国产亚洲欧美精品永久| 国产高清视频在线播放一区| 99riav亚洲国产免费| 欧美乱色亚洲激情| 国产av不卡久久| www.精华液| 伊人久久大香线蕉亚洲五| 青草久久国产| 国语自产精品视频在线第100页| 婷婷六月久久综合丁香| 制服丝袜大香蕉在线| 日本一区二区免费在线视频| 精品国产亚洲在线| 亚洲黑人精品在线| 久久精品亚洲精品国产色婷小说| 香蕉av资源在线| 国产激情久久老熟女| 亚洲自拍偷在线| 国产精品美女特级片免费视频播放器 | 男女午夜视频在线观看| 国产av又大| 黑人巨大精品欧美一区二区mp4| 亚洲黑人精品在线| 女人爽到高潮嗷嗷叫在线视频| 日本精品一区二区三区蜜桃| 国产av一区在线观看免费| 男人舔女人下体高潮全视频| 国内精品久久久久久久电影| 男人舔女人的私密视频| 最近在线观看免费完整版| 一二三四在线观看免费中文在| 两个人视频免费观看高清| 久久精品国产综合久久久| 日韩三级视频一区二区三区| 人人澡人人妻人| 欧美午夜高清在线| 色尼玛亚洲综合影院| 国产三级在线视频| 黄色片一级片一级黄色片| 岛国在线观看网站| 免费观看人在逋| 亚洲专区字幕在线| 女生性感内裤真人,穿戴方法视频| 国产高清videossex| 在线观看66精品国产| 两个人视频免费观看高清| 成人永久免费在线观看视频| 天天添夜夜摸| 亚洲中文日韩欧美视频| АⅤ资源中文在线天堂| 国产精品一区二区免费欧美| 国产私拍福利视频在线观看| 午夜久久久在线观看| 午夜影院日韩av| 国产亚洲欧美98| 精品人妻1区二区| 99热6这里只有精品| 好男人电影高清在线观看| 黄色片一级片一级黄色片| 亚洲国产欧美日韩在线播放| 国产色视频综合| 亚洲三区欧美一区| 麻豆av在线久日| 亚洲av五月六月丁香网| 免费看美女性在线毛片视频| 韩国av一区二区三区四区| 露出奶头的视频| 亚洲av成人一区二区三| 久久久水蜜桃国产精品网| 午夜久久久在线观看| 一a级毛片在线观看| 精品日产1卡2卡| 日韩精品免费视频一区二区三区| 老司机午夜福利在线观看视频| 亚洲美女黄片视频| 国产97色在线日韩免费| 色在线成人网| 久久精品国产综合久久久| 久久亚洲精品不卡| 亚洲精品国产一区二区精华液| 国产精品久久久av美女十八| 精品久久蜜臀av无| 亚洲精品一区av在线观看| 美女 人体艺术 gogo| 欧美一区二区精品小视频在线| 亚洲一卡2卡3卡4卡5卡精品中文| 999久久久精品免费观看国产| 99riav亚洲国产免费| 久久精品91蜜桃| 男人操女人黄网站| 悠悠久久av| 国产真实乱freesex| 丝袜在线中文字幕| 午夜精品久久久久久毛片777| 久久精品91蜜桃| 久久九九热精品免费| 国产单亲对白刺激| 国产男靠女视频免费网站| xxxwww97欧美| 久久久国产欧美日韩av| e午夜精品久久久久久久| 2021天堂中文幕一二区在线观 | 一区二区三区高清视频在线| 可以在线观看的亚洲视频| 欧美+亚洲+日韩+国产| 欧美中文日本在线观看视频| 国产在线观看jvid| 亚洲精品国产一区二区精华液| 一区二区日韩欧美中文字幕| 欧美中文综合在线视频| 在线永久观看黄色视频| 国产aⅴ精品一区二区三区波| 久久中文看片网| 日韩大尺度精品在线看网址| 国产欧美日韩精品亚洲av| 脱女人内裤的视频| 90打野战视频偷拍视频| 99久久99久久久精品蜜桃| 日韩视频一区二区在线观看| 国产亚洲精品综合一区在线观看 | 欧美成人性av电影在线观看| 日韩精品中文字幕看吧| a级毛片在线看网站| 亚洲五月色婷婷综合| 好男人在线观看高清免费视频 | 欧美日韩福利视频一区二区| 国产爱豆传媒在线观看 | 色av中文字幕| 在线永久观看黄色视频| 黄片播放在线免费| 国产精品九九99| 黄片播放在线免费| 哪里可以看免费的av片| 一级毛片高清免费大全| 免费观看人在逋| 色老头精品视频在线观看| 免费看十八禁软件| 国产午夜福利久久久久久| 一进一出抽搐动态| 国产精品影院久久| 午夜福利成人在线免费观看| 色综合亚洲欧美另类图片| 男女那种视频在线观看| 色综合站精品国产| 中文字幕精品免费在线观看视频| 人成视频在线观看免费观看| 女同久久另类99精品国产91| 国产欧美日韩一区二区精品| 欧美又色又爽又黄视频| 久久久久久亚洲精品国产蜜桃av| 国产又色又爽无遮挡免费看| 伦理电影免费视频| 欧美黑人欧美精品刺激| 久久 成人 亚洲| av在线天堂中文字幕| 欧美激情久久久久久爽电影| 日韩大码丰满熟妇| 少妇粗大呻吟视频| 女人被狂操c到高潮| 久久99热这里只有精品18| 国内精品久久久久精免费|