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

    基于能量守恒原理的金屬壓縮流動(dòng)應(yīng)力修正方法

    2017-10-13 05:55:42黃西成胡文軍
    中國有色金屬學(xué)報(bào) 2017年7期
    關(guān)鍵詞:實(shí)驗(yàn)

    賈 東,黃西成,莫 軍,胡文軍

    ?

    基于能量守恒原理的金屬壓縮流動(dòng)應(yīng)力修正方法

    賈 東,黃西成,莫 軍,胡文軍

    (中國工程物理研究院總體工程研究所,綿陽 621900)

    通過軸對(duì)稱缺口圓柱試件設(shè)計(jì),消除了端面摩擦對(duì)金屬(特別是延性金屬)試件壓縮變形行為的影響?;趬嚎s過程中外力功與試件內(nèi)能之間的能量守恒原理,推導(dǎo)了壓縮流動(dòng)應(yīng)力的近似修正公式,并給出了相應(yīng)的迭代修正過程。針對(duì)鋁合金2024,利用修正后的流動(dòng)應(yīng)力?應(yīng)變關(guān)系對(duì)壓縮過程中試件的最小橫截面直徑和內(nèi)能進(jìn)行了預(yù)測,并與實(shí)驗(yàn)結(jié)果進(jìn)行了比較。研究結(jié)果表明:在整個(gè)壓縮過程中,試件的最小橫截面直徑的預(yù)測值與實(shí)驗(yàn)值吻合,且試件內(nèi)能與實(shí)驗(yàn)外力功始終保持一致,從而驗(yàn)證了該修正方法的有效性。

    流動(dòng)應(yīng)力;端面摩擦;缺口圓柱試件;能量守恒;迭代修正

    有限元數(shù)值模擬是金屬塑性成形過程分析的重要研究手段,而準(zhǔn)確的流動(dòng)應(yīng)力?應(yīng)變關(guān)系是對(duì)該過程進(jìn)行數(shù)值分析的重要前提。由于金屬成形過程主要處于壓應(yīng)力狀態(tài),且壓縮實(shí)驗(yàn)可以獲取更大塑性變形量的應(yīng)力?應(yīng)變曲線,所以如何獲取準(zhǔn)確的壓縮流動(dòng)應(yīng)力?應(yīng)變曲線是非常重要的。

    目前,金屬材料壓縮應(yīng)力?應(yīng)變關(guān)系主要通過軸對(duì)稱圓柱試件壓縮實(shí)驗(yàn)獲取,但由于壓頭與試件端面之間的摩擦約束,試件在整個(gè)壓縮過程中無法保持均勻應(yīng)力狀態(tài),這對(duì)金屬材料(特別是延性金屬)壓縮流動(dòng)應(yīng)力?應(yīng)變曲線的準(zhǔn)確性影響很大[1]。消除或修正這種端面摩擦約束影響的方法主要有以下幾種:一是通過在接觸面上添加潤滑劑或在試件端面設(shè)計(jì)凹槽等方式來直接減小摩擦作用[2],這種方法對(duì)于大變形的延性金屬效果并不理想,當(dāng)壓縮載荷逐漸增大時(shí),兩接觸面間的潤滑作用會(huì)逐漸減弱,凹槽面會(huì)產(chǎn)生明顯的塑性變形,從而影響對(duì)摩擦的消除效果[3];二是COOK等[4]提出的外推法,通過多組不同長徑比的圓柱試件壓縮實(shí)驗(yàn)找到應(yīng)力隨長徑比變化的規(guī)律,再外推得到長徑比無窮大(無摩擦狀態(tài))時(shí)的應(yīng)力,該方法避免了對(duì)壓縮端面真實(shí)摩擦約束情況的分析,但其外推結(jié)果精度受限于試件的長徑比范圍,當(dāng)長徑比較大時(shí)試件易發(fā)生屈曲,從而導(dǎo)致實(shí)驗(yàn)無效;三是利用試件變形特征求解端面摩擦因數(shù),進(jìn)而對(duì)流動(dòng)應(yīng)力進(jìn)行修正,這其中以圓柱試件鼓度法[5?10]以及圓環(huán)鐓粗法[11?15]最為典型,在這些方法中,端面摩擦行為常常被假設(shè)為特定的線性摩擦模型,并以平均摩擦因數(shù)表征整個(gè)接觸面的摩擦特性。但在壓縮過程中,隨著載荷的增大,試件端面的不同區(qū)域常常處于不同的摩擦類型,且不同類型之間還會(huì)發(fā)生相互轉(zhuǎn)換,摩擦常表現(xiàn)出非線 性[16],所以端面摩擦行為描述的準(zhǔn)確性將會(huì)直接影響到流動(dòng)應(yīng)力的修正計(jì)算。除此以外,目前最為常用的是實(shí)驗(yàn)和數(shù)值模擬相結(jié)合的各種逆向求解方法[17?19],但這些方法需要給出應(yīng)力、幾何參數(shù)以及摩擦因數(shù)等之間準(zhǔn)確的數(shù)學(xué)函數(shù)關(guān)系??偟膩碚f,上述各種實(shí)驗(yàn)及有限元方法對(duì)端面摩擦約束的修正都有其自身的缺點(diǎn)或局限性,這是由于端面摩擦行為的復(fù)雜性所造成的。

    本文作者為了避免對(duì)試件端面復(fù)雜摩擦行為的分析,采用新型的缺口圓柱試件設(shè)計(jì)來消除端面摩擦約束對(duì)試件壓縮變形的影響。基于缺口試件壓縮過程中的能量守恒原理推導(dǎo)壓縮流動(dòng)應(yīng)力的近似修正公式,同時(shí)結(jié)合有限元分析,對(duì)流動(dòng)應(yīng)力進(jìn)行迭代修正,獲取金屬材料準(zhǔn)確的壓縮流動(dòng)應(yīng)力?應(yīng)變關(guān)系。

    1 流動(dòng)應(yīng)力修正方法

    1.1 近似修正公式

    在端面摩擦等約束存在時(shí),壓縮外力所做的功F將轉(zhuǎn)化為兩部分能量:一部分轉(zhuǎn)化為試件的塑性應(yīng)變能P儲(chǔ)存在試件中;另一部分則轉(zhuǎn)化為試件端面與壓頭之間的摩擦功f。所以外力功可以表示為

    式中:為壓縮過程中試件所受的外力,為試件的位移。而塑性變形能和壓縮過程中摩擦作用所耗散的能量可以分別表示為

    (2)

    式中:0為試件的應(yīng)變能密度,是試件的體積,為試件端面的摩擦因數(shù),Δ為試件端面的滑移速度。為了避免對(duì)端面復(fù)雜的摩擦行為進(jìn)行描述,這里采用缺口圓柱試件設(shè)計(jì)來消除端面摩擦的影響[20],如圖1所示。在壓縮過程中,試件端部不會(huì)發(fā)生變形,壓頭和試件端面之間沒有相對(duì)滑動(dòng),從而也沒有磨擦作用所帶來的能量耗散。

    圖1 缺口圓柱試件壓縮

    Fig. 1 Compression of notched cylindrical specimen

    所以,通過這樣的試件設(shè)計(jì)就可以使試件在壓縮過程中外力所做的功將完全轉(zhuǎn)化為試件的塑性應(yīng)變能:

    由于外力功完全轉(zhuǎn)化為試件的變形能,所以這里將缺口試件模型近似等效為一個(gè)無端面摩擦作用的理想圓柱試件壓縮模型。對(duì)于理想圓柱壓縮模型,試件在壓縮過程中應(yīng)力、應(yīng)變狀態(tài)均勻,所以塑性變形能P可以表示為

    (5)

    在壓縮過程中,試件在任一時(shí)刻的軸向位移()和真實(shí)應(yīng)變()可以表示為

    (7)

    (8)

    式中:U為=1時(shí)的軸向位移,為一無量綱時(shí)間參數(shù),0≤≤1;0為理想圓柱試件的初始高度;為=1時(shí)的真實(shí)應(yīng)變。所以,在U/0較小時(shí),忽略二階及其以上的高階項(xiàng),等式(8)可以表示為

    將式(6)和式(9)分別代入式(2)和(5)中,可以得到外力功和塑性變性能以及Δ的表達(dá)式為

    (10)

    (12)

    由于在壓縮過程中能量守恒,所以對(duì)于任一時(shí)刻,上述等式均成立,所以在時(shí)刻,即=1時(shí),由上述方程可以得到

    假設(shè)此時(shí)等效圓柱試件的軸向壓縮位移為,高度為,真實(shí)應(yīng)變?yōu)椋鎸?shí)壓縮流動(dòng)應(yīng)力為1,待修正的壓縮流動(dòng)應(yīng)力為2,于是可以得到

    (14)

    (16)

    式中:1為缺口試件壓縮實(shí)驗(yàn)的真實(shí)載荷,2為以流動(dòng)應(yīng)力2作為輸入進(jìn)行缺口試件壓縮實(shí)驗(yàn)?zāi)M得到的載荷。可以看出,在小應(yīng)變范圍內(nèi),上述修正關(guān)系能較好地成立;而當(dāng)試件塑性變形較大時(shí),公式(8)中忽略的高階項(xiàng)所引進(jìn)的誤差將會(huì)影響流動(dòng)應(yīng)力的修正結(jié)果。所以對(duì)于塑性變形較大的延性金屬材料而言,單次修正的效果是不夠理想的。

    1.2 迭代修正步驟

    為了減小近似修正公式在大變形條件下產(chǎn)生的誤差,得到更加準(zhǔn)確的材料應(yīng)力?應(yīng)變曲線,這里需要結(jié)合有限元分析對(duì)流動(dòng)應(yīng)力進(jìn)行迭代修正。在迭代修正過程中,和的值會(huì)逐漸減小,流動(dòng)應(yīng)力?應(yīng)變曲線也會(huì)逐漸收斂。基于上述分析,這里給出有限元迭代修正方法的具體修正過程:

    1) 通過光滑圓柱試件壓縮實(shí)驗(yàn)獲取材料的初始流動(dòng)應(yīng)力?應(yīng)變曲線;

    2) 利用設(shè)計(jì)的缺口試件進(jìn)行壓縮實(shí)驗(yàn)獲取其真實(shí)的載荷?位移曲線;

    3) 將獲取的流動(dòng)應(yīng)力?應(yīng)變曲線作為本構(gòu)輸入對(duì)缺口試件壓縮過程進(jìn)行有限元模擬,獲取材料的載 荷?位移曲線;

    4) 將步驟(3)得到的載荷?位移曲線與缺口試件壓縮實(shí)驗(yàn)得到的載荷?位移曲線進(jìn)行比較,如果載荷的相對(duì)誤差值大于規(guī)定值,則利用公式(16)對(duì)真實(shí)流動(dòng)應(yīng)力值進(jìn)行計(jì)算,得到修正后的流動(dòng)應(yīng)力?應(yīng)變曲線;

    5) 重復(fù)步驟3)和4),直至載荷的相對(duì)誤差值收斂到規(guī)定的誤差值范圍內(nèi)。

    2 鋁合金2024軸對(duì)稱圓柱試件壓縮

    2.1 圓柱試件壓縮實(shí)驗(yàn)

    本實(shí)驗(yàn)中所用材料為鋁合金2024,其主要力學(xué)性能參數(shù)如下:彈性模量71 GPa,泊松比0.3,密度2.78 g/cm3;所用試件為兩種軸對(duì)稱圓柱試件,分別為12 mm×12 mm的光滑圓柱試件以及整體尺寸為16 mm×24 mm的缺口圓柱試件。在常溫、準(zhǔn)靜態(tài)條件下對(duì)兩種試件進(jìn)行單軸壓縮加載,并記錄整個(gè)加載過程的載荷?位移曲線以及試件的變形過程。

    在光滑圓柱試件壓縮實(shí)驗(yàn)中,試件會(huì)發(fā)生整體變形,當(dāng)壓頭與試件端面之間存在摩擦作用時(shí),試件的變形會(huì)變得不均勻,在試件中部會(huì)出現(xiàn)鼓脹效應(yīng);而對(duì)于缺口圓柱,無論端面摩擦是否存在,其全部變形都集中在試件的中部區(qū)域,試件的端部不會(huì)發(fā)生變形,如圖2所示。

    圖2 不同類型圓柱試件的壓縮

    從圖2還可以看出,無論是光滑圓柱試件還是缺口圓柱試件,其最終的破壞形式均為沿45°方向的剪切破壞,并且裂紋均是從試件中部表面開始產(chǎn)生,這是由于試件端部約束(光滑圓柱試件主要為端面摩擦,缺口圓柱試件可以視為端部固約束)作用下,試件中部的鼓脹區(qū)域表面會(huì)出現(xiàn)拉應(yīng)力,從而導(dǎo)致此處的剪切應(yīng)力增大,最終使試件從這里開始發(fā)生破壞。而由于缺口試件的特殊構(gòu)型,在壓縮加載初期,試件缺口部位處于多向壓縮應(yīng)力狀態(tài),所以其壓縮過程中的鼓脹效應(yīng)要明顯弱于光滑圓柱試件,試件中部的拉應(yīng)力也相對(duì)更小,從而“減緩”了試件的破壞,可以使試件獲得更大范圍的塑性應(yīng)變,這也說明利用缺口圓柱試件壓縮實(shí)驗(yàn)對(duì)流動(dòng)應(yīng)力?應(yīng)變曲線進(jìn)行修正是可行的。

    通過壓縮實(shí)驗(yàn)得到的光滑圓柱試件的真實(shí)應(yīng)力?應(yīng)變曲線和缺口圓柱試件的載荷?位移曲線分別如圖3和圖4所示。

    2.2 缺口試件壓縮過程數(shù)值模擬

    本研究利用ABAQUS/Standard對(duì)缺口試件的壓縮過程進(jìn)行數(shù)值模擬分析,鑒于缺口試件的幾何特征,這里采用軸對(duì)稱模型,如圖5所示。試件網(wǎng)格單元類型采用4節(jié)點(diǎn)軸對(duì)稱單元,兩端壓頭采用剛性單元,壓頭與試件端面之間接觸類型為面面接觸。在壓縮過程中,模型下端軸向位移固定不變,上端以恒定速度進(jìn)行壓縮加載,數(shù)值模擬加載位移與真實(shí)實(shí)驗(yàn)位移保持一致。

    圖3 光滑圓柱試件的真實(shí)應(yīng)力?應(yīng)變曲線

    圖4 缺口圓柱試件的載荷?位移曲線

    為保證有限元模型單元尺寸不對(duì)修正結(jié)果造成影響,這里以4種不同的單元尺寸(1.6 mm、0.8 mm、0.4 mm以及0.2 mm)對(duì)試件模型進(jìn)行網(wǎng)格劃分,并以不同單元尺寸下的載荷?位移關(guān)系響應(yīng)對(duì)網(wǎng)格密度敏感性進(jìn)行了分析,如圖6所示。

    從比較結(jié)果可以看出,試件有限元模型對(duì)網(wǎng)格密度存在一定的敏感性,當(dāng)單元尺寸逐漸減小時(shí),試件的載荷?位移曲線逐漸收斂,直至單元尺寸減小到0.4 mm以內(nèi)時(shí),載荷?位移曲線基本趨于一致。所以,為了避免網(wǎng)格密度對(duì)修正結(jié)果的影響,這里取0.2 mm為試件有限元模型的單元尺寸。

    圖5 缺口試件模型

    圖6 缺口圓柱試件在不同單元尺寸下的載荷?位移曲線

    3 結(jié)果與討論

    3.1 端面摩擦對(duì)缺口試件壓縮變形行為的影響

    為了進(jìn)一步研究端面摩擦對(duì)缺口試件壓縮過程的影響,這里利用光滑圓柱試件壓縮實(shí)驗(yàn)得到的流動(dòng)應(yīng)力本構(gòu)關(guān)系對(duì)不同摩擦情況(無摩擦、=0.2、=0.4、=0.8以及粗糙表面)下缺口試件壓縮變形進(jìn)行數(shù)值模擬,得到了試件軸向壓縮載荷隨位移的變化,如圖7所示。

    從圖7可以看出,不同摩擦情況下,試件在相同壓縮位移量下所受的載荷是相同的。所以對(duì)于同一材料和幾何尺寸的缺口試件來說,摩擦約束存在與不存在時(shí)試件的載荷?位移曲線都是相同的,進(jìn)而說明壓縮載荷所做的功也是一樣的。

    數(shù)值模擬結(jié)果表明,不同摩擦條件下試件端部幾乎沒有變形,所有變形均集中在試件中部。以無摩擦和粗糙表面情況下的變形結(jié)果進(jìn)行比較知道,在相同壓縮位移下,試件徑向位移和等效塑性應(yīng)變分布一致,都分布在試件缺口區(qū)域,且最大偏差均不超過0.5%,如圖8所示,這與缺口試件壓縮實(shí)驗(yàn)現(xiàn)象也是吻合的。所以,對(duì)于缺口試件,壓頭與試件端面之間摩擦力的影響可以忽略不計(jì)。

    圖7 缺口圓柱試件在不同端面摩擦條件下的載荷?位移曲線

    3.2 流動(dòng)應(yīng)力迭代修正

    按照修正方法中的迭代修正步驟對(duì)初始流動(dòng)應(yīng)力進(jìn)行修正,得到了不同迭代次數(shù)下的鋁合金2024的流動(dòng)應(yīng)力?塑性應(yīng)變曲線以及對(duì)應(yīng)載荷?位移曲線分別如圖9和10所示。

    從圖9和10中可以看出,在迭代修正過程中,鋁合金2024的流動(dòng)應(yīng)力逐漸趨于一致,且經(jīng)過4次的迭代修正后,其流動(dòng)應(yīng)力基本保持穩(wěn)定,對(duì)應(yīng)的載荷?位移曲線與實(shí)驗(yàn)結(jié)果吻合很好。在修正過程中,載荷相對(duì)誤差的絕對(duì)值隨著迭代次數(shù)的增加而不斷減小,如圖11所示。

    從圖11比較結(jié)果可以看出,利用修正公式進(jìn)行單次修正后,相對(duì)載荷誤差值依舊較大,且隨著應(yīng)變的增加誤差逐漸增大,但經(jīng)過4次迭代修正后,載荷的最大相對(duì)誤差值已從初始狀態(tài)的11%減小到1%左右,這說明迭代過程是收斂的。為了更進(jìn)一步驗(yàn)證修正結(jié)果的準(zhǔn)確性,這里利用光滑圓柱試件壓縮實(shí)驗(yàn)得到的流動(dòng)應(yīng)力本構(gòu)關(guān)系以及經(jīng)過4次迭代修正后的流動(dòng)應(yīng)力本構(gòu)關(guān)系對(duì)缺口試件壓縮過程中缺口根部橫截面的直徑變化進(jìn)行模擬預(yù)測,并與不同載荷水平下的實(shí)驗(yàn)測量結(jié)果進(jìn)行比較,其結(jié)果如圖12所示。

    圖8 缺口圓柱試件的壓縮變形數(shù)值模擬

    從圖12對(duì)比結(jié)果可以看出,經(jīng)過迭代修正后的流動(dòng)應(yīng)力本構(gòu)關(guān)系能非常好地模擬出缺口試件壓縮實(shí)驗(yàn)過程中試件根部橫截面直徑的變化情況,而未經(jīng)修正的本構(gòu)關(guān)系預(yù)測結(jié)果則隨著壓縮載荷的增加偏差越來越大。同時(shí),從功能轉(zhuǎn)化的角度來說,本修正方法成立的前提是外力功與試件內(nèi)能之間的能量守恒,所以這里采用公式(1)對(duì)缺口試件壓縮實(shí)驗(yàn)獲取的載荷?位移曲線進(jìn)行積分計(jì)算實(shí)驗(yàn)過程中的外力功,然后在相同壓縮位移下分別采用初始流動(dòng)應(yīng)力本構(gòu)以及經(jīng)過4次迭代修正后的流動(dòng)應(yīng)力本構(gòu)對(duì)缺口試件壓縮過程進(jìn)行模擬,得到兩種本構(gòu)輸入下試件的內(nèi)能,并將其與真實(shí)實(shí)驗(yàn)過程的外力功進(jìn)行比較,其結(jié)果如圖13所示。

    圖9 不同迭代次數(shù)下的流動(dòng)應(yīng)力?塑性應(yīng)變曲線

    圖10 載荷?位移曲線模擬結(jié)果與實(shí)驗(yàn)結(jié)果的比較

    圖11 不同迭代次數(shù)下的相對(duì)載荷誤差值

    從圖13可以看出,在不同壓縮位移量下,經(jīng)過迭代修正后的流動(dòng)應(yīng)力本構(gòu)關(guān)系模擬得到的試件內(nèi)能與實(shí)驗(yàn)載荷?位移曲線計(jì)算得到的外力功始終保持一致,這也從能量守恒的角度證明了迭代修正結(jié)果的準(zhǔn)確性。

    圖12 缺口圓柱試件根部橫截面直徑模擬預(yù)測與實(shí)驗(yàn)結(jié)果的比較

    圖13 缺口圓柱試件內(nèi)能與實(shí)驗(yàn)外力功的比較

    4 結(jié)論

    1) 缺口圓柱試件的設(shè)計(jì)有效地避免了對(duì)壓縮實(shí)驗(yàn)中壓頭與試件端面之間復(fù)雜摩擦行為的分析,基于此設(shè)計(jì)提出了以能量守恒原則為基礎(chǔ)的流動(dòng)應(yīng)力近似修正公式。

    2) 借助于有限元數(shù)值模擬分析,結(jié)合流動(dòng)應(yīng)力近似修正公式可以對(duì)金屬材料壓縮流動(dòng)應(yīng)力進(jìn)行迭代修正,得到收斂后的流動(dòng)應(yīng)力值。

    3) 以鋁合金2024為例,對(duì)該修正方法進(jìn)行了驗(yàn)證,其結(jié)果證明經(jīng)該修正方法得到的流動(dòng)應(yīng)力?應(yīng)變關(guān)系可以很好地預(yù)測鋁合金2024在壓縮過程中的變形行為,并使試件內(nèi)能與實(shí)驗(yàn)得到的外力功在整個(gè)壓縮過程中保持一致。

    REFERENCES

    [1] EVANS R W, SCHARNING P J. Strain inhomogeneity in hot axisymmetric compression test[J]. Materials Science and Technology, 2002, 18(11): 1389?1398.

    [2] JOHNSON W. Experiments in plane-strain extrusion[J]. Journal of the Mechanics and Physics of Solids, 1956, 4(4): 269?282.

    [3] BAO Y B. Prediction of ductile crack formation in uncracked bodies[D]. Cambridge: Massachusetts Institute of Technology, 2003.

    [4] COOK M, LARKE E C. Resistance of Copper and Copper Alloys to Homogeneous Deformation in Compression[J]. Journal of the Institute of Metals, 1945, 71(12): 371?390.

    [5] EBRAHIMI R, NAJAFIZADEH A. A new method for evaluation of friction in bulk metal forming[J]. Journal of Materials Processing Technology, 2004, 152(20): 136?143.

    [6] LIN J F, FENG L I, ZHANG J F. A new approach to investigate real flow stress in micro-extrusion[J]. Transactions of Nonferrous Metals Society of China, 2012, 22(S2): s232?s238.

    [7] YAO Z H, MEI D Q, SHEN H, CHEN Z C. A friction evaluation method based on barrel compression test[J]. Tribology Letters, 2013, 51(3): 525?535.

    [8] 段園培, 黃仲佳, 余小魯, 邢 昌, 孫宇峰. 基于摩擦修正的TB6合金流變應(yīng)力行為研究及本構(gòu)模型建立[J]. 稀有金屬, 2014, 38(2): 202?209. DUAN Yuan-pei, HUANG Zhong-jia, YU Xiao-lu, XIN Chang, SUN Yu-feng. Flow stress behavior and constitutive model of As-Cast TB6 titanium alloy based on friction correction. Chinese Journal of Rare Metals, 2014, 38(2): 202?209.

    [9] 陳 前, 王 巖.相時(shí)效態(tài)GH4169合金的熱加工行為[J]. 中國有色金屬學(xué)報(bào), 2015, 25(10): 2727?2737. CHEN Qian, WANG Yan. Hot working behavior of delta-processed GH4169 alloy[J]. The Chinese Journal of Nonferrous Metals, 2015, 25(10): 2727?2737.

    [10] 司家勇, 廖曉航, 黃再旺, 江 亮. DP工藝inconel718合金高溫流變曲線修正及應(yīng)變本構(gòu)模型[J]. 中國有色金屬學(xué)報(bào), 2015, 25(12): 3361?3369. SI Jia-yong, LIAO Xiao-hang, HUANG Zai-wang, JIANG Liang. Revised hot derofmation behavior and strain constitutive modeling of delta-processed Inconel718 alloy[J]. The Chinese Journal of Nonferrous Metals, 2015, 25(12): 3361?3369.

    [11] WANG J P, LIN F L, HUANG B C, YUN C C. A new experimental approach to evaluate friction in ring test[J]. Journal of Materials Processing Technology, 2008, 197(1): 68?76.

    [12] 邰清安, 李曉光, 李治華, 劉 艷, 朱艷春, 曾衛(wèi)東. 應(yīng)用FEM與圓環(huán)鐓粗實(shí)驗(yàn)測定TC4鈦合金高溫變形時(shí)的摩擦因數(shù)[J]. 材料工程, 2011(6): 23?26, 31. TAI Qing-an,LI Xiao-guang,LI Zhi-hua,LIU Yan,ZHU Yan-chun,ZENG Wei-dong. On the measurement of friction coefficient of Ti-6Al-4V titanium alloy utilizing the hot compression test of ring and FE simulation[J]. Journal of Materials Engineering, 2011(6): 23?26, 31.

    [13] 陶 歡, 曾衛(wèi)東, 朱艷春, 邰清安, 李曉光, 李治華. 熱傳導(dǎo)系數(shù)對(duì)TC4鈦合金熱變形過程摩擦因數(shù)測定的影響[J]. 中國有色金屬學(xué)報(bào), 2012, 22(10): 2762?2767. TAO Huan, ZENG Wei-dong, ZHU Yan-chun, TAI Qing-an, LI Xiao-guang, LI Zhi-hua. Effect of heat transfer coefficient on measurement of friction factor in hot deformation of Ti-6Al-4V alloy[J]. The Chinese Journal of Nonferrous Metals, 2012, 22 (10): 2762?2767.

    [14] FERESHTEH-SANIEE F, NOURBAKHSH S H, PEZESHKI S M. Estimation of flow curve and friction coefficient by means of a one-step ring test using a neural network coupled with FE simulations[J]. Journal of Mechanical Science and Technology, 2012, 26(1): 153?160.

    [15] DESALEGN W, DAVIDSON M J, KHANRA A K. Friction factor evaluation using experimental and finite element methods for Al-4%Cu Preforms[J]. Journal of Materials Engineering and Performance, 2014, 23(8): 3058?3068.

    [16] BARRAU O, BOHER C, GRAS R, REZAI-ARIA F. Analysis of the friction and wear behavior of hot work tool steel for forging[J]. Wear, 2003, 255(7): 1444?1454.

    [17] HAN H. Determination of mean flow stress and friction coefficient by the modified two-specimen method in cold rolling[J]. Journal of Materials Processing Technology, 2005, 159(3): 401?408.

    [18] SHIN K, CHI S, KIM N. Prediction of flow stress of metallic material and interfacial friction condition at high temperature using inverse analysis[J]. Journal of Mechanical Science and Technology, 2010, 24(3): 639?648.

    [19] CHRISTIANSEN P, MARTINS P A F, BAY N. Friction compensation in the upsetting of cylindrical test specimens[J]. Experimental Mechanics, 2016, 56(7): 1271?1279.

    [20] BAO Y B, WIERZBICKI T. On fracture locus in the equivalent strain and stress triaxiality space[J]. International Journal of Mechanical Sciences, 2004, 46(1): 81?98.

    (編輯 何學(xué)鋒)

    Correction method of compression flow stress of metals based on energy conservation principle

    JIA Dong, HUANG Xi-cheng, MO Jun, HU Wen-jun

    (Institute of Systems Engineering, China Academy of Engineering Physics, Mianyang 621900, China)

    The influence of interfacial friction on the compressive deformation behavior of metals (especially ductile metals) was eliminated by design of axisymmetric notched cylindrical specimen. Based on the energy conservation between external work and internal energy of specimen in compression test, an approximate correction formula of compression flow stress was derived and the corresponding iterative correction procedure was proposed. For aluminum alloy 2024, both the diameter of minimum cross-section and the internal energy of specimen in compression process were predicted by the corrected flow stress-strain relationships and compared with the experimental results. The results show that the predicted diameters of minimum cross-section are agree with the experimental data and the internal energy of specimen is always consistent with the experimental external work in the whole compression process, which verify the effectiveness of the correction method.

    flow stress; interfacial friction; notched cylindrical specimen; energy conservation; iterative correction

    Project(11472257) supported by the National Natural Science Foundation of China; Project(13CXJ09) supported by Innovation and Development Foundation of Institute of Systems Engineering, China Academy of Engineering Physics

    2016-05-26; Accepted date: 2016-11-02

    JIA Dong; Tel: +86-816-2484030; E-mail: jiadong@mail.ustc.edu.cn

    10.19476/j.ysxb.1004.0609.2017.07.15

    1004-0609(2017)-07-1433-08

    TG115

    A

    國家自然科學(xué)基金資助項(xiàng)目(11472257);中國工程物理研究院總體工程研究所創(chuàng)新與發(fā)展基金資助項(xiàng)目(13CXJ09)

    2016-05-26;

    2016-11-02

    賈 東,助理研究員;電話:0816-2484030;E-mail:jiadong@mail.ustc.edu.cn

    猜你喜歡
    實(shí)驗(yàn)
    我做了一項(xiàng)小實(shí)驗(yàn)
    記住“三個(gè)字”,寫好小實(shí)驗(yàn)
    我做了一項(xiàng)小實(shí)驗(yàn)
    我做了一項(xiàng)小實(shí)驗(yàn)
    記一次有趣的實(shí)驗(yàn)
    有趣的實(shí)驗(yàn)
    微型實(shí)驗(yàn)里看“燃燒”
    做個(gè)怪怪長實(shí)驗(yàn)
    NO與NO2相互轉(zhuǎn)化實(shí)驗(yàn)的改進(jìn)
    實(shí)踐十號(hào)上的19項(xiàng)實(shí)驗(yàn)
    太空探索(2016年5期)2016-07-12 15:17:55
    成人无遮挡网站| 国国产精品蜜臀av免费| freevideosex欧美| 美女脱内裤让男人舔精品视频| 久久久久久久久久人人人人人人| 欧美 日韩 精品 国产| 又粗又硬又长又爽又黄的视频| 国产成人一区二区在线| 国产 精品1| 欧美国产精品va在线观看不卡| 男女国产视频网站| 美女国产视频在线观看| 日本色播在线视频| 一区二区三区精品91| 欧美xxⅹ黑人| 免费人成在线观看视频色| 新久久久久国产一级毛片| 丝瓜视频免费看黄片| 亚洲av综合色区一区| 亚洲欧美日韩另类电影网站| 男人爽女人下面视频在线观看| av卡一久久| 午夜福利视频精品| 黄片播放在线免费| 亚洲精品色激情综合| 久久精品国产亚洲av涩爱| 久久精品久久久久久久性| 人人澡人人妻人| 国产麻豆69| 国产有黄有色有爽视频| 欧美日韩亚洲高清精品| 春色校园在线视频观看| 久久ye,这里只有精品| 亚洲国产av新网站| 国产极品天堂在线| 国产在线免费精品| 精品一区二区三卡| av卡一久久| 欧美最新免费一区二区三区| 亚洲成人一二三区av| 日韩一本色道免费dvd| 国产 精品1| 国产精品成人在线| av在线老鸭窝| 久久青草综合色| 国产精品无大码| 尾随美女入室| 青春草国产在线视频| 菩萨蛮人人尽说江南好唐韦庄| 国产探花极品一区二区| 激情视频va一区二区三区| 亚洲精品色激情综合| 五月开心婷婷网| 午夜老司机福利剧场| 99久久综合免费| 国产片内射在线| 中国美白少妇内射xxxbb| 97精品久久久久久久久久精品| 男女国产视频网站| 国产精品秋霞免费鲁丝片| 免费高清在线观看日韩| 免费少妇av软件| 欧美xxxx性猛交bbbb| 99久久综合免费| 我的女老师完整版在线观看| av有码第一页| 国产xxxxx性猛交| 水蜜桃什么品种好| 国产一区二区在线观看av| 99精国产麻豆久久婷婷| 国产亚洲精品第一综合不卡 | 嫩草影院入口| 国产一区二区在线观看av| 中文欧美无线码| 亚洲av国产av综合av卡| 伊人亚洲综合成人网| 精品福利永久在线观看| 在线看a的网站| 国产成人精品久久久久久| 欧美精品人与动牲交sv欧美| 天天操日日干夜夜撸| 亚洲国产精品一区二区三区在线| 男女下面插进去视频免费观看 | 人人妻人人添人人爽欧美一区卜| 视频中文字幕在线观看| 亚洲欧美日韩卡通动漫| 久久97久久精品| 91精品伊人久久大香线蕉| 嫩草影院入口| 中文字幕av电影在线播放| 黑人欧美特级aaaaaa片| 亚洲色图综合在线观看| 欧美精品国产亚洲| 色吧在线观看| 精品久久久久久电影网| 国产精品一国产av| 国产精品 国内视频| 精品国产一区二区三区久久久樱花| 欧美激情国产日韩精品一区| 日本欧美视频一区| 亚洲国产精品国产精品| 妹子高潮喷水视频| 亚洲精品第二区| 日本欧美国产在线视频| 国产成人免费无遮挡视频| xxxhd国产人妻xxx| 丰满迷人的少妇在线观看| 中文精品一卡2卡3卡4更新| 少妇人妻 视频| 亚洲av中文av极速乱| 欧美xxⅹ黑人| 色网站视频免费| 岛国毛片在线播放| 久久毛片免费看一区二区三区| 巨乳人妻的诱惑在线观看| 99久久人妻综合| 成年av动漫网址| 国产精品女同一区二区软件| 女人精品久久久久毛片| 国产熟女午夜一区二区三区| 免费人妻精品一区二区三区视频| 亚洲激情五月婷婷啪啪| 亚洲国产精品一区二区三区在线| 精品国产露脸久久av麻豆| 91精品三级在线观看| 欧美老熟妇乱子伦牲交| 男女无遮挡免费网站观看| 精品人妻一区二区三区麻豆| 欧美亚洲日本最大视频资源| 久久精品夜色国产| 免费日韩欧美在线观看| 王馨瑶露胸无遮挡在线观看| 制服人妻中文乱码| 久久久久久久精品精品| 最后的刺客免费高清国语| 老司机亚洲免费影院| 97人妻天天添夜夜摸| 成年人午夜在线观看视频| 王馨瑶露胸无遮挡在线观看| 美女内射精品一级片tv| 成人漫画全彩无遮挡| 久久久精品94久久精品| 日本猛色少妇xxxxx猛交久久| 国产精品人妻久久久影院| 精品久久久久久电影网| 男女边摸边吃奶| 精品熟女少妇av免费看| av有码第一页| 久久久国产精品麻豆| 精品亚洲乱码少妇综合久久| 欧美 日韩 精品 国产| 久久精品国产a三级三级三级| 99久久人妻综合| 亚洲欧美色中文字幕在线| 国产成人精品一,二区| tube8黄色片| 久久久精品免费免费高清| 欧美激情极品国产一区二区三区 | www.色视频.com| 桃花免费在线播放| 一本—道久久a久久精品蜜桃钙片| 久久久久精品久久久久真实原创| 日韩伦理黄色片| 亚洲成国产人片在线观看| 97在线视频观看| 欧美少妇被猛烈插入视频| 人人澡人人妻人| 亚洲激情五月婷婷啪啪| 成年美女黄网站色视频大全免费| 精品国产一区二区三区久久久樱花| 成人亚洲欧美一区二区av| 夫妻性生交免费视频一级片| 欧美日韩一区二区视频在线观看视频在线| 看十八女毛片水多多多| 国产精品一区二区在线观看99| 亚洲成色77777| h视频一区二区三区| av视频免费观看在线观看| 久久av网站| 捣出白浆h1v1| 亚洲精品,欧美精品| 精品久久蜜臀av无| 咕卡用的链子| 女的被弄到高潮叫床怎么办| 久久婷婷青草| 黄色一级大片看看| 亚洲第一区二区三区不卡| 全区人妻精品视频| 国产永久视频网站| av福利片在线| 亚洲av电影在线观看一区二区三区| 亚洲av免费高清在线观看| 狠狠精品人妻久久久久久综合| 秋霞伦理黄片| 全区人妻精品视频| 一区在线观看完整版| 国产激情久久老熟女| 91在线精品国自产拍蜜月| 美女国产高潮福利片在线看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲精品自拍成人| 亚洲美女视频黄频| 国产精品欧美亚洲77777| 一级毛片电影观看| 黑人欧美特级aaaaaa片| av不卡在线播放| 亚洲国产成人一精品久久久| 韩国高清视频一区二区三区| 久久久国产欧美日韩av| 国产日韩欧美在线精品| 热99久久久久精品小说推荐| 国产熟女欧美一区二区| 如日韩欧美国产精品一区二区三区| 婷婷色综合大香蕉| 最新中文字幕久久久久| 免费久久久久久久精品成人欧美视频 | 水蜜桃什么品种好| 婷婷色综合大香蕉| 自拍欧美九色日韩亚洲蝌蚪91| av卡一久久| 欧美丝袜亚洲另类| 日韩成人伦理影院| 99国产综合亚洲精品| 男人操女人黄网站| 男女高潮啪啪啪动态图| av女优亚洲男人天堂| 久久精品国产a三级三级三级| 欧美日本中文国产一区发布| 精品熟女少妇av免费看| 两性夫妻黄色片 | 日本av免费视频播放| 黑人高潮一二区| 久久97久久精品| 九色亚洲精品在线播放| 乱码一卡2卡4卡精品| 一二三四中文在线观看免费高清| 人妻少妇偷人精品九色| 精品亚洲成国产av| 亚洲,欧美精品.| 天天操日日干夜夜撸| 国产成人免费观看mmmm| 一二三四中文在线观看免费高清| 丰满饥渴人妻一区二区三| 日韩三级伦理在线观看| 久久精品夜色国产| 视频在线观看一区二区三区| 高清毛片免费看| 亚洲色图 男人天堂 中文字幕 | 成年人免费黄色播放视频| 飞空精品影院首页| 秋霞在线观看毛片| 天堂8中文在线网| 亚洲av.av天堂| 亚洲av成人精品一二三区| 欧美激情国产日韩精品一区| 大香蕉久久网| 成年人免费黄色播放视频| 天堂中文最新版在线下载| 伦精品一区二区三区| 18禁国产床啪视频网站| 日韩一区二区视频免费看| 国产免费又黄又爽又色| 丁香六月天网| 丝袜喷水一区| 老熟女久久久| 久久久久久久国产电影| 人妻系列 视频| 亚洲av免费高清在线观看| 欧美日韩国产mv在线观看视频| 欧美精品av麻豆av| 777米奇影视久久| www.色视频.com| av福利片在线| 少妇被粗大猛烈的视频| 欧美人与性动交α欧美精品济南到 | 亚洲久久久国产精品| 亚洲精品日本国产第一区| 一本色道久久久久久精品综合| 国产精品一区二区在线观看99| 成人二区视频| 日日爽夜夜爽网站| 少妇人妻久久综合中文| 蜜桃在线观看..| 国产成人aa在线观看| 丰满乱子伦码专区| 亚洲av国产av综合av卡| kizo精华| 久久久久久久大尺度免费视频| 日韩三级伦理在线观看| 亚洲av国产av综合av卡| 亚洲性久久影院| 麻豆精品久久久久久蜜桃| 亚洲欧美清纯卡通| 国产精品久久久久成人av| 热99国产精品久久久久久7| 亚洲欧美一区二区三区国产| 亚洲国产欧美日韩在线播放| 国产亚洲av片在线观看秒播厂| 男人添女人高潮全过程视频| 色5月婷婷丁香| 人体艺术视频欧美日本| 久久免费观看电影| 观看av在线不卡| 免费不卡的大黄色大毛片视频在线观看| 在线观看国产h片| 亚洲国产精品一区三区| 久久97久久精品| 成年人午夜在线观看视频| 飞空精品影院首页| 精品一区二区三卡| 人人妻人人爽人人添夜夜欢视频| 国产探花极品一区二区| 亚洲欧美一区二区三区黑人 | 黄色配什么色好看| 国产有黄有色有爽视频| 蜜臀久久99精品久久宅男| 色吧在线观看| 亚洲情色 制服丝袜| 免费日韩欧美在线观看| 国产精品蜜桃在线观看| 欧美精品高潮呻吟av久久| 国产精品不卡视频一区二区| 亚洲国产精品一区二区三区在线| 在线观看免费日韩欧美大片| 菩萨蛮人人尽说江南好唐韦庄| 一边亲一边摸免费视频| 一区二区三区四区激情视频| 亚洲人与动物交配视频| 精品亚洲成国产av| 人妻一区二区av| 精品福利永久在线观看| 亚洲国产av新网站| 亚洲精品国产色婷婷电影| 最近手机中文字幕大全| 亚洲av综合色区一区| 人妻 亚洲 视频| 夫妻性生交免费视频一级片| 欧美亚洲 丝袜 人妻 在线| 青春草亚洲视频在线观看| 女人被躁到高潮嗷嗷叫费观| 国产又爽黄色视频| 亚洲av电影在线进入| 91午夜精品亚洲一区二区三区| 各种免费的搞黄视频| 观看av在线不卡| 久久久精品94久久精品| a级毛片黄视频| 午夜福利在线观看免费完整高清在| av免费在线看不卡| 中文字幕精品免费在线观看视频 | 香蕉丝袜av| 大片电影免费在线观看免费| 亚洲精品第二区| 99re6热这里在线精品视频| 中文字幕人妻熟女乱码| 99九九在线精品视频| 在线观看人妻少妇| 国产精品久久久久久久电影| 欧美精品国产亚洲| 亚洲国产色片| 热99久久久久精品小说推荐| 午夜av观看不卡| 午夜福利视频精品| 国产亚洲av片在线观看秒播厂| 美女大奶头黄色视频| 日日摸夜夜添夜夜爱| 97在线人人人人妻| 亚洲av日韩在线播放| 日韩av免费高清视频| 欧美3d第一页| 亚洲人成77777在线视频| 爱豆传媒免费全集在线观看| 97人妻天天添夜夜摸| 日韩在线高清观看一区二区三区| 26uuu在线亚洲综合色| 中文乱码字字幕精品一区二区三区| 久久久欧美国产精品| 国产欧美日韩综合在线一区二区| 国产片特级美女逼逼视频| 精品一品国产午夜福利视频| 日本与韩国留学比较| 一二三四中文在线观看免费高清| 欧美成人午夜免费资源| 天天操日日干夜夜撸| 丰满乱子伦码专区| 欧美日韩亚洲高清精品| 在线天堂最新版资源| 欧美精品av麻豆av| 自线自在国产av| 中文字幕人妻丝袜制服| 伦理电影免费视频| 久久久久久久大尺度免费视频| 全区人妻精品视频| 色哟哟·www| 免费看av在线观看网站| 亚洲中文av在线| 国产一区亚洲一区在线观看| 国产精品国产三级国产专区5o| 波野结衣二区三区在线| 男女高潮啪啪啪动态图| 亚洲国产成人一精品久久久| 免费观看av网站的网址| 久久久国产一区二区| 飞空精品影院首页| 免费女性裸体啪啪无遮挡网站| 大片电影免费在线观看免费| 麻豆精品久久久久久蜜桃| 90打野战视频偷拍视频| 久久人人爽人人爽人人片va| 国产亚洲一区二区精品| 日韩熟女老妇一区二区性免费视频| 久久久久视频综合| 国产精品国产三级国产专区5o| 久久女婷五月综合色啪小说| 日韩av免费高清视频| 飞空精品影院首页| 亚洲欧美一区二区三区黑人 | 欧美少妇被猛烈插入视频| 狠狠精品人妻久久久久久综合| 亚洲欧美精品自产自拍| 日韩成人av中文字幕在线观看| 十八禁高潮呻吟视频| 亚洲精品日韩在线中文字幕| 综合色丁香网| 日韩视频在线欧美| 亚洲美女黄色视频免费看| 卡戴珊不雅视频在线播放| 久久狼人影院| 欧美亚洲 丝袜 人妻 在线| av不卡在线播放| 在线精品无人区一区二区三| 精品久久蜜臀av无| 建设人人有责人人尽责人人享有的| 免费av不卡在线播放| 亚洲欧洲日产国产| 久久精品久久久久久久性| 亚洲第一区二区三区不卡| 国产国拍精品亚洲av在线观看| 亚洲综合精品二区| 国产探花极品一区二区| 精品午夜福利在线看| 久久人人97超碰香蕉20202| 国产成人精品久久久久久| 精品少妇久久久久久888优播| 亚洲第一av免费看| 欧美人与性动交α欧美精品济南到 | 全区人妻精品视频| 9色porny在线观看| 国产成人精品一,二区| 91国产中文字幕| 亚洲四区av| 亚洲中文av在线| 人人澡人人妻人| 精品一区二区三卡| 少妇被粗大猛烈的视频| 最近2019中文字幕mv第一页| 在线天堂最新版资源| 亚洲精品久久午夜乱码| av视频免费观看在线观看| 日本vs欧美在线观看视频| 久久精品久久久久久久性| 欧美xxⅹ黑人| 国产精品久久久久久av不卡| 在线观看免费视频网站a站| 亚洲av男天堂| 性色av一级| 日韩人妻精品一区2区三区| 91精品伊人久久大香线蕉| 人人妻人人爽人人添夜夜欢视频| 国产男女内射视频| 一个人免费看片子| 精品国产一区二区久久| 国产亚洲一区二区精品| 成人毛片a级毛片在线播放| 精品国产国语对白av| 亚洲精品美女久久av网站| 在线观看国产h片| 国产又爽黄色视频| 国产日韩欧美亚洲二区| 一区二区日韩欧美中文字幕 | 日日摸夜夜添夜夜爱| 亚洲,欧美精品.| 久久久久视频综合| 伊人久久国产一区二区| 国产熟女午夜一区二区三区| 日韩 亚洲 欧美在线| 免费观看av网站的网址| 美女内射精品一级片tv| 国产一区二区三区综合在线观看 | 日韩精品免费视频一区二区三区 | 亚洲四区av| 国产熟女欧美一区二区| 中文字幕av电影在线播放| 伦精品一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 国产成人午夜福利电影在线观看| 十分钟在线观看高清视频www| 婷婷成人精品国产| 亚洲av国产av综合av卡| 国产 精品1| 最近最新中文字幕大全免费视频 | 国产69精品久久久久777片| 狠狠婷婷综合久久久久久88av| 在线观看美女被高潮喷水网站| 97超碰精品成人国产| 亚洲内射少妇av| 午夜影院在线不卡| 久久久久久伊人网av| 高清毛片免费看| 中文字幕人妻熟女乱码| 丝袜在线中文字幕| 欧美日韩成人在线一区二区| 热99国产精品久久久久久7| 美女脱内裤让男人舔精品视频| 国产成人a∨麻豆精品| 看免费成人av毛片| 国产精品女同一区二区软件| 免费人成在线观看视频色| 爱豆传媒免费全集在线观看| 伊人久久国产一区二区| 最新中文字幕久久久久| 中文字幕亚洲精品专区| 国产精品一国产av| 少妇熟女欧美另类| 亚洲精品日韩在线中文字幕| 一二三四在线观看免费中文在 | 中国三级夫妇交换| 高清在线视频一区二区三区| 日韩制服骚丝袜av| 18禁动态无遮挡网站| 69精品国产乱码久久久| 亚洲精品456在线播放app| 亚洲一区二区三区欧美精品| 曰老女人黄片| 一级片免费观看大全| 亚洲av免费高清在线观看| 在线天堂最新版资源| 中文字幕精品免费在线观看视频 | 久久韩国三级中文字幕| 久久精品人人爽人人爽视色| 久久 成人 亚洲| 伦理电影大哥的女人| 校园人妻丝袜中文字幕| 久久久久国产网址| 国精品久久久久久国模美| 九草在线视频观看| 亚洲 欧美一区二区三区| 国产一区二区三区综合在线观看 | 伦精品一区二区三区| 青青草视频在线视频观看| 宅男免费午夜| 青青草视频在线视频观看| 久久久欧美国产精品| 国产高清国产精品国产三级| 18+在线观看网站| 精品一区在线观看国产| 日韩在线高清观看一区二区三区| 国产日韩一区二区三区精品不卡| 国产av一区二区精品久久| 九色成人免费人妻av| 欧美变态另类bdsm刘玥| 精品一区在线观看国产| 午夜免费男女啪啪视频观看| 欧美精品av麻豆av| 国产午夜精品一二区理论片| 一边亲一边摸免费视频| 久久久精品免费免费高清| 久久ye,这里只有精品| 久久亚洲国产成人精品v| 成人影院久久| 久久av网站| 人人澡人人妻人| 国产欧美亚洲国产| 亚洲精品av麻豆狂野| 久久久国产精品麻豆| 大片免费播放器 马上看| 人人妻人人澡人人爽人人夜夜| 最近中文字幕高清免费大全6| 国产成人精品福利久久| 中文字幕av电影在线播放| 亚洲成色77777| av国产久精品久网站免费入址| 国产精品久久久久久久电影| 高清欧美精品videossex| 亚洲人成77777在线视频| 国产免费一区二区三区四区乱码| 国产有黄有色有爽视频| 国产精品嫩草影院av在线观看| 最近最新中文字幕大全免费视频 | 少妇被粗大猛烈的视频| 99热全是精品| 国产精品无大码| 飞空精品影院首页| 亚洲国产av新网站| 亚洲国产毛片av蜜桃av| 久久久久视频综合| 国产av一区二区精品久久| 精品第一国产精品| 黄色配什么色好看| 亚洲欧美清纯卡通| 男的添女的下面高潮视频| 少妇 在线观看| 亚洲天堂av无毛| 中文字幕制服av| 国产成人免费观看mmmm| 欧美精品高潮呻吟av久久| 欧美精品国产亚洲| 97超碰精品成人国产| 日日爽夜夜爽网站| 在线天堂最新版资源| 欧美日韩av久久| 久久99精品国语久久久| 国产一级毛片在线| 日本wwww免费看| 巨乳人妻的诱惑在线观看|