收稿日期:2022-08-11ꎮ
基金项目:国家自然科学基金项目(51378112)ꎻ福建省建设科技项目(2017-K-59)ꎮ
作者简介:林友勤(1973 )ꎬ男ꎬ高级实验师ꎬ博士ꎬ研究方向为结构健康监测㊁损伤诊断与结构异常识别㊁桥梁检测和加
固ꎮE ̄mail:lyq@fzu.edu.cnꎮ
林友勤ꎬ郑学善ꎬ余印根ꎬ等.基于BP神经网络的斜拉索损伤识别方法[J].南昌大学学报(工科版)ꎬ2023ꎬ45(4):370-376.LINYQꎬZHENGXSꎬYUYGꎬetal.DamageidentificationmethodforstaycablebasedonBPneuralnetwork[J].JournalofNanchang
University(Engineering&Technology)ꎬ2023ꎬ45(4):370-376.
基于BP神经网络的斜拉索损伤识别方法
林友勤1ꎬ郑学善2ꎬ余印根3ꎬ王志俸1
(1.福州大学土木工程学院ꎬ福建福州350116ꎻ2.福建省高华建设工程有限公司ꎬ福建福清350301ꎻ
3.福建省永正工程质量检测有限公司ꎬ福建福州350001)
㊀㊀摘要:斜拉桥拉索损伤识别中不仅需要拉索发生损伤时自身索力的变化ꎬ还需要考虑拉索发生损伤引起的其他拉索索力变化ꎬ以单索索力变化为指标的方法难以适应斜拉桥的拉索损伤识别ꎮ以不同损伤下全桥索力变化率作为输入向量ꎬ拉索损伤位置和损伤程度工况为输出向量ꎬ建立了较为精准的BP神经网络模型ꎮ对斜拉桥模型在发生单根和2根拉索损伤工况时进行损伤识别ꎬ并随机选取不同损伤工况对BP神经网络模型的准确性进行验证ꎬ拉索损伤识别误差在6.0%以内ꎮ结果表明:以拉索索力变化率为指标ꎬ基于BP神经网络模型能够有效地识别斜拉桥拉索的损伤位置和损伤程度ꎬ为桥梁安全运营提供保障ꎮ
关键词:斜拉桥ꎻ拉索ꎻ神经网络ꎻ索力ꎻ损伤识别
中图分类号:U466.2㊀㊀㊀㊀文献标志码:A㊀㊀㊀㊀文章编号:1006-0456(2023)04-0370-07
DamageidentificationmethodforstaycablebasedonBPneuralnetwork
LINYouqin1ꎬZHENGXueshan2ꎬYUYingen3ꎬWANGZhifeng1
(1.CollegeofCivilEngineeringꎬFuzhouUniversityꎬFuzhou350116ꎬChinaꎻ
2.FujianGaohuaConstructionEngineeringCo.ꎬLtd.ꎬFuqing350301ꎬFujianꎬChinaꎻ
3.FujianYongzhenConstructionQualityInspectionCo.ꎬLtd.ꎬFuzhou350001ꎬChina)
Abstract:Forcable ̄stayedbridgesꎬthedamageofasinglecablewillcausethestressredistributionofthewholecables.Althoughu ̄
singfrequencyofsinglecableasthedamageindexwilldetectabnormalitiesinmultiplecordsꎬitisimpossibletodeterminethelocationanddegreeofdamagea
ccurately.ToaddressthischallengeꎬaneffectivemethodwasexploredtorecognizecabledamagelocationanditsdegreethroughBPneuralnetworkinthisstudy.Intotalꎬabenchmarkfiniteelementmodelofadouble ̄towerdouble ̄cablesemi ̄buoycontinuouscompositebeamcable ̄stayedbridgewasestablished.Thedamageconditionsofsinglecableandtwocablesweresimulatedthroughthechangeoftheelasticmodulus.ABPneuralnetworkmodelwasestablishedbasedoncableforcegradientasidentificationin ̄dexofcabledamageanditsaccuracywasverified.Adamageidentificationmethodforlong ̄spancable ̄stayedbridgeswasestablished.Theresultsshowedthattheestablishedneuralnetworkmodelidentifythedamagelocationanddamagedegreeofstaycableswithhighaccuracy.Theidentificationerrorwaswithin6.0%ꎬsotheBPneuralnetworkcanbeusedforcabledamageidentificationofcable ̄stayed
bridges.
KeyWords:cable ̄stayedbridgeꎻcableꎻneuralnetworkꎻcableforceꎻdamageidentification
㊀㊀改革开放以来ꎬ我国经济快速发展ꎬ桥梁网络建设不断扩大ꎮ然而桥梁在使用的过程中由于外界环境㊁交通荷载㊁疲劳等因素的作用ꎬ会引起材料的老化和结构承载能力的降低ꎬ这一直都是影响结构的安全性和适用性的重要问题ꎮ对于斜拉桥而言ꎬ拉
索在服役过程中是相互协同工作的ꎬ当某一根拉索发生损伤时ꎬ会引起整个拉索的索力发生一定的变化ꎮ那如何利用智能算法将索力变化与结构损伤建立起某种关系ꎬ从而实现通过索力的变化来对结构的损伤进行识别并作出全面综合判断和评估值得进
第45卷第4期2023年12月
㊀㊀㊀㊀㊀㊀
南昌大学学报(工科版)
JournalofNanchangUniversity(Engineering&Technology)
Vol.45No.4Dec.2023
一步探索ꎮ
人工神经网络(artificialneuralnetworkꎬANN)
是从信息处理角度对人脑神经元网络进行抽象ꎬ建
立某种简单模型ꎬ主要的用途是做数据处理和各种
参数拟合ꎮ根据环境的变化ꎬ对权值进行调整ꎬ改善
系统的行为ꎬ建立起大数据之间的相互关系ꎬ最终达
到数据处理的要求ꎮ神经网络具有强大的计算能
力㊁联想能力㊁适应性㊁容错能力和自组织能力ꎬ因而
它在结构损伤识别领域受到了广泛的关注和研究ꎮ
在众多神经网络中ꎬBP神经网络应用最为广泛[1]ꎮBP神经网络具有非线性映射㊁自我组织和并行处理等特点[2-3]ꎬ特别是在大数据处理和解决复杂问题能力方面与其他传统诊断方法相比具有很大优势ꎬ
在结构损伤识别中得到了有效的应用[4-6]ꎮ
李忠献等[7]用自组织竞争神经网络迅速准确地识别出子结构发生损伤后ꎬ再通过BP神经网络所建立的结构损伤识别方法ꎬ能够进一步准确识别子结构中结构损伤的位置和程度ꎮ杨杰等[8]通过固有频率与拉索损伤的对应模式ꎬ建立BP神经网络对斜拉桥拉索损伤识别定位ꎬ但对于损伤程度的识别不敏感ꎮ孙宗光等[9]在已经确定了结构损伤区域的情况下ꎬ应用BP神经网络实现了结构的损伤位置及其损伤程度识别的方法ꎮ谭冬梅等[10]将能量累积变异值作为特征值ꎬ再将AdaBoost算法和BP神经网络相结合ꎬ在有噪声的情况下ꎬ仍可有效地识别出大跨斜拉索的损伤程度ꎮ
合理的损伤指标是进行拉索损伤识别的关键ꎬ
这方面也开展了相关研究ꎮ闫维明等[11]通过计算索力在最不利荷载效应时相对恒载索力值的变化ꎬ确立吊索与主缆的索力预警阈值ꎬ结合对索力异常值的成因分析ꎬ确立悬索桥的损伤预警流程ꎮ李延强等[12]通过斜拉索张力指标实现了对斜拉桥主梁不同位置㊁不同程度单点损伤的识别ꎮHua等[13]通过索力变化来检测桥梁主梁是否发生损伤ꎬ结果表明当使用无噪声索力时ꎬ可以正确识别桥梁主梁中的损坏
位置和损坏程度ꎬ在测量噪声低的情况下也能对结构的损伤进行良好的识别ꎮ
可以看出ꎬBP神经网络可以对结构损伤进行有
效的识别ꎬ索力是拉索损伤识别的常用指标ꎬ但其作
为一个绝对量ꎬ通常只针对所研究的某个桥梁ꎬ不具
有通用性ꎮ此外ꎬ拉索损伤还会引起其他拉索索力
的变化ꎬ多根拉索的损伤则会引起更复杂的变化ꎮ
针对这些问题ꎬ本文以1座斜拉桥单根和2根拉索
损伤识别为研究对象ꎬ尝试以无量纲的索力变化率作为拉索损伤识别指标ꎬ结合BP神经网络模型ꎬ进行大跨度斜拉桥的拉索损伤识别方法研究ꎮ
正则化损伤识别matlab1㊀BP神经网络
㊀㊀神经网络的分类方式有学习方式㊁时间特性㊁模型结构等ꎬ从学习方式可分为监督学习㊁无监督学习以及半监督学习ꎻ从时间特性可分离散型与连续型ꎻ从模型结构角度分类是最为常见的分类形式ꎬ将神
经网络分为前向网络和反馈网络(图1)ꎮBP神经网络属于一种按照误差逆向传播算法训练的多层前馈神经网络ꎬ它包括了输入层㊁隐含层㊁输出层ꎬ通过正向传播和误差反向传播ꎬ对各层的权值进行不断的修正ꎬ直至网络误差平方小于预先设定的阈值ꎬ是一种有导师学习的网络模型[14]
隐含层
输入层输出层
(a)前向网络
隐含层
输入层输出层
(b)反馈网络
图1㊀神经网络
Fig.1㊀Neuralnetwork
1986年ꎬMcClelland和Rumelhaart定义了δ规则ꎬ其核心内容就是计算目标值与计算值的差值ꎬ通过不断地调整函数单元之间的连接强度来减小这个差值ꎮ
δ=[dj-f(WTjX)]fᶄ(WTjX)(1)式中:dj为期望输出值ꎻW为神经元的连接权值矩阵ꎻX为神经元矩阵ꎻfᶄ(WTjX)是传递函数的导数ꎬ显然可以看出δ学习规则需要激活函数可导ꎮ首先定义输出与期望输出之间的平方误差函数E为
E=(1/2)[dj-f(WTjX)]2(2)㊀㊀如有要以最快的速度降低误差的话就是Wj应该与误差的负梯度成正比
173
第4期㊀㊀㊀㊀㊀林友勤等:基于BP神经网络的斜拉索损伤识别方法
ΔWj=-ηÑE
(3)㊀㊀其中误差梯度为
ÑE=-[dj-f(WTjX)]fᶄ(WT
jX)X
(4)
㊀㊀代入可得
ΔWj=η[dj-f(WTjX)]fᶄ(WT
jX)X=ηδX(5)进而得到ΔWj中的分量调整公式
Δwij=η[dj-f(WT
X)]fᶄ
(WTj
X)xi=ηδxi(6)
2㊀斜拉桥模型建立与验证
2.1㊀工程背景
模拟1座双塔双索面半漂浮连续组合梁斜拉桥(图2)ꎮ桥梁的主跨布置为(32.9+115.4+340.0+
115.4+32.9)mꎬ全长为636.6mꎮ斜拉索为双索面
结构ꎬ共104根拉索ꎬ拉索均采用OVM250系列环氧涂层钢绞线ꎬ标准强度为1860MPaꎬ塔上索距0.6mꎬ梁上索距3.5mꎬ与索塔和主梁均采用钢锚箱锚固ꎮ主梁采用的是Q345C的单箱三室扁平流线型栓焊钢箱梁ꎬ桥面宽为32.4mꎬ桥面坡度为2%ꎬ纵向梁和横梁是由接头㊁高强度螺栓相
连接ꎮ索塔采用的C50混凝土ꎬ塔柱的截面形式为单箱单室的空心截面ꎬ塔底左右塔柱中心间距35.6mꎬ总高为
96.54mꎮ桥面与钢梁采用剪力螺栓连接
图2㊀某斜拉桥现场照片
Fig.2㊀Sitephotosofacable ̄stayedbridge
2.2㊀有限元模型建立
利用有限元软件ANASYS建立斜拉桥的实体模型ꎬ模型索力采用竣工后成桥索力定义ꎮ其中ꎬ
1)采用Beam188单元对有限元模型中的纵梁㊁横梁和索塔进行模拟ꎮ2)采用Shell63板壳单元对桥面板进行模拟ꎮ3)桥梁的二期恒载直接转换成质量块ꎬ采用Mass21单元进行模拟ꎬ设置每个质量元在XYZ3个平动方向的质量为8013.25kgꎮ4)支座用Combin14弹簧单元模拟ꎮ5)承台采用的是Solid45实体单元进行模拟ꎬ将墩底的每个节点的6个自由度全部约束ꎮ
全桥斜拉索共有104根ꎬ单侧52根ꎬ从南至北方向定义拉索编号ꎬ右侧依次为R1~R52ꎬ左侧依次
为L1~L52ꎬ成桥左右对称拉索的索力是相同的ꎮ在ANASYS斜拉桥模型中的拉索初始应力是通过定义拉索的初应变所得到[15]ꎮ斜拉桥模型单元数
13060个ꎬ节点数总计58221个ꎮ基本材料参数见表1ꎮ
表1㊀结构模型材料参数
Tab.1㊀Materialparametersofstructuralmodel结构部位
材料
弹性模量/
GPa密度/(kg m-3
)泊松
比混凝土桥塔混凝土C5034.525000.2钢主梁钢材21078000.3桥面板混凝土C5034.525000.2斜拉索钢绞线19578000.3边墩
混凝土C40
32.5
2500
0.2
2.3㊀模型验证
为验证斜拉桥ANASYS有限元模型(图3)的准确性ꎬ通过对有限元模型计算频率和桥梁的实测频率值进行比较ꎬ从表2可知ꎬ在竖向㊁横向㊁扭转这3个振动的方向的最大频率值误差为11%ꎬ在可允许误差范围之内ꎬ斜拉桥ANASYS有限元模型与实际结构相逼近ꎮ
表2㊀大桥模态参数
Tab.2㊀Modalparametersofbridge
振动方向
频率阶次f1/Hzf2/HzE/%竖向
一阶
0.380.380二阶0.500.492三阶0.770.735四阶
0.880.880横向一阶0.390.428二阶0.500.4410扭转
一阶0.630.653二阶
0.81
0.72
11
㊀㊀注:f1为实测值ꎬf2为计算值ꎬE为实测值与计算值的误差
图3㊀斜拉桥有限元模型
Fig.3㊀Finiteelementmodelofcable ̄stayedbridge
273 南昌大学学报(工科版)2023年㊀
3㊀斜拉索损伤识别
3.1㊀损伤指标的定义
在大跨度斜拉桥中ꎬ拉索发生损伤的同时所引起索力的变化也是十分显著的ꎬ将索力的变化与智能算法相结合运用于斜拉桥拉索的损伤识别具有更加实际的工程意义ꎮ定义拉索损伤指标索力变化率[16]见式(7)ꎮ
RCOF=T-Tᶄ
(7)式中:RCOF㊁T㊁Tᶄ分别表示为索力变化率㊁未发生损伤时的拉索索力㊁发生损伤时的拉索索力
通过拉索发生损伤前后的索力的变化计算得出全桥拉索的索力变化率ꎬ建立索力变化率与结构损伤的BP神经网络模型ꎬ进而可根据拉索的索力的变化来判断拉索损伤的位置和损伤程度ꎮ
3.2㊀斜拉索损伤识别流程
基于BP神经网络的拉索损伤识别过程可分为以下3个部分:
1)获取样本数据库ꎮ建立斜拉桥基准有限元模型ꎬ通过模拟拉索在发生不同损伤情况下的索力变化率以及拉索的损伤位置和损伤程度ꎬ以此作为样本数据库ꎮ
2)建立BP神经神经网络ꎮ以索力变化率为输入向量ꎬ相对应的损伤位置和损伤程度为输出向量ꎬ对样本进行训练ꎬ通过不断地迭代和修正得到了较为精准的BP神经网络系统ꎮ
3)实际桥梁的验证ꎮ在不同的外界激励下对实际桥梁的拉索进行振动测试ꎬ获取到的振动信号进行频谱分析ꎬ计算出拉索的实际索力ꎬ进而得到作为BP神经网络的输入向量的索力变化率ꎬ最后得到桥梁发生损伤及损伤的位置和程度ꎬ具体流程见图4ꎮ
3.3㊀BP神经网络的输入向量
在斜拉桥拉索损伤识别当中ꎬ由于该斜拉桥拉索是相互对称ꎬ故拉索损伤的模拟方式选取26根拉索进行分析ꎬ输入向量和输出向量均为26维ꎬ输入向量代表的是26根拉索的索力变化率ꎬ输出向量为26根斜拉索相对应的损伤程度和损伤位置参数ꎮ现模拟单根拉索损伤程度依次为弹性模量损失为0㊁10%㊁20%㊁30%㊁40%㊁50%㊁60%㊁70%㊁80%㊁90%这10种损伤ꎬ共有26ˑ10=260个工况ꎮ2根拉索损伤是从拉索L1~L13中选取2根拉索ꎬ共78组ꎬ每组拉索的损伤程度分别为弹性模量损失40%㊁40%ꎬ
建立基准桥梁有限元模型
不同损伤参数的索力变化率
训练样本的输入
网络权值和阈值的确定
隐含层输出层误差计算
权值和阈值的更新
是否满足要求
开始
实际桥梁的拉索
振动测试
振动信号的频谱
分析
根据索力公式
计算索力
BP神经网络
系统
拉索损伤识别否是
图4㊀斜拉桥神经网络损伤识别流程
Fig.4㊀Damageidentificationprocessof
neuralnetworkforcable ̄stayedbridge
40%㊁80%ꎬ80%㊁40%和80%㊁80%ꎬ总共有312个工况ꎮ总计572个拉索损伤工况数据作为神经网络模型的输入向量ꎬ图5列出了斜拉桥在单根拉索L1㊁L2㊁L5㊁L6在损伤程度为0㊁10%㊁20%㊁30%㊁40%㊁50%㊁60%㊁70%㊁80%㊁90%时的索力变化率数值和2根拉索L11㊁L12分别发生80%㊁40%与L9㊁L12分别发生80%㊁80%损伤时拉索的索力变化率ꎬX轴㊁Y轴㊁Z轴分别表示的是损伤程度D㊁拉索编号n㊁索力变化率RCOFꎮ
3.4㊀BP神经网络的建立
通过对有限元的模型修正ꎬ得到较为精准的有限元模型ꎬ为获取拉索发生损伤时的样本数据库ꎬ则通过降低拉索的弹性模量来模拟拉索损伤ꎬ通过上述式(7)所定义的索力变化率ꎬ构建出单根拉索发生损伤和2根拉索发生损伤时的索力变化率数据库ꎬ以编号为L1~L26的26根拉索的索力变化率作为神经网络的输入向量Iꎬ见式(8)ꎬ拉索损伤位置和损伤程度则作为输出向量Oꎬ见式(9)ꎬ最终建立起较为精准的BP神经网络模型(图6)ꎮ
I={RCOF1ꎬRCOF2ꎬ ꎬRCOF25ꎬRCOF26}(8)
O={k1ꎬk2ꎬ ꎬk25ꎬk26}(9)㊀㊀通过多次训练发现ꎬ当神经网络隐含层的节点数量为100层时该模型运行效果最佳ꎬ建立26ˑ100ˑ26的BP神经网络模型ꎮ首先按照LM最优化算法的网络对已构建的样本数据库进行训练ꎬ之后再做进一步的损伤位置和损伤程度的识别ꎮ本文BP网络结构中采用MSE函数来对网络的性能进行
373
第4期㊀㊀㊀㊀㊀林友勤等:基于BP神经网络的斜拉索损伤识别方法
R C O F /%
D /
%n
(a)L1拉索损伤
R C O F /%
D /
%
n
(b)L2拉索损伤
R C O F /%
D /
%
n
(c)L5拉索损伤
R C O F /%
D /
%
n
(d)L6
拉索损伤
6040200-20-40
R C O F
/%
20
15
105
25
(e)L11索损伤80%㊁L12索损伤40%
6040200-20-40
R C O F /%
2015105
25
(f)L9索损伤80%㊁L12索损伤80%
图5㊀神经网络部分输入向量
Fig.5㊀Partialinputvectorofneuralnetwork
输入层
隐含层
输出层
输出层
图6㊀神经网络的建立Fig.6㊀Establishedneuralnetwork
473 南昌大学学报(工科版)2023年㊀

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。