1. 引言
模态分析常用来计算结构振动特性的数值 [1]。通过对双螺杆膨胀机主箱体进行模态分析,可以得到结构的固有频率和危险振动范围,从而可以避免结构在特定频率下的共振或振动 [2],以提高机械加工中的精度。
目前国内外学者采用模态分析的方法对如何优化机器振动进行了大量研究。李红华对双螺杆挤出机传动箱进行有限元模态分析,根据这些模态参数对箱体性能的影响,指出了传动箱体结构改进的方向 [3]。赵海涛 [1] 通过对主轴箱分别进行热模态和振动模态分析,处理结果提出了两种模态分析不同的优化设计方向。刘碧俊通过对YND100减速器箱体进行有限元模态分析,为减速器在工作过程中的共振及噪音产生分析提供了数据支持 [4]。
本文主要通过ANSYS Workbench针对双螺杆膨胀机的重要零件——主箱体进行模态分析,根据分析得出的数据和结果,对主箱体的优化提出相关方案,以期为螺杆膨胀机的相关研究提供参考。
2. 模态分析理论基础
无阻尼模态分析在特征值问题中很典型,其动力学运动方程为 [5]:
(1)
其中M为质量;K为刚度;
为加速度;
为位移。
其解的形式为:
(2)
代入式(1)中得:
(3)
式(3)中:特征值为
,自振圆频率为
。自振频率为
。
3. 双螺杆膨胀机主箱体模型建立
双螺杆膨胀机主箱体结构如图1所示。外壳材料选用45钢,其材料特性为:密度为7890 kg/m3,弹性模量为2.09 × 1011 Pa,泊松比为0.31。
本文利用SolidWorks软件对双螺杆膨胀机主箱体进行三维建模后,再导入到ANSYS Workbench软件中对其进行网格划分 [6]。划分网格时为方便计算和保证结果准确性,可使用四面体网格划分法,且在关键区域可使用曲度和近似尺寸功能自动细化网格,网格密度应适中 [7]。经过网格划分,共得到69046个单元。网格划分结果如图1所示。
![](//html.hanspub.org/file/8-2340652x17_hanspub.png)
Figure 1. Grid division diagram of main box
图1. 主箱体网格划分图
4. 双螺杆膨胀机主箱体模态分析
4.1. 约束条件
根据螺杆膨胀机实际的工作过程及原理,将四个支撑体的下平面进行固定约束,作为整个模型的边界条件。具体如图2所示:
![](//html.hanspub.org/file/8-2340652x18_hanspub.png)
Figure 2. Apply fixed constraints schematic
图2. 施加固定约束示意图
4.2. 求解结果
结构的振动可以用多阶固有振型的线性组合来表示,其中低阶固有振型会对结构振动产生较大影响,振型阶数越低就越会影响结构的振动,因此低阶振型对结构的动态特性起着决定性的作用 [8]。所以在分析时常选取固有振型的前4阶。
经过模态分析得到主箱体的前四阶固有频率如表1所示:
![](Images/Table_Tmp.jpg)
Table 1. The first four natural frequencies
表1. 前四阶固有频率
由表1可以看出,双螺杆膨胀机主箱体的固有频率分布在739.6~1443.7 Hz之间,频率分布较为集中,其中三、四阶固有频率较接近,该现象和箱体结构有一定关系。
主箱体前四阶振动形态 [9] 如图3~6所示,在上述图形中,坐标系XY平面平行于螺杆径向,Z轴正方向平行与螺杆轴向并指向外。
(1) 第一阶振型的振动形态:主箱体主要是沿X向振动,振幅较大;Y轴振动幅度较小;Z轴基本不变。中部壳体发生变形,上下出入口变形不明显,四个支撑基本不动。变形主要集中在壳体上表面部分。最大变形量为0.23877 mm,在四个振型中,一阶振型的变形量最小。
![](//html.hanspub.org/file/8-2340652x19_hanspub.png)
Figure 3. First order mode diagram of main box
图3. 主箱体一阶振型图
![](//html.hanspub.org/file/8-2340652x20_hanspub.png)
Figure 4. Second order mode diagram of main box
图4. 主箱体二阶振型图
![](//html.hanspub.org/file/8-2340652x21_hanspub.png)
Figure 5. Third order mode diagram of main box
图5. 主箱体三阶振型图
![](//html.hanspub.org/file/8-2340652x22_hanspub.png)
Figure 6. Fourth order mode diagram of main box
图6. 主箱体四阶振型图
(2) 第二阶振型的振动形态:主箱体主要是沿Z向振动,振幅较大;Y轴振动幅度较小;X轴基本不变。整体变形较小,四个支撑体变形较小。变形主要集中在出入口的上下表面。最大变形量为0.24693 mm,比一阶振型的变形量要略大。
(3) 第三阶振型的振动形态:主箱体主要是沿X和Z向振动,振幅较大;Y轴振动幅度较小。中部壳体和上下出入口变形较大,四个支撑体变形较小。变形主要集中在壳体外表面的两端面。最大变形量为0.32408 mm,是四个振型图中变形量最大的。
(4) 第四阶振型的振动形态:主箱体主要是沿X和Y向振动,振幅较大;Z轴振动幅度较小。中部壳体和上下出入口产生明显变形。四个支撑体基本不动。变形主要集中在出入口的外表面。最大变形量为0.27348 mm,比第一阶、二阶振型的变形量都要大,小于第三阶振型的变形量。
四个支撑体的使用增加了箱体的强度,施加固定约束的位置是支撑体的下表面,这样施加固定约束可降低该区域的振幅,从而使结构趋于稳定,有利于工作的稳定性 [10]。
根据公式:
(4)
得出不同固有频率下的螺杆膨胀机转速n在44374~86622 r/min间变化,因在实际工作中考虑到装置的安全性和能源的有限性,输入的蒸汽压力不会太高 [11],所以工作转速不可能接近上述数值,即产生共振的可能性非常小。
5. 结论
通过对双螺杆膨胀机主箱体进行模态分析,得到了其前四阶振型图及固有频率,进而计算出双螺杆膨胀机危险转速在44374~86622 r/min间变化。得出了影响主箱体变形的因素,并对如何使结构趋于稳定进行了简要说明;同时通过控制输入转速范围,有效避免了共振现象的产生,为双螺杆膨胀机主箱体的分析及共振研究提供了一定的数据支持,也进一步分析奠定了基础。
致谢
本文的研究工作得到浙江省基础公益研究计划项目(No.LGG18E050014)与浙江省大学生创新创业训练计划项目资助(项目编号:2019R435014)。
参考文献