CN105512427A - 考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法 - Google Patents
考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法 Download PDFInfo
- Publication number
- CN105512427A CN105512427A CN201511018869.XA CN201511018869A CN105512427A CN 105512427 A CN105512427 A CN 105512427A CN 201511018869 A CN201511018869 A CN 201511018869A CN 105512427 A CN105512427 A CN 105512427A
- Authority
- CN
- China
- Prior art keywords
- pantograph
- osculatory
- contact
- catenary
- speed
- 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
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法,在弓网发生离线时,通过弓网离线结束瞬间的弓头速度与接触线模态速度计算再接触时刻的速度增量,从而可以考虑考虑弓网再接触时受电弓弓头对接触线的动量冲击作用,使弓网动力学仿真更加精确,避免了传统仿真结果过于保守的缺点,该发明弥补了传统方法的不足,具有考虑因素更加全面,更安全的优点。
Description
技术领域
本发明涉及电动车辆的电源线路或沿路轨的仿真技术,尤其是一种考虑高速铁路受电弓与接触网离线再接触时的动量冲击的弓网动态仿真方法。
背景技术
高速铁路接触网与受电弓动力学行为十分复杂,数值仿真计算是在研究弓网动态行为、设计其关键参数的重要途径,目前,有大量的文章都采用弓网数值仿真模型研究弓网的动态性能,为工程实际提供了诸多有益借鉴。比如期刊《VehicleSystemDynamics》2015年的第53卷第10期1455-1479页YSong,ZLiu,HWang,XLu,JZhang等人的论文《Nonlinearmodellingofhigh-speedcatenarybasedonanalyticalexpressionsofcableandtrusselements》提出了高速铁路接触网的非线性动态仿真建模方法,并研究了接触网不同参数对弓网接触力的影响。然而目前所有的弓网数值仿真方法中,并未考虑受电弓与接触网脱离后,再接触的动量冲击,这种忽略会影响仿真计算的精度,导致仿真计算结果的保守,不利于弓网系统的后期设计。本发明的目的在于提出一种考虑弓网再接触动量冲击的动力学仿真方法,采用弓网接触前瞬间时刻的弓头速度与接触线各阶模态速度,确定接触线在受电弓再接触时刻由动量冲击引起的速度增量,并将其带入提高仿真的精确性和后续设计的安全性。
发明内容
本发明的目的在于提出一种考虑弓网再接触动量冲击的动力学仿真方法,采用弓网接触前瞬间时刻的弓头速度与接触线各阶模态速度,确定接触线在受电弓再接触时刻由动量冲击引起的速度增量,并将其带入提高仿真的精确性和后续设计的安全性。
本发明实现发明目的的技术手段为:
考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法,在考虑受电弓脱离接触线后与接触线再接触瞬间弓头对接触线的动量冲击影响下进行高铁接触网与受电弓动力学数值仿真,包含以下具体步骤:
1)、采用传统模态分解法与多体动力学方法分别构建接触网和受电弓的动力学模型,建立弓网耦合运动方程,
其中,M、C、K分别为弓网整体质量、阻尼和刚度矩阵,X(t)分别为弓网整体加速度、速度、位移矩阵;并根据此方程进行弓网动力学仿真计算,F(t)为外界击扰;
2)、在1)弓网动力学数值仿真计算中,若弓网离线发生,将受电弓与接触线再接触瞬间的时刻t2分解为受电弓与接触线离线结束时刻和弓头与接触线刚接触瞬间的时刻
3)、有1)和2)的结果,确定在时刻受电弓弓头速度和接触线各阶模态速度其中,n为接触线的模态阶数;
4)、通过以下公式计算接触线在受电弓再接触冲击下的各阶模态速度增量:
Bn(t)为接触线第n阶模态的广义位移,y1(t)为受电弓弓头位移,ψn(t)为接触线第n阶模态的振形函数,v为列车行驶速度,ρb为接触线的线密度,L为接触线的锚段长度,m1为受电弓弓头质量;
5)、确定接触线在接触点vt2处,时刻的速度增量:
6)、将动量冲击产生的附加速度带入弓网动力学方程中,进行后续仿真。
这样,在弓网动态仿真过程中,一旦弓网发生离线与再接触,可以通过接触前时刻的受电弓弓头速度与接触线各阶模态速度计算出再接触时刻接触线由于弓头冲击产生的速度增量。为更为清楚地描述本方法的,对上述技术手段的获取基础作以下展开叙述:
1、建立高速铁路弓网动力学模型,如图1所示,分别采用模态分析法和多体动力学方法建立高速弓网动力学模型。承力索和接触线采用欧拉伯努利梁单元,吊弦、定位器以及支撑杆等悬挂装置看作集中质量、刚度点。承力索的运动方程如公式(1)所示。接触线方程如公式(2)所示。
其中,Ta和Tb分别为承力索和接触线的张力。ρa和ρb分别为承力索和接触线的线密度。EIa和EIb分别为承力索和接触线的弯曲刚度。Ca和Cb分别为承力索和接触线的阻尼,L为锚段长度。δ为狄力克雷函数符号。wa(x,t)和wb(x,t)分别为承力索和接触线位移,和分别为第i根吊弦对承力索和接触线的作用力,其作用点为 是第j根支撑杆对承力索的作用力,其作用点为为第j根定位器对接触线的作用力,其作用点为 和为第i根吊弦的刚度,和分别是第j根支撑杆的质量和刚度。是第j根定位器的质量。nd、nt和nr分别是吊弦、支撑杆、定位器的数量。
通过模态分解方法,承力索和接触线的位移可以写成以下模态叠加形式。
其中,An(t)和Bn(t)分别为第n阶模态的承力索和接触线广义位移。ψn(x)和ψn(x)分别是承力索和接触线的第n阶模态振形函数。通过将式(3)带入式(1-2)则可以构建接触网的动力学方程。受电弓的动力学运动方程采用传统的三元受电弓建模方法,如式(4)所示,
其中,m1、m2、m3分别为受电弓弓头、上下框架的归算质量,c1、c2、c3分别为对应的阻尼,k1、k2、k3分别为对应的刚度。f0为静态抬升力,fc(vt,t)为弓网接触力,y1(t)、y2(t)、y3(t)分别为弓头、上下框架位移。v为列车行进速度。弓网之间接触采用拉格朗日乘子法实现,从而式(1-5)构建了高速铁路接触网-受电弓动力学方程,式(1-5)可以写成矩阵形式为:
其中,M、C、K分别为弓网整体质量、阻尼和刚度矩阵,X(t)分别为弓网整体加速度、速度、位移矩阵。并根据此方程进行弓网动力学仿真计算,F(t)为外界击扰;
2、对式(1-4)的求解采用Newmark法实现,当wb(vt,t)-y1(t)>0时,认为接触线与受电弓分离。捕捉弓头与接触线再次接触的时刻t2,即下一个满足wb(vt2,t2)=y1(t2)关系的时刻。
3、通过求解弓网动力学方程(1-4),确定在时刻受电弓弓头速度和接触线各阶模态速度其中,n为接触线的模态阶数;
4、通过以下公式计算接触线在受电弓再接触冲击下的各阶模态速度增量:
5、确定接触线在接触点vt2处,时刻的速度增量:
6、将动量冲击产生的附加速度带入弓网动力学方程(1-4)中,进行后续仿真。
采用本发明的方法,在弓网动态仿真过程中,一旦弓网发生离线与再接触,可以通过接触前时刻的受电弓弓头速度与接触线各阶模态速度计算出再接触时刻接触线由于弓头冲击产生的速度增量。与传统的仿真方法相比,本方法考虑了弓网离线再接触时的动量冲击,提高了仿真计算的准确性。
附图说明
图1为弓网模型示意图
图2为弓网接触力仿真结果图
图3为接触线接触点位移图
图4为接触线接触点速度图
具体实施方式
下面结合附图对本发明的实施例做详细说明:本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施过程,但本发明的保护范围不限于下述的实施例。
本实例以设计列车运行速度为350Km/h的京津铁路接触网和DSA380型受电弓为例,接触网与受电弓参数均源于文献[Dynamicperformanceofpantograph/overheadlineinteractionfor4spanoverlaps–TPS/OCSportion.SIEMENS,2006:4-21.]。采用本发明提出的仿真方法对高速铁路弓网的动态响应进行仿真计算,并于传统未考虑再接触动量冲击方法的计算结果进行对比。
首先按照步骤1建立接触网-受电弓动力学模型,并捕捉受电弓离线时刻,
本算例中,受电弓与接触线发生离线的时刻为0.437s(42.5m处),二者再次接触的时刻为0.448s(43.6m)。随后按照2、3、4、5等步骤确定再接触时刻的由于动量冲击所产生的附加速度,并按照步骤6带入到动力学方程(1-4)中继续后续仿真。接触力的仿真结果如图2所示,接触线接触点位移和速度仿真结果如图3、图4所示,其中红色线条为本方法计算结果,黑色线条为未考虑弓网再接触动量冲击传统方法计算结果。从两种方法的计算结果可以看到,由于动量冲击的作用,弓网接触力、接触线位移、速度在接触点均有一个明显的突升。而传统方法并不能够体现这种突升,对弓网系统的安全运行评估过于保守。
Claims (1)
1.考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法,在考虑受电弓脱离接触线后与接触线再接触瞬间弓头对接触线的动量冲击影响下进行高铁接触网与受电弓动力学数值仿真,包含以下具体步骤:
1)、采用传统模态分解法与多体动力学方法分别构建接触网和受电弓的动力学模型,建立弓网耦合运动方程,
其中,M、C、K分别为弓网整体质量、阻尼和刚度矩阵,X(t)分别为弓网整体加速度、速度、位移矩阵;并根据此方程进行弓网动力学仿真计算,F(t)为外界击扰;
2)、在1)弓网动力学数值仿真计算中,若弓网离线发生,将受电弓与接触线再接触瞬间的时刻t2分解为受电弓与接触线离线结束时刻和弓头与接触线刚接触瞬间的时刻
3)、有1)和2)的结果,确定在时刻受电弓弓头速度和接触线各阶模态速度其中,n为接触线的模态阶数;
4)、通过以下公式计算接触线在受电弓再接触冲击下的各阶模态速度增量:
Bn(t)为接触线第n阶模态的广义位移,y1(t)为受电弓弓头位移,ψn(t)为接触线第n阶模态的振形函数,v为列车行驶速度,ρb为接触线的线密度,L为接触线的锚段长度,m1为受电弓弓头质量;
5)、确定接触线在接触点vt2处,时刻的速度增量:
6)、将动量冲击产生的附加速度带入弓网动力学方程中,进行后续仿真。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511018869.XA CN105512427B (zh) | 2015-12-29 | 2015-12-29 | 考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511018869.XA CN105512427B (zh) | 2015-12-29 | 2015-12-29 | 考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105512427A true CN105512427A (zh) | 2016-04-20 |
CN105512427B CN105512427B (zh) | 2018-09-28 |
Family
ID=55720407
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511018869.XA Expired - Fee Related CN105512427B (zh) | 2015-12-29 | 2015-12-29 | 考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105512427B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107621319A (zh) * | 2017-09-29 | 2018-01-23 | 西南交通大学 | 一种弓网接触力测量方法及其测量装置 |
CN110362786A (zh) * | 2019-06-14 | 2019-10-22 | 西南交通大学 | 一种高速铁路弓网耦合系统动力学响应分析的解耦建模方法 |
CN110457822A (zh) * | 2019-08-13 | 2019-11-15 | 西南交通大学 | 一种判定弓网接触电弧发生的接触力阈值模型构建方法 |
CN113255195A (zh) * | 2021-07-02 | 2021-08-13 | 成都国佳电气工程有限公司 | 一种考虑实测分布式阻尼的接触网有限元建模方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010122102A (ja) * | 2008-11-20 | 2010-06-03 | East Japan Railway Co | 紫外線検出装置及び離線検出記録計 |
CN104361196A (zh) * | 2014-09-22 | 2015-02-18 | 西南交通大学 | 一种计及列车速度的弓网离线电弧数学模型 |
CN104597354A (zh) * | 2015-01-27 | 2015-05-06 | 南车青岛四方机车车辆股份有限公司 | 基于列车供电系统的弓网电弧的检测方法及系统 |
CN104765916A (zh) * | 2015-03-31 | 2015-07-08 | 西南交通大学 | 一种高速列车动力学性能参数优化方法 |
-
2015
- 2015-12-29 CN CN201511018869.XA patent/CN105512427B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010122102A (ja) * | 2008-11-20 | 2010-06-03 | East Japan Railway Co | 紫外線検出装置及び離線検出記録計 |
CN104361196A (zh) * | 2014-09-22 | 2015-02-18 | 西南交通大学 | 一种计及列车速度的弓网离线电弧数学模型 |
CN104597354A (zh) * | 2015-01-27 | 2015-05-06 | 南车青岛四方机车车辆股份有限公司 | 基于列车供电系统的弓网电弧的检测方法及系统 |
CN104765916A (zh) * | 2015-03-31 | 2015-07-08 | 西南交通大学 | 一种高速列车动力学性能参数优化方法 |
Non-Patent Citations (2)
Title |
---|
JORGE AMBROSIO等: ""Optimization of high-speed railway pantographs for improving pantograph-catenary contact"", 《THEORETICAL & APPLIED MECHANICS LETTERS》 * |
姜静等: ""考虑受电弓非线性的弓网动态特性仿真研究"", 《计算机仿真》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107621319A (zh) * | 2017-09-29 | 2018-01-23 | 西南交通大学 | 一种弓网接触力测量方法及其测量装置 |
CN107621319B (zh) * | 2017-09-29 | 2023-10-13 | 西南交通大学 | 一种弓网接触力测量方法及其测量装置 |
CN110362786A (zh) * | 2019-06-14 | 2019-10-22 | 西南交通大学 | 一种高速铁路弓网耦合系统动力学响应分析的解耦建模方法 |
CN110362786B (zh) * | 2019-06-14 | 2023-03-14 | 西南交通大学 | 一种高速铁路弓网耦合系统动力学响应分析的解耦建模方法 |
CN110457822A (zh) * | 2019-08-13 | 2019-11-15 | 西南交通大学 | 一种判定弓网接触电弧发生的接触力阈值模型构建方法 |
CN110457822B (zh) * | 2019-08-13 | 2022-04-29 | 西南交通大学 | 一种判定弓网接触电弧发生的接触力阈值模型构建方法 |
CN113255195A (zh) * | 2021-07-02 | 2021-08-13 | 成都国佳电气工程有限公司 | 一种考虑实测分布式阻尼的接触网有限元建模方法 |
CN113255195B (zh) * | 2021-07-02 | 2021-10-08 | 成都国佳电气工程有限公司 | 一种考虑实测分布式阻尼的接触网有限元建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105512427B (zh) | 2018-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105512427A (zh) | 考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法 | |
CN105631128A (zh) | 一种高速铁路弓-网-车-轨垂向耦合大系统动力学建模仿真方法 | |
CN106873369B (zh) | 针对列车输入受限和执行器故障的自适应容错控制方法 | |
CN102339350B (zh) | 一种基于参数化设计的整车碰撞仿真方法 | |
CN106503362A (zh) | 一种高速受电弓多目标鲁棒h∞控制器设计方法 | |
CN107958098A (zh) | 一种基于频谱分析的高速列车倾覆危险性评价方法 | |
CN107202584A (zh) | 一种行星精确着陆抗扰制导方法 | |
CN104217253A (zh) | 台风天气下配电线路荷载可靠度分析方法 | |
CN107102542A (zh) | 一种列车自动运行的鲁棒自适应非奇异终端滑模控制方法 | |
CN101733750B (zh) | 基于动态响应谱的关节含间隙机器人轨迹规划方法 | |
CN109375510A (zh) | 一种用于高速列车的自适应滑模容错控制方法 | |
CN105607485A (zh) | 基于故障特征模型的挠性充液卫星姿态自适应容错控制方法 | |
CN105487384A (zh) | 一种基于事件触发机制的整车悬架控制系统及其设计方法 | |
CN107817683A (zh) | 一种动车组精准停车方法 | |
CN103010268A (zh) | 一种动车邻域子系统的动力分布式优化调度方法 | |
CN104155572A (zh) | 一种同塔双回直流输电线路故障选线方法 | |
CN103123668B (zh) | 一种基于混合单元法的空间绳系机器人系统的仿真方法 | |
Peng et al. | A new macro model of traffic flow by incorporating both timid and aggressive driving behaviors | |
CN102890737B (zh) | 一种着陆冲击下的机械结构累积损伤计算方法 | |
CN103678781A (zh) | 一种阀门动部件运动与碰撞仿真分析方法 | |
Ma et al. | Solution algorithm and experimental verification on behavior of crossing frame impacted by conductor breaking | |
Zhengbao et al. | The response of collision speed caused by the large bus to new flexible barrier' | |
Zhang et al. | Effect of current vehicle’s interruption on traffic stability in cooperative car-following theory | |
CN110285979A (zh) | 一种侧面柱碰撞台车试验方法及其参数提取方法 | |
CN110705164B (zh) | 冰—船相互作用模拟方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180928 Termination date: 20181229 |
|
CF01 | Termination of patent right due to non-payment of annual fee |