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

    翼柱型與環(huán)向開槽型燃燒室聲學(xué)特性對比

    2021-09-07 12:03:46王大鵬趙靜陳林君楊奔馬瀟健中國運(yùn)載火箭技術(shù)研究院研究發(fā)展部北京00076哈爾濱工程大學(xué)航天與建筑工程學(xué)院黑龍江哈爾濱5000
    關(guān)鍵詞:聲頻環(huán)槽環(huán)向

    王大鵬, 趙靜, 陳林君, 楊奔, 馬瀟健 (.中國運(yùn)載火箭技術(shù)研究院研究發(fā)展部, 北京 00076;.哈爾濱工程大學(xué) 航天與建筑工程學(xué)院,黑龍江 哈爾濱 5000)

    火箭發(fā)動機(jī)中的不穩(wěn)定燃燒現(xiàn)象自得到觀測起便引起了學(xué)者們的廣泛關(guān)注,周期性的壓力振蕩是不穩(wěn)定燃燒現(xiàn)象的重要表現(xiàn)。即使在發(fā)動機(jī)正常工作時,燃燒室內(nèi)仍會存在壓力振蕩,而大振幅壓力振蕩很有可能引起嚴(yán)重的結(jié)構(gòu)損壞和災(zāi)難性的發(fā)動機(jī)故障[1]。因此不穩(wěn)定燃燒一直以來都是火箭發(fā)動機(jī)領(lǐng)域最具挑戰(zhàn)性的問題之一[2]。不穩(wěn)定燃燒按照其壓力振蕩頻率與燃燒室聲頻的關(guān)系可以分為聲不穩(wěn)定和非聲不穩(wěn)定。我國在對大長徑比翼柱裝藥固體火箭發(fā)動機(jī)的研制中出現(xiàn)了聲不穩(wěn)定現(xiàn)象[3]。美國末段高空防御系統(tǒng)(terminal high altitude area defense,THAAD)反導(dǎo)固體發(fā)動機(jī)采用了一種特殊的環(huán)向開槽型裝藥[4]結(jié)構(gòu),并且未見有關(guān)其不穩(wěn)定燃燒的報(bào)道。對以上2種典型裝藥結(jié)構(gòu)展開對比研究,或可從裝藥結(jié)構(gòu)角度為我國戰(zhàn)術(shù)火箭發(fā)動機(jī)不穩(wěn)定燃燒問題找到解決途徑。

    發(fā)動機(jī)燃燒室常作為聲腔進(jìn)行分析,通過有限元技術(shù)可以獲得燃燒室的聲模態(tài)[5],也有學(xué)者通過揚(yáng)聲器激勵測量燃燒室內(nèi)的聲學(xué)響應(yīng)特性[5-6]。在美國海軍與空軍共同開發(fā)的聲不穩(wěn)定性預(yù)估軟件SSP[7-9]中,也包含發(fā)動機(jī)的聲學(xué)特性分析。我國學(xué)者也在燃燒不穩(wěn)定性的聲學(xué)問題上取得了一定成果:樊曉波[10]研究了聲抑制裝置對燃燒不穩(wěn)定的抑制機(jī)理;楊向明等[11]對翼柱型固體發(fā)動機(jī)燃燒室聲場進(jìn)行分析,得出了噴管潛入段空腔的阻尼作用;張嶠等[12]對H型聲腔的固有頻率進(jìn)行了計(jì)算,得到了該類型聲腔縱向聲頻的變化規(guī)律,并且在對發(fā)動機(jī)頭部空腔的研究中,使用了有限元手段與CFD仿真結(jié)果進(jìn)行了相互驗(yàn)證[13];蘇萬興[14]在對大長徑比固體火箭發(fā)動機(jī)的渦-聲耦合研究中,將燃燒室聲頻與渦脫落頻率進(jìn)行了對比,表明某翼柱型發(fā)動機(jī)燃燒末段的燃燒不穩(wěn)定不是由渦-聲耦合所導(dǎo)致的。

    對通過聲模態(tài)分析可迅速獲得確定聲壓波形和頻率[1-2],但所得聲振幅是相對值,不能得到發(fā)動機(jī)燃燒室在激勵下的聲學(xué)響應(yīng)特性。推進(jìn)劑的壓力耦合響應(yīng)、速度耦合響應(yīng)[15]以及發(fā)動機(jī)結(jié)構(gòu)振動[16-17]會使燃速發(fā)生變化,漩渦脫落[18]與內(nèi)流道微小彈性形變[19]也會導(dǎo)致燃燒室壓力的變化,最終影響燃燒室聲腔響應(yīng)。Wang等[20]通過聲學(xué)有限元方法對燃燒室的響應(yīng)聲場進(jìn)行了研究,獲得了聲壓在時間-空間的分布規(guī)律。

    為解決在大長徑比翼柱裝藥發(fā)動機(jī)研制過程中出現(xiàn)的不穩(wěn)定燃燒問題,對燃燒室在特定激勵下的聲學(xué)響應(yīng)展開研究是迫切需要的。本文從聲學(xué)角度出發(fā),對比分析了典型翼柱型和環(huán)向開槽型裝藥發(fā)動機(jī)在燃面退移過程中燃燒室聲頻和聲學(xué)響應(yīng)特性變化規(guī)律,并給出了裝藥結(jié)構(gòu)設(shè)計(jì)的一些建議。

    1 聲學(xué)有限元方法

    理論方法計(jì)算聲場的固有頻率和聲振型一般只應(yīng)用于簡單的聲腔形狀,發(fā)動機(jī)燃燒室聲腔形狀多為較復(fù)雜的結(jié)構(gòu),解析法求解三維波動方程極為復(fù)雜,不同結(jié)構(gòu)需采用不同的修正系數(shù),準(zhǔn)確性也不能保證,而低維的求解不能表征其物理特性[21]。本文采用有限元方法,對不同內(nèi)流道形狀聲腔進(jìn)行模擬。該方法的優(yōu)點(diǎn)是可以求解復(fù)雜幾何形狀聲腔的固有頻率和聲振型,且可以模擬燃燒室中不同位置在頻域的響應(yīng)特性。

    1.1 控制方程

    根據(jù)質(zhì)量守恒可以得到流體的連續(xù)性方程:

    (1)

    式中:v是速度向量;ρ是密度;Q是質(zhì)量源項(xiàng)。

    根據(jù)動量守恒,納維-斯托克斯方程為:

    (2)

    式中:b是體積力;S是粘性應(yīng)力張量;p是壓力。

    使用下列假設(shè)可以將方程化簡并得到聲波方程:1)流體是可壓縮的(密度是壓力的函數(shù));2)流體是無旋的;3)無體積力;4)流體的壓力擾動很??;5)在流體中沒有平均流動;6)流體為理想氣體,絕熱并且可逆。

    線性化的連續(xù)性方程為:

    (3)

    線性化的納維-斯托克斯方程為:

    (4)

    式中:va是聲學(xué)速度;pa是聲壓;ρ0是流體的平均密度。

    聲波在流體中的有損傳播方程可以表示為:

    (5)

    使用伽遼金法對式(5)進(jìn)行處理,與測試函數(shù)w相乘并且在體積上積分可以獲得[22]:

    (6)

    式中:dV是聲學(xué)域ΩF的微分體積;dS是聲學(xué)域邊界ΓF的微分面積;n是聲學(xué)域邊界ΓF的單位外法向量。

    從動量守恒方程可以得到聲學(xué)域邊界的法向速度:

    (7)

    將流體的法向加速度使用位移表示,可得:

    (8)

    將式(7)、(8)代入到式(6)中可以得到式(5)的弱形式:

    (9)

    式中:流體壓力p和位移uF作為待求解量。

    使用有限元近似得到壓力和位移的形函數(shù)為:

    p=NTpe

    (10)

    u=N′Tue

    (11)

    式中:N是壓力的單元形函數(shù);N′是位移的單元形函數(shù);pe是節(jié)點(diǎn)壓力向量;ue是節(jié)點(diǎn)位移向量。

    根據(jù)式(10)、(11)待求解量對時間的二階導(dǎo)數(shù)以及壓力的變化為:

    (12)

    (13)

    δp=NTδpe

    (14)

    將式(10)、(11)代入到式(9)中可以得到有限元形式的聲波方程:

    (15)

    式中:n是流體邊界的外法向量;q是節(jié)點(diǎn)質(zhì)量源向量。

    (16)

    式中:MF是聲質(zhì)量矩陣;CF是聲流體阻尼矩陣;KF是聲流體剛度矩陣;RT是聲流體邊界矩陣;fF是聲流體載荷向量。

    1.2 幾何模型及網(wǎng)格劃分

    文中的研究對象是2種典型裝藥的固體火箭發(fā)動機(jī)內(nèi)流場。根據(jù)燃面的平行退移假設(shè)分別建立了2種內(nèi)流道結(jié)構(gòu)從初始時刻到結(jié)束工作時刻共計(jì)22個幾何模型。為方便后文討論,將發(fā)動機(jī)的工作時刻做歸一化處理:

    Tn=ta/td

    (17)

    式中:Tn為歸一化工作時刻;ta為目前工作時刻;td為設(shè)計(jì)工作時長。

    Tn=0時的內(nèi)流道幾何模型如圖1所示,喉部直徑L0=78 mm。由于發(fā)動機(jī)工作時喉部為聲速,喉部下游的聲學(xué)振蕩不會傳遞回燃燒室,采用文獻(xiàn)[20]對模型的處理方法移除了噴管擴(kuò)張段。并且由于發(fā)動機(jī)頭部位置和后封蓋位置在多種軸向振型下始終為聲壓波腹[20],位置對聲不穩(wěn)定下的壓力振蕩更為敏感,定義位于發(fā)動機(jī)頭部附近(0.013, 0.04, 0 m)為監(jiān)測點(diǎn)P1;P2位于發(fā)動機(jī)后封蓋附近(3.238, 0.125, 0 m)為監(jiān)測點(diǎn)P2。

    圖1 Tn=0時內(nèi)流道幾何模型Fig.1 Geometric model of internal flow-field when Tn=0

    1.3 燃?xì)鈪?shù)及邊界條件

    表1 典型燃?xì)馕镄詤?shù)Table 1 Gas property parameters

    1.4 聲學(xué)有限元方法驗(yàn)證

    為驗(yàn)證本文方法,對馮卡門流體力學(xué)研究所(Von Karman institute for fluid dynamics,VKI)采用的P230發(fā)動機(jī)縮比模型[23-24]進(jìn)行了有限元仿真,并與結(jié)果進(jìn)行了對比。為與文獻(xiàn)保持一致,驗(yàn)證中使用溫度為285 K的空氣作為聲流體。模型喉部截面定義為聲輻射阻抗邊界,其余外表面定義為剛性壁面,所有表面無位移[20]。如圖2所示,模型燃燒室2 000 Hz以下共有4階聲模態(tài),并且都為軸向聲模態(tài)。模型燃燒室頭部位置與潛入式噴管空前尾端始終是軸向聲模態(tài)聲壓的波腹。驗(yàn)證工況結(jié)果與文獻(xiàn)[20]中的聲壓振型進(jìn)行了理論計(jì)算聲壓分布規(guī)律一致。如表2所示,驗(yàn)證工況的聲頻計(jì)算結(jié)果與對比組[23-24]結(jié)果的一致性良好。前4階聲頻的數(shù)值驗(yàn)證結(jié)果與實(shí)驗(yàn)值依次相差1.9%、3.9%、1.0%和2.5%,現(xiàn)有模型計(jì)算結(jié)果是合理并且有效的。因此,文中使用的聲學(xué)有限元方法能夠?qū)θ紵衣晥鲞M(jìn)行合理且有效地預(yù)測。

    圖2 模態(tài)相位角φ=0°時,模型燃燒室的前4階聲模態(tài)振型Fig.2 First four acoustic mode shapes of the computational model, when modal phase angle φ=0°

    表2 模型燃燒室前4階聲模態(tài)頻率與文獻(xiàn)對比

    2 聲學(xué)特性仿真結(jié)果與討論

    2.1 聲模態(tài)頻率隨燃面退移的變化

    對2種結(jié)構(gòu)發(fā)動機(jī)不同時刻下的燃燒室做聲模態(tài)分析,如圖3所示。工作時間內(nèi)翼柱型發(fā)動機(jī)與環(huán)向開槽型發(fā)動機(jī)前5階聲模態(tài)均為軸向模態(tài)。對于一定幅值不同頻率的擾動對燃燒室的激勵,當(dāng)擾動頻率接近燃燒室聲頻時燃燒室內(nèi)壓力振幅顯著增大[20]。因此,在工程設(shè)計(jì)角度上,燃燒室聲頻變化范圍不應(yīng)太大。小的聲頻變化范圍更有助于通過結(jié)構(gòu)設(shè)計(jì)避免聲頻與結(jié)構(gòu)頻率接近所引起的共振現(xiàn)象。因此,本節(jié)重點(diǎn)討論2種結(jié)構(gòu)發(fā)動機(jī)隨工作時間的聲頻變化。

    圖3中顯示了翼柱型發(fā)動機(jī)和環(huán)向開槽型發(fā)動機(jī)隨工作時間的第1階聲模態(tài)振型。翼柱型發(fā)動機(jī)隨工作時間翼結(jié)構(gòu)不斷變寬加深,并形成一個整體突擴(kuò)腔,并且腔的容積逐漸增大。環(huán)向開槽型發(fā)動機(jī)槽的深度逐漸變淺,寬度逐漸變大。處在尾部的槽結(jié)構(gòu)在工作中期已經(jīng)變得不明顯,并逐漸消失。燃燒室內(nèi)尾端的擴(kuò)張結(jié)構(gòu)會使得發(fā)動機(jī)的聲學(xué)穩(wěn)定性變差[25]。而環(huán)槽裝藥發(fā)動機(jī)中,環(huán)槽特征結(jié)構(gòu)在工作中期的逐漸退化可能有聲學(xué)穩(wěn)定性的考量。

    圖3 相位角φ=0°時,翼柱裝藥發(fā)動機(jī)聲模態(tài)基振型隨工作時間變化Fig.3 First order acoustic mode shape of finocyl motor over time when phase angle φ=0°

    對2種結(jié)構(gòu)發(fā)動機(jī)不同時刻下的燃燒室做聲模態(tài)分析可以得到其隨時間變化的前5階聲頻,如圖4。本文中計(jì)算得到的2種發(fā)動機(jī)各時刻前5階聲模態(tài)均為軸向模態(tài),文中所討論的所有聲頻均為軸向聲模態(tài)頻率。

    如圖4,隨著工作過程燃面退移,翼柱型發(fā)動機(jī)的前2階聲頻先降低再升高,聲頻最低值出現(xiàn)在Tn=0.3時刻。翼柱型發(fā)動機(jī)的第3~5階聲頻均存在著隨工作時間先升高再降低的規(guī)律,聲頻最高值出現(xiàn)在Tn=0.3附近。環(huán)向開槽型發(fā)動機(jī)的前5階聲頻均具有隨工作時間先降低再升高的變化規(guī)律,聲基頻最低值均出現(xiàn)在Tn=0.3,第2~5階聲頻最低值出現(xiàn)在Tn=0.4時刻附近。2種裝藥發(fā)動機(jī)由于特征結(jié)構(gòu)不同,導(dǎo)致其高階聲頻的變化規(guī)律也不一致。2種裝藥發(fā)動機(jī)的前2階聲頻最低值均出現(xiàn)在工作時間的中前段。圖4中,翼柱型發(fā)動機(jī)的聲頻高于環(huán)向開槽型發(fā)動機(jī)。2種發(fā)動機(jī)聲頻在相同階數(shù)下均有差別,前5階聲頻依次相差7.9、18.3、36.1、76.4、79.3 Hz。該結(jié)果表明,在相同長度不同特征結(jié)構(gòu)的燃燒室中,特征結(jié)構(gòu)對聲頻有影響。此外,聲頻相差的最大值都出現(xiàn)在工作時間的中前段。這是由于在燃燒室燃面退移的過程中,如翼、環(huán)槽的特征結(jié)構(gòu)先逐步加深、擴(kuò)大,而后隨著裝藥燃盡逐漸消失。翼型結(jié)構(gòu)在發(fā)動機(jī)工作中前期已不明顯,而槽型結(jié)構(gòu)依舊顯著。較大的特征結(jié)構(gòu)差別導(dǎo)致了較大的聲頻差別。如圖4所示,翼柱型發(fā)動機(jī)聲基頻在179.2~186.8 Hz變化;第4階聲頻在759.3~789.1 Hz變化,隨工作時間改變最顯著。環(huán)向開槽型發(fā)動機(jī)聲基頻在171.5~183.3 Hz變化;第4階聲頻在712.4~754.4 Hz變化,隨工作時間變化范圍最大。

    圖4 聲頻隨歸一化工作時刻的變化Fig.4 Acoustic frequency variation with normalized operating time

    綜上,特征結(jié)構(gòu)不同將導(dǎo)致主尺寸相同的燃燒室聲頻變化規(guī)律具有明顯差異。相比于環(huán)向開槽型發(fā)動機(jī)翼柱型發(fā)動機(jī)隨工作時間聲頻變換范圍更小。雖然較小的聲頻變化范圍更有利于通過結(jié)構(gòu)設(shè)計(jì)避免結(jié)構(gòu)與聲腔的共振,但是對于學(xué)者們廣泛重視的聲基頻變化范圍兩型發(fā)動機(jī)差別不明顯。

    2.2 聲學(xué)響應(yīng)特性隨燃面退移分析

    最大響應(yīng)幅值實(shí)際上出現(xiàn)在聲模態(tài)頻率處,而由于在諧響應(yīng)分析中,頻率間隔的設(shè)定導(dǎo)致了最大響應(yīng)幅值的頻率出現(xiàn)在實(shí)際聲頻的附近。為方便討論將最大響應(yīng)幅值的頻率直接稱為聲頻。

    如圖5(a)所示,P1在翼柱型發(fā)動機(jī)中最大聲振幅為2 823.5 Pa,其出現(xiàn)在Tn=0.8的模型中,為第2階聲頻下的振蕩。在翼柱型發(fā)動機(jī)中,P1前5階聲頻下的最大振幅依次為1 992.6、2 823.5、2 660.7、1 465.9、1 566.1 Pa。前3階聲頻下的最大響應(yīng)幅值依次出現(xiàn)在Tn=0.5,0.8,0.6時刻模型中;第4階和第5階聲頻下的最大響應(yīng)幅值出現(xiàn)在Tn=1時刻的模型中。將第5階聲頻下的壓力振蕩幅值選做基準(zhǔn),第1~4階聲頻下的振蕩幅值在同一時刻下依次為第5階聲頻下振蕩幅值284%、452%、364%、221%。而且P1點(diǎn)在大部分工作時間內(nèi),第1到3階聲頻下的振蕩比第4、第5階聲頻下的振蕩更為明顯。如圖5(a)所示,P1在翼柱型發(fā)動機(jī)中的主要振蕩頻率為第1~3階聲頻下的振蕩。

    圖5 P1與P2處聲壓時-頻幅值Fig.5 Acoustic pressure amplitude at P1 and P2 in time and frequency domain

    如圖5(b)所示,P1在環(huán)向開槽發(fā)動機(jī)中,最大聲振幅為2 946.2 Pa,其出現(xiàn)在Tn=1的模型中,為第4階聲頻下的振蕩。在環(huán)向開槽型發(fā)動機(jī)中,P1前5階聲頻下的最大振幅依次為1 418.0、854.1、1 430.0、2 946.2、1 749.1 Pa。第1階和第4階聲頻下的最大響應(yīng)幅值出現(xiàn)在Tn=1時刻的模型中;第2階和第3階聲頻下的最大響應(yīng)幅值出現(xiàn)在Tn=0.8時刻的模型中;第5階聲頻下的最大響應(yīng)幅值出現(xiàn)在Tn=0.7時刻的模型中。將第2階聲頻下的壓力振蕩幅值選做基準(zhǔn),第1、3、4、5階聲頻下的振蕩幅值在同一時刻下依次為第2階聲頻下振蕩幅值224%、200%、381%、302%。而且P1點(diǎn)在大部分工作時間內(nèi),第4、第5階聲頻下的振蕩比第1到第3階聲頻下的振蕩更為明顯。如圖5(b)所示,P1在環(huán)向開槽型發(fā)動機(jī)中的主要振蕩頻率為第4、第5階聲頻下的振蕩。

    如圖5(c)所示,P2在翼柱型發(fā)動機(jī)中最大聲振幅為1 378.1 Pa,其出現(xiàn)在Tn=0.9時刻,為第2階聲頻下的振蕩。在翼柱型發(fā)動機(jī)中,P2前5階聲頻下的最大振幅依次為1 243.4、1 378.1、1 094.1、698.3、889.2 Pa。第1、第2和第4階聲頻下的最大響應(yīng)幅值出現(xiàn)在Tn=0.9時刻的模型中;第3階聲頻下的最大響應(yīng)幅值出現(xiàn)在Tn=0.6時刻的模型中;第5階聲頻下的最大響應(yīng)幅值出現(xiàn)在Tn=0.8時刻的模型中。將第5階聲頻下的壓力振蕩幅值選做基準(zhǔn),第1到第4階聲頻下的振蕩幅值在同一時刻下依次為第5階聲頻下振蕩幅值159%、203%、185%、148%。而且P2點(diǎn)在大部分工作時間內(nèi),第1到3階聲頻下的振蕩比第4、第5階聲頻下的振蕩更為明顯。如圖5(c)所示,P2在翼柱型發(fā)動機(jī)中的主要振蕩頻率為第1~3階聲頻下的振蕩。

    如圖5(d)所示,P2在環(huán)向開槽發(fā)動機(jī)中,最大聲振幅為1 515.8 Pa,其出現(xiàn)在Tn=1時刻,為第4階聲頻下的振蕩。在環(huán)向開槽型發(fā)動機(jī)中P2前5階聲頻下的最大振幅依次為1 032.0 Pa、700.2 Pa、914.8 Pa、1 515.8 Pa、976.6 Pa。第1、第4和第5階聲頻下的最大響應(yīng)幅值出現(xiàn)在Tn=1時刻的模型中;第2階和第3階聲頻下的最大響應(yīng)幅值出現(xiàn)在Tn=0.8時刻的模型中。將第2階聲頻下的壓力振蕩幅值選做基準(zhǔn),第1和第3~5階聲頻下的振蕩幅值在同一時刻下依次為第2階聲頻下振蕩幅值189%、147%、277%、198%。而且P2點(diǎn)在大部分工作時間內(nèi),第3~5階聲頻下的振蕩比第1和第2階聲頻下的振蕩更為明顯。如圖5(d)所示,P2在環(huán)向開槽型發(fā)動機(jī)中的主要振蕩頻率為第3~5階聲頻下的振蕩。

    圖5中,2種發(fā)動機(jī)聲頻隨工作時間有不同程度的變化,翼柱型發(fā)動機(jī)第1和第2階聲頻隨工作時間先降低再升高,第3~5階聲頻先升高再下降。環(huán)槽型發(fā)動機(jī)的聲頻隨工作時間按先降低再升高。聲頻隨工作時間的變化規(guī)律與聲模態(tài)分析中得到的結(jié)果一致。2種發(fā)動機(jī)隨著工作時間響應(yīng)幅值在各階聲頻下都有變大趨勢,也就是說2種發(fā)動機(jī)的聲學(xué)穩(wěn)定性都會隨著工作時間逐漸變差。

    由于2種類型發(fā)動機(jī)聲壓響應(yīng)幅值隨時間不斷變化,為了比較2種發(fā)動機(jī)的主要聲振型,將2個發(fā)動機(jī)的聲壓響應(yīng)幅值在工作時間內(nèi)做加權(quán)平均如表3所示。翼柱發(fā)動機(jī)中,發(fā)動機(jī)頭部P1處響應(yīng)幅值最大的3階聲頻,由大幅值到小幅值依次為第2、3、1階聲頻。發(fā)動機(jī)尾部P2處響應(yīng)幅值最大的3階聲頻,按大幅值到小幅值依次為第2、1、3階聲頻。可以得出,翼柱型發(fā)動機(jī)的主要響應(yīng)頻率為前3階聲頻,第2階聲頻處的響應(yīng)能力最強(qiáng)。環(huán)槽型發(fā)動機(jī)中,發(fā)動機(jī)頭部P1與發(fā)動機(jī)尾部P2處響應(yīng)幅值最大的3階聲頻,由大幅值到小幅值依次為第4、5、3階聲頻。環(huán)槽型發(fā)動機(jī)的主要響應(yīng)頻率為第3~5階聲頻,第4階聲頻處的響應(yīng)能力最強(qiáng)。翼柱型2個發(fā)動機(jī)在監(jiān)測點(diǎn)處的響應(yīng)幅值均大于環(huán)槽型發(fā)動機(jī),因此環(huán)槽發(fā)動機(jī)在前3階聲頻下更穩(wěn)定。

    在某一階聲頻響應(yīng)能力強(qiáng),則說明受到該聲頻激勵情況下壓力振幅大。因此,翼柱燃燒室聲場在受到第1~3階聲頻激勵下,振幅放大最為明顯。在翼柱型發(fā)動機(jī)的穩(wěn)定性設(shè)計(jì)中,需要著重考慮避免外界激勵頻率與燃燒室第1~3階聲頻相接近。同理,環(huán)槽型發(fā)動機(jī)應(yīng)著重避免激勵頻率與第3~5階聲頻接近。

    近年來,固體發(fā)動機(jī)中多發(fā)的不穩(wěn)定燃燒為軸向一階不穩(wěn)定。因此,對于聲基頻處的響應(yīng)能力是工程上最應(yīng)關(guān)注的。如表 3,環(huán)槽發(fā)動機(jī)P1處在第1階聲頻的響應(yīng)幅值是翼柱發(fā)動機(jī)的54.8%,環(huán)槽發(fā)動機(jī)P2處在第1階聲頻的響應(yīng)幅值是翼柱發(fā)動機(jī)的75.4%。環(huán)向開槽型發(fā)動機(jī)在工作時間內(nèi)聲基頻下的響應(yīng)幅值顯著低于翼柱型發(fā)動機(jī)。

    表3 時間加權(quán)平均響應(yīng)幅值Table 3 Time-weighted average response amplitude kPa

    尤其在發(fā)動機(jī)工作中期,Tn=0.5時,環(huán)向開槽型發(fā)動機(jī)在P1與P2處的一階聲振幅依次為翼柱型發(fā)動機(jī)的34.5%與48.2%,如圖5。環(huán)槽型發(fā)動機(jī)在一階聲頻下的聲學(xué)穩(wěn)定性優(yōu)勢最為顯著,即便發(fā)生一階聲頻下的壓力振蕩,環(huán)槽發(fā)動機(jī)也比翼柱發(fā)動機(jī)產(chǎn)生的壓力振幅更小。

    綜上,特征結(jié)構(gòu)不同將導(dǎo)致主尺寸相同的燃燒室聲學(xué)響應(yīng)特性具有明顯差異。1 000 Hz以內(nèi),翼柱型發(fā)動機(jī)的主要聲振頻率為第1~3階聲頻;環(huán)向開槽型發(fā)動機(jī)的主要聲振蕩頻率為第4、第5階聲頻。環(huán)槽發(fā)動機(jī)P1處在一階聲頻的時間加權(quán)平均響應(yīng)幅值是翼柱發(fā)動機(jī)的54.8%,環(huán)槽發(fā)動機(jī)P2處在一階聲頻的時間加權(quán)平均響應(yīng)幅值是翼柱發(fā)動機(jī)的75.4%。在大部分工作時間中,環(huán)向開槽型發(fā)動機(jī)在一階聲頻下的聲壓響應(yīng)幅值低于翼柱型裝藥發(fā)動機(jī)。特別是在發(fā)動機(jī)工作中段,環(huán)向開槽型燃燒室的穩(wěn)定性優(yōu)勢最為顯著,環(huán)向開槽型發(fā)動機(jī)P1和P2處的一階聲振幅分別是翼柱型發(fā)動機(jī)的34.5%和48.2%。

    3 結(jié)論

    1)特征結(jié)構(gòu)不同,導(dǎo)致2種等長度發(fā)動機(jī)的聲頻特性有顯著差異。隨著發(fā)動機(jī)工作,燃燒室的各階聲頻隨在一定范圍內(nèi)變化,低階聲頻變化范圍小,高階聲頻變化范圍大。

    2)2種典型發(fā)動機(jī)在相同激勵下,聲學(xué)響應(yīng)隨著工作過程整體呈增大趨勢。

    3)低階聲頻下,尤其在聲基頻下,環(huán)向開槽型發(fā)動機(jī)相比翼柱型發(fā)動機(jī)聲學(xué)穩(wěn)定性更強(qiáng)??梢酝ㄟ^更改特征結(jié)構(gòu),增強(qiáng)發(fā)動機(jī)聲學(xué)穩(wěn)定性。

    猜你喜歡
    聲頻環(huán)槽環(huán)向
    金屬管材力學(xué)性能多角度環(huán)向拉伸測試方法研究
    自承式鋼管跨越結(jié)構(gòu)鞍式支承處管壁環(huán)向彎曲應(yīng)力分析
    汽車零部件(2022年9期)2022-09-30 00:54:20
    一種新型蒸汽聲頻清灰裝置在鍋爐吹灰上的應(yīng)用
    旋轉(zhuǎn)式壓縮機(jī)主軸承減摩技術(shù)研究
    低合金鋼7 倍徑環(huán)槽精銑加工方案
    城市供水管網(wǎng)中鋼筋混凝土岔管受力分析
    初冬游河套
    英國MACAW公司依據(jù)CEPA 2015提出管道環(huán)向應(yīng)力腐蝕開裂預(yù)防處理改進(jìn)方法
    YSZ-50型聲頻振動鉆機(jī)的工程應(yīng)用及技術(shù)改進(jìn)
    亚洲精华国产精华精| 免费观看精品视频网站| 亚洲午夜理论影院| 亚洲中文字幕一区二区三区有码在线看 | 欧美在线黄色| 观看美女的网站| 成人性生交大片免费视频hd| 久久人妻av系列| 日本熟妇午夜| 99re在线观看精品视频| 1024手机看黄色片| 欧美丝袜亚洲另类 | a级毛片a级免费在线| 国产精品精品国产色婷婷| 久久亚洲精品不卡| 黄色 视频免费看| 动漫黄色视频在线观看| av福利片在线观看| 日日摸夜夜添夜夜添小说| 99热精品在线国产| 一级毛片高清免费大全| 欧美日韩国产亚洲二区| 欧美日韩精品网址| 黄色女人牲交| 1024手机看黄色片| 黄色片一级片一级黄色片| 成年免费大片在线观看| 欧美丝袜亚洲另类 | 欧美日韩中文字幕国产精品一区二区三区| 国产精品1区2区在线观看.| 高清在线国产一区| 日本在线视频免费播放| 成人性生交大片免费视频hd| 俺也久久电影网| 久久午夜亚洲精品久久| 久99久视频精品免费| 99精品久久久久人妻精品| 色尼玛亚洲综合影院| 757午夜福利合集在线观看| 亚洲美女黄片视频| 精品国产乱码久久久久久男人| 亚洲成av人片免费观看| 18禁黄网站禁片午夜丰满| 99riav亚洲国产免费| 97超级碰碰碰精品色视频在线观看| 亚洲欧洲精品一区二区精品久久久| 亚洲第一电影网av| xxx96com| 亚洲av片天天在线观看| 免费无遮挡裸体视频| 国产一级毛片七仙女欲春2| 午夜福利18| 欧美日本亚洲视频在线播放| 精华霜和精华液先用哪个| 久久精品影院6| 亚洲av成人不卡在线观看播放网| 欧美又色又爽又黄视频| 亚洲国产看品久久| 国产精品av视频在线免费观看| 19禁男女啪啪无遮挡网站| 俄罗斯特黄特色一大片| 嫩草影院精品99| 夜夜爽天天搞| 精品欧美国产一区二区三| 免费在线观看亚洲国产| 亚洲第一欧美日韩一区二区三区| 黑人欧美特级aaaaaa片| 人人妻人人看人人澡| 白带黄色成豆腐渣| 网址你懂的国产日韩在线| 三级男女做爰猛烈吃奶摸视频| 99久久久亚洲精品蜜臀av| 久久人人精品亚洲av| 久久99热这里只有精品18| 国产三级黄色录像| 18禁美女被吸乳视频| 一级作爱视频免费观看| 高清在线国产一区| 国产精品电影一区二区三区| 国产亚洲精品一区二区www| 1000部很黄的大片| 中文字幕人成人乱码亚洲影| 国内揄拍国产精品人妻在线| 男女午夜视频在线观看| 亚洲专区字幕在线| 麻豆一二三区av精品| 国产高清videossex| 韩国av一区二区三区四区| 国产精品一区二区三区四区免费观看 | 熟女电影av网| 国产精品免费一区二区三区在线| 狂野欧美白嫩少妇大欣赏| 99国产精品99久久久久| 成人国产一区最新在线观看| 少妇裸体淫交视频免费看高清| 色视频www国产| 一本综合久久免费| 亚洲美女视频黄频| 十八禁人妻一区二区| 久久午夜综合久久蜜桃| 久久久久久九九精品二区国产| 久久久久久大精品| 久久久久国产一级毛片高清牌| 欧美不卡视频在线免费观看| 草草在线视频免费看| 日本 欧美在线| 国产aⅴ精品一区二区三区波| 亚洲中文字幕日韩| 舔av片在线| 中出人妻视频一区二区| 国产精品亚洲一级av第二区| 久久香蕉国产精品| e午夜精品久久久久久久| 欧美乱色亚洲激情| av欧美777| 色精品久久人妻99蜜桃| 成人一区二区视频在线观看| 日韩欧美精品v在线| 黄片小视频在线播放| 午夜日韩欧美国产| 999精品在线视频| 少妇的逼水好多| 一级黄色大片毛片| 此物有八面人人有两片| 免费电影在线观看免费观看| 国产精品野战在线观看| 最近最新中文字幕大全电影3| 看黄色毛片网站| 国产99白浆流出| 精品国产三级普通话版| 欧美性猛交黑人性爽| 高清毛片免费观看视频网站| 成年女人毛片免费观看观看9| 国产视频一区二区在线看| 国产精品久久久久久精品电影| 丝袜人妻中文字幕| 一级毛片高清免费大全| 草草在线视频免费看| 亚洲成av人片免费观看| 日本熟妇午夜| 国产精品免费一区二区三区在线| 一个人免费在线观看的高清视频| 国模一区二区三区四区视频 | 99视频精品全部免费 在线 | 午夜福利成人在线免费观看| 亚洲 欧美 日韩 在线 免费| 国产黄a三级三级三级人| 精品无人区乱码1区二区| 国产免费男女视频| 国产av不卡久久| x7x7x7水蜜桃| 超碰成人久久| 天堂动漫精品| e午夜精品久久久久久久| 国产乱人伦免费视频| 亚洲18禁久久av| 欧美乱色亚洲激情| 午夜福利18| 久久精品91无色码中文字幕| 日韩大尺度精品在线看网址| 日韩 欧美 亚洲 中文字幕| 国产精品久久视频播放| 99国产精品99久久久久| 精品国产美女av久久久久小说| 国产亚洲av嫩草精品影院| 国产伦精品一区二区三区四那| 天天躁狠狠躁夜夜躁狠狠躁| 两性夫妻黄色片| 少妇的逼水好多| 国产一区二区在线观看日韩 | 成年免费大片在线观看| 午夜福利在线观看免费完整高清在 | 欧美日韩黄片免| 九九热线精品视视频播放| 国产亚洲欧美98| 久久久久国产精品人妻aⅴ院| 国产私拍福利视频在线观看| 国产亚洲欧美98| 变态另类成人亚洲欧美熟女| 成人特级av手机在线观看| 欧美另类亚洲清纯唯美| 精品一区二区三区视频在线观看免费| 嫩草影院精品99| 久久伊人香网站| 真人做人爱边吃奶动态| 18禁观看日本| 男人舔女人的私密视频| 国产精品国产高清国产av| 日本一本二区三区精品| 国产又黄又爽又无遮挡在线| 一个人看视频在线观看www免费 | 久久国产精品影院| 一a级毛片在线观看| 免费观看人在逋| 狂野欧美白嫩少妇大欣赏| 欧美激情久久久久久爽电影| 国产成年人精品一区二区| 非洲黑人性xxxx精品又粗又长| 国产乱人伦免费视频| 19禁男女啪啪无遮挡网站| 成人亚洲精品av一区二区| 人妻夜夜爽99麻豆av| 51午夜福利影视在线观看| 一本精品99久久精品77| 国产精品久久久av美女十八| 一进一出抽搐动态| 在线a可以看的网站| 19禁男女啪啪无遮挡网站| 成熟少妇高潮喷水视频| 嫩草影视91久久| 亚洲黑人精品在线| 亚洲av片天天在线观看| 亚洲精品美女久久久久99蜜臀| 亚洲精品456在线播放app | 国产欧美日韩一区二区精品| 国产精品98久久久久久宅男小说| 91av网站免费观看| 99热这里只有是精品50| 99re在线观看精品视频| 国产日本99.免费观看| h日本视频在线播放| 19禁男女啪啪无遮挡网站| 日日夜夜操网爽| 欧美日韩精品网址| 听说在线观看完整版免费高清| 国产久久久一区二区三区| 亚洲专区中文字幕在线| 麻豆av在线久日| 黑人操中国人逼视频| 日韩欧美一区二区三区在线观看| 亚洲成a人片在线一区二区| 麻豆成人av在线观看| 色综合站精品国产| 18禁黄网站禁片午夜丰满| 亚洲国产精品合色在线| 丰满的人妻完整版| 99精品久久久久人妻精品| 亚洲午夜理论影院| 中文字幕人成人乱码亚洲影| 一边摸一边抽搐一进一小说| 真实男女啪啪啪动态图| 两个人的视频大全免费| 黄色女人牲交| 欧美色视频一区免费| 国产又黄又爽又无遮挡在线| 欧美三级亚洲精品| 色吧在线观看| 一本一本综合久久| 一级毛片女人18水好多| 伊人久久大香线蕉亚洲五| 久久久久久久午夜电影| 最近在线观看免费完整版| 中文字幕高清在线视频| 欧美+亚洲+日韩+国产| 午夜亚洲福利在线播放| 欧美日韩亚洲国产一区二区在线观看| 欧美日韩综合久久久久久 | 免费高清视频大片| 久9热在线精品视频| 精品久久久久久久人妻蜜臀av| 神马国产精品三级电影在线观看| 制服丝袜大香蕉在线| 国产精品久久电影中文字幕| 午夜福利在线观看免费完整高清在 | 禁无遮挡网站| 国产乱人伦免费视频| 亚洲性夜色夜夜综合| 亚洲在线自拍视频| 久久99热这里只有精品18| 欧美3d第一页| 亚洲精品一区av在线观看| 国产亚洲av嫩草精品影院| 亚洲午夜精品一区,二区,三区| 亚洲国产欧美人成| 最近最新中文字幕大全电影3| 美女午夜性视频免费| 熟女人妻精品中文字幕| 999久久久精品免费观看国产| 免费人成视频x8x8入口观看| 国产成+人综合+亚洲专区| 国产一区二区在线av高清观看| 中文字幕高清在线视频| 色综合亚洲欧美另类图片| 国产97色在线日韩免费| 伊人久久大香线蕉亚洲五| 美女cb高潮喷水在线观看 | 午夜福利视频1000在线观看| 国产乱人伦免费视频| 俄罗斯特黄特色一大片| 国产欧美日韩精品一区二区| 国产av不卡久久| 亚洲天堂国产精品一区在线| avwww免费| 首页视频小说图片口味搜索| 2021天堂中文幕一二区在线观| 1024手机看黄色片| 国内精品美女久久久久久| 久久亚洲真实| 欧美日韩精品网址| 免费搜索国产男女视频| 曰老女人黄片| 男女视频在线观看网站免费| 超碰成人久久| а√天堂www在线а√下载| 成人一区二区视频在线观看| 国产av不卡久久| 久久热在线av| 亚洲av美国av| 国产 一区 欧美 日韩| 日韩 欧美 亚洲 中文字幕| 国产亚洲欧美98| 男人舔女人的私密视频| 狂野欧美白嫩少妇大欣赏| av在线蜜桃| 免费看十八禁软件| 老熟妇乱子伦视频在线观看| 亚洲国产精品成人综合色| 老司机深夜福利视频在线观看| 亚洲中文av在线| 一级毛片女人18水好多| 哪里可以看免费的av片| 欧美一级a爱片免费观看看| 国产主播在线观看一区二区| 校园春色视频在线观看| 国产一区在线观看成人免费| 亚洲精品456在线播放app | 伦理电影免费视频| 婷婷六月久久综合丁香| 99久久成人亚洲精品观看| 国产一区在线观看成人免费| 成人亚洲精品av一区二区| 欧美av亚洲av综合av国产av| 国产亚洲欧美98| 亚洲精品在线观看二区| x7x7x7水蜜桃| 男女午夜视频在线观看| 亚洲一区高清亚洲精品| 丁香欧美五月| 校园春色视频在线观看| 久久亚洲精品不卡| 一卡2卡三卡四卡精品乱码亚洲| 久久人人精品亚洲av| 精品一区二区三区四区五区乱码| 好男人电影高清在线观看| 一本一本综合久久| 一区二区三区高清视频在线| 国产美女午夜福利| 欧美日韩黄片免| 禁无遮挡网站| 1000部很黄的大片| 亚洲成人久久爱视频| 国产高清视频在线播放一区| 一进一出抽搐动态| 亚洲精品在线观看二区| 欧美丝袜亚洲另类 | 久久伊人香网站| 男女下面进入的视频免费午夜| 性欧美人与动物交配| 午夜福利18| 国产视频内射| 国产男靠女视频免费网站| 免费在线观看亚洲国产| 亚洲成人精品中文字幕电影| 国产单亲对白刺激| 一级a爱片免费观看的视频| 狂野欧美激情性xxxx| 少妇人妻一区二区三区视频| 久久精品国产99精品国产亚洲性色| 久久久国产欧美日韩av| 欧美成人免费av一区二区三区| 18禁黄网站禁片免费观看直播| 欧美黄色淫秽网站| 日韩 欧美 亚洲 中文字幕| 欧美日韩亚洲国产一区二区在线观看| 久久亚洲精品不卡| 伦理电影免费视频| 五月玫瑰六月丁香| 国产激情偷乱视频一区二区| 九九热线精品视视频播放| 国产成人精品久久二区二区免费| 最近视频中文字幕2019在线8| 欧美高清成人免费视频www| 噜噜噜噜噜久久久久久91| 国产成人精品久久二区二区免费| 久久中文字幕人妻熟女| 特级一级黄色大片| 搡老妇女老女人老熟妇| 精品福利观看| 国产精华一区二区三区| 黄色 视频免费看| 欧美日韩一级在线毛片| 久久精品人妻少妇| 午夜激情福利司机影院| 日韩成人在线观看一区二区三区| 欧美另类亚洲清纯唯美| 宅男免费午夜| 亚洲欧美日韩高清在线视频| 国产成人福利小说| 国产熟女xx| 91av网站免费观看| 成年版毛片免费区| 亚洲五月婷婷丁香| 午夜福利高清视频| 久久久久性生活片| 日韩 欧美 亚洲 中文字幕| 国产精品久久电影中文字幕| 18禁国产床啪视频网站| 国产精品电影一区二区三区| 成在线人永久免费视频| 久久久久国产精品人妻aⅴ院| 母亲3免费完整高清在线观看| 麻豆成人av在线观看| 禁无遮挡网站| 成人国产一区最新在线观看| www国产在线视频色| 国产精品久久久av美女十八| 在线播放国产精品三级| 久久中文看片网| 国产又黄又爽又无遮挡在线| 国产成人福利小说| 亚洲一区二区三区不卡视频| 国产黄色小视频在线观看| 婷婷亚洲欧美| 又大又爽又粗| 日韩免费av在线播放| 俺也久久电影网| 国产精品国产高清国产av| 亚洲狠狠婷婷综合久久图片| 亚洲欧美日韩高清在线视频| 久久天堂一区二区三区四区| 老司机福利观看| 色在线成人网| 两人在一起打扑克的视频| 日韩免费av在线播放| 久久久色成人| 国产精品一区二区精品视频观看| 少妇人妻一区二区三区视频| 国产亚洲精品综合一区在线观看| 国产欧美日韩一区二区三| 亚洲一区高清亚洲精品| 日本免费一区二区三区高清不卡| 久久天躁狠狠躁夜夜2o2o| 欧美日韩瑟瑟在线播放| 久久久久久大精品| 在线十欧美十亚洲十日本专区| 色精品久久人妻99蜜桃| 99久久久亚洲精品蜜臀av| 在线看三级毛片| 一级a爱片免费观看的视频| 丰满人妻熟妇乱又伦精品不卡| 男女床上黄色一级片免费看| 国产亚洲精品综合一区在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品亚洲一级av第二区| av女优亚洲男人天堂 | 亚洲精品久久国产高清桃花| av在线蜜桃| 中文字幕精品亚洲无线码一区| 美女大奶头视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲,欧美精品.| 真人做人爱边吃奶动态| 亚洲av熟女| 老司机在亚洲福利影院| 99在线视频只有这里精品首页| 日韩精品青青久久久久久| 久久久久国产一级毛片高清牌| 免费av毛片视频| 高潮久久久久久久久久久不卡| 久久久精品欧美日韩精品| 超碰成人久久| a在线观看视频网站| 午夜两性在线视频| 很黄的视频免费| 国产精品香港三级国产av潘金莲| 级片在线观看| 国产精品九九99| 99国产精品一区二区三区| 久久亚洲真实| 成人av一区二区三区在线看| 欧美黄色淫秽网站| 免费电影在线观看免费观看| 亚洲午夜理论影院| 精品久久久久久久末码| 老熟妇仑乱视频hdxx| 国内精品久久久久久久电影| 欧美在线黄色| 淫妇啪啪啪对白视频| 听说在线观看完整版免费高清| 色老头精品视频在线观看| 日韩高清综合在线| 午夜福利在线在线| 国产av一区在线观看免费| АⅤ资源中文在线天堂| 一本综合久久免费| 日韩高清综合在线| 午夜免费成人在线视频| 亚洲欧洲精品一区二区精品久久久| 国产午夜精品久久久久久| 黑人操中国人逼视频| 母亲3免费完整高清在线观看| 我的老师免费观看完整版| 日本撒尿小便嘘嘘汇集6| 少妇丰满av| 看免费av毛片| 日本a在线网址| 精品国产超薄肉色丝袜足j| 97超级碰碰碰精品色视频在线观看| 亚洲欧美日韩无卡精品| 亚洲成a人片在线一区二区| 99久久成人亚洲精品观看| 国产精品亚洲一级av第二区| 欧美zozozo另类| 又黄又爽又免费观看的视频| 一夜夜www| 欧美日韩瑟瑟在线播放| 亚洲精品中文字幕一二三四区| 无限看片的www在线观看| 一个人免费在线观看电影 | 丰满的人妻完整版| 亚洲天堂国产精品一区在线| 国产91精品成人一区二区三区| 国产aⅴ精品一区二区三区波| 美女黄网站色视频| 亚洲人成网站高清观看| 国产精品一区二区三区四区久久| 真人做人爱边吃奶动态| 久久亚洲真实| 国产三级中文精品| av黄色大香蕉| 女生性感内裤真人,穿戴方法视频| 成人国产综合亚洲| 熟女人妻精品中文字幕| 国产精品久久久av美女十八| 国产激情久久老熟女| 在线永久观看黄色视频| 国产精品野战在线观看| 国产精品 国内视频| 国产亚洲精品久久久久久毛片| 两性夫妻黄色片| 此物有八面人人有两片| 级片在线观看| 亚洲九九香蕉| 不卡一级毛片| 校园春色视频在线观看| 成年女人永久免费观看视频| 黄色成人免费大全| 午夜激情福利司机影院| 欧美午夜高清在线| 亚洲成av人片在线播放无| 99riav亚洲国产免费| 久久人人精品亚洲av| 91老司机精品| 国产黄色小视频在线观看| 999精品在线视频| 亚洲av电影不卡..在线观看| 琪琪午夜伦伦电影理论片6080| 国产精品98久久久久久宅男小说| 亚洲人成伊人成综合网2020| 色吧在线观看| 日本熟妇午夜| 亚洲成人久久爱视频| 99久久精品热视频| 久久久久免费精品人妻一区二区| 欧美乱妇无乱码| av国产免费在线观看| 日韩大尺度精品在线看网址| 夜夜爽天天搞| 国产av在哪里看| 久久久久久久久免费视频了| 观看美女的网站| 日韩欧美 国产精品| 午夜成年电影在线免费观看| 久久久国产精品麻豆| 九九在线视频观看精品| 少妇的丰满在线观看| 国产高清视频在线观看网站| 免费在线观看视频国产中文字幕亚洲| av女优亚洲男人天堂 | 日日摸夜夜添夜夜添小说| 国产精品久久久久久人妻精品电影| 日韩免费av在线播放| 久久久久精品国产欧美久久久| 中文字幕熟女人妻在线| 国产极品精品免费视频能看的| 精品久久蜜臀av无| 九色成人免费人妻av| 老司机在亚洲福利影院| xxxwww97欧美| 亚洲av免费在线观看| 婷婷亚洲欧美| 99久久精品国产亚洲精品| 日韩欧美三级三区| 国产精品一区二区三区四区免费观看 | 日本撒尿小便嘘嘘汇集6| 深夜精品福利| 岛国在线观看网站| 免费在线观看成人毛片| 97超视频在线观看视频| 亚洲欧美日韩东京热| 午夜激情欧美在线| 亚洲精品色激情综合| 亚洲av五月六月丁香网| 国产视频内射| 99在线人妻在线中文字幕| 在线a可以看的网站| 亚洲第一欧美日韩一区二区三区| 亚洲av日韩精品久久久久久密| 亚洲 欧美 日韩 在线 免费| 成熟少妇高潮喷水视频| 成在线人永久免费视频| 麻豆久久精品国产亚洲av| 亚洲男人的天堂狠狠| 日韩欧美免费精品| 男人舔女人的私密视频|