桥梁混凝土结构健康诊断的理论和方法-代东辉

1 第一章第一章健康诊断体系健康诊断体系 2 1 1 基于原位监测的健康诊断体系基于原位监测的健康诊断体系 2 第二章第二章监控模型监控模型 3 2 1 监控模型监控模型 3 2 2 因果分析模型因果分析模型 基于积分概念的健康诊断的因果分析模型基于积分概念的健康诊断的因果分析模型 3 2 2 1车辆荷载 3 2 2 2温度荷载 5 2 2 3时效分量 6 2 2 4效应量 监控模型 6 2 3 因果分析模型因果分析模型 基于遗传模拟退火算法的健康诊断的因果分析模型基于遗传模拟退火算法的健康诊断的因果分析模型 7 2 3 1遗传模拟退火算法的实现 7 2 4 因果分析模型因果分析模型 基于偏最小二乘回归法的健康诊断的因果分析模型基于偏最小二乘回归法的健康诊断的因果分析模型 8 2 5 因果分析模型因果分析模型 桥梁健康诊断的组合因果分析模型桥梁健康诊断的组合因果分析模型 8 第三章第三章诊断标准诊断标准 9 3 1 监控模型的诊断标准监控模型的诊断标准 9 3 2 结构健康控制指标的诊断标准结构健康控制指标的诊断标准 10 3 2 1力学解释 10 3 2 2变形控制指标的诊断标准 11 3 2 3应力应变控制指标的诊断标准 12 3 3 时效分量的诊断标准时效分量的诊断标准 12 3 4 桥梁结构整体健康等级及标准桥梁结构整体健康等级及标准 13 3 4 1整体健康等级 13 3 4 2整体健康等级的标准 14 2 第一章第一章 健康诊断体系健康诊断体系 1 1 基于原位监测的健康诊断体系基于原位监测的健康诊断体系 原位监测资料及其分析成果是重大桥梁混凝土结构健康诊断的定量依据 构造的基于原位监测 的健康诊断体系如下图所示 原位监测的健康诊断体系 基于原位监测的健康诊断体系包括七层 第一层为需诊断的桥梁结构的健康状态 第二层为反 映结构健康的变形 应力等监测类别 第三层为第二层的效应量 第四层为第三层的监测方法 第 五层为监测方法的测点 第六层为依据测点的测量值 并依据集成和融合的原则 建立桥梁结构健 康诊断的因果分析模型 健康标准和诊断模型 第七层为依据第六层的诊断结果 建立桥梁结构健 康综合诊断 3 第二章第二章 监控模型监控模型 2 1 监控模型监控模型 采用监控模型对桥梁结构健康进行定量分析与评价 即监测效应量一般分解为车辆荷载分量 温度荷载分量以及实效分量 即 式 1 VTg YYYYY 式中 Y 为监测的效应量 如变形 应力 变 为车辆荷载分量 即车辆荷载作用对效应 V Y 量的影响量 为温度分量 即温度荷载作用对效应量的影响量 为时效分量 即因徐变 蠕 T Y Y 变 或结构形态 或监测系统 变化等引起的影响量 为其他分量 如结构裂缝引起结构变形 g Y 的影响量等 2 2 因果分析模型因果分析模型 基于积分概念的健康诊断的因果分析模型基于积分概念的健康诊断的因果分析模型 2 2 1车辆荷载车辆荷载 由于作用在桥梁结构上的车辆荷载大小不确定 位置不确定 所以无法用准确的公式描述车辆 荷载分量引起的效应量值 这里将车辆荷载从监测数据的效应量中滤掉 那么剩下的效应量值就是 由温度分量和实效分量组成的 滤掉车辆荷载是从统计学的角度进行的 下面以效应量是挠度值为 例进行说明 监测到的挠度值是包括车辆荷载分量 温度分量 实效分量三者组成的 监测数据内 容为挠度值 传感器采集相应挠度值时对应的时间 传感器采集相应挠度值时对应的温度 下面一 一个传感器采集到的数据为例画出挠度随时间变化的曲线 如下图所示 4 050100150200250300350400450500 20 0 20 40 60 80 100 如果没有车辆荷载的影响 由于温度的因素 桥梁挠度曲线在一年当中应该呈现周期性的变化 大致应呈正弦曲线 上图可能由于采集挠度值的传感器由于一些因素引起位置的变化 所以在第 100 天和第 200 天左右时挠度值产生阶跃 分析要从 200 天之后的数据进行 因为监测的挠度值都 是以一个基准值为标准的 传感器位置的变化就会引起基准值的改变 基准值一旦改变就会使监测 到的挠度不准确 为方便观察取出上面图形当中的一个时间段 蓝线表示总的挠度值与时间的关系 曲线 从图形可以看出挠度曲线是由无数个小的突变图形组成的 这里认为这些小的突变就对应车 辆荷载引起的挠度效应量 滤掉车辆荷载采用的方法是用 matlab 软件编写程序取平均值使挠度曲 线光滑 例如第 i 个数据点的挠度值 取其前面一点的挠度值和后面一点的平均值 进行多次这样 的循环操作得到的值作为该点滤掉车辆荷载的挠度值 经过多次这样的操作 挠度曲线就去除了突 变值变得光滑 如下图所示 蓝线表示监测到的总的挠度值与时间的关系曲线 红线表示滤掉车辆 荷载后的挠度值与时间之间的关系曲线 5 299 8300300 2300 4300 6300 8301 87 88 89 90 91 92 93 94 95 2 2 2温度荷载温度荷载 温度荷载变化规律较复杂 对桥梁结构的影响结果又往往与其他荷载共同作用 温度荷载在一 天 一年中作用的周期性 对桥梁结构的变形 应力产生不同的影响结果 所以用简单的线性表达 式很难准确的表示温度对桥梁结构的影响 所以使用积分回归的的概念对温度分量进行修正 i i aT 从而建立优化模型 应用积分概念反映温度对变形的影响为 式 4 0 Tii Aa t T t dt 式中 为位移温度荷载分量 A 为积分常数 为温度因子的系数 表示某一时刻 T i a t 温度因子每变化一个单位时对位移的影响 非常数 tt 应用 式 3 对统计模型中的温度分量的因子选择进行改进 将 式 3 处理为与某一温度因 素 T t 和某一时间函数 f t 的二元函数 F f T 将此函数按二元多项式展开 则考虑不同 时间温度的影响因子形式为 式 5 11 lm ij Tij ij a ft T 式中 为位移温度荷载分量 为时间函数 为回归系数 l m 为各影响因子个数 T f t ij a 6 根据温度的周期性 建立温度的如下因子 1 将设为周期函数 用以反映温度因子在一天 年 内对变形和应力的不同影响 f t 2 将设为时间函数 用以反映温度因子对变形和应力影响的长期变化 f t 综上所述 可表示为 T 式 6 其实这个式子也是由拟合确定 8760 也是选 j m j ijijT T it b it a 1 8760 2 cos 8760 2 sin 择较好的拟合模型所得的数值 即 24365 而且通过拟合发现的最高次项也只能是 1 T 式中 为积分回归系数 t 为观测日至始测日的累积天数 T 为观测日相对 ij a ij b ij c ij d 于始测日的变温 S 为季节系数 为各影响因子个数 m 1 l 2 l 3 l 2 2 3时效分量时效分量 定义时效分量表达式为 式 7 12ln cc 式中 为回归系数 为时效因子 取值待定 1 c 2 c 2 2 4效应量效应量 监控模型监控模型 综上效应量的温度分量和实效分量优化监控模型表达式为 式 8 ln 8760 2 cos 8760 2 sin 21 1 ccT it b it a j m j ijijT 用 matlab 软件程序建立效应量的温度分量与实效分量之和的拟合曲线 如下图所示 其中蓝 线代表监测到的挠度值与时间的关系 红线代表滤掉车辆荷载后的挠度值与时间的关系 绿线代表 根据拟合公式画出的滤掉车辆荷载后的挠度值与时间的关系 7 0100200300400500600 20 0 20 40 60 80 100 120 显然绿线与红线偏离较大 说明曲线拟合的效果不好 需要重新选定拟合公式对滤掉车辆荷载 后的挠度曲线进行重新拟合 拟合公式正在改进当中 直至所使用的拟合公式可以使得绿线与红线 较为接近时方可 经过多次改变拟合公式 效果都不是很理想 拟合公式有待进一步尝试 2 3 因果分析模型因果分析模型 基于遗传模拟退火算法的健康诊断的因果分析模型基于遗传模拟退火算法的健康诊断的因果分析模型 遗传算法和模拟退火算法是模拟生物遗传和进化过程中的选择 交叉 变异机理而建立起的一 种自适应全局优化概率搜索方法 它能够以随机搜索技术从概率意义上找出目标函数的全局最优点 遗传算法把握总体的能力较强 但局部搜索能力较差 模拟退火算法具有较强的局部搜索能力 能 避免搜索过程陷入局部最优 遗传模拟退火算法即是将遗传算法与模拟退火算法相结合 互相取长 补短 2 3 1遗传模拟退火算法的实现遗传模拟退火算法的实现 应用遗传模拟退火算法 建立因果分析模型包括一下几个因素 1 参数编码 遗传算法不对所求解的实际决策变量直接进行操作 因此需要把一个问题的可行解从解空间转换 到遗传算法所能处理的搜索空间 二进制编码是遗传算法中最常用的一种编码方法对于桥梁健康检 测和监测中的建模的变量筛选问题 用二进制编码也是较合理的方法 2 群体初始化 8 初始群体中各个个体的基因值可以随机或根据一定的限制条件生成 群体的大小即群体中所 有个体的数量一般为 20 100 对于桥梁健康监测中群体数的具体数值还有待进一步确定 3 个体适应度函数的设计 遗传算法中使用适应度来度量群体中各个个体在优化计算中有可能达到或接近于或有助于找 到最优解的优良程度 并且按与个体适应度成正比的概率来决定当前群体中每个个体遗传到下一代 群体中的机会多少 4 遗传算子的设计 1 选择算子 其作用是从当前群体中选择一些比较优良的个体 并将其复制到下一代群体中 2 交叉算子 其作用是产生一些较好的新个体模式 寻优的搜索过程主要通过它来实现 3 变异算子 其目的是改善遗传的局部搜索能力 维持群体的多样性 防止出现早熟现象 5 遗传模拟退火算法的设计 首先通过选择 交叉 变异等遗传操作产生一组新个体 然后独立地随机选择各个个体中的 两个基因作为饶动后的个体所得的复相关系数 R 值增大则一定接受这个新个体 而新个体所得的复 相关 R 值减小时 一某一概率接受这个新个体 6 算法的终止条件 遗传算法的判据一般都是启发式的 通常以进化代数 计算时间 解的质量等确定判据 2 4 因果分析模型因果分析模型 基于偏最小二乘回归法的健康诊断的因果分析模型基于偏最小二乘回归法的健康诊断的因果分析模型 1 基本原理 偏最小二乘回归法是一种新颖的多元分析方法 该方法融合多元线性回归分析 典型相关分 析和主成份分析的基本功能 将建模预测类型的数据分析方法与非模型式的数据认识性分析有机地 结合起来 2 5 因果分析模型因果分析模型 桥梁健康诊断的组合因果分析模型桥梁健康诊断的组合因果分析模型 桥梁工作条件及其复杂 用单一的数学模型进行诊断和分析还存在一定的局限性 为此 将组 合模型引入桥梁健康诊断中 基本思路是 选择几个 较优 的监控模型 尤其是检测发现隐患病 害处的相应监控模型 然后采用某种方式将它们组合起来 以得到更好的监控模型 一下讨论线性 和非线性两种组合方法 1 线性组合模型 线性组合模型也被称为最有权组合模型 其基本思路是 先优选几个监控模型 然后选取一组最优 9 的权重 加大隐患病害出相应监控模型的权 将它们线性组合起来 从而得到组合模型 而建立组 合模型的关键就是合理地确定权系数向量 2 非线性组合模型 采用小波神经网络来实现多模型的非线性组合 以上的几种因果分析模型对于第一种虽然公式没有拟合的很理想 但思路还是比较清晰的 后几种 模型的建立有一定的难度 看了一些资料较难理解 第三章第三章 诊断标准诊断标准 3 1 监控模型的诊断标准监控模型的诊断标准 由 式 1 结合实测资料 用回归分析法和范数最小原理建立安全监控模型 即 式 9 YYYTT 式中 分别为滤掉车辆荷载后模型的计算值 温度荷载计算值 时效分量计 TY T Y Y 算值和其他分量计算值 某一时刻 将相应的温度荷载 F 时间等参数输入 式 9 中 即可求得滤掉车辆荷载后模 型的计算值 则可用以下公式判别反映桥梁结构健康因素的某一效应量 Y 是否正常 TY 当 式 10 SYY T T2 则为正常 当 式 11 3 2 SSYY T T 则跟踪监测 2 3 次 当无趋势性变化或检测未发现该部位有隐患病害 则为正常 若有趋势 性变化或检测发现该部位有隐患病害 则为异常 10 当 式 12 SYY T T3 则为异常 在 式 10 式 11 式 12 中 效应量的实测值滤掉车辆荷载所得值 为模型 T Y TY 计算值 S 为模型的标准差 用健康模型的上述标准 仅能判断某一检测效应量的测量值是否正常 还难以对结构的健康做 出判断 因此 尚需研究用健康诊断的控制值 实效分量等综合诊断桥梁混凝土结构的健康状况 3 2 结构健康控制指标的诊断标准结构健康控制指标的诊断标准 控制指标对桥梁结构的健康诊断具有特殊的意义 对监测桥梁的安全运行起着重要的作用 此 外 依据控制指标可及时 简单地判断桥梁结构的健康状况 控制指标主要包括变形和应力 变 而变形是监控和诊断桥梁结构健康的最主要监测量 其 是桥梁结构在各种荷载组合作用的综合反映 但其影响因素也十分复杂 3 2 1力学解释力学解释 变形控制指标拟定的方法 及应用控制指标诊断桥梁结构健康的标准 主要是依据材料的力学 性能 利用小变形的粘弹性 粘塑性理论以及大变形弹塑性理论 提出拟定多级变形的控制指标 桥梁混凝土结构的破坏主要经历四个阶段 1 线弹性或粘弹性阶段 在荷载作用下 混凝土中任一部位的应力均小于材料强度的比例界 限 此时 结构总体上处于线弹性工作阶段 应力应变呈线性关系 荷载与位移关系基本呈线性 此阶段对应的控制指标为一级 2 准线弹性阶段 随荷载的增加和结构参数的变化 桥梁结构受拉区出现裂缝 拉应力转由 受拉区钢筋承受 虽然结构出现裂缝 但结构变形与荷载仍呈线性关系 此阶段对应控制指标为二 级 3 屈服变形阶段 随荷载的继续增加或结构参数的继续变化 结构出现拉应力 尤其在梁底 拉应力出现最为明显 受拉区钢筋进入极限状态 结构变形显著增加 与荷载呈非线性关系 此时 结构处于弹塑性工作状态 此阶段对应的控制指标为三级 4 失稳破坏阶段 随荷载进一步增加或参数进一步变化 当超过一定值后 结构变形急剧增 11 加 结构出现超过规范规定的裂缝 屈服区急剧扩展 发生大变形 当荷载达到临界值时 结构完 全破坏 从而三式继续承载的能力 此阶段对应的控制指标为四级 3 2 2变形控制指标的诊断标准变形控制指标的诊断标准 根据控制指标 诊断的控制准则如下 一级 若 式 13 01 m 桥梁结构未出现裂缝 变形正常 即为弹性工作工作阶段 则结构健康 二级 若 式 14 102 mm 桥梁结构出现裂缝 开始带裂缝工作 变形较正常 即准弹性工作阶段 则结构正常 三级 若 式 15 203 mm 桥梁结构变形异常 裂缝有开展趋势 即为弹塑性阶段 并检测发现有隐患病害 有使桥梁结 构恶化的趋势 结构为病态 四级 若 式 16 03 m 桥梁结构处于大变形阶段 裂缝急剧开裂 结构失事破坏 在 式 13 式 14 式 15 式 16 中 分别为变形的一级 二级 三 1 m2 m3 m 级监控指标 其值根据结构破坏的理论与方法确定 为桥梁结构的实测变形 变形控制指标的 0 拟定方法如下 采用设计荷载组合工况 由对应的允许应力和允许的抗滑稳定安全系数作为约束条件 1 m 应用粘弹性有限元分析计算 选择出现裂缝是的最不利荷载组合工况 有对应强度 稳定和抗裂等指标 采用小变形 2 m 的粘弹塑性有限元程序拟定二级变形控制指标 用大变形的粘弹塑形有限元分析计算 3 m 对于这些指标的具体拟定会是一个复杂的工作 12 3 2 3应力应变控制指标的诊断标准应力应变控制指标的诊断标准 对于桥梁结构的控制部位的应力是否在允许范围内 控制准则如下 1 当控制部位 1 4 跨 1 2 跨 3 4 跨 以及支座部位 在静荷载作用时 实测值的最大主 压应力小于等于允许压应力 实测的最大拉应力小于或等于允许拉应力 则结构处于线弹性阶段 为健康 2 当控制部位在动荷载作用时 实测的最大主压应力小于等于规范规定的允许压应力 实测的最 大拉应力小于或等于规范规定的允许拉应力 则结构处于准弹性阶段 为健康 3 当控制部位的实测值应力大于规范或设计允许值 而小于或等于混凝土极限强度是 则强 度超过规范标准 结构处于弹塑性阶段 并检测发现有隐患病害 则结构局部有病变 4 当控制部位的实测值大于极限强度 则结构出现险情 结构失事 3 3 时效分量的诊断标准时效分量的诊断标准 1 变形实效分量的诊断标准 在 式 1 中 时效分量反映了桥梁结构健康的变化趋势 Y 变形的时效位移 综合反映了桥梁结构的徐变 蠕变 和结构裂缝以及隐患等不可逆变 形 可以反映结构的趋势性变化 是衡量结构的重要指标 因此对时效分量进行进行分离 并对变 化趋势进行分析 对桥梁的健康诊断很有帮助 t a t b t c 时效分量变化规律 1 时效分量变化平稳 时效分量变化过程如上图 线所示 则有 式 17 0 d dt 13 说明桥梁结构变形的变化规律正常 2 时效分量呈收敛趋势 如图 时效分量变化过程如上图 线所示 则有 但或 式 18 0 0 2 2 dt d dt d m 0m10 式中 分别为变形的一级控制值和历史极值 为观测的变形值 为时效分量 1 m m 0 这说明时效分量减速率增长 但实测位移小于等于一级控制值或历史极值 若检测又尚未 0 发现有危机结构的隐患病害 则变化规律基本正常 3 时效分量呈线性增大趋势 如上图 b 所示 则 但 式 19 0 2 2 dt d mm201 式中 为变形的二级控制值 2 m 说明时效分量等速增长 但实测位移值小于二级控制值 则变形规律异常 0 4 时效分量加速增长 如上图 c 所示 则 但 式 20 0 2 2 dt d 0 3 3 dt d mm302 式中 为变形的二级控制值 3 m 这说明变形规律处于险情 2 神经网络分离实效分量的诊断标准 3 尖点突变模型诊断 第二种和第三种诊断模型不太明白 有待学习 14 3 4 桥梁结构整体健康等级及标准桥梁结构整体健康等级及标准 3 4 1整体健康等级整体健康等级 依据混凝土结构设计规范 监测原理 安全检查方法的规定 结合医学界对人体健康诊断的等 级等 对桥梁结构的健康状况划分为五个等级 健康 基本健康 病变 病重 病危 式 21 12345 NNNNNN 3 4 2整体健康等级的标准整体健康等级的标准 1 健康标准 1