CN103852282A - 索力监测问题索载荷递进式识别方法 - Google Patents

索力监测问题索载荷递进式识别方法 Download PDF

Info

Publication number
CN103852282A
CN103852282A CN201410085128.2A CN201410085128A CN103852282A CN 103852282 A CN103852282 A CN 103852282A CN 201410085128 A CN201410085128 A CN 201410085128A CN 103852282 A CN103852282 A CN 103852282A
Authority
CN
China
Prior art keywords
cable
temperature
data
vector
load
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201410085128.2A
Other languages
English (en)
Inventor
韩玉林
韩佳邑
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201410085128.2A priority Critical patent/CN103852282A/zh
Publication of CN103852282A publication Critical patent/CN103852282A/zh
Pending legal-status Critical Current

Links

Abstract

索力监测问题索载荷递进式识别方法基于索力监测、通过监测索结构温度、环境温度、载荷变化程度和支承索健康程度来决定是否需要更新索结构的力学计算基准模型,得到新的计入载荷变化程度、支承索健康程度、温度的索结构的力学计算基准模型,在此模型的基础上依据被监测量的当前数值向量同被监测量当前初始数值向量、单位损伤被监测量数值变化矩阵和待求的当前名义损伤向量间存在的近似线性关系,据此可以在有温度变化时,能够剔除干扰因素的影响,准确地识别问题索和载荷变化量。

Description

索力监测问题索载荷递进式识别方法
技术领域
[0001] 斜拉桥、悬索桥、桁架结构等结构有一个共同点,就是它们有许多承受拉伸载荷的部件,如斜拉索、主缆、吊索、拉杆等等,该类结构的共同点是以索、缆或仅承受拉伸载荷的杆件为支承部件,为方便起见,本方法将该类结构表述为“索结构”,并将索结构的所有承载索、承载缆,及所有仅承受轴向拉伸或轴向压缩载荷的杆件(又称为二力杆件),为方便起见统一称为“索系统”,本方法中用“支承索”这一名词指称承载索、承载缆及仅承受轴向拉伸或轴向压缩载荷的杆件,有时简称为“索”,所以在后面使用“索”这个字的时候,对桁架结构实际就是指二力杆件。支承索的受损和松弛对索结构安全是一项重大威胁,本方法将受损索和松弛索统称为有健康问题的支承索,简称为问题索。在结构服役过程中,对支承索或索系统的健康状态的正确识别关系到整个索结构的安全。在环境温度发生变化时,索结构的温度一般也会随着发生变化,在索结构温度发生变化时,索结构承受的载荷也可能发生变化,实际上即使索结构的温度不发生变化,索结构承受的载荷也可能单独发生变化,同时索结构的健康状态也可能在发生变化,在这种复杂条件下,本方法基于索力监测(本方法将被监测的索力称为“被监测量”)来识别问题索和索结构承受的载荷的变化量,属工程结构健康监测领域。
背景技术
[0002] 剔除载荷变化和结构温度变化对索结构健康状态识别结果的影响,从而准确地识别结构的健康状态的变化,是目前迫切需要解决的问题;同样的,剔除结构温度变化和结构健康状态变化对结构承受的载荷的变化量的识别结果的影响,对结构安全同样具有重要意义;目前还没有一种有效的方法能够同时解决这两个问题,本方法公开了解决这两个问题的一种有效方法。也就是说,本方法实现了已有方法不可能具备的两种功能。
发明内容
[0003] 技术问题:本方法公开了一种方法,实现了已有方法不可能具备的两种功能,分别是,一、在结构承受的载荷和结构(环境)温度变化时,能够剔除载荷变化和结构温度变化对索结构健康状态识别结果的影响,从而准确地识别出问题索;二、本方法在识别出问题索的同时,还能同时识别出载荷的变化,即本方法能够剔除结构温度变化和支承索健康状态变化的影响,实现载荷变化程度的正确识别。
[0004] 在索结构服役过程中,支承索自由状态(此时索张力也称索力为O)下的索长度(称为自由长度,本方法专指支承索两支承端点间的那段索的自由长度)会发生变化,本方法的目的之一就是要识别出自由长度发生了变化的支承索,并识别出它们的自由长度的改变量,此改变量为该索的索力调整提供了直接依据,为了方便,本方法将自由长度发生变化的支承索统称为松弛索。
[0005] 技术方案:物体、结构承受的外力可称为载荷,载荷包括面载荷和体积载荷。面载荷又称表面载荷,是作用于物体表面的载荷,包括集中载荷和分布载荷两种。体积载荷是连续分布于物体内部各点的载荷,如物体的自重和惯性力。
[0006] 集中载荷分为集中力和集中力偶两种,在坐标系中,例如在笛卡尔直角坐标系中,一个集中力可以分解成三个分量,同样的,一个集中力偶也可以分解成三个分量,如果载荷实际上是集中载荷,在本方法中将一个集中力分量或一个集中力偶分量称为一个载荷,此时载荷的变化具体化为一个集中力分量或一个集中力偶分量的变化。
[0007] 分布载荷分为线分布载荷和面分布载荷,分布载荷的描述至少包括分布载荷的作用区域和分布载荷的大小,分布载荷的大小用分布集度来表达,分布集度用分布特征(例如均布、正弦函数等分布特征)和幅值来表达(例如两个分布载荷都是均布,但其幅值不同,可以均布压力为例来说明幅值的概念:同一个结构承受两个不同的均布压力,两个分布载荷都是均布载荷,但一个分布载荷的幅值是IOMPa,另一个分布载荷的幅值是50MPa)。如果载荷实际上是分布载荷,本方法谈论载荷的变化时,实际上是指分布载荷分布集度的幅值的改变,而分布载荷的作用区域和分布集度的分布特征是不变的。在坐标系中,一个分布载荷可以分解成若干个分量,如果这分布载荷的若干个分量的各自的分布集度的幅值发生变化,且变化的比率不全部相同,那么在本方法中把这若干个分布载荷的分量看成同样数量的独立的分布载荷,此时一个载荷就代表一个分布载荷的分量,也可以将其中分布集度的幅值变化比率相同的分量合成为一个分布载荷或称为一个载荷。
[0008] 体积载荷是连续分布于物体内部各点的载荷,如物体的自重和惯性力,体积载荷的描述至少包括体积载荷的作用区域和体积载荷的大小,体积载荷的大小用分布集度来表达,分布集度用分布特征(例如均布、线性函数等分布特征)和幅值来表达(例如两个体积载荷都是均布,但其幅值不同,可以自重为例来说明幅值的概念:同一个结构的两个部分的材料不同,故密度不同,所以虽然这两个部分所受的体积载荷都是均布的,但一个部分所受的体积载荷的幅值可能是10kN/m3,另一个部分所受的体积载荷的幅值是50kN/m3)。如果载荷实际上是体积载荷,在本方法中实际处理的是体积载荷分布集度的幅值的改变,而体积载荷的作用区域和分布集度的分布特征是不变的,此时在本方法中提到载荷的改变时实际上是指体积载荷的分布集度的幅值的改变,此时,发生变化的载荷是指那些分布集度的幅值发生变化的体积载荷。在坐标系中,一个体积载荷可以分解成若干个分量(例如在笛卡尔直角坐标系中,体积载荷可以分解成关于坐标系的三个轴的分量,也就是说,在笛卡尔直角坐标系中体积载荷可以分解成三个分量),如果这体积载荷的若干个分量的各自的分布集度的幅值发生变化,且变化的比率不全部相同,那么在本方法中把这若干个体积载荷的分量看成同样数量的独立的载荷,也可以将其中分布集度的幅值变化比率相同的体积载荷分量合成为一个体积载荷或称为一个载荷。
[0009]当载荷具体化为集中载荷时,在本方法中,“载荷单位变化”实际上是指“集中载荷的单位变化”,类似的,“载荷变化”具体指“集中载荷的大小的变化”,“载荷变化量”具体指“集中载荷的大小的变化量”,“载荷变化程度”具体指“集中载荷的大小的变化程度”,“载荷的实际变化量”是指“集中载荷的大小的实际变化量”,“发生变化的载荷”是指“大小发生变化的集中载荷”,简单地说,此时“某某载荷的某某变化”是指“某某集中载荷的大小的某某变化”。
[0010] 当载荷具体化为分布载荷时,在本方法中,“载荷单位变化”实际上是指“分布载荷的分布集度的幅值的单位变化”,而分布载荷的分布特征是不变的,类似的,“载荷变化”具体指“分布载荷的分布集度的幅值的变化”,而分布载荷的分布特征是不变的,“载荷变化量”具体指“分布载荷的分布集度的幅值的变化量”,“载荷变化程度”具体指“分布载荷的分布集度的幅值的变化程度”,“载荷的实际变化量”具体指“分布载荷的分布集度的幅值的实际变化量”,“发生变化的载荷”是指“分布集度的幅值发生变化的分布载荷”,简单地说,此时“某某载荷的某某变化”是指“某某分布载荷的分布集度的幅值的某某变化”,而所有分布载荷的作用区域和分布集度的分布特征是不变的。
[0011] 当载荷具体化为体积载荷时,在本方法中,“载荷单位变化”实际上是指“体积载荷的分布集度的幅值的单位变化”,类似的,“载荷变化”是指“体积载荷的分布集度的幅值的变化”,“载荷变化量”是指“体积载荷的分布集度的幅值的变化量”,“载荷变化程度”是指“体积载荷的分布集度的幅值的变化程度”,“载荷的实际变化量”是指“体积载荷的分布集度的幅值的实际变化量”,“发生变化的载荷”是指“分布集度的幅值发生变化的体积载荷”,简单地说,“某某载荷的某某变化”是指“某某体积载荷的分布集度的幅值的某某变化”,而所有体积载荷的作用区域和分布集度的分布特征是不变的。
[0012] 本方法具体包括:
[0013] a.为叙述方便起见,本方法统一称被评估的支承索和载荷为被评估对象,设被评估的支承索的数量和载荷的数量之和为N,即被评估对象的数量为N ;确定被评估对象的编号规则,按此规则将索结构中所有的被评估对象编号,该编号在后续步骤中将用于生成向量和矩阵;本方法用变量k表示这一编号,k=l, 2,3,…,N ;设索系统中共有M1根支承索,索结构索力数据包括这M1根支承索的索力,显然M1小于被评估对象的数量N ;仅仅通过M1个支承索的M1个索力数据来求解未知的N个被评估对象的状态是不可能的,本方法在监测全部札根支承索索力的基础上,在索结构上人为增加M2根索,称为传感索,在索结构健康监测过程中将监测这新增加的M2根传感索的索力;综合上述被监测量,整个索结构共有M根索的M个索力被监测,即有M个被监测量,其中M为M1与M2之和;M应当大于被评估对象的数量N ;新增加的M2根传感索的刚度同索结构的任意一根支承索的刚度相比,应当小得多;新增加的M2根传感索的各传感索的索力应当比索结构的任意一根支承索的索力小得多,这样可以保证即使这新增加的M2根传感索出现了损伤或松弛,对索结构其他构件的应力、应变、变形的影响微乎其微;新增加的M2根传感索的横截面上正应力应当小于其疲劳极限,这些要求可以保证新增加的M2根传感索不会发生疲劳损伤;新增加的M2根传感索的两端应当充分锚固,保证不会出现松弛;新增加的M2根传感索应当得到充分的防腐蚀保护,保证新增加的M2根传感索不会发生损伤和松弛;为方便起见,在本方法中将“索结构的被监测的所有参量”简称为“被监测量”;给M个被监测量连续编号,本方法用用变量j表示这一编号,j=l,2,3,".,Μ,该编号在后续步骤中将用于生成向量和矩阵;在本方法中新增加的M2根传感索作为索结构的一部分,后文再提到索结构时,索结构包括增加M2根传感索前的索结构和新增加的M2根传感索,也就是说后文提到索结构时指包括新增加的M2根传感索的索结构;因此后文提到按照“本方法的索结构的温度测量计算方法”测量计算得到“索结构稳态温度数据”时,其中的索结构包括新增加的M2根传感索,得到的“索结构稳态温度数据”包括新增加的M2根传感索的稳态温度数据,获得新增加的M2根传感索的稳态温度数据的方法同于索结构的M1根支承索的稳态温度数据的获得方法,在后文不再一一交代;测量得到新增加的M2根传感索的索力的方法同于索结构的M1根支承索的索力的测量方法,在后文不再一一交代;对索结构的支承索进行任何测量时,同时对新增加的M2根传感索进行同样的测量,在后文不再一一交代;新增加的仏根传感索除了不发生损伤和松弛外,对新增加的仏根传感索的信息量的要求和获得方法与索结构的支承索的信息量的要求和获得方法相同,在
后文不再--交代;在后文建立索结构的各种力学模型时,将新增加的M2根传感索视同索
结构的支承索对待;在后文中,除了提到支承索的损伤和松弛的场合外,当提到支承索时所说的支承索包括索结构的支承索和新增加的M2根传感索;本方法中对同一个量实时监测的任何两次测量之间的时间间隔不得大于30分钟,测量记录数据的时刻称为实际记录数据时刻;物体、结构承受的外力可称为载荷,载荷包括面载荷和体积载荷;面载荷又称表面载荷,是作用于物体表面的载荷,包括集中载荷和分布载荷两种;体积载荷是连续分布于物体内部各点的载荷,包括物体的自重和惯性力在内;集中载荷分为集中力和集中力偶两种,在包括笛卡尔直角坐标系在内的坐标系中,一个集中力可以分解成三个分量,同样的,一个集中力偶也可以分解成三个分量,如果载荷实际上是集中载荷,在本方法中将一个集中力分量或一个集中力偶分量计为或统计为一个载荷,此时载荷的变化具体化为一个集中力分量或一个集中力偶分量的变化;分布载荷分为线分布载荷和面分布载荷,分布载荷的描述至少包括分布载荷的作用区域和分布载荷的大小,分布载荷的大小用分布集度来表达,分布集度用分布特征和幅值来表达;如果载荷实际上是分布载荷,本方法谈论载荷的变化时,实际上是指分布载荷分布集度的幅值的改变,而所有分布载荷的作用区域和分布集度的分布特征是不变的;在包括笛卡尔直角坐标系在内的坐标系中,一个分布载荷可以分解成三个分量,如果这分布载荷的三个分量的各自的分布集度的幅值发生变化,且变化的比率不全部相同,那么在本方法中把这分布载荷的三个分量计为或统计为三个分布载荷,此时一个载荷就代表分布载荷的一个分量;体积载荷是连续分布于物体内部各点的载荷,体积载荷的描述至少包括体积载荷的作用区域和体积载荷的大小,体积载荷的大小用分布集度来表达,分布集度用分布特征和幅值来表达;如果载荷实际上是体积载荷,在本方法中实际处理的是体积载荷分布集度的幅值的改变,而所有体积载荷的作用区域和分布集度的分布特征是不变的,此时在本方法中提到载荷的改变时实际上是指体积载荷的分布集度的幅值的改变,此时,发生变化的载荷是指那些分布集度的幅值发生变化的体积载荷;在包括笛卡尔直角坐标系在内的坐标系中,一个体积载荷可以分解成三个分量,如果这体积载荷的三个分量的各自的分布集度的幅值发生变化,且变化的比率不全部相同,那么在本方法中把这体积载荷的三个分量计为或统计为三个分布载荷;
[0014] b.本方法定义“ 本方法的索结构的温度测量计算方法”按步骤bl至b3进行;
[0015] bl:查询或实测得到索结构组成材料及索结构所处环境的随温度变化的传热学参数,利用索结构的设计图、竣工图和索结构的几何实测数据,利用这些数据和参数建立索结构的传热学计算模型;查询索结构所在地不少于2年的近年来的气象资料,统计得到这段时间内的阴天数量记为T个阴天,在本方法中将白天不能见到太阳的一整日称为阴天,统计得到T个阴天中每一个阴天的O时至次日日出时刻后30分钟之间的最高气温与最低气温,日出时刻是指根据地球自转和公转规律确定的气象学上的日出时刻,不表示当天一定可以看见太阳,可以查询资料或通过常规气象学计算得到所需的每一日的日出时刻,每一个阴天的O时至次日日出时刻后30分钟之间的最高气温减去最低气温称为该阴天的日气温的最大温差,有T个阴天,就有T个阴天的日气温的最大温差,取T个阴天的日气温的最大温差中的最大值为参考日温差,参考日温差记为ΛΊ;;查询索结构所在地和所在海拔区间不少于2年的近年来的气象资料或实测得到索结构所处环境的温度随时间和海拔高度的变化数据和变化规律,计算得到索结构所在地和所在海拔区间不少于2年的近年来的索结构所处环境的温度关于海拔高度的最大变化率ATh,为方便叙述取ATh的单位为。C/m;在索结构的表面上取“R个索结构表面点”,取“R个索结构表面点”的具体原则在步骤b3中叙述,后面将通过实测得到这R个索结构表面点的温度,称实测得到的温度数据为“R个索结构表面温度实测数据”,如果是利用索结构的传热学计算模型,通过传热计算得到这R个索结构表面点的温度,就称计算得到的温度数据为“R个索结构表面温度计算数据”;从索结构所处的最低海拔到最高海拔之间,在索结构上均布选取不少于三个不同的海拔高度,在每一个选取的海拔高度处、在水平面与索结构表面的交线处至少选取两个点,从选取点处引索结构表面的外法线,所有选取的外法线方向称为“测量索结构沿壁厚的温度分布的方向”,测量索结构沿壁厚的温度分布的方向与“水平面与索结构表面的交线”相交,在选取的测量索结构沿壁厚的温度分布的方向中必须包括索结构的向阳面外法线方向和索结构的背阴面外法线方向,沿每一个测量索结构沿壁厚的温度分布的方向在索结构中均布选取不少于三个点,测量所有被选取点的温度,测得的温度称为“索结构沿厚度的温度分布数据”,其中沿与同一“水平面与索结构表面的交线”相交的、“测量索结构沿壁厚的温度分布的方向”测量获得的“索结构沿厚度的温度分布数据”,在本方法中称为“相同海拔高度索结构沿厚度的温度分布数据”,设选取了 H个不同的海拔高度,在每一个海拔高度处,选取了 B个测量索结构沿壁厚的温度分布的方向,沿每个测量索结构沿壁厚的温度分布的方向在索结构中选取了 E个点,其中H和E都不小于3,B不小于2,设HBE为H与B和E的乘积,对应的共有HBE个“测量索结构沿厚度的温度分布数据的点”,后面将通过实测得到这HBE个“测量索结构沿厚度的温度分布数据的点”的温度,称实测得到的温度数据为“HBE个索结构沿厚度温度实测数据”,如果是利用索结构的传热学计算模型,通过传热计算得到这HBE个测量索结构沿厚度的温度分布数据的点的温度,就称计算得到的温度数据为“HBE个索结构沿厚度温度计算数据”;设BE为B和E的乘积,本方法中在每一个选取的海拔高度处共有BE个“相同海拔高度索结构沿厚度的温度分布数据”;在索结构所在地按照气象学测量气温要求选取一个位置,将在此位置实测得到符合气象学测量气温要求的索结构所在环境的气温;在索结构所在地的空旷无遮挡处选取一个位置,该位置应当在全年的每一日都能得到该地所能得到的该日的最充分的`日照,在该位置安放一块碳钢材质的平板,称为参考平板,参考平板与地面不可接触,参考平板离地面距离不小于1.5米,该参考平板的一面向阳,称为向阳面,参考平板的向阳面是粗糙的和深色的,参考平板的向阳面应当在全年的每一日都能得到一块平板在该地所能得到的该日的最充分的日照,参考平板的非向阳面覆有保温材料,将实时监测得到参考平板的向阳面的温度;
[0016] b2:实时监测得到上述R个索结构表面点的R个索结构表面温度实测数据,同时实时监测得到前面定义的索结构沿厚度的温度分布数据,同时实时监测得到符合气象学测量气温要求的索结构所在环境的气温数据;通过实时监测得到当日日出时刻到次日日出时刻后30分钟之间的索结构所在环境的气温实测数据序列,索结构所在环境的气温实测数据序列由当日日出时刻到次日日出时刻后30分钟之间的索结构所在环境的气温实测数据按照时间先后顺序排列,找到索结构所在环境的气温实测数据序列中的最高温度和最低温度,用索结构所在环境的气温实测数据序列中的最高温度减去最低温度得到索结构所在环境的当日日出时刻到次日日出时刻后30分钟之间的最大温差,称为环境最大温差,记为Λ Traiax ;由索结构所在环境的气温实测数据序列通过常规数学计算得到索结构所在环境的气温关于时间的变化率,该变化率也随着时间变化;通过实时监测得到当日日出时刻到次日日出时刻后30分钟之间的参考平板的向阳面的温度的实测数据序列,参考平板的向阳面的温度的实测数据序列由当日日出时刻到次日日出时刻后30分钟之间的参考平板的向阳面的温度的实测数据按照时间先后顺序排列,找到参考平板的向阳面的温度的实测数据序列中的最高温度和最低温度,用参考平板的向阳面的温度的实测数据序列中的最高温度减去最低温度得到参考平板的向阳面的温度的当日日出时刻到次日日出时刻后30分钟之间的最大温差,称为参考平板最大温差,记为ATpmax ;通过实时监测得到当日日出时刻到次日日出时刻后30分钟之间的所有R个索结构表面点的索结构表面温度实测数据序列,有R个索结构表面点就有R个索结构表面温度实测数据序列,每一个索结构表面温度实测数据序列由一个索结构表面点的当日日出时刻到次日日出时刻后30分钟之间的索结构表面温度实测数据按照时间先后顺序排列,找到每一个索结构表面温度实测数据序列中的最高温度和最低温度,用每一个索结构表面温度实测数据序列中的最高温度减去最低温度得到每一个索结构表面点的温度的当日日出时刻到次日日出时刻后30分钟之间的最大温差,有R个索结构表面点就有R个当日日出时刻到次日日出时刻后30分钟之间的最大温差数值,其中的最大值称为索结构表面最大温差,记为ATsmax ;由每一索结构表面温度实测数据序列通过常规数学计算得到每一个索结构表面点的温度关于时间的变化率,每一个索结构表面点的温度关于时间的变化率也随着时间变化;通过实时监测得到当日日出时刻到次日日出时刻后30分钟之间的、在同一时刻、HBE个“索结构沿厚度的温度分布数据”后,计算在每一个选取的海拔高度处共计BE个“相同海拔高度索结构沿厚度的温度分布数据”中的最高温度与最低温度的差值,这个差值的绝对值称为“相同海拔高度处索结构厚度方向最大温差”,选取了 H个不同的海拔高度就有H个“相同海拔高度处索结构厚度方向最大温差”,称这H个“相同海拔高度处索结构厚度方向最大温差”中的最大值为“索结构厚度方向最大温差”,记为 ATtniax ;
[0017] b3:测量计算获得索`结构稳态温度数据;首先,确定获得索结构稳态温度数据的时刻,与决定获得索结构稳态温度数据的时刻相关的条件有六项,第一项条件是获得索结构稳态温度数据的时刻介于当日日落时刻到次日日出时刻后30分钟之间,日落时刻是指根据地球自转和公转规律确定的气象学上的日落时刻,可以查询资料或通过常规气象学计算得到所需的每一日的日落时刻;第二项条件的a条件是在当日日出时刻到次日日出时刻后30分钟之间的这段时间内,参考平板最大温差Λ Tpmax和索结构表面最大温差八1;_都不大于5摄氏度;第二项条件的b条件是在当日日出时刻到次日日出时刻后30分钟之间的这段时间内,在前面测量计算得到的环境最大误差不大于参考日温差ΛΊ;,且参考平板最大温差Λ Tpmax减去2摄氏度后不大于Δ Temax,且索结构表面最大温差Δ Tsmax不大于ATpmax ;只需满足第二项的a条件和b条件中的一项就称为满足第二项条件;第三项条件是在获得索结构稳态温度数据的时刻,索结构所在环境的气温关于时间的变化率的绝对值不大于每小时0.1摄氏度;第四项条件是在获得索结构稳态温度数据的时刻,R个索结构表面点中的每一个索结构表面点的温度关于时间的变化率的绝对值不大于每小时0.1摄氏度;第五项条件是在获得索结构稳态温度数据的时刻,R个索结构表面点中的每一个索结构表面点的索结构表面温度实测数据为当日日出时刻到次日日出时刻后30分钟之间的极小值;第六项条件是在获得索结构稳态温度数据的时刻,“索结构厚度方向最大温差” ATtmax不大于I摄氏度;本方法利用上述六项条件,将下列三种时刻中的任意一种称为“获得索结构稳态温度数据的数学时刻”,第一种时刻是满足上述“与决定获得索结构稳态温度数据的时刻相关的条件”中的第一项至第五项条件的时刻,第二种时刻是仅仅满足上述“与决定获得索结构稳态温度数据的时刻相关的条件”中的第六项条件的时刻,第三种时刻是同时满足上述“与决定获得索结构稳态温度数据的时刻相关的条件”中的第一项至第六项条件的时刻;当获得索结构稳态温度数据的数学时刻就是本方法中实际记录数据时刻中的一个时,获得索结构稳态温度数据的时刻就是获得索结构稳态温度数据的数学时刻;如果获得索结构稳态温度数据的数学时刻不是本方法中实际记录数据时刻中的任一个时刻,则取本方法最接近于获得索结构稳态温度数据的数学时刻的那个实际记录数据的时刻为获得索结构稳态温度数据的时刻;本方法将使用在获得索结构稳态温度数据的时刻测量记录的量进行索结构相关健康监测分析;本方法近似认为获得索结构稳态温度数据的时刻的索结构温度场处于稳态,即此时刻的索结构温度不随时间变化,此时刻就是本方法的“获得索结构稳态温度数据的时刻”;然后,根据索结构传热特性,利用获得索结构稳态温度数据的时刻的“R个索结构表面温度实测数据”和“HBE个索结构沿厚度温度实测数据”,利用索结构的传热学计算模型,通过常规传热计算得到在获得索结构稳态温度数据的时刻的索结构的温度分布,此时索结构的温度场按稳态进行计算,计算得到的在获得索结构稳态温度数据的时刻的索结构的温度分布数据包括索结构上R个索结构表面点的计算温度,R个索结构表面点的计算温度称为R个索结构稳态表面温度计算数据,还包括索结构在前面选定的HBE个“测量索结构沿厚度的温度分布数据的点”的计算温度,HBE个“测量索结构沿厚度的温度分布数据的点”的计算温度称为“HBE个索结构沿厚度温度计算数据”,当R个索结构表面温度实测数据与R个索结构稳态表面温度计算数据对应相等时,且“HBE个索结构沿厚度温度实测数据”与“HBE个索结构沿厚度温度计算数据”对应相等时,计算得到的在获得索结构稳态温度数据的时刻的索结构的温度分布数据在本方法中称为“索结构稳态温度数据”,此时的“R个索结构表面温度实测数据”称为“R个索结构稳态表面温度实测数据”,“HBE个索结构沿厚度温度实测数据”称为“HBE个索结构沿厚度稳态温度实测数据”;在索结构的表面上取“R个索结构表面点”时,“R个索结构表面点”的数量与分布必须满足三个条件,第一个条件是当索结构温度场处于稳态时,当索结构表面上任意一点的温度是通过“R个索结构表面点”中与索结构表面上该任意点相邻的点的实测温度线性插值得到时,线性插值得到的索结构表面上该任意点的温度与索结构表面上该任意点的实际温度的误差不大于5% ;索结构表面包括支承索表面;第二个条件是“R个索结构表面点”中在同一海拔高度的点的数量不小于4,且“R个索结构表面点”中在同一海拔高度的点沿着索结构表面均布;“R个索结构表面点”沿海拔高度的所有两两相邻索结构表面点的海拔高度之差的绝对值中的最大值Ah不大于0.2°C除以Λ Th得到的数值,为方便叙述取Λ Th的单位为。C/m,为方便叙述取Ah的单位为m ;“R个索结构表面点”沿海拔高度的两两相邻索结构表面点的定义是指只考虑海拔高度时,在“R个索结构表面点”中不存在一个索结构表面点,该索结构表面点的海拔高度数值介于两两相邻索结构表面点的海拔高度数值之间;第三个条件是查询或按气象学常规计算得到索结构所在地和所在海拔区间的日照规律,再根据索结构的几何特征及方位数据,在索结构上找到全年受日照时间最充分的那些表面点的位置,“R个索结构表面点”中至少有一个索结构表面点是索结构上全年受日照时间最充分的那些表面点中的一个占.[0018] c.按照“本方法的索结构的温度测量计算方法”直接测量计算得到初始状态下的索结构稳态温度数据,初始状态下的索结构稳态温度数据称为初始索结构稳态温度数据,记为“初始索结构稳态温度数据向量T。”;实测或查资料得到索结构所使用的各种材料的随温度变化的物理和力学性能参数;在实测得到初始索结构稳态温度数据向量T。的同一时刻,直接测量计算得到所有支承索的初始索力,组成初始索力向量F。;依据包括索结构设计数据、竣工数据在内的数据得到所有支承索在自由状态即索力为O时的长度、在自由状态时的横截面面积和在自由状态时的单位长度的重量,以及获得这三种数据时所有支承索的温度,在此基础上利用所有支承索的随温度变化的物理性能参数和力学性能参数,按照常规物理计算得到所有支承索在初始索结构稳态温度数据向量T。条件下的索力为O时所有支承索的长度、索力为O时所有支承索的横截面面积以及索力为O时所有支承索的单位长度的重量,依次组成支承索的初始自由长度向量、初始自由横截面面积向量和初始自由单位长度的重量向量,支承索的初始自由长度向量、初始自由横截面面积向量和初始自由单位长度的重量向量的元素的编号规则与初始索力向量F。的元素的编号规则相同;在实测得到T。的同时,也就是在获得初始索结构稳态温度数据向量T。的时刻的同一时刻,直接测量计算得到初始索结构的实测数据,初始索结构的实测数据是包括索结构集中载荷测量数据、索结构分布载荷测量数据、索结构体积载荷测量数据、所有被监测量的初始数值、所有支承索的初始索力数据、初始索结构模态数据、初始索结构应变数据、初始索结构几何数据、初始索结构支座广义坐标数据、初始索结构角度数据、初始索结构空间坐标数据在内的实测数据,在得到初始索结构的实测数据的同时,测量计算得到包括支承索的无损检测数据在内的能够表达支承索的健康状态的数据,此时的能够表达支承索的健康状态的数据称为支承索初始健康状态数据;所有被监测量的初始数值组成被监测量初始数值向量C。,被监测量初始数值向量C。的编号规则与M个被监测量的编号规则相同;利用支承索初始健康状态数据以及索结构载荷测量数据建立被评估对象初始损伤向量d。,向量d。表示用初始力学计算基准模型A。表示的索结构的被评估对象的初始健康状态;被评估对象初始损伤向量
d。的元素个数等于N,d。的元素与被评估对象是一一对应关系,向量d。的元素的编号规则与被评估对象的编号规则相同;如果d。的某一个元素对应的被评估对象是索系统中的一根支承索,那么d。的该元素的数值代表对应支承索的初始损伤程度,若该元素的数值为0,表示该元素所对应的支承索是完好的,没有损伤的,若其数值为100%,则表示该元素所对应的支承索已经完全丧失承载能力,若其数值介于O和100%之间,则表示该支承索丧失了相应比例的承载能力;如果d。的某一个元素对应的被评估对象是某一个载荷,本方法中取d。的该元素数值为0,代表这个载荷的变化的初始数值为O ;如果没有支承索的无损检测数据及其他能够表达支承索的健康状态的数据时,或者可以认为结构初始状态为无损伤无松弛状态时,向量d。中与支承索相关的各元素数值取O ;
[0019] d.根据索结构的设计图、竣工图和初始索结构的实测数据、支承索初始健康状态数据、索结构集中载荷测量数据、索结构分布载荷测量数据、索结构体积载荷测量数据、索结构所使用的各种材料的随温度变化的物理和力学性能参数、初始索结构稳态温度数据向量T。和和前面步骤得到的所有的索结构数据,建立计入“索结构稳态温度数据”的索结构的初始力学计算基准模型A。,基于A。计算得到的索结构计算数据必须非常接近其实测数据,其间的差异不得大于5% ;对应于A。的“索结构稳态温度数据”就是“初始索结构稳态温度数据向量T。” ;对应于A。的被评估对象健康状态用被评估对象初始损伤向量d。表示;对应于A。的所有被监测量的初始数值用被监测量初始数值向量C。表示;T。和d。是A。的参数,由A。的力学计算结果得到的所有被监测量的初始数值与C。表示的所有被监测量的初始数值相同,因此也可以说C。由A。的力学计算结果组成,在本方法中A。、C。、d。和T。是不变的;
[0020] e.在本方法中,字母i除了明显地表示步骤编号的地方外,字母i仅表示循环次数,即第i次循环;第i次循环开始时需要建立的或已建立的索结构的当前初始力学计算基准模型记为当前初始力学计算基准模型A^A。和Ai。计入了温度参数,可以计算温度变化对索结构的力学性能影响;第i次循环开始时,对应于Aitj的“索结构稳态温度数据”用当前初始索结构稳态温度数据向量Titj表示,向量Titj的定义方式与向量T。的定义方式相同,Τ1。的元素与T。的元素一一对应;第i次循环开始时需要的被评估对象当前初始损伤向量记为Cli0, Cli0表示该次循环开始时索结构Ai0的被评估对象的健康状态,Cli0的定义方式与d。的定义方式相同,dt的元素与d。的元素一一对应;第i次循环开始时,所有被监测量的初始值,用被监测量当前初始数值向量C1。表示,向量C1。的定义方式与向量C。的定义方式相同,Ci0
的元素与C。的元素--对应,被监测量当前初始数值向量C1。表示对应于A1。的所有被监测
量的具体数值Ji。和d、是Ai0的特性参数,Ci0由Aitj的力学计算结果组成;第一次循环开始时,Ai0记为A1。,建立A1。的方法为使A1。等于A。;第一次循环开始时,Ti0记为T1。,建立T1。的方法为使T1。等于T。;第一次循环开始时,Cli0记为d1。,建立d1。的方法为使d1。等于d。;第一次循环开始时,Ci0记为C1。,建立C1。的方法为使C1。等于C。;
[0021] f.从这里进入由第f步到第q步的循环;在结构服役过程中,按照“本方法的索结构的温度测量计算方法”不断实测计算获得索结构稳态温度数据的当前数据,所有“索结构稳态温度数据”的当前数据组成当前索结构稳态温度数据向量Ti,向量Ti的定义方式与向量T。的定义方式相同,Ti的元素与T。的元素一一对应;在实测得到向量!^的同时,实测得到在获得当前索结构稳态温度数据向量Ti的时刻的同一时刻的索结构中所有被监测量的当前值,所有这些数值组成被监测量当前数值向量Ci,向量Ci的定义方式与向量C。的定义方式相同,Ci的元素与C。的元素一一对应,表示相同被监测量在不同时刻的数值;在实测得到当前索结构稳态温度数据向量Ti的同时,对新增加的M2根传感索进行无损检测,从中鉴别出出现损伤或松弛的传感索,依据被监测量编号规则,从本方法之前出现的按照被监测量编号规则编号的各向量中去除与鉴别出的出现损伤或松弛的传感索对应的元素,在本方法之后出现的各向量和矩阵中也不再出现与鉴别出的出现损伤或松弛的传感索对应的元素,在本方法之后提到传感索时不再包括这里被鉴别出出现损伤或松弛的传感索,在本方法之后提到被监测量时不再包括这里被鉴别出出现损伤或松弛的传感索的索力;从索结构上鉴别出几根出现损伤或松弛的传感索,就将M2和M减小同样的数量;在实测得到当前索结构稳态温度数据向量Ti的同一时刻,实测得到索结构中所有M1根支承索的索力数据,所有这些索力数据组成当前索力向量Fi,向量Fi的元素与向量F。的元素的编号规则相同;在实测得到当前索结构稳态温度数据向量Ti的同一时刻,实测计算得到所有M1根支承索的两个支承端点的空间坐标,两个支承端点的空间坐标在水平方向分量的差就是两个支承端点水平距离,所有支承索的两个支承端点水平距离数据组成当前支承索两支承端点水平距离向量,当前支承索两支承端点水平距离向量的元素的编号规则与初始索力向量F。的元素的编号规则相同;
[0022] g.根据当前索结构稳态温度数据向量Ti,按照步骤gl至g3更新当前初始力学计算基准模型A1。、被监测量当前初始数值向量C1。和当前初始索结构稳态温度数据向量f。,而被评估对象当前初始损伤向量f。保持不变;
[0023] gl.比较Ti与T:如果Ti等于f。,则A1。、Ci0和Titj保持不变;否则需要按下列步骤对Ai。和ή进行更新;
[0024] g2.计算Ti与T。的差,Ti与T。的差就是当前索结构稳态温度数据关于初始索结构稳态温度数据的变化,Ti与T。的差用稳态温度变化向量S表示,S等于Ti减去T。,S表示索结构稳态温度数据的变化;
[0025] g3.对A。中的索结构施加温度变化,施加的温度变化的数值就取自稳态温度变化向量S,对A。中的索结构施加的温度变化后得到更新的当前初始力学计算基准模型Ai。,更新Aitj的同时,Ti。所有元素数值也用Ti的所有元素数值对应代替,即更新了 Τ1。,这样就得到了正确地对应于Aitj的Titj ;此时(保持不变;当更新Aitj后,Aitj的索的健康状况用被评估对象当前初始损伤向量f。表示,Aitj的索结构稳态温度用当前索结构稳态温度数据向量Ti表示,更新C1。的方法是:当更新(后,通过力学计算得到Aitj中所有被监测量的、当前的具体数值,这些具体数值组成C1。;
[0026] h.在当前初始力学计算基准模型Aitj的基础上,按照步骤hi至步骤h4进行若干次力学计算,通过计算建立单位损伤被监测量数值变化矩阵ACi和被评估对象单位变化向
量 Diu;
[0027] h1.在第i次循环开始时,直接按步骤h2至步骤h4所列方法获得ACi和Diu;在其它时刻,当在步骤g中对Ai0进行更新后,必须按步骤h2至步骤h4所列方法重新获得Λ Ci和Diu,如果在步骤g中没有对Ai0进行更新,则在此处直接转入步骤i进行后续工作;
[0028] h2.在当前初始力学计算基准模型Aitj的基础上进行若干次力学计算,计算次数数值上等于所有被评估对象的数量N,有N个评估对象就有N次计算;依据被评估对象的编号规则,依次进行计算;每一次计算假设只有一个被评估对象在原有损伤或载荷的基础上再增加单位损伤或载荷单位变化,具体的,如果该被评估对象是索系统中的一根支承索,那么就假设该支承索再增加单位损伤,如果该被评估对象是一个载荷,就假设该载荷再增加载荷单位变化,用Diuk记录这一增加的单位损伤或载荷单位变化,其中k表示增加单位损伤或载荷单位变化的被评估对象的编号,Diuk是被评估对象单位变化向量Diu的一个元素,被评估对象单位变化向量Diu的元素的编号规则与向量d。的元素的编号规则相同;每一次计算中再增加单位损伤或载荷单位变化的被评估对象不同于其它次计算中再增加单位损伤或载荷单位变化的被评估对象,每一次计算都利用力学方法计算索结构的所有被监测量的当前计算值,每一次计算得到的所有被监测量的当前计算值组成一个被监测量计算当前向量;当假设第k个被评估对象再增加单位损伤或载荷单位变化时,用Cidt表示对应的“被监测量计算当前向量”;在本步骤中给各向量的元素编号时,应同本方法中其它向量使用同一编号规则,以保证本步骤中各向量中的任意一个元素,同其它向量中的、编号相同的元素,表达了同一被监测量或同一对象的相关信息K1d5的定义方式与向量C。的定义方式相同,C1ck的兀素与C。的兀素对应;
[0029] h3.每一次计算得到的向量Cidt减去向量C1。得到一个向量,再将该向量的每一个元素都除以本次计算所假设的单位损伤或载荷单位变化数值后得到一个“被监测量的数值变化向量δ CY ;有N个被评估对象就有N个“被监测量的数值变化向量”;
[0030] h4.由这N个“被监测量的数值变化向量”按照N个被评估对象的编号规则,依次组成有N列的“单位损伤被监测量数值变化矩阵单位损伤被监测量数值变化矩阵ACi的每一列对应于一个被监测量单位变化向量;单位损伤被监测量数值变化矩阵ACi的每一行对应于同一个被监测量在不同被评估对象增加单位损伤或载荷单位变化时的不同的单位变化幅度;单位损伤被监测量数值变化矩阵ACi的列的编号规则与向量d。的元素的编号规则相同,单位损伤被监测量数值变化矩阵ACi的行的编号规则与M个被监测量的编号规则相同;
[0031] 1.定义当前名义损伤向量d、和当前实际损伤向量ClSdic^Pdi的元素个数等于被评估对象的数量,Clic和Cli的元素和被评估对象之间是一一对应关系,f。的元素数值代表对应被评估对象的名义损伤程度或名义载荷变化量,Clic和Cli与被评估对象初始损伤向量d。的元素编号规则相同,Clic的元素、Cli的元素与d。的元素是一一对应关系;
[0032] j.依据被监测量当前数值向量Ci同“被监测量当前初始数值向量 <”、“单位损伤被监测量数值变化矩阵AC1”和“当前名义损伤向量屮。”间存在的近似线性关系,该近似线性关系可表达为式1,式I中除(外的其它量均为已知,求解式I就可以算出当前名义损伤向量d1。;
Figure CN103852282AD00221
[0034] k.利用式2表达的当`前实际损伤向量Cli的第k个元素Clik同被评估对象当前初始损伤向量必的第k个元素d、和当前名义损伤向量 < 的第k个元素Clidt间的关系,计算得到当前实际损伤向量Cli的所有元素;
[0035]
Figure CN103852282AD00222
[0036] 式2中k=l, 2, 3,……,N;向量Cli的元素的编号规则与式I中向量d。的元素的编号规则相同表示第i次循环中第k个被评估对象的当前实际健康状态,如果该被评估对象是索系统中的一根支承索,那么Clik表示其当前健康问题的严重程度,有健康问题的支承索可能是松弛索、也可能是受损索,Clik数值反应了该支承索的松弛或损伤的程度;如果该被评估对象是一个载荷,那么Clik表示该载荷的实际变化量;将被评估对象当前实际损伤向量Cli中与M1根支承索相关的M1个元素取出,组成支承索当前实际损伤向量Cfi,支承索当前实际损伤向量Cfi的元素的编号规则与初始索力向量F。的元素的编号规则相同;支承索当前实际损伤向量Cfi的第h个元素表示索结构中第h根支承索的当前实际损伤量,h=l, 2,3,…….,M1 ;支承索当前实际损伤向量Cfi中数值不为O的元素对应于有健康问题的支承索,从这些有健康问题的支承索中鉴别出受损索,剩下的就是松弛索;与受损索对应的支承索当前实际损伤向量Cfi中的元素的数值表达的是该受损索的当前实际损伤,元素数值为100%时表示该支承索彻底丧失承载能力,介于O与100%之间时表示丧失相应比例的承载能力;利用在当前索结构稳态温度数据向量Ti条件下的、在第I步鉴别出的松弛索及用支承索当前实际损伤向量Cfi表达的这些松弛索的、与其松弛程度力学等效的当前实际等效损伤程度,利用在第f步获得的在当前索结构稳态温度数据向量Ti条件下的当前索力向量Fi和当前支承索两支承端点水平距离向量,利用在第C步获得的在初始索结构稳态温度数据向量T。条件下的支承索的初始自由长度向量、初始自由横截面面积向量和初始自由单位长度的重量向量、初始索力向量F。,利用当前索结构稳态温度数据向量Ti表示的支承索当前稳态温度数据,利用在第c步获得的在初始索结构稳态温度数据向量T。表示的支承索初始稳态温度数据,利用在第c步获得的索结构所使用的各种材料的随温度变化的物理和力学性能参数,计入温度变化对支承索物理、力学和几何参数的影响,通过将松弛索同受损索进行力学等效来计算松弛索的、与当前实际等效损伤程度等效的松弛程度,力学等效条件是:一、两等效的索的无松弛和无损伤时的初始自由长度、几何特性参数、密度及材料的力学特性参数相同;二、松弛或损伤后,两等效的松弛索和损伤索的索力和变形后的总长相同;满足上述两个力学等效条件时,这样的两根支承索在索结构中的力学功能就是完全相同的,即如果用等效的松弛索代替受损索后,索结构不会发生任何变化,反之亦然;依据前述力学等效条件求得那些被判定为松弛索的松弛程度,松弛程度就是支承索自由长度的改变量,也就是确定了那些需调整索力的支承索的索长调整量;这样就实现了支承索的松弛识别和损伤识别;计算时所需索力由当前索力向量Fi对应元素给出;本方法将受损索和松弛索统称为有健康问题的支承索,简称为问题索,至此本方法实现了剔除载荷变化和结构温度变化的影响的、索结构的问题索识别,同时实现了剔除结构温度变化和支承索健康状态变化影响的、载荷变化量的识别;
[0037] 1.在求得当前名义损伤向量后,按照式3建立标识向量Bi,式4给出了标识向量#的第k个元素的定义;
Figure CN103852282AD00231
[0040] 式4中元素Bik是标识向量Bi的第k个元素,Diuk是被评估对象单位变化向量Diu的第k个元素,Click是被评估对象当前名义损伤向量的第k个元素,它们都表示第k个被评估对象的相关信息,式4中k=l,2,3,......,N ;
[0041] m.如果标识向量Bi的元素全为0,则回到步骤f继续本次循环;如果标识向量Bi的元素不全为0,则进入下一步、即步骤η ;
[0042] η.根据式5计算得到下一次、即第i+Ι次循环所需的被评估对象当前初始损伤向量di+1。的每一个元素;
[0043]
Figure CN103852282AD00232
[0044] 式5中(T1。,是下一次、即第i+Ι次循环所需的被评估对象当前初始损伤向量di+1。的第k个元素,Cliok是本次、即第i次循环的被评估对象当前初始损伤向量f。的第k个元素,Diuk是第i次循环的被评估对象单位变化向量Diu的第k个元素,Bik是第i次循环的标识向量Bi的第k个元素,式5中k=l,2,3,……,N ;
[0045] ο.取下一次、即第i+Ι次循环所需的当前初始索结构稳态温度数据向量Ti+1。等于第i次循环的当前初始索结构稳态温度数据向量Titj ;
[0046] p.在初始力学计算基准模型A。的基础上,对A。中的索结构施加温度变化,施加的温度变化的数值就取自稳态温度变化向量S,再令索的健康状况为di+1。后得到的就是下一次、即第i+Ι次循环所需的力学计算基准模型Ai+1 ;得到Ai+1后,通过力学计算得到Ai+1中所有被监测量的、当前的具体数值,这些具体数值组成下一次、即第i+Ι次循环所需的被监测量当前初始数值向量Ci+1。;
[0047] q.回到步骤f,开始下一次循环。
[0048] 有益效果:本方法实现了已有方法不可能具备的两种功能,分别是:一、在结构承受的载荷和结构(环境)温度变化时,能够剔除载荷变化和结构温度变化对索结构健康状态识别结果的影响,从而准确地识别出问题索的结构健康监测方法;二、本方法在识别出问题索的同时,还能同时识别出载荷的变化,即本方法能够剔除结构温度变化和支承索健康状态变化的影响,实现载荷变化程度的正确识别。
具体实施方式
[0049] 本方法的实施例的下面说明实质上仅仅是示例性的,并且目的绝不在于限制本方法的应用或使用。具体实施时,下列步骤是可采取的各种步骤中的一种。
[0050] 第一步:首先确认索结构承受的可能发生变化的载荷的数量。根据索结构所承受的载荷的特点,确认其中“所有可能发生变化的载荷”,或者将所有的载荷视为“所有可能发生变化的载荷”,设共有JZW个可能发生变化的载荷,本方法通过识别这JZW个“所有可能发生变化的载荷”的变化程度来表达“所有可能发生变化的载荷”的变化量。
[0051] 设索结构的支承索的数量和JZW个“所有可能发生变化的载荷”的数量之和为N。为叙述方便起见,本方法统一称被评估的支承索和“所有可能发生变化的载荷”为“被评估对象”,共有N个被评估对象。给被评估对象连续编号,该编号在后续步骤中将用于生成向量和矩阵。
[0052] 设索系统中共有M1根支承索,结构索力数据包括这M1根支承索的索力,显然M1小于被评估对象的数量N。仅仅通过M1根支承索的M1个索力数据来求解未知的N个被评估对象的状态是不可能的,本方法在监测全部M1根支承索索力的基础上,增加对不少于(N-M1)个其他被监测量。 [0053] 增加的不少于(N-M1)个的其他被监测量仍然是索力,叙述如下:
[0054] 在结构健康检测系统开始工作前,先在索结构上人为增加M2 (M2不小于N-M1)根索,称为传感索,新增加的M2根传感索的刚度同索结构的任意一根支承索的刚度相比,应当小很多,例如小20倍,新增加的M2根传感索的索力应当较小,例如其横截面正应力应当小于其疲劳极限,这些要求可以保证新增加的M2根传感索不会发生疲劳损伤,新增加的M2根传感索的两端应当充分锚固,保证不会出现松弛,新增加的M2根传感索应当得到充分的防腐蚀保护,保证新增加的M2根传感索不会发生损伤和松弛,在结构健康监测过程中将监测这新增加的M2根传感索的索力。
[0055] 还可以采用多增加传感索的方式来保证健康监测的可靠性,例如使仏不小于N-M1的2倍,在结构健康监测过程中只挑选其中的完好的传感索的索力数据(称为实际可以使用的被监测量,记录其数量为K,K不得小于N)和对应的索结构被监测量单位变化矩阵AC进行健康状态评估,由于M2不小于N-M1的2倍,可以保证实际可以使用的有效传感索的数量加上M1不小于N。在结构健康监测过程中将监测这新增加的M2根传感索的索力。新增加的M2根传感索应当安装在结构上、人员易于到达的部位,便于人员对其进行无损检测。
[0056] 在本方法中新增加的M2根传感索作为索结构的一部分,后文再提到索结构时,索结构包括增加M2根传感索前的索结构和新增加的M2根传感索,也就是说后文提到索结构时指包括新增加的M2根传感索的索结构。因此后文提到按照“本方法的索结构的温度测量计算方法”测量计算得到“索结构稳态温度数据”时,其中的索结构包括新增加的M2根传感索,得到的“索结构稳态温度数据”包括新增加的M2根传感索的稳态温度数据,获得新增加的M2根传感索的稳态温度数据的方法同于索结构的M1根支承索的稳态温度数据的获得方
法,在后文不再--交代;测量得到新增加的M2根传感索的索力的方法同于索结构的M1根
支承索的索力的测量方法,在后文不再一一交代;对索结构的支承索进行任何测量时,同时对新增加的M2根传感索进行同样的测量,在后文不再一一交代;新增加的M2根传感索除了不发生损伤和松弛外,新增加的M2根索的信息量与索结构的支承索的信息量相同,在后文不再一一交代;新增加的M2根传感索的索力就是增加的不少于(N-M1)个的其他被监测量。在后文建立索结构的各种力学模型时,将新增加的M2根传感索视同索结构的M1根支承索对待,除了提到支承索的损伤和松弛的场合,在其他场合提到支承索时包括新增加的M2根索。
[0057] 综合上述被监测 量,整个索结构共有M (M=MJM2)根索的M个被监测量,M不得小于被评估对象的数量N。
[0058] 为方便起见,在本方法中将“索结构的被监测的所有参量”简称为“被监测量”。给M个被监测量连续编号,该编号在后续步骤中将用于生成向量和矩阵。本方法用用变量j表示这一编号,j=l, 2, 3,…,Μ。
[0059] 确定“本方法的索结构的温度测量计算方法”,该方法具体步骤同于本方法“权利要求书”中“本方法的索结构的温度测量计算方法”的具体步骤。
[0060] 第二步:建立初始力学计算基准模型Α。。
[0061] 在索结构竣工之时,或者在建立健康监测系统前,按照“本方法的索结构的温度测量计算方法”测量计算得到“索结构稳态温度数据”(可以用常规温度测量方法测量,例如使用热电阻测量),此时的“索结构稳态温度数据”用向量Τ。表示,称为初始索结构稳态温度数据向量Τ。。在实测得到Τ。的同时,也就是在获得初始索结构稳态温度数据向量的时刻的同一时刻,使用常规方法直接测量计算得到索结构的所有被监测量的初始数值,组成被监测量初始数值向量C。。
[0062] 本方法中可以具体按照下列方法在获得某某(例如初始或当前等)索结构稳态温度数据向量的时刻的同一时刻,使用某某方法测量计算得到某某被测量量被监测量(例如索结构的所有被监测量)的数据:在测量记录温度(包括索结构所在环境的气温、参考平板的向阳面的温度和索结构表面温度)的同时,例如每隔10分钟测量记录一次温度,那么同时同样也每隔10分钟测量记录某某被测量量被监测量(例如索结构的所有被监测量)的数据。一旦确定了获得索结构稳态温度数据的时刻,那么与获得索结构稳态温度数据的时刻同一时刻的某某被测量量被监测量(例如索结构的所有被监测量)的数据就称为在获得索结构稳态温度数据的时刻的同一时刻,使用某某方法测量计算方法得到的某某被测量量被监测量的数据。
[0063] 使用常规方法(查资料或实测)得到索结构所使用的各种材料的随温度变化的物理参数(例如热膨胀系数)和力学性能参数(例如弹性模量、泊松比)。
[0064] 在实测得到初始索结构稳态温度数据向量T。的同一时刻,直接测量计算得到所有支承索的初始索力,组成初始索力向量F。;依据索结构设计数据、竣工数据得到所有支承索在自由状态即索力为O时的长度、在自由状态时的横截面面积和在自由状态时的单位长度的重量,以及获得这三种数据时所有支承索的温度,在此基础上利用所有支承索的随温度变化的物理性能参数和力学性能参数,按照常规物理计算得到所有支承索在初始索结构稳态温度数据向量T。条件下的索力为O时所有支承索的长度、索力为O时所有支承索的横截面面积以及索力为O时所有支承索的单位长度的重量,依次组成支承索的初始自由长度向量10、初始自由横截面面积向量A。和初始自由单位长度的重量向量ω。,支承索的初始自由长度向量I。、初始自由横截面面积向量Α。和初始自由单位长度的重量向量ω。的元素的编号规则与初始索力向量F。的元素的编号规则相同。
[0065] 在实测计算得到初始索结构稳态温度数据向量Τ。的同时,使用常规方法实测计算得到索结构的实测计算数据。索结构的实测计算数据包括支承索的无损检测数据等能够表达索的健康状态的数据、索结构初始几何数据、索力数据、拉杆拉力数据、索结构支座广义坐标数据、索结构集中载荷测量数据、索结构分布载荷测量数据、索结构体积载荷测量数据、索结构模态数据、结构应变数据、结构角度测量数据、结构空间坐标测量数据等实测数据。索结构的初始几何数据可以是所有索的端点的空间坐标数据加上结构上一系列的点的空间坐标数据,目的在于根据这些坐标数据确定索结构的几何特征。对斜拉桥而言,初始几何数据可以是所有索的端点的空间坐标数据加上桥梁两端上若干点的空间坐标数据,这就是所谓的桥型数据。利用支承索的无损检测数据等能够表达支承索的健康状态的数据以及索结构载荷测量数据建立被评估对象初始损伤向量d。,用d。表示索结构(用初始力学计算基准模型A。表示)的被评估对象的初始健康状态。如果没有支承索的无损检测数据及其他能够表达支承索的健康状态的数据时,或者可以认为结构初始状态为无损伤无松弛状态时,向量d。的中与支承索相关的各元素数值取0,如果d。的某一个元素对应的被评估对象是某一个载荷,本方法中取d。的该元素数值为0,代表这个载荷的变化的初始数值为O。利用索结构的设计图、竣工图和初始索结构的实测数据、支承索的无损检测数据、索结构所使用的各种材料的随温度变化的物理和力学性能参数和初始索结构稳态温度数据向量T。,利用力学方法(例如有限元法)计入“索结构稳态温度数据”建立初始力学计算基准模型A。。
[0066] 不论用何种方法获得初始力学计算基准模型A。,计入“索结构稳态温度数据”(即初始索结构稳态温度数据向量T。)、基于A。计算得到的索结构计算数据必须非常接近其实测数据,误差一般不得大于5%。这样可保证利用A。计算所得的模拟情况下的索力计算数据、应变计算数据、索结构形状计算数据和位移计算数据、索结构角度数据、索结构空间坐标数据等,可靠地接近所模拟情况真实发生时的实测数据。模型A。中支承索的健康状态用被评估对象初始损伤向量d。表示,索结构稳态温度数据用初始索结构稳态温度数据向量T。表示。由于基于A。计算得到所有被监测量的计算数值非常接近所有被监测量的初始数值(实测得到),所以也可以用在A。的基础上、进行力学计算得到的、A0的每一个被监测量的计算数值组成被监测量初始数值向量C。。对应于A。的“索结构稳态温度数据”就是“初始索结构稳态温度数据向量T。” ;对应于A。的被评估对象健康状态用被评估对象初始损伤向量d。表示;对应于A。的所有被监测量的初始数值用被监测量初始数值向量C。表示。T。和d。是A0的参数,C0由A。的力学计算结果组成。
[0067] 第三步:在本方法中,字母i除了明显地表示步骤编号的地方外,字母i仅表示循环次数,即第i次循环;第i次循环开始时需要建立的或已建立的索结构的当前初始力学计算基准模型记为当前初始力学计算基准模型A^A。和Ai。计入了温度参数,可以计算温度变化对索结构的力学性能影响;第i次循环开始时,对应于Aitj的“索结构稳态温度数据”用当前初始索结构稳态温度数据向量Titj表示,向量Titj的定义方式与向量T。的定义方式相同,Ti0的元素与T。的元素一一对应;第i次循环开始时需要的被评估对象当前初始损伤向量记为d1。,Cli0表示该次循环开始时索结构Aitj的被评估对象的健康状态,Cli0的定义方式与d。的定义方式相同,d、的元素与d。的元素一一对应;第i次循环开始时,所有被监测量的初始值,用被监测量当前初始数值向量C1。表示,向量C1。的定义方式与向量C。的定义方式相同,Citj的元素与C。的元素一一对应,被监测量当前初始数值向量Citj表示对应于Ai0的所有被监测量的具体数值和必是Aitj的特性参数由Ai0的力学计算结果组成;第一次循环开始时,Ai0记为A1。,建立A1。的方法为使A1。等于A。;第一次循环开始时,Ti0记为T1。,建立T1。的方法为使T10等于T。;第一次循环开始时,Cli0记为d1。,建立d1。的方法为使d1。等于d。;第一次循环开始时,Ci0记为C1。,建立C1。的方法为使C1。等于C。。
[0068] 第四步:安装索结构健康监测系统的硬件部分。硬件部分至少包括:被监测量监测系统(例如含索力测量系统、信号调理器等)、索结构温度监测系统(含温度传感器、信号调理器等)和索结构环境温度测量系统(含温度传感器、信号调理器等)、支承索的支承端点的空间坐标监测系统、信号(数据)采集器、计算机和通信报警设备。每一个被监测量、每一个温度、每一根支承索的索力、每一根支承索的支承端点的空间坐标都必须被监测系统监测到,监测系统将监测到的信号传输到信号(数据)采集器;信号经信号采集器传递到计算机;计算机则负责运行索结构的被评估对象的健康监测软件,包括记录信号采集器传递来的信号;当监测到被评估对象健康状态有变化时,计算机控制通信报警设备向监控人员、业主和(或)指定的人员报警。
[0069] 第五步:编制并在计算机上安装运行本方法的系统软件,该软件将完成本方法任务所需要的监测、记录、控制、存储、计算、通知、报警等功能(即本具体实施方法中所有可以用计算机完成的工作)。
[0070]第六步:由此步开始循环运作,在结构服役过程中,按照“本方法的索结构的温度测量计算方法”不断实测计算获得索结构稳态温度数据的当前数据,所有“索结构稳态温度数据”的当前数据组成当前索结构稳态温度数据向量Ti,向量Ti的定义方式与向量T。的定义方式相同,Ti的元素与T。的元素一一对应;在实测向量Ti的同时,也就是在获得当前索结构稳态温度数据向量Ti的时刻的同一时刻,实测得到索结构中所有被监测量的当前值,所有这些数值组成被监测量当前数值向量Ci,向量Ci的定义方式与向量C。的定义方式相同,Ci的元素与C。的元素一一对应,表示相同被监测量在不同时刻的数值。[0071] 在实测得到当前索结构稳态温度数据向量Ti的同一时刻,实测得到索结构中所有M1根支承索的索力数据,所有这些索力数据组成当前索力向量F1,向量F1的元素与向量F。的元素的编号规则相同;在实测得到当前索结构稳态温度数据向量Ti的同一时刻,实测计算得到所有M1根支承索的两个支承端点的空间坐标,两个支承端点的空间坐标在水平方向分量的差就是两个支承端点水平距离,所有M1根支承索的两个支承端点水平距离数据组成当前支承索两支承端点水平距离向量Iix,当前支承索两支承端点水平距离向量Iix的元素的编号规则与初始索力向量F。的元素的编号规则相同。
[0072] 在实测得到当前索结构稳态温度数据向量Ti的同时,对新增加的M2根传感索进行无损检测,例如超声波探伤、目视检查、红外成像检查,从中鉴别出出现损伤或松弛的传感索,依据被监测量编号规则,从本方法之前出现的按照被监测量编号规则编号的各向量中去除与鉴别出的出现损伤或松弛的传感索对应的元素,在本方法之后出现的各向量和矩阵中也不再出现与鉴别出的出现损伤或松弛的传感索对应的元素,在本方法之后提到传感索时不再包括这里被鉴别出出现损伤或松弛的传感索,在本方法之后提到被监测量时不再包括这里被鉴别出出现损伤或松弛的传感索的索力;从索结构上鉴别出几根出现损伤或松弛的传感索,就将M2和M减小同样的数量。
[0073] 第七步:在得到当前索结构稳态温度数据向量疒后,比较疒和砍,如果Ti等于f。,则不需要对Ai0和Titj进行更新,否则需要对当前初始力学计算基准模型Α1。、当前初始索结构稳态温度数据向量Titj和被监测量当前初始数值向量C1。进行更新,而被评估对象当前初始损伤向量f。保持不变,更新方法按技术方案和权利要求书中的规定进行。
[0074] 第八步:在当前初始力学计算基准模型Aitj的基础上,按照技术方案和权利要求书中的规定步骤进行若干次力学计算,通过计算建立单位损伤被监测量数值变化矩阵ACi和被评估对象单位变化向量Diutl其中可以取5%、10%,20%或30%等损伤为单位损伤,载荷单位变化根据载荷的不同情况可以取不同值,如果该载荷是分布载荷,且该分布载荷是线分布载荷,载荷单位变化可以取lkN/m、2kN/m、3kN/m或lkNm/m、2kNm/m、3kNm/m等为单位变化;如果该载荷是分布载荷,且该分布载荷是是面分布载荷,载荷单位变化可以取lMPa、2MPa、3MPa或lkNm/m2、2kNm/m2、3kNm/m2等为单位变化;如果该载荷是集中载荷,且该集中载荷是力偶,载荷单位变化可以取lkNm、2kNm、3kNm等为单位变化;如果该载荷是集中载荷,且该集中载荷是集中力,载荷单位变化可以取lkN、2kN、3kN等为单位变化;如果该载荷是体积载荷,载荷单位变化可以取lkN/m3、2kN/m3、3kN/m3等为单位变化。
[0075] 第九步:建立线性关系误差向量¥和向量g、利用前面的数据(“被监测量当前初始数值向量C1。”、“单位损伤被监测量数值变化矩阵ACU’),在第八步进行每一次计算的同时,即在每一次计算假设被评估对象中只有一个被评估对象的增加单位损伤或载荷单位变
化的同时,当假设第k (k=l, 2,3,......,N)个被评估对象增加单位损伤或载荷单位变化时,
每一次计算组成一个损伤向量,用Clitk表示该损伤向量,对应的被监测量计算当前向量为Citk (参见第八步或 技术方案和权利要求书中的相关内容),损伤向量Clitk的元素个数等于被评估对象的数量,向量Clitk的所有元素中只有一个元素的数值取每一次计算中假设增加单位损伤或载荷单位变化的被评估对象的单位损伤或载荷单位变化值,Clitk的其它元素的数值取0,那个不为O的元素的编号与假定增加单位损伤或载荷单位变化的被评估对象的对应关系、同其他向量的同编号的元素同该被评估对象的对应关系是相同的Mitk与被评估对象初始损伤向量d。的元素编号规则相同,Clitk的元素与d。的元素是--对应关系。将Cιtk、
Cι0, Δ Ci, Clitk带入式(1 ),得到一个线性关系误差向量eιk,每一次计算得到一个线性关系误差向量eιk;eιk的下标k表示第k (k=l,2,3,……,N)个被评估对象增加单位损伤或载荷单位变化。有N个被评估对象就有N次计算,就有N个线性关系误差向量θ\,将这N个线性关系误差向量相加后得到一个向量,将此向量的每一个元素除以N后得到的新向量就是最终的线性关系误差向量ei。向量f等于最终的误差向量将向量¥保存在运行健康监测系统软件的计算机硬盘上,供健康监测系统软件使用。
[0076]
Figure CN103852282AD00291
[0077] 式(I)中abs()是取绝对值函数,对括号内求得的向量的每一个元素取绝对值。
[0078] 第十步:定义当前名义损伤向量dic和当前实际损伤向量di, dic和di的元素个数等于被评估对象的数量, dic和di的元素和被评估对象之间是一一对应关系,dic和di的元素数值代表对应被评估对象的损伤程度或载荷变化程度,di。和di与被评估对象初始损伤向量d。的元素编号规则相同,dic的元素、di的元素与d。的元素是一一对应关系。
[0079] 第十一步:依据被监测量当前数值向量Ci同“被监测量当前初始数值向量Ci。”、“单位损伤被监测量数值变化矩阵AC1”和“当前名义损伤向量屮。”间存在的近似线性关系,该近似线性关系可表达为式(2),按照多目标优化算法计算当前名义损伤向量f。的非劣解,也就是带有合理误差、但可以比较准确地从所有索中确定受损索的位置及其名义损伤程度的解。
[0080]
Figure CN103852282AD00292
[0081] 按照多目标优化算法中的目标规划法(Goal Attainment Method),式(2)可以转化成式(3)和式(4)所示的多目标优化问题,式(3)中Y是一个实数,R是实数域,空间区域Ω限制了向量dic的每一个元素的取值范围。式(3)的意思是寻找一个最小的实数Y,使得式(4)得到满足。式(4)中G(dic)由式(5)定义,式(4)中加权向量W与Y的积表示式(4)中G(Cic)与向量¥之间允许的偏差,其值已在第九步计算得到。实际计算时向量W可以与向量f相同。目标规划法的具体编程实现已经有通用程序可以直接采用。使用目标规划法就可以求得当前名义损伤向量dιc。
Figure CN103852282AD00293
[0085] 第十二步:依据当前实际损伤向量士的定义和其元素的定义计算得到当前实际损伤向量di的每一个元素,从而可由di确定被评估对象的健康状态。被评估对象当前实际损伤向量di的第k个元素dik表示第i次循环中第k个被评估对象的当前实际健康状态。
[0086]dik表示第i次循环中第k个被评估对象的当前实际健康状态,如果该被评估对象是索系统中的一根支承索,那么dik表示其当前实际损伤,dik为O时表示其对应的支承索无健康问题,dik数值不为O时表示其对应的支承索是有健康问题的支承索,有健康问题的支承索可能是松弛索、也可能是受损索,其数值反应了松弛或损伤的程度。[0087] Clik表示第i次循环中第k个被评估对象的当前实际健康状态,如果该被评估对象是一个载荷,那么Clik表示其当前实际载荷变化数值。
[0088] 将被评估对象当前实际损伤向量Cli中与支承索相关的M1个元素取出,组成支承索当前实际损伤向量d'支承索当前实际损伤向量Cfi的元素的编号规则与初始索力向量F。的元素的编号规则相同。支承索当前实际损伤向量Cfi的第h个元素表示索结构中第h根支承索的当前实际损伤量,h=l, 2,3,…….,M1 ;支承索当前实际损伤向量Cfi中数值不为O的元素对应于有健康问题的支承索,采用无损检测方法从这些有健康问题的支承索中鉴别出受损索,剩下的就是就是需调整索力的索,就是松弛索,这些需调整索力的索在支承索当前实际损伤向量Cfi中所对应的元素数值(例如其中一个元素可用Cfih表示)表示与这些支承索的松弛程度力学等效的损伤程度,由此就确定了松弛索,具体松弛量的计算方法在下面说明。与受损索对应的支承索当前实际损伤向量Cfi中的元素的数值表达的是该受损索的当前实际损伤,元素数值为100%时表示该支承索彻底丧失承载能力,介于O与100%之间时表示丧失相应比例的承载能力;受损索在支承索当前实际损伤向量0-中对应的元素的数值就表示其损伤程度,对应元素的数值为100%时表示该支承索彻底丧失承载能力,介于O与100%之间时表示该支承索丧失相应比例的承载能力,至此便识别出了受损索及其损伤程度。
[0089] 通过将松弛索同受损索进行力学等效来计算松弛索的、与当前实际等效损伤程度等效的松弛程度。具体地可以依据式(6)可以求得这些索的松弛程度(即索长调整量)。这样就实现了支承索的松弛识别。至此便全部识别了受损索和松弛索。
Figure CN103852282AD00301
[0091] 式(6)中中Λ Iih是第h根支承索当前松弛程度,Eih是在索结构的稳态温度数据用当前初始索结构稳态温度数据向量Titj表示时,第h根支承索的弹性模量,Aih是在索结构的稳态温度数据用当前初始索结构稳态温度数据向量Titj表示时,第h根支承索的横截面面积,Fih是在索结构的稳态温度数据用当前初始索结构稳态温度数据向量Titj表示时,第h根支承索的当前索力,Cfih是第h根支承索的当前实际损伤程度,ω 是在索结构的稳态温度数据用当前初始索结构稳态温度数据向量Titj表示时,第h根支承索的单位长度的重量,Iixh是在索结构的稳态温度数据用当前初始索结构稳态温度数据向量Titj表示时,第h根支承索的两个支承端点的水平距离,Iixh是当前支承索两支承端点水平距离向量Iix的一个元素,当前支承索两支承端点水平距离向量Iix的元素的编号规则与初始自由长度向量I。的元素的编号规则相同,Eih可以根据查或实测第h根支承索的材料特性数据得到,Aih和ω \可以根据第h根支承索的热膨胀系数、Aa、Ti。通过常规物理和力学计算得到。式
(7)中[]内的项是该支承索的Ernst等效弹性模量。
[0092] 第十三步:健康监测系统中的计算机定期自动或由人员操作健康监测系统生成索系统健康情况报表。
[0093] 第十四步:在指定条件下,健康监测系统中的计算机自动操作通信报警设备向监控人员、业主和(或)指定的人员报警。
[0094] 第十五步:建立标识向量Bi,如果标识向量Bi的元素全为0,则回到第六步继续进行对索系统的健康监测和计算;如果标识向量Bi的元素不全为0,则完成后续步骤后,进入下一次循环。
[0095] 第十六步:计算得到下一次(即第i+Ι次,i=l, 2,3,4,…)循环所需的初始损伤向量di+1。的每一个元素cT1。,(k=l, 2,3,……,N);第二,在初始力学计算基准模型A。的基础上,对A。中的索结构施加温度变化,施加的温度变化的数值就取自稳态温度变化向量S,再令索的健康状况为di+1。后得到的就是下一次、即第i+Ι次(i=l,2,3,4,…)循环所需的力学计算基准模型Ai+1 ;下一次(即第i+Ι次,i=l, 2,3,4,…)循环所需的当前初始索结构稳态温度数据向量Ti+1。等于f。。得到Ai+1、di+1。和Ti+1。后,通过力学计算得到Ai+1中所有被监测量的、当前的具体数值,这些具体数值组成下一次、即第i+Ι次循环所需的被监测量当前初始数值向量Ci+1。。
[0096] 第十七步:回到第六步,开始由第六步到第十七步的循环。

Claims (1)

1.索力监测问题索载荷递进式识别方法,其特征在于所述方法包括: a.为叙述方便起见,本方法统一称被评估的支承索和载荷为被评估对象,设被评估的支承索的数量和载荷的数量之和为N,即被评估对象的数量为N ;确定被评估对象的编号规贝U,按此规则将索结构中所有的被评估对象编号,该编号在后续步骤中将用于生成向量和矩阵;本方法用变量k表示这一编号,k=l, 2,3,…,N ;设索系统中共有M1根支承索,索结构索力数据包括这M1根支承索的索力,显然M1小于被评估对象的数量N ;仅仅通过M1个支承索的M1个索力数据来求解未知的N个被评估对象的状态是不可能的,本方法在监测全部M1根支承索索力的基础上,在索结构上人为增加M2根索,称为传感索,在索结构健康监测过程中将监测这新增加的M2根传感索的索力;综合上述被监测量,整个索结构共有M根索的M个索力被监测,即有M个被监测量,其中M为M1与M2之和;M应当大于被评估对象的数量N ;新增加的M2根传感索的刚度同索结构的任意一根支承索的刚度相比,应当小得多;新增加的M2根传感索的各传感索的索力应当比索结构的任意一根支承索的索力小得多,这样可以保证即使这新增加的M2根传感索出现了损伤或松弛,对索结构其他构件的应力、应变、变形的影响微乎其微;新增加的M2根传感索的横截面上正应力应当小于其疲劳极限,这些要求可以保证新增加的M2根传感索不会发生疲劳损伤;新增加的M2根传感索的两端应当充分锚固,保证不会出现松弛;新增加的M2根传感索应当得到充分的防腐蚀保护,保证新增加的M2根传感索不会发生损伤和松弛;为方便起见,在本方法中将“索结构的被监测的所有参量”简称为“被监测量”;给M个被监测量连续编号,本方法用用变量j表示这一编号,j=l, 2,3,…,M,该编号在后续步骤中将用于生成向量和矩阵;在本方法中新增加的M2根传感索作为索结构的一部分,后文再提到索结构时,索结构包括增加M2根传感索前的索结构和新增加的M2根传感索,也就是说后文提到索结构时指包括新增加的M2根传感索的索结构;因此后文提到按照“本方法的索结构的温度测量计算方法”测量计算得到“索结构稳态温度数据”时,其中的索结构包括新增加的M2根传感索,得到的“索结构稳态温度数据”包括新增加的M2根传感索的稳态温度数据,获得新增加的M2根传感索的稳态温度数据的方法同于索结构的M1根支承索的稳态温度数据的获得方法,在后文不再一一交代;测量得到新增加的M2根传感索的索力的方法同于索结构的M1根支承索的索力的测量方法,在后文不再一一交代;对索结构的支承索进行任何测量时`,同时对新增加的M2根传感索进行同样的测量,在后文不再一一交代;新增加的仏根传感索除了不发生损伤和松弛外,对新增加的仏根传感索的信息量的要求和获得方法与索结构的支承索的信息量的要求和获得方法相同,在后文不再--交代;在后文建立索结构的各种力学模型时,将新增加的M2根传感索视同索结构的支承索对待;在后文中,除了提到支承索的损伤和松弛的场合外,当提到支承索时所说的支承索包括索结构的支承索和新增加的M2根传感索;本方法中对同一个量实时监测的任何两次测量之间的时间间隔不得大于30分钟,测量记录数据的时刻称为实际记录数据时刻;物体、结构承受的外力可称为载荷,载荷包括面载荷和体积载荷;面载荷又称表面载荷,是作用于物体表面的载荷,包括集中载荷和分布载荷两种;体积载荷是连续分布于物体内部各点的载荷,包括物体的自重和惯性力在内;集中载荷分为集中力和集中力偶两种,在包括笛卡尔直角坐标系在内的坐标系中,一个集中力可以分解成三个分量,同样的,一个集中力偶也可以分解成三个分量,如果载荷实际上是集中载荷,在本方法中将一个集中力分量或一个集中力偶分量计为或统计为一个载荷,此时载荷的变化具体化为一个集中力分量或一个集中力偶分量的变化;分布载荷分为线分布载荷和面分布载荷,分布载荷的描述至少包括分布载荷的作用区域和分布载荷的大小,分布载荷的大小用分布集度来表达,分布集度用分布特征和幅值来表达;如果载荷实际上是分布载荷,本方法谈论载荷的变化时,实际上是指分布载荷分布集度的幅值的改变,而所有分布载荷的作用区域和分布集度的分布特征是不变的;在包括笛卡尔直角坐标系在内的坐标系中,一个分布载荷可以分解成三个分量,如果这分布载荷的三个分量的各自的分布集度的幅值发生变化,且变化的比率不全部相同,那么在本方法中把这分布载荷的三个分量计为或统计为三个分布载荷,此时一个载荷就代表分布载荷的一个分量;体积载荷是连续分布于物体内部各点的载荷,体积载荷的描述至少包括体积载荷的作用区域和体积载荷的大小,体积载荷的大小用分布集度来表达,分布集度用分布特征和幅值来表达;如果载荷实际上是体积载荷,在本方法中实际处理的是体积载荷分布集度的幅值的改变,而所有体积载荷的作用区域和分布集度的分布特征是不变的,此时在本方法中提到载荷的改变时实际上是指体积载荷的分布集度的幅值的改变,此时,发生变化的载荷是指那些分布集度的幅值发生变化的体积载荷;在包括笛卡尔直角坐标系在内的坐标系中,一个体积载荷可以分解成三个分量,如果这体积载荷的三个分量的各自的分布集度的幅值发生变化,且变化的比率不全部相同,那么在本方法中把这体积载荷的三个分量计为或统计为三个分布载荷; b.本方法定义“本方法的索结构的温度测量计算方法”按步骤bl至b3进行;b 1:查询或实测得到索结构组成材料及索结构所处环境的随温度变化的传热学参数,利用索结构的设计图、竣工图和索结构的几何实测数据,利用这些数据和参数建立索结构的传热学计算模型;查询索结构所在地不少于2年的近年来的气象资料,统计得到这段时间内的阴天数量记为T个阴天,在本方法中将白天不能见到太阳的一整日称为阴天,统计得到T个阴天中每一个阴天的O时至次日日出时刻后30分钟之间的最高气温与最低气温,日出时刻是指根据地球自转和公转规律确定的气象学上的日出时刻,不表示当天一定可以看见太阳,可以查询资料或通过常规气象学计算得到所需的每一日的日出时刻,每一个阴天的O时至次日日出时刻后30分钟之间的最高气温减去最低气温称为该阴天的日气温的最大温差,有T个阴天,就有T个阴天的日气温的最大温差,取T个阴天的日气温的最大温差中的最大值为参考日温差,参考日温差记为△ I;;查询索结构所在地和所在海拔区间不少于2年的近年来的气象资料或实测得到索结构所处环境的温度随时间和海拔高度的变化数据和变化规律,计算得到索结构所在地`和所在海拔区间不少于2年的近年来的索结构所处环境的温度关于海拔高度的最大变化率ATh,为方便叙述取ATh的单位为。C/m;在索结构的表面上取“R个索结构表面点”,取“R个索结构表面点”的具体原则在步骤b3中叙述,后面将通过实测得到这R个索结构表面点的温度,称实测得到的温度数据为“R个索结构表面温度实测数据”,如果是利用索结构的传热学计算模型,通过传热计算得到这R个索结构表面点的温度,就称计算得到的温度数据为“R个索结构表面温度计算数据”;从索结构所处的最低海拔到最高海拔之间,在索结构上均布选取不少于三个不同的海拔高度,在每一个选取的海拔高度处、在水平面与索结构表面的交线处至少选取两个点,从选取点处引索结构表面的外法线,所有选取的外法线方向称为“测量索结构沿壁厚的温度分布的方向”,测量索结构沿壁厚的温度分布的方向与“水平面与索结构表面的交线”相交,在选取的测量索结构沿壁厚的温度分布的方向中必须包括索结构的向阳面外法线方向和索结构的背阴面外法线方向,沿每一个测量索结构沿壁厚的温度分布的方向在索结构中均布选取不少于三个点,测量所有被选取点的温度,测得的温度称为“索结构沿厚度的温度分布数据”,其中沿与同一“水平面与索结构表面的交线”相交的、“测量索结构沿壁厚的温度分布的方向”测量获得的“索结构沿厚度的温度分布数据”,在本方法中称为“相同海拔高度索结构沿厚度的温度分布数据”,设选取了 H个不同的海拔高度,在每一个海拔高度处,选取了 B个测量索结构沿壁厚的温度分布的方向,沿每个测量索结构沿壁厚的温度分布的方向在索结构中选取了 E个点,其中H和E都不小于3,B不小于2,设HBE为H与B和E的乘积,对应的共有HBE个“测量索结构沿厚度的温度分布数据的点”,后面将通过实测得到这HBE个“测量索结构沿厚度的温度分布数据的点”的温度,称实测得到的温度数据为“HBE个索结构沿厚度温度实测数据”,如果是利用索结构的传热学计算模型,通过传热计算得到这HBE个测量索结构沿厚度的温度分布数据的点的温度,就称计算得到的温度数据为“HBE个索结构沿厚度温度计算数据” 为B和E的乘积,本方法中在每一个选取的海拔高度处共有BE个“相同海拔高度索结构沿厚度的温度分布数据”;在索结构所在地按照气象学测量气温要求选取一个位置,将在此位置实测得到符合气象学测量气温要求的索结构所在环境的气温;在索结构所在地的空旷无遮挡处选取一个位置,该位置应当在全年的每一日都能得到该地所能得到的该日的最充分的日照,在该位置安放一块碳钢材质的平板,称为参考平板,参考平板与地面不可接触,参考平板离地面距离不小于1.5米,该参考平板的一面向阳,称为向阳面,参考平板的向阳面是粗糙的和深色的,参考平板的向阳面应当在全年的每一日都能得到一块平板在该地所能得到的该日的最充分的日照,参考平板的非向阳面覆有保温材料,将实时监测得到参考平板的向阳面的温度; b2:实时监测得到上述R个索结构表面点的R个索结构表面温度实测数据,同时实时监测得到前面定义的索结构沿厚度的温度分布数据,同时实时监测得到符合气象学测量气温要求的索结构所在环境的气温数据;通过实时监测得到当日日出时刻到次日日出时刻后30分钟之间的索结构所在环境的气温实测数据序列,索结构所在环境的气温实测数据序列由当日日出时刻到次日日出时刻后30分钟之间的索结构所在环境的气温实测数据按照时间先后顺序排列,找到索结构所在环境的气温实测数据序列中的最高温度和最低温度,用索结构所在环境的气温实测数据序列中的最高温度减去最低温度得到索结构所在环境的当日日出时刻到次日日出时刻后·30分钟之间的最大温差,称为环境最大温差,记为ATemax ;由索结构所在环境的气温实测数据序列通过常规数学计算得到索结构所在环境的气温关于时间的变化率,该变化率也随着时间变化;通过实时监测得到当日日出时刻到次日日出时刻后30分钟之间的参考平板的向阳面的温度的实测数据序列,参考平板的向阳面的温度的实测数据序列由当日日出时刻到次日日出时刻后30分钟之间的参考平板的向阳面的温度的实测数据按照时间先后顺序排列,找到参考平板的向阳面的温度的实测数据序列中的最高温度和最低温度,用参考平板的向阳面的温度的实测数据序列中的最高温度减去最低温度得到参考平板的向阳面的温度的当日日出时刻到次日日出时刻后30分钟之间的最大温差,称为参考平板最大温差,记为ATpmax ;通过实时监测得到当日日出时刻到次日日出时刻后30分钟之间的所有R个索结构表面点的索结构表面温度实测数据序列,有R个索结构表面点就有R个索结构表面温度实测数据序列,每一个索结构表面温度实测数据序列由一个索结构表面点的当日日出时刻到次日日出时刻后30分钟之间的索结构表面温度实测数据按照时间先后顺序排列,找到每一个索结构表面温度实测数据序列中的最高温度和最低温度,用每一个索结构表面温度实测数据序列中的最高温度减去最低温度得到每一个索结构表面点的温度的当日日出时刻到次日日出时刻后30分钟之间的最大温差,有R个索结构表面点就有R个当日日出时刻到次日日出时刻后30分钟之间的最大温差数值,其中的最大值称为索结构表面最大温差,记为ATsmax ;由每一索结构表面温度实测数据序列通过常规数学计算得到每一个索结构表面点的温度关于时间的变化率,每一个索结构表面点的温度关于时间的变化率也随着时间变化;通过实时监测得到当日日出时刻到次日日出时刻后30分钟之间的、在同一时刻、HBE个“索结构沿厚度的温度分布数据”后,计算在每一个选取的海拔高度处共计BE个“相同海拔高度索结构沿厚度的温度分布数据”中的最高温度与最低温度的差值,这个差值的绝对值称为“相同海拔高度处索结构厚度方向最大温差”,选取了 H个不同的海拔高度就有H个“相同海拔高度处索结构厚度方向最大温差”,称这H个“相同海拔高度处索结构厚度方向最大温差”中的最大值为“索结构厚度方向最大温差”,记为Δ Ttmax ; b3:测量计算获得索结构稳态温度数据;首先,确定获得索结构稳态温度数据的时刻,与决定获得索结构稳态温度数据的时刻相关的条件有六项,第一项条件是获得索结构稳态温度数据的时刻介于当日日落时刻到次日日出时刻后30分钟之间,日落时刻是指根据地球自转和公转规律确定的气象学上的日落时刻,可以查询资料或通过常规气象学计算得到所需的每一日的日落时刻;第二项条件的a条件是在当日日出时刻到次日日出时刻后30分钟之间的这段时间内,参考平板最大温差Λ Tpmax和索结构表面最大温差ATsmax都不大于5摄氏度;第二项条件的b条件是在当日日出时刻到次日日出时刻后30分钟之间的这段时间内,在前面测量计算得到的环境最大误差Λ Traiax不大于参考日温差Λ Ί;,且参考平板最大温差ATpmax减去2摄氏度后不大于Λ Temax,且索结构表面最大温差ATsmax不大于ATpmax;只需满足第二项的a条件和b条件中的一项就称为满足第二项条件;第三项条件是在获得索结构稳态温度数据的时刻,索结构所在环境的气温关于时间的变化率的绝对值不大于每小时0.1摄氏度;第四项条件是在获得索结构稳态温度数据的时刻,R个索结构表面点中的每一个索结构表面点的温度关于时间的变化率的绝对值不大于每小时0.1摄氏度;第五项条件是在获得索结构稳态温度数据的时刻,R个索结构表面点中的每一个索结构表面点的索结构表面温度实测数据为当日日出时`刻到次日日出时刻后30分钟之间的极小值;第六项条件是在获得索结构稳态温度数据的时刻,“索结构厚度方向最大温差” Λ Ttmax不大于I摄氏度;本方法利用上述六项条件,将下列三种时刻中的任意一种称为“获得索结构稳态温度数据的数学时刻”,第一种时刻是满足上述“与决定获得索结构稳态温度数据的时刻相关的条件”中的第一项至第五项条件的时刻,第二种时刻是仅仅满足上述“与决定获得索结构稳态温度数据的时刻相关的条件”中的第六项条件的时刻,第三种时刻是同时满足上述“与决定获得索结构稳态温度数据的时刻相关的条件”中的第一项至第六项条件的时刻;当获得索结构稳态温度数据的数学时刻就是本方法中实际记录数据时刻中的一个时,获得索结构稳态温度数据的时刻就是获得索结构稳态温度数据的数学时刻;如果获得索结构稳态温度数据的数学时刻不是本方法中实际记录数据时刻中的任一个时刻,则取本方法最接近于获得索结构稳态温度数据的数学时刻的那个实际记录数据的时刻为获得索结构稳态温度数据的时刻;本方法将使用在获得索结构稳态温度数据的时刻测量记录的量进行索结构相关健康监测分析;本方法近似认为获得索结构稳态温度数据的时刻的索结构温度场处于稳态,即此时刻的索结构温度不随时间变化,此时刻就是本方法的“获得索结构稳态温度数据的时刻”;然后,根据索结构传热特性,利用获得索结构稳态温度数据的时刻的“R个索结构表面温度实测数据”和“HBE个索结构沿厚度温度实测数据”,利用索结构的传热学计算模型,通过常规传热计算得到在获得索结构稳态温度数据的时刻的索结构的温度分布,此时索结构的温度场按稳态进行计算,计算得到的在获得索结构稳态温度数据的时刻的索结构的温度分布数据包括索结构上R个索结构表面点的计算温度,R个索结构表面点的计算温度称为R个索结构稳态表面温度计算数据,还包括索结构在前面选定的HBE个“测量索结构沿厚度的温度分布数据的点”的计算温度,HBE个“测量索结构沿厚度的温度分布数据的点”的计算温度称为“HBE个索结构沿厚度温度计算数据”,当R个索结构表面温度实测数据与R个索结构稳态表面温度计算数据对应相等时,且“HBE个索结构沿厚度温度实测数据”与“HBE个索结构沿厚度温度计算数据”对应相等时,计算得到的在获得索结构稳态温度数据的时刻的索结构的温度分布数据在本方法中称为“索结构稳态温度数据”,此时的“R个索结构表面温度实测数据”称为“R个索结构稳态表面温度实测数据”,“HBE个索结构沿厚度温度实测数据”称为“HBE个索结构沿厚度稳态温度实测数据”;在索结构的表面上取“R个索结构表面点”时,“R个索结构表面点”的数量与分布必须满足三个条件,第一个条件是当索结构温度场处于稳态时,当索结构表面上任意一点的温度是通过“R个索结构表面点”中与索结构表面上该任意点相邻的点的实测温度线性插值得到时,线性插值得到的索结构表面上该任意点的温度与索结构表面上该任意点的实际温度的误差不大于5% ;索结构表面包括支承索表面;第二个条件是“R个索结构表面点”中在同一海拔高度的点的数量不小于4,且“R个索结构表面点”中在同一海拔高度的点沿着索结构表面均布;“R个索结构表面点”沿海拔高度的所有两两相邻索结构表面点的海拔高度之差的绝对值中的最大值Ah不大于0.2°C除以Λ Th得到的数值,为方便叙述取Λ Th的单位为。C/m,为方便叙述取Ah的单位为m ;“R个索结构表面点”沿海拔高度的两两相邻索结构表面点的定义是指只考虑海拔高度时,在“R个索结构表面点”中不存在一个索结构表面点,该索结构表面点的海拔高度数值介于两两相邻索结构表面点的海拔高度数值之间;第三个条件是查询或按气象学常规计算得到索结构所在地和所在海拔区间的日照规律,再根据索结构的几何特征及方位数据,在索结构上找到全年 受日照时间最充分的那些表面点的位置,“R个索结构表面点”中至少有一个索结构表面点是索结构上全年受日照时间最充分的那些表面点中的一个点; c.按照“本方法的索结构的温度测量计算方法”直接测量计算得到初始状态下的索结构稳态温度数据,初始状态下的索结构稳态温度数据称为初始索结构稳态温度数据,记为“初始索结构稳态温度数据向量T。” ;实测或查资料得到索结构所使用的各种材料的随温度变化的物理和力学性能参数;在实测得到初始索结构稳态温度数据向量T。的同一时刻,直接测量计算得到所有支承索的初始索力,组成初始索力向量F。;依据包括索结构设计数据、竣工数据在内的数据得到所有支承索在自由状态即索力为O时的长度、在自由状态时的横截面面积和在自由状态时的单位长度的重量,以及获得这三种数据时所有支承索的温度,在此基础上利用所有支承索的随温度变化的物理性能参数和力学性能参数,按照常规物理计算得到所有支承索在初始索结构稳态温度数据向量T。条件下的索力为O时所有支承索的长度、索力为O时所有支承索的横截面面积以及索力为O时所有支承索的单位长度的重量,依次组成支承索的初始自由长度向量、初始自由横截面面积向量和初始自由单位长度的重量向量,支承索的初始自由长度向量、初始自由横截面面积向量和初始自由单位长度的重量向量的元素的编号规则与初始索力向量F。的元素的编号规则相同;在实测得到T。的同时,也就是在获得初始索结构稳态温度数据向量T。的时刻的同一时刻,直接测量计算得到初始索结构的实测数据,初始索结构的实测数据是包括索结构集中载荷测量数据、索结构分布载荷测量数据、索结构体积载荷测量数据、所有被监测量的初始数值、所有支承索的初始索力数据、初始索结构模态数据、初始索结构应变数据、初始索结构几何数据、初始索结构支座广义坐标数据、初始索结构角度数据、初始索结构空间坐标数据在内的实测数据,在得到初始索结构的实测数据的同时,测量计算得到包括支承索的无损检测数据在内的能够表达支承索的健康状态的数据,此时的能够表达支承索的健康状态的数据称为支承索初始健康状态数据;所有被监测量的初始数值组成被监测量初始数值向量C。,被监测量初始数值向量C。的编号规则与M个被监测量的编号规则相同;利用支承索初始健康状态数据以及索结构载荷测量数据建立被评估对象初始损伤向量d。,向量d。表示用初始力学计算基准模型A。表示的索结构的被评估对象的初始健康状态;被评估对象初始损伤向量d。的元素个数等于N,d。的元素与被评估对象是--对应关系,向量d。的元素的编号规则与被评估对象的编号规则相同;如果d。的某一个元素对应的被评估对象是索系统中的一根支承索,那么d。的该元素的数值代表对应支承索的初始损伤程度,若该元素的数值为O,表示该元素所对应的支承索是完好的,没有损伤的,若其数值为100%,则表示该元素所对应的支承索已经完全丧失承载能力,若其数值介于O和100%之间,则表示该支承索丧失了相应比例的承载能力;如果d。的某一个元素对应的被评估对象是某一个载荷,本方法中取d。的该元素数值为0,代表这个载荷的变化的初始数值为O ;如果没有支承索的无损检测数据及其他能够表达支承索的健康状态的数据时,或者可以认为结构初始状态为无损伤无松弛状态时,向量d。中与支承索相关的各元素数值取O ; d.根据索结构的设计图、竣工图和初始索结构的实测数据、支承索初始健康状态数据、索结构集中载荷测量数据、索结·构分布载荷测量数据、索结构体积载荷测量数据、索结构所使用的各种材料的随温度变化的物理和力学性能参数、初始索结构稳态温度数据向量T。和和前面步骤得到的所有的索结构数据,建立计入“索结构稳态温度数据”的索结构的初始力学计算基准模型A。,基于A。计算得到的索结构计算数据必须非常接近其实测数据,其间的差异不得大于5% ;对应于A。的“索结构稳态温度数据”就是“初始索结构稳态温度数据向量T。” ;对应于A。的被评估对象健康状态用被评估对象初始损伤向量d。表示;对应于A。的所有被监测量的初始数值用被监测量初始数值向量C。表示;T。和d。是A。的参数,由A。的力学计算结果得到的所有被监测量的初始数值与C。表示的所有被监测量的初始数值相同,因此也可以说C。由A。的力学计算结果组成,在本方法中A。、C。、d。和T。是不变的; e.在本方法中,字母i除了明显地表示步骤编号的地方外,字母i仅表示循环次数,即第i次循环;第i次循环开始时需要建立的或已建立的索结构的当前初始力学计算基准模型记为当前初始力学计算基准模型A^A。和Ai。计入了温度参数,可以计算温度变化对索结构的力学性能影响;第i次循环开始时,对应于Aitj的“索结构稳态温度数据”用当前初始索结构稳态温度数据向量Titj表示,向量Titj的定义方式与向量T。的定义方式相同,Titj的元素与T。的元素一一对应;第i次循环开始时需要的被评估对象当前初始损伤向量记为(^。,屮。表示该次循环开始时索结构Ai0的被评估对象的健康状态,Cli0的定义方式与d。的定义方式相同,dt的元素与d。的元素一一对应;第i次循环开始时,所有被监测量的初始值,用被监测量当前初始数值向量C1。表示,向量C1。的定义方式与向量C。的定义方式相同,C1。的元素与C。的元素一一对应,被监测量当前初始数值向量Ci。表示对应于Aitj的所有被监测量的具体数值;1"。和f。是Aitj的特性参数,C1。由Ai0的力学计算结果组成;第一次循环开始时,Aitj记为A1。,建立A1。的方法为使A1。等于A。;第一次循环开始时,Titj记为T1。,建立T1。的方法为使T1。等于T。;第一次循环开始时,Cli0记为d1。,建立d1。的方法为使d1。等于d。;第一次循环开始时,Ci0记为C10,建立C1。的方法为使C1。等于C0 ; f.从这里进入由第f步到第q步的循环;在结构服役过程中,按照“本方法的索结构的温度测量计算方法”不断实测计算获得索结构稳态温度数据的当前数据,所有“索结构稳态温度数据”的当前数据组成当前索结构稳态温度数据向量Ti,向量Ti的定义方式与向量T0的定义方式相同,Ti的元素与T。的元素一一对应;在实测得到向量Ti的同时,实测得到在获得当前索结构稳态温度数据向量Ti的时刻的同一时刻的索结构中所有被监测量的当前值,所有这些数值组成被监测量当前数值向量Ci,向量Ci的定义方式与向量C。的定义方式相同,Ci的元素与C。的元素一一对应,表示相同被监测量在不同时刻的数值;在实测得到当前索结构稳态温度数据向量Ti的同时,对新增加的M2根传感索进行无损检测,从中鉴别出出现损伤或松弛的传感索,依据被监测量编号规则,从本方法之前出现的按照被监测量编号规则编号的各向量中去除与鉴别出的出现损伤或松弛的传感索对应的元素,在本方法之后出现的各向量和矩阵中也不再出现与鉴别出的出现损伤或松弛的传感索对应的元素,在本方法之后提到传感索时不再包括这里被鉴别出出现损伤或松弛的传感索,在本方法之后提到被监测量时不再包括这里被鉴别出出现损伤或松弛的传感索的索力;从索结构上鉴别出几根出现损伤或松弛的传感索,就将M2和M减小同样的数量;在实测得到当前索结构稳态温度数据向量Ti的同一时刻,实测得到索结构中所有M1根支承索的索力数据,所有这些索力数据组成当前索力向量Fi,向量Fi的元素与向量F。的元素的编号规则相同;在实测得到当前索结构稳态温度数据向量Ti的同一时刻,实测计算得到所有M1根支承索的两个支承端点的空间坐标,两个支承端点的空间坐标在水平方向分量的差就是两个支承端点水平距离,所有支承索的两个支承端点水平距离数据组成当前支承索两支承端点水平距离向量,当前支承索两支承端点水平距离向量的元素的编号规则与初始索力向量F。的元素的编号规则相同;` g.根据当前索结构稳态温度数据向量T1,按照步骤gl至g3更新当前初始力学计算基准模型A1。、被监测量当前初始数值向量C1。和当前初始索结构稳态温度数据向量Τ1。,而被评估对象当前初始损伤向量f。保持不变; gl.比较Ti与f。,如果Ti等于f。,则A1。、Ci0和Titj保持不变;否则需要按下列步骤对Ai0和Ti0进行更新; g2.计算Ti与T。的差,Ti与T。的差就是当前索结构稳态温度数据关于初始索结构稳态温度数据的变化,Ti与T。的差用稳态温度变化向量S表示,S等于Ti减去T。,S表示索结构稳态温度数据的变化; g3.对A。中的索结构施加温度变化,施加的温度变化的数值就取自稳态温度变化向量S,对A。中的索结构施加的温度变化后得到更新的当前初始力学计算基准模型A1。,更新A1。的同时,Ti。所有元素数值也用Ti的所有元素数值对应代替,即更新了 Ti。,这样就得到了正确地对应于Aitj的ή ;此时(保持不变;当更新Aitj后,Aitj的索的健康状况用被评估对象当前初始损伤向量(表示,Aitj的索结构稳态温度用当前索结构稳态温度数据向量Ti表示,更新C1。的方法是:当更新(后,通过力学计算得到Aitj中所有被监测量的、当前的具体数值,这些具体数值组成C、; h.在当前初始力学计算基准模型Ai。的基础上,按照步骤hi至步骤h4进行若干次力学计算,通过计算建立单位损伤被监测量数值变化矩阵Λ Ci和被评估对象单位变化向量Diu ; h1.在第i次循环开始时,直接按步骤h2至步骤h4所列方法获得Λ Ci和Diu ;在其它时刻,当在步骤g中对Ai0进行更新后,必须按步骤h2至步骤h4所列方法重新获得Λ Ci和Diu,如果在步骤g中没有对Aitj进行更新,则在此处直接转入步骤i进行后续工作; h2.在当前初始力学计算基准模型Aitj的基础上进行若干次力学计算,计算次数数值上等于所有被评估对象的数量N,有N个评估对象就有N次计算;依据被评估对象的编号规贝1J,依次进行计算;每一次计算假设只有一个被评估对象在原有损伤或载荷的基础上再增加单位损伤或载荷单位变化,具体的,如果该被评估对象是索系统中的一根支承索,那么就假设该支承索再增加单位损伤,如果该被评估对象是一个载荷,就假设该载荷再增加载荷单位变化,用Diuk记录这一增加的单位损伤或载荷单位变化,其中k表示增加单位损伤或载荷单位变化的被评估对象的编号,Diuk是被评估对象单位变化向量Diu的一个元素,被评估对象单位变化向量 Diu的元素的编号规则与向量d。的元素的编号规则相同;每一次计算中再增加单位损伤或载荷单位变化的被评估对象不同于其它次计算中再增加单位损伤或载荷单位变化的被评估对象,每一次计算都利用力学方法计算索结构的所有被监测量的当前计算值,每一次计算得到的所有被监测量的当前计算值组成一个被监测量计算当前向量;当假设第k个被评估对象再增加单位损伤或载荷单位变化时,用Cidt表示对应的“被监测量计算当前向量”;在本步骤中给各向量的元素编号时,应同本方法中其它向量使用同一编号规则,以保证本步骤中各向量中的任意一个元素,同其它向量中的、编号相同的元素,表达了同一被监测量或同一对象的相关信息K1d5的定义方式与向量C。的定义方式相同,C1ck的元素与C。的元素一一对应; h3.每一次计算得到的向量Cidt减去向量C1。得到一个向量,再将该向量的每一个元素都除以本次计算所假设的单位损伤或载荷单位变化数值后得到一个“被监测量的数值变化向量δ CY ;有N个被评估对象就有N个“被监测量的数值变化向量”; h4.由这N个“被监测量的数值变化向量”按照N个被评估对象的编号规则,依次组成有N列的“单位损伤被监测量数值变化矩阵△ Ci”;单位损伤被监测量数值变化矩阵八^勺每一列对应于一个被监测量单位变化向量;单位损伤被监测量数值变化矩阵ACi的每一行对应于同一个被监测量在不同被评估对象增加单位损伤或载荷单位变化时的不同的单位变化幅度;单位损伤被监测量数值变化矩阵ACi的列的编号规则与向量d。的元素的编号规则相同,单位损伤被监测量数值变化矩阵ACi的行的编号规则与M个被监测量的编号规则相同;i.定义当前名义损伤向量(和当前实际损伤向量dSd1。和Cli的元素个数等于被评估对象的数量,Clic和Cli的元素和被评估对象之间是--对应关系,f。的元素数值代表对应被评估对象的名义损伤程度或名义载荷变化量,Clic和Cli与被评估对象初始损伤向量d。的元素编号规则相同,Clic的元素、Cli的元素与d。的元素是一一对应关系; j.依据被监测量当前数值向量Ci同“被监测量当前初始数值向量C1。”、“单位损伤被监测量数值变化矩阵AC1”和“当前名义损伤向量屮。”间存在的近似线性关系,该近似线性关系可表达为式1,式I中除f。外的其它量均为已知,求解式I就可以算出当前名义损伤向量札;
Figure CN103852282AC00101
k.利用式2表达的当前实际损伤向量Cli的第k个元素Clik同被评估对象当前初始损伤向量必的第k个元素d、和当前名义损伤向量的第k个元素Clidt间的关系,计算得到当前实际损伤向量Cli的所有元素;
Figure CN103852282AC00102
式2中k=l,2,3,……,N ;向量Cli的元素的编号规则与式I中向量d。的元素的编号规则相同表示第i次循环中第k个被评估对象的当前实际健康状态,如果该被评估对象是索系统中的一根支承索,那么Clik表示其当前健康问题的严重程度,有健康问题的支承索可能是松弛索、也可能是受损索,Clik数值反应了该支承索的松弛或损伤的程度;如果该被评估对象是一个载荷,那么Clik表示该载荷的实际变化量;将被评估对象当前实际损伤向量Cli中与M1根支承索相关的M1个元素取出,组成支承索当前实际损伤向量Cfi,支承索当前实际损伤向量Cfi的元素的编号规则与初始索力向量F。的元素的编号规则相同;支承索当前实际损伤向量Cfi的第h个元素表示索结构中第h根支承索的当前实际损伤量,h=l, 2,3,…….,M1 ;支承索当前实际损伤向量Cfi中数值不为O的元素对应于有健康问题的支承索,从这些有健康问题的支承索中鉴别出受损索,剩下的就是松弛索;与受损索对应的支承索当前实际损伤向量Cfi中的元素的数值表达的是该受损索的当前实际损伤,元素数值为100%时表示该支承索彻底丧失承载能力,介于O与100%之间时表示丧失相应比例的承载能力;利用在当前索结构稳态温度数据向量Ti条件下的、在第I步鉴别出的松弛索及用支承索当前实际损伤向量Cfi表达的这些松弛索的、与其松弛程度力学等效的当前实际等效损伤程度,利用在第f步获得的在当前索结构稳态温度数据向量Ti条件下的当前索力向量Fi和当前支承索两支承端点水平距离向量,利用在第c步获得的在初始索结构稳态温度数据向量T。条件下的支承索的初始自由长度向量、初始自由横截面面积向量和初始自由单位长度的重量向量、初始索力向量F。,利用当前索结构稳态温度数据向量Ti表示的支承索当前稳态温度数据,利用在第c步获得的在初始索结构稳态温度数据向量T。表示的支承索初始稳态温度数据,利用在第c步获得的索结构所使用的各种材料的随温度变化的物理和力学性能参数,计入温度变化对支承索物理、力学和几何参数的影响,通过将松弛索同受损索进行力学等效来计算松弛索的、与当前实际等效损伤程度等效的松弛程度,力学等效条件是:一、两等效的索的无松弛和无损伤时的初始自由长度、几何特性参数、密度及材料的力学特性参数相同;二、松弛或损伤后,两等效的松弛索和损伤索的索力和变形后的总长相同;满足上述两个力学等效条件时,这样的两根支承索在索结构中的力学功能就是完全相同的,即如果用等效的松弛索代替受损索后,索结构不会发生任何变化,反之亦然;依据前述力学等效条件求得那些被判定为松弛索的松弛程度,松弛程度就是支承索自由长度的改变量,也就是确定了那些需调整索力的支承索的索长调整量;这样就实现了支承索的松弛识别和损伤识别;计算时所需索力由当前索力向量Fi对应元素给出;本方法将受损索和松弛索统称为有健康问题的支承索,简称为问题索,至此本方法实现了剔除载荷变化和结构温度变化的影响的、索结构的问题索识别,同时实现了剔除结构温度变化和支承索健康状态变化影响的、载荷变化量的识别;l.在求得当前名义损伤向量t后,按照式3建立标识向量Bi,式4给出了标识向量Bi的第k个元素的定义;
Figure CN103852282AC00111
式4中元素Bik是标识向量Bi的第k个元素,Diuk是被评估对象单位变化向量Diu的第k个元素,Click是被评估对象当前名义损伤向量f。的第k个元素,它们都表示第k个被评估对象的相关信息,式4中k=l,2,3,……,N ; m.如果标识向量B1的元素全为O,则回到步骤f继续本次循环;如果标识向量B1的元素不全为O,则进入下一步、即步骤η ; η.根据式5计算得到下一次、即第i+1次循环所需的被评估对象当前初始损伤向量di+1。的每一个元素;
Figure CN103852282AC00112
式5中(T1。,是下一次、即第i+Ι次循环所需的被评估对象当前初始损伤向量di+1。的第k个元素,d、是本次、即第i次循环的被评估对象当前初始损伤向量必的第k个元素,Diuk是第i次循环的被评估对象单位变化向量Diu的第k个元素,Bik是第i次循环的标识向量Bi的第k个元素,式5中k=l,2,3,……,N ; ο.取下一次、即第i+Ι次循环所需的当前初始索结构稳态温度数据向量Ti+1。等于第i次循环的当前初始索结构稳态温度数据向量Titj ; P.在初始力学计算基准模型A。的基础上,对A。中的索结构施加温度变化,施加的温度变化的数值就取自稳态温度变化向量S,再令索的健康状况为di+1。后得到的就是下一次、SP第i+Ι次循环所需的力学计算基准模型Ai+1 ;得到Ai+1后,通过力学计算得到Ai+1中所有被监测量的、当前的具体数值,这些具体数值组成下一次、即第i+Ι次循环所需的被监测量当前初始数值向量卜+1。; q.回到步骤f,开始下一次循环。
CN201410085128.2A 2014-03-10 2014-03-10 索力监测问题索载荷递进式识别方法 Pending CN103852282A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410085128.2A CN103852282A (zh) 2014-03-10 2014-03-10 索力监测问题索载荷递进式识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410085128.2A CN103852282A (zh) 2014-03-10 2014-03-10 索力监测问题索载荷递进式识别方法

Publications (1)

Publication Number Publication Date
CN103852282A true CN103852282A (zh) 2014-06-11

Family

ID=50860200

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410085128.2A Pending CN103852282A (zh) 2014-03-10 2014-03-10 索力监测问题索载荷递进式识别方法

Country Status (1)

Country Link
CN (1) CN103852282A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104990588A (zh) * 2015-07-23 2015-10-21 东南大学 精简索力监测受损索载荷递进式识别方法
CN104990737A (zh) * 2015-07-23 2015-10-21 东南大学 精简索力监测问题索载荷递进式识别方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102768123A (zh) * 2012-05-30 2012-11-07 东南大学 温度变化索力监测的问题索和支座角位移递进式识别方法
CN103604631A (zh) * 2013-12-09 2014-02-26 东南大学 线位移时索力监测的问题索集中载荷递进式识别方法
CN103616230A (zh) * 2013-12-09 2014-03-05 东南大学 索力监测的受损索和集中载荷递进式识别方法
CN103616197A (zh) * 2013-12-09 2014-03-05 东南大学 混合监测的问题索和集中载荷递进式识别方法
CN103616114A (zh) * 2013-12-09 2014-03-05 东南大学 角位移时索力监测的问题索集中载荷递进式识别方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102768123A (zh) * 2012-05-30 2012-11-07 东南大学 温度变化索力监测的问题索和支座角位移递进式识别方法
CN103604631A (zh) * 2013-12-09 2014-02-26 东南大学 线位移时索力监测的问题索集中载荷递进式识别方法
CN103616230A (zh) * 2013-12-09 2014-03-05 东南大学 索力监测的受损索和集中载荷递进式识别方法
CN103616197A (zh) * 2013-12-09 2014-03-05 东南大学 混合监测的问题索和集中载荷递进式识别方法
CN103616114A (zh) * 2013-12-09 2014-03-05 东南大学 角位移时索力监测的问题索集中载荷递进式识别方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104990588A (zh) * 2015-07-23 2015-10-21 东南大学 精简索力监测受损索载荷递进式识别方法
CN104990737A (zh) * 2015-07-23 2015-10-21 东南大学 精简索力监测问题索载荷递进式识别方法

Similar Documents

Publication Publication Date Title
CN103913342A (zh) 角度监测问题索载荷广义位移递进式识别方法
CN103852304A (zh) 索力监测问题索载荷线位移识别方法
CN103852309A (zh) 混合监测问题索载荷线位移递进式识别方法
CN103868745A (zh) 索力监测受损索载荷广义位移递进式识别方法
CN103852291A (zh) 广义位移空间坐标监测问题索载荷递进式识别方法
CN103852317A (zh) 角位移混合监测问题索载荷递进式识别方法
CN103868732A (zh) 线位移应变监测问题索载荷递进式识别方法
CN103913327A (zh) 线位移空间坐标监测问题索载荷递进式识别方法
CN103868727A (zh) 线位移角度监测问题索载荷递进式识别方法
CN103852312A (zh) 角位移角度监测问题索载荷递进式识别方法
CN103868743A (zh) 广义位移应变监测问题索载荷递进式识别方法
CN103913323A (zh) 线位移索力监测问题索载荷递进式识别方法
CN103852319A (zh) 索力监测问题索载荷线位移递进式识别方法
CN103852288A (zh) 索力监测受损索载荷递进式识别方法
CN103868728A (zh) 角位移索力监测问题索载荷递进式识别方法
CN103868707A (zh) 应变监测问题索载荷递进式识别方法
CN103852321A (zh) 角度监测问题索载荷递进式识别方法
CN103868715A (zh) 索力监测问题索载荷角位移递进式识别方法
CN103868738A (zh) 广义位移索力监测受损索载荷递进式识别方法
CN103868744A (zh) 索力监测问题索载荷广义位移递进式识别方法
CN103913336A (zh) 索力监测受损索载荷角位移递进式识别方法
CN103868742A (zh) 角位移应变监测问题索载荷递进式识别方法
CN103868722A (zh) 角位移空间坐标监测问题索载荷递进式识别方法
CN103852282A (zh) 索力监测问题索载荷递进式识别方法
CN103868736A (zh) 线位移索力监测受损索载荷递进式识别方法

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20140611

RJ01 Rejection of invention patent application after publication