999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

汽車轉向節臂端部冷精整工藝優化

2022-09-07 08:22:48宛加雄武建祥晏洋丁華鋒
精密成形工程 2022年9期
關鍵詞:變形工藝

宛加雄,武建祥,晏洋,丁華鋒,

汽車轉向節臂端部冷精整工藝優化

宛加雄1,武建祥2,晏洋2,丁華鋒1,2

(1.純電動汽車電力系統設計與測試湖北省重點實驗室,湖北 襄陽 441053;2.湖北三環鍛造有限公司,湖北 襄陽 441700)

得到單向壓縮冷精整的精整量與預變形孔半徑和深度之間的變化規律,以及摩擦因數大小同回彈量之間的關系。對轉向節臂端部鍛件上下表面進行預變形凹槽孔處理,并建立鍛件單向壓縮冷精整的彈塑性有限元模型,利用ABAQUS有限元軟件進行數值模擬仿真,分析摩擦因數對表面質量的影響,得到優化后的工藝參數,并進行實驗驗證模擬結果的準確性。摩擦因數越大,冷精整后的轉向節臂端部鍛件的鼓形越明顯。在精整量為1 mm的情況下,當預變形凹槽孔的半徑為4 mm、深度為0.3 mm、單向壓縮量為1.1 mm、摩擦因數為0.2時,能夠得到最好的表面質量。當摩擦因數相同時,預變形凹槽孔的半徑越大,徑向位移越小,而回彈量基本保持不變。

汽車轉向節臂;冷精整;工藝優化;回彈量

轉向節臂在汽車行駛過程中起重要作用,是汽車上重要的安保件,其幾何精度直接影響裝配效果。隨著有限元技術的發展,學者們利用CAE軟件進行了大量研究[1-5]。章志兵等[6]基于CAE模型船舶艙室的相關研究,提出了一種半封閉空間的新概念,并設計了全船CAE模型的識別算法,實現了半面擴展算法的高效性和準確性。劉江等[7]針對轉向節在實際生產過程中出現的塌角等缺陷問題,通過數值模擬和實驗驗證提出了復合成形工藝,使材料利用率達到了80.8%,成品合格率達到了97.9%。趙毅等[8]分析了某皮卡的轉向節熱鍛成形過程,采用大小臺階式阻力模具設計方法,解決了鍛件填充不滿的問題。在實際生產中,通常采用冷精壓的方法來提高轉向節臂端部的幾何精度。王德林等[9]分析了三銷軸叉冷精整工藝過程中型腔底部應力、應變的變化情況,優化設計了不同的入模角,使冷精整后的型腔底部范圍區域沒有產生裂紋,顯著提高了冷精整鍛件的產品質量。Stone等[10]通過ABAQUS軟件建立了厚壁圓筒的有限元模型,研究了精整量對表面質量的影響,結果表明,冷精整能夠改善工件的表面質量且不會對工件外部形狀產生影響。朱懷沈等[11]研究了冷精整工藝過程中直齒輪的齒面冷精整量與齒面之間回彈量的變化關系,得到了精整量與壁厚之間的優化選擇。郭嘉晨等[12]用Kriging模型和粒子群算法改善了冷精整工藝中萬向節滑套杯壁的表面質量。Kwanghyun等[13]整合了歐拉光束理論并預測了螺旋管蒸汽發生器的彎曲回彈情況,將流變應力和彈性模量的數值應用到分析模型當中,精確預測了最終尺寸。Farhad等[14]建立了非線性運動模型,發現改變載荷的順序能夠明顯減小彎曲回彈角。Soheil等[15]和Uten等[16]研究發現,鍛造時,在低摩擦條件、不考慮折疊的情況下,當最大摩擦因數≤0.5時,Avitzur模型具有最高的準確性。Fras等[17]研究發現,高載荷下表面之間的直接接觸會導致摩擦力顯著增加,對工件會造成不可預估的缺陷。在實際的冷精整過程中,壓頭與工件表面存在摩擦,這會導致工件的精整表面粗糙,而摩擦因數的大小對精整表面的影響很大。

目前針對汽車轉向節臂冷精整工藝優化的研究還比較少,上述研究忽略了摩擦因數對表面質量的影響。文中采用在轉向節臂端部上下表面設置預變形凹槽的方法來提高表面質量,將壓縮量作為關鍵參數,對轉向節臂端部圓棒進行單向壓縮冷精整加工,通過ABAQUS有限元仿真軟件分析孔的半徑和深度對回彈量的影響,同時研究摩擦因數對表面質量的影響規律,表面質量通過冷精整回彈后轉向節臂端部上表面節點的軸向位移表現。

1 轉向節臂端部冷精壓工藝

1.1 轉向節臂結構

轉向節臂鍛件圖如圖1所示。轉向節臂是汽車轉向系統的最后一級傳力部件,轉向節臂一端安裝在上下轉向節之間,另外一端則用球銷與橫拉桿連接。因此端部的表面質量將直接影響轉向節臂的力學性能。

圖1 轉向節臂鍛件圖

1.2 塑性應變和摩擦影響

對鍛件進行單向壓縮,由于外端表面受摩擦力的影響,變形難度增加,因此會發生凸緣現象。單向壓縮的初始和最終形狀如圖2所示,得到的曲率半徑表達式如式(1)所示。

式中:H1為最終高度;DM為最大直徑;D1為頂部直徑。

摩擦因數的測定對塑性成形至關重要[18-21]。在轉向節臂端部單向壓縮冷精整過程中,當不考慮摩擦條件時,冷精整為理想單向壓縮模型,此時轉向節臂端部軸向尺寸變化相同。隨著摩擦因數的增大,接觸面外端變形困難,使上下接觸表面變形小于中間部分的變形,進而產生鼓形。隨著摩擦因數的進一步增大,鼓形越加明顯。下面將通過有限元數值模擬得到摩擦因數大小對轉向節臂回彈量的定量影響規律。

2 冷精壓有限元仿真

在轉向節臂端部冷精壓過程中,在摩擦力的影響下,精整表面質量較為粗糙。為得到冷精整時的變形規律及較高的表面質量,在轉向節臂端部先加工一個預變形凹槽孔,研究預變形凹槽孔的開口大小、深度以及表面摩擦因數對回彈量的影響,以實現冷精整量的優化選擇。文中轉向節臂端部外徑為57 mm,軸向長度為35 mm,冷精整后的軸向長度為34 mm,端部外圓直徑增大不超過1%。設計預變形凹槽孔的半徑為2、3、4 mm,深度為0.5、0.3 mm,研究摩擦因數分別為0、0.1、0.2、0.3、0.4時轉向節臂端部的回彈大小。文獻[22]表明,摩擦因數通常小于0.4,因此文中的摩擦因數最大取值為0.4。

圖3 轉向節臂端部精整模型

2.1 彈塑性有限元精整模型的建立

在冷精整成形過程中,轉向節臂的內部發生了彈塑性變形,積攢了彈性勢能,當卸除載荷后,彈性勢能釋放從而發生回彈,因此用彈塑性模型分析。材料為42CrMo高強度合金鋼,其屈服強度為930 MPa、抗拉強度為1 080 MPa、密度為7.85 g/cm3、彈性模量為2.02×105MPa、泊松比為0.28,材料為均質且各向同性,在整個過程中不計算體積力與慣性力,遵循米塞斯屈服準則。由于轉向節臂端部的圓柱體是軸對稱零件,因此在模擬時選擇二維坯料,只分析1/2模型以減少計算量[23-25],再對中間預變形凹槽孔進行局部網格的細化處理。冷精整開始前先設置一步靜接觸,使接觸關系平穩建立,并設置不同的摩擦因數來進行對比分析。

在冷精整過程中,金屬一部分沿徑向方向向內流動,填充預變形凹槽孔,另一部分沿著徑向方向向外流動,使徑向尺寸增大。轉向節臂端部精整過程位移場分布情況如圖4所示,其內部金屬組織以軸向方向流動為主,以徑向方向流動為輔。

2.2 工藝參數的影響規律研究

為控制轉向節臂端部冷精整的尺寸精度,研究轉向節臂端部的回彈規律。根據仿真分析結果,得到不同預變形孔的大小和深度,以及在表面摩擦因數不同的條件下,單向壓縮冷精整后徑向位移1和軸向位移2的回彈量數據。

轉向節臂端部預變形凹槽孔半徑為0.5 mm和0.3 mm時,不同半徑下的回彈量大小分別見圖5和圖6。根據回彈量數據可知,當預變形凹槽孔的半徑相同時,預變形凹槽孔深度的影響可以忽略不計,徑向位移和回彈量相等。在單向壓縮精整過程中,由于精整量為1 mm,深度對轉向節臂回彈量影響較小,仿真模擬預變形孔深度最小為0.3 mm,預變形凹槽孔能夠有效提高精整后轉向節臂端部上下表面質量。預變形凹槽孔的半徑越大,徑向位移越小,而回彈量基本不變。當摩擦因數為0.2時,徑向位移的回彈量最大。摩擦因數越大,徑向位移也就越大。因此,由數據對比得出,在半徑為4 mm、深度為0.3 mm時,徑向位移1和軸向位移2的綜合效果最好。

圖4 轉向節臂端部精整過程位移場分布圖

圖5 預變形孔深度為0.5 mm時的回彈量

圖6 預變形孔深度為0.3 mm時的回彈量

在相同的精整量下,隨著摩擦因數的增大,上下表面變形困難,彈塑性變形過程中的彈性變形量先增大后減小,因而相應的回彈量也就呈先增大后減小的變化趨勢。在深度為0.3 mm的精整條件下,3種不同半徑大小的預變形凹槽孔在摩擦因數為0~0.4時的徑向回彈量和軸向回彈量分別見圖7和圖8。由圖7和圖8可知,回彈量隨著摩擦因數的增大,總體呈現先增大后減小的趨勢,當摩擦因數為0.2時,回彈量最大。通過對比3種不同預變形凹槽孔半徑的大小,得到在預變形凹槽孔半徑為4 mm時,摩擦因數對回彈量影響最小。

考慮到摩擦對表面粗糙度的影響,仿真精整過程中上表面的軸向位移隨時間的變化情況如圖9所示。由圖9可以看出,在前1 s的時間內,壓頭與轉向節臂端部平穩接觸,為降低速度對其產生的影響,將壓縮精整的過程控制在3 s內,隨后卸載壓頭,出現了回彈現象。

當精整量為1 mm時,壓力卸載后出現了回彈現象,由于摩擦的影響,回彈結束后轉向節端部上表面的位移并沒有重合。為了準確地控制回彈量,將精整量確定為1.1 mm,預變形凹槽孔的深度為0.3 mm,得到的表面節點的軸向位移如圖10所示。

圖7 不同摩擦因數下的徑向回彈量

圖8 不同摩擦因數下的軸向回彈量

圖9 精整過程中軸向位移隨時間變化圖形

圖10是摩擦因數不同時,轉向節臂端部圓棒上表面軸向位移值。由圖10可知,在摩擦因數為0的理想狀態下,整個上表面位移近似相同。當摩擦因數為0.2時,得到的位移值最接近1 mm。當摩擦因數不同時,轉向節臂端部圓棒徑向位移值如圖11所示,根據圖11可知,隨著摩擦因數的增大,轉向節臂端部在中間位置的位移增大,當摩擦因數為0時則不會產生鼓形,隨著摩擦因數的增大,出現的鼓形越明顯。由于摩擦的存在,轉向節臂端部圓棒的上下表面與壓頭接觸之間的摩擦切應力會導致靠近上下表面的位移很小。

圖10 轉向節臂端部圓棒上表面軸向位移

圖11 轉向節臂端部圓棒徑向位移

3 實驗結果驗證

對轉向節臂零件進行實物實驗,以驗證上述的仿真結果。首先在轉向節臂的端部兩端加工預變形凹槽孔,然后沿著厚度的方向對轉向節臂端部進行冷精整,在冷精整過程中金屬沿著端部外圓和預變形凹槽孔內徑的方向流動,最后得到的汽車轉向節臂如圖12所示,測量其端部尺寸。

圖12 冷精整后的轉向節臂

考慮到實驗成本,為減少實驗時間,根據上述的模擬值,確定轉向節臂端部預變形凹槽孔的半徑為4 mm、深度為0.3 mm,對多個轉向節臂試件進行實驗,得到的實驗測量數據如下:徑向回彈量的模擬值為0.024 mm,測量值為0.023 1 mm,相對誤差為3.8%;軸向回彈量的模擬值為0.096 3 mm,測量值為0.098 4 mm,相對誤差為2.2%。由此能夠看出,轉向節臂的回彈測量值與實驗模擬值的相對誤差皆不超過5%,加工預變形凹槽孔能夠有效保證轉向節臂端部的表面質量,通過實驗驗證了模擬結果的準確性,實驗的結果能夠指導實際工業生產。

4 結論

通過彈塑性有限元分析,得到了摩擦因數與回彈量之間的變化關系以及預變形孔的大小對回彈的影響規律。

1)當摩擦因數相同時,預變形凹槽孔的半徑越大,徑向位移越小,而回彈量基本保持不變。在冷精整過程中,當預變形凹槽孔的半徑相同時,預變形孔的深度影響可以忽略不計。

2)當摩擦因數為0.2時,徑向位移的回彈量最大。隨著摩擦因數的增大,轉向節臂端部鍛件鼓形越明顯。

3)當預變形凹槽孔的半徑為4 mm、深度為0.3 mm、摩擦因數為0.2時,轉向節臂端部鍛件的冷精整達到了工藝要求。基于所得到的結論,可以對精整量進行優化選擇,并應用于實際的工藝工程上。

[1] 郭晶玉, 鄧小虎, 鄭寶星, 等. Ti55531鈦合金扭力臂熱鍛成形工藝設計及優化[J]. 精密成形工程, 2021, 13(2): 96-104.

GUO Jing-yu, DENG Xiao-hu, ZHENG Bao-xing, et al. Design and Optimization of Hot Forging Forming Process of Ti55531 Titanium Alloy Torsion Arm[J]. Journal of Netshape Forming Engineering, 2021, 13(2): 96-104.

[2] 楊天云, 張晴朗, 楊兵, 等. CAE在汽車座椅輕量化設計中的應用[J]. 精密成形工程, 2012, 4(1): 73-77.

YANG Tian-yun, ZHANG Qing-lang, YANG Bing, et al. Application of CAE in Lightweight Design of Car Seats[J]. Journal of Netshape Forming Engineering, 2012, 4(1): 73-77.

[3] 王江波, 強寶民. 基于ABAQUS圓柱殼減振裝置振動特性研究[J]. 包裝工程, 2017, 38(17): 112-117.

WANG Jiang-bo, QIANG Bao-min. Research on Vibration Characteristics of Cylindrical Shell Vibration Reduction Device Based on ABAQUS[J]. Packaging Engineering, 2017, 38(17): 112-117.

[4] 李炎輝, 胡青春. PET瓶基于ABAQUS的軸壓分析[J]. 包裝工程, 2009, 30(10): 18-20.

LI Yan-hui, HU Qing-chun. Axial Compression Analysis of PET Bottle Based on ABAQUS[J]. Packaging Engineering, 2009, 30(10): 18-20.

[5] 趙寧, 張運軍, 陳天賦, 等. V5型轉向節臂部鍛造開裂研究[J]. 塑性工程學報, 2017, 24(3): 7-12.

ZHAO Ning, ZHANG Yun-jun, CHEN Tian-fu, et al. Research on Forging Cracking of V5 Steering Knuckle Arm[J]. Chinese Journal of Plastic Engineering, 2017, 24(3): 7-12.

[6] 章志兵, 楊明潤, 王麗榮, 等. 基于CAE模型的船舶艙室自動識別高效算法[J]. 精密成形工程, 2020, 12(4): 160-165.

ZHANG Zhi-bing, YANG Ming-run, WANG Li-rong, et al. High-Efficiency Algorithm for Automatic Identification of Ship Cabins Based on CAE Model[J]. Journal of Netshape Forming Engineering, 2020, 12(4): 160-165.

[7] 劉江, 徐皓. 基于Deform 3D的長桿類汽車轉向節鍛模設計及鍛造工藝生產驗證[J]. 鍛壓技術, 2021, 46(2): 9-13.

LIU Jiang, XU Hao. Design of Forging Die and Forging Process Production Verification of Long-Rod Automobile Steering Knuckle Based on Deform 3D[J]. Forging Technology, 2021, 46(2): 9-13.

[8] 趙毅, 劉淑梅, 何文濤, 等. 汽車轉向節熱鍛工藝分析及模具結構優化設計[J]. 精密成形工程, 2016, 8(3): 40-44.

ZHAO Yi, LIU Shu-mei, HE Wen-tao, et al. Hot Forging Process Analysis and Die Structure Optimization Design of Automobile Steering Knuckles[J]. Journal of Netshape Forming Engineering, 2016, 8(3): 40-44.

[9] 王德林, 陸有根, 吳兵. 三銷軸叉精整底部裂紋產生和防護[J]. 精密成形工程, 2018, 10(1): 177-180.

WANG De-lin, LU You-gen, WU Bing. The Occurrence and Protection of Bottom Cracks in the Finishing of Three-Pin Shaft Fork[J]. Journal of Netshape Forming Engineering, 2018, 10(1): 177-180.

[10] STONE E, CAI J, HU Z M, et al. An Exercise in Cold Ironing as the Post-Forging Operation for Net-Shape Manufacture[J]. Journal of Materials Processing Technology, 2003, 135(2/3): 278-283.

[11] 朱懷沈, 夏巨諶, 金俊松, 等. 大模數直齒輪溫冷鍛精整量的優化選擇[J]. 塑性工程學報, 2011, 18(1): 53-57.

ZHU Huai-shen, XIA Ju-chen, JIN Jun-song, et al. Optimization of Sizing Amount during Warm-Cold Compound Forging for Supr Gear With Big Modulus[J]. Journal of Plasticity Engineering, 2011, 18(1): 53-57.

[12] 郭嘉晨, 汪建敏, 潘小迎, 等. 基于Kriging模型的等速萬向節滑套冷精整工藝的優化[J]. 鍛壓技術, 2017, 42(5): 176-180.

GUO Jia-chen, WANG Jian-min, PAN Xiao-ying, et al. Optimization of Cold Finishing Process of Constant Velocity Universal Joint Sliding Sleeve Based on Kriging Model[J]. Forging Technology, 2017, 42(5): 176-180.

[13] KWANGHYUN A, KANG-HEON L, JAE-SEON L, et al. Analytic Springback Prediction in Cylindrical Tube Bending for Helical Tube Steam Generator[J]. Nuclear Engineering and Technology, 2020, 52(9): 2100-2106.

[14] FARHAD A, NAYEBI A. Springback Analysis of Thick-Walled Tubes under Combined Bending-Torsion Loading with Consideration of Nonlinear Kinematic Hardening[J]. Production Engineering, 2020, 14(8): 135-145.

[15] SOHEIL Solhjoo, SHAHIN Khoddam. Evaluation of Barreling and Friction in Uniaxial Compression Test: A Kinematic Analysis[J]. International Journal of Mechanical Sciences, 2019, 156: 486-493.

[16] UTEN Khanawapee, SUTHEP Butdee. A Study of Barreling and DEFORM 3D Simulation in Cold Upsetting of Bi-Material[J]. Materials Today: Proceedings, 2020, 26(2): 1262-1270.

[17] FRAS T, RUSINEK A, PECHERSKI R B, et al. Analysis of Friction Influence on Material Deformation under Biaxial Compression State[J]. Tribology International, 2014, 80: 14-24.

[18] PERSSON B, SIVEBAEK I M, SAMOILOV V N, et al. On the Origin of Amonton's Friction Law[J]. Journal of Physics Condensed Matter, 2008, 20(39): 1-11.

[19] 秦芳, 李子亮. TC18鈦合金機械支座模鍛成形特性的研究[J]. 熱加工工藝, 2022, 51(9): 110-113.

QIN Fang, LI Zi-liang. Research on Die Forging Characteristics of TC18 Titanium Alloy Mechanical Support[J]. Hot Working Technology, 2022, 51(9): 110-113.

[20] 孫世仁, 劉虹, 陳文琳, 等. 基于有限元軟件的鍛造工藝參數對牽引拉桿成形的影響分析[J]. 熱加工工藝, 2020, 49(15): 68-72.

SUN Shi-ren, LIU Hong, CHEN Wen-lin, et al. Analysis of the Influence of Forging Process Parameters on the Forming of Traction Rod Based on Finite Element Software[J]. Hot Working Technology, 2020, 49(15): 68-72.

[21] 張慶旭, 劉飛, 周偉. 基于Deform的鋁合金盤體件鍛造工藝參數對成形性能影響的分析[J]. 熱加工工藝, 2020, 49(13): 112-114.

ZHANG Qing-xu, LIU Fei, ZHOU Wei. Analysis of the Effect of Forging Process Parameters on the Formability of Aluminum Alloy Disc Body Parts Based on Deform[J]. Hot Working Technology, 2020, 49(13): 112-114.

[22] YE Xian-jue, GONG Xiao-juan, YANG Biao-biao, et al. Deformation Inhomogeneity Due to Sample-Anvil Friction in Cylindrical Compression Test[J]. Transactions of Nonferrous Metals Society of China, 2019, 29(2): 279-286.

[23] 郭宇飛, 趙康, 海永清. 面向有限元分析的三角網格布爾運算方法[J]. 系統仿真學報, 2022, 34(5): 1003-1014.

GUO Yu-fei, ZHAO Kang, HAI Yong-qing. A Triangular Mesh Boolean Operation Method for Finite Element Analysis[J]. Journal of System Simulation, 2022, 34(5): 1003-1014.

[24] JIANG Xiao-tong, PENG Qing-jin, CHENG Xiao-sheng, et al. Efficient Booleans Algorithms for Triangulated Meshes of Geometric Modeling[J]. Computer-Aided Design and Applications, 2016, 13(4): 419-430.

[25] DONNIE Curington. Real Time Data Verification of Load and Pressure Testing Using von Mises Yield Criterion for Thick Walled Tubular Products[J]. International Journal of Mechanical Engineering and Robotics Research, 2019, 8(2): 167-172.

Optimization of Cold Finishing Process for the End of Automobile Knuckle Arm

WAN Jia-xiong1,WU Jian-xiang2,YAN Yang2,DING Hua-feng1,2

(1. Hubei Key Laboratory of Power System Design and Test for Electrical Vehicle, Hubei Xiangyang 441053, China; 2. Hubei Sanhuan Forging Co., Ltd., Hubei Xiangyang 441700, China)

The work aims to obtain the relationship between the change rule of the one-way compression cold finishing and the radius and depth of the pre-deformed hole and between the friction coefficient and the rebound. The upper and lower surfaces of the forgings at the end of the knuckle arm were treated with pre-deformed grooves, and the elastoplastic finite element model of the forgings with unidirectional compression and cold finishing was established. The ABAQUS finite element software was used for numerical simulation to analyze the effects of the friction coefficient on the surface quality. The optimized process parameters were obtained. The accuracy of the simulation result was verified by experiment. The results showed that the greater the friction coefficient, the more obvious the drum shape of the knuckle arm end forging after cold finishing. When the finishing amount was 1 mm, the radius of the pre-deformed groove hole was 4 mm and the depth was 0.3 mm, the unidirectional compression was 1.1 mm and the friction coefficient was 0.2, the best surface quality can be obtained. When the friction coefficient is the same, the larger the radius of the pre-deformed groove hole, the smaller the radial displacement, and the rebound remains basically unchanged.

automobile knuckle arm; cold finishing; process optimization; rebound

10.3969/j.issn.1674-6457.2022.09.009

TG316;U463.46

A

1674-6457(2022)09-0066-07

2021–12–22

湖北省教育廳科學技術研究計劃青年人才項目(Q20202602)

宛加雄(1996—),男,碩士生,主要研究方向為汽車輕量化及CAE分析。

丁華鋒(1986—),男,博士,副教授,主要研究方向為汽車輕量化設計及制造技術。

責任編輯:蔣紅晨

猜你喜歡
變形工藝
談詩的變形
中華詩詞(2020年1期)2020-09-21 09:24:52
轉爐高效復合吹煉工藝的開發與應用
山東冶金(2019年6期)2020-01-06 07:45:54
5-氯-1-茚酮合成工藝改進
世界農藥(2019年2期)2019-07-13 05:55:12
“我”的變形計
變形巧算
例談拼圖與整式變形
會變形的餅
一段鋅氧壓浸出與焙燒浸出工藝的比較
銅業工程(2015年4期)2015-12-29 02:48:39
FINEX工藝與高爐工藝的比較
新疆鋼鐵(2015年3期)2015-11-08 01:59:52
絡合鐵脫硫工藝在CK1井的應用
主站蜘蛛池模板: 欧美翘臀一区二区三区| 亚洲国产精品一区二区高清无码久久 | 国产区免费精品视频| 怡春院欧美一区二区三区免费| 91免费观看视频| 久热中文字幕在线观看| 婷婷亚洲天堂| 免费人成视网站在线不卡 | 国产精品污污在线观看网站| 影音先锋丝袜制服| 亚洲无码91视频| 99久久免费精品特色大片| 久久免费精品琪琪| 亚洲欧洲日韩综合色天使| 2020最新国产精品视频| 无码高清专区| 欧美一道本| 青青青视频蜜桃一区二区| av一区二区三区高清久久| 找国产毛片看| 久久人搡人人玩人妻精品| 黑色丝袜高跟国产在线91| 呦女亚洲一区精品| 波多野结衣亚洲一区| 不卡午夜视频| 露脸国产精品自产在线播| 国产毛片一区| 日韩毛片在线播放| 国产视频自拍一区| 久久一色本道亚洲| yjizz国产在线视频网| 蜜臀AV在线播放| Aⅴ无码专区在线观看| 国产91无码福利在线| 青青青亚洲精品国产| 国产精品19p| 国产经典三级在线| 久久久四虎成人永久免费网站| 99偷拍视频精品一区二区| 中文字幕在线播放不卡| 婷婷亚洲最大| 好紧好深好大乳无码中文字幕| 最新国产高清在线| 91国语视频| 欧美精品v| 97久久精品人人| 69av在线| 免费毛片全部不收费的| 免费一级无码在线网站 | 欧美精品在线视频观看| 人妻无码中文字幕一区二区三区| 四虎影视8848永久精品| 夜夜操狠狠操| 综合色88| 欧美日在线观看| 午夜国产精品视频| 亚洲中文字幕精品| 国产丝袜丝视频在线观看| 六月婷婷激情综合| 国产日韩欧美精品区性色| 国产亚洲欧美在线中文bt天堂| 综合久久五月天| 小13箩利洗澡无码视频免费网站| 99re66精品视频在线观看| 亚洲aaa视频| 四虎影院国产| 国产精品免费p区| 最新国产在线| 国产另类视频| 亚洲视频免费在线看| 久无码久无码av无码| 人妻中文字幕无码久久一区| 操美女免费网站| 午夜日b视频| 97国产精品视频人人做人人爱| 最新国产高清在线| 午夜日b视频| 国产一区免费在线观看| 国产亚洲精品97AA片在线播放| 亚洲国产中文欧美在线人成大黄瓜| 国产福利小视频在线播放观看| 亚洲欧美成人影院|