第43卷第7期中南大學學報(自然科學版)Vol, 43 No. 72012年7月Jourmal of Central South University(Science and Technology)July 2012熱聲熱機結(jié)構(gòu)動力學分析劉益才,顏鵬,馬衛(wèi)武,陳思明,陳麗新,李照龍(中南大學能源科學與工程學院,湖南長沙,410083)摘要:利用商業(yè)軟件 ANSYS建立熱聲熱機模型,從結(jié)構(gòu)動力學角度進行模擬研究,對其進行結(jié)構(gòu)靜力、模態(tài)和諧響應分析。在熱聲系統(tǒng)的關(guān)鍵位置布置監(jiān)測點,得到熱聲系統(tǒng)內(nèi)部位移、應力和應變等參數(shù)隨壓力和頻率的變化關(guān)系。研究結(jié)果表明:隨著充氣壓力的提升,應力應變基本呈現(xiàn)線性增大的趨勢,且峰值點一一對應。隨著結(jié)構(gòu)模態(tài)的提高,熱聲熱機的應力和應變因位置的不同而存在差異;熱聲系統(tǒng)的3階和4階模態(tài)下固有頻率接近氣體介質(zhì)的振蕩頻率:由諧振管內(nèi)表面布置的監(jiān)測點所得結(jié)果可知,高密度的氣體介質(zhì)的振蕩會在回熱器中產(chǎn)生振動激勵和聲激勵,并作用于諧振管,使聲波和固體介質(zhì)產(chǎn)生共振和耦合關(guān)鍵詞:熱聲熱機;結(jié)構(gòu)動力學;模態(tài);諧響應中圖分類號:TB65文獻標志碼:A文章編號:1672-7207(2012)07-2844-06Analysis of structural dynamics in thermoacoustic engineLIU Yi-cai, YAN Peng, MA Wei-wu, CHEN Si-ming, CHEN Li-xin, LI Zhao-longSchool of Energy Science and Engineering, Central South University, Changsha 410083, China)Abstract: Mechanism of thermoacoustic engine was studied using ANSYS software. Simulation was made from theperspective of structural dynamics with the analysis of static force of structure and modal and harmonic responsMonitoring points were arranged at the crucial positon ofvariation of the displacement, stress and strain versuscncy or pressure. The results show that with thepressure, the trends of stress and strain generally linearly increase, and the peaks are one-to-one corrWithimprovement of modal, thermoacoustic stress and deformation are dissimilar in different locations; the natural frequencyof thermoacoustic system is close to the oscillation frequency of the gas medium at modal 3 and 4; the result ofmonitoring point on the inner surface of the resonator shows that high-density gas medium would produce vibrationexcitation and sound excitation, which makes the acoustic and structural modes couple with each other.Key word: thermoacoustic engine; structural dynamics; modal; harmonic response近年來國內(nèi)外學者從多個角度討論熱聲振蕩的產(chǎn)構(gòu)聲產(chǎn)生的一個重要原因即振動激勵1。根據(jù)結(jié)構(gòu)生機理研究表明熱聲熱機的聲學特性和頻率存在聲理論,結(jié)構(gòu)聲是熱聲熱機所固有的聲學特性,在發(fā)定關(guān)系,在某一頻率下聲學特性表現(xiàn)最佳,當某生結(jié)構(gòu)振動時機體是以多個模態(tài)振動。因此,當機體個力作用于結(jié)構(gòu)時,會引起結(jié)構(gòu)自身的振動,這是結(jié)振動頻率與熱聲系統(tǒng)中氣體介質(zhì)的振動發(fā)生耦合時,收稿日期:201-0916:修回日期:2012-0108中國煤化工基金項目:國家自然科學基金資助項目(50676110):中南大學博士后基金資助項目(2010年)岸改革課題資助項目(2011jg48)CNMH通信作者:馬衛(wèi)武(1974)男,湖南隆回人,博士副教授,從事制冷及空調(diào)系統(tǒng)、熱聲熱機等研究;電話:13974871366:E-mail:mweru@csu.edu.cn劉益才熱聲熱機結(jié)構(gòu)動力學分析845勢必會對系統(tǒng)效率造成影響。黃謙等21利用Ma+ Cc +Ku=F((4)Solidworks和 ANSYS軟件對平板和回熱器的熱聲機理從結(jié)構(gòu)動力學的角度進行了模擬。結(jié)果顯示,回熱其中:M為結(jié)構(gòu)質(zhì)量矩陣;a為節(jié)點加速度矢量:C器的振動中存在多種模態(tài),其中在低頻段各點應力應為結(jié)構(gòu)阻尼矩陣;c為節(jié)點速度矢量:K為結(jié)構(gòu)剛度變基本吻合,而在高頻段各點差異逐漸增大;在一階矩陣;"為節(jié)點位移矢量,當進行模態(tài)分析時,F(固有頻率3300H附近,回熱器內(nèi)各點的應力和應變?yōu)榱?包含流-固耦合自由振動,有無阻尼自由振動達到了最大。然而將平板或回熱器從熱聲系統(tǒng)中獨立根據(jù)具體情況而定:當進行諧響應分析時,FO)和u出來孤立分析,會產(chǎn)生諸多誤差。如回熱器的套筒在都為諧函數(shù)壓力載荷下發(fā)生了較大應變,而實際上,當回熱器在取初始頻率G100Hz時,氮氣的半波長17m作諧振管的約束下發(fā)生應變的幅度很小;且就一般熱聲為諧振管的長度。固體介質(zhì)為鍛造不銹鋼,其密度p=系統(tǒng)而言系統(tǒng)工作時的諧振頻率無法達到如此之高。7800kg/m3,彈性模量E=200GNm2,泊松比=028本文作者對熱聲整機進行了三維建模,利用商業(yè)軟件板疊回熱器長度為60mm,外徑為44mm,平板厚度ANSYS對其進行了結(jié)構(gòu)動力學分析。和板間距均為06mm,諧振管厚度為2mm。諧振管各封閉端和連接處利用Glυe命令進行焊接。由于平板模型概述是直接焊接在諧振管的內(nèi)表面,因此在平板和諧振管的相交面設置了DOF面約束,并在諧振管內(nèi)表面施加了內(nèi)壓力。結(jié)構(gòu)中的聲波指的是某一物理量(力、位移、速度根據(jù)熱聲熱機結(jié)構(gòu)特性和諧振管內(nèi)部的受力情或能量)以結(jié)構(gòu)為媒體向某一方向傳播,結(jié)構(gòu)本身作為況,單元類型采用 PLANE82和 SOLID95。模型總共一個整體并不運動,在一維情況下,假設結(jié)構(gòu)的質(zhì)包含31607個單元及70889個節(jié)點該模型基于以下量元僅在x方向分布,且波動沿x方向以速度c傳播,幾種假設:(1)諧振管內(nèi)均勻分布氣體介質(zhì),不考慮則結(jié)構(gòu)聲可用下式表示:周圍流體介質(zhì)的黏性,較之在實際流體中板疊的振動E(x, r=Ae(of+p-bar)(1)幅度衰減較慢;(2)固體材料的物性參數(shù)如彈性模量、密度等不隨外界因素的變化而變化;(3)設熱聲熱機式中:為時間:為角頻率:A為常數(shù),k0為波處于恒壓環(huán)境中數(shù),表示聲波相位沿空間的變化規(guī)律。結(jié)構(gòu)中聲波的表達式與流體中聲波的表達式一致,各物理量意義也相似。結(jié)構(gòu)中的聲波存在縱波、橫波和彎曲波,縱波2數(shù)值研究是指結(jié)構(gòu)中微小質(zhì)量單元的運動方向和聲波的傳播方向一致,表達式為21結(jié)構(gòu)靜力分析e分別選取了03,05和07MPa3種充氣壓力分ax B at析熱聲熱機內(nèi)部的結(jié)構(gòu)響應。靜力分析可以模擬充氣式中:B為材料剛度的相關(guān)參數(shù),與材料彈性模量E開始時到熱聲熱機起振前系統(tǒng)的結(jié)構(gòu)參數(shù),此時假定及泊松比a有關(guān),B=E(1-)載荷和響應是固定不變的,即載荷和結(jié)構(gòu)的響應隨時(1+a(1-2a)間的變化非常緩慢。圖1所示為系統(tǒng)中最大應力和最固體中,質(zhì)量微元振動的方向可以和聲波傳播方大應變隨充氣壓力的變化關(guān)系。隨著壓力的提升,應向垂直,即橫波。橫波導致質(zhì)量做元發(fā)生形變導致切力應變基本呈現(xiàn)線性增大的趨勢,且峰值點一一對應應力,故被稱為橫向切變波:圖2和3所示為05MPa壓力下熱聲熱機不同角度的應力云圖。分析結(jié)果表明系統(tǒng)內(nèi)部的應力、應變ax G at分布是基本吻合的,總體看來系統(tǒng)中外表面受力微小,式中:G為剪切應變和剪切應力的比值,稱之為切變應力和應變主要八個產(chǎn)中國煤工以看到:由于填充同性多孔器所處位置模量,G=ECNMH2(1+σ)的諧振管外部應和m又共,而之前單獨對回動力學基本控制方程如下:熱器的結(jié)構(gòu)動力分析中,該位置請振管的應變不容忽中南大學學報(自然科學版)第43卷諧振管中段呈現(xiàn)了上凸下凹的彎曲狀態(tài),且整個諧振3.6管徑向應力由里至外逐漸加強。由圖3可知:從宏觀上看板疊整體各個面的中部位置幾乎沒有受到應力作用,應力載荷集中分布在四周邊緣處,各個面應力等28值圖呈現(xiàn)出由里及外逐漸增大的特點;而對單獨一塊20平板而言,應力載荷在軸向和徑向均呈現(xiàn)出“中間小兩側(cè)大”的特點,且隨著平板表面積的增大,應力載荷的極大值與極小值之比也逐漸增大,這與之前的研10究相符[2固體介質(zhì)應力和應變最終影響回熱器和諧振管的0.50.60.7聲學特性,進一步影響了回熱器中氣體微團的運動。壓力MPa由結(jié)構(gòu)聲輻射理論,在熱聲熱機的固體介質(zhì)(如平板圖13種不同壓力下熱聲熱機最大應力和最大應變的分布中彎曲波所引起的振動幅度最大,而彎曲波的相速度Fig-1 Distributions of maximum stress point and strain point與頻率有關(guān),因此都存在色散現(xiàn)象。任意形狀的彎曲in thermoacoustic engine under three different pressures波經(jīng)過一段時間的傳播后,該聲波的各個頻率分量的相對相位的變化有可能導致其形狀的變化,從而進一步影響到結(jié)構(gòu)和流體之間的聲耦合。因而分析系統(tǒng)的模態(tài)和諧響應尤為重要。2.2模態(tài)以及諧響應分析模態(tài)分析可以確定熱聲熱機固體介質(zhì)的固有頻率和振型,同時也是諧響應分析的出發(fā)點。根據(jù)熱聲系統(tǒng)情況,選取充氣壓力為05MPa,相位角為0°。表1列出了系統(tǒng)1-8階模態(tài)的頻率值。在不同振動模態(tài)時,熱聲系統(tǒng)各位置的應變不同,說明靜力下熱聲熱應力(N-m2機的應變起主要作用,而動力下回熱器的應變是靜力63818813.xa3x×m和波動兩者的疊加。可以看到:熱聲系統(tǒng)的3階和4圖20.5MPa時熱聲熱機的應力云圖階模態(tài)下固有頻率接近氣體介質(zhì)的振蕩頻率,因而可Fig2 Stress contour of thermoacoustic engine under 0.5MPa以重點觀察3階和4階模態(tài)時熱聲系統(tǒng)的結(jié)構(gòu)動力參數(shù)表1熱聲熱機各階模態(tài)頻率Table 1 Frequencies of each modal in thermoacoustic engine模態(tài)頻率/Hz載荷步數(shù)子步數(shù)累積部數(shù)17.37431077106o883755圖30.5MPa時回熱器的應力云圖6300.3006Fig 3 Stress contour of regenerator under 0.5 MPa5096008576.960視13;諧振管其他位置所受應力較大,且同樣分布均勻:管左邊封閉端的不銹鋼圓板上應力和應變呈同心諧響應分YH史8中國煤化工NMH承受隨時間按圓狀由外及內(nèi)逐漸增大,并在圓板中心位置達到峰值。正弦(簡諧)規(guī)律變化的載荷時穩(wěn)態(tài)響應的一種算法劉益才,等:熱聲熱機結(jié)構(gòu)動力學分析分析的目的是計算一定材料和結(jié)構(gòu)在一定頻率范圍內(nèi)的周期性載荷作用下的響應情況。該技術(shù)計算結(jié)構(gòu)的穩(wěn)態(tài)受迫振動時,不考慮結(jié)構(gòu)發(fā)在激勵開始時的瞬態(tài)振動。諧響應分析是一種線性分析,非線性特性被忽210略。在熱聲熱機內(nèi)部重點研究位置選取了監(jiān)測點來研究結(jié)構(gòu)振動的響應情況,其中:點a設定在板疊內(nèi)中間平板的上端面幾何中心;點b位于該板的幾何中心;點c位于該板下端面的幾何中心;點d位于諧振管封閉端端面中心,也就是系統(tǒng)最大應力位置處;點e位于諧振管上部內(nèi)表面。圖48所示為上述點處應變隨頻率的變化曲線,其中對點d和點e處曲線進行了擬合。模態(tài)圖46分別代表了板疊中間板xy和z方向上應1-x方向;2-y方向;3-z方向變隨頻率的變化關(guān)系。由圖46可以看到:位于平板圖6c點處應變隨頻率的變化上、中和下部的3個監(jiān)測點的應變變化高度一致:xFlg. Curves of strain versus frequency at point和z方向的應變均非常微小,且隨著頻率的增大亦無十分明顯的變化,只是當頻率超過500H以后,上端面和下端面分別出現(xiàn)了應變減小和應變增大的趨勢,不過變化幅度不大;反觀監(jiān)測點y方向的應變在300H以后明顯上升,直至500H附近達到峰值,隨后隨著頻率的增大應變不斷減小直至變?yōu)樨撝??？梢钥吹?隨著結(jié)構(gòu)頻率的上升,平板的應變變化幅度總體是不斷增大的,且各處受振動的影響相似,在垂直于平板的y方向應變最大而其他2個方向在低頻段幾乎沒有發(fā)生應變,這也與實際情況相吻合。在模態(tài)頻率為500H時平板的應力應變達到極大值,且隨著x方向;2-y方向;3-z方向頻率的進一步提高而增大,說明回熱器作用是激發(fā)和圖4a點處應變隨頻率的變化維持高階熱聲振蕩Fig 4 Curves of strain versus frequency at point a圖7和8所示為諧振管內(nèi)表面監(jiān)測點的應變隨頻率變化的擬合曲線。在熱聲系統(tǒng)這一有限結(jié)構(gòu)內(nèi),結(jié)構(gòu)聲的傳播在遇到內(nèi)表面時會發(fā)生折射、反射和衍射等現(xiàn)象,其中邊界的反射現(xiàn)象尤為突出,這些邊界反射導致結(jié)構(gòu)中的聲波在穩(wěn)態(tài)條件下僅在某些頻率處存在較大的響應。圖中曲線以一種類似正弦或余弦的形式發(fā)展,因而,應變在某些頻率處達到極值,這些頻率就是熱聲系統(tǒng)的本征頻率。點d位于諧振管封閉端,由之前的應力云圖可知此處諧振管的應變達到系統(tǒng)峰值,尤其是z方向也就是軸向應變最為強烈。該點處各個方向的應變在5階、6階模態(tài)即300Hz附近同時達到最大值點。值得注意的是位于諧振管內(nèi)表面的e模態(tài)點,在3階、4中國煤化工變幅度最大x方向;2-y方向:3-z方向這與氣體振蕩日5b點處應變隨頻率的變化CNMH度的氣體介質(zhì)的振動在固體,xJ中介質(zhì)的結(jié)構(gòu)振Flg.5 Curves of strain versus frequency at point b動,結(jié)構(gòu)振動會反作用于氣體介質(zhì)。當固體介質(zhì)在某2848中南大學學報(自然科學版第43卷一模態(tài)中達到其本征頻率,且該頻率與聲激勵頻率相應力由里至外逐漸加強,而封閉端在圓板中心位置達近時,聲波與固體介質(zhì)會在結(jié)構(gòu)中產(chǎn)生以此階模態(tài)為到峰值并呈同心圓狀由內(nèi)及外逐漸減小:板疊處各個基礎(chǔ)的相互作用,從而產(chǎn)生共振和耦合41。平板前后兩端面幾乎沒有受到應力作用,上下端面發(fā)生了較大應變(3)隨著結(jié)構(gòu)頻率的上升,平板的應變變化幅度總體是不斷加強的,且各處受振動的影響相似,在垂1.0直于平板的y方向應變最大而其他2個方向在低頻段幾乎沒有發(fā)生應變,與實際情況相吻合。在頻率為5000.5時平板的應力應變達到極大值,且隨著頻率的進步提高而增大,說明回熱器激發(fā)和維持高階熱聲振蕩。(4)熱聲系統(tǒng)的3階和4階模態(tài)下固有頻率接近氣體介質(zhì)的振蕩頻率,根據(jù)檢測點d的結(jié)果驗證了高密度的氣體介質(zhì)的振動會在回熱器中產(chǎn)生振動激勵和1.0聲激勵并作用于諧振管,使聲波和固體介質(zhì)產(chǎn)生共振模態(tài)和耦合。x方向;2-y方向;3z方向圖7d點處應變隨頻率變化擬合曲線參考文獻:Fig7 Fitting curves of strain versus frequency at point d[〕張曉青,郭方中,熱聲熱機中回熱器的數(shù)值計算與分析團門0.8華中科技大學學報,2001,29(3):87-89ZHANG Xiao-qing, GUO Fangand analysis of regenerator in thermoacoustic heat engines[J]Journal of Huazhong University of Science and Technology2001,29(3)87-89[2] HU Zhong-jun,g. LI Qiang, et al. A high frequency04cascade thermogine ] Cryogenics, 2006, 46(11)771-7773]陳煕,李青,李正宇,等.熱聲熱機網(wǎng)絡模型演化與機理初探[低溫工程,2004(1)27-31CHEN Xi, LI Qing, LI Zheng-yu, et al. Network model of模態(tài)d study on the mechanism[J]1-x方向:2-y方向;3-z方向Cryogenics,2004(1):27-31.圖8e點處應變隨頻率變化擬合曲線[4]羅二倉,戴巍,RayR交變流動回熱器的熱聲功能和回熱功能[工程熱物理學報,2006,27(1)1-4Fig& Fitting curves of strain versus frequency at point eLUO Er-cang, DAI Wei, Ray R Thermoacousticrecuperative functions of a cyclic-flow regenerator]. Journal of3結(jié)論engineering thermophysics, 2006, 27(1): 1-4.5]劉益才,張明研,黃謙,等.熱聲熱機板疊式回熱器結(jié)構(gòu)數(shù)值計算門中南大學學報:自然科學版2010,41(3):186-119(1)利用商業(yè)軟件 ANSYS建立了熱聲熱機模型LIU Yi-cai, ZHANG Ming-yan, HUANG Qian, et al. Numerical并對其進行了結(jié)構(gòu)靜力、模態(tài)和諧響應分析。在分別施加了03,05和0.7MPa3種壓力載荷后得到了熱t(yī)hermoacoustic engine J]. Journal of Central South University聲系統(tǒng)內(nèi)部位移、應力和應變等參數(shù)的隨壓力和頻率Science and Technology, 2010, 41(3): 1186-1189的變化關(guān)系。隨著壓力的提升,應力應變基本呈現(xiàn)線DAYU Guo-yBo, LUO Er-cang, e al.300 Hz中國煤化工 T for temperature性增大的趨勢,且峰值點一一對應。CNMHG 902): 104-124(2)在板疊回熱器位置處,諧振管外部應力和應⑦謝秀嫻,李,禹廠m凹器聲學特性的理變極其微弱,而其他位置所受應力較大:諧振管徑向論及實驗研究低溫工程,20046)33-37劉益才,等」熱聲熱機結(jié)構(gòu)動力學分析ⅫEⅫ