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

    考慮寬應(yīng)變率的修正Johnson-Cook模型及蜂窩板動(dòng)態(tài)力學(xué)響應(yīng)模擬

    2021-07-08 06:52:32李雙蓓常巖軍
    關(guān)鍵詞:芯層蒙皮金屬材料

    李雙蓓,趙 璇,吳 園,常巖軍

    ( 1. 廣西大學(xué)土木建筑工程學(xué)院,廣西 南寧,530004;2. 廣西大學(xué)工程防災(zāi)與結(jié)構(gòu)安全教育部重點(diǎn)實(shí)驗(yàn)室,廣西 南寧,530004;3.廣西大學(xué)廣西防災(zāi)減災(zāi)與工程安全重點(diǎn)實(shí)驗(yàn)室,廣西 南寧,530004)

    金屬蜂窩板是由蒙皮、膠層和蜂窩芯層組成的夾層結(jié)構(gòu),憑借其比強(qiáng)度高、抗沖擊性能好等優(yōu)點(diǎn)而廣泛應(yīng)用于建筑幕墻、室內(nèi)裝飾、車輛船舶等領(lǐng)域,近年來(lái),該結(jié)構(gòu)在動(dòng)荷載下的力學(xué)研究已成為相關(guān)領(lǐng)域的熱點(diǎn)[1-2],尤其數(shù)值仿真技術(shù)的不斷進(jìn)步,不僅突破了傳統(tǒng)實(shí)驗(yàn)研究手段的局限[3],還能明顯提高分析效率,而仿真結(jié)果的精確程度則直接取決于材料本構(gòu)模型的正確與否。

    通常,在沖擊荷載作用下,金屬材料的應(yīng)變率變化范圍較寬,低、中、高3種應(yīng)變率階段均會(huì)出現(xiàn)。Johnson等[4]通過(guò)分析實(shí)驗(yàn)數(shù)據(jù),總結(jié)出應(yīng)變率效應(yīng)項(xiàng)、應(yīng)變硬化項(xiàng)及溫度項(xiàng)乘積形式的Johnson-Cook經(jīng)驗(yàn)型本構(gòu)模型(JC模型),該模型因具有物理意義清晰、參數(shù)易標(biāo)定等優(yōu)點(diǎn)而成為分析金屬材料動(dòng)態(tài)響應(yīng)的主流模型,但其缺陷在于不能同時(shí)滿足金屬寬應(yīng)變率下的本構(gòu)關(guān)系,因此需對(duì)模型進(jìn)行修正[5-6]。在修正JC模型應(yīng)變率效應(yīng)項(xiàng)方面,孟利平等[7]使用Cowper-Symonds(CS)應(yīng)變率強(qiáng)化模型[8]替換JC模型中的應(yīng)變率項(xiàng)從而獲得了能夠更好表達(dá)應(yīng)變率強(qiáng)化效應(yīng)的JC-CS模型;Kang等[9]將JC模型中的1次應(yīng)變率效應(yīng)項(xiàng)修正為高階應(yīng)變率效應(yīng)項(xiàng),所得修正模型更適用于描述車身在高速碰撞下的動(dòng)態(tài)響應(yīng);Huh等[10]在JC模型中添加2次應(yīng)變率效應(yīng)項(xiàng)后,模型預(yù)測(cè)值與實(shí)驗(yàn)數(shù)據(jù)雖仍有一定差距,但二者偏差較之前明顯減小。此外,因?yàn)楦邞?yīng)變率下金屬的應(yīng)變硬化特征會(huì)發(fā)生顯著變化,所以金屬在不同應(yīng)變率階段所遵循的應(yīng)變硬化準(zhǔn)則也不盡相同[5]。在對(duì)JC模型應(yīng)變硬化項(xiàng)進(jìn)行修正時(shí),Kim等[11]提出的H/V-R模型線性組合了Hollomon模型[12]和Voce準(zhǔn)則[13],能夠較為準(zhǔn)確地反映鋼材動(dòng)態(tài)行為隨應(yīng)變率提高而變化的現(xiàn)象,但由于該模型具有高度非線性性質(zhì),加之材料參數(shù)設(shè)置也存在一定局限性,因此難以準(zhǔn)確估計(jì)不同應(yīng)變率條件下的材料應(yīng)力-應(yīng)變關(guān)系,不利于實(shí)際工程應(yīng)用;李營(yíng)等[14]修正后的線性疊加模型以應(yīng)變1次項(xiàng)和2次項(xiàng)的和為應(yīng)變硬化項(xiàng),研究表明,該模型以加和形式雖然提高了對(duì)金屬行為預(yù)測(cè)的能力,但擬合所得曲線接近于一次函數(shù),與實(shí)測(cè)值變化的趨勢(shì)差異較大。綜上可知,僅僅通過(guò)修正JC模型中的應(yīng)變率效應(yīng)項(xiàng)或應(yīng)變硬化項(xiàng)均不能有效提高模型在寬應(yīng)變率條件下的適用性。有鑒于此,本文擬對(duì)常溫下的JC本構(gòu)模型應(yīng)變硬化項(xiàng)和應(yīng)變率效應(yīng)項(xiàng)進(jìn)行修正,提出適用于寬應(yīng)變率條件的H/V-CS本構(gòu)模型并對(duì)模型的正確性進(jìn)行驗(yàn)證,再通過(guò)顯式VUMAT接口將H/V-CS本構(gòu)模型導(dǎo)入ABAQUS有限元分析軟件中,建立沖擊荷載作用下的蜂窩板結(jié)構(gòu)模型,進(jìn)而基于該模型分析蒙皮厚度、芯層壁厚和芯層高度等因素對(duì)蜂窩板吸能效果的影響,以期為建立更符合實(shí)際的寬應(yīng)變率金屬動(dòng)態(tài)本構(gòu)關(guān)系及蜂窩板在實(shí)際工程中的防護(hù)和優(yōu)化提供參考。

    1 寬應(yīng)變率下H/V-CS金屬動(dòng)態(tài)本構(gòu)模型

    1.1 H/V-CS本構(gòu)模型的建立

    1.1.1 JC本構(gòu)模型

    JC本構(gòu)模型是被廣泛應(yīng)用的一種唯象型動(dòng)態(tài)本構(gòu)模型,將金屬材料的真實(shí)應(yīng)力直觀地表示為應(yīng)變硬化、應(yīng)變率修正項(xiàng)和溫度項(xiàng)的乘積,當(dāng)不考慮溫度效應(yīng)時(shí),其基本形式為:

    (1)

    1.1.2 應(yīng)變硬化項(xiàng)修正

    大多數(shù)金屬材料在不同應(yīng)變率下所表現(xiàn)出的硬化規(guī)律并不一致,所遵循的準(zhǔn)則往往隨著應(yīng)變率的增大由Ludwik準(zhǔn)則[15]轉(zhuǎn)變?yōu)閂oce準(zhǔn)則,這主要是因?yàn)長(zhǎng)udwik準(zhǔn)則為冪函數(shù)形式,真實(shí)應(yīng)力隨真實(shí)應(yīng)變的增加而不斷增大,而Voce準(zhǔn)則為指數(shù)型,應(yīng)力在大應(yīng)變條件下會(huì)不斷接近飽和狀態(tài),即應(yīng)力-應(yīng)變曲線會(huì)出現(xiàn)水平漸近線[5],這與實(shí)際高應(yīng)變率條件下金屬材料真實(shí)應(yīng)力會(huì)隨應(yīng)變的增大而趨于飽和的現(xiàn)象[16]一致,故Voce準(zhǔn)則更適于展現(xiàn)金屬材料在高應(yīng)變率下的本構(gòu)關(guān)系。

    Kim等[11]和李營(yíng)等[14]針對(duì)硬化規(guī)律對(duì)JC本構(gòu)模型的應(yīng)變硬化項(xiàng)進(jìn)行改進(jìn),相應(yīng)得到H/V-R模型和線性疊加模型分別如式(2)和式(3)所示:

    (2)

    (3)

    式(2)~式(3)中,α1、α2、m1、m2、a1、a2、P及q均為材料常數(shù);σ0為靜態(tài)真實(shí)應(yīng)力;fH表示Hollomon準(zhǔn)則,引入初始屈服強(qiáng)度值后發(fā)展成為L(zhǎng)udwik準(zhǔn)則,fV表示Voce準(zhǔn)則,二者表達(dá)式分別為:

    fH=Bεn

    (4)

    fV=A2-Dexp(-Eε)

    (5)

    式(5)中,A2為材料常數(shù);D為金屬?gòu)?qiáng)度系數(shù);E為硬化指數(shù)(Voce硬化準(zhǔn)則),可通過(guò)擬合金屬真實(shí)應(yīng)力-應(yīng)變關(guān)系而得到。

    為了準(zhǔn)確描述金屬材料在不同應(yīng)變率階段的應(yīng)變硬化性質(zhì),本文參考H/V-R模型(式(2))及線性疊加模型(式(3))對(duì)應(yīng)變硬化項(xiàng)的處理方式,根據(jù)真實(shí)應(yīng)力-應(yīng)變曲線特征,將Ludwik準(zhǔn)則與Voce準(zhǔn)則疊加獲得新的應(yīng)變硬化項(xiàng),表達(dá)式為:

    σ=A+Bεn-Dexp(-Eε)

    (6)

    式中,A為材料常數(shù)。

    1.1.3 應(yīng)變率項(xiàng)修正

    利用CS模型[8]可獲得材料動(dòng)態(tài)真實(shí)應(yīng)力與應(yīng)變率之間的關(guān)系表達(dá)式為:

    (7)

    式中,σd為動(dòng)態(tài)真實(shí)應(yīng)力;C和m為材料常數(shù)。

    以文獻(xiàn)[14]所報(bào)道的船用低碳鋼、文獻(xiàn)[17]所報(bào)道的921A鋼以及文獻(xiàn)[18]所報(bào)道的EH36鋼為研究對(duì)象,分別利用JC模型應(yīng)變率效應(yīng)項(xiàng)和CS應(yīng)變率強(qiáng)化模型獲得上述鋼種的屈服強(qiáng)度-應(yīng)變率關(guān)系擬合曲線,并將擬合值與相應(yīng)文獻(xiàn)所提供的實(shí)驗(yàn)值進(jìn)行比較,結(jié)果如圖1所示。從圖1

    中可以看出,利用JC本構(gòu)模型應(yīng)變率效應(yīng)項(xiàng)對(duì)船用低碳鋼及921A鋼力學(xué)行為進(jìn)行擬合時(shí),在中、高應(yīng)變率階段所得數(shù)據(jù)與相應(yīng)實(shí)驗(yàn)值差異較大,而利用CS模型所得3種鋼的擬合結(jié)果均與相應(yīng)實(shí)驗(yàn)值吻合度較高,兼顧了材料在低、中、高應(yīng)變率階段的力學(xué)行為變化,滿足寬應(yīng)變率的要求。因此,本文選擇CS模型作為金屬材料本構(gòu)模型中的應(yīng)變率效應(yīng)項(xiàng)。

    1.1.4 H/V-CS本構(gòu)模型

    將修正后的應(yīng)變硬化項(xiàng)與CS模型應(yīng)變率效應(yīng)項(xiàng)耦合,獲得適用于寬應(yīng)變率條件的H/V-CS本構(gòu)模型,表達(dá)式為:

    (8)

    由式(8)可見(jiàn),該本構(gòu)模型表達(dá)式中的參數(shù)較H/V-R模型(式(2))更少,模型應(yīng)變硬化項(xiàng)考慮了金屬在不同應(yīng)變率范圍內(nèi)的真實(shí)應(yīng)力-應(yīng)變關(guān)系,包容性較強(qiáng),擬合方式也較為靈活。同時(shí),模型選用更能準(zhǔn)確反映應(yīng)變率敏感性的CS模型作為應(yīng)變率強(qiáng)化項(xiàng),考慮了金屬材料的應(yīng)變率效應(yīng),使得整個(gè)模型能夠充分體現(xiàn)金屬材料的力學(xué)性能,物理意義清晰,具有通用性。

    1.2 H/V-CS本構(gòu)模型參數(shù)的確定

    本文基于文獻(xiàn)[5]所報(bào)道的Q420鋼拉伸實(shí)驗(yàn)和文獻(xiàn)[19]所報(bào)道的2024-T3鋁的霍普金森拉桿實(shí)驗(yàn),對(duì)JC模型及本文所提H/V-CS模型在寬應(yīng)變率條件下的適用性進(jìn)行討論,這2種金屬材料分別廣泛應(yīng)用于建筑及航空領(lǐng)域,應(yīng)變率變化涵蓋低、中、高3個(gè)階段,具有一定的普適性。

    首先選用數(shù)值優(yōu)化仿真計(jì)算軟件平臺(tái)First Optimization(1stOpt)中的全局優(yōu)化算法,再編寫(xiě)程序輸入待擬合參數(shù)的模型表達(dá)式以及不同應(yīng)變率下真實(shí)應(yīng)力-應(yīng)變實(shí)驗(yàn)數(shù)據(jù),無(wú)需依賴初值,即可得到2種材料的H/V-CS模型參數(shù)如表1所示。

    表1 H/V-CS本構(gòu)模型參數(shù)

    待獲得H/V-CS模型參數(shù)后,再分別使用JC模型和H/V-CS模型對(duì)Q420鋼及2024-T3鋁在多種應(yīng)變率條件下的應(yīng)力-應(yīng)變關(guān)系進(jìn)行預(yù)測(cè),將所得理論模型預(yù)測(cè)結(jié)果與相應(yīng)文獻(xiàn)所提供的實(shí)驗(yàn)值進(jìn)行對(duì)比,結(jié)果如圖2所示。由圖2可知,2種金屬材料在不同應(yīng)變率條件下都表現(xiàn)出明顯的應(yīng)變率強(qiáng)化效應(yīng)。當(dāng)使用JC模型對(duì)金屬材料的應(yīng)力-應(yīng)變關(guān)系進(jìn)行預(yù)測(cè)時(shí),高應(yīng)變率條件下的模型預(yù)測(cè)值與實(shí)驗(yàn)實(shí)測(cè)值差異較大,如圖2(a)中應(yīng)變率分別為133、288 s-1條件下的2條模型預(yù)測(cè)曲線大部分重疊,未表現(xiàn)出應(yīng)變率敏感性,且預(yù)測(cè)值與相應(yīng)的實(shí)驗(yàn)值存在明顯偏差,這種現(xiàn)象在圖2(b)中(應(yīng)變率為0.01、500、1300 s-1)表現(xiàn)得更為明顯,應(yīng)歸因于JC模型所采用的應(yīng)變硬化準(zhǔn)則和應(yīng)變率強(qiáng)化項(xiàng)更適用于中、低應(yīng)變率條件,對(duì)金屬材料高應(yīng)變率階段的力學(xué)行為變化不夠敏感,難以真實(shí)反映金屬材料的應(yīng)力-應(yīng)變關(guān)系,所以在寬應(yīng)變率條件下使用JC模型不能準(zhǔn)確預(yù)測(cè)金屬材料的應(yīng)力-應(yīng)變關(guān)系;當(dāng)使用H/V-CS模型對(duì)金屬材料應(yīng)力-應(yīng)變關(guān)系進(jìn)行預(yù)測(cè)時(shí)(圖2(c)~圖2(d)),2種金屬材料在低、中、高應(yīng)變率條件下的模型預(yù)測(cè)結(jié)果均與相應(yīng)實(shí)驗(yàn)值高度吻合,體現(xiàn)了該模型在寬應(yīng)變率條件下的適用性與準(zhǔn)確性,由此可證明在寬應(yīng)變率條件下,將Voce模型和Ludwik模型組合的應(yīng)變強(qiáng)化模型耦合CS應(yīng)變率強(qiáng)化模型后,所得H/V-CS模型可以更好地表達(dá)金屬材料的應(yīng)變率敏感性,對(duì)金屬材料應(yīng)力-應(yīng)變關(guān)系的預(yù)測(cè)結(jié)果符合相應(yīng)的真實(shí)發(fā)展趨勢(shì)。

    (a) 船用低碳鋼 (b) 921A鋼 (c) EH36鋼

    (a)Q420鋼,JC本構(gòu)模型 (b)2024-T3鋁,JC本構(gòu)模型

    (c)Q420鋼,H/V-CS模型 (d)2024-T3鋁,H/V-CS模型

    1.3 H/V-CS本構(gòu)模型VUMAT開(kāi)發(fā)

    (9)

    (10)

    圖3 VUMAT子程序流程圖

    2 基于H/V-CS本構(gòu)模型的蜂窩板動(dòng)態(tài)力學(xué)響應(yīng)模擬

    2.1 蜂窩板模型的建立與驗(yàn)證

    為驗(yàn)證VUMAT子程序的正確性,本文以蜂窩夾層板為研究對(duì)象,參照文獻(xiàn)[2]建立模型,該結(jié)構(gòu)如圖4所示,其中蒙皮厚度Tf為1 mm,蜂窩單元直徑R為6.35 mm,蜂窩芯層高度h為18mm,蜂窩芯層厚度tc為0.07 mm,蒙皮(2024-T3鋁材質(zhì))的H/V-CS模型參數(shù)見(jiàn)表1,其余部件屬性均參照文獻(xiàn)[2]設(shè)定。采用ABAQUS有限元軟件構(gòu)建蜂窩板有限元模型(見(jiàn)圖5)并模擬蜂窩板沖擊試驗(yàn)(沖擊能為81.07 J),獲得沖擊中心點(diǎn)處沿X、Y、Z方向的應(yīng)變率時(shí)程曲線如圖6所示。由圖6可見(jiàn),沖擊點(diǎn)處沿不同方向的應(yīng)變率均變化復(fù)雜且無(wú)規(guī)律,應(yīng)變率值覆蓋了低、中、高3個(gè)階段,故采用傳統(tǒng)的動(dòng)態(tài)本構(gòu)模型難以滿足精確計(jì)算的要求,在此基于本文所提H/V-CS寬應(yīng)變率動(dòng)態(tài)本構(gòu)模型,借助ABAQUS軟件對(duì)蜂窩板在不同沖擊能下的沖擊過(guò)程進(jìn)行模擬,將蜂窩板損傷深度模擬計(jì)算值與文獻(xiàn)[2]所提供的相應(yīng)的實(shí)驗(yàn)值進(jìn)行對(duì)比,結(jié)果如圖7所示。由圖7可見(jiàn),在不同沖擊能量條件下,由ABAQUS模擬計(jì)算所得蜂窩板損傷深度均與相應(yīng)實(shí)驗(yàn)結(jié)果基本一致,最大誤差不超過(guò)5%,從而驗(yàn)證了H/V-CS本構(gòu)模型與有限元仿真模型的準(zhǔn)確性,因此可利用這兩個(gè)模型進(jìn)一步探究蜂窩板結(jié)構(gòu)參數(shù)對(duì)蜂窩板在沖擊荷載下吸能效果的影響。

    圖4 模型結(jié)構(gòu)

    圖5 蜂窩板有限元模型

    圖6 應(yīng)變率時(shí)程曲線

    圖7 蜂窩板損傷深度

    2.2 蜂窩板結(jié)構(gòu)參數(shù)對(duì)其吸能效果的影響

    設(shè)定圓柱體平頭彈底面半徑為10 mm,高為56.7 mm,保持沖擊能量為81.07 J,利用本文所構(gòu)建的H/V-CS本構(gòu)模型與有限元仿真模型研究了蜂窩板遭受平頭彈沖擊時(shí),蒙皮厚度(0.5、1.0、1.5 mm)、蜂窩芯壁厚(0.035、0.07、0.105 mm)以及蜂窩芯高度(9、18、27 mm)等結(jié)構(gòu)參數(shù)對(duì)其吸能效果的影響,分析結(jié)果見(jiàn)圖8,相應(yīng)條件下蜂窩板的吸收能量穩(wěn)定值如表2所示。由圖8(a)可見(jiàn),當(dāng)蒙皮厚度分別為0.5、1.0、1.5 mm時(shí),相應(yīng)的蜂窩板吸收全部能量的時(shí)長(zhǎng)依次為4.2、3.2、2.8 ms,吸能速率隨蒙皮厚度的增加而增大,因此,增加蒙皮厚度可有效提高蜂窩板吸能效果;由圖8(b)可見(jiàn),當(dāng)蜂窩芯壁厚分別為0.035、0.07、0.105 mm時(shí),相應(yīng)的蜂窩板吸收全部能量的時(shí)長(zhǎng)依次為8.7、3.2、2.5 ms,吸能速率隨蜂窩芯層壁厚的增加而增大,并且還需指出的是,較薄的蜂窩芯層將導(dǎo)致蜂窩板剛度嚴(yán)重降低;由圖8(c)可見(jiàn),當(dāng)蜂窩芯高度分別為9、18、27 mm時(shí),相應(yīng)的蜂窩板吸收全部能量的時(shí)長(zhǎng)依次為4.0、3.2、2.9 ms,吸能速率隨蜂窩芯高度的增加而增大,蜂窩芯的高度變化對(duì)其吸能效率也有較明顯的影響;從表2中可以看出,蜂窩板吸收能量穩(wěn)定值隨蒙皮厚度或蜂窩芯壁厚的增加而減小,并隨蜂窩芯高度的增加而增大。這表明蜂窩板結(jié)構(gòu)參數(shù)對(duì)蜂窩板力學(xué)行為的影響在變化程度和變動(dòng)方向上均有體現(xiàn),因此實(shí)際工程中在對(duì)蜂窩結(jié)構(gòu)進(jìn)行優(yōu)化和使用時(shí),需謹(jǐn)慎考量材料參數(shù)、結(jié)構(gòu)性能和功能需求之間的關(guān)系。

    (b)蜂窩芯壁厚(蒙皮厚度、蜂窩芯高度分別為1.0、18 mm)

    (c)蜂窩芯高度(蒙皮厚度、蜂窩芯壁厚分別為1.0、0.07 mm)

    表2 吸收能量穩(wěn)定值

    3 結(jié)語(yǔ)

    在沖擊荷載作用下,金屬材料的應(yīng)變率范圍包含了低、中、高三種階段,而傳統(tǒng)的Johnson-Cook模型不適于體現(xiàn)這種寬應(yīng)變率條件下的金屬材料本構(gòu)關(guān)系,故本文以Ludwik準(zhǔn)則和Voce準(zhǔn)則線性組合作為新應(yīng)變硬化項(xiàng)、CS模型作為新應(yīng)變率效應(yīng)項(xiàng)對(duì)Johnson-Cook本構(gòu)模型進(jìn)行改造,獲得了能適用于寬應(yīng)變率條件的H/V-CS本構(gòu)模型?;诮饘僬鎽?yīng)力-應(yīng)變實(shí)驗(yàn),采用全局優(yōu)化算法確定了金屬材料的H/V-CS本構(gòu)模型參數(shù)并對(duì)模型的正確性進(jìn)行了驗(yàn)證,所得模型參數(shù)不依賴于初值,是全局最優(yōu)組合,具有準(zhǔn)確性高、方便易得和應(yīng)用便捷的優(yōu)點(diǎn),且不受參數(shù)數(shù)量限制。此外,還利用所建H/V-CS本構(gòu)模型并結(jié)合有限元分析研究了蜂窩板結(jié)構(gòu)參數(shù)對(duì)其吸能效果的影響。結(jié)果表明,改變蜂窩板的蒙皮厚度、芯層厚度和高度均會(huì)在不同程度上影響蜂窩板結(jié)構(gòu)的吸能效果,蜂窩板吸收能量穩(wěn)定值隨蒙皮或芯層厚度的增加而增大,并隨芯層高度的增加而減小,因此在實(shí)際工程中進(jìn)行多方案比較和優(yōu)化分析時(shí)應(yīng)綜合考慮材料使用功能及多種結(jié)構(gòu)參數(shù)的影響。

    猜你喜歡
    芯層蒙皮金屬材料
    一種適用于變彎度機(jī)翼后緣的蒙皮設(shè)計(jì)方法
    不同材質(zhì)客車頂蒙皮與前后圍蒙皮接縫處理方案
    一種耐高溫能力強(qiáng)的橡膠輸送帶
    運(yùn)載火箭框桁蒙皮結(jié)構(gòu)鉚接殼段多余物分析與控制
    把握考查角度 學(xué)好金屬材料
    飛機(jī)蒙皮上的幽默
    航空世界(2018年12期)2018-07-16 08:34:46
    空中爆炸載荷下梯度波紋夾層板抗爆性能仿真研究
    車用新型結(jié)構(gòu)的夾層構(gòu)件
    汽車文摘(2017年4期)2017-04-27 02:05:07
    負(fù)載型納米金屬材料的最新研究進(jìn)展
    為金屬材料披上“防護(hù)衣”——金屬材料的腐蝕與防護(hù)研究
    国产免费又黄又爽又色| 国产亚洲最大av| 国产免费一区二区三区四区乱码| 三上悠亚av全集在线观看| 日韩视频在线欧美| 精品国产露脸久久av麻豆| 婷婷色综合大香蕉| 国产精品女同一区二区软件| 麻豆乱淫一区二区| 人人妻人人爽人人添夜夜欢视频| 一级,二级,三级黄色视频| 女人精品久久久久毛片| av在线播放精品| 人妻系列 视频| 秋霞伦理黄片| 亚洲欧美一区二区三区黑人 | 丝袜喷水一区| 久久久久久大精品| 日本vs欧美在线观看视频| 国产激情欧美一区二区| 伊人久久大香线蕉亚洲五| 长腿黑丝高跟| 日韩有码中文字幕| 在线国产一区二区在线| 99国产极品粉嫩在线观看| 亚洲中文av在线| 午夜精品久久久久久毛片777| 精品久久久久久久久久免费视频 | 51午夜福利影视在线观看| 亚洲欧美激情综合另类| 免费在线观看视频国产中文字幕亚洲| 三上悠亚av全集在线观看| 欧美日本中文国产一区发布| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看免费高清a一片| 免费在线观看影片大全网站| 欧美黑人欧美精品刺激| 午夜免费激情av| 国内久久婷婷六月综合欲色啪| 亚洲国产精品合色在线| 国产精品久久久av美女十八| 一个人观看的视频www高清免费观看 | 90打野战视频偷拍视频| x7x7x7水蜜桃| 国产亚洲精品综合一区在线观看 | 午夜久久久在线观看| 日韩国内少妇激情av| 91成人精品电影| 亚洲狠狠婷婷综合久久图片| 十分钟在线观看高清视频www| 午夜福利影视在线免费观看| 一边摸一边做爽爽视频免费| 欧美激情 高清一区二区三区| 视频在线观看一区二区三区| 黄色怎么调成土黄色| 在线观看www视频免费| 一进一出好大好爽视频| a在线观看视频网站| 国内久久婷婷六月综合欲色啪| 精品福利观看| 亚洲精品久久成人aⅴ小说| 满18在线观看网站| 欧美中文综合在线视频| 夜夜看夜夜爽夜夜摸 | 亚洲国产中文字幕在线视频| 亚洲成人免费电影在线观看| 99国产精品免费福利视频| 欧美成狂野欧美在线观看| 精品国产乱子伦一区二区三区| 久久香蕉精品热| avwww免费| cao死你这个sao货| 成人18禁高潮啪啪吃奶动态图| 女同久久另类99精品国产91| 久久精品影院6| 亚洲人成电影免费在线| 亚洲专区字幕在线| 精品久久久久久,| 国产伦人伦偷精品视频| 高清在线国产一区| 女人被躁到高潮嗷嗷叫费观| 成在线人永久免费视频| 操美女的视频在线观看| 无限看片的www在线观看| 无遮挡黄片免费观看| 国产有黄有色有爽视频| 多毛熟女@视频| 久久香蕉精品热| 88av欧美| 国产成人一区二区三区免费视频网站| 香蕉久久夜色| 免费av毛片视频| 大香蕉久久成人网| 黄色怎么调成土黄色| 午夜精品国产一区二区电影| 91av网站免费观看| 高清在线国产一区| 91字幕亚洲| 国产一区二区三区在线臀色熟女 | 日韩欧美一区二区三区在线观看| 99在线视频只有这里精品首页| 精品国内亚洲2022精品成人| 久久国产乱子伦精品免费另类| 亚洲国产精品一区二区三区在线| 国产aⅴ精品一区二区三区波| 激情视频va一区二区三区| 免费高清视频大片| 亚洲一区高清亚洲精品| 少妇 在线观看| 黄片小视频在线播放| 久久久久国内视频| 99国产精品99久久久久| 99re在线观看精品视频| 久久性视频一级片| 男男h啪啪无遮挡| av中文乱码字幕在线| 欧美激情高清一区二区三区| 正在播放国产对白刺激| 精品电影一区二区在线| 在线看a的网站| 日韩欧美在线二视频| 亚洲欧美激情在线| 久久久久久久午夜电影 | 少妇被粗大的猛进出69影院| 在线观看日韩欧美| 日韩国内少妇激情av| 日本 av在线| a在线观看视频网站| 99精品在免费线老司机午夜| 男女下面进入的视频免费午夜 | 亚洲,欧美精品.| 欧美色视频一区免费| 精品高清国产在线一区| 国产精品香港三级国产av潘金莲| 后天国语完整版免费观看| 一级毛片高清免费大全| 一区福利在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩精品中文字幕看吧| 在线永久观看黄色视频| 美女 人体艺术 gogo| 中文字幕精品免费在线观看视频| 男女做爰动态图高潮gif福利片 | 人妻久久中文字幕网| 国产精品免费视频内射| 在线天堂中文资源库| 日韩一卡2卡3卡4卡2021年| 久久久久国产精品人妻aⅴ院| 80岁老熟妇乱子伦牲交| 黄色片一级片一级黄色片| 女人爽到高潮嗷嗷叫在线视频| 国产又色又爽无遮挡免费看| 午夜免费激情av| 一级毛片高清免费大全| 亚洲视频免费观看视频| 亚洲成a人片在线一区二区| 手机成人av网站| 亚洲男人天堂网一区| 久久久精品国产亚洲av高清涩受| 黑人巨大精品欧美一区二区mp4| 免费高清视频大片| 午夜久久久在线观看| 欧美亚洲日本最大视频资源| 日日干狠狠操夜夜爽| 国产高清videossex| 亚洲专区中文字幕在线| 男女下面插进去视频免费观看| 久久午夜综合久久蜜桃| 男女午夜视频在线观看| 18美女黄网站色大片免费观看| 亚洲 欧美一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 后天国语完整版免费观看| 搡老熟女国产l中国老女人| 国产高清激情床上av| 久久午夜亚洲精品久久| netflix在线观看网站| 美女午夜性视频免费| 国产成人影院久久av| 免费在线观看日本一区| 亚洲片人在线观看| 99在线视频只有这里精品首页| 欧美乱色亚洲激情| 午夜福利,免费看| 国产精品秋霞免费鲁丝片| 久久久久国产精品人妻aⅴ院| 男女床上黄色一级片免费看| 免费久久久久久久精品成人欧美视频| 亚洲全国av大片| 欧美成人性av电影在线观看| 人人澡人人妻人| 91成人精品电影| 高清毛片免费观看视频网站 | 国产精品日韩av在线免费观看 | 高清欧美精品videossex| 1024香蕉在线观看| 老司机午夜十八禁免费视频| 777久久人妻少妇嫩草av网站| 亚洲精品粉嫩美女一区| 亚洲成国产人片在线观看| 午夜两性在线视频| cao死你这个sao货| 国产成人精品久久二区二区免费| 1024香蕉在线观看| 亚洲成a人片在线一区二区| 欧美成人午夜精品| 亚洲欧美一区二区三区黑人| 淫妇啪啪啪对白视频| 久久人人精品亚洲av| 中国美女看黄片| 亚洲中文日韩欧美视频| 欧美黄色片欧美黄色片| 很黄的视频免费| 国产精品1区2区在线观看.| 午夜福利,免费看| 久久久国产成人精品二区 | 桃色一区二区三区在线观看| 日本黄色视频三级网站网址| 国产成人精品无人区| 夜夜躁狠狠躁天天躁| 在线视频色国产色| 日本一区二区免费在线视频| 麻豆成人av在线观看| 嫩草影院精品99| 亚洲精品国产色婷婷电影| 国产深夜福利视频在线观看| 91国产中文字幕| 久久精品亚洲精品国产色婷小说| 中文字幕人妻丝袜一区二区| 大码成人一级视频| 日韩欧美国产一区二区入口| 在线观看免费日韩欧美大片| 亚洲午夜精品一区,二区,三区| 国产三级黄色录像| 夜夜躁狠狠躁天天躁| 欧美在线黄色| 夜夜夜夜夜久久久久| 夜夜看夜夜爽夜夜摸 | 色综合婷婷激情| av视频免费观看在线观看| 日本撒尿小便嘘嘘汇集6| 黑人巨大精品欧美一区二区mp4| 国产av精品麻豆| 丰满迷人的少妇在线观看| 美女大奶头视频| 麻豆国产av国片精品| 老司机在亚洲福利影院| 日本免费a在线| 黄片播放在线免费| 国产一区在线观看成人免费| 国产亚洲欧美98| 亚洲精品在线美女| 一级,二级,三级黄色视频| 久久久久久久久中文| 少妇被粗大的猛进出69影院| 美女扒开内裤让男人捅视频| 免费看a级黄色片| 国产精品成人在线| 亚洲人成网站在线播放欧美日韩| 一级片免费观看大全| 亚洲午夜理论影院| 国产精品九九99| 黄色丝袜av网址大全| 亚洲欧美日韩无卡精品| 美女 人体艺术 gogo| 日本黄色视频三级网站网址| 一边摸一边抽搐一进一小说| 美女国产高潮福利片在线看| 琪琪午夜伦伦电影理论片6080| www.熟女人妻精品国产| 国产精品香港三级国产av潘金莲| 少妇的丰满在线观看| 99久久久亚洲精品蜜臀av| 老鸭窝网址在线观看| 精品卡一卡二卡四卡免费| 日本精品一区二区三区蜜桃| 在线播放国产精品三级| 国产主播在线观看一区二区| 性欧美人与动物交配| 国产欧美日韩一区二区三区在线| 91国产中文字幕| 国产精品一区二区精品视频观看| 亚洲狠狠婷婷综合久久图片| 国产精品电影一区二区三区| 99热只有精品国产| 国产人伦9x9x在线观看| 在线永久观看黄色视频| 免费av中文字幕在线| 亚洲成人免费av在线播放| 久9热在线精品视频| 又紧又爽又黄一区二区| a级毛片在线看网站| 高清欧美精品videossex| 成人免费观看视频高清| av国产精品久久久久影院| 国产蜜桃级精品一区二区三区| 成人国语在线视频| 乱人伦中国视频| tocl精华| 精品国产超薄肉色丝袜足j| 夜夜爽天天搞| 国产精品国产高清国产av| 丝袜美腿诱惑在线| 国产亚洲精品综合一区在线观看 | 高清欧美精品videossex| 琪琪午夜伦伦电影理论片6080| 人人澡人人妻人| 国产精品国产av在线观看| av网站在线播放免费| 婷婷精品国产亚洲av在线| 欧美日韩福利视频一区二区| 99久久久亚洲精品蜜臀av| 嫁个100分男人电影在线观看| 亚洲国产精品合色在线| 国产成人精品久久二区二区91| 亚洲午夜理论影院| 黑丝袜美女国产一区| 久久青草综合色| 亚洲第一青青草原| av天堂久久9| av欧美777| 婷婷丁香在线五月| 亚洲第一av免费看| 久久草成人影院| 日韩欧美国产一区二区入口| 很黄的视频免费| 久久人人97超碰香蕉20202| 久久久国产成人免费| 在线观看www视频免费| 视频在线观看一区二区三区| 亚洲五月色婷婷综合| 99国产精品一区二区蜜桃av| 欧美精品啪啪一区二区三区| 成在线人永久免费视频| 女人爽到高潮嗷嗷叫在线视频| 久久久国产精品麻豆| 中文字幕av电影在线播放| 亚洲精品美女久久久久99蜜臀| 久久久久精品国产欧美久久久| 美女午夜性视频免费| 中文字幕色久视频| 亚洲国产看品久久| www.精华液| 国产又色又爽无遮挡免费看| 亚洲国产中文字幕在线视频| av视频免费观看在线观看| 亚洲精品国产一区二区精华液| 亚洲第一av免费看| 精品久久久久久成人av| 精品一品国产午夜福利视频| 成熟少妇高潮喷水视频| 久久精品aⅴ一区二区三区四区| 精品久久蜜臀av无| 免费少妇av软件| 成年女人毛片免费观看观看9| 婷婷精品国产亚洲av在线| 日本wwww免费看| 精品电影一区二区在线| 美女 人体艺术 gogo| 伊人久久大香线蕉亚洲五| 一进一出抽搐动态| 伊人久久大香线蕉亚洲五| bbb黄色大片| 成人av一区二区三区在线看| 精品熟女少妇八av免费久了| 999久久久国产精品视频| 欧美人与性动交α欧美精品济南到| 精品一区二区三区视频在线观看免费 | 国产成人av激情在线播放| 18禁国产床啪视频网站| 亚洲在线自拍视频| 亚洲精品美女久久av网站| 黄色a级毛片大全视频| 亚洲成av片中文字幕在线观看| 精品一区二区三区av网在线观看| 免费高清视频大片| 看黄色毛片网站| 久久久久久久午夜电影 | av网站免费在线观看视频| 亚洲色图av天堂| 亚洲aⅴ乱码一区二区在线播放 | 18禁黄网站禁片午夜丰满| 亚洲专区中文字幕在线| 三级毛片av免费| 波多野结衣高清无吗| 国产成人av教育| 99国产极品粉嫩在线观看| 97超级碰碰碰精品色视频在线观看| 亚洲av成人不卡在线观看播放网| 桃色一区二区三区在线观看| 美女午夜性视频免费| 国产激情久久老熟女| 在线观看免费视频网站a站| 亚洲男人的天堂狠狠| 男女下面进入的视频免费午夜 | 又大又爽又粗| 久久婷婷成人综合色麻豆| 在线观看日韩欧美| 亚洲少妇的诱惑av| 色综合欧美亚洲国产小说| 免费少妇av软件| 免费久久久久久久精品成人欧美视频| 久久精品国产清高在天天线| 免费观看精品视频网站| 午夜福利免费观看在线| 美国免费a级毛片| 欧美激情高清一区二区三区| 9热在线视频观看99| 母亲3免费完整高清在线观看| 久久国产精品影院| 麻豆国产av国片精品| 亚洲全国av大片| 性少妇av在线| 亚洲国产看品久久| 成人18禁高潮啪啪吃奶动态图| 天堂动漫精品| 色播在线永久视频| 免费女性裸体啪啪无遮挡网站| www日本在线高清视频| 亚洲国产欧美一区二区综合| 999久久久精品免费观看国产| 久9热在线精品视频| 欧美+亚洲+日韩+国产| 两人在一起打扑克的视频| 成人三级黄色视频| 日韩欧美在线二视频| 18禁国产床啪视频网站| 中文字幕av电影在线播放| 亚洲狠狠婷婷综合久久图片| 欧美激情久久久久久爽电影 | xxx96com| 久久香蕉国产精品| 国产伦一二天堂av在线观看| 亚洲黑人精品在线| 亚洲中文字幕日韩| 在线观看免费午夜福利视频| 日韩成人在线观看一区二区三区| 黑人猛操日本美女一级片| 日韩欧美三级三区| 一级毛片高清免费大全| 日本免费a在线| 欧美不卡视频在线免费观看 | 天天躁夜夜躁狠狠躁躁| 少妇的丰满在线观看| 男女做爰动态图高潮gif福利片 | 两性夫妻黄色片| 激情视频va一区二区三区| 久99久视频精品免费| 久久久久亚洲av毛片大全| 可以免费在线观看a视频的电影网站| 99re在线观看精品视频| 丰满迷人的少妇在线观看| 黄色视频不卡| 欧美性长视频在线观看| 中国美女看黄片| 悠悠久久av| 两人在一起打扑克的视频| 国产不卡一卡二| 最新美女视频免费是黄的| 亚洲精品粉嫩美女一区| 亚洲自偷自拍图片 自拍| 色精品久久人妻99蜜桃| 亚洲成人久久性| 亚洲人成伊人成综合网2020| 美女午夜性视频免费| 制服人妻中文乱码| videosex国产| 久久久久国产精品人妻aⅴ院| 日韩大码丰满熟妇| 亚洲 国产 在线| 在线av久久热| 国产免费av片在线观看野外av| 两个人看的免费小视频| 久久精品人人爽人人爽视色| 757午夜福利合集在线观看| 久久精品91无色码中文字幕| 三上悠亚av全集在线观看| 中文字幕色久视频| 欧美黄色片欧美黄色片| 91九色精品人成在线观看| 国产成人精品无人区| 黄网站色视频无遮挡免费观看| 黄色片一级片一级黄色片| 啦啦啦免费观看视频1| 最近最新中文字幕大全电影3 | www.自偷自拍.com| 成人亚洲精品av一区二区 | 亚洲一区二区三区不卡视频| 一级片免费观看大全| 免费在线观看影片大全网站| 成年女人毛片免费观看观看9| 91麻豆精品激情在线观看国产 | 亚洲一码二码三码区别大吗| 成人特级黄色片久久久久久久| 久久这里只有精品19| 一区福利在线观看| 国产精品电影一区二区三区| 欧美成人性av电影在线观看| 欧美亚洲日本最大视频资源| 久久精品国产99精品国产亚洲性色 | 最好的美女福利视频网| 多毛熟女@视频| 国产一区二区在线av高清观看| 亚洲va日本ⅴa欧美va伊人久久| 桃红色精品国产亚洲av| 神马国产精品三级电影在线观看 | 亚洲熟妇熟女久久| 色婷婷久久久亚洲欧美| 成人手机av| 亚洲美女黄片视频| 精品第一国产精品| 成在线人永久免费视频| 欧美成人午夜精品| e午夜精品久久久久久久| 久久精品国产亚洲av香蕉五月| 久久青草综合色| 亚洲欧美一区二区三区黑人| 妹子高潮喷水视频| 十分钟在线观看高清视频www| 1024香蕉在线观看| 国产熟女xx| 日韩高清综合在线| 国产av在哪里看| 免费av中文字幕在线| 色综合站精品国产| 一级a爱视频在线免费观看| 中文字幕人妻熟女乱码| 男女下面进入的视频免费午夜 | 美女 人体艺术 gogo| 久久精品aⅴ一区二区三区四区| 一级毛片精品| 久久久久久久久久久久大奶| 嫩草影视91久久| 免费av中文字幕在线| 国产免费现黄频在线看| 亚洲欧美日韩另类电影网站| 久久久精品国产亚洲av高清涩受| 在线永久观看黄色视频| 韩国av一区二区三区四区| 99riav亚洲国产免费| 新久久久久国产一级毛片| 制服人妻中文乱码| 欧美日韩乱码在线| 色哟哟哟哟哟哟| av网站在线播放免费| 国产午夜精品久久久久久| 亚洲一区中文字幕在线| 精品熟女少妇八av免费久了| 国产主播在线观看一区二区| 久久久久久久久中文| 视频在线观看一区二区三区| 亚洲精品一区av在线观看| 国产高清videossex| 国产av精品麻豆| 日本黄色视频三级网站网址| 可以免费在线观看a视频的电影网站| 亚洲久久久国产精品| 老司机深夜福利视频在线观看| 久久精品亚洲av国产电影网| 亚洲成人免费av在线播放| 50天的宝宝边吃奶边哭怎么回事| 最近最新中文字幕大全免费视频| 欧美一区二区精品小视频在线| 法律面前人人平等表现在哪些方面| 国产免费现黄频在线看| 亚洲第一青青草原| 欧美性长视频在线观看| 操出白浆在线播放| 少妇被粗大的猛进出69影院| 香蕉久久夜色| 真人一进一出gif抽搐免费| 国产不卡一卡二| 一区二区三区国产精品乱码| 色精品久久人妻99蜜桃| 午夜影院日韩av| 高潮久久久久久久久久久不卡| 最近最新中文字幕大全免费视频| 精品电影一区二区在线| 国产av一区在线观看免费| 亚洲九九香蕉| 丝袜在线中文字幕| 欧美日韩亚洲国产一区二区在线观看| 亚洲av日韩精品久久久久久密| 在线观看免费视频日本深夜| 国产亚洲欧美98| 亚洲av成人av| av免费在线观看网站| 夜夜夜夜夜久久久久| 一二三四社区在线视频社区8| 在线十欧美十亚洲十日本专区| 中文字幕人妻熟女乱码| 丰满的人妻完整版| 国产免费男女视频| 欧美人与性动交α欧美精品济南到| 久久精品aⅴ一区二区三区四区| 黑人操中国人逼视频| tocl精华| 国产av一区二区精品久久| 一本综合久久免费| 欧美日韩一级在线毛片| 婷婷六月久久综合丁香| 亚洲午夜理论影院| 如日韩欧美国产精品一区二区三区| 国产一区在线观看成人免费| 久久久久久久久中文| 久久九九热精品免费| 精品一品国产午夜福利视频| 韩国av一区二区三区四区| 麻豆国产av国片精品| 另类亚洲欧美激情| 国产成+人综合+亚洲专区| 在线观看免费午夜福利视频| 18禁国产床啪视频网站| 色老头精品视频在线观看| 97超级碰碰碰精品色视频在线观看|