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

        ?

        基于損傷等效原理的商用車蓄電池箱耐久性分析

        2022-01-05 11:21:22任全陳文斐周肖飛杜嘉峰郭寧
        汽車實用技術(shù) 2021年23期
        關(guān)鍵詞:信號處理箱體蓄電池

        任全,陳文斐,周肖飛,杜嘉峰,郭寧

        基于損傷等效原理的商用車蓄電池箱耐久性分析

        任全,陳文斐,周肖飛,杜嘉峰,郭寧

        (陜西汽車控股集團有限公司,陜西 西安 710200)

        文章通過加速度傳感器采集蓄電池箱載荷信號,基于損傷等效原理對加速度載荷譜進行信號處理,輸出與原始信號損傷一致的PSD信號。建立蓄電池箱總成有限元模型,對其進行頻率響應(yīng)分析,得到輸入激勵與響應(yīng)應(yīng)力的傳遞函數(shù)。對蓄電池箱進行隨機振動疲勞仿真分析,得出箱體疲勞損傷結(jié)果,并與強化路試驗結(jié)果進行對比,驗證該方法的可行性。最后根據(jù)仿真結(jié)果對原方案進行結(jié)構(gòu)優(yōu)化,對新方案進行第二輪疲勞分析與強化路試驗驗證,并成功達到試驗?zāi)繕死锍?,從而完成蓄電池箱體方案整改。

        蓄電池箱;損傷等效;PSD信號;疲勞分析

        前言

        蓄電池是商用車供電系統(tǒng)的核心部件,作為電池組的載體,蓄電池箱對電池組起著支撐與防護作用,若其結(jié)構(gòu)失效將引起整車供電異常等危險情況,時刻影響著整車的安全性,因此在商用車正向研發(fā)過程中,對蓄電池箱體進行準確的疲勞壽命預(yù)測顯得尤為重要[1-2]。

        現(xiàn)階段對蓄電池箱體的耐久性開發(fā)主要通過隨機振動仿真及臺架振動試驗的方式進行。彭倩等結(jié)合《GB/T 31467— 2015 電動汽車用鋰離子動力蓄電池箱和系統(tǒng)》中的加速度功率譜PSD(Power Spectrum Density)對箱體進行了隨機振動疲勞分析,預(yù)估了箱體的疲勞損傷情況[3-4];吳長德等以定頻振動試驗臺的振動波形為載荷譜,結(jié)合Miner線性累積損傷法則計算得出了焊點的疲勞壽命[5]。以上文獻雖對箱體結(jié)構(gòu)進行了疲勞耐久計算,但其輸入載荷譜與實際路面載荷譜存在一定差異,未充分考慮實際路面激勵對箱體的損傷分布影響,從而導致了仿真結(jié)果與路試試驗結(jié)果對應(yīng)性不高。

        針對上述問題,本文基于損傷等效原理將路面時域信號轉(zhuǎn)換為加速度PSD信號,充分保證二者損傷的一致性,通過調(diào)節(jié)等效作用時間避免頻域加速現(xiàn)象,并以PSD隨機振動分析方法對蓄電池箱體進行疲勞計算,將仿真結(jié)果與強化路試驗情況進行對比,驗證該方法的正確性,同時根據(jù)仿真結(jié)果對箱體進行結(jié)構(gòu)優(yōu)化,并對新方案進行疲勞仿真與二次試驗驗證,最終新方案達到試驗?zāi)繕死锍虝r箱體無損壞,以此證明方案優(yōu)化成功,體現(xiàn)出基于損傷等效原理的耐久性分析在項目正向研發(fā)中的可行性與重要性,為后續(xù)此類產(chǎn)品的耐久性開發(fā)及驗證提供理論依據(jù)。

        1 理論背景

        1.1 沖擊響應(yīng)譜

        損傷等效理論是基于單自由度系統(tǒng)(SDOF),解決單一激勵振動部件的疲勞耐久問題,在汽車耐久性開發(fā)中有著較為廣泛的應(yīng)用[6]。

        沖擊響應(yīng)譜(Shock Response Spectrum, SRS)通常稱為“沖擊譜”,是當一系列線性的、相互獨立的單自由度系統(tǒng)受到?jīng)_擊激勵時,將各個單自由度系統(tǒng)的最大響應(yīng)值作為對應(yīng)于系統(tǒng)固有頻率的函數(shù)響應(yīng)曲線[7]。隨著系統(tǒng)固有頻率的不斷變化,輸入信號施加于不同固有頻率的單自由度系統(tǒng)后得到相應(yīng)的最大響應(yīng)值,可繪制出一條隨頻率變化的包絡(luò)線,該包絡(luò)線即為沖擊響應(yīng)譜曲線。沖擊響應(yīng)譜的生成流程如圖1所示。

        圖1 沖擊響應(yīng)譜生成流程圖

        1.2 疲勞損傷譜

        疲勞損傷譜(Fatigue Damage Spectrum, FDS)表示一系列單自由度系統(tǒng)在受到外部激勵時所產(chǎn)生疲勞損傷值與頻率變化的曲線,可反映外部激勵在某個頻率下對結(jié)構(gòu)產(chǎn)生損傷值的大小[8]。工程中的頻域加速試驗都是基于損傷等效原則進行的,故需對原始采集信號的疲勞損傷進行計算,其生成流程與沖擊響應(yīng)譜類似,如圖2所示。

        圖2 疲勞損傷譜生成流程圖

        疲勞損傷譜也可通過基于功率譜密度計算的方式進行獲取。在得到某原始時域信號的功率譜密度曲線(Power Spec- trum Density, PSD)后,通過式(1)可計算出某一頻率的損傷值,隨后對該系統(tǒng)下不同固有頻率的損傷值進行計算,將各頻率所對應(yīng)的損傷值用包絡(luò)線進行連接,即疲勞損傷譜,Lalanne給出了用于計算疲勞損傷譜的計算公式[9]。

        式中,(f)為頻率為f的疲勞損傷值;GZ?(f)為頻率為f的加速度PSD幅值;為單自由度系統(tǒng)彈簧剛度;、為疲勞曲線參數(shù);(1+?2)為(g)=0x (g?1)的Gamma函數(shù);為動態(tài)放大因子;為加速度等效時間,s;f為固有頻率,Hz。

        1.3 極端響應(yīng)譜

        極端響應(yīng)譜(Extreme Response Spectrum, ERS)表示持續(xù)時間為的PSD信號通過單自由度系統(tǒng)時所得到的最大概率的極限幅值所形成的譜圖[10]。在工程領(lǐng)域中,對于一條隨機振動的激勵曲線,一般更適合于用功率譜密度來表示數(shù)據(jù)振動的概率分布情況,因此通過極端響應(yīng)譜可考察該信號在不同頻率下的輸出響應(yīng)情況。極端響應(yīng)譜的計算公式為:

        式中,(f)為頻率為f的極端響應(yīng)譜幅值;GZ?(f)為頻率為f的加速度PSD幅值;為動態(tài)放大因子;為加速度等效時間;f為頻率,Hz。

        2.4 PSD信號生成方法

        基于上述理論,在得到激勵信號所對應(yīng)的總疲勞損傷譜后,可將其分解為、、三個方向的疲勞損傷譜,基于損傷譜可以分解為對應(yīng)方向的PSD激勵信號,并保證三個方向的PSD激勵信號損傷值之和與原始信號總損傷值一致。

        基于損傷等效理論的PSD信號獲取方法為:對疲勞損傷譜的推導公式(2)進行求逆,即可得到激勵時間為T的合成PSD激勵信號,該信號可用于疲勞耐久仿真計算。PSD信號的推導公式為:

        2 載荷譜采集及信號處理

        2.1 載荷譜采集

        載荷譜采集是指在路試試驗中,通過對行駛車輛的目標部位安裝各類傳感器,經(jīng)數(shù)據(jù)采集儀接收傳感器信號,得到試驗路面對目標部位的激勵信息,激勵信息通常包括位移、速度、加速度、力、力矩、應(yīng)變等[11]。本次試驗需用到的載荷譜采集設(shè)備包括:eDAQ數(shù)據(jù)采集器、Kistler加速度傳感器等,如圖3所示。

        圖3 數(shù)據(jù)采集儀與加速度傳感器

        試驗車輛為某款牽引車,樣車狀體為滿載狀態(tài),試驗環(huán)境為某試驗基地,其中強化路跑道包括坑洼路、比利時路、卵石路、井蓋路、魚鱗坑等。其中蓄電池箱固定于右側(cè)車架中部,為采集路面載荷對電瓶箱的振動信號,將加速度傳感器安裝在箱體所對應(yīng)的車架端,安裝位置如圖4所示。

        圖4 實車狀態(tài)下電瓶箱根部傳感器布置位置

        圖5 三向加速度載荷采集信號

        圖6 三向加速度載荷雨流計數(shù)圖

        通過對采集數(shù)據(jù)的初步檢查,判斷采集過程中試驗設(shè)備均正常工作,試驗采集的電池箱根部三向加速度信號如圖5所示,對其進行雨流計數(shù)統(tǒng)計可得,三向加速度信號均符合高斯分布,采集信號無明顯漂移現(xiàn)象,雨流統(tǒng)計譜如圖6所示,載荷數(shù)據(jù)滿足后續(xù)疲勞分析需求。

        2.2 信號處理

        2.2.1信號處理流程

        在得到試驗載荷譜后,需基于損傷等效原理將原始時域信號轉(zhuǎn)換為PSD信號,該方法既保證PSD信號損傷的一致性,避免出現(xiàn)幅值平均現(xiàn)象,并且無需分路段計算,有效地提高了仿真精度與計算效率。

        原始信號中通常包含較多毛刺或噪聲信號,首先需對其進行去毛刺與濾波處理,去除干擾信號后再計算時域信號的沖擊響應(yīng)譜(SRS)與疲勞損傷譜(FDS),通過公式(3)可基于疲勞損傷譜反推出PSD信號,將PSD信號的極端損傷譜(ERS)與沖擊響應(yīng)譜進行幅值對比后用于確定等效時間,信號處理流程如圖7所示。

        圖7 信號處理流程

        為避免過度加速,需將PSD信號的極端響應(yīng)譜(ERS)與原始信號的沖擊響應(yīng)譜(SRS)進行比較,若極端響應(yīng)譜幅值不超過時域信號沖擊響應(yīng)譜幅值的1.5倍,則該生成信號可用于疲勞計算。

        2.2.2信號處理模塊搭建

        本文基于損傷等效原理對所采集的載荷譜進行信號處理。根據(jù)信號處理流程,分別調(diào)用濾波、去毛刺、PSD合成器等模塊,并設(shè)置相關(guān)參數(shù),根據(jù)處理流程搭建信號處理模塊,如圖8所示。

        圖8 信號處理模塊搭建圖

        首先,采用濾波模塊進行帶寬濾波處理,以過濾原始信號中的噪聲干擾信號;其次調(diào)用重采樣模塊進行重新采樣,采樣算法為線性插值法,采樣頻率設(shè)為512 Hz;最后,調(diào)用PSD生成器模塊生成疲勞計算所需的PSD信號,通過調(diào)節(jié)等效作用時間保證與原始信號的損傷一致性,若極端響應(yīng)譜的幅值過高,需手動減小PSD生成器的信號等效時間以降低幅值。

        將等效時間調(diào)為50 s后,PSD信號的ERS峰值未超過原始時域信號的SRS峰值,說明未出現(xiàn)加速過渡現(xiàn)象,該信號可用于疲勞耐久仿真計算,新生成的PSD信號如圖9所示。

        圖9 加速度PSD信號圖

        3 蓄電池箱振動疲勞耐久分析

        本文采用PSD隨機振動分析方法對蓄電池箱體進行疲勞計算,該方法需將加速度PSD功率譜與有限元頻響應(yīng)力傳遞函數(shù)作為輸入條件,其中PSD功率譜通過基于損傷等效原理的信號處理方法獲得,應(yīng)力傳遞函數(shù)通過有限元仿真計算獲得。

        3.1 有限元模型搭建及計算

        3.1.1有限元模型搭建

        蓄電池箱體主要由外蒙皮、連接支架、底板、電池彈簧壓板、固定螺桿組成,將壓板壓至電池上表面,并通過兩根螺桿將電池組固定于箱體底板,箱體通過固定支架與車架進行連接,本文采用CATIA建立蓄電池及箱體幾何模型,如圖10所示。

        圖10 蓄電池箱體總成幾何模型圖

        在模型前處理過程中,箱體底板、側(cè)板、壓板等沖壓件采用殼單元模擬,根部連接支架為鑄件,采用四面體單元模擬,連接螺栓采用Rbe2+Bar單元進行模擬,中部連接螺桿采用Beam單元模擬,將電池芯體簡化為質(zhì)點并用集中質(zhì)量單元CONM2模擬。蓄電池箱體側(cè)板及加強板材料為510L,底板及電池壓板材料為Q235,連接支架材料為ZL201B,電池總質(zhì)量為113 kg,材料屬性如表1所示。

        表1 蓄電池箱體材料參數(shù)表

        材料名稱楊氏模量/MPa泊松比密度/(t/m-3)屈服強度/MPa抗拉強度/MPa Q235210 0000.37.85×10-9235370 510L210 0000.37.85×10-9355510 ZL201B70 0000.342.7×10-9245295

        為保證蓄電池箱有限元模型與實際裝車狀態(tài)一致,需同時建立氣瓶模型,單個氣瓶重量為4.7 kg,箱體總成有限元模型如圖11所示;該模型共計12.9萬個單元、41.8萬個節(jié)點。經(jīng)查網(wǎng)格質(zhì)量良好,可用于后續(xù)分析計算。

        圖11 蓄電池箱體總成有限元模型圖

        3.1.2頻率響應(yīng)仿真計算

        為得到PSD隨機振動疲勞耐久分析的應(yīng)力傳遞函數(shù),需對該箱體進行頻率響應(yīng)分析。頻率響應(yīng)分析主要用于計算結(jié)構(gòu)在周期載荷作用下的動態(tài)響應(yīng)結(jié)果,是對汽車支架類結(jié)構(gòu)進行強度分析的常規(guī)分析方法[12]。

        因路面激勵頻率較低,頻響分析頻率范圍選取為0~50 Hz,激勵源為箱體根部連接支架螺栓孔處,并對箱體激勵源、、三個方向分別施加1加速度激勵,經(jīng)過頻率響應(yīng)分析計算箱體的響應(yīng)應(yīng)力,其頻率響應(yīng)應(yīng)力曲線如圖12所示。

        圖12 頻率響應(yīng)應(yīng)力曲線圖

        由圖12可得,向最大應(yīng)力為210 MPa,出現(xiàn)在18 Hz時,向最大應(yīng)力為513 MPa,出現(xiàn)在10 Hz時,向最大應(yīng)力為152 MPa,出現(xiàn)在11 Hz時。

        3.2 隨機振動疲勞耐久仿真計算

        在獲得三向PSD激勵載荷與應(yīng)力傳遞函數(shù)后,將其同時導入疲勞軟件進行仿真計算,其中計算模塊可分為兩層:信號處理層與疲勞分析層,原始信號經(jīng)處理后從PSD生成器模塊輸出激勵信號,與頻響應(yīng)力傳遞函數(shù)共同作為輸入信息,導入疲勞分析模塊中進行仿真計算,計算模塊搭建示意圖如圖13所示。

        圖13 疲勞分析模塊搭建示意圖

        本文采用剛性單元模擬焊縫,無法對其進行疲勞壽命計算,文中計算過程主要針對箱體側(cè)板、底板及連接支架等主體部分。在軟件計算中,材料存活率設(shè)置為93%,表2為材料S-N曲線擬合所需參數(shù)。

        表2 S-N曲線材料參數(shù)設(shè)置表

        名稱參數(shù)值 應(yīng)力比?1 存活率系數(shù)93% 缺口敏感系數(shù)1 應(yīng)力集中系數(shù)1.2

        經(jīng)計算,箱體底板損傷分布明顯高于箱體側(cè)板,疲勞周次最小部位位于箱體底板與螺桿連接處,最小疲勞周次為802次,強化路試驗場一圈為6 km,換算實際斷裂里程為4 812 km,蓄電池箱疲勞仿真云圖如圖14所示。

        圖14 疲勞壽命仿真云圖

        在強化路試驗中,經(jīng)試驗人員查驗,箱體斷裂位置為底板螺栓孔處,斷裂里程為6 035 km,未通過強化路試驗?zāi)繕? 000 km,箱體斷裂位置放大圖如圖15所示,斷裂部位疲勞壽命對比結(jié)果如表3所示。

        圖15 箱體底板實際斷裂圖

        表3 仿真與試驗結(jié)果對比表

        名稱斷裂位置斷裂里程/km 疲勞仿真底板中部螺栓孔處4 812 路試試驗底板中部螺栓孔處6 035

        通過對比仿真與試驗結(jié)果可得,該方法基于損傷等效原理成功復現(xiàn)強化路試驗箱體斷裂現(xiàn)象,在信號激勵作用下,疲勞損傷均集中在箱體底板中部螺栓孔附近,與實際斷裂處一致,仿真斷裂里程為4 812 km,與實際斷裂里程較為接近。

        4 蓄電池箱方案優(yōu)化及驗證

        根據(jù)仿真及試驗結(jié)果,初步判斷造成箱體底板斷裂原因為雙螺桿加彈簧壓板結(jié)構(gòu)無法緊固蓄電池,試驗過程中電池出現(xiàn)上下跳動情況,導致箱體底板中部螺栓孔處頻繁出現(xiàn)應(yīng)力集中,從而導致局部斷裂情況,現(xiàn)針對此問題進行方案優(yōu)化。

        根據(jù)分析判斷,更改原方案的箱體底板結(jié)構(gòu)并進行輕量化設(shè)計,材料為Q235,厚度為3 mm;將兩點固定形式更改為六點固定形式;取消電池彈簧壓板,更改為剛性壓板結(jié)構(gòu),其材料為Q345,壓板厚度為2 mm;將蓄電池組更改為縱置形式,通過六根螺桿及壓板將電池組固定于箱體底板上,新方案幾何模型如圖16所示。

        圖16 蓄電池箱新方案幾何模型圖

        為驗證方案整改效果,建立新方案有限元模型,并采用相同載荷譜及仿真方法對新方案進行疲勞仿真分析,箱體損傷結(jié)果如圖17所示。由圖17可得,箱體底板損傷集中區(qū)域依舊為中部螺栓孔處,最小疲勞周次較原方案上升為1 472次,換算為實際斷裂里程為8 832 km。

        圖17 新方案疲勞壽命仿真云圖

        根據(jù)試驗要求,強化路試驗?zāi)繕死锍虨? 000 km,完成新方案樣件試制及裝車后,需對其進行第二輪強化路試驗驗證,新方案裝車狀態(tài)如圖18所示。

        圖18 蓄電池箱新方案裝車實物圖

        經(jīng)試驗,箱體各部件在8 000 km時無損傷斷裂情況,電池壓板無變形,蓄電池組工作正常,試驗過程中無供電異常現(xiàn)象,新方案達成強化路試驗?zāi)繕死锍?,蓄電池箱體方案整改成功。

        5 結(jié)論

        本文基于損傷等效原理將時域信號轉(zhuǎn)換為PSD信號,保證了二者的損傷一致性,對于復雜路面無需分路段計算,保證了損傷精度;并以某牽引車蓄電池箱為研究對象,完成箱體的疲勞仿真與強化路試驗;通過對比可得,該方法較準靜態(tài)法等傳統(tǒng)疲勞分析方法可更好地復現(xiàn)箱體損傷部位,同時破壞里程也較為相近;根據(jù)該仿真結(jié)果對箱體進行方案優(yōu)化,對新方案進行疲勞分析及二次試驗驗證,最終新方案完成輕量化設(shè)計的同時成功通過了強化路試驗?zāi)繕?,有效地減少了項目支出與研發(fā)周期,可將基于損傷等效原理的仿真方法作為此類產(chǎn)品疲勞分析的通用方法,并為其耐久性開發(fā)及驗證提供了理論依據(jù)。

        [1] 黃培鑫,藍鳳崇,陳吉清.隨機振動與沖擊條件下電動車電池包結(jié)構(gòu)響應(yīng)分析[J].汽車工程,2017(9):1087-1099.

        [2] 王文偉,程雨婷,姜衛(wèi)遠,等.電動汽車電池箱結(jié)構(gòu)隨機振動疲勞分析[J].汽車工程學報,2016(6):10-14.

        [3] 彭倩,王繼輝,林有淮,等.基于頻率響應(yīng)分析的車用電池箱結(jié)構(gòu)疲勞壽命預(yù)估[J].機械工程與自動化,2016(5):31-36.

        [4] 冷曉偉,戴作強,鄭莉莉,等.基于nCodeDesignlife的電池箱疲勞壽命研究[J].客車技術(shù)與研究,2018(17):17-19.

        [5] 吳長德,戴江梁,唐煒,等.基于某電動汽車電池箱焊點的疲勞壽命預(yù)測與優(yōu)化[J].機械強度, 2013(3): 663-667.

        [6] Hooper J M, Marco J. Characterising the in-vehicle vibration inputs to the high voltage battery of an electric vehicle[J]. Journal of power sources, 2014, 245:510-519.

        [7] Lang J F, Kjell G. Comparing vibration measurements in an electric vehicle with standard vibration requirements for Li-ion batteries using power spectral density analysis[J].International Journal of Electric and Hybrid Vehicles, 2015,7(3):272-286.

        [8] 唐利濤,楊舟,李剛,等.基于疲勞損傷譜的隨機振動試驗方法在智能電表模擬公路運輸中的研究[J].裝備環(huán)境工程,2019(5):38-42.

        [9] 張勇,董釗志,侯之超,等.基于疲勞損傷譜的動力電池包振動標準分析[J].汽車技術(shù),2020(4):27-32.

        [10] Lalanne C.Mechanical vibration and shock analysis, specification development[M]. John Wiley & Sons,2013.

        [11] 劉金,李偉,劉春立,等.基于載荷譜的主減速器耐久性分析與預(yù)測[J].汽車技術(shù), 2019(10):19-22.

        [12] 吳趙佳,侯永平,張建文.隨機振動條件下的橡膠襯套疲勞壽命預(yù)測[J].汽車技術(shù), 2017(3): 24-28.

        Durability Analysis of Accumulator Case of Commercial Vehicle Based on Principle of Damage Equivalence

        REN Quan, CHEN Wenfei, ZHOU Xiaofei, DU Jiafeng, GUO Ning

        (Shaanxi Automobile Holding Group Co., Ltd., Shaanxi Xi’an 710200 )

        In this research,the load signal of accumulator case was collected by acceleration sensor and the PSD signal consistent with the value of the original signal was processed based on the principle of damage equivalence. The finite element model of accumulator case was established, then the frequency response was analyzed and the transfer function of input excitation and response stress was obtained by frequency response analysis. The damage result of the accumulator case was acquired through the random vibration fatigue analysis, and the feasibility of this method was verified by compared the simulation results with the experimental results. Finally, according to the result of simulation, the structure of the original scheme was optimized, the new scheme was verified by the second round of the fatigue analysis and road test, the target mileage of test was reached successfully, and the new scheme was reformed.

        Accumulator case;Damage equivalence;PSD signal; Fatigue analysis

        U469.72

        A

        1671-7988(2021)23-30-06

        U469.72

        A

        1671-7988(2021)23-30-06

        10.16638/j.cnki.1671-7988.2021.023.009

        任全,碩士研究生,就職于陜西汽車控股集團有限公司,研究方向為疲勞耐久性分析。

        猜你喜歡
        信號處理箱體蓄電池
        《信號處理》征稿簡則
        信號處理(2018年5期)2018-08-20 06:16:02
        《信號處理》第九屆編委會
        信號處理(2018年5期)2018-08-20 06:16:00
        《信號處理》征稿簡則
        信號處理(2018年8期)2018-07-25 12:25:42
        《信號處理》第九屆編委會
        信號處理(2018年8期)2018-07-25 12:24:56
        高牌號灰鐵前端箱體質(zhì)量提升
        聊聊蓄電池的那點事兒(1) 汽車蓄電池的前世
        蓄電池去哪兒了
        超大型冷剪箱體加工難點分析
        基于ANSYS Workbench 的ATB260 減速器箱體模態(tài)分析
        蓄電池去哪兒了?VOL12.雷克薩斯RX450h
        国产一区二区三区免费av| 国产免费专区| 国产一区二区三区免费精品| 国产精品自拍视频在线| 日产精品99久久久久久| 中文无码乱人伦中文视频在线v| 日本免费一区二区三区在线看| 麻豆av在线免费观看精品| 日本道免费一区二区三区日韩精品| 51国产黑色丝袜高跟鞋| 五月激情婷婷丁香| 免费女同毛片在线不卡| 女色av少妇一区二区三区| 天干天干天啪啪夜爽爽av| 国产自产精品露脸刺激91在线| 一区二区日本影院在线观看| 人妻精品视频一区二区三区| 在线综合亚洲欧洲综合网站| 亚洲熟妇网| 国产成人亚洲综合二区| 亚洲另类无码专区首页| 欧美性xxxx狂欢老少配| 四虎国产精品成人影院| 人妻少妇中文字幕,久久精品| 国产一区二区三区乱码| 国内免费AV网站在线观看| 亚洲最黄视频一区二区| 日本道色综合久久影院| 亚洲人成色777777老人头| 高清高速无码一区二区| av在线天堂国产一区| 亚洲裸男gv网站| 久久精品无码一区二区2020| 精品日韩一区二区三区av| 看久久久久久a级毛片| 自拍偷自拍亚洲精品播放| 国产精品一区二区三区蜜臀| 天堂网av一区二区三区在线观看| 久久亚洲精品成人av| 欧美高h视频| 亚洲av高清天堂网站在线观看|