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

    高速列車動車轉(zhuǎn)向架氣動噪聲數(shù)值分析

    2016-10-21 01:13:38西南交通大學(xué)牽引動力國家重點實驗室四川成都610031
    西南交通大學(xué)學(xué)報 2016年5期
    關(guān)鍵詞:遠場聲壓級動車

    (西南交通大學(xué)牽引動力國家重點實驗室,四川成都610031)

    (西南交通大學(xué)牽引動力國家重點實驗室,四川成都610031)

    為研究高速列車動車轉(zhuǎn)向架氣動噪聲特性,建立了動車轉(zhuǎn)向架空氣動力學(xué)模型,采用定常RNG k-ε湍流模型與寬頻帶噪聲源模型對其氣動噪聲聲源進行初步探討,并結(jié)合非定常LES大渦模擬與Lighthill聲學(xué)比擬理論進行了遠場氣動噪聲分析.研究結(jié)果表明:動車轉(zhuǎn)向架氣動噪聲源為輪對、構(gòu)架、牽引電機1、枕梁、垂向減振器、抗側(cè)滾扭桿等結(jié)構(gòu)的迎風側(cè)凸起部位,且構(gòu)架對動車轉(zhuǎn)向架遠場氣動噪聲的貢獻最大,其次為輪對和抗側(cè)滾扭桿,然后為垂向減振器和枕梁,牽引電機1、牽引電機2、空氣彈簧和橫向減振器對遠場氣動噪聲的貢獻較小.動車轉(zhuǎn)向架遠場氣動噪聲是寬頻噪聲,具有衰減特性、幅值特性和氣動噪聲指向性.在低頻部分能量較大,中心頻率為25、50 Hz,且分布規(guī)律不隨運行速度的改變而變化.

    高速列車;動車轉(zhuǎn)向架;氣動噪聲;大渦模擬;Lighthill聲學(xué)比擬理論;噪聲貢獻量

    為滿足交通運輸體系安全、高效、可靠和舒適性的需要,我國高速鐵路發(fā)展迅速.但與此同時產(chǎn)生了諸多空氣動力學(xué)問題,包括氣動阻力、列車風、橫風安全性能、列車交會壓力波、氣動噪聲等,這些問題均與高速列車轉(zhuǎn)向架密切相關(guān).而轉(zhuǎn)向架是高速列車結(jié)構(gòu)中最重要的部件之一,位于車輛最下部車體與軌道之間.它牽引和引導(dǎo)車輛沿著軌道行駛,并承受和傳遞來自車體及線路的各種載荷,同時緩和其作用力,它是保證高速列車以200 km/h及以上速度安全平穩(wěn)運行品質(zhì)的關(guān)鍵部件之一[1].研究表明,低速運行時,列車阻力中的氣動阻力所占比例很小,但當速度達到200和300 km/h時,氣動阻力在總阻力中所占的比例將上升到70%和80%左右[2],且轉(zhuǎn)向架阻力占整車氣動阻力的20%以上.高速帶來的噪聲問題更為嚴重,當列車的運行速度超過300 km/h時,氣動噪聲超過輪軌噪聲成為主要噪聲聲源,且轉(zhuǎn)向架是高速列車氣動噪聲的主要聲源部位[3-6].所以,研究高速列車轉(zhuǎn)向架的氣動噪聲特性、主要聲源對總噪聲的貢獻量等對高速列車減阻降噪具有重要的工程參考意義.

    當前多數(shù)對高速列車轉(zhuǎn)向架空氣動力學(xué)性能的研究采用試驗研究和數(shù)值模擬方法.數(shù)值模擬主要針對轉(zhuǎn)向架與整車阻力關(guān)系及整車遠場氣動噪聲特性.由于轉(zhuǎn)向架結(jié)構(gòu)的復(fù)雜性,大多數(shù)數(shù)值計算過多地簡化了轉(zhuǎn)向架結(jié)構(gòu),或者采用縮比模型研究轉(zhuǎn)向架氣動噪聲特性,但對組成轉(zhuǎn)向架各部件的氣動噪聲問題研究甚少.文獻[7]通過對動車轉(zhuǎn)向架的氣動阻力分析,得出頭車第一個轉(zhuǎn)向架阻力是中間車第二個轉(zhuǎn)向架阻力的4倍以上,在側(cè)風作用下轉(zhuǎn)向架阻力占列車總阻力的40%左右.文獻[8]采用分離渦數(shù)值計算方法,研究了橫風中高速列車轉(zhuǎn)向架在時域、頻域的非定常氣動力特性,得到基本氣動力在時域內(nèi)具有隨機波動性,在頻域內(nèi)存在明顯的主頻段1.664~12.990 Hz.文獻[9]通過在頭尾車第一個轉(zhuǎn)向架處安裝裙板并進行數(shù)值研究,結(jié)果表明整車阻力減少較多.文獻[10]重點研究了以簡單轉(zhuǎn)向架(只考慮輪對和枕梁結(jié)構(gòu)的轉(zhuǎn)向架)為噪聲源的高速列車氣動噪聲問題,分析了優(yōu)化轉(zhuǎn)向架裙板后,轉(zhuǎn)向架對遠場噪聲評估點的氣動噪聲降噪效果.文獻[11]基于延遲獨立渦模擬方法,對只包括輪對和構(gòu)架結(jié)構(gòu)的1∶10縮比簡化轉(zhuǎn)向架流場特性和偶極子分布規(guī)律進行了預(yù)測,并通過風洞試驗驗證了數(shù)值模擬結(jié)果的正確性.

    本文進行數(shù)值計算時,同時考慮高速列車動車轉(zhuǎn)向架的細微結(jié)構(gòu),例如橫向減振器、垂向減振器、抗蛇形減振器、空氣彈簧和抗側(cè)滾扭桿等結(jié)構(gòu)的空氣動力學(xué)模型.研究得到動車轉(zhuǎn)向架的漩渦分布特性、氣動噪聲聲源分布特點,以及遠場氣動噪聲的衰減特性、幅值特性、氣動噪聲指向性等規(guī)律,同時得到各噪聲源對動車轉(zhuǎn)向架氣動噪聲的貢獻量.研究成果可為研制設(shè)計高速動車轉(zhuǎn)向架初期的結(jié)構(gòu)優(yōu)化、氣動噪聲改進等提供工程參考依據(jù).

    1 動車轉(zhuǎn)向架氣動噪聲數(shù)值分析理論

    1.1 大渦模擬理論計算方法

    大渦模擬(large eddy simulation,LES)控制方程[12]為不可壓縮粘性流體的N-S方程,LES中大尺度的速度為濾波速度,其定義為

    G(y,y′,Δ)為濾波函數(shù),描述過濾網(wǎng)格的大小.

    假定過濾過程和求導(dǎo)過程可交換,把該函數(shù)用于不可壓縮粘性流體的N-S方程,不管其形式如何,總可得到:

    為使方程組封閉,根據(jù)Smagorinsky的基本亞格子尺度應(yīng)力(sub-grid-scale streese,SGS)模型,假定SGS雷諾應(yīng)力具有如下形式

    式中:δij為單位張量的分量;

    μi為亞格子湍流粘性系數(shù);

    ˉSij為求解尺度下的應(yīng)變張量的分量,

    1.2 計算氣動聲學(xué)方法

    計算氣動聲學(xué)(computational aeroacoustics,CAA)中普遍使用聲類比法,F(xiàn)fowcs Williams-Hawking方程[13](簡稱FW-H方程)的微分形式為

    式中:p′為氣體壓強;

    ni為法向方向;

    a0為聲速;

    vn為法向速度;

    p為靜壓強;

    Tij為Lighthill壓力張量的分量,

    δ(f)為Dirac delta函數(shù);

    H(f)為Heaviside函數(shù).

    FW-H模型表達式表示聲壓是由于瞬態(tài)變化的質(zhì)點力和加速度引起,式(6)右邊分別為單極子聲源、偶極子聲源和四極子聲源.高速列車以一定的速度行駛時,車身表面可以看作是剛性的,體積脈動量幾乎為0,所以可不考慮單極子聲源項[14].文獻[15]指出,流場中四極子聲源與偶極子聲源強度之比正比于馬赫數(shù)的平方,而高速列車的運動仍屬于低速運動,其四極子聲源噪聲強度遠小于偶極子聲源,故四極子聲源項可忽略.即本文只考慮偶極子聲源引起的高速列車動車轉(zhuǎn)向架氣動噪聲問題.

    2 動車轉(zhuǎn)向架氣動噪聲計算模型

    2.1 計算模型

    通常高速列車轉(zhuǎn)向架可分為動力轉(zhuǎn)向架(又稱動車轉(zhuǎn)向架)和非動力轉(zhuǎn)向架(拖車轉(zhuǎn)向架),本文主要以動車轉(zhuǎn)向架為分析對象.動車轉(zhuǎn)向架主要包括輪對、軸箱、一系懸掛,構(gòu)架、二系懸掛、牽引電機和枕梁7部分,結(jié)構(gòu)如圖1所示.

    圖1 動車轉(zhuǎn)向架簡化模型Fig.1 Simplified model of motor car bogie

    2.2 計算區(qū)域、邊界條件及網(wǎng)格劃分

    計算區(qū)域如圖2所示,流場計算區(qū)域長度、寬度和高度分別為25L、12W和6H,動車轉(zhuǎn)向架長度L=3.477 m,寬度W=3.068 m,高度H=1.115 m,輪對踏面與軌道所處地面之間的距離為0.2 m.

    動車轉(zhuǎn)向架來流方向截面為速度入口邊界,正后方截面為壓力出口邊界,左側(cè)、右側(cè)和正上方截面設(shè)置為對稱邊界,動車轉(zhuǎn)向架表面設(shè)置為無滑移壁面的wall邊界.為了模擬地面效應(yīng),地面設(shè)置為滑移地面,其滑移速度為動車轉(zhuǎn)向架的運行速度.

    由于動車轉(zhuǎn)向架結(jié)構(gòu)比較復(fù)雜,因此采用非結(jié)構(gòu)化網(wǎng)格,如圖3所示.計算過程中采用網(wǎng)格自適應(yīng)技術(shù),不斷調(diào)整網(wǎng)格密度,以保證計算精度.為更加精確的考慮動車轉(zhuǎn)向架表面對流體流動的影響,在其表面進行邊界層網(wǎng)格劃分,邊界層增長率為1.2、總厚度為10 mm的5層三棱柱網(wǎng)格,計算域網(wǎng)格量約為2 515萬.

    圖2 計算區(qū)域及邊界設(shè)置Fig.2 Computational domain and boundary setting

    圖3 動車轉(zhuǎn)向架表面網(wǎng)格Fig.3 Surface meshes of motor car bogie

    3 動車轉(zhuǎn)向架氣動噪聲特性

    3.1 動車轉(zhuǎn)向架流場特性

    圖4給出垂直于z軸的縱向速度流線圖.由圖4可以看出,轉(zhuǎn)向架區(qū)域的氣流流速變慢,在有的區(qū)域產(chǎn)生了一定的漩渦.迎風側(cè)軸箱與構(gòu)架之間的外側(cè)、牽引電機后側(cè)、空氣彈簧與抗蛇形減振器之間的臨近區(qū)域均存在不同尺度不同方向的漩渦,尤其是在迎風側(cè)軸箱與構(gòu)架之間、空氣彈簧與抗蛇形減振器的鄰近區(qū)域形成了較大的漩渦,此位置位于構(gòu)架側(cè)梁外側(cè).在進行動車轉(zhuǎn)向架流線型優(yōu)化設(shè)計及減阻降噪時,需要考慮該區(qū)域的流場分布特性.

    3.2 動車轉(zhuǎn)向架氣動噪聲聲源特性

    圖5為動車轉(zhuǎn)向架以350 km/h運行時的聲功率級分布云圖.由圖5可見,迎風側(cè)大部分表面,即輪對、構(gòu)架、牽引電機1、枕梁、垂向減振器和抗側(cè)滾扭桿等聲功率級均達到112 dB以上.其中迎風側(cè)的垂向減振器、抗側(cè)滾扭桿和枕梁的表面聲功率最高.可見,迎風側(cè)凸起部位為動車轉(zhuǎn)向架的氣動噪聲源,且此氣動噪聲源區(qū)域發(fā)生在氣流易分離、湍流運動較劇烈處.

    圖4 動車轉(zhuǎn)向架截面流線圖Fig.4 Streamline diagram in longitudinal sections of the motor car bogie

    圖5 動車轉(zhuǎn)向架聲功率分布云圖Fig.5 Sound power contours of the motor car bogie

    [4]關(guān)于列車表面最大聲功率級Pm與列車運行速度之間的關(guān)系

    式中:a、b為常系數(shù);v0=200 km/h.

    進一步尋求動車轉(zhuǎn)向架最大聲功率與運行速度的對應(yīng)函數(shù)關(guān)系.

    圖6給出動車轉(zhuǎn)向架表面最大聲功率級與運行速度的關(guān)系.由圖6可知,動車轉(zhuǎn)向架的最大聲功率級隨著運行速度的增加而顯著增大,且動車轉(zhuǎn)向架最大聲功率級與運行速度滿足以下函數(shù)關(guān)系:

    動車轉(zhuǎn)向架最大表面聲功率級與運行速度滿足以下函數(shù)關(guān)系

    圖6 最大聲功率級與轉(zhuǎn)向架速度的函數(shù)關(guān)系(v0=200 km/h)Fig.6 Function relationship of the maximum sound power level and motor bogie running speed(v0=200 km/h)

    3.3 動車轉(zhuǎn)向架遠場氣動噪聲特性

    為了研究動車轉(zhuǎn)向架遠場氣動噪聲分布特性,以動車轉(zhuǎn)向架前后輪對中心(距地面高度為0.688 m)為原點,以25 m為半徑,在xz平面以5°間隔布置一周噪聲評估點,在xy平面和yz平面以5°間隔分別布置半周噪聲評估點.

    動車轉(zhuǎn)向架遠場氣動噪聲計算的噪聲評估點布置與坐標平面的對應(yīng)關(guān)系如圖7所示.本文規(guī)定xz平面為垂向平面,xy平面為縱向平面,yz平面為橫向平面,動車轉(zhuǎn)向架正前方、正左方、正后方、正右方測點編號依次為y1、y2、y3、y4,動車轉(zhuǎn)向架正上方測點編號為xz19(橫向平面與縱向平面噪聲評估點的交匯處).

    通常以聲壓級形式對高速鐵路噪聲進行評價,本文采用等效連續(xù)A計權(quán)聲壓級,根據(jù)ISO3095—2005的定義,等效連續(xù)A計權(quán)聲壓級LA的計算式[16]為

    式中:Δt為測量時間間隔,

    Δt=0.5 s;

    pA(t)為瞬時A計權(quán)聲壓,Pa;

    p0為基準聲壓,p0=20 μPa.

    圖7 遠場氣動噪聲評估點分布Fig.7 Distribution of far-field aerodynamic noise evaluation points

    利用快速傅里葉變換將遠場測點的聲壓轉(zhuǎn)換到頻域,并采用A頻域計權(quán)表[17]對頻域進行修正,然后再利用傅里葉逆變換將頻域聲壓轉(zhuǎn)換到時域,即可得到測點的A計權(quán)聲壓pA(t),進而利用式(10)計算出測點的等效連續(xù)A計權(quán)聲壓級.

    圖8為高速列車以200、250、300、350和400 km/h運行時,垂向、縱向和橫向平面噪聲評估點的等效連續(xù)A計權(quán)聲壓級分布對比圖.

    由圖8可見以下特性:

    (1)衰減特性

    運行速度分別為200、250、300、350和400 km/h時,垂向噪聲評估點的等效連續(xù)A計權(quán)聲壓級平均增加幅度為2.96 dBA→2.56 dBA→2.22 dBA→1.94 dBA.說明隨著動車轉(zhuǎn)向架運行速度的增大,垂向平面的噪聲評估點等效連續(xù)A計權(quán)聲壓級增加幅度逐漸減小.其余兩平面均有同樣的分布規(guī)律.

    (2)幅值特性

    各平面噪聲評估點的聲壓級在θ=0°、90°、180°、270°時均出現(xiàn)最大值;垂向平面聲壓級最大值位于θ=90°(噪聲評估點y2)、θ=270°(噪聲評估點y4)處;縱向平面和橫向平面聲壓級最大值均位于θ=90°(噪聲評估點xz19)處,為113.1 dBA(動車轉(zhuǎn)向架以350 km/h運行);不同運行速度下,垂向噪聲評估點最大值和最小值聲壓級相差18.8~21.2 dBA,縱向噪聲評估點最大值和最小值聲壓級相差19.3~21.8 dBA,橫向噪聲評估點最大值和最小值聲壓級相差22.0~26.2 dBA.

    圖8 遠場氣動噪聲評估點的聲壓級分布圖Fig.8 Sound pressure level(SPL)distribution of aerodynamic noise evaluation points

    (3)氣動噪聲指向特性

    動車轉(zhuǎn)向架氣動噪聲是典型的偶極子噪聲,在圓形噪聲評估點中,噪聲輻射表現(xiàn)出很強的對稱性;聲源主要輻射方向為來流枕梁正上方,其次為橫向正左方和正右方,最后為來流正前方和正后方.表明動車轉(zhuǎn)向架的氣動噪聲在垂向、縱向和橫向平面內(nèi)具有明顯的氣動噪聲指向性.

    圖9為動車轉(zhuǎn)向架以350 km/h運行時,噪聲評估點xz19的1/3倍頻程分布圖.由圖9可見,當中心頻率為25和50 Hz時,主頻能量較大.動車轉(zhuǎn)向架的氣動噪聲能量主要集中在高頻部分,隨著運行速度的增大,其主要能量向高頻部分轉(zhuǎn)移且呈現(xiàn)增大趨勢.此外,分析動車轉(zhuǎn)向架以200、250、300和400 km/h運行時的1/3倍頻程可見,中心頻率25和50 Hz為動車轉(zhuǎn)向架的固有中心頻率且在低頻部分能量較大,分布規(guī)律不隨運行速度的改變而變化.

    圖10為動車轉(zhuǎn)向架以350 km/h運行時噪聲評估點xz19的功率譜密度圖.由圖10可見,動車轉(zhuǎn)向架遠場氣動噪聲是寬頻噪聲,在很寬的頻率內(nèi)存在,主要能量集中在28~44.7 kHz頻率范圍內(nèi).當頻率小于11.2 kHz時,功率譜密度的值很小;當頻率在11.2~28 kHz時,功率譜密度隨頻率的增加顯著增大;當頻率大于44.7 kHz時,功率譜密度迅速減小.

    圖9 遠場氣動噪聲評估點的1/3倍頻程Fig.9 1/3 octave band centre frequency of far-field aerodynamic noise evaluation points

    分析動車轉(zhuǎn)向架不同運行速度下的功率譜密度可見,功率譜密度的分布規(guī)律不隨運行速度的改變而變化,其幅值隨運行速度的增加而增大.

    圖10 遠場氣動噪聲評估點的功率譜密度Fig.10 Power spectral density of far-field aerodynamic noise evaluation points

    4 動車轉(zhuǎn)向架遠場氣動噪聲聲源貢獻量分析

    圖11所示為動車轉(zhuǎn)向架各部件分別作為噪聲源以350 km/h運行時得到的等效連續(xù)A計權(quán)聲壓級對比曲線.

    圖11 各部件作為噪聲源的遠場氣動噪聲在垂向平面的聲壓級對比圖Fig.11 SPL comparison of the far-field aerodynamic noise on vertical plane for components of noise sources

    由圖11可看出,采用1個噪聲評估點的聲壓 級評價動車轉(zhuǎn)向架遠場氣動噪聲無意義.因此,可根據(jù)能量疊加原理,采用聲壓級的平均值

    式中:ηi為第i個噪聲評估點所得到的等效連續(xù)A計權(quán)聲壓級;

    I為測點總數(shù).

    在垂向平面中,I=71.由圖11可見:

    (1)構(gòu)架(ηI=85.45 dBA)對動車轉(zhuǎn)向架遠場總噪聲的貢獻最大,其次為輪對(ηI= 79.74 dBA),然后為抗側(cè)滾扭桿(ηI= 79.31 dBA)、垂向減振器(ηI=78.97 dBA)和枕梁(ηI=78.52 dBA).牽引電機1(ηI=77.42 dBA)、牽引電機2(ηI=77.18 dBA)、空氣彈簧(ηI= 75.19 dBA)和橫向減振器(ηI=64.52 dBA)對遠場總噪聲的貢獻量較小.

    (2)對比構(gòu)架與動車轉(zhuǎn)向架的聲壓級可見,在θ=0°和180°處最大聲壓級相差9.9 dBA;在θ= 90°和270°處最大聲壓級相差3.5 dBA.可見構(gòu)架在θ=90°和270°處的遠場噪聲輻射能量最多,指向性較明顯.

    (3)輪對在垂向平面θ=90°和270°處的氣動噪聲指向性很明顯,最大聲壓級達到101.6 dBA,與動車轉(zhuǎn)向架最大聲壓級相差4.3 dBA.

    (4)枕梁的氣動噪聲指向性主要在垂向平面的θ=180°處,在θ=180°處枕梁聲壓級與轉(zhuǎn)向架聲壓級相差12.2 dBA,θ=0°,90°,270°處的最大聲壓級分別相差23.6、13.3和13.3 dBA.

    (5)抗側(cè)滾扭桿在θ=0°,40°,90°,130°,180°,230°,270°,320°處的遠場噪聲輻射較強,且氣動噪聲指向性明顯,平均聲壓級達到79.3 dBA,最大聲壓級與動車轉(zhuǎn)向架聲壓級相差23.5 dBA.

    (6)垂向減振器在垂向平面內(nèi)的氣動噪聲指向性不明顯,其與轉(zhuǎn)向架聲壓級比較,最大聲壓級相差25.3 dBA.

    (7)牽引電機1主要對動車轉(zhuǎn)向架在垂向平面θ=0°,85°,180°,275°處的遠場氣動噪聲輻射有貢獻,牽引電機2主要對動車轉(zhuǎn)向架在垂向平面θ=0°,95°,180°,265°處的遠場氣動噪聲輻射有貢獻,其中在θ=0°和180°處對遠場氣動噪聲輻射的貢獻最大,氣動噪聲指向性很明顯.

    (8)空氣彈簧在垂向平面內(nèi)的遠場氣動噪聲輻射指向性不明確,其與轉(zhuǎn)向架聲壓級比較,最大聲壓級相差29.4 dBA.橫向減振器在垂向平面θ= 0°,90°,180°,270°處的遠場氣動噪聲指向性很明顯,最大聲壓級達到76.5 dBA,與動車轉(zhuǎn)向架最大聲壓級相差32.7 dBA.

    5 結(jié)束語

    本文主要進行了某型動車轉(zhuǎn)向架的氣動性能及氣動噪聲研究,探討了動車轉(zhuǎn)向架的流場特性、氣動噪聲聲源分布、遠場噪聲傳播及動車轉(zhuǎn)向架各部件對遠場總噪聲的貢獻量等問題,對動車轉(zhuǎn)向架氣動噪聲機理、衰減特性、幅值特性及氣動噪聲指向性等均進行了較深入的研究.關(guān)于車體底板等結(jié)構(gòu)對動車轉(zhuǎn)向架氣動噪聲的影響需要進一步研究.

    (1)輪對、構(gòu)架、牽引電機1、枕梁、垂向減振器、抗側(cè)滾扭桿等部件的迎風側(cè)凸起位置為動車轉(zhuǎn)向架的氣動噪聲源,且迎風側(cè)的垂向減振器、抗側(cè)滾扭桿和枕梁對氣動聲源的貢獻量最大.因此,對氣動聲源貢獻較大的部位進行結(jié)構(gòu)設(shè)計,降噪效果會很顯著.

    (2)動車轉(zhuǎn)向架在各平面內(nèi)具有衰減特性、幅值特性和氣動噪聲指向特性等.動車轉(zhuǎn)向架遠場氣動噪聲是寬頻噪聲,主要能量集中在28~44.7 kHz范圍內(nèi).低頻部分的主要能量集中在中心頻率為25和50 Hz附近,且分布規(guī)律不隨運行速度的改變而變化.

    (3)對動車轉(zhuǎn)向架遠場氣動噪聲總噪聲貢獻量最多的部件依次為構(gòu)架、輪對、抗側(cè)滾扭桿、垂向減振器和枕梁.牽引電機1、牽引電機2、空氣彈簧和橫向減振器對遠場總噪聲的貢獻量較小.

    參考文獻:

    [1] 王伯銘.高速動車組總體及轉(zhuǎn)向架[M].成都:西南交通大學(xué)出版社,2014:28-32.

    [2] GAWTHORPE R G.Train drag reduction from simple design changes[J].International Journal of Vehicle Design,1982,3(3):263-274.

    [3] KITAGAWA T,NAGAKURA K.Aerodynamic noise generated by Shinkansen cars[J].Journal of Sound and Vibration,2000,231(5):913-924.

    [4] MELLET C,L TOURNEAUX F,POISSON F,et al. High speed train noise emission:latest investigation of the aerodynamic/rolling noise contribution[J].Journal of Sound and Vibration,2006,293(3):535-546.

    [5] NAGAKURA K.Localization of aerodynamic noise sources of Shinkansen train[J].Journal of Sound and Vibration,2006,293(3):547-556.

    [6] 高陽,王毅剛,王金田,等.聲學(xué)風洞中的高速列車模型氣動噪聲試驗研究[J].聲學(xué)技術(shù),2013,32(6):506-510.

    GAO Yang,WANG Yigang,WANG Jintian,et al. Testing study of aerodynamic noise for high-speed train model in aero-acoustic wind tunnel[J]. Technical Acoustic,2013,32(6):506-510.

    [7] 鄭循皓,張繼業(yè),張衛(wèi)華.高速列車轉(zhuǎn)向架空氣阻力的數(shù)值模擬[J].交通運輸工程學(xué)報,2011,11(2):45-51.

    ZHENG Xunhao,ZHANGJiye,ZHANGWeihua. Numerical simulation of aerodynamic drag for high-speed train bogie[J].Journal of Traffic and Transportation Engineering,2011,11(2):45-51.

    [8] 郗艷紅,毛軍,高亮,等.橫風作用下高速列車轉(zhuǎn)向架非定??諝鈩恿μ匦裕跩].中南大學(xué)學(xué)報:自然科學(xué)版,2014,45(5):1705-1714.

    XIYanhong, MAO Jun, GAO Liang, etal. Aerodynamic force/moment for high-speed train bogie in crosswind field[J].Journal of Central South University:Science and Technology,2014,45(5):1705-1714.

    [9] 楊志剛,高喆,陳羽,等.裙板安裝對高速列車氣動性能影響的數(shù)值分析[J].計算機輔助工程,2010,19(3):16-21.

    YANG Zhigang,GAO Zhe,CHEN Yu,etal. Numericalanalysis on influence on aerodynamic performance of high-speed train caused by installation of skirt plates[J].Computer Aided Engineering,2010,19(3):16-21.

    [10] 黃莎,楊明智,李志偉,等.高速列車轉(zhuǎn)向架部位氣動噪聲數(shù)值模擬及降噪研究[J].中南大學(xué)學(xué)報:自然科學(xué)版,2011,42(12):3899-3904.

    HUANG Sha,YANG Mingzhi,LI Zhiwei,et al. Aerodynamic noise numerical simulation and noise reduction of high-speed train bogie section[J].Journal of Central South University:Science and Technology,2011,42(12):3899-3904.

    [11] ZHU J Y,HU Z W,THOMPSON D J.Analysis of aerodynamic and aeroacoustic behaviour of a simplified high-speed train bogie[J]. Noise and Vibration Mitigation for Rail Transportation Systems,2015,126:489-496.

    [12] GARNIER E,SAGAUT P,ADAMS N.Large eddy simulation for compressible flows[M].[S. l.]:Springer,2009:172-184.

    [13] FFOWCS-WILLIAMS J E,HAWKINGS D L.Sound generation by turbulence and surfaces in arbitrary motion[J].Philosophical Transactions for the Royal Society ofLondon,SeriesA,Mathematicaland Physical Sciences,1969,264(1151):321-342.

    [14] 鄭拯宇,李人憲.高速列車表面氣動噪聲偶極子聲源分布數(shù)值分析[J].西南交通大學(xué)學(xué)報,2011,46(6):996-1002.

    ZHENG Zhengyu,LI Renxian.Numerical analysis of aerodynamic dipole source on high-speed train surface[J].Journal of Southwest Jiaotong University,2011,46(6):996-1002.

    [15] LIGHTHILL M J. On sound generated aerodynamically: part 1: general theory[J]. Proceedings of the Royal Society of London,Series A,Mathematical and Physical Sciences,1952,211:564-587.

    [16] ISO.ISO 3095-2005 Railway application-acousticsmeasurement of noise emitted by railbound vehicles[S].[S.l.]:ISO,2005.

    [17] 馬大猷.現(xiàn)代聲學(xué)理論基礎(chǔ)[M].北京:科學(xué)出版社,2004:57-59.

    高速列車動車轉(zhuǎn)向架氣動噪聲數(shù)值分析

    張亞東, 張繼業(yè), 張 亮, 李 田

    Numerical Analysis of Aerodynamic Noise of Motor Car Bogie for High-Speed Trains

    ZHANG Yadong, ZHANG Jiye, ZHANG Liang, LI Tian
    (State Key Laboratory of Traction Power,Southwest Jiaotong University,Chengdu 610031,China)

    In order to study the aerodynamic noise characteristics of motor car bogie in high-speed trains,an aerodynamic model of motor car bogie was established,a preliminary study on the aerodynamic noise source of the motor car bogie was made using the steady-state RNG k-ε turbulence model and the broadband noise source model,and the far-field aerodynamic noise was analyzed by combined use of the transient-state large eddy simulation(LES)and Lighthill's acoustic analogue theory.The results show that the wheel set,bogie frame,traction motor-1,bolster,vertical shock absorber,anti-rolling torsion bar and other protruding parts on the windward side are the aerodynamic noise source of the bogie.Among them,bogie frame has the largest contribution to the total far-field aerodynamic noise of the bogie,the wheel set and anti-rolling torsion bar have the secondary contribution,and the vertical shock absorber and bolster have the third contribution amount.Compared to the above components,the traction motor-1,traction motor-2,air spring and lateral shock absorber have less contribution to the total noise.In addition,the far-field noise of the motor car bogie is a broadband noise,characterized by attenuation,amplitude and aerodynamic noise directivity.The main energy of the noise in the low frequency band is concentrated at the centre frequencies 25 and 50 Hz,and the power spectral density distribution does not change with the train speed.

    high-speed train;motor car bogie;aerodynamic noise;large eddy simulation;Lighthill's acoustic analogue theory;noise contribution

    張亞東,張繼業(yè),張亮,等.高速列車動車轉(zhuǎn)向架氣動噪聲數(shù)值分析[J].西南交通大學(xué)學(xué)報,2016,51(5):870-877.

    0258-2724(2016)05-870-08

    10.3969/j.issn.0258-2724.2016.05.008

    U270.2;TB115

    A

    2015-03-19

    高速鐵路基礎(chǔ)研究聯(lián)合基金資助項目(U1234208);中央高校基本科研業(yè)務(wù)費資金資助項目(2682014CX042)

    張亞東(1987—),男,博士研究生,研究方向為高速轉(zhuǎn)向架氣動噪聲特性分析及降噪,E-mail:aliyzyd@163.com

    (中文編輯:秦萍玲 英文編輯:蘭俊思)

    猜你喜歡
    遠場聲壓級動車
    機器噪聲平均聲壓級計算方法差異性實證研究
    電動工具(2024年1期)2024-02-29 01:40:24
    一種計算消聲室聲壓級的新方法
    動車西行記
    全新DXR mkll有源揚聲器
    演藝科技(2019年4期)2019-03-30 03:21:46
    基于仿真與實測的列車遠場氣動噪聲分析
    樂!乘動車,看桂林
    某種陣列雷達發(fā)射通道遠場校準簡易方法
    第一次坐動車
    讀寫算(中)(2015年25期)2015-11-07 07:25:43
    Diodes1.9W D類音頻放大器提供高聲壓級水平并延長電池壽命
    戰(zhàn)斗部遠場水下爆炸對艦船沖擊損傷評估
    亚洲人成77777在线视频| 婷婷六月久久综合丁香| 高清在线国产一区| 亚洲av电影不卡..在线观看| 欧洲精品卡2卡3卡4卡5卡区| 脱女人内裤的视频| 中文字幕高清在线视频| 国产伦一二天堂av在线观看| 91大片在线观看| 国产精品 欧美亚洲| 国产亚洲欧美精品永久| 日本在线视频免费播放| 午夜成年电影在线免费观看| 久久久久精品国产欧美久久久| 久久影院123| 亚洲欧洲精品一区二区精品久久久| svipshipincom国产片| 91老司机精品| 香蕉久久夜色| 国产又爽黄色视频| 免费搜索国产男女视频| 亚洲欧美精品综合久久99| 好看av亚洲va欧美ⅴa在| 91字幕亚洲| 国产一区二区三区综合在线观看| 久久这里只有精品19| 国产精品野战在线观看| 一区二区三区高清视频在线| 国产不卡一卡二| 此物有八面人人有两片| 精品久久蜜臀av无| 啦啦啦观看免费观看视频高清 | 欧美日本亚洲视频在线播放| 99在线人妻在线中文字幕| 久久久久久大精品| 日本 欧美在线| 黄片大片在线免费观看| 中亚洲国语对白在线视频| 国产精品久久久久久人妻精品电影| 老鸭窝网址在线观看| 后天国语完整版免费观看| 女警被强在线播放| 亚洲欧美精品综合久久99| 91大片在线观看| 精品国产乱子伦一区二区三区| 窝窝影院91人妻| 成人永久免费在线观看视频| 久久国产精品影院| 大型黄色视频在线免费观看| 久久久精品欧美日韩精品| 国产国语露脸激情在线看| 日韩大尺度精品在线看网址 | 99久久国产精品久久久| 一本大道久久a久久精品| 91成人精品电影| 久久中文字幕一级| 美女扒开内裤让男人捅视频| 少妇熟女aⅴ在线视频| 国产精品久久久av美女十八| 久久久久久久精品吃奶| 久久人妻av系列| 亚洲欧美激情在线| 99精品欧美一区二区三区四区| 久久性视频一级片| 中文字幕高清在线视频| 涩涩av久久男人的天堂| 男女之事视频高清在线观看| 国内久久婷婷六月综合欲色啪| 国产精品1区2区在线观看.| 精品国产美女av久久久久小说| 色精品久久人妻99蜜桃| 黄片播放在线免费| 亚洲专区中文字幕在线| 嫩草影院精品99| e午夜精品久久久久久久| 91九色精品人成在线观看| 日本免费一区二区三区高清不卡 | 日韩国内少妇激情av| 美女午夜性视频免费| 大型av网站在线播放| 午夜成年电影在线免费观看| 老司机靠b影院| 黄色视频不卡| x7x7x7水蜜桃| 欧美激情久久久久久爽电影 | 国内精品久久久久精免费| 女人被躁到高潮嗷嗷叫费观| 一个人免费在线观看的高清视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品国产高清国产av| 男女下面进入的视频免费午夜 | 一级作爱视频免费观看| 精品久久久精品久久久| 欧美精品亚洲一区二区| 欧美日本中文国产一区发布| 十分钟在线观看高清视频www| 国产精品亚洲美女久久久| 搡老岳熟女国产| 日韩av在线大香蕉| 国内精品久久久久精免费| 好男人电影高清在线观看| 国产精品爽爽va在线观看网站 | 国产人伦9x9x在线观看| 国产欧美日韩一区二区三区在线| 久久人妻熟女aⅴ| 午夜久久久在线观看| av免费在线观看网站| 午夜免费观看网址| 国产麻豆69| 日韩国内少妇激情av| 国产高清激情床上av| 91国产中文字幕| 欧美最黄视频在线播放免费| 成人三级黄色视频| 又大又爽又粗| 精品久久久久久久人妻蜜臀av | 别揉我奶头~嗯~啊~动态视频| 看黄色毛片网站| 色在线成人网| 一本久久中文字幕| 免费一级毛片在线播放高清视频 | 成熟少妇高潮喷水视频| 一级a爱视频在线免费观看| 亚洲男人天堂网一区| 欧美日韩福利视频一区二区| 亚洲va日本ⅴa欧美va伊人久久| 国产不卡一卡二| 99riav亚洲国产免费| 一级毛片女人18水好多| 天天躁狠狠躁夜夜躁狠狠躁| 久久午夜亚洲精品久久| 色婷婷久久久亚洲欧美| 岛国视频午夜一区免费看| 久99久视频精品免费| 搡老熟女国产l中国老女人| 国产真人三级小视频在线观看| 免费看十八禁软件| 91在线观看av| 国产精品久久久久久精品电影 | 可以在线观看毛片的网站| av天堂在线播放| 满18在线观看网站| 精品国内亚洲2022精品成人| 国产三级在线视频| 一区二区三区精品91| 窝窝影院91人妻| 午夜影院日韩av| 亚洲五月天丁香| 国产在线精品亚洲第一网站| 波多野结衣高清无吗| 一夜夜www| 欧美大码av| 男人舔女人的私密视频| 一夜夜www| 国产精品一区二区在线不卡| 最好的美女福利视频网| 在线观看免费日韩欧美大片| 成人永久免费在线观看视频| 99热只有精品国产| 久久草成人影院| 女同久久另类99精品国产91| 国产精品自产拍在线观看55亚洲| 中文字幕最新亚洲高清| 涩涩av久久男人的天堂| 欧美成狂野欧美在线观看| 成人精品一区二区免费| 久久狼人影院| 久久午夜亚洲精品久久| 岛国在线观看网站| 久久香蕉国产精品| 成熟少妇高潮喷水视频| 自线自在国产av| 男女之事视频高清在线观看| 窝窝影院91人妻| 巨乳人妻的诱惑在线观看| 99riav亚洲国产免费| 亚洲中文日韩欧美视频| 婷婷六月久久综合丁香| 国产区一区二久久| 国产一区二区三区视频了| 大陆偷拍与自拍| 丰满人妻熟妇乱又伦精品不卡| 12—13女人毛片做爰片一| 国内毛片毛片毛片毛片毛片| 精品国产乱子伦一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 精品熟女少妇八av免费久了| 搞女人的毛片| 午夜福利,免费看| 日本撒尿小便嘘嘘汇集6| 亚洲av五月六月丁香网| 久久久久久久午夜电影| 正在播放国产对白刺激| 免费少妇av软件| 午夜福利视频1000在线观看 | 欧美成人一区二区免费高清观看 | 男女下面插进去视频免费观看| 国产精品一区二区免费欧美| 亚洲精品一区av在线观看| 中亚洲国语对白在线视频| 波多野结衣巨乳人妻| 一a级毛片在线观看| 真人一进一出gif抽搐免费| 一本综合久久免费| 国产一区二区在线av高清观看| 大码成人一级视频| 久久欧美精品欧美久久欧美| 91国产中文字幕| 欧美乱色亚洲激情| 欧美日本中文国产一区发布| 亚洲天堂国产精品一区在线| 在线视频色国产色| 三级毛片av免费| 国产乱人伦免费视频| 久久久久久国产a免费观看| 午夜成年电影在线免费观看| 韩国精品一区二区三区| 视频在线观看一区二区三区| 一本大道久久a久久精品| 国产精品二区激情视频| 亚洲av电影不卡..在线观看| 亚洲伊人色综图| 一边摸一边抽搐一进一小说| 搞女人的毛片| 90打野战视频偷拍视频| 久久青草综合色| 午夜福利欧美成人| 国产亚洲欧美98| 黄网站色视频无遮挡免费观看| 国内毛片毛片毛片毛片毛片| 午夜免费鲁丝| 99re在线观看精品视频| 国产亚洲欧美98| 色哟哟哟哟哟哟| 国产成人精品无人区| 日日爽夜夜爽网站| 午夜福利,免费看| 久久草成人影院| ponron亚洲| 一个人观看的视频www高清免费观看 | 国产精品二区激情视频| 女人被狂操c到高潮| 不卡一级毛片| www.www免费av| 一本大道久久a久久精品| 国产免费av片在线观看野外av| 欧美色欧美亚洲另类二区 | 精品一区二区三区四区五区乱码| 国产99白浆流出| 搡老妇女老女人老熟妇| 精品卡一卡二卡四卡免费| 亚洲视频免费观看视频| 婷婷丁香在线五月| 村上凉子中文字幕在线| 久久久精品国产亚洲av高清涩受| 黄网站色视频无遮挡免费观看| 亚洲七黄色美女视频| 真人做人爱边吃奶动态| 久久久水蜜桃国产精品网| 日韩欧美在线二视频| 波多野结衣高清无吗| 91国产中文字幕| 国产精品国产高清国产av| 露出奶头的视频| 50天的宝宝边吃奶边哭怎么回事| 免费av毛片视频| 亚洲一区高清亚洲精品| 亚洲av成人不卡在线观看播放网| 一区在线观看完整版| 老司机在亚洲福利影院| 两个人视频免费观看高清| 一级毛片高清免费大全| 国产91精品成人一区二区三区| 女人被狂操c到高潮| 999久久久国产精品视频| 亚洲精品久久国产高清桃花| 欧美日本视频| 啦啦啦 在线观看视频| 国产精品永久免费网站| av欧美777| 在线av久久热| 99国产极品粉嫩在线观看| 亚洲色图 男人天堂 中文字幕| 免费无遮挡裸体视频| 大陆偷拍与自拍| 国产精品电影一区二区三区| 免费高清视频大片| 高潮久久久久久久久久久不卡| 男女下面进入的视频免费午夜 | 中文字幕最新亚洲高清| 在线永久观看黄色视频| 大码成人一级视频| 国产精品久久久av美女十八| 久久这里只有精品19| 麻豆久久精品国产亚洲av| 亚洲精品久久国产高清桃花| 久久国产精品影院| 18美女黄网站色大片免费观看| 久久精品亚洲精品国产色婷小说| 黑丝袜美女国产一区| 亚洲精品av麻豆狂野| 婷婷六月久久综合丁香| 日本三级黄在线观看| 亚洲国产精品999在线| 极品人妻少妇av视频| 少妇粗大呻吟视频| 国产熟女午夜一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 久久影院123| 国产高清激情床上av| 大码成人一级视频| 欧美在线黄色| 欧美黄色淫秽网站| 欧美日本亚洲视频在线播放| 久久久精品国产亚洲av高清涩受| 色哟哟哟哟哟哟| 高清毛片免费观看视频网站| 午夜两性在线视频| 久久性视频一级片| 亚洲一区二区三区色噜噜| 18禁裸乳无遮挡免费网站照片 | 亚洲 国产 在线| 亚洲色图 男人天堂 中文字幕| 黑人巨大精品欧美一区二区mp4| 69av精品久久久久久| 日韩国内少妇激情av| 丝袜美腿诱惑在线| 国产成+人综合+亚洲专区| 国产欧美日韩精品亚洲av| 亚洲精品国产色婷婷电影| 久久亚洲精品不卡| 日本在线视频免费播放| 色综合亚洲欧美另类图片| 成人亚洲精品av一区二区| 91成人精品电影| 老司机午夜十八禁免费视频| 日本五十路高清| 非洲黑人性xxxx精品又粗又长| 亚洲国产欧美日韩在线播放| 欧美成狂野欧美在线观看| 亚洲三区欧美一区| 亚洲五月婷婷丁香| 欧美亚洲日本最大视频资源| 精品熟女少妇八av免费久了| 18禁观看日本| 中文字幕人妻熟女乱码| 成人国产一区最新在线观看| 热99re8久久精品国产| 亚洲全国av大片| 久9热在线精品视频| 日本黄色视频三级网站网址| 青草久久国产| 麻豆久久精品国产亚洲av| 如日韩欧美国产精品一区二区三区| 欧美成人午夜精品| 欧美乱妇无乱码| 在线观看免费视频日本深夜| 高清毛片免费观看视频网站| 国产一区二区三区视频了| 国产色视频综合| 国产午夜福利久久久久久| 成人手机av| 亚洲熟女毛片儿| 在线观看午夜福利视频| 国产午夜福利久久久久久| 欧美乱妇无乱码| 久久香蕉精品热| 久久久久九九精品影院| 亚洲国产精品合色在线| 国产精品九九99| e午夜精品久久久久久久| 禁无遮挡网站| 真人做人爱边吃奶动态| www.自偷自拍.com| 亚洲av成人一区二区三| 亚洲欧美日韩无卡精品| 亚洲人成电影观看| 99精品在免费线老司机午夜| 巨乳人妻的诱惑在线观看| 亚洲国产欧美一区二区综合| 波多野结衣巨乳人妻| 搡老岳熟女国产| 国产欧美日韩一区二区精品| 欧美日本视频| 在线免费观看的www视频| 久久这里只有精品19| 国产熟女xx| 久久久久久人人人人人| 免费av毛片视频| 中国美女看黄片| 亚洲 欧美一区二区三区| 亚洲国产毛片av蜜桃av| 中国美女看黄片| 国产日韩一区二区三区精品不卡| 婷婷精品国产亚洲av在线| 99久久国产精品久久久| 午夜影院日韩av| 亚洲成国产人片在线观看| 美女免费视频网站| 亚洲男人的天堂狠狠| 国产蜜桃级精品一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 国产一卡二卡三卡精品| 夜夜躁狠狠躁天天躁| 天堂动漫精品| 国产成人精品久久二区二区免费| 国内毛片毛片毛片毛片毛片| 欧美日韩精品网址| 欧美乱妇无乱码| 麻豆久久精品国产亚洲av| 免费在线观看黄色视频的| 天天躁夜夜躁狠狠躁躁| 女人被狂操c到高潮| 又紧又爽又黄一区二区| 这个男人来自地球电影免费观看| 欧美亚洲日本最大视频资源| 精品一区二区三区四区五区乱码| 国产成人系列免费观看| 国产成人啪精品午夜网站| 精品国产美女av久久久久小说| 欧美久久黑人一区二区| 19禁男女啪啪无遮挡网站| 久久久久久久久中文| 欧美激情极品国产一区二区三区| 国产精品一区二区免费欧美| 一区福利在线观看| 狂野欧美激情性xxxx| 老汉色∧v一级毛片| 狠狠狠狠99中文字幕| 婷婷丁香在线五月| 婷婷丁香在线五月| 中文亚洲av片在线观看爽| 久久久久久久久久久久大奶| 欧美绝顶高潮抽搐喷水| 18禁国产床啪视频网站| 国产精品美女特级片免费视频播放器 | 国内精品久久久久久久电影| 午夜视频精品福利| 伦理电影免费视频| 久久香蕉精品热| 在线天堂中文资源库| 极品人妻少妇av视频| 搞女人的毛片| 黄色 视频免费看| 色综合站精品国产| 美女扒开内裤让男人捅视频| 老司机深夜福利视频在线观看| 亚洲精品一区av在线观看| 好男人在线观看高清免费视频 | 久久精品aⅴ一区二区三区四区| 亚洲人成电影观看| 女同久久另类99精品国产91| 国产av又大| ponron亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 成人av一区二区三区在线看| 精品无人区乱码1区二区| 国产野战对白在线观看| 亚洲avbb在线观看| 日韩欧美三级三区| 中文字幕久久专区| 久久久精品欧美日韩精品| 又黄又粗又硬又大视频| 一级片免费观看大全| 亚洲人成网站在线播放欧美日韩| 亚洲成av人片免费观看| 自线自在国产av| 久久这里只有精品19| 日韩成人在线观看一区二区三区| 国产精品久久电影中文字幕| 精品乱码久久久久久99久播| 日日摸夜夜添夜夜添小说| 国产精品自产拍在线观看55亚洲| 在线十欧美十亚洲十日本专区| a级毛片在线看网站| 日韩中文字幕欧美一区二区| 欧美日韩福利视频一区二区| 久久人人精品亚洲av| 亚洲av电影不卡..在线观看| 国产精品,欧美在线| 欧美最黄视频在线播放免费| 99国产综合亚洲精品| 女性被躁到高潮视频| 中国美女看黄片| 亚洲成人久久性| 亚洲人成伊人成综合网2020| 巨乳人妻的诱惑在线观看| 正在播放国产对白刺激| 黄色视频不卡| 国产免费男女视频| 日韩av在线大香蕉| 亚洲中文字幕一区二区三区有码在线看 | 欧美大码av| 久久香蕉国产精品| 不卡av一区二区三区| 一级,二级,三级黄色视频| 亚洲电影在线观看av| 十八禁人妻一区二区| 老熟妇仑乱视频hdxx| 国产97色在线日韩免费| 在线永久观看黄色视频| 可以在线观看毛片的网站| 国产精品av久久久久免费| 免费无遮挡裸体视频| 禁无遮挡网站| 国产三级在线视频| 色尼玛亚洲综合影院| 十八禁网站免费在线| 黑人操中国人逼视频| 母亲3免费完整高清在线观看| 国产午夜福利久久久久久| 国产精品秋霞免费鲁丝片| 国产又爽黄色视频| 黄色a级毛片大全视频| or卡值多少钱| 欧美日韩精品网址| 久久久国产成人免费| 村上凉子中文字幕在线| 正在播放国产对白刺激| 精品久久久久久久毛片微露脸| 欧美性长视频在线观看| 一区二区三区精品91| 嫩草影视91久久| 午夜福利成人在线免费观看| 99精品欧美一区二区三区四区| 看黄色毛片网站| 男女下面插进去视频免费观看| 激情视频va一区二区三区| 97超级碰碰碰精品色视频在线观看| 午夜a级毛片| 国内精品久久久久精免费| 国内毛片毛片毛片毛片毛片| 精品少妇一区二区三区视频日本电影| 黄色片一级片一级黄色片| 国产真人三级小视频在线观看| 亚洲情色 制服丝袜| 日韩视频一区二区在线观看| av免费在线观看网站| 99精品久久久久人妻精品| 成人三级黄色视频| 亚洲情色 制服丝袜| 91av网站免费观看| 亚洲欧美激情在线| 日本vs欧美在线观看视频| 男女做爰动态图高潮gif福利片 | 午夜精品国产一区二区电影| 最近最新免费中文字幕在线| 亚洲精品av麻豆狂野| 免费看a级黄色片| 男人的好看免费观看在线视频 | 亚洲无线在线观看| 国产精品永久免费网站| 久久久久久国产a免费观看| 久久久久国产精品人妻aⅴ院| 两性夫妻黄色片| 精品久久久久久,| 久久国产精品人妻蜜桃| 国产成人精品在线电影| 侵犯人妻中文字幕一二三四区| 国产熟女午夜一区二区三区| 亚洲欧美激情在线| 亚洲黑人精品在线| 无限看片的www在线观看| 国产av精品麻豆| 99国产精品99久久久久| 国产精品香港三级国产av潘金莲| 18禁国产床啪视频网站| 级片在线观看| 99精品欧美一区二区三区四区| 国产在线精品亚洲第一网站| 少妇被粗大的猛进出69影院| 欧美中文日本在线观看视频| 国产精品一区二区精品视频观看| 女同久久另类99精品国产91| 在线免费观看的www视频| 国产高清videossex| 国产高清激情床上av| 中文字幕色久视频| 日韩欧美国产一区二区入口| 别揉我奶头~嗯~啊~动态视频| 老司机在亚洲福利影院| 免费女性裸体啪啪无遮挡网站| 精品久久久久久久人妻蜜臀av | 老汉色av国产亚洲站长工具| 黄色毛片三级朝国网站| 国产精品精品国产色婷婷| 欧美精品啪啪一区二区三区| 午夜亚洲福利在线播放| 少妇熟女aⅴ在线视频| 亚洲国产精品成人综合色| 成年版毛片免费区| 久久精品91无色码中文字幕| 又黄又粗又硬又大视频| 欧美日韩亚洲综合一区二区三区_| 亚洲国产欧美网| 人人妻人人爽人人添夜夜欢视频| 最近最新中文字幕大全电影3 | 久久人人97超碰香蕉20202| 国产麻豆69| 久久人人97超碰香蕉20202| 亚洲人成伊人成综合网2020| 妹子高潮喷水视频| 视频在线观看一区二区三区| 久热爱精品视频在线9| 国产成人精品无人区| 老司机深夜福利视频在线观看| 露出奶头的视频| 最新美女视频免费是黄的| 午夜a级毛片| 高清在线国产一区| 国产精品av久久久久免费| 成人av一区二区三区在线看| 久久久久精品国产欧美久久久| 99国产综合亚洲精品| 少妇粗大呻吟视频|