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一区二区_国产不卡人人_一区二区欧美日韩视频
欧美精品在线免费| 欧美一区成人| 国产精品久久久久久福利一牛影视| 美女日韩在线中文字幕| 美女精品在线| 欧美日韩在线一区| 国产综合色一区二区三区| 1000部精品久久久久久久久| 亚洲国产成人在线播放| 99视频精品免费观看| 久久精彩视频| 国产精品xnxxcom| 亚洲成色www8888| 性视频1819p久久| 欧美日韩亚洲一区二区三区| 国产一区二区剧情av在线| 日韩午夜精品| 欧美精品日本| 亚洲国产精品嫩草影院| 久久成人精品| 国产精品五月天| 一本色道久久88精品综合| 欧美成人精品1314www| 国产一区二区三区在线观看免费视频| 亚洲一区中文| 国产精品久久波多野结衣| 99视频在线精品国自产拍免费观看 | 久久av一区二区三区| 欧美交受高潮1| 国产亚洲欧美激情| 欧美在线91| 国产综合久久| 久久久久久久欧美精品| 激情综合五月天| 久久躁日日躁aaaaxxxx| 国产亚洲精品aa午夜观看| 欧美一区二区三区四区在线| 国产日韩视频一区二区三区| 久久精品72免费观看| 国产一区二区三区久久悠悠色av | 亚洲精品日韩在线观看| 欧美电影免费观看高清完整版| 亚洲精品美女| 国产精品久久久久一区二区三区共 | 欧美日韩一区二区在线视频 | 欧美日韩中文在线观看| 夜夜爽www精品| 国产精品久久久久久亚洲毛片 | 国产日韩精品视频一区| 小黄鸭精品aⅴ导航网站入口| 欧美日韩国产首页在线观看| 亚洲国产欧美不卡在线观看| 国内精品久久久久影院 日本资源| 在线综合+亚洲+欧美中文字幕| 欧美激情亚洲一区| 一本久久青青| 国产精品日韩精品欧美精品| 久久精品国产综合精品| 亚洲国产精品视频| 欧美系列一区| 久久久久久久久久久一区 | 久久国产精品高清| 国产一区二区久久久| 欧美**人妖| 亚洲综合社区| 在线欧美三区| 国产精品系列在线播放| 老司机凹凸av亚洲导航| 亚洲午夜精品久久| 亚洲电影免费在线| 国产精品高清在线观看| 蜜桃久久av一区| 亚洲女女女同性video| 亚洲国产裸拍裸体视频在线观看乱了中文 | 欧美日韩国产一级片| 性色av一区二区三区在线观看| 亚洲国产91色在线| 国产精品综合久久久| 欧美日韩国产免费| 久久亚洲一区二区| 香蕉成人伊视频在线观看 | 国产精品久久久久久户外露出 | 亚洲精品一区在线观看| 国产日产欧美精品| 久久精彩视频| 99精品免费| 欧美激情综合亚洲一二区| 国产日韩欧美不卡| 老司机精品久久| 一本大道久久a久久精品综合| 欧美国产三区| 久久国产黑丝| 激情五月综合色婷婷一区二区| 欧美激情综合网| 久久在线91| 久久精品中文| 久久国产高清| 久久av在线看| 欧美在线黄色| 久久婷婷成人综合色| 在线观看日韩专区| 欧美国产综合视频| 亚洲一区国产一区| 国产模特精品视频久久久久| 欧美国产日本高清在线| 国产亚洲欧美在线| 男人的天堂成人在线| 久久国产88| 久久久噜噜噜久久人人看| 午夜精品视频在线观看一区二区| 亚洲视频香蕉人妖| 在线视频精品一| 亚洲校园激情| 亚洲免费在线观看视频| 欧美亚洲免费| 性色av一区二区三区| 国产精品久久久久久一区二区三区| 香蕉成人啪国产精品视频综合网| 亚洲欧美一区二区三区极速播放| 久久精品综合网| 国产精品欧美精品| 国产精品毛片大码女人| 国产精品爽黄69| 国模私拍一区二区三区| 国产欧美一区二区在线观看| 狠狠色狠狠色综合日日五| 亚洲国产成人精品视频| 99在线精品视频在线观看| 亚洲欧美另类中文字幕| 久久久一二三| 欧美午夜精品久久久久久人妖 | 亚洲精品在线电影| 最新高清无码专区| 久久嫩草精品久久久精品| 久久全国免费视频| 欧美激情精品久久久久久久变态| 欧美激情第3页| 国产精品自拍小视频| 亚洲国产成人精品久久久国产成人一区 | 欧美成人免费网| 国产精品成人观看视频国产奇米| 国产精品有限公司| 亚洲国产小视频| 欧美与欧洲交xxxx免费观看| 欧美成人tv| 韩日成人在线| 亚洲一区二区三区久久 | 一区在线视频| 亚洲已满18点击进入久久| 久久全球大尺度高清视频| 国产精品www.| 亚洲人成网站色ww在线| 久久精品国产96久久久香蕉| 国产精品成人午夜| 亚洲精品久久| 久久一二三四| 国产伦精品一区二区三区照片91| 日韩视频三区| 欧美金8天国| 在线电影院国产精品| 欧美一区二区三区四区夜夜大片| 欧美午夜激情视频| 在线欧美日韩| 久久综合中文字幕| 国产一区二区三区丝袜| 午夜视频在线观看一区二区三区| 久久偷看各类wc女厕嘘嘘偷窃| 国产精品久久久久国产a级| 亚洲国产高清一区| 免费在线亚洲| 亚洲国产成人在线| 欧美电影免费| 亚洲精品一区二区在线| 男人的天堂亚洲在线| 国产亚洲视频在线观看| 欧美一级夜夜爽| 国产精品区一区二区三| 亚洲视频香蕉人妖| 欧美视频精品一区| 日韩一级二级三级| 欧美成人一区在线| 亚洲国产一区二区视频| 免费成人黄色av| 亚洲美女啪啪| 欧美理论电影在线播放| 亚洲精品一区二区三区在线观看 | 亚洲一区二区三区免费观看| 欧美午夜不卡| 亚洲影视综合| 国产精品嫩草影院av蜜臀| 亚洲一区久久久| 国模精品娜娜一二三区| 美国成人毛片| 亚洲人成高清| 国产精品欧美日韩一区| 久久九九免费| 一区二区精品在线观看| 国产精品试看| 麻豆国产精品一区二区三区| 日韩一级精品视频在线观看|