您好,欢迎来到划驼旅游。
搜索
您的当前位置:首页悬索桥结构基于敏感性分析的动力有限元模型修正

悬索桥结构基于敏感性分析的动力有限元模型修正

来源:划驼旅游
第33卷第l期 2000年2月 土木工程学报 Vd.33 No.1 Feb 2o0o CHINA CIVIL ENGINEERING JOURNAL 悬索桥结构基于敏感性分析的动力有限元模型修正 ,一、、/ (3 范立础袁万城(同济大学) 张启伟Lt f; f 摘要提出一种悬索桥结构基于特征值敏感性分析的有限元模型修正方法 由结构特征值关于模型参数 的灵敏度矩阵建立理论/试验频率差与模型参数摄动之间的关系式:根据待修正参数的不确定性,将参数摄动限 制干预设的范围,并将参数修正转化为不等式约束优化问题迭代求解。通过一1/|5 ̄比例的悬索桥模型的试验分 析,验证了这一模型修正方法的可行性: 关键词模型修正约束优化敏感性分析悬索桥爱{ 内 结构特征值敏感性分析基础上,提出一种悬索桥结构 1 引 言 有限元模型修正的方法。 对大跨度桥梁进行动力响应预测、振动控制和结 构状态评估,首先必需详细了解结构的动力性能。结 构动力性能基本上取决于其模态特性,如自振频率、 振型和阻尼。这些模态特性可以通过有限元分析得到 其理论值,也可由试验模态分析得到其实测值。但 是,结构自振频率的实测值与有限元理论值之间常常 存在较大的差异。产生理论与试验频率差的原因主要 包括:(1)连续系统离散化引起的误差、(2)结构几 何与边界条件的不确定性、(3)材料特性的变异性以 及(4)试验实测和试验信号处理过程中的误差。其 中前三个原因与结构有限元建模的假定和输人数据有 2理式 2.1 基于敏感性分析的参数修正方法 将一连续结构体系离散为n个自由度的有限元 模型,由下列特征方程求解结构第i阶特征值 和 相应的特征向量 : ^ .= 啦 (1) 其中刚度矩阵置和质量矩阵M由结构参数决定。结 构参数包括结构各构件材料、几何特性、构件联结条 件以及结构的边界条件。 关。如果假定实测模态特性与真实值非常接近,那么 需要解决的问题是如何修正结构的有限元模型,使得 结构模态特性的理论值趋近于实测值。 对于建立的有限元离散模型,假设估计的结构参 数集为P.。通过有限元分析,可得结构理论特征值向 量^。和特征向量矩阵(振型矩阵) (下标a表示有 限元模型初始估计值或理论分析值)。对结构进行模 态试验分析,可得结构实测的特征值(圆频率平方)向 量^ 和振型矩阵 (下标e表示试验实测值)。 由于有限元模型的不精确和试验测量值的误差, 结构模态参数的理论值(^., 。)和试验实测值 (^ , )常常存在一定的差异,这一差异称为理 现有的模型修正方法主要有矩阵优化修正” 、 基于敏感性分析的参数修正 j、特征结构分配算 法 以及利用神经网络原理的修正方法 J 这些模型 修正技术已在简支梁、悬臂桁架等简单结构上得到验 证。但是。对于大型复杂结构(如索支承桥梁),由 于有限元模型中存在大量的材料、几何特性及联结条 件的不确定性,模型修正变得困难而复杂。 悬索桥和斜拉桥的振动特性分析中必需考虑悬 (拉)索轴向拉力对结构总体刚度的贡献。由于加劲 梁质量的变化将通过悬(拉)索轴向拉力的变化而影 响结构刚度,对这类非线性结构有限元模型的修正, 必需考虑结构质量变化对刚度变化的影响。本文在 收稿日期:1999-02-24 国家重大基础性研究攀登计划B(批准号:85-40)及国家自然科 学基金(批准号 59878037)联台资助 论/实测值余量。有限元模型的不精确因素主要来自 三个方面 ]:模型结构误差、模型阶次误差和模型参 数误差。假定对结构构造了一个精细的有限元模型, 并且描述结构运动的物理方程是精确的,那么模型参 数误差便是有限元模型误差的最主要因素。 在桥梁结构模态试验中,振型坐标的测量值常混 有较大的测量噪声而精度不足。如果不考虑振型项而 只计人测量精度高的频率项,通过特征值对结构参数 的Taylor级数展开并略去高阶项,可得结构参数修正 10· 土木工程学报 20O0正 与特征值余量之间的关系: (2) 其中 A 一A (3) :P—P. (4) 其中,aA为结构特征值余量;aP为结构参数的修正 值向量;P为修正后的参数向量; 为结构特征值的 阶灵敏度矩阵。 从式(2)出发,由目标函数 J= P+( 一 ) ( 一 ) (5) 的最小化可得到参数修正值 。式(5)中w 和 w 为加权矩阵。 值得注意的是,当有参数对测量值(频率)影响 非常小,或几个参数对测量值影响相近,那么灵敏度 矩阵 常常会是病态的l8 J。另外,当待修正参数较 多或实铡数据存在较大误差时,修正的参数往往失去 其物理意义。 2.2不等式约束优化 首先,为使修正后参数具有物理意义,引人参数 修正向量 的上、下限值: bJ≤ ≤b (6) 于是,结构参数修正值的求解可以表示为式(5)和 (6)构成的约束优化问题。这一约束优化问题可以转 化为以下二次规划问题来求解: minimize.,( ):去 矸 (7) subjectto ≤d { ) (8) 】 (9) ^=【一,o】 (10) f 0 (“) 0 w 其中,J为单位矩阵。 2.3非线性模型参数修正方法 对于非线性结构理论模型的修正,必需考虑结构 参数对振动的非线性影响。这里提出的方法是:利用 前面提出的不等式约束优化方法,采用多次迭代进行 模型修正。假定每步迭代修正中对结构振动呈非线性 影响的参数发生微小的摄动,这时这些参数对振动的 影响可近似作线性处理。非线性参数的摄动量可以通 过设定 的上下限值或权矩阵 来。图1为 非线性模型参数修正的流程图。 臣豆五匾 亟匿困 叵互五互 1  图1非线性模型参数修正藏程图 参数上、下限的设定影响模型修正的结果,可根 据结构的材料、几何特点结合工程经验加以确定。 加权矩阵w 的选择同样影响修正结果。建议根 据对待修正参数初始估计值的自信程度来构造加权矩 阵1 ,同时考虑使得对结构频率呈非线性影响的参 数在每次迭代修正过程不发生较大的变化。另一加权 矩阵’ 用来强调修正后理论频率与实铡频率在指定 频率点的吻合程度。 迭代修正的收敛准则定义为: ”, i≤tolerance (12) 或 J — J≤t≤ ole… u 其中,岫为实测频率向量:,d^ 为 步迭代后的 理论频率向量。 3实验室悬索桥模型的理论与试验模态分析 3.1悬索桥模型简介 实验室悬索桥模型的原型是主跨888 m的虎门大 桥。桥模根据动力模型模拟准则,以实桥1/150的几 何比例制做。模型总长1025.3 cm,中跨长592.0 cm, 边跨分别为201.3伽l和232.0 cm。箱形桥道梁宽 23.7 cm,纵向主跨与边跨在桥塔处断开。桥塔从承 台顶到塔顶高98.4 crⅡ;两主索间距21.3 cm,索的矢 跨比为1:10.5(参见图2) 第站卷第1期 范立础等-悬索桥结构基于敏感性分析的动力有限元模型修正 ·l1 粗,塔柱与横梁单元问由刚臂联接。桥道粱与吊杆单 元问也用刚臂联接。构造的有限元模型共有222个节 点,208个粱单元和39个刚臂单元,全桥共有1211 T__T_ 。。————。 r 个自由度。 假定桥梁围绕其静力平衡位置作微幅振动,运用 \ = / 兰 (a) 总体布置与测点位置 (b) 桥道粱与传感器布置 线弹性理论进行结构特征值分析。表1(第2列)给 出了有限元分析得到的自振频率(0—20H卫)。从振型 分析可见,低阶侧弯和扭转振动有较大的耦合。图3 绘出了一些振型曲线。 图2虎门大桥模型结构与测点布置(咖) 3.3试验模态分析 模型试验所需测点数量多,而加速度传感器和数 据采集通道有限,因此所有的测点被分成几组并分别 3.2有限元分析 测试。每一组测点测试时取同一个参考点。测点和参 构造的三维有限元模型所有单元均采用三维梁单 元,每个节点6个自由度。桥塔简化为三维多层门式 框架,塔柱下端固结于基底。由于塔柱与横梁相对较 考点布置见图2。振动试验中,借用电扇吹风来模拟 环境激励。表1(第3列)给出了振动模态试验分析 得到的模型桥自振频率。  :、 第一阶侧弯,f-3】7Hz 、 图3有限元理论振型曲线示例 衰1暑素桥模型理论与实测自摄频率 注:V=竖弯;L=捌弯;T=扭弯;PL=桥堵蜊弯;叩口=对称;^叩n:反对称(由于桥梁模型中两侧背索并非完全对称,振型描连中的 对 称 与“反对称”非严格定义) l2· 土木工程学报 20130 对所考虑的各阶振型的影响几乎一致。为避免由此带 4有限元模型修正 4.1模态参数理论与实测结果比较 首先将理论振型和实测振型“配对”。利用MAC 准则来定量地检验实测与理论模态参数的相关性 。 MAC为1或接近1,说明两振型相关性好,或者说是 同一振型。通过振型配对得到相应的理论和实测振型 对(在0—16Hz频率范围),包括8个竖向弯曲、2个 横向弯曲和1个扭转振动振型(见表1)。由表1第4 列中的MAc值可以看出, 这些理论与实测振型对具 有较好的相关性,但理论频率与实测频率具有较大的 差异。 来的病态问题,下面的模型修正中将不考虑对截面积 l4 的修正 表2待修正的不确定参数 4.2待修正结构参数的选择 模型参数误差一般由不精确的材料、几何参数和 联结、边界条件估计引起。对于这里的悬索桥试验模 型,部分结构参数(如各杆件的长度等)可以较精确 地量测得到;粱与吊杆、主缆与塔顶、粱与桥塔的联 结以及结构边界条件也可以较好地控制在理想化的状 态。因而初始有限元模型中的部分几何参数与边界、 联结条件可以认为无需修正 另外,加劲粱、主缆及 吊杆的截面均匀性也由精细的模型制作得到保证,因 而可以假定加劲粱、主缆及吊杆各单元的误差是分别 *主缆索和吊杆的轴向刚度的试验标定值为2 83×104N。为验证 后面参鼓修正的正确性,在初始有限元模型中故意将E A。和 ^h 错误值输人。 迭代收敛的准则设为:理论计算频率与实测频率 误差小于5%,而第一阶频率误差不超过1%。 经过三次迭代,理论,实测频率差迅速收敛到预设 的误差范围。每次迭代后有限元模型计算频率列于表 3。从修正的结果可以看出,首次迭代后,理论,实测频 率误差便缩减至10%的范围内;而经过三次迭代修正 后,理论,实测频率最大误差为4.8%(振型T1),最小 误差为1.0%(振型V1)。迭代各步结构参数的修正量 列于表4。加劲粱弹性模量 和竖弯惯矩L 主缆 致的。最后确定的待修正参数及其初始估计值列于 表2。 4.3模型参数修正 由初始有限元模型参数的特征值灵敏度分析可 知,某些参数(如加劲梁的弹性模量 、质量密度 及截面积A )对所有振型都有影响;而另外一些 和吊杆的轴向刚度E A 与E A 有较大的修正,而加 劲粱的剪切模量G 、桥塔的弹性模量E.和质量密度 参数(如加劲粱截面惯矩,正. 、 , 和 、主缆索轴 向刚度E A 与塔的弹性模量E )则分别对特定的振 型有影响。而且,加劲粱质量密度J"d和截面积A 基本保持不变 值得注意的是,修正后的主缆和吊 杆的轴向刚度与真实值已非常接近。 表3迭代修正模型理论频率与实翻频率比较 *三次造代后有限元理论频率与实测额率误差 第孙卷第1期 范立础等·悬索桥结构基于敏感性分析的动力有限元模型修正 ·13 参数单位同表2 别.同济大学学报,1998,26(5):528~532 5结 论 3 C.F“ and F.M.Hemez.u 1993,31(9):1702—1711 4 K.K.Denoyef and L.D.Pete ̄ing dy ̄mic011y me. ̄ured 唱 ̄,ite dement曲呻 cⅡ M_ib咿 , usirIg址l element-brelemem ̄msiifvity meII10dDI。 针对必需考虑结构几何非线性的悬索桥模型修 正,提出了一种线性化的迭代修正方法,并利用一实验 室悬索桥模型和试验模态分析结果,对建立的悬索桥 动力有限元模型进行了修正。由于以特征值敏感分析 .Me for m da · c flexibility m c∞.删 w , 199"7,35(2):362—368 5 D.C.Zimmealr ̄and M. de r∞.C ̄'Teeting finite element 为基础,因此模型修正过程的中间结果可以清楚地显 示出对结构各阶振动的敏感参数。结果显示,该方法 可以较好的修正初始建立的悬索桥理论模型。使得理 论模态分析结果大大接近于实测值,从而为全面了解 结构的动力性能、确定控制策略以及结构状态评估提 供更精确的理论模型。 参考文献 c直lⅡ ro ̄aels u. ̄iug a 删m曲 西群∞B 眦h垃 曲ent · nique.AⅢ ,1990,28(9)1670~1676 I L嵋i。1g hemal net. 6 M.J.Ats ̄.且nd D.J.Immn.OnⅡ w池M 跏 妇 A ,1998,12(1) 135~161 7 J.E I鼬d且nd M.I. 讪.M0dd up dIlg in咖柚啊l dy· ,1993,167:347—375 Mode/ l— rⅡr1i :A Su ey.J.岛加d d阳 8 M I Ffi ̄eU and J.E.M曲 1e8d.Fin 凡n妇剧删1 A.B m且Il且nd E.J.Nagy. 1pr帆el嘟ll 0f4la 学越 d口痂 Sm ̄turo/功仉删 .Klw ̄er:Ae ̄mie,Pree.1995 叫j啦砌idetlt for modal Ilgtest .删/o ̄rr,o/,1982,21(8):I168~1173 9 11 J..AlIdTmr《且nd D.L.BrowrI.A c 2张启伟,范立础.利用动静力测量数据的桥梁结构损伤识 vector 出 1 LI.fAC}110—116,1982 目 册,Ⅱ BA: D FE M0DEL U )』 Ⅲ G oF A sI 『s_0N BRⅡIGE 厶 hm Wancheng撇(1嘶i University) Abstract Q TNs 衄pI鹧曲t8帅eigenv ̄lue semidvi旷based ifmte e]ement(FE)model uN.ting ap Bacll for suspension d es. 111e FEmodeli8坤d in aIlitemtlvefashion 80 ast0 reducethe difea'enees betweerIthe肌al cal 8Ildtheme ̄uled natural nthefirst order"r ̄ylor series麟pans|m1 ofthe elgenvalue ̄ f 哪|et e8.An eige ̄alue鲫怛jt ty matrix isfirst ot ̄ined Mth re吕pectto B0fne pE吸吡-e 啪t0 be adj峭t耐.A set oflinear equations rda 增the perturbations of p日瑚In芒船惜tothe difer- 朗ce8 between the ̄'Mtieal and the measured呻tu商f 。 曲ci髑are then established.The l ̄.1llmeters to be叩dated are as一  ̄lmed t0 be bounded within 80rrle prescribed re西ons ̄ecording to the degrees of uncertainty and vari ̄iom that ted in the I工毛哪e erB.The changes of these pa皿meter8 a聪found by solving a quadratic matrix eq ̄mtion thaI c0m签from a eomtroined 14· 土木工程学报 200O越 optimizati ̄m problem.1he eur ̄nt improved model updating metodolbogy is applied on a 1/150 scald esuspension d萨rood- e1.1tIefmal uI 曲ed FEmodelforthe suspension b ̄asemodelis ableto produceIISIL!I ̄fl'eque ̄lciesin good agreementwith the measured one8. Key words:model updating;constrained optimizaticm;sensitivity analysis;suspension d . 范立础教授,博导。现任土木工程防灾国家重点实验室学术委员会常任副主任,并兼任中土桥梁与结构工程学会理事长 等国内外学术团体五个职务。我国著名的桥梁与结构工程专家,桥梁抗震学科的主要学术带头人。近十年来承担 “七五”公关项目,“八五”国家攀登B计划,国家自然科学基金会重大项目的多项课题,主持国家自然科学基金 会重点项目及其它2o多个重大工程科研项目,完成了绝大部分大跨桥梁的抗震研究。获省部级一等奖(参与)一 项,二等奖四项,三等奖一项,为我国大桥建设作出了重要贡献。1980年 来主编著作5部,发表论文120余篇, 培养了近40名研究生。通讯地址:2(XE02上海市四平路同济大学土木工程学院 寰万壤张启伟教授,博士,博导。现任同济大学桥梁工程系主任。主要从事桥梁结构的振动与抗震研究。近年来共发表论文30 余篇。1995年度获上海市优秀青年教师称号。 博士,副教授。主要从事桥梁可靠度设计标准与方法、桥梁结构空间应力分析、大跨度桥梁结构状态识别和健康 监测的理论研究以及桥梁评估与管理系统的开发和桥梁动静力试验与损伤检测等方面的工作。. (上接第8页) 所介绍的方法是有效的,算法是正确的,且具有较高 的求解效率与精度。 参考文献 ed and Prestressed Co/ ̄rete Beams.Proe.dthe 2nd Easec ( n ∞e .1989 2潘家英,程庆国.大跨度悬索桥有限位移分析.土木工程 学报,1994(1) 3张国政.铁路悬索桥非线性分析及其极限承载力的研究. [博士学位论文]北京:铁道部科学研究院,1994 1 Pan Jiaylns YangⅪn , 舭0f Strain S ̄fterdngfor Rein- GE0ⅧI1虹CAL AND MATERIAI,No Ⅱ胍^R^N^I Y圊瞎FOR DI ̄I'ERMINING III Tlm^ IE LIDAD CAPAaⅡ_Y OF LONG-SPAN BRIDGI ̄S Pan^咖 删 ( 档 (China Academy of Railway Sciences) Abstract An analytical pxoeedure considering both ge ̄netifcal and materil nonlaineariites for detemdning ultimate load o ̄iyt d l。 span bdd is describedinthis paper.upc edk afIgl岫fommlatin combionedwith rered Finiet ElementMethodis doeptdfeor i wiht geometrical aswell as rmterialnonlnieairties 0f brid8e stt'de ̄"ahiehmade it possibleto 鲫弛the accuracy of caleulati ̄andto save computationaltime and storage simultaneously.The procedure has been verifid eby expefi- mental results and employedto analyseanalternativeof railload suspensionHd suecessfuUy. Keywold:nonli ̄analysis,ultimateload capaciy,ltong-span grids ̄,raiUoad sspuausin d .o 潘家英研究员、博士生导师。中国土木工程学会桥粱与结构工程学会剐理事长。长期从事桥粱结构分析、研究。曾主持 多项部级 上重点课题。负责多座大跨度桥粱的施工监理工作。从事各类桥梁的可行性研究、方案设计及优化结 构分析、静动力测试、设计审查、施工监理、技术咨询等。通讯地址:1OO ̄l北京大柳树路2号铁道部科学研究 铁建所 张置政博士。现在美国休斯敦大学任教。 f疆i 司中科院院士、曾任研究员、中国土木工程学会及中国铁道学会剐理事长。曾主持研究开发一系列桥粱新结构及施 工新工艺。实现预应力混凝土粱工业化生产。在成昆线创议并建成申联梁、悬臂施工大跨径桥粱。曾主持我国第 座铁路斜拉桥——红水河桥的设计、施工及试验研究,并主持完成多项国家及部级重点课题。 (本文责任编辑刘芳春) 

因篇幅问题不能全部显示,请点此查看更多更全内容

Copyright © 2019- huatuo6.com 版权所有 湘ICP备2023023988号-11

违法及侵权请联系:TEL:199 1889 7713 E-MAIL:2724546146@qq.com

本站由北京市万商天勤律师事务所王兴未律师提供法律服务