货物在船舱存储时,需要进行水平搬运和竖直搬运,升降装置是该环节的重要设备之一。某船用货物升降装置具有提升载荷大、运行速度高的特点,使用过程中存在振动和噪声辐射等问题,对船舱内工作环境和产品本身的性能产生影响。当结构发生共振时,设备的动应力和噪声会显著增加,影响设备使用寿命。为保证设备长期可靠工作,有必要对该升降装置工作时的振动特性进行测试和分析。
振动烈度是能够反映设备振动状态的实用有效的特征量,目前在大、中型机电设备的振动检测和状态评估研究中已长期应用[1 – 4],其理论基础比较成熟,利用该理论开展升降装置状态评估是完全可行的。本文在工作状态下对升降装置的振动加速度进行测试和分析,并使用频域算法对振动烈度进行分析,旨在掌握该升降装置的振动特性和振动烈度等级。在Ansys Workbench平台研究该升降装置的振动模态,对比了振动测试和模态分析结果,用于指导设备振动特性改进设计。
1 计算和分析方法 1.1 振动烈度计算方法振动烈度是表示机电设备振动强烈程度的量,通常采用振动速度有效值的合成值表征振动烈度。试验中,测得升降装置基座各测点在各方向的振动加速度离散信号,通过时域积分的方法得到各测点在各方向的振动速度信号,然后通过频域计算方法计算各测点在各方向的振动烈度,最后合成得到各测点的振动烈度。振动烈度的频域计算方法[5 – 6]利用离散傅里叶变换(DFT)理论,可得频域上信号
$X(k) = \sum\limits_{n = 0}^{N - 1} {x(n){e^{ - j2nk{\text{π}} /N}}} ,k = 0,1, \cdots N\text{。}$ | (1) |
式中:
${A_k} = \frac{2}{N}\left| {X(k)} \right|,k = 0,1, \cdots N/2\text{,}$ | (2) |
${f_k} = \frac{{k{f_s}}}{N},k = 0,1, \cdots N/2\text{。}$ | (3) |
计算频率范围为
$ \begin{split} {V_{rms}} & = \sqrt {\frac{1}{2}\sum\limits_{k = {k_a}}^{{k_b}} {A_k^2} } = \sqrt {\frac{1}{2}\sum\limits_{k = {k_a}}^{{k_b}} {{{\left( {\frac{2}{N}\left| {X(k)} \right|} \right)}^2}} }= \\ & \frac{{\sqrt 2 }}{N}\sqrt {\sum\limits_{k = {k_a}}^{{k_b}} {{{\left| {X(k)} \right|}^2}} } \text{。} \end{split} $ | (4) |
式中:
振动烈度的频域计算方法相对于时域计算方法适应性强,能够提取振动位移、速度和加速度的烈度特征,可以根据需要选取一定的频域范围计算烈度值。
根据国家标准GB/T16301《船舶机舱辅机振动烈度的测量和评价》,当测点的振动按照3个相互垂直方向测试时,对各测点的振动烈度按照下式合成:
${V_S} = \sqrt {{{\left( {\frac{{\displaystyle\sum {{V_X}} }}{{{N_X}}}} \right)}^2} + {{\left( {\frac{{\displaystyle\sum {{V_Y}} }}{{{N_Y}}}} \right)}^2} + {{\left( {\frac{{\displaystyle\sum {{V_Z}} }}{{{N_Z}}}} \right)}^2}}\text{。} $ | (5) |
式中:
设备的模态由结构本身特性和材料的特性所决定,与外部载荷大小无关。模态分析利用系统固有频率的正交性,以系统的各阶模态向量所组成的模态矩阵作为变换矩阵,对选取的物理坐标进行线性变换,使得振动系统以物理坐标和物理参数所描述的、相互耦合的运动方程组能够变为一组独立的模态方程。模态分析主要有试验模态分析和有限元模态分析2种方法。试验模态分析是将试验与理论分析结合起来辨识系统的固有频率和振型;有限元模态分析主要运用有限元法对振动结构进行离散,求取振动结构固有频率及其相应振型。由于升降装置体积较大、结构较复杂,本文采用有限元模态分析方法进行研究。升降装置主体为钢结构体,可视作小阻尼多自由度系统[7 – 8],其振动微分方程为:
${ M}\ddot X(t) + { C}\dot X(t) + { K}X(t) = P(t)\text{。}$ | (6) |
式中:
式(6)为有阻尼的多自由度系统的强迫振动微分方程,而固有频率分析一般可看作无阻尼自由振动,则式(6)变为:
${ M}\ddot X(t) + { K}X(t) = 0\text{,}$ | (7) |
有非零解的充分必要条件是系数矩阵行列式等于0,即特征方程为:
$\left| {K - {\omega ^2}M} \right| = 0\text{。}$ | (8) |
求解特征方程可得到升降装置的固有频率。目前有较多的有限元模态分析的软件,本文选用Ansys Workbench软件进行升降装置的模态分析研究。
2 试验情况根据升降装置的电机工作转速为1 700 r/min,确定试验系统组成如图1所示。其中使用1台西门子SCM05型动态信号分析仪,12个美国PCB 333B32型加速度计,测量范围为0.5~3 000 Hz,灵敏度为100 mV/g,12个333B32型加速度计合成4组三向加速度计,能够满足测量要求。
在电机额定工作转速1 700 r/min,提升载荷2 t的条件下进行测试,升降装置的顶部和底部基座分别固定在船体结构上,在升降装置的顶部和底部基座上各布置2个测点,对每个测点采集X,Y,Z三向的振动加速度,采样频率为20 480 Hz,其中1个加速度计测量方向垂直基座安装面,升降装置的顶部基座具体测点位置和方向如图2所示。
试验中,升降装置稳定运行过程中利用加速度计对4个测点的振动加速度信号进行采集,每个测点有X,Y,Z三向加速度信息,选取稳定运行阶段的5 s的信息进行分析,下面列出2号测点和4号测点分析结果,2号测点和4号测点加速度时域图如图3所示。通过信号处理和傅里叶变换得到加速度1~1 000 Hz频谱图如图4所示。通过加速度信息时域积分的方法得到测点各方向的振动速度信号[9 – 11],速度时域图如图5所示。试验中,振动测量误差主要来源于加速度计线性误差、频率响应误差、动态信号分析仪误差和粘贴时传感器轴向与测点轴线偏差所引起误差。振动加速度计线性误差为1%;振动加速度计频率响应误差为5%;动态信号分析仪误差为0.5%;粘贴引入误差为1%;按均匀分布计算,则振动烈度试验系统的扩展不确定度为:U=1.8×10–2(k=2),满足GB/t 16301的测量要求。
统计2号、4号测点在0.5~1 000 Hz范围内振动加速度比较大时的频率,如表1所示。电机额定工作转速1 700 r/min,对应转动频率为28.33 Hz,从表1和图4可以看出,振动加速度比较大时均出现在电机转动频率的倍频处,且大部分出现在电机转动频率的10倍频(283.3 Hz)以上,说明电机转动频率低倍频工作时升降装置振动力较小。
按照式(1)~式(5),根据各测点的振动速度计算各测点的振动烈度,结果如表2所示。
可以看出1号、2号测点在Y向振动烈度最大,3号、4号测点在Z向振动烈度最大,根据GB/T 16301对升降装置整体进行评价,升降装置整体在Y向振动烈度最大,振动烈度合成值为6.74E–4 m/s,符合第一类辅机A级(优良工作状态)的要求,说明升降装置工作时振动强度小,对周围环境的设备影响较小。
3.2 模态分析建立正确的模型是有限元模态分析的基础。本文使用Pro/E三维软件对升降装置进行建模。Pro/E是目前先进的计算机辅助设计、分析与制造软件之一,具有强大的实体造型、虚拟装配和产生工程图等设计功能,同时Pro/E可以与Ansys Workbench联合仿真,便于模型处理和提高仿真效率。为提高仿真单元质量、减少仿真分析耗时,考虑到升降装置传动丝杠为受拉状态,不影响整体结构的模态振型,对其和升降装置模型中的孔、倒角等结构适当的简化。在Pro/E中完成三维模型建立后导入Ansys Workbench中,选用modal模块进行分析。首先进行网格划分,在门架主体采用大尺寸网格划分,其余的零部件采用确定尺寸进行网格划分,如图6所示。
按照实际使用工况对升降装置的顶部基座和底部基座采用固定约束。以升降装置中心位置为基准建立了有限元网格模型,该模型中总节点数为2 288 131个,总单元格数为1 154 788个。通过对升降装置模型进行有限元模态分析,可得到前8阶模态振型图和固有频率,如图7和表3所示。
图7中,1阶模态图表示门架中部沿X轴方向弯曲变形,工作中为减小门架的变形,需采用刚性较强的结构,2阶模态图表示门架中部整体沿Z轴方向的单方向弯曲变形,3阶模态图表示门架顶部和中部沿X轴方向的弯曲变形,4阶模态图表示门架中部在ZY平面的两方向弯曲变形,5阶模态和6阶模态图表示升降台前沿Y方向的弯曲变形,7阶模态图表示升降台绕X轴的扭曲变形和门架沿X轴方向的弯曲变形,8阶模态图表示门架中下部在ZY平面的两方向弯曲变形。通过对整个升降装置进行模态分析,前8阶固有频率集中在50~300 Hz的范围内。当升降装置空载或满载工作时,避免激振频率与固有频率相同或相近,防止升降装置零部件发生较大振动而失效。
对比表2和表3,可以看出在升降装置电机转动频率的10倍频(283.3 Hz)内,4号测点振动加速度在转动频率的5倍频(141.7 Hz)、6倍频(170 Hz)和3倍频(85 Hz)附近出现较大值,分别对应升降装置的3阶、5阶和2阶模态固有频率,说明在固有频率与转动频率各倍频接近时,易产生较大的振动力,在设计时应尽量将两者分开。
4 结 语根据振动烈度测试结果对升降装置振动强度进行评价,升降装置在Y向振动烈度最大,整体振动烈度合成值为6.74E–4 m/s,符合第一类辅机A级(优良工作状态)的要求,说明升降装置工作时振动强度小,对周围环境的设备影响较小。通过对比振动加速度测试结果和升降装置固有频率,说明在固有频率与转动频率各倍频接近时,更易产生较大的振动力,在设计时应尽量将两者分开,对于设计过程中降低设备振动提供指导和参考。
[1] |
韩伟, 贾启芬, 邱家俊. 异步电机定子的振动与模态分析[J]. 振动与冲击, 2012, 31(17): 91-94. HAN Wei, JIA Qi-fen, QIU Jia-jun. Vibration and modal analysis for stator of an induction motor[J]. JOURNAL OF VIBRATION AND SHOCK, 2012, 31(17): 91-94. DOI:10.3969/j.issn.1000-3835.2012.17.015 |
[2] |
王富民, 李捷, 杨建伟, 等. 地铁齿轮箱箱体模态及谐响应分析[J]. 机械传动, 2015, 39(9): 146-150. WANG Fumin, LI Jie, YANG Jianwei, et al. Modal and harmonic response analysis of subway gearbox housing[J]. Journal of Mechanical Transmission, 2015, 39(9): 146-150. |
[3] |
张顺风, 吴怀超, 李哲. 一种可调式高速精密轧辊磨头的模态及谐响应分析[J]. 机床与液压, 2018, 46(1): 115-120. ZHANG Shunfeng, WU Huaichao, LI Zhe, et al. Modal and harmonic response analysis of one adjustable high-precision roll grinding head[J]. Machine Tool and Hydraulics, 2018, 46(1): 115-120. DOI:10.3969/j.issn.1001-3881.2018.01.024 |
[4] |
马维金, 熊诗波. 自动化立体仓库巷道堆垛机的振动测试与工作模态分析[J]. 中国机械工程, 2004, 15(4): 287-289. MA Wei-jin, XIONG Shi-bo. Vibration testing and operational modal analysis of automated warehouse lane stacker[J]. China Mechanical Engineering, 2004, 15(4): 287-289. DOI:10.3321/j.issn:1004-132X.2004.04.002 |
[5] |
樊新海, 赵智勇, 安钢, 等. 机械振动烈度的频域算法研究[J]. 装甲兵工程学院学报, 2008, 22(1): 42-45. FAN Xin-hai, ZHAO Zhi-yong, AN Gang, et al. Research on the arithmetic of mechanical vibration severity in frequency domain[J]. Journal of Academy of Armored Force Engineering, 2008, 22(1): 42-45. DOI:10.3969/j.issn.1672-1497.2008.01.010 |
[6] |
班希翼, 余建勇, 贺小龙. 有源设备对车辆舒适度指标和设备振动烈度影响分析[J]. 噪声与振动控制, 2017, 37(4): 95-99. BAN Xi-yi, YU Jian-yong, HE Xiao-long. Influence of equipment excitation on vehicle comfort index and equipment vibration severity[J]. Noise and Vibration Control, 2017, 37(4): 95-99. DOI:10.3969/j.issn.1006-1355.2017.04.019 |
[7] |
李连生, 唐斌, 王乐. 往复压缩机曲轴变工况条件下有限元分析[J]. 流体机械, 2011, 39(12): 14-17. LI Lian-sheng, TANG Bin, WANG Le. FEM analysis of reciprocating compressor crankshaft at variable working condition[J]. Fluid Machinery, 2011, 39(12): 14-17. DOI:10.3969/j.issn.1005-0329.2011.12.004 |
[8] |
吴文轩, 张纪平, 马维金, 等. 汽车变速箱的振动特性与模态联合分析研究[J]. 机械传动, 2017, 41(8): 137-140. WU Wen-xuan, ZHANG Ji-ping, MA Wei-jin, et al, Study on the vibration characteristic and modal analysis of automobile transmission shell[J]. Journal of Mechanical Transmission, 2017, 41(8): 137-140. |
[9] |
刘思远, 李晓明, 刘建勋. 基于振动烈度的液压泵故障多信息特征提取方法研究[J]. 振动与冲击, 2018, 37(14): 269-275. LIU Si-yuan1, LI Xiao-ming, LIU Jian-xun. Multi-information fault feature extraction method for hydraulic pumps based on the vibration intensity[J]. Journal of Vibration and Shock, 2018, 37(14): 269-275. |
[10] |
于金朋, 黄冠华, 张卫华, 等. 基于振动能量的内燃动车组柴油机组振动烈度分析[J]. 北京交通大学学报, 2014, 38(6): 47-51. YU Jin-peng, HUANG Guan-hua, ZHANG Wei-hua, et al. Analysis of vibration intensity for diesel group based on vibration energy for DMU[J]. Journal of Beijing jiaotong university, 2014, 38(6): 47-51. DOI:10.11860/j.issn.1673-0291.2014.06.008 |
[11] |
顾名坤, 吕振华. 基于振动加速度测量的振动速度和位移信号识别方法探讨[J]. 机械科学与技术, 2011, 30(4): 522-526. GU Ming-kun, LV Zhen-hua. Identification of a mechanism's vibration velocity and displacement based on the acceleration measurement[J]. Mechanical Science and Technology for Aerospace Engineering, 2011, 30(4): 522-526. |