文章编号： (200)0-000-0

Size: px

Start display at page:

Download "文章编号： (200*)0*-000*-0*"

征任
7 years ago
Views:

1 偶页湖南大学学报 ( 自然科学版 ) Journal of Hunan University (Natural Sciences) 2016/3/7. 文章编号 : (200*)0*-000*-0* 大跨度拱桥静动力试验与结构识别的实践周云 1,2 张军凯 1 陈松柏 1 易伟建 1,2 (1. 湖南大学土木工程学院, 湖南省长沙市 ; 2. 工程结构损伤诊断湖南省重点实验室, 湖南大学, 湖南省长沙市 ) 摘要 : 由于有限元建模中存在不确定性因素的影响, 大跨度拱桥初始设计模型的理论分析结果与现场测试数据之间存在一定差异为了消除模型校验技术中的认知不确定性和随机不确定性, 从而更好地进行结构识别, 论文以来华大桥和巴溪洲大桥为例, 对其中若干关键问题进行了讨论通过合理地现场勘察图纸资料和与精确的有限元模拟相结合, 基于灵敏度分析方法进行不确定性参数筛选, 并最终采用最小二乘法进行参数优化, 实现了可靠的结构工作状况评估与响应预测研究包括现场静动力试验有限元模拟模型校验响应预测和校验合理性评估等基于精确建模和基于模态参数的校验技术可以有效地减小有限元模型与现场静动力试验结果之间的误差, 使用未经试验结果校验的设计模型, 进行两座桥梁校验系数的评估分别存在着 21% 和 23% 的误差, 从而将导致桥梁的成桥评估的技术参数不准确关键词 : 大跨度拱桥 ; 静动力试验 ; 工作模态分析 ; 结构识别 ; 模型校验中图分类号 :U446.1 文献标识码 :A Practice of Static and Dynamic Experiments and Structural Identification of Long Span Arch Bridges Zhou Yun 1,2 Zhang Junkai 1 Chen Songbai 1 Yi Weijian 1,2 (1.Hunan University, Department of civil engineering, Changsha, Hunan province, , China China) 2.Hunan Provincial Key Laboratory of Damage Detection, Hunan University, Changsha, Hunan province, , Abstract: Various uncertainties existing in the FE modeling may have significant impact on structural identification and could result in obvious discrepancy between FE model analysis and field test results. In order to eliminate epistemic uncertainty and aleatory uncertainty during model calibration procedure, which contributes to achieve reliable condition assessment and response prediction, two long-span bridges were taken for example, and several critical issues were discussed. By virtue of field inspection, review of blueprints and precise modeling, parameter calibration were carried out by sensitivity-based uncertain parameter selection combined with least square method. Research in this paper is consisted of field tests, precise finite element simulation, model calibration and model admissibility check. Discrepancies between FE model and field test results can be mitigated by efficient model calibration method, and initial design model generally shows 20 percent relative difference considering the calibration coefficients. * 收稿日期 : 基金项目 : 国家自然科学基金青年基金 ( ), 高校博士点基金新教师基金项目资助 ( ), 湖南省优秀青年骨干教师培养计划作者简介 : 周云 (1979 ), 男, 湖南长沙人, 湖南大学副教授, 博士通讯联系人, jkzhang1119@126.com

2 Key words: Long-span bridge; Static and dynamic test; Operational modal analysis; Structural identification; Model calibration; 解析解和实验数据之间差异建立目标函数, 通过目引言标函数的优化实现参数识别由于不确定性参数和模态参数之间的非线性关系, 参数修正通常采用迭代算法, 识别结果可以保证模型单元的连续性和系统矩阵的物理意义 [15][16] 本文以两座新建的大跨度拱桥为研究对象, 依据工作模态分析测试结果和有限元计算模态分析结果建立目标函数通过模型校验技术, 最终得到能够真实准确反映结构状况的有限元校验模型静载试验结果和校验系数表明, 经过现场观测环境振动试验和模型校验, 校验模型能够更准确地模拟实际结构, 其响应预测结果更加准确可靠在大跨度桥梁结构中, 拱桥具有跨越能力大自重轻外形美观等优点 [1] 然而, 由于拱桥是以推力为主的结构, 其施工质量难以控制, 且构件的性能变化难以观测 [2] 大跨度拱桥的构件参数存在一定的不确定性, 对其进行状态评估的难度较大通过结构识别技术确定其真实工作状况, 对桥梁的使用维护评估以及决策支持具有重要意义随着结构识别理论在工程领域的发展, 国内外许多学者对大跨度拱桥的静动力试验和结构识别进行了大量的研究工作 [3] 1997 年,Nazmy A S 等探讨了不同设计参数下钢拱桥结构性能的变化, 分析参数变化对拱桥的强度和稳定性的影响, 并提出了拱桥的结构设计优化 [4] 建议 2004 年,Ren W X 等对田纳西河桥进行了环境振动测试, 并分别使用峰值提取法 (PP) 和随机子空间识别法 (SSI) 进行了参数识别 2006 年, [5] Yong X 等对比进行了 MRWA bridge no 使用碳纤维复合材料加固前后的模态测试, 在工作模态分析的基础上通过模型修正对结构工作状况进行了评估 2008 年,Filipe M [6] 研究了一种多通道动态监测系统, 该系统能够实现桥梁的固有频率与频率域分解方法的在线自动识别 2008 年,Song Y 等 [7] 使用密布的无线传感器完成了厦门五缘大桥的模 [8] 态试验和工作模态分析 2011 年,Wei L 等对一座 5 跨连续拱桥进行了环境振动测试, 结合车桥耦合振动试验以及静载测试的结果, 对桥梁的健康状 [9] 况进行了评估 2011 年,Qin P 等研究了结构识别中各种不确定性因素的来源, 并且通过对 Henry Hudson 桥的结构识别和模型校验实例说明了如何减小不确定性因素对识别结果的影响 2013 年, 李晰 [10] 等利用设置在一座钢管混凝土拱桥上的健康监测系统, 对实际结构的工作模态进行识别并采用零阶近似法对初始有限元模型进行了修正 2014 年, 胡 [11] 志坚等针对九江长江大桥主桥三联拱桥部分开展大缩尺比的静动力模型试验研究, 通过各类构件不同的损伤模拟方法研究了结构动力特性的变化 [12] 2014 年,Costa B 等对比了路易斯一世桥加固前后随机振动试验的模态分析结果结合有限元模型结构识别, 论证了桥梁刚度模态参数和边界条件之间的关系 2014 年, 陈清军以一座经历地震后的框架结构为研究对象, 通过环境振动测试, 采用小波变换的方法进行了损伤识别研究 2015 年, 周云等通过对一块钢 - 混凝土组合板进行静动力试验, 通过基于脉冲锤击测试的模态柔度理论进行了结构损伤的识别基于灵敏度分析的模型校验是依据有限元模型 1 桥梁结构识别结构识别理论的提出始于上世纪 60 年代末, 发展至今已有 50 多年的历史结构识别的范式最早于 1977 年被 Hart G 和 Yao J [17] 提出自 1978 年,Liu S 和 Yao J [18] 将其引入到土木工程领域后, 结构识别成为了工程力学界和土木工程界的一个富有活力的研究方向 2005 年, 结构识别的概念被 ASCE St-Id of Constructed Systems Committee 所接受, 其范式被描述为理论 - 试验 - 决策的六步综合圆 [19], 如图 1 所示图 1 结构识别六步圆 Fig.1 Six steps of Structural Identification circle 为保证研究工作的缜密和完善, 本文的大跨度桥梁结构识别严格按照六步圆的顺序开展具体工作内容划分如下 : (1) 观察和概念详细查阅结构的图纸记录, 并结合现场观测结果确定结构的建筑材料结构形式等, 明确结构当前状况和可能出现的损伤 (2) 先验模型基于己知信息建立初始的有限元模型先验模型相对于实际结构的准确程度可以通过对比其预测结果和现场测试数据来判定通常可取用设计单位出具的设计模型作为先验模型 (3) 控制试验根据校验目标选定合适的测试方案对于初始建成的大跨度拱桥, 在通车之前对其进行静载试验和随机振动模态试验是有效评定其

3 工作状况的方法之一静载卡车试验按照预先设定的加载工况使用卡车施加荷载, 待加载稳定, 观测桥面和拱肋控制截面的应力与变形来判断桥梁结构的受力性能是否满足要求随机振动试验以环境振动作为激励, 通过安放于结构上的单轴加速度传感器拾取结构振动的时域信号试验测点分为参考点与移动测点, 并假定参考点输入信号为高斯白噪声 [18] (4) 数据分析与阐述静载试验主要测试桥梁控制截面的应变应力和挠度, 通过结构在外荷载作用下的反应来评估结构能否达到正常使用的要求模态试验中, 根据各测点的加速度时程信号, 经由快速傅立叶变换转换为频域信号分别计算参考点的自功率谱以及参考点关于各移动测点的互功率谱, 最终由 Poly-Max 方法求解结构的模态频率阻尼与振型 (5) 模型校验与参数识别, 也称模型试验匹配在结构识别的第 5 步模型校验前, 首先要判断有限元模型的准确性, 在精确建模并消除有限元模型不合理的模拟方式引入的误差的基础上进行参数识别 [20] 基于灵敏度分析的参数校验有助于明确对结构响应较为敏感的建模不确定参数并通过目标函数迭代寻优的方法实现参数校验 [21] 其中, 目标函数和不确定参数的选取决定了修正的方向, 本文选取了基于模态频率误差的目标函数 : n fai ( x) fei ( x) F ( x, i) = 100% (1) f ( x) i= 1 f ai : 有限元模型计算所得到的模态频率 ; f ei : 实际测试分析得到的模态频率 ; x: 有限元模型中的建模不确定参数 ; i: 所选取的参与目标函数的第 i 阶模态其中, 试验模态分析与理论模态分析通过模态置信准则 (MAC) 匹配 ( 如公式 (2) 所示 ) MAC 值接近于 1 说明两种振型匹配良好 T 2 ( ϕaiφei ) MAC i = (2) T T ( ϕ ϕ )( φ φ ) ai 其中, φ ai 和 φ ei 分别表示有限元理论模态分析和环境振动测试工作模态分析得到的振型向量 ;i 指代进行分析的第 i 阶模态 (6) 结构响应预测在结构识别第 6 步, 校验系数 η [22] 是评定结构工作状态和变形性能的一个重要指标, 可以由其判定有限元模型的准确程度 η 由静载试验数据 ( 应力挠度 ) 与有限元模型预测结果 ( 应力挠度 ) 的比值得到, 如公式 (3) 所示, De η = (3) Dr 其中,D e,d r 分别表示有限元模型预测的理论挠度与实际结构现场测试挠度通常,η 1 说明有限元计算结果偏于安全, 实际结构工作状况良好并尚有一定的安全储备 η=1 则说明有限元计算结果与试验 ai ai ei ei 结果十分接近, 模型较为准确地模拟了真实结构 2 例 1: 来华大桥结构识别的研究 2.1 桥梁概况广西来华大桥位于广西省来宾市, 是一座中承式钢管混凝土拱桥大桥全长 465m, 桥面宽 36m, 主跨拱肋为双肋悬链线无铰拱, 计算跨径为 210m, 计算拱高 60m, 矢跨比为 1/3.5 每片拱肋由 4 根直径 750mm, 厚度为 20mm 的 Q345C 钢管组成, 内部灌注 C50 微膨胀混凝土拱肋之间由水平联通钢管和斜腹钢管连接, 整体结构形成桁架式拱肋桥面板与钢管混凝土拱由立柱和吊杆连接, 桥面上下游方向各有 8 根立柱,21 根吊杆, 实际结构照片见图 2 湖南大学土木工程结构健康监测研究团队 ( 于 2013 年 6 月对该桥进行了系统的静动力试验图 2 来华大桥照片 Fig.2 Picture of Laihua Bridge 2.2 静载试验在各种不同的荷载组合作用下, 钢管混凝土拱和主跨桥面竖向变形及应变是本次静载试验的主要内容加载的位置主要集中在主跨 1/2 位置和主跨 1/4 位置按照不同的车辆数分别设置沿桥梁纵轴对称加载和偏心加载的两种加载形式, 每台加载车辆的重量均为 300kN 来华大桥静载试验共预设了 20 个加载工况, 其中主跨 1/4 位置加载 10 台车的工况如图 3 所示北岸正载 1/4 跨 10 台车 220m 南岸图 3 静载试验加载布置图 Fig.3 Overview of static loading test 桥面测点布置在沿桥面均匀分布的 8 等分点上, 位于桥面防撞护栏内侧平坦处测试过程采用电子水准仪和光学人工水准仪分别进行测量读数, 并相互校核测试结果在进行钢管拱的静载卡车试验前,

4 预先将全站仪棱镜安放在钢管拱上各个 8 等分点位置, 分别在上下游河滩安置全站仪测量各个测点的位移结果 2.3 模态试验本文使用 LMS Cadax-8 系统进行数据采集工作,3 个参考点 ( 第 4# 11# 和 16# 测点, 可见图 4) 的传感器固定布置, 每次测试移动其余 5 个传感器所测试的位置测试过程首先将传感器竖直放置拾取竖向的加速度时程响应, 然后将传感器水平放置, 完成桥梁横向振动测试, 数据采集时间为 10~20 分钟, 采样频率为 512Hz, 由随机子空间识别法获得的桥面竖向振动稳态图可见于图 5 南岸 220 北岸使用两节点的 beam 单元模拟 ; 拉索的模拟使用仅承受拉力的三维 cable 单元 ; 桥面板人行道板和人行道扶手由 plate 单元模拟桥面板两端与引桥部分通过伸缩缝隔开, 其边界条件在有限元模型中被模拟为理想铰支座 ( 限制相应节点的纵横向和竖向平动位移 ) 钢拱的两端边界条件假定为固定约束 ( 限制相应位置节点的全部平动及转动位移 ) 图 6 初始设计模型 Fig.6 Initial design model DAQ DAQ 图 4 来华大桥模态测试测点布置图 Fig.4 Layout of sensors in the ambient vibration test 5.4 e-8 g 2 Log do sf df vo svd fs d s ds ds s sd d s f s s ss d s d s f s s s d s s d f ds sd s s s d d s sdsf 7 fo vf oo fo do f fd df df ff f f f o d o f f f f fd f f d o d f f d f f f o f f f f d f f o f f f f 7 d os o f d df sf d fs sd df ds fs f s f s o fd f f os d s o d s d d sd sd f f ds dd s d s f s s ss 7 s s d sd ds s so ds s ss s ds s s s s s ss ss s f s s d s s d dd s s d s s sd d d s fd d dds 7 o f d so fd f f o f f sf d ds d f d d o fd ff f o df d d f d o ff d f o d df f f d d d fd f d f f 7 d s d s vs sd o ds sd osf d ds d s s s ds d dd s o sd s d s d f d d d d df d dd d d d fd f d ff 7 d s d s s sd o ds sd s s ds s s s s o ds o df s d s d f d o f f f d f f d df d f s ff f d fd 7 d s o d fs sfd fds sd sd ds ds d d d do ds o ff f df d d f df f o f d ff df d f d f d f f df f 7 o s o o fd sfd fd o sf df d fs f d f s odf f d o df d d d d f dd d fd d dd d f d ff f f f s 7 d s d s ds sfd ds d s ds ds s s d f s sdf s f d o df d f d f f o d f f d f df f f d f f s f f d 6 o s f f ofd sff ds sd sd fd ds d d f o s osd f f f fd s f d f f fff f f df f d d f d f fd f 6 d s d s vd svd s f s s ds sd s d s o ss s s s ds d f d d d f f dd d sd d d d d s f d d d ff 6 d s o v s svs ds fs s s ds sf s s s o ss d f d ds s f s d df d ff f ds f f f f s o f f f df 6 d s of fs svs s s s s ds s o s d s o ss vd s ds s d s d f ddf d fds d d d d d d f d d d 6 d s f vd sff ds f sd s fs d f s f s o ds d f d fs d f f d d f ff f f d d f f d d d f f d 6 d s o od sf fs fd s ds fs f o s f s f ss d f d fs d f f f f f f o f odf f d f f f d d f ff 6 d s d d o sfd f df sd d fs f fd f d o f f o d fd d f d f d f f f f ds f d f d d f d d fd 6 o s f fs sff f ds sd sf d ds d dd fd fs f f f fd s o f d d df o f df f f f f f f f ff f 6 o f v fd fd f fd d sf d ds f s f od fd o f d o f s o f d f f f f f ff f f f f f f f f d 6 s d fd dffoo fd d sf d dv f d fo d ds ff f fd d d f f f f f f o d f f f f d f d f d 5 f d d dfd dfo d sf d d d d f o v dd ff d ff s d d o f df ff df od d f f f f f f 5 s d fd sfd fs df sd d fs f s f o s dd d f d f d d d d f f f f df f f d f d f d ff f 5 s s fs svd ds o sd sd d ds s s f f s ds f d d fd s d d f f f d f df d f f d f f f f f 5 s o df sfd fd o sd df d ds f s f d s dd f f f ds d d f f f f ff fff d d f f o f f f 5 s f f dfd fv f fdf d f f d f o d ff f o f f f f f f o ff o ff f d f f f f f f d 5 v d df dvf d fd sd d fd d fd f d fd f f fd d d f f d f f fd f f f d f f f f 5 s d df svd fd sd d d fs ff d f d o ff o f f o f d f f f f f f f f f f f f f f d 5 s d df sfd ds sf sf d ds df s f s o ds d df s d f f d f d fff d f f f f f f d 5 s f sd svf df df df df ds do s f d o df f df d d f f f f d f fd d f f d f f f 5 v d df sf f df f d d o df f d f d o df f f d f f o o f f f ff f d o f f f f 4 s d fd sf d fs sf dd ds ds ff df s fds d d f d f f o f d f o ff f f f o f f 4 s s sf svd ds f ss s ds ds sd s df f f fd d o o f f f dff f f d f f 4 v o fd sf d fd f df d fs ff ff d dff f f o f d o ff f f dd f f f f f f 4 s d df sf d d d df ds ds d d d o df d d d o f o f f d fd df f f f f 4 f d df sf d df d sd d o ds s dd d fd f f d d f f o f f f ff f o f f f 4 f o d ff o f df df df fs d d d od f s f f d o d f f ff f o f o f 4 f f fd df d f fd sd d fs f ff d f f f f f d f o f f f fd f o f f f e Linear 10.0 Hz 图 5 基于 SSI 方法获得桥面振动的稳态图 Fig.5 SSI Stabilization diagram for the deck 2.4 有限元建模及分析根据来华大桥在设计和成桥试验的两个阶段的研究目的, 笔者分别使用有限元软件 Midas Civil 和 Strand7 建立了来华大桥的空间有限元模型 ( 后文分别以设计模型和校验模型指代 ), 从模型校验前后所得到的不同结果可以预测实际结构和计算模型在性能方面的差异来华大桥设计模型通过 Midas 建立 ( 图 6), 根据其进行全桥施工过程控制, 各个施工阶段的线型和应力都在合理的范围内主要控制参数能直观地反映实际工作状况, 相应的有限元计算结果可以作为成桥静动力试验的有力参考在完成静载控制试验和模态测试后, 笔者使用 Strand7 建立了来华大桥校验模型 ( 见于图 7) 所有单元的几何特征和物理参数都详细按照现场勘测结果确定拱肋纵梁和横梁 K 型横撑立柱等均图 7 经过校验的有限元模型 Fig.7 Calibrated finite element model 灵敏度分析反映模型输出的分析结果对于各建模参数的敏感程度通常, 在有限元模型中选取主要单元的弹性模量几何尺寸以及边界条件为不确定性参数不确定参数是模型校验中的自变量, 校验结果通过改变参数的数值实现, 每个不确定性参数的上下限定义了修正参数的变化范围, 即模型校验的约束条件校验结果在缩小模型与真实结构之间误差的前提下, 首先应保证校验后的不确定性参数具备现实物理意义 [23] 另外参数变化范围也根据工程经验施工工艺和相关参考文献的常用取值范围选取 [24], 本文对每个不确定性参数选取的合理约束界限统计于表 1 来华大桥部分不确定性参数的基于模态参数的灵敏度分析结果可见于图 8 其中, 横轴为不确定性参数的标准化弹性模量, 竖轴为基于 7 阶模态频率的目标函数值, 见公式 (2) 根据灵敏度分析结果, 选取 4 个对模态分析结果有显著影响的关键不确定性参数为目标函数中待修正的自变量, 将工作模态分析结果最可靠的前 7 阶模态频率选为修正参数 ( 即 n=7), 通过调整不确定参数的取值寻找有限元模态分析结果与工作模态分析结果前 7 阶频率相对误差的迭代寻优目标函数迭代优化过程中, 通过公式 (3) 实现有限元计算模态与试验模态振型的匹配, 并定义自变量的约束条件为 : l x (4) b u b 其中,x, l b, u b 分别表示不确定参数变量下限和变量上限本文采用有限元软件 Strand7 的 API 模块实现其与 Matlab 的交互访问, 并使用最小二乘

5 迭代优化算法 Isqnonlin 进行校验, 收敛准则为 : 3 abs ( f i + 1 f i ) 或 50 iter (4) abs(f i+1 -f i ) 与 iter 分别表示相邻两次迭代运算目标函数的差值和当前迭代次数, 所有参数的校验结果统计于表目标函数值标准化弹性模量图.8 灵敏度分析曲线图 Fig.8 Sensitivity analysis of modulus 表 1 基于模态频率的灵敏度分析所选取参数及其校验结果 Tab.1 Sensitivity analysis and parameters selections based on modal frequency 参数参数设计值变量下界变量上界校验结果混凝土拱弹性模量 /MPa E E 钢拱弹性模量 /MPa E E 人行道板密度 /(kg/m 3 ) 支座竖向刚度 /(N/mm) 校验结果合理性评估本文选取加载挠曲变形最大的主跨 1/4 位置加载 10 台车的工况, 对现场试验结果设计模型和校验模型预测的结果进行对比桥梁钢拱和桥面板的竖向变形可见于图 9~ 图 10 从静载变形曲线可以看出, 两种模型均能正确反映大跨度拱桥在荷载作用下的变形状况, 对比于设计模型, 校验模型的变形结果与试验测试数据更为一致设计模型中各个测点的校验系数值较小, 说明由其模拟的桥梁结构具有较高的安全储备校验模型的结果则分布集中, 所有校验系数都稳定于 0.9~1.0 之间, 模型评估结果较为准确可靠 ( 见于表 2) 此外, 来华大桥的工作模态分析和计算模态分析结果列于表 3 对比两种模型的模态分析结果, 来华大桥校验模型的理论模态分析与现场试验工作模态分析的频率误差相较于设计模型大幅减小位移 ( m) 位移 ( m) 实际测量设计模型校验模型测点沿桥纵向坐标 (m) 图 9 1/4 跨加载 10 台车作用下的钢拱变形 Fig.9 Displacement shape of steel arch under 10 cars loaded on the 1/4 span 实际测量设计模型校验模型测点沿桥纵向坐标 (m) 图 10 1/4 跨加载 10 台车作用下的桥面变形 Fig.10 Displacement shape of deck under 10 cars loaded on the 1/4 span 表 2 挠度校验系数对比 Tab.2 Comparison of deflection calibration coefficient 位置拱肋桥面主跨 1/4 主跨 3/4 主跨 1/4 主跨 3/4 实测挠度 /mm 设计模型挠度 /mm 校验模型挠度 /mm 校验前 η 校验后 η 表 3 前三阶有限元软件计算模态振型与实测结果对比 Tab.3 Comparison of modal analysis between ambient vibration testing and 3D FE model Mode 1 Mode 2 Mode 3 实测模态 f 1=0.69Hz f 2=1.04Hz f 3=1.53Hz

6 附表 3 前三阶有限元软件计算模态振型与实测结果对比 Tab.3 Comparison of modal analysis between ambient vibration testing and 3D FE model Mode 1 Mode 2 Mode 3 设计模型模态 f 1=0.46Hz f 2=0.93Hz f 3=1.28Hz 校验模型模态 f 1=0.70Hz f 2=1.07Hz f 3=1.42Hz 3 例 2: 巴溪州大桥结构识别的研究 3.1 桥梁概述巴溪州景观桥位于长沙市坪塘地区, 全桥长 280m, 主跨部分为 180m 拱肋结构为全焊提篮式钢箱拱, 钢拱肋的拱轴线在空间上为悬链线形式, 拱轴系数 1.8, 跨径 142m, 矢高 24.89m, 拱肋截面内倾角为 8 箱型截面尺寸为 1.0m 1.8m, 两片拱肋通过 5 道钢箱横撑连为整体拱肋沿纵向分三种节段类型 : 钢筋混凝土结合段钢拱肋 A 段和钢拱肋 B 段吊索直接承受来自桥面主梁的荷载, 间距为 6.25m, 沿桥梁纵向布设, 全桥共设 23 对, 桥梁照片可见于图 11 正载跨中 4 台车图 12 静载试验加载布置图 Fig.12 Overview of static loading test 3.3 模态试验为准确模拟桥梁实际工作状况, 巴溪洲桥模态测试通过环境振动 ( 风水波冲击行人荷载和地脉动等 ) 激励结构的响应工作模态分析结果通过随机子空间识别法提取, 测点布设如图 13 所示图 11 巴溪州桥照片 Fig.11 Picture of Ba Xi Zhou Bridge 3.2 静载试验静载试验记录了巴溪州桥在载重卡车荷载作用下各个控制截面的应力应变等数据, 全面反映了结构的工作状况试验共设置 14 个工况, 分为桥梁纵向 1/2 位置和 1/4 位置加载加载车辆平均重量为 270kN, 加载稳定后通过全站仪读取桥面和钢拱的变形在后续的校验模型评估阶段, 本文选取了桥面挠曲变形最大的 1/2 位置加载 4 台车工况, 如图 12 所示图 13 模态测试测点布置图 Fig.13 Instrumentation layout in modal test 其中, 第 11# 点和 44# 点被设为参考点, 其余为移动测点, 数据采集工作由 8 通道 LMS Cadax 系统完成测试过程先后测取桥梁的竖向模态以及横向弯曲模态, 每次采集时间为 15 分钟, 采样频率为 256Hz 竖向模态的测试中, 经过信号处理和频谱分析所得到的模态提取图如图 14 所示

7 4.3 e-5 g 2 Log 1.8 e Linear 10.0 Hz 图 14 SSI 方法计算桥面模态稳态图 Fig.14 SSI Stabilization diagram of the deck 3.4 有限元建模及分析巴溪州桥的设计模型通过有限元软件 Midas 建立 ( 如图 15 所示 ), 所有构件都使用空间梁单元模拟, 设计模型采用的建模方法基本满足施工控制和结构工作状况评估的要求校验模型的建立综合采用了梁单元桁架单元 [25] 实体单元等桥面系通过梁格法建立, 利用等效的梁格代替桥梁上部结构, 将分散在梁板每一区段内的弯曲刚度和扭转刚度集中于最邻近的等效梁格内, 每一梁格内的弯矩剪力和扭矩均等于该梁格所代表的实际结构的内力校验模型中, 拱肋两端的大体积混凝土由实体单元模拟 ; 吊索采用仅承受拉力的桁架单元模拟 ; 刚性吊杆采用梁单元模拟, 并且释放梁端的约束桥头支座被模拟为简单支撑, 支座节点沿 x 轴 z 轴的转动位移和沿 x y z 三个方向平动位移均被限制, 模型示意可见于图 16 图 15 巴溪州桥设计模型 Fig.15 Design model of Ba Xi Zhou Bridge 图 16 巴溪州桥校验模型 Fig.16 Calibrated model of Ba Xi Zhou Bridge 对于校验模型, 根据环境振动测试的工作模态分析和有限元模型理论模态分析的频率相对误差建立目标函数根据工程经验选取部分不确定性较高的灵敏度参数及其相应的变化范围 C40 混凝土容重的灵敏度分析曲线如图 17 所示其中, 横轴表示混凝土容重相对于初始设计值的变化过程,E 为初始设计值, 竖轴为前 3 阶模态频率的相对误差之和参数校验基于最小二乘法原理, 所选取的 3 个待修正的不确定性参数可见于表 4 目标函数选取如公式 (2) 所示其中,n=3, 即模型校验的目标是实现有限元模态分析前 3 阶模态频率与工作模态分析结果的契合修正过程的自变量个数为 3, 每次修正均固定其余两个自变量, 在当前修正自变量的界限内寻找满足目标函数极小值的最优点, 并将其确定为校验结果 14 频率相对误差 (%) E 1.10E 混凝土容重图.17 混凝土容重灵敏度分析 Fig.17 Sensitivity analysis of density of the concrete 表 4 基于模态频率的灵敏度分析所选取参数及其校验结果 Tab.4 Sensitivity analysis and parameters selections on modal frequency 参数参数设计值变量下界变量上界校验结果钢拱弹性模量 /MPa E E C40 混凝土容重 /(kn/m 3 ) 桥面铺装重量 /(kn/m) 校验结果合理性评估为了评估模型校验结果的合理性, 本文选取加载变形最大的主跨 1/2 位置加载 4 台车工况下结构的变形进行研究钢拱测点布置在其四等分点上, 桥面板的边跨测点布设在其跨中, 主跨测点布置在所有的八等分点上, 挠曲变形数据结果统计于图 18~ 图 19 对比试验和两种有限元模型的计算结果, 钢拱和桥面测点的变形与模型预测结果一致, 且大部分测点位移小于相应的有限元解析解为了更加清晰地反映巴溪洲大桥在加载时的工作状况, 部分关键位置的校验系数列于表 5 对比两种模型的校验系数, 设计模型能够完成结构在外部荷载作用下承载变形性能的初步预测, 而经过精确的现场勘验和单元模拟, 以及由模态参数目标函数完成参数校验的校验模型对结构工作状况的评估和响应预测更为准确可靠, 其控制测点的校验系数均高于设计模型的相应结果

8 Fig.19 Displacement shape of deck under 4 cars loaded on the mid-span 巴溪洲桥的工作模态分析和计算模态分析结果列于表 6 相比于设计模型, 校验模型的静载变形结果和模态分析结果均更加符合巴溪洲桥的现场试验数据, 真实全面地模拟了实际结构, 验证了有限元模型校验的合理性图 18 跨中加载 4 台车作用下的钢拱变形 Fig.18 Displacement shape of steel arch under 4 cars loaded on the mid-span 表 5 挠度校验系数对比 Tab.5 Comparison of deflection calibration coefficient 拱肋桥面位置拱跨 4/8 拱跨 6/8 主跨 4/8 主跨 6/8 实测挠度 /mm 设计模型挠度 /mm 校验模型挠度 /mm 校验前 η 校验后 η 图 19 跨中加载 4 台车作用下的桥面变形表 6 前三阶有限元软件计算模态振型与实测结果对比 Tab.6 Comparison of modal analysis between ambient vibration testing and 3D FE model Mode 1 Mode 2 Mode 3 实测模态 f 1=0.61Hz f 2=1.00Hz f 3=1.83Hz 设计模型模态 f 1=0.58Hz f 2=0.98Hz f 3=1.71Hz 校验模型模态 f 1=0.60Hz f 2=0.99Hz f 3=1.81Hz 4 结论大跨度拱桥是以推力为主的结构, 其施工质量控制难度大, 且构件的性能变化难以观测构件的物理参数存在一定的不确定性, 传统方法难以实现准确的状态评估本文以团队测试过的两座大跨度拱桥为研究对象, 按照初始设计模型现场试验模型校验校验结果评估等步骤完成了两座大跨度拱桥的结构识别 (1) 大跨度拱桥结构中构件与材料性质复杂, 使用有限元方法对其模拟的过程引入了大量的物理参数和几何参数不确定性通过对比两座大桥的初始设计模型和其工作模态分析的结果, 使用仅基于

9 设计图纸建立的设计模型进行两座桥梁校验系数的评估分别存在着 21% 和 23% 的误差, 设计模型仅能作为实际结构的先验与近似 (2) 综合考虑大跨度拱桥模态频率与模态振型匹配的目标函数配合最小二乘优化算法可以有效地减小有限元模型与真实结构的参数误差其中, 灵敏度分析可以迅速筛选重要的不确定性参数完成参数的迭代优化计算后, 来华大桥与巴溪洲桥有限元模型的计算模态分析误差分别缩小了 16.24% 和 3.25% (3) 基于现场勘测的真实结构细节模拟和基于模态参数目标函数的重要不确定性参数校验过程, 大幅降低了初始有限元模型中的几何参数与物理参数不确定性校验模型的变形数据和现场静载试验结果更为一致, 为进行准确可靠的大跨度拱桥工作状况评估和响应预测奠定了基础参考文献 [1] 钟善桐. 钢管混凝土结构 [M]. 北京 : 清华大学出版社, 2003: ZHONG Shan-tong. The concrete-filled steel tubular bridge structures [M]. Beijing: Tsinghua University Press, 2003: (In Chinese) [2] 张跃红. 钢管混凝土拱桥的损伤识别研究 [D]. 大连 : 大连理工大学土木工程学院, 2006:1-3. ZHANG Yue-hong. Study on damage detection of the concrete-filled steel tube arch bridge [D]. Dalian: Dalian University of Technology, 2006:1-3. (In Chinese) [3] Nazmy A S. Stability and load-carrying capacity of three-dimensional long-span steel arch bridges [J]. Journal of Computers and Structures, 1997,65(6): [4] Ren W, Zhao T, and Harik I. Experimental and analytical modal analysis of steel arch bridge [J]. Journal of Structural Engineering, 2004, 130(7): [5] Giovanna Z, Hong H, Yong X, et al. Stiffness assessment through modal analysis of an RC slab bridge before and after strengthening [J]. Journal of Bridge Engineering, 2006, 11(5): [6] Magalhaes F, Cunha A, Caetano E. Dynamic monitoring of a long span arch bridge [J]. Engineering Structures, 2008, 30(11): [7] Song Y. and Lei Y. Application of wireless sensing system for the wuyuan steel arch bridge [C]// Earth and Space, Long Beach, California, United States:ISBN,2008:1-6. [8] Wei L. and He Q. Health state evaluation for multi-span double-curvature arch bridge based on static-dynamic load testing [C]// Condition, Reliability, and Resilience Assessment of Tunnels and Bridges.Zhangjiajie,China,GeoHunan International Conference 2011: Curran Associates Inc,2011: [9] Qin P, Grimmelsman K, Moon F, et al, Mitigating epistemic uncertainty in structural identification:case study for a long-span steel arch bridge [J]. Journal of Structural Engineering,2011,137(1):1-13. [10] 李晰, 张德义, 闰维明, 等. 基于环境激励的钢管混凝土拱桥工作模态识别及修正 [J]. 工程力学,2013,30(9): LI Xi, ZHANG De-yi, YAN Wei-ming, et al. Operational Model identification and calibration of cfst arch bridge based on ambient excitation [J]. Engineering Mechanics,2013,30(9): (In Chinese) [11] 胡志坚, 张明辉, 孔祥韶, 等. 大跨度钢桁架拱桥静动力相似模型 [J]. 中国公路学报,2014,9(27): HU Zhi-jian, ZHANG Ming-hui, KONG Xiang-shao, et al. Static and dynamic experimental model for long-span steel-truss arch bridges[j]. China Journal Of Highway And Transport, 2014, 9(27): (In Chinese) [12] Costa B, Magalhães F, Cunha Á, et al. Modal analysis for the rehabilitation assessment of the Luiz I bridge [J]. Journal of Bridge Engineering.2014,19(12) [13] 陈清军, 周成杰, 杨永胜. 基于环境振动信号的框架结构震后损伤识别 [J]. 湖南大学学报 : 自然科学版,2014,41(9): CHEN Qing-jun, ZHOU Cheng-jie, YANG Yong-sheng. Damage identification of frame structure after earthquakes based on environmental vibration records [J]. Journal of Hunan University: Natural Sciences. 2014,41(9): (In Chinese) [14] 周云, 蒋运忠, 易伟建, 等. 基于模态柔度理论的结构损伤诊断试验研究 [J]. 湖南大学学报 : 自然科学版,2015,42(5): ZHOU yun, JIANG Yun-zhong, YI Wei-jian, et al. Experimental research on structural damage detection based on modal flexibility theory [J]. Journal of Hunan University: Natural Sciences. 2015,42(5): (In Chinese) [15] Friswell M I, Mottershead J E. Finite element model updating in structural dynamics [M]. Kluwer Academic Publishers, Dordrecht, Netherlands: Kluwer Academic Publishers,1995: [16] Mottershead J E, Friswell M I. Model updating in structural dynamics: a survey [J]. Journal of Sound and Vibration, 1993,167(2): [17] Hart G C, Yao J T P. System identification in structural dynamics [J]. Journal of the Engineering Mechanics Division,1977,103(6): [18] Liu S C. and Yao J T P.Structural identification concept [J]. Journal of Structural Division,1978,104(10): [19] Catbas F N, Kijewski-Correa T, Aktan A E. Structural identification (St-Id) of constructed facilities: Approaches, methods and technologies for effective practice of St-Id [R]. A State-of-the-Art Report by ASCE SEI Committee on Structural Identification of Constructed Systems, Reston, United States, American Society of Civil Engineers,2011: [20] Smith I, Saitta S. Improving Knowledge of Structural System Behavior through Multiple Models [J]. Journal of Structural Engineering. 2008,134(4): [21] Jaishi B, Weixin R. Structural finite element model updating using ambient vibration test results [J]. Journal of Structural Engineering, 2005,131(4): [22] 谌润水, 胡钊芳. 公路桥梁荷载试验与结构评定 [M]. 北京 : 人民交通出版社, 2003:8-12. CHEN Run-shui, HU Zhao-fang. Highway bridge load test [M]. Beijing: China Communication Press, 2003:8-12. (In Chinese) [23] 董小军. 基于灵敏度分析的结构动力模型修正及损伤识别研究 [D]. 兰州 : 兰州理工大学土木工程学院,2009: DONG Xiao-jun. Study on structural dynamic model updating and damage detecting based on sensitivity analysis [D]. Lanzhou: Lanzhou University of Technology,2009: (In Chinese) [24] Aktan E, Çatbaş N, Türer A, et al. Structural identification: analytical aspects [J]. Journal of Structural Engineering, 1998,124(7): [25] ( 英 ) 汉勃利 (Hambly E C) 著 ; 郭文辉译. 桥梁上部构造性能 [M]. 北京 : 人民交通出版社,1982: Hambly E C. Bridge deck behavior [M]. Beijing: China Communication Press, 1982: (In Chinese)

中国科技论文在线中文稿件模板

中国科技论文在线中文稿件模板 5 10 15 20 25 30 35 40 基于静动力试验的大跨度拱桥结构识别与 # 状况评估 * 周云, 张军凯, 陈松柏 ( 湖南大学土木工程学院 ) 摘要 : 大跨度拱桥是连接相邻交通的枢纽, 在公路桥梁现代化建设事业中占据重要地位结构识别是在现场观测先验模型静动力试验和模型校验的基础上进行的结构评估与预测工作, 通过静动力试验校验过的有限元模型能够真实可靠地反映大型结构的工作状况本文以一座大跨度拱桥为研究对象,