[发明专利]基于索力监测识别受损索松弛索支座广义位移的递进方法

摘要

基于索力监测识别受损索松弛索支座广义位移的递进方法，即对全部支承索和人为增加的索的索力进行监测，考虑到了被监测量的当前数值向量同被监测量的初始数值向量、单位损伤被监测量变化矩阵和当前名义健康状态向量间的线性关系是近似的，为克服此缺陷，本发明给出了使用线性关系分段逼近非线性关系的方法，将大区间分割成连续的一个个小区间，在每一个小区间内上述线性关系都是足够准确的，在每一个小区间内可以利用多目标优化算法等合适的算法快速识别出支座广义位移、受损索和松弛索。

主权项

一种基于索力监测识别受损索松弛索支座广义位移的递进方法，其特征在于所述方法包括：a.为叙述方便起见，统一称被评估的支承索和支座广义位移分量为被评估对象，设被评估的支承索的数量和支座广义位移分量的数量之和为N，即被评估对象的数量为N；确定被评估对象的编号规则，按此规则将索结构中所有的被评估对象编号，该编号在后续步骤中将用于生成向量和矩阵；用变量j表示这一编号，j＝1，2，3，...，N；b.设索系统中共有Q根支承索，结构索力数据包括这Q根支承索的索力，显然Q小于被评估对象的数量N；仅仅通过Q个支承索的Q个索力数据来求解未知的N个被评估对象的状态是不可能的，在监测全部Q根支承索索力的基础上，在结构上人为增加M2根索，在结构健康监测过程中将监测这新增加的M2根索的索力；综合上述被监测量，整个结构共有M根索的M个索力被监测，即有M个被监测量，其中M为Q与M2之和；M不得小于被评估对象的数量N；新增加的M2根索的刚度同索结构的任意一根支承索的刚度相比，应当小得多；新增加的M2根索的索力应当比索结构的任意一根支承索的索力小得多，这样可以保证即使这新增加的M2根索出现了损伤或松弛，对索结构其他构件的应力、应变、变形的影响微乎其微；新增加的M2根索的横截面上正应力应当小于其疲劳极限，这些要求可以保证新增加的M2根索不会发生疲劳损伤；新增加的M2根索的两端应当充分锚固，保证不会出现松弛；新增加的M2根索应当得到充分的防腐蚀保护，保证新增加的M2根索不会发生损伤和松弛；为方便起见，将“结构的被监测的所有参量”简称为“被监测量”；给M个被监测量连续编号，该编号在后续步骤中将用于生成向量和矩阵；c.利用被评估对象的无损检测数据等能够表达被评估对象的健康状态的数据建立被评估对象初始健康状态向量dio；如果没有被评估对象的无损检测数据时，向量dio的各元素数值取0；向量dio的元素的编号规则和被评估对象的编号规则相同；本发明用i表示循环次数，i＝1，2，3，......；这里是第一次循环，i取1，即这里建立的初始健康状态向量dio可以具体化为d1o；d.在建立初始健康状态向量d1o的同时，直接测量计算得到索结构的所有被监测量的初始数值，组成被监测量的初始数值向量Cio；这里是第一次循环，i取1，即这里建立的被监测量的初始数值向量Cio可以具体化为C1o；在实测得到被监测量初始数值向量C1o的同时，实测得到索结构的初始几何数据和初始索结构支座广义坐标数据；直接测量计算得到所有支承索的初始索力，组成初始索力向量Fo；同时，依据结构设计数据、竣工数据得到所有支承索的初始自由长度，组成初始自由长度向量Io；向量Fo和向量Io是不变的；同时，实测或根据结构设计、竣工资料得到所有索的弹性模量、密度、初始横截面面积；支座广义坐标包括线量和角量两种；e.根据索结构的设计图、竣工图和索结构的实测数据、索的无损检测数据和初始索结构支座广义坐标数据建立索结构的力学计算基准模型Ai；这里是第一次循环，i取1，即这里建立的索结构的力学计算基准模型Ai可以具体化为A1；f.在力学计算基准模型Ai的基础上进行若干次力学计算，通过计算获得“单位损伤被监测量数值变化矩阵ΔCi”和“名义单位损伤向量Diu”；g.实测得到索结构的所有指定被监测量的当前实测数值，组成“被监测量的当前数值向量Ci”；给本步及本步之前出现的所有向量的元素编号时，应使用同一编号规则，这样可以保证本步及本步之前出现的各向量的、编号相同的元素，表示同一被监测量的、对应于该元素所属向量所定义的相关信息；实测得到索结构的所有支承索的当前索力，组成当前索力向量Fi；实测计算得到所有支承索的两个支承端点的空间坐标，两个支承端点的空间坐标在水平方向分量的差就是两个支承端点水平距离；h.在结构健康监测过程中，对新增加的M2根索进行无损检测，从中鉴别出出现损伤或松弛的索；i.依据被监测量编号规则，从被监测量的初始数值向量Cio中去除步骤h中鉴别出的出现损伤或松弛的索对应的元素；依据被监测量编号规则，从单位损伤被监测量数值变化矩阵ΔCi中去除步骤h中鉴别出的出现损伤或松弛的索对应的行；j.定义当前名义健康状态向量dic和当前实际健康状态向量di，两个损伤向量的元素个数等于被评估对象的数量，当前名义健康状态向量dic的元素数值代表对应被评估对象的当前名义损伤程度或支座广义位移，当前实际健康状态向量di的元素数值代表对应被评估对象的当前实际损伤程度或支座广义位移，两个损伤向量的元素的元素个数等于被评估对象的数量，两个损伤向量的元素和被评估对象之间是一一对应关系，两个损伤向量的元素的编号规则和被评估对象的编号规则相同；k.依据“被监测量的当前数值向量Ci”同“被监测量的初始数值向量Cio”、“单位损伤被监测量数值变化矩阵ΔCi”和“当前名义健康状态向量dic”间存在的近似线性关系，该近似线性关系可表达为式1，式1中除dic外的其它量均为已知，求解式1就可以算出当前名义健康状态向量dic； C i = C o i + Δ C i · d c i 式1l.利用式2表达的当前实际健康状态向量di同初始损伤向量dio和当前名义健康状态向量dic的元素间的关系，计算得到当前实际健康状态向量di的所有元素； d j i = 1 - ( 1 - d oj i ) ( 1 - d ci i ) 式2式2中j＝1，2，3，……，N；当前实际健康状态向量di的元素数值代表对应被评估对象的实际损伤程度或实际支座广义位移，根据当前实际健康状态向量di就能确定有哪些索受损及其损伤程度，就能确定实际支座广义位移；若当前实际健康状态向量的某一元素对应于是索系统中的一根索，且其数值为0，表示该元素所对应的索是完好的，没有损伤或松弛的的，若其数值为100％，则表示该元素所对应的索已经完全丧失承载能力，若其数值介于0和100％之间，则表示该索丧失了相应比例的承载能力；如果当前实际健康状态向量的某一元素对应于一个支座的一个广义位移分量，那么dij表示其当前广义位移数值；m.从第l步中识别出的有问题的支承索中鉴别出受损索，剩下的就是松弛索；n.利用在第l步获得的当前实际虚拟损伤向量di得到松弛索的当前实际虚拟损伤程度，利用在第g步获得的当前索力向量Fi，利用在第g步获得的所有支承索的两个支承端点的水平距离，利用在第d步获得的初始自由长度向量Io，利用在第d步获得的所有索的弹性模量、密度、初始横截面面积数据，通过将松弛索同受损索进行力学等效来计算松弛索的、与当前实际虚拟损伤程度等效的松弛程度，等效的力学条件是：一、两等效的索的无松弛和无损伤时的初始自由长度、几何特性参数、密度及材料的力学特性参数相同；二、松弛或损伤后，两等效的松弛索和损伤索的索力和变形后的总长相同；满足上述两个等效条件时，这样的两根支承索在结构中的力学功能就是完全相同的，即如果用等效的松弛索代替受损索后，索结构不会发生任何变化，反之亦然；依据前述力学等效条件求得那些被判定为松弛索的松弛程度，松弛程度就是支承索自由长度的改变量，也就是确定了那些需调整索力的支承索的索长调整量；这样就实现了支承索的松弛识别；计算时所需索力由当前索力向量Fi对应元素给出；o.在求得当前名义健康状态向量dic后，按照式3建立标识向量Bi，式4给出了标识向量Bi的第j个元素的定义； B i = B 1 i B 2 i . . . B j i . . . B N i T 式3 B j i = 0 , if d cj i < D uj i 1 , if d cj i ≥ D uj i 式4式4中元素Bij是标识向量Bi的第j个元素，Diuj是名义单位损伤向量Diu的第j个元素，dicj是当前名义健康状态向量dic的第j个元素，它们都表示第j个被评估对象的相关信息，式4中j＝1，2，3，……，N；p.如果标识向量Bi的元素全为0，则回到第g步继续本次循环；如果标识向量Fi的元素不全为0，则进入下一步、即第q步；q.根据式5计算得到下一次、即第i+1次循环所需的初始损伤向量di+1o的每一个元素di+1oj； d oj i + 1 = 1 - ( 1 - d oj i ) ( 1 - D uj i F j i ) 式5式5中Diuj是名义单位损伤向量Diu的第j个元素，dicj是当前名义健康状态向量dic的第j个元素，Fij是标识向量Fi的第j个元素，式5中j＝1，2，3，……，N；向量di+1o的元素的编号规则和被评估对象的编号规则相同；r.在力学计算基准模型Ai的基础上，令被评估对象的健康状况为di+1o后更新得到下一次、即第i+1次循环所需的力学计算基准模型Ai+1；s.通过对力学计算基准模型Ai+1的计算得到对应于模型Ai+1的结构的所有被监测应变的点的、将被监测的应变方向的应变数值，这些数值组成下一次、即第i+1次循环所需的被监测量的初始数值向量Ci+1o；t.回到第f步，开始下一次循环。