![](https://lxxb.cstam.org.cn/fileLXXB/richhtml/0459-1879/richHtml_jats1_1/images/REemail.gif)
STUDY ON COUPLED VIBRATIONS OF ROTOR WITH AN ARBITRARY SLANT CRACK 1)
Jiao Weidong![](https://lxxb.cstam.org.cn/fileLXXB/richhtml/0459-1879/richHtml_jats1_1/images/REemail.gif)
通讯作者: 2)焦卫东, 教授, 主要研究方向: 转子动力学、智能检测与故障诊断. E-mail:jiaowd1970@zjnu.cn
收稿日期:2019-10-21接受日期:2020-02-15网络出版日期:2020-03-18
基金资助: |
Received:2019-10-21Accepted:2020-02-15Online:2020-03-18
作者简介 About authors
![](https://lxxb.cstam.org.cn/fileLXXB/richhtml/0459-1879/richHtml_jats1_1/images/down.png)
摘要
对包含不同类型裂纹(横裂纹、横-斜裂纹以及任意斜裂纹)的转子的耦合振动进行研究,以揭示裂纹转子在不同方向上刚度参数的变化规律及其交叉耦合机理,特别是由此引发的振动特征. 对于包含不同类型裂纹的转子轴段,采用六自由度Timoshenko梁单元模型对其进行单元建模,并基于应变能理论推导计算柔度参数和刚度矩阵. 在此基础上, 采用纽马克-$\beta$数值算法求解裂纹转子的运动方程,获得裂纹转子在单故障或多故障激励(不平衡激励、扭转激励或不平衡激励加扭转激励)作用下的耦合振动响应,进而分析耦合振动谱特征. 与横裂纹和横-斜裂纹相比,任意斜裂纹使转子刚度矩阵的交叉耦合效应更显著,导致转子发生更强烈的弯-扭耦合甚至是纵-弯-扭耦合振动.无论是在不平衡激励还是扭转激励作用下, 弯曲振动与扭转振动幅度都更大. 而且,包含不同类型裂纹的转子的耦合振动特征频率,例如旋转基频与二倍频、扭转激励频率及其边带成分的幅值,对裂纹面方向角具有不同的敏感性. 所得的这些研究结果,可以为转子裂纹的特征参数辨识与诊断提供理论依据.
关键词:
Abstract
The coupled vibrations of rotor with different type of cracks were studied, including transverse crack, transverse-slant crack and arbitrary slant crack. It aims to reveal the variation rule of stiffness parameters in different directions, their cross-coupling mechanism and especially the feature of resulted vibrations of the cracked rotor. The rotor segment with different type of cracks, including transverse crack, transverse-slant crack and arbitrary slant crack, was respectively modeled using the Timoshenko beam element, considering all six degrees of freedom including longitudinal, bending and torsional directions. Flexibility coefficients and then stiffness matrix was derived, based on strain energy theory. On this basis, motion equation of the cracked rotor was solved by numerical algorithm `Newmark-beta' to obtain the time-domain response of coupled vibrations of the cracked rotor under different excitations by single or multiple faults such as unbalanced excitation, torsional excitation and combined excitation of unbalanced with torsional. The spectral characteristics of coupled vibrations were then analyzed. Compared to either transverse crack or transverse-slant crack, the stiffness matrix governed by arbitrary slant crack is more populated with additional cross coupled coefficients, and its cross coupling effect is more significant, which contributes to stronger coupled vibrations in bending-torsional direction or even longitudinal-bending-torsional direction. The amplitudes of both bending vibration and torsional vibration are larger than that at the situation of either transverse crack or transverse-slant crack, under the action of either unbalanced excitation or torsional excitation. Moreover, the characteristic frequencies of coupled vibrations of rotor with different type of cracks including transverse crack, transverse-slant crack and arbitrary slant crack, are of different sensitivities to oriented angle of crack surface, for example the rotating fundamental frequency and its second harmonic frequency, the torsional excitation frequency and its sideband components. These research results can provide theoretical basis for both characteristic parameters identification and diagnosis of rotor cracks.
Keywords:
PDF (5953KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
焦卫东, 蒋永华, 李刚, 蔡建程. 任意斜裂纹转子的耦合振动研究 1). 力学学报[J], 2020, 52(2): 533-543 DOI:10.6052/0459-1879-19-292
Jiao Weidong, Jiang Yonghua, Li Gang, Cai Jiancheng.
引言
裂纹是旋转机械的常见故障, 对设备运行安全的潜在危害很大. 转子在旋转过程中, 裂纹面承受拉、压应力的交替作用,裂纹开/合(或呼吸)行为引起转轴刚度的周期性变化, 导致转子复杂的耦合振动[1-6]. 对裂纹转子的耦合振动机理与特性进行研究,是转子裂纹诊断的基础.对于最简单的横裂纹, 即裂纹面同时垂直于转轴轴线与基面, Darpe 等[7]深入研究了裂纹转子刚度变化的机理,分析了裂纹转子的纵向、弯曲与扭转耦合振动的特征. 随后,又将横裂纹推广至更一般的横-斜裂纹情形, 即裂纹面垂直于基面但不垂直于转轴轴线.从刚度变化机理与耦合振动特性两方面, 与横裂纹进行了对比研究. 相比于横裂纹,横-斜裂纹导致转子更多刚度参数发生耦合,引起弯曲、扭转甚至纵向耦合振动[8]. 纵观国内、外现有的研究,主要集中于横裂纹或横-斜裂纹转子的振动机理与特性问题[9-16].但在某些特殊工况下, 转子裂纹会呈现出更加复杂的几何形态,例如裂纹面既不垂直于基面也不垂直于转轴轴线, 即任意斜裂纹. 例如,在大扭矩和强弯矩载荷作用下裂纹会沿着螺旋方向扩展,从而形成螺旋裂纹或斜裂纹[17]; 由齿轮啮合力所导致的大扭矩,引发裂纹的斜向扩展[18]. 转子旋转过程中,张开型裂纹承受恒定方向的拉应力作用, 导致转子刚度发生定值削弱.不同于张开型裂纹, 呼吸型裂纹承受拉、压应力的交替作用,激起转轴刚度的周期性变化[1,19]. 因此当转子包含有呼吸型任意斜裂纹时,由于裂纹面两个方向角的交互作用,其刚度参数的变化规律及其交叉耦合机理明显不相同于横裂纹与横-斜裂纹转子,不同方向上的刚度参数发生广泛、强烈且复杂的耦合,引起具有不同特征的多自由度耦合振动[20]. 目前, 这方面的研究还不够深入.
本工作采用计算机仿真方法, 基于力学以及数值分析理论,研究包含不同类型裂纹特别是任意斜裂纹转子的变刚度机理,以及由此引发的不同故障激励作用下的耦合振动特性,以期为转子裂纹诊断提供参考数据.
1 裂纹转子的耦合振动机理
1.1 裂纹转子的刚度[7-8,20]
横裂纹、横-斜裂纹以及任意斜裂纹的坐标系统及其空间变换关系如图1 所示.图1中, xoy面为基面. 图1(d)描述了各坐标系统的空间变换关系. 显然,从最基础的横裂纹(坐标系xyz)出发, 其裂纹面绕 $z$ 轴旋转一定角度 $\theta _{1}$ $(0^{\circ}<\theta_{1}< 90^{\circ})$,即可得到横-斜裂纹(坐标系 $x' y' z' )$;横-斜裂纹的裂纹面绕其坐标轴 $y' $再旋转一定角度$\theta_{2}$ $(0^{\circ}<\theta_{2} <90^{\circ})$, 即可得到任意斜裂纹(坐标系 $x''y''z'')$.
图1
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/html/0459-1879-52-2-533/thumbnail/img_1.png)
图1裂纹转子的坐标系统
Fig. 1The coordinate systems of cracked rotor
在材料疲劳裂纹与断裂性能分析中,基于应变能理论的有限元方法得到广泛应用[21-24]. 例如,文龙飞等[21]重点研究了动载荷作用下扩展裂纹尖端应力强度因子的求解方法;王晓明等[22]将表征能量耗散的变量引入到应变能函数中, 形成新的弹性势的显式表达,从而得到精确匹配实验数据的数值模拟结果.本文采用Timoshenko梁单元对转子单元进行建模, 考虑如图1所示的纵向、弯曲及扭转所有6个方向自由度. 根据卡斯蒂利亚诺定理,裂纹单元的柔度参数表达为[25-26]
式中, $u_{i}$和 $P_{i}$ 分别为沿着第 $i$ 个坐标方向的节点位移与节点力. $U^{0}$ 为无裂纹单元的应变能, $U^{\rm c}$为裂纹导致的外加应变能.
无裂纹单元的弹性应变能表达式为
式中, $A = \pi R^{2}$ 为转子轴横截面积, $E$ 为杨氏模量, $G$ 为刚性模量, $I$为转轴截面面积矩, $I_{0}$ 为截面极惯矩, $\alpha_{\rm s}$ 为Timoshenko 梁剪切系数.
由裂纹导致的外加应变能为
式中, $E'=E/(1 - v)$, $m_{\rm s}= 1 + v$, $v$ 为泊松比.
裂纹面的位移可以用张开、滑移与剪开3种模式来描述. 基于三向应力分析,分别推导3种位移模式的应力密度因子 $K_{i}^{\rm I}$, $K_{i}^{\rm II}$ 与$K_{i}^{\rm III}$, 进而通过面积积分计算 $U^{\rm c}$, 最终得到裂纹单元的各个柔度参数$g_{ij}$, $i,j = 1,2,\cdots, 6$. 考虑裂纹单元各节点位移 $q_{i}$, $i =1,2,\cdots, 12$ 的静平衡条件, 有 $[q_{1-12}]^{\rm T} = T[q_{1-6}]^{\rm T}$, $T$为变换矩阵. 从而, 裂纹单元的刚度矩阵为 $K^{\rm c} = TG^{-1}T^{\rm T}$, 其中柔度矩阵 $G= [g_{ij}]$.
1.2 裂纹转子的振动响应
在全局静态坐标系 $q$下, 裂纹转子的运动学方程以矩阵向量形式表达为式中, $M^{\rm s}$, $C^{\rm s}$以及$K^{\rm s}$ 分别为总体质量、阻尼以及刚度矩阵, $f^{\rm s}$ 则为总体激励力向量.
总体质量矩阵 $M^{\rm s}$ 由单元质量矩阵组装而成, 单元质量矩阵的计算采用一致性方法[27], 考虑六自由度 Timoshenko梁的剪切变形与转动惯量效应. 同理, 总体刚度矩阵 $K^{\rm s}$ 也需要由所有单元(包括无裂纹单元与裂纹单元)的刚度矩阵借助合适的方法组装生成[28].对于总体阻尼矩阵 $C^{\rm s}$, 采用文献[7] 或 [8]所建议的比例阻尼进行估算, 估计式为 $C^{\rm s} =\alpha_{d}M^{\rm s} +\beta _{d}K^{\rm s}$. 此外, 在构建总体激励力向量 $f^{\rm s}$ 时, 需要全面考虑转子有限元模型中各个节点所受到的外部激励作用,例如转子质量不平衡离心力、动-静碰摩力、非线性油膜力或外加的扭转激励等[29-34].
相比于无裂纹转子, 裂纹转子的振动响应直接受裂纹单元刚度值的影响, 而由裂纹引起的刚度参数变化则需要借助总体应力密度因子值的符号由振动响应进行估算.实际仿真过程中, 假设总体质量与阻尼矩阵 $M^{\rm s}$ 与 $C^{\rm s}$保持不变, 只有 $K^{\rm s}$ 因裂纹的呼吸行为而不断变化,其值在转子每旋转一度后被更新一次. 具体实施时,以节点力估计总体的应力密度因子值, 其符号用于确定裂纹闭合线的位置,进而确定裂纹柔度系数积分运算的积分限, 获得裂纹单元柔度矩阵与刚度矩阵的估计.再经过静态坐标系变换, 即可组装为总体刚度矩阵$K^{\rm s}$. $K^{\rm s}$ 连同根据转子最新位置更新的 $f^{\rm s}$一起用于估算下一个旋转角度的转子振动响应. 如此不断重复, 即可获得裂纹转子的振动响应.
2 仿真分析
考虑一个两端支撑单盘转子系统, 裂纹位于转子盘右端靠近盘的位置. 转轴长度 $L =0.7$ m, 直径 $D = 0.015$ m, 转子盘质量 $m = 1$ kg. 整个转轴共划分为 14 个单元, 包含 15个节点. 裂纹几何参数包括裂纹深度 $a$ 与裂纹面方向角($\theta_{1}$ 与$\theta_{2})$, 裂纹的类型主要取决于后者.在随后的仿真计算中, 所用到的仿真参数如裂纹单元柔度系数的积分计算参数$\alpha_{i}$ 与$\delta_{i}$ 以及纽马克-$\beta $ 数值算法中的时间步长 $\Delta t$ 等,参照相关文献进行选取.
2.1 裂纹单元的变刚度特性分析
裂纹单元长度 $l = 0.7 / 14 = 0.05$ m, 裂纹深度比$\bar{a}/D = 0.3$, 裂纹面方位角设置为$\theta_{1} = 45^{\circ}$, $\theta_{2} = 60^{\circ}$. 如图2 所示,包含横裂纹(T)、横-斜裂纹(TS)以及任意斜裂纹(AS)单元的刚度参数分别以线型\mbox{"- - - - -''}、 "- - -''和"--''表示. 图2 中横坐标"CCL Position''意为裂纹闭合线(crack closure line, CCL)位置. CCL 概念由 Darpe等提出, 用于精确求解式(3)所示的面积积分[7].图2
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/html/0459-1879-52-2-533/thumbnail/img_2.png)
图2三类不同裂纹的交叉耦合刚度参数($k_{ij}$, $i, j = 1,2,\cdots,6)$
Fig. 2The cross-coupled stiffness coefficients $k_{ij}$, $i, j = 1,2,\cdots,6$ of three different types of crack
交叉耦合刚度参数 $k_{ij}$, $i\ne j$, 是造成裂纹转子不同方向振动耦合的内因[35]. 对于 AS 型裂纹转子,刚度参数交叉耦合的现象更明显. 如图2, 在水平剪切-垂直剪切($k_{23})$、水平剪切-扭转($k_{24})$、水平剪切-垂直弯曲($k_{25})$以及垂直剪切-水平弯曲($k_{36})$等方向均出现强烈的交叉耦合现象.
保持裂纹面方位角$\theta_{1} = 45^{\circ}$ 不变, $\theta_{2}$ 在 30$^{\circ}$到 90$^{\circ}$ 之间变化. 具有不同方位角的 AS 型裂纹单元的刚度特性曲线如图3 所示. 图中分别以线型\mbox{"- - - - -''、}"- - -''、"--''、"-+-''以及"-o-''线型按照$\theta_{2}$ 递增的顺序加以对比描述.
由图3 可见, 随着方向角$\theta_{2}$ 的增大, 在纵向($k_{11})$、垂直剪切方向($k_{33})$、扭转方向($k_{44})$以及垂直弯曲方向($k_{55})$的刚度值单调下降.但是, 在水平弯曲方向($k_{66})$以及水平剪切方向($k_{22})$则不存在这种变化趋势, 刚度参数曲线彼此出现了明显交叉, 表明 AS 型裂纹的裂纹面方向角$\theta_{1}$ 与 $\theta_{2}$ 之间的交互作用效应. 显然, 在整个旋转周期的不同位置, 这种交互效应对刚度的影响是不同的,特别是在水平弯曲($k_{66})$和水平剪切($k_{22})$方向, 呈现出明显的非线性特征.
图3
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/html/0459-1879-52-2-533/thumbnail/img_3.png)
图3方位角($\theta_{2})$对包含任意斜裂纹的转子单元的变刚度特性影响
Fig. 3The influence of oriented angle ($\theta_{2})$ on stiffness variation of the rotor element including an arbitrary slant crack
接下来, 考虑一个包含 AS 型裂纹的转子轴段, 保持裂纹面方位角$\theta_{1} =45^{\circ}$, $\theta_{2} = 60^{\circ}$ 不变, 裂纹深度 $a$ 以 1/10倍的转轴直径 $D$ (0.015 m)为增量从 0.001 5 m 均匀增大到 0.007 5 m, 即裂纹深度比$(\bar{a}/D)$从 0.1 均匀增大到 0.5, 裂纹转子的刚度特性如图4 所示. 对应于不同裂纹深度的刚度特性曲线, 采用与图3 相同的线型按照 $a$ 或 $\bar{a}$递增的顺序进行描述.
由图4 可以看到, 随着裂纹深度 $a$ 或 $\bar{a}$ 的增大, AS 型裂纹转子在所有 6个自由度方向上的正刚度参数 $k_{ii}$, $i = 1,2,\cdots, 6$ 均呈现出单调减小的趋势, 即裂纹深度越大, 裂纹转子的刚度值就越小. 而且,这种刚度值单调变化趋势是整体性的、连续性的, 发生在转子整个旋转周期的不同转角位置(以裂纹面闭合线位置"CCL position"来确定, 见图4 横坐标). 虽然裂纹转子的刚度参数值相对于转角位置的变化曲线明显是非线性的, 但是对应于不同裂纹深度的转子变刚度特性曲线几乎具有相同的走向或变化方向.
2.2 裂纹转子的耦合振动特征分析
考虑无外加的扭转激励情况. 不平衡质量偏心距为$1.6\times 10^{-5}$ m, 转子转速 22 rad/s (3.5 Hz), 约等于 1/10 倍的弯曲自然频率(35.2 Hz). 当无裂纹转子只受到不平衡激励作用时, 振动谱图中以旋转频率(或基频)为主,水平与垂直两个方向的弯曲振动基频的幅值水平差不多; 当转子包含不同空间方向($\theta_{1} = 45^{\circ}$ 且$\theta_{2} =30^{\circ}$, $60^{\circ}$ 和 90${^\circ})$的任意斜裂纹时, 裂纹导致转子的柔度增大, 基频分量的幅值水平相比于无裂纹情况也显著增加.在水平与垂直两个方向的弯曲振动谱中, 还包括二倍频和三倍频谐波分量. 而且, 从纵向与扭转方向的振动谱中可以发现明显的基频、二倍频以及微弱的四倍频谐波成分, 而此时裂纹转子系统只受到弯曲方向的不平衡激励, 并无外加的纵向与扭转方向激励, 说明此时系统的振动行为主要受弯曲与扭转耦合机理的支配. 对$\theta_{1} =45^{\circ}$, $\theta_{2}$ 分别为 $30^{\circ}$, $60^{\circ}$ 和 $90^{\circ}$ 三种情况进行对比研究, 发现水平、垂直弯曲方向振动谱的基频、二倍频、三倍频以及扭转方向振动谱的基频、二倍频、四倍频谐波分量幅值水平与裂纹面方向角$\theta_{2}$ 呈负相关关系, 即$\theta_{2}$ 越大, 幅值水平越低; 纵向振动谱则相反,其基频、二倍频、四倍频谐波分量幅值水平与$\theta_{2}$ 呈正相关关系, 即$\theta_{2}$越大, 幅值水平越高. 由于篇幅限制, 无外加扭转激励情况的仿真图未给出. 文献[8]尽管以 TS 型裂纹转子作为研究对象, 但在只有不平衡激励的情况下,转子的振动响应特性与这里研究的 AS 型裂纹转子具有较大的相似性, 可以作为参考,特别是该文献给出的图8 (无裂纹转子的不平衡振动响应)和图9 (裂纹转子的不平衡振动响应).
图4
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/html/0459-1879-52-2-533/thumbnail/img_4.png)
图4裂纹深度 $a$ 对裂纹转子变刚度特性的影响
Fig. 4The influence of crack depth $a$ on stiffness variation characteristics of cracked rotor
考虑存在外加的扭转激励情况. 此时, 不平衡激励与扭转激励同时作用. 采用谐波扭转激励 $T\sin(\omega_{\rm e}t)$, 其中$T = 10$ N$\cdot $m, 扭转频率$\omega _{\rm e} = 35$ Hz (约等于弯曲自然频率$\omega _{0})$. 在不同裂纹面方向角($\theta_{1} = 45^{\circ}$, $\theta_{2}=30^{\circ}$, $60^{\circ}$, 90${^\circ})$下弯曲 ($y$与$z)$、纵向($u)$与扭转($\theta$)方向振动响应的时域与频域波形如图5 $\sim\!$图7所示. 需要注意的是, 当$\theta_{2} = 90^{\circ}$ 时实际上对应的是横-斜裂纹.
从图5 $\sim\!$图7 可以观察到, 不平衡激励叠加扭转激励下裂纹转子的变刚度特性明显不同于无外加扭转激励情况. 当只有不平衡激励作用时, 裂纹转子的刚度参数与转轴旋转频率同步变化且趋势平缓; 不平衡激励与扭转激励共同作用时, 裂纹转子的刚度参数与扭转激励频率的变化趋势剧烈.
扭转激励作用下裂纹转子刚度的突变特性, 是导致拍振与振动调制的根本原因[7]. 其中, 振动调制表现为以扭转激励频率$\omega_{\rm e}$ 为中心、转轴旋转频率$\omega $ 为半带宽的对称边带调制现象.在水平和垂直两个方向上, 转轴旋转频率$\omega $ 以及扭转激励频率$\omega _{\rm e}$两侧的边带频率的幅值近似相等; 弯曲与扭转振动的耦合,使扭转激励频率出现在弯曲振动谱中, 而且被旋转频率及其高次谐波分量所调制.扭转激励频率$\omega _{\rm e}$ 与转轴旋转频率$\omega $ 及其高次谐波分量 $m\omega$ 之间的相互作用,导致围绕$\omega_{\rm e}$ 两侧的边带频率$\omega_{\rm e}\pm m\omega$ 的出现. 在图5 $\sim\!$图7的原始振动响应谱图中, 由于各特征频率的幅值彼此差异较大, 导致一些围绕扭转激励频率$\omega_{\rm e}$ 的高级调制边带成分无法被观察到.通过调整原始谱图的纵坐标对其进行局部放大, 可以清楚地显示出这些特征频率的存在. 例如图5 右上角的局部放大图所示的$\omega_{\rm e} \pm m\omega $ 成分. 在纵向上, 与只有不平衡激励情况相比, 不平衡激励叠加扭转激励下裂纹转子的振动谱出现明显变化. 扭转激励的引入, 导致谱图中高频分量的产生, 并在时域波形中出现显著的削波现象, 其机理见文献[7], 不再赘述.
图5
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/html/0459-1879-52-2-533/thumbnail/img_5.png)
图5不平衡激励与扭转激励共同作用下任意斜裂纹转子的不平衡振动响应($\theta_{1} = 45^{{\circ}}$, $\theta_{2} = 30^{{\circ}})$
Fig. 5Unbalance response of the arbitrary slant crack rotor with torsional excitation: $\theta_{1} = 45^{{\circ}}$, $\theta_{2} = 30^{{\circ}}$
图6
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/html/0459-1879-52-2-533/thumbnail/img_6.png)
图6不平衡激励与扭转激励共同作用下任意斜裂纹转子的不平衡振动响应($\theta_{1} =45^{{\circ}}$, $\theta_{2} = 60^{{\circ}})$
Fig. 6Unbalance response of the arbitrary slant crack rotor with torsional excitation: $\theta_{1} = 45^{{\circ}}$, $\theta_{2} = 60^{{\circ}}$
图7
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/html/0459-1879-52-2-533/thumbnail/img_7.png)
图7不平衡激励与扭转激励共同作用下任意斜裂纹转子的不平衡振动响应($\theta_{1} =45^{{\circ}}$, $\theta_{2} = 90^{{\circ}})$
Fig. 7Unbalance response of the arbitrary slant crack rotor with torsional excitation: $\theta_{1} = 45^{{\circ}}$, $\theta_{2} = 90^{{\circ}}$
最后, 将不平衡激励与谐波扭转激励联合作用下任意斜裂纹转子的弯曲振动谱特征参数($\omega$, $2\omega$, $\omega_{\rm e}$与$\omega _{\rm e}\pm \omega $)的幅值对方向角$\theta_{2}$ 的敏感性综合于图8. 横-斜裂纹转子的弯曲振动谱特征参数对方向角$\theta_{1}$ 的敏感性在图9中也对比给出.
从图9 可以清楚地看到, 横-斜裂纹转子的弯曲振动谱中旋转频率$\omega $ 以及扭转激励频率$\omega_{\rm e}$ 的第一级边带分量$\omega_{\rm e}\pm \omega $ 的幅值对方向角$\theta_{1}$ 相当敏感. 随着$\theta_{1}$ 从 $10^{\circ}$ 增大到 $90^{\circ}$, 这些分量的幅值迅速减小;而对于其他的谱分量如旋转频率$\omega $ 的二次谐波 $2\omega$、扭转激励频率$\omega_{\rm e}$及其第一级边带分量$\omega_{\rm e}\pm \omega $, 情况则有所不同. 例如, 在水平方向, 随着$\theta_{1}$ 增大, 旋转频率$\omega$ 的二次谐波 $2\omega$ 的幅值增大, 扭转激励频率$\omega_{\rm e}$的幅值减小; 而在垂直方向, 两者均增大. 当$\theta_{1} = 90^{\circ}$ 时,情况出现了变化, 即扭转激励频率$\omega_{\rm e}$的幅值在水平方向上增大而在垂直方向上减小. 实际上,此时的裂纹类型已经由横-斜裂纹转变为横裂纹了.
对于如图8 所示的任意斜裂纹转子振动, 这些特征频率的幅值与方向角$\theta_{1}$ 与$\theta_{2}$ 两者均有着密切的关系, 呈现出不同的变化模式. 当$\theta_{1}$ 很小时, 例如$\theta_{1} = 30^{\circ}$ 的情况, 随着$\theta_{2}$ 从 $10^{\circ}$ 增大到$90^{\circ}$, 弯曲振动谱中旋转频率$\omega $、扭转激励频率$\omega_{\rm e}$的第一级边带频率$\omega_{\rm e}\pm \omega $ 的幅值均单调减小; 当$\theta_{1}$ 中等大小时,例如$\theta_{1} = 45^{\circ}$ 的情况, 这些特征频率幅值先增大然后减小,趋势变化拐点近似在$\theta_{2} = 36^{\circ}$ 左右; 当$\theta_{1}$ 较大时, 例如$\theta_{1} = 60^{\circ}$ 或 $75^{\circ}$ 的情况, 这些特征频率幅值先减小后增大, 然后又减小, 变化拐点分别近似在$\theta_{2} =[25.71^{\circ}$, $45^{\circ}$] 和$\theta_{2} = [30^{\circ}$, $45^{\circ}$] 左右.
图8
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/html/0459-1879-52-2-533/thumbnail/img_8.png)
图8任意斜裂纹转子振动谱特征参数对方向角的敏感性: $\theta_{1} =[30^{{\circ}}$, 45$^{{\circ}}$, 60$^{{\circ}}$, 75$^{{\circ}}$],$\theta_{2} = 10^{{\circ}}\!\sim\! 90^{{\circ}}$
Fig. 8Sensitivity of spectral parameters of the arbitrary slant crack to orientation angle: $\theta_{1} = [30^{{\circ}}$, 45$^{{\circ}}$, 60$^{{\circ}}$, 75$^{{\circ}}$], $\theta_{2} = 10^{{\circ}}\!\sim\! 90^{{\circ}}$
图9
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/html/0459-1879-52-2-533/thumbnail/img_9.png)
图9横-斜裂纹转子振动谱特征参数对方向角的敏感性: $\theta_{1} =10^{{\circ}}\sim 90^{{\circ}}$, $\theta_{2} = 90^{{\circ}}$
Fig. 9Sensitivity of spectral parameters of the transverse slant crack to orientation angle: $\theta_{1} = 10^{{\circ}}\sim 90^{{\circ}}$, $\theta_{2} =90^{{\circ}}$
横-斜裂纹与任意斜裂纹转子弯曲振动谱中特征频率幅值的不同变化特性,体现了两类裂纹对裂纹面方向角参数敏感性的差异. 应该注意的是, 在$\theta_{1}$较大的情况, 任意斜裂纹转子弯曲振动谱中$\omega $ 和$\omega_{\rm e}\pm \omega $两者的幅值在第一个减小阶段的变化趋势相比于图9 所示的横-斜裂纹情况,突变性更加强烈, 这说明任意斜裂纹比横-斜裂纹对裂纹面方向角的敏感度更高. 此外,与横-斜裂纹相比, 如图8 所示的任意斜裂纹具有明显的方向敏感性,暗示了任意斜裂纹转子所具有的不同振动特征. 文献[8]已导出一些用于裂纹类型(横或横-斜)辨识的有效准则. 本仿真分析所获得的一些结果,则有助于揭示任意斜裂纹转子的振动特征,便于对振动检测所获取的谱特征甚至裂纹参数进行评估,例如弯曲振动谱中旋转基频与二倍频分量、扭转激励频率及其边带成分以及裂纹面方向角等.
3 结论
由于裂纹面两个方向角的交互作用,任意斜裂纹转子在不同方向上刚度参数的变化规律及其交叉耦合机理明显不同于横裂纹与横-斜裂纹.这种交互效应不仅体现在刚度参数的量值上, 还体现在刚度曲线的变化趋势上,具有明显的非线性特征. 任意斜裂纹所具有的复杂变刚度特性,导致裂纹转子发生复杂的非线性振动,并在纵向、弯曲与扭转多个自由度方向上发生强烈耦合.无论是在不平衡激励还是在扭转激励的作用下, 弯曲振动与扭转振动的幅度都更大,危害也更严重.任意斜裂纹转子相比于其他两类裂纹转子, 其振动谱特征也不同.从弯曲振动谱特征参数(如旋转基频与二倍频、扭转激励频率及其边带成分的幅值)的敏感性分析结果来看,这些参数对裂纹面方向参数相当敏感. 从而,通过比较弯曲振动谱中垂直与水平方向的扭转激励频率的幅值大小,可以对裂纹面的方向参数进行估计. 此外, 在只有不平衡激励情况下,扭转振动的时域与频域波形显示: 当裂纹类型从任意斜裂纹转变为横-斜裂纹时,特征频率的幅值呈现单调变化, 这一结果对于裂纹类型的辨识也是有帮助的.
参考文献 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子
[本文引用: 2]
[本文引用: 1]
[本文引用: 1]
[本文引用: 6]
[本文引用: 5]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 2]
[本文引用: 2]
[本文引用: 2]
[本文引用: 2]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]