CN107609317B - 一种钢桥疲劳评估荷载建模及疲劳评估方法 - Google Patents
一种钢桥疲劳评估荷载建模及疲劳评估方法 Download PDFInfo
- Publication number
- CN107609317B CN107609317B CN201711009863.5A CN201711009863A CN107609317B CN 107609317 B CN107609317 B CN 107609317B CN 201711009863 A CN201711009863 A CN 201711009863A CN 107609317 B CN107609317 B CN 107609317B
- Authority
- CN
- China
- Prior art keywords
- passenger
- train
- goods
- load
- variable
- 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.)
- Active
Links
- 229910000831 Steel Inorganic materials 0.000 title claims abstract description 62
- 239000010959 steel Substances 0.000 title claims abstract description 62
- 238000011156 evaluation Methods 0.000 title claims abstract description 46
- 238000001228 spectrum Methods 0.000 claims abstract description 40
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 23
- 238000000034 method Methods 0.000 claims abstract description 23
- 238000004458 analytical method Methods 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims abstract description 7
- 238000009826 distribution Methods 0.000 claims description 100
- 230000003137 locomotive effect Effects 0.000 claims description 73
- 239000013256 coordination polymer Substances 0.000 claims description 27
- 238000007619 statistical method Methods 0.000 claims description 16
- 150000001875 compounds Chemical class 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 6
- 238000009827 uniform distribution Methods 0.000 claims description 5
- 230000004044 response Effects 0.000 claims description 4
- 238000012544 monitoring process Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 230000006872 improvement Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000002902 bimodal effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种钢桥疲劳评估荷载建模及疲劳评估方法,包括:步骤1,建立列车荷载模型;步骤2,整理实测数据,并进行分类;步骤3,根据实测数据,对各参数建立概率模型;步骤4,建立列车荷载概率模型;步骤5,采用Monte‑Carlo法生成随机荷载谱;步骤6,采用有限元分析,得到钢桥疲劳细节的应力影响线;步骤7,采用随机荷载谱进行影响线加载,得到应力历程;步骤8,对应力历程进行分析,评估钢桥疲劳。本发明的有益效果:以实测数据为基础,建立疲劳评估荷载模型,描述当前通行列车荷载状况,实现对钢桥任意构件的疲劳评估,建立了完善的列车荷载建模方法,提高了实测数据的利用价值,使得钢桥疲劳评估更为准确,更接近于实际情况。
Description
技术领域
本发明涉及铁路桥梁技术领域,具体而言,涉及一种钢桥疲劳评估荷载建模及疲劳评估方法。
背景技术
随着铁路运营时间的增长,如何对铁路结构进行评估,其主要方面是确定评估荷载。特别对于铁路钢桥而言,在进行疲劳评估时,疲劳破坏特性不仅与荷载作用下的应力变化大小有关,而且还与应力变化的次数有关,因此,钢桥在进行设计时,通常采用相关规范规定的疲劳荷载谱进行验算,以控制疲劳应。
在现有的钢桥疲劳评估荷载模型中,钢构件疲劳损伤度的检算一般根据现场的实测数据推算或通过典型列车荷载计算得到,但与实际情况存在一定差别,且随着运营列车荷载和数量的增加,桥梁通行的荷载越来越复杂,已有的典型列车荷载模型已很难反映实际情况。另外,疲劳荷载一般包括实测荷载谱、标准荷载谱和疲劳荷载模型等,为得到更为准确的评价结果,多采用实测数据进行分析,但由于实测数据量庞大,应用复杂,一般都将其简化为标准的疲劳荷载谱,再采用相适应的方法生成标准疲劳荷载模型,然而如何根据实测数据建立相应的列车荷载模型,暂无统一标准,这使得依据实测数据的疲劳荷载模型建模方法不完善。
发明内容
为解决上述问题,本发明的目的在于提供一种钢桥疲劳评估荷载建模及疲劳评估方法,以实测数据为基础,建立相关疲劳评估荷载模型,描述当前通行列车荷载状况,实现对钢桥的疲劳评估,该方法建立了完善的列车荷载建模方法,提高了实测数据的利用价值,使得钢桥疲劳评估更为准确,更接近于实际情况。
本发明提供了一种钢桥疲劳评估荷载建模及疲劳评估方法,包括:
步骤1,建立列车荷载模型:
通过布设在钢桥上的传感器监测列车荷载响应及轮轨力,将钢桥上通行列车的荷载模型表示为:F=f(c,l,n,v);
式中,
参数c表示列车类型,参数l表示机车类型,参数n表示车辆编组,参数v表示列车运行速度;
步骤2,整理实测数据,并进行分类:
根据列车荷载模型包含的参数,将实测数据中各参数所包含的变量进行分类、整理;
步骤3,根据实测数据,对各参数中各变量分别建立概率模型:
对步骤2中整理分类的变量的实测数据进行统计分析,针对各个变量分别采用随机过程来拟合各变量的分布,并建立各变量的概率模型;
步骤4,建立列车荷载概率模型:
将步骤3中建立的各变量的概率模型进行组合,形成完整的列车荷载概率模型;
步骤5,采用Monte-Carlo法生成随机荷载谱:
根据所建立的列车荷载概率模型,采用Monte-Carlo法从已知的变量分布中进行抽样,生成一列列的随机列车荷载;
步骤6,采用有限元分析,得到钢桥疲劳细节的应力影响线:
在已知列车荷载的情况下,针对钢桥某部分的应力历程,采用影响线进行加载,并通过有限元分析,求得疲劳细节的应力影响线;
步骤7,采用随机荷载谱进行影响线加载,得到应力历程:
将上述生成的一列列随机列车荷载在应力影响线上加载,得到该疲劳细节的应力历程;
步骤8,对应力历程进行分析,评估钢桥疲劳:
采用雨流计算法对生成的应力历程进行处理,得到疲劳应力谱及循环次数,再根据相应的疲劳分析准则来评估钢桥疲劳状况。
作为本发明进一步的改进,参数c包含的变量为列车数量T或列车数量T及客货比α;
当通行线路为客运线路或货运线路时,参数c包含的变量为一个周期内通行的客车或货车的列车数量T;
当通行线路为客货共线线路时,参数c包含的变量为一个周期内通行的客车和货车的列车数量T和客货比α两个变量,其中,客货比α表示一个周期内通行的客车与货车的数量比。
作为本发明进一步的改进,参数l包含的变量为机车轴数β和机车轴重L;
当通行线路为客运线路时,机车轴数β包含的变量为客车机车轴数β客,机车轴重L包含的变量为客车机车轴重L客;
当通行线路为货运线路时,机车轴数β包含的变量为货车机车轴数β货,机车轴重L包含的变量为货车机车轴重L货;
当通行线路为客货共线线路时,机车轴数β包含的变量为客车机车轴数β客和货车机车轴数β货,机车轴重L包含的变量为客车机车轴重L客和货车机车轴重L货。
作为本发明进一步的改进,参数n包含的变量为车辆编组数量TC和车辆轴重LC;
当通行线路为客运线路时,车辆编组数量TC包含的变量为客车编组数量TCP,车辆轴重LC包含的变量为客车车辆轴重LCP;
当通行线路为货运线路时,车辆编组数量TC包含的变量为货车编组数量TCF,车辆轴重LC包含的变量为货车车辆轴重LCF;
当通行线路为客货共线线路时,车辆编组数量TC包含的变量为客车编组数量TCP和货车编组数量TCF,车辆轴重LC包含的变量为客车车辆轴重LCP和货车车辆轴重LCF。
作为本发明进一步的改进,参数v包含的变量为客车运行速度v客和/或货车运行速度v货;
当通行线路为客运线路时,参数v包含的变量为客车运行速度v客;
当通行线路为货运线路时,参数v包含的变量为货车运行速度v货;
当通行线路为客货共线线路时,参数v包含的变量为客车运行速度v客和货车运行速度v货。
作为本发明进一步的改进,步骤5具体包括:
步骤501,根据所建立的列车荷载概率模型中各变量分布的概率模型,采用Monte-Carlo法随机生成钢桥上一个周期内通行的列车数量T或列车数量T及客货比α;
步骤502,根据所建立的列车荷载概率模型中各变量分布的概率模型,采用Monte-Carlo法随机生成每列客车的客车机车轴数β客与客车机车轴重L客或每列货车的货车机车轴数β货与货车机车轴重L货或每列客车的客车机车轴数β客与客车机车轴重L客及每列货车的货车机车轴数β货与货车机车轴重L货;
步骤503,根据所建立的列车荷载概率模型中各变量分布的概率模型,采用Monte-Carlo法随机生成每列客车的客车编组数量TCP与客车车辆轴重LCP或每列货车的货车编组数量TCF与货车车辆轴重LCF或每列客车的客车编组数量TCP与客车车辆轴重LCP及每列货车的货车编组数量TCF与货车车辆轴重LCF;
步骤504,根据所建立的列车荷载概率模型中各变量分布的概率模型,采用Monte-Carlo法随机生成每列客车的客车运行速度v客或每列货车的货车运行速度v货或每列客车的客车运行速度v客及每列货车的货车运行速度v货;
步骤505,将上述步骤501-步骤504中随机生成的变量组合在一起,形成一列列完整的随机列车荷载。
作为本发明进一步的改进,在对各变量的分布进行拟合时,采用正态分布、0-1分布、Weibull分布、均匀分布中的一种分布或几种分布的加权组合。
本发明的有益效果为:
1、解决了钢桥疲劳评估荷载的问题,以实际通行的列车荷载为基础,根据采集到的实测数据,建立了疲劳评估荷载模型,生成的随机荷载谱更接近于实际情况,且能较好的描述当前列车荷载的通行情况,使得评估结果更为准确,能够较好的对钢桥进行疲劳评估,预测其疲劳寿命;
2、解决了列车荷载建模方法的问题,以采集到的实测数据为基础,对实测数据进行处理,建立了列车荷载概率模型,生成了随机荷载谱,形成了一套完整的基于实测数据的列车荷载模型建模方法,大大提高了实测数据的利用价值。
附图说明
图1为本发明实施例所述的一种钢桥疲劳评估荷载建模及疲劳评估方法的流程示意图;
图2下承式钢桁梁的应变片布设位置示意图;
图3为机车轴重的统计结果示意图;
图4为客车编组数量的统计结果示意图;
图5为货车编组数量的统计结果示意图;
图6为客车车辆轴重的统计结果示意图;
图7为货车车辆轴重的统计结果示意图;
图8为图1中步骤5的流程示意图;
图9为图2中下弦杆和斜腹杆上应变测试点的应力影响线;
图10为加载后的模拟应力历程示意图;
图11为下弦杆的实测日应力谱示意图;
图12为斜腹杆的实测日应力谱示意图;
图13为下弦杆的模拟日应力谱示意图;
图14为斜腹杆的模拟日应力谱示意图。
具体实施方式
下面通过具体的实施例并结合附图对本发明做进一步的详细描述。
本发明实施例所述的一种钢桥疲劳评估荷载建模及疲劳评估方法,该方法以通行车辆的实测数据为基础,实现列车荷载建模。以桥跨为75m的某下承式钢桁梁桥为例说明本发明的钢桥疲劳评估荷载建模及疲劳评估方法。
如图1所示,具体包括:
步骤1,建立列车荷载模型:
通过布设在钢桥上的传感器监测列车荷载响应及轮轨力,将钢桥上通行列车的荷载模型表示为:F=f(c,l,n,v);
式中,参数c表示列车类型,参数l表示机车类型,参数n表示车辆编组,参数v表示列车运行速度。
如图2所示,分别在一下弦杆和斜腹杆上布设应变片,监测杆件列车荷载响应,并在桥头某轨道截面位置布设轮轨力测试,监测列车通行状况。
步骤2,整理实测数据,并进行分类:
根据列车荷载模型包含的参数,将实测数据中各参数所包含的变量进行分类、整理。
钢桥上通行的列车的构成主要包括机车和车辆两部分,对桥梁结构产生影响的主要因素包括:机车轴重、机车轴数、列车类型,车辆轴重、车辆编组,以及列车运行速度等等。在铁路行车过程中,机车与车辆类型较为固定。
列车类型c这个参数主要用来表示通行的列车种类,即客车和/或货车。参数c包含的变量为列车数量T或列车数量T及客货比α。
当通行线路为客运线路或货运线路时,参数c包含的变量为一个周期内通行的客车或货车的列车数量T;
当通行线路为客货共线线路时,参数c包含的变量为一个周期内通行的客车和货车的列车数量T和客货比α两个变量,其中,客货比α表示一个周期内通行的客车与货车的数量比。
本实施例中,由于该桥为客货共线线路,且该桥每日通行的列车数量基本一致,因此,列车类型c这个参数中,包括钢桥上通行的客车与货车的客货比α和钢桥上通行的客车与货车的列车数量T两个变量。
机车类型l这个参数主要用来表示机车轴数和机车轴重,参数l包含的变量为机车轴数β和机车轴重L。
当通行线路为客运线路时,机车轴数β包含的变量为客车机车轴数β客,机车轴重L包含的变量为客车机车轴重L客;
当通行线路为货运线路时,机车轴数β包含的变量为货车机车轴数β货,机车轴重L包含的变量为货车机车轴重L货;
当通行线路为客货共线线路时,机车轴数β包含的变量为客车机车轴数β客和货车机车轴数β货,机车轴重L包含的变量为客车机车轴重L客和货车机车轴重L货。
本实施例中,由于该桥为客货共线线路,因此,需分别对客车与货车的机车进行统计分析,机车类型l包含的变量为客车机车轴数β客、货车机车轴数β货、客车机车轴重L客和货车机车轴重L货。
车辆编组n这个参数主要用来表示车辆编组数量和车辆轴重,参数n包含的变量为车辆编组数量TC和车辆轴重LC。
当通行线路为客运线路时,车辆编组数量TC包含的变量为客车编组数量TCP,车辆轴重LC包含的变量为客车车辆轴重LCP;
当通行线路为货运线路时,车辆编组数量TC包含的变量为货车编组数量TCF,车辆轴重LC包含的变量为货车车辆轴重LCF;
当通行线路为客货共线线路时,车辆编组数量TC包含的变量为客车编组数量TCP和货车编组数量TCF,车辆轴重LC包含的变量为客车车辆轴重LCP和货车车辆轴重LCF。
本实施例中,由于该桥为客货共线线路,因此,需分别对客车与货车的编组数量和轴重进行统计分析,参数n包含的变量为客车编组数量TCP、货车编组数量TCF、客车车辆轴重LCP和货车车辆轴重LCF。
列车运行速度v表示客车或货车的运行速度,参数v包含的变量为客车运行速度v客和/或货车运行速度v货。
当通行线路为客运线路时,参数v包含的变量为客车运行速度v客;
当通行线路为货运线路时,参数v包含的变量为货车运行速度v货;
当通行线路为客货共线线路时,参数v包含的变量为客车运行速度v客和货车运行速度v货。
本实施例中,由于该桥为客货共线线路,因此,需分别对客车与货车的运行速度进行分析,参数v包含的变量为客车运行速度v客和货车运行速度v货。
步骤3,根据实测数据,对各参数建立概率模型:
对步骤2中整理分类的变量的实测数据进行统计分析,针对各个变量分别采用随机过程来拟合各变量的分布,并建立各变量的概率模型。
具体的:
对该钢桥一段时间内的实测数据进行统计分析,结果显示列车数量与客车与货车的数量比均呈现单峰分布,因此,可以分别采用正态分布来拟合客车与货车的客货比α和列车总数量T的分布。
客货比α和列车总数量T的概率模型分别为:
本实施例中,根据统计分析的结果,98≤T≤124,0.9≤α≤1.59,μT=111,σT=7,μα=1.21,σα=0.13。
该钢桥客车机车以6轴机车为主,货车机车以8轴机车为主,因此,采用0-1分布拟合机车轴数β的分布,其中,客车机车轴数β客的概率模型用f(β客)表示,货车机车轴数β货的概率模型用f(β货)表示。
本实施例中,根据统计分析的结果,6轴机车的f(β客)=94.9%,8轴机车的f(β客)=95.1%;8轴机车的f(β货)=87.9%,6轴机车的f(β货)=12.1%。
客车的机车轴重和货车的机车轴重相差不大,因此,可以不区分,从图3机车轴重统计分析的结果中可以看出,机车轴重呈现三峰分布的概率模型,机车各轴的轴重主要集中在三个区域:20.05t、22.15t以及23.35t。因此,不能采用常用的单峰概率密度函数进行拟合。
本实施例中,采用2个正态分布和1个Weibull分布的加权来拟合机车轴重L的分布,其中,机车轴重L的前部分布(三峰分布前部)和尾部分布(三峰分布尾部)均分别采用1个正态分布拟合,机车轴重L的中部分布(三峰分布中部)采用Weibull分布拟合,对于3个分布的重叠处设置变量范围来加以区分。机车轴重L的概率模型为:
式中,pL1、pL2、1-pL1-pL2分别表示前部分布、中部分布和尾部分布的加权系数,μL1表示前部分布中机车轴重L的均值,σL1表示前部分布中机车轴重L的标准差,表示前部分布中机车轴重L的方差;μL2表示尾部分布中机车轴重L的均值,σL2表示尾部分布中机车轴重L的标准差,表示尾部分布中机车轴重L的方差;bL表示中部分布的形状参数,θL表示中部分布的比例参数。
本实施例中,根据统计分析的结果,pL1=0.1219,pL2=0.3022,μL1=20.085,σL1=0.358,bL=0.058,θL=22.176,μL2=23.490,σL2=0.469。其中,三个加权分布的取值区间分别为[19.02,21]、[21,22.65]以及[22.65,26.55]。
其中,由于可以不区分客车的机车轴重和货车的机车轴重,因此,客车机车轴重L客和货车机车轴重L货的分布均可以采用上述概率模型来表示。
如图4所示,为客车编组数量的统计结果,可以看出,客车编组数量主要集中在16~19之间,呈现单峰分布。因此,采用1个正态分布拟合客车编组数量TCP的分布,客车编组数量TCP的概率模型为:
如图5所示,为货车编组数量的统计结果,可以看出,货车编组数量主要集中在30和50附近,呈现双峰分布。采用2个正态分布的加权来拟合货车编组数量TCF的分布,2个分布的重叠处设置变量范围,货车编组数量TCF的概率模型为:
式中,分别表示两个正态分布的加权系数;表示第一正态分布中货车编组数量TCF的均值,表示第一正态分布中货车编组数量TCF的标准差,表示第一正态分布中货车编组数量TCF的方差;表示第二正态分布中货车编组数量TCF的均值,表示第二正态分布中货车编组数量TCF的标准差,表示第二正态分布中货车编组数量TCF的方差。
如图6所示,为客车车辆轴重的统计结果,可以看出,客车车辆轴重呈现单峰分布,主要集中在12.56t附近。因此,采用1个正态分布拟合客车车辆轴重LCP的分布,客车车辆轴重LCP的概率模型为:
如图7所示,为货车车辆轴重的统计结果,可以看出,货车荷载分布较为复杂,呈现多峰分布,前半部分各峰值相差较小,后半部分两个峰值较为突出,主要集中在19.84t以及22.19t附近。因此,采用1个均匀分布、1个Weibull分布以及1个正态分布的加权来拟合货车车辆轴重LCF的分布,3个分布的重叠处设置变量范围,货车车辆轴重LCF的概率模型为:
式中,分别表示均匀分布、Weibull分布以及正态分布的加权系数;表示Weibull分布的形状参数,表示Weibull分布的比例参数;表示货车车辆轴重LCF的均值,表示货车车辆轴重LCF的标准差,表示货车车辆轴重LCF的方差。
统计分析结果显示,客车和货车运行速度均集中在35km/h附近,都呈现单峰分布。采用正态分布拟合客车运行速度v客的分布,客车运行速度v客的概率模型为:
采用Weibull分布拟合货车运行速度v货的分布,货车运行速度v货的概率模型为:
需要说明的是,在对各变量的分布进行拟合时,并不仅限于正态分布、0-1分布、Weibull分布、均匀分布中的一种分布或几种分布的加权组合,可以根据实测数据的统计分析后,对不同的变量各自选择合适的分布进行拟合。
步骤4,建立列车荷载概率模型:
将步骤3中建立的各变量的概率模型进行组合,形成完整的列车荷载概率模型。
步骤5,采用Monte-Carlo法生成随机荷载谱:
根据所建立的列车荷载概率模型,采用Monte-Carlo法从已知的变量分布中进行抽样,生成一列列的随机列车荷载。
Monte-Carlo法称为统计模拟法,随机抽样技术,是以概率和统计理论方法为基础的一种计算方法。其基本思想是在样本量足够大的情况下事件发生的“频率”可等效为“概率”,基于Monte-Carlo法的随机列车荷载谱生成主要分为两个步骤:1)建立各参数分布的概率模型,如上所述;2)从所建立的概率模型中进行抽样,先产生一系列随机数,再根据随机数从已知分布进行抽样。
如图8所示,具体包括:
步骤501,根据所建立的列车荷载概率模型中各变量分布的概率模型,采用Monte-Carlo法随机生成钢桥上一个周期内通行的列车数量T或列车数量T及客货比α;
步骤502,根据所建立的列车荷载概率模型中各变量分布的概率模型,采用Monte-Carlo法随机生成每列客车的客车机车轴数β客与客车机车轴重L客或每列货车的货车机车轴数β货与货车机车轴重L货或每列客车的客车机车轴数β客与客车机车轴重L客及每列货车的货车机车轴数β货与货车机车轴重L货;
步骤503,根据所建立的列车荷载概率模型中各变量分布的概率模型,采用Monte-Carlo法随机生成每列客车的客车编组数量TCP与客车车辆轴重LCP或每列货车的货车编组数量TCF与货车车辆轴重LCF或每列客车的客车编组数量TCP与客车车辆轴重LCP及每列货车的货车编组数量TCF与货车车辆轴重LCF;
步骤504,根据所建立的列车荷载概率模型中各变量分布的概率模型,采用Monte-Carlo法随机生成每列客车的客车运行速度v客或每列货车的货车运行速度v货或每列客车的客车运行速度v客及每列货车的货车运行速度v货;
步骤505,将上述步骤501-步骤504中随机生成的变量组合在一起,形成一列列完整的随机列车荷载。
本实施例中,由于该桥为客货共线线路,因此,需要采用Monte-Carlo法随机生成以下变量:
钢桥上一个周期内通行的客车和货车列车数量T及客货比α;
每列客车的客车机车轴数β客与客车机车轴重L客及每列货车的货车机车轴数β货与货车机车轴重L货;
每列客车的客车编组数量TCP与客车车辆轴重LCP及每列货车的货车编组数量TCF与货车车辆轴重LCF;
每列客车的客车运行速度v客及每列货车的货车运行速度v货。
步骤6,采用有限元分析,得到钢桥疲劳细节的应力影响线:
在已知列车荷载的情况下,针对钢桥某部分的应力历程,采用影响线进行加载,并通过MIDAS/Civil建模进行有限元分析,求得疲劳细节的应力影响线。
图9为本实施例中下弦杆和斜腹杆上应变测试点的应力影响线。
步骤7,采用随机荷载谱进行影响线加载,得到应力历程:
将上述生成的一列列随机列车荷载在应力影响线上加载,得到该疲劳细节的应力历程。
图10为本实施例中加载后的模拟应力历程示意图。
步骤8,对应力历程进行分析,评估钢桥疲劳:
采用雨流计算法对生成的应力历程进行处理,得到疲劳应力谱及循环次数,再根据相应的疲劳分析准则来评估钢桥疲劳状况。
雨流计算法的计算规则为:
(1)雨流依次从载荷时间历程的峰值位置的内侧沿着斜坡往下流;
(2)雨流从某一个峰值点开始流动,当遇到比其起始峰值更大的峰值时要停止流动;
(3)雨流遇到上面流下的雨流时,必须停止流动;
(4)取出所有的全循环,记下每个循环的幅度;
(5)将第一阶段计数后剩下的发散收敛载荷时间历程等效为一个收敛发散型的载荷时间历程,进行第二阶段的雨流计数。计数循环的总数等于两个计数阶段的计数循环之和。
雨流计算法的要点是载荷-时间历程的每一部分都参与计数,且只计数一次,一个大的幅值所引起的损伤不受截断它的小循环的影响,截出的小循环迭加到较大的循环和半循环上去。再根据累积损伤理论结合相应的S-N曲线就可对钢构件的疲劳寿命进行评估。
如图11~图14所示,分别为本实施例一天内的实测应力谱与模拟应力谱。实测应力谱根据实测数据进行统计得到,模拟应力谱采用雨流计算法对步骤7中加载的模拟应力历程进行计算得到。可以看出,下弦杆实测应力谱主要集中于10MPa以下以及20MPa左右的应力循环,模拟应力谱主要集中于10MPa以下、25MPa以及35MPa左右的应力循环;斜杆实测应力谱主要集中于15MPa以下、25MPa以及40MPa左右的应力循环,模拟应力谱主要于15MPa以下、25MPa、40MPa以及50MPa左右的应力循环;对于低应力循环,实测谱循环次数大于模拟谱,主要是由于未能充分剔除干扰对应力历程的干扰,导致低应力循环次数偏大,对于中高幅应力循环,实测谱循环次数与模拟谱循环次数接近,总的来看,模拟应力谱与实测应力谱吻合良好。
从上述对荷载的疲劳效应分析可以看出,采用该方法建立的列车荷载模型能很好的反映当前列车荷载的通行状况,模拟应力谱与实测应力谱吻合良好,说明基于概率论的列车荷载建模方法具有良好的适用性,可用作钢桥疲劳分析,提高了实测数据的利用价值,使得钢桥疲劳评估更为准确,更接近于实际情况。
本发明的方法将列车荷载模型分成四个参数来表示,分别包括列车类型,机车类型,车辆编组以及列车运行速度,每个参数根据线路情况的不同包含不同的变量,这些变量的运用使得建立的车辆荷载模型与实际情况一致,提高了建模的准确性。列车荷载模型中各参数都根据实测数据以概率论为基础建立参数中各变量的概率分布,采用一随机过程来描述通行的列车荷载,能更好的描述当前铁路通行列车荷载的复杂性,更加反映列车通行的实际情况。采用Monte-Carlo法根据列车荷载模型生成随机列车荷载谱,形成了一套完整的基于实测数据的列车荷载模型建模方法及列车荷载谱模拟方法,大大提高了实测数据的利用价值,使得钢桥疲劳评估更为准确。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种钢桥疲劳评估荷载建模及疲劳评估方法,其特征在于,包括:
步骤1,建立列车荷载模型:
通过布设在钢桥上的传感器监测列车荷载响应及轮轨力,将钢桥上通行列车的荷载模型表示为:F=f(c,l,n,v);
式中,
参数c表示列车类型,参数l表示机车类型,参数n表示车辆编组,参数v表示列车运行速度;
步骤2,整理实测数据,并进行分类:
根据列车荷载模型包含的参数,将实测数据中各参数所包含的变量进行分类、整理;
步骤3,根据实测数据,对各参数中各变量分别建立概率模型:
对步骤2中整理分类的变量的实测数据进行统计分析,针对各个变量分别采用随机过程来拟合各变量的分布,并建立各变量的概率模型;
步骤4,建立列车荷载概率模型:
将步骤3中建立的各变量的概率模型进行组合,形成完整的列车荷载概率模型;
步骤5,采用Monte-Carlo法生成随机荷载谱:
根据所建立的列车荷载概率模型,采用Monte-Carlo法从已知的变量分布中进行抽样,生成一列列的随机列车荷载;
步骤6,采用有限元分析,得到钢桥疲劳细节的应力影响线:
在已知列车荷载的情况下,针对钢桥的应力历程,采用影响线进行加载,并通过有限元分析,求得疲劳细节的应力影响线;
步骤7,采用随机荷载谱进行影响线加载,得到应力历程:
将上述生成的一列列随机列车荷载在应力影响线上加载,得到该疲劳细节的应力历程;
步骤8,对应力历程进行分析,评估钢桥疲劳:
采用雨流计算法对生成的应力历程进行处理,得到疲劳应力谱及循环次数,再根据相应的疲劳分析准则来评估钢桥疲劳状况。
2.根据权利要求1所述的钢桥疲劳评估荷载建模及疲劳评估方法,其特征在于,
参数c包含的变量为列车数量T或列车数量T及客货比α;
当通行线路为客运线路或货运线路时,参数c包含的变量为一个周期内通行的客车或货车的列车数量T;
当通行线路为客货共线线路时,参数c包含的变量为一个周期内通行的客车和货车的列车数量T和客货比α两个变量,其中,客货比α表示一个周期内通行的客车与货车的数量比。
3.根据权利要求1所述的钢桥疲劳评估荷载建模及疲劳评估方法,其特征在于,
参数l包含的变量为机车轴数β和机车轴重L;
当通行线路为客运线路时,机车轴数β包含的变量为客车机车轴数β客,机车轴重L包含的变量为客车机车轴重L客;
当通行线路为货运线路时,机车轴数β包含的变量为货车机车轴数β货,机车轴重L包含的变量为货车机车轴重L货;
当通行线路为客货共线线路时,机车轴数β包含的变量为客车机车轴数β客和货车机车轴数β货,机车轴重L包含的变量为客车机车轴重L客和货车机车轴重L货。
4.根据权利要求1所述的钢桥疲劳评估荷载建模及疲劳评估方法,其特征在于,
参数n包含的变量为车辆编组数量TC和车辆轴重LC;
当通行线路为客运线路时,车辆编组数量TC包含的变量为客车编组数量TCP,车辆轴重LC包含的变量为客车车辆轴重LCP;
当通行线路为货运线路时,车辆编组数量TC包含的变量为货车编组数量TCF,车辆轴重LC包含的变量为货车车辆轴重LCF;
当通行线路为客货共线线路时,车辆编组数量TC包含的变量为客车编组数量TCP和货车编组数量TCF,车辆轴重LC包含的变量为客车车辆轴重LCP和货车车辆轴重LCF。
5.根据权利要求1所述的钢桥疲劳评估荷载建模及疲劳评估方法,其特征在于,
参数v包含的变量为客车运行速度v客和/或货车运行速度v货;
当通行线路为客运线路时,参数v包含的变量为客车运行速度v客;
当通行线路为货运线路时,参数v包含的变量为货车运行速度v货;
当通行线路为客货共线线路时,参数v包含的变量为客车运行速度v客和货车运行速度v货。
6.根据权利要求1所述的钢桥疲劳评估荷载建模及疲劳评估方法,其特征在于,步骤5具体包括:
步骤501,根据所建立的列车荷载概率模型中各变量分布的概率模型,采用Monte-Carlo法随机生成钢桥上一个周期内通行的列车数量T或列车数量T及客货比α;
步骤502,根据所建立的列车荷载概率模型中各变量分布的概率模型,采用Monte-Carlo法随机生成每列客车的客车机车轴数β客与客车机车轴重L客或每列货车的货车机车轴数β货与货车机车轴重L货或每列客车的客车机车轴数β客与客车机车轴重L客及每列货车的货车机车轴数β货与货车机车轴重L货;
步骤503,根据所建立的列车荷载概率模型中各变量分布的概率模型,采用Monte-Carlo法随机生成每列客车的客车编组数量TCP与客车车辆轴重LCP或每列货车的货车编组数量TCF与货车车辆轴重LCF或每列客车的客车编组数量TCP与客车车辆轴重LCP及每列货车的货车编组数量TCF与货车车辆轴重LCF;
步骤504,根据所建立的列车荷载概率模型中各变量分布的概率模型,采用Monte-Carlo法随机生成每列客车的客车运行速度v客或每列货车的货车运行速度v货或每列客车的客车运行速度v客及每列货车的货车运行速度v货;
步骤505,将上述步骤501-步骤504中随机生成的变量组合在一起,形成一列列完整的随机列车荷载。
7.根据权利要求1所述的钢桥疲劳评估荷载建模及疲劳评估方法,其特征在于,在对各变量的分布进行拟合时,采用正态分布、0-1分布、Weibull分布、均匀分布中的一种分布或几种分布的加权组合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711009863.5A CN107609317B (zh) | 2017-10-25 | 2017-10-25 | 一种钢桥疲劳评估荷载建模及疲劳评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711009863.5A CN107609317B (zh) | 2017-10-25 | 2017-10-25 | 一种钢桥疲劳评估荷载建模及疲劳评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107609317A CN107609317A (zh) | 2018-01-19 |
CN107609317B true CN107609317B (zh) | 2020-11-17 |
Family
ID=61079842
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711009863.5A Active CN107609317B (zh) | 2017-10-25 | 2017-10-25 | 一种钢桥疲劳评估荷载建模及疲劳评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107609317B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108170991B (zh) * | 2018-01-24 | 2021-03-16 | 西南交通大学 | 基于Weibull分布的综合应力恒加试验方案可靠性评价方法 |
CN108363892A (zh) * | 2018-04-26 | 2018-08-03 | 扬州大学 | 随机车辆荷载作用下应力谱编制方法 |
CN109991092B (zh) * | 2019-03-18 | 2022-08-19 | 华南理工大学 | 变幅荷载下构件的自然暴露实验装置、系统及方法 |
CN111985084B (zh) * | 2020-07-24 | 2023-10-27 | 华设设计集团股份有限公司 | 基于有限变量的桥梁力学实时分析方法 |
CN112347668B (zh) * | 2020-09-29 | 2022-04-12 | 华东交通大学 | 基于概率断裂力学的钢桥面板疲劳可靠度评估方法 |
CN115408754B (zh) * | 2022-09-05 | 2024-03-12 | 北京交通大学 | 一种考虑时变效应的组合梁桥动力疲劳性能分析方法 |
CN116011195B (zh) * | 2022-12-19 | 2024-03-15 | 广东水电二局股份有限公司 | 一种基于数据驱动的风机基础疲劳荷载谱建模方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103279588A (zh) * | 2013-04-09 | 2013-09-04 | 东南大学 | 车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法 |
CN106055784A (zh) * | 2016-05-30 | 2016-10-26 | 东南大学 | 一种钢桥细节疲劳裂纹扩展评估方法 |
CN107092735A (zh) * | 2017-04-07 | 2017-08-25 | 东南大学 | 一种桥梁疲劳状态评估方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8761000B2 (en) * | 2010-11-29 | 2014-06-24 | Edge Velocity Corporation | Router and rapid response network |
-
2017
- 2017-10-25 CN CN201711009863.5A patent/CN107609317B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103279588A (zh) * | 2013-04-09 | 2013-09-04 | 东南大学 | 车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法 |
CN106055784A (zh) * | 2016-05-30 | 2016-10-26 | 东南大学 | 一种钢桥细节疲劳裂纹扩展评估方法 |
CN107092735A (zh) * | 2017-04-07 | 2017-08-25 | 东南大学 | 一种桥梁疲劳状态评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107609317A (zh) | 2018-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107609317B (zh) | 一种钢桥疲劳评估荷载建模及疲劳评估方法 | |
Milković et al. | Wayside system for wheel–rail contact forces measurements | |
CN105718633B (zh) | 一种底盘件的载荷分析方法 | |
Lu et al. | Fatigue reliability evaluation of orthotropic steel bridge decks based on site-specific weigh-in-motion measurements | |
CN109002622B (zh) | 一种随机车流作用下大跨径桥梁总体荷载响应估算方法 | |
CN107107927A (zh) | 铁道车辆状态监视装置 | |
Funfschilling et al. | Uncertainty quantification in vehicle dynamics | |
Yang et al. | Methodology to evaluate fatigue damage of high-speed train welded bogie frames based on on-track dynamic stress test data | |
Banerji et al. | Structural health monitoring of a steel railway bridge for increased axle loads | |
Wang et al. | A diagnostic method of freight wagons hunting performance based on wayside hunting detection system | |
Brzozowski et al. | A weigh-in-motion system with automatic data reliability estimation | |
Kim et al. | Calculation of resistance to motion of a high-speed train using acceleration measurements in irregular coasting conditions | |
Manouchehrynia et al. | Reliability-based fatigue life assessment using random road load condition and local damage criteria | |
CN113515810A (zh) | 基于可靠性和安全性分析的动车组转向架设计开发方法 | |
Figlus et al. | Selection of a semi-trailer for the haulage of long oversize loads, taking into account an analysis of operational damage | |
CN116046302B (zh) | 基于应变时程曲线识别装配式梁桥损伤的方法 | |
CN110853164A (zh) | 基于路网损伤的通行收费方法、系统、介质及收费设备 | |
Spiryagin et al. | Comparison of simplified and complete contact modelling approaches in simulations of high adhesion locomotives | |
AT518759A1 (de) | Verfahren und Vorrichtung zur Überwachung zumindest einer im Bahnbau verlegten Fahrwegkomponente | |
Luo et al. | Fatigue damage evaluation for a railway vehicle bogie using appropriate sampling frequencies | |
Milković et al. | Experimental approach to assessment of safety against derailment of freight wagons | |
Gao | Load rating highway bridges in the United States: The state of practice | |
Nurmikolu et al. | Statistical analysis of wheel impact load data and review for Finnish impact load limits | |
CN115791033B (zh) | 基于通行车辆荷载及动挠度监测的桥梁承载能力评定方法 | |
Parkinson et al. | An intelligent track monitoring system |
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 |