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

    Ti6554鈦合金脈沖電流輔助壓縮本構(gòu)模型的建立及應(yīng)用

    2022-05-16 03:09:08劉斌武川周宇杰
    精密成形工程 2022年5期
    關(guān)鍵詞:脈沖電流本構(gòu)鈦合金

    劉斌,武川,周宇杰

    Ti6554鈦合金脈沖電流輔助壓縮本構(gòu)模型的建立及應(yīng)用

    劉斌,武川,周宇杰

    (天津職業(yè)技術(shù)師范大學(xué) 汽車(chē)模具智能制造國(guó)家地方聯(lián)合工程實(shí)驗(yàn)室,天津 300222)

    基于電輔助壓縮試驗(yàn),建立Ti6554鈦合金的電輔助壓縮本構(gòu)模型,并對(duì)本構(gòu)模型進(jìn)行驗(yàn)證。針對(duì)Ti6554鈦合金材料,在一定脈沖電流參數(shù)(電流密度為40、32、24 A/mm2,頻率為1、5、10 Hz,占空比為0.7、0.5、0.3)以及應(yīng)變速率為0.002、0.005、0.01 s?1、下壓量為50%的條件下,進(jìn)行電輔助壓縮試驗(yàn),并基于唯象和物理機(jī)制,建立電輔助壓縮的本構(gòu)模型,采用回歸擬合的方法確定本構(gòu)模型參數(shù),得到具體方程,并基于電輔助壓縮試驗(yàn)的應(yīng)力–應(yīng)變曲線,分析脈沖電流參數(shù)對(duì)應(yīng)力–應(yīng)變的影響規(guī)律。將電輔助壓縮試驗(yàn)測(cè)得的試驗(yàn)值與本構(gòu)模型得到的預(yù)測(cè)值進(jìn)行相對(duì)誤差分析,得到最大相對(duì)誤差小于7%。采用回歸擬合的方法建立的本構(gòu)模型相對(duì)誤差較小,能夠準(zhǔn)確預(yù)測(cè)材料在不同脈沖電流參數(shù)下的流動(dòng)應(yīng)力。

    電致塑性;Ti6554;本構(gòu)方程;脈沖電流

    Ti6554(Ti–6Mo–5Cr–5V–4Al)鈦合金是一種亞穩(wěn)定高強(qiáng)韌β型鈦合金,因其具有密度小、比強(qiáng)度大、無(wú)磁性、強(qiáng)度高等特點(diǎn),成為制造航空器結(jié)構(gòu)件的重要材料[1-2]。鈦合金擁有較高的屈服強(qiáng)度和較低的彈性模量,室溫條件下塑性變形范圍小,在傳統(tǒng)成形過(guò)程中變形阻力大、成形性能差,如采用冷沖壓法沖壓,所需的沖壓噸位大,容易開(kāi)裂,一般采用高溫成形和熱成形,然而在高溫成形和熱成形過(guò)程中能量損耗大,并且會(huì)對(duì)材料表面造成氧化進(jìn)而導(dǎo)致材料利用率低,對(duì)模具性能影響較大[3]。鑒于以上因素,發(fā)展了一種將脈沖電流引入金屬塑性成形過(guò)程的成形方式,其原理是在金屬成形過(guò)程中,利用漂流電子促進(jìn)位錯(cuò)運(yùn)動(dòng),降低流動(dòng)應(yīng)力,提高成形性能,并在一定程度上細(xì)化晶粒,從而改善鈦合金成形性能,是一種很有前景的成形方法[4-6]。

    過(guò)去大量科研人員對(duì)鈦合金熱變形過(guò)程中的變形行為、組織演化和本構(gòu)建模進(jìn)行了系統(tǒng)深入的研究,相關(guān)結(jié)論對(duì)材料變形機(jī)理的理解和實(shí)際熱變形工藝制定起到一定作用。Wu等[7]在Gleeble–3800D模擬機(jī)上進(jìn)行了溫度為860~900 ℃、應(yīng)變速率為0.001~ 10 s?1的熱壓縮試驗(yàn),研究了Ti55531合金在β單相區(qū)的變形和顯微組織演變。王立穎等[8]通過(guò)熱模擬試驗(yàn)機(jī)對(duì)Ti–6Al–4V進(jìn)行熱壓剪處理,對(duì)比分析了不同熱機(jī)械加工參數(shù)下的流變行為和微觀組織演變規(guī)律。Wu等[9]在回火狀態(tài)、應(yīng)變速率為0.001~10 s?1的條件下對(duì)Ti55531進(jìn)行熱壓縮試驗(yàn),根據(jù)熱變形對(duì)實(shí)測(cè)應(yīng)力進(jìn)行修正,建立了整個(gè)應(yīng)變范圍內(nèi)的Arrhenius型本構(gòu)模型。肖寒等[10]通過(guò)使用雙曲正弦形式修正的Arrhenius本構(gòu)方程描述了TC4鈦合金高溫壓縮變形時(shí)的峰值應(yīng)力。朱鴻昌等[11]通過(guò)材料參數(shù)與真應(yīng)變之間的關(guān)系,利用Arrhenius本構(gòu)方程和參數(shù)建立了流動(dòng)應(yīng)力和變形溫度、應(yīng)變速率和真應(yīng)變之間的本構(gòu)關(guān)系。目前國(guó)內(nèi)外研究人員主要針對(duì)鈦合金熱變形宏觀和微觀變化進(jìn)行了大量研究,但是針對(duì)材料電輔助變形的研究還比較少。

    杜引等[12]基于鋁鋰合金電輔助單向拉伸試驗(yàn)的流動(dòng)應(yīng)力曲線,建立了2060–T8鋁鋰合金電塑性本構(gòu)方程,進(jìn)而為制定電輔助鋁鋰合金成形工藝提供理論依據(jù)。張寧等[13]建立了耦合脈沖電流密度和脈沖頻率的超塑性本構(gòu)方程,預(yù)測(cè)了1420鋁鋰合金在電致超塑性變形中流動(dòng)應(yīng)力的變化規(guī)律。彭書(shū)華[14]通過(guò)分析脈沖電流對(duì)金屬應(yīng)力–應(yīng)變的影響,選擇合理的經(jīng)驗(yàn)方程以簡(jiǎn)化金屬材料的應(yīng)力–應(yīng)變模型,得到了脈沖電流對(duì)金屬應(yīng)力–應(yīng)變的影響函數(shù),并結(jié)合不同加工硬化模型,提出了新的電致塑性效應(yīng)本構(gòu)方程—— H–V模型。盡管針對(duì)鋁鋰合金和鋁合金的本構(gòu)方程的研究較多,但是鈦合金與這些材料的變形機(jī)理有本質(zhì)區(qū)別,而針對(duì)Ti6554鈦合金電塑性的變形行為以及本構(gòu)方程的建立研究較少,因此建立鈦合金的電塑性本構(gòu)方程具有十分重要的意義。文中基于此背景開(kāi)展了鈦合金電塑性壓縮本構(gòu)模型的研究,為后續(xù)有限元模擬提供理論和數(shù)據(jù)支撐。

    1 試驗(yàn)

    1.1 材料

    試驗(yàn)材料采用Ti6554鈦合金,其化學(xué)成分如表1所示,采用4 mm×8 mm的圓柱體試樣,采用矩形方波直流脈沖電源對(duì)壓縮試樣施加電脈沖,方波電源可調(diào)節(jié)脈沖電流的電流密度、占空比、頻率。試驗(yàn)裝置及原理如圖1所示,由壓縮試驗(yàn)機(jī)的上下壓塊為試樣提供下壓力,電源電極為試樣提供脈沖電流,并通過(guò)FLIR熱成像儀檢測(cè)其溫度變化。

    表1 Ti6554鈦合金化學(xué)成分

    Tab.1 Chemical composition of Ti6554 titanium alloy wt.%

    圖1 電塑性壓縮試驗(yàn)裝置及原理

    1.2 方法

    在萬(wàn)能壓縮試驗(yàn)機(jī)上進(jìn)行試驗(yàn)。首先在常溫不施加脈沖電流時(shí),測(cè)得壓縮變形量為50%時(shí)Ti6554鈦合金圓柱試樣的屈服極限強(qiáng)度為1 277 MPa,然后進(jìn)行電塑性壓縮試驗(yàn),在應(yīng)力達(dá)到屈服極限、開(kāi)始轉(zhuǎn)向塑性階段時(shí)施加脈沖電流,分別測(cè)得電流密度為40、32、24 A/mm2,占空比為0.7、0.5、0.3,頻率為1、5、10 Hz,測(cè)量壓縮過(guò)程中的應(yīng)力–應(yīng)變曲線,具體試驗(yàn)方案如表2所示。

    表2 電塑性壓縮試驗(yàn)方案

    Tab.2 Electro-plastic compression experiment plan

    2 結(jié)果與分析

    Ti6554鈦合金電輔助壓縮在不同電流密度、占空比、脈沖頻率和應(yīng)變速率下的真實(shí)應(yīng)力–應(yīng)變曲線如圖2所示??梢钥闯?,在塑性階段加入脈沖電流后,流動(dòng)應(yīng)力迅速下降。由圖2a可以得出,在相同試驗(yàn)條件下,隨著電流密度的增大,塑性階段的流動(dòng)應(yīng)力減小,電流密度為40 A/mm2時(shí)的應(yīng)力最??;從圖2b可以發(fā)現(xiàn),脈沖頻率對(duì)金屬的流動(dòng)應(yīng)力影響較大,當(dāng)頻率為10 Hz時(shí)幾乎不產(chǎn)生應(yīng)力降,當(dāng)脈沖頻率為1 Hz時(shí),塑性階段流動(dòng)應(yīng)力最小,隨著頻率減小,塑性階段產(chǎn)生的流動(dòng)應(yīng)力減??;從圖2c可以得出,當(dāng)應(yīng)變速率增大時(shí),塑性階段的流動(dòng)應(yīng)力減??;從圖2d可以得出,隨著占空比的增大,金屬在塑性階段產(chǎn)生的流動(dòng)應(yīng)力減小,當(dāng)占空比為0.7時(shí),塑性階段的流動(dòng)應(yīng)力最小。

    圖2 不同脈沖電流因素下的應(yīng)力–應(yīng)變曲線

    Ti6554的壓縮變形進(jìn)入塑性階段后,試樣的流動(dòng)應(yīng)力在通電后迅速下降,隨著對(duì)試樣載荷的加大,應(yīng)力逐漸升高。這種現(xiàn)象可以通過(guò)鈦合金塑性變形時(shí)內(nèi)部發(fā)生的位錯(cuò)進(jìn)行解釋:壓縮過(guò)程的彈性階段沒(méi)有施加脈沖電流,此時(shí)可以發(fā)生滑移的位錯(cuò)密度較低,滑移時(shí)的動(dòng)力來(lái)源于試樣的壓縮力,此階段金屬的位錯(cuò)增殖還未發(fā)生;當(dāng)金屬的變形進(jìn)入塑性階段后,位錯(cuò)密度成倍增殖,脈沖電流可作用的晶界大幅增加,因此在通電的瞬間迅速促進(jìn)金屬塑性流動(dòng),宏觀上表現(xiàn)為“軟化”現(xiàn)象,進(jìn)而流動(dòng)應(yīng)力隨之降低[15]。金屬通電后,流動(dòng)應(yīng)力先迅速降低后逐漸上升,這是因?yàn)閷?duì)試樣持續(xù)通電會(huì)使試樣產(chǎn)生溫升,發(fā)生熱膨脹[16],進(jìn)而表現(xiàn)為金屬的變形抗力升高,流動(dòng)應(yīng)力升高。

    3 流動(dòng)應(yīng)力本構(gòu)模型建立

    3.1 唯象本構(gòu)模型

    唯象模型本構(gòu)是基于試驗(yàn)數(shù)據(jù)建立的流動(dòng)應(yīng)力與工藝參數(shù)之間的數(shù)學(xué)模型,該模型形式較為簡(jiǎn)單,不能反映材料變形的物理機(jī)制和物理意義。目前常用的本構(gòu)模型一般不考慮電流密度、占空比、頻率這些電參數(shù)對(duì)變形的影響,文中將電流密度、占空比、頻率參數(shù)引入唯象模型中,建立鈦合金電輔助壓縮的變形本構(gòu)模型。

    在金屬材料的研究過(guò)程中建立本構(gòu)模型預(yù)測(cè)應(yīng)力–應(yīng)變的變化規(guī)律尤為重要,在眾多金屬塑性研究中,有很多不同形式的本構(gòu)方程,其中大部分包含金屬的變形溫度、應(yīng)變速率、晶粒尺寸等影響因素,而脈沖電流與應(yīng)力、應(yīng)變之間的函數(shù)關(guān)系仍然處于探索階段。文中基于電輔助壓縮實(shí)驗(yàn)脈沖電流對(duì)應(yīng)力–應(yīng)變曲線的影響規(guī)律,建立唯象的電塑性本構(gòu)模型,該本構(gòu)模型是在Johnson–Cook模型的基礎(chǔ)上添加了多個(gè)函數(shù)來(lái)考慮電塑性效應(yīng),其中基礎(chǔ)的Johnson–Cook模型如式(1)所示。

    式中:1為電流密度;2為占空比;3為頻率;1為電流密度影響下的應(yīng)變強(qiáng)化系數(shù);2為占空比影響下的應(yīng)變強(qiáng)化系數(shù);3為頻率影響下的應(yīng)變強(qiáng)化系數(shù);1為電流密度影響下的應(yīng)變強(qiáng)化常數(shù);2為占空比影響下的應(yīng)變強(qiáng)化常數(shù);3為頻率影響下的應(yīng)變強(qiáng)化常數(shù);1為電流密度影響下的應(yīng)變硬化系數(shù);2為占空比影響下的應(yīng)變硬化系數(shù);3為頻率影響下的應(yīng)變硬化系數(shù);1為電流密度影響下的應(yīng)變硬化常數(shù);2為占空比影響下的應(yīng)變硬化常數(shù);3為頻率影響下的應(yīng)變硬化常數(shù);為流動(dòng)應(yīng)力。

    對(duì)求解時(shí),不考慮e對(duì)的影響,因此將式(5)取對(duì)數(shù)得到式(6)。

    ln()=ln+ln(6)

    針對(duì)40、32、24 A/mm2電流密度下–曲線,對(duì)ln–ln()擬合如圖3所示,由式(6)可知,ln為ln與ln()擬合曲線的截距,為ln與ln()的斜率。

    由截距l(xiāng)n和斜率得到電流密度1為40、32、24 A/mm2時(shí)的應(yīng)變強(qiáng)化系數(shù)和應(yīng)變硬化指數(shù),根據(jù)式(7)在不同電流密度下對(duì)應(yīng)變強(qiáng)化系數(shù)與電流密度1進(jìn)行擬合,如圖4a所示,得到1=1 860.71,1=?42 926.72。

    根據(jù)式(8)在不同電流密度下對(duì)應(yīng)變硬化指數(shù)n與電流密度I1進(jìn)行擬合,其擬合曲線如圖4b所示,得到k1=0.231,n1=?4.86。

    式(5)兩邊取對(duì)數(shù)后用lnln(A+Bε)除以電參數(shù)的乘積,得到電流系數(shù),因此電流密度影響下的電流系數(shù)=0.000 023,故得到電流密度的電塑性本構(gòu)模型如式(9)所示。

    在電輔助壓縮試驗(yàn)中發(fā)現(xiàn)占空比對(duì)流動(dòng)應(yīng)力的影響較大,占空比為一個(gè)工作周期內(nèi),脈沖所占用的時(shí)間與總時(shí)間的比值。試驗(yàn)方案具體參數(shù)如表2所示,采用占空比大小為0.3、0.5、0.7,其余各項(xiàng)參數(shù)不變。因此式(4)中當(dāng)2時(shí),式(4)為針對(duì)占空比的電塑性本構(gòu)模型表達(dá)式。對(duì)占空比影響的電塑性本構(gòu)參數(shù)進(jìn)行求解時(shí),用占空比為0.3、0.5、0.7下的曲線,在不考慮e對(duì)的影響下,對(duì)式(5)中的應(yīng)變強(qiáng)化系數(shù)和應(yīng)變硬化指數(shù)進(jìn)行擬合,得到lnln(?)擬合曲線如圖5所示,因此ln為式(7)中l(wèi)n與ln(?)的截距,為ln與ln(?)的斜率。

    圖5 不同占空比下的ln ε與ln(σ–A)擬合曲線

    由截距l(xiāng)n和斜率得到占空比2為0.3、0.5、0.7時(shí)對(duì)應(yīng)的應(yīng)變強(qiáng)化系數(shù)和應(yīng)變硬化指數(shù),根據(jù)式(10)對(duì)不同占空比下的應(yīng)變強(qiáng)化系數(shù)與占空比2進(jìn)行擬合,如圖6a所示,得到2=15 969.95,2=?1 397.90。

    根據(jù)式(11)對(duì)不同占空比下應(yīng)變硬化指數(shù)與占空比2進(jìn)行擬合,如圖6b所示,得到2=4.26,2=?0.85。

    式(5)兩邊取對(duì)數(shù)后用lnln(A+Bε)除以電參數(shù)的乘積,得到電流系數(shù),因此占空比影響下的電流系數(shù)=0.003 9,故得到占空比的電塑性本構(gòu)模型如式(12)所示。

    為建立頻率的電輔助壓縮本構(gòu)方程,試驗(yàn)方案采用的具體參數(shù)如表2所示,頻率大小為1、5、10 Hz,其余各項(xiàng)參數(shù)不變。因此式(4)中當(dāng)=3時(shí),式(4)為針對(duì)頻率的電塑性本構(gòu)模型表達(dá)式。用1、5、10 Hz頻率下的曲線對(duì)頻率的電塑性本構(gòu)參數(shù)求解,在不考慮e對(duì)的影響時(shí),對(duì)式(5)中的應(yīng)變強(qiáng)化系數(shù)和應(yīng)變硬化指數(shù)進(jìn)行擬合,得到參數(shù)和。對(duì)lnln(?)進(jìn)行擬合如圖7所示,ln為式(7)中l(wèi)n與ln(?)的截距,為ln與ln(?)的斜率。

    由截距l(xiāng)n和斜率得到頻率3為1、5、10 Hz時(shí)的應(yīng)變強(qiáng)化系數(shù)和應(yīng)變硬化指數(shù),根據(jù)式(13)對(duì)不同頻率下的應(yīng)變強(qiáng)化系數(shù)與頻率3進(jìn)行擬合,如圖8a所示,得到3=?267.77,3=6 491.19。

    圖6 不同占空比下的材料參數(shù)和值的確定

    圖7 針對(duì)于占空比的ln ε與ln(σ?A)擬合曲線

    根據(jù)式(14)對(duì)不同頻率下的應(yīng)變硬化指數(shù)與頻率3進(jìn)行擬合,如圖8b所示,得到3=?0.143,3=1.942。

    式(5)兩邊取對(duì)數(shù)后用lnln(A+Bε)除以電參數(shù)的乘積,得到電流系數(shù),因此頻率影響下的電流系數(shù)=0.004 22,故得到頻率的電塑性本構(gòu)模型,如式(15)所示。

    3.2 基于位錯(cuò)密度的物理機(jī)制本構(gòu)模型

    鑒于唯象本構(gòu)模型不能反映相應(yīng)的物理意義,文中以位錯(cuò)密度為中間變量,建立Ti6554鈦合金材料電塑性壓縮的物理機(jī)制本構(gòu)模型,以反映鈦合金電輔助壓縮變形流動(dòng)應(yīng)力變化的物理機(jī)制。其物理機(jī)制是金屬在施加脈沖電流后,定向漂流電子加快位錯(cuò)滑移速率,在外力作用下降低了應(yīng)變速率[17]。施加脈沖電流后的應(yīng)變速率是由外力壓縮產(chǎn)生的應(yīng)變速率和電塑性效應(yīng)引起的應(yīng)變速率同時(shí)作用而產(chǎn)生的,根據(jù)式(16)可得,在脈沖電流對(duì)試樣進(jìn)行壓縮試驗(yàn)的過(guò)程中,流動(dòng)應(yīng)力降低是因?yàn)槊}沖電流引起金屬的應(yīng)變速率發(fā)生變化,在試樣壓縮過(guò)程中應(yīng)變速率始終不變,在不通電的情況下壓縮材料產(chǎn)生的應(yīng)變速率全部源于外力,在有脈沖電流參與壓縮時(shí),應(yīng)變速率由外力和脈沖電流同時(shí)作用而產(chǎn)生變化的,從而宏觀上表現(xiàn)為金屬的電塑性效應(yīng)導(dǎo)致流動(dòng)應(yīng)力下降。由此可知,位錯(cuò)密度可以從微觀角度解釋脈沖電流對(duì)應(yīng)變速率的影響,因此文中基于位錯(cuò)密度的物理機(jī)制建立應(yīng)變速率的電輔助壓縮的本構(gòu)模型。

    文中對(duì)Ti6554在0.002、0.005、0.01 s?1不同應(yīng)變速率條件下進(jìn)行脈沖電流的壓縮試驗(yàn),采用的K–M位錯(cuò)密度模型如式(17)所示,設(shè)定初始位錯(cuò)密度為6×105m?2,隨著應(yīng)變?cè)黾樱诲e(cuò)密度也隨之遞增,每迭代一步應(yīng)變?cè)黾恿繛?.01,位錯(cuò)密度迭代改變量見(jiàn)式(17),位錯(cuò)密度與流動(dòng)應(yīng)力之間的關(guān)系見(jiàn)式(18)[18]。

    式中:Δρ為位錯(cuò)密度增量;Δ為應(yīng)變?cè)黾恿?;為泰勒因子?3.1;為柏氏矢量,=2.86×10?6;ρ為當(dāng)前位錯(cuò)密度;1為加工硬化系數(shù);2為動(dòng)態(tài)軟化系數(shù);1為初始位錯(cuò)密度,1=10?7m?2;為流動(dòng)應(yīng)力;為位錯(cuò)間強(qiáng)度,=0.3;為剪切模量,=67 kPa。

    在位錯(cuò)密度迭代過(guò)程中對(duì)參數(shù)進(jìn)行擬合,得到式(19)中1、1、2、2共4個(gè)參數(shù)的最優(yōu)解如下:1=424 015.3,1=?348.587,2=1 401,2=?29.918。

    3.3 模型準(zhǔn)確性驗(yàn)證

    3.3.1 唯象本構(gòu)模型驗(yàn)證

    文中所建立的本構(gòu)方程可以準(zhǔn)確預(yù)測(cè)電流密度為40、32、24 A/mm2,占空比為0.7、0.5、0.3,頻率為1、5、10 Hz條件下的流動(dòng)應(yīng)力,電流密度本構(gòu)模型的相對(duì)誤差=7.3%,頻率本構(gòu)模型的相對(duì)誤差=5.6%,占空比的本構(gòu)模型相對(duì)誤差=6.4%,總體誤差小于7%,相對(duì)誤差計(jì)算見(jiàn)式(20),因此唯象本構(gòu)模型能夠準(zhǔn)確預(yù)測(cè)不同電參數(shù)下材料的流動(dòng)應(yīng)力。為進(jìn)一步驗(yàn)證本構(gòu)模型的準(zhǔn)確性,對(duì)電流密度為28 A/mm2、頻率為0.1 Hz、占空比為0.9時(shí)的流動(dòng)應(yīng)力做出預(yù)測(cè),如圖9所示,隨著電流密度增大,金屬的流動(dòng)應(yīng)力減小,隨著頻率的減小,金屬的流動(dòng)應(yīng)力也減小,隨著占空比的增加,金屬的流動(dòng)應(yīng)力隨著減小,由圖9可知,基于文中所建立的本構(gòu)模型預(yù)測(cè)的流動(dòng)應(yīng)力趨勢(shì)與試驗(yàn)數(shù)據(jù)相同,符合文中所建立的唯象本構(gòu)模型規(guī)律,能夠?qū)A柱試樣電塑性壓縮做出較好的預(yù)測(cè)。

    式中:e為峰值應(yīng)力的預(yù)測(cè)值;p為峰值應(yīng)力的實(shí)驗(yàn)值;為實(shí)驗(yàn)次數(shù)。

    3.3.2 位錯(cuò)密度本構(gòu)模型驗(yàn)證

    文中基于位錯(cuò)密度所建立的本構(gòu)模型,可以準(zhǔn)確預(yù)測(cè)應(yīng)變速率為0.002、0.005、0.01 s?1時(shí)的流動(dòng)應(yīng)力,總體誤差=7.4%,如圖10所示,位錯(cuò)密度本構(gòu)模型具有較高的準(zhǔn)確度。為進(jìn)一步驗(yàn)證位錯(cuò)密度本構(gòu)模型的準(zhǔn)確性,文中在應(yīng)變速率為0.05 s?1下進(jìn)行預(yù)測(cè),可以得到當(dāng)應(yīng)變速率為0.05 s?1時(shí),與電輔助壓縮試驗(yàn)具有相同變化趨勢(shì),因此基于位錯(cuò)密度的本構(gòu)模型,能夠?qū)﹄娸o助壓縮的流動(dòng)應(yīng)力進(jìn)行準(zhǔn)確預(yù)測(cè)。

    圖10 基于所建位錯(cuò)密度模型對(duì)不同應(yīng)變速率下的流動(dòng)應(yīng)力預(yù)測(cè)

    4 結(jié)論

    1)針對(duì)Ti6554鈦合金在壓縮過(guò)程中的塑性階段對(duì)試樣施加脈沖電流,使材料在塑性階段出現(xiàn)較大的應(yīng)力降。其中電流密度越大、頻率越小、占空比越大以及應(yīng)變速率越小,所產(chǎn)生的應(yīng)力降越大。

    2)采用數(shù)學(xué)回歸擬合的方法建立了電流密度、頻率和占空比的唯象本構(gòu)模型,以及基于位錯(cuò)密度的物理機(jī)制本構(gòu)模型。

    3)對(duì)建立的本構(gòu)模型進(jìn)行驗(yàn)證,其中唯象本構(gòu)模型誤差為6.4%,位錯(cuò)密度本構(gòu)模型誤差為7.4%,總體誤差小于7%,因此建立的本構(gòu)模型對(duì)電輔助壓縮的流動(dòng)應(yīng)力能進(jìn)行準(zhǔn)確的預(yù)測(cè)。

    [1] LI Chang-min, HUANG Liang, ZHAO Ming-jie, et al. Hot Deformation Behavior and Mechanism of a New Metastable β Titanium Alloy Ti-6Cr-5Mo-5V-4Al in Single Phase Region[J]. Materials Science & Engineering A, 2021, 814: 141231.

    [2] 武川, 劉斌, 周宇杰, 等. Ti6554鈦合金高溫變形行為與微觀組織演化機(jī)制研究[J]. 精密成形工程, 2022, 14(1): 114-125.

    WU Chuan, LIU Bin, ZHOU Yu-jie, et al. Investigation on Hot Deformation Behavior and Microstructural Evolution of Ti6554 Titanium Alloy[J]. Journal of Netshape Forming Engineering, 2022, 14(1): 114-125.

    [3] LI Xi-feng, JI Bo-yu, ZHOU Qiang, et al. Influence of Grain Size on Electrically Assisted Tensile Behavior of Ti-6Al-4V Alloy[J]. Journal of Materials Engineering and Performance, 2016, 25(10): 4514-4520.

    [4] 張昕, 李智軍, 李宏偉, 等. 鎳基高溫合金初始織構(gòu)對(duì)電輔助拉伸宏微觀行為的影響[J]. 精密成形工程, 2019, 11(5): 32-36.

    ZHANG Xin, LI Zhi-jun, LI Hong-wei, et al. Effects of Initial Texture on Macro-Micro Behaviors of Ni-Based Superalloy during the Electrically-Assisted Tensile[J]. Journal of Netshape Forming Engineering, 2019, 11(5): 32-36.

    [5] 江雙雙, 湯澤軍, 杜浩, 等. 鈦合金電流輔助成形工藝研究進(jìn)展[J]. 精密成形工程, 2017, 9(2): 7-13.

    JIANG Shuang-shuang, TANG Ze-jun, DU Hao, et al. Research Progress of the Forming Process of Titanium Alloy Assisted by Current[J]. Journal of Netshape Forming Engineering, 2017, 9(2): 7-13.

    [6] RUSZKIEWICZ B J, RAGAI I, ROTH J T, et al. A Review of Electrically-Assisted Manufacturing with Emphasis on Modeling and Understanding of the Electroplastic Effect[J]. Journal of Manufacturing Science and Engineering, 2017, 139(11): 110801.

    [7] WU Chuan, HUANG Liang. Hot Deformation and Dynamic Recrystallization of a Near-Beta Titanium Alloy in the β Single Phase Region[J]. Vacuum, 2018, 156: 384-401.

    [8] 王立穎, 楊友, 劉春蘭. 熱機(jī)械加工Ti-6Al-4V鈦合金的流變行為和顯微組織演變[J]. 鍛壓技術(shù), 2020, 45(12): 183-190.

    WANG Li-ying, YANG You, LIU Chun-lan. Rheological Behavior and Microstructure Evolution on Ti-6Al-4V Titanium Alloy by Hot Machining[J]. Forging & Stamping Technology, 2020, 45(12): 183-190.

    [9] WU Chuan, HUANG Liang, LI Chang-min. Experimental Investigation on Dynamic Phase Transformation and Texture Evolution of Ti55531 High Strength Titanium Alloy during Hot Compression in the α+β Region[J]. Materials Science & Engineering A, 2020, 773: 138851.

    [10] 肖寒, 丁平, 張宏宇, 等. TC4鈦合金EB爐扁錠高溫壓縮變形行為和熱加工圖[J]. 材料熱處理學(xué)報(bào), 2021, 42(4): 175-182.

    XIAO Han, DING Ping, ZHANG Hong-yu, et al. High Temperature Compression Deformation Behavior and Hot Processing Map of TC4 Titanium Alloy Slab Smelted by Electron Beam Cold Bed Furnace[J]. Transactions of Materials and Heat Treatment, 2021, 42(4): 175-182.

    [11] 朱鴻昌, 羅軍明, 朱知壽. TB17鈦合金高溫壓縮變形行為[J]. 航空材料學(xué)報(bào), 2019, 39(3): 44-52.

    ZHU Hong-chang, LUO Jun-ming, ZHU Zhi-shou. High Temperature Compression Deformation Behavior of TB17 Titanium Alloy[J]. Journal of Aeronautical Materials, 2019, 39(3): 44-52.

    [12] 杜引, 趙亦希, 于忠奇, 等. 2060-T8鋁鋰合金電致塑性本構(gòu)方程[J]. 塑性工程學(xué)報(bào), 2017, 24(1): 133-139.

    DU Yin, ZHAO Yi-xi, YU Zhong-qi, et al. Constitutive Equation of Electroplastic Effect for 2060-T8 Al-Li Alloy[J]. Journal of Plasticity Engineering, 2017, 24(1): 133-139.

    [13] 張寧, 張艷苓, 畢靜, 等. 1420鋁鋰合金電致超塑性本構(gòu)方程[J]. 鍛壓技術(shù), 2015, 40(5): 63-68.

    ZHANG Ning, ZHANG Yan-ling, BI Jing, et al. Constitutive Equation of Electro-Superplastic for 1420 Al- Li Alloy[J]. Forging & Stamping Technology, 2015, 40(5): 63-68.

    [14] 彭書(shū)華. 電致塑性效應(yīng)本構(gòu)方程的研究[J]. 塑性工程學(xué)報(bào), 2018, 25(2): 202-206.

    PENG Shu-hua. Study on Constitutive Equation of Electroplastic Effect[J]. Journal of Plasticity Engineering, 2018, 25(2): 202-206.

    [15] 武偉超, 王永軍, 孫寶龍, 等. Al-Cu合金脈沖電塑性效應(yīng)研究[J]. 熱加工工藝, 2014, 43(2): 82-85.

    WU Wei-chao, WANG Yong-jun, SUN Bao-long, et al. Study on Electroplastic Effect of Al-Cu Alloy[J]. Hot Working Technology, 2014, 43(2): 82-85.

    [16] 楊以鵬. TC4曲面構(gòu)件自阻加熱電塑性成形工藝研究[D]. 秦皇島: 燕山大學(xué), 2019.

    YANG Yi-peng. Study on Resistance Heating and Electro Plastic Forming Process of TC4 Curved Surfaces[D]. Qinhuangdao: Yanshan University, 2019.

    [17] 李大龍. 電塑性效應(yīng)中金屬流動(dòng)應(yīng)力的理論及實(shí)驗(yàn)研究[D]. 秦皇島: 燕山大學(xué), 2014.

    LI Da-long. Theoretical and Experimental Study on Metal's Flow Stress in Electroplastic Effect[D]. Qinhuangdao: Yanshan University, 2014.

    [18] YANAGIMOTO J, KARHAUSEN K, BRAND A J, et al. Incremental Formulation for the Prediction of Flow Stress and Microstructural Change in Hot Forming[J]. Journal of Manufacturing Science and Engineering, Transactions of the ASME, 1998, 120(2): 318-322.

    Establishment and Application of Ti6554 Titanium Alloy Pulse Current Assisted Compression Constitutive Model

    LIU Bin, WU Chuan, ZHOU Yu-jie

    (National-local Joint Engineering Laboratory of Intelligent Manufacturing Oriented Automobile Die & Mould, Tianjin University of Technology and Education, Tianjin 300222, China)

    This paper is to establish the electro-assisted compression constitutive model of Ti6554 titanium alloy and verify the constitutive model based on the electro-assisted compression experiment. For Ti6554 titanium alloy material, under certain pulse current parameters (current density is 40, 32, 24 A/mm2, frequency is 1, 5, 10 Hz, duty cycle is 0.7, 0.5, 0.3) and strain rate is 0.002, 0.005, 0.01 s?1and the depression amount is 50%, the electro-assisted compression experiment is carried out, and the constitutive model of electro-assisted compression is established based on the phenomenological and physical mechanisms. The constitutive model parameters are fitted with experiment data by the regression fitting method and the specific equation is obtained. Based on the stress-strain curve of the electro-assisted compression experiment, the influence of pulse current parameters on the stress-strain is analyzed. The result shows that the relative error analysis is carried out between the experimental value measured by the electro-assisted compression experiment and the predicted value obtained by the constitutive model, and the maximum relative error is less than 7%.The constitutive model established by the regression fitting method has a small relative error and can accurately predict the flow stress of the material under different pulse current parameters.

    electro-plasticity; Ti6554; constitutive equation; pulse current

    10.3969/j.issn.1674-6457.2022.05.005

    TG306

    A

    1674-6457(2022)05-0027-09

    2021–07–17

    國(guó)家自然科學(xué)基金面上項(xiàng)目(52075386);中國(guó)博士后科學(xué)基金面上項(xiàng)目(2020M672309);天津市教委科研項(xiàng)目(2020KJ107)

    劉斌(1996—),男,碩士生,主要研究方向?yàn)榻饘俨牧纤苄猿尚喂に嚒?/p>

    武川(1981—),男,博士,講師,主要研究方向?yàn)榻饘俨牧纤苄猿尚喂に嚒?/p>

    責(zé)任編輯:蔣紅晨

    猜你喜歡
    脈沖電流本構(gòu)鈦合金
    離心SC柱混凝土本構(gòu)模型比較研究
    鋸齒形結(jié)構(gòu)面剪切流變及非線性本構(gòu)模型分析
    “神的金屬”鈦合金SHINE YOUR LIFE
    鈦合金板鍛造的工藝實(shí)踐
    四川冶金(2017年6期)2017-09-21 00:52:30
    一種新型超固結(jié)土三維本構(gòu)模型
    脈沖電流密度對(duì)納米晶鎳鍍層結(jié)構(gòu)及性能的影響
    醫(yī)用鈦合金的研究與應(yīng)用
    低頻脈沖電流經(jīng)皮刺激足三里穴對(duì)高強(qiáng)度軍訓(xùn)士兵抗疲勞作用觀察
    一種新型的雙極性脈沖電流源
    TC17鈦合金超大規(guī)格棒材的制備
    99热这里只有精品一区| 免费看光身美女| 免费观看在线日韩| 日韩中字成人| 青春草亚洲视频在线观看| 日韩av在线大香蕉| 日韩欧美精品免费久久| 国产成人a∨麻豆精品| 一级片'在线观看视频| 99久久精品热视频| 免费观看a级毛片全部| 午夜免费观看性视频| 色视频www国产| 成人av在线播放网站| 能在线免费看毛片的网站| 麻豆成人午夜福利视频| 亚洲天堂国产精品一区在线| 99视频精品全部免费 在线| 街头女战士在线观看网站| 亚洲国产欧美人成| 噜噜噜噜噜久久久久久91| 最近视频中文字幕2019在线8| 亚洲乱码一区二区免费版| 黄色日韩在线| 亚洲在线自拍视频| 天堂俺去俺来也www色官网 | 久久亚洲国产成人精品v| 身体一侧抽搐| 国产成人精品婷婷| 又粗又硬又长又爽又黄的视频| 又大又黄又爽视频免费| 你懂的网址亚洲精品在线观看| 又爽又黄a免费视频| 蜜桃亚洲精品一区二区三区| 久99久视频精品免费| 亚洲精品国产av成人精品| 青春草国产在线视频| 中文字幕久久专区| 国产午夜精品一二区理论片| 久久亚洲国产成人精品v| 国产精品久久久久久久久免| 国产精品嫩草影院av在线观看| 综合色av麻豆| 国产精品日韩av在线免费观看| 亚洲熟妇中文字幕五十中出| 最近视频中文字幕2019在线8| 国产亚洲精品久久久com| 国产精品一区www在线观看| 国产精品美女特级片免费视频播放器| 99九九线精品视频在线观看视频| 又大又黄又爽视频免费| 国产黄频视频在线观看| 亚洲欧美精品专区久久| 国产真实伦视频高清在线观看| 777米奇影视久久| 六月丁香七月| 国产乱人偷精品视频| 成人高潮视频无遮挡免费网站| 国产白丝娇喘喷水9色精品| 亚洲国产高清在线一区二区三| 国产一区二区三区综合在线观看 | 午夜老司机福利剧场| 午夜福利视频精品| 成人鲁丝片一二三区免费| 极品少妇高潮喷水抽搐| 成年女人在线观看亚洲视频 | 噜噜噜噜噜久久久久久91| 国产欧美另类精品又又久久亚洲欧美| 亚洲怡红院男人天堂| 亚州av有码| 欧美成人a在线观看| 精品一区二区三区视频在线| 韩国高清视频一区二区三区| 国产老妇伦熟女老妇高清| 丰满少妇做爰视频| 亚洲av电影在线观看一区二区三区 | 天堂俺去俺来也www色官网 | 伦精品一区二区三区| 午夜福利在线观看免费完整高清在| 国产爱豆传媒在线观看| 婷婷色麻豆天堂久久| 三级经典国产精品| 高清毛片免费看| 亚洲av不卡在线观看| 99久久人妻综合| 青春草视频在线免费观看| 精品酒店卫生间| 精品一区二区三区人妻视频| 国产精品1区2区在线观看.| 一级毛片aaaaaa免费看小| 国产久久久一区二区三区| 精品亚洲乱码少妇综合久久| 国产91av在线免费观看| 精品久久久久久成人av| 免费av观看视频| 肉色欧美久久久久久久蜜桃 | 联通29元200g的流量卡| 国产精品伦人一区二区| 日韩 亚洲 欧美在线| 久热久热在线精品观看| 欧美极品一区二区三区四区| 老女人水多毛片| 精品久久久噜噜| 成人一区二区视频在线观看| 久久精品国产亚洲av涩爱| 免费看a级黄色片| 国产精品爽爽va在线观看网站| 精品欧美国产一区二区三| 十八禁国产超污无遮挡网站| 熟女电影av网| 国产精品国产三级国产av玫瑰| 一区二区三区乱码不卡18| 色尼玛亚洲综合影院| 高清视频免费观看一区二区 | 国产极品天堂在线| 两个人的视频大全免费| 波野结衣二区三区在线| 白带黄色成豆腐渣| 成人亚洲精品一区在线观看 | 国产片特级美女逼逼视频| 成人亚洲精品av一区二区| 最近最新中文字幕大全电影3| 成年免费大片在线观看| 亚洲成人一二三区av| 中文字幕人妻熟人妻熟丝袜美| 国产成人午夜福利电影在线观看| 午夜免费男女啪啪视频观看| 看免费成人av毛片| 亚洲国产精品专区欧美| 能在线免费观看的黄片| 国产伦在线观看视频一区| 精华霜和精华液先用哪个| 午夜福利视频精品| 免费少妇av软件| 欧美区成人在线视频| 欧美日韩精品成人综合77777| 日韩制服骚丝袜av| 在线观看一区二区三区| av国产免费在线观看| 亚洲成人中文字幕在线播放| 亚洲av.av天堂| 亚洲最大成人av| 婷婷色av中文字幕| 国产黄片美女视频| 国产不卡一卡二| 国产亚洲最大av| 亚洲精品色激情综合| 亚洲人成网站在线观看播放| 干丝袜人妻中文字幕| 99久国产av精品国产电影| 我要看日韩黄色一级片| 久久久久九九精品影院| 蜜桃亚洲精品一区二区三区| 亚洲,欧美,日韩| 国产高清不卡午夜福利| 国产老妇伦熟女老妇高清| 女人被狂操c到高潮| 日韩一区二区视频免费看| 熟女电影av网| 久久99精品国语久久久| or卡值多少钱| 男人爽女人下面视频在线观看| 成人午夜精彩视频在线观看| 两个人的视频大全免费| 我的女老师完整版在线观看| 国产精品久久久久久av不卡| av黄色大香蕉| h日本视频在线播放| 国产成人免费观看mmmm| 亚洲欧美成人精品一区二区| 亚洲乱码一区二区免费版| 日韩av免费高清视频| 亚洲一区高清亚洲精品| 亚洲国产日韩欧美精品在线观看| 在线观看一区二区三区| 日本av手机在线免费观看| 国产精品国产三级国产专区5o| 国内精品一区二区在线观看| 亚洲欧美成人综合另类久久久| 毛片女人毛片| 国产一区有黄有色的免费视频 | 国产不卡一卡二| 超碰av人人做人人爽久久| av女优亚洲男人天堂| 中文资源天堂在线| 亚洲激情五月婷婷啪啪| 夫妻午夜视频| 亚洲欧美一区二区三区国产| 久久久精品免费免费高清| 免费观看a级毛片全部| 国产综合精华液| 亚洲精品成人av观看孕妇| 国产淫语在线视频| 午夜福利高清视频| 精品午夜福利在线看| 韩国高清视频一区二区三区| 搞女人的毛片| 搡女人真爽免费视频火全软件| 成人漫画全彩无遮挡| 成人av在线播放网站| 亚洲国产精品国产精品| 麻豆久久精品国产亚洲av| 国产精品人妻久久久影院| 97热精品久久久久久| 亚洲久久久久久中文字幕| 伊人久久精品亚洲午夜| 国内揄拍国产精品人妻在线| 乱系列少妇在线播放| 欧美97在线视频| 你懂的网址亚洲精品在线观看| 少妇被粗大猛烈的视频| 搡老妇女老女人老熟妇| 如何舔出高潮| 大香蕉久久网| 搡老乐熟女国产| 亚洲精品第二区| 亚洲丝袜综合中文字幕| 中文字幕av成人在线电影| 国产精品女同一区二区软件| 成人av在线播放网站| 联通29元200g的流量卡| 纵有疾风起免费观看全集完整版 | 成人性生交大片免费视频hd| 天堂√8在线中文| 天天躁日日操中文字幕| 亚洲天堂国产精品一区在线| av天堂中文字幕网| 青春草亚洲视频在线观看| 免费看av在线观看网站| 看免费成人av毛片| 久久精品久久久久久噜噜老黄| 丝袜美腿在线中文| 99久久精品国产国产毛片| 成人美女网站在线观看视频| av.在线天堂| 97超视频在线观看视频| 如何舔出高潮| 小蜜桃在线观看免费完整版高清| 久久热精品热| av在线天堂中文字幕| 成人鲁丝片一二三区免费| kizo精华| 久久国产乱子免费精品| 日韩亚洲欧美综合| 亚洲丝袜综合中文字幕| 国产av在哪里看| 日韩强制内射视频| 久久人人爽人人片av| 99久久精品一区二区三区| 成人亚洲欧美一区二区av| 国产成年人精品一区二区| 成年人午夜在线观看视频 | 国产有黄有色有爽视频| 色吧在线观看| 成年人午夜在线观看视频 | 久久精品久久久久久噜噜老黄| 成人特级av手机在线观看| 18禁动态无遮挡网站| 黄色配什么色好看| 又粗又硬又长又爽又黄的视频| 国产乱来视频区| 人人妻人人看人人澡| 国产精品一及| 国产精品国产三级国产专区5o| 国产精品不卡视频一区二区| 少妇人妻一区二区三区视频| 国产亚洲av嫩草精品影院| 久久久欧美国产精品| 黄片无遮挡物在线观看| 性插视频无遮挡在线免费观看| av专区在线播放| 午夜日本视频在线| 国产av国产精品国产| 午夜久久久久精精品| 国产成人精品久久久久久| 视频中文字幕在线观看| 日韩欧美国产在线观看| 看非洲黑人一级黄片| 一区二区三区乱码不卡18| 国产精品不卡视频一区二区| 亚洲欧洲国产日韩| 非洲黑人性xxxx精品又粗又长| 国产毛片a区久久久久| 亚洲内射少妇av| 亚洲一区高清亚洲精品| 女人久久www免费人成看片| 亚洲最大成人av| 久久草成人影院| 中文字幕av成人在线电影| 美女xxoo啪啪120秒动态图| 免费看av在线观看网站| 日韩av不卡免费在线播放| 超碰av人人做人人爽久久| 久久精品熟女亚洲av麻豆精品 | 视频中文字幕在线观看| 丝瓜视频免费看黄片| 最近中文字幕2019免费版| 亚洲av一区综合| 日韩欧美一区视频在线观看 | 69人妻影院| 国产精品一区www在线观看| 少妇熟女aⅴ在线视频| 激情五月婷婷亚洲| 高清毛片免费看| 国产午夜福利久久久久久| 在线观看免费高清a一片| 少妇的逼水好多| 欧美日韩在线观看h| 能在线免费观看的黄片| 伊人久久精品亚洲午夜| 在线免费观看的www视频| 特大巨黑吊av在线直播| 久久亚洲国产成人精品v| 在线免费观看的www视频| 99热6这里只有精品| 99re6热这里在线精品视频| 哪个播放器可以免费观看大片| 精品久久久久久久末码| 国产成人午夜福利电影在线观看| 五月玫瑰六月丁香| 韩国高清视频一区二区三区| 天堂中文最新版在线下载 | 国产成人免费观看mmmm| 麻豆乱淫一区二区| 一个人看视频在线观看www免费| 久久99热这里只有精品18| 久热久热在线精品观看| 三级国产精品欧美在线观看| 国产免费视频播放在线视频 | 日韩欧美三级三区| 日韩一区二区三区影片| 久久久午夜欧美精品| 亚洲高清免费不卡视频| 99久久九九国产精品国产免费| 国产精品一及| 久久人人爽人人片av| 国产av码专区亚洲av| av一本久久久久| 韩国av在线不卡| 亚洲,欧美,日韩| 久久久久久久久大av| 在现免费观看毛片| 在线天堂最新版资源| 狂野欧美激情性xxxx在线观看| 成人鲁丝片一二三区免费| 九九在线视频观看精品| 乱人视频在线观看| av专区在线播放| 男人舔女人下体高潮全视频| 五月天丁香电影| 精品熟女少妇av免费看| 国产精品精品国产色婷婷| 亚洲三级黄色毛片| 亚洲国产av新网站| 亚洲av不卡在线观看| 纵有疾风起免费观看全集完整版 | 亚洲成人一二三区av| 国产av在哪里看| 久久99热这里只频精品6学生| 久久久a久久爽久久v久久| 寂寞人妻少妇视频99o| 日韩欧美精品免费久久| 特级一级黄色大片| av在线天堂中文字幕| 丝袜喷水一区| 国产 亚洲一区二区三区 | 日本欧美国产在线视频| 国产有黄有色有爽视频| av在线播放精品| 国产精品麻豆人妻色哟哟久久 | 在线免费观看不下载黄p国产| 三级国产精品欧美在线观看| 国产亚洲av嫩草精品影院| 成人午夜精彩视频在线观看| 亚洲av电影在线观看一区二区三区 | 欧美激情国产日韩精品一区| 欧美bdsm另类| 国产中年淑女户外野战色| 中文字幕av成人在线电影| 日韩av在线大香蕉| av福利片在线观看| 最近中文字幕高清免费大全6| 搞女人的毛片| 日韩欧美三级三区| 亚洲三级黄色毛片| 国产视频首页在线观看| 精品99又大又爽又粗少妇毛片| 欧美xxⅹ黑人| 日本欧美国产在线视频| 99久久九九国产精品国产免费| 免费观看无遮挡的男女| 成人综合一区亚洲| 亚洲国产欧美人成| 中文在线观看免费www的网站| 99久久精品热视频| 秋霞在线观看毛片| 少妇人妻一区二区三区视频| 成人av在线播放网站| .国产精品久久| 99九九线精品视频在线观看视频| 天堂网av新在线| 精品国产三级普通话版| 男女那种视频在线观看| 高清日韩中文字幕在线| 一二三四中文在线观看免费高清| 免费人成在线观看视频色| 26uuu在线亚洲综合色| 欧美日韩亚洲高清精品| ponron亚洲| 最近2019中文字幕mv第一页| 高清午夜精品一区二区三区| 嫩草影院入口| 国产伦在线观看视频一区| 久久久久久伊人网av| 国产大屁股一区二区在线视频| 日本午夜av视频| 欧美极品一区二区三区四区| 日日啪夜夜撸| 一级av片app| av网站免费在线观看视频 | 国产探花在线观看一区二区| 亚洲内射少妇av| 乱码一卡2卡4卡精品| 色吧在线观看| 亚洲精品乱码久久久v下载方式| 七月丁香在线播放| 日本三级黄在线观看| 亚洲av男天堂| 如何舔出高潮| 午夜福利网站1000一区二区三区| 欧美日韩国产mv在线观看视频 | 久久精品国产亚洲av天美| 人妻夜夜爽99麻豆av| 日韩 亚洲 欧美在线| 波野结衣二区三区在线| 午夜精品在线福利| 大香蕉久久网| 一本久久精品| 久久鲁丝午夜福利片| 午夜福利视频精品| 国产伦在线观看视频一区| 一区二区三区四区激情视频| 免费看不卡的av| 久久99蜜桃精品久久| 日韩欧美精品免费久久| h日本视频在线播放| 毛片女人毛片| 在线天堂最新版资源| 乱系列少妇在线播放| 亚洲第一区二区三区不卡| 午夜爱爱视频在线播放| 日韩一本色道免费dvd| 午夜福利在线观看吧| 欧美xxⅹ黑人| 亚洲aⅴ乱码一区二区在线播放| av天堂中文字幕网| 舔av片在线| 日本熟妇午夜| 插阴视频在线观看视频| 人人妻人人澡人人爽人人夜夜 | 干丝袜人妻中文字幕| 国内少妇人妻偷人精品xxx网站| 国产高清国产精品国产三级 | 欧美精品国产亚洲| 小蜜桃在线观看免费完整版高清| 天堂√8在线中文| 免费看不卡的av| 亚洲熟妇中文字幕五十中出| 在线播放无遮挡| 免费大片18禁| 天堂影院成人在线观看| 一级片'在线观看视频| 亚洲欧美成人精品一区二区| 精品一区在线观看国产| 91狼人影院| 久久精品综合一区二区三区| 不卡视频在线观看欧美| 一个人免费在线观看电影| 三级经典国产精品| 国产免费视频播放在线视频 | 国产精品精品国产色婷婷| 免费大片黄手机在线观看| av在线天堂中文字幕| 岛国毛片在线播放| 免费观看无遮挡的男女| 国产一级毛片七仙女欲春2| 国产在视频线精品| 亚洲精品自拍成人| 国产成人aa在线观看| 亚洲自拍偷在线| 国产中年淑女户外野战色| 在线观看一区二区三区| 天堂网av新在线| 菩萨蛮人人尽说江南好唐韦庄| 日韩视频在线欧美| av在线天堂中文字幕| 国产成人午夜福利电影在线观看| 欧美成人a在线观看| 免费看光身美女| 18禁裸乳无遮挡免费网站照片| 国产成人freesex在线| 一本久久精品| av天堂中文字幕网| 最近2019中文字幕mv第一页| 欧美日韩一区二区视频在线观看视频在线 | 91aial.com中文字幕在线观看| 青春草国产在线视频| 亚洲av成人精品一二三区| 日韩视频在线欧美| 国产精品人妻久久久影院| 欧美日本视频| 毛片一级片免费看久久久久| 欧美日韩综合久久久久久| 国产黄色视频一区二区在线观看| 亚洲精品亚洲一区二区| 精品久久久久久成人av| 男人舔女人下体高潮全视频| 日韩视频在线欧美| 日韩av免费高清视频| 久久久久国产网址| 少妇裸体淫交视频免费看高清| 国产有黄有色有爽视频| 人妻系列 视频| 亚洲成人精品中文字幕电影| 在线观看人妻少妇| 久久久精品欧美日韩精品| 久久久久久久久中文| 热99在线观看视频| 色视频www国产| ponron亚洲| 久久久久久久大尺度免费视频| 亚洲av成人av| 99久久精品一区二区三区| 日日干狠狠操夜夜爽| 内射极品少妇av片p| 婷婷色综合www| 国产午夜精品久久久久久一区二区三区| 三级国产精品片| 最近最新中文字幕免费大全7| 日本一本二区三区精品| 一夜夜www| 精品少妇黑人巨大在线播放| 大香蕉97超碰在线| 日韩欧美精品v在线| 婷婷色av中文字幕| 亚洲自偷自拍三级| 国产久久久一区二区三区| 色吧在线观看| 色5月婷婷丁香| 国产精品女同一区二区软件| 亚洲第一区二区三区不卡| 日日干狠狠操夜夜爽| 国产黄片美女视频| 亚洲欧美成人综合另类久久久| 美女黄网站色视频| 日本一二三区视频观看| 麻豆成人av视频| 成人特级av手机在线观看| 十八禁网站网址无遮挡 | 久久久久久久亚洲中文字幕| 精品午夜福利在线看| 午夜久久久久精精品| 国产精品福利在线免费观看| 久久久久久久大尺度免费视频| 欧美日韩在线观看h| 精品欧美国产一区二区三| 免费看光身美女| 狂野欧美激情性xxxx在线观看| 直男gayav资源| 亚洲自偷自拍三级| 欧美成人午夜免费资源| 91久久精品电影网| 非洲黑人性xxxx精品又粗又长| 国产高潮美女av| 欧美丝袜亚洲另类| 你懂的网址亚洲精品在线观看| 青春草亚洲视频在线观看| 看免费成人av毛片| 狂野欧美激情性xxxx在线观看| 亚洲性久久影院| 国产高清国产精品国产三级 | 色5月婷婷丁香| 日韩欧美三级三区| 超碰97精品在线观看| 日韩人妻高清精品专区| 欧美日韩亚洲高清精品| 美女国产视频在线观看| 亚洲精华国产精华液的使用体验| 搡老妇女老女人老熟妇| 成人毛片a级毛片在线播放| 日韩一区二区三区影片| 久久久久久久久久黄片| 亚洲精品乱码久久久v下载方式| 青春草亚洲视频在线观看| 欧美人与善性xxx| 五月天丁香电影| 午夜福利在线观看免费完整高清在| 99久国产av精品国产电影| 日韩精品有码人妻一区| 成人午夜精彩视频在线观看| 久久久精品欧美日韩精品| 精品熟女少妇av免费看| 男插女下体视频免费在线播放| 街头女战士在线观看网站| 少妇熟女欧美另类| 免费电影在线观看免费观看| 亚洲在久久综合| av国产久精品久网站免费入址| 街头女战士在线观看网站| 成人亚洲精品av一区二区| 久久久久久国产a免费观看| 亚洲国产日韩欧美精品在线观看| 亚洲最大成人手机在线| 国内少妇人妻偷人精品xxx网站| 丝瓜视频免费看黄片|