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

    泵噴分布式脈動壓力激勵下泵噴艇體耦合系統(tǒng)振動聲輻射

    2021-10-11 09:49:22黃修長蘇智偉師帥康饒志強華宏星
    振動與沖擊 2021年18期

    黃修長, 蘇智偉, 師帥康, 饒志強, 華宏星

    (1.上海交通大學 振動、沖擊、噪聲實驗室,上海 200240;2.上海交通大學 機械系統(tǒng)與振動國家重點實驗室,上海 200240)

    螺旋槳、泵噴等推進器在表面脈動壓力激勵下,一方面會產(chǎn)生結構振動直發(fā)聲以及流噪聲直發(fā)聲;另一方面會通過軸系、定子等傳遞激勵艇體,導致艇體產(chǎn)生強烈的聲輻射。因此,推進器-軸系-艇體的耦合振動聲輻射研究受到廣泛關注[1-3]。

    目前已針對螺旋槳-軸系-艇體耦合系統(tǒng)的振動噪聲特性開展豐富的研究,明確了螺旋槳激勵下螺旋槳-軸系-艇體耦合系統(tǒng)的聲振響應特性[4-5]。但上述分析均采用在螺旋槳槳轂或0.7R處施加集中力(R為螺旋槳半徑),無法考慮分布式載荷的相位不同以及螺旋槳空間模態(tài)特性對響應的影響。為考慮上述問題,部分學者在分析螺旋槳-軸系耦合系統(tǒng)的響應或螺旋槳的流固耦合響應時,采用條帶法將隨機寬帶力映射到螺旋槳結構網(wǎng)格上進行施加[6],或利用面元法和有限元法(finite element method,F(xiàn)EM)耦合將計算得到的脈動壓力施加到結構網(wǎng)格[7-8]。

    泵噴推進器作為一種組合推進器,由于具有較螺旋槳更優(yōu)越的推進及靜音性能而被廣泛應用。泵噴推進器的激勵力更為復雜,其脈動壓力是在導管、定子、轉子吸力面和壓力面上呈復雜空間分布的力,在該空間分布力作用下會激勵起不同模態(tài)產(chǎn)生復雜的響應,其響應特性和響應機理更加復雜。目前針對泵噴推進器的研究主要集中在3個方面:一是水動力性能、表面脈動壓力激勵特性、流場特性等的研究[9-10];二是泵噴參數(shù)(如間隙、導管幾何參數(shù)和定轉子幾何參數(shù)等)對水動力學性能、脈動壓力等的影響[11-12];三是推進器脈動壓力激勵下艇體的聲輻射特性以及新型結構(如鋸齒后緣導管)的降噪效果[13]。但對于泵噴激勵下泵噴-軸系-艇體耦合系統(tǒng)振動聲輻射的研究鮮見公開報道。對于泵噴推進器激勵下的艇體振動聲輻射響應,泵噴推進器同時通過推進軸系激勵艇體以及通過導管-定子激勵艇體,這兩種途徑的傳遞特性、傳遞幅值相對大小等仍不明確。

    本論文考慮suboff艇后泵噴推進器的激勵力為轉子、導管和定子上的分布式脈動壓力分布,通過插值的方法把流體脈動壓力映射到推進器結構濕表面,建立了泵噴-軸系-艇體的耦合動力學模型,獲得了其結構振動聲輻射響應,并分析了通過軸系激勵艇體、通過導管-定子激勵艇體這兩條傳遞途徑的振動聲輻射響應特性。

    1 理論基礎

    1.1 suboff艇后泵噴表面脈動壓力計算

    不考慮推進器內(nèi)部流場空化影響,可假設流體為不可壓縮黏性單相流,基于雷諾平均納維-斯托克斯(reynolds-averaged navier-stokes,RANS)的控制方程可寫為[14]

    式中:i,j=1,2,3;ρ為流體密度;xi和xj分別為笛卡爾坐標系下3個方向的坐標分量;P為雷諾平均壓力;Uj和u′i分別為平均速度和脈動速度;μ為流體動力黏度系數(shù);為雷諾應力項,為了使式(2)封閉,引入渦流黏度vt,雷諾應力可寫為

    式中:k為雷諾時間平均下湍流動能;Cμ=0.09。采用剪切應力輸運(shear stress transfer,SST)k-ω湍流模型進行模擬。

    1.2 泵噴-軸系-艇體耦合系統(tǒng)振動聲輻射模型

    采用耦合有限元法開展泵噴-軸系-艇體耦合系統(tǒng)的動力學建模

    式中:M和K分別為泵噴-軸系-艇體耦合系統(tǒng)的質量和剛度矩陣;C為阻尼矩陣,C=αM+βK;X為結構的位移向量;p為流場中的聲壓向量;F為泵噴的分布式脈動激勵力;Ma和Ka分別為流體的質量矩陣和剛度矩陣;H為流固耦合矩陣;∏建立了結構濕表面上的節(jié)點法向位移與結構振動位移向量X之間的關系。Ma,Ka和H的表達式分別為

    式中:N和Nψ分別為結構和流體的插值形函數(shù);ρa為流體密度;c為流體中波速。

    由結構濕表面的法向振速viα,利用直接邊界元法(boundary element method,BEM)求得濕表面上的聲壓piα以及場點聲壓pj[15]

    式中:J(ξ)為Jacobian矩陣;Ai為第i個單元的面積;(ξ)為第j個場點與泵噴-軸系-艇體耦合系統(tǒng)結構濕表面上第i個單元的第α個節(jié)點之間的距離;N(ξ)為邊界元中聲壓的插值形函數(shù)。

    1.3 表面分布式脈動壓力和結構濕表面節(jié)點力映射

    將計算流體動力學(computational fluid dynamics,CFD)模型中轉子、定子、導管表面單元的壓力以及泵噴結構濕表面的節(jié)點導出,通過徑向基函數(shù)插值將流體單元節(jié)點的壓力映射到泵噴結構濕表面節(jié)點上。假設空間存在一組點群xj及點群上的函數(shù)值fj(j=1,2,…,N),可構建各點之間的距離函數(shù)?(‖x-xj‖)與函數(shù)值之間的關系

    式中:?(‖x-xj‖)為徑向基函數(shù),它是插值點與目標點之間距離的函數(shù);λj為徑向基函數(shù)的系數(shù)。通過xj可建立點群之間的λj系數(shù)求解線性方程組??臻g任意點上的函數(shù)值f(x)則可通過式(8)求解。具體計算時首先找到泵噴結構濕表面單元中心附近的一組流體節(jié)點及壓力值,并求解徑向基函數(shù)插值系數(shù),然后將泵噴結構濕表面單元中心代入式(8)中求解泵噴結構濕表面單元上的壓力分布。在求得泵噴結構濕表面單元的壓力分布后,將壓力乘以面積得到面單元上的集中力,將集中力按節(jié)點數(shù)平均分配到各個節(jié)點上,即可得到每個節(jié)點上的集中力分布。在每個時刻進行徑向基函數(shù)插值即可得到時域中的結構濕表面節(jié)點的三向集中力。利用傅里葉變換可得到頻率域結構濕表面節(jié)點的三向集中力。

    2 數(shù)值模型

    研究對象為suboff艇后泵噴推進器模型,如圖1所示。泵噴-軸系-suboff艇參數(shù),如表1所示。殼體厚度為6 mm,并采用T型加強筋進行加強,縱向方向有4個艙壁;泵噴推進器通過定子支撐在殼體艉部。泵噴-軸系模型見圖1,推進軸系的長度為0.89 m、直徑為0.022 m,通過艉后軸承支撐在殼體艉部以及通過推力軸承支撐在最后一個隔艙壁的中心。艉后軸承的橫向和垂向剛度分別為4.5×109N/m,1.7×108N/m;推力軸承的縱向、橫向和垂向剛度分別為6.6×108N/m,2.3×108N/m和6.7×108N/m。殼體、加強筋、軸系、泵噴定子和導管的材料為鋼,其密度、彈性模量和泊松比分別為7 850 kg/m3,2.1×1011Pa,0.3。 轉子的材料為青銅,其密度、彈性模量和泊松比分別為8 790 kg/m3,1.1×1011Pa,0.3。 結構阻尼均取為0.001。

    圖1 泵噴-軸系-艇體耦合系統(tǒng)Fig.1 Pump-jet-shafting-suboff system

    表1 泵噴-軸系-suboff艇參數(shù)Tab.1 Parameters of pump-jet-shafting-suboff system

    2.1 CFD模型和計算驗證

    創(chuàng)建泵噴-軸系-suboff的CFD計算區(qū)域,如圖2(a)所示,將計算域劃分成靜止域和旋轉域,靜止域包含定子、艇體計算域和導管計算域,將轉子所在旋轉域設置為滑移網(wǎng)格,兩個計算域之間的交界面設置為interface邊界條件,進行數(shù)據(jù)傳遞。計算中遠前方入口設置為速度入口條件,遠后方出口設置為壓力出口條件,外圍遠場邊界與入口一樣為速度入口條件,轉子、定子及導管等固體表面采用無滑移固體壁面條件。部分進行流體計算的面網(wǎng)格,如圖2(b)所示。結合網(wǎng)格無關性分析結果,選擇合適的網(wǎng)格密度進行流域網(wǎng)格劃分。計算域網(wǎng)格共計2 760萬,其中:艇體域網(wǎng)格為1 000萬;泵噴定子域網(wǎng)格數(shù)量為660萬;泵噴轉子域網(wǎng)格為1 100萬。艇體表面和泵噴推進器均劃分精細邊界層網(wǎng)格,目標Y+值為1。

    圖2 suboff艇后泵噴的計算域和CFD網(wǎng)格Fig.2 Computational domain and CFD model of the behind-hull suboff

    計算工況為進速9.42 m/s,轉子轉速1 824 r/min,進行非定常計算,設置時間步長為4.568 7×10-5s,保證轉子每個時間步旋轉0.5°,計算的最高頻率為10 944 Hz,采樣頻率為21 888 Hz。數(shù)值計算的收斂性殘差滿足:連續(xù)性殘差收斂在10-5以下,3個方向速度殘差均在10-6以下,湍動能k殘差低于10-7。

    CFD數(shù)值求解采用商業(yè)軟件ANSYS Fluent進行,采用SST k-ω湍流模型,壓力、速度耦合求解采用SIMPLEC(semi-implicit method for pressure-linked equations consistent)方法,其他項均采用二階迎風格式。

    采用對suboff艇和泵噴單獨進行CFD計算驗證的方法對建立的模型進行驗證。對suboff艇體進行驗證時,保留圖2中的suboff網(wǎng)格,流域和網(wǎng)格密度與圖2相同。采用上述網(wǎng)格模型對進速3.051 m/s的工況進行計算,得到全附體suboff艇阻力(103.2 N)與實測阻力(102.3 N)誤差為0.9%[16],同時得到的槳盤面處的軸向伴流結果與試驗結果相似,如圖3所示。對泵噴進行驗證時,采用與圖2相同的流域、網(wǎng)格密度(其中:Grid2與圖2中網(wǎng)格數(shù)量相同;Grid3將轉子域中網(wǎng)格增加40%)和計算設置,對一五葉前置定子七葉轉子泵噴進行計算,如圖4所示。將計算結果與空泡水筒中的試驗結果進行對比,采用Grid2推力最大誤差為3.59%,最小誤差為0.23%,扭矩最大誤差為5.5%,效率的最大誤差為5.17%,精度較高,增加網(wǎng)格數(shù)量后(Grid3)增加了計算成本,精度提高不明顯,故采用Grid2滿足分析要求。

    圖3 槳盤面處軸向伴流對比Fig.3 Comparison of axial velocity in the rotor plane

    圖4 泵噴水動力性能計算模型和試驗驗證Fig.4 CFD model and experimental validation for propulsion performance of the pump-jet

    2.2 結構動力學模型

    在ABAQUS軟件中建立泵噴-軸系-艇體耦合系統(tǒng)的結構有限元模型,如圖5所示。殼體采用33 861個S4R單元進行模擬,加強筋采用4 354個B31單元進行模擬,艙壁采用4 362個S4R單元進行模擬。定子、轉子和導管采用22 641個C3D8實體單元進行模擬。推進軸系采用91個B31梁單元進行模擬。為考慮流固耦合,在泵噴-軸系-艇體耦合系統(tǒng)的結構有限元模型外建立圓柱形的水體,在結構有限元模型的濕表面和水體模型的表面之間建立壓力和速度耦合條件。其中圓柱形水體的長度和直徑分別為6 600 mm,2 500 mm(直徑為6倍導管最大直徑)。水體采用2 601 967個AC3D8聲學單元進行模擬。水的密度和體積模量分別為1 000 kg/m3和2×109Pa。在水體外表面施加無反射聲學邊界條件。將泵噴導管、定子和轉子濕表面以及艇體濕表面的速度導出到LMS VirtualLab中,利用直接邊界元法進行聲輻射的計算。結構單元和聲輻射單元的網(wǎng)格大小為5 mm,根據(jù)拇指法則,結果在600 Hz以內(nèi)具有可信度[17]。

    圖5 泵噴-軸系-艇體耦合系統(tǒng)結構振動聲輻射模型Fig.5 Vibro-acoustic model of coupled pump-jet-shafting-suboff system

    3 仿真結果與分析

    3.1 表面脈動壓力

    取出某個轉子葉片表面0.7R,0.8R,0.9R和相應轉子區(qū)域對應的導管內(nèi)壁面某點的脈動壓力,其頻譜結果如圖6所示。由圖6可知,葉片的脈動壓力以軸頻30.4 Hz及其倍頻為主,導管上的脈動壓力以葉頻273.7 Hz(9倍軸頻)和2倍葉頻547.4 Hz為主??傮w而言,導管上的脈動壓力幅值比轉子葉片上監(jiān)測點的脈動壓力幅值大1~2個數(shù)量級。

    圖6 轉子不同半徑處和導管的表面脈動壓力頻譜Fig.6 Spectrum of pressure on the different radius of rotor surface and the duct

    轉子縱向、橫向和垂向合力的時程曲線和頻譜,如圖7所示。由圖7可知,由于不同轉子葉片上的脈動壓力相位不同,疊加后各個方向的脈動力在葉頻及其倍頻處出現(xiàn)較大的峰值,在軸頻及其倍頻處也表現(xiàn)為峰值,縱向脈動力的幅值大于垂向和橫向脈動力。隨著頻率的增加,脈動力的幅值呈衰減趨勢。

    圖7 轉子三向脈動力的時程曲線和頻譜Fig.7 Time-history,spectrum of pulsation force in three directions for the rotor

    圖8給出了CFD計算得到的表面脈動壓力結果和插值后導入到ABAQUS中結構濕表面脈動壓力結果的對比,給出了轉子壓力面和導管內(nèi)壁面的結果。由圖8可知,采用徑向基函數(shù)插值后,脈動壓力的分布和幅值均能夠保證足夠的精度。

    圖8 分布式脈動壓力的CFD結果和FEM結果比較Fig.8 Comparison of fluctuation pressure distribution on surface of propeller,duct between CFD and mapped result in FEM

    3.2 模態(tài)結果與分析

    對泵噴-軸系-艇體耦合系統(tǒng)動力學模型開展模態(tài)分析,得到典型模態(tài)結果如圖9、表2所示。

    表2 耦合系統(tǒng)模態(tài)頻率列表 Tab.2 List of modal frequencies of the coupled system Hz

    圖9 泵噴-軸系-艇體耦合系統(tǒng)典型模態(tài)Fig.9 Typical modal shapes of the coupled pump-jet-shafting-suboff system

    下面計算泵噴-軸系-艇體耦合系統(tǒng)的振動聲輻射響應??紤]4種載荷施加方法:一是將導管、轉子和定子上的脈動壓力施加到結構網(wǎng)格上(‘Pumpjet’);二是僅將轉子上的脈動壓力施加到結構網(wǎng)格上(‘Ro-tor’);三是對轉子葉片上的脈動壓力進行求和,將每個葉片上的三向集中力施加到轉子葉片0.7R處(‘0.7R’);四是對整個轉子的脈動壓力求和,將整個轉子的三向集中力施加到槳轂中心(‘Hub’)。

    3.3 通過軸承的傳遞力結果與分析

    對通過各個軸承的傳遞力進行分析。推力軸承和艉后軸承在不同載荷施加方式下的多向傳遞力,如圖10所示。對于推力軸承,通過其傳遞的三向力對于‘Pumpjet’和‘Rotor’基本重合,‘0.7R’將使結果偏大;‘Hub’將使212.9 Hz以前的結果偏小。推力軸承縱向傳遞力的最大響應峰值為軸頻30.4 Hz。5個顯著的峰值分別為209.9 Hz(軸系一階縱振),212.9 Hz(7倍軸頻),243.3 Hz(8倍軸頻)和273.7 Hz(9倍軸頻,即葉頻)以及547.5 Hz(18倍軸頻,即2倍葉頻)。軸系一階縱振的峰值頻率顯著和相鄰的峰值具有可比性。一些小的響應峰值出現(xiàn)在334.6 Hz(艇體一階縱振),360.4 Hz,391.3 Hz(軸系二階縱振),486.6 Hz和577.9 Hz。對于推力軸承垂向傳遞力,‘0.7R’和‘Hub’施加方式在葉頻273.7 Hz處出現(xiàn)峰值,但是在‘Pumpjet’和‘Rotor’情況下并未出現(xiàn)。主要的峰值頻率為30.4 Hz,118.6 Hz(軸系一階彎曲),121.7 Hz(4倍軸頻),243.3 Hz(8倍軸頻),365.0 Hz(12倍軸頻),415.2 Hz(在泵噴處出現(xiàn)最大變形的軸系彎曲模態(tài)),517.0 Hz(17倍軸頻)和547.5 Hz(18倍軸頻)。推力軸承橫向傳遞力和垂向傳遞力具有類似的規(guī)律,主要的峰值頻率為30.4 Hz,118.6 Hz,121.7 Hz,152.1 Hz(5倍軸頻),243.3 Hz,304.1 Hz(10倍軸頻),365.0 Hz,416.7 Hz(在泵噴處出現(xiàn)最大變形的軸系彎曲模態(tài))和517.0 Hz。除了軸頻,4個最大的響應峰值出現(xiàn)在121.7 Hz,243.3 Hz,365.0 Hz和415.2 Hz;小的峰值出現(xiàn)在273.7 Hz,334.6 Hz(11倍軸頻)。推力軸承縱向傳遞力的幅值比橫向和垂向傳遞力大一個數(shù)量級。

    圖10 推力軸承和艉后軸承的力傳遞Fig.10 Force transmission by thrust bearing and stern tube bearing

    對于艉后軸承傳遞力,可得到類似結論。垂向傳遞力的顯著峰值出現(xiàn)在30.4 Hz,118.6 Hz,121.7 Hz,152.1 Hz,182.5 Hz(5倍軸頻),212.9 Hz(7倍軸頻),243.3 Hz,304.1 Hz,365.0 Hz,416.7 Hz和547.5 Hz。橫向傳遞力的顯著峰值出現(xiàn)在30.4 Hz,118.6 Hz,121.7 Hz,152.1 Hz,182.9 Hz,212.9 Hz,243.3 Hz,304.1 Hz,365.0 Hz,416.7 Hz和517.0 Hz。在艉后軸承處的垂向和橫向傳遞力比推力軸承處的垂向和橫向傳遞力大一個數(shù)量級。

    在上述峰值頻率中,大部分為泵噴推進器脈動壓力的軸頻及其倍頻,以及結構的模態(tài),如艇體一階縱振模態(tài)、軸系的一階和二階縱振模態(tài),但是‘Hub’無法激勵起豐富的脈動壓力峰值頻率和模態(tài)頻率。

    3.4 水下聲輻射結果與分析

    泵噴-軸系-艇體耦合系統(tǒng)的水下輻射聲功率(sound radiation power,SRP)和均方振速(mean-square velocity,MSV)級,如圖11所示。不同的載荷施加方式下輻射聲功率和均方振速的變化規(guī)律和傳遞力類似??梢娹D子上分布式脈動壓力引起的振動聲輻射將占主導作用,主要的峰值頻率出現(xiàn)在30.4 Hz,91.3 Hz(3倍軸頻),121.7 Hz,152.1 Hz,182.5 Hz(6倍軸頻),209.9 Hz,212.9 Hz,243.3 Hz,273.7 Hz,304.1 Hz,334.6 Hz(艇體一階縱振),365.0 Hz,395.4 Hz(13倍軸頻),416.7 Hz,456.2 Hz(15倍軸頻),486.6 Hz(16倍軸頻),517.0 Hz,547.5 Hz和577.9 Hz。 一些小的峰值出現(xiàn)在209.9 Hz,235.7 Hz和252.4 Hz。并且在艇體一階縱振模態(tài)334.6 Hz處的響應比相鄰峰值均大。傳遞力中在軸系一階彎曲模態(tài)118.6 Hz出現(xiàn)的峰值在輻射聲功率和均方振速中未得到明顯體現(xiàn),這是由于軸系彎曲模態(tài)通過艉后軸承傳遞的力以及通過推力軸承傳遞到艙壁的力都難以有效激勵起艇體的振動聲輻射。推力軸承傳遞力在軸系一階縱振模態(tài)209.9 Hz出現(xiàn)的峰值在輻射聲功率和均方振速中仍具有較大的峰值。為進一步研究不同方向脈動壓力激勵起的艇體聲輻射,取出泵噴的縱向脈動壓力施加到轉子、定子和泵噴表面,將得到的輻射聲功率和同時在泵噴表面施加三向脈動壓力的結果進行比較,如圖11(c)所示,可見在2倍軸頻以上,聲輻射主要由縱向脈動壓力主導。

    圖11 輻射聲功率和均方振速Fig.11 Sound radiation power level and mean-square velocity level

    給出4個典型頻率30.4 Hz(軸頻),243.3 Hz(8倍軸頻),273.7 Hz(葉頻),334.6 Hz(11倍軸頻)的聲輻射指向性,如圖12所示。場點位于水平面內(nèi),是以原點為中心直徑50 m的圓。軸頻處指向性為橢圓,首尾方向的輻射噪聲較大;8倍軸頻處為蝴蝶形,首尾方向和正橫方向的輻射聲壓較小。除葉頻外,其他3個頻率處‘Pumpjet’和‘Rotor’基本重合。葉頻處各種載荷施加方式的指向性差異較大,這是由于聲輻射同時包含了艇體、泵噴轉子、泵噴定子和導管的輻射聲壓,導管本身脈動壓力在葉頻處較大,也產(chǎn)生直接聲輻射,并且和轉子、定子等輻射聲壓疊加,而其他施加方式下導管的振動較小。

    圖12 4個頻率處的指向性圖Fig.12 Sound radiation directivity pattern at four frequencies

    為衡量導管和定子上的脈動壓力對水下聲輻射的貢獻,僅將導管和定子上的脈動壓力施加到結構網(wǎng)格上,計算輻射聲功率級和均方振速級,結果如圖13所示。由圖13可知,導管和定子上的脈動壓力對艇體水下聲輻射的主要貢獻體現(xiàn)在葉頻和2倍葉頻處,這是由于導管上的脈動壓力在這兩個頻率處為顯著峰值。另一方面,結合圖9的模態(tài)結果可知,與這兩個激勵頻率接近的結構模態(tài)也被激發(fā),出現(xiàn)峰值。分別在整個泵噴推進器表面、只有轉子表面、只有導管-定子表面施加脈動激勵力,計算得到各自的輻射聲功率和均方振速,并進行比較。可知,在單獨導管-定子脈動壓力激勵下,存在一些突出線譜峰值,量級和轉子脈動壓力激勵下的輻射聲功率和均方振速量級可比擬,但它們對整體聲輻射的貢獻可以忽略。在兩者同時激勵時,均方振速和輻射聲功率級并沒有顯著增大。這是由于在導管-定子和轉子脈動壓力同時激勵下,相位差會導致相互抵消。為驗證抵消效應,取出艇體和導管上某監(jiān)測點的法向速度響應,如圖14所示。由圖14可知,在整個泵噴推進器脈動壓力激勵時,法向速度會比導管-定子、轉子分別施加脈動壓力時小,這表明,同時施加脈動壓力會由于相位差相互抵消,不同位置處的相互抵消情況會有所差別。

    圖13 導管-定子、轉子分別施加脈動壓力時響應Fig.13 SRPand MSV under distributed pressure applied on duct-stator,rotor separately

    圖14 導管-定子、轉子分別施加脈動壓力時法向振速響應Fig.14 Normal velocity under distributed pressure applied on duct-stator,rotor separately

    4 結 論

    開展了suboff艇后泵噴推進器脈動壓力激勵下泵噴-軸系-艇體耦合系統(tǒng)的振動聲輻射響應計算和分析。利用CFD模型獲得了艇體泵噴推進器的脈動壓力。采用耦合有限元法和邊界元法獲得了分布式脈動壓力激勵下泵噴-軸系-艇體耦合系統(tǒng)的傳遞力、聲輻射響應。得到結論如下:

    (1)泵噴-軸系-艇體耦合系統(tǒng)的振動聲輻射以轉子脈動壓力激勵下的響應占主導,而其中2倍軸頻以上以縱向脈動壓力的貢獻占主導。

    (2)泵噴-軸系-艇體耦合系統(tǒng)振動聲輻射響應中,顯著的峰值包含泵噴脈動壓力的軸頻及其倍頻成分(以軸頻、葉頻最大),也包含軸系一階縱振和二階縱振、艇體縱向振動等結構模態(tài)。

    (3)導管的脈動壓力以葉頻和2倍葉頻為主,導管-定子脈動壓力激勵下振動聲輻射響應峰值出現(xiàn)在上述頻率處;且響應量級和轉子表面脈動壓力激勵處于同一量級;但導管-定子、轉子表面同時施加脈動壓力時,由于相位相互抵消,均方振速和輻射聲功率級并沒有顯著增大。

    真人做人爱边吃奶动态| 成年免费大片在线观看| 99热这里只有是精品50| 欧美日韩乱码在线| 少妇丰满av| 观看免费一级毛片| 欧美性感艳星| 欧美黄色片欧美黄色片| 日韩欧美在线乱码| x7x7x7水蜜桃| 99久久久亚洲精品蜜臀av| 狠狠狠狠99中文字幕| 日韩大尺度精品在线看网址| 亚洲av免费在线观看| 国产精品久久久久久亚洲av鲁大| 91九色精品人成在线观看| 欧美日韩福利视频一区二区| 久久婷婷人人爽人人干人人爱| 精品久久久久久久久久免费视频| 亚洲乱码一区二区免费版| 亚洲精品一卡2卡三卡4卡5卡| 99久久精品国产亚洲精品| 别揉我奶头~嗯~啊~动态视频| 欧美一级a爱片免费观看看| 无人区码免费观看不卡| 我要搜黄色片| 香蕉久久夜色| 99久久精品国产亚洲精品| 极品教师在线免费播放| 在线观看av片永久免费下载| 国产久久久一区二区三区| 99精品久久久久人妻精品| 久久久久久久精品吃奶| 成年女人永久免费观看视频| 亚洲国产精品999在线| 日本黄大片高清| 免费看光身美女| 亚洲精品影视一区二区三区av| 91九色精品人成在线观看| 一级黄片播放器| 天天添夜夜摸| 97碰自拍视频| 成人欧美大片| 真人一进一出gif抽搐免费| 久久久国产精品麻豆| 午夜福利18| 亚洲午夜理论影院| 国产麻豆成人av免费视频| 少妇裸体淫交视频免费看高清| 91字幕亚洲| 免费在线观看日本一区| svipshipincom国产片| www国产在线视频色| 国产亚洲欧美98| 毛片女人毛片| 女人十人毛片免费观看3o分钟| 丁香欧美五月| 九色国产91popny在线| 18美女黄网站色大片免费观看| 亚洲成人免费电影在线观看| 高清在线国产一区| 嫁个100分男人电影在线观看| 欧美午夜高清在线| 又粗又爽又猛毛片免费看| 欧美+日韩+精品| 制服人妻中文乱码| 国产单亲对白刺激| 国产精品嫩草影院av在线观看 | 国产午夜精品久久久久久一区二区三区 | 黄色成人免费大全| 又黄又粗又硬又大视频| 免费在线观看日本一区| 国产精品久久久久久久久免 | 在线十欧美十亚洲十日本专区| 欧美成狂野欧美在线观看| 18禁国产床啪视频网站| 嫩草影视91久久| 国内少妇人妻偷人精品xxx网站| 久久国产精品影院| 色老头精品视频在线观看| 成人精品一区二区免费| 熟女少妇亚洲综合色aaa.| 久久精品91无色码中文字幕| 久久久精品欧美日韩精品| 色噜噜av男人的天堂激情| 伊人久久大香线蕉亚洲五| 脱女人内裤的视频| 内地一区二区视频在线| 99久久综合精品五月天人人| 99久国产av精品| 色综合欧美亚洲国产小说| 亚洲第一欧美日韩一区二区三区| tocl精华| 男女视频在线观看网站免费| a级毛片a级免费在线| 国产乱人视频| 天天一区二区日本电影三级| 人人妻人人看人人澡| 一级黄片播放器| 狠狠狠狠99中文字幕| 变态另类丝袜制服| av中文乱码字幕在线| 久久精品国产亚洲av涩爱 | 老司机在亚洲福利影院| 蜜桃久久精品国产亚洲av| 性色avwww在线观看| 久久久久久国产a免费观看| 国产不卡一卡二| 日本黄大片高清| 国产视频一区二区在线看| 亚洲中文字幕一区二区三区有码在线看| 无限看片的www在线观看| 日韩国内少妇激情av| 国产精品国产高清国产av| 国产毛片a区久久久久| 又紧又爽又黄一区二区| 一个人看视频在线观看www免费 | 国产精品98久久久久久宅男小说| 日本 欧美在线| 精品免费久久久久久久清纯| 午夜福利欧美成人| 日韩欧美国产一区二区入口| 日韩中文字幕欧美一区二区| 国产精品1区2区在线观看.| 一个人看的www免费观看视频| 天堂av国产一区二区熟女人妻| 成人18禁在线播放| 老司机福利观看| 久久久久性生活片| 黄片大片在线免费观看| 十八禁网站免费在线| 欧美黄色片欧美黄色片| 男人舔女人下体高潮全视频| 免费看a级黄色片| 亚洲精品成人久久久久久| 国产蜜桃级精品一区二区三区| 国产午夜精品久久久久久一区二区三区 | 中文资源天堂在线| 中文字幕人妻熟人妻熟丝袜美 | 一区二区三区激情视频| 欧美绝顶高潮抽搐喷水| 窝窝影院91人妻| 久久久久久久久中文| 亚洲av中文字字幕乱码综合| 国产色婷婷99| 亚洲中文日韩欧美视频| 伊人久久精品亚洲午夜| 国产精品一区二区三区四区免费观看 | 成人特级黄色片久久久久久久| 国产成人a区在线观看| 国内久久婷婷六月综合欲色啪| 性色avwww在线观看| or卡值多少钱| 高清日韩中文字幕在线| 国产精品av视频在线免费观看| 亚洲人成网站在线播放欧美日韩| av专区在线播放| 99久国产av精品| 欧美色欧美亚洲另类二区| 69av精品久久久久久| 少妇裸体淫交视频免费看高清| 欧美黄色片欧美黄色片| 国产真人三级小视频在线观看| 中文字幕人成人乱码亚洲影| 欧美bdsm另类| 88av欧美| 亚洲,欧美精品.| 琪琪午夜伦伦电影理论片6080| 深夜精品福利| 精华霜和精华液先用哪个| 国产99白浆流出| 亚洲专区国产一区二区| 午夜日韩欧美国产| 国产成人av激情在线播放| 免费看日本二区| 国产美女午夜福利| 18禁裸乳无遮挡免费网站照片| 亚洲av熟女| 欧美xxxx黑人xx丫x性爽| 亚洲五月婷婷丁香| 亚洲专区国产一区二区| 成人亚洲精品av一区二区| 欧美大码av| 免费一级毛片在线播放高清视频| 国产成人a区在线观看| 成人精品一区二区免费| 人妻丰满熟妇av一区二区三区| 午夜久久久久精精品| 看片在线看免费视频| 欧美日本视频| 看片在线看免费视频| 国产高清videossex| 久久国产精品人妻蜜桃| 免费av观看视频| 变态另类成人亚洲欧美熟女| 大型黄色视频在线免费观看| 国产一区在线观看成人免费| 国产成年人精品一区二区| 国产精品香港三级国产av潘金莲| 一区二区三区高清视频在线| 99国产精品一区二区三区| 欧美日本视频| 日韩国内少妇激情av| 天天一区二区日本电影三级| 婷婷精品国产亚洲av| av在线蜜桃| 国产三级在线视频| 热99在线观看视频| 亚洲激情在线av| 久久久国产精品麻豆| 我要搜黄色片| 久久午夜亚洲精品久久| 桃色一区二区三区在线观看| 一级毛片女人18水好多| 免费在线观看成人毛片| 一进一出抽搐动态| 男人和女人高潮做爰伦理| 免费人成视频x8x8入口观看| 男女视频在线观看网站免费| а√天堂www在线а√下载| 一个人免费在线观看的高清视频| 亚洲国产精品sss在线观看| 精华霜和精华液先用哪个| 2021天堂中文幕一二区在线观| 我的老师免费观看完整版| 91久久精品电影网| 麻豆久久精品国产亚洲av| 国产一级毛片七仙女欲春2| 国产精品亚洲一级av第二区| 老司机在亚洲福利影院| 老熟妇仑乱视频hdxx| 一级a爱片免费观看的视频| 亚洲激情在线av| 51午夜福利影视在线观看| 欧美成狂野欧美在线观看| 色av中文字幕| 国产探花在线观看一区二区| 在线观看免费视频日本深夜| 欧美最新免费一区二区三区 | 香蕉av资源在线| 丰满的人妻完整版| 欧美黑人欧美精品刺激| 国产欧美日韩一区二区三| 亚洲精品一区av在线观看| 一区二区三区激情视频| 亚洲av成人不卡在线观看播放网| 成年女人看的毛片在线观看| 少妇熟女aⅴ在线视频| 欧美激情久久久久久爽电影| 99国产精品一区二区三区| 国产一级毛片七仙女欲春2| 精品久久久久久,| 亚洲av免费在线观看| 最近最新中文字幕大全免费视频| 国模一区二区三区四区视频| 麻豆成人午夜福利视频| 人妻久久中文字幕网| 天堂影院成人在线观看| 脱女人内裤的视频| 国产毛片a区久久久久| 一个人看的www免费观看视频| 三级国产精品欧美在线观看| 一区福利在线观看| 听说在线观看完整版免费高清| 久久久久久久精品吃奶| 无人区码免费观看不卡| 男女那种视频在线观看| 韩国av一区二区三区四区| 少妇熟女aⅴ在线视频| 免费在线观看影片大全网站| 欧美日本亚洲视频在线播放| 岛国视频午夜一区免费看| 搡老妇女老女人老熟妇| 免费电影在线观看免费观看| 久久久精品欧美日韩精品| 嫁个100分男人电影在线观看| 亚洲狠狠婷婷综合久久图片| 一个人观看的视频www高清免费观看| 久久久久久久精品吃奶| 国产69精品久久久久777片| 国产精品一区二区免费欧美| 伊人久久精品亚洲午夜| 亚洲 国产 在线| 国产精品香港三级国产av潘金莲| 99久久综合精品五月天人人| 亚洲不卡免费看| 日本一本二区三区精品| 亚洲av成人精品一区久久| netflix在线观看网站| 亚洲中文日韩欧美视频| 亚洲自拍偷在线| 欧美成人性av电影在线观看| 欧美成人一区二区免费高清观看| 亚洲aⅴ乱码一区二区在线播放| 久久中文看片网| www.999成人在线观看| 国产精品自产拍在线观看55亚洲| 亚洲激情在线av| 日本免费a在线| 天天躁日日操中文字幕| 成年人黄色毛片网站| 国产aⅴ精品一区二区三区波| 在线播放国产精品三级| 波多野结衣巨乳人妻| 动漫黄色视频在线观看| 一区福利在线观看| 久久久精品欧美日韩精品| 国产真实乱freesex| 最新美女视频免费是黄的| 婷婷亚洲欧美| 日韩成人在线观看一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 波多野结衣高清无吗| 香蕉av资源在线| 午夜老司机福利剧场| 狂野欧美激情性xxxx| 啦啦啦韩国在线观看视频| 国产精品亚洲一级av第二区| 3wmmmm亚洲av在线观看| 99久久精品热视频| 嫩草影院精品99| 中亚洲国语对白在线视频| 一级黄片播放器| 夜夜躁狠狠躁天天躁| 小说图片视频综合网站| 国产真人三级小视频在线观看| 日韩欧美精品v在线| 亚洲国产欧美人成| 嫁个100分男人电影在线观看| 床上黄色一级片| 亚洲av第一区精品v没综合| tocl精华| 蜜桃久久精品国产亚洲av| 久久久久免费精品人妻一区二区| bbb黄色大片| 国产单亲对白刺激| xxx96com| 亚洲av二区三区四区| 久久国产乱子伦精品免费另类| 婷婷精品国产亚洲av在线| 黄片小视频在线播放| 亚洲av熟女| 国产精华一区二区三区| 桃红色精品国产亚洲av| 18禁黄网站禁片午夜丰满| 蜜桃亚洲精品一区二区三区| 免费人成在线观看视频色| 欧美三级亚洲精品| 国产真实伦视频高清在线观看 | 色在线成人网| 亚洲熟妇中文字幕五十中出| 精品久久久久久久末码| 国产精品亚洲av一区麻豆| 中文字幕高清在线视频| 国产99白浆流出| 国产精品,欧美在线| 亚洲精品影视一区二区三区av| 午夜福利视频1000在线观看| av女优亚洲男人天堂| 首页视频小说图片口味搜索| 亚洲人成网站高清观看| 国产伦在线观看视频一区| h日本视频在线播放| 97人妻精品一区二区三区麻豆| 此物有八面人人有两片| 麻豆国产97在线/欧美| 欧美成人性av电影在线观看| svipshipincom国产片| 看免费av毛片| 日韩免费av在线播放| 色噜噜av男人的天堂激情| 久久6这里有精品| 免费观看的影片在线观看| 亚洲人成网站高清观看| 内射极品少妇av片p| 人人妻,人人澡人人爽秒播| 久久精品综合一区二区三区| 久久九九热精品免费| 国产高清videossex| 日本成人三级电影网站| 香蕉av资源在线| 久久性视频一级片| 淫妇啪啪啪对白视频| 亚洲 欧美 日韩 在线 免费| 香蕉av资源在线| 麻豆国产av国片精品| 成人精品一区二区免费| 国产精品影院久久| 超碰av人人做人人爽久久 | 亚洲五月天丁香| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品影院久久| 91在线观看av| 国内久久婷婷六月综合欲色啪| 亚洲欧美日韩高清在线视频| 中文亚洲av片在线观看爽| 草草在线视频免费看| 老鸭窝网址在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产精品1区2区在线观看.| 精品国产美女av久久久久小说| 日日干狠狠操夜夜爽| 天堂网av新在线| 午夜福利高清视频| 变态另类丝袜制服| 女警被强在线播放| a级一级毛片免费在线观看| 12—13女人毛片做爰片一| 欧美性猛交黑人性爽| 三级毛片av免费| 国产探花在线观看一区二区| 麻豆久久精品国产亚洲av| 国产av不卡久久| 男插女下体视频免费在线播放| 午夜老司机福利剧场| 校园春色视频在线观看| 少妇熟女aⅴ在线视频| 国产一区二区三区视频了| 久久久久久久久久黄片| 男人和女人高潮做爰伦理| 一区二区三区国产精品乱码| 日韩欧美精品v在线| 欧美激情在线99| 99久久成人亚洲精品观看| 亚洲成人久久性| 亚洲精品色激情综合| 桃红色精品国产亚洲av| 级片在线观看| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利高清视频| 国产一区二区在线av高清观看| 一个人免费在线观看电影| 国产美女午夜福利| 成人精品一区二区免费| 免费看a级黄色片| 欧美日韩乱码在线| 精品午夜福利视频在线观看一区| 国产精品一区二区免费欧美| 三级国产精品欧美在线观看| 国产一区二区在线av高清观看| 国产精品久久久久久人妻精品电影| netflix在线观看网站| 久久久国产成人精品二区| 在线播放无遮挡| 日本黄色视频三级网站网址| 波多野结衣高清作品| 日韩亚洲欧美综合| 99久久99久久久精品蜜桃| 91字幕亚洲| 制服丝袜大香蕉在线| 亚洲不卡免费看| 亚洲中文字幕一区二区三区有码在线看| 日韩成人在线观看一区二区三区| 国产免费av片在线观看野外av| 国产精品日韩av在线免费观看| 国产高清激情床上av| 国产国拍精品亚洲av在线观看 | 成人欧美大片| 亚洲一区二区三区不卡视频| 精品福利观看| 99久久九九国产精品国产免费| 黑人欧美特级aaaaaa片| 99精品在免费线老司机午夜| 91在线观看av| 少妇熟女aⅴ在线视频| 国产视频一区二区在线看| 日韩欧美 国产精品| 午夜免费成人在线视频| 欧美在线黄色| 日韩欧美三级三区| 90打野战视频偷拍视频| 日本 欧美在线| 国产淫片久久久久久久久 | 成人无遮挡网站| 久久精品国产99精品国产亚洲性色| 亚洲欧美精品综合久久99| 午夜视频国产福利| 精品国产亚洲在线| 亚洲无线在线观看| svipshipincom国产片| 国内精品美女久久久久久| 狠狠狠狠99中文字幕| 两个人看的免费小视频| 无遮挡黄片免费观看| 久久久久国内视频| 亚洲成av人片在线播放无| 亚洲美女视频黄频| 欧美日韩福利视频一区二区| 老司机在亚洲福利影院| 99国产综合亚洲精品| 亚洲天堂国产精品一区在线| 老汉色∧v一级毛片| 国产三级在线视频| 亚洲黑人精品在线| 日日摸夜夜添夜夜添小说| eeuss影院久久| 亚洲人与动物交配视频| 国产精品亚洲一级av第二区| 国产精品自产拍在线观看55亚洲| 别揉我奶头~嗯~啊~动态视频| 亚洲在线自拍视频| 色噜噜av男人的天堂激情| 别揉我奶头~嗯~啊~动态视频| 精品国内亚洲2022精品成人| eeuss影院久久| 国产真实乱freesex| 可以在线观看毛片的网站| 日本 av在线| 亚洲人成电影免费在线| 亚洲无线观看免费| 久久中文看片网| 91麻豆av在线| av天堂在线播放| 99riav亚洲国产免费| 日韩国内少妇激情av| 亚洲国产欧美人成| 天堂影院成人在线观看| 在线视频色国产色| 51国产日韩欧美| av天堂中文字幕网| 男女之事视频高清在线观看| 国产 一区 欧美 日韩| 国产一区在线观看成人免费| 一级作爱视频免费观看| 欧美区成人在线视频| 国产精品精品国产色婷婷| 淫妇啪啪啪对白视频| 国产一区在线观看成人免费| 国产成人av教育| 夜夜看夜夜爽夜夜摸| av视频在线观看入口| 小蜜桃在线观看免费完整版高清| 国产精品亚洲一级av第二区| 亚洲国产欧美人成| 蜜桃亚洲精品一区二区三区| 女同久久另类99精品国产91| 久久久国产精品麻豆| 欧美一级毛片孕妇| 我的老师免费观看完整版| 偷拍熟女少妇极品色| av黄色大香蕉| 国产精品影院久久| 久久精品91无色码中文字幕| 一区二区三区国产精品乱码| 少妇丰满av| 久久久国产成人免费| 激情在线观看视频在线高清| 五月玫瑰六月丁香| 欧美成人免费av一区二区三区| 怎么达到女性高潮| 婷婷六月久久综合丁香| 久久性视频一级片| 无限看片的www在线观看| 久久久久久久久久黄片| 法律面前人人平等表现在哪些方面| 少妇高潮的动态图| 亚洲av电影在线进入| 一级毛片女人18水好多| 亚洲片人在线观看| e午夜精品久久久久久久| 身体一侧抽搐| 久久天躁狠狠躁夜夜2o2o| 久久久久免费精品人妻一区二区| 国产三级中文精品| 欧美在线一区亚洲| 亚洲国产欧美网| 国产成人福利小说| 国产又黄又爽又无遮挡在线| 最新中文字幕久久久久| 国产亚洲欧美98| 午夜影院日韩av| bbb黄色大片| 日韩免费av在线播放| 毛片女人毛片| 手机成人av网站| 中文字幕人妻丝袜一区二区| 美女高潮的动态| 亚洲国产精品合色在线| 国产精品免费一区二区三区在线| 黄色视频,在线免费观看| 国产精品亚洲av一区麻豆| 国产黄片美女视频| 床上黄色一级片| 亚洲第一欧美日韩一区二区三区| 国产综合懂色| 亚洲真实伦在线观看| 18美女黄网站色大片免费观看| 无人区码免费观看不卡| 手机成人av网站| 麻豆国产97在线/欧美| 色av中文字幕| 亚洲人成网站在线播放欧美日韩| 婷婷精品国产亚洲av在线| 男女视频在线观看网站免费| 日韩大尺度精品在线看网址| 免费在线观看影片大全网站| 听说在线观看完整版免费高清| 99热6这里只有精品| 成人国产一区最新在线观看| 给我免费播放毛片高清在线观看| 一区二区三区国产精品乱码| 国产成年人精品一区二区| 久久国产精品人妻蜜桃| 美女黄网站色视频| 精品久久久久久久久久久久久| 日韩欧美 国产精品| 久久九九热精品免费| 全区人妻精品视频| 俺也久久电影网| 亚洲电影在线观看av| 亚洲成av人片免费观看| 国产高潮美女av| 在线观看免费午夜福利视频| 中文字幕av成人在线电影| 三级毛片av免费|