发明名称 基于空间坐标监测的识别受损索和支座位移的递进式方法
摘要 基于空间坐标监测的识别受损索和支座位移的递进式方法基于空间坐标监测,考虑到了被监测量的当前数值向量同被监测量的初始数值向量、单位损伤被监测量变化矩阵和当前名义损伤向量间的线性关系是近似的,为克服此缺陷,本发明给出了使用线性关系分段逼近非线性关系的方法,将大区间分割成连续的一个个小区间,在每一个小区间内上述线性关系都是足够准确的,在每一个小区间内可以利用多目标优化算法等合适的算法快速识别出支座位移和受损索。
申请公布号 CN101832872A 申请公布日期 2010.09.15
申请号 CN201010137562.2 申请日期 2010.03.31
申请人 东南大学 发明人 韩玉林;王芳
分类号 G01M19/00(2006.01)I;G01B21/32(2006.01)I;G01B21/02(2006.01)I;G01B21/00(2006.01)I;G01N19/08(2006.01)I;G01L5/00(2006.01)I 主分类号 G01M19/00(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 一种基于空间坐标监测的识别受损索和支座位移的递进式方法,其特征在于所述方法包括:a.称被评估的支承索和支座位移分量为被评估对象,设被评估的支承索的数量和支座位移分量的数量之和为N,即被评估对象的数量为N;确定被评估对象的编号规则,按此规则将索结构中所有的被评估对象编号,该编号在后续步骤中将用于生成向量和矩阵;本发明用变量j表示这一编号,j=1,2,3,...,N;b.确定指定的将被监测空间坐标的被测量点,给所有指定点编号;确定过每一测量点的将被监测的空间坐标分量,给所有被测量空间坐标分量编号;上述编号在后续步骤中将用于生成向量和矩阵;“结构的全部被监测的空间坐标数据”由上述所有被测量空间坐标分量组成;为方便起见,在本发明中将“结构的被监测的空间坐标数据”称为“被监测量”;所有被测量空间坐标分量的数量之和不得小于N;c.利用被评估对象的无损检测数据等能够表达被评估对象的健康状态的数据建立被评估对象初始健康状态向量dio;如果没有被评估对象的无损检测数据时,向量dio的各元素数值取0;向量dio的元素的编号规则和被评估对象的编号规则相同;本发明用i表示循环次数,i=1,2,3,......;这里是第一次循环,i取1,即这里建立的初始健康状态向量dio可以具体化为d1o;d.在建立初始健康状态向量d1o的同时,直接测量计算得到索结构的所有被监测量的初始数值,组成被监测量的初始数值向量Cio;这里是第一次循环,i取1,即这里建立的被监测量的初始数值向量Cio可以具体化为C1o;在实测得到被监测量初始数值向量C1o的同时,实测得到索结构的所有索的初始索力数据、结构的初始几何数据和初始索结构支座坐标数据;e.根据索结构的设计图、竣工图和索结构的实测数据、索的无损检测数据和初始索结构支座坐标数据建立索结构的力学计算基准模型Ai;这里是第一次循环,i取1,即这里建立的索结构的力学计算基准模型Ai可以具体化为A1;f.在力学计算基准模型Ai的基础上进行若干次力学计算,通过计算获得“单位损伤被监测量数值变化矩阵ΔCi”和“名义单位损伤向量Diu”;g.实测得到索结构的所有指定被监测量的当前实测数值,组成“被监测量的当前数值向量Ci”;给本步及本步之前出现的所有向量的元素编号时,应使用同一编号规则,这样可以保证本步及本步之前出现的各向量的、编号相同的元素,表示同一被监测量的、对应于该元素所属向量所定义的相关信息;h.定义当前名义损伤向量dic和当前实际损伤向量di,两个损伤向量的元素个数等于被评估对象的数量,当前名义损伤向量dic的元素数值代表对应被评估对象的当前名义损伤程度或支座位移,当前实际损伤向量di的元素数值代表对应被评估对象的当前实际损伤程度或支座位移,两个损伤向量的元素的元素个数等于被评估对象的数量,两个损伤向量的元素和被评估对象之间是一一对应关系,两个损伤向量的元素的编号规则和被评估对象的编号规则相同;i.依据“被监测量的当前数值向量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>式1j.利用式2表达的当前实际损伤向量di同初始损伤向量dio和当前名义损伤向量dic的元素间的关系,计算得到当前实际损伤向量di的所有元素; <mrow> <msub> <msup> <mi>d</mi> <mi>i</mi> </msup> <mi>j</mi> </msub> <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,表示该元素所对应的索是完好的,没有损伤的,若其数值为100%,则表示该元素所对应的索已经完全丧失承载能力,若其数值介于0和100%之间,则表示该索丧失了相应比例的承载能力;如果当前实际损伤向量的某一元素对应于一个支座的一个位移分量,那么dij表示其当前位移数值;k.在求得当前名义损伤向量dic后,按照式3建立标识向量Fi,式4给出了标识向量Fi的第j个元素的定义; <mrow> <msup> <mi>F</mi> <mi>i</mi> </msup> <mo>=</mo> <msup> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msubsup> <mi>F</mi> <mn>1</mn> <mi>i</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>F</mi> <mn>2</mn> <mi>i</mi> </msubsup> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <msubsup> <mi>F</mi> <mi>j</mi> <mi>i</mi> </msubsup> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <msubsup> <mi>F</mi> <mi>N</mi> <mi>i</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> </mrow>式3 <mrow> <msubsup> <mi>F</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式4中元素Fij是标识向量Fi的第j个元素,Diuj是名义单位损伤向量Diu的第j个元素,dicj是当前名义损伤向量dic的第j个元素,它们都表示第j个被评估对象的相关信息,式4中j=1,2,3,……,N;l.如果标识向量Fi的元素全为0,则回到第g步继续本次循环;如果标识向量Fi的元素不全为0,则进入下一步、即第m步;m.根据式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>F</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> </mrow>式5式5中Diuj是名义单位损伤向量Diu的第j个元素,dicj是当前名义损伤向量dic的第j个元素,Fij是标识向量Fi的第j个元素,式5中j=1,2,3,……,N;向量di+1o的元素的编号规则和被评估对象的编号规则相同;n.在力学计算基准模型Ai的基础上,令被评估对象的健康状况为di+1o后更新得到下一次、即第i+1次循环所需的力学计算基准模型Ai+1;o.通过对力学计算基准模型Ai+1的计算得到对应于模型Ai+1的结构的所有被监测应变的点的、将被监测的应变方向的应变数值,这些数值组成下一次、即第i+1次循环所需的被监测量的初始数值向量Ci+1o;p.回到第f步,开始下一次循环。
地址 210009 江苏省南京市江宁开发区东南大学路2号