您的当前位置:首页正文

基于振动模态分量的结构损伤定位方法

2021-03-30 来源:个人技术集锦
第39卷第5期 2013年1O月 兰州理工大学学报 Vo1.39 No.5 OeL 2013 Journal of1.anzhouUniversity ofTechnology 文章编号:1673-5196(2013)05-0116-05 基于振动模态分量的结构损伤定位方法 方、治华 ,杜晓旭 (1.内蒙古科技大学建筑与土木工程学院,内蒙古包头014010;2.平顶山煤矿机械有限责任公司,河南平顶山467001) 摘要:三维实体结构振动时,结构上的任一点在三维空间各个方向均产生位移分量,振动模态也包含各个方向的模 态分量.当结构损伤时这些模态分量会发生变化,但对损伤的敏感程度是不同的.基于这一特点,提出位移模态分 量变化率、应变模态分量变化率和第一阶模态下的轴向位移差变化率3种针对三维实体结、构的损伤标识量,并通 过算例对损伤位置和损伤程度与3种标识量的关系进行计算分析,研究3种标识量对损伤位置的敏感程度.计算 实例表明,3种结构损伤标识量都在损伤部位出现尖峰,可以用于损伤定位。3种损伤标识量中,位移模态分量变化 率较小,对损伤的敏感程度不高;应变模态分量变化率较大,对损伤的敏感性较高;第一阶模态下的轴向位移差变 化率最大,对损伤最敏感. 关键词:振动模态分量;结构损伤;定位方法;敏感性 中图分类号:TU312.3 文献标识码:A Structural damage positioning method based on vibration modal component FANG Zhi-hua .DU Xiao-xu (1.School of Civil Engineering,Inner Mongolia University of Science&Technology,Baotou 014010,China;2.Pigdinngshan Coal Mining Machinery Co.,Ltd,Pigdiongshan,467001,China) Abstract:When a three dimensional structure vibrates.every point in the structure will be displaced three-dimensionally.Vibration mode also includes modal components in itself in all directions.These mo— dal components will change when the structure is damaged.But their sensitivity to the damage are differ— ent.According to this characteristics,three identification indexes of structure damage,such as the change rate of variation displacement modal component,the change rate of variation strain modal component,and the first-order change rate of axial displacement difference were presented.The calculation and analysis of the relationship between damage loeation and damage extent were carried out with these identification inde— xes by means of a calculation example.The sensitivity degree of three identification indexes to the damage location was studied.The calculation results showed that the peak appeared in these identification indexes spikes at the damage location.which could be used to determine the damage location.Among these three identification indexes,the change rate of variation displacement modal component was smaller at the dam— age location,SO sensitivity degree of damage was lower.The change rate of variation strain modal compo— nent was 1arger at the damage 1ocation than the change rate of variation displacement modal component, but the sensitivity degree of damage was higher.The first-order change rate of axial displacement differ— ence was the largest at the damage location.being the most sensitive damage. Key words:vibration modal component;structural damage;positioning method;sensitivity 当结构存在损伤时,会引起结构本身动力参数 以及模态参数的改变,在物理状态空间中表现为刚 度下降,柔度增大;在模态状态空间中表现为固有频 收稿日期:2012—06—22 率下降,阻尼比增大,振动模态发生变化.因此可以 通过结构动态参数变化对损伤做出早期诊断.基于 动态特性的结构损伤识别方法的研究始于2O世纪 7O年代.近4O年来,人们为了解决结构中的早期损 基金项目:内蒙古自然科学基金(20080404MS0810) 作者简介:方治华(1962-),男,内蒙古包头人,教授. 伤识别问题,进行了很多的研究工作,出现了许多种 方法[1.引.这些方法选取能够表征结构损伤出现的标 第5期 方治华等:基于振动模态分量的结构损伤定位方法 识量作为观测值,通过观测这些标识量的变化来判 断结构损伤的位置和程度,损伤标识量的选取有固 有频率、位移模态、应变模态、曲率模态及由此构造 出的频响函数、柔度矩阵、刚度矩阵等.目前大量的 动态损伤识别研究的领域是关于简单的平面杆、梁、 板结构,对于空间结构的研究有针对于海洋石油平 台[引,桥梁[ ]和航空航天结构[ ]的,其组成成分仍然 是杆、梁、板.三维实体结构损伤识别的研究比较少, 研究最多的是转子,如曲轴[6],也有学者研究了裂纹 齿轮E ]和裂纹叶片[8],但其标识量的选取仅限于固 有频率.固有频率作为结构整体特性参数,对结构局 部损伤往往不敏感,因此适用范围有限;用其它标识 量进行损伤识别虽然可以,但敏感程度不高. 由于三维结构振动在空间各方向上产生位移和 应变,而损伤可能只出现在某一方向,对该方向的振 动模态分量影响较大,其他方向影响较小,因此用对 损伤敏感的模态参数分量构造损伤标识量可以提高 敏感程度.本文根据三维实体结构振动的特点,提出 以模态在空间某一方向的分量构造损伤标识量的思 路,建立了位移模态分量变化率、应变模态分量变化 率和第一阶模态下的轴向位移差变化率3种损伤标 识量.计算实例表明,3种损伤标识量都能够较准确 地识别出三维实体结构损伤位置,但对损伤的敏感 程度相差较大. 1基于空间振动模态分量的结构损伤 定位标识量 三维实体结构振动时,结构上的任一点在三维 空间各个方向均产生位移分量.模态分析时,无论是 位移模态还是应变模态,也会有3个方向的模态分 量.当结构损伤时这些模态分量会发生变化,但对损 伤的敏感程度是不同的.在目前研究较多的平面结 构中,由于振动形式较简单,这一影响并不突出,但 是在三维空间结构中,这一影响需要考虑;选取较敏 感的参数(分量)构造损伤标识量是提高三维空间结 构损伤标识量敏感性的有效途径.本文提出以模态 在空间某一方向的分量构造出位移模态分量变化 率、应变模态分量变化率和第一阶模态下的轴向位 移差变化率3个损伤标识量. 1.1位移模态分量变化率 设三维空间实体结构损伤前后第i阶 方向上 的位移模态分量为 和 ,定义第 阶模态在第 个节点处 方向的位移模态分量变化率为 RD}z’( )一 尊 ×lOO %,O ( 一1,2,…,s) 考虑实际测量时只能测得实体结构表面点的振 动形态,因此在分析计算时只讨论结构表面上点的 位移模态分量变化率.具体操作时,可以在实体结构 表面选取一条直线路径(例如轴向母线),提取该路 径上等间距节点的位移模态分量,计算各点的位移 模态分量变化率. 1.2应变模态分量变化率 设三维空间实体结构损伤前后第i阶 方向上 的应变模态分量为e ’和e ,,定义第i阶模态在第 个节点处z方向的应变模态分量变化率为 Y (. )一 ×lOO% £ ) ( 一1,2,…, ) 同样考虑实际测量时只能测得表面点的振动形 态,因此在计算中,根据振动特点通过定义表面路径 (例如轴向)的方法提取实体表面一条线段上的应变 模态,计算应变模态分量变化率. 1.3第一阶模态下的轴向位移差变化率 当三维实体结构按第一阶模态做主振动时,如 果结构存在裂纹损伤,必然引起局部刚度的降低,裂 纹处在振动时处于张开、闭合的运动中,使裂纹两边 节点之间相对位移增大.基于这一分析,提出以第一 阶模态下的轴向位移差作为标识量.具体操作时,在 实体结构表面选定一条平行轴线的母线,从母线上 等间距取若干点作为取值的节点,提取其在第一阶 模态下的轴向模态位移分量,计算相邻两节点之间 的轴向模态位移差.定义损伤前后的第一阶模态下 的轴向位移差变化率A 为 觚= 等 ×lOO% (忌一1,2,…,s) 式中:X 一嬲和础 一 分别表示结构损伤前 后第k个节点和第k+1个节点之间在第一阶模态 下的轴向位移差. 要根据裂纹的宽度和计算工作量选取节点数 量,既不能节点间距过于大,使轴向模态位移差对裂 纹敏感性降低,也不能节点间距过小,使计算工作量 太大. 2算例 以包钢薄板厂的圆柱形轧辊为例计算含表面裂 纹损伤时,轧辊表面一条母线上的位移模态分量变 化率、应变模态分量变化率和第一阶模态下的轴向 位移差变化率3个损伤标识量.轧辊可分成工作段 兰州理工大学学报 第39卷 和支撑段两个部分.中间工作段的长度为2 ITI,直径 为1 m;两端支承段均为长1 13"1,直径0.5 m;轴向设 为 方向,垂直轴线分别建立正交的Y和 轴.利用 ANSYS软件建立无裂纹轧辊的三维有限元模型; 在轧辊上设定一条圆环形裂纹,建立含裂纹轧辊的 三维有限元模型,见图1.通过模态分析,分别获得 无裂纹和有裂纹轧辊的位移模态和应变模态参数. 图1含裂纹轧辊的有限元模型 Fig-1 Finite element model of cracked roller 2.1裂纹损伤位置和损伤程度对轧辊位移模态分 量变化率的影响规律 2.1.1裂纹损伤位置对轧辊位移模态分量变化率 的影响规律 在轧辊上设置深度为z=0.08 rn,圆心角为6O。 的裂纹.分别选取裂纹位置z一0.1 ITI等7个裂纹 位置,在轧辊工作段表面取一条母线,计算第一阶 方向位移模态分量变化率,其变化曲线如图2所示. 5 褥 4 删 3 桶 登z 鑫 I l 0 l 2 轴向坐b ̄,x/m 图2第一阶z方向位移模态分量对于裂纹损伤位置的 变化率 F.晷2 Variational rate of f'wst-order z-directional dis- plao既netIt modal component’vitll crack locations 从图2可以看出,在相同的裂纹深度下,当裂纹 位置靠近轧辊中部时,第一阶 方向位移模态分量 变化率在裂纹处出现明显尖峰,最大变化率在4 左右.当裂纹在靠近两端时,这种尖峰的现象不太明 显. 2.1.2裂纹损伤程度对轧辊位移模态分量变化率 的影响规律 分别选取裂纹深度 一0.02 m等6个裂纹深 度,计算在 =0.7 1TI处的第一阶 方向位移模态 分量变化率,其变化曲线如图3所示. 7 董6 鎏。4 鬟z 鎏1 —1 0 l 2 轴向坐b ̄x/m 图3第一阶z方向位移模态分量对于裂纹深度的变化率 Fig.3 Variational rate of first-order z-directional dis- placementmodal componentwith crack depth 从图3可以看出,在轧辊的相同位置,随着裂纹 深度的增加,轧辊第一阶 方向位移模态分量变化 率在裂纹处也增加,因此,可以根据第一阶位移模态 分量变化率大小大致判断损伤程度. 2.2裂纹损伤位置和损伤程度对轧辊应变模态分 量变化率的影响规律 2.2.1裂纹损伤位置对轧辊应变模态分量变化率 的影响规律 在轧辊上设置深度为z=0.08 ITI,圆心角为6O。 的裂纹.分别选取裂纹位置 一0.1 ITI等7个裂纹 位置,在轧辊工作段表面取一条母线,计算第一阶z 方向应变模态分量变化率,其变化曲线如图4所示. \ 褥 揪 嘲 鑫 1 轴向坐Jb ̄x/m 图4第一阶 方向应变模态分量对于裂纹位置的变化率 Fig.4 Variational rate of f'wst-order x-directional strain modal component with crack locations 从图4可以看出,在相同的裂纹深度下,第一阶 方向应变模态分量变化率曲线在裂纹位置产生尖 峰,且峰值接近150 .但是无论裂纹在什么地方, 在轧辊两端处都会出现较大的尖峰,会对裂纹位置 的判断造成影响. 2.2.2裂纹损伤程度对轧辊应变模态分量变化率 的影响规律 , 分别选取裂纹深度 一0.02 ITI等6个裂纹深 第5期 方治华等:基于振动模态分量的结构损伤定位方法 \诗 砸 i}}锹趟鑫 度,计算在x=0.7 ITI处的第一阶z方向应变模态 ∞如0如∞如加∞如 分量变化率,其变化曲线如图5所示. 0 l 2 轴向坐标 /m 图5第一阶 方向应变模态分■对于裂纹深度的变化率 5 Variation rate of rn'st-order x-directional strain modal campmlentwith crackdepth 从图5可以看出,在裂纹深度较浅时,随着裂纹 深度的增加,轧辊第一阶X方向应变模态分量变化 \褂 楼覃l星暴 鑫 捺 率在裂纹处逐渐增加,但是当裂纹增大到一定程度 时其裂纹处的应变模态分量变化率虽然也出现尖 峰,但其值反而减小了,同时在靠近裂纹一定范围内 应变模态分量变化率也明显增大了.按照这一规律 也可大致判断裂纹深度. 2.3裂纹损伤位置和损伤程度对轧辊第一阶模态 轴向位移差变化率的影响规律 2.3.1裂纹损伤位置对轧辊第一阶模态轴向位移 差变化率的影响规律 在轧辊上设置深度为z=0.08 ITI,圆心角为60。 的裂纹.分别选取裂纹位置z一0.1 m等7个裂纹 位置,计算第一阶模态下的轴向位移差变化率,其变 化曲线如图6所示. 0。 900 捻 600 墨 3oo j霉5 0 鑫 2.3.2裂纹损伤深度对轧辊第一阶模态轴向位移 差变化率的影响规律 分别选取裂纹深度 一0.02 rfl等6个裂纹深 度,计算在x=O.7 133处的第一阶模态下的轴向位 移差变化率,其变化曲线如图7所示. 轴向坐b ̄,'x/m 图7第一阶模态下的轴向位移差变化率对于裂纹深度 的变化率 F唔7 Variational rate of fast-order modal axial dis- placementdiferencewitII crack depth 从图7可以看出,在轧辊的相同位置,随着裂纹 深度的增加,轧辊第一阶模态下的轴向位移差变化 率在裂纹处也明显增加.该标识量在裂纹较小时也 变化明显,对裂纹最敏感. 3结论 三维实体结构振动时,结构上的任一点在三维 空间各个方向均产生模态分量.当结构损伤发生时, 这些模态分量会发生变化,但对损伤的敏感程度不 同.本文根据三维实体结构振动模态的这一特点,从 提高损伤标识量对损伤的敏感性出发,提出用对损 伤敏感的模态分量构造损伤标识量的思路,构造了 位移模态分量变化率、应变模态分量变化率和第一 阶模态下的轴向位移差变化率3个结构损伤标识 量,并对含表面裂纹轧辊的裂纹位置和裂纹深度与 3种标识量的关系进行了计算分析.研究不同位置 损伤对上述3种标识量的敏感程度和变化规律. 从计算结果可以看出,在存在裂纹的轧辊上,在 各种裂纹损伤程度下,第一阶 方向位移模态分量 变化率在裂纹处出现尖峰,能反映出裂纹的位置和 损伤程度,但它的变化率较小,对裂纹的敏感程度不 高;应变模态分量变化率也在裂纹处出现尖峰,且对 裂纹位置和损伤程度的敏感性要强于位移模态分量 变化率,但在两端变截面处也出现尖峰,会对裂纹位 置判断造成影响;第一阶模态下的轴向位移差变化 率在裂纹处出现尖峰,并且尖峰值明显大于前两种 ・120・ 兰州理工大学学报 第39卷 变化率标识量,因此,第一阶模态下的轴向位移差变 化率是几种参数中最敏感的损伤标识量. 致谢:本文得到内蒙古科技大学创新基金项目 (2009NC040)的资助,在此表示感谢. 参考文献: [1]孙晓丹,欧进萍.基于动力检测的损伤指标评价方法[J].振动 与冲击,2009,28(1):9-13. [2]顾培英,陈厚群。李同春,等.基于损伤应变模态的结构损伤识 别直接指标法口].自然科学进展,2007,17(2):240-247. [3]王俊荣.海洋平台结构物损伤检测与模型修正方法研究[D]. 北京:中国海洋大学,2009. [4]宋雨。项贻强.基于结构振动的桥梁损伤识别[J].振动、测 试与诊断,2005,25(3):222—226. [5] ROB1Ns0N N A,PEHE sON L D,JAMES G H,et a1.Dam— age detection in aircraft structures using dynamically measored static flexibility matrices[C]//Proc of the 14th International Modal Analysis Conference.1 ̄rbom,USA:SEM,1996:857- 865. . [63赵国文,赵启元,蔡振雄,等.船用柴油机曲轴裂纹的有限元模 态分析[J].舰船科学技术,2008,30(3):148-151. [7]阳建宏,罗桂荣,杨德斌,等.裂纹齿轮的模态分析及诊断[J]. 新技术新工艺,2008(3):50-52. [8]丁镇军,纪国宜.基于有限元法的疲劳裂纹叶片振动特性研究 [J].机械科学与技术,2010,29(10):1385-1388. 

因篇幅问题不能全部显示,请点此查看更多更全内容