CN103279588B - 车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法 - Google Patents
车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法 Download PDFInfo
- Publication number
- CN103279588B CN103279588B CN201310121658.3A CN201310121658A CN103279588B CN 103279588 B CN103279588 B CN 103279588B CN 201310121658 A CN201310121658 A CN 201310121658A CN 103279588 B CN103279588 B CN 103279588B
- Authority
- CN
- China
- Prior art keywords
- fatigue
- time
- fatigue stress
- load
- bridge deck
- 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.)
- Expired - Fee Related
Links
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法,包括如下步骤:建立车辆荷载与温度共同作用下的钢桥面板-铺装层一体化疲劳分析模型;选取单辆标准疲劳车,加载于上述的一体化疲劳分析模型,计算预定工况下焊接细节的疲劳应力时程曲线,在此基础上建立任一工况下焊接细节疲劳应力时程计算的分析模型,确定车辆荷载与温度共同作用下的焊接细节疲劳应力时程。该计算方法可实现任一时间段内的钢桥面板疲劳应力时程模拟,且模拟结果能准确反映实际运营状态下的钢桥面板疲劳应力变化规律,填补了现有技术中没有合适的车辆荷载与温度共同作用下钢桥面板疲劳应力计算方法的空白。
Description
技术领域
本发明属于大跨缆索承重桥梁结构钢桥面板的疲劳分析与设计领域,具体涉及一种车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法。
背景技术
钢桥面板疲劳性能的影响因素有桥梁结构所处的环境以及外界荷载作用,其中车辆荷载作为主要影响因素,在钢桥面板疲劳问题的研究领域中受到重点关注,例如林日常等人基于实测车辆荷载数据建立了全国各地的标准疲劳车模型,在此基础上经分析得到了标准疲劳车辆模型的整车车重、轴重分配率、轴长、轴距分配等参数对钢桥面板疲劳损伤和疲劳寿命的影响规律;此外,研究表明环境温度对钢桥面板的疲劳性能也会产生影响,其影响源于沥青铺装层具有温度敏感性,铺装层温度的变化改变了铺装层的弹性模量,进而改变了车辆荷载在铺装层中的传递压力,从而间接影响了钢桥面板的疲劳性能。然而,目前国内外学者一直致力于解决钢桥面板在车辆荷载持续作用下的疲劳问题,而忽略温度作用的影响;同时,大多数钢桥面板疲劳问题的研究仅局限在有限时间长度内,考虑到钢桥面板的疲劳损伤是一个不断累积的过程,有限时间长度内的疲劳计算结果是难以满足分析要求的。
因此,需要一种新的车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法以解决上述问题。
发明内容
发明目的:本发明针对现有技术中钢桥面板在车辆荷载持续作用下的疲劳问题未考虑温度作用的影响的缺陷,提供一种车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法。
技术方案:为解决上述技术问题,本发明的车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法采用如下技术方案:
一种车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法,包括以下步骤:
1)、建立车辆荷载与温度共同作用下的钢桥面板-铺装层一体化疲劳分析模型;
2)、选取单辆标准疲劳车,加载于步骤1)建立的所述钢桥面板-铺装层一体化疲劳分析模型,计算预定工况下焊接细节的疲劳应力时程曲线,在此基础上建立任一工况下焊接细节疲劳应力时程计算的分析模型;
3)、选取k辆标准疲劳车进行随机模拟,其中,k≥2,确定各辆标准疲劳车在时间轴上的起始时刻t1,m、结束时刻t2,m以及在时间段[t1,m,t2,m]内的行驶车速vm,其中m=1,……,k;
4)、对沥青铺装层的温度值进行随机模拟,确定步骤3)时间段[t1,m,t2,m]内对应的铺装层温度值T;
5)、将时间段[t1,m,t2,m]内的行驶车速vm以及铺装层温度值T代入步骤2)的分析模型,计算每个时间段各个时刻对应的疲劳应力值,并将不同时间段同一时刻的疲劳应力值线性叠加,确定车辆荷载与温度共同作用下的焊接细节疲劳应力时程。
更进一步的,步骤1)所述的钢桥面板-铺装层一体化疲劳分析模型通过以下方法得到:
基于ANSYS有限元分析软件,建立车辆荷载与温度共同作用下的钢桥面板-铺装层一体化疲劳分析模型:
①钢桥面板部分由顶板和纵肋组成,采用壳单元进行模拟;
②沥青铺装层采用实体单元模拟,其中对于沥青铺装层弹性模量的模拟,利用沥青铺装层的弹性模量E与铺装层温度T之间的相关性,采用下式表示:
E=1.2×100.01693(20-T)+3。
更进一步的,钢桥面板的横桥向尺寸取6个U形肋,顺桥向尺寸取5个横隔板间距,四周边界约束模拟为铰接支撑,中间横隔板对钢桥面板的边界约束模拟为刚度大小为1.88×107N/m的竖向弹簧支撑。
更进一步的,步骤2)中选取单辆标准疲劳车,加载于步骤1)建立的所述钢桥面板-铺装层一体化疲劳分析模型,计算预定工况下焊接细节的疲劳应力时程曲线,在此基础上建立任一工况下焊接细节疲劳应力时程计算的分析模型,具体包括以下步骤:
(a)取单辆标准疲劳车荷载模型,以行驶车速v0沿车道中心线驶过步骤1)建立的一体化疲劳分析模型,根据沥青铺装层温度T的不同分为n种不同的计算工况,计算钢桥面板焊接细节的疲劳应力时程曲线S1(t)~Sn(t),用集合表示为S(t)={S1(t),……,Sn(t)},其中t为时间,Sr(t)为第r种计算工况下的疲劳应力时程曲线,r=1,……,n;
(b)根据焊接细节的n条疲劳应力时程曲线S(t),选取其在t0时刻的疲劳应力值S(t0)及其对应的沥青铺装层温度值T,绘制S(t0)-T散点图并利用最小二乘法进行二次多项式函数拟合,建立S(t0)-T相关模型S(t0,T),据此可得到t0时刻任一铺装层温度T时的疲劳应力值,根据此方法对S(t)的所有时刻遍历可得到任一时刻t、任一铺装层温度T时的疲劳应力值,用函数S(t,T)表示;
(c)对函数S(t,T)的时间变量t乘以比例系数v/v0,用函数S(vt/v0,T)表示,据此可得到自变量v、t、T在给定值时的焊接细节疲劳应力值,其中函数S(tv/v0,T)可作为焊接细节疲劳应力时程计算的分析模型。
更进一步的,步骤(a)中沥青铺装层温度T分为-10℃、0℃、10℃、20℃、30℃、40℃、50℃和60℃八种不同情况。
更进一步的,步骤(a)中所述行驶车速v0=60km/h。
更进一步的,步骤3)中确定各辆标准疲劳车在时间轴上的起始时刻t1,m、结束时刻t2,m以及在时间段[t1,m,t2,m]内的行驶车速vm,具体包括以下步骤:
(a1)、确定标准疲劳车的行车间距s和行驶车速v的概率密度统计模型f(s)和g(v):
标准疲劳车行车间距s的概率密度统计模型f(s)服从对数正态分布,用公式表达为:
式中,η、γ为待估参数,利用Matlab统计工具箱中的对数正态分布函数,对行车间距采集样本进行概率密度统计特性分析,确定η、γ参数估计值,
标准疲劳车行驶车速v的概率密度统计模型g(v)采用下式表达:
g(v)=f(s)p(tp)
式中,tp为某辆标准疲劳车最后一车轴经过观测点至后一辆标准疲劳车第一车轴经过观测点的所用时间,p(tp)为时间tp的概率密度统计模型,服从指数分布,其中ρ为p(tp)的待估参数,利用Matlab统计工具箱中的指数分布函数,对时间tp的采集样本进行概率密度统计特性分析,确定ρ参数估计值;
(b1)、对概率密度统计模型f(s)和g(v)进行蒙特卡罗抽样,确定行车间距s和行驶车速v在时间轴上的样本序列S和V:
利用隐式蒙特卡罗抽样公式,对标准疲劳车行车间距s的概率密度统计模型f(s)进行均匀随机抽样,得到标准疲劳车行车间距s的样本序列S:
F(S)=F(Sa)+[F(Sb)-F(Sa)]·rand(0,1)
式中,Sa、Sb分别为样本序列S的上下限;F(S)、F(Sa)、F(Sb)分别为f(s)在区间(-∞,S]、(-∞,Sa]、(-∞,Sb]上的积分值;[F(Sb)-F(Sa)]·rand(0,1)为在区间(0,F(Sb)-F(Sa)]上的均匀随机抽样,其中k-1为样本序列S的容量;
利用隐式蒙特卡罗抽样公式,对标准疲劳车行驶车速v的概率密度统计模型g(v)进行均匀随机抽样,得到标准疲劳车行驶车速v的样本序列V:
F(V)=F(Va)+[F(Vb)-F(Va)]·rand(0,1)
式中,Va、Vb分别为样本序列V的上下限;F(V)、F(Va)、F(Vb)分别为f(v)在区间(-∞,V]、(-∞,Va]、(-∞,Vb]上的积分值;[F(Vb)-F(Va)]·rand(0,1)为在区间(0,F(Vb)-F(Va)]上的均匀随机抽样,其中k为样本序列V的容量;
(c1)基于样本序列S和V,对k辆标准疲劳车进行随机模拟:假定单辆标准疲劳车的各轴距之和为l1,一体化分析模型的顺桥向长度为l2,若在t1,m时刻,第m辆标准疲劳车第一车轴以行驶车速vm(位于样本序列V的第m列)驶入一体化分析模型,则经历时间l1/vm后,第m辆标准疲劳车的最后一轴会经过观测点,此后再等待时间间隔△tm(△tm=Sm/vm+1,Sm位于样本序列S的第m列)后,第m+1辆标准疲劳车第一车轴会以行驶车速vm+1驶入一体化分析模型,据此分析可得到第m辆标准疲劳车在时间轴上的结束时刻t2,m=t1,m+l1/vm+l2/vm,第m+1辆标准疲劳车在时间轴上的起始时刻t1,m+1=t1, m+l1/vm+△tm;依次类推,便可根据抽样序列确定k辆标准疲劳车在时间轴上的起始时刻t1,m和结束时刻t2,m(m=1,……,k),以及在时间段[t1,m,t2,m]内的行驶车速vm。
模拟结果能准确反映实际运营状态下的钢桥面板疲劳应力变化规律:①在对随机变量s、v的概率密度统计模型进行蒙特卡罗抽样时,通过调整抽样个数k-1、k控制最后一辆标准疲劳车在时间轴上的结束时刻t2,k,可以实现任一时间长度内的钢桥面板疲劳应力时程模拟;②随机变量s、v、T的概率密度统计模型是基于实测样本的概率密度统计特性得到的,因此其模拟样本序列也具有实测样本的概率密度统计特征。
更进一步的,步骤5)中将时间段[t1,m,t2,m]内的行驶车速vm以及铺装层温度值T代入步骤2)的分析模型,计算每个时间段各个时刻对应的疲劳应力值,并将不同时间段同一时刻的疲劳应力值线性叠加,确定车辆荷载与温度共同作用下的焊接细节疲劳应力时程,具体包括以下步骤:
(a2)、基于步骤3)和步骤4),根据第m个时间段[t1,m,t2,m]中tm时刻对应的行驶车速vm(tm)以及铺装层温度值T(tm),令t=tm-t1,m、v=vm(tm)、T=T(tm),代入步骤2)的分析模型,确定第m辆标准疲劳车在时间轴上tm时刻的疲劳应力值,对所有时间段的时刻遍历确定k辆标准疲劳车在时间轴上的疲劳应力时程;
(b2)、将k辆标准疲劳车的疲劳应力时程在时间轴上同一时刻的疲劳应力值线性叠加,得到疲劳应力值在时间轴上的模拟序列,从而确定车辆荷载与温度共同作用下的钢桥面板疲劳应力时程。
有益效果:本发明的车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法,可以实现任一时间长度内的钢桥面板疲劳应力时程模拟,基于模拟样本序列得到的疲劳应力时程可以准确反映实际运营状态下的钢桥面板疲劳应力变化规律。该计算方法填补了现有技术中没有合适的车辆荷载与温度共同作用下钢桥面板疲劳应力计算方法的空白,可得到广泛推广和应用。
附图说明
图1为本发明实施例润扬大桥斜拉桥钢桥面板-铺装层一体化分析模型;
图2为本发明实施例标准疲劳车荷载模型的主视图;
图3为本发明实施例标准疲劳车荷载模型的俯视图;
图4为本发明实施例单侧轮压顺桥向布置及疲劳应力监测部位;
图5为本发明实施例单侧轮压横桥向布置及疲劳应力监测部位;
图6为本发明实施例不同铺装层温度时测点A的疲劳应力时程;
图7为本发明实施例不同铺装层温度时测点B的疲劳应力时程;
图8为本发明实施例测点A、B的S(t0)-T相关模型散点图及估计曲线;
图9为本发明实施例分析模型S(vt/60,T)在v、T给定取值时的疲劳应力时程;
图10为本发明实施例随机变量s模拟样本的概率密度统计特性;
图11为本发明实施例随机变量t模拟样本的概率密度统计特性;
图12为本发明实施例温度模拟样本Tyear调整后的时程序列;
图13为本发明实施例温度模拟样本Tyear的概率密度统计特性;
图14为本发明实施例车辆荷载与温度共同作用下钢桥面板疲劳应力时程。
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
本发明的车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法,该计算方法包括如下步骤:
步骤1)建立车辆荷载与温度共同作用下的钢桥面板-铺装层一体化疲劳分析模型:
基于ANSYS有限元分析软件,根据大跨缆索承重桥梁结构体系的设计图纸和资料,建立车辆荷载与温度共同作用下的钢桥面板-铺装层一体化疲劳分析模型。①其中钢桥面板部分由顶板和纵肋组成,采用壳单元进行模拟,钢桥面板的横桥向尺寸取6个U形肋,顺桥向尺寸取5个横隔板间距,四周边界约束模拟为铰接支撑,中间横隔板对钢桥面板的边界约束模拟为刚度大小为1.88×107N/m的竖向弹簧支撑;②沥青铺装层采用实体单元模拟,其中对于沥青铺装层弹性模量的模拟,利用沥青铺装层的弹性模量E与铺装层温度T之间的相关性,采用式(1)表示:
E=1.2×100.01693(20-T)+3(1)
步骤2)选取单辆标准疲劳车,加载于步骤1)建立的一体化疲劳分析模型,计算特定工况(给定行驶车速v和铺装层温度T)下焊接细节的疲劳应力时程曲线,在此基础上建立任一工况下焊接细节疲劳应力时程计算的分析模型:
(a)取单辆标准疲劳车荷载模型,以行驶车速v=60km/h沿车道中心线驶过步骤10)建立的一体化疲劳分析模型,根据沥青铺装层温度T(-10℃、0℃、10℃、20℃、30℃、40℃、50℃、60℃)的不同分为八种计算工况,计算钢桥面板焊接细节的疲劳应力时程曲线S1(t)~S8(t),用集合表示为S(t)={S1(t),……,S8(t)},其中t为时间,Sr(t)为第r种计算工况下的疲劳应力时程曲线,r=1,……,8;
(b)根据焊接细节的八条疲劳应力时程曲线S(t),选取其在t0时刻的疲劳应力值S(t0)及其对应的沥青铺装层温度值T,绘制S(t0)-T散点图并利用最小二乘法进行二次多项式函数拟合,建立S(t0)-T相关模型S(t0,T),据此可得到t0时刻任一铺装层温度T时的疲劳应力值,根据此方法对S(t)的所有时刻遍历可得到任一时刻t、任一铺装层温度T时的疲劳应力值,用函数S(t,T)表示;
(c)对函数S(t,T)的时间变量t乘以比例系数v/60,用函数S(vt/60,T)表示,据此可得到自变量v、t、T在给定值时的焊接细节疲劳应力值,其中函数S(tv/60,T)可作为焊接细节疲劳应力时程计算的分析模型;
步骤3)选取k(k≥2)辆标准疲劳车进行随机模拟,确定各辆标准疲劳车在时间轴上的起始时刻t1,m和结束时刻t2,m(m=1,……,k),以及在时间段[t1,m,t2,m]内的行驶车速vm:
(a1)确定标准疲劳车的行车间距s和行驶车速v的概率密度统计模型f(s)和g(v):
标准疲劳车行车间距s的概率密度统计模型f(s)服从对数正态分布,用公式表达为:
式中,η、γ为待估参数,利用Matlab统计工具箱中的对数正态分布函数,对行车间距采集样本进行概率密度统计特性分析,确定η、γ参数估计值。标准疲劳车行驶车速v的概率密度统计模型g(v)采用下式表达:
g(v)=f(s)p(tp)(3a)
式中,tp为某辆标准疲劳车最后一车轴经过观测点至后一辆标准疲劳车第一车轴经过观测点的所用时间,p(tp)为时间tp的概率密度统计模型,服从指数分布,其中ρ为p(tp)的待估参数,利用Matlab统计工具箱中的指数分布函数,对时间tp的采集样本进行概率密度统计特性分析,确定ρ参数估计值;
(b1)对概率密度统计模型f(s)和g(v)进行蒙特卡罗抽样,确定行车间距s和行驶车速v在时间轴上的样本序列S和V:
利用隐式蒙特卡罗抽样公式(4),对式(2)的概率密度统计模型f(s)进行均匀随机抽样,得到标准疲劳车行车间距s的样本序列S:
F(S)=F(Sa)+[F(Sb)-F(Sa)]·rand(0,1)(4)
式(4)中,Sa、Sb分别为样本序列S的上下限,分别取为13.2m和1000m;F(S)、F(Sa)、F(Sb)分别为f(s)在区间(-∞,S]、(-∞,Sa]、(-∞,Sb]上的积分值,利用Matlab数学软件中的积分函数命令vpa(int(f(s),-∞,S))、vpa(int(f(s),-∞,Sa))、vpa(int(f(s),-∞,Sb))得到;[F(Sb)-F(Sa)]·rand(0,1)为在区间(0,F(Sb)-F(Sa)]上的均匀随机抽样,利用Matlab数学软件中的抽样函数命令(F(Sb)-F(Sa))*rand(1,k-1)得到,其中k-1为样本序列S的容量;
利用隐式蒙特卡罗抽样公式(5),对式(3a)的概率密度统计模型g(v)进行均匀随机抽样,得到标准疲劳车行驶车速v的样本序列V:
F(V)=F(Va)+[F(Vb)-F(Va)]·rand(0,1)(5)
式(5)中,Va、Vb分别为样本序列V的上下限,分别取为13.8m/s和30.6m/s;F(V)、F(Va)、F(Vb)分别为f(v)在区间(-∞,V]、(-∞,Va]、(-∞,Vb]上的积分值,利用Matlab数学软件中的积分函数命令vpa(int(f(v),-∞,V))、vpa(int(f(v),-∞,Va))、vpa(int(f(v),-∞,Vb))得到;[F(Vb)-F(Va)]·rand(0,1)为在区间(0,F(Vb)-F(Va)]上的均匀随机抽样,利用Matlab数学软件中的抽样函数命令(F(Vb)-F(Va))*rand(1,k)实现,其中k为样本序列V的容量;
(c1)基于样本序列S和V,对k辆标准疲劳车进行随机模拟:
假定单辆标准疲劳车的各轴距之和为l1,一体化分析模型的顺桥向长度为l2,若在t1,m时刻,第m辆标准疲劳车第一车轴以行驶车速vm(位于样本序列V的第m列)驶入一体化分析模型,则经历时间l1/vm后,第m辆标准疲劳车的最后一轴会经过观测点,此后再等待时间间隔△tm(△tm=Sm/vm+1,Sm位于样本序列S的第m列)后,第m+1辆标准疲劳车第一车轴会以行驶车速vm+1驶入一体化分析模型,据此分析可得到第m辆标准疲劳车在时间轴上的结束时刻t2,m=t1,m+l1/vm+l2/vm,第m+1辆标准疲劳车在时间轴上的起始时刻t1,m+1=t1,m+l1/vm+△tm;依次类推,便可根据抽样序列确定k辆标准疲劳车在时间轴上的起始时刻t1,m和结束时刻t2,m(m=1,……,k),以及在时间段[t1, m,t2,m]内的行驶车速vm;
步骤4)对沥青铺装层的温度值进行随机模拟,确定步骤3)各时间段[t1,m,t2,m]内对应的铺装层温度值T:
参照相关专利,申请号:201110195454.5,公开号:102393877A,名称为:一种桥梁结构钢箱梁随机温度场的模拟方法,对时间段[t1,1,t2,k]内的铺装层温度值T进行随机模拟,并从中筛选出步骤3)各时间段[t1,m,t2,m]内对应的铺装层温度值T;
步骤5)将各时间段[t1,m,t2,m]内的行驶车速vm以及铺装层温度值T代入步骤2)的分析模型,计算每个时间段各个时刻对应的疲劳应力值,并将不同时间段同一时刻的疲劳应力值线性叠加,确定车辆荷载与温度共同作用下的焊接细节疲劳应力时程:
(a2)基于步骤3)和步骤4),根据第m个时间段[t1,m,t2,m]中tm时刻对应的行驶车速vm(tm)以及铺装层温度值T(tm),令t=tm-t1,m、v=vm(tm)、T=T(tm),代入步骤2)的疲劳荷载效应分析模型S(tv/60,T),确定第m辆标准疲劳车在时间轴上tm时刻的疲劳应力值,对所有时间段的时刻遍历确定k辆标准疲劳车在时间轴上的疲劳应力时程;
(b2)将k辆标准疲劳车的疲劳应力时程在时间轴上同一时刻的疲劳应力值线性叠加,得到疲劳应力值在时间轴上的模拟序列,从而确定车辆荷载与温度共同作用下的钢桥面板疲劳应力时程;
实施例:
下面以润扬大桥北汊斜拉桥钢桥面板在车辆荷载和温度共同作用下的疲劳应力计算分析为例,说明本发明的具体实施过程:
(1)根据步骤1)建立的润扬大桥北汊斜拉桥的钢桥面板-铺装层一体化分析模型如图1所示,在此基础上建立步骤2)任一工况下钢桥面板焊接细节疲劳应力时程计算的分析模型:
①确定标准疲劳车及加载方式:步骤2)所用的标准疲劳车荷载模型如图2和图3所示,考虑到标准疲劳车同一轴的两侧轮压对疲劳应力的计算结果互不影响,采用单侧轮压加载方式,加载部位和疲劳应力监测部位分别如图4和图5所示,其中图5中测点A为顶板-纵肋焊接细节监测部位,测点B为纵肋对接焊接细节监测部位;
②按照步骤2)中的第(a)步进行加载,得到两类焊接细节的八条疲劳应力时程曲线S1(t)~S8(t),部分曲线分别如图6、7所示;进一步根据步骤2)中的第(b)步建立两类焊接细节的S(t)-T相关模型S(t,T),用二次多项式函数表示为S(t,T)=a(t)T2+b(t)T+c(t),其中a(t)、b(t)、c(t)为待估参数,采用最小二乘法和二次非线性回归方法可确定参数估计值,表1给出了相关模型S(t,T)在时刻t部分取值时的参数估计值,图8进一步给出了时刻t=0.768s时的S(t)-T散点图及估计曲线;
③基于相关模型S(t,T),根据步骤2)中的第(c)步确定任一工况下焊接细节疲劳应力时程计算的分析模型S(vt/60,T),用二次多项式函数表示为S(vt/60,T)=a(vt/60)T2+b(vt/60)T+c(vt/60),图9给出了分析模型S(vt/60,T)在自变量v、T部分取值(v=50km/h、70km/h、90km/h、110km/h,T=10℃)时顶板-纵肋焊接细节的疲劳应力时程;
表1
(2)选取400万辆标准疲劳车在时间轴上进行模拟,时间长度约为一年。基于步骤3)第(a1)步的概率密度统计模型f(s)和p(tp),对行车间距s和时间tp的采集样本进行概率密度统计特性分析,得到f(s)服从η=2.2451、γ=1.3497的对数正态分布,p(tp)服从ρ=0.0436的指数分布;进一步根据步骤3)第(b1)步对确定的概率密度统计模型f(s)和g(v)进行蒙特卡罗抽样,确定行车间距s和行驶车速v在时间轴上的样本序列S和V,S和V的概率密度统计特性如图10和图11所示,可以看出与概率密度统计模型f(s)和p(tp)的估计曲线吻合效果良好;基于样本序列S和V,根据步骤3)第(c1)步方法实现对400万辆标准疲劳车的随机模拟,确定各辆标准疲劳车在时间轴上的起始时刻t1,m和结束时刻t2,m(m=1,……,k),以及在时间段[t1,m,t2,m]内的行驶车速vm;
(3)参照相关专利,申请号:201110195454.5,公开号:102393877A,名称为:一种桥梁结构钢箱梁随机温度场的模拟方法,对一年的沥青铺装层温度值进行随机模拟,确定步骤3)各时间段[t1,m,t2,m]内对应的铺装层温度值T:
①确定一年温度模拟样本的最大值和最小值:以10min为基本时距生成1年的温度模拟样本Tyear,其中模拟样本Tyear的温度值个数N=144×365=52560。若Tyear的最大值和最小值分别用Tyear,maxTyear,min表示,则Tyear,max和Tyear,min的超越概率pmax、pmin均为1/N,利用下式可求得Tyear,max和Tyear,min:
式中,Q(T)为温度T的概率累加分布函数,q(T)为温度T的概率密度函数,其中q(T)服从n个正态分布的加权和,采用下式表示为:
式中,αj为权重,且表示正态分布函数的均值,σj表示正态分布函数的方差,n为≥2的整数,j为整数,且j=1、2、…、n,αj、μj和σj均为待估参数;取n=2并对铺装层温度的采集样本进行概率密度统计特性分析,得到参数估计值如表2所示,进一步根据式(7)、(8)计算得到Tyear,max=51.3℃,Tyear,max=-9.8℃,由此可知温度模拟样本的区间为[-9.8℃,51.3℃]。
表2
待估参数 | α1 | μ1 | σ1 | α2 | μ2 | σ2 |
估计值 | 0.28 | 5.63 | 5.52 | 0.72 | 25.66 | 9.13 |
②在温度模拟样本的区间内进行抽样:利用概率密度函数q(T)在温度模拟样本的区间[-9.8℃,51.3℃]内进行抽样,确定一年内温度模拟样本的温度值。已计算出Tyear的最大值Tyear,max=51.3℃,最小值Tyear,max=-9.8℃,则在区间(-9.8℃,51.3℃)内的剩余抽样个数为N-2。将温度区间(-9.8℃,51.3℃)等间隔划分为100个子区间,则每个子区间的增量ΔT为0.61℃,第i个子区间的范围为[(i-1)×0.61-9.8,i×0.61-9.8],其中i=1,2,...,100。设每个子区间的温度值的个数为Ni,则Ni为:
式中,表示向下取整。由于式(8)向下取整,导致实际生成的温度样本个数少于要求的样本数,其差值ΔN为:
将多余的温度样本数随机分配到100个子区间中,若第i个子区间分配的多余样本数为ΔNi,则第i个子区间最终的温度值的个数为:
根据第i个子区间的温度样本数生成[0,1]之间均匀分布的随机数rand则第i个子区间的温度样本序列为:
对每个子区间进行遍历即可生成满足目标概率密度函数的温度模拟样本序列Tyear。
③对温度模拟样本进行季节变化特征的调整:对第②步得到的温度模拟样本序列Tyear进行重新排列,使其具有实测温度样本的日昼夜变化特征和年季节变化特征。首先根据实测温度样本的年变化特征从Tyear中选出月波峰最大值、月波峰最小值、月波谷最大值、月波谷最小值(Tmax,peak、Tmin,peak、Tmax,trough、Tmin,trough),在此基础上从Tyear中选出合适的2个温度值作为某1天的日波峰值Tpeak和日波谷值Ttrough,使得Tpeak、Ttrough在这1天所在的月份内满足Tmin,peak<Tpeak<Tmax,peak,Tmin,trough<Ttrough<Tmax,trough,对365天遍历便得到1年的日波峰值Tpeak和日波谷值Ttrough,进而根据日波峰值Tpeak和日波谷值Ttrough确定日温度变化的理想正弦曲线。利用非线性最小二乘法从Tyear选取与某时刻理想值最接近的样本,便可得到调整后的温度模拟样本序列Tyear,如图12所示,可以看出调整后的温度模拟样本序列Tyear反映出实测温度的年变化和日变化规律。图13进一步给出了温度模拟样本序列Tyear与实测温度样本的概率密度统计特性对比图,可以发现两者的概率密度曲线(ProbabilityDensityFunction,简称PDF)基本吻合,说明温度模拟样本可以准确反映实测样本的概率密度统计特性;
(4)基于步骤5),确定车辆荷载与温度共同作用下的钢桥面板疲劳应力时程:
①根据第m个时间段[t1,m,t2,m]中tm时刻对应的行驶车速vm(tm)以及铺装层温度值T(tm),令t=tm-t1,m、v=vm(tm)、T=T(tm),代入步骤2)的疲劳荷载效应分析模型S(tv/60,T),确定第m辆标准疲劳车在时间轴上tm时刻的疲劳应力值,对所有时间段的时刻遍历确定400万辆标准疲劳车在时间轴上的疲劳应力时程;
②将400万辆标准疲劳车的疲劳应力时程在时间轴上同一时刻的疲劳应力值线性叠加,得到疲劳应力值在时间轴上的模拟序列,从而确定车辆荷载与温度共同作用下的钢桥面板疲劳应力时程,如图14所示。
Claims (6)
1.一种车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法,其特征在于,包括以下步骤:
1)、建立车辆荷载与温度共同作用下的钢桥面板-铺装层一体化疲劳分析模型;
2)、选取单辆标准疲劳车,加载于步骤1)建立的所述钢桥面板-铺装层一体化疲劳分析模型,计算预定工况下焊接细节的疲劳应力时程曲线,在此基础上建立任一工况下焊接细节疲劳应力时程计算的分析模型,其中,步骤2)中选取单辆标准疲劳车,加载于步骤1)建立的所述钢桥面板-铺装层一体化疲劳分析模型,计算预定工况下焊接细节的疲劳应力时程曲线,在此基础上建立任一工况下焊接细节疲劳应力时程计算的分析模型,具体包括以下步骤:
(a)取单辆标准疲劳车荷载模型,以行驶车速v0沿车道中心线驶过步骤1)建立的一体化疲劳分析模型,根据沥青铺装层温度T的不同分为n种不同的计算工况,计算钢桥面板焊接细节的疲劳应力时程曲线S1(t)~Sn(t),用集合表示为S(t)={S1(t),……,Sn(t)},其中t为时间,Sr(t)为第r种计算工况下的疲劳应力时程曲线,r=1,……,n;
(b)根据焊接细节的n条疲劳应力时程曲线S(t),选取其在t0时刻的疲劳应力值S(t0)及其对应的沥青铺装层温度值T,绘制S(t0)-T散点图并利用最小二乘法进行二次多项式函数拟合,建立S(t0)-T相关模型S(t0,T),据此可得到t0时刻任一铺装层温度T时的疲劳应力值,根据此方法对S(t)的所有时刻遍历可得到任一时刻t、任一铺装层温度T时的疲劳应力值,用函数S(t,T)表示;
(c)对函数S(t,T)的时间变量t乘以比例系数v/v0,用函数S(vt/v0,T)表示,据此可得到自变量v、t、T在给定值时的焊接细节疲劳应力值,其中函数S(tv/v0,T)可作为焊接细节疲劳应力时程计算的分析模型;
3)、选取k辆标准疲劳车进行随机模拟,其中,k≥2,确定各辆标准疲劳车在时间轴上的起始时刻t1,m、结束时刻t2,m以及在时间段[t1,m,t2,m]内的行驶车速vm,其中m=1,……,k,步骤3)中确定各辆标准疲劳车在时间轴上的起始时刻t1,m、结束时刻t2,m以及在时间段[t1,m,t2,m]内的行驶车速vm,具体包括以下步骤:
(a1)、确定标准疲劳车的行车间距s和行驶车速v的概率密度统计模型f(s)和g(v):
标准疲劳车行车间距s的概率密度统计模型f(s)服从对数正态分布,用公式表达为:
式中,η、γ为待估参数,利用Matlab统计工具箱中的对数正态分布函数,对行车间距采集样本进行概率密度统计特性分析,确定η、γ参数估计值,
标准疲劳车行驶车速v的概率密度统计模型g(v)采用下式表达:
g(v)=f(s)/p(tp)
式中,tp为某辆标准疲劳车最后一车轴经过观测点至后一辆标准疲劳车第一车轴经过观测点的所用时间,p(tp)为时间tp的概率密度统计模型,服从指数分布,其中ρ为p(tp)的待估参数,利用Matlab统计工具箱中的指数分布函数,对时间tp的采集样本进行概率密度统计特性分析,确定ρ参数估计值;
(b1)、对概率密度统计模型f(s)和g(v)进行蒙特卡罗抽样,确定行车间距s和行驶车速v在时间轴上的样本序列S和V:
利用隐式蒙特卡罗抽样公式,对标准疲劳车行车间距s的概率密度统计模型f(s)进行均匀随机抽样,得到标准疲劳车行车间距s的样本序列S:
F(S)=F(Sa)+[F(Sb)-F(Sa)]·rand(0,1)
式中,Sa、Sb分别为样本序列S的上下限;F(S)、F(Sa)、F(Sb)分别为f(s)在区间(-∞,S]、(-∞,Sa]、(-∞,Sb]上的积分值;[F(Sb)-F(Sa)]·rand(0,1)为在区间(0,F(Sb)-F(Sa)]上的均匀随机抽样,其中k-1为样本序列S的容量;
利用隐式蒙特卡罗抽样公式,对标准疲劳车行驶车速v的概率密度统计模型g(v)进行均匀随机抽样,得到标准疲劳车行驶车速v的样本序列V:
F(V)=F(Va)+[F(Vb)-F(Va)]·rand(0,1)
式中,Va、Vb分别为样本序列V的上下限;F(V)、F(Va)、F(Vb)分别为f(v)在区间(-∞,V]、(-∞,Va]、(-∞,Vb]上的积分值;[F(Vb)-F(Va)]·rand(0,1)为在区间(0,F(Vb)-F(Va)]上的均匀随机抽样,其中k为样本序列V的容量;
(c1)基于样本序列S和V,对k辆标准疲劳车进行随机模拟:假定单辆标准疲劳车的各轴距之和为l1,一体化分析模型的顺桥向长度为l2,若在t1,m时刻,第m辆标准疲劳车第一车轴以行驶车速vm驶入一体化分析模型,其中行驶车速vm位于样本序列V的第m列,则经历时间l1/vm后,第m辆标准疲劳车的最后一轴会经过观测点,此后再等待时间间隔△tm后,其中△tm=Sm/vm+1,Sm位于样本序列S的第m列,第m+1辆标准疲劳车第一车轴会以行驶车速vm+1驶入一体化分析模型,据此分析可得到第m辆标准疲劳车在时间轴上的结束时刻t2,m=t1,m+l1/vm+l2/vm,第m+1辆标准疲劳车在时间轴上的起始时刻t1,m+1=t1,m+l1/vm+△tm;依次类推,便可根据抽样序列确定k辆标准疲劳车在时间轴上的起始时刻t1,m和结束时刻t2,m(m=1,……,k),以及在时间段[t1,m,t2,m]内的行驶车速vm;
4)、对沥青铺装层的温度值进行随机模拟,确定步骤3)时间段[t1,m,t2,m]内对应的铺装层温度值T;
5)、将时间段[t1,m,t2,m]内的行驶车速vm以及铺装层温度值T代入步骤2)的分析模型,计算每个时间段各个时刻对应的疲劳应力值,并将不同时间段同一时刻的疲劳应力值线性叠加,确定车辆荷载与温度共同作用下的焊接细节疲劳应力时程。
2.如权利要求1所述的车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法,其特征在于,步骤1)所述的钢桥面板-铺装层一体化疲劳分析模型通过以下方法得到:
基于ANSYS有限元分析软件,建立车辆荷载与温度共同作用下的钢桥面板-铺装层一体化疲劳分析模型:
①钢桥面板部分由顶板和纵肋组成,采用壳单元进行模拟;
②沥青铺装层采用实体单元模拟,其中对于沥青铺装层弹性模量的模拟,利用沥青铺装层的弹性模量E与铺装层温度T之间的相关性,采用下式表示:
E=1.2×100.01693(20-T)+3。
3.如权利要求2所述的车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法,其特征在于,钢桥面板的横桥向尺寸取6个U形肋,顺桥向尺寸取5个横隔板间距,四周边界约束模拟为铰接支撑,中间横隔板对钢桥面板的边界约束模拟为刚度大小为1.88×107N/m的竖向弹簧支撑。
4.如权利要求1所述的车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法,其特征在于,步骤(a)中沥青铺装层温度T分为-10℃、0℃、10℃、20℃、30℃、40℃、50℃和60℃八种不同情况。
5.如权利要求1所述的车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法,其特征在于,步骤(a)中所述行驶车速v0=60km/h。
6.如权利要求1所述的车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法,其特征在于,步骤5)中将时间段[t1,m,t2,m]内的行驶车速vm以及铺装层温度值T代入步骤2)的分析模型,计算每个时间段各个时刻对应的疲劳应力值,并将不同时间段同一时刻的疲劳应力值线性叠加,确定车辆荷载与温度共同作用下的焊接细节疲劳应力时程,具体包括以下步骤:
(a2)、基于步骤3)和步骤4),根据第m个时间段[t1,m,t2,m]中tm时刻对应的行驶车速vm(tm)以及铺装层温度值T(tm),令t=tm-t1,m、v=vm(tm)、T=T(tm),代入步骤2)的分析模型,确定第m辆标准疲劳车在时间轴上tm时刻的疲劳应力值,对所有时间段的时刻遍历确定k辆标准疲劳车在时间轴上的疲劳应力时程;
(b2)、将k辆标准疲劳车的疲劳应力时程在时间轴上同一时刻的疲劳应力值线性叠加,得到疲劳应力值在时间轴上的模拟序列,从而确定车辆荷载与温度共同作用下的钢桥面板疲劳应力时程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310121658.3A CN103279588B (zh) | 2013-04-09 | 2013-04-09 | 车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310121658.3A CN103279588B (zh) | 2013-04-09 | 2013-04-09 | 车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103279588A CN103279588A (zh) | 2013-09-04 |
CN103279588B true CN103279588B (zh) | 2016-02-10 |
Family
ID=49062105
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310121658.3A Expired - Fee Related CN103279588B (zh) | 2013-04-09 | 2013-04-09 | 车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103279588B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104268415B (zh) * | 2014-09-30 | 2017-03-08 | 北京交通发展研究中心 | 出行拥堵概率的估计方法 |
CN105005694B (zh) * | 2015-07-13 | 2018-02-13 | 天津大学 | 一种基于动态称重***的桥梁疲劳寿命频域分析方法 |
CN105628373B (zh) * | 2016-03-21 | 2019-05-14 | 金陵科技学院 | 一种焊接钢结构的疲劳损伤计算方法 |
CN105868493B (zh) * | 2016-04-14 | 2018-10-09 | 中铁大桥勘测设计院集团有限公司 | 连续钢桁拱桥盆式橡胶支座的损伤诊断和定位方法 |
CN107220219B (zh) * | 2017-07-10 | 2020-11-20 | 长安大学 | 铁路钢箱梁桥温度梯度模式评价方法 |
CN107609317B (zh) * | 2017-10-25 | 2020-11-17 | 中国铁道科学研究院铁道建筑研究所 | 一种钢桥疲劳评估荷载建模及疲劳评估方法 |
CN108363892A (zh) * | 2018-04-26 | 2018-08-03 | 扬州大学 | 随机车辆荷载作用下应力谱编制方法 |
CN113868749B (zh) * | 2021-10-19 | 2024-03-26 | 大连理工大学 | 一种基于车辆动态称重数据的车致桥梁疲劳损伤分析方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393877A (zh) * | 2011-07-13 | 2012-03-28 | 东南大学 | 一种桥梁结构钢箱梁随机温度场的模拟方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5402274B2 (ja) * | 2009-06-11 | 2014-01-29 | 横浜ゴム株式会社 | 回転体のシミュレーション方法、装置及びプログラム |
-
2013
- 2013-04-09 CN CN201310121658.3A patent/CN103279588B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393877A (zh) * | 2011-07-13 | 2012-03-28 | 东南大学 | 一种桥梁结构钢箱梁随机温度场的模拟方法 |
Non-Patent Citations (2)
Title |
---|
基于长期监测数据的润扬大桥斜拉桥钢箱梁横向温差特性研究;王高新等;《工程力学》;20130131;第30卷(第1期);第163-167,197页 * |
钢桥面铺装在行车荷载与环境温度共同作用下疲劳损伤分析;徐勋倩等;《工程力学》;20101130;第27卷(第11期);第76-81页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103279588A (zh) | 2013-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103279588B (zh) | 车辆荷载与温度共同作用下的钢桥面板疲劳应力计算方法 | |
CN102023317B (zh) | 高速铁路大风监测点布设方法 | |
CN102384856B (zh) | 基于概率有限元法的钢桥疲劳可靠度评估方法 | |
CN103246943B (zh) | 基于马尔可夫链的汽车运行工况多尺度预测方法 | |
CN109815523B (zh) | 基于分解的列车运行多目标差分进化算法 | |
CN104008229A (zh) | 一种街区污染物浓度分布模型建立方法 | |
CN104408924B (zh) | 一种基于耦合隐马尔可夫模型的城市道路异常交通流检测方法 | |
Proppe et al. | A probabilistic approach for assessing the crosswind stability of ground vehicles | |
Solhmirzaei et al. | Road profile estimation using wavelet neural network and 7-DOF vehicle dynamic systems | |
CN105046956A (zh) | 一种基于转向系数的交通流模拟及预测方法 | |
CN106021839B (zh) | 一种地铁隧道电缆支架疲劳可靠性的预测方法 | |
Guo et al. | A novel energy consumption prediction model with combination of road information and driving style of BEVs | |
CN114997027B (zh) | 一种智能求解车桥***随机信号的方法 | |
CN104951607A (zh) | 一种基于流固耦合模拟的光伏支撑***风振计算方法 | |
CN103344395A (zh) | 一种桥梁加固目标承载力的确认方法及装置 | |
CN102141385B (zh) | 沥青路面表面形貌曲面的测试方法 | |
CN109389310B (zh) | 基于蒙特卡洛模拟的电动汽车充电设施成熟度评价方法 | |
CN104166778B (zh) | 汽车座椅导轨接触部位刚度和阻尼优化方法 | |
CN116702096B (zh) | 车辆高原环境道路滑行阻力测算方法及装置 | |
CN108846571A (zh) | 一种网联化电动汽车宏观能量消耗估计方法 | |
CN103547900B (zh) | 在高效室内耐久性试验期间应当施加于轮胎的应力的确定方法 | |
EP2480440B1 (en) | Method for forecasting the evolution of the magnitude of a data for a vehicle journey | |
KR100941968B1 (ko) | 인공윤하중을 이용한 교량의 동적거동 해석방법 | |
Parker et al. | Verification and validation of a winter driving simulator | |
Xu et al. | Neural network predictive control of vehicle suspension |
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: 20160210 Termination date: 20190409 |
|
CF01 | Termination of patent right due to non-payment of annual fee |