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

    高速列車轉(zhuǎn)向架艙對(duì)轉(zhuǎn)向架區(qū)域流場(chǎng)與氣動(dòng)噪聲影響

    2018-12-04 08:27:00朱劍月任利惠雷震宇
    關(guān)鍵詞:偶極子車軸聲源

    朱劍月, 任利惠, 雷震宇

    (1.同濟(jì)大學(xué) 鐵道與城市軌道交通研究院,上海 201804;2.同濟(jì)大學(xué) 上海地面交通工具風(fēng)洞中心,上海 201804)

    當(dāng)高速列車運(yùn)行速度超過(guò)300 km·h-1時(shí),氣動(dòng)噪聲將占據(jù)主導(dǎo)地位[1-2].研究表明,轉(zhuǎn)向架部位存在著表面結(jié)構(gòu)變化的不連續(xù)性,氣流流經(jīng)時(shí)在該部位產(chǎn)生流動(dòng)分離和較強(qiáng)的流體相互作用,從而形成氣動(dòng)噪聲[2-3].已有高速列車轉(zhuǎn)向架部位氣動(dòng)噪聲研究多采用現(xiàn)場(chǎng)和風(fēng)洞測(cè)試方法.Fremion等[4]對(duì)法國(guó)高速列車(TGV)轉(zhuǎn)向架和車輛連接部位的氣動(dòng)噪聲源特性進(jìn)行了現(xiàn)場(chǎng)測(cè)試,發(fā)現(xiàn)轉(zhuǎn)向架區(qū)域噪聲由多種非相干聲源組成,形成機(jī)理復(fù)雜.Noh等[5]對(duì)不同在線運(yùn)行速度的高速列車采用聲陣列與波束成形技術(shù)進(jìn)行測(cè)試,獲得了各主要?dú)鈩?dòng)噪聲源的貢獻(xiàn)量及頻率特性.隨著計(jì)算能力的提高,很多研究采用數(shù)值方法模擬分析了高速列車整車和各關(guān)鍵部位的氣動(dòng)噪聲特性.杜健等[6]對(duì)高速列車受電弓氣動(dòng)噪聲的計(jì)算結(jié)果表明,聲輻射具有指向性并存在主頻,并認(rèn)為受電弓頂部橫梁是誘發(fā)氣動(dòng)噪聲的主要因素.張亞?wèn)|等[7]建立了細(xì)化轉(zhuǎn)向架結(jié)構(gòu)的列車模型,數(shù)值計(jì)算了遠(yuǎn)場(chǎng)氣動(dòng)噪聲,并分析了轉(zhuǎn)向架艙內(nèi)各主要部件的氣動(dòng)力性能以及轉(zhuǎn)向架的速度場(chǎng)與氣動(dòng)噪聲源特性.Zhu等[8]考慮轉(zhuǎn)向架艙外的高速列車車輪部分直接受到來(lái)流沖擊,分析了單獨(dú)輪對(duì)氣動(dòng)噪聲的形成機(jī)理及輻射特性.由于高速列車結(jié)構(gòu)復(fù)雜,已有研究多采用非結(jié)構(gòu)化網(wǎng)格進(jìn)行流場(chǎng)計(jì)算,然后根據(jù)聲比擬方法計(jì)算遠(yuǎn)場(chǎng)氣動(dòng)噪聲,探討時(shí)域與頻域內(nèi)氣動(dòng)噪聲輻射特性及變化規(guī)律.聲比擬方法進(jìn)行氣動(dòng)噪聲預(yù)測(cè)時(shí),無(wú)法深入理解流體流動(dòng)發(fā)聲機(jī)理.為了進(jìn)一步分析轉(zhuǎn)向架區(qū)域各結(jié)構(gòu)部件間流體相互作用與氣動(dòng)噪聲產(chǎn)生機(jī)理,本文結(jié)合渦聲理論和 Ffowcs Williams-Hawkings(FW-H)聲比擬方法,分析高速列車轉(zhuǎn)向架艙結(jié)構(gòu)對(duì)轉(zhuǎn)向架區(qū)域流場(chǎng)與氣動(dòng)噪聲特性的影響,為研究該區(qū)域氣動(dòng)噪聲控制的有效措施建立理論基礎(chǔ).

    1 氣動(dòng)噪聲預(yù)測(cè)

    高速列車周圍流動(dòng)處于低馬赫數(shù)區(qū)域,幾何體結(jié)構(gòu)表面壓力脈動(dòng)產(chǎn)生的面偶極子聲源是流動(dòng)誘發(fā)噪聲的主要聲源,可以根據(jù)不可壓縮流體進(jìn)行計(jì)算[9].采用開源軟件OpenFOAM求解流體動(dòng)力學(xué)控制方程,空間和時(shí)間離散均采用二階精度格式.數(shù)值模擬采用基于Spalart-Allmaras湍流模型的延遲分離渦(DDES)模型.對(duì)于DDES模型,在固體表面附近的邊界層區(qū)域內(nèi)采用雷諾平均模型求解流場(chǎng),而在流動(dòng)分離區(qū)域內(nèi)應(yīng)用大渦模擬模型進(jìn)行計(jì)算流體力學(xué)(CFD)數(shù)值計(jì)算[10].

    根據(jù)聲比擬方法,近場(chǎng)流體計(jì)算結(jié)果將為遠(yuǎn)場(chǎng)聲學(xué)預(yù)測(cè)提供聲源數(shù)據(jù).假設(shè)流場(chǎng)參數(shù)為廣義函數(shù),F(xiàn)fowcs-Williams等[11]得出了廣義波傳播方程,如下所示:

    (1)

    式中:c0為周圍無(wú)擾動(dòng)時(shí)聲在空氣中的傳播速度;p′為遠(yuǎn)場(chǎng)聲壓;ρ為空氣密度,ρ0為靜止流場(chǎng)的空氣密度;un和vn為垂直于固體表面的流體和固體表面速度分量;ui為i方向上的流體表面速度分量;δ(f)為Dirac delta函數(shù);xi或xj為笛卡爾坐標(biāo)系的3個(gè)方向;Pij為可壓縮應(yīng)力張量;nj為對(duì)應(yīng)于單位法向量的分量;Tij為L(zhǎng)ighthill湍流應(yīng)力張量;H(f)為Heaviside函數(shù).Brentner等[12]對(duì)方程(1)進(jìn)行了數(shù)值求解.

    聲比擬方法無(wú)法深入理解流動(dòng)發(fā)聲機(jī)理,渦聲理論則揭示了流體自由湍流渦運(yùn)動(dòng)誘導(dǎo)噪聲的形成機(jī)理[13-14],對(duì)于等熵低速流動(dòng),渦聲方程可以表示為

    (2)

    式中:ω、v分別為渦矢量與速度矢量,(ω×v)為L(zhǎng)amb矢量.方程右端第1項(xiàng)表示聲源來(lái)自于渦線在速度場(chǎng)中的拉伸變形,為體偶極子聲源;方程右端第2項(xiàng)為流體動(dòng)能分布不均勻形成的聲源,通常為體四極子聲源.對(duì)于低馬赫數(shù)下流動(dòng),渦聲方程中由體偶極子和體四極子組成的聲源項(xiàng)在數(shù)學(xué)上近似于聲比擬方法中Lighthill四極子聲源項(xiàng).

    2 轉(zhuǎn)向架區(qū)域計(jì)算模型

    高速列車轉(zhuǎn)向架區(qū)域1∶10縮比模型如圖1所示,圖1a為單獨(dú)轉(zhuǎn)向架模型(工況1),圖1b為轉(zhuǎn)向架位于轉(zhuǎn)向架艙內(nèi)模型(工況2).一般認(rèn)為轉(zhuǎn)向架結(jié)構(gòu)位于轉(zhuǎn)向架艙外部分對(duì)轉(zhuǎn)向架區(qū)域流場(chǎng)影響較大[2],故對(duì)轉(zhuǎn)向架內(nèi)部結(jié)構(gòu)進(jìn)行了簡(jiǎn)化,忽略兩輪軸間各部件.由于轉(zhuǎn)向架區(qū)域結(jié)構(gòu)模型沿展向?qū)ΨQ,故車體縱向?qū)ΨQ面上流場(chǎng)作用的隨機(jī)性影響較小,為緩解對(duì)計(jì)算資源的要求,采用展向半車模型進(jìn)行數(shù)值模擬[15].計(jì)算區(qū)域采用全結(jié)構(gòu)化網(wǎng)格劃分,邊界層內(nèi)第1層網(wǎng)格至幾何體固體表面距離為1×10-5m,并沿著壁面法向量方向以1.1的增長(zhǎng)率向外拓展,使得y+(第1層網(wǎng)格質(zhì)心到壁面的量綱一距離)的最大值不超過(guò)1.各計(jì)算區(qū)域內(nèi)共生成1.46×107個(gè)六面體網(wǎng)格(工況1)和3.65×107個(gè)六面體網(wǎng)格(工況2).計(jì)算區(qū)域的邊界條件設(shè)置如下:來(lái)流設(shè)置為低湍流度均勻流(U∞=30 m·s-1),頂面、底面與兩側(cè)面均設(shè)為對(duì)稱面(相當(dāng)于剪切為零的滑移壁面),出口設(shè)為壓力出口,幾何體固體表面定義為靜態(tài)無(wú)滑移壁面.數(shù)值仿真時(shí),所取時(shí)間步長(zhǎng)確保Courant-Friedrichs-Lewy數(shù)在大部分計(jì)算區(qū)域內(nèi)小于1而在整個(gè)計(jì)算區(qū)域內(nèi)的最大值不超過(guò)2.

    a 工況1

    b 工況2圖1 轉(zhuǎn)向架區(qū)域縮比模型Fig.1 Scaled model of the bogie region

    3 瞬態(tài)流場(chǎng)分析

    為理解轉(zhuǎn)向架周圍流場(chǎng)特性及其對(duì)氣動(dòng)噪聲的影響,對(duì)上述2種工況下計(jì)算結(jié)果中轉(zhuǎn)向架物面附近流場(chǎng)的Q準(zhǔn)則等值面與渦量場(chǎng)進(jìn)行分析,普遍認(rèn)為湍流渦生成集中與發(fā)展劇烈部位會(huì)成為氣動(dòng)噪聲源.

    單獨(dú)轉(zhuǎn)向架工況流場(chǎng)中物面附近的湍流渦結(jié)構(gòu)如圖2所示.采用Q準(zhǔn)則顯示α=25的等值面(α=Q/(U∞/D)2,其中Q為速度梯度張量第二不變量,D為車輪直徑,U∞為來(lái)流速度).從圖2a可以發(fā)現(xiàn):流體在前輪對(duì)車輪前端邊緣分離后又黏附在車輪外側(cè)表面,形成沿車輪前緣呈月牙形狀的分離泡,隨后馬蹄渦在車輪表面上產(chǎn)生并沿下游傳播;同樣,流動(dòng)在轉(zhuǎn)向架構(gòu)架端部產(chǎn)生分離后,不同尺度的湍流渦沿著構(gòu)架表面向下游傳播,最后均合并入后輪對(duì)尾流中.由圖2b可見:流動(dòng)分離發(fā)生在前輪對(duì)和構(gòu)架前緣,前輪對(duì)車軸頂部生成了準(zhǔn)二維展向渦,隨著渦脫落沿下游發(fā)展成為流向渦;在輪對(duì)與構(gòu)架聯(lián)結(jié)區(qū)域,車軸尾流渦與車輪以及構(gòu)架表面產(chǎn)生的渦相互作用,形成了具有較強(qiáng)湍流度的湍流渦結(jié)構(gòu),使得轉(zhuǎn)向架上游輪對(duì)尾流由不同尺度與強(qiáng)度的湍流渦組成,這些湍流渦向下游傳播,與轉(zhuǎn)向架下游幾何體產(chǎn)生流動(dòng)沖擊作用,導(dǎo)致后輪對(duì)車軸尾流內(nèi)大尺度相干渦結(jié)構(gòu)明顯減少而發(fā)展成大量具有不同方向的小尺度渦結(jié)構(gòu).圖3為幾何體周圍的局部流線圖.從圖3a可以發(fā)現(xiàn),車輪兩側(cè)邊緣發(fā)生流動(dòng)分離后,剪切層離開車輪表面,轉(zhuǎn)捩后的湍流克服逆壓梯度的影響,發(fā)生流動(dòng)再附,在車輪側(cè)面形成分離泡.從圖3b 可以看出,流體在車軸頂部與底部分離點(diǎn)下游鄰近區(qū)域生成了沿輪軸軸向的準(zhǔn)二維展向渦,并逐漸發(fā)展成規(guī)則旋渦脫落后向下游傳播.

    a 側(cè)視圖

    b 仰視圖圖2 Q準(zhǔn)則顯示的轉(zhuǎn)向架物面附近渦結(jié)構(gòu)(工況1)Fig.2 Iso-surface of instantaneous normalized Q- criterion (case 1)

    轉(zhuǎn)向架位于轉(zhuǎn)向架艙內(nèi)的工況中轉(zhuǎn)向架物面附近的湍流渦結(jié)構(gòu)如圖4所示.采用Q準(zhǔn)則顯示α=25的等值面.從圖4a可以發(fā)現(xiàn),由于受轉(zhuǎn)向架艙前端分離的剪切層流動(dòng)的影響,前輪對(duì)和轉(zhuǎn)向架艙頂壁之間產(chǎn)生了較強(qiáng)的流體相互作用,形成大量不規(guī)則湍流渦并沿著轉(zhuǎn)向架艙頂壁向下游傳播.由圖4b可見,流體在轉(zhuǎn)向架艙端部分離后,與轉(zhuǎn)向架前輪對(duì)和構(gòu)架發(fā)生較強(qiáng)的相互作用,形成不同尺度的渦結(jié)構(gòu),由于該區(qū)域來(lái)流湍流度較高,與單獨(dú)轉(zhuǎn)向架工況相比,前輪對(duì)尾流內(nèi)大尺度相干渦結(jié)構(gòu)已不出現(xiàn).轉(zhuǎn)向架艙上游部位產(chǎn)生的渦向下游傳播,與轉(zhuǎn)向架后輪對(duì)和構(gòu)架后端進(jìn)行流動(dòng)作用后產(chǎn)生尾渦,進(jìn)而撞擊轉(zhuǎn)向架艙后壁面,形成較多高湍流度渦結(jié)構(gòu)并向車體下游傳播.

    a 前輪對(duì)車輪前端

    b 前輪對(duì)車軸后端圖3 局部流線(工況1)Fig.3 Local streamline (case 1)

    a 側(cè)視圖

    b 仰視圖圖4 Q準(zhǔn)則顯示的轉(zhuǎn)向架物面附近渦結(jié)構(gòu)(工況2)Fig.4 Iso-surface of instantaneous normalized Q- criterion (case 2)

    對(duì)于工況2,車軸與車輪中截面的渦量場(chǎng)輪廓(ωz)如圖6所示.可以發(fā)現(xiàn),一剪切層在轉(zhuǎn)向架艙前緣脫落并沿流向傳播,與前輪對(duì)車軸上端的分離渦相互作用后,迅速向上彎曲,撞擊在轉(zhuǎn)向架艙的頂壁上,導(dǎo)致該區(qū)域渦結(jié)構(gòu)復(fù)雜;另外,由于該剪切層的直接作用,前車軸上流動(dòng)提前分離,破壞了車軸后規(guī)則渦脫落的形成,使得前車軸尾流內(nèi)均形成了湍流度較高的不規(guī)則渦,所發(fā)展的渦向下游傳播并與后車軸周圍流體發(fā)生相互作用,最終后車軸尾渦與轉(zhuǎn)向架艙內(nèi)渦混合在一起,作用在轉(zhuǎn)向架艙后壁與后緣上,使得該區(qū)域流動(dòng)作用劇烈,渦生成集中.

    圖5 車軸截面瞬態(tài)渦量場(chǎng)輪廓(工況1)

    Fig.5Contourofinstantaneousspanwisevorticityfield(case1)

    圖6 車軸截面瞬態(tài)渦量場(chǎng)輪廓(工況2)

    Fig.6Contourofinstantaneousspanwisevorticityfield(case2)

    4 近場(chǎng)四極子聲源特性

    渦聲方程(2)右端兩聲源項(xiàng)中密度為常量,為便于數(shù)據(jù)處理,省略該密度項(xiàng),由此繪出了2種工況下轉(zhuǎn)向架水平面中心截面內(nèi)作為L(zhǎng)amb矢量散度的體偶極子聲源分布(見圖7)和流體動(dòng)能分布不均勻形成的體四極子聲源分布(見圖8),其中體偶極子和體四極子聲源強(qiáng)度在云圖顯示時(shí)采用了相同數(shù)量級(jí),以便于比較2類不同聲源的分布范圍.從圖7和圖8可見:前后車輪前緣、車軸尾流渦脫落以及構(gòu)架外側(cè)流體分離為轉(zhuǎn)向架周圍主要噪聲源;單獨(dú)轉(zhuǎn)向架結(jié)構(gòu)中(見圖7a和圖8a),上游車軸流動(dòng)分離產(chǎn)生的展向交替分離渦形成了較強(qiáng)的氣動(dòng)聲源,所產(chǎn)生的渦向下游傳播,與后輪對(duì)車軸進(jìn)行流動(dòng)作用后,在轉(zhuǎn)向架尾流中產(chǎn)生了氣動(dòng)噪聲源.當(dāng)轉(zhuǎn)向架位于轉(zhuǎn)向架艙時(shí)(見圖7b和圖8b),與單獨(dú)轉(zhuǎn)向架工況相比,轉(zhuǎn)向架結(jié)構(gòu)周圍形成的氣動(dòng)噪聲源產(chǎn)生了一定程度的削弱,但轉(zhuǎn)向架艙后壁部位由于劇烈的流動(dòng)沖擊作用,發(fā)展成為強(qiáng)氣動(dòng)噪聲源.為進(jìn)行聲源強(qiáng)度對(duì)比,圖9繪出了2種工況下轉(zhuǎn)向架前輪對(duì)車軸尾流中沿近壁流線方向的噪聲源分布.圖9表明,體偶極子聲源的強(qiáng)度遠(yuǎn)大于體四極子聲源的強(qiáng)度.因此,與體四極子聲源相比,體偶極子聲源中由流動(dòng)能量轉(zhuǎn)化為聲能的效率較高.另外,對(duì)于高速列車轉(zhuǎn)向架區(qū)域,各部件表面流動(dòng)分離產(chǎn)生較強(qiáng)的流體相互作用,使得靠近幾何體近壁區(qū)流場(chǎng)內(nèi)的渦運(yùn)動(dòng)劇烈,所形成的體偶極子聲源發(fā)展為近場(chǎng)體四極子聲源的主要組成部分.

    a 工況1

    b 工況2圖7 轉(zhuǎn)向架水平面中心截面體偶極子聲源分布

    Fig.7Contourofvolumedipolesourcedistributiononwheelsethorizontalmid-plane

    a 工況1

    b 工況2圖8 轉(zhuǎn)向架水平面中心截面體四極子聲源分布

    Fig.8Contourofvolumequadrpolesourcedistributiononwheelsethorizontalmid-plane

    5 遠(yuǎn)場(chǎng)氣動(dòng)噪聲特性

    當(dāng)瞬態(tài)流場(chǎng)達(dá)到統(tǒng)計(jì)意義上穩(wěn)定后,運(yùn)用FW-H聲比擬方法,以近場(chǎng)流場(chǎng)數(shù)據(jù)作為源項(xiàng)對(duì)遠(yuǎn)場(chǎng)輻射噪聲進(jìn)行預(yù)測(cè),數(shù)值仿真結(jié)果通過(guò)低噪聲風(fēng)洞進(jìn)行試驗(yàn)驗(yàn)證.采用Welch方法對(duì)計(jì)算所得的遠(yuǎn)場(chǎng)噪聲時(shí)間序列進(jìn)行功率譜密度分析[15],然后進(jìn)行總聲壓級(jí)的計(jì)算,所用頻率分辨率為6 Hz.

    a 工況1

    b 工況2圖9 轉(zhuǎn)向架前輪對(duì)車軸尾流中沿流線方向四極子聲源分布

    Fig.9Quadrpolenoisesourcedistributionalongthestreamlineinthefrontaxlewakeofbogie

    5.1 聲學(xué)風(fēng)洞試驗(yàn)驗(yàn)證

    在英國(guó)南安普頓大學(xué)聲與振動(dòng)研究所(ISVR)的低噪聲風(fēng)洞內(nèi)進(jìn)行單獨(dú)轉(zhuǎn)向架模型氣動(dòng)噪聲數(shù)值模擬的試驗(yàn)驗(yàn)證[15-16].轉(zhuǎn)向架模型被固定在噴口側(cè)向支承板上(見圖10),遠(yuǎn)場(chǎng)測(cè)點(diǎn)位于輪對(duì)頂部;來(lái)流速度為30 m·s-1,湍流度低于0.3%.由于100 Hz以下頻域內(nèi)噴口噴流產(chǎn)生較大的背景噪聲而實(shí)際轉(zhuǎn)向架結(jié)構(gòu)產(chǎn)生的氣動(dòng)噪聲幅值很小,故該頻域內(nèi)噪聲輻射數(shù)據(jù)將被忽略.從圖11可以看出,轉(zhuǎn)向架模型遠(yuǎn)場(chǎng)測(cè)點(diǎn)氣動(dòng)噪聲的數(shù)值模擬與聲學(xué)風(fēng)洞測(cè)試結(jié)果吻合得很好.對(duì)于主頻所對(duì)應(yīng)的幅值,計(jì)算值高于測(cè)試值,這是由于數(shù)值模擬時(shí)采用對(duì)稱面的邊界條件,而在試驗(yàn)中使用支承板固定輪對(duì),削弱了車軸尾流中規(guī)則的相干渦旋脫落的產(chǎn)生與其所誘發(fā)的氣動(dòng)噪聲;同時(shí),數(shù)值計(jì)算所得主頻值(324 Hz)略高于測(cè)試值(314 Hz),誤差約為3%.究其原因,主要是由于轉(zhuǎn)向架結(jié)構(gòu)沿縱向較長(zhǎng),風(fēng)洞測(cè)試時(shí)可能沒(méi)有完全處于噴口自由射流的勢(shì)流核心區(qū)內(nèi).

    圖10 低噪聲風(fēng)洞內(nèi)轉(zhuǎn)向架測(cè)試裝置Fig.10 Experimental setup in anechoic chamber

    圖11 遠(yuǎn)場(chǎng)氣動(dòng)噪聲預(yù)測(cè)的試驗(yàn)驗(yàn)證Fig.11 Experimental verification of far-field noise prediction

    5.2 轉(zhuǎn)向架氣動(dòng)噪聲頻譜分析

    對(duì)于單獨(dú)轉(zhuǎn)向架和轉(zhuǎn)向架位于轉(zhuǎn)向架艙內(nèi)2種工況,采用流場(chǎng)穩(wěn)定后繼續(xù)運(yùn)行0.52 s進(jìn)行物面壓力脈動(dòng)時(shí)間序列信號(hào)的采集,以進(jìn)行遠(yuǎn)場(chǎng)噪聲預(yù)測(cè).近場(chǎng)噪聲源為轉(zhuǎn)向架半模型的固體表面,遠(yuǎn)場(chǎng)測(cè)點(diǎn)位于以轉(zhuǎn)向架中心為原點(diǎn)的展向平面內(nèi).圖12a比較了2種工況下轉(zhuǎn)向架頂部接收點(diǎn)輻射噪聲的功率譜密度.對(duì)于工況2,假設(shè)轉(zhuǎn)向架固面產(chǎn)生的氣動(dòng)噪聲直接輻射至遠(yuǎn)場(chǎng)聲接收點(diǎn).結(jié)果表明,與單獨(dú)轉(zhuǎn)向架相比,位于轉(zhuǎn)向架艙內(nèi)轉(zhuǎn)向架的遠(yuǎn)場(chǎng)輻射噪聲級(jí)顯著降低,在2 kHz以下頻域范圍內(nèi)從72.7 dB降至55.2 dB,下降了17.5 dB.這主要是由于當(dāng)轉(zhuǎn)向架位于轉(zhuǎn)向架艙內(nèi)時(shí),來(lái)流速度降低的同時(shí)轉(zhuǎn)向架周圍不再產(chǎn)生大尺度的規(guī)則渦脫落.圖12b比較了2種工況下轉(zhuǎn)向架側(cè)向接收點(diǎn)輻射噪聲的功率譜密度.在2 kHz以下頻域內(nèi)轉(zhuǎn)向架艙內(nèi)轉(zhuǎn)向架的遠(yuǎn)場(chǎng)輻射噪聲較單獨(dú)轉(zhuǎn)向架下降了約3 dB.與單獨(dú)轉(zhuǎn)向架相比,轉(zhuǎn)向架位于轉(zhuǎn)向架艙內(nèi)后,不同方向的氣動(dòng)噪聲級(jí)下降幅度不同,顯然受到了結(jié)構(gòu)聲輻射指向性的影響.

    a 頂部接收點(diǎn)

    b 側(cè)向接收點(diǎn)圖12 遠(yuǎn)場(chǎng)測(cè)點(diǎn)輻射噪聲頻譜比較Fig.12 Comparisons of noise spectra of far-field receivers

    5.3 轉(zhuǎn)向架聲輻射指向性分析

    根據(jù)展向半轉(zhuǎn)向架模型所得流場(chǎng)數(shù)據(jù),整個(gè)轉(zhuǎn)向架的聲壓級(jí)由沿著展向?qū)ΨQ面對(duì)稱布置的2個(gè)接收點(diǎn)聲壓級(jí)(Lp1和Lp2)計(jì)算得出.假定兩對(duì)稱半轉(zhuǎn)向架聲源互不相關(guān),則Lp=10log(10Lp1/10+10Lp2/10).進(jìn)行聲輻射指向性計(jì)算時(shí),遠(yuǎn)場(chǎng)1 946個(gè)聲接收點(diǎn)均勻分布在以轉(zhuǎn)向架中心為球心、半徑為2.5 m的球面上,每測(cè)點(diǎn)根據(jù)聲壓時(shí)間序列計(jì)算出總聲壓級(jí)值(低于5 kHz頻域內(nèi)).對(duì)于轉(zhuǎn)向架艙內(nèi)轉(zhuǎn)向架工況,假設(shè)氣動(dòng)噪聲由轉(zhuǎn)向架表面向遠(yuǎn)場(chǎng)輻射,不考慮車體對(duì)聲波的反射與遮擋作用,以便于與單獨(dú)轉(zhuǎn)向架工況進(jìn)行比較.圖13a顯示了單獨(dú)轉(zhuǎn)向架的聲輻射指向性為垂向偶極子形狀,這主要是由于單獨(dú)轉(zhuǎn)向架的輪對(duì)車軸表面上較大的壓力脈動(dòng)產(chǎn)生了交替分離渦,從而形成了較強(qiáng)的面偶極子聲源.圖13b表明,轉(zhuǎn)向架艙內(nèi)轉(zhuǎn)向架的聲輻射指向性為橫向偶極子形狀,噪聲級(jí)幅值較單獨(dú)轉(zhuǎn)向架工況下降約12 dB.究其原因,轉(zhuǎn)向架艙的存在,改變了轉(zhuǎn)向架周圍的流動(dòng)特性,抑制了輪對(duì)車軸表面規(guī)則交替分離渦的形成,削弱了偶極子氣動(dòng)噪聲源產(chǎn)生的強(qiáng)度.與單獨(dú)轉(zhuǎn)向架相比,轉(zhuǎn)向架艙改變了轉(zhuǎn)向架氣動(dòng)噪聲的聲輻射指向性,同時(shí)降低了其噪聲級(jí);轉(zhuǎn)向架艙后緣壁面由于受到轉(zhuǎn)向架尾流的流動(dòng)沖擊作用,會(huì)產(chǎn)生較大氣動(dòng)噪聲.

    a 工況1b 工況2

    圖13轉(zhuǎn)向架氣動(dòng)噪聲聲輻射指向性(中心位于50dB)

    Fig.13Noisedirectivityforbogie(startingat50dB)

    6 結(jié)語(yǔ)

    本文根據(jù)渦聲理論和聲比擬方法,計(jì)算分析了轉(zhuǎn)向架艙對(duì)轉(zhuǎn)向架結(jié)構(gòu)流動(dòng)與氣動(dòng)噪聲性能的影響.結(jié)果表明:當(dāng)轉(zhuǎn)向架位于轉(zhuǎn)向架艙內(nèi)時(shí),流體在轉(zhuǎn)向架艙前緣端部分離并與轉(zhuǎn)向架前端結(jié)構(gòu)發(fā)生較強(qiáng)的相互作用,與單獨(dú)轉(zhuǎn)向架工況相比,轉(zhuǎn)向架前輪對(duì)尾流內(nèi)大尺度相干渦結(jié)構(gòu)已不出現(xiàn);2種工況下,靠近幾何體近壁區(qū)流場(chǎng)內(nèi),渦運(yùn)動(dòng)劇烈,所形成的體偶極子聲源為近場(chǎng)四極子噪聲的主要聲源.對(duì)于遠(yuǎn)場(chǎng)輻射噪聲,轉(zhuǎn)向架各部件表面壓力脈動(dòng)形成的面偶極子聲源成為主要噪聲源;與單獨(dú)轉(zhuǎn)向架結(jié)構(gòu)相比,轉(zhuǎn)向架艙改變了轉(zhuǎn)向架周圍流體的流動(dòng)特性與轉(zhuǎn)向架氣動(dòng)噪聲的聲輻射指向性,削弱了轉(zhuǎn)向架結(jié)構(gòu)表面氣動(dòng)噪聲產(chǎn)生的強(qiáng)度,而車體轉(zhuǎn)向架艙后緣區(qū)域由于受到流動(dòng)沖擊作用,成為強(qiáng)氣動(dòng)噪聲源.對(duì)于實(shí)尺轉(zhuǎn)向架,輪軸間布置有牽引與制動(dòng)系統(tǒng)的相關(guān)結(jié)構(gòu)部件,將會(huì)影響單純輪軸尾流內(nèi)規(guī)則渦的形成與發(fā)展,從而影響氣動(dòng)噪聲的產(chǎn)生.本文基于簡(jiǎn)化比例模型的研究,有助于理解實(shí)際運(yùn)行車輛轉(zhuǎn)向架艙對(duì)于轉(zhuǎn)向架區(qū)域復(fù)雜流動(dòng)流場(chǎng)與氣動(dòng)噪聲特性的影響,為高速列車轉(zhuǎn)向架區(qū)域結(jié)構(gòu)優(yōu)化與噪聲控制提供理論依據(jù).

    猜你喜歡
    偶極子車軸聲源
    車軸軸徑微差測(cè)量?jī)x測(cè)量不確定度評(píng)定
    虛擬聲源定位的等效源近場(chǎng)聲全息算法
    高品質(zhì)生產(chǎn)
    基于DDS的正交偶極子聲波測(cè)井儀快檢裝置研究
    化工管理(2021年7期)2021-05-13 00:46:12
    弧形寬帶印刷偶極子5G天線的設(shè)計(jì)
    基于GCC-nearest時(shí)延估計(jì)的室內(nèi)聲源定位
    電子制作(2019年23期)2019-02-23 13:21:12
    運(yùn)用內(nèi)積相關(guān)性結(jié)合迭代相減識(shí)別兩點(diǎn)聲源
    高速列車空心車軸國(guó)產(chǎn)化的選材和試制
    力-聲互易在水下聲源強(qiáng)度測(cè)量中的應(yīng)用
    25CrMo鋼高速車軸坯的試制
    a级毛片a级免费在线| 超碰av人人做人人爽久久| 久久久久久九九精品二区国产| 三级男女做爰猛烈吃奶摸视频| 国产午夜精品论理片| 特级一级黄色大片| 人人妻人人看人人澡| 99精品在免费线老司机午夜| 免费人成在线观看视频色| 精品一区二区三区av网在线观看| 99热6这里只有精品| 日本撒尿小便嘘嘘汇集6| bbb黄色大片| 久久久久性生活片| 一级黄片播放器| 99热6这里只有精品| 91麻豆av在线| 精品一区二区三区av网在线观看| 99热精品在线国产| 亚洲欧美日韩高清在线视频| 亚洲不卡免费看| 成人亚洲精品av一区二区| 欧美国产日韩亚洲一区| 亚洲av不卡在线观看| 91精品国产九色| 国产精品永久免费网站| 欧美国产日韩亚洲一区| 亚洲av成人精品一区久久| 人妻丰满熟妇av一区二区三区| 欧美一级a爱片免费观看看| 精品一区二区三区av网在线观看| 大又大粗又爽又黄少妇毛片口| 久久久久免费精品人妻一区二区| 伊人久久精品亚洲午夜| 亚洲精品一卡2卡三卡4卡5卡| 成人高潮视频无遮挡免费网站| 精品不卡国产一区二区三区| 成人特级黄色片久久久久久久| 久久香蕉精品热| 人妻夜夜爽99麻豆av| 日韩一本色道免费dvd| 五月伊人婷婷丁香| 俄罗斯特黄特色一大片| 国产精品自产拍在线观看55亚洲| 18禁在线播放成人免费| 91在线精品国自产拍蜜月| 一区二区三区免费毛片| 午夜老司机福利剧场| 人妻少妇偷人精品九色| 成人二区视频| 尤物成人国产欧美一区二区三区| 级片在线观看| 国产亚洲精品av在线| 中文资源天堂在线| 免费看光身美女| 免费av不卡在线播放| 又紧又爽又黄一区二区| 两个人的视频大全免费| 欧美区成人在线视频| 三级国产精品欧美在线观看| 麻豆久久精品国产亚洲av| 国产三级在线视频| 热99在线观看视频| 亚洲熟妇中文字幕五十中出| 丰满乱子伦码专区| 搡老妇女老女人老熟妇| 久久久久国内视频| 给我免费播放毛片高清在线观看| 亚洲成人久久爱视频| 一级黄片播放器| 一级av片app| 午夜福利成人在线免费观看| 亚洲第一电影网av| 国产探花在线观看一区二区| 国产精品99久久久久久久久| 日日摸夜夜添夜夜添av毛片 | 欧美日韩精品成人综合77777| 黄片wwwwww| 91麻豆精品激情在线观看国产| 88av欧美| 我要看日韩黄色一级片| 一夜夜www| 禁无遮挡网站| 日本撒尿小便嘘嘘汇集6| 亚洲人成网站在线播| 日韩大尺度精品在线看网址| 国产精品亚洲美女久久久| 伦理电影大哥的女人| 999久久久精品免费观看国产| 国产精品久久久久久亚洲av鲁大| 成人av一区二区三区在线看| 免费av不卡在线播放| 看十八女毛片水多多多| 在线观看免费视频日本深夜| 成人av一区二区三区在线看| 日日啪夜夜撸| 国产精品av视频在线免费观看| 久久中文看片网| 91久久精品国产一区二区三区| 99riav亚洲国产免费| 狂野欧美白嫩少妇大欣赏| 一卡2卡三卡四卡精品乱码亚洲| 最近最新中文字幕大全电影3| 欧美最黄视频在线播放免费| 国产激情偷乱视频一区二区| 国产乱人伦免费视频| 又粗又爽又猛毛片免费看| 91av网一区二区| 亚洲av免费在线观看| 国产精品女同一区二区软件 | 久久精品国产亚洲av涩爱 | 国产人妻一区二区三区在| 桃色一区二区三区在线观看| 久久精品国产亚洲av天美| a级一级毛片免费在线观看| 亚洲欧美日韩东京热| 不卡一级毛片| 亚洲成人精品中文字幕电影| www.色视频.com| 琪琪午夜伦伦电影理论片6080| 国内毛片毛片毛片毛片毛片| 乱人视频在线观看| 亚洲精品久久国产高清桃花| 99热6这里只有精品| 99久久成人亚洲精品观看| 欧美潮喷喷水| 国产白丝娇喘喷水9色精品| 此物有八面人人有两片| 亚洲美女视频黄频| 最近视频中文字幕2019在线8| 69av精品久久久久久| 久99久视频精品免费| 欧美高清性xxxxhd video| 国产探花在线观看一区二区| .国产精品久久| 日本五十路高清| 色吧在线观看| 国产伦人伦偷精品视频| 日韩 亚洲 欧美在线| 日韩欧美国产在线观看| 国内精品一区二区在线观看| 久久精品人妻少妇| 久久久久久久久大av| 亚洲狠狠婷婷综合久久图片| 久久午夜福利片| 国产中年淑女户外野战色| 级片在线观看| 伊人久久精品亚洲午夜| 国内精品美女久久久久久| 国产蜜桃级精品一区二区三区| 91久久精品国产一区二区三区| 亚洲va在线va天堂va国产| 欧美中文日本在线观看视频| 国产精品电影一区二区三区| 在线观看免费视频日本深夜| 国产男靠女视频免费网站| 久久精品综合一区二区三区| 成年人黄色毛片网站| 亚洲国产精品久久男人天堂| 久久久久久久久中文| 久久香蕉精品热| 九九久久精品国产亚洲av麻豆| 欧美性猛交黑人性爽| 美女xxoo啪啪120秒动态图| 日本一二三区视频观看| 婷婷丁香在线五月| 国产av在哪里看| 日本色播在线视频| 精品久久久久久久久av| 69人妻影院| 两人在一起打扑克的视频| 国产 一区 欧美 日韩| 欧美日韩瑟瑟在线播放| 男女啪啪激烈高潮av片| 日韩欧美一区二区三区在线观看| 成年女人永久免费观看视频| 欧美在线一区亚洲| 99热6这里只有精品| 成人美女网站在线观看视频| 免费看光身美女| 亚洲真实伦在线观看| 久久久精品欧美日韩精品| 春色校园在线视频观看| 中出人妻视频一区二区| 赤兔流量卡办理| 999久久久精品免费观看国产| 久久婷婷人人爽人人干人人爱| 欧美日韩乱码在线| 热99re8久久精品国产| 色av中文字幕| 日日摸夜夜添夜夜添av毛片 | 日韩欧美精品免费久久| 国产v大片淫在线免费观看| 中文字幕高清在线视频| 日韩中文字幕欧美一区二区| 欧美bdsm另类| 亚洲狠狠婷婷综合久久图片| 成人av在线播放网站| 九色国产91popny在线| 国产精品亚洲美女久久久| 亚洲专区中文字幕在线| 亚洲精品色激情综合| 69人妻影院| 美女高潮的动态| 人人妻人人看人人澡| 不卡一级毛片| 久久久久久久久大av| 日韩av在线大香蕉| 麻豆精品久久久久久蜜桃| 波多野结衣巨乳人妻| 亚洲第一区二区三区不卡| 久久精品国产清高在天天线| 女人十人毛片免费观看3o分钟| 桃色一区二区三区在线观看| 欧美最黄视频在线播放免费| 欧美三级亚洲精品| 久久久久久久久大av| 色精品久久人妻99蜜桃| 一级av片app| av在线亚洲专区| 日韩欧美国产一区二区入口| 国产欧美日韩精品亚洲av| 最后的刺客免费高清国语| 久久久久久久久久久丰满 | 精品人妻视频免费看| 亚洲最大成人av| bbb黄色大片| 窝窝影院91人妻| 22中文网久久字幕| АⅤ资源中文在线天堂| 69人妻影院| 亚洲av不卡在线观看| 亚洲不卡免费看| 精品福利观看| 久久精品夜夜夜夜夜久久蜜豆| 69人妻影院| 亚洲美女搞黄在线观看 | 三级男女做爰猛烈吃奶摸视频| 在线观看一区二区三区| 国产伦人伦偷精品视频| 国产亚洲欧美98| 日韩在线高清观看一区二区三区 | 欧美一区二区国产精品久久精品| 男人的好看免费观看在线视频| 老师上课跳d突然被开到最大视频| 精品不卡国产一区二区三区| 日韩精品有码人妻一区| 亚洲欧美日韩卡通动漫| 蜜桃久久精品国产亚洲av| 国产精品98久久久久久宅男小说| 久久精品夜夜夜夜夜久久蜜豆| 国产精品久久久久久av不卡| 国产精品永久免费网站| 白带黄色成豆腐渣| 3wmmmm亚洲av在线观看| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩亚洲国产一区二区在线观看| 国产毛片a区久久久久| 国产真实伦视频高清在线观看 | 国产伦精品一区二区三区四那| 欧美日韩综合久久久久久 | 天美传媒精品一区二区| 极品教师在线视频| 伊人久久精品亚洲午夜| 日韩精品有码人妻一区| 午夜福利在线观看吧| 尤物成人国产欧美一区二区三区| 在线观看66精品国产| 色综合亚洲欧美另类图片| 久久午夜福利片| 久久国产精品人妻蜜桃| 国产女主播在线喷水免费视频网站 | 一级毛片久久久久久久久女| 一进一出抽搐动态| 少妇的逼水好多| 国产三级在线视频| 国内毛片毛片毛片毛片毛片| 成年女人毛片免费观看观看9| 亚洲国产精品合色在线| 日本五十路高清| 性插视频无遮挡在线免费观看| 九九爱精品视频在线观看| 内射极品少妇av片p| 少妇丰满av| 欧美人与善性xxx| 深夜a级毛片| 久久6这里有精品| 日本三级黄在线观看| 在线观看66精品国产| 久久精品国产清高在天天线| 午夜免费激情av| 黄色欧美视频在线观看| 一进一出抽搐gif免费好疼| АⅤ资源中文在线天堂| 欧美xxxx性猛交bbbb| 欧美+日韩+精品| 一区福利在线观看| 精品久久久久久久末码| 老女人水多毛片| 免费在线观看影片大全网站| 丰满乱子伦码专区| 国产亚洲欧美98| h日本视频在线播放| 国产精品一及| 在现免费观看毛片| 成人一区二区视频在线观看| 国产爱豆传媒在线观看| 免费大片18禁| 极品教师在线视频| 老熟妇乱子伦视频在线观看| 99热只有精品国产| 18禁黄网站禁片午夜丰满| a级一级毛片免费在线观看| 亚洲美女黄片视频| 亚洲欧美日韩卡通动漫| 日韩欧美国产在线观看| 麻豆久久精品国产亚洲av| 美女xxoo啪啪120秒动态图| 蜜桃久久精品国产亚洲av| 999久久久精品免费观看国产| 淫秽高清视频在线观看| 99久国产av精品| 在线观看66精品国产| 床上黄色一级片| 夜夜夜夜夜久久久久| 欧美激情久久久久久爽电影| www.色视频.com| 别揉我奶头~嗯~啊~动态视频| 午夜日韩欧美国产| 搡老岳熟女国产| 99riav亚洲国产免费| 成人精品一区二区免费| 亚洲自拍偷在线| 69av精品久久久久久| 国产欧美日韩精品一区二区| 国产男靠女视频免费网站| 男女视频在线观看网站免费| 啦啦啦观看免费观看视频高清| 国产精品无大码| 美女cb高潮喷水在线观看| 国产激情偷乱视频一区二区| 婷婷色综合大香蕉| 麻豆国产av国片精品| 亚洲av二区三区四区| 国产精品一区二区三区四区久久| 麻豆成人午夜福利视频| 免费观看精品视频网站| 久久国产乱子免费精品| 亚洲国产精品成人综合色| 男女那种视频在线观看| 日本与韩国留学比较| 毛片一级片免费看久久久久 | 日本黄大片高清| or卡值多少钱| 直男gayav资源| 大型黄色视频在线免费观看| 国产高清有码在线观看视频| 国产综合懂色| 大型黄色视频在线免费观看| 两个人视频免费观看高清| 久久国产精品人妻蜜桃| 欧美又色又爽又黄视频| 69av精品久久久久久| 国产精品电影一区二区三区| 久久国产乱子免费精品| 高清毛片免费观看视频网站| 中文资源天堂在线| 免费av毛片视频| 在线看三级毛片| 看十八女毛片水多多多| 国产成年人精品一区二区| xxxwww97欧美| 日本成人三级电影网站| 丰满乱子伦码专区| 日本在线视频免费播放| 又粗又爽又猛毛片免费看| 国产精品久久视频播放| 在线国产一区二区在线| 久久99热这里只有精品18| 亚洲人与动物交配视频| 欧美不卡视频在线免费观看| 免费高清视频大片| 久久久国产成人精品二区| 不卡一级毛片| 天美传媒精品一区二区| 亚洲第一电影网av| 97碰自拍视频| 最后的刺客免费高清国语| 伦精品一区二区三区| 高清日韩中文字幕在线| 好男人在线观看高清免费视频| 亚洲性久久影院| 色吧在线观看| 成年免费大片在线观看| 亚洲成人免费电影在线观看| 少妇裸体淫交视频免费看高清| 亚洲av在线观看美女高潮| 午夜福利在线在线| 国产精品久久久久成人av| 久久av网站| 亚洲av.av天堂| 亚洲精品成人av观看孕妇| 一级毛片黄色毛片免费观看视频| 综合色丁香网| 精品久久久噜噜| 在线观看一区二区三区| 国产老妇伦熟女老妇高清| 久久久久国产网址| 国产精品欧美亚洲77777| 高清在线视频一区二区三区| 亚洲一区二区三区欧美精品| 一区二区三区免费毛片| 久热这里只有精品99| 一级毛片 在线播放| 男人狂女人下面高潮的视频| 人妻夜夜爽99麻豆av| 久久99热这里只频精品6学生| 男女无遮挡免费网站观看| 观看av在线不卡| 乱系列少妇在线播放| 性高湖久久久久久久久免费观看| 国产男女内射视频| 亚洲不卡免费看| 久久久亚洲精品成人影院| 精品久久久久久久末码| 亚洲精华国产精华液的使用体验| 爱豆传媒免费全集在线观看| 亚洲欧美一区二区三区国产| 久久鲁丝午夜福利片| 国产亚洲欧美精品永久| 国产一级毛片在线| 永久免费av网站大全| 美女福利国产在线 | 久久人人爽av亚洲精品天堂 | 亚洲国产av新网站| 性色av一级| 边亲边吃奶的免费视频| 亚洲欧美日韩卡通动漫| 久久毛片免费看一区二区三区| 日韩大片免费观看网站| 中文字幕制服av| 免费人妻精品一区二区三区视频| 日本欧美国产在线视频| 看十八女毛片水多多多| 久久久久人妻精品一区果冻| 国产人妻一区二区三区在| 美女xxoo啪啪120秒动态图| 丝袜喷水一区| 精品久久久久久久久亚洲| 欧美xxxx性猛交bbbb| 你懂的网址亚洲精品在线观看| 国产视频首页在线观看| 国产大屁股一区二区在线视频| 欧美日韩在线观看h| 最近的中文字幕免费完整| 日本黄色日本黄色录像| 菩萨蛮人人尽说江南好唐韦庄| 精品久久久久久久末码| 建设人人有责人人尽责人人享有的 | 寂寞人妻少妇视频99o| 久久久久久久久久人人人人人人| 亚洲av二区三区四区| 久久99蜜桃精品久久| 春色校园在线视频观看| 国产成人精品福利久久| 国产男女内射视频| 国产亚洲av片在线观看秒播厂| 欧美bdsm另类| 亚洲精品中文字幕在线视频 | 十八禁网站网址无遮挡 | 免费观看无遮挡的男女| 久久av网站| 丰满乱子伦码专区| 精品一区二区三卡| 亚洲av.av天堂| 亚洲美女搞黄在线观看| 亚洲精品日韩av片在线观看| 色视频www国产| 日韩av不卡免费在线播放| 男人爽女人下面视频在线观看| 99视频精品全部免费 在线| 一二三四中文在线观看免费高清| 亚洲最大成人中文| 少妇裸体淫交视频免费看高清| 97超碰精品成人国产| 97在线视频观看| av专区在线播放| av又黄又爽大尺度在线免费看| av在线老鸭窝| 亚洲国产av新网站| 亚洲图色成人| 黄片wwwwww| 熟妇人妻不卡中文字幕| 黄片wwwwww| 欧美 日韩 精品 国产| av播播在线观看一区| 欧美xxxx性猛交bbbb| 国产高潮美女av| 在线免费十八禁| 亚洲av日韩在线播放| 一区二区三区乱码不卡18| 国产男女超爽视频在线观看| 日日摸夜夜添夜夜爱| 免费大片黄手机在线观看| 高清不卡的av网站| 日韩,欧美,国产一区二区三区| 成人18禁高潮啪啪吃奶动态图 | 十分钟在线观看高清视频www | 国产人妻一区二区三区在| av在线蜜桃| 女性被躁到高潮视频| 插逼视频在线观看| 午夜老司机福利剧场| 黑丝袜美女国产一区| 国产精品久久久久久av不卡| 九草在线视频观看| 亚洲欧美日韩卡通动漫| 人妻制服诱惑在线中文字幕| 欧美日韩视频高清一区二区三区二| 97超碰精品成人国产| 国产亚洲最大av| 国产av一区二区精品久久 | 国产黄色视频一区二区在线观看| 国产免费一级a男人的天堂| 在线 av 中文字幕| 我的老师免费观看完整版| 在线 av 中文字幕| 交换朋友夫妻互换小说| 成年人午夜在线观看视频| 国产熟女欧美一区二区| 亚洲国产最新在线播放| 大陆偷拍与自拍| 国语对白做爰xxxⅹ性视频网站| 一个人看的www免费观看视频| 久久精品久久精品一区二区三区| 精品少妇黑人巨大在线播放| 丰满乱子伦码专区| 亚洲欧美成人精品一区二区| 欧美少妇被猛烈插入视频| 久久久精品免费免费高清| 777米奇影视久久| 精品视频人人做人人爽| 亚洲自偷自拍三级| 国产日韩欧美在线精品| 午夜激情久久久久久久| 一级毛片aaaaaa免费看小| 亚洲精品国产成人久久av| 久久毛片免费看一区二区三区| 激情五月婷婷亚洲| 国产毛片在线视频| 亚洲色图综合在线观看| 女人久久www免费人成看片| av卡一久久| 国产有黄有色有爽视频| 99久久精品国产国产毛片| 国产成人精品久久久久久| 在线免费十八禁| 欧美精品一区二区大全| 丝袜脚勾引网站| 天天躁夜夜躁狠狠久久av| 国产乱来视频区| 国产精品无大码| 女的被弄到高潮叫床怎么办| 99热国产这里只有精品6| 久久精品夜色国产| 蜜臀久久99精品久久宅男| 国产视频内射| 亚洲欧美日韩无卡精品| 国产精品伦人一区二区| 国产精品无大码| 国产久久久一区二区三区| 成人一区二区视频在线观看| 高清av免费在线| 欧美日本视频| 国产av国产精品国产| 欧美区成人在线视频| 女性被躁到高潮视频| freevideosex欧美| 中文字幕免费在线视频6| 99热网站在线观看| 欧美97在线视频| 国产精品一区二区在线不卡| 国产成人精品福利久久| 一本久久精品| 男女啪啪激烈高潮av片| 嘟嘟电影网在线观看| 色婷婷av一区二区三区视频| 成年女人在线观看亚洲视频| 中国国产av一级| 成人免费观看视频高清| 热99国产精品久久久久久7| 午夜福利在线在线| 精品99又大又爽又粗少妇毛片| 亚洲一级一片aⅴ在线观看| 久久婷婷青草| 色视频在线一区二区三区| 国产精品一区二区在线观看99| 亚洲av免费高清在线观看| 亚洲精品视频女| 国产v大片淫在线免费观看| 精品久久久久久电影网| 国产免费又黄又爽又色| 一本色道久久久久久精品综合| 午夜免费观看性视频| 我要看日韩黄色一级片| 十八禁网站网址无遮挡 | 我要看日韩黄色一级片| 高清视频免费观看一区二区| 91精品国产国语对白视频| av.在线天堂| 亚洲精品成人av观看孕妇| 国产精品无大码| 高清在线视频一区二区三区| 国产精品三级大全| 精品少妇黑人巨大在线播放|