本发明专利技术提供了一种基于节点曲率和小波分析的变截面梁损伤识别方法，该方法基于结构微损伤对其应力重分布影响很小的假定，首先测量在集中荷载匀速通过有损变截面梁时的位移信号，并计算基于无损单一尺度有限元模型求解变截面梁损伤前后的节点曲率矩阵实现损伤位置识别，进而基于时空多尺度模型对结构损伤程度展开精确计算，从而最终得到变截面梁的损伤状态。
【技术实现步骤摘要】
本专利技术属于土木工程领域，涉及一种基于节点曲率和小波分析的变截面梁损伤识 别方法。
技术介绍
随着社会的进步和土木工程技术的发展，大型土木工程结构，尤其是大型桥梁结 构的健康监测与安全状态评估目前已成为国内外学术界与工程界关注的热点。如何根据采 集的数据和信号反演出桥梁结构的工作状态和健康状况，建立精确的有限元模型，准确识 别出结构损伤部位及其损伤程度，并为桥梁结构的损伤及安全预后服务，是桥梁结构健康 监测、损伤诊断与安全评估最核心和最关键的问题之一。目前，国内外关于损伤识别方法的研宄主要集中在基于模型（model-based)和数 据驱动（data-driven)两方面，其中数据驱动的损伤识别方法即是基于结构健康监测系统 采集的海量时程数据，通过数据本身或相关变化后导出的结构单元非线性特征来识别损 伤，主要方法有时序模型方法、小波分析方法、Hilbert-Huang变换（HHT)、Kalmanfilter方 法等，该方法能快速识别损伤发生的时刻和位置，甚至能做到实时损伤识别，对基于在线健 康监测系统的实时预警具有较积极的作用；而基于模型的损伤识别方法主要是指利用结构 的实测信息通过不断修正有限元模型来进行结构损伤识别，该方法能从整体上认识被监测 的结构，且能有效进行预测，但其计算量和数据存储量较大，往往需要很繁琐的计算过程和 程序，几乎不可能做到实时预测。 随着多尺度数学理论模型的建立，引发了结构多尺度物理和数值模型的广泛讨 论，其中尺度耦合的有效模拟解决了多尺度有限元建模过程中的关键问题，使得建立既高 效又不失精度的数值模型成为可能，并为耦合基于模型与数据驱动的损伤识别方法提供了 一种途径。
技术实现思路
技术问题：本专利技术的目的是提供一种基于节点曲率和小波分析的变截面梁损伤识 别方法。 技术方案：本专利技术的的具体过 程包括如下步骤： (1)采集在集中荷载匀速通过有损变截面梁时的位移时程数据，并将采集到的位 移时程数据作为初始值通过小波去噪去除有损变截面梁实测位移时程数据中的干扰噪声， 得到优化后有损变截面梁节点k在t时刻的位移时程yk(tr，将其对时间t进行2次微分 得节点k的曲率时程y"k(tr;其中节点k为第i和i-1号单元中间的节点； (2)在单一尺度下建立无损变截面梁的有限元模型，模拟并采集在集中荷载匀速 通过无损变截面梁时的位移时程数据，节点k的弯曲刚度为EIk，t时刻的节点k的位移时 程为yk⑴； (3)综合有损变截面梁和无损变截面梁各节点的位移时程y"k(tr和yk(t)，计算 变截面梁损伤位置，具体如下： 定义参数：有损变截面梁与无损变截面梁的长度为1，均划分为n-1个单元，共n 个节点； EIdPElH分别代表无损变截面梁第i和i-1号单元的弯曲刚度，无损变截面梁 节点k的弯曲刚度为EIk;i?/厂和分别代表有损变截面梁第i和i-1号单元的弯曲刚度， 有损变截面梁节点k的弯曲刚度为; 无损变截面梁的节点k与梁端左侧距离为lk处的单元i在外荷载和环境作用下 的弯矩为M队），根据弯矩的连续性可得： M(lk) =EIiy// =EIi_1y//r(lk) 其中，y"i(lk)和y"Jlk)分别为无损变截面梁距梁端lk处节点k的左侧和右 侧的曲率； 节点k的节点曲率为：【主权项】1. 一种，其特征在于，该损伤识别 方法包括以下步骤： (1) 采集在集中荷载匀速通过有损变截面梁时的位移时程数据，并将采集到的位移时 程数据作为初始值通过小波去噪去除有损变截面梁实测位移时程数据中的干扰噪声，得到 优化后有损变截面梁节点k在t时刻的位移时程y k(tr，将其对时间t进行2次微分得节 点k的曲率时程y" k(tr;其中节点k为两个连续单元中间的节点； (2) 在单一尺度下建立无损变截面梁的有限元模型，模拟并采集在集中荷载匀速通过 无损变截面梁时的位移时程数据，节点k的弯曲刚度为EI k，t时刻的节点k的位移时程为 yk⑴； (3) 综合有损变截面梁和无损变截面梁各节点的位移时程y" k(tr和yk(t)，计算变截 面梁损伤位置，具体如下： 定义参数：有损变截面梁与无损变截面梁的长度为1，均划分为n-1个单元，共η个节 占. EIjP EI η分别代表无损变截面梁第i和i-Ι号单元的弯曲刚度，i = 2, 3, 4,. . .，η ; 无损变截面梁节点k的弯曲刚度为EIk; ￡/]和分别代表有损变截面梁第i和i-1号单 元的弯曲刚度，有损变截面梁节点k的弯曲刚度为； 无损变截面梁的节点k与梁端左侧距离为Ik处的单元i在外荷载和环境作用下的弯 矩为M (Ii)，根据弯矩的连续性可得： M(Ik) =EIiyi" (Ik) =EIiJ" r (Ik) 其中，y" Jlk)和y" r(lk)分别为无损变截面梁距梁端Ik处节点k的左侧和右侧的 曲率； 节点k的节点曲率为：故而得tt;同理，有损截面梁也可得出：假定结构微损伤对其应力重分布影响很小，节点k的弯矩时程积分方程：其中，k = 2, 3…n-1 ;对于有η个节点和n-1个单元组成的结构可得n-2个线性方程， 整合成矩阵可得：根据奇异值分解方法求解得到矩阵D的逆ΙΓ1，从而可得矩阵A = D^1Bi 结合MATLAB程序可计算矩阵A和得到损伤的指标其中，表示损伤前后结构的刚度变化； 其中％即为变截面梁的损伤指标，&1中高的跳跃预示着损伤的发生，并指明了损伤的 位置；定义代表单元i的损伤程度为Iii= 1-a i，计算得到有损变截面梁的有损单元的损伤 程度矩阵 H = [Ii1 h2 h3 . . · hn_J ; (4)将步骤（3)中已经识别出的有损单元进一步平均细分为若干个小单元，建立时空 多尺度有限元模型，分别模拟并采集有损变截面梁与无损变截面梁的各个节点在集中荷载 通过变截面梁的时间内的位移时程数据； 重复步骤（3)，计算得到已经识别出的有损单元细分的若干个小单元的表征损伤指标。2.根据权利要求1所述的，其特征 在于，所述步骤（4)中对已经识别出的有损单元为s个，进一步细分为m+s个；其中，m>2s。【专利摘要】本专利技术提供了一种，该方法基于结构微损伤对其应力重分布影响很小的假定，首先测量在集中荷载匀速通过有损变截面梁时的位移信号，并计算基于无损单一尺度有限元模型求解变截面梁损伤前后的节点曲率矩阵实现损伤位置识别，进而基于时空多尺度模型对结构损伤程度展开精确计算，从而最终得到变截面梁的损伤状态。【IPC分类】G06F17-50【公开号】CN104750926【申请号】CN201510141171【专利技术人】宗周红, 钟儒勉, 杨泽刚, 周海飞, 郑沛娟, 臧建斌 【申请人】东南大学【公开日】2015年7月1日【申请日】2015年3月27日本文档来自技高网...
【技术保护点】
一种基于节点曲率和小波分析的变截面梁损伤识别方法，其特征在于，该损伤识别方法包括以下步骤：(1)采集在集中荷载匀速通过有损变截面梁时的位移时程数据，并将采集到的位移时程数据作为初始值通过小波去噪去除有损变截面梁实测位移时程数据中的干扰噪声，得到优化后有损变截面梁节点k在t时刻的位移时程yk(t)*，将其对时间t进行2次微分得节点k的曲率时程y″k(t)*；其中节点k为两个连续单元中间的节点；(2)在单一尺度下建立无损变截面梁的有限元模型，模拟并采集在集中荷载匀速通过无损变截面梁时的位移时程数据，节点k的弯曲刚度为EIk，t时刻的节点k的位移时程为yk(t)；(3)综合有损变截面梁和无损变截面梁各节点的位移时程y″k(t)*和yk(t)，计算变截面梁损伤位置，具体如下：定义参数：有损变截面梁与无损变截面梁的长度为l，均划分为n‑1个单元，共n个节点；EIi和EIi‑1分别代表无损变截面梁第i和i‑1号单元的弯曲刚度，i＝2,3,4,...,n；无损变截面梁节点k的弯曲刚度为EIk；和分别代表有损变截面梁第i和i‑1号单元的弯曲刚度，有损变截面梁节点k的弯曲刚度为EIk*;]]>无损变截面梁的节点k与梁端左侧距离为lk处的单元i在外荷载和环境作用下的弯矩为M(li)，根据弯矩的连续性可得：M(lk)＝EIiyl″(lk)＝EIi‑1y″r(lk)其中，y″l(lk)和y″r(lk)分别为无损变截面梁距梁端lk处节点k的左侧和右侧的曲率；节点k的节点曲率为：y′′(lk)=12(yl′′(lk)+yr′′(lk))]]>综合可得：M(lk)=EIky′′(lk)=12EIk(M(lk)EIi+M(lk)EIi-1)]]>故而得出：EIk=2(EIiEIi-1EIi+EIi-1);]]>同理，有损截面梁也可得出：EIk*=2(EIi*EIi-1*EIi*+EIi-1*);]]>假定结构微损伤对其应力重分布影响很小，节点k的弯矩时程积分方程：∫t=t1t=t2M(lk)dt=∫t=t1t=t2M*(lk)dt]]>可得：(EIiEIi+EIi-1EIi-1*EIi-1+EIi-1EIi+EIi-1EIi*EIi)∫t=t1t=t2yk′′(t)dt=∫t=t1t=t2yk′′(t)*dt]]>令ai=EIEIi*,,bk=∫t=t1t=t2yk′′(t)dt,bk*=∫t=t1t=t2yk′′(t)*dt,ci=EIiEIi+EIi-1]]>得：[aici+ai-1(1-ci)]bk=bk*]]>类推有：[ai+1ci+1+ai(1-ci+1)]bk+1=bk+1*]]>其中，k＝2,3…n‑1；对于有n个节点和n‑1个单元组成的结构可得n‑2个线性方程，整合成矩阵可得：D(n-2)×(n-1)A(n-1)×1=B(n-2)×1*]]>其中，A＝[a1 a2 a3 … an‑1]B*＝[b2 b3 … bn‑1]根据奇异值分解方法求解得到矩阵D的逆D‑1，从而可得矩阵A＝D‑1B*结合MATLAB程序可计算矩阵A和得到损伤的指标其中，表示损伤前后结构的刚度变化；其中ai即为变截面梁的损伤指标，ai中高的跳跃预示着损伤的发生，并指明了损伤的位置；定义代表单元i的损伤程度为hi＝1‑ai，计算得到有损变截面梁的有损单元的损伤程度矩阵H＝[h1 h2 h3 ... hn‑1]；(4)将步骤(3)中已经识别出的有损单元进一步平均细分为若干个小单元，建立时空多尺度有限元模型，分别模拟并采集有损变截面梁与无损变截面梁的各个节点在集中荷载通过变截面梁的时间内的位移时程数据；重复步骤(3)，计算得到已经识别出的有损单元细分的若干个小单元的表征损伤指标。...
【技术特征摘要】
【专利技术属性】
技术研发人员：宗周红，钟儒勉，杨泽刚，周海飞，郑沛娟，臧建斌，
申请(专利权)人：东南大学，
类型：发明
国别省市：江苏;32
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。