船舶中梁的裂纹振动诊断方法.docVIP

船舶中梁的裂纹振动诊断方法 居桦桦 赵德有 李志强 （大连理工大学船舶工程学院，辽宁 大连 116023） 摘 要：板、梁是船舶的基本构件，如果梁和板产生裂纹，将会导致船舶失事。当梁出现裂纹时，其刚度会随之发生改变。为此，基于振动诊断方法，本文首次将修正的线弹簧模型，引入到铁木辛柯梁传递矩阵法中，计算出梁损伤前后的振动信息。利用振动信息，本文提出振型差值曲率方法，来诊断梁的裂纹位置和裂纹深度，经研究证明，此方法可供梁的裂纹诊断参考。 关键词：结构振动；传递矩阵；裂纹诊断 1 引言　 梁是船舶的基本组成部分之一, 在外激励作用下, 易产生裂纹。裂纹的存在给结构造成了潜在的危害,甚至可能导致灾难性事件发生。因此, 在发生海损事故前, 及时发现梁存在的微裂纹是十分重要的。近三十年来,利用结构的动态特性进行损伤诊断一直是国际学术界和工程界关注的热点之一。本文基于铁木辛柯梁理论，采用传递矩阵法计算梁损伤前后的振动信息，并且引入修正后的线弹簧模型来模拟裂纹。相对于传统的诊断指标，本文提出的新诊断指标――振型差值曲率对微小裂纹更为敏感。在对工程结构裂纹损伤故障分析中，一方面要用结构动力学的理论对结构进行分析，另一方面对结构系统进行试验研究。为了探讨梁的裂纹位置和裂纹深度，需要大量的试验模型，试验造价昂贵，时间长，本文采用数值仿真技术，即通过结构振动计算，获得梁的振动信息，进行裂纹诊断。 梁的裂纹损伤振动信息获得，常用的数值计算方法有有限元法和传递矩阵法，对于梁来说，由于有限元法处理梁裂纹处需细化单元，工作量是较大的，传递矩阵法就更为方便。为了适应长梁和短梁及变端面梁，本文采用铁木辛柯梁传递矩阵法，不仅考虑弯曲变形，同时还考虑了剪切变形影响，因此计算得到的振动信息更为准确。 2 裂纹梁的传递矩阵 2.1 传递矩阵法 把梁分为若干段，利用铁木辛柯梁的振动微分方程式，可以求出每一段梁挠度、转角、弯矩和剪力的状态矢量，然后利用各段梁连接处的连接条件和力的平衡条件求出迁移矩阵。 铁木辛柯振动微分方程式为： （1）式中：E—材料弹性模量，G—材料剪切弹性模数，—第j段梁的截面惯性矩，—第段梁的有效剪切面积，—第段梁的剖面回转半径，—第段梁单位长度上的质量，k—剪切滞后影响系数。 对于自由振动时，运动微分方程的解为： （2） （3） 式中：，， 其中：—第节点振动的圆频率，—第段梁的长度。 方程式（3）的通解为： (4) 根据（4）式，可以导出梁段的场传递矩阵 （5） （6） 式中： , , 其中：和由式（3）导出的特征方程求得 （7） 通过上面的分析,推导出梁最左端和最右端的状态参量之间的传递关系为: （8） 式中： ， 由式（8）式出发，利用两端的边界条件，消去首尾某些已知边界条件所含的物理量，然后利用非零条件即可得到频率方程式，求得频率后亦可求出各点处的状态矢量。由于频率方程式是超越方程，故有无穷多个固有频率和固有振型。 2.2 修正的线弹簧模型 文献[2]中首次提出线弹簧模型用于模拟裂纹，文献[3]对文献[2]的线弹簧模型进行修正，并且利用修正的模型, 求解欧拉裂纹梁振动的解析解。传统的弹性铰模型，只涉及到单一荷载—弯矩，修正后的线弹簧模型包括弯距和剪力两种载荷。图1是文献[3]中的修正线弹簧模型基本构成及其相关作用力和方向。 图1 修正的线弹簧模型 本文首次将文献[3]中的线弹簧模型引入到铁木辛柯梁传递矩阵法中。将梁分若干段时，裂纹处必须看作一个节点。裂纹处连续条件为： 其中: 这里，，分别是与弯矩和剪力相对应的裂纹梁柔度系数。其中： (9) (10) 式中：为弹性模量, 为裂纹梁宽度, 为裂纹的长度。是由弯矩产生的应力强度因子, 为剪力引起的应力强度因子, 它们分别与裂纹结构的几何形状系数有关,写成： 其中：,是裂纹梁的高度。 3 裂纹诊断指标的选取 结构产生裂纹会引起结构某些固有参数发生改变, 在裂纹诊断中, 就利用这些变化的参数进行诊断。