鄭賢中，於潛軍，徐小明

(武漢工程大學機電工程學院，湖北 武漢 430205)

0 引 言

熱風爐是鋼鐵冶煉設備高爐的一個重要組成部分[1-2].熱風爐的主要作用是把鼓風加熱到要求的溫度，用以提高高爐的效益和效率.但是在冶煉過程中經(jīng)常要調(diào)節(jié)氣體流量的大小，以配合各種工況下鋼鐵的冶煉.為了實現(xiàn)自由調(diào)節(jié)氣流目的，需要用閥或閥組的開啟或關閉.其中閥在啟閉瞬間產(chǎn)生巨大的噪音，而噪聲又已成為威脅人類生存的三大公害之一，噪聲對人們的心理和生理都有嚴重的影響.故需要對閥組進行聲學仿真分析，以期能得到閥組的噪音值，從而判定噪聲危害的程度，而該噪聲產(chǎn)生的主要原因是振動，又由于閥或閥組在此過程中的使用頻率高，屬于易發(fā)生損壞失效的部件，故對熱風爐閥做聲固耦合分析[3-5]對減噪和強度校核有著重要的實際意義.

筆者以某熱風爐閥為研究對象進行聲固耦合分析，主要從模態(tài)振型和聲固耦合[6]等方面對閥體和閥周圍區(qū)域流體進行分析.重點計算了閥體及閥體周圍的聲場，得到了閥體周圍一定區(qū)域內(nèi)噪音值，以期通過改變結構以達到減噪的效果，并期望此計算結果能對熱風爐閥減噪有一定的參考意義.

1 聲振耦合有限元分析理論基礎

1.1 聲學Helmholtz波動方程

解決聲學外聲場問題的最終目的就是要求解某個特定問題的Helmholtz(人名：亥姆霍茲)波動方程的解，而聲學Helmholtz方程是結合了聲波的連續(xù)方程、運動方程和物態(tài)方程推導得來的[7].

▽2p(x，y，z)-k2p(x，y，z)=-jρ0ωq(x,y,z)

(1)

1.2 耦合聲學有限元方程

聲音作用于結構上的聲壓載荷可以看作是附加的法線載荷，可得到動力學方程如下：

(Ks+jωCs-ω2Ms)·{ui}+Kc{pi}={Fsi}

(2)

在流體和結構耦合的位置處，結構法線方向的振動速度與流體法線方向的振動速度相同，在流固耦合交界面處，結構的振動速度可以看作是聲音的附加速度輸入，這時聲學方程變化為如下方程(3)所示：

(Ka+jωCa-ω2Ma)·{pi}-ω2Mc{ui}={Fai}

(3)

將方程(2)、(3)寫成一個矩陣的形式，并進行耦合處理，得到耦合聲學方程：

(4)

2 微穿孔板消聲閥體聲場計算

2.1 工程閥體概述

圖1所示即為某熱風爐管道的一個閥體，當氣態(tài)流體從左端入口流入，隨活塞推移而經(jīng)穿孔柱板流入腔體，該過程中噪音主要源自氣體流動時對固體管壁產(chǎn)生沖擊而產(chǎn)生的機械振動噪音和氣體湍流振動的氣動噪音.該閥結構是依據(jù)小孔噴流消音機理，小孔噴注消音的設計機理是根據(jù)科學院聲學研究所馬大猷教授等人提出的小孔噴注噪聲極其控制理論[8-9]，從發(fā)聲機理上使它的干擾噪聲減少，由于噴注噪聲峰值頻率與噴口直徑成反比，若噴口直徑變小，噴口輻射的噪聲能量將從低頻移向高頻，于是低頻噪聲被降低，高頻噪聲反而增高.經(jīng)試驗表明，當孔徑d≤4 mm時，人耳能聽到的低頻噪音能被降低，以達到減噪的作用；同時若在中間部分加入一個活塞來回的滑動，以達到控制流量的效果，故該閥體兼有消音減壓的效果.

故選用孔徑為4 mm的穿孔板焊接在熱風爐管道中，左端是管道入口，右端是管道出口，中間部分為帶孔柱板，氣體經(jīng)左端通過柱板再進入大容積腔體，最后再通過管道出口流出，整個過程既實現(xiàn)了閥體減壓的效果，也實現(xiàn)了減噪的效果.

圖1 閥體結構示意圖Fig.1 Schematic drawing of the valve structure

2.2 有限元計算

2 2.1建立幾何模型 建立如圖2所示的帶孔柱板，其具體尺寸為內(nèi)徑150 mm，外徑160 mm，底板厚度為10 mm，然后再在其外圍建立一個半徑為200 mm的實心圓柱體，將實心圓柱體減去帶孔柱板，得到如圖3的實體流體模型.再在外部建立一個半徑為600 mm的球形區(qū)域用來模擬周圍的聲場，也即是周圍的流體.結構模型如圖4所示，其模型透視如圖5所示.

圖2 帶孔柱板圖Fig.2 The chart of cylindrical sheet with holes

2.2.2 建立有限元模型 對于三維聲學問題，ANSYS(有限元分析軟件)聲場分析指定了兩種單元類型：三維模型的流體部分分別使用Fluid30(ANSYS軟件中的一種單元)單元和Fluid130(ANSYS軟件中的另一種流體單元)，且Fluid130單元必須和Fluid30一起使用，用來構造包圍Fluid30單元的無限外殼[10].利用兩種單元類型可以構造流體部分的模型，然后利用相應的結構單元Solid185構造固體模型.只有Fluid30單元才能與結構單元相接觸；Fluid130單元只能與Fluid30單元相接觸，而不能直接與結構單元接觸.由于ANSYS的以上功能，故分別選用Solid185單元來模擬結構模型，選用Fluid30來模擬與固體接觸的流體部分，選用Fluid130單元用來構造無限外殼.

結構模型采用普通低碳鋼制作，其材料屬性分別為彈性模型Ex為2.1×105MPa，泊松比為0.3，密度為7.85×10-9t/mm3，流體模型需要定義速度和密度，速度為3.44×105mm/s，密度為1.21×10-12t/mm3，無限外殼的聲速為2×103mm/s.

圖3 內(nèi)部流體Fig.3 Inner fluid

圖4 周圍聲場模型圖Fig.4 The model figure of Acoustic field surrounding

圖5 線可視模型圖 Fig.5 The model figure of line visibility

由于柱板有很多小孔，若是不作處理自由劃分網(wǎng)格，小孔邊緣位置網(wǎng)格質量很差，若是要保證小孔邊緣的網(wǎng)格質量，又導致整個柱板的單元數(shù)過多而加大計算負擔甚至無法計算.故采用Solid185單元對柱板上小孔進行切分之后再進行網(wǎng)格劃分，切分之后的柱板結構如圖6所示,切分之后對柱板兩端和底板進行掃掠劃分，然后對有小孔的中間柱板部分采用智能劃分，劃分網(wǎng)格之后的網(wǎng)格如圖7所示.對中間部分的流體采用Fluid30單元智能方式劃分，劃分之后的網(wǎng)格如圖8所示，中間網(wǎng)格局部視圖如圖9所示.對外圍球體部分采用Fluid30單元自由劃分，單元尺寸設置為40 mm,劃分網(wǎng)格圖如圖10所示.

圖6 切分之后模型圖Fig.6 The model figure after segmentation

圖7 切分之后網(wǎng)格圖Fig.7 The mesh chart after segmentation

圖8 中間流體網(wǎng)格圖Fig.8 The intermediate fluid mesh chart

圖9 中間流體網(wǎng)格正視圖Fig.9 The front view chart of intermediate fluid mesh

2.2.3加載與求解 劃分網(wǎng)格后，由于流體和固體是兩個單獨的部分，只有通過耦合才能使流體和固體進行有效關聯(lián)從而實現(xiàn)計算，故在流體和固體交界面處定義流固耦合面，即FSI(流固耦合)如圖11和圖12所示.然后在柱板入口處施加0.25 MPa的入口壓力.再在無限外殼面定義吸聲屬性，此處吸聲屬性定義為1，即全吸聲屬性，如圖13所示.最后設置計算參數(shù)并計算，設置步長為20，范圍從500至1 000.

圖10 中間流體網(wǎng)格正視圖Fig.10 The front view chart of intermediate fluid grid

圖11 流固耦合交界面(整體視圖)Fig.11 Interface of fluid-structure interaction

圖12 流固耦合交界面(局部視圖)Fig.12 Interface of fluid-structure interaction

圖13 無限外殼定義吸聲屬性Fig.13 The sound absorption properties chart of infinite shell

3 聲場分析與結論

隨著人們對其所居住環(huán)境越來越重視，環(huán)保部門開始慢慢地關注和管制噪音的危害，而且出臺了相關限制噪音的標準和規(guī)范.工業(yè)企業(yè)環(huán)境標準所允許的最大噪音值為90 dB,實際普通熱風爐閥在啟閉時產(chǎn)生的噪音高達120 dB,嚴重地污染了周圍的環(huán)境.對該有限元模型求解計算成功后通過后處理，得到如圖14所示聲壓級云圖，從該圖中可看出最大的聲壓級值為79.2 dB.由此可知采用該新型減壓閥能從模擬仿真方面得到的噪音值滿足相關法律法規(guī)所要求控制的噪音范圍，可推知小孔噴注是一種行之有效的減噪方法，但必須要選用合理的小孔孔徑及排列方式.

圖14 聲壓級分布云圖Fig.14 Distributed cloud chart of sound pressure level

4 結 語

通過對某熱風爐系統(tǒng)中閥做聲振耦合有限元計算，并分析聲場結果，最后可得出以下幾點結論：

b.僅當小孔的直徑為4 mm左右時，消音減噪的效果較佳.若是孔徑太小，小孔容易被氣體中的雜質堵塞，若是孔徑過大，又不能將高頻噪音濾去，達不到消音減噪的效果.

c.計算過程及計算結果表明，ANSYS是一款比較實用的有限元處理軟件,其用戶界面良好,前、后處理功能強大,計算精確度高,是求解聲場分布行之有效的有限元處理軟件.

致 謝

本文研究工作得到武漢工程大學機電工程學院的鼎力支持，在此表示衷心的感謝！

參考文獻:

[1] 銀漢.現(xiàn)代熱風爐設計的若干問題[J].煉鐵學報,2002, 2(4):27-30.

YIN Han.Some questions on design of high temperature hot stove[J].Ironmaking,2002，2(4):27-30.(in Chinese)

[2] 楊俊,杜濤,蔡九菊，等.國外熱風爐發(fā)展綜述[J].沈陽工程學院學報，2005,4(12):18-20.

剛開始自駕游時，大孫女李子宜還在上小學，現(xiàn)在她已經(jīng)上高三了。這些年的出游大大鍛煉了她的各種能力。從新疆回來后，她將自己的感受寫成了14頁紙的《天路行》，獲得“國家新聞獎”。在學校，她制作的旅游幻燈片圖文并茂，很受同學們歡迎，為此，她成為學校的小記者。她的地理成績在全年級也是名列前茅。

YANG Jun，DU Tao，CAI Jiu-ju,et al.Summarization on the development of overseas hot blast store[J].Journal of Shenyang Institute of Engineering:Natural Science,2005,4(12):18-20.(in Chinese)

[3] 張軍,兆文忠,張維英.聲場-結構耦合系統(tǒng)聲壓約束下板重量優(yōu)化設計研究[J].應用力學學報，2006，23(4)：568-571．

ZHANG Jun,ZHAO Wen-zhong,ZHANG Wei-ying.Optimum weight design of plate with acoustic pressure restrict for coupled acoustic-structure systems[J].Chinese Journal of Applied Mechanics,2006，23(4):568-571.(in Chinese)

[4] 王耀前，陸森林. ANSYS在抗性消聲器分析中的應用[J]. 江蘇大學學報,2003, 24(3):53-56.

WANG Yao-qian,LU Sen-lin.Application of ANSYS software in analyzing reactive muffler[J].Journal of Jiangsu University:Natural Science Edition,2003,24(3):53-56.(in Chinese)

[5] 鄭賢中，於潛軍，周寧波，等.脈動流技術在管殼式換熱器振動分析中的應用[J].武漢工程大學學報,2013,35(8)：68-73.

ZHENG Xian-zhong,YU Qian-jun，ZHOU Ning-bo, et al.Application of pulsating flow technology in vibration analysis of shell-and tube heat exchanger[J].Journal of Wuhan Institute of Technology,2013,35(8):68-73.(in Chinese)

[6] 楊永謙，肖金生. 實用有限元分析技術[M]. 北京：機械工業(yè)出版社，2010:215-225.

YANG Yong-qian,XIAO Jin-sheng.Practical Finite Element Analysis Technology[M].Beijing:Mechanicalindustry Press,2010:215-225.(in Chinese)

[7] 李增剛，詹福良. Virtual.Lab Acoustics 聲學仿真計算高級應用實例[M]. 北京：國防工業(yè)出版社，2010:4-17.

LI Zeng-gang,ZHAN Fu-liang.Senior acoustic simulation application example of Virtual.Lab Acoustics[M].Beijing:National Defence Industry Press,2010:4-17.(in Chinese)

[8] 馬大猷,李沛滋.小孔噴注噪音和小孔消聲器[J].中國科學,1975(1):38-50.

MA Da-you,LI Pei-zi.Small hole jet noise and Small hole muffler[J].China Science,1975(1):38-50.(in Chinese)

[9] 馬大猷.現(xiàn)代聲學基礎[M].北京：科學出版社，2004：310-312.

MA Da-you.Foundation of modern acoustic[M].Beijin:Science Press,2004:310-312.(in Chinese)

[10] 張乃龍,楊文通，費仁元，等.基于ANSYS的抗性消聲器性能仿真分析[J].計算仿真,2006,23(8):306-309.

ZHANG Nai-long,YANG Wen-tong,FEI Ren-yuan,et al.Simulation and analysis of reactive muffler based on ANSYS software[J].Simulation of Computer,2006,23(8):306-309.(in Chinese)