应变能方法进行结构损伤识与模拟计算.docVIP

应变能方法进行结构损伤识与模拟计算 　　结构的无损检测是在不损伤结构中材料的前提下，将结构中得损伤或缺陷等检测出来。目前，对结构的无损检测方法主要有： 　　频率相关分析 　　振型模态相关分析 　　响应函数模态曲率 　　柔度改变 　　应变能 　　这些无损检测方法是基于对结构的动力学特性（质量、刚度、阻尼）的变化跟踪分析，与在没有损伤的理论值进行相关分析，由此来判断结构的损伤部位及程度。由于结构的动力学特性是结构的整体特性的集中反映，因此，结构损伤的无损检测方法是一种结构整体损伤检测方法。无损检测方法一般是从结构的频率响应或者结构的模态响应进行分析的。尽管人们在结构的频率响应对结构的损伤进行检测方面进行了大量的工作，但是由于结构中的缺陷或损伤对结构的频率影响很小，因此实际效果不佳。 　　Stubbs于1995年提出了利用结构的应变能变化来确定结构的损伤，由于其方程简练、物理意义明确，受到广泛的关注。本文根据Stubbs的理论，给出了一个结构损伤判断依据，并对一个简支梁在不同工况下进行了数值模拟计算。 应变能理论 　　对于平面梁，如果梁的挠曲线函数为 　　对于第阶模态的振型，梁内相应的应变能为： 　　当梁被划分为N个单元时（如图 1），则第j个单元对应于第i阶模态，相应的应变能为： 　　图 1 梁单元划分示意图 　　因而，对于第i阶模态，j单元的应变能与梁的应变能的比值（应变能分数）为： 　　显然 　　同理，对于受损的梁结构，如第i阶模态的振型为，则梁的应变能： 　　对于第i阶模态，单元j内相应的应变能： 　　对于受损结构的第i阶模态，单元j的应变能与梁的应变能比值（应变能分数）为： 　　显然 　　如果选取的梁单元相对较小，则第j个单元的弯曲刚度可以近似认为是常量，可以提到积分号外，即： 　　根据损伤局部化的概念，对于没受损的单元j，其应变能分数应相等，即： 　　这样，对于梁内受损的单元k来说，由上面两个可以得出受损单元k的应变能分数也相等，即： 　　从而： 　　由于损伤往往是发生在结构中的某一局部区域，可以以假定损伤前和损伤后，梁总得抗弯刚度不变，则可以用下式来表达受损单元k的刚度变化情况： 　　上式称为结构单元的损伤判据，它反映了单元受损前和受损后，梁单元的刚度变化情况。是对于第i阶模态，单元k的损伤指数。 　　如果考虑前m阶模态，结构的单元的损伤指数： 　　当，表示单元k无损，当，表示单元k受损。 　　2 数值结果分析 　　取一矩形截面简支梁，跨径为，横截面尺寸为，梁被划分为10个单元，本文选择四种工况分别进行分析比较： 工况一：结构无损; 工况二：单元2刚度受损50%; 工况三：单元5刚度受损50%; 工况四：单元5刚度受损80%。 　　对以上个工况计算其固有频率、振型以及各单元的损伤指数，各工况的固有频率如表 1所示;振型计算结果如图 2所示; 　　表 1各工况固有频率 　　阶数 　　无损 　　单元2 　　受损50% 　　单元5 　　受损50% 　　单元5 　　受损80% 　　1 　　0.6981 　　0.6837 　　0.6388 　　0.5232 　　2 　　2.7920 　　2.6272 　　2.7622 　　2.7017 　　3 　　6.2824 　　5.8330 　　5.9142 　　5.4048 　　4 　　11.1760 　　10.5690 　　10.8470 　　10.3350 　　5 　　17.4900 　　16.9070 　　16.8420 　　15.9060 　　6 　　25.2660 　　24.6060 　　24.2740 　　22.9690 　　7 　　28.8970 　　26.4670 　　27.3370 　　23.7160 　　8 　　34.5700 　　33.4400 　　33.3990 　　31.0510 　　9 　　45.4960 　　43.8280 　　43.9290 　　42.3930 　　10 　　58.1160 　　56.4520 　　55.3860 　　50.9420 　　图 2各阶振型图 　　选取前5阶振型计算各单元的损伤指数计算结果如表 2所示，损伤指数变化图如图 3所示。 　　表 2各工况损伤指数计算结果 　　单元 　　工况二 　　损伤指数 　　工况三 　　损伤指数 　　工况是 　　损伤指数 　　1 　　1.0254 　　0.9845 　　0.9954 　　2 　　1.2854 　　0.9825 　　0.9965 　　3 　　1.1251 　　0