发明名称 基于混合监测的识别松弛的支承索的递进式方法
摘要 基于混合监测的识别松弛的支承索的递进式方法用于识别松弛的支承索时,考虑到了被监测量的当前数值向量同被监测量的初始数值向量、虚拟单位损伤被监测量数值变化矩阵和当前名义虚拟损伤向量间的线性关系是近似的,本发明基于混合监测,给出了使用线性关系分段逼近非线性关系的方法,将大区间分割成连续的一个个小区间,在每一个小区间内上述线性关系都是足够准确的,在每一个小区间内可以利用多目标优化算法等合适的算法算出当前索虚拟损伤向量的非劣解,据此可识别出虚拟受损索,在使用无损检测等方法从中鉴别出真实受损索后,剩下的虚拟受损索就是松弛的支承索,依据松弛程度同虚拟损伤程度间的关系就可确定需调整的索长。
申请公布号 CN101714189B 申请公布日期 2012.05.09
申请号 CN200910185353.2 申请日期 2009.11.05
申请人 东南大学 发明人 韩玉林
分类号 G06F19/00(2011.01)I;G01D21/02(2006.01)I;E01D19/16(2006.01)I 主分类号 G06F19/00(2011.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 叶连生
主权项 1.一种基于混合监测的识别松弛的支承索的递进式方法,其特征是该方法包括:a.设共有N根索,首先确定索的编号规则,按此规则将索结构中所有的索编号,该编号在后续步骤中将用于生成向量和矩阵;b.确定混合监测时指定的将被监测索力的支承索,设索系统中共有N根索,结构的被监测的索力数据由结构上M<sub>1</sub>个指定索的M<sub>1</sub>个索力数据来描述,结构索力的变化就是所有指定索的索力的变化,每次共有M<sub>1</sub>个索力测量值或计算值来表征结构的索力信息,M<sub>1</sub>是一个不小于0的整数,确定混合监测时指定的将被监测应变的被测量点,结构的被监测的应变数据可由结构上K<sub>2</sub>个指定点的、及每个指定点的L<sub>2</sub>个指定方向的应变来描述,结构应变数据的变化就是K<sub>2</sub>个指定点的所有被测应变的变化,每次共有M<sub>2</sub>个应变测量值或计算值来表征结构应变,M<sub>2</sub>为K<sub>2</sub>和L<sub>2</sub>之积,M<sub>2</sub>是不小于0的整数,确定混合监测时指定的将被监测角度的被测量点,结构的被监测的角度数据由结构上K<sub>3</sub>个指定点的、过每个指定点的L<sub>3</sub>个指定直线的、每个指定直线的H<sub>3</sub>个角度坐标分量来描述,结构角度的变化就是所有指定点的、所有指定直线的、所有指定的角度坐标分量的变化,每次共有M<sub>3</sub>个角度坐标分量测量值或计算值来表征结构的角度信息,M<sub>3</sub>为K<sub>3</sub>、L<sub>3</sub>和H<sub>3</sub>之积,M<sub>3</sub>是一个不小于0的整数,确定混合监测时指定的将被监测的形状数据,结构的被监测的形状数据由结构上K<sub>4</sub>个指定点的、及每个指定点的L<sub>4</sub>个指定方向的空间坐标来描述,结构形状数据的变化就是K<sub>4</sub>个指定点的所有坐标分量的变化,每次共有M<sub>4</sub>个坐标测量值或计算值来表征结构形状,M<sub>4</sub>为K<sub>4</sub>和L<sub>4</sub>之积,M<sub>4</sub>是一个不小于0的整数,综合上述混合监测的被监测量,整个结构共有M个被监测量,M为M<sub>1</sub>、M<sub>2</sub>、M<sub>3</sub>和M<sub>4</sub>之和,定义参量K,K为M<sub>1</sub>、K<sub>2</sub>、K<sub>3</sub>和K<sub>4</sub>之和,K和M不得小于索的数量N,由于M个被监测量是不同类型的,称为“基于混合监测的识别松弛的支承索的递进式方法”,将本步所列出的“混合监测时结构的被监测的所有参量”简称为“被监测量”;c.利用索的无损检测数据等能够表达索的健康状态的数据建立初始虚拟损伤向量d<sup>i</sup><sub>o</sub>,其中i表示循环次数,后面i及上标i都表示循环次数,i=1,2,3,......;第一次循环时d<sup>i</sup><sub>o</sub>记为d<sup>1</sup><sub>o</sub>;如果没有索的无损检测数据及其他能够表达索的健康状态的数据时,或者可以认为结构初始状态为无松弛、无损伤状态时,向量d<sup>1</sup><sub>o</sub>的各元素数值取0;d.在建立初始虚拟损伤向量d<sup>1</sup><sub>o</sub>的同时,直接测量计算得到索结构的所有被监测量的初始数值,组成被监测量的初始数值向量C<sup>1</sup><sub>o</sub>;e.在建立初始虚拟损伤向量d<sup>1</sup><sub>o</sub>和被监测量的初始数值向量C<sup>1</sup><sub>o</sub>的同时,直接测量计算得到所有支承索的初始索力,组成初始索力向量F<sub>o</sub>;同时,依据结构设计数据、竣工数据得到所有支承索的初始自由长度,组成初始自由长度向量l<sub>o</sub>;同时,依据结构设计数据、竣工数据或实测得到索结构的初始几何数据;同时,实测或根 据结构设计、竣工资料得到所有索的弹性模量、密度、初始横截面面积;f.根据索结构的设计图、竣工图和索结构的上述实测数据,建立索结构的力学计算模型,基于该模型计算所得的计算数据同上述实测数据越接近越好,此时该模型被称为结构的力学计算基准模型A<sup>1</sup>,A<sup>1</sup>对应的索的健康状态由d<sup>1</sup><sub>o</sub>描述;g.在力学计算基准模型A<sup>i</sup>的基础上进行若干次力学计算,通过计算获得“名义虚拟单位损伤向量D<sup>i</sup><sub>u</sub>”和“虚拟单位损伤被监测量数值变化矩阵ΔC<sup>i</sup>”,其中i表示循环次数,后面i及上标i都表示循环次数,i=1,2,3,......,A<sup>i</sup>对应的索的健康状态由d<sup>i</sup><sub>o</sub>描述;h.实测得到索结构的所有支承索的当前索力,组成当前索力向量F<sup>i</sup>;同时,实测得到索结构的所有指定被监测量的当前实测数值,组成“被监测量的当前数值向量C<sup>i</sup>”,实测计算得到所有支承索的两个支承端点的水平距离,给本步及本步之前出现的所有向量的元素编号时,应使用同一编号规则,这样可以保证本步及本步之前和之后出现的各向量的、编号相同的元素,表示同一被监测量的、对应于该元素所属向量所定义的相关信息;i.定义待求的当前名义虚拟损伤向量d<sup>i</sup><sub>c</sub>和当前实际虚拟损伤向量d<sup>i</sup>,损伤向量d<sup>1</sup><sub>o</sub>、d<sup>i</sup><sub>c</sub>和d<sup>i</sup>的元素个数等于索的数量,损伤向量的元素和索之间是一一对应关系,损伤向量的元素数值代表对应索的虚拟损伤程度或健康状态;j.依据“被监测量的当前数值向量C<sup>i</sup>”同“被监测量的初始数值向量C<sup>i</sup><sub>o</sub>”、“虚拟单位损伤被监测量数值变化矩阵ΔC<sup>i</sup>”和“当前名义虚拟损伤向量d<sup>i</sup><sub>c</sub>”间存在的近似线性关系,该近似线性关系可表达为式1,式1中除d<sup>i</sup><sub>c</sub>外的其它量均为已知,求解式1就可以算出当前名义虚拟损伤向量d<sup>i</sup><sub>c</sub>;<img file="FSB00000476571800021.GIF" wi="357" he="58" />式1k.利用式2表达的当前实际虚拟损伤向量d<sup>i</sup>的元素d<sup>i</sup><sub>j</sub>同初始虚拟损伤向量d<sup>i</sup><sub>o</sub>的元素d<sup>i</sup><sub>oj</sub>和当前名义虚拟损伤向量d<sup>i</sup><sub>c</sub>的元素d<sup>i</sup><sub>cj</sub>间的关系,计算得到当前实际虚拟损伤向量d<sup>i</sup>的所有元素,<img file="FSB00000476571800022.GIF" wi="492" he="70" />式2式2中j=1,2,3,......,N,由于当前实际虚拟损伤向量d<sup>i</sup>的元素数值代表对应索的当前实际虚拟损伤程度,即实际松弛程度或实际损伤程度,当前实际虚拟损伤向量d<sup>i</sup>中数值不为0的元素对应的支承索就是有问题的支承索,有问题的支承索可能是松弛索、也可能是受损索,其数值反应了松弛或损伤的程度;l.从第k步中识别出的有问题的支承索中鉴别出受损索,剩下的就是松弛索;m.利用在第k步获得的当前实际虚拟损伤向量d<sup>i</sup>得到松弛索的当前实际虚拟损伤程度,利用在第h步获得的当前索力向量F<sup>i</sup>,利用在第h步获得的所有支承索的两个支承端点的水平距离,利用在第e步获得的初始自由长度向量l<sub>o</sub>,利用在第e步获得的所有索的弹性模量、密度、初始横截面面积数据,通过将松弛索同受损索进行力学等效来计算松弛索的、与当前实际虚拟损伤程度等效的松弛程度,等效的力学条件是:一、两等效的索的无松弛和无损伤时的初始自由长度、几何特性参数、密度及材料的力学特性参数相同;二、松弛或损伤后,两等效的松弛索和损伤 索的索力和变形后的总长相同,满足上述两个等效条件时,这样的两根支承索在结构中的力学功能就是完全相同的,即如果用等效的松弛索代替受损索后,索结构不会发生任何变化,反之亦然,依据前述力学等效条件求得那些被判定为松弛索的松弛程度,松弛程度就是支承索自由长度的改变量,也就是确定了那些需调整索力的支承索的索长调整量,这样就实现了支承索的松弛识别和损伤识别,计算时所需索力由当前索力向量F<sup>i</sup>对应元素给出,n.在求得当前名义虚拟损伤向量d<sup>i</sup><sub>c</sub>后,按照式3建立标识向量B<sup>i</sup>,式4给出了标识向量B<sup>i</sup>的第j个元素的定义;<img file="FSB00000476571800031.GIF" wi="902" he="65" />式3<img file="FSB00000476571800032.GIF" wi="457" he="150" />式4式3、式4中元素B<sup>i</sup><sub>j</sub>是标识向量B<sup>i</sup>的第j个元素,D<sup>i</sup><sub>uj</sub>是名义虚拟单位损伤向量D<sup>i</sup><sub>u</sub>的第j个元素,d<sup>i</sup><sub>cj</sub>是当前名义虚拟损伤向量d<sup>i</sup><sub>c</sub>的第j个元素,它们都表示第j根索的相关信息,式4中j=1,2,3,......,N;o.如果标识向量B<sup>i</sup>的元素全为0,则回到第h步继续本次循环;如果标识向量B<sup>i</sup>的元素不全为0,则进入下一步、即第p步;p.根据式5计算得到下一次、即第i+1次循环所需的初始虚拟损伤向量d<sup>i+1</sup><sub>o</sub>的每一个元素d<sup>i+1</sup><sub>oj</sub>;<img file="FSB00000476571800033.GIF" wi="539" he="66" />式5式5中D<sup>i</sup><sub>uj</sub>是第i次循环名义虚拟单位损伤向量D<sup>i</sup><sub>u</sub>的第j个元素,d<sup>i</sup><sub>cj</sub>是第i次循环当前名义虚拟损伤向量d<sup>i</sup><sub>c</sub>的第j个元素,B<sup>i</sup><sub>j</sub>是第i次循环标识向量B<sup>i</sup>的第j个元素,式5中j=1,2,3,......,N;q.在力学计算基准模型A<sup>i</sup>的基础上,令索的健康状况为d<sup>i+1</sup><sub>o</sub>后更新得到下一次、即第i+1次循环所需的力学计算基准模型A<sup>i+1</sup>;r.通过对力学计算基准模型A<sup>i+1</sup>的计算得到对应于模型A<sup>i+1</sup>的结构的所有被监测量的数值,这些数值组成下一次、即第i+1次循环所需的被监测量的初始数值向量C<sup>i+1</sup><sub>o</sub>;s.回到第g步,开始下一次循环;在步骤g中,获得“名义虚拟单位损伤向量D<sup>i</sup><sub>u</sub>”和“虚拟单位损伤被监测量数值变化矩阵ΔC<sup>i</sup>”的方法为:g1.在结构的力学计算基准模型A<sup>i</sup>的基础上进行若干次力学计算,计算次数数值上等于所有索的数量,有N根索就有N次计算,每一次计算假设索系统中只有一根索在原有虚拟损伤的基础上再增加虚拟单位损伤,每一次计算中出现虚拟单位损伤的索不同于其它次计算中出现虚拟单位损伤的索,并且每一次假定有虚拟单位损伤的索的虚拟单位损伤值可以不同于其他索的虚拟单位损伤值,用“名义虚拟单位损伤向量D<sup>i</sup><sub>u</sub>”记录所有索的假定的单位损伤,每一次计算得到所有被监测量的当前数值,每一次计算得到的所有被监测量的当前数值组成一个“被监测量的计算当前数值向量”;当假设第j根索有单位损伤时,可用C<sup>i</sup><sub>tj</sub>表示对应的“被监测量的当 前计算数值向量C<sup>i</sup><sub>tj</sub>”,在本步骤中给各向量的元素编号时,应同本方法中其它向量使用同一编号规则,这样可以保证本步骤中各向量中的任意一个元素,同其它向量中的、编号相同的元素,表达了同一被监测量或同一对象的相关信息;g2.每一次计算得到的那个“被监测量的当前计算数值向量C<sup>i</sup><sub>tj</sub>”减去“被监测量的初始数值向量C<sup>i</sup><sub>o</sub>”得到一个向量,再将该向量的每一个元素都除以本次计算中假定的虚拟单位损伤值后得到一个“被监测量的数值变化向量”;有N根索就有N个“被监测量的数值变化向量”;g3.由这N个“被监测量的数值变化向量”依次组成有N列的“虚拟单位损伤被监测量数值变化矩阵ΔC<sup>i</sup>”;“虚拟单位损伤被监测量数值变化矩阵”的列的编号规则与当前名义虚拟损伤向量d<sup>i</sup><sub>c</sub>和当前实际虚拟损伤向量d<sup>i</sup>的元素编号规则相同。 
地址 211109 江苏省南京市江宁开发区东南大学路2号