综合曲率和协整识别环境温度影响下的结构损伤.pdf
《综合曲率和协整识别环境温度影响下的结构损伤.pdf》由会员分享,可在线阅读,更多相关《综合曲率和协整识别环境温度影响下的结构损伤.pdf(9页珍藏版)》请在咨信网上搜索。
1、综合曲率和协整识别环境温度影响下的结构损伤周翔袁 周泽文袁 常军渊苏州科技大学 土木工程学院袁江苏 苏州 215011冤摘要院 结构损伤引起的损伤特征参数的变化往往被变化的环境温度引起的损伤特征参数的变化所掩盖袁从而导致基于静载荷的损伤识别方法失效遥 论文提出了综合曲率模态和协整的方法袁在变化的环境温度影响下进行简支梁结构损伤识别研究遥 环境温度变化下袁首先袁对无损结构采集挠度数据并计算成曲率袁选择曲率为协整变量曰其次袁对不同节点的协整变量进行协整曰然后袁将协整所得的协整残差作为损伤指标曰最后袁利用 X-bar 控制图建立置信区间用于评定损伤状态下的协整残差袁识别结构损伤遥 混凝土与工字钢简支
2、梁数值模拟结果表明该方法能够有效地识别环境温度影响下的结构损伤遥关键词院 协整曰曲率曰环境温度曰损伤识别中图分类号院 TU375.1文献标识码院 A文章编号院 2096-3270渊2023冤02-0008-09近 30 年来袁基于静载荷的结构损伤识别方法发展迅速1-3遥 在实际工程中袁环境温度引起的损伤特征参数的变化往往会掩盖由结构损伤引起的损伤特征参数的改变4袁导致损伤识别结果出现偏差袁从而影响桥梁维护策略的正确性遥 因此袁变化的环境温度对结构损伤识别的影响已成为国内外学者日渐研究的热点5遥目前袁相关研究方法大致可以分为两类院一是根据桥梁健康监测信号建立环境温度与结构响性的对应关系袁例如院插
3、值分析尧随机子空间等方法6-7袁该方法的缺点是需要记录环境变量的值袁实际使用中存在困难曰二是根据监测信号识别结构损伤袁例如院主成分分析法尧自适应神经网络和协整法等8-10袁这种方法的优点是不需要记录环境变量的值袁使用简单且适用性更高遥源自于计量经济学的协整被引入结构健康监测领域后袁可以有效地去除环境温度的影响11遥Cross 等人12以频谱幅值作为协整变量用于复合板结构损伤识别袁有效地消除了环境温度对损伤识别的影响且能准确地识别损伤曰Liu 等人13以梁伸长率作为协整变量用于钢架桥的损伤识别袁有效地消除了环境温度对损伤识别的影响并能较好地识别出损伤曰周翠14以挠度为协整变量用于识别简支梁结构的
4、损伤袁有效地消除了温度对损伤识别的影响且能定位损伤单元遥 梁亚斌15将频率作为协整变量袁以简支梁结构为研究对象进行损伤识别袁剥离出了温度效应对损伤识别的影响袁并能较快地识别出损伤遥 曹亚东等人16以 AR 模型系数为协整变量用于海洋平台结构损伤识别袁有效地消除了环境温度对损伤识别的影响且能很好地定位损伤单元遥 这些研究大大推进了协整在温度影响下结构损伤识别的发展袁但已有的研究成果仍不能满足现阶段对于损伤识别的需求袁且无法在实际工程应用中很好地发挥作用遥 例如袁以挠度为协整变量会引起无损梁单元节点的损伤指标发生较大的变化袁导致损伤识别步骤繁琐且时间较长曰以 AR 模型系数为协整变量识别结构损伤时
5、存在一定的局限性遥 本文提出以曲率作为协整变量袁将曲率与协整相结合袁研究了环境温度对简支梁结构损伤识别的影响袁通过数值模拟的方式袁探讨并提出的曲率模态和协整的综合方法的可行性袁为在实际工程应用中提供参考遥要要要要要要要要要要要要要要要要要眼收稿日期演 圆园22原05-17眼作者简介演 周翔渊1996-冤袁男袁江苏宿迁人袁硕士研究生遥眼基金项目演 国家自然科学基金项目渊51508368冤曰江苏省研究生科研与实践创新计划项目渊SJCX20-1112冤通信联系人院常军渊1973-冤袁男袁教授袁博士袁主要从事健康监测与振动控制研究袁Email院遥第 36 卷第 2 期苏州科技大学学报渊工程技术版冤灾燥
6、造援 36晕燥援 2圆园23 年 6 月允燥怎则灶葬造 燥枣 Suzhou 哉灶蚤增藻则泽蚤贼赠 燥枣 杂糟蚤藻灶糟藻 葬灶凿 栽藻糟澡灶燥造燥早赠渊耘灶早蚤灶藻藻则蚤灶早 葬灶凿 栽藻糟澡灶燥造燥早赠 Edition冤Jun援 圆园23第 2 期1基本原理1.1协整理论1978 年 Engle 和 Granger 提出协整理论及其方法17袁为非平稳时间序列变成平稳时间序列提供了途径遥如果序列 y 通过d 次差分变成平稳序列袁而d-1 次差分却不平稳袁那么称序列y 为d 阶单整序列袁可为记yI渊d冤遥 在 k 维向量时间序列yt=y1t袁y2t袁噎袁yktT渊t=1袁2袁噎袁T冤中袁其分量序列
7、皆为非平稳时间序列遥 若存在一个非零向量 茁 使得序列 yt的线性组合茁Tyt=ut渊1冤满足平稳性要求袁即残差 ut符合零均值尧同方差的白噪声序列袁那么称向量 茁T为序列 yt的协整向量遥 本文采用 Johansen 协整检验确定 茁18遥1.2ADF 单位根检验Johansen 协整检验的前提条件是协整变量序列为一阶单整序列遥 序列阶数平稳性的检验方法通常采用Augmented Dickey-Fuller渊ADF冤检验遥 检验模型驻yt=a+浊yt-1+p-1i=1移茁i驻yt-i+ut渊2冤浊=p-1i=1移准i-1袁茁i=-pj=i+1移驻准j渊3冤式中袁驻 为差分算子袁定义为 驻yt
8、-i=yt-i-yt-i-1曰浊 和为 茁i参数曰a 为常数项曰p 为滞后阶数遥 参数 浊 的计算通常采用最小二乘法曰p 值通常由 AIC渊Akaike Information Criterion冤准则确定遥平稳性检验如下t浊=浊赞/滓浊渊4冤浊赞是通过 浊 的最小二乘法估计得到曰滓浊是 浊 的标准差估计遥将计算得到的 t浊与 DF 表19中临界值 t琢进行比较袁如果 t浊t琢袁则表明原序列为非平稳的曰再对一阶差分后的序列 驻yt进行平稳性检验袁此时 t浊t琢袁则说明原序列 yt为一阶单整平稳序列遥1.3Johansen 协整检验协整检验旨在找到一个合适的协整向量将非平稳时间序列变成平稳时间序
9、列遥 检验步骤如下袁首先建立一个 VAR渊p冤模型袁变量为 k渊k逸2冤个情形院yt=椎1yt-1+噎+椎pyt-p+Hxt+ut渊t=1袁2袁噎袁T冤渊5冤式中袁yt的各分量都是非平稳 I渊1冤的变量曰xt是一个确定的 d 维外生向量袁代表趋势项尧常数项等确定项曰ut是 k 维扰动向量遥 式渊5冤经变换可得驻yt=装yt-1+p-1i=1移祝i驻yt-i+Hxt+ut渊6冤装=p-1i=1移椎i-I袁祝i=-pj=i+1移椎i渊7冤将 yt的协整检验转变为对矩阵 装 的分析问题袁这也是 Johansen 协整检验的基本原理遥最后讨论 装 的秩为 r渊0臆r臆k冤的情形渊其它无意义冤袁可以将
10、装 分解为 装=琢茁T袁其中 琢 和 茁 都是 k 行 r 列的矩阵袁茁 的列向量即为目标协整向量袁通常采用 茁 第一列向量以取得最好的协整效果遥1.4X-bar 控制图当基准样本的残差序列符合正态分布时袁可通过准则得到残差序列的控制限值袁以此建立置信区间4CL=NI渊中心线冤渊8冤UCL=CL+3滓渊控制上限冤渊9冤周翔袁等院综合曲率和协整识别环境温度影响下的结构损伤92023苏州科技大学学报渊工程技术版冤LCT=CL-3滓渊控制下限冤渊10冤式中院CL 为残差的均值曰滓 为残差的标准差值遥 对于检测样本袁当残差序列处于控制限之内时说明该节点没有发生损伤袁当残差序列超出控制限值后则判断为该节
11、点发生了损伤遥2综合曲率和协整的损伤识别曲率的物理意义表示曲线在某一点的弯曲程度袁由材料力学可知其公式表达为籽=M/渊EI冤渊11冤式中袁籽 为截面的曲率曰M 为的弯矩曰E 为弹性模量遥 在截面弯矩 M 不变情况下袁截面刚度 EI 决定了曲率 籽 的大小袁而影响 EI 的因素有单元损伤的设定与环境温度的影响遥 因此在考虑为曲率识别简支梁损伤时受到环境温度的影响下袁利用协整法去除环境温度的影响遥 本文采用曲率为协整变量的原因是在满足作为协整变量的条件下袁相对于位移袁曲率对局部损伤更加敏感遥基于环境温度的影响袁将结构采集的挠度数据计算成曲率袁选择曲率为协整变量袁并对曲率进行协整处理袁以协整残差作为
12、损伤指标袁对去除温度影响的理论与定位损伤理论公式进行推导说明袁并将其应用于简支梁在环境温度影响下的损伤识别中遥2.1简支梁推导混凝土简支梁模型如图 1 所示袁在跨中集中荷载 P作用下袁损伤发生在左半跨处遥 则简支梁损伤前后的曲率20可表达为式渊12冤尧渊13冤遥准u渊x冤=Px2EI渊0臆x臆L2冤渊12冤准d渊x冤=1EI1+琢啄渊x-x0冤P渊x-xp冤U渊x-xp冤-P渊L-xp冤Lx渊0臆x臆L2冤渊13冤式中袁m 表示第 m 个测点曰xp表示荷载作用位置曰L 表示简支梁跨度曰琢 为附加柔度系数袁表示损伤程度遥U渊x冤为阶跃函数袁当 xxp冤时袁U渊x-xp冤=0曰当 x逸xp时袁U渊
13、x-xp冤=1遥 啄渊x-x0冤为狄拉克函数21袁当 x屹x0时袁啄渊x-x0冤=0曰当 x=x0时袁啄渊0冤=肄遥 EI 表示简支梁的截面抗弯刚度遥 由文献1可知曲率模态反映了结构的固有特性遥协整法中袁在无损状态下利用参考节点与待测节点之间的曲率关系建立协整方程遥 假设与为参考测点和检测测点的曲率袁则有茁m准m-茁n准n=ud渊14冤式中袁u 表示为符合零均值尧同方差的白噪声序列遥 式渊13冤对应于论文 1.1 节协整理论中的式渊1冤袁茁m尧茁n表示为协整系数遥 由此可知袁简支梁任意两个节点曲率序列满足协整关系袁当对应节点发生损伤时袁其相应的协整关系也会遭到破坏而发生改变袁此时袁原来的协整关
14、系不再满足现在的要求袁转而将损伤后的简支梁结构第 m尧n 节点的曲率序列代入原协整方程时袁u 将变化为 ud遥茁m准m袁d-茁n准n袁d=ud渊15冤当损伤发生时袁曲率变化可表达为驻准=准d-准u=Px2EI琢啄渊x-x0冤x渊0臆x臆L2冤渊16冤由式渊16冤可知袁当 x屹x0时袁啄渊x-x0冤=0曰所以 驻准 几乎无变化曰当 x=x0时袁啄渊0冤=肄袁所以 驻准 发生较大变化遥 此时袁对式渊14冤尧式渊15冤作差得到测点曲率的协整残差的变化为驻ud=茁m驻准m-茁n驻准n渊17冤此处的协整残差 驻ud=u-ud袁且 驻ud屹0遥 从而可得袁协整残差与温度变化表现无关袁可以有效地去除温度的影
15、图 1两端简支梁模型10第 2 期图 2北京该地的室外温度图 3混凝土和钢的弹性模量随温度的变化关系渊a冤混凝土渊b冤钢图 4温度与结构响应关系图图 5温度影响下曲率幅值变化图样本点T/益T/益0200400600800100012005.0E-035.2E-035.4E-035.6E-035.8E-036.0E-036.2E-036.4E-036.6E-036.8E-037.0E-03挠度曲率(1/m)样本点基准损伤MAXMIN样本点响袁当简支梁结构发生损伤时袁原来的协整关系将会遭到破坏袁不再适应现在的协整要求袁进而使得该节点的协整残差序列 驻ud发生大幅度变化袁因此可以判断为该点发生了损伤
16、遥2.2损伤识别步骤由以上理论推导可以看出院当梁单元发生损伤时袁单元相邻节点的检测样本的残差幅值会同时超出控制限值袁而无损单元相邻节点的检测样本的残差幅值则不会出现此种状况遥 具体识别步骤如下院渊1冤无损状态下袁在变化的环境温度中采集简支梁单元节点的挠度数据曰渊2冤选择挠度数据计算出曲率并生成曲率序列曰渊3冤采用 ADF 单位根检验方法袁检验各个节点的曲率序列是否满足一阶单整序列条件曰渊4冤对满足一阶单整的曲率数列袁随机选择一个节点曲率序列作为参考序列袁并依次对每个需要检测的节点曲率序列进行 Johansen 协整检验以确定协整向量 茁曰渊5冤使用协整向量对各个节点曲率序列建立协整方程袁生成协
17、整残差序列曰渊6冤用 X-bar 控制图为协整残差序列建立置信区间袁用于评定损伤状态下的协整残差序列袁以此识别结构是否发生损伤遥3数值模拟3.1环境因素及其影响简支梁数值模拟的环境温度由两部分组成袁一是基准工况的温度样本袁二是待检工况的温度样本遥 基准工况温度样本采用了 2019 年 1 月 1 日到 2019 年 12 月 31 日北京某地的室外温度数据袁该温度数据记录为每天 8 点尧14 点和 24 点的温度袁共计 1 095 个样本点遥 待检工况温度样本由-1434 益每隔 3 益取样 8 次袁共计 136 个样本点袁如图 2 所示遥 混凝土和钢的弹性模量与温度的关系如图 3 所示4遥当
18、环境温度不变时袁结构响应在简支梁单元损伤后发生了改变袁如图 4 中 1尧3 点所示22袁所以袁结构响应可以评估简支梁单元的安全状态遥 但在环境温度发生变化时袁即使简支梁单元发生了损伤袁结构响应也可能不会发生改变袁如图 4 中 2尧3 点所示袁此时袁结构响应将无法评估简支梁单元的安全状态遥 因此袁环境温度的改变可能会掩盖结构的损伤状态遥 在本文简支梁数值模拟中袁随机提取一组单元节点的基准数据袁如图 5 所示袁随着样本点的变化袁曲率幅值发生了较大变化遥 因此袁直接使用曲率作为损伤指标将无法正确评估简支梁结构的安全性遥周翔袁等院综合曲率和协整识别环境温度影响下的结构损伤112023苏州科技大学学报渊
19、工程技术版冤图 6混凝土简支梁数值模型渊单位院mm冤表 1数值模拟损伤工况损伤工况 损伤梁单元 弹性模量折减/%一 9 10 二 6、9 10 图 7无噪声下工况一识别结果渊a院节点 1曰b院节点 2曰c院节点 3曰d院节点 4冤020040060080010001200-0.00048-0.00040-0.00032-0.00024-0.00016-0.000080.000000.000080.000160.000240.000320.000400.00048残差样本点JD7020040060080010001200-0.00048-0.00040-0.00032-0.00024-0.000
20、16-0.000080.000000.000080.000160.000240.000320.000400.00048残差样本点JD8020040060080010001200-0.00048-0.00040-0.00032-0.00024-0.00016-0.000080.000000.000080.000160.000240.000320.000400.00048残差样本点JD9020040060080010001200-0.00048-0.00040-0.00032-0.00024-0.00016-0.000080.000000.000080.000160.000240.000320.0
21、00400.00048残差样本点JD10渊a冤渊b冤渊c冤渊d冤020040060080010001200-0.0006-0.0005-0.0004-0.0003-0.0002-0.00010.00000.00010.00020.00030.00040.00050.0006残差样本点JD7020040060080010001200-0.0006-0.0005-0.0004-0.0003-0.0002-0.00010.00000.00010.00020.00030.00040.00050.0006残差样本点JD8020040060080010001200-0.0006-0.0005-0.0004
22、-0.0003-0.0002-0.00010.00000.00010.00020.00030.00040.00050.0006残差样本点JD9020040060080010001200-0.0006-0.0005-0.0004-0.0003-0.0002-0.00010.00000.00010.00020.00030.00040.00050.0006残差样本点JD10渊a冤渊b冤渊c冤渊d冤图 820%噪声水平下工况一识别结果渊a院节点 7曰b院节点 8曰c院节点 9曰d院节点 10冤3.2混凝土梁模型为了解决提出的综合曲率和协整方法适用于变化的环境温度和运行状况的损伤问题袁在通用有限元软件
23、ABAQUS 中建立了混凝土简支梁数值模型袁如图 6 所示遥混凝土简支梁模型基本参数为密度 2 300 kg/m3袁泊松比 滋=0.2袁总长为 1.8 m尧单元长度 0.09 m袁截面高宽尺寸为 0.16 m伊0.2 m袁跨中沿梁顶竖直向下施加集中荷载 P=30 kN遥 采用了 2 种损伤工况袁工况一为单点损伤袁工况二为多点损伤遥 具体损伤工况设定见表 1遥渊1冤工况一遥 工况一为单点损伤袁考虑为验证所提方法在变化的环境温度下和跨中集中荷载作用下对混泥土梁识别损伤与定位损伤问题的有效性遥 损伤的模拟在 9 单元袁即节点 8 与 9 之间遥 混泥土梁单点损伤识别结果如图 7 所示遥由图 7 可知
24、袁节点 7尧10 检测样本残差序列均处于控制限值内袁节点 8 与节点 9 检测样本残差序列同时超出了控制限值袁由本文所提方法可知袁单元 9 发生了损伤袁识别结果正确遥 其它单元节点检测样本残差序列均处于控制限值内袁单元无损袁限于篇幅限制及重复性已省略遥上面的模拟没有考虑系统噪声的影响遥 为了验证该方法的抗噪性袁添加一定量的噪声遥 噪声模型参考文献14中公式袁引用公式Sn=S+琢伊A伊R渊18冤A=rms渊准max冤-rms渊准min冤2渊19冤其中袁Sn添加噪声后的信号曰S 为真实测量信号曰琢 为噪声水平曰rms 为均方根曰准max和 准min分布表示曲率序列的最大和最小包络曰R 表示独立正态
25、分布变量袁满足零均值尧单位标准差遥在相同条件下袁 分别将 10%尧15%和 20%的噪声加入到曲率序列中袁 所有测试都表现良好遥 此处给出20%噪声水平的检测结果袁如图 8 所示遥 与没有噪声的结果相比袁协整的残差序列幅值波动更大袁控制限值也相应增大袁但噪声并没有掩盖损伤信息袁表明该方法具有较强的抗噪性遥渊2冤工况二遥工况二为了验证所提方法对于混凝土梁模型同时识别多处损伤的准确性袁设定为多点损伤袁损伤模拟在 6 单元和 9 单元袁即节点 5 与 6 和节点 8 与 9 之间袁且在提取的结构位移序列中加入 20%高斯白噪声后计算为曲率变量袁然后根据所提方法识别混凝土梁单元的损伤袁识别结果如图 9
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 综合 曲率 和协 识别 环境温度 影响 结构 损伤
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。