基于Riccati传递矩阵法的桥梁结构概率损伤识别方法分析-结构工程专业论文.docxVIP

郑州大学硕士学位论文——————_————————-————●——_—_—————————_____———————————————●__-__————————————————一一 郑州大学硕士学位论文 ——————_————————-————●——_—_—————————_____———————————————●__-__————————————————一一 摘 要 桥梁是交通运输线路中的重要组成部分，桥梁的健康状态对确保道路畅通至关重 要，近年来国内外频繁出现桥塌人亡的重大事故，使得桥梁健康检测与损伤识别成为 众多学者和工程师们研究的热点问题之一。由于桥梁健康检测与损伤识别是一个相当 复杂的问题，国内外对此问题仍未很好解决。 本文在回顾结构近年来健康检测与损伤识别方法的研究现状和存在问题的基础 上，确定本文采用基于摄动Kiccati传递矩阵的概率损伤识别方法对桥梁进行损伤识别 研究，研究内容分为3大部分： 第l部分：根据桥梁结构具有一维特征的实际特点，利用一般传递矩阵法的单元 矩阵传递公式，给出了适用于桥梁结构在不同情况下的一般传递矩阵的零阶(不摄动)、 一阶、二阶摄动单元矩阵公式，再按照Riccati传递矩阵法的基本思想对该单元矩阵传 递公式重新组合，就可以得到各种情况下PJccati传递矩阵的各阶摄动单元矩阵传递公 式；同时，又由Pdccati传递矩阵法，给出摄动Kiccati传递矩阵法中的单元传递矩阵 以及各截面状态向量的零阶、一阶、二阶摄动递推关系式，根据此递推关系式和边界 条件，可以推出具有不确定参数的桥梁结构特征值的零阶、一阶、二阶摄动频率方程： lsl。．。=” §；，杪r西h，和}N+I,0I=0 ‰肚=胪Ns】Ⅳ+l’小h。+(1一颤)‰“。{e}Ⅳ。，+‰扎业{8)㈣)}，=o 采用数值方法求解摄动频率方程，可以得到特征值的零阶、一阶、二阶摄动结果。再 利用矩阵的奇异值分解理论，对末端截面的传递矩阵进行奇异值分解，推导出求解末 端截面状态向量的零阶、一阶、二阶的计算公式： {8}Ⅳ+1．o=IV忙jⅣ+l。={Vrj ‰哪=矿弦}Ⅳ+1，， 叠}。．庸=旷弦}Ⅳ卅，业 根据末端截面状态向量的计算结果，由截面状态向量递推关系式可以解出各截面状态 向量数值，由此得到特征向量的零阶、一阶、二阶摄动结果。 依据上述计算公式，在MATLAB环境下编制了计算机分析程序，以简支梁为例， 将各单元抗弯刚度和集中质量分别作为不定参数进行摄动计算，通过数值解和理论解 的比较，证明了所采用的计算公式和编制计算机分析程序的正确性。这些计算公式可 用于第2部分计算具有随机参数结构特征值和特征向量的统计特性。 第2部分：根据具有随机参数结构固有频率的Taylor展开式，给出了求解具有随 机参数结构固有频率和振型的均值和方差计算公式，采用摄动Kiccati传递矩阵法，编 NT求解结构特征值和特征向量统计特性(均值和方差)的计算机应用程序，仍以简 基于Riccati传递矩阵法的桥梁结构概率损伤识别方法研究支梁为例，将质量密度作为随机参数进行统计特性计算，通过数值解和理论解的比较， 基于Riccati传递矩阵法的桥梁结构概率损伤识别方法研究 支梁为例，将质量密度作为随机参数进行统计特性计算，通过数值解和理论解的比较， 证明了公式的正确性，该计算公式为桥梁结构概率损伤识别奠定了理论基础。 第3部分：依照确定的结构损伤识别方法，假定桥梁结构各个单元待识别的随机 参数是相互独立，且服从正态分布，提出了使用概率损伤识别方法识别桥梁结构损伤 的基本方法：将桥梁结构固有频率和振型的零阶、一阶、二阶摄动以及其统计特性作 为已知数据，代入频率和振型统计特性的计算公式，建立含有所有待识别参数均值和 方差的目标函数，以单元无损伤时参数的均值和方差为初值，利用MATLAB中优化 程序对该目标函数进行优化计算，求得待识别参数的最优解，由于每个待识别参数在 桥梁发生损伤前后均服从正态分布，且一般都满足方差相等的条件，根据数理统计中 两个正态总体均值的检验方法，查t双侧分布表，则可以计算得到各个待识别单元的 损伤概率，对于判断为有损伤的单元，由待识别参数均值的最优解与其初值的比值， 可进一步得到各损伤单元的近似损伤量。 文中利用摄动Riccati传递矩阵法，对于2种不同形式的桥梁结构(预应力混凝土 简支梁桥和预应力混凝土5跨连续梁桥)，分别进行概率损伤识别的数值仿真分析。在 仿真计算过程中，按照损伤位置和损伤量的不同，对于预应力混凝土简支梁桥，分4 种工况进行数值仿真计算；对于预应力混凝土5跨连续梁桥，分3种工况进行数值仿 真计算，得到了不同工况下的结构损伤识别结果： 1．对于有损伤桥梁，采用基于摄动Riccati传递矩阵法的结构概率损伤识别方法 进行桥梁结构损伤