91麻豆成人久久精品二区三区_亚洲成人av一区二区_国产不卡人人_一区二区欧美日韩视频

  1. 首頁
  2. 噴丸處理

激光噴丸參數對 7075 鋁合金殘余應力分布影響規律的數值分析

為了研究不同激光噴丸參數對殘余應力分布的影響規律, 較快獲取期望的殘余應力分布, 通過利用有限元軟件 ABAQUS 對激光噴丸過程進行數值模擬, 完成了 7075 航空鋁合金靶材在不同沖擊波峰值壓力、 沖擊次數、 光斑大小及搭接率 等工藝參數下的數值模擬。 模擬結果表明: 隨著峰值壓力的增加, 材料表面殘余壓應力也增加, 當峰值壓力增大到一定值時, 材料表面殘余壓應力達到最大; 當峰值壓力一定時, 材料表面及內部殘余壓應力隨光斑大小及沖擊次數的增加而增加; 隨著 光斑搭接率的增加, 材料表面殘余壓應力增大, 且隨著搭接率的增加, 表面殘余壓應力的波動逐漸減小。

飛機服役一定年限后, 其結構件存在大量的腐 蝕損傷, 時刻威脅著飛機的飛行安全 。 因此, 解 決材料表面腐蝕損傷問題、 提高金屬材料耐應力腐 蝕性能顯得尤為迫切。 目前, 國內外通常采用表面 強化技術來改善材料的耐應力腐蝕性能。 與其他 表面強化技術相比, 激光噴丸具有靈活可控、 應變 率高、 效果顯著等特點, 因而在航空航天領域具有 較為廣闊的應用前景 。 掌握各種不同激光噴丸工藝參數對材料表面殘余應力分布的影響規律, 將有 助于更加快捷地獲取所期望的殘余應力分布, 為實 際飛機維修工作提供一定的理論指導。

激光噴丸參數對殘余應力分布影響規律方面, 周建忠 、 胡永祥 等分別以 2024 – T3 鋁合金和 35CD4 50HRC 鋼為研究對象, 通過有限元模擬發 現, 材料表面殘余壓應力隨著沖擊次數的增加而增 加, 壓應力影響層深度也隨之增加。 陳瑞芳等[6] 借 助 ANSYS 軟件對激光噴丸 40Cr 鋼試樣分析發現,材料表面殘余壓應力隨沖擊波峰值壓力的增加而增 加, 當峰值壓力增加至一定值時, 材料表面殘余壓 應力達到最大并不再增加。 黃舒等[7] 利用 ABAQUS軟件對 6061 – T6 鋁合金試樣模擬分析發現, 試樣殘余應力場隨光斑尺寸的增加而增加。 錢紹祥等通過激光沖擊 AISI202 不銹鋼焊接接頭試驗發現, 提高光斑搭接率可以優化激光噴丸后焊接接頭的應 力分布, 提高材料表面壓應力水平。 掌握不同激光 噴丸參數對材料表面殘余應力分布的影響規律有助 于激光噴丸在實際操作過程中更好地運用。

本文以飛機常用 7075 航空鋁合金靶材為研究對 象, 利用有限元軟件 ABAQUS 對激光噴丸過程進行 數值模擬, 分析 7075 鋁合金靶材在不同峰值壓力、 沖擊次數、 光斑大小及搭接率等條件下材料表面殘 余應力分布情況, 獲取各不同工藝參數對材料表面 及內部殘余應力分布的影響規律, 為實際的飛機維 修工作奠定一定的理論基礎。

1  沖擊波壓力模型

由于激光噴丸過程中, 沖擊波對材料的作用過程極為復雜, 因而難以運用有限元軟件 ABAQUS 對 其進行精確的數值模擬。 目前, 國內外學者一般選 擇將沖擊波對材料作用等效為隨時間變化的壓力對 材料的作用。 1990 年, Fabbro R 等通過對毫米尺度 約束模型下沖擊波的半理論研究, 提出沖擊波峰值 壓力與激光功率密度的關系如下 :

公式
公式

式中: Pmax 為峰值壓力; α 為等離子體的比熱比; Z 為固體靶及約束層對沖擊波的聲阻抗; I0 為激光功 率密度。

Zhang W W 等沿用 Fabbro R 等的理論, 但其認 為激光能量隨其與光斑中心距離呈準高斯分布趨勢, 并得到沖擊波峰值壓力隨空間分布關系如下 :

公式
公式

式中: r 為某點與光斑中心距離; R 為光斑半徑。 本文采用式 (2) 對沖擊波峰值壓力隨空間分 布關系進行計算, 采用的激光器基本參數為光斑能 量2J,脈寬≤10ns,焦點直徑Ф2mm,通過計算 可得, 光斑中心處沖擊波峰值壓力最大可達到3. 4 GPa。

Peyre P 等通過試驗研究發現, 沖擊波壓力作用時間為激光脈寬的 3 ~ 4 倍。 本文所采用激光 器脈寬為 10 ns, 故而沖擊波壓力作用時間可取30ns模擬采用的沖擊波壓力隨時間變化曲線如圖1所示。

沖擊壓力時間曲線
沖擊壓力時間曲線

2  材料的本構模型

激光噴丸是利用短脈沖?(幾十?ns)、?高峰值功 率密度?( > 109?W·cm?– 2?)?的激光束輻射金屬材料 表面,?其應變速率極高?( >10的6次方?/S)。?激光噴丸過 程中,?由于金屬材料表面吸收層的隔熱作用,?金屬 材料的塑性變形可看做冷變形,?材料表面只受到沖 壓力的作用。?材料屈服強度和流體應力隨著應變率 的提高而提高,?而硬化率則會隨之降低,?因此,?可 近似采用理想彈塑性材料模型。?在一維應變條件下,?材料的初始屈服極限應力可定義為?Hugoniot?彈性極 限?(σHEL)。?當材料中應力波壓力小于?σHEL?時,?材 料處于彈性極限范圍內,?只會發生彈性變形;?當壓 力大于?σHEL?時,?材料將產生無法恢復的塑性變形;?當材料中應力波壓力達到?σHEL?時,?材料產生屈服。?材料的動態屈服強度?σr?計算式如下?:

公式
公式

式中: v 為材料的泊松比。 7075 航空鋁合金機械性能如表 1 所示 。

鋁合金參數
鋁合金參數

3  有限元模型的建立

激光沖擊處理 7075 鋁合金靶材, 由于激光光斑為圓形, 其在材料表面的壓力作用區域亦為圓形, 且沖擊波壓力成中心對稱, 為提高計算機運算速率, 僅對模型1/2進行建模。 本文在建模時, 將模型簡 化為20 ×8 ×4 的立方體, 如圖2 所示, 其中1/2 圓 形為沖擊波壓力作用范圍。 胡永祥[13] 通過研究發 現, 激光噴丸沖擊波主要沿縱向傳播, 因而對縱向 網格尺寸要求遠遠高于橫向網格。 在橫向上, 對激 光噴丸區域, 采用R/Lxy =20 (Lxy為激光噴丸模型 橫向尺寸大小) 進行網格劃分。 縱向上, 越靠近激 光沖擊區域表面需要越細的網格, 而對靠近下表面 的網格尺寸要求較低。 因此, 本文采用 single 的布 種方法, 從模型上表面到下表面, 其種子分布間 距呈等比數列遞增。 為提高運算效率, 對模型非 關鍵部分選用較為粗大的網格進行劃分。 得到的 激光噴丸有限元模型如圖 2 所示, 其單元類型為 C3D8R, 單元數量為 137025, 模型的剖面設置為 對稱邊界。

有限元分析
有限元分析

4  模擬結果及分析

利用有限元軟件 ABAQUS 對激光噴丸過程進行 數值模擬, 當沖擊波峰值壓力為 2. 04 GPa、 光斑直 徑為 Ф4 mm、 脈寬為 10 ns 時, 最終得到的穩定殘 余應力場如圖 3 所示。

有限元分析圖
有限元分析圖

圖 4 為激光噴丸后 7075 鋁合金沿表面方向及深度方向的殘余應力分布曲線。 由圖 3 及圖 4a 可知,當沖擊波峰值壓力為 2. 04 GPa 時, 材料表面殘余應力最大值出現在光斑中心處, 即圖 4a 中的 A 點處,為 – 315. 9 MPa。 隨著距光斑中心距離的增加, 材料表面殘余壓應力逐漸減小, 其原因主要是光斑中心處激光功率密度較高, 產生沖擊波峰值壓力較大, 而邊緣處激光功率密度較低, 產生的沖擊波壓力也較小, 造成光斑中心殘余壓應力較大

激光噴丸分析
激光噴丸分析

而邊緣處殘余應力較小。 在距光斑中心 1 mm 處, 即圖 4a 中 B 點處, 材料表面殘余壓應力減小為 -237.7MPa。 在光斑邊緣處, 如圖4a中C點所 示, 材料表面殘余壓應力進一步減小為 – 30. 9 MPa。 由圖 4b 可知, 材料內部殘余壓應力隨著深 度的加深而逐漸減小,約在0.37mm處(圖4b中 D 點) 減小為 0, 即殘余壓應力影響層深度為 0. 37mm。 當深度進一步加深時材料內部出現殘余拉應力并逐漸較小為 0。

5  激光噴丸參數對殘余應力分布的 影響規律

5. 1  沖擊波峰值壓力對殘余應力分布的影響規律 沖擊波峰值壓力大小是影響殘余應力分布的一 個關鍵性因素。 由激光器基本參數并根據式 (2) 計算可知, 激光器所能提供的最大沖擊波峰值壓力 為 3. 4 GPa, 在模擬時選取沖擊波峰值壓力應小于

該值。?選取沖擊波峰值壓力?P?分別為?0. 5Pmax, 0. 6Pmax, 0. 7Pmax, 0. 8Pmax, 0. 9Pmax,?即?1. 7, 2. 04, 2. 38, 2. 72?和3. 06 GPa?進行有限元模擬,?其 他激光噴丸參數為:?沖擊波壓力作用時間?30 ns、?光 斑大小?Ф4 mm、?沖擊?1?次,?得到不同峰值壓力下殘余應力?σx?沿表面方向及深度方向分布曲線,?如圖?5所示,?由圖?5?可以給出以下結果。

不同峰值下壓應力分布
不同峰值下壓應力分布

(1)?沖擊波峰值壓力的大小對材料表面及內部 殘余應力分布有較大的影響。?當峰值壓力為1. 7 GPa時,?材料表面殘余壓應力最大值為?– 196. 3 MPa,?位于光斑中心處; 材料表面最大殘余壓應力隨沖擊波 峰值壓力的增加而增加。 當沖擊波峰值壓力為 2. 72 GPa 時, 材料表面峰值壓力達到最大。 進一步 增加峰值壓力時, 材料表面殘余壓應力幾乎不再增 加。 其原因主要是沖擊區域材料塑性變形趨于飽和, 使殘余壓應力難以提高。

(2) 材料內部最大殘余壓應力及殘余壓應力影 響層深度隨著沖擊波峰值壓力的增加而增加, 當峰 值壓力為 1. 7 GPa 時, 殘余壓應力影響層深度為 0. 37 mm。 當峰值壓力增加至 3. 06 GPa, 殘余應力 影響層深度增加至 0. 49 mm。

(3) 當峰值壓力為2. 72 和3. 06 GPa 時, 靶材 表面均出現了明顯的殘余應力洞現象, 即殘余壓 應力峰值并不出現在沖擊區域的幾何中心, 而是 出現在距中心一定距離的環形區域。 且隨著沖擊 波峰值壓力的增加, 殘余應力洞現象越明顯。 這 主要是由于沖擊波峰值壓力越大, 沖擊區域邊緣 邊界效應產生的諸如 Rayleigh 波和切應變波等強度越 大, 造成沖擊區域中心材料反向塑性變形越顯著[14]。

5. 2??沖擊次數對殘余應力分布的影響規律?

在某些情況下,?為了使靶材表面及內部達到較 大的殘余壓應力值,?但由于設備功率限制,?無法在 一次沖擊后實現時,?可采用單點多次沖擊[5]?。?為了研究激光沖擊次數?T?對殘余應力分布的影響規律,?分別對靶材進行單點?1 ~ 5?次沖擊進行數值分析,?其 他激光噴丸參數為峰值壓力?2. 2 GPa、?光斑直徑?Ф4 mm、?沖擊波壓力持續時間?30 ns。?得到不同沖擊 次數下殘余應力?σx?沿表面方向及深度方向的分布曲 線,?如圖?6?所示,?可以得出以下結果。

(1) 隨著沖擊次數的增加, 材料表面殘余壓應 力增加。 當沖擊次數T由1增加至5時, 材料表面 最大殘余壓應力值分別為 – 369.3, – 404.2,-424.5, -443.1和-454.7MPa,沖擊2,3,4 和 5 次后殘余壓應力最大值分別相對前一次沖擊后 殘余壓應力最大值增加 9. 45% , 5. 02% , 4. 38% 和2. 61% 。 由此可見, 在沖擊第 2 次時, 材料表面最 大殘余應力增長幅度最大, 為9.45%, 再增加沖擊 次數時, 其增益效果減弱。 其原因是由于多次沖擊 后, 材料產生硬化現象, 再次產生塑性變形必須在 材料表面施加更大的沖擊壓力。

(2)?材料內部最大殘余壓應力及殘余壓應力影 響層深度隨著沖擊次數的增加而增加。?當沖擊次數 從?1?增加至?5?時,?殘余壓應力影響層深度分別為?0.39, 0.47, 0.52, 0.56?和?0.58 mm。?由此可見,隨著沖擊次數的增加,?殘余壓應力影響層深度增益 效果減弱。

分布圖
分布圖

5. 3?光斑大小對殘余應力分布的影響規律

為了探索光斑大小對激光噴丸殘余應力分布的 影響規律, 選取直徑 D 分別為 Ф2, Ф4, Ф6 和 Ф8 mm 的光斑進行數值分析。 其他激光噴丸工藝參 數為: 峰值壓力 2. 38 GPa、 沖擊波壓力持續時間 30 ns、 沖擊 1 次, 得到的不同光斑大小下殘余應力 σx沿表面方向及深度方向的分布曲線, 如圖7所示, 由圖 7 可以得出以下結果。

(1)?當光斑直徑分別為?Ф2,?Ф4,?Ф6?和?Ф8 mm?時,?材料表面最大殘余應力分別為?– 363.2, -421.2, – 464.9?和?– 479.3 MPa。?光斑直徑為Ф4,?Ф6,?和?Ф8 mm?分別相對光斑半徑為?R2 mm?時表面最大殘余壓應力增加?15.9%, 28.0%?和31. 97% 。?由此可知,?材料表面殘余壓應力隨著光斑尺寸的增加而增加。?當光斑尺寸增大到一定值時,其表面最大殘余壓應力增益效果減弱。?其原因主要是小直徑光斑下沖擊波在材料內部傳播類似球面波,大光斑直徑下沖擊波在材料內部傳播類似平面波,球面波的衰減速率遠大于平面波。

2) 當光斑尺寸分別為 Φ2, Φ4, Φ6 和Φ8 mm 時, 殘余壓應力影響層深度分別為 0. 36, 0. 40, 0. 43, 和 0. 46 mm。 由此可見, 在峰值壓力一定時,殘余應力影響層深度隨著光斑尺寸的增加而增加。

(3) 當光斑直徑為 Φ2 mm 時, 靶材表面出現 了非常明顯的殘余應力洞現象, 且殘余應力洞現象 隨著光斑尺寸的增大而逐漸消失。 這主要是由于光 斑直徑越大, 光斑邊界效應越小, 產生的稀疏波的 強度越低, 且大光斑下稀疏波傳播行程長, 匯聚到 光斑中心時強度進一步減弱, 造成沖擊區域中心反 向塑性變形較小。

5. 4??光斑搭接率對殘余應力分布的影響規律?

在激光噴丸過程中,?為了使靶材表面沖擊波峰 值壓力達到某一大小,?激光器需保證一定的功率密 度。?由于激光器本身功率限制,?在激光噴丸時選用 的光斑不宜過大,?單個光斑不可能對整個待沖擊區 域一次性沖擊強化處理。?因而,?在實際激光噴丸過 程中需采用一定的光斑搭接處理工藝。?光斑搭接率計算公式為[8]?:

p =?(1?– L / D)?×?100%?(4)

式中:?L?為兩相鄰光斑間距離;?D?為光斑直徑。

為研究光斑搭接率對殘余應力分布的影響規律,?分別采用搭接率為?0% , 25% , 50%?和?75%?進行模 擬,?激光噴丸參數選擇沖擊波峰值壓力?2. 04 GPa、?光斑直徑?Ф?4 mm、?沖擊波壓力作用時間?30 ns、?沖擊?1?次,?得到的不同搭接率時材料表面殘余應力?σx?沿表面?X?方向分布規律,?如圖?8?所示,?圖中圓形區 域為激光沖擊區域。?圖?9?為不同搭接率下殘余應力?σx?沿材料表面?X?方向分布曲線對比。

分布規律
分布規律

由圖8a可知,?當光斑搭接率為0%時,?由于相鄰光斑之間不存在任何搭接,?因而各次激光噴丸過程是相對獨立的過程,?材料表面殘余壓應力極大值出現在光斑中心處,?極小值出現在各激光沖擊區域邊緣處,?且材料表面殘余壓應力波動較明顯。?由圖8b可知,?當光斑搭接率為25%時,?由于相鄰光斑之 間存在少量搭接,?且各激光沖擊區域中心不存在搭 接,?因而材料表面殘余壓應力極大值出現在激光沖 擊區域中心處。?而在沖擊邊緣處,?由于與后一沖擊 區域光斑半徑中點搭接,?使該處殘余壓應力得到較 大改善,?殘余壓應力極小值并未出現在光斑邊緣處,?而是出現在相鄰兩光斑搭接區域中點處。?且材料表 面殘余壓應力波動值較搭接率為?0%?時減小很多。?由圖?8c?可知,?當搭接率為?50%?時,?前一沖擊區域 中心正好與下一沖擊區域邊緣搭接,?前一沖擊區 域邊緣恰好與后一沖擊區域中心搭接,?由于沖擊 區域中心殘余壓應力較大,?沖擊區域邊緣壓應力 較小,?搭接后材料表面殘余應力波動值較搭接率 為25%時有明顯改善。?殘余壓應力極大值出現在 沖擊區域中心,?極小值出現在光斑半徑中點處。?由圖8d可知,?當光斑搭接率為75%時,?前1?個沖 擊區域中心正好與后?1?個沖擊區域半徑中點處及 其后第?2?個沖擊區域邊緣處搭接,?由于沖擊區域中心殘余應力值較大, 搭接后材料表面殘余應力極 大值出現在該區域。 材料表面殘余應力波動值在搭 接率為 75% 時更小。

分布圖
分布圖

由圖 9 可知, 材料表面整體殘余壓應力隨著光 斑搭接率的增加而增加。 其原因主要是光斑搭接率 越大, 同一區域激光沖擊次數越多, 因而其殘余壓 應力越大。

6  結論

(1) 材料表面殘余壓應力隨沖擊波峰值壓力的 增加而增加。 當峰值壓力增加至一定值 (如7075 航 空鋁合金靶材當光斑半徑為 R4 mm 時、 峰值壓力為 2. 72 GPa) 時, 材料表面殘余壓應力達到最大。 過 大的沖擊波峰值壓力作用在材料表面會在沖擊區域 中心處產生極為明顯的殘余應力洞現象, 在實際激 光噴丸過程中應盡量避免。

(2) 隨著沖擊次數的增加, 材料表面殘余壓應 力增加。 但當沖擊次數達到一定值 (本文中對應沖 擊次數為 3 次) 時, 其增益效果減弱。 在實際激光 噴丸過程中, 為提高生產效率, 應盡量避免采用過 多的沖擊次數。

(3)?光斑搭接率是影響材料表面殘余應力分布 的一個重要參數。?當搭接率為?75%?時,?材料表面殘 余應力較大,?且沿表面方向波動較小,?在實際操作 中應優先選取。

江蘇激光噴丸設備制造商,激光噴丸代工服務商,竭誠為您服務

噴拋丸設備、噴丸強化技術服務,182 0189 8806 文章鏈接:http://www.8beet.com/1752/

發表評論

登錄后才能評論

聯系我們

131-6246-0999

在線咨詢:點擊這里給我發消息

微 ? ? ? ?信:DS2060798

工作時間:24小時服務在線

QR code
error: Content is protected !!
91麻豆成人久久精品二区三区_亚洲成人av一区二区_国产不卡人人_一区二区欧美日韩视频
一区二区三区欧美激情| 天堂蜜桃一区二区三区| 性久久久久久久久久久久| 91免费国产在线| 亚洲卡通欧美制服中文| 欧美日本韩国一区二区三区视频| 五月婷婷色综合| 久久嫩草精品久久久精品| 国产盗摄视频一区二区三区| 国产精品久久久久影院色老大 | 成人激情综合网站| 亚洲人成7777| 日韩欧美在线一区二区三区| 国产91精品精华液一区二区三区 | 欧美亚洲国产bt| 蜜臀91精品一区二区三区| 国产精品久久久久一区 | 在线不卡一区二区| 国产精品123| 亚洲成av人片在线观看无码| 久久综合色婷婷| 欧美日韩一区小说| 国产sm精品调教视频网站| 日韩成人一区二区| 亚洲免费观看高清在线观看| 欧美成人性战久久| 在线观看日韩精品| 成人av在线电影| 精一区二区三区| 亚洲制服丝袜一区| 国产精品成人免费在线| 欧美精品一区二区在线播放| 欧美日韩在线一区二区| zzijzzij亚洲日本少妇熟睡| 国产一区二区日韩精品| 日韩电影在线一区| 亚洲成av人片在www色猫咪| 国产精品乱码久久久久久| 欧美一区二区福利视频| 欧美日韩国产在线观看| 色哟哟一区二区| 91亚洲国产成人精品一区二区三| 国产精一品亚洲二区在线视频| 日韩中文字幕麻豆| 亚洲第一成年网| 伊人色综合久久天天| 中文字幕一区在线观看| 日本一区二区久久| 国产女人18水真多18精品一级做 | caoporn国产一区二区| 国产一区二区电影| 精品影院一区二区久久久| 日本特黄久久久高潮| 日本美女一区二区| 日本vs亚洲vs韩国一区三区| 视频一区视频二区在线观看| 午夜久久久久久电影| 午夜国产精品影院在线观看| 丝袜a∨在线一区二区三区不卡| 亚洲成人先锋电影| 亚洲成人av中文| 青青青伊人色综合久久| 奇米色777欧美一区二区| 麻豆一区二区在线| 国产一区二区三区综合| 成人午夜私人影院| 91日韩在线专区| 欧美日韩亚洲综合在线| 欧美一区二区精品在线| 精品国产区一区| 国产欧美精品日韩区二区麻豆天美| 国产人伦精品一区二区| 日韩毛片一二三区| 午夜精品在线视频一区| 国产麻豆精品一区二区| 99re这里都是精品| 7777精品伊人久久久大香线蕉的 | 不卡的av电影| 91激情五月电影| 在线不卡欧美精品一区二区三区| 日韩视频免费直播| 国产精品色哟哟| 午夜精彩视频在线观看不卡| 韩国精品在线观看| 91浏览器在线视频| 日韩视频中午一区| 亚洲精品自拍动漫在线| 麻豆一区二区三| 色欧美日韩亚洲| 久久精品一二三| 亚洲一二三专区| 国产精品夜夜嗨| 欧美精品日日鲁夜夜添| 日本一区二区三区免费乱视频 | 国产乱对白刺激视频不卡| 91在线国内视频| 久久在线观看免费| 偷窥国产亚洲免费视频 | 成人国产一区二区三区精品| 欧美日韩高清在线| 成人欧美一区二区三区视频网页| 日本怡春院一区二区| 色网综合在线观看| 国产亚洲一区二区三区| 奇米一区二区三区| 欧美网站大全在线观看| 日本一区二区三区dvd视频在线| 首页综合国产亚洲丝袜| 一本大道综合伊人精品热热| 国产丝袜在线精品| 麻豆成人久久精品二区三区红| 在线亚洲高清视频| 国产精品久久精品日日| 国产一区 二区 三区一级| 日韩一区二区三区四区五区六区| 亚洲精品网站在线观看| av电影在线观看一区| 中文字幕av不卡| 成人综合婷婷国产精品久久 | 亚洲国产成人91porn| 91小视频在线| 亚洲丝袜制服诱惑| 99精品视频一区二区三区| 中文无字幕一区二区三区 | 欧美精品丝袜中出| 亚洲一区二区三区四区五区黄| av高清不卡在线| 一区二区三区中文字幕| 色嗨嗨av一区二区三区| 亚洲综合精品久久| 在线观看视频一区二区欧美日韩| 伊人婷婷欧美激情| 欧美肥大bbwbbw高潮| 日本欧美久久久久免费播放网| 91精选在线观看| 蜜臂av日日欢夜夜爽一区| 日韩欧美不卡一区| 国产精品自在在线| 国产精品女主播在线观看| 9久草视频在线视频精品| 亚洲乱码国产乱码精品精可以看| 在线观看亚洲一区| 蜜臀av亚洲一区中文字幕| 精品国产一区二区三区四区四| 国产精品一区不卡| 亚洲日本在线a| 欧美日韩国产三级| 激情成人综合网| 一区精品在线播放| 欧美日韩免费不卡视频一区二区三区| 亚洲va欧美va天堂v国产综合| 91精品欧美福利在线观看| 国产精品18久久久| 亚洲国产一区视频| 精品国产乱码久久久久久老虎 | 成人免费av资源| 亚洲国产欧美在线| 精品国产区一区| 97se亚洲国产综合自在线| 天天影视涩香欲综合网| 久久综合五月天婷婷伊人| 91欧美一区二区| 久久电影网站中文字幕| 中文字幕一区二区三区色视频| 欧美专区在线观看一区| 国产美女精品在线| 亚洲地区一二三色| 中文字幕一区二区三区乱码在线| 91精品国产一区二区三区蜜臀 | 亚洲男同性视频| 26uuu色噜噜精品一区二区| 日本二三区不卡| 国产成人免费视| 日本特黄久久久高潮| 一区二区三区欧美在线观看| 久久欧美一区二区| 91精品国产综合久久久久久| 99九九99九九九视频精品| 国产一区不卡精品| 五月婷婷久久丁香| 亚洲欧美一区二区久久 | 99re成人精品视频| 精品中文字幕一区二区| 午夜一区二区三区在线观看| 亚洲欧洲精品一区二区三区| 久久伊人蜜桃av一区二区| 91麻豆精品91久久久久同性| 欧美最猛性xxxxx直播| 99re免费视频精品全部| 国产成a人亚洲精品| 激情六月婷婷综合| 极品少妇一区二区三区精品视频| 亚洲国产精品视频| 国产精品水嫩水嫩| 久久久99久久| 久久色中文字幕| 久久久久综合网| 久久―日本道色综合久久| 91精品国产综合久久精品图片 | 91啪九色porn原创视频在线观看|