余剑武,仝瑞庆†,罗红,陆岳托,胡俊之
(1.湖南大学国家高效磨削工程技术研究中心,湖南长沙 410082;2.湖南大学 机械与运载工程学院,湖南 长沙 410082)
目前超声振动辅助技术在机械加工、材料成形和材料焊接领域得到了十分广泛的应用[1-2].变幅杆作为超声振动系统中必不可少的一个部件,作用是将超声波能量从源元件(换能器或另一变幅杆)高效传输到工具头、耦合器、或负载元件,以达到振幅最大化[3].变幅杆与工具头之间多采用螺栓连接,但在超声频振动工况下,变幅杆与工具头之间的螺栓连接损耗[4]以及材料内部阻尼[5]等能量损耗形式均会导致变幅杆产生明显的温升,严重影响变幅杆振动性能的稳定性.
Rani 等[6]模拟并测量了超声波焊接中不同阻尼材料的变幅杆的温升,结果表明,温升大小主要取决于变幅杆材料的阻尼.Yu 等[7]研究了超声振动玻璃模压中高温环境对变幅杆振动特性的影响,发现高温致使变幅杆谐振频率降低,且谐振频率的降低主要与变幅杆弹性模量降低和热膨胀相关.崔方方等[8]研究了超声辅助磨削中不同材质砂轮基体的振动功耗与振幅,发现连续振动状态下砂轮基体的温升越大,其输出振幅就越小,且该温升与材料阻尼高度相关.
上述研究表明,变幅杆温度升高将引起其谐振频率降低,进而影响其输出特性.考虑到目前工业超声波发生器的频率追踪范围有限,变幅杆谐振频率降低过多时将造成整个超声振动系统失谐[9].因此,分析变幅杆的发热原因并减少其温升对谐振特性的影响是超声辅助加工和成形领域共同关注的问题.另一方面,课题组在前期研究中发现碳化钨模具摩擦系数对微结构光学玻璃元件模压成形精度有较大的影响[10-11],因此有必要进行模具抛光.考虑到碳化钨材料的难加工性和模具结构的复杂性,将采用超声抛光方法来降低模具表面粗糙度.这对变幅杆在连续工况下的振动稳定性提出了更高要求,其中控制变幅杆的发热是关键.
本文以模具抛光用钛合金超声纵振变幅杆为例,对变幅杆的振动自热现象进行数值模拟和实验测定,研究其能量耗散机理和温升特性,并分析不同螺栓材料和工具头材料对温升特性的影响.结合有限元分析和实验,进一步分析变幅杆温升对其谐振频率和振幅的影响,从而为连续工况下变幅杆的振动稳定性设计优化提供参考和指导.
对于由均匀、各向同性材料所构成的变截面杆,在发生自由振动的情况下,变幅杆的设计通常基于经典的一维波动方程[3]:
式中:S 为变幅杆截面积;函数ξ 为质点位移函数;k=ω/c 为圆波数,c 为纵波在变幅杆中的传播速度.
考虑到变幅杆内部谐振损耗机制,纵振波动方程可表示为[12]:
式中:Rm为材料内部损耗因子;cd为考虑阻尼作用时均匀等截面细杆中的波速.cd可表示为:
式中:E 为材料弹性模量;ρ 为材料密度.由于Rm<<ρ,可知阻尼对波速的影响较小,在变幅杆的理论设计中可以忽略.
因此,基于半波长设计方法,可推导出变幅杆的长度l 与振动频率f 及波长λ 的关系:
变幅杆结构设计为典型的两段式,在第1 段l1大端处连接换能器,在节点位置设置法兰;第2 段l2为工具头,两段之间采用M8×1×20 的螺栓连接.变幅杆材料选择TC4 钛合金,中心工作频率f 设定为35 kHz.结合TC4 钛合金材料参数可确定l1=l2=71.28 mm.初步选定螺栓材料为45 钢,得到变幅杆的几何模型如图1(a)所示.
通过有限元分析软件对变幅杆进行模态分析,搜索35 kHz 附近的谐振频率,得到其纵振模态如图1(b)所示.该纵振频率为34 976 Hz,与设定频率35 kHz 接近.初步验证了理论设计的准确性.
图1 变幅杆几何模型与纵振模态Fig.1 Geometrical model and longitudinal mode of horn
为分析变幅杆的振动自热特性,对变幅杆进行热力耦合数值模拟,并采用红外热成像仪对不同材料螺栓连接的变幅杆以及不同阻尼材料的工具头进行温升监测.
变幅杆振动自热模拟基于结构频域分析与传热分析的耦合,考虑了材料与螺栓处的阻尼,模拟变幅杆在35 kHz 振动时的温升与能量损耗.
在频域分析中,需要为材料添加各向同性内部损耗因子.Mason 等[13]测得退火TC4 钛合金在低应变条件下的品质因子Qm=20 000,则各向同性内部损耗因子取1/Qm=5×10-5.相对于钛合金材料的自身阻尼,当变幅杆与工具头之间采用螺栓连接时,螺栓处阻尼水平要大得多,可等效为等效黏性阻尼[14],其阻尼水平多处于0.01 数量级.
在法兰处设定固定约束,变幅杆大端处承受周期位移激励载荷.先要测量实际工作中变幅杆的激励振幅,但在换能器与变幅杆装配后,不方便直接测量换能器的输出振幅,即变幅杆的激励振幅,因此,首先使用KEYENCE LK-H020 激光位移传感器测量小端端面的振幅为13.67 μm,通过谐响应分析获得变幅杆轴线振幅分布曲线,如图2 所示.由于放大系数为3.19,因此设置变幅杆大端端面的激励振幅为4.28 μm.
在传热分析模块中,温度升高由传热方程给出:
式中:k 为导热系数;ρCp为体积比热,与温度无关,符合Dulong-Petit 定律;T 为时间段2π/ω 的平均温度;Qh为热源,是非弹性力在这段时间内的内部功,在频域结构分析中计算得出,并在结构分析与传热分析之间传递.Qh可表示为[6]:
式中:η 为损耗因子;ω 为角激励频率;ε 和C 分别为应变张量和弹性矩阵.
图2 变幅杆轴线振幅分布曲线Fig.2 Axial amplitude distribution curve of horn
对于其他边界使用对流冷却条件:
式中:T 为温度;h 为传热系数;Text为远离边界的外部流体的温度.在仿真过程中,节点固定端和周围环境保持20 ℃的恒定温度,变幅杆在空气中自然对流换热,传热系数取h=20 W/(m2·K),外部温度Text=20 ℃.
为深入了解变幅杆发热的主要原因,对变幅杆温度与总功耗密度进行分析,其中总功耗密度为单位体积内的功率损耗.在30 min 时变幅杆轴线温度与总功耗密度曲线见图3,可知总功耗密度在两个螺栓端面出现峰值,达到3.88×107W/m3,远大于由变幅杆材料自身阻尼引起的最大总功耗密度,而位于变幅杆的1/4 波长附近,只有1.59×105W/m3.分析螺栓处功耗与整体功耗的比值,得到螺栓处发热量占变幅杆发热总量的97.7%,因而在螺栓处有最高温度117.14 ℃.值得注意的是,变幅杆大端端面的温升几乎可以忽略,这对保持换能器的热稳定性起着良好作用.
为了分析变幅杆的温升变化规律,分别取螺栓连接处与小端端面一点,其温度随时间变化规律见图4.由图4 可知,螺栓处温度在开始振动时急剧上升,然后温度升速减缓,并在20 min 左右达到平衡,螺栓处温度稳定在117.14 ℃,小端端面处温度稳定在64.91 ℃.
图3 变幅杆轴线温度与总功耗密度分布曲线Fig.3 Distribution curve of axial temperature and total power consumption density of horn
图4 螺栓连接处与小端端面处的温度变化趋势Fig.4 Trend of temperature change at bolted joint and the end of tool head
由仿真结果可知,螺栓处能量损耗为发热首要原因.为对比分析不同螺栓材料对变幅杆发热的影响,使用工程中常用的45 钢螺栓以及与变幅杆相同材料的TC4 钛合金螺栓分别连接变幅杆与工具头,并利用FLIR SC7000 红外热成像仪对变幅杆温度进行测量.测量用变幅杆成品以及测量区域如图5(a)所示,测量区域为法兰以下部分,获得其温度分布分别如图5(b)和5(c)所示,绘制变幅杆的轴线温度曲线见图6.由图5 和图6 可知,最高温度均分布于螺栓处,变幅杆轴线温度变化规律与仿真相似.使用45钢螺栓连接时,最高温度达到了118.52 ℃;而使用TC4 钛合金螺栓连接时,最高温度明显降低,只有84.85 ℃.在采用45 钢螺栓时,45 钢的阻尼较TC4 钛合金更大[15],并且机械波从钛合金材料入射到45 钢时,会有一部分超声能量反射,变成热能损耗在器件中,因而温升更快.螺栓材料的选择在变幅杆实际应用中容易被忽视.在满足加工强度的情况下,应选用与变幅杆材料相同的材料.
变幅杆材料的自身阻尼也是发热不可或缺的一个原因.为了分析不同阻尼的材料对变幅杆最终温度的影响,分别对高阻尼TU2 紫铜工具头以及钛合金工具头进行测量.测试时连接螺栓为TC4 钛合金螺栓,变幅杆温度分布与轴线温度分布曲线分别见图5(d)和图6.与低阻尼材料不同的是,使用高阻尼TU2 紫铜工具头时的温升首先遍布在整个工具头与螺栓处,工具头温度均匀分布,最终达到141.99 ℃,比TC4 钛合金工具头的温升高出57.14 ℃.
图5 变幅杆测量区域与最终温度分布Fig.5 Measuring area and final temperature of horn
图6 不同材料的螺栓/工具头连接时幅杆轴线温度分布Fig.6 Temperature distribution of the axis of the horn under connecting bolt and tool head with different materials
通过对仿真值与实验测量值的分析可以发现,温度上升最终会达到某一平衡点,如图4 与图7 所示.不过对于高阻尼材料的工具头,温度上升更快,在5 min 便达到最高温度.
温度达到平衡的主要原因是:由公式(7)可知,随变幅杆温度上升,变幅杆与周围环境的热量交换不断增大;另一方面,由文献[14]可知,能量损耗与振幅成正比,温度上升导致振幅减小,进而致使能量损耗降低.两者共同作用使温度达到平衡,这对变幅杆的振动稳定性是有利的.
图7 不同材料的螺栓/工具头连接时螺栓连接处的温度变化趋势Fig.7 Temperature change trend at bolted joint under connecting bolt and tool head with different materials
在仿真后处理中取变幅杆小端端面一点,绘制振幅随时间的变化曲线,如图8 所示.使用激光位移传感器分别测得初始时刻与振动30 min 时的振幅如图9 所示.仿真结果显示振幅从13.67 μm 下降到11.85 μm,实验结果显示振幅从13.67 μm 下降到10.9 μm,从图8 和图9 可知,振动发热后的振幅均比初始振幅要小.同时由图8 可以得出,振幅的降速与温升速度相关.
图8 小端端面处位移随时间变化曲线Fig.8 Displacement curve with time at the end of tool head
振幅下降的原因主要是:随温度的升高,变幅杆的谐振频率不断下降,变幅杆大端处的激励频率35 000 Hz 将不断偏离变幅杆的谐振频率.在共振的情况下振幅最大,在变幅杆的谐振频率不断偏离共振频率之后,振幅将有所减小.
图9 变幅杆振幅测试Fig.9 Amplitude test of horn
对于1/2 波长变幅杆,当温度变化时,假设变幅杆温度均匀分布,其纵向谐振频率满足:
式中:L0为变幅杆初始长度;α 为热膨胀系数;E0为初始弹性模量;β 为弹性模量的相对变化;ρ0为初始密度;r 为工具头的半径.
由式(8)可得,漂移后谐振频率f 与初始谐振频率f0的比值为:
在变幅杆的实际设计中,为避免横向振动,变幅杆中直径小于1/4 波长,同时又由于系数α 与系数β的值非常小,可以近似地认为式(9)中忽略(1+πr2αΔT)项,其一阶泰勒级数展开是准确的,则
由式(10)可知,该比值始终小于1,即频率漂移与温升呈线性负相关,随温度的升高,变幅杆的频率漂移将不断增加.此外,还与材料热膨胀系数和弹性模量的相对变化有关.
对于TC4 钛合金材料,取α=9.0×10-6m/(m·℃),β=4.602×10-4K-1[16-17],结合公式(10)可知,温升下的弹性模量相对变化为变幅杆谐振频率降低的主要原因.对于初始谐振频率为34 976 Hz 变幅杆,在T=120 ℃时,即由室温20 ℃升温100 ℃时,可得到f/f0=0.976 1,其谐振频率变化至34 140 Hz,频率下降了836 Hz.
将振动发热仿真中瞬态分析的结果,即变幅杆的温度分布作为模态分析的初始值,并考虑到不同温度下的材料参数变化,对变幅杆进行振动发热后模态分析,得到振动发热前后变幅杆谐振频率值,如表1 所示.
表1 振动发热前后变幅杆谐振频率值Tab.1 Resonant frequency values of horns before and after vibration heating
频率漂移的仿真值要比其计算值小,主要原因是变幅杆的温度分布不均,在轴向有较大的温度梯度.
本文结合数值模拟和实验验证对变幅杆振动发热问题进行了分析,并研究了振动发热对超声振动稳定性的影响,得到以下结论:
1)变幅杆材料与连接螺栓材料对温升有较大影响,而最终温度会达到平衡.在变幅杆阻尼较小时,以螺栓处发热为主,45 钢螺栓处发热量占变幅杆发热总量的97.7%,稳定后温度达到118.52°C,大于TC4 钛合金螺栓连接时的84.85°C.采用较高阻尼的TU2 紫铜工具头时,温升更快,达到141.99°C,并且高温存在于螺栓与整个工具头区域.
2)变幅杆振幅与谐振频率均随温度的升高而衰减,并最终保持稳定.振幅的降速与温度的升速相关,并且当变幅杆温度均匀分布时,谐振频率与温升呈线性负相关.
3)变幅杆材料宜采用低阻尼、高强度、热稳定性好的材料,螺栓材料应与变幅杆相同,在长时间工作下有必要采取通风等冷却方式以减小其频率漂移.