Processing math: 100%
  • ISSN 0258-2724
  • CN 51-1277/U
  • EI Compendex
  • Scopus 收录
  • 全国中文核心期刊
  • 中国科技论文统计源期刊
  • 中国科学引文数据库来源期刊

高铁声屏障连接螺栓松弛对疲劳寿命的影响

卫星 汪蓉蓉 温宗意 戴李俊 胡喆

张露, 李冰, 王少华, 李怀仙. 载重车辆-伸缩缝耦合系统的垂向振动数值模拟方法[J]. 西南交通大学学报, 2022, 57(5): 1032-1039. doi: 10.3969/j.issn.0258-2724.20200712
引用本文: 卫星, 汪蓉蓉, 温宗意, 戴李俊, 胡喆. 高铁声屏障连接螺栓松弛对疲劳寿命的影响[J]. 西南交通大学学报, 2023, 58(2): 373-380. doi: 10.3969/j.issn.0258-2724.20210060
ZHANG Lu, LI Bing, WANG Shaohua, LI Huaixian. Numerical Simulation Method for Vertical Vibration of Heavy Vehicle-Expansion Joint Coupled System[J]. Journal of Southwest Jiaotong University, 2022, 57(5): 1032-1039. doi: 10.3969/j.issn.0258-2724.20200712
Citation: WEI Xing, WANG Rongrong, WEN Zongyi, DAI Lijun, HU Zhe. Influence of Bolt Relaxation of High-Speed Railway Sound Barrier on Fatigue Life[J]. Journal of Southwest Jiaotong University, 2023, 58(2): 373-380. doi: 10.3969/j.issn.0258-2724.20210060

高铁声屏障连接螺栓松弛对疲劳寿命的影响

doi: 10.3969/j.issn.0258-2724.20210060
基金项目: 四川省科技创新人才(2020JDRC0009)
详细信息
    作者简介:

    卫星(1976—),男,教授,研究方向为钢结构桥梁行为,E-mail: we_star@swjtu.cn

  • 中图分类号: U213.8

Influence of Bolt Relaxation of High-Speed Railway Sound Barrier on Fatigue Life

  • 摘要:

    螺栓的疲劳寿命和松弛寿命影响着螺栓的使用寿命,在疲劳和松弛的共同作用下,柱脚处连接不断退化,为探究高铁声屏障连接螺栓松弛对疲劳寿命的影响,以某速度为400 km/h高铁声屏障非对称排布和对称排布螺栓为研究对象,利用ANSYS建立柱脚螺栓有限元模型,通过降温法施加预紧力,并施加正负单位弯矩荷载,计算柱脚最不利螺栓在不同预紧力作用下的应力幅,提出了应力幅随预紧力变化的拟合关系式;利用Midas建立声屏障整体模型,分析列车在400 km/h行驶速度下结构动力响应特性,提取柱脚螺栓弯矩时程结果,对仅考虑疲劳失效的螺栓寿命和考虑松弛疲劳共同影响下的螺栓寿命进行比较. 研究结果表明:螺栓松弛会使预紧力下降,导致两种柱脚模型的螺栓应力幅增大;在已有的柱脚螺栓时程计算中,考虑松弛疲劳共同影响下的疲劳寿命比仅考虑疲劳作用时大大降低,当松弛导致预紧力下降至55%以后,将会产生疲劳效应,该结果可为连接结构领域设计人员定量评估螺栓寿命以及对螺栓的维修养护方面提供参考依据.

     

  • 大型桥梁上使用的大位移伸缩缝属于定制产品. 受研发成本的限制,目前大位移伸缩缝的设计和检验均以静态性能指标为依据,对产品动态性能的试验研究应用较少,产品性能无法满足国内桥梁重载和大车流量的需求,使用寿命大幅缩减 [1]. 同时由于缺乏大位移伸缩缝动态特性的研究和高效、精确的计算方法,现有的动态监测数据多用于滑动部件的磨损预测和位移监测[2-3],基于振动特性的故障监测难于实现. 因此,急需对大位移伸缩缝动态特性数值模拟方法展开研究.

    目前,国内外学者的研究多把大位移伸缩缝作为单独的系统进行研究,研究中使用正弦荷载模拟车辆载荷代替车辆的作用,分析了车速、缝宽、弹性元件刚度等参数对伸缩缝性能的影响规律[4-6]. 这类方法只适用于研究伸缩缝的振动规律,由于缺少车辆模型,往往不能准确地反映车辆荷载的冲击效应[7]. 因此,车辆与大位移伸缩缝的耦合作用问题值得考虑.

    Sun等[8]采用轮胎点接触模型建立一种考虑车辆作用的大位移伸缩缝实体单元有限元模型,分析了车速影响下大位移伸缩缝的振动特性,研究结果表明车辆荷载会引起大位移伸缩缝的显著动力响应. Friedl等[9-11]提出了一种基于印迹恒定的平面弹簧轮胎模型,并基于该模型建立了车辆与大位移伸缩缝耦合系统的动力学方程,研究了车速和横梁位置对大位移伸缩缝垂向振动冲击系数的影响,但由于轮胎模型过于简单,需要考虑对轮胎力的修正问题. 丁勇等[12-13]提出了一种基于分布弹簧阻尼单元的车-桥-伸缩缝耦合系统动力学模型建模方法,模型中采用迭代方法求解轮胎的动态印迹长度,对车辆和桥梁的动态响应分析结果表明车辆轴载和桥梁内力的冲击效应高于我国现行桥梁规范的建议值. 吴延平等[14]基于刚性车轮建立了车辆-简支钢箱梁数值模型,验证了伸缩缝的存在会显著增大车辆的冲击效应. 以上的研究结果表明了车辆对伸缩缝有明显的冲击效应. 研究中有限元模型仿真求解方法精度高,但对计算资源的需要过高、计算效率低、难于推广. 数值模型求解速度快,但恒定印迹假设会造成耦合系统中轮胎力的冲击过大,迭代求解动态印迹的方法求解又效率偏低、且采用的轮胎力传递模型没有轮胎的刚性位移对印迹长度的影响,因此,需要对车辆-大位移伸缩缝耦合系统寻求更加细致化合理的方法.

    本文以ZL1600模数式大位移伸缩缝为研究对象,首先在分析轮胎载重车辆-伸缩缝耦合系统受力特点的基础上提出了一种考虑轮胎力真实激励特性的耦合动力学数值模型,并采用新型积分算法对数值模型进行求解,然后进行了轮胎载重车辆-伸缩缝耦合系统动态试验测试,通过仿真结果与试验测试结果对比验证动力学模型的有效性,最后分析了车辆对伸缩缝的冲击效应,为耦合系统性能的进一步研究提供基础.

    ZL1600模数式大位移伸缩缝主要由边梁、中梁、横梁、中梁支承组件、横梁支承组件、位移箱锚固系统和防水密封条等组成. 车辆选用国产某型四轴轮胎载重车辆为计算实例. 车辆通过车轮与伸缩缝中梁顶面接触,中梁通过中梁弹性元件将垂向荷载传递给横梁,横梁通过横梁弹性元件将垂向荷载传递给桥梁梁体,耦合系统受力如图1所示. 图中:F011F021分别为梁体对第1号横梁1端和2端的支反力;FHjij号横梁对i号中梁的支反力;FIiji号中梁对j号横梁的支反力;FTki为第k轴车轮作用于i号中梁的垂向载荷;kS1cS1分别为1轴悬架刚度和阻尼;kS2cS2分别为后悬架刚度和阻尼;mC为车体质量;l1l2l3分别为1轴、2轴和后悬架中心至车体质心的距离;yHyZyTyC分别为横梁、中梁、轮胎轮心和车体质心的垂向位移;zC为车体质心的水平向位移;v为速度.

    图  1  车辆-伸缩缝耦合系统受力分析
    Figure  1.  Force analysis of the vehicle-expansion joint coupled system

    四轴轮胎载重车辆采用十一自由度整车模型[15],其中车体具有3个自由度(垂向振动、侧滚运动和俯仰运动). 每个平衡悬架具有2个自由度(垂向振动和俯仰运动). 前轴左、右车轮分别具有一个垂向振动自由度. 采用D’Alembert原理建立车辆振动微分方程为

    MVAV+CVVV+KVXV=RV (1)

    式中:MVCVKV分别为四轴车辆的质量矩阵、阻尼矩阵和刚度矩阵[15]AVVVXV分别为四轴车辆的加速度矢量、速度矢量和位移矢量;RV为车辆垂向轮胎力矢量.

    大位移伸缩缝的中梁、横梁、中梁弹性元件、横梁弹性元件和剪切弹簧为重要的传力部件. 由于伸缩缝的剪切弹簧对伸缩缝的垂向振动特性影响很小[4],因此以下的研究内容均忽略剪切弹簧. 大位移伸缩缝主要的结构参数如表1所示.

    表  1  大位移伸缩缝主要结构参数
    Table  1.  Main structural parameters of large displacement expansion joint
    项目数值项目数值
    弹性模量 E/GPa205泊松比 ν0.3
    材料密度 ρ/(kg·m−37800缝宽 B/mm0~80
    中梁单位长度质量 mZ/(kg·m−149.33横梁单位长度质量 mH/(kg·m−1315.9
    中梁长度 LZ/m10.2横梁长度 LH/m4.2
    中梁截面惯性矩 IZ/m41.24 × 10−5横梁截面惯性矩 IH/m42.46 × 10−4
    中梁、横梁弹性元件刚度/(kN·mm−180中梁、横梁弹性元件阻尼/(N·s·mm−15
    下载: 导出CSV 
    | 显示表格

    i号中梁(i = 1 ~ NINI为中梁的根数)的振动位移变量为$y_{{\text{Z}}i}^{}(x,t)$,则其振动微分方程为

    EZIZ4yZi(x,t)x4+mZ2yZi(x,t)t2=NJj=1FHji(t)δ(xxj)NKk=1[FTLki(t)δ(xxLk)+FTRki(t)δ(xxRk)], (2)

    式中:EZ为中梁的弹性模量,下标Z、H分别表示中梁和横梁;t为时间;δ(·)为Dirac函数;FTLkiFTRki为第k轴左、右侧车轮作用于i号中梁的垂向载荷;NJ为横梁的根数;NK为车轴的根数;xjj号横梁中心轴的横向坐标;xLk、xRk为第k根车轴左、右车轮的横向坐标.

    中梁与横梁通过弹性元件连接,j号横梁对i号中梁的支反力为

    FHji(t)δ(xxj)=k1[yZi(xi,t)yHj(zi,t)]+c1[˙yZi(xi,t)˙yHj(zi,t)], (3)

    式中:k1为中梁弹性元件的刚度;c1为中梁弹性元件的阻尼;$y_{{\text{H}}j}^{}({\textit{z}}_i^{},t)$j号横梁在坐标zi处的振动位移;$y_{{\text{Z}}i}^{}(x_i^{},t)$i号中梁在坐标xi处的振动位移;zii号中梁中心轴的纵向坐标.

    同理,可得j号横梁(j=1~NJ)振动位移量为$y_{{\text{H}}j}^{}({\textit{z}},t)$,振动微分方程为

    EHIH4yHj(z,t)z4+mH2yHj(z,t)t2=NIi=1FIij(t)δ(zzi)2h=1F0hj(t)δ(zz0h), (4)

    式中:EH为横梁的弹性模量;h为横梁端部弹性元件编号;z0hh端横梁弹性元件的纵向坐标.

    i号中梁对第j号横梁的支反力为

    FIij(t)δ(zzi)=k1[yZi(xj,t)yHj(zi,t)]c1[˙yZi(xj,t)˙yHj(zi,t)]. (5)

    横梁与桥梁梁体通过横梁弹性元件连接,梁体对第j号横梁h端的支反力为

    F0hj(t)δ(zz0h)=k2yHj(z0h,t)+c2˙yHj(z0h,t), (6)

    式中:k2为横梁弹性元件的刚度;c2为横梁弹性元件的阻尼.

    采用Ritz法[16],引入中梁和横梁正则振型坐标qZie(t)和qHjf(t),式(2)和式(4)的解为

    yZi(x,t)=NEe=1YZie(x)qZie(t) (7)
    yHj(z,t)=NFf=1YHjf(z)qHjf(t) (8)

    式中:YZie(x)为i号中梁的e阶振型函数,${Y_{{\text{Z}}ie}}(x) = \sin \left(\dfrac{{e{\text{π }}}}{{L_{\text{Z}}^{}}}x\right)$e =1~NEYHjf(z)为第j号横梁的f阶振型函数,${Y_{{\text{H}}jf}}({\textit{z}}) = \sin \left(\dfrac{{f{\text{π }}}}{{{L_{\text{H}}}}}{\textit{z}}\right)$f为横梁计算阶数,f=1~NF.

    根据模态正交性和Dirac函数的性质,式(2)两边同时乘以YZie(x)后沿中梁长度L1取积分得

    mZLZ2¨qZie(t)+NJj=1c1YZie(xj)NEe=1YZie(xj)˙qZie(t)NJj=1c1YZie(xj)NFf=1YHjf(zi)˙qHjf(t)+EZIZLZ2(eπLZ)4qZie(t)+NJj=1k1YZie(xj)NEe=1YZie(xj)qZie(t)NJj=1k1YZie(xj)NFf=1YHjf(zi)qHjf(t)=NKk=1[YZie(xLk)FTLki+YZie(xRk)FTRki]. (9)

    同理,式(4)两边同时乘以$Y_{\mathrm{H} j f}$(z)后,沿横梁长度LH取积分得

    mHLH2¨qHjf(t)+NIi=1c1YHjf(zi)NFf=1YHjf(zi)˙qHjf(t)NIi=1c1YHjf(zi)NEe=1YZie(xj)˙qZie(t)+2h=1c2YHjf(z0h)NFf=1YHjf(z0h)˙qHjf(t)+EHIHLH2(fπLH)4qHjf(t)+NIi=1k1YHjf(zi)NFf=1YHjf(zi)qHjf(t)NIi=1k1YHjf(zi)NEe=1YZie(xj)qZie(t)+2h=1k2YHjf(z0h)NFf=1YHjf(z0h)qHjf(t)=0. (10)

    大位移伸缩缝的中梁垂直于路面中轴线布置,因此,轮胎与中梁的接触面具有非连续性的特点,需考虑轮胎印迹范围内局部脱空的情况. 为了准确地表达轮胎与中梁的接触关系,轮胎模型在弹性滚子模型[17-18]中加入了伸缩缝中梁顶面轮廓函数,轮胎与中梁接触关系如图2所示,图中:$y_{0}\left(Z_{{\rm{T}}k} + {Z_{{\rm{A}}k}}\right)$为伸缩缝中梁顶面轮廓函数,ZTk为第k号轮胎轮心在总体坐标系中的纵向坐标,ZTk=vt,${Z_{{\rm{A}}k}}$为第k号轮胎接触印内的局部纵向坐标,−r0${Z_{{\rm{A}}k}}$r0r0为轮胎自由半径.

    图  2  轮胎与中梁接触关系
    Figure  2.  The contact relationship between the center beam and tire

    弹性滚子模型考虑了轮胎的包容性[17](即几何滤波和弹性滤波特性)、轮胎轴心振动和路面振动条件,能较好地模拟轮胎过缝的过程. 模型中第k号轮胎与伸缩缝第i号中梁各接触点的压缩量为

    Δy(ZTk,ZAk)={yZi(xk,t)+y0(ZTk+ZAk)+r20Z2Akr0yTk(ZTk)yTks,Δy>0,0,Δy0, (11)

    式中:yTks为车辆静止时第k号车轮的下压量,$ y_{{\text{T}}ks}^{} = - G_{{\text{T}}k}^{}/k_{{\text{T}}ks}^{} $GTk为第k号轮胎分配的车辆重力,kTks为第k号轮胎静垂直刚度.

    根据伸缩缝中梁顶面的非连续性特点,$y_{0}\left(Z_{{\rm{T}}k} + {Z_{{\rm{A}}k}}\right)$表示为

    y0(ZTk+ZAk+t)={0,ziBZ2ZTk+ZAkzi+BZ2,HZ,ziBZ2BF<ZTk+ZAk<ziBZ2zi+BZ2<ZTk+ZAk<zi+BZ2+BF, (12)

    式中:BZHZ分别为中梁顶面宽度和中梁高度;BF为缝宽.

    k号轮胎与伸缩缝中梁之间的垂向轮胎合力为

    FTk(t)=NIi=1FTki(t) (13)

    式中:FTki(t)为k号轮胎作用于i号中梁的轮胎力,如式(14).

    FTkit=aki[kTkΔy(ZTk,ZAk,t)+cTkΔ˙y(ZTk,ZAk,t)]dZAk, (14)

    式中:kTk为轮胎垂直分布刚度;cTk为轮胎垂直分布阻尼;aki为第k号轮胎与伸缩缝第i号中梁的接触印迹长度.

    由于轮胎的分布垂直刚度、阻尼测定较为困难,因此计算时采用简化公式[17],轮胎的分布垂直刚度kT、阻尼cT

    {kT=kTs/2ak,cT=cTs/2ak, (15)

    式中:kTscTs分别为轮胎静垂直刚度、阻尼;ak为轮胎与接触面的接触印迹半长度,如式(16).

    ak=r20f2a (16)

    式中:fa为轮胎在荷载作用下的胎面位移.

    子午线轮胎在荷载作用下变形时存在带束层刚性位移现象[19-20]. 带束层刚性位移即滚动轴心偏离原轮胎中心的距离为fb,则轮胎在荷载作用下的下沉量ΔyA

    ΔyA=fb+fa=r0r (17)

    式中:r为轮胎滚动半径.

    联立式(1)、(9)、(10)、(11)、(12)得到车辆-伸缩缝耦合系统的运动微分方程组. 采用求解效率高、稳定性好的新型快速积分法[16],在MATLAB软件中编制计算程序求解方程组.

    为了得到测试结果来验证耦合系统数值模型仿真结果的有效性,对ZL1600大位移伸缩缝展开了试验测试.

    在成都某公司试验场完成试验. 大位移伸缩缝缝宽B为40 mm. 车辆选用国产某型四轴轮胎载重车辆,车辆参数见表2所示. 伸缩缝由6根横梁(H1~H6)和19根中梁(Z1~Z19)组成. 在实际应用中,载重车辆一般行走于外侧货车道,因此试验中车辆右前轮从伸缩缝第5根横梁上方位置驶过. 在Z1、Z10、Z13、Z19号中梁位于H5、H6号横梁和H4、H5号横梁跨中处分别布置A、B两组速度传感器,如图3所示. 传感器频率范围10~1000 Hz,灵敏度30 mV/(mm·s−1).

    表  2  四轴载重车辆动力学参数
    Table  2.  Dynamic parameters of four-axle vehicle
    参数数值参数数值
    车体质量 mC/kg64501 轴至车体质心
    的距离 l1/m
    3.327
    车体俯仰惯量JCx/(kg·m2883272 轴至车体质心
    的距离 l2/m
    1.227
    车体侧滚惯量JCz/(kg·m217665后悬架中心至车体质心的距离 l3/m3.948
    后悬架平衡
    杆质量 mS/kg
    2003、4 轴之间的
    距离 l4/m
    1.35
    后悬架平衡杆俯仰惯量 JS/(kg·m23801、2 轴的轮距
    一半 b1/m
    1.0135
    1、2 轴非簧载
    的质量 mT1 /kg
    5003、4 轴的轮距
    一半 b2/m
    0.93
    3、4 轴非簧载
    的质量 mT2/kg
    8001、2 悬架阻尼
    cS1/(kN·s·m−1
    25.32
    1、2 轴悬架刚
    kS1 /(kN·m−1
    284后悬架阻尼cS2/(kN·s·m−150.636
    后悬架刚度kS2/(kN·m−120641、2 轴轮胎阻尼
    cT1/(kN·s·m−1
    3.5
    1、2 轴轮胎刚
    kT1 /(kN·m−1
    14023、4 轴轮胎阻尼 cT2/(kN·s·m−17
    3、4 轴轮胎刚
    kT2 /(kN·m−1
    2804轮胎型号11.00R20-18RP
    下载: 导出CSV 
    | 显示表格

    图4为车辆以30 km/h时车速驶过伸缩缝时Z1号中梁A组测点的速度时程曲线对比结果. 由图可知,仿真曲线与试验测试曲线振动趋势基本一致. 由于试验测试时实际车速低于试验设计车速30 km/h,因此,试验测试曲线中轮胎作用在1号中梁的时间较仿真曲线有滞后.

    图  3  传感器布置
    Figure  3.  Layout of the sensors
    图  4  1号中梁A组测点的速度时程曲线(v = 30 km/h)
    Figure  4.  Speed time history curves of group A measuring points of No.1 middle beam(v = 30 km/h)

    图5为车辆以不同车速通过伸缩缝时Z1号中梁A组测点振动速度最大值对比结果. 由图可知:同一车速下1轴和4轴分别通过Z1号中梁时测点的振动速度最大值仿真结果与试验测试结果基本一致,测点的振动速度最大值随速度的变化趋势也一致.

    图  5  1号中梁A组测点振动速度最大值对比
    Figure  5.  Comparison results of maximum speed at group A measuring points of No.1 middle beam

    对中梁测点的速度时程曲线进行积分[21],可以得到测点的位移时程曲线. 试验测试和仿真计算分别得到车辆以不同车速通过伸缩缝时Z1、Z10、Z13、Z19号中梁测点最大下沉量对比结果见表3. 由表可知:同一车速下各中梁A组、B组测点的最大下沉量仿真结果与试验测试结果最大偏差分别为7.8%、8.7%.

    表  3  中梁测点最大下沉量对比结果
    Table  3.  Comparison results of maximum sinking displacement of center beam test points
    传感器
    位置
    车速/
    (km·h−1
    1 号中梁10 号中梁13 号中梁19 号中梁
    试验
    值/mm
    仿真
    值/mm
    偏差/
    %
    试验
    值/mm
    仿真
    值/mm
    偏差/
    %
    试验
    值/mm
    仿真
    值/mm
    偏差/
    %
    试验
    值/mm
    仿真
    值/mm
    偏差/
    %
    A 组 30 0.42 0.41 2.4 0.49 0.51 4.1 0.54 0.56 3.7 0.55 0.58 5.5
    40 0.50 0.48 4.0 0.56 0.60 7.1 0.58 0.61 5.2 0.64 0.69 7.8
    60 0.54 0.51 5.6 0.58 0.59 1.7 0.62 0.65 4.8 0.60 0.64 6.7
    B 组 30 0.14 0.13 7.1 0.23 0.24 4.3 0.20 0.21 5.0 0.22 0.21 4.5
    40 0.18 0.17 5.6 0.25 0.27 8.0 0.29 0.31 6.9 0.26 0.28 7.7
    60 0.20 0.19 5.0 0.29 0.31 6.9 0.31 0.30 3.2 0.23 0.25 8.7
    下载: 导出CSV 
    | 显示表格

    用冲击系数μ来表示车辆轮胎力对伸缩缝冲击效应的大小,如式(18).

    μ=(FTmax/FTs)/FTs, (18)

    式中:FTmax为轮胎动态力最大值;FTs为静轮胎力.

    图6为车辆以80 km/h驶过伸缩缝各中梁(Z1~Z19)及上桥(B2)时车辆轮胎力冲击系数. 由图可知:车辆驶过伸缩缝各中梁时,1、2轴轮胎力冲击系数小于3、4轴;车轮驶过伸缩缝上桥时轮胎力冲击系数最大,此处边梁与桥梁梁体端部连接处混凝土最容易发生破坏,与实际应用中梁体端部混凝土病害高发的情况一致.

    图  6  车辆驶过伸缩缝时轮胎力冲击系数(v = 80 km/h)
    Figure  6.  Impact factors of tire force when vehicle passing through expansion joint (v = 80 km/h)

    图7为车辆驶上伸缩缝时轮胎力最大冲击系数随车速变化的关系. 由图可知:车辆驶上伸缩缝时各轴轮胎力最大冲击系数随着车速的增大而增大,且3、4轴轮胎力最大冲击系数增速高于1、2轴;当车速为120 km/h时,1、2轴轮胎力的冲击系数最大值均为0.31,小于我国现行规范[21]的推荐值0.45,3、4轴轮胎力的冲击系数最大值分别为0.62、0.67,均大于规范[21]的推荐值0.45.

    图  7  车辆驶上伸缩缝时的最大冲击系数与车速的关系
    Figure  7.  Relationships between maximum impact factor and speed when vehicle getting on the expansion joint

    图8为车辆驶过伸缩缝上桥时轮胎力最大冲击系数随车速变化的关系. 由图可知:车辆上桥时各轴轮胎力最大冲击系数随着车速的增大而增大;当车速为120 km/h时,轮胎力的冲击系数最大值均为0.82,大于规范[21]的推荐值0.45.

    图  8  车辆驶上桥梁时的最大冲击系数与车速的关系
    Figure  8.  Relationship between maximum impact factor and speed when vehicle getting on the bridge

    本文提出了一种轮胎载重车辆-大位移伸缩缝耦合系统动力学模型,并通过试验测试结果对动力学模型的有效性进行了验证,在此基础上研究了车辆对大位移伸缩缝的冲击效应,研究结果表明:

    1) 车辆-大位移伸缩缝耦合系统动力学模型能有效地仿真耦合系统动态响应,仿真得到的中梁振动速度时程曲线与试验测试值趋势相同,中梁测点的振动速度最大值仿真结果与试验测试值随速度的变化趋势表现出较好的一致性,中梁测点最大下沉量与试验测试值的最大偏差小于10.0%,可满足工程精度.

    2) 车辆后双轴轮胎力的冲击系数大于前轴,且车辆对大位移伸缩缝末端边梁与桥梁连接处混凝土的冲击大于对伸缩缝中梁的冲击,实际应用中应考虑对此处桥面进行加强.

    3) 车辆对大位移伸缩缝中梁和伸缩缝末端桥梁的轮胎力冲击系数均随车速的增大为增大,当车速为120 km/h时,对两者的冲击系数最大值分别为0.67和0.82,大于我国现行规范中的局部冲击系数取0.45的建议值. 建议在伸缩缝和桥梁端部局部结构的设计中引入车辆冲击系数时应在现有规范推荐值的基础上进行适当的放大. 同时建议进一步优化大位移伸缩缝中梁与车道线的夹角以降低车辆的冲击效应.

    本文所建立的轮胎载重车辆-大位移伸缩缝耦合系统动力学模型能可靠地求解耦合系统动态响应,为大位移伸缩缝的设计和动态特性研究提供了一种计算方法,后续可以针对各参数和病害状态对大位移伸缩缝的性能影响进行进一步的研究.

  • 图 1  螺栓柱脚细部构造

    Figure 1.  Detailed structure of bolt base

    图 2  建立接触单元

    Figure 2.  Establish contact unit

    图 3  螺栓柱脚整体模型

    Figure 3.  Whole model of bolt base

    图 4  荷载施加模型

    Figure 4.  Model with applied load

    图 5  螺栓应力变化曲线

    Figure 5.  Variation of bolt stress with preload

    图 6  螺栓应力幅变化曲线

    Figure 6.  Variation of bolt stress amplitude with preload

    图 7  螺栓应力变化曲线

    Figure 7.  Variation of bolt stress with preload

    图 8  螺栓应力幅变化曲线

    Figure 8.  Variation of bolt stress amplitude with preload

    图 9  连接部分

    Figure 9.  Connection parts

    图 10  桩基础边界条件

    Figure 10.  Boundary conditions of pile foundation

    图 11  声屏障整体模型

    Figure 11.  Whole model of sound barrier

    图 12  脉动风压下立柱根部弯矩

    Figure 12.  Bending moment of column root under fluctuating wind pressure

    图 13  底板螺栓应力时程变化曲线

    Figure 13.  Time history curve of stress of base plate bolt

    表  1  材料参数

    Table  1.   Material parameters

    材料弹性模量/Pa泊松比切线模量/GPa
    Q2352.06 × 10110.32.06
    C303.00 × 10100.2
    下载: 导出CSV

    表  2  预紧力大小与对应应力幅

    Table  2.   Preload value and the corresponding stress amplitude

    预紧力
    大小/kN
    单位弯矩/
    (MN·mm)
    预紧力/
    MPa
    应力幅/
    MPa
    0 1 5.02 5.20
    −1 −0.26
    16.0 1 2.18 2.32
    −1 −0.20
    48.0 1 0.94 1.06
    −1 −0.17
    80.0 1 0.65 0.76
    −1 −0.15
    160.0 1 0.45 0.55
    −1 −0.14
    下载: 导出CSV

    表  3  预紧力大小与对应应力幅

    Table  3.   Preload value and the corresponding stress amplitude

    预紧力
    大小/kN
    单位弯矩/
    (MN·mm)
    应力/
    MPa
    应力幅/
    MPa
    0 1 6.63 6.85
    −1 −0.31
    7.5 1 4.92 5.05
    −1 −0.18
    22.5 1 2.27 2.38
    −1 −0.15
    37.5 1 1.42 1.50
    −1 −0.12
    75.0 1 0.80 0.83
    −1 −0.04
    下载: 导出CSV

    表  4  材料参数

    Table  4.   Material parameters

    部件材料密度/
    (kg·m−3
    弹性模量/
    GPa
    泊松比
    钢构件Q23578502060.30
    底梁和桩基础C302360300.20
    铝合金复合
    吸声板
    复合
    材料
    250710.33
    下载: 导出CSV

    表  5  桩基主要控制参数

    Table  5.   Main control parameters of pile foundation

    名称控制值
    桩侧水平抗力系数的比例系数/(MN·m−4≥60
    桩侧土土壤容重/(kN·m−317.5~19.5
    桩侧土土壤内摩擦角/(°)≥30
    桩侧限端阻力标准值/kPa≥850
    桩侧土的阻力/kPa≥70
    下载: 导出CSV

    表  6  雨流计数结果

    Table  6.   Rain flow counting results

    应力幅值/MPa作用次数
    12.19 1
    4.69
    4.14
    3.39
    3.08
    下载: 导出CSV

    表  7  雨流计数结果

    Table  7.   Rain flow counting results

    应力幅值/MPa作用次数
    35.111
    13.50
    11.93
    9.75
    8.86
    下载: 导出CSV

    表  8  400辆/d疲劳损伤分析

    Table  8.   Fatigue damage analysis for 400 veh/d

    循环
    次数/次
    应力幅值/
    MPa
    n/次疲劳
    损伤
    T/d
    40035.116.41×1066.24×10−516037.23
    下载: 导出CSV
  • [1] 孔繁晓,言婷,周海波. 预紧力对风电叶片根部螺栓疲劳寿命的影响分析[J]. 风机技术,2017,59(6): 49-52. doi: 10.16492/j.fjjs.2017.06.0008

    KONG Fanxiao, YAN Ting, ZHOU Haibo. Influence of the pretightening stress on the fatigue life of a bolt used for wind turbine blades[J]. Chinese Journal of Turbomachinery, 2017, 59(6): 49-52. doi: 10.16492/j.fjjs.2017.06.0008
    [2] 朱若燕,李厚民. 高强度螺栓的预紧力及疲劳寿命[J]. 湖北工学院学报,2004,19(3): 135-136,141.

    ZHU Ruoyan, LI Houmin. Preload and fatigue life of high strength bolt[J]. Journal of Hubei University of Technology, 2004, 19(3): 135-136,141.
    [3] 吴勇,陈琴珠,邹慧君. 摩擦离合器螺栓联接预紧力对疲劳寿命的影响[J]. 机械设计与研究,2012,28(6): 67-69.

    WU Yong, CHEN Qinzhu, ZOU Huijun. Effect of preload on fatigue life of friction clutch bolt[J]. Machine Design & Research, 2012, 28(6): 67-69.
    [4] 刘嘉慧,林腾蛟,吕和生,等. 多工况下风电齿轮箱联接螺栓疲劳寿命分析[J]. 机械研究与应用,2020,33(3): 95-101. doi: 10.16576/j.cnki.1007-4414.2020.03.027

    LIU Jiahui, LIN Tengjiao, LV Hesheng, et al. Fatigue life analysis of the bolt connections for wind turbine gearbox under multiple working conditions[J]. Mechanical Research & Application, 2020, 33(3): 95-101. doi: 10.16576/j.cnki.1007-4414.2020.03.027
    [5] 祝雨. 螺栓联接的预紧力与疲劳强度的讨论[J]. 化工管理,2016(35): 256.
    [6] 郭卫凡,唐文良. 螺栓联接的预紧力与疲劳强度的讨论[J]. 科技视界,2013(23): 65-66.

    GUO Weifan, TANG Wenliang. A discussion of effect of preload on bolt joint fatigue strength[J]. Science & Technology Vision, 2013(23): 65-66.
    [7] LIU Y Z, CHEN J, ZHANG X F, et al. Fatigue behaviour of blind bolts under tensile cyclic loads[J]. Journal of Constructional Steel Research, 2018, 148: 16-27. doi: 10.1016/j.jcsr.2018.05.019
    [8] PENNEC F, DURIF S, CAMARA A B, et al. 01.16: Fatigue behaviour analysis of bolts in tee-stub steel connections[J]. ce/papers, 2017, 1(2/3): 298-307.
    [9] 曹罚君. 多轴非比例载荷谱条件下高强度螺栓疲劳强度分析[J]. 机械强度,2019,41(1): 232-237.

    CAO Fajun. Fatigue strength analysis of high-strength joint bolts under multiaxial nonproportional loading spectrum[J]. Journal of Mechanical Strength, 2019, 41(1): 232-237.
    [10] 刘育,晋健,王勇飞,等. 基于PD-STFA的水轮机顶盖联接螺栓疲劳寿命计算方法研究[J]. 机械强度,2020,42(6): 1459-1465.

    LIU Yu, JIN Jian, WANG Yongfei, et al. Research on the fatigue life calculation method of hydro-turbine head cover connecting bolts based on pd-stfa[J]. Journal of Mechanical Strength, 2020, 42(6): 1459-1465.
    [11] 汤春球,周科帆,莫易敏,等. 基于有限元仿真的螺栓疲劳性能分析[J]. 数字制造科学,2017,15(增1): 30-35.

    TANG Chunqiu, ZHOU Kefan, MO Yimin, et al. FEA fatigue simulation on the bolt performance[J]. Digital Manufacture Science, 2017, 15(S1): 30-35.
    [12] 彭飞. 螺栓预紧力对发动机气缸盖疲劳特性的影响[J]. 唐山学院学报,2020,33(3): 56-59.

    PENG Fei. Effect of the bolt pretightening force on the fatigue properties of engine cylinder head[J]. Journal of Tangshan University, 2020, 33(3): 56-59.
    [13] 杜静,黄文,王磊,等. 基于接触分析的高强度螺栓疲劳寿命分析[J]. 现代科学仪器,2013(1): 73-77,90.

    DU Jing, HUANG Wen, WANG Lei, et al. Fatigue life analysis of the high intention bolt based on contact analysis[J]. Modern Scientific Instruments, 2013(1): 73-77,90.
    [14] MAJZOOBI G H, FARRAHI G H, HABIBI N. Experimental evaluation of the effect of thread pitch on fatigue life of bolts[J]. International Journal of Fatigue, 2005, 27(2): 189-196. doi: 10.1016/j.ijfatigue.2004.06.011
    [15] RAHMAN N A, TIZANI W. Fatigue performance of blind bolt in concrete-filled hollow section[C]// Proceedings of 2013 4th International Conference on Information Technology for Manufacturing Systems (ITMS 2013 Ⅳ). Auckland: Trans Tech Publications Ltd., 2013: 5-10.
    [16] WANG H L, YIN H J, LIU K. Fatigue performance analysis of frictional type high strength bolts of overlapped joints[C]//13th International Conference on Fracture (ICF13). Beijing: Curran Associates, Inc., 2013: 2132-2137.
    [17] 张新鹏,张广泰,张辉亮,等. 螺旋流道水冷IGBT散热器数值模拟及试验研究[J]. 电力电子技术,2014,48(2): 71-73. doi: 10.3969/j.issn.1000-100X.2014.02.023

    ZHANG Xinpeng, ZHANG Guangtai, ZHANG Huiliang, et al. Numerical simulation and experimental research on spiral flow channel water-cooling IGBT radiator[J]. Power Electronics, 2014, 48(2): 71-73. doi: 10.3969/j.issn.1000-100X.2014.02.023
    [18] 中华人民共和国国家标准. 六角头螺栓: GB/T 5782—2016[S]. 北京: 中国标准出版社, 2016.
    [19] 中华人民共和国国家标准. 1型六角螺母: GB//T 6170—2015[S]. 北京: 中国标准出版社, 2016.
    [20] 中华人民共和国国家标准. 钢结构设计标准: GB 50017—2017[S]. 北京: 中国建筑工业出版社, 2017.
    [21] 赵丽滨,龙丽平,蔡庆云. 列车风致脉动力下声屏障的动力学性能[J]. 北京航空航天大学学报,2009,35(4): 505-508. doi: 10.13700/j.bh.1001-5965.2009.04.012

    ZHAO Libin, LONG Liping, CAI Qingyun. Dynamic properties of noise barrier structure subjected to train-induced impulsive wind pressure[J]. Journal of Beijing University of Aeronautics and Astronautics, 2009, 35(4): 505-508. doi: 10.13700/j.bh.1001-5965.2009.04.012
    [22] 屈爱平,高淑英. 梁-墩-桩基的动力特性研究[J]. 西南交通大学学报,2001,36(6): 641-644. doi: 10.3969/j.issn.0258-2724.2001.06.021

    QU Aiping, GAO Shuying. Vibration characteristics of a 3-dimensional beam-pier-pile system[J]. Journal of Southwest Jiaotong University, 2001, 36(6): 641-644. doi: 10.3969/j.issn.0258-2724.2001.06.021
  • 加载中
图(13) / 表(8)
计量
  • 文章访问数:  430
  • HTML全文浏览量:  188
  • PDF下载量:  50
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-01-21
  • 修回日期:  2021-07-24
  • 网络出版日期:  2022-11-05
  • 刊出日期:  2021-08-05

目录

/

返回文章
返回