CN103226626A - 一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法 - Google Patents
一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法 Download PDFInfo
- Publication number
- CN103226626A CN103226626A CN2013100629678A CN201310062967A CN103226626A CN 103226626 A CN103226626 A CN 103226626A CN 2013100629678 A CN2013100629678 A CN 2013100629678A CN 201310062967 A CN201310062967 A CN 201310062967A CN 103226626 A CN103226626 A CN 103226626A
- Authority
- CN
- China
- Prior art keywords
- bridge deck
- model
- steel bridge
- steel
- value
- 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
Images
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法,包括如下步骤:步骤10):对大跨缆索承重桥梁结构体系进行有限元模拟,建立整体结构尺度模型;步骤20):利用子模型法,建立整体结构尺度模型的钢箱梁局部精细模型,并计算其钢桥面板部位的等效应力幅精确解;步骤30):将钢箱梁局部精细模型简化为钢桥面板局部受力结构体系,并基于步骤20)的等效应力幅精确解对该结构体系的参数取值进行优化分析,确定钢桥面板疲劳荷载效应分析的简化模型。该简化模型可在保证计算结果准确性的同时,大幅提高正交异性钢桥面板疲劳荷载效应分析的计算效率。
Description
技术领域
本发明属于大跨缆索承重桥梁结构钢桥面板的疲劳分析与设计领域,具体来说,涉及一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法。
背景技术
正交异性钢桥面板由顶板、纵肋和横肋(扁平钢箱梁中称为横隔板)三部分组成,其中纵肋分布一般较密,而横肋分布则相对疏松,导致钢桥面板的纵横方向刚度不同,进而使得钢桥面板的受力特性表现为正交异性。正交异性钢桥面板的结构行为十分复杂,其除了具有桥面板和桥面系的作用外,还作为主梁的一部分发挥作用,例如当正交异性钢桥面板作为承受汽车荷载的局部受力构件时,顶板可视为四周边弹性支承在纵肋上,纵肋是弹性支承在横肋上,横肋则是弹性支承在主梁上。
因此,在对正交异性钢桥面板的疲劳荷载效应进行有限元分析时,通常需要建立其大跨缆索承重桥梁结构体系的整体结构尺度模型,然而此模型的单元网格划分较粗会导致计算结果的不够精确,单元网格划分较细又会导致庞大的计算量,因此在分析时往往利用子模型法,建立整体结构尺度模型的钢箱梁局部精细模型,采用钢箱梁局部精细模型计算正交异性钢桥面板的疲劳荷载效应精确解。然而子模型法也存在一定缺陷:①其边界约束条件会因疲劳荷载的不断变化而改变,使得计算过程变得复杂;②其单元网格的精细划分也会导致单元个数达到几十万个,使得计算量仍旧庞大。因此,在对正交异性钢桥面板进行疲劳荷载效应的有限元分析时,有必要找出一种计算效率高、计算结果精确的合理简化模型。
发明内容
发明目的:针对上述现有技术存在的问题和不足,本发明的目的是提供一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法,得到的简化模型可以高效准确地计算正交异性钢桥面板的疲劳荷载效应。
技术方案:为实现上述发明目的,本发明采用的技术方案为一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法,包括如下步骤:
步骤10):对大跨缆索承重桥梁结构体系进行有限元模拟,建立整体结构尺 度模型;
步骤20):利用子模型法,建立整体结构尺度模型的钢箱梁局部精细模型,并计算其钢桥面板部位的等效应力幅精确解Sa;
步骤30):将钢箱梁局部精细模型简化为钢桥面板局部受力结构体系,并基于步骤20)的等效应力幅精确解Sa对该结构体系的参数取值进行优化分析,确定钢桥面板疲劳荷载效应分析的简化模型。
进一步的,所述步骤10)中,根据大跨缆索承重桥梁结构体系的设计图纸和资料,对所述大跨缆索承重桥梁结构体系的桥塔、桥面系、边墩、斜拉索、主缆和吊杆进行整体结构尺度模型的三维有限单元模拟,其中桥塔采用三维梁单元模拟,桥面系中钢箱梁部分采用壳单元进行模拟,桥面系中铺装层部分采用实体单元进行模拟,边墩采用实体单元模拟,斜拉索、主缆和吊杆采用三维桁架单元模拟。
进一步的,所述整体结构尺度模型的边界连接条件处理为:①主梁与索塔的连接:主梁与索塔之间设有横向抗风支座和竖向抗风支座;②索塔与边墩底部固结;以及③主梁在各边墩处,横桥向、竖向、扭转及绕竖向的转动自由度与墩顶节点建立主从关系。
进一步的,所述步骤20)中,包括如下步骤:
(a)从整体结构尺度模型中切割出一段钢箱梁,并进行精细网格划分,网格划分水平取为5,得到钢箱梁局部精细模型;
(b)对步骤10)得到的整体结构尺度模型进行较粗网格划分,网格划分水平取为1,施加车辆荷载并计算求解,得出钢箱梁局部精细模型的边界部位的位移响应;
(c)将第(b)步所得的位移响应作为边界条件,采用线性插值法施加到钢箱梁局部精细模型的切割边界;
(d)保持车辆荷载加载不变,对钢箱梁局部精细模型进行计算求解,得出其钢桥面板部位的等效应力幅精确解Sa。
进一步的,所述步骤30)中,包括如下步骤:
a)将钢箱梁局部精细模型简化为仅由顶板和纵肋组成的钢桥面板局部受 力结构体系;
b)将该结构体系的顺桥向长度和横桥向长度分别采用横隔板间距个数W和纵肋个数L表示;将该结构体系的边界约束条件分为3类:①边界I:周边横隔板的约束条件;②边界II:周边纵肋的约束条件;③边界III:中间横隔板的约束条件;每类边界约束条件包括沿节点坐标系X、Y、Z轴的平动约束刚度Cx、Cy、Cz和绕节点坐标系X、Y、Z轴的转动约束刚度Cyz、Cxz、Cxy;因此,钢桥面板局部受力结构体系的待确定参数包括几何尺寸参数W、L以及三类边界约束刚度参数Cx、Cy、Cz、Cyz、Cxz、Cxy,共计20个计算参数;
c)保持车辆荷载加载,以等效应力幅S为分析指标对20个计算参数进行敏感度分析,认为等效应力幅S是以20个计算参数为自变量的函数,其中20个计算参数采用集合{V}表示,若集合{V}中第i个计算参数用Vi表示,Vi的第j、j+1次取值用Vi,j、Vi,j+1表示,则利用式(1)可计算Vi在区间(Vi,j,Vi,j+1)内的敏感系数F(Vi,j,Vi,j+1):
式中,分别为集合{V}中Vi取值为Vi,j、Vi,j+1时钢桥面板局部受力结构体系的等效应力幅计算值,且i=1,2,…,20,j为正整数;若根据式(1)计算得到的F(Vi,j,Vi,j+1)<0.2,且与Sa的相对误差均小于2.0%,则将Vi视为在该区间(Vi,j,Vi,j+1)内的非敏感参数,从区间(Vi,j,Vi,j+1)中直接选取任一值作为该参数的最优值;若根据式(1)计算得到的F(Vi,j,Vi,j+1)≥0.2,则将Vi视为在该区间(Vi,j,Vi,j+1)内的敏感参数,需要对其取值作进一步优化分析;
其中[Fk]+=[Fk]T([Fk][Fk]T)-1 (2b)
式中,和分别为向量的第k次和第k+1次迭代值,Sk为基 于向量得到的钢桥面板局部受力结构体系的等效应力幅计算值,[Fk]为在区间内的敏感系数向量,是式(1)的向量表示形式,[Fk]+是[Fk]的可逆矩阵,k=0、1、…、m,m的最小取值应满足Sk与Sa之间的相对误差小于2.0%;
e)利用集合{V}中非敏感参数和敏感参数的最优值,最终确定钢桥面板疲劳荷载效应分析的简化模型。施加其它车辆荷载加载方式,求解步骤20)钢箱梁局部精细模型的精确解和简化模型的计算解,对比验证简化模型的合理性。
有益效果:与现有技术相比,本发明具有以下有益效果:
该简化模型在用于正交异性钢桥面板的疲劳荷载效应分析时,可在保证计算结果准确性的同时,大幅提高计算效率:①简化模型是以等效应力幅精确解Sa为分析指标对20个计算参数进行优化分析得到的,因此可保证简化模型计算结果的准确性;②与整体结构尺度模型和钢箱梁局部精细模型相比,简化模型仅由顶板和纵肋组成,大大减少了单元网格划分个数;此外,简化模型的边界约束条件是固定不变的,不会因疲劳荷载的变化而改变,降低了疲劳荷载效应计算的复杂性。因此,该简化模型在用于正交异性钢桥面板的疲劳荷载效应分析时,建模简单、计算高效、结果准确,具有很好的实用价值,可得到广泛推广和应用。
附图说明
图1为本发明实施例钢桥面板疲劳荷载效应分析模型简化流程图;
图2为本发明实施例润扬长江大桥斜拉桥整体结构尺度模型;
图3为本发明实施例润扬长江大桥斜拉桥钢箱梁局部精细模型;
图4(a)为本发明实施例标准疲劳车荷载模型,图4(b)为图4(a)中的标准疲劳车各轴轮胎着地面示意图;
图5(a)、图5(b)和图5(c)分别为本发明实施例的三种车辆荷载加载工况及疲劳应力监测位置示意图;
图6为本发明实施例钢桥面板局部受力结构体系图;
图7(a)为本发明实施例钢桥面板局部受力结构体系的横桥向横截面边界条件示意图,图7(b)为本发明实施例钢桥面板局部受力结构体系的顺桥向横截面边界条件示意图;
图8(a)为本发明实施例顶板-纵肋焊接几何尺寸参数W、L的敏感度分析图,图8(b)为本发明实施例纵肋对接焊接几何尺寸参数W、L的敏感度分析图;
图9(a)为本发明实施例顶板-纵肋焊接边界条件I的敏感度分析图,图9(b)为本发明实施例纵肋对接焊接边界条件I的敏感度分析图;
图10(a)为本发明实施例顶板-纵肋焊接边界条件II的敏感度分析图,图10(b)为本发明实施例纵肋对接焊接边界条件II的敏感度分析图;
图11(a)为本发明实施例顶板-纵肋焊接边界条件III的敏感度分析图,图11(b)为本发明实施例纵肋对接焊接边界条件III的敏感度分析图;
图12(a)为本发明实施例顶板-纵肋焊接边界条件III中3个平动约束刚度的敏感度分析图,图12(b)为本发明实施例纵肋对接焊接边界条件III中3个平动约束刚度的敏感度分析图。
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
本发明的一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法,该简化方法包括如下步骤,如图1所示:
步骤10):对大跨缆索承重桥梁结构体系进行有限元模拟,建立整体结构尺度模型:
基于ANSYS有限元分析软件,根据大跨缆索承重桥梁结构体系的设计图纸和资料,对其桥塔、桥面系、斜拉索(或主缆及吊杆)、边墩和支承连接装置等进行整体结构尺度模型的三维有限单元模拟,其中桥塔采用三维梁单元模拟,钢箱梁的顶板、纵肋、横隔板和纵隔板等采用壳单元进行模拟,桥面铺装层采用实体单元进行模拟,斜拉索或主缆及吊杆采用三维桁架单元模拟,模拟时考虑缆索垂度效应引起的弹性模量折减以及恒载作用下的几何刚度;
此外,整体结构尺度模型的边界连接条件处理为:①主梁与索塔的连接:主梁与索塔之间设有横向抗风支座和竖向抗风支座;②索塔与边墩底部固结;③主梁在各边墩处,横桥向、竖向、扭转及绕竖向的转动自由度与墩顶节点建立主从 关系;
步骤20):利用子模型法,建立整体结构尺度模型的钢箱梁局部精细模型,并计算其钢桥面板部位的等效应力幅精确解Sa:
(a)从整体结构尺度模型中切割出一段钢箱梁,并进行精细网格划分,得到钢箱梁局部精细模型;
(b)对步骤10)得到的整体结构尺度模型进行较粗网格划分,施加车辆荷载并计算求解,得出其特定部位(钢箱梁局部精细模型的边界部位)的位移响应;
(c)将第(b)步所得的位移响应作为边界条件,采用线性插值法施加到钢箱梁局部精细模型的切割边界;
(d)保持车辆荷载加载不变,对钢箱梁局部精细模型进行计算求解,得出其钢桥面板部位的等效应力幅精确解Sa;
步骤30):将钢箱梁局部精细模型简化为钢桥面板局部受力结构体系,并基于步骤20)的等效应力幅精确解Sa对此结构体系的参数取值进行优化分析,确定钢桥面板疲劳荷载效应分析的简化模型:
a)将钢箱梁局部精细模型简化为仅由顶板和纵肋组成的钢桥面板局部受力结构体系;
b)将该结构体系的顺桥向长度和横桥向长度分别采用横隔板间距个数W和纵肋个数L表示。此外,将该结构体系的边界约束条件分为3类:①边界I:周边横隔板的约束条件;②边界II:周边纵肋的约束条件;③边界III:中间横隔板的约束条件。每类边界约束条件包括沿节点坐标系X、Y、Z轴的平动约束刚度Cx、Cy、Cz和绕节点坐标系X、Y、Z轴的转动约束刚度Cyz、Cxz、Cxy。因此,钢桥面板局部受力结构体系的待确定参数包括几何尺寸参数W、L以及三类边界约束刚度参数Cx、Cy、Cz、Cyz、Cxz、Cxy,共计20个计算参数。
c)保持车辆荷载加载与步骤20)相同,以等效应力幅S为分析指标对20个计算参数进行敏感度分析,认为等效应力幅S是以20个计算参数为自变量的函数,其中20个计算参数采用集合{V}表示,若集合{V}中第i个计算参数用Vi表示,Vi的第j、j+1次取值用Vi,j、Vi,j+1表示,则利用式(1)可计算Vi在区间(Vi,j, Vi,j+1)内的敏感系数F(Vi,j,Vi,j+1):
式中,分别为集合{V}中Vi取值为Vi,j、Vi,j+1时钢桥面板局部受力结构体系的等效应力幅计算值,且i=1,2,…,20,j=1,2,…。若根据式(1)计算得到的F(Vi,j,Vi,j+1)较小且接近Sa,则将Vi视为在该区间(Vi,j,Vi,j+1)内的非敏感参数,且可从区间(Vi,j,Vi,j+1)中直接选取任一值作为该参数的最优值;若根据式(1)计算得到的F(Vi,j,Vi,j+1)较大,则将Vi视为在该区间(Vi,j,Vi,j+1)内的敏感参数,需要对其取值作进一步优化分析;
其中[Fk]+=[Fk]T([Fk][Fk]T)-1 (2b)
式中,和分别为向量的第k、k+1次迭代值,Sk为基于向量 得到的钢桥面板局部受力结构体系的等效应力幅计算值,[Fk]为在区间 内的敏感系数向量,是式(1)的向量表示形式,k=0、1、…、m,m的最小取值应满足Sk与Sa之间的相对误差小于2.0%。
e)利用集合{V}中非敏感参数和敏感参数的最优值,可最终确定钢桥面板疲劳荷载效应分析的简化模型。施加其它车辆荷载加载方式,求解步骤20)钢箱梁局部精细模型的精确解和简化模型的计算解,对比验证简化模型的合理性。
本发明的一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法,首先对大跨缆索承重桥梁结构体系进行有限元模拟,建立整体结构尺度模型,然后利用子模型法,建立整体结构尺度模型的钢箱梁局部精细模型,并计算其钢桥面板部位的等效应力幅精确解,最后将钢箱梁局部精细模型简化为钢桥面板局部受力结构体系,并基于步骤20)的等效应力幅精确解Sa对此结构体系的参数取值进行 优化分析,从而确定钢桥面板疲劳荷载效应分析的简化模型。
实施例
下面以润扬大桥北汊斜拉桥正交异性钢桥面板的疲劳荷载效应分析为例,说明本发明的具体实施过程:
(1)基于步骤10)建立的润扬大桥北汊斜拉桥整体结构尺度模型,如图2所示,基于步骤20)中的第(a)步建立的钢箱梁局部精细模型,如图3所示,选取图4(a)和图4(b)中的标准疲劳车作为加载荷载,其横桥向加载位置如图5(a)所示,按照步骤20)中的第(b)~(d)步计算钢桥面板顶板-纵肋焊接细节(图5(a)中A点所示)和纵肋对接焊接细节(图5(a)中B点所示)的等效应力幅精确解Sa,计算结果分别为29.58MPa和38.15MPa。
(2)基于步骤20)中的第(a)步建立的钢桥面板局部受力结构体系如图6所示,其边界条件示意图如图7(a)和图7(b)所示,按照步骤30)中的第b)~c)步对钢桥面板局部受力结构体系的20个计算参数进行敏感度分析:
1)首先对横隔板间距个数W和纵肋个数L进行敏感度分析。横隔板间距个数W分别取为1、3、5、7;纵肋个数L分别取为2、4、6、8。计算与步骤20)相同的加载荷载下顶板-纵肋焊接细节和纵肋对接焊接细节的等效应力幅Seq。此时,简化分析模型中的边界条件I、II和III均为铰接。为方便等效应力幅计算值和准确值的比较,在给出等效应力幅随参数变化的趋势线的同时,在图中用水准线标识出等效应力幅准确值Seq的位置。敏感度分析结果如图8(a)和图8(b)所示,利用趋势线的变化率反映参数不同取值的敏感系数大小。
从图8(a)和图8(b)中可以看出:几何尺寸参数W的敏感系数变化以W=3为分界点,取值小于3视为敏感参数,大于3视为非敏感参数;几何尺寸参数L的趋势线以L=4为分界点,分为斜坡段和平直段,在斜坡段取值的L视为敏感参数,在平直段取值的L视为非敏感参数。当横隔板间距个数W取值为5、纵肋个数L取值为6时,几何尺寸参数对等效应力幅计算误差的影响已经很小。根据上述分析,钢桥面板简化分析模型的横隔板间距个数W取为5,纵肋个数L取为6。
2)下面对边界约束刚度参数Cx、Cy、Cz、Cyz、Cxz、Cxy进行敏感度分析。为了简化分析,首先将每一边界约束简化为自由、铰结和固结三种约束形式,除 三类边界条件全部自由外共存在26种可能的组合,计算每种组合在与步骤20)相同的加载荷载下两类焊接细节的等效应力幅。图9~图11给出了边界条件I、II、III分别按自由、铰结和固结三种约束时的计算结果。每一边界约束可得到多条趋势线。为便于分析,利用均值法绘出反映趋势线整体走势的拟合趋势线,并进行敏感度分析。
由图9~图11中可以看出:①三类边界条件下顶板-纵肋焊接细节和纵肋对接焊接细节的等效应力幅拟合趋势线均以铰接为分界点,分为斜坡段和平直段,斜坡段内的边界约束取值视为敏感参数,平直段内的边界约束取值视为非敏感参数。当边界条件为铰接和固结时,等效应力幅基本相同,因此,三类边界条件的约束刚度参数Cyz、Cxz、Cxy均为非敏感参数,并且对等效应力幅计算误差的影响很小;②对图11(a)和图11(b)的分析结果表明,各趋势线的走势趋于一致且均逼近准确值,并且拟合趋势线在斜坡段具有较大的斜率。这说明边界条件I和II中约束刚度的取值对第3类边界的趋势线干扰性很小,因此,边界条件I和II的约束刚度参数Cx、Cy、Cz为非敏感参数,等效应力幅的计算精度主要受第III类边界约束在斜坡段取值的影响,需要对第3类边界约束在斜坡段的取值作进一步分析。
3)第III类边界条件在斜坡段的约束刚度可分以解为x、y、z三个方向的平动弹性刚度约束,当约束刚度均为0时为自由边界,无穷大时为铰结边界。分别计算加载工况1(如图5(a)所示)下不同平动弹性刚度取值时两类焊接细节的等效应力幅,并利用其趋势线作敏感度分析,如图12(a)和图12(b)所示。从中可以看出,x、y向平动约束刚度的取值对等效应力幅的影响很小,等效应力幅主要受z向平动约束刚度的影响。因此,边界条件III的约束刚度参数中Cx、Cy为非敏感参数,Cz为敏感参数,需要通过优化分析确定。
(3)基于步骤30)中的第d)步对第III类边界约束条件的计算参数Cz进行优化分析。首先建立钢桥面板疲劳荷载效应分析的初始模型,其横隔板间距个数W和纵肋个数L分别取为5和6,第I类和第II类边界约束条件均处理为铰接,第III类边界约束条件中计算参数Cz为敏感参数,初始迭代值取为106N/m,其余边界约束均为铰接。迭代优化分析时取与步骤20)相同的加载荷载作为目标工况,采用式(2a)、(2b)和(2c)对计算参数Cz进行迭代,当参数迭代值使 等效应力幅的计算值与步骤20)的等效应力幅精确值Sa之间的相对误差小于2.0%时即停止迭代,此时的参数迭代值作为参数优化值,计算参数Cz的迭代过程及结果如表1所示,由表1可知约束刚度Cz的第5次迭代值为1.388×107N/m时,两类焊接细节的等效应力幅相对误差均小于2.0%。基于步骤30)中的第e)步,以迭代结果为参数值建立钢桥面板疲劳分析模型,并进一步计算图5(b)和图5(c)所示加载工况下两类焊接细节的等效应力幅,表2给出了等效应力幅的对比结果。从表2中可以看出,等效应力幅计算值与准确值的最大相对误差仅为3.18%。
综上分析,通过对模型计算参数的敏感度分析以及敏感参数的优化结果,润扬大桥斜拉桥钢桥面板疲劳荷载效应的最优简化分析模型为:横桥向尺寸取6个U形纵肋,顺桥向尺寸取5个横隔板间距,第I类边界铰接,第II类边界自由,第III类边界沿Z向弹性约束刚度为1.388×107N/m,其余为铰接。
表1
表2
Claims (5)
1.一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法,其特征在于,包括如下步骤:
步骤10):对大跨缆索承重桥梁结构体系进行有限元模拟,建立整体结构尺度模型;
步骤20):利用子模型法,建立整体结构尺度模型的钢箱梁局部精细模型,并计算其钢桥面板部位的等效应力幅精确解Sa;
步骤30):将钢箱梁局部精细模型简化为钢桥面板局部受力结构体系,并基于步骤20)的等效应力幅精确解Sa对该结构体系的参数取值进行优化分析,确定钢桥面板疲劳荷载效应分析的简化模型。
2.根据权利要求1所述一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法,其特征在于:所述步骤10)中,对所述大跨缆索承重桥梁结构体系的桥塔、桥面系、边墩、斜拉索、主缆和吊杆进行整体结构尺度模型的三维有限单元模拟,其中桥塔采用三维梁单元模拟,桥面系中钢箱梁部分采用壳单元进行模拟,桥面系中铺装层部分采用实体单元进行模拟,边墩采用实体单元模拟,斜拉索、主缆和吊杆采用三维桁架单元模拟。
3.根据权利要求1所述一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法,其特征在于:所述整体结构尺度模型的边界连接条件处理为:①主梁与索塔的连接:主梁与索塔之间设有横向抗风支座和竖向抗风支座;②索塔与边墩底部固结;以及③主梁在各边墩处,横桥向、竖向、扭转及绕竖向的转动自由度与墩顶节点建立主从关系。
4.根据权利要求1所述一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法,其特征在于:所述步骤20)中,包括如下步骤:
(a)从整体结构尺度模型中切割出一段钢箱梁,并进行精细网格划分,网格划分水平取为5,得到钢箱梁局部精细模型;
(b)对步骤10)得到的整体结构尺度模型进行较粗网格划分,网格划分水平取为1,施加车辆荷载并计算求解,得出钢箱梁局部精细模型的边界部位的位移响应;
(c)将第(b)步所得的位移响应作为边界条件,采用线性插值法施加到钢箱梁局部精细模型的切割边界;
(d)保持车辆荷载加载不变,对钢箱梁局部精细模型进行计算求解,得出其钢桥面板部位的等效应力幅精确解Sa。
5.根据权利要求1所述一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法,其特征在于:所述步骤30)中,包括如下步骤:
a)将钢箱梁局部精细模型简化为仅由顶板和纵肋组成的钢桥面板局部受力结构体系;
b)将该结构体系的顺桥向长度和横桥向长度分别采用横隔板间距个数W和纵肋个数L表示;将该结构体系的边界约束条件分为3类:①边界I:周边横隔板的约束条件;②边界II:周边纵肋的约束条件;③边界III:中间横隔板的约束条件;每类边界约束条件包括沿节点坐标系X、Y、Z轴的平动约束刚度Cx、Cy、Cz和绕节点坐标系X、Y、Z轴的转动约束刚度Cyz、Cxz、Cxy;因此,钢桥面板局部受力结构体系的待确定参数包括几何尺寸参数W、L以及三类边界约束刚度参数Cx、Cy、Cz、Cyz、Cxz、Cxy,共计20个计算参数;
c)保持车辆荷载加载,以等效应力幅S为分析指标对20个计算参数进行敏感度分析,认为等效应力幅S是以20个计算参数为自变量的函数,其中20个计算参数采用集合{V}表示,若集合{V}中第i个计算参数用Vi表示,Vi的第j、j+1次取值用Vi,j、Vi,j+1表示,则利用式(1)可计算Vi在区间(Vi,j,Vi,j+1)内的敏感系数F(Vi,j,Vi,j+1):
式中,分别为集合{V}中Vi取值为Vi,j、Vi,j+1时钢桥面板局部受力结构体系的等效应力幅计算值,且i=1,2,…,20,j为正整数;若根据式(1)计算得到的F(Vi,j,Vi,j+1)<0.2,且与Sa的相对误差均小于2.0%,则将Vi视为在该区间(Vi,j,Vi,j+1)内的非敏感参数,从区间(Vi,j,Vi,j+1)中直接选取任一值作为该参数的最优值;若根据式(1)计算得到的F(Vi,j,Vi,j+1)≥0.2,则将Vi视为在该区间(Vi,j,Vi,j+1)内的敏感参数;
其中[Fk]+=[Fk]T([Fk][Fk]T)-1 (2b)
式中,和分别为向量的第k次和第k+1次迭代值,Sk为基于向量得到的钢桥面板局部受力结构体系的等效应力幅计算值,[Fk]为在区间内的敏感系数向量,是式(1)的向量表示形式,[Fk]+是[Fk]的可逆矩阵,k=0、1、…、m,m的最小取值应满足Sk与Sa之间的相对误差小于2.0%;
e)利用集合{V}中非敏感参数和敏感参数的最优值,最终确定钢桥面板疲劳荷载效应分析的简化模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310062967.8A CN103226626B (zh) | 2013-02-27 | 2013-02-27 | 一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310062967.8A CN103226626B (zh) | 2013-02-27 | 2013-02-27 | 一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103226626A true CN103226626A (zh) | 2013-07-31 |
CN103226626B CN103226626B (zh) | 2016-06-29 |
Family
ID=48837071
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310062967.8A Expired - Fee Related CN103226626B (zh) | 2013-02-27 | 2013-02-27 | 一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103226626B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104991991A (zh) * | 2015-05-22 | 2015-10-21 | 南通大学 | 一种切条机切割过程数值模拟方法 |
CN105160142A (zh) * | 2014-06-03 | 2015-12-16 | 中铁大桥科学研究院有限公司 | 用于延长正交异性桥面板钢箱梁桥梁寿命的管理养护方法 |
CN105716965A (zh) * | 2016-03-29 | 2016-06-29 | 同济大学 | 正交异性钢桥面混凝土铺装疲劳性能评价装置 |
CN105808876A (zh) * | 2016-03-21 | 2016-07-27 | 潘燕 | 钢桁梁斜拉桥主塔墩顶钢梁架设方法 |
CN106326600A (zh) * | 2016-09-18 | 2017-01-11 | 中国铁建大桥工程局集团有限公司 | 基于张拉力或强制位移值的主缆敏感性工程化分析方法 |
CN106991255A (zh) * | 2017-05-04 | 2017-07-28 | 中铁二院工程集团有限责任公司 | 一种通过桥面系控制铁路大跨度斜拉桥刚度的方法 |
CN106991254A (zh) * | 2017-04-28 | 2017-07-28 | 武汉理工大学 | 一种公铁两用大桥钢桥面铺装的多尺度设计方法 |
CN107016194A (zh) * | 2017-04-07 | 2017-08-04 | 东南大学 | 基于多尺度有限元建模的桥梁缆索体系劣化性能分析方法 |
CN107895094A (zh) * | 2017-12-14 | 2018-04-10 | 河南省交通规划设计研究院股份有限公司 | 用于正交异性钢‑混组合桥面板界面抗剪的设计方法 |
CN108959725A (zh) * | 2018-06-11 | 2018-12-07 | 东南大学 | 钢箱梁焊接残余应力和结构应力耦合计算模型的建立方法 |
CN111274646A (zh) * | 2020-03-12 | 2020-06-12 | 北京建筑大学 | 钢桥疲劳敏感区的合成疲劳应力获取方法及装置 |
CN111797447A (zh) * | 2019-12-11 | 2020-10-20 | 中国电建集团成都勘测设计研究院有限公司 | 一种复杂地基混凝土坝孔道配筋方法 |
CN113295688A (zh) * | 2021-04-16 | 2021-08-24 | 开封大学 | 一种焊接部位疲劳强度的检测系统及检测方法 |
CN113722794A (zh) * | 2021-08-16 | 2021-11-30 | 中交第二航务工程局有限公司 | 一种正交异性钢桥面系板桁组合结构有限元模拟方法 |
CN114406552A (zh) * | 2022-03-04 | 2022-04-29 | 中国石油大学(华东) | 一种用于转鼓式压滤机转鼓的抗疲劳焊接制造方法 |
CN114840897A (zh) * | 2022-05-16 | 2022-08-02 | 金陵科技学院 | 钢箱梁焊接节点残余应力与车致结构应力耦合计算方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013008223A (ja) * | 2011-06-24 | 2013-01-10 | Toyota Motor Corp | スポット溶接部応力算出方法 |
-
2013
- 2013-02-27 CN CN201310062967.8A patent/CN103226626B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013008223A (ja) * | 2011-06-24 | 2013-01-10 | Toyota Motor Corp | スポット溶接部応力算出方法 |
Non-Patent Citations (4)
Title |
---|
丁幼亮等: "润扬大桥斜拉桥结构安全评估的有限元建模与修正", 《东南大学学报(自然科学版)》 * |
李爱群等: "子模型法在超大跨悬索桥钢箱梁应力分析中的应用", 《工程力学》 * |
王华毅等: "拱塔斜拉桥正交异性钢桥面板疲劳荷载应力幅分析", 《中国水运》 * |
邓强民等: "大跨钢桥桥面铺装有限元分析合理简化模型", 《交通运输工程学报》 * |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105160142A (zh) * | 2014-06-03 | 2015-12-16 | 中铁大桥科学研究院有限公司 | 用于延长正交异性桥面板钢箱梁桥梁寿命的管理养护方法 |
CN105160142B (zh) * | 2014-06-03 | 2018-06-26 | 中铁大桥科学研究院有限公司 | 用于延长正交异性桥面板钢箱梁桥梁寿命的管理养护方法 |
CN104991991B (zh) * | 2015-05-22 | 2018-05-04 | 南通大学 | 一种切条机切割过程数值模拟方法 |
CN104991991A (zh) * | 2015-05-22 | 2015-10-21 | 南通大学 | 一种切条机切割过程数值模拟方法 |
CN105808876A (zh) * | 2016-03-21 | 2016-07-27 | 潘燕 | 钢桁梁斜拉桥主塔墩顶钢梁架设方法 |
CN105716965B (zh) * | 2016-03-29 | 2018-06-26 | 同济大学 | 正交异性钢桥面混凝土铺装疲劳性能评价装置 |
CN105716965A (zh) * | 2016-03-29 | 2016-06-29 | 同济大学 | 正交异性钢桥面混凝土铺装疲劳性能评价装置 |
CN106326600A (zh) * | 2016-09-18 | 2017-01-11 | 中国铁建大桥工程局集团有限公司 | 基于张拉力或强制位移值的主缆敏感性工程化分析方法 |
CN107016194A (zh) * | 2017-04-07 | 2017-08-04 | 东南大学 | 基于多尺度有限元建模的桥梁缆索体系劣化性能分析方法 |
CN106991254A (zh) * | 2017-04-28 | 2017-07-28 | 武汉理工大学 | 一种公铁两用大桥钢桥面铺装的多尺度设计方法 |
CN106991254B (zh) * | 2017-04-28 | 2020-03-10 | 武汉理工大学 | 一种公铁两用大桥钢桥面铺装的多尺度设计方法 |
CN106991255A (zh) * | 2017-05-04 | 2017-07-28 | 中铁二院工程集团有限责任公司 | 一种通过桥面系控制铁路大跨度斜拉桥刚度的方法 |
CN107895094A (zh) * | 2017-12-14 | 2018-04-10 | 河南省交通规划设计研究院股份有限公司 | 用于正交异性钢‑混组合桥面板界面抗剪的设计方法 |
CN107895094B (zh) * | 2017-12-14 | 2020-03-10 | 河南省交通规划设计研究院股份有限公司 | 用于正交异性钢-混组合桥面板界面抗剪的设计方法 |
CN108959725A (zh) * | 2018-06-11 | 2018-12-07 | 东南大学 | 钢箱梁焊接残余应力和结构应力耦合计算模型的建立方法 |
CN111797447A (zh) * | 2019-12-11 | 2020-10-20 | 中国电建集团成都勘测设计研究院有限公司 | 一种复杂地基混凝土坝孔道配筋方法 |
CN111797447B (zh) * | 2019-12-11 | 2022-05-31 | 中国电建集团成都勘测设计研究院有限公司 | 一种复杂地基混凝土坝孔道配筋方法 |
CN111274646A (zh) * | 2020-03-12 | 2020-06-12 | 北京建筑大学 | 钢桥疲劳敏感区的合成疲劳应力获取方法及装置 |
CN113295688A (zh) * | 2021-04-16 | 2021-08-24 | 开封大学 | 一种焊接部位疲劳强度的检测系统及检测方法 |
CN113722794A (zh) * | 2021-08-16 | 2021-11-30 | 中交第二航务工程局有限公司 | 一种正交异性钢桥面系板桁组合结构有限元模拟方法 |
CN113722794B (zh) * | 2021-08-16 | 2023-09-19 | 中交第二航务工程局有限公司 | 一种正交异性钢桥面系板桁组合结构有限元模拟方法 |
CN114406552A (zh) * | 2022-03-04 | 2022-04-29 | 中国石油大学(华东) | 一种用于转鼓式压滤机转鼓的抗疲劳焊接制造方法 |
CN114406552B (zh) * | 2022-03-04 | 2023-10-10 | 中国石油大学(华东) | 一种用于转鼓式压滤机转鼓的抗疲劳焊接制造方法 |
CN114840897A (zh) * | 2022-05-16 | 2022-08-02 | 金陵科技学院 | 钢箱梁焊接节点残余应力与车致结构应力耦合计算方法 |
CN114840897B (zh) * | 2022-05-16 | 2024-02-02 | 金陵科技学院 | 钢箱梁焊接节点残余应力与车致结构应力耦合计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103226626B (zh) | 2016-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103226626A (zh) | 一种正交异性钢桥面板的疲劳荷载效应分析模型简化方法 | |
CN102767133B (zh) | 一种钢箱梁桥疲劳应力监测的传感器布置方法 | |
CN110377943A (zh) | 一种基于移动荷载试验的桥梁承载能力评估方法 | |
Bell et al. | Objective load rating of a steel-girder bridge using structural modeling and health monitoring | |
CN102254068A (zh) | 大跨度桥梁抖振响应的多尺度分析方法 | |
CN105092261A (zh) | 道路载荷测试方法及系统 | |
CN104133959A (zh) | 桥梁有限元模型的修正方法 | |
CN104217094A (zh) | 用于计算结构的疲劳以及疲劳破坏的方法 | |
CN102384856A (zh) | 基于概率有限元法的钢桥疲劳可靠度评估方法 | |
CN102243671A (zh) | 大跨钢桥扁平钢箱梁的温度梯度效应分析方法 | |
CN103279611A (zh) | 一种应变传感器的优化布置方法 | |
Lin et al. | Model updating using bridge influence lines based on an adaptive metamodel global optimization method | |
CN108090283A (zh) | 一种列车动载下隧道—车辆耦合振动的有限元分析方法 | |
CN107817038A (zh) | 一种轨道车辆轮重通过式智能化检测装置及方法 | |
Sadeghi et al. | Damage detection of composite beams via variational mode decomposition of shear-slip data | |
KR101094221B1 (ko) | 수문 및 수리분석 통합시스템 | |
CN106354929A (zh) | 基于刚度变化原理的承载结构载荷传递路径可视化方法 | |
Calçada et al. | Analysis of traffic-induced vibrations in a cable-stayed bridge. Part I: Experimental assessment | |
Akintunde et al. | Singular value decomposition and unsupervised machine learning for virtual strain sensing: Application to an operational railway bridge | |
Tarboton et al. | Taudem 5.1: guide to using the Taudem command line functions | |
CN107886573A (zh) | 一种复杂地质条件下边坡三维有限元网格生成方法 | |
CN103399234B (zh) | 一种电网功角暂态稳定在线预测方法 | |
Lin et al. | A multiscale modeling and updating framework for suspension bridges based on modal frequencies and influence lines | |
CN114386137B (zh) | 基于倾斜摄影技术的道路地形曲面优化设计方法 | |
Chen et al. | Recent Research and Applications of Numerical Simulation for Dynamic Response of Long‐Span Bridges Subjected to Multiple Loads |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160629 Termination date: 20190227 |
|
CF01 | Termination of patent right due to non-payment of annual fee |