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

        ?

        基于試驗對比的汽車輪轂動特性及疲勞壽命研究

        2021-12-21 10:24:28唐貴基喻自力王曉龍高會超
        關(guān)鍵詞:輪轂壽命模態(tài)

        唐貴基,喻自力,王曉龍,白 潔,高會超

        基于試驗對比的汽車輪轂動特性及疲勞壽命研究

        唐貴基1,喻自力1,王曉龍1,白 潔1,高會超2

        (1. 華北電力大學(xué) 機械工程系,河北 保定 071000;2. 保定市立中車輪制造有限公司,河北 保定 071000)

        為研究汽車輪轂承擔(dān)動態(tài)載荷下的性能,以鋁合金輪轂為研究對象,基于有限元建模仿真的方法對輪轂結(jié)構(gòu)的固有頻率及振型等固有特性進行分析,并通過現(xiàn)場模態(tài)試驗測試對有限元分析結(jié)果進行了驗證。此外,對輪轂做進一步靜力學(xué)仿真分析,在獲取輪轂受載狀態(tài)下應(yīng)力云圖基礎(chǔ)上,分析了輪轂的疲勞壽命及最易失效部位,并通過輪轂全壽命周期加速疲勞試驗對仿真分析結(jié)果進行驗證。通過有限元仿真分析結(jié)果與實際試驗結(jié)果的對比,驗證了所述建模仿真方法的準(zhǔn)確性及可靠性,為研究輪轂固有特性、提高疲勞壽命及結(jié)構(gòu)優(yōu)化設(shè)計提供理論依據(jù)。

        輪轂;固有特性;模態(tài)試驗;疲勞壽命分析;疲勞試驗

        輪轂作為汽車的重要安全部件,行駛時承受各種交變載荷,因此要求輪轂具有良好的動特性與足夠的疲勞壽命。隨著有限元計算與試驗設(shè)備、技術(shù)的不斷發(fā)展,可以更好地研究結(jié)構(gòu)動態(tài)性能,為結(jié)構(gòu)優(yōu)化提供依據(jù)。文獻[1]基于ANSYS建立參數(shù)化有限元模型,再通過試驗?zāi)B(tài)進行驗證,在確定模型的可靠性基礎(chǔ)上逐漸優(yōu)化實體結(jié)構(gòu)模型。文獻[2]基于有限元方法模擬物體被接觸擠壓中的內(nèi)部受力,從而分析內(nèi)部結(jié)構(gòu)損傷規(guī)律,為結(jié)構(gòu)設(shè)計與控制方法提供數(shù)據(jù)參考。文獻[3]基于ANSYS對結(jié)構(gòu)進行模態(tài)分析與預(yù)測,理論分析上提高動態(tài)性能,模型修改前后結(jié)構(gòu)動特性與模態(tài)特性發(fā)生變化。文獻[4]基于有限元理論計算某圓管結(jié)構(gòu)承受動態(tài)軸向載荷下的沖擊能量,修改的模型在實際試驗中較之前的結(jié)構(gòu)在能量吸收特性方面有明顯的優(yōu)化。文獻[5]通過處理實驗平臺下的加速度信號,將試驗獲取結(jié)構(gòu)的前四階固有頻率與振動模式與數(shù)值計算相驗證。文獻[6]通過有限元數(shù)值計算方法對于集成動力學(xué)模型的動態(tài)特性建立評估函數(shù),并基于評估函數(shù)簡化結(jié)構(gòu)的等效模型。文獻[7]將有限元分析引入傳統(tǒng)輪轂設(shè)計流程,增強產(chǎn)品造型創(chuàng)新與設(shè)計效率。為優(yōu)化輪轂結(jié)構(gòu)設(shè)計,分析輪轂結(jié)構(gòu)動特性與疲勞壽命性能,本文通過有限元仿真與試驗結(jié)合對比分析,獲得輪轂的動特性,驗證仿真分析的可行性;并依據(jù)疲勞壽命曲線進行疲勞壽命分析,再通過彎曲疲勞試驗發(fā)現(xiàn)疲勞裂紋處與仿真分析一致,驗證了該方法的準(zhǔn)確性,為結(jié)構(gòu)優(yōu)化提供理論與試驗依據(jù)。

        1 有限元模型

        1.1 計算原理

        有限元方法基本思想將結(jié)構(gòu)進行離散,將載荷離散分布到各節(jié)點,分別計算節(jié)點位移與應(yīng)力分布。為分析單元應(yīng)力特性,利用多種方程求解應(yīng)力與應(yīng)變關(guān)系,通過幾何方程建立單元體應(yīng)變矩陣:

        為分析結(jié)構(gòu)內(nèi)的小變形與線性彈性,獲得其單元體內(nèi)的應(yīng)力矩陣,根據(jù)物理方程進行變換:

        式中:[]為幾何矩陣,[]為彈性矩陣。

        通過虛功原理得出各單元體中的力:

        式中:[]為整體剛度矩陣。

        對于整體結(jié)構(gòu)上的任意一點,建立平衡方程:

        式中:{R}為節(jié)點上的外載荷。

        結(jié)合單元體受力情況與邊界條件,得出所結(jié)構(gòu)體的整體平衡方程:

        式中:[]為整體剛度矩陣,{}為整體結(jié)構(gòu)節(jié)點位移矩陣,由式(1)、式(2)求出各單元的應(yīng)變和應(yīng)力,式(5)求出節(jié)點位移。

        1.2 模型建立

        考慮輪轂的強度、質(zhì)量以及散熱等性能,輪轂?zāi)P筒馁|(zhì)為鋁合金A356,參數(shù)如表1所示。

        表1 輪轂材料參數(shù)

        材料類型彈性模量/MPa泊松比密度/(kg·m-3) 鋁合金710000.332770

        由于輪轂?zāi)P椭杏新菟?、氣門孔以及連接中心孔,因此通過Space Claim對于模型中的線以及圓面進行合并以及刪除處理,所建立的有限元模型如圖1所示。

        圖1 汽車輪轂?zāi)P?/p>

        1.3 網(wǎng)格劃分與約束施加

        基于有限元思想進行網(wǎng)格劃分以及模態(tài)分析,由于有限元模型網(wǎng)格質(zhì)量關(guān)系著是否能夠得到高精度、高吻合的分析結(jié)果[8],因此,對輪轂進行網(wǎng)格劃分時,全局網(wǎng)格控制規(guī)則采取Curvature,對于螺栓孔、中心孔附近曲率變化處網(wǎng)格自動加密[9],通過ANSYS workbench對模型進行前處理,模型采用四面體單元,獲得網(wǎng)格如圖2所示。

        圖2 網(wǎng)格劃分

        模型包括69 726個節(jié)點,37 238個單元??紤]到輪轂實際使用中與汽車車軸通過螺栓緊固連接,因此,對4個螺栓孔施加全約束進行模態(tài)分析。

        2 固有特性分析

        2.1 模態(tài)分析

        模態(tài)分析是一種確定結(jié)構(gòu)振動特性的技術(shù),是結(jié)構(gòu)運動學(xué)的基礎(chǔ)[10]。結(jié)構(gòu)設(shè)計過程中,利用有限元仿真進行模態(tài)計算,并據(jù)此進行優(yōu)化設(shè)計[11]。經(jīng)過有限元分析獲得固有頻率、模態(tài)振型等,來解決復(fù)雜的結(jié)構(gòu)振動問題、振動測量與結(jié)構(gòu)的動力學(xué)分析,表2為通過有限元仿真分析獲得的動特性參數(shù),對應(yīng)的四階模態(tài)振型如圖3所示。

        表2 模態(tài)分析結(jié)果

        階數(shù)模態(tài)頻率/Hz模態(tài)振型 1365.09輪輞彎曲 2787.4輪輞偏轉(zhuǎn) 31276.9三角形振動 41834.8四邊形彎曲

        圖3 四階振型圖

        固有特性分析結(jié)果通??捎糜诮Y(jié)構(gòu)優(yōu)化中的數(shù)據(jù)指導(dǎo),來盡量避免結(jié)構(gòu)出現(xiàn)共振導(dǎo)致功能失效。考慮汽車實際行駛中的路況,城市路況較好,平坦路面激勵頻率或者石子路等激勵頻率一般低于15 Hz,發(fā)動機的轉(zhuǎn)速對應(yīng)頻率如公式(6)所示。

        式中:為內(nèi)燃機轉(zhuǎn)速,r/min;為氣缸數(shù)目;為沖程數(shù)目,四行程2,兩行程1。

        一般發(fā)動機轉(zhuǎn)速可達8 000r/min,則0=267 Hz,汽車行業(yè)中的共振帶范圍一般取固有頻率的13%到20%的范圍,當(dāng)激勵頻率和結(jié)構(gòu)固有頻率滿足式(7)時,會發(fā)生共振。

        式中:為固有頻率;0為激勵頻率。

        輪轂一階固有頻率為365.09 Hz,若按照20%的頻率范圍,共振帶范圍為292~438 Hz,當(dāng)四缸發(fā)動機加速到最高轉(zhuǎn)速,對應(yīng)的激勵頻率僅為0=267 Hz,未達到輪轂的共振頻帶覆蓋范圍內(nèi),因此正常行駛狀態(tài)下輪轂不會出現(xiàn)共振現(xiàn)象,能夠保障車輛行駛的安全性。為驗證所構(gòu)建輪轂有限元模型分析結(jié)果的準(zhǔn)確性,進行輪轂?zāi)B(tài)測試試驗,用于對比有限元分析結(jié)果及現(xiàn)場試驗測試結(jié)果。

        2.2 模態(tài)試驗驗證

        模態(tài)試驗采用北京東方振動和噪聲研究所的成套設(shè)備對于輪轂進行分析,配套設(shè)備的具體型號如表3所示.

        表3 試驗儀器設(shè)備

        設(shè)備名稱型號 北京東方所數(shù)據(jù)采集系統(tǒng)INV3018A 模態(tài)試驗分析系統(tǒng)DASP-V10 加速度傳感器INV9824 模態(tài)力錘INV9310

        現(xiàn)場模態(tài)試驗設(shè)備如圖4所示,圖4中(a)為模態(tài)力錘,(b)為加速度傳感器,(c)為試驗輪轂,使用DASP-V10系統(tǒng)進行加速度數(shù)據(jù)的分析處理。

        圖4 模態(tài)試驗設(shè)備

        試驗中采用螺栓孔螺栓固定約束,對于四個螺栓孔施加120 N·m的扭矩進行固定;激勵方式為力錘激勵法,試驗共進行20次激勵,激勵點分布在傳感器附近的2點以及對稱的2點,以及輪轂背面與正面激勵點平行的4點;2傳感器安裝位置呈90°分布,在傳感器底座上涂抹熱熔膠,如圖(b)所示垂直安裝在輪緣上;在INV3018A采集系統(tǒng)中設(shè)置采樣頻率為20 kHz,激勵時輪轂正面4個激勵點各激勵3次,背面4個激勵點各激勵2次,為保證加速度傳感器的數(shù)據(jù)的準(zhǔn)確性,每次使用力錘快速敲擊激勵點并采樣,等采集系統(tǒng)中的數(shù)據(jù)波形完全振蕩衰減后再進行下1次激勵。

        為避免每次激勵力度不同可能對加速度傳感器中數(shù)據(jù)帶來的影響,在計算頻響函數(shù)時,加窗函數(shù)對采集到的波形進行截斷處理,通過模態(tài)試驗分析系統(tǒng)對20次加速度數(shù)據(jù)進行加窗函數(shù),并處理擬合得到如圖5所示頻響函數(shù)曲線圖。

        圖5 頻響函數(shù)曲線

        對圖5中的頻響函數(shù)圖進行標(biāo)定得出如表4所示輪轂前四階固有頻率,將之與仿真分析所得固有頻率進行對比,計算誤差。

        表4 有限元分析與試驗結(jié)果對比

        實測頻率/Hz分析頻率/Hz誤差/% 357.8365.092.0 813.4787.43.2 1358.11276.95.9 1980.91834.87.3

        通過表4中有限元仿真分析與模態(tài)試驗結(jié)果對比發(fā)現(xiàn),試驗分析中的頻響函數(shù)曲線可以明顯地得到輪轂的四階固有頻率,但是仿真分析中含有多階頻率,存在同階不同方向上的振型,且模態(tài)試驗中傳感器會受環(huán)境影響,但兩者結(jié)果對比得出誤差較小,說明仿真模型與實際模型誤差在允許范圍內(nèi)。

        3 靜力學(xué)分析

        3.1 載荷分析

        通過仿真分析與模態(tài)試驗的結(jié)果對比驗證了有限元模型的可靠性,再對承載實際載荷下的汽車輪轂進行靜力學(xué)分析,分析其動態(tài)特性。車輛行駛時力與力矩通過輪胎傳遞到輪輞,在傾斜路面行駛轉(zhuǎn)向時,車軸還將承受彎矩,因此依據(jù)動態(tài)彎曲疲勞試驗進行載荷的計算與施加。

        根據(jù)文獻[12],輪轂所受最大載荷可按照公式(8)確定。

        輪轂在汽車行駛中,支撐汽車載重的同時,承受彎曲載荷,根據(jù)GB/T 5334-2005《轎車鋼制車輪性能要求和試驗方法》中動態(tài)彎曲疲勞試驗要求,彎矩按照公式(9)計算。

        經(jīng)計算,=3 046 N·m。

        3.2 應(yīng)力分布

        基于如上計算得出輪轂所承受的載荷,再進行載荷施加及仿真動力學(xué)分析得到如圖6所示等效應(yīng)力圖,分析圖6可以得出應(yīng)力主要集中在輪輻與窗口上,其中最大等效應(yīng)力為117.65 MPa,分布在輪輻中間的窗口處,窗口附近的輪輻部位承受應(yīng)力也達到了104.58 MPa,由于鋁合金A356的屈服強度值為240 MPa,因此,應(yīng)力最大處有較高的安全系數(shù)和較高的強度儲備。

        圖6 等效應(yīng)力圖

        3.3 疲勞壽命分析

        承受長期交變載荷后,應(yīng)力最大處可能最先出現(xiàn)疲勞失效,因此需進行疲勞壽命分析。金屬構(gòu)件的壽命預(yù)測方法主要有基于應(yīng)力的壽命預(yù)測方法和基于應(yīng)變的壽命預(yù)測方法[13]。應(yīng)力疲勞壽命預(yù)測中,基于S-N曲線進行預(yù)測,工程中用Basquin方程表示橫幅載荷下疲勞壽命應(yīng)力的關(guān)系[14],如式(10)所示。

        零件或材料的疲勞曲線是零件疲勞壽命預(yù)測的基礎(chǔ),材料曲線一般通過試驗得到,也可以通過估算繪制出簡單的曲線[15],基于鋁合金A365的一些試驗參數(shù)繪制如圖7所示的疲勞壽命曲線。

        圖7 材料疲勞壽命曲線

        通過Fatigue Tool對輪轂進行疲勞壽命分析[16],得到圖8所示疲勞壽命云圖,輪轂疲勞壽命最小值出現(xiàn)在輪輻窗口處,最小疲勞壽命能承受2.343×106次周期的循環(huán)載荷。國家GB/T 5334—2005《乘用車車輪性能要求和試驗方法》中疲勞試驗要求輕合金車輪最低循環(huán)次數(shù)為5×105次,因此符合國家標(biāo)準(zhǔn)。

        圖8 疲勞壽命分析

        3.4 疲勞試驗驗證

        通過分析疲勞壽命云圖得出疲勞壽命最低處,為驗證其分析結(jié)果,使用CFT-3轎車車輪彎曲疲勞試驗機對輪轂進行彎曲疲勞試驗,設(shè)備如圖9所示。

        試驗中法蘭盤為4孔均勻分布法蘭盤,孔內(nèi)徑

        圖9 疲勞試驗臺

        為100 mm,使用120 N·m的扭矩將4個螺栓進行緊固,而后通過8個壓板將輪轂的下邊緣固定,試驗機所給彎矩為3 046 N·m,設(shè)置自動停機的主軸偏移量為初始的10%,主軸轉(zhuǎn)動直至主軸偏移量達到設(shè)定偏移比,得到如圖10所示疲勞失效裂紋。由于輪輻窗口屬于應(yīng)力最集中處,在承受交變載荷后最先出現(xiàn)疲勞失效,圖10試驗結(jié)果與仿真分析結(jié)果一致,通過有限元仿真分析結(jié)果與實際試驗結(jié)果的對比,驗證了所述建模仿真方法的準(zhǔn)確性,為研究同類汽車輪轂的固有特性及結(jié)構(gòu)優(yōu)化設(shè)計提供理論依據(jù)與技術(shù)支持。

        圖10 疲勞失效部位

        4 結(jié)論

        本文基于有限元仿真進行汽車輪轂的動特性及疲勞壽命分析,并進行模態(tài)試驗以及疲勞試驗得到輪轂的應(yīng)力最大以及疲勞壽命最小處對仿真結(jié)果進行驗證。為研究輪轂固有特性、提高強度及結(jié)構(gòu)優(yōu)化設(shè)計提供理論依據(jù)。

        (1)通過模態(tài)試驗的頻響函數(shù)曲線,計算輪轂一階固有頻率與仿真分析結(jié)果誤差較小,都避開一階共振頻率,不會產(chǎn)生共振。

        (2)進行靜力學(xué)分析,應(yīng)力集中在輪輻窗口處,應(yīng)力最高為117.65 MPa,該應(yīng)力小于材料的屈服強度240 MPa,具有較高的安全系數(shù)與強度儲備。

        (3)繪制疲勞壽命曲線,進行疲勞壽命分析,得出輪轂疲勞壽命最小值在輪輻窗口處,最小疲勞壽命能承受2.343×106次周期的循環(huán)載荷,完全滿足國家標(biāo)準(zhǔn),失效部位與疲勞試驗相符合。

        (4)將有限元分析與相關(guān)試驗對比分析,并依據(jù)疲勞壽命曲線進行分析與彎曲疲勞試驗驗證疲勞薄弱處,為結(jié)構(gòu)優(yōu)化提供理論基礎(chǔ)與依據(jù)。

        [1] 邱白晶, 何耀杰, 盛云輝, 等.噴霧機噴桿有限元模態(tài)分析與結(jié)構(gòu)優(yōu)化[J]. 農(nóng)業(yè)機械學(xué)報, 2014, 45(8): 112-116, 105.

        [2] 姬偉, 李俊樂, 楊俊, 等. 機器手采摘蘋果抓取損傷機理有限元分析及驗證[J]. 農(nóng)業(yè)工程學(xué)報, 2015, 31(5): 17-22.

        [3] 李天箭, 吳晨帆, 沈磊, 等. 基于模態(tài)預(yù)測及敏度分析的機床動特性設(shè)計方法[J]. 機械工程學(xué)報, 2019, 55(7): 178-186.

        [4] Ahmad M, Ismail K A, Hanid M H M, et al. Modification of the design of circular thin-walled tubes to enhance dynamic energy absorption characteristics: Experimental and finite element analysis[J]. IOP Conference Series: Materials Science and Engineering, 2020, 917(1): 1-13.

        [5] Fan W, Yang Y, Su X. Dynamic Modeling and Vibration Characteristics Analysis of Transmission Process for Dual-Motor Coupling Drive System[J]. Symmetry, 2020, 12(7): 1171-1190.

        [6] Xu M T, Cai B, Li C Y, et al. Dynamic characteristics and reliability analysis of ball screw feed system on a lathe[J]. Mechanism and Machine Theory, 2020, 150:1-21.

        [7] 閆勝昝, 劉偉鵬. 基于有限元分析的鋁合金輪轂造型設(shè)計[J]. 機械設(shè)計, 2014, 31(10):109-112.

        [8] 張穎, 王元清, 張俊光, 等. 鋁合金網(wǎng)殼結(jié)構(gòu)箱型-工字型盤式節(jié)點單肢受力性能有限元分析[J]. 工程力學(xué), 2020, 37(S1): 130-138.

        [9] 李海峰, 吳冀川, 劉建波, 等. 有限元網(wǎng)格剖分與網(wǎng)格質(zhì)量判定指標(biāo)[J]. 中國機械工程, 2012, 23(3): 368-377.

        [10] 代穎, 崔淑梅, 宋立偉. 車用電機的有限元模態(tài)分析[J]. 中國電機工程學(xué)報, 2011,31(9): 100-104.

        [11] 馬驪凕, 朱智民, 安占飛, 等. 轎車白車身模態(tài)分析[J]. 振動與沖擊, 2013, 32(21): 214-218.

        [12] 王明明. 鋁合金汽車輪轂結(jié)構(gòu)設(shè)計及優(yōu)化[D]. 長春: 吉林大學(xué), 2011.

        [13] 朱峰, 楊宏偉, 王本勁, 等. 航空導(dǎo)管彎曲試驗的疲勞裂紋萌生壽命[J]. 力學(xué)季刊, 2020, 41(3): 519-527.

        [14] 左芳君. 機械結(jié)構(gòu)的疲勞壽命預(yù)測與可靠性方法研究[D]. 成都: 電子科技大學(xué), 2016.

        [15] 黃寧. 大型結(jié)構(gòu)件的疲勞壽命預(yù)測方法研究[D]. 長沙: 中南大學(xué), 2013.

        [16] 沈明明, 李榮, 劉祖國, 等. 數(shù)控機床軸承結(jié)構(gòu)振動諧響應(yīng)及疲勞壽命研究[J]. 組合機床與自動化加工技術(shù), 2019(11): 143-147.

        Research On Dynamic Characteristics and Fatigue Life of Automobile Wheel Hub Based On Test Comparison

        TANG Gui-ji1, YU Zi-li1, WANG Xiao-long1, BAI Jie1, GAO Hui-chao2

        (1. Department of Mechanical Engineering, North China Electric Power University, Baoding 071000, China; 2. Baoding Lizhong Wheel Manufacturing Co., Ltd., Baoding, 071000, China)

        In order to study the performance of automobile wheel hub under dynamic load, aluminum alloy wheel hub taken as the research object, the natural frequency and vibration mode of wheel hub structure were analyzed based on the method of finite element modeling and simulation. The results of finite element analysis were verified by field modal test. In addition, a further static simulation analysis was performed on the hub, and on the basis of obtaining the stress cloud diagram of the hub under load, the fatigue life and the most prone failure position of the hub were analyzed. The simulation results were verified by accelerated fatigue test of the hub life cycle. The accuracy and reliability of the modeling and simulation method are verified by comparing the finite element simulation results with the actual test results, which provides a theoretical basis for studying the inherent characteristics of the hub, the fatigue life and structural optimization design improved.

        Wheel hub; inherent characteristics; modal test; fatigue life analysis; fatigue test

        10.15916/j.issn1674-3261.2021.06.001

        TH133

        A

        1674-3261(2021)06-0351-06

        2021-02-03

        國家自然科學(xué)基金項目(52005180);河北省自然科學(xué)基金項目(E2019502047);中央高?;究蒲袠I(yè)務(wù)費專項資金項目(2021MS069)

        唐貴基(1962-),男,山東龍口人,教授,博士,博士生導(dǎo)師。

        王曉龍(1989-),男,黑龍江大慶人,副教授,博士。

        責(zé)任編輯:陳 明

        猜你喜歡
        輪轂壽命模態(tài)
        抗壓痕透明粉在精車鋁輪轂上的應(yīng)用研究
        上海涂料(2021年5期)2022-01-15 06:09:26
        人類壽命極限應(yīng)在120~150歲之間
        中老年保健(2021年8期)2021-12-02 23:55:49
        倉鼠的壽命知多少
        馬烈光養(yǎng)生之悟 自靜其心延壽命
        華人時刊(2018年17期)2018-12-07 01:02:20
        基于CPS 的汽車輪轂制造系統(tǒng)設(shè)計
        人類正常壽命為175歲
        奧秘(2017年12期)2017-07-04 11:37:14
        基于ANSYS的輪轂支架結(jié)構(gòu)設(shè)計
        風(fēng)電輪轂鑄造技術(shù)研究
        大型鑄鍛件(2015年1期)2016-01-12 06:33:29
        國內(nèi)多模態(tài)教學(xué)研究回顧與展望
        基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識別
        91麻豆精品国产91久久麻豆| 成人无码视频在线观看网站| 国产综合一区二区三区av| 亚洲精品中文字幕一二三四| 国产播放隔着超薄丝袜进入| 在线亚洲人成电影网站色www| 中出高潮了中文字幕| 国产大屁股熟女流白浆一区二区 | 亚洲人成网站色在线入口口| 欧美真人性野外做爰| 99久久婷婷国产综合亚洲91| 中文字幕一区二区三区97| 在线观看视频国产一区二区三区 | 精品区2区3区4区产品乱码9| 又大又粗弄得我出好多水| 国产激情一区二区三区在线蜜臀 | 亚洲av永久青草无码精品| 久久亚洲网站中文字幕| 日韩人妻无码一区二区三区久久 | 亚洲熟妇无码av不卡在线播放| 国产av一区二区凹凸精品| 精品国产一区二区三区av免费| 久久久无码人妻精品无码| 日本亚洲欧美高清专区| WWW拍拍拍| 狠狠久久av一区二区三区| 风韵丰满熟妇啪啪区老老熟妇| 亚洲日韩v无码中文字幕| 亚洲AⅤ男人的天堂在线观看| 人妻av不卡一区二区三区| 亚洲综合网国产精品一区| 免费特级毛片| 国产一区二区欧美丝袜| 青青草免费在线视频久草| 久久久久av无码免费网| 亚洲无码一二专区| 精品国产精品久久一区免费| 国产a在亚洲线播放| 色多多a级毛片免费看| 精品国产亚洲一区二区三区演员表| 国产丝袜一区丝袜高跟美腿|