没有合适的资源?快使用搜索试试~ 我知道了~
列车荷载下冻土路基非平稳随机振动分析的时域显式法.docx
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 42 浏览量
2023-02-23
20:19:44
上传
评论
收藏 368KB DOCX 举报
温馨提示
试读
23页
列车荷载下冻土路基非平稳随机振动分析的时域显式法.docx
资源推荐
资源详情
资源评论
0 引言
我国为世界第三冻土大国,多年冻土面积占国土总面积的 21.5%,其中高
原冻土面积更是居于世界之最
[ 1]
。随着“一带一路”倡议深入推进实施,冻土区
的铁路建设与维护得到高度重视。同时,国民经济的发展使得列车行驶频次有
所提升,加之受全球变暖
[ 2]
等气候因素的影响,冻土区铁路病害问题越来越
突出。为保障冻土区铁路的安全运营,亟须发展高效准确的计算方法,分析列
车行驶对冻土路基动力响应的影响。
冻土路基在列车荷载作用下的动力响应是当前的研究热点之一。马巍等
[ 3 ]
和朱占元等
[ 4 ]
开展了列车荷载作用下冻土路基动力响应的现场监测研究,
为列车通过冻土路段时路基稳定性的预测分析提供了实测资料依据。Stevens
[ 5]
、孔祥兵等
[ 6]
和李双洋等
[ 7-9]
采用数值分析方法,研究了一次行车荷载作
用下路基的位移和应力等动力响应问题。此外,还有部分学者
[ 10-13 ]
进行了列
车荷载作用下冻土路基的稳定性、蠕变、轨枕响应和车速对路基响应影响等方
面的数值模拟研究。
上述研究中,都采用了确定性的列车荷载,选取某一次列车驶过产生的荷
载样本计算路基响应,其结果难以反映铁路实际运行中由轨道不平顺及列车类
型和轴重差异导致的荷载随机性
[ 14]
。对于路基的特定横断面,列车驶入和驶
出的过程导致其顶部荷载的统计特征具有显著的时变特性。因此,在研究路基
动力响应时应将其顶部列车荷载按非平稳随机荷载考虑
[ 7, 15]
。由于冻土独特
的物理、热力学性质以及非平稳随机振动分析的复杂性,目前关于随机列车荷
载作用下冻土路基动力响应的研究较少。
目前,解决结构非平稳随机振动问题的主要方法有蒙特卡罗模拟法、虚拟
激励法
[ 16]
、概率密度演化法
[ 17]
以及时域显式法
[ 18]
等。其中,时域显式法通
过构建系统动力响应的时域显式表达式,在随机分析中可以针对任意关键响应
进行降维,具有理想的计算效率,已成功应用于高层建筑
[ 19 ]
、大跨度桥梁
[ 20]
和车桥耦合系统
[ 21]
等大型复杂结构系统的随机振动分析。在随机振动分
析中,若已知随机激励的均值和互相关函数等统计特征,则可按统计矩的运算
规律根据时域显式表达式直接求得响应的前两阶统计矩,这一方法称为时域显
式直接法。若通过数值模拟手段生成大量随机激励样本,则可结合时域显式表
达式进行高效蒙特卡罗模拟,这一方法称为时域显式蒙特卡罗模拟法。
为了解决随机列车荷载作用下冻土路基动力响应分析这一工程问题,本文
首先提出一种随机列车荷载的简化模拟方法,以获取随机分析所需的荷载样本。
然后,通过有限元脉冲响应分析构建路基结构动力响应的时域显式表达式。在
此基础上,进一步利用时域显式蒙特卡罗模拟法,对随机列车荷载作用下冻土
路基的关键响应开展统计分析,可大幅提高冻土路基结构随机振动的计算效率。
最后,通过数值算例展示时域显式蒙特卡罗模拟法在冻土路基非平稳随机振动
分析中的计算精度与效率。
1 随机列车荷载模拟方法
在冻土路基随机振动问题中,需要列车驶过引起的随机荷载的统计特性或
大量样本作为输入。本文采用常见的二维路基模型
[ 7]
进行计算,此类模型需
要列车行驶引起的随机轨枕竖向作用力作为荷载输入
[ 15]
。目前,相关实测数
据较少,难以满足冻土路基随机振动分析的需求。为了获取分析所需输入,可
以对确定性列车荷载的数值模拟方法
[ 22-24 ]
加以改进,在其中引入随机参数,
进而生成大量随机列车荷载样本。
列车在轨道上行驶的简图如图 1 所示,其中列车行驶速度为 vv。为了描
述列车在沿轨道方向上的位置,建立如图 1 所示的直线坐标系 Ox。列车前端
的初始坐标假定为−L'0-L0'。
图 1
图 1 列车行驶简图
Fig. 1 Moving train schematic
考虑如图 1 所示的列车模型,其车厢的数量为 NN,包括机车和普通两类
车厢。在列车行驶过程中,每个轮对的荷载通过轨道及轨枕传递到路基顶面。
对于某一关注轨枕,列车对路基的竖向作用力 F(t)F(t)的表达式为
F(t)=∑i=1Nfi(t) (i=1,2,⋯,N)F(t)=∑i=1Nfi(t) (i=1,2,⋯,N)
(1)
式中:fi(t)fi(t)表示第 ii 节车厢的荷载传递到关注轨枕下方的部分,可进一
步表示为该车厢所有轮对产生的竖向轨枕力之和。若该车厢为机车车厢,则其
表达式为
fi(t)=∑j=1wPij(t)K(xij)fi(t)=∑j=1wPij(t)K(xij)
(2)
式中:ww 为车厢轮对数量,其中机车和普通两类车厢的轮对数分别为 6
和 4;Pij(t)Pij(t)为第 ii 节车厢第 jj 组轮对在轨道上产生的荷载;K(xij)K(xij)为
关注轨枕对轮对作用力的分担系数,按日本规范的经验假设轮对荷载由附近 5
根轨枕按 0.1∶0.2∶0.4∶0.2∶0.1 分担
[ 25]
,可进一步表示为
K(xij)=⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪0.4,0.2,0.1,0,∣∣xij−5Ls2∣∣≤Ls2Ls2<∣∣xij
−5Ls2∣∣≤3Ls23Ls2<∣∣xij−5Ls2∣∣≤5Ls2∣∣xij−5Ls2∣∣>5Ls2K(xij)=0.4,xij-
5Ls2≤Ls20.2,Ls2<xij-5Ls2≤3Ls20.1,3Ls2<xij-5Ls2≤5Ls20,xij-5Ls2>5Ls2
(3)
式中:LsLs 为轨枕间距,如图 2 所示;xijxij 为轮对在 Ox 坐标系中的坐标。
当第 ii 节车厢为机车车厢时,xijxij 具体表示为
图 2
图 2 道床顶面荷载沿线路方向分布图
Fig. 2 Distribution of loading on top surface of road bed along the road
⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪
⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪xi1=−L'0+vt−∑k=1iLd,k+L0+4ai+bixi2=−L'
0+vt−∑k=1iLd,k+L0+3ai+bixi3=−L'0+vt−∑k=1iLd,k+L0+2ai+bixi4=−L'0+vt−∑k=1
iLd,k+L0+2aixi5=−L'0+vt−∑k=1iLd,k+L0+aixi6=−L'0+vt−∑k=1iLd,k+L0xi1=-
L0'+vt-∑k=1iLd,k+L0+4ai+bixi2=-L0'+vt-∑k=1iLd,k+L0+3ai+bixi3=-L0'+vt-
∑k=1iLd,k+L0+2ai+bixi4=-L0'+vt-∑k=1iLd,k+L0+2aixi5=-L0'+vt-
∑k=1iLd,k+L0+aixi6=-L0'+vt-∑k=1iLd,k+L0
(4)
式中:aiai 第 ii 节车厢同一个转向架内相邻两个轮对的轴距;bibi 为第 ii
节车厢不同转向架相邻两个轮对的轴距;Ld,kLd,k 为第 kk 节车厢长度;L0L0
为车厢后端与该车厢最后一组轮对水平方向上的距离;−L'0-L0'为列车前端在
该坐标系内的初始坐标。当第 ii 节车厢为普通车厢时,也可以参考图 1 采用类
似的方式计算轮对坐标。
列车第 ii 节车厢第 jj 组轮对产生的作用力 Pij(t)Pij(t)可以采用反映行车平
稳性、附加动载和轨面波形磨耗效应三个控制条件的激励力来模拟
[ 24]
,其表
达式为
Pij(t)=P(0)ij+P(1)ijsinω1t+P(2)ijsinω2t+P(3)ijsinω3tPij(t)=Pij(0)+Pij(1)sinω1
t+Pij(2)sinω2t+Pij(3)sinω3t
(5)
式中:P(0)ijPij(0)为轮对静载;P(1)ijPij(1)、P(2)ijPij(2)和 P(3)ijPij(3)分别为反
映前述三个控制条件的某一典型值的振动荷载幅值;ω1、ω2ω1、ω2 和 ω3ω3
分别为反映控制条件的轨道几何不平顺曲线波长下的振动圆频率;tt 为荷载作
用时间。其中 P(k)ijPij(k)和 ωkωk 的表达式分别为
{P(k)ij=M(0)ijhkω2kωk=2πv/lk (k=1, 2, 3)Pij(k)=Mij(0)hkωk2ωk=2πv/lk
(k=1, 2, 3)
(6)
剩余22页未读,继续阅读
资源评论
罗伯特之技术屋
- 粉丝: 3694
- 资源: 1万+
下载权益
C知道特权
VIP文章
课程特权
开通VIP
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- 5G RAN NR TDD单用户峰值吞吐率计算过程.pptx
- 基于Objective-C开发的手势密码解锁和面容解锁(代码简洁高效通俗易懂)+源码+开发文档+源码解析(毕业设计&课程设计)
- LangChain学习之 Question And Answer的操作附带数据文件
- linux常用命令大全.zip
- CATIA入门操作案例-正八边异形带孔凸台绘制,等距点绘制正多边形,凸台绘制
- 大型代码语言模型的项目级提示生成pdf
- 纸片战争The War of paper.sb3
- testedtestedtested
- 基于C语言+STM32开发的智能门锁优秀项目+包含指纹识别、人脸识别、RFID解锁、密码解锁、蓝牙解锁功能+源码+项目解析
- 目标检测-智能零售柜商品检测数据集-5000张图-+对应VOC-COCO-YOLO三种格式标签+数据集划分脚本
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功