发明名称 精简广义位移空间坐标监测载荷受损索识别方法
摘要 精简广义位移空间坐标监测载荷受损索识别方法基于空间坐标监测,通过监测支座广义位移来决定是否需要更新索结构的力学计算基准模型,得到计入支座广义位移的索结构的力学计算基准模型,在此模型的基础上计算获得单位损伤被监测量数值变化矩阵。依据被监测量当前数值向量同被监测量当前初始数值向量、单位损伤被监测量数值变化矩阵和待求的被评估对象当前名义损伤向量间存在的近似线性关系算出被评估对象当前名义损伤向量的非劣解,据此可以识别出核心被评估对象的健康状态。
申请公布号 CN105115752A 申请公布日期 2015.12.02
申请号 CN201510438060.6 申请日期 2015.07.23
申请人 东南大学 发明人 韩玉林;韩佳邑
分类号 G01M99/00(2011.01)I;G01K13/00(2006.01)I;G01L5/00(2006.01)I;G01B21/00(2006.01)I;G01B21/02(2006.01)I 主分类号 G01M99/00(2011.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 杨晓玲
主权项 精简广义位移空间坐标监测载荷受损索识别方法,其特征在于所述方法包括:a.当索结构承受的载荷虽有变化,但索结构正在承受的载荷没有超出索结构初始许用载荷时,本方法适用;索结构初始许用载荷指索结构在竣工时的许用载荷,能够通过常规力学计算获得;本方法统一称被评估的支承索和载荷为“被评估对象”,设被评估的支承索的数量和载荷的数量之和为N,即“被评估对象”的数量为N;本方法用名称“核心被评估对象”专指“被评估对象”中的被评估的支承索,本方法用名称“次要被评估对象”专指“被评估对象”中的被评估的载荷;确定被评估对象的编号规则,按此规则将索结构中所有的被评估对象编号,该编号在后续步骤中将用于生成向量和矩阵;本方法用变量k表示这一编号,k=1,2,3,…,N;设索系统中共有M<sub>1</sub>根支承索,显然核心被评估对象的数量就是M<sub>1</sub>;确定指定的将被监测空间坐标的被测量点,给所有指定点编号;确定过每一测量点的将被监测的空间坐标分量,给所有被测量空间坐标分量编号;上述编号在后续步骤中将用于生成向量和矩阵;“索结构的全部被监测的空间坐标数据”由上述所有被测量空间坐标分量组成;为方便起见,在本方法中将“索结构的被监测的空间坐标数据”简称为“被监测量”;所有被监测量的数量之和记为M,M必须大于核心被评估对象的数量,M小于被评估对象的数量;物体、结构承受的外力可称为载荷,载荷包括面载荷和体积载荷;面载荷又称表面载荷,是作用于物体表面的载荷,包括集中载荷和分布载荷两种;体积载荷是连续分布于物体内部各点的载荷,包括物体的自重和惯性力在内;集中载荷分为集中力和集中力偶两种,在包括笛卡尔直角坐标系在内的坐标系中,一个集中力可以分解成三个分量,同样的,一个集中力偶也可以分解成三个分量,如果载荷实际上是集中载荷,在本方法中将一个集中力分量或一个集中力偶分量计为或统计为一个载荷,此时载荷的变化具体化为一个集中力分量或一个集中力偶分量的变化;分布载荷分为线分布载荷和面分布载荷,分布载荷的描述至少包括分布载荷的作用区域和分布载荷的大小,分布载荷的大小用分布集度来表达,分布集度用分布特征和幅值来表达;如果载荷实际上是分布载荷,本方法谈论载荷的变化时,实际上是指分布载荷分布集度的幅值的改变,而所有分布载荷的作用区域和分布集度的分布特征是不变的;在包括笛卡尔直角坐标系在内的坐标系中,一个分布载荷可以分解成三个分量,如果这分布载荷的三个分量的各自的分布集度的幅值发生变化,且变化的比率不全部相同,那么在本方法中把这分布载荷的三个分量计为或统计为三个分布载荷,此时一个载荷就代表分布载荷的一个分量;体积载荷是连续分布于物体内部各点的载荷,体积载荷的描述至少包括体积载荷的作用区域和体积载荷的大小,体积载荷的大小用分布集度来表达,分布集度用分布特征和幅值来表达;如果载荷实际上是体积载荷,在本方法中实际处理的是体积载荷分布集度的幅值的改变,而所有体积载荷的作用区域和分布集度的分布特征是不变的,此时在本方法中提到载荷的改变时实际上是指体积载荷的分布集度的幅值的改变,此时,发生变化的载荷是指那些分布集度的幅值发生变化的体积载荷;在包括笛卡尔直角坐标系在内的坐标系中,一个体积载荷可以分解成三个分量,如果这体积载荷的三个分量的各自的分布集度的幅值发生变化,且变化的比率不全部相同,那么在本方法中把这体积载荷的三个分量计为或统计为三个分布载荷;b.实测或查资料得到索结构所使用的各种材料的物理和力学性能参数;c.在实测或查资料得到索结构所使用的各种材料的物理和力学性能参数的同时,直接测量计算得到初始索结构的实测数据,初始索结构的实测数据是包括索结构集中载荷测量数据、索结构分布载荷测量数据、索结构体积载荷测量数据、所有被监测量的初始数值、所有支承索的初始索力数据、初始索结构模态数据、初始索结构应变数据、初始索结构几何数据、初始索结构支座广义坐标数据、初始索结构角度数据、初始索结构空间坐标数据在内的实测数据,初始索结构支座广义坐标数据包括初始索结构支座空间坐标数据和初始索结构支座角坐标数据,在得到初始索结构的实测数据的同时,测量计算得到包括支承索的无损检测数据在内的能够表达支承索的健康状态的数据,此时的能够表达支承索的健康状态的数据称为支承索初始健康状态数据;所有被监测量的初始数值组成被监测量初始数值向量C<sub>o</sub>,被监测量初始数值向量C<sub>o</sub>的编号规则与M个被监测量的编号规则相同;利用支承索初始健康状态数据以及索结构载荷测量数据建立被评估对象初始损伤向量d<sub>o</sub>,向量d<sub>o</sub>表示用初始力学计算基准模型A<sub>o</sub>表示的索结构的被评估对象的初始健康状态;被评估对象初始损伤向量d<sub>o</sub>的元素个数等于N,d<sub>o</sub>的元素与被评估对象是一一对应关系,向量d<sub>o</sub>的元素的编号规则与被评估对象的编号规则相同;如果d<sub>o</sub>的某一个元素对应的被评估对象是索系统中的一根支承索,那么d<sub>o</sub>的该元素的数值代表对应支承索的初始损伤程度,若该元素的数值为0,表示该元素所对应的支承索是完好的,没有损伤的,若其数值为100%,则表示该元素所对应的支承索已经完全丧失承载能力,若其数值介于0和100%之间,则表示该支承索丧失了相应比例的承载能力;如果d<sub>o</sub>的某一个元素对应的被评估对象是某一个载荷,本方法中取d<sub>o</sub>的该元素数值为0,代表这个载荷的变化的初始数值为0;如果没有支承索的无损检测数据及其他能够表达支承索的健康状态的数据时,或者可以认为结构初始状态为无损伤无松弛状态时,向量d<sub>o</sub>中与支承索相关的各元素数值取0;初始索结构支座广义坐标数据组成初始索结构支座广义坐标向量U<sub>o</sub>;d.根据索结构的设计图、竣工图和初始索结构的实测数据、支承索初始健康状态数据、索结构集中载荷测量数据、索结构分布载荷测量数据、索结构体积载荷测量数据、索结构所使用的各种材料的物理和力学性能参数、初始索结构支座广义坐标向量U<sub>o</sub>和前面步骤得到的所有的索结构数据,建立索结构的初始力学计算基准模型A<sub>o</sub>,基于A<sub>o</sub>计算得到的索结构计算数据必须非常接近其实测数据,其间的差异不得大于5%;对应于A<sub>o</sub>的被评估对象健康状态用被评估对象初始损伤向量d<sub>o</sub>表示;对应于A<sub>o</sub>的所有被监测量的初始数值用被监测量初始数值向量C<sub>o</sub>表示;第一次建立索结构的当前初始力学计算基准模型A<sup>t</sup><sub>o</sub>和被监测量当前初始数值向量C<sup>t</sup><sub>o</sub>;第一次建立索结构的当前初始力学计算基准模型A<sup>t</sup><sub>o</sub>和被监测量当前初始数值向量C<sup>t</sup><sub>o</sub>时,索结构的当前初始力学计算基准模型A<sup>t</sup><sub>o</sub>就等于索结构的初始力学计算基准模型A<sub>o</sub>,被监测量当前初始数值向量C<sup>t</sup><sub>o</sub>就等于被监测量初始数值向量C<sub>o</sub>;对应于索结构的当前初始力学计算基准模型A<sup>t</sup><sub>o</sub>的索结构支座广义坐标数据组成当前初始索结构支座广义坐标向量U<sup>t</sup><sub>o</sub>,第一次建立索结构的当前初始力学计算基准模型A<sup>t</sup><sub>o</sub>时,U<sup>t</sup><sub>o</sub>就等于U<sub>o</sub>;A<sup>t</sup><sub>o</sub>的被评估对象的初始健康状态与A<sub>o</sub>的被评估对象的健康状态相同,也用被评估对象初始损伤向量d<sub>o</sub>表示,在后面的循环过程中A<sup>t</sup><sub>o</sub>的被评估对象的初始健康状态始终用被评估对象初始损伤向量d<sub>o</sub>表示;U<sub>o</sub>和d<sub>o</sub>是A<sub>o</sub>的参数,由A<sub>o</sub>的力学计算结果得到的所有被监测量的初始数值与C<sub>o</sub>表示的所有被监测量的初始数值相同,因此也可以说C<sub>o</sub>由A<sub>o</sub>的力学计算结果组成;U<sup>t</sup><sub>o</sub>和d<sub>o</sub>是A<sup>t</sup><sub>o</sub>的参数,C<sup>t</sup><sub>o</sub>由A<sup>t</sup><sub>o</sub>的力学计算结果组成;e.从这里进入由第e步到第m步的循环;在结构服役过程中,不断实测得到索结构支座广义坐标当前数据,所有索结构支座广义坐标当前数据组成当前索结构实测支座广义坐标向量U<sup>t</sup>,向量U<sup>t</sup>的定义方式与向量U<sub>o</sub>的定义方式相同;f.根据当前索结构实测支座广义坐标向量U<sup>t</sup>,按照步骤f1至f3更新当前初始力学计算基准模型A<sup>t</sup><sub>o</sub>、当前初始索结构支座广义坐标向量U<sup>t</sup><sub>o</sub>和被监测量当前初始数值向量C<sup>t</sup><sub>o</sub>;f1.比较U<sup>t</sup>与U<sup>t</sup><sub>o</sub>,如果U<sup>t</sup>等于U<sup>t</sup><sub>o</sub>,则A<sup>t</sup><sub>o</sub>、U<sup>t</sup><sub>o</sub>和C<sup>t</sup><sub>o</sub>保持不变,否则需要按下列步骤对A<sup>t</sup><sub>o</sub>、U<sup>t</sup><sub>o</sub>和C<sup>t</sup><sub>o</sub>进行更新;f2.计算U<sup>t</sup>与U<sub>o</sub>的差,U<sup>t</sup>与U<sub>o</sub>的差就是索结构支座关于初始位置的支座广义位移,用支座广义位移向量V表示支座广义位移,V等于U<sup>t</sup>减去U<sub>o</sub>,支座广义位移向量V中的元素与支座广义位移分量之间是一一对应关系,支座广义位移向量V中一个元素的数值对应于一个指定支座的一个指定方向的广义位移;f3.对A<sub>o</sub>中的索结构支座施加支座广义位移约束,支座广义位移约束的数值就取自支座广义位移向量V中对应元素的数值,对A<sub>o</sub>中索结构支座施加支座广义位移约束后得到更新的当前初始力学计算基准模型A<sup>t</sup><sub>o</sub>,更新A<sup>t</sup><sub>o</sub>的同时,U<sup>t</sup><sub>o</sub>所有元素数值也用U<sup>t</sup>所有元素数值对应代替,即更新了U<sup>t</sup><sub>o</sub>,这样就得到了正确地对应于A<sup>t</sup><sub>o</sub>的U<sup>t</sup><sub>o</sub>;更新C<sup>t</sup><sub>o</sub>的方法是:当更新A<sup>t</sup><sub>o</sub>后,通过力学计算得到A<sup>t</sup><sub>o</sub>中所有被监测量的、当前的具体数值,这些具体数值组成C<sup>t</sup><sub>o</sub>;A<sup>t</sup><sub>o</sub>的支承索的初始健康状态始终用被评估对象初始损伤向量d<sub>o</sub>表示;g.在当前初始力学计算基准模型A<sup>t</sup><sub>o</sub>的基础上按照步骤g1至g4进行若干次力学计算,通过计算获得索结构单位损伤被监测量数值变化矩阵ΔC和被评估对象单位变化向量D<sub>u</sub>;g1.索结构单位损伤被监测量数值变化矩阵ΔC是不断更新的,即在更新当前初始力学计算基准模型A<sup>t</sup><sub>o</sub>、当前初始索结构支座广义坐标向量U<sup>t</sup><sub>o</sub>和被监测量当前初始数值向量C<sup>t</sup><sub>o</sub>之后,必须接着更新索结构单位损伤被监测量数值变化矩阵ΔC和被评估对象单位变化向量D<sub>u</sub>;g2.在索结构的当前初始力学计算基准模型A<sup>t</sup><sub>o</sub>的基础上进行若干次力学计算,计算次数数值上等于所有被评估对象的数量N,有N个评估对象就有N次计算;依据被评估对象的编号规则,依次进行计算;每一次计算假设只有一个被评估对象在原有损伤或载荷的基础上再增加单位损伤或载荷单位变化,具体的,如果该被评估对象是索系统中的一根支承索,那么就假设该支承索在向量d<sub>o</sub>表示的该支承索已有损伤的基础上再增加单位损伤,如果该被评估对象是一个载荷,就假设该载荷在向量d<sub>o</sub>表示的该载荷已有变化量的基础上再增加载荷单位变化,用D<sub>uk</sub>记录这一增加的单位损伤或载荷单位变化,其中k表示增加单位损伤或载荷单位变化的被评估对象的编号,D<sub>uk</sub>是被评估对象单位变化向量D<sub>u</sub>的一个元素,被评估对象单位变化向量D<sub>u</sub>的元素的编号规则与向量d<sub>o</sub>的元素的编号规则相同;每一次计算中增加单位损伤或载荷单位变化的被评估对象不同于其它次计算中增加单位损伤或载荷单位变化的被评估对象,每一次计算都利用力学方法计算索结构的所有被监测量的当前计算值,每一次计算得到的所有被监测量的当前计算值组成一个被监测量计算当前向量,被监测量计算当前向量的元素编号规则与被监测量初始数值向量C<sub>o</sub>的元素编号规则相同;g3.每一次计算得到的被监测量计算当前向量减去被监测量当前初始数值向量C<sup>t</sup><sub>o</sub>得到一个向量,再将该向量的每一个元素都除以该次计算所假设的单位损伤或载荷单位变化数值,得到一个被监测量单位变化向量,有N个被评估对象就有N个被监测量单位变化向量;g4.由这N个被监测量单位变化向量按照N个被评估对象的编号规则,依次组成有N列的索结构单位损伤被监测量数值变化矩阵ΔC;索结构单位损伤被监测量数值变化矩阵ΔC的每一列对应于一个被监测量单位变化向量;索结构单位损伤被监测量数值变化矩阵ΔC的每一行对应于同一个被监测量在不同被评估对象增加单位损伤或载荷单位变化时的不同的单位变化幅度;索结构单位损伤被监测量数值变化矩阵ΔC的列的编号规则与向量d<sub>o</sub>的元素的编号规则相同,索结构单位损伤被监测量数值变化矩阵ΔC的行的编号规则与M个被监测量的编号规则相同;h.在实测得到当前索结构实测支座广义坐标向量U<sup>t</sup>的同时,实测得到索结构的所有被监测量的当前实测数值,组成被监测量当前数值向量C;被监测量当前数值向量C和被监测量当前初始数值向量C<sup>t</sup><sub>o</sub>与被监测量初始数值向量C<sub>o</sub>的定义方式相同,三个向量的相同编号的元素表示同一被监测量在不同时刻的具体数值;i.定义被评估对象当前名义损伤向量d,被评估对象当前名义损伤向量d的元素个数等于被评估对象的数量,被评估对象当前名义损伤向量d的元素和被评估对象之间是一一对应关系,被评估对象当前名义损伤向量d的元素数值代表对应被评估对象的名义损伤程度或名义载荷变化量;向量d的元素的编号规则与向量d<sub>o</sub>的元素的编号规则相同;j.依据被监测量当前数值向量C同被监测量当前初始数值向量C<sup>t</sup><sub>o</sub>、索结构单位损伤被监测量数值变化矩阵ΔC和待求的被评估对象当前名义损伤向量d间存在的近似线性关系,该近似线性关系可表达为式1,式1中除d外的其它量均为已知,求解式1就可以算出被评估对象当前名义损伤向量d;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>C</mi><mo>=</mo><msubsup><mi>C</mi><mi>o</mi><mi>t</mi></msubsup><mo>+</mo><mi>&Delta;</mi><mi>C</mi><mo>&CenterDot;</mo><mi>d</mi></mrow>]]></math><img file="FDA0000765553310000051.GIF" wi="332" he="73" /></maths>   式1k.定义被评估对象当前实际损伤向量d<sup>a</sup>,被评估对象当前实际损伤向量d<sup>a</sup>的元素个数等于被评估对象的数量,被评估对象当前实际损伤向量d<sup>a</sup>的元素和被评估对象之间是一一对应关系,被评估对象当前实际损伤向量d<sup>a</sup>的元素数值代表对应被评估对象的实际损伤程度或实际载荷变化量;向量d<sup>a</sup>的元素的编号规则与向量d<sub>o</sub>的元素的编号规则相同;l.利用式2表达的被评估对象当前实际损伤向量d<sup>a</sup>的第k个元素d<sup>a</sup><sub>k</sub>同被评估对象初始损伤向量d<sub>o</sub>的第k个元素d<sub>ok</sub>和被评估对象当前名义损伤向量d的第k个元素d<sub>k</sub>间的关系,计算得到被评估对象当前实际损伤向量d<sup>a</sup>的所有元素;<img file="FDA0000765553310000061.GIF" wi="1318" he="172" />式2式2中k=1,2,3,…….,N,d<sup>a</sup><sub>k</sub>表示第k个被评估对象的当前实际健康状态,如果该被评估对象是索系统中的一根支承索,那么d<sup>a</sup><sub>k</sub>表示其当前实际损伤,d<sup>a</sup><sub>k</sub>为0时表示无损伤,为100%时表示该支承索彻底丧失承载能力,介于0与100%之间时表示丧失相应比例的承载能力;所以根据被评估对象当前实际损伤向量d<sup>a</sup>能够确定核心被评估对象的健康状态;m.回到第e步,开始由第e步到第m步的下一次循环。
地址 211189 江苏省南京市江宁区东南大学路2号