亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        鋁合金厚板殘余應(yīng)力測試方法有效性分析

        2010-11-01 14:08:30吳運(yùn)新
        材料工程 2010年1期
        關(guān)鍵詞:心部柔度厚板

        龔 海,吳運(yùn)新,廖 凱

        (中南大學(xué)機(jī)電工程學(xué)院,長沙410083)

        鋁合金厚板殘余應(yīng)力測試方法有效性分析

        龔 海,吳運(yùn)新,廖 凱

        (中南大學(xué)機(jī)電工程學(xué)院,長沙410083)

        通過分析試件應(yīng)力分布的均勻性,對裂紋柔度法和層削法測試結(jié)果的有效性進(jìn)行了評價(jià)。結(jié)果表明:深層應(yīng)力測試方法能有效描述應(yīng)力分布均勻的試件內(nèi)應(yīng)力的真實(shí)分布;對于應(yīng)力分布均勻性差的試件,測試結(jié)果反映的是綜合應(yīng)力,不能有效描述內(nèi)應(yīng)力的真實(shí)分布,內(nèi)應(yīng)力分布的準(zhǔn)確描述需要結(jié)合有限元仿真、表面測試和深層應(yīng)力測試進(jìn)行綜合評價(jià)。

        殘余應(yīng)力;鋁合金;裂紋柔度法;層削法;X射線衍射;有限元法

        Abstract:The validity of both crack compliance method and layer removal method measurement results was evaluated by analyzing the stress distribution uniformity of the test coupons.The results show that the two methods can describe the real stress distribution effectively for test coupon with mean stress distribution;while for test coupon with unmean stress distribution,the measurement results are synthetic stresses,and can not describe a real stress distribution.In order to describe the real residual stress distribution of test coupon accurately,it is necessary to combine the use of FEM,surface stress measurement and deep floor stress measurement,and to give an integrate evaluation.

        Key words:residual stress;aluminum alloy;crack compliance method;layer removal method;X-ray diffraction;finite element method

        提高航空結(jié)構(gòu)件殘余應(yīng)力的測試精度,意味著可提高其安全可靠性,減輕質(zhì)量及降低成本[1]。到目前為止,厚板內(nèi)部殘余應(yīng)力的測量和評估還沒有一個(gè)成熟、有效的方法。中子衍射法作為無損檢測技術(shù)能夠相對有效地檢測厚板內(nèi)應(yīng)力,但目前國內(nèi)尚不具備相應(yīng)的實(shí)驗(yàn)條件,因此對于材料內(nèi)部宏觀殘余應(yīng)力的檢測,通常采用有損的力學(xué)測量方法,其中最常用的方法是層削法和裂紋柔度法[1-3]。

        PRIME等[2]最先采用裂紋柔度法對美國鋁業(yè)公司生產(chǎn)的7050-T74軋制板進(jìn)行了殘余應(yīng)力測試,獲得了不同狀態(tài)的鋁厚板內(nèi)部殘余應(yīng)力分布,之后王秋成、王樹宏等[1,3]分別采用裂紋柔度法和層削法對鋁厚板內(nèi)部殘余應(yīng)力分布進(jìn)行了測試。然而,上述研究僅僅將方法應(yīng)用于鋁厚板的殘余應(yīng)力測試,對于測試結(jié)果的有效性及測試方法的適用性沒有進(jìn)行深入的分析。本工作在實(shí)驗(yàn)研究的基礎(chǔ)上,采用有限元法對淬火鋁厚板內(nèi)應(yīng)力場的分布狀況進(jìn)行了分析,深入探討了應(yīng)力分布均勻性對測試結(jié)果有效性的影響。

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

        1.1 測試方法原理

        裂紋柔度法和層削法均假設(shè)板為均勻連續(xù)材料,殘余應(yīng)力沿板的長度、寬度方向的分布基本不變,只沿厚度方向(z)存在較大的應(yīng)力梯度。

        裂紋柔度法的測定原理是在被測物體表面引入一條深度逐漸增加的裂紋來釋放殘余應(yīng)力,通過測定工件表面特定點(diǎn)的應(yīng)變隨裂紋深度的變化來計(jì)算殘余應(yīng)力[2],如圖1所示。方法中進(jìn)一步假設(shè),工件內(nèi)部應(yīng)力梯度可以用多項(xiàng)式連續(xù)函數(shù)描述為

        式中:Ai為待定系數(shù);Pi(z)為插值函數(shù),Pi(z)的選擇必須滿足板內(nèi)殘余應(yīng)力的力和力矩自平衡條件,本工作中采用2~9階Legendre級數(shù)。為了確定待定系數(shù)Ai,首先要計(jì)算圖1中所示應(yīng)變片位置上的應(yīng)變隨裂紋深度aj變化的響應(yīng)值Cij(稱為柔度函數(shù))。根據(jù)疊加原理,該應(yīng)變可以表示為

        圖1 裂紋柔度法的應(yīng)變測量Fig.1 Strain measurement of crack compliance method

        在計(jì)算Ai時(shí),為了減少計(jì)算誤差,根據(jù)柔度函數(shù)計(jì)算得到的應(yīng)變值與實(shí)驗(yàn)得到的m次應(yīng)變讀數(shù),采用最小二乘法擬合,得到

        通過求解上述的n+1個(gè)方程,可以得到方程(1)的系數(shù)Ai,即

        層削法測試原理如圖2所示,假設(shè)試件處于內(nèi)應(yīng)力平衡的初始狀態(tài),當(dāng)圖中標(biāo)識為A的部分被剝除后,引起試件剩余部分B的應(yīng)力重新分布而產(chǎn)生變形,通過測定試件底部應(yīng)變的大小,即可計(jì)算出殘余應(yīng)力的分布,其常用公式為[4]

        式中:ε為測量得到的試件下表面應(yīng)變對深度a的函數(shù)ε(a)。將逐層銑削中獲得的ε-a數(shù)據(jù)擬合成函數(shù)關(guān)系ε(a)或ε(z),代入式(5)中,就能求得應(yīng)力σx,y(z)。

        1.2 實(shí)驗(yàn)條件

        研究對象為國產(chǎn)7075鋁合金軋制板,厚板尺寸為1200mm×220mm×50mm,固溶處理后,分別采用浸沒淬火和噴淋淬火,淬火溫度為475℃,淬火介質(zhì)為室溫自來水。浸沒淬火采用直接將鋁板投入水池的方法;噴淋淬火采取鋁板側(cè)立放置,從兩邊噴水的方法,在水池兩側(cè)壁各安置一組噴嘴,每組上下各一排,每排9個(gè)呈對稱分布,噴水壓力為 0.34MPa,水流量為136L/min,噴射面形狀為圓錐形。

        采用裂紋柔度法和層削法分別對兩種淬火板進(jìn)行內(nèi)部殘余應(yīng)力測試,試樣從淬火板的長度和寬度中間部位截取,裂紋柔度法和層削法試件尺寸分別為120mm×120mm×50mm和 160mm×160mm× 50mm。層削法實(shí)驗(yàn)在L EADWELL V-60A數(shù)控加工中心進(jìn)行,適當(dāng)調(diào)整銑削參數(shù)以減小加工應(yīng)力:銑刀轉(zhuǎn)速1000r/min,進(jìn)給速度50mm/min,每剝除2mm后記錄一次應(yīng)變數(shù)據(jù),每層分兩刀進(jìn)行,第一刀剝除1.5mm,第二刀剝除0.5mm。裂紋柔度法實(shí)驗(yàn)采用DK77線切割機(jī)來引入裂紋,鉬絲直徑φ0.18mm,加工后切縫尺寸約為0.5mm,進(jìn)給速度1mm/min,切割深度每增加1mm記錄一次應(yīng)變數(shù)據(jù)。選用BX120-5AA應(yīng)變片,其柵長×柵寬尺寸為5mm×3mm,用 YE2533靜態(tài)應(yīng)變儀記錄數(shù)據(jù)。

        圖2 層削法原理圖Fig.2 Principal diagram of layer removal method

        1.3 測試結(jié)果

        測試結(jié)果如圖3所示,圖3(a),(b)分別對應(yīng)浸沒淬火板和噴淋淬火板。可以看出,淬火應(yīng)力沿板厚度方向呈外壓內(nèi)拉的拱形分布。對比圖3(a),(b)可以看出,由于淬火冷卻速率的差異,兩種淬火板內(nèi)的應(yīng)力存在較大的差異,浸沒淬火板內(nèi)應(yīng)力明顯大于噴淋淬火板。

        對于同種淬火板,總體上兩種測試結(jié)果一致,顯示了相同的應(yīng)力分布規(guī)律,但在板厚度的表層和心部存在不同程度的差異:對于浸沒淬火板,層削法計(jì)算結(jié)果為表層-216MPa,心部143MPa,明顯高于裂紋柔度法的表層-130MPa,心部90MPa;而噴淋淬火板,兩種測試結(jié)果一致性較好,最大差異出現(xiàn)在厚板表層,層削法和裂紋柔度法結(jié)果分別為-77,-52MPa。

        圖3 淬火鋁厚板內(nèi)部軋向殘余應(yīng)力 (a)浸沒淬火板;(b)噴淋淬火板Fig.3 Rolling direction residual stresses in quenched plates (a)immerging quenched plate;(b)spraying quenched plate

        2 有效性分析

        2.1 有限元仿真

        運(yùn)用有限元軟件MSC.Marc,對鋁厚板淬火過程進(jìn)行數(shù)值仿真,模擬實(shí)驗(yàn)中兩種不同的淬火方式:浸沒淬火和噴淋淬火,通過熱力準(zhǔn)耦合計(jì)算淬火殘余應(yīng)力。模型中假設(shè)試樣為連續(xù)均勻、各向同性固體材料,初始溫度均勻且無內(nèi)應(yīng)力場,淬火過程為熱應(yīng)力分析,材料的彈性模量取E=71.7GPa,泊松比取μ=0.33,其他熱物理參數(shù)和力學(xué)性能參數(shù)(如比熱、熱膨脹系數(shù)和表面換熱系數(shù)等)均隨溫度變化[5,6],通過實(shí)測計(jì)算得到。

        考慮到厚板具有幾何對稱性和邊界條件對稱性,模型取實(shí)際板尺寸的 1/8(即 600mm×110mm× 25mm)如圖4所示,坐標(biāo)軸x,y,z分別表示厚板軋制(長度)方向、寬度方向和厚度方向,坐標(biāo)原點(diǎn)o設(shè)在實(shí)際厚板的正心部,順著坐標(biāo)軸的正方向表示從心部往板表面延伸。

        圖4 厚板1/8模型及分析截面示意圖Fig.4 Sketch-map of 1/8 plate model and analyzed sections

        為了研究鋁厚板淬火后不同厚度上的應(yīng)力分布狀態(tài),同時(shí)考慮到實(shí)驗(yàn)試件從厚板的長度和寬度中心截取,取模型心部長度為200mm的一段上的四個(gè)截面進(jìn)行分析,分別以對角線字母標(biāo)記為截面ao,ae,bf和cg,其中ao為垂直截面,ae,bf,cg為水平截面,如圖4中虛線所示。浸沒淬火板和噴淋淬火板四個(gè)截面上的軋向應(yīng)力分布分別如圖5,6所示。

        從圖5(a),6(a)可以看出,兩種淬火板內(nèi)的軋向殘余應(yīng)力隨厚度的變化都表現(xiàn)為明顯的外壓內(nèi)拉分布,對應(yīng)不同的寬度位置,這種趨勢都一致;而軋向殘余應(yīng)力沿板寬度方向的變化存在較大的差異,同一厚度上噴淋淬火板內(nèi)應(yīng)力隨寬度變化較小,浸沒淬火板內(nèi)應(yīng)力隨寬度變化較大。

        從圖5(c),(d),(e)和圖6(c),(d),(e)可以看出,兩種淬火板內(nèi)的軋向殘余應(yīng)力沿板長度方向的變化都很小,對應(yīng)不同的寬度、厚度位置這種趨勢一致。而軋向殘余應(yīng)力沿板寬度方向的變化非常明顯:對于浸沒淬火板,厚板心部cg截面、表面ae截面上軋向殘余應(yīng)力沿寬度變化很大,差異分別達(dá)到40MPa和30MPa,而中間截面bf上殘余應(yīng)力沿寬度變化相對較小,只在寬度的邊緣出現(xiàn)了一定的起伏;對于噴淋淬火板,總體上軋向殘余應(yīng)力沿寬度變化都比浸沒淬火板小,心部cg截面、表面ae截面上的變化略大,而中間截面bf上殘余應(yīng)力沿寬度變化很小。

        可見,兩種淬火板中,軋向殘余應(yīng)力沿板長度方向的分布都是比較均勻的;沿板寬度方向的分布,整體上噴淋淬火板比較均勻,浸沒淬火板則均勻性較差;從全厚度來看,厚板心部和表層應(yīng)力分布均勻性較差,厚度中心面兩側(cè)的對稱區(qū)域應(yīng)力分布均勻性較好。

        考慮到厚板淬火的邊緣效應(yīng)[7],板的長度和寬度邊緣部分應(yīng)力的分布比較復(fù)雜,應(yīng)力測試中也是從板的長度、寬度中間部位截取試件,因此可以認(rèn)為,試樣上板厚的表層和心部區(qū)域軋向應(yīng)力沿寬度分布的均勻性較差,而表層和心部之間區(qū)域則比較均勻。

        2.2 表面應(yīng)力測試

        對截取上述深層應(yīng)力測試試件的鋁板母板,采用Proto X射線衍射儀進(jìn)行表面應(yīng)力測試,在每塊板的寬度中心測試3點(diǎn),寬度邊緣測試2點(diǎn),測試結(jié)果如表1所示。可以看出,浸沒淬火板寬度中心位置和邊緣位置的表面軋向應(yīng)力表現(xiàn)為中心小邊緣大,差值達(dá)到60MPa,可見其表面應(yīng)力分布很不均勻;而噴淋淬火板寬度中心位置和邊緣位置的表面軋向應(yīng)力值則相差很小,可以認(rèn)為其表面應(yīng)力分布比較均勻。

        表1 表面應(yīng)力X射線衍射測試結(jié)果Table 1 Surface residual stresses by X-ray diffraction

        2.3 實(shí)驗(yàn)結(jié)果討論

        裂紋柔度法和層削法測試結(jié)果的差異,需要從方法原理上進(jìn)行分析。兩種方法均假設(shè)板內(nèi)殘余應(yīng)力沿長度和寬度方向基本不變,只隨板厚度方向變化。層削法實(shí)驗(yàn)中每銑削一層(2mm左右)記錄一次應(yīng)變數(shù)據(jù),這個(gè)應(yīng)變是整個(gè)剝除層內(nèi)應(yīng)力的綜合結(jié)果,因而層削法求解出的是整個(gè)剝除層內(nèi)的綜合應(yīng)力;裂紋柔度法則是采取線切割的方法在試件上引入裂紋,切深每增加1mm(或更小)后記錄一次應(yīng)變數(shù)據(jù),這個(gè)應(yīng)變是作用在整個(gè)裂紋面上的應(yīng)力的綜合結(jié)果,因而裂紋柔度法求解出的是作用在整個(gè)裂紋面上的綜合應(yīng)力。相比層削法,裂紋柔度法測得的仍然是厚板內(nèi)部的局部應(yīng)力,因此,它受應(yīng)力分布狀態(tài)的影響要小得多。

        理論上,一塊應(yīng)力分布均勻的試樣,兩種方法的測試結(jié)果會很接近。上述測試結(jié)果的差異,應(yīng)當(dāng)從以下方面來分析:①測試方法固有的計(jì)算誤差;②從大板上鋸切測試用小試塊導(dǎo)致應(yīng)力釋放引起的誤差;③實(shí)驗(yàn)過程中的加工及讀數(shù)誤差;④厚板本身內(nèi)應(yīng)力分布不均勻引起的誤差。

        首先,采用層削法和裂紋柔度法對無應(yīng)力板進(jìn)行測試(尺寸分別為 160mm×160mm×50mm和120mm×120mm×50mm,鋸切好后固溶處理再隨爐冷卻,板內(nèi)應(yīng)力可視為0),層削法的結(jié)果在±15MPa以內(nèi),而裂紋柔度法的結(jié)果則在±3MPa以內(nèi)。銑削過程中,加工應(yīng)力會對層削法的測試帶來一定的影響,而裂紋柔度法中線切割帶來的加工應(yīng)力可以忽略不計(jì)??梢?在存在加工應(yīng)力和計(jì)算誤差的情況下,兩種方法的誤差較小,方法本身計(jì)算誤差和實(shí)驗(yàn)誤差不是引起兩種測試結(jié)果差異(浸沒淬火板表層相差86MPa,心部相差53MPa)的主要原因。其次,鋸切試塊會導(dǎo)致應(yīng)力釋放,改變厚板內(nèi)部初始應(yīng)力狀態(tài),引起一定的測試誤差,參閱文獻(xiàn)[2]可知,鋸切過程中應(yīng)力釋放引起的誤差小于1%。可見,試件制備過程中的應(yīng)力損失也不是導(dǎo)致測試結(jié)果差異的主要原因,兩種測試結(jié)果的差異應(yīng)當(dāng)歸因于淬火板內(nèi)應(yīng)力分布的不均勻性。

        仿真結(jié)果中,在ao垂直截面的表層和心部位置上,軋向應(yīng)力沿寬度分布不均勻,這就必然導(dǎo)致裂紋柔度法測試結(jié)果的誤差;在表面ae水平截面、心部cg水平截面上,軋向應(yīng)力沿寬度分布也不均勻,且這種不均勻性在整個(gè)長度上都存在,這就必然導(dǎo)致層削法測試結(jié)果的誤差。單個(gè)垂直截面上應(yīng)力沿寬度分布不均勻,與試件整個(gè)長度上應(yīng)力沿寬度分布不均勻,其綜合結(jié)果不同,這將導(dǎo)致層削法測試結(jié)果明顯不同于裂紋柔度法測試結(jié)果??梢?正是淬火鋁厚板內(nèi)部殘余應(yīng)力分布的不均勻性,造成了力學(xué)測試方法的測量誤差。

        表面應(yīng)力測試結(jié)果表明,淬火鋁厚板表面應(yīng)力分布的不均勻性真實(shí)存在,浸沒淬火板表層應(yīng)力沿寬度方向中心小邊緣大、并在整個(gè)長度上保持此趨勢的分布特點(diǎn),使得整個(gè)表層的綜合應(yīng)力大于表層裂紋面上的綜合應(yīng)力,即層削法測試結(jié)果大于裂紋柔度法測試結(jié)果。

        綜上所述,從國內(nèi)現(xiàn)有的技術(shù)條件和測試手段來看,裂紋柔度法和層削法是深層殘余應(yīng)力測試的最佳選擇,但二者都建立在平面應(yīng)力假設(shè)的基礎(chǔ)上,因此,在應(yīng)用中需要對測試結(jié)果有效性進(jìn)行評價(jià)。對于應(yīng)力分布均勻的工件,滿足測試方法的假設(shè),其測試結(jié)果是有效的;但對應(yīng)力分布均勻性較差的工件,如上述淬火鋁厚板,測試結(jié)果反映的是一個(gè)綜合力,不能準(zhǔn)確描述厚板的應(yīng)力狀態(tài),需要借助準(zhǔn)確的數(shù)值模擬和精度高的表面應(yīng)力測試來綜合評定試件內(nèi)部的應(yīng)力分布。

        3 結(jié)論

        (1)鋁厚板長寬尺寸越大、淬火冷卻速率越慢,內(nèi)部淬火應(yīng)力分布越均勻;厚板全厚度上,表層和心部的應(yīng)力分布均勻性較差,表層和心部之間區(qū)域的應(yīng)力分布均勻性較好。

        (2)裂紋柔度法和層削法都是建立在平面應(yīng)力假設(shè)的基礎(chǔ)上,測得的是應(yīng)力釋放區(qū)域的綜合應(yīng)力,因此在實(shí)際運(yùn)用中需要分析測試結(jié)果的有效性。

        (3)對于應(yīng)力分布均勻的工件,裂紋柔度法和層削法的測試結(jié)果是有效的,噴淋淬火鋁厚板內(nèi)部殘余應(yīng)力約為表層-60MPa,心部30MPa;對于應(yīng)力分布不均勻的工件,裂紋柔度法和層削法測試結(jié)果不能有效描述厚板內(nèi)應(yīng)力的真實(shí)分布,必須借助其他研究手段作出綜合評價(jià)。

        [1]王秋成,柯映林,章巧芳.7075鋁合金板材殘余應(yīng)力深度梯度的評估[J].航空學(xué)報(bào),2003,24(4):336-338.

        [2]PRIME MICHAEL B,HILL MICHAEL R.Residual stress,stress relief,and inhomogeneity in aluminum plate[J].Scripta Materialia,2002,46(1):77-82.

        [3]王樹宏,左敦穩(wěn),潤長生,等.L Y12,B95和7050鋁合金預(yù)拉伸厚板內(nèi)部殘余應(yīng)力分布特征評估與分析[J].材料工程,2004,(10): 32-35.

        [4]米谷茂.殘余應(yīng)力的產(chǎn)生和對策[M].北京:機(jī)械工業(yè)出版社, 1983.60-67.

        [5]姚燦陽.7050鋁合金厚板淬火溫度場及內(nèi)應(yīng)力場的數(shù)值模擬研究[D].長沙:中南大學(xué)機(jī)電工程學(xué)院,2007.

        [6]廖凱,吳運(yùn)新,張園園.鋁合金厚板淬火表面換熱系數(shù)的離散解析求法[J].材料熱處理學(xué)報(bào),2009,39(1):192-196.

        [7]胡少虬,曾蘇民.無相變合金淬火熱應(yīng)力演變機(jī)理的理論模型-(1)角端和邊緣淬火熱應(yīng)力模型與淬火角端效應(yīng)[J].稀有金屬材料與工程,2006,35(4):538-541.

        Analysis on Validity of Residual Stress Measurement Methods for Aluminum Alloy Thick-plate

        GONG Hai,WU Yun-xin,LIAO Kai
        (College of Mechanical and Electrical Engineering, Central South University,Changsha 410083,China)

        TG146.2

        A

        1001-4381(2010)01-0042-05

        國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃資助項(xiàng)目(2005CB623708)

        2009-01-12;

        2009-11-17

        龔海(1982—),男,博士研究生,從事鋁合金厚板殘余應(yīng)力檢測和消減的研究,聯(lián)系地址:湖南長沙中南大學(xué)機(jī)電工程學(xué)院(410083),E-mail:gonghai88@163.com

        猜你喜歡
        心部柔度厚板
        軋制壓縮比對特厚板心部組織影響分析
        鋼結(jié)構(gòu)厚板及中厚板的焊接
        重言“憂”義字略辨
        薄壁齒套的心部硬度控制
        基于模態(tài)柔度矩陣識別結(jié)構(gòu)損傷方法研究
        鎳基合金UNS N08810厚板埋弧焊工藝
        基于柔度比優(yōu)化設(shè)計(jì)杠桿式柔性鉸鏈放大機(jī)構(gòu)
        厚板鋁合金MIG+TIG雙槍焊接工藝
        焊接(2016年9期)2016-02-27 13:05:25
        厚板、單向板和雙向板的工程量統(tǒng)計(jì)
        基于模態(tài)柔度矩陣的結(jié)構(gòu)損傷識別
        天堂69亚洲精品中文字幕| 国产无套中出学生姝| 国产夫妇肉麻对白| 欧美性开放bbw| 性做久久久久久久| 久久偷拍国内亚洲青青草| 国产精品一区二区三区在线免费| 最新中文字幕av无码不卡| 亚洲av成人精品日韩一区| 久久久久久免费播放一级毛片| 国产中文字幕亚洲国产| 大地资源高清在线视频播放| 久久久噜噜噜www成人网| 中文字幕在线观看国产双飞高清| 91久久精品一二三区色| 久久99精品久久久大学生| 无码av免费一区二区三区试看| 亚洲伊人久久成人综合网| 日韩中文字幕一区在线| 天天综合网网欲色| 在教室伦流澡到高潮hgl视频| 亚洲先锋影院一区二区| 亚洲一区二区三区色偷偷| 夹得好湿真拔不出来了动态图 | 一区二区三区蜜桃av| 性无码专区无码| 国产第一草草影院| 精品人妻一区二区三区蜜臀在线| 日本一二三区在线观看视频| 婷婷亚洲久悠悠色悠在线播放| 人妻丰满熟妇AV无码片| 国产麻豆国精精品久久毛片| 亚洲自偷精品视频自拍| 国产va免费精品观看| 亚洲AV秘 无码一区二区在线| 99久久婷婷国产一区| 色噜噜av亚洲色一区二区| 无码片久久久天堂中文字幕| 超碰青青草手机在线免费观看 | 午夜免费观看国产视频| 色偷偷偷久久伊人大杳蕉|