基于遗传算法的梁类构件损伤识别

1、基于遗传算法的梁类构件动力损伤基于遗传算法的梁类构件动力损伤识别方法研究识别方法研究Civil Engineering Department, College of Engineering, Ocean University of China2汇报内容汇报内容n1.1.选题的意义与研究现状选题的意义与研究现状n2.2.论文主要研究内容与创新点论文主要研究内容与创新点n3.3.梁类构件的动力学损伤识别模型梁类构件的动力学损伤识别模型n4.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究n5.5.结论与展望结论与展望Civil Engineering Department, College

2、 of Engineering, Ocean University of China31.1.选题的意义与研究现状选题的意义与研究现状1.1 1.1 选题的意义选题的意义 图1.1莫斯科水上乐园屋顶坍塌图1.2宜宾南门大桥桥面断裂 结构的损伤不能及时被发现并进行加固处理，则可能威结构的损伤不能及时被发现并进行加固处理，则可能威胁到结构的安全，甚至会使结构失效，带来不可估量的经胁到结构的安全，甚至会使结构失效，带来不可估量的经济损失和人员伤亡。济损失和人员伤亡。Civil Engineering Department, College of Engineering, Ocean Universi

3、ty of China41.1.选题的意义与研究现状选题的意义与研究现状1.2 1.2 研究现状研究现状动力检测方法动力检测方法模态参数模态参数直接比较直接比较模型修正模型修正技术技术模态参数模态参数衍生指标衍生指标固有固有频率频率阻阻尼尼振振型型计算计算智能智能神经神经网络网络小波小波 技术技术曲率曲率振型振型柔度柔度矩阵矩阵模态模态应变能应变能 遗传遗传算法算法图1.3常用的动力损伤识别方法整体检测与局整体检测与局部检测相结合。部检测相结合。缺点是：缺点是：不能够准确得到不能够准确得到 损伤位置和程度。损伤位置和程度。本文方法：本文方法：Civil Engineering Departme

4、nt, College of Engineering, Ocean University of China52.2.论文主要研究内容与创新点论文主要研究内容与创新点2.1 2.1 研究主要内容研究主要内容梁类构件两端弹性约束与中间单元。梁类构件两端弹性约束与中间单元。（1）研究对象）研究对象（2）研究目标）研究目标能够定位定量识别梁类构件两端弹性约束和中间单元的损伤情况。能够定位定量识别梁类构件两端弹性约束和中间单元的损伤情况。（3）研究理论基础）研究理论基础 Timoshenko梁模型梁模型 Timoshenko梁自由振动动力学方程梁自由振动动力学方程 遗传算法优化理论遗传算法优化理论 Ci

5、vil Engineering Department, College of Engineering, Ocean University of China62.2.论文主要研究内容与创新点论文主要研究内容与创新点图2.1利用遗传算法梁类构件损伤识别的基本结构图本本文文研研究究方方法法（4）Civil Engineering Department, College of Engineering, Ocean University of China72.2.论文主要研究内容与创新点论文主要研究内容与创新点2.2 2.2 论文的创新点论文的创新点 建立建立TimoshenkoTimoshenko梁类

6、构件两端弹性约束损伤识别理论模型。梁类构件两端弹性约束损伤识别理论模型。 推导出梁类构件中间某单元损伤时，在某阶振型下引起梁上各个节推导出梁类构件中间某单元损伤时，在某阶振型下引起梁上各个节 点应变模态的变化量的理论公式。点应变模态的变化量的理论公式。 利用一阶和二阶单元平均应变模态差的绝对值之和最小、三利用一阶和二阶单元平均应变模态差的绝对值之和最小、三 阶固有频率之差绝对值之和最小分别作为适应度函数，对弹阶固有频率之差绝对值之和最小分别作为适应度函数，对弹 性约束和中间单元不同位置、不同程度进行了损伤系数的优性约束和中间单元不同位置、不同程度进行了损伤系数的优 化识别识别结果较为准确化识别

7、识别结果较为准确。 Civil Engineering Department, College of Engineering, Ocean University of China83.3.梁类构件的动力学损伤识别模型梁类构件的动力学损伤识别模型3.1 Timoshenko3.1 Timoshenko梁损伤识别的动力学模型梁损伤识别的动力学模型图3.1 梁构件损伤识别的动力学模型（1 1）TimoshenkoTimoshenko梁两端梁两端 弹性约束损伤识别模型如弹性约束损伤识别模型如图（图（3.13.1）1234,k k k k 为梁两端的弹性约为梁两端的弹性约束，模拟临近构件结束，模拟临近构

8、件结点约束对本梁的影响点约束对本梁的影响情况情况120,0kk两端简支梁两端简支梁12,kk两端完全固结两端完全固结 123,4, 弹性约束损伤的程度弹性约束损伤的程度 1 12 23 34 4,kkkk表达结点的损伤状况表达结点的损伤状况 Civil Engineering Department, College of Engineering, Ocean University of China93.梁类构件的动力学损伤识别模型梁类构件的动力学损伤识别模型（2 2）弹性约束梁损伤识别模型的振型方程）弹性约束梁损伤识别模型的振型方程0sinhcoshcoshsinhsincoscossinco

9、shsinhsinhcoshcossinsincos344344344344222222222222333333211211DCBAEIkEIkEIkEIkEIkEIkEIkEIkkEIkEIEIkEIk0; 0:0; 0:03344222233332211dxdEIkdxdEIdxdklxdxdEIkdxdEIdxdkxxDxCxBxAxcoshsinhcossin)(2)4(22144ggh2)4(22144gghAGkmrag22424224)1 (aAGkrmh自由振动形状函数表达式自由振动形状函数表达式: 弹性约束边界条件弹性约束边界条件: （3-1） （3-2） （3-3） Civ

10、il Engineering Department, College of Engineering, Ocean University of China103.3.梁类构件的动力学损伤识别模型梁类构件的动力学损伤识别模型2212341234221231232212412422134134222342342(1c o s c o sh )() sin sin h ()(sin c o shc o s sin h )()(sin c o shc o s sin h )()(sin c o shc o s sin h )()(sik k k kk k kE Ik k kE Ik k kE Ik k

11、kE I 222222121224433222213132222141422222323244332424n c o shc o s sin h )()()sin sin h() ()2c o s c o sh() sin sin h ()()c o s c o sh()()c o s c o sh() ()2c o sk kE Ik kE Ik kE Ik kE Ik kE I 2222222234343332211333222232222333222244c o sh() sin sin h () ()sin sin h() ()(sin c o shc o s sin h )() ()(

12、sin c o shc o s sin h )() ()(c o s sin hsin c o sh )() (k kE IkE IkE IkE IkE I 4554422)(c o s sin hsin c o sh )() 2(1c o s c o sh )() sin sin h 0E Ilsinsin lcoscos lsinhsinh lcoshcosh (3) 弹性约束梁损伤识别模型频率方程弹性约束梁损伤识别模型频率方程如果（如果（3-3）有非零解则）有非零解则其系数矩阵的行列式为零。其系数矩阵的行列式为零。其中：其中：（3-4） Civil Engineering Departm

13、ent, College of Engineering, Ocean University of China113.3.梁类构件的动力学损伤识别模型梁类构件的动力学损伤识别模型由频率方程（3-4）24maEI由振型方程（3-3）表（3-1）梁的物理性质及其支座刚度表（3-2）求得固有频率和A，B，C，D的值Civil Engineering Department, College of Engineering, Ocean University of China123.3.梁类构件的动力学损伤识别模型梁类构件的动力学损伤识别模型1)1) 频率方程得到的固有频率系数的值频率方程得到的固有频率系数

14、的值2)2) 振型方程得到的振型曲线振型方程得到的振型曲线图3.3 函数的零点放大图( )yf a图3.2 函数 的图象 图3.4 第一阶理论解振型图3.5 第二阶理论解振型图3.6 第三阶理论解振型Civil Engineering Department, College of Engineering, Ocean University of China133.3.梁类构件的动力学损伤识别模型梁类构件的动力学损伤识别模型3.2 ANSYS建立弹性约束梁模型建立弹性约束梁模型 采用采用BEAM188单元单元进行划分进行划分20等分；一维等分；一维弹簧单元弹簧单元COMBIN14对两端弹性约束进

15、行模对两端弹性约束进行模拟。拟。图3.7 ANSYS模型第一阶振型表（3-3）对称约束梁的物理性质及其支座刚度图3.8 ANSYS模型第二阶振型 图3.9 ANSYS模型第三阶振型Civil Engineering Department, College of Engineering, Ocean University of China143.梁类构件的动力学损伤识别模型 说明了用两端弹性约束说明了用两端弹性约束TimoshenkoTimoshenko梁模型对结构中梁模型对结构中的梁构件进行假设建模具有合理性和可行性的梁构件进行假设建模具有合理性和可行性 ANSYS模型和理论解析解的振型曲线模

16、型和理论解析解的振型曲线图3.10 ANSYS模型和理论模型解析解第一阶振型图图3.11 ANSYS模型和理论模型解析解第二阶振型图图3.12 ANSYS模型和理论模型解析解第三阶振型图表（3-4）理论模型与ANSYS模型的频率值Civil Engineering Department, College of Engineering, Ocean University of China153.3.梁类构件的动力学损伤识别模型梁类构件的动力学损伤识别模型3.3 中间单元损伤识别理论公式中间单元损伤识别理论公式,(1),(1),11,(1),(1)21()()(2)TTNNjudtrtitrtij

17、i tji tji tr jr jr jr jrrjjrijriTNtrtir jr jr jrjrihKKlllhKil；（）mttjiji1, 单元发生损伤在第单元发生损伤在第 阶振型下阶振型下 节点处的应变模态差节点处的应变模态差 :tij对于多处损伤可由单处损伤进行线性叠加求得对于多处损伤可由单处损伤进行线性叠加求得: （3-5） （3-6） Civil Engineering Department, College of Engineering, Ocean University of China164.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究4.1 弹性约束损伤系数

18、的多元多峰值优化目标函数弹性约束损伤系数的多元多峰值优化目标函数min( ( )min( )JE011234Test为损伤前后所有单元的测量平均应变模态值的向量为损伤前后所有单元的测量平均应变模态值的向量12,.TestTestTestTestTest Tkn即：即：( ) |( )|TestcalE其中：其中：( )cal为损伤前后所有单元的理论平均应变模态值的向量为损伤前后所有单元的理论平均应变模态值的向量12( )( ),( ).( ).( )calcalcalcalcalTkn 即：即： 为损伤单元的刚度影响系数为损伤单元的刚度影响系数， 为将梁构件划分的单元数。为将梁构件划分的单元数

19、。 nCivil Engineering Department, College of Engineering, Ocean University of China174.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究4.2 中间损伤单元应变残差改变量的多目标优中间损伤单元应变残差改变量的多目标优 化目标函数化目标函数min( )min( )kkJE1,2,.,kn01( ) |TestcalkkkE 为优化子目标为优化子目标 。,Testcalkk为损伤前后单元为损伤前后单元 的测量平均应变改变量的测量平均应变改变量k与理论平均应变改变量与理论平均应变改变量 。 为损伤单元的刚度影

20、响系数为损伤单元的刚度影响系数 。 为梁划分的单元数。为梁划分的单元数。 123,.,.Tkn nCivil Engineering Department, College of Engineering, Ocean University of China184.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究4.3 遗传算法设置遗传算法设置（1）编码：本文采用二进制编码）编码：本文采用二进制编码 ，损伤识别的计算精度为，损伤识别的计算精度为0.001 。（2）初始种群：）初始种群：80120不等。不等。（3）适应度函数）适应度函数：20, ,( )|( )|TestCaljj ij

21、iif1j 123420, ,( )|( )|TestCaljj ij iif2j 12343,1( )|( )|testcaljjjfff1234（4）交叉概率：（根据具体情况，有时候需要改变）本文用多点交叉）交叉概率：（根据具体情况，有时候需要改变）本文用多点交叉 （5）终止准则：）终止准则： 达到规定迭代的最大数目。达到规定迭代的最大数目。 最优个体的适应度值在一个位置停留次数最优个体的适应度值在一个位置停留次数50代或相差很小代或相差很小Civil Engineering Department, College of Engineering, Ocean University of C

22、hina194.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究（1）function sresult=myGAs(n,a,b,pc,pm,e,z,sumstep)（2）function sresult= plmyGAs (n,a,b,pc,pm,e,z,sumstep) myGAs为编制的matlab程序函数名；是适用于第一、第二种目标函 数遗传算法程序。 plmyGAs为编制的matlab程序函数名；是适用于第三种目标函数遗传算法程序。 其中n群体规模；a搜索上限；b搜索下限；pc交叉概率；pm变异概率；e计算精度；z-表示前z阶数之和；sumstep表示计算的终止代数。 Civi

23、l Engineering Department, College of Engineering, Ocean University of China204.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究4.4 4.4 两端弹性约束不同损伤情况下的识别结果两端弹性约束不同损伤情况下的识别结果 没有发生损伤没有发生损伤Civil Engineering Department, College of Engineering, Ocean University of China214.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究发生一处损伤采用适发生一处损伤采用适应度函数一识别

24、结果应度函数一识别结果Civil Engineering Department, College of Engineering, Ocean University of China224.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究发生一处损伤采用适发生一处损伤采用适应度函数二识别结果应度函数二识别结果Civil Engineering Department, College of Engineering, Ocean University of China234.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究发生一处损伤采用适发生一处损伤采用适应度函数三识别结果应度函数

25、三识别结果Civil Engineering Department, College of Engineering, Ocean University of China244.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究（1）一处损伤识别结果总结）一处损伤识别结果总结 采用一、二阶单元平均应变模态差的绝对值之和做为采用一、二阶单元平均应变模态差的绝对值之和做为 适应度函数，能够很精确的定位定量。适应度函数，能够很精确的定位定量。 采用三阶固有频率之差绝对值之和作为适应度函数，能采用三阶固有频率之差绝对值之和作为适应度函数，能 够精确的定位，大致能够定量。够精确的定位，大致能够定量。

26、 采用单元平均应变模态差的绝对值之和的适应度函数，采用单元平均应变模态差的绝对值之和的适应度函数， 比采用固有频率作为适应的函数优化识别效果好比采用固有频率作为适应的函数优化识别效果好Civil Engineering Department, College of Engineering, Ocean University of China254.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究发生两处损伤采用适发生两处损伤采用适应度函数一识别结果应度函数一识别结果Civil Engineering Department, College of Engineering, Ocean U

27、niversity of China264.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究发生两处损伤采用适发生两处损伤采用适应度函数二识别结果应度函数二识别结果Civil Engineering Department, College of Engineering, Ocean University of China274.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究发生两处损伤采用适发生两处损伤采用适应度函数三识别结果应度函数三识别结果Civil Engineering Department, College of Engineering, Ocean Univers

28、ity of China284.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究（2）二处损伤识别结果总结）二处损伤识别结果总结 采用一、二阶单元平均应变模态差的绝对值之和做为采用一、二阶单元平均应变模态差的绝对值之和做为 适应度函数，能够很精确的定位定量。适应度函数，能够很精确的定位定量。 采用三阶固有频率之差绝对值之和作为适应度函数，能采用三阶固有频率之差绝对值之和作为适应度函数，能 够精确的定位，大致能够定量。够精确的定位，大致能够定量。 采用单元平均应变模态差的绝对值之和的适应度函数，采用单元平均应变模态差的绝对值之和的适应度函数， 比采用固有频率作为适应的函数优化识别效果好比

29、采用固有频率作为适应的函数优化识别效果好Civil Engineering Department, College of Engineering, Ocean University of China294.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究发生三处损伤采用适发生三处损伤采用适应度函数一识别结果应度函数一识别结果Civil Engineering Department, College of Engineering, Ocean University of China304.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究发生三处损伤采用适发生三处损伤采用适应度函数

30、二识别结果应度函数二识别结果Civil Engineering Department, College of Engineering, Ocean University of China314.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究发生三处损伤采用适发生三处损伤采用适应度函数三识别结果应度函数三识别结果Civil Engineering Department, College of Engineering, Ocean University of China324.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究（3）三处损伤识别结果总结）三处损伤识别结果总结 采用一

31、、二阶单元平均应变模态差的绝对值之和做为采用一、二阶单元平均应变模态差的绝对值之和做为 适应度函数，基本能够进行定位定量，但精度不高。适应度函数，基本能够进行定位定量，但精度不高。 采用三阶固有频率之差绝对值之和作为适应度函数，采用三阶固有频率之差绝对值之和作为适应度函数，大大 致能够定位，定量，但是精度不高。致能够定位，定量，但是精度不高。Civil Engineering Department, College of Engineering, Ocean University of China334.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究4.5 中间单元的不同损伤情况下的

32、识别结果中间单元的不同损伤情况下的识别结果 发生两处损伤采用适发生两处损伤采用适应度函数一识别结果应度函数一识别结果Civil Engineering Department, College of Engineering, Ocean University of China344.4.基于遗传算法的损伤识别研究基于遗传算法的损伤识别研究发生两处损伤采用适发生两处损伤采用适应度函数二识别结果应度函数二识别结果Civil Engineering Department, College of Engineering, Ocean University of China354.4.基于遗传算法的损伤识

33、别研究基于遗传算法的损伤识别研究4.5 4.5 中间单元的不同损伤情况下的识别结果总结中间单元的不同损伤情况下的识别结果总结 以第一、二阶单元应变模态值改变量之差的绝对值为最小以第一、二阶单元应变模态值改变量之差的绝对值为最小 的多目标优化遗传算法，对一、二、三处损伤能够较好识的多目标优化遗传算法，对一、二、三处损伤能够较好识 别出损伤的位置与程度。别出损伤的位置与程度。 在损伤识别的时候，要注意检测单元避开，各阶应变模态节点。在损伤识别的时候，要注意检测单元避开，各阶应变模态节点。Civil Engineering Department, College of Engineering, Ocean University of China365.5.结论与展望结论与展望5.1 5.1 结论结论 通过弹性约束损伤识别模型，可以明确的得到，损伤的通过弹性约束损伤识别模型，可以明确的得到，损伤的 弹性约束对固有频率、振型的影响。弹性约束对固有频率、振型的影响。 得到梁类构件，中间