发明名称 基于空间坐标监测的识别松弛的支承索的递进式方法
摘要 基于空间坐标监测的识别松弛的支承索的递进式方法,考虑到了被监测量的当前数值向量同被监测量的初始数值向量、虚拟单位损伤被监测量数值变化矩阵和当前名义虚拟损伤向量间的线性关系是近似的,本发明基于空间坐标监测,给出了使用线性关系分段逼近非线性关系的方法,将大区间分割成连续的一个个小区间,在每一个小区间内上述线性关系都是足够准确的,在每一个小区间内可以利用多目标优化算法等合适的算法算出当前索虚拟损伤向量的非劣解,据此可识别出虚拟受损索,在使用无损检测等方法从中鉴别出真实受损索后,剩下的虚拟受损索就是松弛的支承索,依据松弛程度同虚拟损伤程度间的关系就可确定需调整的索长。
申请公布号 CN101696910B 申请公布日期 2011.05.25
申请号 CN200910185352.8 申请日期 2009.11.05
申请人 东南大学 发明人 韩玉林
分类号 G01M99/00(2011.01)I;G01B21/00(2006.01)I;G01B21/20(2006.01)I;G01N19/08(2006.01)I;G01L5/04(2006.01)I;G06F19/00(2006.01)I 主分类号 G01M99/00(2011.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 叶连生
主权项 一种基于空间坐标监测的识别松弛的支承索的递进式健康监测方法,其特征是该方法包括:a.设共有N根索,首先确定索的编号规则,按此规则将索结构中所有的索编号,该编号在后续步骤中将用于生成向量和矩阵;b.确定指定的将被监测空间坐标的被监测点,给所有指定点编号;确定过每一监测点的将被监测的空间坐标分量,给所有指定的被监测空间坐标分量编号,“被监测空间坐标分量编号”在后续步骤中将用于生成向量和矩阵,“结构的全部被监测的空间坐标数据”由上述所有被监测空间坐标分量组成,将“结构的被监测的空间坐标数据”简称为“被监测量”,被监测点的数量不得小于索的数量;所有被监测空间坐标分量的数量之和不得小于索的数量;c.利用索的无损检测数据能够表达索的健康状态的数据建立初始虚拟损伤向量dio,其中i表示循环次数,后面i及上标i都表示循环次数,i=1,2,3,......;第一次循环时dio记为d1o;如果没有索的无损检测数据时,或者可以认为结构初始状态为无松弛、无损伤状态时,向量d1o的各元素数值取0;d.在建立初始虚拟损伤向量d1o的同时,直接测量计算得到索结构的所有被监测量的初始数值,组成被监测量的初始数值向量C1o;e.在建立初始虚拟损伤向量d1o和被监测量的初始数值向量C1o的同时,直接测量计算得到所有支承索的初始索力,组成初始索力向量Fo;同时,依据结构设计数据、竣工数据得到所有支承索的初始自由长度,组成初始自由长度向量lo;同时,依据结构设计数据、竣工数据或实测得到索结构的初始几何数据;同时,实测或根据结构设计、竣工资料得到所有索的弹性模量、密度、初始横截面面积;f.根据索结构的设计图、竣工图和索结构的上述实测数据,建立索结构的力学计算模型,基于该模型计算所得的计算数据同上述实测数据越接近越好,此时该模型被称为结构的力学计算基准模型A1,A1对应的索的健康状态由d1o描述;g.在力学计算基准模型Ai的基础上进行若干次力学计算,通过计算获得“名义虚拟单位损伤向量Diu”和“虚拟单位损伤被监测量数值变化矩阵ΔCi”;其中i表示循环次数,后面i及上标i都表示循环次数,i=1,2,3,......;Ai对应的索的健康状态由dio描述;h.实测得到索结构的所有支承索的当前索力,组成当前索力向量Fi;同时,实测得到索结构的所有指定被监测量的当前实测数值,组成“被监测量的当前数值向量Ci”;实测计算得到所有支承索的两个支承端点的水平距离,给本步及本步之前出现的所有向量的元素编号时,应使用同一编号规则,这样可以保证本步及本步之前和之后出现的各向量的、编号相同的元素,表示同一被监测量的、对应于该元素所属向量所定义的相关信息;i.定义待求的当前名义虚拟损伤向量dic和当前实际虚拟损伤向量di,损伤向量d1o、dic和di的元素个数等于索的数量,损伤向量的元素和索之间是一一对应关系,损伤向量的元素数值代表对应索的虚拟损伤程度或健康状态;j.依据“被监测量的当前数值向量Ci”同“被监测量的初始数值向量Cio”、“虚拟单位损伤被监测量数值变化矩阵ΔCi”和“当前名义虚拟损伤向量dic”间存在的近似线性关系,该近似线性关系可表达为式1,式1中除dic外的其它量均为已知,求解式1就可以算出当前名义虚拟损伤向量dic; <mrow> <msup> <mi>C</mi> <mi>i</mi> </msup> <mo>=</mo> <msubsup> <mi>C</mi> <mi>o</mi> <mi>i</mi> </msubsup> <mo>+</mo> <mi>&Delta;</mi> <msup> <mi>C</mi> <mi>i</mi> </msup> <mo>&CenterDot;</mo> <msubsup> <mi>d</mi> <mi>c</mi> <mi>i</mi> </msubsup> </mrow>式1k.利用式2表达的当前实际虚拟损伤向量di的元素dij同初始虚拟损伤向量dio的元素dioj和当前名义虚拟损伤向量dic的元素dicj间的关系,计算得到当前实际虚拟损伤向量di的所有元素, <mrow> <msubsup> <mi>d</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msubsup> <mi>d</mi> <mi>oj</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msubsup> <mi>d</mi> <mi>cj</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> </mrow>式2式2中j=1,2,3,……,N,由于当前实际虚拟损伤向量di的元素数值代表对应索的当前实际虚拟损伤程度,即实际松弛程度或实际损伤程度,当前实际虚拟损伤向量di中数值不为0的元素对应的支承索就是有问题的支承索,有问题的支承索可能是松弛索、也可能是受损索,其数值反应了松弛或损伤的程度;l.从第k步中识别出的有问题的支承索中鉴别出受损索,剩下的就是松弛索;m.利用在第k步获得的当前实际虚拟损伤向量di得到松弛索的当前实际虚拟损伤程度,利用在第h步获得的当前索力向量Fi,利用在第h步获得的所有支承索的两个支承端点的水平距离,利用在第e步获得的初始自由长度向量lo,利用在第e步获得的所有索的弹性模量、密度、初始横截面面积数据,通过将松弛索同受损索进行力学等效来计算松弛索的、与当前实际虚拟损伤程度等效的松弛程度,等效的力学条件是:一、两等效的索的无松弛和无损伤时的初始自由长度、几何特性参数、密度及材料的力学特性参数相同;二、松弛或损伤后,两等效的松弛索和损伤索的索力和变形后的总长相同,满足上述两个等效条件时,这样的两根支承索在结构中的力学功能就是完全相同的,即如果用等效的松弛索代替受损索后,索结构不会发生任何变化,反之亦然;依据前述力学等效条件求得那些被判定为松弛索的松弛程度,松弛程度就是支承索自由长度的改变量,也就是确定了那些需调整索力的支承索的索长调整量;这样就实现了支承索的松弛识别和损伤识别,计算时所需索力由当前索力向量Fi对应元素给出;n.在求得当前名义虚拟损伤向量dic后,按照式3建立标识向量Bi,式4给出了标识向量Bi的第j个元素的定义; <mrow> <msup> <mi>B</mi> <mi>i</mi> </msup> <mo>=</mo> <msup> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msubsup> <mi>B</mi> <mn>1</mn> <mi>i</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>B</mi> <mn>2</mn> <mi>i</mi> </msubsup> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msubsup> <mi>B</mi> <mi>j</mi> <mi>i</mi> </msubsup> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msubsup> <mi>B</mi> <mi>N</mi> <mi>i</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> </mrow>式3 <mrow> <msubsup> <mi>B</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> <mo>,</mo> </mtd> <mtd> <mi>if</mi> </mtd> <mtd> <msubsup> <mi>d</mi> <mi>cj</mi> <mi>i</mi> </msubsup> <mo>&lt;</mo> <msubsup> <mi>D</mi> <mi>uj</mi> <mi>i</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> <mo>,</mo> </mtd> <mtd> <mi>if</mi> </mtd> <mtd> <msubsup> <mi>d</mi> <mi>cj</mi> <mi>i</mi> </msubsup> <mo>&GreaterEqual;</mo> <msubsup> <mi>D</mi> <mi>uj</mi> <mi>i</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> </mrow>式4式3、式4中元素Bij是标识向量Bi的第j个元素,Diuj是名义虚拟单位损伤向量Diu的第j个元素,dicj是当前名义虚拟损伤向量dic的第j个元素,它们都表示第j根索的相关信息,式4中j=1,2,3,……,N;o.如果标识向量Bi的元素全为0,则回到第h步继续本次循环;如果标识向量Bi的元素不全为0,则进入下一步、即第p步;p.根据式5计算得到下一次、即第i+1次循环所需的初始虚拟损伤向量di+1o的每一个元素di+1oj; <mrow> <msubsup> <mi>d</mi> <mi>oj</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msubsup> <mi>d</mi> <mi>oj</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msubsup> <mi>D</mi> <mi>uj</mi> <mi>i</mi> </msubsup> <msubsup> <mi>B</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> </mrow>式5式5中Diuj是第i次循环名义虚拟单位损伤向量Diu的第j个元素,dioj是第i次循环当前名义虚拟损伤向量dio的第j个元素,Bij是第i次循环标识向量Bi的第j个元素;式5中j=1,2,3,……,N;q.在力学计算基准模型Ai的基础上,令索的健康状况为di+1o后更新得到下一次、即第i+1次循环所需的力学计算基准模型Ai+1;r.通过对力学计算基准模型Ai+1的计算得到对应于模型Ai+1的结构的所有被监测量的数值,这些数值组成下一次、即第i+1次循环所需的被监测量的初始数值向量Ci+1o;s.回到第g步,开始下一次循环。
地址 211109 江苏省南京市江宁开发区东南大学路2号