• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    考慮輪胎載荷轉(zhuǎn)移影響的復(fù)合式胎壓監(jiān)測(cè)補(bǔ)償修正方法

    2023-12-17 07:28:30唐先智張曉壯郝少澎王波張宇
    中國(guó)機(jī)械工程 2023年23期

    唐先智 張曉壯 郝少澎 王波 張宇

    摘要:車輛在加速和爬坡時(shí)會(huì)產(chǎn)生載荷轉(zhuǎn)移,導(dǎo)致胎壓監(jiān)測(cè)結(jié)果受到干擾,針對(duì)該問(wèn)題,提出了基于輪胎載荷的胎壓監(jiān)測(cè)結(jié)果補(bǔ)償修正方法。基于雙重?zé)o跡卡爾曼濾波對(duì)車輛質(zhì)量和質(zhì)心位置進(jìn)行了雙重估計(jì),采用遞推最小二乘法(RLS)進(jìn)行了坡度識(shí)別。利用輪胎轉(zhuǎn)轂試驗(yàn)臺(tái)進(jìn)行試驗(yàn),得到載荷、車速、胎壓和輪胎下沉量的數(shù)據(jù),擬合了載荷、車速、胎壓與輪胎下沉量的關(guān)系表達(dá)式。利用坡度識(shí)別,通過(guò)車輛質(zhì)量和質(zhì)心位置的估計(jì),計(jì)算四輪獨(dú)立載荷,修正相應(yīng)的輪速脈沖。試驗(yàn)結(jié)果表明:沒(méi)有載荷修正的胎壓計(jì)算值存在20%的穩(wěn)定誤差和大幅度波動(dòng);加入載荷修正,雖有一定的遲滯性和波動(dòng),但會(huì)很快收斂到真實(shí)值,穩(wěn)定誤差在5%以內(nèi),解決了因載荷變化導(dǎo)致胎壓估算不準(zhǔn)確的問(wèn)題。

    關(guān)鍵詞:輪胎載荷轉(zhuǎn)移;復(fù)合式胎壓監(jiān)測(cè);雙重?zé)o跡卡爾曼濾波;輪胎轉(zhuǎn)轂試驗(yàn)臺(tái)

    中圖分類號(hào):U463.6

    DOI:10.3969/j.issn.1004-132X.2023.23.012

    開放科學(xué)(資源服務(wù))標(biāo)識(shí)碼(OSID):

    Compensation and Correction Method of Compound Tire Pressure Monitoring

    Considering Influences of Tire Load Transfer

    TANG Xianzhi1ZHANG Xiaozhuang1HAO Shaopeng1WANG Bo1ZHANG Yu2

    1.School of Vehicle and Energy,Yanshan University,Qinhuangdao,Hebei,066004

    2.Luan Automobile Technology Co.,Ltd.,Changzhou,Jiangsu,213000

    Abstract:The vehicle would produce load transfer during acceleration and hill climbing,which led to the interference of tire pressure monitoring results. To address this problem,a tire pressure monitoring result compensation correction method was proposed based on tire loads herein. Dual estimation of vehicle mass and center of mass position was carried out based on Dual-UKF,and recursive least square (RLS) was used for slope identification. The data of load,vehicle speed,tire pressure and tire sinkage were collected by using a tire rotation hub test bench. The relational expression of load, speed, tire pressure and tire sinking were obtained by fitting the data. Using slope identification,the four-wheel independent load was calculated by estimating the vehicle mass and center of mass position to correct the corresponding wheel speed pulse. The test results show that the calculated values of tire pressure without load correction have 20% stable error and large fluctuation. with load correction,although there is a certain hysteresis and fluctuation,it will eventually converge to the true value quickly and the stable error is within 5%,which solves the problem of inaccurate tire pressure estimation due to load changes.

    Key words:tire load transfer; composite tire pressure monitoring; dual unscented Kalman filtering(Dual-UKF); tire hub test bench

    收稿日期:2022-12-14

    基金項(xiàng)目:河北省高等學(xué)??茖W(xué)技術(shù)研究重點(diǎn)項(xiàng)目(ZD2019063);國(guó)家自然科學(xué)基金(51505414)

    0 引言

    在汽車行駛安全方面,制動(dòng)防抱死系統(tǒng)(ABS)、胎壓監(jiān)測(cè)系統(tǒng)(TPMS)和安全氣囊被稱為保證汽車安全行駛的三大重要配置。對(duì)胎壓監(jiān)測(cè)系數(shù)的研究主要分為三類:直接式胎壓監(jiān)測(cè)、間接式胎壓監(jiān)測(cè)和復(fù)合式胎壓監(jiān)測(cè)[1。

    直接式胎壓監(jiān)測(cè)通過(guò)安裝在輪胎上的傳感器實(shí)時(shí)更新胎壓胎溫,并且可以檢測(cè)出輪胎是否出現(xiàn)漏氣。張威等[2采用單片機(jī),通過(guò)傳感器采集壓力和溫度數(shù)據(jù),由液晶顯示器顯示胎壓和溫度,并且增加了儲(chǔ)存器,能夠在系統(tǒng)斷電時(shí)保存數(shù)據(jù)。 郭玉言[3基于ARM處理器對(duì)直接式胎壓監(jiān)測(cè)系統(tǒng)進(jìn)行了設(shè)計(jì),利用WIFI模塊與手機(jī)APP進(jìn)行了連接。韓志嶸等[4對(duì)大型汽車直接式胎壓監(jiān)測(cè)進(jìn)行了設(shè)計(jì),不僅可以實(shí)時(shí)監(jiān)測(cè)胎壓和胎溫,還可以自行調(diào)節(jié)胎壓。雖然直接式胎壓監(jiān)測(cè)對(duì)胎壓有著較高的測(cè)量精度,但是成本較高,使用壽命有限,并且容易被外界信號(hào)干擾造成數(shù)據(jù)丟失。

    間接式胎壓監(jiān)測(cè)系統(tǒng)主要是基于ABS的輪速信號(hào),估算方法一般為輪速脈沖法、頻域法、動(dòng)力學(xué)模型法和參數(shù)辨別法。周福強(qiáng)等[5對(duì)實(shí)車采集的數(shù)據(jù)通過(guò)快速傅里葉變換進(jìn)行頻域分析,采用最高峰值法獲得輪胎扭轉(zhuǎn)頻率,從而獲得胎壓。單經(jīng)緯[6首先對(duì)輪速信號(hào)進(jìn)行處理,采用分段線性差值法對(duì)輪速信號(hào)進(jìn)行了重構(gòu),通過(guò)快速傅里葉變換進(jìn)行頻譜分析,提取輪胎共振頻率判斷輪胎是否缺氣。彭加耕等7通過(guò)里程修正輪胎磨損量,基于BP神經(jīng)網(wǎng)絡(luò)對(duì)四輪同時(shí)缺氣進(jìn)行辨識(shí)。王寶琳等[8采用BP神經(jīng)網(wǎng)絡(luò)修正車輛轉(zhuǎn)彎所造成的內(nèi)外側(cè)輪胎脈沖存在差異的問(wèn)題,利用三均值比較法對(duì)右后輪胎壓進(jìn)行判斷。羅新帥[9、董倩倩10基于輪速脈沖法和頻域法設(shè)計(jì)了間接式胎壓監(jiān)測(cè)算法,對(duì)兩種方法進(jìn)行融合實(shí)現(xiàn)了四輪同程度欠壓識(shí)別。尹龍川11通過(guò)輪速分析,融合頻率法和半徑法對(duì)胎壓進(jìn)行監(jiān)測(cè)。間接式胎壓監(jiān)測(cè)成本低,但是無(wú)法顯示胎壓值,只能在胎壓達(dá)到報(bào)警線時(shí)提示駕乘人員。

    復(fù)合式胎壓監(jiān)測(cè)安裝一個(gè)或兩個(gè)直接式胎壓壓力傳感器,使用間接式的方法實(shí)現(xiàn)對(duì)其他輪胎胎壓的估算。劉澍[12嵌入了基于同側(cè)輪輪速脈沖比較法的間接式胎壓監(jiān)測(cè)單輪算法,并在直線工況下實(shí)現(xiàn)了單輪缺氣定位報(bào)警。吳飛龍等13設(shè)計(jì)的復(fù)合式胎壓監(jiān)測(cè)方法在輪胎只有一個(gè)輪胎單獨(dú)欠壓或是兩個(gè)輪胎欠壓情況下使用半徑法,其余情況使用頻率法。復(fù)合式胎壓監(jiān)測(cè)系統(tǒng)成本比直接式胎壓監(jiān)測(cè)系統(tǒng)低,并且還可以實(shí)時(shí)監(jiān)測(cè)輪胎胎壓。然而,車輛在實(shí)際行駛過(guò)程中,由于加速度、道路坡度的存在,車輛會(huì)產(chǎn)生載荷轉(zhuǎn)移,導(dǎo)致前軸載荷與后軸載荷分配情況發(fā)生變化;并且油箱內(nèi)現(xiàn)有的油量、載客量、載貨量等變化,也會(huì)導(dǎo)致軸荷分配發(fā)生變化;車輛人員的數(shù)量和分布情況不同,也會(huì)在一定程度上影響軸荷分配。載荷的變化會(huì)影響胎壓監(jiān)測(cè)的準(zhǔn)確性。

    基于以上問(wèn)題,本文采用輪速比較法設(shè)計(jì)復(fù)合式胎壓監(jiān)測(cè)系統(tǒng),并且對(duì)車輛質(zhì)量、質(zhì)心位置、道路坡度進(jìn)行識(shí)別,對(duì)不同載荷下胎壓與輪胎下沉量進(jìn)行分析,驗(yàn)證載荷和胎壓之間存在一定的關(guān)系式,提出了載荷補(bǔ)償修正胎壓的方法,并且通過(guò)試驗(yàn)驗(yàn)證了方法的正確性。

    1 車輛動(dòng)力學(xué)模型搭建

    本文以四輪獨(dú)立懸掛、前置前驅(qū)、前輪轉(zhuǎn)向車輛為研究對(duì)象,測(cè)量值信號(hào)來(lái)自于車載傳感器和CAN總線。

    根據(jù)牛頓第二定律建立質(zhì)量與坡度估計(jì)所需的車輛縱向動(dòng)力學(xué)模型,如圖1所示。

    由牛頓第二定律進(jìn)行受力分析可得到車輛在坡道上的縱向動(dòng)力學(xué)表達(dá)式:

    式中,δm為車輛換算質(zhì)量;v縱向速度;Ft為車輛驅(qū)動(dòng)力;Fi為坡道阻力;Fw為空氣阻力;Ff為滾動(dòng)阻力;Tm為滾動(dòng)阻力矩;ig為變速器傳動(dòng)比;i0為主減速器傳動(dòng)比;ηt為傳動(dòng)效率;r0為車輪半徑;Cd為空氣阻力系數(shù);A為車輛迎風(fēng)面積;ρ為空氣密度;f為滾動(dòng)阻力系數(shù);m為車輛質(zhì)量;g為重力加速度,取9.8 m/s2;α為道路坡度。

    2 基于Dual-UKF的車輛參數(shù)估計(jì)

    選擇合適的濾波算法,利用汽車縱向動(dòng)力學(xué)對(duì)車輛參數(shù)進(jìn)行估計(jì)是核心環(huán)節(jié)。在車輛使用過(guò)程中,車輛上發(fā)生的加載或卸載會(huì)導(dǎo)致整車質(zhì)量和質(zhì)心位置變化,這些變化會(huì)導(dǎo)致前后軸載荷不同,從而使后輪車輪輪速脈沖發(fā)生變化,導(dǎo)致胎壓的計(jì)算估計(jì)值增大或減小,所以,為了減小這些影響,在傳統(tǒng)估計(jì)方法中,會(huì)將質(zhì)心位置和整車質(zhì)量等參數(shù)直接作為常數(shù)來(lái)使用,沒(méi)有考慮車輛參數(shù)的改變對(duì)結(jié)果造成的影響。雙重?zé)o跡卡爾曼濾波算法(Dual-UKF)的優(yōu)點(diǎn)在于可將可能變化的車輛參數(shù)量考慮進(jìn)去,并且考慮了車輛參數(shù)量與狀態(tài)量之間的互相影響關(guān)系(互為輸入輸出),提高估計(jì)準(zhǔn)確率,因此本文基于雙重?zé)o跡卡爾曼濾波算法實(shí)現(xiàn)對(duì)車輛質(zhì)量與車輛質(zhì)心位置的高效追蹤估計(jì),從而對(duì)胎壓估計(jì)值進(jìn)行修正。

    經(jīng)典卡爾曼濾波的系統(tǒng)狀態(tài)方程是線性的。對(duì)于線性系統(tǒng),其狀態(tài)空間表達(dá)式如下:

    式中,N為狀態(tài)轉(zhuǎn)移矩陣;B、H為轉(zhuǎn)換矩陣;xk為狀態(tài)量;uk-1為控制量;wk-1為過(guò)程噪聲;vk為測(cè)量噪聲;k為離散時(shí)間。

    先驗(yàn)估計(jì)誤差、協(xié)方差計(jì)算為

    式中,X^-k為先驗(yàn)狀態(tài)估計(jì)值;P-k為先驗(yàn)狀態(tài)協(xié)方差;Q為過(guò)程噪聲的協(xié)方差;Pk-1為后驗(yàn)狀態(tài)協(xié)方差。

    后驗(yàn)估計(jì)更新計(jì)算誤差協(xié)方差為

    式中,X^k為當(dāng)前狀態(tài)的估計(jì)值;Zk為觀測(cè)量;Kk為卡爾曼增益。

    建立考慮參數(shù)變化的車輛狀態(tài)空間方程與觀測(cè)方程如下:

    式中,θ~k為狀態(tài)量;f(xk,θ~k,uk)、h(xk,uk)為多元函數(shù)式;ws,k為車輛狀態(tài)估計(jì)時(shí)的系統(tǒng)過(guò)程噪聲,其協(xié)方差為Qs;vs,k為車輛狀態(tài)估計(jì)時(shí)的測(cè)量噪聲,其協(xié)方差為Rs;qk是被估計(jì)的車輛參數(shù)量與狀態(tài)量。

    觀測(cè)量Zk表示為

    式中,ax為縱向加速度;ay為側(cè)向加速度;ωr為橫擺角速度。

    xk為車輛的參數(shù)量,x~k表示在狀態(tài)估計(jì)時(shí)將參數(shù)視作定值,要估計(jì)的車輛參數(shù)量與狀態(tài)量分別如下:

    式中,m為整車質(zhì)量;a0為質(zhì)心至前軸距離;vx為縱向車速;vy為側(cè)向車速;β為轉(zhuǎn)向盤轉(zhuǎn)角。

    輸入量uk

    式中,δ為前輪轉(zhuǎn)角;ωij為四輪輪速信號(hào);Tij為四輪扭矩。

    輸入量離散方程表示如下:

    其中,Ts為最小仿真步長(zhǎng),輸入量基于前項(xiàng)歐拉法進(jìn)行離散化處理,使其可以正常輸出。

    由于車輛參數(shù)量xk的變化速率緩慢,因此在第一層估計(jì)結(jié)果出來(lái)后,將車輛參數(shù)量認(rèn)作常數(shù),得到參數(shù)方程與觀測(cè)方程如下:

    式中,wvar,k為車輛狀態(tài)估計(jì)時(shí)的系統(tǒng)過(guò)程噪聲,其協(xié)方差為Qvar;vvar,k為車輛狀態(tài)估計(jì)時(shí)的測(cè)量噪聲,其協(xié)方差為Rvar。

    利用CarSim創(chuàng)建測(cè)試工況,輸出所需的控制和測(cè)量參數(shù)量,利用MATLAB/Simulink搭建Dual-UKF濾波器,對(duì)比CarSim輸出的車輛行駛狀態(tài)以驗(yàn)證估計(jì)結(jié)果的正確性。

    車型選擇CarSim中的C級(jí)車,設(shè)置駕駛員模型中的控制量。仿真路面設(shè)置為附著系數(shù)μ=0.9的平坦路面,以60 km/h的初速度進(jìn)行仿真,輪缸壓力設(shè)置為pb=0.5 MPa,以0.5 Hz頻率的正弦波信號(hào)作為轉(zhuǎn)向盤輸入,角度范圍為[-90°,90°]。

    通過(guò)參數(shù)調(diào)試得到雙重?zé)o跡卡爾曼濾波參數(shù)估計(jì)器中的過(guò)程噪聲Qvarc=diag(1,0.1,1),測(cè)量噪聲Rvarc=0.01I3×3;狀態(tài)估計(jì)器的過(guò)程噪聲Qvars=10-8I3×3,測(cè)量噪聲Rvars=I3×3。

    CarSim的C級(jí)車的整車質(zhì)量m=1655 kg,質(zhì)心至前軸的距離a0=1.30 m。整車質(zhì)量估計(jì)如圖2所示,質(zhì)心位置至前軸的距離如圖3所示。

    由圖2與圖3可知,隨著仿真的進(jìn)行,整車質(zhì)量m、質(zhì)心至前軸的距離a0估計(jì)量有所波動(dòng),但最終都收斂到對(duì)應(yīng)的真實(shí)值附近,因此該觀測(cè)器可以實(shí)現(xiàn)對(duì)車輛質(zhì)心及質(zhì)量預(yù)測(cè)。

    3 基于RLS算法的道路坡度識(shí)別

    由于遞推最小二乘(RLS)估計(jì)算法相比其他方法需求的參量更少、實(shí)時(shí)性更好,且對(duì)單輸入單輸出系統(tǒng)而言,無(wú)需額外先驗(yàn)條件,因此本文選用RLS估計(jì)算法實(shí)現(xiàn)對(duì)路面坡度的估計(jì)。

    聯(lián)立式(1)~式(5),得到車輛縱向動(dòng)力學(xué)公式表達(dá):

    建立n階自回歸(AR)模型:

    式中,y(t)為系統(tǒng)輸出;a1,a2,…,an為待估計(jì)參數(shù);y(t-1),…,y(t-n)為測(cè)量數(shù)據(jù);ε(t)為系統(tǒng)白噪聲,均值為零且方差大于零。

    遞推最小二乘法參數(shù)估計(jì)表達(dá)式如下:

    式中,θ^(t+1)為待估計(jì)參數(shù);K(t)為增益矩陣;P(t)為協(xié)方差矩陣;φ(t)為輸入向量;λ為遺忘因子。

    由于新舊觀測(cè)值相互影響較小,故令遺忘因子λ=0.95。根據(jù)RLS算法,由式(24)~式(26)可得路面坡度估計(jì)值α為

    為驗(yàn)證坡度估計(jì)算法模塊有效性,同樣利用MATLAB/Simulink搭建狀態(tài)估計(jì)觀測(cè)器模型,并與CarSim中設(shè)定好的車輛模型進(jìn)行聯(lián)合仿真測(cè)試,仿真流程如圖4所示。圖4中,T_總表示發(fā)動(dòng)機(jī)輸出扭矩,a_x表示縱向加速度,V_x表示縱向速度,m_0表示汽車質(zhì)量初始值。CarSim輸出信號(hào)至MATLAB/Simulink所搭建好的模型中,輸出估計(jì)的坡度值。

    中道路條件設(shè)置為30°的坡度,仿真車輛速度為40 km/h勻速,仿真結(jié)果如圖5所示。

    由圖5可知,真實(shí)路面坡度值在26%~32%之間,觀測(cè)結(jié)果在5s內(nèi)快速收斂到真實(shí)坡度值附近,估計(jì)誤差在5%之內(nèi),該觀測(cè)器具有較高的估計(jì)精度且穩(wěn)定性較好,可以實(shí)現(xiàn)對(duì)道路坡度的預(yù)測(cè)。

    4 載荷對(duì)輪速的影響分析

    為進(jìn)一步探究載荷對(duì)輪胎滾動(dòng)的影響機(jī)制,本文利用輪胎轉(zhuǎn)轂試驗(yàn)臺(tái)采集輪胎在不同垂向載荷、胎壓狀態(tài)下的輪胎變形量信息,根據(jù)試驗(yàn)結(jié)果建立各參量之間的關(guān)系。

    搭建的單輪轉(zhuǎn)轂試驗(yàn)臺(tái)如圖6所示,該試驗(yàn)臺(tái)由驅(qū)動(dòng)電機(jī)、轉(zhuǎn)轂、氣泵、壓力缸、單邊懸架、輪胎以及相應(yīng)的傳感器組成,轉(zhuǎn)轂試驗(yàn)臺(tái)采用SC160復(fù)動(dòng)氣缸進(jìn)行加載,通過(guò)打氣泵使氣體經(jīng)過(guò)油水分離器的過(guò)濾作用于加載氣缸上,控制復(fù)動(dòng)氣缸輸出壓力模擬載荷,并將載荷施加在麥弗遜懸架的減振彈簧及擺臂上,最終實(shí)現(xiàn)控制輪胎所受到的垂直載荷,調(diào)節(jié)范圍為0~10 kN。同時(shí)利用輪胎壓力傳感器、激光位移傳感器、氣缸壓力傳感器來(lái)獲取對(duì)應(yīng)的載荷、胎壓、輪胎變形量等信息。

    當(dāng)輪胎壓力固定時(shí),隨著輪胎加載力的增大,輪胎的變形量逐漸增大。根據(jù)日本學(xué)者FUJIMOTO經(jīng)驗(yàn)結(jié)果[14,輪胎下沉量經(jīng)驗(yàn)公式可表示如下:

    式中,F(xiàn)z為載荷,kN;p為胎壓,kPa;σ為下沉量,mm;a1、b1、c1、d1、e1為待標(biāo)定系數(shù)。

    本文選擇的輪胎型號(hào)為Cross Wind HP品牌的205 60 R16 92H;胎壓的選擇為標(biāo)準(zhǔn)胎壓、85%標(biāo)準(zhǔn)胎壓、70%標(biāo)準(zhǔn)胎壓;載荷從2.0 kN開始,以0.5 kN為間隔,逐漸增大載荷至6.0 kN;載荷5 kN,胎壓分別為150 kPa(1.5 bar)、240 kPa(2.4 bar)時(shí),車速以10 km/h為間隔,從40 km/h依次增大至120 km/h。建立輪胎承受載荷、輪胎氣壓和車速與輪胎變形量對(duì)應(yīng)關(guān)系,得到胎壓、載荷與輪胎下沉量的變化關(guān)系,如圖7所示,車速與輪胎下沉量的關(guān)系如圖8所示。

    由圖7可以看出,隨著對(duì)輪胎施加的載荷增大,輪胎的下沉量也增大,并且在同一載荷的情況下,胎壓越高,輪胎下沉量越小。載荷2~6 kN處于線性區(qū)域,并且在胎壓較低的情況下,變形量變化范圍也會(huì)增加。

    由圖8可以看出,車速對(duì)輪胎變形量影響較小,當(dāng)車速?gòu)?0 km/h增加至120 km/h時(shí),在輪胎低壓的情況下,輪胎下沉量變化為2.4 mm,相對(duì)于輪胎低壓狀態(tài),標(biāo)準(zhǔn)胎壓時(shí),從低速至高速,輪胎下沉量不到2 mm。

    利用式(28)對(duì)圖7和圖8中的試驗(yàn)數(shù)據(jù)進(jìn)行擬合,得到未知系數(shù)a1、b1、c1、d1、e1,擬合后的下沉量經(jīng)驗(yàn)公式為

    式中,E1、F1為整理后的擬合系數(shù),E1為1.771、1.747、1.476,F(xiàn)1為10.43、8.594、7.524。

    為驗(yàn)證載荷是否對(duì)輪速脈沖信號(hào)有影響,將測(cè)試輪胎安裝在試驗(yàn)車上,通過(guò)在車內(nèi)加載或卸載,模擬后軸車輪分別在2.0 kN、4.0 kN、6.0 kN的載荷工況,通過(guò)使用10 kg的沙袋改變后軸輪胎垂直載荷,車輛以60 km/h速度做勻速直線行駛運(yùn)動(dòng);通過(guò)試驗(yàn)車的車載CAN通信網(wǎng)絡(luò),利用Kvaser leaf v2 CAN總線分析儀,從總線上采集所需輪速脈沖等信號(hào),并將所采集到的信號(hào)以.log文件形式儲(chǔ)存至指定文件夾12,之后利用MATLAB建立函數(shù)腳本,按照試驗(yàn)車型的CAN網(wǎng)絡(luò)通信協(xié)議,針對(duì)所需信號(hào)進(jìn)行解析、換算、濾波、修正[9,得到各組試驗(yàn)下的輪速信號(hào),結(jié)果如圖9a、圖9b、圖9c所示,各組載荷速度差值對(duì)比如圖9d所示。

    將三種不同載荷下的前后輪速濾波處理后做均差,結(jié)果如圖9d所示,結(jié)合圖7可以看出,隨著載荷的增大,后輪輪胎下沉量增大,輪胎滾動(dòng)半徑減小,使車輛前后軸之間的輪速脈沖差值增大,后續(xù)會(huì)影響輪速法計(jì)算得到的實(shí)際胎壓值。

    5 基于載荷修正的復(fù)合式胎壓監(jiān)測(cè)算法

    復(fù)合式胎壓監(jiān)測(cè)是在間接式胎壓監(jiān)測(cè)算法基礎(chǔ)上,引入一個(gè)或兩個(gè)胎壓壓力傳感器的信息用以估算未安裝輪胎壓力傳感器車輪的胎壓。現(xiàn)有的復(fù)合式胎壓監(jiān)測(cè)均以輪速法為基本方法。

    輪速脈沖比較法的原理是利用ABS系統(tǒng)中輪速傳感器的脈沖信號(hào)計(jì)算各個(gè)車輪轉(zhuǎn)速,由于輪胎缺氣時(shí)其滾動(dòng)半徑會(huì)減小轉(zhuǎn)速會(huì)增加,故通過(guò)對(duì)比四輪輪速脈沖信號(hào)之間的差異就可以判斷出輪胎是否出現(xiàn)缺氣狀況,復(fù)合式胎壓監(jiān)測(cè)系統(tǒng)主要流程如圖10所示。

    脈沖數(shù)組的累計(jì)是通過(guò)累計(jì)輪速脈沖數(shù)據(jù)的方法放大各輪速信號(hào)Ni之間的差異。為方便數(shù)據(jù)后期處理,更新之前首先對(duì)輪速信號(hào)進(jìn)行歸一化處理。同時(shí)也利用輪速信號(hào)累計(jì)的過(guò)程計(jì)算車輛速度、加速度、方向轉(zhuǎn)角等狀態(tài)值在累計(jì)期間的均值,用于判斷樣本可靠性,決定本組樣本的剔除或修正。本文選取固定長(zhǎng)度Nstd=5000,即每累計(jì)5000個(gè)輪速脈沖樣本進(jìn)行一次更新。其中脈沖數(shù)組的累計(jì)具體步驟與部分表達(dá)式如圖11所示。其中,Ni為輪速信號(hào)(i=0,1,2,3),Acc±為采樣周期內(nèi)的加速度、減速度均值,v為車速信號(hào)值,str為轉(zhuǎn)向盤轉(zhuǎn)角信號(hào)值,Degr,l為采樣周期內(nèi)右轉(zhuǎn)角、左轉(zhuǎn)角均值,Ni為輪速脈沖樣本數(shù)據(jù),v為采樣周期內(nèi)的車速均值。

    為提高更新速率,本文采用滾動(dòng)時(shí)域算法數(shù)據(jù)結(jié)構(gòu),數(shù)據(jù)結(jié)構(gòu)如圖12所示,以Nstd為固定窗口長(zhǎng)度,按時(shí)域順序依次存放樣本數(shù)據(jù),數(shù)據(jù)隊(duì)列的方式保證了隊(duì)列中的樣本數(shù)據(jù)維度不變。

    利用自適應(yīng)雙重疊加濾波器模塊對(duì)輪速信號(hào)進(jìn)行預(yù)處理。之后使用Busmaster中的Replay功能,對(duì)采集到的不同階梯載荷下的勻速直線行駛工況樣本進(jìn)行回放仿真,利用標(biāo)準(zhǔn)輪速脈沖值與胎壓計(jì)算公式進(jìn)行擬合標(biāo)定。胎壓計(jì)算擬合公式如下:

    式中,pi為各非基準(zhǔn)輪胎壓,kPa;p0為基準(zhǔn)輪胎壓,kPa;Ai、Bi為胎壓計(jì)算公式標(biāo)定系數(shù);N^i為各非基準(zhǔn)輪修正并濾波后輪速脈沖樣本值,i=1,2,3;N^0為基準(zhǔn)輪濾波后輪速脈沖樣本值。

    Ai、Bi標(biāo)定后的胎壓擬合系數(shù)數(shù)值為

    對(duì)標(biāo)定后的胎壓擬合公式設(shè)計(jì)測(cè)試工況,測(cè)試方法按照GB/T 26149《乘用車輪胎氣壓監(jiān)測(cè)系統(tǒng)的性能要求和試驗(yàn)方法》進(jìn)行,測(cè)試內(nèi)容為直線行駛工況胎壓計(jì)算顯示精度測(cè)試。設(shè)計(jì)階梯載荷試驗(yàn)工況,以150 kg載荷工況為基準(zhǔn)進(jìn)行胎壓計(jì)算公式參數(shù)擬合,待參數(shù)擬合完成后依次進(jìn)行一次卸載與兩次階梯加載過(guò)程。為保證試驗(yàn)結(jié)果的準(zhǔn)確性,測(cè)試以60 km/h速度進(jìn)行勻速直線行駛,并保證每個(gè)載荷區(qū)間段均有足夠的有效樣本數(shù)量進(jìn)入。

    試驗(yàn)車輛為前置前驅(qū)車型,因?yàn)榘l(fā)動(dòng)機(jī)艙、變速箱等因素,導(dǎo)致車輛的簧上載荷大部分都分布于前軸,在車內(nèi)載荷發(fā)生變化時(shí),前軸載荷發(fā)生的變化相對(duì)后軸載荷變化較小,本文所設(shè)計(jì)的復(fù)合式胎壓監(jiān)測(cè),選擇的基準(zhǔn)輪為左前輪和右前輪。圖13所示是車輛載荷發(fā)生變化,但沒(méi)有對(duì)載荷進(jìn)行修正所計(jì)算得到的兩后輪胎壓結(jié)果。

    由圖13可以看出,由于標(biāo)定過(guò)程中未考慮載荷變化帶來(lái)的影響,所以當(dāng)載荷差值發(fā)生較大變化時(shí),計(jì)算結(jié)果出現(xiàn)較大的穩(wěn)態(tài)誤差。以上述工況為例,標(biāo)準(zhǔn)胎壓為230 kPa時(shí),最大計(jì)算值為246.3 kPa,最小計(jì)算值為198.6 kPa,計(jì)算誤差波動(dòng)范圍超過(guò)20%。

    車輛在行駛過(guò)程中,輪胎的垂向載荷會(huì)根據(jù)不同的車輛工況和路況的改變而發(fā)生改變。軸荷會(huì)根據(jù)加速度的變化而發(fā)生載荷轉(zhuǎn)移,利用車輛動(dòng)力學(xué)對(duì)四輪獨(dú)立載荷進(jìn)行計(jì)算,實(shí)現(xiàn)實(shí)時(shí)更新,具體計(jì)算公式為

    式中,a為質(zhì)心至前軸的距離;b為質(zhì)心至后軸的距離;lf為前軸輪距;lr為后軸輪距;h為車輛質(zhì)心高度。

    結(jié)合第二節(jié)對(duì)m0、a0的估計(jì)結(jié)果可以得到整車質(zhì)量分布在車輛后軸的載荷得到修正量,具體表達(dá)式為

    結(jié)合式(30),可以得到修正后的胎壓計(jì)算公式:

    式中,Δ為載荷影響修正量;τ為擬合換算系數(shù);a1為車輛軸距。

    第三節(jié)對(duì)坡度識(shí)別的結(jié)果顯示,坡度識(shí)別具有一定的遲滯性,并且僅通過(guò)車輛動(dòng)力學(xué)模型實(shí)現(xiàn)坡度高精度識(shí)別或?qū)崿F(xiàn)變坡度高精度識(shí)別較困難,所以本文利用坡度識(shí)別對(duì)路面進(jìn)行篩選,若路面坡度過(guò)大,四輪獨(dú)立載荷變化較大,則此段的數(shù)據(jù)會(huì)被過(guò)濾,以保證數(shù)據(jù)的準(zhǔn)確性,實(shí)現(xiàn)對(duì)胎壓的準(zhǔn)確識(shí)別。圖14所示為車輛載荷發(fā)生變化且對(duì)載荷進(jìn)行修正后所計(jì)算得到的兩后輪胎壓結(jié)果。

    分析上述試驗(yàn)結(jié)果,由于濾波算法估計(jì)參數(shù)的滯后性,當(dāng)載荷差值發(fā)生較大變化時(shí),計(jì)算依然出現(xiàn)小波動(dòng)的偏移誤差,但隨著估計(jì)值的逐步趨穩(wěn),胎壓計(jì)算值收斂至胎壓真實(shí)值附近。如圖14所示,標(biāo)準(zhǔn)胎壓為230 kPa時(shí),最大計(jì)算值為238.1 kPa,最小計(jì)算值為225.6 kPa,計(jì)算誤差波動(dòng)范圍被控制在了5%以內(nèi),穩(wěn)態(tài)誤差得到了很好的控制。

    6 結(jié)論

    (1)基于輪速比較法采用車輛載荷修正復(fù)合式胎壓監(jiān)測(cè)系統(tǒng),使用車輛縱向動(dòng)力學(xué)模型基于雙重?zé)o跡卡爾曼濾波(Dual-UKF)對(duì)車輛質(zhì)量與質(zhì)心位置進(jìn)行了雙重估計(jì),從CarSim仿真試驗(yàn)中可以看出,該方法有著較高的精度與收斂速度,實(shí)現(xiàn)了對(duì)目標(biāo)的跟蹤估計(jì)。利用遞推最小二乘法(RLS)對(duì)道路坡度進(jìn)行了識(shí)別,并在CarSim軟件中進(jìn)行了仿真試驗(yàn),誤差在5%左右,驗(yàn)證了算法的可行性。

    (2)搭建了輪胎轉(zhuǎn)轂試驗(yàn)臺(tái),研究了載荷、車速、胎壓與輪胎變形量之間的關(guān)系,擬合出輪胎下沉量經(jīng)驗(yàn)公式。將測(cè)試輪胎裝入實(shí)車,進(jìn)行了在后軸車輪施加2 kN、4 kN和6 kN不同載荷的試驗(yàn),將前后輪速濾波處理后取平均作差,發(fā)現(xiàn)載荷會(huì)導(dǎo)致輪胎輪速發(fā)生變化,后續(xù)會(huì)影響基于輪速法的復(fù)合式胎壓監(jiān)測(cè)的胎壓計(jì)算。

    (3)通過(guò)坡度識(shí)別過(guò)濾車輛在坡道上行駛的數(shù)據(jù),使用Busmaster軟件對(duì)采集到的不同階梯載荷下的勻速直線行駛工況樣本進(jìn)行回放仿真,擬合標(biāo)準(zhǔn)輪速脈沖值對(duì)胎壓的計(jì)算。仿真試驗(yàn)表明,沒(méi)有加入載荷修正時(shí),存在20%的穩(wěn)定誤差和大幅度波動(dòng);加入載荷修正,雖有一定的遲滯性和波動(dòng),但是最終會(huì)很快收斂到真實(shí)值,穩(wěn)定誤差在5%以內(nèi),解決了因坡度和載荷導(dǎo)致胎壓值計(jì)算不準(zhǔn)確的問(wèn)題。

    參考文獻(xiàn):

    [1]王勇. 淺談胎壓監(jiān)測(cè)系統(tǒng)TPMS的類型和標(biāo)準(zhǔn)實(shí)施現(xiàn)狀[J]. 科技資訊,2017,15(9):26-27.

    WANG Yong. Discuss the Types and Implementation Status of TPMS for Tire Pressure Monitoring Systems[J]. Science & Technology Information,2017,15(9):26-27.

    [2]張威,王寧波,陳美玲. 基于單片機(jī)的胎壓監(jiān)測(cè)裝置設(shè)計(jì)[J]. 現(xiàn)代信息科技,2022,6(2):70-74.

    ZHANG Wei,WANG Ningbo,CHEN Meiling. Design on Tire Pressure Monitoring Device Based on Single Chip Microcomputer[J]. Modern Information Technology,2022,6(2):70-74.

    [3]郭玉言. 基于ARM的汽車胎壓監(jiān)測(cè)系統(tǒng)[J]. 自動(dòng)化技術(shù)與應(yīng)用,2020,39(12):173-178.

    GUO Yuyan. ARM-based Automotive Tire Pressure Monitoring System[J]. Techniques of Automation and Applications,2020,39(12):173-178.

    [4]韓志嶸,林蕓,劉潔. 基于ARM的汽車胎壓監(jiān)測(cè)系統(tǒng)[J]. 汽車實(shí)用技術(shù),2019(1):114-115.

    HAN Zhirong,LIN Yun,LIU Jie. ARM Based Automotive Tire Pressure Monitoring System[J]. Automobile Applied Technology,2019(1):114-115.

    [5]周福強(qiáng),王少紅,蘭寧,等. 基于輪胎扭轉(zhuǎn)高階頻點(diǎn)的汽車胎壓監(jiān)測(cè)信號(hào)分析[J]. 汽車科技,2021(6):7-11.

    ZHOU Fuqiang,WANG Shaohong,LAN Ning,et al. Analysis of Automobile Tire Pressure Monitoring Signal Based on High Order Frequency Point of Tire Torsion[J]. Auto Sci-Tech. ,2021(6):7-11.

    [6]單經(jīng)緯. 基于頻域特征分析的間接式胎壓監(jiān)測(cè)算法研究[D]. 長(zhǎng)春:吉林大學(xué),2016.

    SHAN Jingwei. The Research of Indirect Tire Pressure Monitoring Algorithm Based on Frequency Domain Feature Analysis[D]. Changchun:Jilin University,2016.

    [7]彭加耕,王宇岱,韓宗奇,等. 高速工況下間接式胎壓監(jiān)測(cè)的影響因素研究與優(yōu)化[J]. 汽車技術(shù),2020(9):45-51.

    PENG Jiageng,WANG Yudai,HAN Zongqi,et al. Research and Optimization on Influencing Factors on Indirect Tire Pressure Monitoring in High Speed Driving Conditions[J]. Automobile Technology,2020(9):45-51.

    [8]王寶琳,夏懷成,董倩倩. 車輛轉(zhuǎn)彎工況下間接式胎壓監(jiān)測(cè)系統(tǒng)脈沖數(shù)的修正[J]. 汽車安全與節(jié)能學(xué)報(bào),2021,12(4):507-515.

    WANG Baolin,XIA Huaicheng,DONG Qianqian. Correction for Pulse Number of Indirect Tire Pressure Monitoring System under Vehicle Turning Condition[J]. Journal of Automotive Safety and Energy,2021,12(4):507-515.

    [9]羅新帥. 基于輪胎滾動(dòng)半徑和頻域特征的間接式胎壓監(jiān)測(cè)算法研究[D]. 長(zhǎng)春:吉林大學(xué),2021.

    LUO Xinshuai. Research on Indirect Tire Pressure Monitoring Algorithm Based on Tire Rolling Radius and Frequency Domain Feature[D]. Changchun:Jilin University,2021.

    [10]董倩倩. 基于輪速脈沖數(shù)和頻率法融合的胎壓監(jiān)測(cè)算法研究[D]. 秦皇島:燕山大學(xué),2021.

    DONG Qianqian. Research on Tire Pressure Monitoring Algorithm Based on the Fusion of Wheel Speed Pulse Number and Frequency Method[D]. Qinhuangdao:Yanshan University,2021.

    [11]尹龍川. 基于輪速分析的間接式胎壓監(jiān)測(cè)算法研究[D]. 武漢:湖北工業(yè)大學(xué),2021.

    YIN Longchuan. Research on Indirect Tire Pressure Monitoring Algorithm Based on Wheel Speed Analysis[D]. Wuhan:Hubei University of Technology,2021.

    [12]劉澍. 直線行駛工況下的復(fù)合式胎壓監(jiān)測(cè)算法研究[D]. 秦皇島:燕山大學(xué),2021.

    LIU Shu. Research on Compound Tire Pressure Monitoring Algorithm under Straight Line Driving Condition[D]. Qinhuangdao:Yanshan University,2021.

    [13]吳飛龍,郭世永. 復(fù)合間接式胎壓監(jiān)測(cè)機(jī)制研究[J]. 青島理工大學(xué)學(xué)報(bào),2022,43(1):106-114.

    WU Feilong,GUO Shiyong. Research on Monitoring Mechanism of Compound Indirect Tire Pressure[J]. Journal of Qingdao University of Technology,2022,43(1):106-114.

    [14]莊繼德. 汽車輪胎學(xué)[M]. 北京:化學(xué)工業(yè)出版社,1988.

    ZHUANG Jide. Automotive Tyres[M]. Beijing:Chemical Industry Press,1988.

    (編輯 袁興玲)

    作者簡(jiǎn)介:唐先智,男,1982年生,副教授。研究方向?yàn)橹悄茌喬ァ⑵囅到y(tǒng)動(dòng)力學(xué)及控制、無(wú)人駕駛汽車駕駛意圖識(shí)別及駕乘感受識(shí)別、混合動(dòng)力/電動(dòng)汽車能量管理策略。E-mail:tangxz@ysu.edu.cn。王 波(通信作者),女,1983年生,博士,講師。研究方向?yàn)楣?jié)能與新能源汽車和液壓混合動(dòng)力控制策略。E-mail:wangbo2015@ysu.edu.cn。

    国产97色在线日韩免费| 国产私拍福利视频在线观看| 欧美日韩黄片免| 亚洲人成网站在线播放欧美日韩| 在线观看一区二区三区| 欧美性猛交╳xxx乱大交人| 高清毛片免费观看视频网站| 正在播放国产对白刺激| 国产精品98久久久久久宅男小说| 国产不卡一卡二| 舔av片在线| 成年人黄色毛片网站| www.精华液| 黑人操中国人逼视频| 麻豆久久精品国产亚洲av| 免费av毛片视频| 丝袜美腿诱惑在线| 又爽又黄无遮挡网站| 一区二区三区国产精品乱码| 97碰自拍视频| 波多野结衣高清作品| 这个男人来自地球电影免费观看| 亚洲精品国产一区二区精华液| 国产伦在线观看视频一区| ponron亚洲| 亚洲,欧美精品.| 日本a在线网址| 国内精品久久久久久久电影| 韩国av一区二区三区四区| 少妇粗大呻吟视频| 久久久久久国产a免费观看| 久久性视频一级片| 色尼玛亚洲综合影院| 男女床上黄色一级片免费看| 日韩欧美国产一区二区入口| 黄色视频,在线免费观看| 在线a可以看的网站| 变态另类成人亚洲欧美熟女| 亚洲色图av天堂| 欧美日韩瑟瑟在线播放| 成人高潮视频无遮挡免费网站| 中文字幕熟女人妻在线| 亚洲成人中文字幕在线播放| 日韩中文字幕欧美一区二区| 日本熟妇午夜| 国产成+人综合+亚洲专区| 国产精品永久免费网站| 欧洲精品卡2卡3卡4卡5卡区| 成人亚洲精品av一区二区| 国产成人欧美在线观看| 欧美黄色淫秽网站| 老司机福利观看| 亚洲男人的天堂狠狠| 国内精品久久久久精免费| 国产精品亚洲一级av第二区| 深夜精品福利| 日韩有码中文字幕| 黑人操中国人逼视频| 欧美又色又爽又黄视频| 少妇裸体淫交视频免费看高清 | 丰满人妻熟妇乱又伦精品不卡| 亚洲一区中文字幕在线| 在线观看一区二区三区| 久久久久国内视频| 色尼玛亚洲综合影院| 亚洲专区国产一区二区| 中文字幕人妻丝袜一区二区| 亚洲成人久久爱视频| 国产精品野战在线观看| 亚洲 国产 在线| 哪里可以看免费的av片| 老司机午夜十八禁免费视频| 老司机深夜福利视频在线观看| 亚洲真实伦在线观看| 久久久精品欧美日韩精品| 99热6这里只有精品| 中文字幕av在线有码专区| 18禁观看日本| 久久久久久大精品| 国产成人av教育| 精品高清国产在线一区| www.自偷自拍.com| 极品教师在线免费播放| 欧美日韩亚洲国产一区二区在线观看| 在线观看免费视频日本深夜| 国产精品久久久久久精品电影| 精品久久久久久久末码| 精品国产乱子伦一区二区三区| 全区人妻精品视频| 88av欧美| 女警被强在线播放| 久久精品亚洲精品国产色婷小说| 丁香欧美五月| 久久精品成人免费网站| 看免费av毛片| 国产亚洲精品第一综合不卡| 亚洲aⅴ乱码一区二区在线播放 | 亚洲av电影在线进入| 精品一区二区三区四区五区乱码| 老司机深夜福利视频在线观看| 精华霜和精华液先用哪个| 亚洲精品中文字幕一二三四区| 精品午夜福利视频在线观看一区| 亚洲国产高清在线一区二区三| 18禁国产床啪视频网站| 午夜精品在线福利| 精品国产乱子伦一区二区三区| 一进一出抽搐gif免费好疼| 亚洲熟妇中文字幕五十中出| 桃色一区二区三区在线观看| 国产精品国产高清国产av| 色噜噜av男人的天堂激情| 99国产精品一区二区蜜桃av| 久久久久久久久免费视频了| www.自偷自拍.com| 人成视频在线观看免费观看| √禁漫天堂资源中文www| 久热爱精品视频在线9| 亚洲激情在线av| 在线观看免费日韩欧美大片| 国产亚洲精品综合一区在线观看 | 正在播放国产对白刺激| 久久精品人妻少妇| 国产高清有码在线观看视频 | 午夜福利在线在线| 一个人观看的视频www高清免费观看 | 国产麻豆成人av免费视频| svipshipincom国产片| 两性夫妻黄色片| 国产不卡一卡二| 一级a爱片免费观看的视频| 在线免费观看的www视频| 午夜亚洲福利在线播放| 好男人在线观看高清免费视频| 久久性视频一级片| 日本三级黄在线观看| 日韩大尺度精品在线看网址| 一级毛片高清免费大全| 成人18禁高潮啪啪吃奶动态图| 日日摸夜夜添夜夜添小说| 国内精品一区二区在线观看| 亚洲午夜精品一区,二区,三区| 999久久久精品免费观看国产| 亚洲一区二区三区色噜噜| 免费看a级黄色片| 草草在线视频免费看| www.999成人在线观看| 久久久国产精品麻豆| 中国美女看黄片| 欧美色欧美亚洲另类二区| 免费在线观看成人毛片| 精品国产乱子伦一区二区三区| 亚洲熟女毛片儿| 黄片小视频在线播放| 免费av毛片视频| 91大片在线观看| 国产久久久一区二区三区| www.熟女人妻精品国产| 老司机福利观看| 成人欧美大片| 欧美日韩精品网址| 婷婷精品国产亚洲av在线| 两性午夜刺激爽爽歪歪视频在线观看 | 国产av又大| 女警被强在线播放| 亚洲avbb在线观看| 99精品在免费线老司机午夜| 精品高清国产在线一区| 亚洲人成网站在线播放欧美日韩| 午夜福利在线在线| av国产免费在线观看| 黄片小视频在线播放| 国产黄色小视频在线观看| 久99久视频精品免费| 欧美乱码精品一区二区三区| 国产蜜桃级精品一区二区三区| 制服诱惑二区| 日韩精品中文字幕看吧| 波多野结衣高清作品| 夜夜夜夜夜久久久久| 激情在线观看视频在线高清| 哪里可以看免费的av片| 国产精品 欧美亚洲| 亚洲中文字幕一区二区三区有码在线看 | 最近视频中文字幕2019在线8| 熟女少妇亚洲综合色aaa.| 动漫黄色视频在线观看| 亚洲欧美一区二区三区黑人| 亚洲精品av麻豆狂野| 久久国产乱子伦精品免费另类| 一级片免费观看大全| 中文字幕久久专区| av欧美777| 国产欧美日韩一区二区精品| 成年版毛片免费区| 淫秽高清视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 日本黄色视频三级网站网址| 国产伦人伦偷精品视频| 国产探花在线观看一区二区| 欧美在线黄色| 黄色片一级片一级黄色片| 国内精品一区二区在线观看| 久久久国产精品麻豆| 亚洲精品国产一区二区精华液| 身体一侧抽搐| 色播亚洲综合网| 国产99久久九九免费精品| 久久久国产欧美日韩av| 亚洲人成电影免费在线| 757午夜福利合集在线观看| 国产精品一区二区三区四区久久| 亚洲avbb在线观看| 国产男靠女视频免费网站| 国模一区二区三区四区视频 | 波多野结衣高清无吗| 人成视频在线观看免费观看| 精品久久久久久久人妻蜜臀av| www日本黄色视频网| 亚洲色图 男人天堂 中文字幕| 女人爽到高潮嗷嗷叫在线视频| 婷婷六月久久综合丁香| 国产成人aa在线观看| 岛国在线观看网站| 嫩草影院精品99| 国产片内射在线| 波多野结衣高清无吗| 久久天躁狠狠躁夜夜2o2o| 日日摸夜夜添夜夜添小说| 一区福利在线观看| 国产精品美女特级片免费视频播放器 | 国产亚洲精品久久久久久毛片| 淫秽高清视频在线观看| 午夜精品一区二区三区免费看| 少妇人妻一区二区三区视频| 一进一出抽搐动态| 一区二区三区激情视频| 国产成人精品无人区| 国产激情欧美一区二区| 香蕉av资源在线| 很黄的视频免费| av福利片在线观看| 精品免费久久久久久久清纯| 国产成人av激情在线播放| 搡老岳熟女国产| 后天国语完整版免费观看| 99国产精品一区二区三区| 日本精品一区二区三区蜜桃| 特大巨黑吊av在线直播| 神马国产精品三级电影在线观看 | 波多野结衣高清无吗| 久久久久国内视频| 久久久久国产精品人妻aⅴ院| 国产av不卡久久| 熟女少妇亚洲综合色aaa.| 91字幕亚洲| 免费看十八禁软件| 亚洲最大成人中文| 极品教师在线免费播放| 国产99久久九九免费精品| 欧美黑人精品巨大| 欧美一级毛片孕妇| 欧美又色又爽又黄视频| 99国产精品一区二区蜜桃av| 中文字幕最新亚洲高清| 欧美成人午夜精品| 午夜免费成人在线视频| 久久香蕉国产精品| 给我免费播放毛片高清在线观看| 成人欧美大片| 搡老岳熟女国产| 亚洲国产欧美一区二区综合| 伦理电影免费视频| 精品国产超薄肉色丝袜足j| 国产欧美日韩精品亚洲av| 桃红色精品国产亚洲av| 亚洲国产中文字幕在线视频| 啦啦啦韩国在线观看视频| 午夜激情av网站| 亚洲国产精品999在线| 丰满人妻一区二区三区视频av | 亚洲va日本ⅴa欧美va伊人久久| 国产一级毛片七仙女欲春2| 欧美不卡视频在线免费观看 | 国产v大片淫在线免费观看| 午夜激情av网站| 他把我摸到了高潮在线观看| 日韩有码中文字幕| 一区二区三区激情视频| 十八禁人妻一区二区| 好看av亚洲va欧美ⅴa在| 亚洲av成人av| 亚洲欧美日韩高清在线视频| 国产免费av片在线观看野外av| 亚洲一区高清亚洲精品| 一区二区三区高清视频在线| 真人一进一出gif抽搐免费| 一级黄色大片毛片| 中国美女看黄片| 日韩国内少妇激情av| 精品久久久久久久毛片微露脸| 精品一区二区三区av网在线观看| 亚洲av美国av| 国产在线精品亚洲第一网站| 99热只有精品国产| 最近最新中文字幕大全电影3| 国产精品香港三级国产av潘金莲| 午夜福利成人在线免费观看| 午夜福利欧美成人| 熟女少妇亚洲综合色aaa.| 村上凉子中文字幕在线| 两人在一起打扑克的视频| 12—13女人毛片做爰片一| 青草久久国产| 在线a可以看的网站| 亚洲国产精品成人综合色| 欧美乱色亚洲激情| 免费观看人在逋| 久久久久久久久中文| 黄色视频不卡| 精品国产乱子伦一区二区三区| 1024香蕉在线观看| www.熟女人妻精品国产| 欧美午夜高清在线| 91大片在线观看| 男女做爰动态图高潮gif福利片| 中出人妻视频一区二区| 啦啦啦免费观看视频1| 国产97色在线日韩免费| 在线视频色国产色| 国产v大片淫在线免费观看| 色在线成人网| 精品欧美一区二区三区在线| 国产一区二区在线观看日韩 | 国产成人精品无人区| 一本大道久久a久久精品| 国产精品自产拍在线观看55亚洲| 国产av麻豆久久久久久久| 长腿黑丝高跟| 99精品在免费线老司机午夜| 亚洲中文av在线| 亚洲av成人精品一区久久| 亚洲激情在线av| 国产av一区在线观看免费| 亚洲美女视频黄频| 亚洲第一欧美日韩一区二区三区| 成人三级黄色视频| 亚洲成人久久爱视频| 久热爱精品视频在线9| cao死你这个sao货| 好看av亚洲va欧美ⅴa在| 日韩欧美 国产精品| 岛国在线免费视频观看| 在线国产一区二区在线| 俄罗斯特黄特色一大片| 亚洲一区二区三区不卡视频| 国产av不卡久久| 看免费av毛片| 亚洲熟妇熟女久久| 久久午夜综合久久蜜桃| 十八禁网站免费在线| 精品久久久久久久久久免费视频| 男女午夜视频在线观看| 欧美中文日本在线观看视频| 无遮挡黄片免费观看| 久久性视频一级片| 一区二区三区高清视频在线| 日本一区二区免费在线视频| 欧美日韩黄片免| 亚洲欧美日韩高清专用| 我的老师免费观看完整版| 99国产精品一区二区蜜桃av| 亚洲欧美精品综合一区二区三区| 午夜福利18| 99精品久久久久人妻精品| 两个人的视频大全免费| 18美女黄网站色大片免费观看| 99国产精品一区二区三区| а√天堂www在线а√下载| 亚洲精品av麻豆狂野| 国产亚洲欧美在线一区二区| 精品高清国产在线一区| 亚洲av成人av| 看黄色毛片网站| 国产一区二区三区在线臀色熟女| 可以在线观看毛片的网站| 哪里可以看免费的av片| xxx96com| 欧美午夜高清在线| 免费观看精品视频网站| 香蕉av资源在线| 黑人欧美特级aaaaaa片| 天堂√8在线中文| 国产人伦9x9x在线观看| www.999成人在线观看| 两人在一起打扑克的视频| 99久久久亚洲精品蜜臀av| 精品熟女少妇八av免费久了| 777久久人妻少妇嫩草av网站| 老鸭窝网址在线观看| 一级黄色大片毛片| 母亲3免费完整高清在线观看| 又紧又爽又黄一区二区| 青草久久国产| 日韩中文字幕欧美一区二区| 90打野战视频偷拍视频| 精品一区二区三区av网在线观看| 欧美黄色片欧美黄色片| 麻豆一二三区av精品| 久久久久久大精品| 91大片在线观看| 一区福利在线观看| 久久人妻av系列| 一进一出好大好爽视频| 级片在线观看| 亚洲av成人av| 国产亚洲精品综合一区在线观看 | 国产在线精品亚洲第一网站| 一二三四社区在线视频社区8| 精品久久久久久久久久久久久| 国产精品一及| 久久人妻av系列| 香蕉av资源在线| 人妻丰满熟妇av一区二区三区| 午夜福利在线在线| 欧美3d第一页| av在线播放免费不卡| 久久久国产成人精品二区| 一a级毛片在线观看| 国产片内射在线| avwww免费| 正在播放国产对白刺激| 波多野结衣高清作品| 黑人欧美特级aaaaaa片| 91成年电影在线观看| 99久久99久久久精品蜜桃| 久久香蕉国产精品| 老司机深夜福利视频在线观看| 亚洲一码二码三码区别大吗| 老熟妇乱子伦视频在线观看| 亚洲成av人片免费观看| 国产高清视频在线播放一区| 制服诱惑二区| 此物有八面人人有两片| 老司机午夜福利在线观看视频| 国产三级在线视频| svipshipincom国产片| 免费电影在线观看免费观看| 中文在线观看免费www的网站 | 香蕉国产在线看| 亚洲自偷自拍图片 自拍| 国内揄拍国产精品人妻在线| 麻豆国产av国片精品| 村上凉子中文字幕在线| 男女视频在线观看网站免费 | 又粗又爽又猛毛片免费看| 男人的好看免费观看在线视频 | 91在线观看av| 国产精品av视频在线免费观看| 国产69精品久久久久777片 | 两个人看的免费小视频| 久久国产精品影院| 日韩三级视频一区二区三区| 97超级碰碰碰精品色视频在线观看| 国产精品国产高清国产av| 成人午夜高清在线视频| 精品一区二区三区视频在线观看免费| 一进一出抽搐动态| 亚洲 国产 在线| 精品欧美一区二区三区在线| 久久久久久久久久黄片| 特大巨黑吊av在线直播| 一a级毛片在线观看| 欧美久久黑人一区二区| 国产精品爽爽va在线观看网站| 免费在线观看日本一区| 国产黄色小视频在线观看| 又紧又爽又黄一区二区| 桃红色精品国产亚洲av| 亚洲精品在线观看二区| 欧美一级毛片孕妇| 国产精品香港三级国产av潘金莲| 欧美一级a爱片免费观看看 | 国产日本99.免费观看| avwww免费| 色av中文字幕| 天堂√8在线中文| 特级一级黄色大片| 两个人视频免费观看高清| 女同久久另类99精品国产91| 免费在线观看黄色视频的| 精品一区二区三区视频在线观看免费| 亚洲欧美激情综合另类| 国产精品av视频在线免费观看| 日韩欧美在线二视频| 国产av不卡久久| 欧美 亚洲 国产 日韩一| 欧美日韩亚洲综合一区二区三区_| 成年版毛片免费区| 99精品久久久久人妻精品| www.自偷自拍.com| svipshipincom国产片| 女生性感内裤真人,穿戴方法视频| 午夜精品一区二区三区免费看| 叶爱在线成人免费视频播放| 中文亚洲av片在线观看爽| 色播亚洲综合网| 亚洲中文字幕一区二区三区有码在线看 | 黄片大片在线免费观看| 精品电影一区二区在线| 18禁黄网站禁片午夜丰满| 午夜免费成人在线视频| 欧美黄色淫秽网站| 亚洲精品在线观看二区| 亚洲欧洲精品一区二区精品久久久| 国产高清视频在线观看网站| 久久伊人香网站| 久久精品国产清高在天天线| 淫妇啪啪啪对白视频| 亚洲国产欧美人成| 免费无遮挡裸体视频| 日韩欧美免费精品| 午夜两性在线视频| 欧美黄色淫秽网站| 宅男免费午夜| 国产成人av激情在线播放| 国内毛片毛片毛片毛片毛片| 亚洲国产高清在线一区二区三| 真人一进一出gif抽搐免费| 国产成人系列免费观看| 啦啦啦观看免费观看视频高清| 色综合亚洲欧美另类图片| 国产激情偷乱视频一区二区| 母亲3免费完整高清在线观看| 亚洲一区二区三区不卡视频| 成人午夜高清在线视频| 久久久久久人人人人人| 午夜福利在线观看吧| 亚洲熟女毛片儿| 国产区一区二久久| 国产精品,欧美在线| 91老司机精品| 日日干狠狠操夜夜爽| 悠悠久久av| 精品一区二区三区四区五区乱码| 久久天堂一区二区三区四区| 亚洲av成人一区二区三| 九九热线精品视视频播放| 777久久人妻少妇嫩草av网站| 国产麻豆成人av免费视频| 亚洲人成网站在线播放欧美日韩| 又大又爽又粗| 女同久久另类99精品国产91| 亚洲欧美一区二区三区黑人| 黑人操中国人逼视频| av在线播放免费不卡| 国产精品亚洲av一区麻豆| 亚洲中文字幕一区二区三区有码在线看 | 欧美性猛交╳xxx乱大交人| 国产亚洲av高清不卡| 制服丝袜大香蕉在线| 日韩欧美在线乱码| 日本a在线网址| 亚洲精品一卡2卡三卡4卡5卡| av福利片在线| 亚洲自拍偷在线| 国产激情久久老熟女| 午夜免费观看网址| 国产免费男女视频| av福利片在线| 一进一出好大好爽视频| 老司机午夜福利在线观看视频| 久久精品影院6| 亚洲色图av天堂| aaaaa片日本免费| 日韩 欧美 亚洲 中文字幕| 国产1区2区3区精品| 欧美成人一区二区免费高清观看 | 每晚都被弄得嗷嗷叫到高潮| 可以免费在线观看a视频的电影网站| 少妇人妻一区二区三区视频| 日日夜夜操网爽| 中国美女看黄片| 免费在线观看视频国产中文字幕亚洲| 精品欧美国产一区二区三| av国产免费在线观看| 成人欧美大片| 久久精品国产亚洲av高清一级| 999久久久国产精品视频| 精品久久久久久久久久久久久| 国产三级黄色录像| 三级男女做爰猛烈吃奶摸视频| aaaaa片日本免费| 午夜影院日韩av| 夜夜躁狠狠躁天天躁| 男女视频在线观看网站免费 | 制服诱惑二区| av在线播放免费不卡| netflix在线观看网站| 99国产综合亚洲精品| 两个人的视频大全免费| www.熟女人妻精品国产| 国产精品av久久久久免费| 女同久久另类99精品国产91| 久久草成人影院| 亚洲精品国产精品久久久不卡| 熟妇人妻久久中文字幕3abv| 国产欧美日韩一区二区精品| 黄色女人牲交| 成人永久免费在线观看视频| 中国美女看黄片| 18禁黄网站禁片免费观看直播| 免费一级毛片在线播放高清视频| 成在线人永久免费视频| 12—13女人毛片做爰片一|