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

    Mg-Gd-Y-Zr合金熱變形本構(gòu)方程

    2011-11-30 02:00:12張新明吳懿萍鄧運(yùn)來(lái)唐昌平
    關(guān)鍵詞:激活能雙曲指數(shù)函數(shù)

    張新明,吳懿萍,鄧運(yùn)來(lái),唐昌平

    (中南大學(xué) 材料科學(xué)與工程學(xué)院,長(zhǎng)沙 410083)

    Mg-Gd-Y-Zr合金熱變形本構(gòu)方程

    張新明,吳懿萍,鄧運(yùn)來(lái),唐昌平

    (中南大學(xué) 材料科學(xué)與工程學(xué)院,長(zhǎng)沙 410083)

    采用Geeble?1500熱模擬實(shí)驗(yàn)機(jī)測(cè)試了高強(qiáng)耐熱Mg-6Gd-3Y-0.5Zr合金在變形溫度為300~500 ℃、應(yīng)變速率為10?3~1 s?1下的流變力學(xué)行為,采用掃描電子顯微鏡對(duì)其微觀組織進(jìn)行觀察,分析了冪函數(shù)(PI)、指數(shù)函數(shù)(EI)和雙曲正弦函數(shù)(SI)半經(jīng)驗(yàn)本構(gòu)方程對(duì)該合金變形行為擬合的適用性。結(jié)果表明:EI擬合實(shí)驗(yàn)結(jié)果精度大于PI和SI的,即使對(duì)SI函數(shù)中材料常數(shù)α進(jìn)行一定優(yōu)化處理得到SIO函數(shù)的擬合精度,也與EI函數(shù)十分接近。合金中高溫耐熱相提高合金高溫下的強(qiáng)度是EI擬合優(yōu)于PI和SI的原因。

    Mg-Gd-Y-Zr合金;熱壓縮變形;本構(gòu)方程;材料常數(shù)優(yōu)化

    高強(qiáng)耐熱鎂合金在航空航天中應(yīng)用越來(lái)越廣泛[1],但鎂合金屬于HCP結(jié)構(gòu)金屬,室溫條件下伸長(zhǎng)率偏低,塑性加工性能差,故實(shí)際應(yīng)用的鎂合金材料大都采用熱變形,而流變力學(xué)規(guī)律與微觀組織形成機(jī)理是熱加工的基礎(chǔ),因而受到廣泛關(guān)注[2?4]。

    研究鎂合金熱變形行為時(shí),一般使用冪函數(shù)、指數(shù)函數(shù)和雙曲正弦函數(shù)這3種本構(gòu)方程對(duì)其流變力學(xué)規(guī)律進(jìn)行表征。其中,在研究ZK60合金的熱壓縮變形行為時(shí),本文作者指出在變形溫度分別為423 K和423~723 K范圍內(nèi)宜對(duì)應(yīng)采用指數(shù)函數(shù)和冪函數(shù)形式本構(gòu)方程進(jìn)行擬合,則兩者都能夠得到最佳結(jié)果[5]。GUO等[6?7]在研究 AZ31合金熱壓縮流變力學(xué)行為時(shí)指出,低于 350 ℃下時(shí)適合用指數(shù)函數(shù)擬合,高于350 ℃則適合用冪函數(shù)擬合。同時(shí),也有作者指出分別使用冪函數(shù)和雙曲正弦函數(shù)擬合AZ31合金熱壓縮變形行為時(shí)得到的擬合結(jié)果相差不大[8]。由上可見(jiàn),擬合同種ZK60合金的熱壓縮變形行為時(shí),在不同的變形條件下宜采用不同類(lèi)型本構(gòu)方程才能得到各自最佳的擬合結(jié)果;在擬合AZ31合金熱壓縮變形行為時(shí)本構(gòu)方程采用方案也有所不同。不同于ZK和AZ系列合金,Mg-Gd-Y-X系合金高溫性能優(yōu)異[9],熱壓縮變形時(shí)比前兩種合金流變應(yīng)力較大,文獻(xiàn)[3, 10]中采用雙曲正弦函數(shù)擬合該系列合金熱變形行為,其中材料常數(shù)(α)設(shè)為冪函數(shù)硬化指數(shù)(n)與指數(shù)函數(shù)硬化系數(shù)(β)的比值(α =β/n ),而至今很少有研究涉及對(duì)比冪函數(shù)、指數(shù)函數(shù)和雙曲正弦函數(shù)3種類(lèi)型本構(gòu)方程擬合Mg-Gd-Y-X系合金熱變形行為時(shí)的優(yōu)劣。

    本文作者以一種 Mg-Gd-Y-Zr合金為試驗(yàn)材料,采用冪函數(shù)、指數(shù)函數(shù)和雙曲正弦函數(shù)3種類(lèi)型本構(gòu)方程擬合了該合金熱壓縮流變行為,綜合分析了這 3種函數(shù)擬合時(shí)的優(yōu)劣,并對(duì)雙曲正弦函數(shù)的材料常數(shù)α進(jìn)行優(yōu)化處理以提高其擬合精度。在此基礎(chǔ)上,結(jié)合組織掃描討論了本構(gòu)方程擬合該合金熱變形行為不同適應(yīng)性的微觀機(jī)理,期望能對(duì)揭示Mg-Gd-Y-X系合金的熱變形規(guī)律提供有益參考。

    1 實(shí)驗(yàn)

    熱壓縮試驗(yàn)材料為自制鑄態(tài)Mg-6Gd-3Y-0.5Zr稀土鎂合金(簡(jiǎn)稱(chēng)GW63)。在熱壓縮前鑄錠經(jīng)過(guò)520 ℃、24 h均勻化處理。熱壓縮試樣(d 10 mm×15 mm)用線(xiàn)切割機(jī)進(jìn)行切取。熱壓縮試驗(yàn)在Gleeble?1500熱模擬試驗(yàn)機(jī)上進(jìn)行:變形溫度為300~500 ℃,應(yīng)變速率為10?3~1 s?1,最大真應(yīng)變?yōu)?0.6。熱壓縮前試樣兩端涂抹含石墨的固體潤(rùn)滑劑以減小變形時(shí)的摩擦。熱壓縮時(shí)1 min內(nèi)升溫到變形溫度并保溫3 min,以盡量保證試樣壓縮時(shí)溫度均勻。所有試樣熱壓縮變形完成后迅速水冷以固定高溫變形結(jié)束時(shí)的顯微組織。熱壓縮變形后的試樣沿垂直于壓縮軸剖切,經(jīng)粗磨、精磨、拋光和腐蝕后在 Sirion200場(chǎng)發(fā)射掃描電子顯微鏡下觀察微觀組織,并分析合金中耐熱相成分。根據(jù)計(jì)算機(jī)采集的載荷—位移數(shù)據(jù)按體積不變?cè)砝L制真應(yīng)力—真應(yīng)變曲線(xiàn)。

    2 結(jié)果及分析

    2.1 真應(yīng)力—真應(yīng)變曲線(xiàn)

    圖1 試驗(yàn)合金在不同條件下壓縮的應(yīng)力—應(yīng)變曲線(xiàn)Fig.1 True stress—strain curves of alloy under thermal compress at different temperatures and strain rates: (a) ε&=10?3 s?1;(b) ε&=10?2 s?1; (c) ε&=10?1 s?1; (d) ε&=1 s?1

    GW63合金圓柱形試樣熱壓縮變形的真應(yīng)力—真應(yīng)變(σ—ε)曲線(xiàn)如圖1所示。從圖1中可以看出,變形開(kāi)始階段(約ε≤0.05),應(yīng)力隨應(yīng)變上升很快,但達(dá)到峰值應(yīng)力(pσ)所需應(yīng)變 (pε)相差很大,各試驗(yàn)條件下pε的范圍為0.05≤pε≤0.35。pσ及其對(duì)應(yīng)的應(yīng)變pε均隨溫度降低和應(yīng)變速率升高而增大。當(dāng)應(yīng)力達(dá)到峰值后,σ—ε曲線(xiàn)分為兩類(lèi),一類(lèi)是下降型(D型),即應(yīng)力開(kāi)始隨應(yīng)變的增加而不斷降低;另一類(lèi)是穩(wěn)定型(S型),即應(yīng)力基本不隨應(yīng)變的增加而增加。σ—ε曲線(xiàn)的D型或S型特征明顯受ε&和T的耦合作用,如在應(yīng)變速率較低(10?3~10?2s?1)的情況下,D 型曲線(xiàn)出現(xiàn)在溫度為300~350 ℃;隨應(yīng)變速率升高到1 s?1,溫度為400~450 ℃時(shí),σ—ε曲線(xiàn)也轉(zhuǎn)變?yōu)镈型。

    值得特別注意的是溫度為 300 ℃、應(yīng)變速率為10?3s?1時(shí)達(dá)到峰值應(yīng)力后,仍能維持穩(wěn)態(tài)變形,但升高應(yīng)變速率,變形達(dá)到峰值應(yīng)力后試樣均失穩(wěn)斷裂。因此,此后計(jì)算熱變形參數(shù)時(shí)不考慮溫度為300 ℃時(shí)的數(shù)據(jù)以便得到更加精確的擬合結(jié)果。

    2.2 熱變形本構(gòu)方程

    一般金屬材料熱變形可視為熱力學(xué)過(guò)程,引入變形表觀激活能(Q)的概念,熱變形的變形溫度與應(yīng)變速率條件可采用帶激活能項(xiàng)的 Zener-Hollomon(Z)參數(shù)描述(見(jiàn)式(1))。在 σ -ε&-T 作為獨(dú)立變量的前提下,熱變形本構(gòu)方程有冪函數(shù)(式(2))(簡(jiǎn)稱(chēng) PI)、指數(shù)函數(shù)(式(3))(簡(jiǎn)稱(chēng)EI)和雙曲正弦函數(shù)(式(4))(簡(jiǎn)稱(chēng)SI) 3種半經(jīng)驗(yàn)類(lèi)型的表征方案,它們均可轉(zhuǎn)化為σ=f( Z)型函數(shù)。

    式中:σ是流變應(yīng)力,MPa;R是氣體常數(shù),8.314 J/(mol·K); API、 AEI、A 是應(yīng)變獨(dú)立因子,s?1;nPI和n是應(yīng)力指數(shù);βEI是硬化系數(shù);α是應(yīng)力因子,MPa?1,且 α = βEI/ nPI。

    一般認(rèn)為,PI形式(式(2))適用于低應(yīng)力水平,EI形式(式(3))適用于高應(yīng)力水平[11],SI形式(式(4))適用于較寬的應(yīng)力與應(yīng)變速率范圍[12?13]。然而,應(yīng)力水平的高低并沒(méi)有嚴(yán)格的衡量指標(biāo),只能依具體情況而定,這就意味著對(duì)于同類(lèi)型合金,有必要對(duì)上述3類(lèi)本構(gòu)方程函數(shù)的適應(yīng)性進(jìn)行全面的驗(yàn)證。以下內(nèi)容將采用上述 3種類(lèi)型函數(shù)分別對(duì)實(shí)驗(yàn)合金流變行為進(jìn)行擬合,求解相應(yīng)的材料常數(shù)并進(jìn)行優(yōu)劣性驗(yàn)證。擬合數(shù)據(jù)的范圍如下:變形溫度為350~500 ℃,應(yīng)變速率為10?3~1 s?1。

    從冪函數(shù)(PI)出發(fā),作 l nε&—lnσ 和lnσ—103(RT)?1關(guān)系圖(見(jiàn)圖2),按 QPI= KTKnPI可得應(yīng)變速率為10?3~1 s?1時(shí)各變形溫度下的表觀激活能QPI,如表1中所列。取表1中 PI下數(shù)據(jù)的算術(shù)平均,得=286.44 kJ/mol。

    從指數(shù)函數(shù)(EI)出發(fā),作 l nε &—σ和σ—103(RT)?1關(guān)系圖(見(jiàn)圖3),按= K K同樣可得 TβEI應(yīng)變速率為10?3~1 s?1時(shí)各變形溫度下的表觀激活能QEI,如表1中所示。取表1中EI下數(shù)據(jù)的算術(shù)平均數(shù),得 QEAI=254.65 kJ/mol。

    從雙曲正弦函數(shù)(SI)出發(fā),按文獻(xiàn)[14]中常用的方法,從圖2(a)和 3(a)中可得到各變形條件下的 nPI和βEI值,進(jìn)一步可得到 αSI=βEI/ nPI=0.009 2 MPa?1。圖4所示為 l nε &— l n[sinh(αSIσ)]和ln[sinh(αSIσ)]—103(RT)?1關(guān)系圖,根據(jù)=同樣可得應(yīng)變速 率為 10?3~1 s?1時(shí)各變形溫度下的表觀激活能 QSI,如表1中所列。取表1中SI下數(shù)據(jù)的算術(shù)平均數(shù)值,得=257.79 kJ/mol。

    圖2 PI下的 ln ε&—lnσ和lnσ—103(RT)?1曲線(xiàn)Fig.2 l nε&—lnσ(a) and lnσ—103(RT)?1 (b) curves under PI

    表1 變形速率為10?3~1 s?1時(shí)各變形溫度下各本構(gòu)方程獲得的合金的表觀激活能QTable 1 Q of alloy at different deformation temperatures and strain rates of 10?3?1 s?1 by using different constitutive equations

    圖3 EI下的 ln ε &—σ和σ—103(RT)?1曲線(xiàn)Fig.3 l n ε &—σ(a) and σ—103(RT)?1 (b) curves under EI

    圖4 SI下 的 l n ε &— l n[sinh(αSIσ)]和ln[sinh(αSIσ)]—103(RT)?1曲線(xiàn)Fig.4 l n ε &— ln [sinh(αSIσ)](a) and ln[sinh(αSIσ)]—103(RT)?1(b) curves under SI

    獲得上述PI、EI、SI 3個(gè)函數(shù)所定義的nPI、βEI、αSI以及變形表觀激活能 QPAI、 QEAI、 QSAI后,依照式(1)~(4)可進(jìn)一步求出上述函數(shù)中的 API、 AEI、 ASI等材料常數(shù),如圖5(a)~(c)所示。圖5(a)~(c)的數(shù)據(jù)點(diǎn)及其線(xiàn)性擬合結(jié)果表明,上述3個(gè)函數(shù)均能建立本研究中GW63鎂合金熱變形的σ=f( Z)型本構(gòu)方程。在本研究中,根據(jù)式(5)用這些函數(shù)預(yù)測(cè)應(yīng)力值(σi)與試驗(yàn)值(σPi)間的標(biāo)準(zhǔn)方差(δ)作為判據(jù)來(lái)檢驗(yàn)上述3個(gè)函數(shù)的擬合精度,計(jì)算得δ由小到大的順序?yàn)镋I(9.83)、SI(16.78)、PI(44.87)。可見(jiàn),擬合 GW63合金變形行為時(shí), EI比SI更加接近試驗(yàn)值。EI函數(shù)適合在高應(yīng)力下使用,而 Mg-Gd-Y-X系合金屬于高強(qiáng)耐熱合金[15?17],因此,這也符合本研究在擬合GW63合金熱變形時(shí)指數(shù)函數(shù)(EI)較雙曲正弦函數(shù)(SI)更加精確的實(shí)驗(yàn)事實(shí)。

    圖5 Z與σ的關(guān)系Fig.5 Relationship between Z and peak stress σ of alloys under hot compression deformation: (a) lnZ—lnσ; (b) lnZ—σ; (c) lnZ—ln[sinh(αSIσ)]; (d) lnZ—ln[sinh(αSIOσ)]

    然而,如圖2(a)和3(a)所示,冪函數(shù)(PI)中材料常數(shù)nPI和指數(shù)函數(shù)(EI)中 βEI是分別通過(guò)lnσ—ln ε &和σ— ln ε &線(xiàn)性擬合得到的,因此,無(wú)法對(duì)其進(jìn)行相關(guān)數(shù)學(xué)處理。而 SI中材料常數(shù)α是通過(guò)關(guān)系式α=βEI/ nPI獲得,有一定的數(shù)學(xué)物理基礎(chǔ),可以對(duì)α進(jìn)行優(yōu)化處理,本研究中采用式(5)對(duì) SI函數(shù)中的材料常數(shù)α值進(jìn)行優(yōu)化,并記優(yōu)化后的SI為SIO函數(shù)。

    式中:σi為按式(2)~(4)計(jì)算的各ε&—T條件下的應(yīng)力值;為對(duì)應(yīng)條件下的峰值應(yīng)力。

    在雙曲正弦函數(shù)中,根據(jù)式(5)基于Matlab平臺(tái)進(jìn)行編程,通過(guò)程序搜索得到最小的δ值為9.17,此時(shí)得到SIO函數(shù)下的αSIO=0.033 MPa?1,把此值代入SIO函數(shù),則類(lèi)似可按 QSIO= KTKnSIO得到應(yīng)變速率為10?3~1 s?1時(shí)各變形溫度下的表觀激活能 QSIO如表1所列,取表中SIO下數(shù)據(jù)的算術(shù)平均值,得252.94 kJ/mol。可由圖5(d)求出為 exp33.89,由圖可直觀看出SIO下的線(xiàn)性擬合誤差小于圖5(a)~(c)的。

    綜上可知,在使用PI、EI和SI函數(shù)擬合實(shí)驗(yàn)合金熱壓縮變形行為時(shí),由δ檢驗(yàn)結(jié)果知 EI擬合優(yōu)于SI。PI、EI與SI對(duì)本研究GW63合金的適應(yīng)性與文獻(xiàn)報(bào)道的 AZ系列合金的不同,出現(xiàn)這種情況可能與Mg-Gd-Y-X系列合金的高溫下強(qiáng)度高有關(guān)。在相同的高溫變形條件下,該系合金的強(qiáng)度遠(yuǎn)高于AZ系列合金的,所以,適用于較高應(yīng)力水平的指數(shù)函數(shù)(EI)本構(gòu)方程更加接近實(shí)驗(yàn)結(jié)果。若對(duì)SI函數(shù)中材料常數(shù)α進(jìn)行優(yōu)化處理,則最終δ由小到大依次為 SIO、EI、SI、PI,可見(jiàn),優(yōu)化處理對(duì)提高雙曲正弦函數(shù)的精度具有非常重要的意義。本文作者提出的方法(式(5))提高預(yù)測(cè)精度的效果十分明顯,如果不進(jìn)行材料常數(shù)優(yōu)化,在本研究的試驗(yàn)條件下,雙曲正弦函數(shù)的精度低于指數(shù)函數(shù)的精度。

    3 討論

    圖6 采用各本構(gòu)方程得到的合金的表觀激活能變化趨勢(shì)Fig.6 Variation trend of apparent activation energies of alloy attained at 350?500 ℃ by using different constitutive equations

    在實(shí)驗(yàn)數(shù)據(jù)擬合范圍設(shè)定為變形溫度 350~500℃、應(yīng)變速率 10?3~1 s?1時(shí),各本構(gòu)方程下得到的表觀激活能變化趨勢(shì)如圖6所示(見(jiàn)表1)。從圖6中可以看出,對(duì)GW63合金熱壓縮變形擬合較精確的EI和SIO函數(shù)下得到的表觀激活能隨變形溫度變化規(guī)律相似,均為先下降后升高,而在擬合性較差的 SI和 PI函數(shù)下表觀激活能數(shù)值則一直呈下降趨勢(shì)。由于GW63合金中稀土元素Gd與Y的原子半徑較大,熱壓縮后組織中仍然存在富Gd和Y的稀土粒子,如圖7及表2所示,使得合金熱壓縮流變應(yīng)力處于高應(yīng)力水平且表觀激活能整體數(shù)值偏大。同時(shí)也看到,EI和SIO函數(shù)在變形溫度為500 ℃時(shí),表觀激活能值升高,這可能是500 ℃變形時(shí)多系滑移啟動(dòng)導(dǎo)致的。

    擬合AZ和ZK系列合金熱壓縮流變行為時(shí),文獻(xiàn)均指出在較低溫度下適合使用指數(shù)函數(shù),溫度升高適合使用冪函數(shù)[5?7]。簽于此,本研究中重新將實(shí)驗(yàn)數(shù)據(jù)擬合范圍設(shè)定為:變形溫度400~500 ℃,應(yīng)變速率10?3~1 s?1,分別使用 PI、EI、SI和 SIO 本構(gòu)方程擬合實(shí)驗(yàn)合金在較高溫度下的熱壓縮流變行為并分析其優(yōu)劣性。采用式(1)~(5)相同的方法計(jì)算得到此時(shí)各擬合函數(shù)下合金的表觀激活能(見(jiàn)表3),δ由小到大的順序?yàn)?SIO(3.58)、EI(5.38)、SI(6.78)、PI(24.05)。由此可見(jiàn),無(wú)論擬合時(shí)變形溫度范圍選為350~500 ℃,還是400~500 ℃,各本構(gòu)方程擬合優(yōu)劣性結(jié)果一致,均是SIO與EI函數(shù)擬合更加接近實(shí)驗(yàn)的σ—ε&—T結(jié)果。

    圖7 變形溫度為350~500 ℃時(shí)合金的SEM像Fig.7 SEM images of alloy compressed at 350?500 ℃: (a) 350 ℃, 10?3 s?1; (b) 400 ℃, 10?2 s?1; (c) 450 ℃, 1 s?1; (d) 500 ℃, 1 s?1

    表2 圖7(d)中A點(diǎn)第二相能譜分析結(jié)果Table 2 EDS results of point A in Fig.7(d)

    表3 變形溫度為 350~500 ℃和 400~500 ℃時(shí)用 PI、EI、SI和SIO得到的合金的表觀激活能Table 3 Apparent activation energy attained by PI, EI, SI and SIO at temperatures of 350?500 ℃ and 400?500 ℃

    本研究中,指數(shù)函數(shù)(EI)形式適用于高應(yīng)力水平的物理意義符合本研究的實(shí)驗(yàn)擬合事實(shí),即EI擬合精度優(yōu)于雙曲正弦函數(shù)(SI)的。雖然對(duì)SI中材料常數(shù)α進(jìn)行優(yōu)化處理后得到了擬合精度略?xún)?yōu)于指數(shù)函數(shù)的SIO函數(shù),然而優(yōu)化式(5)中α值范圍的選定是多次嘗試搜索得到的最佳結(jié)果,計(jì)算過(guò)程相對(duì)于采用指數(shù)函數(shù)擬合時(shí)更加繁瑣。在Matlab編好的程序中,假如輸入的α值搜索范圍內(nèi)沒(méi)有包含能夠得到最小δ時(shí)的αSIO(0.033 MPa?1)值,那么優(yōu)化后的雙曲正弦函數(shù)(SIO)也將得不到比 EI函數(shù)更好的擬合結(jié)果。因此,擬合GW63合金熱壓縮流變行為時(shí)采用指數(shù)函數(shù)(EI)最簡(jiǎn)便且能保證足夠的擬合精度。

    4 結(jié)論

    1) Mg-6Gd-3Y-0.5Zr合金中Gd和Y的原子半徑較大,合金中高溫耐熱富稀土相使得合金高溫強(qiáng)度整體升高,合金熱壓縮流變應(yīng)力較大,得到的表觀激活能較高。

    2) 擬合 Mg-6Gd-3Y-0.5Zr合金熱壓縮變形行為時(shí),采用指數(shù)函數(shù)(EI)本構(gòu)方程擬合的精度高于冪函數(shù)(PI)和雙曲正弦函數(shù)(SI)的,能得到最簡(jiǎn)便并與指數(shù)函數(shù)精度足夠的擬合結(jié)果。

    3) 采用對(duì)雙曲正弦函數(shù)(SI)中材料常數(shù)α進(jìn)行優(yōu)化處理后的SIO函數(shù)擬合實(shí)驗(yàn)數(shù)據(jù)雖然能夠得到比指數(shù)函數(shù)(EI)擬合精度略?xún)?yōu)的實(shí)驗(yàn)結(jié)果,但擬合時(shí)α取值范圍需多次搜索嘗試并設(shè)定為[SIα/10,10SIα]時(shí)才能得到最佳的SIOα值。

    4) 變形溫度為 350~500 ℃、變形速率為 0.001~1 s?1時(shí),求解得到SIOα=0.033 MPa?1,進(jìn)而得到采用SIO擬合時(shí)合金熱壓縮變形表觀激活能和常數(shù)A分別為252.94 kJ/mol和exp33.89。

    REFERENCES

    [1] 張新明, 彭卓凱, 陳健美, 鄧運(yùn)來(lái). 耐熱鎂合金及其研究進(jìn)展[J]. 中國(guó)有色金屬學(xué)報(bào), 2004, 14(9): 1443?1450.ZHANG Xin-ming, PENG Zhuo-kai, CHEN Jian-mei, DENG Yun-lai. Heat-resistant magnesium alloys and their development[J]. The Chinese Journal of Nonferrous Metals,2004, 14(9): 1443?1450.

    [2] PRASAD Y V R K, RAO K P. Processing maps for hot deformation of rolled AZ31 magnesium alloy plate: Anisotropy of hot workability[J]. Materials Science and Engineering A, 2008,487(1/2): 316?327.

    [3] YANG Z, GUO Y C, LI J P, HE F, XIA F, LIANG M X. Plastic deformation and dynamic recrystallization behaviors of Mg-5Gd-4Y-0.5Zn-0.5Zr alloy[J]. Materials Science and Engineering A, 2008, 485(1/2): 487?491.

    [4] YANG Y Q, LI B C, ZHANG Z M. Flow stress of wrought magnesium alloys during hot compression deformation at medium and high temperatures[J]. Materials Science and Engineering A, 2009, 499: 238?241.

    [5] GALIYEV A, KAIBYSHEV R, GOTTSTEIN G. Correlation of plastic deformation and dynamic recrystallization in magnesium alloy ZK60[J]. Acta Materialia, 2001, 49(7): 1199?1207.

    [6] GUO Q, YAN H G, ZHANG H, CHEN Z H, WANG Z F.Behaviour of AZ31 magnesium alloy during compression at elevated temperatures[J]. Materials Science and Technology,2005, 21(11): 1349?1354.

    [7] GUO Q, YAN H G, CHEN Z H, ZHANG H. Elevated temperature compression behaviour of Mg-Al-Zn alloys[J].Materials Science and Technology, 2006, 22(6): 725?729.

    [8] FATEMI-VARZANEH S M, ZAREI-HANZAKI A,HAGHSHENAS M. A study on the effect of thermo-mechanical parameters on the deformation behavior of Mg-3Al-1Zn[J].Materials Science and Engineering A, 2008, 497: 438?444.

    [9] 肖 陽(yáng), 張新明, 陳健美, 蔣 浩, 鄧楨楨. 高強(qiáng)耐熱Mg-9Gd-4Y-0.6Zr合金的性能[J]. 中南大學(xué)學(xué)報(bào): 自然科學(xué)版,2006, 16(5): 850?855.XIAO Yang, ZHANG Xin-ming, CHEN Jian-mei, JIANG Hao,DENG Zhen-zhen. Performance of Mg-9Gd-4Y-0.6Zr alloy with high strength and heat resistance[J]. Journal of Central South University: Science and Technology, 2006, 16(5): 850?855.

    [10] LI L, ZHANG X M. Hot compression deformation behavior and processing parameters of a cast Mg-Gd-Y-Zr alloy[J]. Materials Science and Engineering A,

    [11] MCQUEEN H J, RYAN N D. Constitutive analysis in hot working[J]. Materials Science and Engineering A, 2002,322(1/2): 43?63.

    [12] LI L, ZHOU J, DUSZCZYK J. Determination of a constitutive relationship for AZ31B magnesium alloy and validation through comparison between simulated and real extrusion[J]. Journal of Materials Processing Technology, 2006, 172(3): 372?380.

    [13] SLOOFF F A, ZHOU J, DUSZCZYK J, KATGERMAN L.Constitutive analysis of wrought magnesium alloy Mg-Al4-Zn1[J]. Scripta Materialia, 2007, 57(8): 759?762.

    [14] POIRIER J P. 晶體的高溫塑性變形[M]. 關(guān)德林, 譯. 大連:大連理工大學(xué)出版社, 1988: 42.POIRIER J P. The plastic deformation of crystals at high temperature[M]. GUAN De-lin, transl. Dalian: Dalian University of Science and Technology Press, 1988: 42.

    [15] HE S M, ZENG X Q, PENG L M, GAO X, NIE J F, DING W J.Microstructure and strengthening mechanism of high strength Mg-10Gd-2Y-0.5Zr alloy[J]. Journal of Alloys and Compounds,2007, 427(1/2): 316?323.

    [16] YANG Z, LI J P, GUO Y C, LIU T, XIA F, ZENG Z W, LIANG M X. Precipitation process and effect on mechanical properties of Mg-9Gd-3Y-0.6Zn-0.5Zr alloy[J]. Materials Science and Engineering A, 2007, 454/455: 274?280.

    [17] 張新明, 陳健美, 鄧運(yùn)來(lái), 肖 陽(yáng), 蔣 浩. Mg-Gd-Y-Zr耐熱鎂合金的壓縮變形行為[J]. 中國(guó)有色金屬學(xué)報(bào), 2005, 15(12):1925?1932.ZHANG Xin-ming, CHEN Jian-mei, DENG Yun-lai, XIAO Yang, JIANG Hao. Deformation behavior of Mg-Gd-Y-Zr heat-resistant magnesium alloy during hot-compression[J]. The Chinese Journal of Nonferrous Metals, 2005, 15(12):1925?1932.

    Constitutive equation during hot compression deformation of Mg-Gd-Y-Zr alloy

    ZHANG Xin-ming, WU Yi-ping, DENG Yun-lai, TANG Chang-ping
    (School of Materials Science and Engineering, Central South University, Changsha 410083, China)

    Three kinds of constitutive equations, power law (PI), expression law (EI) and sine hyperbolic law (SI) were used to analyze the flow behavior of the Mg-6Gd-3Y-0.5Zr alloy after compression at temperatures of 300?500 ℃ and strain rates of 10?3?1 s?1by using Geeble?1500 test machine. The microstructures were studied using scanning electron microscope (SEM). The results show that EI fits the experimental results better than PI and SI. And even, the sine hyperbolic law with the optimization of material parameter α (SIO) exhibits the close precision to EI. The precipitated heat-resisting phase after compression is the reason for EI better fitting than PI and SI.

    Mg-Gd-Y-Zr alloy; hot compression deformation; constitutive equation; material parameter optimization

    TG146

    A

    1004-0609(2011)12-2987-08

    國(guó)防預(yù)研基金資助項(xiàng)目(51312010503)

    2010-11-12;

    2011-04-08

    張新明,教授,博士;電話(huà):0731-88830265;E-mail: wuyipingjia@126.com

    (編輯 李艷紅)

    猜你喜歡
    激活能雙曲指數(shù)函數(shù)
    基于加速退化試驗(yàn)的脈寬調(diào)制器貯存壽命預(yù)測(cè)研究
    中國(guó)科學(xué)技術(shù)館之“雙曲隧道”
    軍事文摘(2021年22期)2022-01-18 06:22:48
    冪函數(shù)、指數(shù)函數(shù)、對(duì)數(shù)函數(shù)(2)
    冪函數(shù)、指數(shù)函數(shù)、對(duì)數(shù)函數(shù)(1)
    冪函數(shù)、指數(shù)函數(shù)、對(duì)數(shù)函數(shù)(1)
    冪函數(shù)、指數(shù)函數(shù)、對(duì)數(shù)函數(shù)(2)
    雙曲型交換四元數(shù)的極表示
    一種微合金化鋼的動(dòng)態(tài)再結(jié)晶行為研究
    一階雙曲型偏微分方程的模糊邊界控制
    Zr57.5Cu27.3Al8.5Ni6.7非晶合金的非等溫和等溫晶化動(dòng)力學(xué)研究
    黄网站色视频无遮挡免费观看| 亚洲欧洲精品一区二区精品久久久| 欧美黑人巨大hd| 欧美激情高清一区二区三区| 日本a在线网址| 视频在线观看一区二区三区| 国产片内射在线| 美国免费a级毛片| 国产爱豆传媒在线观看 | av在线天堂中文字幕| av福利片在线| 久久午夜综合久久蜜桃| 久热爱精品视频在线9| 午夜福利18| 欧美激情极品国产一区二区三区| 男人的好看免费观看在线视频 | 老熟妇乱子伦视频在线观看| 波多野结衣av一区二区av| 久久精品影院6| 国产亚洲av嫩草精品影院| 欧美日韩亚洲国产一区二区在线观看| 一本综合久久免费| 国产在线观看jvid| 亚洲九九香蕉| 欧美一级毛片孕妇| 国产伦一二天堂av在线观看| 亚洲国产精品999在线| 国产精品爽爽va在线观看网站 | 亚洲熟女毛片儿| 久久亚洲精品不卡| 亚洲一码二码三码区别大吗| 国产成人啪精品午夜网站| 亚洲自偷自拍图片 自拍| www.自偷自拍.com| 丝袜美腿诱惑在线| 久久久国产成人精品二区| 精品高清国产在线一区| 2021天堂中文幕一二区在线观 | 曰老女人黄片| 欧美在线黄色| 欧美黑人巨大hd| 欧美中文综合在线视频| 国产精品一区二区精品视频观看| 欧美黑人精品巨大| 国产成人精品无人区| 在线av久久热| 不卡av一区二区三区| 精品久久久久久,| 久久中文字幕人妻熟女| 精品午夜福利视频在线观看一区| 成人欧美大片| 国产精品一区二区精品视频观看| av欧美777| 国产91精品成人一区二区三区| 99国产精品一区二区蜜桃av| 丰满人妻熟妇乱又伦精品不卡| 亚洲av五月六月丁香网| 亚洲成人国产一区在线观看| 最近在线观看免费完整版| 国产在线精品亚洲第一网站| 妹子高潮喷水视频| 天堂√8在线中文| 最近在线观看免费完整版| 久久青草综合色| 身体一侧抽搐| 黄色视频不卡| 很黄的视频免费| 9191精品国产免费久久| 亚洲国产精品999在线| 一夜夜www| 特大巨黑吊av在线直播 | 一级片免费观看大全| 欧美中文综合在线视频| 国产精品亚洲一级av第二区| 变态另类丝袜制服| 99热6这里只有精品| 午夜久久久在线观看| 女人爽到高潮嗷嗷叫在线视频| 少妇粗大呻吟视频| 黑丝袜美女国产一区| 精品一区二区三区四区五区乱码| 少妇熟女aⅴ在线视频| 午夜福利欧美成人| 又大又爽又粗| 欧美大码av| 日韩国内少妇激情av| 国产熟女xx| 91av网站免费观看| 欧美黑人欧美精品刺激| 老汉色av国产亚洲站长工具| 男女床上黄色一级片免费看| 在线天堂中文资源库| 久久国产精品男人的天堂亚洲| 久久国产精品男人的天堂亚洲| 最好的美女福利视频网| 成人国产一区最新在线观看| 欧美亚洲日本最大视频资源| 老汉色∧v一级毛片| 999精品在线视频| 不卡一级毛片| 一级a爱片免费观看的视频| 国产激情偷乱视频一区二区| 成人av一区二区三区在线看| 免费看美女性在线毛片视频| 成人国语在线视频| 757午夜福利合集在线观看| 久久狼人影院| 国产色视频综合| 黄片播放在线免费| 99国产极品粉嫩在线观看| 99国产综合亚洲精品| 中出人妻视频一区二区| 黄色成人免费大全| 哪里可以看免费的av片| 成年女人毛片免费观看观看9| 女人爽到高潮嗷嗷叫在线视频| 日韩精品青青久久久久久| 免费高清视频大片| 久久天堂一区二区三区四区| 欧美色欧美亚洲另类二区| 午夜福利在线观看吧| 久久精品91无色码中文字幕| 三级毛片av免费| x7x7x7水蜜桃| 在线av久久热| 国产精品 欧美亚洲| 91麻豆av在线| 桃色一区二区三区在线观看| 女人被狂操c到高潮| 欧美一级a爱片免费观看看 | 美女国产高潮福利片在线看| 久久精品国产亚洲av高清一级| 黄频高清免费视频| 高清毛片免费观看视频网站| av片东京热男人的天堂| 夜夜躁狠狠躁天天躁| 十八禁人妻一区二区| 亚洲国产欧美一区二区综合| 欧美黑人巨大hd| 欧美日韩一级在线毛片| 性欧美人与动物交配| 成人一区二区视频在线观看| 亚洲 欧美 日韩 在线 免费| 一区二区日韩欧美中文字幕| 亚洲性夜色夜夜综合| 免费看十八禁软件| 精品久久久久久成人av| 黄色毛片三级朝国网站| 国产在线精品亚洲第一网站| 50天的宝宝边吃奶边哭怎么回事| 国产精品永久免费网站| 精品一区二区三区av网在线观看| 国产aⅴ精品一区二区三区波| 最新美女视频免费是黄的| 亚洲在线自拍视频| 啦啦啦韩国在线观看视频| 欧美日本视频| 黄色成人免费大全| 久久人人精品亚洲av| 宅男免费午夜| 国产亚洲精品久久久久5区| 香蕉国产在线看| 亚洲成人久久爱视频| 男人舔女人下体高潮全视频| 两个人免费观看高清视频| 国产又爽黄色视频| 性色av乱码一区二区三区2| 99热6这里只有精品| 国产精品亚洲美女久久久| av天堂在线播放| 亚洲精品国产精品久久久不卡| cao死你这个sao货| 亚洲成人久久爱视频| 老司机福利观看| 国产av又大| 最新美女视频免费是黄的| 亚洲 欧美 日韩 在线 免费| 1024香蕉在线观看| 免费看十八禁软件| 国产人伦9x9x在线观看| 麻豆成人av在线观看| 一区二区日韩欧美中文字幕| 国产午夜精品久久久久久| 欧美精品啪啪一区二区三区| 精品少妇一区二区三区视频日本电影| 国产视频一区二区在线看| 久久久久久九九精品二区国产 | 日本成人三级电影网站| 日韩欧美 国产精品| 人人妻人人看人人澡| 曰老女人黄片| 亚洲国产欧美日韩在线播放| 波多野结衣高清无吗| 久久国产乱子伦精品免费另类| 中文字幕高清在线视频| 成人国产综合亚洲| 亚洲精品久久成人aⅴ小说| 性色av乱码一区二区三区2| 黄色片一级片一级黄色片| 亚洲国产欧美网| 18禁观看日本| 国产精品av久久久久免费| 久99久视频精品免费| 亚洲av片天天在线观看| 亚洲三区欧美一区| 亚洲精品在线观看二区| 国产一区二区激情短视频| 桃红色精品国产亚洲av| 精品久久久久久久久久久久久 | 国产真实乱freesex| 国产成人啪精品午夜网站| 亚洲中文日韩欧美视频| 国产99白浆流出| 一级黄色大片毛片| 国产精品亚洲av一区麻豆| 亚洲精品久久国产高清桃花| 美女扒开内裤让男人捅视频| 亚洲av五月六月丁香网| av在线播放免费不卡| 色精品久久人妻99蜜桃| av超薄肉色丝袜交足视频| 国产精品二区激情视频| 久久人妻福利社区极品人妻图片| 欧美乱妇无乱码| 成人三级做爰电影| 男女那种视频在线观看| 国产精品乱码一区二三区的特点| 久久久久久免费高清国产稀缺| 欧美在线黄色| 啦啦啦 在线观看视频| 香蕉国产在线看| 日韩欧美国产在线观看| 不卡av一区二区三区| 国产真人三级小视频在线观看| av在线天堂中文字幕| 久久草成人影院| 人妻丰满熟妇av一区二区三区| 高潮久久久久久久久久久不卡| 在线观看免费日韩欧美大片| 色尼玛亚洲综合影院| 国产亚洲av高清不卡| 免费在线观看完整版高清| 又大又爽又粗| 伦理电影免费视频| 成人亚洲精品一区在线观看| www.www免费av| 精品第一国产精品| 丝袜在线中文字幕| 变态另类丝袜制服| 男人的好看免费观看在线视频 | 亚洲国产精品999在线| 亚洲国产看品久久| 99国产综合亚洲精品| 一二三四社区在线视频社区8| 两个人视频免费观看高清| 黄色丝袜av网址大全| 国产精品98久久久久久宅男小说| 天天一区二区日本电影三级| 精品国内亚洲2022精品成人| av免费在线观看网站| 国产精品二区激情视频| 桃色一区二区三区在线观看| 欧美不卡视频在线免费观看 | 国产在线观看jvid| 免费看十八禁软件| 黑人欧美特级aaaaaa片| 国产精品影院久久| 免费人成视频x8x8入口观看| 国产区一区二久久| 国产成+人综合+亚洲专区| 欧美一区二区精品小视频在线| 午夜免费观看网址| 午夜福利在线观看吧| 国产真实乱freesex| 久久国产精品男人的天堂亚洲| 国产精品av久久久久免费| 国产亚洲欧美精品永久| 丝袜在线中文字幕| 丝袜美腿诱惑在线| 亚洲avbb在线观看| 黑人巨大精品欧美一区二区mp4| 1024香蕉在线观看| 午夜福利在线观看吧| 成熟少妇高潮喷水视频| 成人一区二区视频在线观看| 日本黄色视频三级网站网址| 国产高清激情床上av| 18禁美女被吸乳视频| 99久久综合精品五月天人人| 真人做人爱边吃奶动态| 欧美亚洲日本最大视频资源| 国产欧美日韩精品亚洲av| 最近在线观看免费完整版| 精品久久久久久久久久久久久 | 18禁观看日本| 国产一区二区激情短视频| 在线看三级毛片| 在线十欧美十亚洲十日本专区| 亚洲精品久久成人aⅴ小说| а√天堂www在线а√下载| 亚洲一码二码三码区别大吗| 一级a爱片免费观看的视频| 国产精品自产拍在线观看55亚洲| 别揉我奶头~嗯~啊~动态视频| 亚洲中文日韩欧美视频| 此物有八面人人有两片| 欧美成人一区二区免费高清观看 | 午夜视频精品福利| 欧美黑人欧美精品刺激| 后天国语完整版免费观看| 99精品在免费线老司机午夜| 色综合婷婷激情| 夜夜夜夜夜久久久久| 免费观看精品视频网站| 国产精品综合久久久久久久免费| 久久精品aⅴ一区二区三区四区| 色播亚洲综合网| 九色国产91popny在线| av中文乱码字幕在线| 国产亚洲精品久久久久5区| 色尼玛亚洲综合影院| 天天躁夜夜躁狠狠躁躁| 天天添夜夜摸| 69av精品久久久久久| 最新在线观看一区二区三区| 亚洲精品在线美女| 国产欧美日韩精品亚洲av| 成人手机av| 男女那种视频在线观看| 18禁黄网站禁片免费观看直播| 亚洲色图 男人天堂 中文字幕| 精品高清国产在线一区| 一区福利在线观看| 很黄的视频免费| 亚洲成人免费电影在线观看| 丝袜美腿诱惑在线| 国产久久久一区二区三区| 99久久久亚洲精品蜜臀av| 久久中文字幕一级| 婷婷丁香在线五月| 亚洲三区欧美一区| 亚洲熟女毛片儿| 亚洲av中文字字幕乱码综合 | 欧美黑人欧美精品刺激| 一本一本综合久久| 国产1区2区3区精品| 国产精品久久久久久亚洲av鲁大| 日本a在线网址| 国产在线精品亚洲第一网站| 精品人妻1区二区| 欧美激情高清一区二区三区| 夜夜夜夜夜久久久久| 国产成+人综合+亚洲专区| 男人舔女人下体高潮全视频| 听说在线观看完整版免费高清| 夜夜看夜夜爽夜夜摸| 长腿黑丝高跟| 特大巨黑吊av在线直播 | 精品不卡国产一区二区三区| 亚洲av美国av| 国产精品久久电影中文字幕| 亚洲精品一区av在线观看| 久久九九热精品免费| 99国产精品一区二区三区| 亚洲精品粉嫩美女一区| svipshipincom国产片| 亚洲国产毛片av蜜桃av| 一进一出抽搐动态| 侵犯人妻中文字幕一二三四区| 日韩视频一区二区在线观看| 中文在线观看免费www的网站 | 久久精品国产亚洲av高清一级| 亚洲成国产人片在线观看| 日本成人三级电影网站| 亚洲熟妇熟女久久| 老司机午夜福利在线观看视频| 亚洲人成网站在线播放欧美日韩| 久久久久久大精品| x7x7x7水蜜桃| av电影中文网址| 亚洲国产看品久久| 欧美人与性动交α欧美精品济南到| 最近最新中文字幕大全免费视频| 俺也久久电影网| 日韩大尺度精品在线看网址| 日韩精品青青久久久久久| cao死你这个sao货| 亚洲国产欧美网| 免费在线观看完整版高清| 欧美乱妇无乱码| 女性被躁到高潮视频| 亚洲真实伦在线观看| 日韩av在线大香蕉| 黄片播放在线免费| 国产免费av片在线观看野外av| 最新在线观看一区二区三区| 欧美日韩瑟瑟在线播放| 免费高清视频大片| 国产激情久久老熟女| 免费在线观看影片大全网站| a级毛片a级免费在线| 日韩欧美在线二视频| 国产99白浆流出| 日本 欧美在线| 国产激情偷乱视频一区二区| 国产av一区二区精品久久| 人人妻人人澡人人看| 级片在线观看| 真人一进一出gif抽搐免费| 成人一区二区视频在线观看| 波多野结衣高清无吗| 亚洲精品久久国产高清桃花| 99国产精品99久久久久| 国产又色又爽无遮挡免费看| 国产色视频综合| 美国免费a级毛片| 一二三四社区在线视频社区8| 亚洲人成伊人成综合网2020| 一级毛片高清免费大全| 久久 成人 亚洲| 午夜福利18| 欧美日韩黄片免| 一个人免费在线观看的高清视频| 两个人看的免费小视频| 午夜福利一区二区在线看| 中文字幕精品亚洲无线码一区 | 亚洲中文字幕一区二区三区有码在线看 | 午夜福利免费观看在线| 日本在线视频免费播放| 午夜福利欧美成人| 免费在线观看成人毛片| 精品一区二区三区视频在线观看免费| 亚洲九九香蕉| 黄片播放在线免费| 黄片大片在线免费观看| 两个人看的免费小视频| 国产真人三级小视频在线观看| 男人操女人黄网站| 国产精品乱码一区二三区的特点| 精品久久久久久久末码| 黄频高清免费视频| 久久 成人 亚洲| 99热只有精品国产| 成人亚洲精品av一区二区| 99国产极品粉嫩在线观看| 色播亚洲综合网| 91九色精品人成在线观看| 十八禁网站免费在线| 久久久久久久久免费视频了| 一进一出好大好爽视频| www.精华液| 中文字幕人成人乱码亚洲影| 91成年电影在线观看| 亚洲avbb在线观看| 波多野结衣av一区二区av| 亚洲一区高清亚洲精品| 亚洲美女黄片视频| 最新在线观看一区二区三区| 男女午夜视频在线观看| 麻豆成人午夜福利视频| 国产一区二区三区视频了| 中国美女看黄片| av免费在线观看网站| 亚洲七黄色美女视频| 91在线观看av| 亚洲av美国av| 久久久久国产一级毛片高清牌| 亚洲精品中文字幕在线视频| 老鸭窝网址在线观看| 久久精品国产综合久久久| www.999成人在线观看| 国产av一区在线观看免费| 91成人精品电影| 俺也久久电影网| 国产乱人伦免费视频| 亚洲av熟女| 99国产精品一区二区三区| 国产精品日韩av在线免费观看| 法律面前人人平等表现在哪些方面| 人人妻人人澡欧美一区二区| 亚洲精华国产精华精| 成人永久免费在线观看视频| 亚洲精品粉嫩美女一区| 欧美另类亚洲清纯唯美| 国产精品久久久av美女十八| 99久久精品国产亚洲精品| 婷婷精品国产亚洲av在线| 亚洲人成电影免费在线| 国产一级毛片七仙女欲春2 | 亚洲色图 男人天堂 中文字幕| 亚洲人成电影免费在线| 亚洲电影在线观看av| 国产精品久久久久久人妻精品电影| 亚洲专区字幕在线| 欧美日韩乱码在线| 真人做人爱边吃奶动态| 女性生殖器流出的白浆| 久久 成人 亚洲| 男女床上黄色一级片免费看| 亚洲国产精品久久男人天堂| 午夜福利18| 美女 人体艺术 gogo| 一区二区三区激情视频| 日本免费a在线| 脱女人内裤的视频| 美女高潮喷水抽搐中文字幕| 男人操女人黄网站| 黄片大片在线免费观看| 日日夜夜操网爽| 宅男免费午夜| 在线播放国产精品三级| 日韩精品免费视频一区二区三区| 午夜福利视频1000在线观看| 成人国产综合亚洲| 国产黄色小视频在线观看| 精品久久久久久久人妻蜜臀av| 久久久久久久久久黄片| 正在播放国产对白刺激| 黄色成人免费大全| 侵犯人妻中文字幕一二三四区| av在线天堂中文字幕| 亚洲av片天天在线观看| 亚洲精品美女久久久久99蜜臀| 日本一本二区三区精品| 国产精品久久久久久精品电影 | 国产精品98久久久久久宅男小说| a级毛片在线看网站| 美女国产高潮福利片在线看| 久久人妻福利社区极品人妻图片| 国产激情偷乱视频一区二区| 亚洲熟妇中文字幕五十中出| 91老司机精品| 色综合亚洲欧美另类图片| 欧美国产日韩亚洲一区| 亚洲片人在线观看| 亚洲第一欧美日韩一区二区三区| www.自偷自拍.com| 久久香蕉激情| 免费无遮挡裸体视频| 免费无遮挡裸体视频| 国产亚洲精品久久久久5区| 美国免费a级毛片| 欧美性猛交黑人性爽| 国产亚洲精品av在线| av在线天堂中文字幕| 尾随美女入室| 91在线精品国自产拍蜜月| 精品福利观看| 国产精品久久久久久av不卡| 亚洲色图av天堂| 毛片一级片免费看久久久久| av天堂中文字幕网| 卡戴珊不雅视频在线播放| 可以在线观看的亚洲视频| 观看免费一级毛片| 91麻豆精品激情在线观看国产| 久久久久精品国产欧美久久久| av国产免费在线观看| 国内久久婷婷六月综合欲色啪| 又粗又爽又猛毛片免费看| 桃色一区二区三区在线观看| 91av网一区二区| 国产91av在线免费观看| 久久精品久久久久久噜噜老黄 | 免费在线观看影片大全网站| 午夜精品在线福利| 在线观看66精品国产| 97热精品久久久久久| 亚洲成人av在线免费| 午夜老司机福利剧场| 全区人妻精品视频| 一边摸一边抽搐一进一小说| 三级男女做爰猛烈吃奶摸视频| 亚洲精品影视一区二区三区av| 一级黄色大片毛片| 春色校园在线视频观看| 国产午夜福利久久久久久| 高清毛片免费观看视频网站| 不卡一级毛片| 一进一出抽搐动态| 国内精品美女久久久久久| 成年女人看的毛片在线观看| 亚洲高清免费不卡视频| 熟女电影av网| 久久久国产成人免费| 97碰自拍视频| 91精品国产九色| 日韩制服骚丝袜av| 国产高清三级在线| 久久久久久久午夜电影| 18+在线观看网站| 日本爱情动作片www.在线观看 | 国产私拍福利视频在线观看| 淫妇啪啪啪对白视频| 亚洲熟妇熟女久久| 亚洲欧美成人精品一区二区| 最近中文字幕高清免费大全6| 三级经典国产精品| videossex国产| 国产一级毛片七仙女欲春2| 国产精品美女特级片免费视频播放器| 免费看a级黄色片| 欧美另类亚洲清纯唯美| 在线免费观看的www视频| 免费人成在线观看视频色| 欧美性猛交╳xxx乱大交人| 国产精品国产高清国产av| 亚洲欧美清纯卡通| 大型黄色视频在线免费观看| 亚洲欧美精品综合久久99| www.色视频.com| 国产精品电影一区二区三区| 在线观看美女被高潮喷水网站| 欧美区成人在线视频| 精品人妻熟女av久视频|