余劍武 仝瑞慶 羅紅 陸岳托 胡俊之

摘? ?要：在變幅桿振動中，變幅桿的自身阻尼特性以及與螺栓/工具頭的耦合會導致明顯溫升. 為研究該溫升特性并探討其對變幅桿振動穩(wěn)定性的影響，采用數(shù)值模擬得到了縱振變幅桿的總功耗密度和溫度分布;通過紅外熱成像測試對變幅桿理論溫升進行了驗證，并分析了不同材料的螺栓與工具頭連接對溫升的影響. 結(jié)合數(shù)值模態(tài)分析與實驗分析，進一步探討了溫升對變幅桿諧振頻率與振幅的影響規(guī)律. 結(jié)果表明，變幅桿連續(xù)運行時溫升趨于平衡，該平衡點取決于螺栓材料、變幅桿材料和環(huán)境因素. 對于低阻尼TC4鈦合金變幅桿，螺栓處能量損耗為主要熱源，其中45鋼螺栓發(fā)熱量占97.7%，而采用TC4鈦合金螺栓和TU2紫銅工具頭可顯著減小溫升. 同時發(fā)現(xiàn)，變幅桿諧振頻率與溫升成線性負相關(guān)，振幅因溫升有所降低，二者均穩(wěn)定在平衡溫度點.

關(guān)鍵詞：超聲變幅桿;能量損耗;諧振頻率;振幅;振動發(fā)熱

Abstract：In the vibration of the horn， the damping characteristics of the ultrasonic horn and its coupling with the bolt / tool head will lead to obvious temperature rise. In order to analyze the temperature rise characteristic and explore its influence on the vibration stability of the ultrasonic horn， the total power density and temperature distribution of the longitudinal vibration horn were obtained by numerical simulation. The theoretical temperature rise of the ultrasonic horn was verified by infrared thermal imaging test， and the influence of bolts and tool heads of different materials on the temperature rise was analyzed. Combining the numerical modal analysis and experimental analysis， the influence of temperature rise on the resonant frequency and amplitude of the horn was further discussed. The results show that the temperature rise tends to be balanced during continuous operation of the ultrasonic horn， and the equilibrium point depends on the bolt material， horn material and environmental factors. For low damping TC4 horns， energy loss at bolt is the main heat source， of which 45 steel bolts generate 97.7% of the heat， while TC4 bolts and tool heads can significantly reduce the temperature rise. At the same time， it is found that the resonance frequency of the ultrasonic horn is linearly negatively correlated with the temperature rise. The amplitude is reduced due to the temperature rise， and both of them are stable at the equilibrium temperature point.

Key words：ultrasonic horn;energy loss;resonance frequency;amplitude;vibration heating

目前超聲振動輔助技術(shù)在機械加工、材料成形和材料焊接領(lǐng)域得到了十分廣泛的應(yīng)用[1-2]. 變幅桿作為超聲振動系統(tǒng)中必不可少的一個部件，作用是將超聲波能量從源元件（換能器或另一變幅桿）高效傳輸?shù)焦ぞ哳^、耦合器、或負載元件，以達到振幅最大化[3]. 變幅桿與工具頭之間多采用螺栓連接，但在超聲頻振動工況下，變幅桿與工具頭之間的螺栓連接損耗[4]以及材料內(nèi)部阻尼[5]等能量損耗形式均會導致變幅桿產(chǎn)生明顯的溫升，嚴重影響變幅桿振動性能的穩(wěn)定性.

Rani等[6]模擬并測量了超聲波焊接中不同阻尼材料的變幅桿的溫升，結(jié)果表明，溫升大小主要取決于變幅桿材料的阻尼. Yu等[7]研究了超聲振動玻璃模壓中高溫環(huán)境對變幅桿振動特性的影響，發(fā)現(xiàn)高溫致使變幅桿諧振頻率降低，且諧振頻率的降低主要與變幅桿彈性模量降低和熱膨脹相關(guān). 崔方方等[8]研究了超聲輔助磨削中不同材質(zhì)砂輪基體的振動功耗與振幅，發(fā)現(xiàn)連續(xù)振動狀態(tài)下砂輪基體的溫升越大，其輸出振幅就越小，且該溫升與材料阻尼高度相關(guān).

上述研究表明，變幅桿溫度升高將引起其諧振頻率降低，進而影響其輸出特性. 考慮到目前工業(yè)超聲波發(fā)生器的頻率追蹤范圍有限，變幅桿諧振頻率降低過多時將造成整個超聲振動系統(tǒng)失諧[9].? 因此，分析變幅桿的發(fā)熱原因并減少其溫升對諧振特性的影響是超聲輔助加工和成形領(lǐng)域共同關(guān)注的問題. 另一方面，課題組在前期研究中發(fā)現(xiàn)碳化鎢模具摩擦系數(shù)對微結(jié)構(gòu)光學玻璃元件模壓成形精度有較大的影響[10-11]，因此有必要進行模具拋光. 考慮到碳化鎢材料的難加工性和模具結(jié)構(gòu)的復雜性，將采用超聲拋光方法來降低模具表面粗糙度. 這對變幅桿在連續(xù)工況下的振動穩(wěn)定性提出了更高要求，其中控制變幅桿的發(fā)熱是關(guān)鍵.

本文以模具拋光用鈦合金超聲縱振變幅桿為例，對變幅桿的振動自熱現(xiàn)象進行數(shù)值模擬和實驗測定，研究其能量耗散機理和溫升特性，并分析不同螺栓材料和工具頭材料對溫升特性的影響. 結(jié)合有限元分析和實驗，進一步分析變幅桿溫升對其諧振頻率和振幅的影響，從而為連續(xù)工況下變幅桿的振動穩(wěn)定性設(shè)計優(yōu)化提供參考和指導.

1? ?TC4鈦合金變幅桿設(shè)計

變幅桿結(jié)構(gòu)設(shè)計為典型的兩段式，在第1段l1大端處連接換能器，在節(jié)點位置設(shè)置法蘭;第2段l2為工具頭，兩段之間采用M8×1×20的螺栓連接. 變幅桿材料選擇TC4鈦合金，中心工作頻率f設(shè)定為35 kHz. 結(jié)合TC4鈦合金材料參數(shù)可確定l1 = l2 = 71.28 mm. 初步選定螺栓材料為45鋼，得到變幅桿的幾何模型如圖1（a）所示.

通過有限元分析軟件對變幅桿進行模態(tài)分析，搜索35 kHz附近的諧振頻率，得到其縱振模態(tài)如圖1（b）所示. 該縱振頻率為34 976 Hz，與設(shè)定頻率35 kHz接近. 初步驗證了理論設(shè)計的準確性.

2? ?變幅桿振動自熱分析

為分析變幅桿的振動自熱特性，對變幅桿進行熱力耦合數(shù)值模擬，并采用紅外熱成像儀對不同材料螺栓連接的變幅桿以及不同阻尼材料的工具頭進行溫升監(jiān)測.

2.1? ?振動自熱仿真模型

變幅桿振動自熱模擬基于結(jié)構(gòu)頻域分析與傳熱分析的耦合，考慮了材料與螺栓處的阻尼，模擬變幅桿在35 kHz振動時的溫升與能量損耗.

在頻域分析中，需要為材料添加各向同性內(nèi)部損耗因子. Mason等[13]測得退火TC4鈦合金在低應(yīng)變條件下的品質(zhì)因子Qm = 20 000，則各向同性內(nèi)部損耗因子取1/Qm = 5 × 10-5. 相對于鈦合金材料的自身阻尼，當變幅桿與工具頭之間采用螺栓連接時，螺栓處阻尼水平要大得多，可等效為等效黏性阻尼[14]，其阻尼水平多處于0.01數(shù)量級.

在法蘭處設(shè)定固定約束，變幅桿大端處承受周期位移激勵載荷. 先要測量實際工作中變幅桿的激勵振幅，但在換能器與變幅桿裝配后，不方便直接測量換能器的輸出振幅，即變幅桿的激勵振幅，因此，首先使用 KEYENCE LK-H020 激光位移傳感器測量小端端面的振幅為13.67 μm，通過諧響應(yīng)分析獲得變幅桿軸線振幅分布曲線，如圖2所示. 由于放大系數(shù)為3.19，因此設(shè)置變幅桿大端端面的激勵振幅為4.28 μm.

2.2? ?仿真結(jié)果分析

為深入了解變幅桿發(fā)熱的主要原因，對變幅桿溫度與總功耗密度進行分析，其中總功耗密度為單位體積內(nèi)的功率損耗. 在30 min時變幅桿軸線溫度與總功耗密度曲線見圖3，可知總功耗密度在兩個螺栓端面出現(xiàn)峰值，達到3.88×107 W/m3，遠大于由變幅桿材料自身阻尼引起的最大總功耗密度，而位于變幅桿的1/4波長附近，只有1.59×105 W/m3. 分析螺栓處功耗與整體功耗的比值，得到螺栓處發(fā)熱量占變幅桿發(fā)熱總量的97.7%，因而在螺栓處有最高溫度117.14 ℃. 值得注意的是，變幅桿大端端面的溫升幾乎可以忽略，這對保持換能器的熱穩(wěn)定性起著良好作用.

為了分析變幅桿的溫升變化規(guī)律，分別取螺栓連接處與小端端面一點，其溫度隨時間變化規(guī)律見圖4. 由圖4可知，螺栓處溫度在開始振動時急劇上升，然后溫度升速減緩，并在20 min左右達到平衡，螺栓處溫度穩(wěn)定在117.14 ℃，小端端面處溫度穩(wěn)定在64.91 ℃.

2.3? ?振動溫升實驗分析

由仿真結(jié)果可知，螺栓處能量損耗為發(fā)熱首要原因. 為對比分析不同螺栓材料對變幅桿發(fā)熱的影響，使用工程中常用的45鋼螺栓以及與變幅桿相同材料的TC4鈦合金螺栓分別連接變幅桿與工具頭，并利用 FLIR SC7000紅外熱成像儀對變幅桿溫度進行測量. 測量用變幅桿成品以及測量區(qū)域如圖5（a）所示，測量區(qū)域為法蘭以下部分，獲得其溫度分布分別如圖5（b）和5（c）所示，繪制變幅桿的軸線溫度曲線見圖6. 由圖5和圖6可知，最高溫度均分布于螺栓處，變幅桿軸線溫度變化規(guī)律與仿真相似. 使用45鋼螺栓連接時，最高溫度達到了118.52 ℃;而使用TC4鈦合金螺栓連接時，最高溫度明顯降低，只有84.85 ℃. 在采用45鋼螺栓時，45鋼的阻尼較TC4鈦合金更大[15]，并且機械波從鈦合金材料入射到45鋼時，會有一部分超聲能量反射，變成熱能損耗在器件中，因而溫升更快. 螺栓材料的選擇在變幅桿實際應(yīng)用中容易被忽視. 在滿足加工強度的情況下，應(yīng)選用與變幅桿材料相同的材料.

變幅桿材料的自身阻尼也是發(fā)熱不可或缺的一個原因. 為了分析不同阻尼的材料對變幅桿最終溫度的影響，分別對高阻尼TU2紫銅工具頭以及鈦合金工具頭進行測量. 測試時連接螺栓為TC4鈦合金螺栓，變幅桿溫度分布與軸線溫度分布曲線分別見圖5（d）和圖6. 與低阻尼材料不同的是，使用高阻尼TU2紫銅工具頭時的溫升首先遍布在整個工具頭與螺栓處，工具頭溫度均勻分布，最終達到141.99 ℃，比TC4鈦合金工具頭的溫升高出57.14 ℃.

通過對仿真值與實驗測量值的分析可以發(fā)現(xiàn)，溫度上升最終會達到某一平衡點，如圖4與圖7所示. 不過對于高阻尼材料的工具頭，溫度上升更快，在5 min便達到最高溫度.

溫度達到平衡的主要原因是：由公式 （7） 可知，隨變幅桿溫度上升，變幅桿與周圍環(huán)境的熱量交換不斷增大;另一方面，由文獻[14]可知，能量損耗與振幅成正比，溫度上升導致振幅減小，進而致使能量損耗降低. 兩者共同作用使溫度達到平衡，這對變幅桿的振動穩(wěn)定性是有利的.

3? ?溫升對變幅桿諧振特性的影響

3.1? ?溫升對變幅桿振幅的影響

在仿真后處理中取變幅桿小端端面一點，繪制振幅隨時間的變化曲線，如圖8所示. 使用激光位移傳感器分別測得初始時刻與振動30 min時的振幅如圖9所示. 仿真結(jié)果顯示振幅從13.67 μm下降到11.85 μm，實驗結(jié)果顯示振幅從13.67 μm下降到10.9 μm，從圖8和圖9可知，振動發(fā)熱后的振幅均比初始振幅要小. 同時由圖8可以得出，振幅的降速與溫升速度相關(guān).

振幅下降的原因主要是：隨溫度的升高，變幅桿的諧振頻率不斷下降，變幅桿大端處的激勵頻率35 000 Hz將不斷偏離變幅桿的諧振頻率. 在共振的情況下振幅最大，在變幅桿的諧振頻率不斷偏離共振頻率之后，振幅將有所減小.

3.2? ?溫升對諧振頻率的影響

對于1/2波長變幅桿，當溫度變化時，假設(shè)變幅桿溫度均勻分布，其縱向諧振頻率滿足：

由式（10）可知，該比值始終小于1，即頻率漂移與溫升呈線性負相關(guān)，隨溫度的升高，變幅桿的頻率漂移將不斷增加. 此外，還與材料熱膨脹系數(shù)和彈性模量的相對變化有關(guān).

對于 TC4鈦合金材料，取α=9.0×10-6 m/（m·℃），β=4.602×10-4? K-1[16-17]，結(jié)合公式 （10） 可知，溫升下的彈性模量相對變化為變幅桿諧振頻率降低的主要原因. 對于初始諧振頻率為34 976 Hz變幅桿，在T=120 ℃ 時，即由室溫20 ℃升溫100 ℃時，可得到f/f0=0.976 1，其諧振頻率變化至34 140 Hz，頻率下降了836 Hz.

將振動發(fā)熱仿真中瞬態(tài)分析的結(jié)果，即變幅桿的溫度分布作為模態(tài)分析的初始值，并考慮到不同溫度下的材料參數(shù)變化，對變幅桿進行振動發(fā)熱后模態(tài)分析，得到振動發(fā)熱前后變幅桿諧振頻率值，如表1所示.

4? ?結(jié)? ?論

本文結(jié)合數(shù)值模擬和實驗驗證對變幅桿振動發(fā)熱問題進行了分析，并研究了振動發(fā)熱對超聲振動穩(wěn)定性的影響，得到以下結(jié)論：

1）變幅桿材料與連接螺栓材料對溫升有較大影響，而最終溫度會達到平衡. 在變幅桿阻尼較小時，以螺栓處發(fā)熱為主，45鋼螺栓處發(fā)熱量占變幅桿發(fā)熱總量的97.7%，穩(wěn)定后溫度達到118.52 °C，大于TC4鈦合金螺栓連接時的84.85 °C. 采用較高阻尼的TU2紫銅工具頭時，溫升更快，達到141.99 °C，并且高溫存在于螺栓與整個工具頭區(qū)域.

2）變幅桿振幅與諧振頻率均隨溫度的升高而衰減，并最終保持穩(wěn)定. 振幅的降速與溫度的升速相關(guān)，并且當變幅桿溫度均勻分布時，諧振頻率與溫升呈線性負相關(guān).

3）變幅桿材料宜采用低阻尼、高強度、熱穩(wěn)定性好的材料，螺栓材料應(yīng)與變幅桿相同，在長時間工作下有必要采取通風等冷卻方式以減小其頻率漂移.

參考文獻

[1]? ? 王愛玲，祝錫晶，吳秀玲. 功率超聲振動加工技術(shù)[M]. 北京：國防工業(yè)出版社，2007：258—293.

WANG A L，ZHU X J，WU X L. Power ultrasonic vibration machining technology[M]. Beijing：National Defense Industry Press，2007：258—293. （In Chinese）

[2]? ? 賀地求，趙志峰，賴瑞林，等. 2219-T87超聲輔助攪拌摩擦焊接頭組織與性能[J]. 湖南大學學報（自然科學版），2018，45（4）：41—47.

HE D Q，ZHAO Z F，LAI R L，et al. Research on the microstructure and mechanical properties of ultrasonic assisted friction stir welding joints of 2219-T87 aluminum alloy[J]. Journal of Hunan University （Natural Sciences），2018，45（4）：41—47.（In Chinese）

[3]? ? 林仲茂. 超聲變幅桿的原理和設(shè)計[M]. 北京：科學出版社，1987：65—66.? ? ? LIN Z M. Principle and design of ultrasonic horn[M]. Beijing：Science Press，1987：65—66. （In Chinese）

[4]? ? 傅俊慶，榮見華，張玉萍，等. 螺栓連接接口軸向振動能量耗散特性研究[J]. 振動、測試與診斷，2005，25（3）：42—46.

FU J Q，RONG J H，ZHANG Y P，et al. Study on energy dissipation of axial vibration in bolted joint[J]. Journal of Vibration，Measurement & Diagnosis，2005，25（3）：42—46.（In Chinese）

[5]? ? 鄭冀. 材料物理性能[M]. 天津：天津大學出版社，2008：236—238.ZHENG J. Physical properties of materials[M]. Tianjin：Tianjin University Press，2008：236—238. （In Chinese）

[6]? ? RANI R M，PRAKASAN K，RUDRAMOORTHY R. Studies on thermoelastic heating of horns used in ultrasonic plastic welding[J]. Ultrasonics，2015，55：123—132.

[7]? ? YU J W，LUO H，NGUYEN T V，et al. Eigenfrequency characterization and tuning of Ti-6Al-4V ultrasonic horn at high temperatures for glass molding[J]. Ultrasonics，2020，101：106002.

[8]? ? 崔方方，丁凱，李奇林，等. 超聲輔助磨削砂輪基體材料的振動性能研究[J]. 工具技術(shù)，2017，51（7）：20—23.

CUI F F，DING K，LI Q L，et al. Study on vibration performance of wheel matrix materials for ultrasonic assisted grinding[J]. Tool Engineering，2017，51（7）：20—23.（In Chinese）

[9]? ? 余劍武，柳波，羅紅，等. 玻璃振動模壓中變幅桿高溫模態(tài)分析及實驗研究[J]. 聲學技術(shù)，2018，37（1）：51—56.

YU J W，LIU B，LUO H，et al. High temperature modal analysis and experimental study of the ultrasonic horn for glass vibration molding[J]. Technical Acoustics，2018，37（1）：51—56. （In Chinese）

[10]? 謝望清. 微結(jié)構(gòu)光學元件模壓成形數(shù)值仿真分析[D].長沙：湖南大學機械與運載工程學院，2013：35—36. XIE W Q. The simulation analysis on glass molding process of microstructured optical components[D]. Changsha：College of Mechanical and Vehicle Engineering，Hunan University，2013：35—36. （In Chinese）

[11]? 余劍武，李嬋，尹韶輝，等. 微V形槽玻璃元件模壓成形有限元應(yīng)力分析[J]. 湖南大學學報（自然科學版），2017，44（8）：13—19.YU J W，LI C，YIN S H，et al. Finite element analysis on stress of micro V groove components in GMP[J]. Journal of Hunan University（Natural Sciences），2017，44（8）：13—19. （In Chinese）

[12]? ENSMINGER D，STULEN F B. Ultrasonics：data，equations and their practical uses[M]. Boca Raton：CRC Press，2008：61—62.

[13]? MASON W P，WEHR J. Internal friction and ultrasonic yield? stress of the alloy 90Ti6Al4V[J]. Journal of Physics & Chemistry of Solids，1970，31（8）：1925—1933.

[14]? SINGIRESU S R. 機械振動[M]. 李欣業(yè)，張明路，譯. 北京：清華大學出版社，2007：158—159，495—496.SINGIRESU S R. Mechanical vibration[M]. Translated by LI X Y，ZHANG M L. Beijing：Tsinghua University Press，2007：158—159，495—496. （In Chinese）

[15]? ZHANG J，PEREZ R J，LAVERNIA E J. Documentation of damping capacity of metallic，ceramic and metal matrix composite materials[J]. Journal of Materials Science，1993，28（9）：2395—2404.

[16]? WELSCH G，BOYER R，COLLINGS E W. Materials properties handbook：titanium alloys[M]. Ohio：ASM International，1993：641—649.

[17] MILLS K C. Ti：Ti-6Al-4V（Imi 318），Recommended values of thermophysical properties for selected commercial alloys[M]. Cambridge：Woodhead Publishing Ltd，2002：211—217.