CN106650067A - 一种输出中高频响应的车‑轨‑桥耦合振动频域模型 - Google Patents
一种输出中高频响应的车‑轨‑桥耦合振动频域模型 Download PDFInfo
- Publication number
- CN106650067A CN106650067A CN201611136073.9A CN201611136073A CN106650067A CN 106650067 A CN106650067 A CN 106650067A CN 201611136073 A CN201611136073 A CN 201611136073A CN 106650067 A CN106650067 A CN 106650067A
- Authority
- CN
- China
- Prior art keywords
- wheel
- rail
- track
- contact point
- bridge
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/10—Noise analysis or noise optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Toys (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种输出中高频响应的车‑轨‑桥耦合振动频域模型,采用先选取车‑轨‑桥耦合振动仿真预测参数并输入车‑轨‑桥耦合模型中,根据车辆动力学方程及轨道‑桥梁耦合模型分别求解车轮导纳矩阵、轨道导纳矩阵,根据轮轨相互作用模型求解轮轨接触弹簧导纳矩阵并代入动态轮轨力公式求得动态轮轨力;最后将求得的动态轮轨力代入车辆运动方程、轨道‑桥梁耦合模型运动方程作为车辆模型和轨道‑桥梁耦合模型的激励源,进一步求解动态轮轨力引起的车辆、轨道和桥梁结构频域动力响应;该方法能够揭示车‑轨‑桥结构振动和能量产生与传递的机理及影响规律,从而为有针对性的采取减振降噪措施提供理论依据。
Description
技术领域
本发明涉及一种输出中高频响应的车-轨-桥耦合振动频域模型,可应用于高架轨道交通车辆、轨道及桥梁结构中、高频振动响应的预测,模型输出的上述中高频响应可作为轨道交通环境振动及噪声预测的激励源,进而为研究轨道交通环境振动及噪声问题提供理论支持。
背景技术
随着社会经济的发展和城市化进程的不断推进,交通问题日益突出。轨道交通作为一种大容量交通运输工具,能有效解决客流大密度出行问题,对经济发展的推动作用不容忽视。轨道交通高架线路由于建设周期短、费用低,具有很好的线路平顺性和稳定性,同时可以更为可靠地控制地基长期沉降等优点,在整个轨道交通线路中占有很大比例。但是,高架轨道交通车辆运行引起的环境振动问题及噪声辐射问题却在某种程度上降低了它的优势。随着人们生活水平的提高,公众对振动和噪声问题的容忍度不断降低,并由此引发了诸多社会问题。
轨道交通车辆在运行过程中,钢轨表面和车轮踏面之间的组合不平顺引起动态轮轨力,动态轮轨力向上传递给列车引起车辆振动,向下传递给轨道致使轨道结构振动,轨道结构的振动再经由轨道支承***传递给桥梁结构,进而引发桥梁结构的振动和噪声辐射问题。因而,准确建立车-轨-桥耦合振动分析模型,探明车-轨-桥耦合***中振动能量的产生和传递机理对控制高架轨道交通环境振动及噪声问题具有重要作用。
传统的对车-轨-桥耦合***进行建模和动力分析时通常建立时域模型,但时域模型主要用于分析列车运行时桥梁结构的强度,关注的桥梁结构振动频率通常低于20Hz,然而,高架轨道交通环境振动问题关注的振动频率范围为1~80Hz,而高架轨道交通噪声问题关注的频率范围则更高,若采用时域方法建立车-轨-桥耦合振动分析模型,则输出的响应难以用于对轨道交通环境振动及噪声进行准确预测。
因此,提出一种输出中高频响应的车-轨-桥耦合振动频域模型至关重要,通过该模型计算输出的响应更适合于高架轨道交通环境振动及噪声预测,该方法能够揭示车-轨-桥结构振动和能量产生与传递的机理及影响规律,从而为有针对性的采取减振降噪措施提供理论依据。
发明内容
本发明的目的是提出一种输出中高频响应的车-轨-桥耦合振动频域模型,解决传统的对车-轨-桥耦合***进行建模和动力分析时通常建立时域模型,但时域模型主要用于分析列车运行时桥梁结构的强度,关注的桥梁结构振动频率通常低于20Hz,然而,高架轨道交通环境振动问题关注的振动频率范围为1~80Hz,而高架轨道交通噪声问题关注的频率范围则更高,若采用时域方法建立车-轨-桥耦合振动分析模型,则输出的响应难以用于对轨道交通环境振动及噪声进行准确预测的问题。
本发明解决其技术问题所采用的技术方案为:
一种输出中高频响应的车-轨-桥耦合振动频域模型,包括以下几个步骤:
第一步、选取车-轨-桥耦合振动仿真预测所包含的车辆参数、轨道结构参数、轮轨组合粗糙度谱、桥梁结构参数及行车速度参数;
第二步、根据步骤(一)所选的车辆参数,建立10个自由度的车辆模型,且沿纵向的1/2车体、1/2转向架和车轮均考虑为刚体,考虑车体和转向架的沉浮、点头自由度以及车轮的沉浮自由度;
通过逐一对各刚体应用D′Alembert原理,并根据车辆***运动方程表达式:
其中,M、C和K分别为车辆***的质量、阻尼和刚度矩阵,z为车辆***的自由度向量,F为作用在车辆***上的外力;
在频域求解车辆***的稳态响应,令z=Z(ω)eiωt,可通过车辆运动方程求解车辆***频域响应Z(ω),进而根据位移导纳定义,可求得车轮位移导纳,其数学表达式为
式中Fj表示第j个轮轨接触点位置时的轮轨力,表示当轮轨力Fj作用于第j个轮轨接触点位置时第i个轮轨接触点处车轮的响应,Nw表示轮轨接触点的总数;导纳表示单位力作用于第j个轮轨接触点位置时第i个轮轨接触点处车轮的响应;钢轨上作用Nw个车轮,则车轮导纳可写成一个Nw×Nw的矩阵,其表达式为:
上式中,Nw即为轮轨接触点的总数;
第三步、再根据步骤(一)所选的轨道结构参数和桥梁结构参数,建立四层叠合梁模拟轨-桥耦合子***;
第四步、然后根据步骤(一)所选轮轨组合粗糙度谱及车轮和钢轨的参数,建立轮轨相互作用模型,轮轨接触弹簧即为车轮粗糙度弹簧和轨道粗糙度弹簧串联得到,则其刚度满足
式中KC为轮轨接触弹簧刚度,KCR为轨道粗糙度弹簧刚度,KCW为车轮粗糙度弹簧刚度;
即可得到轮轨接触弹簧导纳为αC=1/KC
同时可得轮轨接触弹簧导纳矩阵为
第五步、根据车-轨-桥耦合模型涉及的多个轮轨耦合作用,建立轮轨耦合模型叠加法,根据叠加法原理,首先设定仅第一车轮(1)与轨道之间存在不平顺激励,产生的轮轨力称为主动力P11,在第二车轮(2)、第三车轮(3)和第四车轮(4)与轨道接触点产生的力分别为被动力P21、P31和P41,将轮轨接触点的第二轮轨接触点(6)、第三轮轨接触点(7)和第四轮轨接触点(8)分别视为主动激励点,类推可得到第k个轮轨接触点激励作用下的第j个接触点处的主动力与被动力,其计算表达式为:
式中Pjk表示当激励作用于k点时j点处产生的力,Nw为轮轨接触点的个数,Fwrj表示第j个轮轨接触点处的动态轮轨力;
因此可以得出轮轨作用力,以数学表达式以向量的形式表示为
Fwr={Fwrj} j=1,2,...,Nw
式中:Nw为轮轨接触点的个数,Fwr表示各轮轨接触点的动态轮轨力向量,Fwrj表示第j个轮轨接触点处的动态轮轨力。
第六步、根据车轮模型轴距和定距的存在,不同轮轨接触点之间的激励出现时间滞后关系,其轮轨接触点的不平顺的数学表达式为
式中,Δ(ω)为Nw个轮轨接触点处的轮轨粗糙度向量,△1(ω)为第1个轮轨接触点处的轮轨粗糙度值,i为虚数单位,ω为圆频率,ti为第i个轮轨接触点与第1个轮轨接触点之间由于距离的存在产生的时间差。
第七步、根据步骤(二)得到的车轮导纳矩阵αw、步骤(三)得到的轨道结构导纳矩阵αR、步骤(四)得到的轮轨接触弹簧导纳矩阵αC、步骤(五)得到多轮轨相互作用的动态轮轨力关系以及步骤(六)中的Nw个轮轨接触点处的不平顺向量Δ(ω),可求得动态轮轨力,其数学表达式为
Fwr(ω)=-(αw+αR+αC)-1Δ(ω)
式中,Fwr(ω)为各轮轨接触点处的动态轮轨力向量;αw为车轮导纳矩阵;αR为轨道的钢轨导纳矩阵;为轮轨接触导纳矩阵;Δ(ω)={△j},j=1,2,...,Nw为轮轨接触点的轨道粗糙度向量。
第八步、将步骤(七)得到的动态轮轨力作为激励力,施加到步骤(二)中的车辆运动方程中及步骤(三)中的轨-桥耦合模型运动方程中进一步求得车辆、轨道结构及桥梁结构的频域位移响应Z(ω)。
第九步、根据速度及加速度与位移之间的关系式,进一步求得车辆、轨道结构及桥梁结构的速度及加速度频域响应;其数学表达式为:
式中,i为虚数单位,ω为圆频率,Z(ω)表示频率位移响应,表示频率速度响应,表示频率加速度响应。
采用本发明的技术方案,采用一种输出中高频响应的车-轨-桥耦合振动频域模型,首先,选取车-轨-桥耦合振动仿真预测参数,将其输入车-轨-桥耦合模型中,根据车辆动力学方程及轨道-桥梁耦合模型分别求解车轮导纳矩阵、轨道导纳矩阵,根据轮轨相互作用模型求解轮轨接触弹簧导纳矩阵;将求得的车轮导纳矩阵、轨道导纳矩阵和轮轨接触弹簧导纳矩阵代入动态轮轨力公式求得动态轮轨力;最后将求得的动态轮轨力代入车辆运动方程、轨道-桥梁耦合模型运动方程作为车辆模型和轨道-桥梁耦合模型的激励源,进一步求解动态轮轨力引起的车辆、轨道和桥梁结构频域动力响应;该方法能够揭示车-轨-桥结构振动和能量产生与传递的机理及影响规律,从而为有针对性的采取减振降噪措施提供理论依据。
与现有技术相比,本发明的有益效果是,从车-轨-桥耦合***振动产生的机理出发,先求解动态轮轨力频谱,再以动态轮轨力作为激励求解车辆、轨道、桥梁结构的频域动力响应;研究关注的频段范围为20~500Hz,属于中高频段,输出的***响应可用于进一步研究高架轨道交通环境振动和噪声问题,与现有常规的车-轨-桥耦合振动时域模型相比,一是更加准确地考虑了动态轮轨力的中高频激励影响;二是避开时域方法所须采用的分步积分,计算效率大大提高,同时频域结果保证了满足工程要求的计算精度;三是该预测模型适用性强,能应用与不同的车辆类型、轨道结构类型、桥梁结构类型。此外,由于该模型计算输出的响应可作为高架轨道交通环境振动和噪声问题的激励源,通过对该模型的参数研究,从而指导轨道形式选择及桥梁结构设计,对高架轨道交通环境振动及噪声控制具有重要理论意义。
综上所述,采用本发明的一种输出中高频响应的车-轨-桥耦合振动频域模型可以解决传统的对车-轨-桥耦合***进行建模和动力分析时通常建立时域模型,但时域模型主要用于分析列车运行时桥梁结构的强度,关注的桥梁结构振动频率通常低于20Hz,然而,高架轨道交通环境振动问题关注的振动频率范围为1~80Hz,而高架轨道交通噪声问题关注的频率范围则更高,若采用时域方法建立车-轨-桥耦合振动分析模型,则输出的响应难以用于对轨道交通环境振动及噪声进行准确预测的问题。
附图说明
图1为本发明的车-轨-桥耦合振动分析模型;
图2为本发明的四层叠合梁模拟的轨-桥耦合模型;
图3为本发明的Remington轮轨相互作用模型;
图4为本发明的多车轮作用情况下的轮轨耦合模型叠加法;
图5为本发明的车-轨-桥耦合模型建模流程图;
图6为本发明的列车过桥照片;
图7为本发明的直线桥梁段CRTSⅡ型板式无砟轨道横截面示意图;
图8为本发明的32m双线混凝土简支箱梁典型横截面图;
图9为本发明的京沪高铁实测轨道不平顺谱;
图10为本发明的第一车轮位置处作用单位力时的车轮导纳;
图11为本发明的第二车轮位置处作用单位力时的车轮导纳;
图12为本发明的第三车轮位置处作用单位力时的车轮导纳;
图13为本发明的第四车轮位置处作用单位力时的车轮导纳;
图14为本发明的第一车轮位置处作用单位力时的轨道结构导纳;
图15为本发明的第二车轮位置处作用单位力时的轨道结构导纳;
图16为本发明的第三车轮位置处作用单位力时的轨道结构导纳;
图17为本发明的第四车轮位置处作用单位力时的轨道结构导纳;
图18为本发明的各轮轨接触点处的动态轮轨力幅值和相位;
图19各车速下车辆加速度幅值对比;
图20各车速下桥梁跨中加速度的幅值对比;
图21为F轴/2示意图。
具体实施方式
下面结合附图和实施例对本发明的技术方案作进一步说明。
如图1~21所示,一种输出中高频响应的车-轨-桥耦合振动频域模型,包括以下几个步骤:
第一步、选取车-轨-桥耦合振动仿真预测所包含的车辆参数、轨道结构参数、轮轨组合粗糙度谱、桥梁结构参数及行车速度参数;
第二步、根据步骤(一)所选的车辆参数,建立10个自由度的车辆模型,且沿纵向的1/2车体、1/2转向架和车轮均考虑为刚体,考虑车体和转向架的沉浮、点头自由度以及车轮的沉浮自由度;
通过逐一对各刚体应用D′Alembert原理,并根据车辆***运动方程表达式:
其中,M、C和K分别为车辆***的质量、阻尼和刚度矩阵,z为车辆***的自由度向量,F为作用在车辆***上的外力;
在频域求解车辆***的稳态响应,令z=Z(ω)eiωt,可通过车辆运动方程求解车辆***频域响应Z(ω),进而根据位移导纳定义,可求得车轮位移导纳,其数学表达式为
式中,Fj表示第j个轮轨接触点位置时的轮轨力,表示当轮轨力Fj作用于第j个轮轨接触点位置时第i个轮轨接触点处车轮的响应;导纳表示单位力作用于第j个轮轨接触点位置时第i个轮轨接触点处车轮的响应;若钢轨上作用Nw个车轮,则车轮导纳可写成一个Nw×Nw的矩阵,其表达式为:
式中,Nw即为轮轨接触点的总数;
第三步、再根据步骤(一)所选的轨道结构参数和桥梁结构参数,建立四层叠合梁模拟轨-桥耦合子***;
第四步、然后根据步骤(一)所选轮轨组合粗糙度谱及车轮和钢轨的参数,建立轮轨相互作用模型,轮轨接触弹簧即为车轮粗糙度弹簧和轨道粗糙度弹簧串联得到,则其刚度满足
式中KC为轮轨接触弹簧刚度,KCR为轨道粗糙度弹簧刚度,KCW为车轮粗糙度弹簧刚度;
即可得到轮轨接触弹簧导纳为αC=1/KC
同时可得轮轨接触弹簧导纳矩阵为
式中,为Nw阶单位阵。
第五步、根据车-轨-桥耦合模型涉及的多个轮轨耦合作用,建立轮轨耦合模型叠加法,根据叠加法原理,首先设定仅第一车轮1与轨道之间存在不平顺激励,产生的轮轨力称为主动力P11,在第二车轮2、第三车轮3和第四车轮4与轨道接触点产生的力分别为被动力P21、P31和P41,将轮轨接触点的第二轮轨接触点(6)、第三轮轨接触点(7)和第四轮轨接触点(8)分别视为主动激励点,类推可得到Nw个轮轨接触点时,第k个轮轨接触点激励作用下的第j个轮轨接触点处的总轮轨力,其计算表达式为:
式中Nw为轮轨接触点个数,Pjk表示当激励作用于k点时j点处产生的相互作用力,Fwrj表示第j个轮轨接触点处的动态轮轨力;
因此可以得出轮轨作用力,以数学表达式以向量的形式表示为
Fwr(ω)={Fwrj} j=1,2,...,Nw
式中Fwr(ω)表示Nw个轮轨接触点的动态轮轨力向量;
第六步、由于车轮模型轴距和定距的存在,不同轮轨接触点之间的激励存在时间滞后关系,各轮轨接触点粗糙度的数学表达式为
式中Δ(ω)表示各个轮轨接触点处的轮轨粗糙度向量,△1(ω)表示第1个轮轨接触点处的轮轨粗糙度值,为虚数单位,ω为圆频率,ti为第i个轮轨接触点与第1个轮轨接触点间距离产生的时间差。
第七步、根据步骤(二)得到的车轮导纳矩阵αw、步骤(三)得到的轨道结构导纳矩阵αR、步骤(四)得到的轮轨接触弹簧导纳矩阵αC、步骤(五)得到多轮轨相互作用的动态轮轨力关系以及步骤(六)中的Nw个轮轨接触点处的不平顺向量Δ(ω),可求得动态轮轨力Fwr(ω),其数学表达式为
Fwr(ω)=-(αw+αR+αC)-1Δ(ω)
式中,Fwr(ω)动态轮轨力向量;αw为车轮导纳矩阵;αR为轨道的钢轨导纳矩阵;αC=INw×Nwαc为轮轨接触导纳矩阵,Δ(ω)={△j},j=1,2,...,Nw为轮轨接触点的轨道不平顺向量。
第八步、将步骤(七)得到的动态轮轨力作为激励力,施加到步骤(二)中的车辆运动方程中及步骤(三)中的轨-桥耦合模型运动方程中进一步求得车辆、轨道结构及桥梁结构的频域位移响应Z(ω)。
第九步、根据速度及加速度与位移之间的关系式,进一步求得车辆、轨道结构及桥梁结构的速度及加速度频域响应。其数学表达式为:
式中,i为虚数单位,ω为圆频率,Z(ω)表示频率位移响应,表示频率速度响应,表示频率加速度响应。
试验时以某客运专线高速列车通过32m混凝土简支箱梁桥为例,高速列车过桥如图6所示,给出本发明方法对车-轨-桥耦合振动仿真预测的具体过程。
(1)确定车-轨-桥耦合振动仿真预测参数,主要为列车类型(CRH380A)及相关参数。如表1所示、轨道结构类型(CRTSⅡ型版式无砟轨道如图7所示)、桥梁结构类型(32m混凝土简支箱梁如图8所示)及相关参数如表2所示、轮轨粗糙度谱参数(采用京沪高铁实测轨道谱如图9所示)、车辆运行速度(350km/h)。
其参数如下表所示:
表1 车辆动力学参数(单根钢轨上支承的车辆)
表2 4层叠合梁模拟的轨道-桥梁耦合模型动力学参数
(2)将车辆模型参数如表1所示,代入车辆运动方程首先求出车轮导纳矩阵,得到车轮导纳如图10~13所示。
(3)将轨-桥耦合模型参数如表2所示,代入轨-桥耦合振动方程求出轨道结构导纳矩阵,得到轨道结构导纳如图14~17所示。
(4)将车轮、钢轨参数代入轮轨接触弹簧导纳计算式,求得轮轨接触弹簧导纳矩阵。
(5)将车轮导纳矩阵、轨道结构导纳矩阵和轮轨接触弹簧导纳矩阵以及轮轨粗糙度谱参数代入动态轮轨力计算式,求出轮轨接触点处的动态轮轨力,如图18所示。
(6)将上一步求得的动态轮轨力作为激励力输入车辆运动方程和轨-桥耦合运动方程,进而求得车辆***动力响应如图19所示和桥梁结构动力响应如图20所示。
其具体计算方法包括如下步骤:
(1)明确车-轨-桥耦合振动频域模型预测参数,所述预测参数包括车辆参数、轨道结构参数、轮轨组合粗糙度谱、桥梁结构参数及行车速度;如下表所示:
表1 车辆动力学参数(单根钢轨上支承的车辆)
表2 4层叠合梁模拟的轨道-桥梁耦合模型动力学参数
(2)根据步骤(1)所选的车辆参数,建立如图21所示单根钢轨上作用如图1所示10个自由度的车辆模型,沿纵向的1/2车体、1/2转向架和车轮均考虑为刚体,考虑了车体和转向架的沉浮、点头自由度以及车轮的沉浮自由度。
通过逐一对各刚体应用D’Alembert原理,可得车辆***运动方程为:
其中,M、C和K分别为车辆***的质量、阻尼和刚度矩阵,z为车辆***的自由度向量,F为作用在车辆***上的外力,其具体表达式如下
z={zc βc zt1 βt1 zt2 βt2 zw1 zw2 zw3 zw4}T (2)
F={0 0 0 0 0 0 -Fwr1 -Fwr2 -Fwr3 -Fwr4}T (3)
M=diag(mc,Ic.mb,Ib,mb,Ib,mw,mw,mw,mw) (4)
以上车辆模型共考虑了10个独立自由度,分别为zc,βc,zt1,βt1,zt2,βt2,zw1,zw2,zw3,zw4,其中z代表沉浮自由度、β代表点头自由度,下标c代表车体,t1和t2分别代表前后两个转向架,w1~w4分别代表4个车轮。mc、mt和mw分别为1/2车体、1/2转向架和车轮质量,Ic、It分别为1/2车体和1/2转向架的点头转动惯量,k1、k2、c1和c2分别为车辆一系、二系悬挂***刚度、阻尼系数,Lt为车辆定距之半,Lw分车辆轴距之半,Lc为车辆钩到钩的距离。
在频域求解车辆***的稳态响应,令z=Z(ω)eiωt,可通过车辆运动方程求解车辆***频域响应Z(ω)。如图1所示有4个轮轨接触点时,根据位移导纳定义,可求得车轮位移导纳为
式中,表示当轮轨力Fj作用于第j个轮轨接触点位置时,第i个轮轨接触点处车轮的响应;导纳表示单位力作用于第j个轮轨接触点位置时,第i个轮轨接触点处车轮的响应。如图1所示有4个轮轨接触点时,则车轮导纳可写成一个4×4的矩阵
(3)根据步骤(1)所选的轨道结构参数和桥梁结构参数,建立四层叠合梁模拟模拟轨-桥耦合子***,如图2所示。四层叠合梁由上到下分别模拟钢轨、轨道板、底座板和桥梁,各层梁之间通过弹性支承进行联结,共包括四层弹性支承,由上到下分别表示扣件层、CA砂浆层、“两布一膜”滑动层和桥梁支座。为方便描述,钢轨、轨道板、底座板和箱梁的单位长度质量、抗弯刚度和损耗因子分别记为m、EI和η,并以下标“r”、“s”、“d”和“b”分别表示钢轨、轨道板、底座板和箱梁;跨度L范围内有ns块长度为Ls的轨道板,不计相邻轨道板间距;弹性支承弹簧刚度、损耗因子和间距分别记为k、η和d,并以下标“f”、“c”、“h”和“bs”分别表示钢轨扣件、CA砂浆层、“两布一膜”滑动层和桥梁支座,需指出的是dbs表示桥梁支座距离箱梁最近梁端距离(沿纵向);每块轨道板上部和下部的离散弹簧数目分别记为Nsf和Nsc,每跨箱梁跨度L范围内的扣件、CA砂浆层、“两布一膜”滑动层和桥梁支座离散弹簧总数分别记为nf、nc、nh和nbs,则应满足nf=ns×Nsf和nc=ns×Nsc。根据轨道结构的几何参数,其应满足关系式L=ns×Ls和Ls=Nsf×df=Nsc×dc。
建立如图2所示的xoz坐标系,坐标原点o位于叠合梁单元左端端点处。钢轨扣件弹簧的坐标记为xfi,i=1,…,nf;CA砂浆层弹簧的坐标记为xci,i=1,…,nc;“两布一膜”滑动层弹簧的坐标记为xhi,i=1,…,nh;桥梁支座弹簧的坐标记为xbsi,i=1,…,nbs。第h块轨道板上部和下部弹簧的局部坐标分为记为和h=1,…,ns;i=1,…,Nsf;j=1,…,Nsc。对于扣件总体坐标系坐标xfm与局部坐标系坐标当m=(h-1)×Nsf+n时,坐标xfm和对应同一个钢轨扣件;对于CA砂浆层弹簧总体坐标系坐标xcm与局部坐标系坐标当m=(h-1)×Nsc+n时,坐标xcm和对应同一个CA砂浆层弹簧。
对轨道结构中各层叠合梁分别取出进行分析,各层叠合梁的垂向振动位移响应与其上作用的荷载激励有关,根据叠加原理,各层梁体的垂向振动响应由各荷载激励引起的振动响应叠加得到。结合导纳法的思想,钢轨、轨道板、底座板和桥梁的垂向位移响应可分别表示为
式(9)~(12)中,Gr(x,xP)、Gs(x,xP)、Gd(x,xP)和Gb(x,xP)分别表示钢轨、轨道板、底座板和桥梁在xP处作用单位1的激励时x处的位移响应,即为各层梁体的位移导纳函数;F为钢轨上作用的外荷载激励;Pfi表示钢轨下部第i个扣件的扣件力;分别表示第h块轨道板上部的第j个扣件力和下部的第k个CA砂浆层弹性支承力;Pci、Phi分别表示底座板上部的第i个CA砂浆层弹性支承力和下部的第i个滑动层弹性支承力;Pbsi表示桥梁下部的第i个支座反力。
钢轨扣件力Pfi或CA砂浆层弹性支承力Pci或滑动层弹性支承力Phi和桥梁支座反力Pbsi可由其自身复刚度与其变形量的乘积求得。
对于钢轨扣件力,当i=(h-1)×Nsf+j时,Pfi和对应同一扣件的扣件力,可表示为
式中,i=1,2,…,nf;h=1,2,…,ns;j=1,2,…,Nsf。
对于CA砂浆层弹性支承力,当i=(h-1)×Nsc+k时,Pci和对应同一CA砂浆层弹性支承力,可表示为
式中,i=1,2,…,nc;h=1,2,…,ns;k=1,2,…,Nsc。
滑动层弹性支承力Phi可表示为
Phi=kh(1+iηh)[Zd(xhi)-Zb(xhi)],i=1,2,…,nh (15)
桥梁支座反力Pbsi可表示为
Pbsi=kbs(1+iηbs)Zb(xbsi),i=1,…,nbs (16)
由式(9)~式(12)求出式(13)~式(16)中钢轨、轨道板、底座板和桥梁在各弹性支承点处的位移响应,可表示为
将式(17)~(23)代入式(13)~(16)后可得到
i=1,2,…,nf且i=(h-1)×Nsf+j,h=1,2,…,ns,j=1,2,…,Nsf
i=1,2,…,nc且i=(h-1)×Nsc+k,h=1,2,…,ns,k=1,2,…,Nsc
将式(24)等号右端第2~4项移到方程左端,同时将式(25)~(27)等号右端项均移到方程左端,经整理将各方程合并写成矩阵形式得
[A(ω)]{P}={U} (28)
式中,[A(ω)]为由钢轨、轨道板、底座板和桥梁的位移导纳以及钢轨扣件、CA砂浆层、滑动层和支座的刚度、阻抗参数等构成的(nf+nc+nh+nbs)×(nf+nc+nh+nbs)维动柔度矩阵;{P}为由钢轨扣件力、CA砂浆层弹性支承力、滑动层弹性支承力和桥梁支座反力等组成的(nf+nc+nh+nbs)维待求列向量;{U}为外荷载激励引起的各层叠合梁各弹性支承处的(nf+nc+nh+nbs)维位移响应,其中仅前nf个元素非0,为FGr(xfi,xF),i=1,2,…,nf,其余元素为0。
由式(28)求出轨道桥梁耦合***中各弹性支承力及支座反力列向量{P},将其代入式(9)计算钢轨上任意位置处由外荷载F作用引起的位移响应Zr(xi)。该响应与外荷载F的比值即为轨-桥耦合***的位移导纳,记为
上式中,表示钢轨上xj处单位力,引起钢轨上xi处的位移响应。当i=j时,为钢轨位移点导纳,当i≠j时,为钢轨位移传递导纳。
如图1所示有4个轮轨接触点时,则轨道结构的钢轨位移导纳矩阵为
(4)根据步骤(1)所选轮轨组合粗糙度谱及车轮和钢轨等的参数,建立轮轨相互作用模型,如图3所示,轮轨接触弹簧即为车轮粗糙度弹簧和轨道粗糙度弹簧串联得到,则其刚度满足
式中,KC为轮轨接触弹簧刚度,KCR为轨道粗糙度弹簧刚度,KCW为车轮粗糙度弹簧刚度。
轮轨接触弹簧导纳即为αC=1/KC
如图1所示有4个轮轨接触点时,则可得轮轨接触弹簧导纳矩阵为
αC=I4×4αC (32)
(5)如图1所建立的车-轨-桥耦合模型涉及多个轮轨耦合作用,建立轮轨耦合模型叠加法示意图,如图4所示。由图4可知,根据叠加法原理,首先设定仅车轮1与轨道之间存在不平顺激励,由其产生的轮轨力称为主动力P11,在其余车轮2、3、4与轨道接触点产生的力称为被动力P21、P31、P41,此处的Pjk表示当激励作用于k点时j点处产生的力。同理,将轮轨接触点2、3、4分别视为主动激励点,可得到各轮轨接触点激励作用下的各接触点处的主动力与被动力。通过叠加,可得到各轮轨接触处总的轮轨作用力
轮轨作用力以向量的形式表示,有
Fwr={Fwrj} j=1,2,3,4 (34)
(6)由于车轮模型轴距和定距的存在,不同轮轨接触点之间的激励出现时间滞后关系,图4模型中4个轮轨接触点的不平顺可表示为
Δ(t)={△1(t-t1) △2(t-t2) △3(t-t3) △4(t-t4)}T (35)
式中,设定t1=0,则t2=2lw/V,t3=2lb/V,t4=2(lb+lw)/V为各车轮与车轮1之间的时间差。其中,V为车速,lw和lb分别为车轮轴距和定距之半。
假设时域的轨道不平顺为△(t)=△(ω)eiωt,则可得
(7)根据步骤(2)得到的车轮导纳矩阵αw、步骤(3)得到的轨道结构导纳矩阵αR、步骤(4)得到的轮轨接触弹簧导纳矩阵αC、步骤(5)得到多轮轨相互作用的动态轮轨力关系以及步骤(6)中的Nw个轮轨接触点处的不平顺向量Δ(ω),可求得动态轮轨力为
Fwr=-(αw+αR+αC)-1Δ (37)
式中,αw为车轮导纳矩阵;αR为轨道的钢轨导纳矩阵;αC=I4×4αc为轮轨接触导纳矩阵,Δ={△j},j=1,2,3,4为轮轨接触点的轨道不平顺向量。
(8)将步骤(7)得到的轮轨力作为激励力,施加到步骤(2)中的车辆运动方程中及步骤(3)中的轨-桥耦合模型运动方程中进一步求得车辆、轨道结构及桥梁结构的频域位移响应Z(ω)。
(9)根据速度及加速度与位移之间的关系式,进一步求得车辆、轨道结构及桥梁结构的速度及加速度频域响应。
其中,i为虚数单位,ω为圆频率,Z(ω)、和分别表示频率位移响应、频率速度响应和频率加速度响应。
与现有技术相比,本发明的有益效果是,从车-轨-桥耦合***振动产生的机理出发,先求解动态轮轨力频谱,再以动态轮轨力作为激励求解车辆、轨道、桥梁结构的频域动力响应;研究关注的频段范围为20~500Hz,属于中高频段,输出的***响应可用于进一步研究高架轨道交通环境振动和噪声问题,与现有常规的车-轨-桥耦合振动时域模型相比,一是更加准确地考虑了动态轮轨力的中高频激励影响;二是避开时域方法所须采用的分步积分,计算效率大大提高,同时频域结果保证了满足工程要求的计算精度;三是该预测模型适用性强,能应用与不同的车辆类型、轨道结构类型、桥梁结构类型。此外,由于该模型计算输出的响应可作为高架轨道交通环境振动和噪声问题的激励源,通过对该模型的参数研究,从而指导轨道形式选择及桥梁结构设计,对高架轨道交通环境振动及噪声控制具有重要理论意义。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (1)
1.一种输出中高频响应的车-轨-桥耦合振动频域模型,其特征在于:包括以下几个步骤:
第一步、选取车-轨-桥耦合振动仿真预测所包含的车辆参数、轨道结构参数、轮轨组合粗糙度谱、桥梁结构参数及行车速度参数;
第二步、根据步骤(一)所选的车辆参数,建立10个自由度的车辆模型,且沿纵向的1/2车体、1/2转向架和车轮均考虑为刚体,考虑车体和转向架的沉浮、点头自由度以及车轮的沉浮自由度;
通过逐一对各刚体应用D′Alembert原理,并根据车辆***运动方程表达式:
其中,M、C和K分别为车辆***的质量、阻尼和刚度矩阵,z为车辆***的自由度向量,F为作用在车辆***上的外力;
在频域求解车辆***的稳态响应,令z=Z(ω)eiωt,可通过车辆运动方程求解车辆***频域响应Z(ω),进而根据位移导纳定义,可求得车轮位移导纳,其数学表达式为
式中Fj表示第j个轮轨接触点位置处的轮轨力,表示当轮轨力Fj作用于第j个轮轨接触点位置时第i个轮轨接触点处车轮的响应;表示单位力作用于第j个轮轨接触点位置时第i个轮轨接触点处车轮的响应;
钢轨上作用Nw个车轮,则车轮导纳可写成一个Nw×Nw的矩阵,其表达式为:
式中,Nw即为轮轨接触点的总数;
第三步、再根据步骤(一)所选的轨道结构参数和桥梁结构参数,建立四层叠合梁模拟轨-桥耦合子***;
第四步、然后根据步骤(一)所选轮轨组合粗糙度谱及车轮和钢轨的参数,建立轮轨相互作用模型,轮轨接触弹簧即为车轮粗糙度弹簧和轨道粗糙度弹簧串联得到,则其刚度满足
式中KC为轮轨接触弹簧刚度,KCR为轨道粗糙度弹簧刚度,KCW为车轮粗糙度弹簧刚度;
可得到轮轨接触弹簧导纳为αC=1/KC
同时可得轮轨接触弹簧导纳矩阵为
式中为Nw阶单位阵;
第五步、假设有4个轮轨接触点,根据车-轨-桥耦合模型涉及的多个轮轨耦合作用,建立轮轨耦合模型叠加法,根据叠加法原理,首先设定仅第一车轮(1)与轨道之间存在不平顺激励,产生的相互作用力为主动力P11,在第二车轮(2)、第三车轮(3)和第四车轮(4)与轨道接触点产生的相互作用力分别为被动力P21、P31和P41,将轮轨接触点的第二轮轨接触点(6)、第三轮轨接触点(7)和第四轮轨接触点(8)分别视为主动激励点,类推可得到有Nw个轮轨接触点时,第k个轮轨接触点激励作用下的第j个接触点处的主动力与被动力,其计算表达式为:
式中Nw表示轮轨接触点个数,Pjk表示当激励作用于k点时j点处产生的相互作用力,Fwrj表示第j个轮轨接触点处的动态轮轨力;
因此可以得出轮轨作用力,数学表达式以向量的形式表示为
Fwr={Fwrj} j=1,2,...,Nw
式中,Fwr表示Nw个轮轨接触点的动态轮轨力向量;
第六步、根据车轮模型轴距和定距的存在,不同轮轨接触点之间的激励出现时间滞后关系,其轮轨接触点的不平顺的数学表达式为
式中,Δ(ω)为Nw个轮轨接触点处的轮轨粗糙度向量,△1(ω)为第1个轮轨接触点处的轮轨粗糙度值,i为虚数单位,ω为圆频率,ti为第i个轮轨接触点与第1个轮轨接触点之间由于距离的存在产生的时间差;
第七步、根据步骤(二)得到的车轮导纳矩阵αw、步骤(三)得到的轨道结构导纳矩阵αR、步骤(四)得到的轮轨接触弹簧导纳矩阵αC、步骤(五)得到多轮轨相互作用的动态轮轨力关系以及步骤(六)中的Nw个轮轨接触点处的不平顺向量Δ(ω),可求得动态轮轨力,其数学表达式为
Fwr=-(αw+αR+αC)-1Δ
式中,Fwr为各轮轨接触点处的动态轮轨力;αw为车轮导纳矩阵;αR为轨道的钢轨导纳矩阵;αC为轮轨接触导纳矩阵,Δ为轮轨接触点的轨道不平顺向量;
第八步、将步骤(七)得到的轮轨力作为激励力,施加到步骤(二)中的车辆运动方程中及步骤(三)中的轨-桥耦合模型运动方程中进一步求得车辆、轨道结构及桥梁结构的频域位移响应Z(ω);
第九步、根据速度及加速度与位移之间的关系式,进一步求得车辆、轨道结构及桥梁结构的速度及加速度频域响应;其数学表达式为:
其中,i为虚数单位,ω为圆频率,Z(ω)表示频率位移响应,表示频率速度响应,表示频率加速度响应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611136073.9A CN106650067B (zh) | 2016-12-12 | 2016-12-12 | 一种输出中高频响应的车-轨-桥耦合振动频域分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611136073.9A CN106650067B (zh) | 2016-12-12 | 2016-12-12 | 一种输出中高频响应的车-轨-桥耦合振动频域分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106650067A true CN106650067A (zh) | 2017-05-10 |
CN106650067B CN106650067B (zh) | 2019-11-22 |
Family
ID=58824179
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611136073.9A Active CN106650067B (zh) | 2016-12-12 | 2016-12-12 | 一种输出中高频响应的车-轨-桥耦合振动频域分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106650067B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107727338A (zh) * | 2017-06-01 | 2018-02-23 | 重庆大学 | 一种基于车桥耦合***的桥梁损伤诊断方法 |
CN108256278A (zh) * | 2018-03-13 | 2018-07-06 | 西南交通大学 | 一种悬挂式单轨车辆耦合动力学仿真***及其仿真方法 |
CN108763646A (zh) * | 2018-04-25 | 2018-11-06 | 华东交通大学 | 基于fe-sea混合法的箱型梁局部振动预测方法 |
CN108763716A (zh) * | 2018-05-22 | 2018-11-06 | 西南交通大学 | 一种基于振型叠加法的过桥列车倾覆危险性评价方法 |
CN109459129A (zh) * | 2018-09-27 | 2019-03-12 | 成都理工大学 | 一种高效求解高速列车过桥振动随机特性的分析方法 |
CN109556812A (zh) * | 2018-12-13 | 2019-04-02 | 佛山科学技术学院 | 一种行人与人行桥横向相互作用动力模拟装置 |
CN110376452A (zh) * | 2019-06-14 | 2019-10-25 | 中国科学院西安光学精密机械研究所 | 基于耦合机电分析的压电陶瓷促动器电噪声指标确定方法 |
CN111651840A (zh) * | 2020-04-09 | 2020-09-11 | 北京交通大学 | 一种基于深度学习技术的轨道板上拱状态检测方法 |
CN113065187A (zh) * | 2021-04-06 | 2021-07-02 | 华东交通大学 | 一种基于解析映射关系的震后高铁连续梁桥墩不均匀变形阈值计算方法 |
CN113432815A (zh) * | 2021-01-26 | 2021-09-24 | 重庆大学 | 一种基于测量车振动响应的桥面响应重构方法 |
CN114186362A (zh) * | 2021-11-30 | 2022-03-15 | 重庆大学 | 一种基于运营列车振动响应预测轨道不平顺及其传感器优化配置的分析方法 |
CN114239347A (zh) * | 2021-11-30 | 2022-03-25 | 重庆大学 | 一种基于运营列车振动信号实时量测铁路轨道支撑刚度的高效方法 |
CN114707203A (zh) * | 2022-03-07 | 2022-07-05 | 东南大学 | 一种桥梁减振降噪快速化预测评估方法 |
CN114997020A (zh) * | 2022-05-20 | 2022-09-02 | 重庆大学 | 一种基于多自由度车辆动态响应的接触点反演算法 |
CN115808280A (zh) * | 2022-11-16 | 2023-03-17 | 中国矿业大学 | 面向轨道减振评估的双轮对-构架多模式激励设备及方法 |
CN117172136A (zh) * | 2023-11-02 | 2023-12-05 | 成都理工大学 | 基于ssa-lstm算法的风车桥***振动响应预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103150458A (zh) * | 2013-04-01 | 2013-06-12 | 中南大学 | 车辆-轨道-桥梁-地基基础耦合***及其动力分析方法 |
CN103853896A (zh) * | 2014-03-20 | 2014-06-11 | 李小珍 | 铁路桥梁结构车致振动响应计算方法 |
CN104036087A (zh) * | 2014-06-24 | 2014-09-10 | 同济大学 | 基于功率流-边界元模型的高架轨道交通振动噪声仿真预测方法 |
CN105631128A (zh) * | 2015-12-29 | 2016-06-01 | 西南交通大学 | 一种高速铁路弓-网-车-轨垂向耦合大***动力学建模仿真方法 |
-
2016
- 2016-12-12 CN CN201611136073.9A patent/CN106650067B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103150458A (zh) * | 2013-04-01 | 2013-06-12 | 中南大学 | 车辆-轨道-桥梁-地基基础耦合***及其动力分析方法 |
CN103853896A (zh) * | 2014-03-20 | 2014-06-11 | 李小珍 | 铁路桥梁结构车致振动响应计算方法 |
CN104036087A (zh) * | 2014-06-24 | 2014-09-10 | 同济大学 | 基于功率流-边界元模型的高架轨道交通振动噪声仿真预测方法 |
CN105631128A (zh) * | 2015-12-29 | 2016-06-01 | 西南交通大学 | 一种高速铁路弓-网-车-轨垂向耦合大***动力学建模仿真方法 |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107727338B (zh) * | 2017-06-01 | 2019-08-23 | 重庆大学 | 一种基于车桥耦合***的桥梁损伤诊断方法 |
CN107727338A (zh) * | 2017-06-01 | 2018-02-23 | 重庆大学 | 一种基于车桥耦合***的桥梁损伤诊断方法 |
CN108256278A (zh) * | 2018-03-13 | 2018-07-06 | 西南交通大学 | 一种悬挂式单轨车辆耦合动力学仿真***及其仿真方法 |
CN108256278B (zh) * | 2018-03-13 | 2023-08-08 | 西南交通大学 | 一种悬挂式单轨车辆耦合动力学仿真***及其仿真方法 |
CN108763646A (zh) * | 2018-04-25 | 2018-11-06 | 华东交通大学 | 基于fe-sea混合法的箱型梁局部振动预测方法 |
CN108763646B (zh) * | 2018-04-25 | 2022-11-15 | 华东交通大学 | 基于fe-sea混合法的箱型梁局部振动预测方法 |
CN108763716A (zh) * | 2018-05-22 | 2018-11-06 | 西南交通大学 | 一种基于振型叠加法的过桥列车倾覆危险性评价方法 |
CN108763716B (zh) * | 2018-05-22 | 2021-06-08 | 西南交通大学 | 一种基于振型叠加法的过桥列车倾覆危险性评价方法 |
CN109459129A (zh) * | 2018-09-27 | 2019-03-12 | 成都理工大学 | 一种高效求解高速列车过桥振动随机特性的分析方法 |
CN109556812A (zh) * | 2018-12-13 | 2019-04-02 | 佛山科学技术学院 | 一种行人与人行桥横向相互作用动力模拟装置 |
CN110376452A (zh) * | 2019-06-14 | 2019-10-25 | 中国科学院西安光学精密机械研究所 | 基于耦合机电分析的压电陶瓷促动器电噪声指标确定方法 |
CN111651840A (zh) * | 2020-04-09 | 2020-09-11 | 北京交通大学 | 一种基于深度学习技术的轨道板上拱状态检测方法 |
CN111651840B (zh) * | 2020-04-09 | 2024-04-02 | 北京交通大学 | 一种基于深度学习技术的轨道板上拱状态检测方法 |
CN113432815A (zh) * | 2021-01-26 | 2021-09-24 | 重庆大学 | 一种基于测量车振动响应的桥面响应重构方法 |
CN113065187B (zh) * | 2021-04-06 | 2023-05-02 | 华东交通大学 | 一种基于解析映射关系的震后高铁连续梁桥墩不均匀变形阈值计算方法 |
CN113065187A (zh) * | 2021-04-06 | 2021-07-02 | 华东交通大学 | 一种基于解析映射关系的震后高铁连续梁桥墩不均匀变形阈值计算方法 |
CN114239347A (zh) * | 2021-11-30 | 2022-03-25 | 重庆大学 | 一种基于运营列车振动信号实时量测铁路轨道支撑刚度的高效方法 |
CN114186362A (zh) * | 2021-11-30 | 2022-03-15 | 重庆大学 | 一种基于运营列车振动响应预测轨道不平顺及其传感器优化配置的分析方法 |
CN114239347B (zh) * | 2021-11-30 | 2024-06-11 | 重庆大学 | 一种基于运营列车振动信号实时量测铁路轨道支撑刚度的高效方法 |
CN114707203A (zh) * | 2022-03-07 | 2022-07-05 | 东南大学 | 一种桥梁减振降噪快速化预测评估方法 |
CN114997020A (zh) * | 2022-05-20 | 2022-09-02 | 重庆大学 | 一种基于多自由度车辆动态响应的接触点反演算法 |
CN114997020B (zh) * | 2022-05-20 | 2024-05-28 | 重庆大学 | 一种基于多自由度车辆动态响应的接触点反演算法 |
CN115808280A (zh) * | 2022-11-16 | 2023-03-17 | 中国矿业大学 | 面向轨道减振评估的双轮对-构架多模式激励设备及方法 |
CN115808280B (zh) * | 2022-11-16 | 2023-08-22 | 中国矿业大学 | 面向轨道减振评估的双轮对-构架多模式激励设备及方法 |
CN117172136A (zh) * | 2023-11-02 | 2023-12-05 | 成都理工大学 | 基于ssa-lstm算法的风车桥***振动响应预测方法 |
CN117172136B (zh) * | 2023-11-02 | 2024-01-05 | 成都理工大学 | 基于ssa-lstm算法的风车桥***振动响应预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106650067B (zh) | 2019-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106650067B (zh) | 一种输出中高频响应的车-轨-桥耦合振动频域分析方法 | |
Lei et al. | Dynamic response analyses of vehicle and track coupled system on track transition of conventional high speed railway | |
Guo et al. | Integral model for train-track-bridge interaction on the Sesia viaduct: Dynamic simulation and critical assessment | |
Yang et al. | Effect of track irregularity on the dynamic response of a slab track under a high-speed train based on the composite track element method | |
Lei et al. | Influence of track stiffness distribution on vehicle and track interactions in track transition | |
CN108256278A (zh) | 一种悬挂式单轨车辆耦合动力学仿真***及其仿真方法 | |
Yang et al. | A novel dynamics model for railway ballastless track with medium-thick slabs | |
Spiryagin et al. | Control system for maximum use of adhesive forces of a railway vehicle in a tractive mode | |
Zakeri et al. | Investigation on dynamic behavior of railway track in transition zone | |
Zhai et al. | Train/track/bridge dynamic interactions: simulation and applications | |
CN106951668A (zh) | 一种铁路大跨度斜拉桥的刚度控制方法及斜拉桥 | |
Zhang et al. | Train-induced vibration and noise radiation of a prestressed concrete box-girder | |
Podworna et al. | Vertical vibrations of composite bridge/track structure/high-speed train systems. Part 2: Physical and mathematical modelling | |
Shi et al. | A theoretical investigation on influences of slab tracks on vertical dynamic responses of railway viaducts | |
Jiang et al. | Running safety and seismic optimization of a fault-crossing simply-supported girder bridge for high-speed railways based on a train-track-bridge coupling system | |
Xiang et al. | Analysis theory of spatial vibration of high-speed train and slab track system | |
Yang et al. | Modeling slab track for vehicle–track-coupled dynamics analysis using spline function method | |
Li et al. | Resonance analysis of cantilever plates subjected to moving forces by a semi-analytical method | |
Li et al. | Influences of soil-structure interaction on coupled vibration of train-bridge system: theoretical and experimental study | |
Xia et al. | Application of train-bridge-interaction analysis to bridge design of high-speed railways in China | |
Yi | Dynamic analysis of high-speed railway alignment: theory and practice | |
CN208673333U (zh) | 一种悬挂式单轨车辆耦合动力学仿真*** | |
Zhouzhou et al. | Research on vehicle-bridge vertical coupling dynamics of monorail based on multiple road excitations | |
Shi et al. | Interaction dynamic response of a high-speed train moving over curved bridges with deficient or surplus superelevation | |
Xia et al. | Dynamic responses of railway suspension bridges under moving trains |
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 |