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

        ?

        基于偏相干分析方法的怠速工況車內(nèi)噪聲源識(shí)別

        2019-07-04 02:12:56姜少瑋李孟宇
        關(guān)鍵詞:噪聲源頻段測點(diǎn)

        陳 克,姜少瑋,李孟宇

        (沈陽理工大學(xué) 汽車與交通學(xué)院,沈陽 110159)

        對(duì)于汽車行業(yè)而言,消費(fèi)者對(duì)整車品質(zhì)的要求嚴(yán)苛,除了汽車的經(jīng)濟(jì)性、動(dòng)力性等,消費(fèi)者更看重的是駕駛和乘坐舒適性.據(jù)統(tǒng)計(jì),汽車故障投訴中有1/3的問題和噪聲、振動(dòng)與聲振粗糙度(Noise,Vibration,Harshness,NVH)控制有關(guān)[1].汽車車內(nèi)噪聲大小作為影響駕乘人員乘坐舒適性的重要指標(biāo),不僅是提高整車品質(zhì)的關(guān)鍵,更是企業(yè)核心競爭力的體現(xiàn)[2].因此,有效控制車內(nèi)噪聲水平勢在必行.

        為了有效控制汽車車內(nèi)噪聲,國內(nèi)外研究人員針對(duì)噪聲源識(shí)別問題提出的理論方法可以歸納為傳統(tǒng)噪聲源識(shí)別方法、基于信號(hào)處理技術(shù)的分析方法以及基于可視化技術(shù)的聲源識(shí)別方法.傳統(tǒng)噪聲源識(shí)別方法和基于可視化技術(shù)的聲源識(shí)別方法,在簡單聲源識(shí)別和表面噪聲識(shí)別上的應(yīng)用較廣.當(dāng)實(shí)際應(yīng)用于汽車這樣一個(gè)存在著多個(gè)噪聲源且包含多個(gè)振動(dòng)發(fā)聲部件、聲固耦合性復(fù)雜的系統(tǒng)時(shí),有識(shí)別精度低、測試成本高的劣勢;而基于信號(hào)處理技術(shù)的識(shí)別方法,在識(shí)別精度、測試手段和普適性上優(yōu)勢明顯[3-4].

        相干分析法識(shí)別車內(nèi)噪聲源屬于數(shù)字信號(hào)處理類.該方法是以相關(guān)理論為基礎(chǔ),以相干函數(shù)為分析工具,本文主要采用偏相干函數(shù),分析頻域內(nèi)各噪聲源的輸入信號(hào)與響應(yīng)點(diǎn)的輸出信號(hào)間的相干系數(shù)大小[5].偏相干函數(shù)是建立在常相干理論基礎(chǔ)上的[6],可視為特殊的常相關(guān)函數(shù),兩者區(qū)別于偏相干分析可以在頻域范圍內(nèi)計(jì)算系數(shù)中各非獨(dú)立輸入對(duì)輸出的線性影響系數(shù),而常相干分析在此方面的應(yīng)用存在局限性[7].在對(duì)振動(dòng)噪聲信號(hào)源因果關(guān)系判斷正確的前提下進(jìn)行有序輸入,偏相干分析能剔除其他相關(guān)輸入的影響,提高在實(shí)際應(yīng)用中車輛噪聲源識(shí)別精確度[8].

        本文在分析汽車振動(dòng)噪聲產(chǎn)生機(jī)理的基礎(chǔ)上,進(jìn)行車內(nèi)噪聲源測試試驗(yàn)以及相干性分析技術(shù)理論研究,結(jié)合譜分析識(shí)別出結(jié)構(gòu)-空氣噪聲源,以噪聲源類別判斷為基礎(chǔ)進(jìn)行信號(hào)源因果關(guān)系判斷,應(yīng)用偏相干分析法準(zhǔn)確識(shí)別引發(fā)車內(nèi)噪聲的噪聲源.

        1 測試目標(biāo)確定及試驗(yàn)設(shè)計(jì)

        1.1 測試目標(biāo)分析

        行駛中的汽車車內(nèi)噪聲源主要來自兩方面,分別是汽車自身工作條件方面和汽車外部環(huán)境方面,而由汽車外部環(huán)境導(dǎo)致的噪聲(如路噪、風(fēng)噪等)在進(jìn)行汽車設(shè)計(jì)生產(chǎn)過程中無法逐一實(shí)地考量進(jìn)而控制.對(duì)于汽車自身工作條件所引起的噪聲問題,為避免外部環(huán)境對(duì)汽車發(fā)動(dòng)機(jī)懸置振動(dòng)的影響,就需要進(jìn)行汽車在定置怠速工況基礎(chǔ)條件下的車內(nèi)噪聲源識(shí)別研究.基于此,本文所研究的影響定置工況下,車內(nèi)噪聲的主要噪聲源有:發(fā)動(dòng)機(jī)懸置系統(tǒng)、進(jìn)氣系統(tǒng)、排氣系統(tǒng)以及發(fā)動(dòng)機(jī)艙的共振.怠速工況下車內(nèi)噪聲源按結(jié)構(gòu)路徑和空氣路徑傳播如圖1所示.

        圖1 定置怠速工況下車內(nèi)噪聲源傳播示意圖Fig.1 Schematic diagram of noise source transmission

        1.2 車內(nèi)噪聲源測試試驗(yàn)設(shè)計(jì)

        對(duì)某國產(chǎn)SUV車進(jìn)行振動(dòng)與噪聲測試,試驗(yàn)方法依據(jù)GB/T 18697—2002《聲學(xué)汽車車內(nèi)噪聲測量方法》[9]進(jìn)行.在定置怠速工況下,以發(fā)動(dòng)機(jī)懸置被動(dòng)側(cè)的振動(dòng)加速度、進(jìn)氣系統(tǒng)噪聲、排氣系統(tǒng)噪聲和發(fā)動(dòng)機(jī)艙噪聲作為系統(tǒng)模型中的輸入信號(hào),駕駛員右耳處聲壓信號(hào)作為系統(tǒng)輸出信號(hào),一共布置7個(gè)振動(dòng)與噪聲測點(diǎn),共需13個(gè)測試通道.

        試驗(yàn)主要儀器和設(shè)備包括:某國產(chǎn)SUV車、32通道LMS振動(dòng)噪聲數(shù)據(jù)采集儀、美國PCB三向加速度傳感器、丹麥GRAS傳聲器等,測試對(duì)象的振動(dòng)與噪聲各測點(diǎn)具體布置如圖2所示.圖2中測點(diǎn)1~3安裝三向加速度傳感器,測點(diǎn)4~7安裝傳聲器.

        圖2 振動(dòng)與噪聲測點(diǎn)布置圖Fig.2 Vibration and noise measuring points layout

        2 車內(nèi)噪聲源測試數(shù)據(jù)分析

        基于試驗(yàn)設(shè)計(jì)方案以及相干性分析的技術(shù)理論,進(jìn)行試驗(yàn)數(shù)據(jù)處理,相干分析的研究是以相關(guān)理論為基礎(chǔ),以相干函數(shù)為分析工具,其流程如圖3所示.

        2.1 時(shí)域數(shù)據(jù)分析

        在測試過程中試驗(yàn)車輛手剎處于制動(dòng)狀態(tài),變速器位于空擋,將試驗(yàn)車輛預(yù)熱5 min,待各系統(tǒng)處于平穩(wěn)運(yùn)行狀態(tài)時(shí)開始進(jìn)行數(shù)據(jù)采集.在定置怠速800 r/min工況下,采樣時(shí)間10 s內(nèi)測得的車內(nèi)噪聲源Overall曲線如圖4所示.

        上述分析從時(shí)域的角度對(duì)懸置點(diǎn)振動(dòng)和空氣噪聲聲壓進(jìn)行了時(shí)域特征分析,但是結(jié)構(gòu)振動(dòng)和空氣噪聲的能量、幅值的大小,不能反映與車內(nèi)噪聲的關(guān)系,只能反映該時(shí)刻下的幅值變化,不能以此確定噪聲源,試驗(yàn)采集的各測點(diǎn)振動(dòng)與噪聲數(shù)據(jù)可以用來轉(zhuǎn)換到頻域下繼續(xù)進(jìn)行分析.

        OA曲線是對(duì)整個(gè)采樣信號(hào)在時(shí)域上的線性平均,可以看到在10 s的采樣時(shí)間內(nèi)怠速800 r/min的車內(nèi)噪聲源輸入端信號(hào)處于穩(wěn)定狀態(tài),符合怠速工況的實(shí)際情況,數(shù)據(jù)采集可靠,整個(gè)時(shí)段的采樣信號(hào)可以作為穩(wěn)態(tài)數(shù)據(jù)轉(zhuǎn)換到頻域下繼續(xù)進(jìn)行分析.

        圖3 試驗(yàn)數(shù)據(jù)處理流程圖Fig.3 Test data processing flow chart

        圖4 車內(nèi)噪聲源Overall Level曲線Fig.4 Interior noise source overall level curve

        2.2 相干性分析技術(shù)理論

        2.2.1多輸入單輸出系統(tǒng)模型

        相干性分析中,假設(shè)噪聲源之間相互獨(dú)立,建立多輸入單輸出系統(tǒng)模型[10],如圖5所示.圖5中:xn(t)為噪聲源輸入信號(hào);Hny(f)為由輸入到輸出的傳遞函數(shù);yn(t)為常參數(shù)理想線性系統(tǒng)的輸出;n(t)為偏離線性系統(tǒng)的所有偏差;y(t)為偏離線性系統(tǒng)的所有偏差n(t)與理想輸出yn(t)之和.

        圖5 多輸入單輸出系統(tǒng)模型Fig.5 Multi-input and single output system model

        常參數(shù)理想線性系統(tǒng)中相干函數(shù)定義為[7]

        (1)

        2.2.2偏相干分析計(jì)算

        偏相干分析由常相干分析的基礎(chǔ)發(fā)展而來,能克服常相干函數(shù)在應(yīng)用條件下的限制.在多輸入單輸出系統(tǒng)中,若輸入信號(hào)間含有相干成分,偏相干分析可以預(yù)先將輸入信號(hào)與其他信號(hào)的相干因素去掉,運(yùn)用條件功率譜計(jì)算來排除信號(hào)源的相干影響,考慮圖5所示模型中n個(gè)任意輸入xi(t),(i=1,2,…,n)和一個(gè)輸出y(t)的情況,所有這(n+1)個(gè)信號(hào)均可測.為了方便描述,現(xiàn)規(guī)定記錄次序如下:去掉前r個(gè)輸入信號(hào),剩下的是后面的(n-r)個(gè)輸入和輸出y(t).條件記錄記為

        以上每個(gè)條件記錄的條件自譜Sii·r!和Syy·r!,條件互譜Siy·r!均可計(jì)算,即

        其中,最優(yōu)系統(tǒng)Liy的計(jì)算公式為

        (6)

        若用y(t)=xn+1(t)來代替,則Liy變成Lij,得

        (7)

        對(duì)任意i,有Lii=1;任意j

        (8)

        (9)

        實(shí)際車輛工作系統(tǒng)振動(dòng)與噪聲的傳遞是耦合的、非獨(dú)立的[11].偏相干計(jì)算引入剩余譜的概念[12],能解決實(shí)際工程情況中多輸入信號(hào)間非獨(dú)立的問題,可以在各輸入信號(hào)間消除相關(guān)影響后計(jì)算各輸入對(duì)輸出的貢獻(xiàn),即在剩余譜應(yīng)用的基礎(chǔ)上建立起的多輸入單輸出系統(tǒng)模型中輸入信號(hào)間等效于滿足獨(dú)立性.模型理論應(yīng)用條件為:假定噪聲源信號(hào)為平穩(wěn)隨機(jī)信號(hào),系統(tǒng)為常參數(shù)線性系統(tǒng).偏相干計(jì)算的系統(tǒng)模型如圖6所示.

        圖6 基于剩余譜的多輸入單輸出系統(tǒng)模型 Fig.6 Multi input and single output system model

        圖6是用一組有序的條件輸入記錄代替原始的已知輸入記錄,按圖6所示次序選擇的條件輸入記錄為{xi·(i-1)!},i=1,2,…,n,對(duì)于任意i,xi·(i-1)!表示前x1,x2,…,xi-1條件下的xi,即從xi中去掉x1到xi-1的線性影響后的測量信號(hào),這些有序的條件輸入是兩兩不相關(guān)的,偏相干的系統(tǒng)模型解決了相干性理論應(yīng)用中對(duì)輸入信號(hào)相關(guān)性的限制.這種有序輸入就是輸入信號(hào)與輸出信號(hào)之間的因果關(guān)系判斷,在國外文獻(xiàn)中提及的有常相干函數(shù)法、發(fā)展的脈沖響應(yīng)函數(shù)法和Hilbert變換法[13-14].文獻(xiàn)[15]中論證了以上方法的缺陷,區(qū)別于以上方法,本文提出基于結(jié)構(gòu)和空氣噪聲源特性分析法的因果關(guān)系判斷:在判斷噪聲源類型的基礎(chǔ)上,若某頻段或工況下結(jié)構(gòu)噪聲為主,則該結(jié)構(gòu)噪聲源優(yōu)先級(jí)越高;若空氣噪聲為主,則空氣噪聲源優(yōu)先級(jí)越高.

        2.3 車內(nèi)噪聲源特性分析

        振動(dòng)與噪聲測點(diǎn)的頻譜分析,可以在不切斷噪聲源與接收者關(guān)聯(lián)的情況下直接反映測試結(jié)果,在定置怠速800 r/min工況下,采用發(fā)動(dòng)機(jī)左、右、后3個(gè)被動(dòng)側(cè)懸置點(diǎn)振動(dòng)加速度頻譜圖(見圖7)和各系統(tǒng)噪聲頻譜圖(見圖8)作為輔助分析圖.

        圖7 被動(dòng)側(cè)懸置點(diǎn)振動(dòng)加速度頻譜圖Fig.7 Passive side mount point acceleration spectrum

        圖8 排氣、進(jìn)氣、發(fā)動(dòng)機(jī)艙噪聲頻譜圖Fig.8 Exhaust,intake,engine room noise spectrum

        由圖7可知,發(fā)動(dòng)機(jī)懸置的振動(dòng)測量點(diǎn)在1 000 Hz以下有明顯的峰值,300 Hz以下的低頻能量最強(qiáng).由圖8可知,空氣系統(tǒng)噪聲在500 Hz以上有明顯的峰值,其中1 000 Hz以上的高頻能量最強(qiáng).圖7中在0~1 000 Hz頻段曲線峰值多于圖8對(duì)應(yīng)頻段峰值,且能量較強(qiáng);圖8中1 000~2 000 Hz頻段曲線峰值多于圖7對(duì)應(yīng)頻段峰值,且能量較強(qiáng).由此判斷汽車怠速工況下車內(nèi)噪聲源產(chǎn)生的原因如下:結(jié)構(gòu)振動(dòng)和空氣噪聲相互耦合共同作用,但0~1 000 Hz頻段處主要是發(fā)動(dòng)機(jī)的激勵(lì)作用經(jīng)懸置點(diǎn)傳遞至車身引起結(jié)構(gòu)壁板振動(dòng),在駕駛室內(nèi)形成輻射噪聲,結(jié)構(gòu)噪聲為主;在1 000~2 000 Hz頻段處,進(jìn)氣噪聲、排氣噪聲和發(fā)動(dòng)機(jī)艙噪聲通過空氣傳播經(jīng)車身孔隙進(jìn)入車內(nèi),空氣噪聲為主.

        3 偏相干分析的車內(nèi)噪聲源識(shí)別

        車內(nèi)噪聲源特性分析識(shí)別出的結(jié)構(gòu)和空氣噪聲源判別結(jié)果,解決了偏相干分析信號(hào)源的因果關(guān)系判斷問題.被測車輛的橡膠懸置實(shí)際為非線性系統(tǒng),但測試條件為定置怠速工況,可以將其視為弱線性系統(tǒng),數(shù)據(jù)采集是在正態(tài)平穩(wěn)隨機(jī)過程下進(jìn)行的,根據(jù)相干性理論,可以進(jìn)行線性系統(tǒng)分析.

        以結(jié)構(gòu)-空氣噪聲源識(shí)別試驗(yàn)的6個(gè)噪聲源測點(diǎn)作為輸入,駕駛員右耳處噪聲作為輸出,建立一個(gè)6輸入單輸出的線性系統(tǒng),如圖9所示.根據(jù)車內(nèi)噪聲源特性分析的因果關(guān)系判斷,當(dāng)研究0~1 000 Hz 頻段時(shí),結(jié)構(gòu)噪聲源因果關(guān)系大、優(yōu)先級(jí)高,則將發(fā)動(dòng)機(jī)振動(dòng)測點(diǎn)信號(hào)放置首位,其余任意排列;當(dāng)研究1 000~2 000 Hz頻段時(shí),空氣噪聲源因果關(guān)系大、優(yōu)先級(jí)高,將頻譜峰值能量高,頻譜峰值密度大的放置首位,其余任意排列;當(dāng)研究其他頻段時(shí),同理分析.

        圖9 6輸入單輸出系統(tǒng)模型Fig.9 Six-input single output system

        將LMS Test.lab中的實(shí)測振動(dòng)與噪聲數(shù)據(jù)處理后應(yīng)用偏相干程序在Matlab軟件中計(jì)算,計(jì)算怠速800 r/min工況下,3個(gè)發(fā)動(dòng)機(jī)懸置被動(dòng)側(cè)振動(dòng)測點(diǎn)和進(jìn)氣系統(tǒng)、排氣系統(tǒng)、發(fā)動(dòng)機(jī)艙的噪聲測點(diǎn)偏相干系數(shù)曲線,如圖10所示.

        由圖10可知:在0~1 000 Hz頻段內(nèi),振動(dòng)與噪聲對(duì)車內(nèi)噪聲貢獻(xiàn)較大,其中,峰值最大的為發(fā)動(dòng)機(jī)右懸置,其偏相干系數(shù)約為0.91,其次是發(fā)動(dòng)機(jī)后懸置和進(jìn)氣系統(tǒng),振動(dòng)測點(diǎn)貢獻(xiàn)明顯偏高.說明怠速工況下發(fā)動(dòng)機(jī)低速運(yùn)轉(zhuǎn),低頻段發(fā)動(dòng)機(jī)振動(dòng)經(jīng)過懸置點(diǎn)傳遞至車身及車架引起車內(nèi)結(jié)構(gòu)壁板振動(dòng)是造成車內(nèi)噪聲的主因,該頻段范圍內(nèi)發(fā)動(dòng)機(jī)振動(dòng)為主要噪聲源,其中經(jīng)發(fā)動(dòng)機(jī)右懸置傳遞的貢獻(xiàn)量最大.

        圖10 振動(dòng)與噪聲測點(diǎn)偏相干曲線圖Fig.10 Partial coherent curve of measuring points of

        在1 000~2 000 Hz頻段內(nèi),峰值最大的為發(fā)動(dòng)機(jī)左懸置,其偏相干系數(shù)約為0.90;其次是發(fā)動(dòng)機(jī)艙和發(fā)動(dòng)機(jī)后懸置,頻段內(nèi)噪聲達(dá)到較高的原因是空氣噪聲與結(jié)構(gòu)振動(dòng)的耦合,兩者共同作用形成的噪聲泄露至駕駛室內(nèi),導(dǎo)致車內(nèi)噪聲壓力值增大,其中經(jīng)發(fā)動(dòng)機(jī)左懸置傳遞的貢獻(xiàn)量最大.

        在2 000~3 000 Hz頻段內(nèi),峰值最大的為發(fā)動(dòng)機(jī)艙,其偏相干系數(shù)約為0.75;其次是排氣系統(tǒng)和發(fā)動(dòng)機(jī)右懸置,振動(dòng)測點(diǎn)的部分貢獻(xiàn)總體上下降,空氣噪聲對(duì)駕駛室內(nèi)噪聲的貢獻(xiàn)增強(qiáng),表明發(fā)動(dòng)機(jī)的燃燒噪聲和排氣系統(tǒng)的噪聲通過空氣泄漏到駕駛室,其中發(fā)動(dòng)機(jī)艙是主要的噪聲源.

        結(jié)合以上分析可知:定置怠速工況下振動(dòng)與噪聲在各頻段內(nèi)均有耦合作用,造成怠速工況下車內(nèi)噪聲的主要噪聲源是發(fā)動(dòng)機(jī)系統(tǒng),發(fā)動(dòng)機(jī)右懸置的振動(dòng)傳遞至車內(nèi)產(chǎn)生的結(jié)構(gòu)噪聲貢獻(xiàn)最大.表明基于結(jié)構(gòu)-空氣噪聲源特性分析的因果關(guān)系判斷下進(jìn)行信號(hào)有序的條件輸入,應(yīng)用偏相干分析在不同頻率上能準(zhǔn)確識(shí)別出主要噪聲源,在解決實(shí)際工程問題應(yīng)用中效果顯著.

        4 結(jié)論

        通過對(duì)被測車輛振動(dòng)與噪聲測點(diǎn)的實(shí)時(shí)數(shù)據(jù)采集并處理,應(yīng)用偏相干性分析得出如下結(jié)論.

        (1) 車內(nèi)噪聲在全頻域內(nèi)均受結(jié)構(gòu)噪聲和空氣噪聲的耦合作用,在基于結(jié)構(gòu)-空氣噪聲源特性分析基礎(chǔ)上解決了信號(hào)源的因果關(guān)系判斷,做到輸入信號(hào)的有序輸入,應(yīng)用相干性分析確定影響車內(nèi)噪聲的主要噪聲源.

        (2) 針對(duì)某國產(chǎn)SUV車型,在定置怠速工況下進(jìn)行振動(dòng)與噪聲測點(diǎn)實(shí)時(shí)采集數(shù)據(jù),應(yīng)用相干性理論識(shí)別噪聲源分析得出:在0~1 000 Hz頻段,發(fā)動(dòng)機(jī)振動(dòng)為主要噪聲源,其中峰值最大的為發(fā)動(dòng)機(jī)右懸置,其偏相干系數(shù)約為0.91,是主要貢獻(xiàn)源;在1 000~2 000 Hz頻段,噪聲達(dá)到較高的原因是空氣噪聲與結(jié)構(gòu)振動(dòng)的耦合,其中峰值最大的為發(fā)動(dòng)機(jī)左懸置,其偏相干系數(shù)約為0.90,是主要貢獻(xiàn)源;在2 000~3 000 Hz頻段,振動(dòng)測點(diǎn)貢獻(xiàn)量下降,空氣噪聲的貢獻(xiàn)增強(qiáng),其中峰值最大的為發(fā)動(dòng)機(jī)艙,其相干系數(shù)約為0.75,是主要貢獻(xiàn)源.根據(jù)偏相干系數(shù)大小可知,0~1 000 Hz頻段的偏相干系數(shù)明顯高于其他頻段的偏相干系數(shù),這是車輛在定置怠速工況下的主要噪聲頻段.由此可以驗(yàn)證應(yīng)用偏相干分析方法可以很好地進(jìn)行車內(nèi)噪聲源識(shí)別.

        猜你喜歡
        噪聲源頻段測點(diǎn)
        液壓支架整機(jī)靜強(qiáng)度試驗(yàn)及等效應(yīng)力分析
        基于CATIA的汽車測點(diǎn)批量開發(fā)的研究與應(yīng)用
        gPhone重力儀的面波頻段響應(yīng)實(shí)測研究
        地震研究(2021年1期)2021-04-13 01:04:56
        汽車后視鏡-A柱區(qū)域氣動(dòng)噪聲源特征識(shí)別
        汽車工程(2021年12期)2021-03-08 02:34:18
        推擠的5GHz頻段
        CHIP新電腦(2016年3期)2016-03-10 14:07:52
        一種基于相位增量隨機(jī)化的寬帶噪聲源產(chǎn)生技術(shù)
        TD—LTE在D頻段和F頻段的覆蓋能力差異
        中國新通信(2015年1期)2015-05-30 10:30:46
        拱壩結(jié)構(gòu)損傷的多測點(diǎn)R/S分析
        利用源強(qiáng)聲輻射模態(tài)識(shí)別噪聲源
        基于電流噪聲源模型的EMI濾波器設(shè)計(jì)
        激情文学婷婷六月开心久久| 中文字幕日韩精品无码内射| 日本丰满人妻xxxxxhd| 精品一精品国产一级毛片| 少妇裸淫交视频免费看| 91九色人妻精品一区二区三区| 欧洲乱码伦视频免费| 久久精品人妻无码一区二区三区| 亚洲一区中文字幕在线电影网 | 热久久网站| 日本人妖一区二区三区| 日本人妻系列中文字幕| 在线播放免费人成毛片乱码| 欧美aaaaaa级午夜福利视频 | 小宝极品内射国产在线| 免费无码中文字幕A级毛片| 四虎在线中文字幕一区| 国产自拍视频在线观看免费 | 麻豆激情视频在线观看| 中文字幕影片免费人妻少妇 | 国产精品美女久久久免费| 国产精品嫩草影院av| 国产在线观看入口| 欧美h久免费女| 视频国产一区二区在线| 精品国产一区二区三区三| 精品国产一区二区三区av片| 亚洲AV激情一区二区二三区| 中文字幕一区二区在线看| 色又黄又爽18禁免费网站现观看| 欧美性猛交内射兽交老熟妇| 免费大学生国产在线观看p| 国内嫩模自拍偷拍视频| 中国老太婆bb无套内射| 波多野结衣aⅴ在线| 精品少妇人妻久久免费| 日本少妇又色又爽又高潮| 国产真实强被迫伦姧女在线观看| 亚洲中文无码成人影院在线播放| 中文字幕日本韩国精品免费观看 | 一本久道综合在线无码人妻|