蘇明珠，孫玉東，王鎖泉，郝夏影

（中國船舶科學(xué)研究中心船舶振動(dòng)噪聲重點(diǎn)實(shí)驗(yàn)室；江蘇綠色船舶技術(shù)重點(diǎn)實(shí)驗(yàn)室，江蘇無錫214082）

0 引 言

柔性接管是用于管路系統(tǒng)中位移補(bǔ)償和減振降噪的重要元件，其動(dòng)態(tài)特性的研究，對(duì)于管系振動(dòng)噪聲傳遞和響應(yīng)計(jì)算分析及減振降噪效果實(shí)現(xiàn)至關(guān)重要。Pinnington（1997）[1]指出，通過降低波速，由編織橡膠構(gòu)成的柔性管和波紋管能同時(shí)提供振動(dòng)和聲波的衰減。阻抗是連接激勵(lì)和響應(yīng)之間關(guān)系的動(dòng)力學(xué)參數(shù)，可以反映元件本身的動(dòng)力學(xué)特性。對(duì)于柔性接管而言，其彈性部分一般為橡膠或橡膠-金屬（非金屬）織物混合結(jié)構(gòu)，其柔性管壁與管內(nèi)液體介質(zhì)存在較強(qiáng)的聲振耦合效應(yīng)。為了考慮柔性接管結(jié)構(gòu)中振動(dòng)波和流體介質(zhì)中聲波的相互耦合，全面反映柔性接管的聲學(xué)特性，從而精確預(yù)報(bào)管路系統(tǒng)的聲振傳遞特性和元器件的減振降噪效果，需開展充液柔性接管的耦合阻抗研究。

近年來，俄羅斯克雷洛夫國家科學(xué)中心Sokolov（2014）[2]和Popkov（2012）[3]分別進(jìn)行了包含耦合阻抗在內(nèi)的柔性接管全阻抗矩陣元素的理論和測試方法研究，其中Sokolov 給出了正交各向異性和粘彈性充液軟管作為梁模型的阻抗矩陣元素的理論公式，Popkov給出了測量機(jī)械阻抗、聲阻抗和兩者耦合阻抗的測試裝置和方法。針對(duì)包含軟管的管路系統(tǒng)，Munjal和Thawani（1997）[4]根據(jù)管段傳遞矩陣，研究了軟管對(duì)入射的彎曲波、縱波和聲波的傳遞損失，以及軟管材料彈性參數(shù)、阻尼參數(shù)和長度、管段數(shù)、壁厚、內(nèi)徑、彎管角度等各種因素對(duì)隔振和聲衰減的影響。Drew 等（1998）[5]描述了包含協(xié)調(diào)器的柔性軟管的理論模型，得到了管端壓力波和流動(dòng)波關(guān)系的頻率相關(guān)阻抗矩陣，計(jì)算了完整系統(tǒng)的流體噪聲特性。在國內(nèi)，哈爾濱工程大學(xué)李帥軍（2015）[6]對(duì)輸流管路的流體-結(jié)構(gòu)相互作用研究進(jìn)展進(jìn)行了較為全面的回顧。西北工業(yè)大學(xué)的樊澤明（2013）[7]回顧了直管和彎管的充液管道線性、非線性流固耦合振動(dòng)模型及分析方法。對(duì)于軟管類管路元件，孫玉東等（2004）[8]以彈性直管為對(duì)象，不考慮泊松耦合影響，給出了軸向、橫向機(jī)械阻抗和聲阻抗的解析結(jié)果。

本文在考慮泊松耦合的基礎(chǔ)上，對(duì)充液柔性接管聲阻抗、軸向機(jī)械阻抗及耦合阻抗進(jìn)行了理論推導(dǎo)，與不考慮泊松耦合的阻抗結(jié)果進(jìn)行了對(duì)比，分析了材料和幾何參數(shù)等相關(guān)因素對(duì)軸向阻抗的影響。研究結(jié)果可作為試驗(yàn)測量的考核依據(jù)，并為試驗(yàn)裝置的設(shè)計(jì)提供參考。

1 充液柔性接管軸向耦合阻抗的理論推導(dǎo)

1.1 充液柔性接管軸向耦合方程及其通解

管壁的縱波和內(nèi)部介質(zhì)的一維聲波的耦合控制方程組為

式中：Fx、Ux、P 和Uf分別是管壁縱向力、管壁縱向運(yùn)動(dòng)速度、管路中液體動(dòng)壓力和管路中流體運(yùn)動(dòng)速度；A 和Ap分別是管道公稱面積和管壁橫截面積；E 和μ 分別是管道材料的楊氏模量和泊松比；x 為沿管路軸向坐標(biāo)；ρp是管壁材料密度；ρf為管內(nèi)流體密度；E'V為考慮管壁彈性修正后的體積模量。對(duì)于薄壁管，

式中，EV為剛性管中流體的體積模量，R為管道的平均半徑，e為管壁厚度。

對(duì)于式（1）～（4）所表示的考慮泊松耦合的聲彈性梁柔性接管軸向方程，可采用以下無量綱形式：

式中，ux為管壁縱向位移，uf為管道流體縱向位移，L為管長。

設(shè)時(shí)間簡諧因子為ejωt，進(jìn)一步采用以下無量綱參數(shù)：

采用式（6）和式（7）中的無量綱變量和參數(shù)，假設(shè)沿軸向坐標(biāo)的解有以下形式：

將式（8）代入式（1）～（4），得到λ的一元四次方程：

即對(duì)應(yīng)于泊松耦合管道中結(jié)構(gòu)和流體的無量綱縱波波數(shù)λ1和λ2，可按下式計(jì)算：

從而得到方程組的線性疊加通解為

式中，

1.2 充液柔性接管軸向全阻抗矩陣及矩陣元素的推導(dǎo)

根據(jù)管路元器件阻抗的定義，兩端狀態(tài)量之間的關(guān)系以阻抗矩陣表示：

阻抗矩陣中各元素定義為

式中，Y為式（15）左端的激勵(lì)力，X為式（15）右端的速度響應(yīng)。

由式（15）和式（16），很容易得到入口端輸入聲阻抗為

出口端輸入聲阻抗為

即Z22= -Z11，表示柔性圓管兩端的對(duì)稱性。

聲傳遞阻抗為

式（20）表示柔性圓管阻抗符合互易原理。

機(jī)械傳遞阻抗為

當(dāng)μ ≠0，即考慮泊松耦合時(shí)，所有耦合阻抗矩陣元素的推導(dǎo)結(jié)果如下。

聲-機(jī)械輸入阻抗為

聲-機(jī)械傳遞阻抗為

機(jī)械-聲輸入阻抗為

機(jī)械-聲傳遞阻抗為

2 幾何、材料參數(shù)對(duì)阻抗結(jié)果的影響

利用上面推導(dǎo)的公式，對(duì)充液柔性接管幾何和材料參數(shù)變化下包含耦合阻抗在內(nèi)的阻抗特性，進(jìn)行MATLAB程序計(jì)算?；居?jì)算參數(shù)如下：

式中：Qm,i為實(shí)測徑流值，為實(shí)測徑流平均值，m/s；Qs,i為模擬徑流值，為模擬徑流平均值，m/s。

水：ρf= 998 kg/m3，EV= 2.18 × 109Pa；

鋼管：ρps= 7 800 kg/m3，Es= 2.1× 1011Pa，μs= 0.29；

橡膠管[8]：ρpR= 1200 kg/m3，ER= 3× 109Pa，μR= 0.48；

幾何參數(shù)：L = 1 m，R = 0.05 m，e = 5 mm。

2.1 考慮耦合作用下橡膠管的各阻抗計(jì)算結(jié)果

由圖1～4可見，考慮泊松耦合作用下的橡膠管的輸入阻抗和傳遞阻抗在低頻時(shí)吻合，在高頻則出現(xiàn)差異。圖1中聲阻抗最小值所在頻率，即反共振頻率大約為100 Hz。圖3中機(jī)械阻抗最小值所在頻率約為450 Hz。圖2中聲-機(jī)械阻抗和圖4中機(jī)械-聲阻抗的輸入和傳遞阻抗同樣在低頻重合，其中輸入阻抗先共振再反共振，傳遞阻抗先反共振繼而出現(xiàn)共振，聲-機(jī)械阻抗和機(jī)械-聲阻抗的共振頻率接近，第一反共振頻率在170 Hz左右。

圖1 聲阻抗Fig.1 Acoustic impedance

圖2 聲-機(jī)械阻抗Fig.2 Acoustic-mechanical impedance

圖3 機(jī)械阻抗Fig.3 Mechanical impedance

圖4 機(jī)械-聲阻抗Fig.4 Mechanical-acoustic impedance

2.2 考慮耦合與不考慮耦合橡膠管的阻抗結(jié)果比較

當(dāng)不考慮泊松耦合時(shí)，柔性接管軸向阻抗矩陣只存在聲阻抗和機(jī)械阻抗，聲-機(jī)械阻抗和機(jī)械-聲阻抗均為零。圖5～8為考慮與不考慮泊松耦合的橡膠管的阻抗結(jié)果比較，從阻抗幅值來看，考慮耦合阻抗計(jì)算得到的阻抗值略高于不考慮耦合阻抗的。同時(shí)，考慮泊松耦合影響計(jì)算得到的聲阻抗幅值-頻率曲線會(huì)出現(xiàn)管壁結(jié)構(gòu)共振頻率，機(jī)械阻抗幅值-頻率曲線會(huì)出現(xiàn)聲共振頻率，體現(xiàn)了二者的耦合影響。

圖5 考慮耦合與不考慮耦合橡膠管的聲輸入阻抗Fig.5 Acoustic input impedance comparison between coupled and uncoupled rubber pipes

圖6 考慮耦合與不考慮耦合橡膠管的聲傳遞阻抗Fig.6 Acoustic transfer impedance comparison between coupled and uncoupled rubber pipes

圖7 考慮耦合與不考慮耦合橡膠管的機(jī)械輸入阻抗Fig.7 Mechanical input impedance comparison between coupled and uncoupled rubber pipes

圖8 考慮耦合與不考慮耦合橡膠管的機(jī)械傳遞阻抗Fig.8 Mechanical transfer impedance comparison between coupled and uncoupled rubber pipes

2.3 幾何參數(shù)對(duì)阻抗結(jié)果的影響

圖9～12為長度變化對(duì)橡膠管阻抗值變化的影響?？梢钥闯觯S著長度的增加，共振頻率降低，反共振頻率也隨之減小。在低頻范圍內(nèi)，充液橡膠管的阻抗值隨長度的增加而減小。圖13～16 為橡膠管口徑變化對(duì)阻抗值的影響，聲阻抗和耦合阻抗的共振與反共振頻率隨口徑的增大而減小，機(jī)械阻抗與之相反。同樣，在低頻范圍內(nèi)，聲阻抗和耦合阻抗值隨口徑的增加而減小，機(jī)械阻抗值隨口徑的增大而增大。

圖9 長度變化對(duì)聲阻抗的影響Fig.9 Influence of length variation on acoustic impedances

圖10 長度變化對(duì)機(jī)械阻抗的影響Fig.10 Influence of length variation on mechanical impedances

圖11 長度變化對(duì)聲-機(jī)械阻抗的影響Fig.11 Influence of length variation on acoustic-mechanical impedances

圖12 長度變化對(duì)機(jī)械-聲阻抗的影響Fig.12 Influence of length variation on mechanical-acoustic impedances

圖13 口徑變化對(duì)聲阻抗的影響Fig.13 Influence of diameter variation on acoustic impedances

圖14 口徑變化對(duì)機(jī)械阻抗的影響Fig.14 Influence of diameter variation on mechanical impedances

圖15 口徑變化對(duì)聲-機(jī)械阻抗的影響Fig.15 Influence of diameter variation on acoustic-mechanical impedances

圖16 口徑變化對(duì)機(jī)械-聲阻抗的影響Fig.16 Influence of diameter variation on mechanical-acoustic impedances

2.4 充液鋼管與充液橡膠管的阻抗結(jié)果比較

圖17～20是橡膠管與鋼管的阻抗結(jié)果對(duì)比。結(jié)果表明，在低頻范圍內(nèi)，鋼管的阻抗值明顯高于橡膠管的阻抗值，其共振頻率也高于橡膠管。鋼管的聲阻抗最小值所在頻率約為340 Hz，機(jī)械阻抗最小值所在頻率超過1 000 Hz，聲-機(jī)械阻抗和機(jī)械-聲阻抗的最小值所在頻率大約為630 Hz。

圖17 鋼管與橡膠管的聲阻抗比較Fig.17 Comparison of acoustic impedances between steel pipe and rubber pipe

圖18 鋼管與橡膠管的機(jī)械阻抗比較Fig.18 Comparison of mechanical impedances between steel pipe and rubber pipe

圖19 鋼管與橡膠管的聲-機(jī)械阻抗比較Fig.19 Comparison of acoustic-mechanical impedances between steel pipe and rubber pipe

圖20 鋼管與橡膠管的機(jī)械-聲阻抗比較Fig.20 Comparison of mechanical-acoustic impedances between steel pipe and rubber pipe

3 結(jié) 論

本文根據(jù)充液柔性接管軸向耦合控制方程，推導(dǎo)獲得了軸向阻抗的解析解，計(jì)算了各種相關(guān)參數(shù)對(duì)阻抗結(jié)果的影響，對(duì)比分析了考慮與不考慮泊松耦合及長度、口徑、橡膠管與鋼管材料不同時(shí)的阻抗計(jì)算結(jié)果，得出以下主要結(jié)論：

（1）充液柔性圓管軸向阻抗矩陣由聲阻抗、機(jī)械阻抗和聲-機(jī)械耦合阻抗、機(jī)械-聲耦合阻抗等子矩陣組成，各子矩陣元素遵循對(duì)稱性和互易性原理。

（2）耦合阻抗的算例表明，輸入阻抗先共振再反共振，傳遞阻抗先反共振繼而出現(xiàn)共振，聲-機(jī)械阻抗和機(jī)械-聲阻抗的共振及反共振頻率接近。

（3）考慮泊松耦合時(shí)，聲阻抗和機(jī)械阻抗會(huì)受到管壁結(jié)構(gòu)共振與聲共振的影響，直觀體現(xiàn)二者的耦合作用。

（4）隨著充液軟管長度的增加，共振和反共振頻率降低。充液軟管口徑的增大，導(dǎo)致聲阻抗與耦合阻抗的共振與反共振頻率減小，機(jī)械阻抗結(jié)果與之相反。鋼管的阻抗值明顯大于橡膠管的阻抗值，其共振頻率也高于橡膠管，符合一般規(guī)律。

后續(xù)將對(duì)本文結(jié)果進(jìn)行試驗(yàn)測量和驗(yàn)證，并對(duì)橫向耦合阻抗、正交各向異性軟管耦合阻抗等開展更多研究工作。