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

冷整形對(duì)大模數(shù)直齒輪精度的影響

2017-03-22 02:56:47梁強(qiáng)周杰朱姍姍曹金豆
關(guān)鍵詞:工藝模型

梁強(qiáng),周杰,朱姍姍,曹金豆

?

冷整形對(duì)大模數(shù)直齒輪精度的影響

梁強(qiáng),周杰,朱姍姍,曹金豆

(重慶大學(xué)材料科學(xué)與工程學(xué)院,重慶,400044)

為提高大模數(shù)圓柱直齒輪冷整形后的精度,以實(shí)測齒輪預(yù)鍛件形狀尺寸為基礎(chǔ),建立基于誤差模型的直齒輪冷整形過程的彈塑性有限元仿真模型,通過數(shù)值模擬得到精整量和精整凹模定徑帶長度對(duì)齒輪齒廓總偏差和螺旋線總偏差的影響規(guī)律。對(duì)直齒輪冷整形工藝進(jìn)行實(shí)際驗(yàn)證,并對(duì)整形后齒輪的檢測數(shù)據(jù)進(jìn)行方差分析。研究結(jié)果表明:精整量對(duì)齒輪精度影響顯著;當(dāng)精整量為0.15~0.20 mm時(shí),通過冷整形可以獲得齒形精度為8級(jí)、齒向精度為7~8級(jí)的齒輪;工藝實(shí)驗(yàn)得到的齒輪偏差與模擬結(jié)果變化趨勢一致。

大模數(shù)直齒輪;冷整形;齒輪精度;數(shù)值模擬

隨著齒輪成形工藝研究的深入和市場需求的不斷增大,現(xiàn)代制造業(yè)對(duì)齒輪產(chǎn)品精度要求不斷提高,齒輪的精密成形一直人們的研究熱點(diǎn)[1?2]。大模數(shù)直齒輪作為大型重載機(jī)械中重要的零部件,其強(qiáng)度與精度成為加工制造的關(guān)鍵。目前,大模數(shù)直齒輪精密成形工藝主要有“溫鍛+冷整形”復(fù)合工藝和“冷擠壓?冷整形”復(fù)合工藝。溫冷復(fù)合工藝具有成形力小、模具壽命長的優(yōu)點(diǎn),但溫鍛齒坯后需要中間幾道工序處理后才能冷整形,工藝較冗雜。隨著新型高強(qiáng)度冷作模具鋼的應(yīng)用,冷擠壓逐漸成為齒輪精密成形的主要工藝,擠齒后齒面光潔,經(jīng)潤滑可以直接冷整形,精整后精度高,表面質(zhì)量好。在冷擠壓?冷整形過程中,冷整形是保證成形精度的關(guān)鍵工藝。朱懷沈等[3?5]對(duì)溫冷復(fù)合工藝下精整量對(duì)齒輪精度的影響進(jìn)行了研究,但均采用無誤差模型進(jìn)行精整模擬及實(shí)驗(yàn)。STONE 等[6?7]對(duì)冷整形后單齒表面質(zhì)量及齒形尺寸進(jìn)行了研究,但未考慮精整量對(duì)齒輪精度的影響。胡成亮等[8?10]則對(duì)直齒輪1次冷精鍛成形進(jìn)行了數(shù)值模擬和實(shí)驗(yàn)研究,鍛后齒輪填充飽滿,但并未對(duì)精度進(jìn)行分析。為此,本文作者采用彈塑性有限元法,以模數(shù)=4,齒數(shù)=16,壓力角=20°,變位系數(shù)=0.45的齒輪為研究對(duì)象,建立直齒輪冷整形有限元模型,研究精整量及精整凹模定徑帶長度對(duì)冷整形后齒輪精度的影響;同時(shí),通過工藝實(shí)驗(yàn)進(jìn)行驗(yàn)證,得出直齒輪冷整形最佳的精整量和精整凹模定徑帶長度。

1 工藝方案

直齒輪成形工藝為熱鍛制坯—球化退火—拋丸—車削制坯—磷化皂化—冷擠壓—冷整形,材料為20CrMnTi。冷擠壓和冷整形均采用通過式的成形工藝,如圖1所示。

冷擠壓后的齒輪如圖2所示,齒輪上端與下端齒頂部分(見圖2(b))均出現(xiàn)填充不飽滿的現(xiàn)象,上端不飽滿長度為1.5~2.0 mm,上端端面(見圖2(a))凹陷約2.0 mm,并且存在縱向毛刺;下端不飽滿長度為2.0~3.0 mm下端端面(見圖2(c))外凸約1.5 mm。導(dǎo)致冷擠壓齒輪上、下端充填不飽滿的原因主要是擠壓時(shí)上、下端均處于無約束狀態(tài),金屬填充模具能力弱。由于在設(shè)計(jì)擠壓凹模時(shí)對(duì)齒形進(jìn)行修形,故擠壓件車外圓后齒頂可滿足設(shè)計(jì)尺寸要求。

在齒輪測量中心檢測冷擠壓齒輪的精度,齒輪精度檢測的2個(gè)重要指標(biāo)為齒廓總偏差和螺旋線總偏差[11],測量結(jié)果見表1。冷擠壓齒輪的精度較低,未能達(dá)到精鍛齒輪設(shè)計(jì)的精度要求,只能作為預(yù)鍛件,故需要增加冷整形工藝來提高冷擠壓齒輪的精度。

(a) 冷擠壓工藝示意圖;(b) 冷整形工藝示意圖

表1 齒輪精度檢測報(bào)告

(a) 上端;(b) 正面;(c) 下端

2 有限元模型建立及工藝實(shí)驗(yàn)

2.1 有限元模型建立

朱懷忱等[3?5]建立的整形模型均為無誤差的理想模型,即整形前齒輪的齒廓總偏差和齒形總偏差均為0 μm,這與實(shí)際情況并不吻合。為了建立與實(shí)際情況相吻合的模型,需要測量擠壓后齒輪的尺寸。隨機(jī)選擇3個(gè)擠壓后的齒輪,使用公法線千分尺跨3齒測量5組沿齒寬不同截面的公法線長度,每個(gè)截面測量3組不同的齒并求平均值,換算后得到沿齒寬方向齒厚度的變化。圖3所示為冷擠壓齒輪沿齒寬方向齒厚度的變化,齒輪中部齒厚度較上、下端均小,下端即先擠出端齒厚度最大。由于齒根部分不影響齒輪精度,齒頂后續(xù)要車外削加工,故冷整形只精整齒面,齒根和齒頂部分不進(jìn)行精整。

以實(shí)測齒厚分別建立三維模型和有限元模型。由于冷整形是小變形,故齒輪預(yù)鍛件定義為彈塑性體,材質(zhì)為20CrMnTi;為提高計(jì)算精度,采用六面體單元對(duì)齒輪預(yù)鍛件進(jìn)行網(wǎng)格劃分,并對(duì)齒形部分進(jìn)行網(wǎng)格局部細(xì)化;模具定義為剛體,忽略模具的彈性變形;設(shè)備選用液壓機(jī),凸模運(yùn)動(dòng)速度設(shè)為20 mm/s;預(yù)鍛件整形前需經(jīng)磷化皂化潤滑處理,故摩擦因數(shù)選擇0.12,服從剪切摩擦模型;為了提高計(jì)算效率,取單齒進(jìn)行模擬,有限元模型如圖4所示。為研究精整量和精整凹模定徑帶長度對(duì)齒輪精度的影響規(guī)律,以實(shí)測的齒輪預(yù)鍛件平均齒厚度作為基準(zhǔn)設(shè)計(jì)整形模具,選擇不同的精整量和定徑帶長度進(jìn)行數(shù)值模擬,模擬實(shí)驗(yàn)方案如表2所示。對(duì)模擬結(jié)果進(jìn)行量化處理,處理方法如下:在漸開線齒面上選取所有節(jié)點(diǎn)進(jìn)行跟蹤,并獲得整形后節(jié)點(diǎn)的坐標(biāo)信息,導(dǎo)入CATIA中進(jìn)行齒面重建,測量其螺旋線總偏差和齒廓總偏差。

圖3 齒輪預(yù)鍛件沿齒寬方向齒厚的變化

圖4 有限元模型

表2 模擬實(shí)驗(yàn)方案

2.2 工藝實(shí)驗(yàn)

工藝實(shí)驗(yàn)在四柱液壓機(jī)上進(jìn)行,精整凹模采用組合式凹模。為節(jié)約實(shí)驗(yàn)成本,先加工單邊精整量為 0.3 mm的整形凹模進(jìn)行實(shí)驗(yàn),整形凹模定徑帶長度為20 mm,實(shí)驗(yàn)完成后在線切割機(jī)床上對(duì)模具精整工作帶進(jìn)行慢走絲加工,加工至單邊精整量分別為0.25,0.20,0.15和0.10 mm后再次進(jìn)行整形實(shí)驗(yàn),整形凹模如圖5所示。由于不同定徑帶長度的凹模加工費(fèi)用高,故未進(jìn)行工藝實(shí)驗(yàn)。

(a) 整體;(b) 細(xì)節(jié)

利用齒輪測量中心對(duì)不同精整量整形后的齒輪進(jìn)行檢測,對(duì)每個(gè)精整量的齒輪分別抽選3個(gè)進(jìn)行測量,每次測量隨機(jī)取3個(gè)齒,并對(duì)測量結(jié)果進(jìn)行方差分析。

3 結(jié)果比較與分析

3.1 精整量對(duì)齒輪精度的影響規(guī)律

精整量是冷整形工藝的重要工藝參數(shù)之一,直接影響整形后齒輪的精度。表3所示為實(shí)測齒廓總偏差和實(shí)測螺旋線總偏差隨精整量變化的方差齊性檢驗(yàn)結(jié)果,可見當(dāng)前自由度下對(duì)應(yīng)的顯著性分別為0.220和0.062,均大于0.050,可以認(rèn)為樣本來自的總體滿足方差齊性的要求[12?13]。表4所示為實(shí)測齒廓總偏差和實(shí)測螺旋線總偏差隨精整量變化的方差分析結(jié)果,顯著性均小于0.001,說明精整量對(duì)齒廓總偏差和螺旋線總偏差有顯著性影響。

表3 方差齊性檢驗(yàn)結(jié)果

由有限元計(jì)算和實(shí)驗(yàn)所得到的齒廓總偏差與螺旋線總偏差隨精整量的變化分布見圖6。從圖6可見:與數(shù)值模擬結(jié)果相比,實(shí)測齒輪偏差偏大,但變化趨勢相吻合,誤差主要源于工藝實(shí)驗(yàn)中模具彈性變形引起的尺寸誤差。

圖6(a)所示為齒廓總偏差隨精整量的變化,可見齒廓總偏差總體上隨著精整量的增大而增大,精整量在0.1~0.2 mm間齒廓總偏差穩(wěn)定,精整量大于0.2 mm后齒廓總偏差變大。圖6(b)所示為螺旋線總偏差隨精整量變化。從圖6(b)可見:螺旋線總偏差隨著精整量的增加先減小后增大,精整量在0.15~0.20 mm區(qū)間內(nèi)冷整形后齒向精度最高;當(dāng)采用較小的精整量整形齒件時(shí),齒件經(jīng)過整形模具發(fā)生微小塑性變形甚至只有彈性變形,當(dāng)模具的約束消失后便發(fā)生彈性回復(fù),齒面不能得到有效整形,故過小的精整量并不能有效提高齒輪精度;當(dāng)采用過大精整量整形齒件時(shí),變形區(qū)域變大,齒面發(fā)生劇烈的塑性變形,達(dá)不到齒形精整的目的。綜合以上分析,齒件單邊精整量為0.15~ 0.20 mm時(shí)為宜。

(a) 齒廓總偏差;(b) 螺旋線總偏差

表4 方差分析結(jié)果

注:括號(hào)內(nèi)數(shù)據(jù)為組間值;括號(hào)外數(shù)據(jù)為組內(nèi)值。

3.2 定徑帶長度對(duì)齒輪精度的影響規(guī)律

凹模定徑帶是影響齒輪成形質(zhì)量的重要工藝參數(shù)之一[14?15]。圖7(a)所示為數(shù)值模擬計(jì)算的齒廓總偏差隨定徑帶長度的變化,可見齒廓總偏差隨著定徑帶長度的增大而減小,定徑帶長度超過20 mm后齒廓總偏差穩(wěn)定在0.02 mm。圖7(b)所示為數(shù)值模擬計(jì)算的螺旋線總偏差隨定徑帶長度的變化,螺旋線總偏差隨著定徑帶長度的增大而減小,定徑帶長度超過20 mm后螺旋線總偏差穩(wěn)定在0.006 mm。

(a) 齒廓總偏差;(b) 螺旋線總偏差

不同定徑帶長度的等效應(yīng)力場見圖8。從圖8可見:當(dāng)定徑帶長度為5 mm時(shí),齒面塑變區(qū)超過定徑帶長度(見圖8(a)),齒件通過模具后仍會(huì)發(fā)生變形,故整形后齒輪精度較低;當(dāng)定徑帶長度分別為20 mm和25 mm時(shí),齒面塑變區(qū)在定徑帶長度范圍內(nèi)(見圖8(b)和8(c)),齒面塑性變形后在定徑帶內(nèi)只進(jìn)行剛性移動(dòng)。故當(dāng)定徑帶達(dá)到20 mm后,齒輪精度便達(dá)穩(wěn)定值。

3.3 冷整形實(shí)驗(yàn)結(jié)果

圖9(a)所示為采用精整量0.2 mm和定徑帶長度20 mm冷整形后的實(shí)驗(yàn)結(jié)果,冷整形后齒件齒寬為61~63 mm。齒件經(jīng)車兩端面、鏜內(nèi)孔、車外圓和倒棱后,機(jī)加后齒件齒寬56.40~56.45 mm,齒形充填飽滿,齒形精度和齒向精度均達(dá)到8級(jí),較冷整形前齒輪精度顯著提高,如圖9(b)所示。

定徑帶長度/mm:(a) 5;(b) 20;(c) 25

(a) 冷整形后齒件;(b) 機(jī)械加工后齒件

4 結(jié)論

1) 以實(shí)際測量的冷擠壓齒輪數(shù)據(jù)建立誤差模型,數(shù)值模擬得到的精整量與齒輪精度間的規(guī)律與工藝實(shí)驗(yàn)結(jié)果趨勢一致,模擬結(jié)果與實(shí)驗(yàn)結(jié)果雖然在數(shù)值上有些差異,但總體來看能反映精整量與成形精度間的規(guī)律,說明建立的基于實(shí)測誤差模型的直齒輪預(yù)鍛件冷整形彈塑性有限元模型是正確的。

2) 齒廓總偏差隨精整量的增大而增大,螺旋線總偏差隨精整量的增大先減小后增大,精整量合理范圍為0.15~0.20 mm。采用精整量為0.2 mm,凹模定徑帶長度20 mm的工藝參數(shù)進(jìn)行工藝實(shí)驗(yàn),冷整形后能夠獲得齒形精度為8級(jí)、齒向精度為7~8級(jí)的齒輪。

3) 齒廓總偏差和螺旋線總偏差隨凹模定徑帶長度增加而減小,并在定徑帶長度大于20 mm時(shí)達(dá)到穩(wěn)定值。

[1] 馮再新, 周賢賓, 張治民, 等. 齒輪精密成形精整毛坯設(shè)計(jì)[J]. 塑性工程學(xué)報(bào), 2006, 13(1): 67?70. FENG Zaixin, ZHOU Xianbin, ZHANG Zhimin, et al. Designing of semifinished product for spur gear precision forming[J]. Journal of Plasticity Engineering, 2006, 13(1): 67?70.

[2] 周瑞, 陳文琳, 李亨, 等. 基于彈性變形行為的冷精整直齒輪尺寸預(yù)測[J]. 塑性工程學(xué)報(bào), 2013, 20(5): 27?32. ZHOU Rui, CHEN Wenlin, LI Heng, et al. Demension prediction of spur gear in cold precision sizing based on elastic deformation behavior[J]. Journal of Plasticity Engineering, 2013, 20(5): 27?32.

[3] 朱懷沈, 夏巨諶, 金俊松, 等. 大模數(shù)直齒輪溫冷鍛精整量的優(yōu)化選擇[J]. 塑性工程學(xué)報(bào), 2011, 18(1): 53?57. ZHU Huaishen, XIA Juchen, JIN Junsong, et al. Optimization of sizing amount during warm-cold compound forging for spur gear with big modulus[J]. Journal of Plasticity Engineering, 2011, 18(1): 53?57.

[4] 謝晉市. 大模數(shù)直齒輪溫冷復(fù)合精密成形工藝研究[D]. 合肥: 合肥工業(yè)大學(xué)材料科學(xué)與工程學(xué)院, 2012: 1?67. XIE Jinshi. Research on the warm-cold composite precision forming process of large module spur gear[D]. Hefei: Hefei University of Technology. School of Materials Science and Engineering, 2012: 1?67.

[5] 張寶紅. 高性能直齒圓柱齒輪開放成形研究[D]. 太原: 中北大學(xué)材料科學(xué)與工程學(xué)院, 2005: 12?17. ZHANG Baohong. Study on opened-die forming for high property spur gear[D]. Taiyuan: North University of China. School of Materials Science and Engineering, 2005: 12?17.

[6] STONE E R H, CAI J, HU Z M. An exercise in cold ironing as the post-forging operation for net-shape manufactures[J]. Journal of Materials Processing Technology, 2002, 135(2): 278?283.

[7] CHANG Y C, HU Z M, KANG B S. A study of cold ironing as a post-process for net-shape manufactures[J]. International Journal of Machine Tools & Manufacture, 2002, 42(8): 945?952.

[8] 胡成亮, 施衛(wèi)兵, 徐翔龍, 等. 帶轂直齒輪冷精鍛新工藝及其數(shù)值模擬[J]. 上海交通大學(xué)學(xué)報(bào), 2009, 43(9): 1494?1497. HU Chengliang, SHI Weibing, XU Xianglong, et al. A novel cold precision forging process of spur gear with boss and its numerical simulation[J]. Journal of Shanghai Jiaotong University, 2009, 43(9): 1494?1497.

[9] CHOI J C, CHOI Y. Precision forging of spur gears with inside relief[J]. International Journal of Machine Tools & Manufacture, 1999, 39(10): 1575?1588.

[10] EYRCIOGLU O, KUTUK M A, YILMAZ N F. Shrink fit design for precision gear forging dies[J]. Journal of Materials Processing Technology, 2009, 209(4): 2186?2194.

[11] 朱孝錄. 齒輪傳動(dòng)設(shè)計(jì)手冊(cè)[M]. 北京: 化學(xué)工業(yè)出版社, 2010: 234?265. ZHU Xiaolu. Handbook of gear design[M]. Beijing: Chemical Industry Press, 2010: 234?265.

[12] 楊虎, 劉瓊蓀, 鐘波. 數(shù)理統(tǒng)計(jì)[M]. 北京: 高等教育出版社, 2004: 153?154. YANG Hu, LIU Qiongsun, ZHONG Bo. Mathematical statistics[M]. Beijing: Higher Education Press, 2004: 153?154.

[13] 劉文卿. 實(shí)驗(yàn)設(shè)計(jì)[M]. 北京: 清華大學(xué)出版社, 2005: 23?26. LIU Wenqing. Design of experiments[M]. Beijing: Tsinghua University Press, 2005: 23?26.

[14] 馮文杰, 辜志勇, 陳瑩瑩. 定徑帶長度對(duì)花鍵軸冷擠壓彎曲變形的影響[J]. 現(xiàn)代制造工程, 2011(10): 97?99. FENG Wenjie, GU Zhiyong, CHEN Yingying. The influence of die bearing length on bending deformation of spline shaft in cold extrusion[J]. Modern Manufacturing Engineering, 2011(10): 97?99.

[15] 辜志勇. 花鍵軸冷擠壓彎曲變形研究[D]. 重慶: 重慶理工大學(xué)重慶汽車學(xué)院, 2012: 51?74. GU Zhiyong. Research on bending deformation of spline shaft in cold extrusion[D]. Chongqing: Chongqing University of Technology. College of Chongqing Automobile, 2012: 51?74.

(編輯 陳燦華)

Influence of gear accuracy in cold precision sizing for spur gear with big modulus

LIANG Qiang, ZHOU Jie, ZHU Shanshan, CAO Jindou

(College of Material Science and Engineering, Chongqing University, Chongqing 400044, China)

In order to improve gear accuracy in cold precision after cold precision sizing, an elastic plastic FE model based on the measured gear shape of pre-forging was built to simulate the cold precision sizing process. The laws of sizing amount and die bearing length on the tooth profile total deviation and the tooth orientation total deviation were obtained by numerical simulation. The experiment of the spur gear cold precision sizing process was made and data were analyzed by analysis of variance (ANOVA). The results show that the sizing amount has significant effects on gear accuracy. When the sizing amount is 0.15?0.20 mm, the gear of the tooth profile accuracy with 8 magnitudes and the gear direction precision with 7?8 magnitudes are obtained through the cold precision sizing. The change tendency of the deviation of gear obtained by experiment is similar to the FEM result.

spur gear with big modulus; cold precision sizing; gear accuracy; numerical simulation

10.11817/j.issn.1672?7207.2017.02.006

TG316

A

1672?7207(2017)02?0310?06

2016?03?05;

2016?05?22

國家自然科學(xué)基金面上項(xiàng)目(51275543);重慶市科委應(yīng)用開發(fā)項(xiàng)目(cstc2014yykfC70003)(Project(51275543) supported by the National Natural Science Foundation of China; Project(cstc2014yykfC70003) supported by the Science and Technology Application Development of Chongqing Municipal)

周杰,博士,教授,從事體積精密成形研究;E-mail:woliangniliu@163.com

猜你喜歡
工藝模型
一半模型
重要模型『一線三等角』
重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
轉(zhuǎn)爐高效復(fù)合吹煉工藝的開發(fā)與應(yīng)用
山東冶金(2019年6期)2020-01-06 07:45:54
5-氯-1-茚酮合成工藝改進(jìn)
3D打印中的模型分割與打包
一段鋅氧壓浸出與焙燒浸出工藝的比較
FINEX工藝與高爐工藝的比較
新疆鋼鐵(2015年3期)2015-11-08 01:59:52
FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
絡(luò)合鐵脫硫工藝在CK1井的應(yīng)用
主站蜘蛛池模板: 欧美成人午夜影院| 欧美精品二区| 伊人五月丁香综合AⅤ| 日韩在线永久免费播放| 九九九九热精品视频| 日本在线国产| 国产老女人精品免费视频| 青青草综合网| 欧美成人综合视频| 性69交片免费看| 精品五夜婷香蕉国产线看观看| 免费中文字幕在在线不卡| 特级精品毛片免费观看| 国产96在线 | 亚洲国产黄色| 91区国产福利在线观看午夜 | 日韩精品无码免费一区二区三区 | 国产精品亚洲αv天堂无码| 91精选国产大片| 久久精品国产一区二区小说| 亚洲综合色在线| 五月天丁香婷婷综合久久| 亚洲日韩精品伊甸| 久久a毛片| 亚洲天堂777| 99在线观看精品视频| 国产成人艳妇AA视频在线| 91探花在线观看国产最新| 欧美成人综合在线| 亚洲伦理一区二区| 亚洲首页国产精品丝袜| 九色视频最新网址| 免费Aⅴ片在线观看蜜芽Tⅴ | 91综合色区亚洲熟妇p| 久久久黄色片| 亚洲资源在线视频| 亚洲swag精品自拍一区| 亚洲网综合| 91久久国产综合精品女同我| 亚洲成肉网| 中文字幕 欧美日韩| 日韩欧美中文| 黑人巨大精品欧美一区二区区| 国产乱子伦精品视频| 人妻一区二区三区无码精品一区| 91福利在线观看视频| 日本不卡在线播放| 欧美一级视频免费| 四虎永久在线视频| 色老二精品视频在线观看| 亚洲精品国产综合99| 亚洲AV无码不卡无码| 精品国产中文一级毛片在线看| 夜夜拍夜夜爽| 国产一级片网址| 国产黄网永久免费| 午夜啪啪网| 日本国产精品一区久久久| 亚洲人免费视频| 亚洲天堂网视频| 色综合久久久久8天国| 2021国产精品自产拍在线| 国产精选自拍| 扒开粉嫩的小缝隙喷白浆视频| 国产一区二区在线视频观看| a级毛片在线免费观看| 少妇露出福利视频| 欧美视频在线第一页| 黄色一及毛片| 国产91色在线| 大陆国产精品视频| 精品91视频| 又黄又湿又爽的视频| 在线观看免费黄色网址| 国产精品人成在线播放| 久久九九热视频| 亚洲成肉网| 91在线播放免费不卡无毒| 亚洲欧洲日产无码AV| 99视频全部免费| 亚洲区一区| 91精品专区|