CN105512427B - 考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法 - Google Patents
考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法 Download PDFInfo
- Publication number
- CN105512427B CN105512427B CN201511018869.XA CN201511018869A CN105512427B CN 105512427 B CN105512427 B CN 105512427B CN 201511018869 A CN201511018869 A CN 201511018869A CN 105512427 B CN105512427 B CN 105512427B
- Authority
- CN
- China
- Prior art keywords
- contact
- bow net
- contact line
- bow
- pantograph
- 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
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000005094 computer simulation Methods 0.000 title claims abstract description 10
- 238000004088 simulation Methods 0.000 claims abstract description 22
- 241001125292 Balaena mysticetus Species 0.000 claims abstract description 11
- 238000006073 displacement reaction Methods 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 7
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 4
- 238000013016 damping Methods 0.000 claims description 4
- 230000033001 locomotion Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000007493 shaping process Methods 0.000 claims description 3
- 229910052742 iron Inorganic materials 0.000 claims description 2
- 230000005611 electricity Effects 0.000 claims 1
- 238000007796 conventional method Methods 0.000 abstract description 3
- 230000007812 deficiency Effects 0.000 abstract 1
- 238000009527 percussion Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005021 gait Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 230000005428 wave function Effects 0.000 description 1
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
技术领域
本发明涉及电动车辆的电源线路或沿路轨的仿真技术,尤其是一种考虑高速铁路受电弓与接触网离线再接触时的动量冲击的弓网动态仿真方法。
背景技术
高速铁路接触网与受电弓动力学行为十分复杂，数值仿真计算是在研究弓网动态行为、设计其关键参数的重要途径，目前，有大量的文章都采用弓网数值仿真模型研究弓网的动态性能，为工程实际提供了诸多有益借鉴。比如期刊《Vehicle System Dynamics》2015年的第53卷第10期1455-1479页Y Song,Z Liu,H Wang,X Lu,J Zhang等人的论文《Nonlinear modelling of high-speed catenary based on analytical expressionsof cable and truss elements》提出了高速铁路接触网的非线性动态仿真建模方法，并研究了接触网不同参数对弓网接触力的影响。然而目前所有的弓网数值仿真方法中，并未考虑受电弓与接触网脱离后，再接触的动量冲击，这种忽略会影响仿真计算的精度，导致仿真计算结果的保守，不利于弓网系统的后期设计。本发明的目的在于提出一种考虑弓网再接触动量冲击的动力学仿真方法，采用弓网接触前瞬间时刻的弓头速度与接触线各阶模态速度，确定接触线在受电弓再接触时刻由动量冲击引起的速度增量，并将其带入提高仿真的精确性和后续设计的安全性。
发明内容
本发明的目的在于提出一种考虑弓网再接触动量冲击的动力学仿真方法，采用弓网接触前瞬间时刻的弓头速度与接触线各阶模态速度，确定接触线在受电弓再接触时刻由动量冲击引起的速度增量，并将其带入提高仿真的精确性和后续设计的安全性。
本发明实现发明目的的技术手段为：
考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法，在考虑受电弓脱离接触线后与接触线再接触瞬间弓头对接触线的动量冲击影响下进行高铁接触网与受电弓动力学数值仿真，包含以下具体步骤：
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分别是吊弦、支撑杆、定位器的数量。
通过模态分解方法，承力索和接触线的位移可以写成以下模态叠加形式。
(3)
其中，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型受电弓为例，接触网与受电弓参数均源于文献[Dynamic performance of pantograph/overheadline interaction for 4span overlaps–TPS/OCS portion.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 CN105512427A (zh) | 2016-04-20 |
| CN105512427B true 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) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107621319B (zh) * | 2017-09-29 | 2023-10-13 | 西南交通大学 | 一种弓网接触力测量方法及其测量装置 |
| CN110362786B (zh) * | 2019-06-14 | 2023-03-14 | 西南交通大学 | 一种高速铁路弓网耦合系统动力学响应分析的解耦建模方法 |
| CN110457822B (zh) * | 2019-08-13 | 2022-04-29 | 西南交通大学 | 一种判定弓网接触电弧发生的接触力阈值模型构建方法 |
| CN113255195B (zh) * | 2021-07-02 | 2021-10-08 | 成都国佳电气工程有限公司 | 一种考虑实测分布式阻尼的接触网有限元建模方法 |
| CN118544903A (zh) * | 2024-07-30 | 2024-08-27 | 四川工程职业技术大学 | 一种抑制接触线振动装置及方法 |
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 |
|---|
| "Optimization of high-speed railway pantographs for improving pantograph-catenary contact";Jorge Ambrosio等;《THEORETICAL & APPLIED MECHANICS LETTERS》;20130110;第3卷(第1期);第013006-1-013006-7页 * |
| "考虑受电弓非线性的弓网动态特性仿真研究";姜静等;《计算机仿真》;20150215;第32卷(第2期);第170-174页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105512427A (zh) | 2016-04-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105512427B (zh) | 考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法 | |
| Song et al. | Active control of contact force for high-speed railway pantograph-catenary based on multi-body pantograph model | |
| CN106503362A (zh) | 一种高速受电弓多目标鲁棒h∞控制器设计方法 | |
| CN107958098B (zh) | 一种基于频谱分析的高速列车倾覆危险性评价方法 | |
| Song et al. | Sliding Mode Control with PD Sliding Surface for High‐Speed Railway Pantograph‐Catenary Contact Force under Strong Stochastic Wind Field | |
| Gu et al. | Direct identification of flutter derivatives and aerodynamic admittances of bridge decks | |
| CN104589359B (zh) | 一种基于振动观测器的柔性机械臂振动控制方法 | |
| CN103676646B (zh) | 一种具有随机发生的不确定性和分布式传感器时滞的网络化控制系统的状态估计方法 | |
| CN110362786A (zh) | 一种高速铁路弓网耦合系统动力学响应分析的解耦建模方法 | |
| CN104299034A (zh) | 基于bp神经网络的三芯电缆导体暂态温度计算方法 | |
| CN104809292A (zh) | 一种高速列车非线性动力学模型参数的在线辨识方法 | |
| CN102393184A (zh) | 电气化铁路接触网风致振动及风偏的动态计算方法 | |
| CN102866635A (zh) | 基于等价模型的高超声速飞行器离散神经网络自适应控制方法 | |
| CN105224763B (zh) | 一种列车的隧道附加空气阻力迭代学习辨识方法 | |
| CN108021747A (zh) | 一种高铁弓网动态行为动网格不平衡力消除的仿真方法 | |
| CN102521651A (zh) | 基于narx神经网络的弓网接触力预测方法 | |
| Ma et al. | Case study of three-dimensional aeroelastic effect on critical flutter wind speed of long-span bridges | |
| Li et al. | Active control strategy for the running attitude of high-speed train under strong crosswind condition | |
| RU2013104504A (ru) | Система и способ тестирования показателя работы паровой турбины | |
| CN115481488B (zh) | 基于机器学习的方锥式吸能结构耐撞性能多目标优化方法 | |
| CN110966144A (zh) | 一种基于bls的风电机组智能恒功率控制方法 | |
| CN102904255B (zh) | 基于解析灵敏度的电力系统低压切负荷优化控制方法 | |
| CN110095984A (zh) | 一种基于状态估计的高速铁路受电弓多目标鲁棒控制方法 | |
| CN103440406B (zh) | 一种离线迭代质点弹道计算方法 | |
| Arsene et al. | Analysis of the wind influence on the aerodynamic drag in the case of a certain emplacement of the pantograph on the electric rail vehicles |
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 |