CN106777549A - 一种面向荷载试验的桥梁多尺度有限元模拟方法 - Google Patents
一种面向荷载试验的桥梁多尺度有限元模拟方法 Download PDFInfo
- Publication number
- CN106777549A CN106777549A CN201611067415.6A CN201611067415A CN106777549A CN 106777549 A CN106777549 A CN 106777549A CN 201611067415 A CN201611067415 A CN 201611067415A CN 106777549 A CN106777549 A CN 106777549A
- Authority
- CN
- China
- Prior art keywords
- model
- bridge
- loading test
- controlling sections
- finite element
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种面向荷载试验的桥梁多尺度有限元模拟方法,其特征在于该方法包括以下步骤:第一步:结合设计竣工图纸建立全桥杆系单元模型,计算全桥结构的内力包络图,得出最不利受力的控制截面;第二步:将所述控制截面用板壳或实体单元模拟得到板壳或实体单元模型,同时删除所述控制截面的杆系单元模型;第三步:耦合为一体组成多尺度有限元模型,并验证其正确性;第四步:在验证后的所述多尺度有限元模型上进行荷载试验模拟。验证后的多尺度有限元模型由于模拟出了整体和局部受力情况,加快了计算速度和提高精度。
Description
技术领域
本发明涉及桥梁模拟,具体涉及一种面向荷载试验的桥梁多尺度有限元模拟方法。
背景技术
经过改革开放30多年的发展,我国桥梁建设取得了举世瞩目的成就。截至2015年底,我国公路桥梁总数已达77.92万座,4592.77万延米。我国已建成的梁桥、拱桥、斜拉桥和悬索桥这四类桥梁的跨径均居世界同类桥梁前列。然而,经过30多年的大规模基础设施建设,当前步入维修期的桥梁数量日益增多。据不完全统计,目前仅公路桥梁中危桥数量高达7.96万座,约占公路桥梁总数的10.5%,直接影响人民的正常安全出行。面向寿命期内桥梁养护的科学决策,需要桥梁技术人员对桥梁结构的技术状况和承载能力进行科学准确地评定。荷载试验是评定桥梁结构承载能力最直接、有效的鉴定方法,其核心思想是将桥梁结构关键参数的现场加载所获得的实测值与有限元模型计算的理论值进行比较,并对桥梁结构的承载能力进行评定。但是传统的杆系单元有限元模型有诸多弊端,如各测点实测值无法直接与理论值比较,不能够反映结构的局部刚度变化和质量分布等。
同时,桥梁荷载试验不仅承担着评定结构承载能力的作用而且为科研提供实测数据,如桥梁结构的剪力滞效应、偏载效应等。一般地,对于桥梁结构的局部分析是根据圣维南原理,采用子模型的方法,即首先进行整体分析,然后再对局部部位进行“二次分析”的方法。子模型方法需要对结构建立整体简化的有限元模型与局部精细化的有限元模型,并利用整体简化的有限元模型对结构进行整体分析,然后将整体分析得出的相应于局部精细有限元模型边界处的结果以边界约束和边界条件的方式施加到局部模型上,最后对局部精细有限元模型进行再分析。但是,该方法的弊端是局部精细化模型的边界条件析取与添加较为困难,因为从整体到局部,施加在局部模型上的边界条件的选取是很复杂的,边界条件选取不当直接会导致局部分析的结果相距甚远,同时子模型方法仍然是在单一尺度上模拟的,只不过整体简化模型和局部精细模型是分开的。
发明内容
本发明要解决的问题是针对现有技术中所存在的上述不足而提供一种面向荷载试验的桥梁多尺度有限元模拟方法可以很好地模拟桥梁的整体和局部受力情况,
为实现上述目的,本发明采用了如下的技术方案:一种面向荷载试验的桥梁多尺度有限元模拟方法,包括以下步骤:
第一步:结合设计竣工图纸建立全桥杆系单元模型(即单梁模型),通过施加设计荷载等级,计算全桥结构的内力包络图,得出最不利受力的截面,即为控制截面;
第二步:将所述控制截面用板壳或实体单元模拟得到板壳或实体单元模型,同时删除所述控制截面的杆系单元模型;
第三步:将已经删除了控制截面杆系单元模型的所述单梁模型和所述板壳或实体单元模型耦合为一体组成多尺度有限元模型;
第四步:在所述多尺度有限元模型上进行荷载试验模拟。
进一步,所述耦合方法有约束方程法。
进一步,通过对已经删除了控制截面杆系单元模型的所述单梁模型进行静力特性与动力特性的对比来验证多尺度模型的准确性。
一种基于所述面向荷载试验的桥梁多尺度有限元模拟方法的控制作用分析,基于桥梁模态分析的运动微分方程组进行求解并得到固有频率,选取结构前几阶自振频率和振型。
进一步,所述桥梁模态分析采用子空间迭代法。
相比于现有技术,本发明具有如下有益效果:验证后的多尺度有限元模型由于模拟出了整体和局部受力情况,在加快计算速度和提高精度的同时上进行荷载试验模拟,不仅可以在验证后的多尺度有限元模型上进行荷载试验方案设计以用于对桥梁结构的承载能力的评定,为桥梁交(竣)工验收及后期养护管理提供科学依据,还可以进行针对性的科研工作,诸如剪力滞效应、偏载效应(偏载系数、扭转变形、横向增大系数等)等。
附图说明
图1为本发明面向荷载试验的桥梁多尺度有限元模拟流程图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与作用更加清楚及易于了解,下面结合附图和具体实施方式对本发明作进一步阐述:
本发明提出了一种面向荷载试验的桥梁多尺度有限元模拟方法,包括以下步骤:
第一步:结合设计竣工图纸建立全桥杆系单元模型(即单梁模型),建模过程中实现建立、复制节点和单元或者输入荷载等,在划分单元时,需要考虑如下因素:a,研究所述单梁模型所需求的计算精度;b,结构的形态变化(包括梁的长度、宽度和厚度变化);c,桥面荷载和车辆活载作用的模拟要求;d,计算求解过程中为避免病态问题所要求的单元形态和单元之间的匹配;e,计算设备的能力;通过参考《公路桥涵设计通用规范》或者《城市桥梁设计规范》施加设计荷载等级,计算全桥结构的内力包络图采用(1)绘制某一截面内力的影响线;(2)确定最不利荷载位置;(3)计算最大和最小内力;(4)在拟作包络图的杆件上等距指定若干截面(称包络图控制截面),按(1)~(3)步,计算各截面在行列荷载作用下的最大和最小内力值;(5)以截面坐标为包络图控制点的横坐标,以上述最大和最小值为纵坐标,在图上标出控制点的位置;用平滑曲线分别连接最大点和最小点,即得所求包络图;然后选出得出最不利受力的截面即最大弯矩截面,也是控制截面;
第二步:采用有限元分析软件将所述控制截面用板壳或实体单元模拟,建立所述控制截面的板壳或实体单元模型;板壳单元的力学模型取为结构单元的中性面,即以各中性面来代表为不同厚度的板或壳单元的组合体,以此来模拟结构体;板壳单元通常是平面任意四边形或者任意三角形单元,或者空间上四个三角形单元构成锥形板所模拟的四边形薄壳单元;考虑到所述控制截面会存在局部细节的受力情况,如模拟梁的端板,螺栓连接的等情形,这时候需用所述实体单元进行模拟;删除所述控制截面的杆系单元模型,那么此时所述单梁模型将不具有所述控制截面的杆系单元。
第三步:将已经删除了控制截面杆系单元模型的所述单梁模型和所述板壳或实体单元模型耦合为一体组成多尺度有限元模型,即将所述板壳或实体单元嵌入到已删除控制截面的梁单元模型的全桥杆系有限元模型中,组合成多尺度有限元模型,此时结构的一部分离散为梁单元,一部分离散为壳单元,可采用约束方程法以及刚度叠加法,使用约束方程法需要逐一建立约束方程式,用约束方程法模拟梁结构,还可以采用伪梁法、刚性区域法(主从节点法)、MPC法(多点约束法)等;
第四步:在所述多尺度有限元模型上进行荷载试验模拟,试验依据参考《公路桥梁荷载试验规程》或者《城市桥梁检测与评定技术规范》实行,桥梁的自重作用在全桥单元上,桥面铺装和二期恒载等荷载采用均布荷载作用在全桥单元上,车辆活载按照静力等效原理作用在桥面板的所述板壳单元上,所述多尺度有限元模型由于模拟出了整体和局部受力情况,在加快计算速度和提高精度的同时上进行荷载试验方案设计以用于对桥梁结构的承载能力的评定,为桥梁交(竣)工验收及后期养护管理提供科学依据,还可以进行针对性的科研工作,诸如剪力滞效应、偏载效应(偏载系数、扭转变形、横向增大系数等)等。
作为具体实施例,所述耦合方法有约束方程法,因为平面单元没有转动自由度,所以通过采用约束方程法可以实现力矩的传递;例如梁单元与实体单元有一个节点位置重合,为使位移和力矩能够传递,则需要耦合两个节点的三个平移自由度,同时还需要用约束方程限制梁的三个转动自由度。
作为具体实施例,通过对已经删除了控制截面杆系单元模型的所述单梁模型进行静力特性(主要包括结构在恒载作用下的内力、支反力与位移响应)与动力特性(主要包括自振频率与模态振型)的对比来验证多尺度模型的准确性。
一种基于所述面向荷载试验的桥梁多尺度有限元模拟方法的控制作用分析,是基于桥梁模态分析的运动微分方程组进行求解并得到固有频率,在结构动力性能分析中选取结构前几阶自振频率和振型,因为一般只有前几阶自振频率和振型起控制作用,运动微分方程组可采用式中[M]、[C]、[K]分别为总质量矩阵、阻尼矩阵、刚度矩阵;{δ}为节点位移向量;***的特征行列式为|[K]-ω2[M]|=0。
作为具体实施例,所述桥梁模态分析采用子空间迭代法。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (5)
1.一种面向荷载试验的桥梁多尺度有限元模拟方法,其特征在于该方法包括以下步骤:
第一步:结合设计竣工图纸建立全桥杆系单元模型(即单梁模型),通过施加设计荷载等级,计算全桥结构的内力包络图,得出最不利受力的截面,即为控制截面;
第二步:将所述控制截面用板壳或实体单元模拟得到板壳或实体单元模型,同时删除所述控制截面的杆系单元模型;
第三步:将已经删除了控制截面杆系单元模型的所述单梁模型和所述板壳或实体单元模型耦合为一体组成多尺度有限元模型;
第四步:在所述多尺度有限元模型上进行荷载试验模拟。
2.根据权利要求1所述一种面向荷载试验的桥梁多尺度有限元模拟方法,其特征在于:所述耦合方法有约束方程法。
3.根据权利要求1所述一种面向荷载试验的桥梁多尺度有限元模拟方法,其特征在于:通过对已经删除了控制截面杆系单元模型的所述单梁模型进行静力特性与动力特性的对比来验证多尺度模型的准确性。
4.一种基于权利要求1所述面向荷载试验的桥梁多尺度有限元模拟方法的控制作用分析方法,其特征在于:基于桥梁模态分析的运动微分方程组进行求解并得到固有频率,选取结构前几阶自振频率和振型。
5.根据权利要求4所述的一种动力特性分析方法,其特征在于:所述桥梁模态分析采用子空间迭代法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611067415.6A CN106777549B (zh) | 2016-11-28 | 2016-11-28 | 一种面向荷载试验的桥梁多尺度有限元模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611067415.6A CN106777549B (zh) | 2016-11-28 | 2016-11-28 | 一种面向荷载试验的桥梁多尺度有限元模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106777549A true CN106777549A (zh) | 2017-05-31 |
CN106777549B CN106777549B (zh) | 2020-09-11 |
Family
ID=58902339
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611067415.6A Active CN106777549B (zh) | 2016-11-28 | 2016-11-28 | 一种面向荷载试验的桥梁多尺度有限元模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106777549B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108595733A (zh) * | 2018-01-30 | 2018-09-28 | 江西理工大学 | 一种间接测量控制荷载所产生变形的静力测试方法 |
CN109858071A (zh) * | 2018-12-06 | 2019-06-07 | 山西大学 | 一种考虑剪力滞后作用的薄壁箱梁结构动力特性分析方法 |
CN109933864A (zh) * | 2019-02-26 | 2019-06-25 | 广州交通投资集团有限公司 | 一种单箱多室混凝土箱形梁桥主梁腹板的设计方法 |
CN113378265A (zh) * | 2021-05-31 | 2021-09-10 | 中铁二院工程集团有限责任公司 | 基于多体动力学的桥梁动力性能的评估方法 |
CN113505405A (zh) * | 2021-04-15 | 2021-10-15 | 中南大学 | 等效荷载获取方法、基于等效荷载的拓扑优化方法及*** |
CN114692284A (zh) * | 2022-05-16 | 2022-07-01 | 西南交通大学 | 宽幅组合-混合梁斜拉桥结合段试验模型设计方法 |
CN115062522A (zh) * | 2022-08-18 | 2022-09-16 | 天河超级计算淮海分中心 | 基于装配式结构的强度确定方法、电子设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102819632A (zh) * | 2012-07-26 | 2012-12-12 | 苏州工业园区设计研究院股份有限公司 | 一种钢框架结构节点多尺度有限元模型建模方法 |
CN104866676A (zh) * | 2015-05-27 | 2015-08-26 | 东南大学 | 一种基于两阶段多尺度模型修正的结合梁斜拉桥传感器布设方法 |
CN106021644A (zh) * | 2016-05-06 | 2016-10-12 | 西北工业大学 | 确定混合维模型界面约束方程系数的方法 |
-
2016
- 2016-11-28 CN CN201611067415.6A patent/CN106777549B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102819632A (zh) * | 2012-07-26 | 2012-12-12 | 苏州工业园区设计研究院股份有限公司 | 一种钢框架结构节点多尺度有限元模型建模方法 |
CN104866676A (zh) * | 2015-05-27 | 2015-08-26 | 东南大学 | 一种基于两阶段多尺度模型修正的结合梁斜拉桥传感器布设方法 |
CN106021644A (zh) * | 2016-05-06 | 2016-10-12 | 西北工业大学 | 确定混合维模型界面约束方程系数的方法 |
Non-Patent Citations (1)
Title |
---|
杨红 等: "钢筋混凝土梁柱节点区非弹性变形的改进模型研究", 《建筑结构学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108595733A (zh) * | 2018-01-30 | 2018-09-28 | 江西理工大学 | 一种间接测量控制荷载所产生变形的静力测试方法 |
CN108595733B (zh) * | 2018-01-30 | 2022-03-11 | 江西理工大学 | 一种间接测量控制荷载所产生变形的静力测试方法 |
CN109858071A (zh) * | 2018-12-06 | 2019-06-07 | 山西大学 | 一种考虑剪力滞后作用的薄壁箱梁结构动力特性分析方法 |
CN109858071B (zh) * | 2018-12-06 | 2023-03-07 | 山西大学 | 一种考虑剪力滞后作用的薄壁箱梁结构动力特性分析方法 |
CN109933864A (zh) * | 2019-02-26 | 2019-06-25 | 广州交通投资集团有限公司 | 一种单箱多室混凝土箱形梁桥主梁腹板的设计方法 |
CN109933864B (zh) * | 2019-02-26 | 2023-05-23 | 广州交通投资集团有限公司 | 一种单箱多室混凝土箱形梁桥主梁腹板的设计方法 |
CN113505405A (zh) * | 2021-04-15 | 2021-10-15 | 中南大学 | 等效荷载获取方法、基于等效荷载的拓扑优化方法及*** |
CN113378265A (zh) * | 2021-05-31 | 2021-09-10 | 中铁二院工程集团有限责任公司 | 基于多体动力学的桥梁动力性能的评估方法 |
CN114692284A (zh) * | 2022-05-16 | 2022-07-01 | 西南交通大学 | 宽幅组合-混合梁斜拉桥结合段试验模型设计方法 |
CN115062522A (zh) * | 2022-08-18 | 2022-09-16 | 天河超级计算淮海分中心 | 基于装配式结构的强度确定方法、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106777549B (zh) | 2020-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106777549A (zh) | 一种面向荷载试验的桥梁多尺度有限元模拟方法 | |
Wang et al. | Fatigue design of steel bridges considering the effect of dynamic vehicle loading and overloaded trucks | |
CN106894328B (zh) | 一种π形结合梁剪力滞的处理方法 | |
CN104615840B (zh) | 一种数字仿真模型的修正方法及*** | |
Chen et al. | Multiscale numerical modeling of steel bridge deck pavements considering vehicle–pavement interaction | |
Xue et al. | Free vibration analysis of a flat stiffened plate with side crack through the Ritz method | |
Jiang et al. | The shear pin strength of friction pendulum bearings (FPB) in simply supported railway bridges | |
CN103678762A (zh) | 优化的复合材料机翼气动弹性风洞模型的缩比建模方法 | |
Beskou et al. | Dynamic analysis of an elastic plate on a cross-anisotropic elastic half-space under a rectangular moving load | |
Lu et al. | Computer simulation of the dynamic layered soil pile structure interaction system | |
Pippi et al. | Dynamic response to different models of adjacent coupled buildings | |
Zhang et al. | Dynamic response of multitower suspension bridge deck pavement under random vehicle load | |
Zhao et al. | Bridge weigh-in-motion algorithms based on the field calibrated simulation model | |
AsadiGorgi et al. | Effects of all-over part-through cracks on the aeroelastic characteristics of rectangular panels | |
Aydınoğlu et al. | Evaluation of analysis procedures for seismic assessment and retrofit design | |
Pu et al. | A two-level nesting smoothed extended meshfree method for static and dynamic fracture mechanics analysis of orthotropic materials | |
Zhang et al. | Simplified calculation method for the fundamental period of floating cable-stayed bridge | |
Aktas et al. | Computer-Aided Mold Design Algorithm for Precast Concrete Elements. | |
Zhou et al. | Damage Identification Method of Tied‐Arch Bridges Based on the Equivalent Thrust‐Influenced Line | |
Kikuchi et al. | Fatigue crack growth simulation using S-version FEM | |
Yuen | Structural identification of Tsing Ma bridge | |
Ng et al. | A new hybrid meshless-differential order reduction (hM-DOR) method with applications to shape control of smart structures via distributed sensors/actuators | |
Xu et al. | Fatigue behaviour study of the shear stud using crack box technology and m integral method | |
Attarchian et al. | Cyclic behavior modeling of rectangular RC bridge piers using OpenSees | |
Heitz et al. | Damping identification and quantification: experimental evidences and first numerical results |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |