• <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福利片在线观看| 乱码一卡2卡4卡精品| 大香蕉97超碰在线| 美女黄网站色视频| 国产精品一区二区在线观看99 | 少妇被粗大猛烈的视频| 人妻制服诱惑在线中文字幕| 亚洲美女搞黄在线观看| 日本av手机在线免费观看| 26uuu在线亚洲综合色| 网址你懂的国产日韩在线| 老女人水多毛片| 亚洲不卡免费看| 在线免费十八禁| 亚洲国产欧美人成| 国产精品久久久久久精品电影| kizo精华| 韩国高清视频一区二区三区| 视频中文字幕在线观看| 欧美日韩精品成人综合77777| 国产精品.久久久| 男女那种视频在线观看| 国产精品美女特级片免费视频播放器| 久久久精品94久久精品| 美女黄网站色视频| 三级国产精品欧美在线观看| 看十八女毛片水多多多| 欧美高清性xxxxhd video| 性色avwww在线观看| 国产成人福利小说| 免费看a级黄色片| 国产精品乱码一区二三区的特点| 国产精品蜜桃在线观看| 日韩高清综合在线| 精品久久久久久久久av| 亚洲成人精品中文字幕电影| 亚洲经典国产精华液单| 欧美一级a爱片免费观看看| 国产亚洲一区二区精品| 少妇猛男粗大的猛烈进出视频 | 最近2019中文字幕mv第一页| 校园人妻丝袜中文字幕| 最新中文字幕久久久久| 日本-黄色视频高清免费观看| 亚洲国产日韩欧美精品在线观看| 欧美高清性xxxxhd video| 99久国产av精品国产电影| 噜噜噜噜噜久久久久久91| 成人国产麻豆网| 我的老师免费观看完整版| 国产精品99久久久久久久久| 级片在线观看| 国产美女午夜福利| 久久久久久国产a免费观看| 亚洲国产精品sss在线观看| 欧美最新免费一区二区三区| 美女被艹到高潮喷水动态| 成人鲁丝片一二三区免费| 男插女下体视频免费在线播放| 精品午夜福利在线看| 亚洲精品亚洲一区二区| 国产精品一区二区三区四区免费观看| 成人美女网站在线观看视频| 久久久久性生活片| 一个人观看的视频www高清免费观看| 国产高清视频在线观看网站| 国产午夜精品久久久久久一区二区三区| 国产精品一区二区三区四区免费观看| 亚洲三级黄色毛片| 老司机福利观看| 一本一本综合久久| 少妇裸体淫交视频免费看高清| 色吧在线观看| 欧美日本亚洲视频在线播放| 中文亚洲av片在线观看爽| 欧美高清成人免费视频www| 日本三级黄在线观看| 中文字幕精品亚洲无线码一区| 国产免费视频播放在线视频 | 91av网一区二区| 一区二区三区高清视频在线| 网址你懂的国产日韩在线| 啦啦啦韩国在线观看视频| 亚洲国产成人一精品久久久| 成人av在线播放网站| 六月丁香七月| 听说在线观看完整版免费高清| 亚洲伊人久久精品综合 | 美女xxoo啪啪120秒动态图| 美女内射精品一级片tv| 国产精品美女特级片免费视频播放器| 天天躁日日操中文字幕| 丰满乱子伦码专区| 我的女老师完整版在线观看| 汤姆久久久久久久影院中文字幕 | 国产探花在线观看一区二区| 青春草视频在线免费观看| 女人被狂操c到高潮| 亚洲国产精品成人久久小说| 国产免费视频播放在线视频 | 久久久精品94久久精品| 亚洲国产欧洲综合997久久,| 99热网站在线观看| 日本wwww免费看| 久久精品夜夜夜夜夜久久蜜豆| av国产免费在线观看| 村上凉子中文字幕在线| 国产午夜精品久久久久久一区二区三区| 亚洲av福利一区| 听说在线观看完整版免费高清| 免费看美女性在线毛片视频| av国产久精品久网站免费入址| 国产亚洲精品久久久com| 亚洲欧美精品综合久久99| 99久久人妻综合| 直男gayav资源| 精品99又大又爽又粗少妇毛片| 女的被弄到高潮叫床怎么办| 亚洲不卡免费看| 欧美3d第一页| 久久精品夜色国产| 日韩欧美精品免费久久| 精品欧美国产一区二区三| 插阴视频在线观看视频| 久久韩国三级中文字幕| 日本免费a在线| 男女国产视频网站| 中文乱码字字幕精品一区二区三区 | 深夜a级毛片| 一个人看的www免费观看视频| 欧美潮喷喷水| 欧美成人午夜免费资源| 国产人妻一区二区三区在| 亚洲精品影视一区二区三区av| 噜噜噜噜噜久久久久久91| 国产精品人妻久久久影院| 婷婷色综合大香蕉| 日本五十路高清| 夜夜看夜夜爽夜夜摸| 国产毛片a区久久久久| 亚洲性久久影院| 免费一级毛片在线播放高清视频| 亚洲av成人精品一区久久| 亚洲天堂国产精品一区在线| 国产91av在线免费观看| 又爽又黄无遮挡网站| 干丝袜人妻中文字幕| 麻豆精品久久久久久蜜桃| 亚洲欧美中文字幕日韩二区| 国产亚洲精品久久久com| av国产久精品久网站免费入址| 欧美zozozo另类| 超碰av人人做人人爽久久| 国产一区二区在线av高清观看| 久久久久久伊人网av| 国产不卡一卡二| 人人妻人人澡欧美一区二区| 日韩欧美 国产精品| 亚洲色图av天堂| 亚洲最大成人中文| 久久久久免费精品人妻一区二区| 国产精品电影一区二区三区| 国产乱人偷精品视频| 亚洲国产精品久久男人天堂| 午夜精品国产一区二区电影 | 99久久成人亚洲精品观看| 亚洲国产精品久久男人天堂| 日韩视频在线欧美| 你懂的网址亚洲精品在线观看 | www.色视频.com| 国产高清不卡午夜福利| 亚洲欧美日韩东京热| 亚洲在久久综合| 国产伦在线观看视频一区| 午夜视频国产福利| 欧美变态另类bdsm刘玥| 久久这里有精品视频免费| 亚洲丝袜综合中文字幕| 国产欧美日韩精品一区二区| 国产激情偷乱视频一区二区| 午夜精品一区二区三区免费看| 国产一区有黄有色的免费视频 | 亚洲av免费在线观看| 免费播放大片免费观看视频在线观看 | 成人高潮视频无遮挡免费网站| 亚洲人成网站在线观看播放| 国产日韩欧美在线精品| 久久久久久久久久久丰满| 成人亚洲精品av一区二区| 国产精品伦人一区二区| 麻豆av噜噜一区二区三区| 国产亚洲一区二区精品| 国产一级毛片在线| 国产单亲对白刺激| 成人鲁丝片一二三区免费| 蜜桃亚洲精品一区二区三区| 欧美极品一区二区三区四区| 亚洲国产精品久久男人天堂| 乱码一卡2卡4卡精品| 夜夜爽夜夜爽视频| 国产伦精品一区二区三区四那| 国产激情偷乱视频一区二区| 久久久久久久久中文| 中国美白少妇内射xxxbb| 久久久久免费精品人妻一区二区| 在线播放无遮挡| 亚洲精华国产精华液的使用体验| 联通29元200g的流量卡| 亚洲在线观看片| 亚洲三级黄色毛片| 精品国内亚洲2022精品成人| 国产爱豆传媒在线观看| 神马国产精品三级电影在线观看| 亚洲精品日韩av片在线观看| 久久精品夜夜夜夜夜久久蜜豆| 菩萨蛮人人尽说江南好唐韦庄 | 有码 亚洲区| 国产亚洲午夜精品一区二区久久 | 亚洲图色成人| 亚洲国产精品国产精品| 亚洲在线观看片| 小蜜桃在线观看免费完整版高清| 国产午夜精品论理片| 成年版毛片免费区| 日本黄色视频三级网站网址| 又粗又爽又猛毛片免费看| 亚洲欧美日韩无卡精品| 婷婷六月久久综合丁香| 日韩三级伦理在线观看| 精品久久久久久久久av| 欧美性猛交黑人性爽| 舔av片在线| 亚洲av一区综合| 久久99热这里只频精品6学生 | 国产精品久久久久久精品电影小说 | 亚洲精品久久久久久婷婷小说 | 日本猛色少妇xxxxx猛交久久| 成人特级av手机在线观看| 深爱激情五月婷婷| 亚洲aⅴ乱码一区二区在线播放| 天天躁日日操中文字幕| 日韩高清综合在线| 亚洲av免费高清在线观看| 麻豆成人午夜福利视频| 日本黄色片子视频| 蜜桃亚洲精品一区二区三区| 精品久久久久久久久久久久久| 99热这里只有是精品50| 久久精品熟女亚洲av麻豆精品 | 成人午夜高清在线视频| 能在线免费看毛片的网站| 欧美成人午夜免费资源| 国产精品久久久久久久电影| 神马国产精品三级电影在线观看| 久久国产乱子免费精品| 亚洲国产欧美在线一区| 偷拍熟女少妇极品色| 日韩亚洲欧美综合| 在线免费观看不下载黄p国产| 青春草国产在线视频| 亚洲欧美精品综合久久99| 久久精品国产亚洲av天美| 直男gayav资源| 精品国内亚洲2022精品成人| 久久亚洲国产成人精品v| 国产精品一区二区三区四区免费观看| av.在线天堂| 久久99蜜桃精品久久| 亚洲av.av天堂| 一级二级三级毛片免费看| 一本久久精品| 麻豆av噜噜一区二区三区| 国产av一区在线观看免费| 青春草视频在线免费观看| 色5月婷婷丁香| 亚洲经典国产精华液单| 久久热精品热| 久久精品夜夜夜夜夜久久蜜豆| 男女边吃奶边做爰视频| 69av精品久久久久久| 国产精品一区www在线观看| 内地一区二区视频在线| 少妇的逼水好多| 亚洲成色77777| 久久热精品热| 在线观看一区二区三区| 国产高清视频在线观看网站| 最近视频中文字幕2019在线8| 91在线精品国自产拍蜜月| 男插女下体视频免费在线播放| 秋霞在线观看毛片| 日本三级黄在线观看| kizo精华| 永久网站在线| 亚洲怡红院男人天堂| or卡值多少钱| 最近手机中文字幕大全| 国产视频内射| 最近最新中文字幕大全电影3| 看黄色毛片网站| 亚洲经典国产精华液单| 少妇人妻一区二区三区视频| 精品午夜福利在线看| av卡一久久| 成人二区视频| 少妇的逼水好多| 久久精品91蜜桃| 久久这里只有精品中国| 在线观看美女被高潮喷水网站| 国语对白做爰xxxⅹ性视频网站| 少妇熟女欧美另类| 级片在线观看| 一个人看的www免费观看视频| 永久网站在线| 免费观看a级毛片全部| 国产高清视频在线观看网站| 狠狠狠狠99中文字幕| 国产美女午夜福利| 高清在线视频一区二区三区 | 熟女电影av网| 少妇丰满av| 嫩草影院精品99| 免费观看的影片在线观看| 久久精品国产99精品国产亚洲性色| 天堂网av新在线| 成年版毛片免费区| 亚洲乱码一区二区免费版| 国产亚洲91精品色在线| 国产精品久久视频播放| 久久草成人影院| 国产伦理片在线播放av一区| 成人毛片a级毛片在线播放| 亚洲美女视频黄频| 欧美性感艳星| 午夜精品在线福利| 精品久久久久久成人av| 国产午夜精品论理片| 中文天堂在线官网| 97超视频在线观看视频| 亚洲五月天丁香| 亚洲av不卡在线观看| kizo精华| 欧美一区二区精品小视频在线| 亚洲国产精品专区欧美| 一级毛片久久久久久久久女| 亚洲国产色片| 国产精品麻豆人妻色哟哟久久 | av专区在线播放| 你懂的网址亚洲精品在线观看 | 色视频www国产| 中文字幕av在线有码专区| 亚洲欧美精品自产自拍| 日本wwww免费看| 日本av手机在线免费观看| 久久精品国产亚洲网站| 黄色欧美视频在线观看| 在线a可以看的网站| 热99在线观看视频| 国产免费视频播放在线视频 | 人妻少妇偷人精品九色| 日本一本二区三区精品| 欧美日韩综合久久久久久| 亚洲欧美成人综合另类久久久 | 国产精品电影一区二区三区| 全区人妻精品视频| 婷婷六月久久综合丁香| 亚洲国产精品成人综合色| 亚洲av不卡在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 自拍偷自拍亚洲精品老妇| 久久精品综合一区二区三区| 国产精华一区二区三区| 亚洲美女搞黄在线观看| 国产精品国产高清国产av| 免费观看a级毛片全部| 嫩草影院新地址| 极品教师在线视频| 久久精品国产亚洲av涩爱| 欧美日韩在线观看h| or卡值多少钱| 能在线免费观看的黄片| 久久久精品大字幕| 国内精品美女久久久久久| 国产精品一区二区三区四区久久| 最近手机中文字幕大全| 国产综合懂色| 国产av码专区亚洲av| 国产黄片视频在线免费观看| 欧美3d第一页| 伊人久久精品亚洲午夜| 亚洲欧美日韩无卡精品| 我要看日韩黄色一级片| 午夜激情欧美在线| 亚洲成av人片在线播放无| 午夜福利在线在线| 久久综合国产亚洲精品| 成人二区视频| 国产成人freesex在线| 亚洲精品日韩av片在线观看| 深爱激情五月婷婷| 桃色一区二区三区在线观看| 欧美激情国产日韩精品一区| 亚洲精品日韩av片在线观看| 国产高清有码在线观看视频| 成年av动漫网址| 欧美日韩综合久久久久久| 高清日韩中文字幕在线| 能在线免费看毛片的网站| 高清在线视频一区二区三区 | 久久99精品国语久久久| 国产精品女同一区二区软件| 如何舔出高潮| 三级男女做爰猛烈吃奶摸视频| 国产精品,欧美在线| 美女cb高潮喷水在线观看| 亚洲在线观看片| 三级国产精品欧美在线观看| 高清av免费在线| 一级黄片播放器| av在线亚洲专区| 男女下面进入的视频免费午夜| 18禁在线无遮挡免费观看视频| 国产在视频线在精品| 亚洲中文字幕一区二区三区有码在线看| 国产一区亚洲一区在线观看| 国语对白做爰xxxⅹ性视频网站| 中文字幕人妻熟人妻熟丝袜美| 永久网站在线| a级毛片免费高清观看在线播放| 成人鲁丝片一二三区免费| 国产精品一区二区三区四区免费观看| www.色视频.com| 亚洲av免费高清在线观看| 精品国内亚洲2022精品成人| 亚洲av免费在线观看| 亚洲国产精品成人久久小说| 日韩三级伦理在线观看| 最近视频中文字幕2019在线8| 国产精品国产三级国产专区5o | 99在线人妻在线中文字幕| 国产淫片久久久久久久久| 欧美变态另类bdsm刘玥| 亚洲成人精品中文字幕电影| 色播亚洲综合网| 又黄又爽又刺激的免费视频.| 午夜福利网站1000一区二区三区| 亚洲国产欧美人成| 免费在线观看成人毛片| 久久久久久久久久成人| 婷婷色综合大香蕉| 国模一区二区三区四区视频| 十八禁国产超污无遮挡网站| 久久久亚洲精品成人影院| 日韩成人av中文字幕在线观看| 久久久久久久国产电影| 精品熟女少妇av免费看| 久久久精品94久久精品| 欧美丝袜亚洲另类| 一个人看视频在线观看www免费| 黄片wwwwww| 一区二区三区乱码不卡18| 亚洲av中文字字幕乱码综合| 亚洲欧美精品专区久久| av卡一久久| 国产欧美日韩精品一区二区| 黄片无遮挡物在线观看| 99久久无色码亚洲精品果冻| 国产色爽女视频免费观看| 中文天堂在线官网| 一级二级三级毛片免费看| 直男gayav资源| 欧美精品一区二区大全| 嫩草影院新地址| 国产 一区 欧美 日韩| 日韩av在线免费看完整版不卡| 久久99热这里只频精品6学生 | 日日摸夜夜添夜夜添av毛片| 在线观看一区二区三区| 2021少妇久久久久久久久久久| 成人午夜精彩视频在线观看| av国产免费在线观看| 2021少妇久久久久久久久久久| 97超碰精品成人国产| 99久久精品热视频| 最近2019中文字幕mv第一页| 国产毛片a区久久久久| 国产一区亚洲一区在线观看| 少妇熟女欧美另类| 亚洲色图av天堂| 久久久久久国产a免费观看| 亚洲国产成人一精品久久久| 人妻夜夜爽99麻豆av| 嫩草影院新地址| 欧美97在线视频| 国产一区二区三区av在线| 成人无遮挡网站| 秋霞伦理黄片| 精品国产三级普通话版| 水蜜桃什么品种好| 成年版毛片免费区| 一级毛片电影观看 | 亚洲人成网站在线观看播放| 久久久久精品久久久久真实原创| 在现免费观看毛片| 日韩国内少妇激情av| 亚洲熟妇中文字幕五十中出| 成人二区视频| 三级国产精品欧美在线观看| 亚洲欧美日韩卡通动漫| 亚洲丝袜综合中文字幕| 男人舔女人下体高潮全视频| 在线免费十八禁| 女人久久www免费人成看片 | 国产色爽女视频免费观看| 久久99蜜桃精品久久| 国产精品嫩草影院av在线观看| 亚洲最大成人中文| 99在线人妻在线中文字幕| 别揉我奶头 嗯啊视频| 久久久午夜欧美精品| 禁无遮挡网站| 亚洲欧美日韩东京热| 精品久久国产蜜桃| 亚洲人成网站在线播| 亚洲精品自拍成人| 久久午夜福利片| 搞女人的毛片| 成年女人永久免费观看视频| 色尼玛亚洲综合影院| 只有这里有精品99| 国产淫片久久久久久久久| h日本视频在线播放| 日韩强制内射视频| 男女那种视频在线观看| 少妇被粗大猛烈的视频| 亚洲人与动物交配视频| 永久网站在线| 欧美xxxx黑人xx丫x性爽| 欧美+日韩+精品| 国产一区二区在线av高清观看| 99久久人妻综合| 美女大奶头视频| 久久久久性生活片| 边亲边吃奶的免费视频| 能在线免费看毛片的网站| 亚洲在久久综合| 一级毛片电影观看 | 亚洲在久久综合| 九九热线精品视视频播放| 男女那种视频在线观看| 国产精品av视频在线免费观看| 美女高潮的动态| 又粗又爽又猛毛片免费看| 插阴视频在线观看视频| 伦理电影大哥的女人| 一级毛片久久久久久久久女| 天堂av国产一区二区熟女人妻| 美女黄网站色视频| 成年av动漫网址| 日韩欧美精品v在线| 热99re8久久精品国产| 精品熟女少妇av免费看| 六月丁香七月| 国内精品一区二区在线观看| 噜噜噜噜噜久久久久久91| 日韩视频在线欧美| 国产精品,欧美在线| 国产精品久久电影中文字幕| 国产一区二区三区av在线| 三级毛片av免费| 国内精品美女久久久久久| 一本久久精品| 人体艺术视频欧美日本| 国产精品久久久久久av不卡| 亚洲,欧美,日韩| 人妻制服诱惑在线中文字幕| 亚洲精品国产成人久久av| 国产精品电影一区二区三区| 亚洲av福利一区| 国产一区有黄有色的免费视频 | 亚洲欧美成人精品一区二区| 亚洲国产欧美在线一区| 国产大屁股一区二区在线视频| 少妇丰满av| 一卡2卡三卡四卡精品乱码亚洲| 日本色播在线视频| 美女高潮的动态| 国产亚洲5aaaaa淫片| 简卡轻食公司| 国内少妇人妻偷人精品xxx网站| 日韩亚洲欧美综合| 亚洲国产高清在线一区二区三| 好男人视频免费观看在线| 天天一区二区日本电影三级| 日本爱情动作片www.在线观看| 国产午夜精品论理片| 成人漫画全彩无遮挡| 在线播放无遮挡| 91久久精品国产一区二区三区| 久久这里有精品视频免费| 亚洲内射少妇av| 成人一区二区视频在线观看| 精品午夜福利在线看| 国产一区有黄有色的免费视频 | 看免费成人av毛片| 国产黄色小视频在线观看| 久久精品国产亚洲av涩爱| 欧美一区二区亚洲| 麻豆成人av视频| 午夜亚洲福利在线播放| 成人高潮视频无遮挡免费网站| 69人妻影院|