CN113255195B - 一种考虑实测分布式阻尼的接触网有限元建模方法 - Google Patents
一种考虑实测分布式阻尼的接触网有限元建模方法 Download PDFInfo
- Publication number
- CN113255195B CN113255195B CN202110747666.3A CN202110747666A CN113255195B CN 113255195 B CN113255195 B CN 113255195B CN 202110747666 A CN202110747666 A CN 202110747666A CN 113255195 B CN113255195 B CN 113255195B
- Authority
- CN
- China
- Prior art keywords
- damping
- matrix
- contact
- carrier cable
- unit
- 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.)
- Active
Links
Images
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
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)
- Current-Collector Devices For Electrically Propelled Vehicles (AREA)
Abstract
本发明公开了一种考虑实测分布式阻尼的接触网有限元建模方法,具体为:通过现场实测获得接触网一跨内关键点的阻尼常数,通过spline样条函数拟合分别对接触线和承力索的阻尼常数进行拟合;采用有限元方法构建接触网模型,通过插值获得接触线和承力索上每一个节点的阻尼常数,生成单元阻尼矩阵;最后根据有限元方法生成接触网总体阻尼矩阵,形成接触网运动方程,完成考虑实测分布式阻尼的接触网有限元建模。本发明克服了传统接触网建模中的恒定阻尼比假设,使弓网动力学仿真更加精确,避免了传统仿真结果过于激进的缺点,具有考虑因素更加全面,计算结果更加符合实际、更安全的优点。
Description
技术领域
本发明属于接触网建模技术领域,尤其涉及一种考虑实测分布式阻尼的接触网有限元建模方法。
背景技术
高速铁路接触网是受电弓通行的滑倒,其结构十分复杂。由于线路试验高额的成本,数值仿真是研究受电弓-接触网动态交互行为主流方法。目前,有大量的文章都采用受电弓-接触网数值仿真模型研究弓网的动态性能,为接触网的设计建造提供了诸多有益借鉴。比如期刊《Vehicle System Dynamics》2015年的第53卷第3期412-435页BruniStefano,Ambrosio Jorge,Carnicero Alberto,Cho Yong Hyeon,Finner Lars,IkedaMitsuru,Kwon Sam Young,Massat Jean Pierre,Stichel Sebastian,Tur Manuel,ZhangWeihua等人的论文《The results of the pantograph-catenary interactionbenchmark》对比了目前主流的接触网仿真软件,并制定了受电弓-接触网仿真结果验证标准。然而目前所有的接触网模型中,整个接触网系统均假定为阻尼比恒定。但是,实际观测发现,接触网各个点的阻尼特性均不相同。恒定阻尼假设可能会引起一定的仿真误差。
发明内容
为解决上述问题,本发明提供一种考虑实测分布式阻尼的接触网有限元建模方法。
本发明的一种考虑实测分布式阻尼的接触网有限元建模方法具体为:
其中,表示接触线质量矩阵对应的阻尼比系数向量,表示承力索质量矩阵对应的阻尼比系数向量,表示接触线刚度矩阵对应的阻尼比系数向量,表示承力索刚度矩阵对应的阻尼比系数向量;关键点包括:吊弦与接触线连接点、吊弦与承力索连接点、定位器点、承力索座点。
步骤2:通过spline样条函数拟合分别对接触线和承力索的阻尼常数进行拟合。
步骤3:采用有限元方法构建接触网模型,通过插值获得接触线和承力索上每一个节点的阻尼常数,生成单元阻尼矩阵。
将接触线和承力索划分为N element个单元,形成总体质量矩阵M和总体刚度矩阵K;根据四条阻尼常数拟合曲线,通过插值获得接触线和承力索上每个单元的阻尼常数,计算单元阻尼矩阵。
步骤4:根据有限元方法生成接触网总体阻尼矩阵,形成接触网运动方程,完成考虑实测分布式阻尼的接触网有限元建模。
根据有限元方法,将单元阻尼矩阵组装为整体阻尼矩阵C,结合总体质量矩阵M和总体刚度矩阵K,获得接触网运动方程:
本发明的有益技术效果为:
本发明在接触网有限元建模中,根据实测接触网关键点的阻尼比生成接触网总体阻尼矩阵,能够在数值仿真中考虑接触网的实测分布式阻尼。与传统的仿真方法相比,本方法舍弃了接触网系统恒定阻尼比假设,考虑了接触网阻尼比沿跨距方向的变化,提高了仿真计算的准确性。在后续的受电弓-接触网交互仿真中能够得到更加符合实际的结果。
附图说明
图1为实测一跨接触线阻尼常数(对应质量)。
图2为实测一跨承力索阻尼常数(对应质量)。
图3为实测一跨接触线阻尼常数(对应刚度)。
图4为实测一跨承力索阻尼常数(对应刚度)。
图5为实测一跨接触线阻尼常数曲线拟合结果(对应质量)。
图6为实测一跨承力索阻尼常数曲线拟合结果(对应质量)。
图7为实测一跨接触线阻尼常数曲线拟合结果(对应刚度)。
图8为实测一跨承力索阻尼常数曲线拟合结果(对应刚度)。
图9为弓网接触力仿真结果对比。
图10为弓网接触力仿真结果对比细节示意图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细说明。
本发明的一种考虑实测分布式阻尼的接触网有限元建模方法具体为:
其中,表示接触线质量矩阵对应的阻尼比系数向量,表示承力索质量矩阵对应的阻尼比系数向量,表示接触线刚度矩阵对应的阻尼比系数向量,表示承力索刚度矩阵对应的阻尼比系数向量;关键点包括:吊弦与接触线连接点、吊弦与承力索连接点、定位器点、承力索座点。
步骤2:通过spline样条函数拟合分别对接触线和承力索的阻尼常数进行拟合。
步骤3:采用有限元方法构建接触网模型,通过插值获得接触线和承力索上每一个节点的阻尼常数,生成单元阻尼矩阵。
将接触线和承力索划分为N element个单元,将每个单元矩阵按如下形式进行组装,形成总体质量矩阵M和总体刚度矩阵K。
根据以下流程生成每个单元的单元阻尼矩阵。
对N element个单元逐个进行判断,即令i = 1→N element
如果单元i在接触线上:
(2)根据拟合曲线计算:
如果单元i在承力索上,
(2)根据拟合曲线计算:
步骤4:根据有限元方法生成接触网总体阻尼矩阵,形成接触网运动方程,完成考虑实测分布式阻尼的接触网有限元建模。
根据有限元方法,将单元阻尼矩阵组装为整体阻尼矩阵C,结合总体质量矩阵M和总体刚度矩阵K,获得接触网运动方程:
实施例
以设计某普速铁路接触网为例,接触网与受电弓参数均源于文献[Identification of short-wavelength contact wire irregularities inelectrified railway pantograph – catenary system, Mech. Mach. Theory. 162(2021) 104338]。采用本发明提出的考虑实测分布式阻尼的接触网模型进行受电弓-接触网动态响应计算,并与传统恒定阻尼比接触网模型计算结果进行对比。
首先按照步骤1获得实测接触网关键点阻尼常数,如图1-图4所示。由于实测阻尼常数呈现出随机性,每个位置的实测阻尼常数用箱线图进行表示,中心圆点为每个位置阻尼常数的均值,上下边缘分别表示阻尼常数的离散程度。然后根据步骤2对接触网阻尼常数进行拟合,分别获得接触线和承力索a和b的拟合曲线,如图5-图8所示。随后按照步骤3、4构建接触网动力学模型。最后,引入归算质量受电弓模型求解受电弓-接触网动态交互响应。弓网接触力的计算结果如图9-图10所示,其中实线为本方法计算结果,虚线为采用恒定阻尼比假设的传统模型计算结果。从两种方法的计算结果可以看到,传统模型计算出的离线间距更小,按照传统模型计算结果进行设计,设计方案可能会更加激进。
Claims (1)
1.一种考虑实测分布式阻尼的接触网有限元建模方法,其特征在于,包括以下步骤:
关键点包括:吊弦与接触线连接点、吊弦与承力索连接点、定位器点、承力索座点;
步骤2:通过spline样条函数拟合分别对接触线和承力索的阻尼常数进行拟合:
步骤3:采用有限元方法构建接触网模型,通过插值获得接触线和承力索上每一个节点的阻尼常数,生成单元阻尼矩阵;
将接触线和承力索划分为N element个单元,将每个单元矩阵按如下形式进行组装,形成总体质量矩阵M和总体刚度矩阵K;
根据以下流程生成每个单元的单元阻尼矩阵;
对N element个单元逐个进行判断,即令i = 1→N element;
如果单元i在接触线上:
(2)根据拟合曲线计算:
如果单元i在承力索上:
(2)根据拟合曲线计算:
步骤4:根据有限元方法生成接触网总体阻尼矩阵,形成接触网运动方程,完成考虑实测分布式阻尼的接触网有限元建模;
根据有限元方法,将单元阻尼矩阵组装为整体阻尼矩阵C,结合总体质量矩阵M和总体刚度矩阵K,获得接触网运动方程:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110747666.3A CN113255195B (zh) | 2021-07-02 | 2021-07-02 | 一种考虑实测分布式阻尼的接触网有限元建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110747666.3A CN113255195B (zh) | 2021-07-02 | 2021-07-02 | 一种考虑实测分布式阻尼的接触网有限元建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113255195A CN113255195A (zh) | 2021-08-13 |
CN113255195B true CN113255195B (zh) | 2021-10-08 |
Family
ID=77190424
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110747666.3A Active CN113255195B (zh) | 2021-07-02 | 2021-07-02 | 一种考虑实测分布式阻尼的接触网有限元建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113255195B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114936487A (zh) * | 2022-04-29 | 2022-08-23 | 中铁第一勘察设计院集团有限公司 | 接触网定位器故障的定位方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011007279A (ja) * | 2009-06-26 | 2011-01-13 | Ihi Corp | 回転体の軸受支持装置 |
CN105512427A (zh) * | 2015-12-29 | 2016-04-20 | 西南交通大学 | 考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法 |
CN105631128A (zh) * | 2015-12-29 | 2016-06-01 | 西南交通大学 | 一种高速铁路弓-网-车-轨垂向耦合大系统动力学建模仿真方法 |
CN106503362A (zh) * | 2016-10-26 | 2017-03-15 | 西南交通大学 | 一种高速受电弓多目标鲁棒h∞控制器设计方法 |
CN108021747A (zh) * | 2017-11-30 | 2018-05-11 | 西南交通大学 | 一种高铁弓网动态行为动网格不平衡力消除的仿真方法 |
CN109766574A (zh) * | 2018-12-03 | 2019-05-17 | 西南交通大学 | 一种基于估计的h∞控制的双受电弓主动控制方法 |
CN111367173A (zh) * | 2020-03-06 | 2020-07-03 | 西南交通大学 | 一种基于状态估计的高速铁路受电弓鲁棒预测控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110850716B (zh) * | 2019-11-13 | 2022-11-22 | 兰州交通大学 | 基于接触网先验信息的受电弓神经网络滑模变结构主动控制方法 |
-
2021
- 2021-07-02 CN CN202110747666.3A patent/CN113255195B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011007279A (ja) * | 2009-06-26 | 2011-01-13 | Ihi Corp | 回転体の軸受支持装置 |
CN105512427A (zh) * | 2015-12-29 | 2016-04-20 | 西南交通大学 | 考虑弓网再接触动量冲击的高速铁路弓网动力学仿真方法 |
CN105631128A (zh) * | 2015-12-29 | 2016-06-01 | 西南交通大学 | 一种高速铁路弓-网-车-轨垂向耦合大系统动力学建模仿真方法 |
CN106503362A (zh) * | 2016-10-26 | 2017-03-15 | 西南交通大学 | 一种高速受电弓多目标鲁棒h∞控制器设计方法 |
CN108021747A (zh) * | 2017-11-30 | 2018-05-11 | 西南交通大学 | 一种高铁弓网动态行为动网格不平衡力消除的仿真方法 |
CN109766574A (zh) * | 2018-12-03 | 2019-05-17 | 西南交通大学 | 一种基于估计的h∞控制的双受电弓主动控制方法 |
CN111367173A (zh) * | 2020-03-06 | 2020-07-03 | 西南交通大学 | 一种基于状态估计的高速铁路受电弓鲁棒预测控制方法 |
Non-Patent Citations (4)
Title |
---|
Effective Measures to Improve Current Collection Quality for Double Pantographs and Catenary Based on Wave Propagation Analysis;Xu, Zhao 等;《IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY 》;20200630;第69卷(第6期);第6299-6309页 * |
有限元分析的开式预弯机轻量化设计;刘志刚;《中国优秀硕士学位论文全文数据库 (工程科技Ⅰ辑)》;20200215(第2期);第217-460页 * |
高速列车的关键力学问题;杨国伟 等;《力学进展》;20150731;第45卷;B022-588 * |
高速铁路接触网研究进展;刘志刚 等;《西南交通大学学报》;20160630;第51卷(第3期);第495-581页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113255195A (zh) | 2021-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107153894B (zh) | 一种风电场的预测风速校正方法及装置 | |
Limebeer et al. | Optimal control of a formula one car on a three-dimensional track—Part 2: Optimal control | |
KR100921575B1 (ko) | 지형해상도 차이를 고려한 대기 유동장 해석 방법 | |
CN100468422C (zh) | 利用交叉模型交叉模态的结构模型修正方法 | |
CN110334371A (zh) | 一种基于有限元模型的车-桥耦合系统振动计算方法 | |
WO2017031718A1 (zh) | 弹性物体变形运动的建模方法 | |
CN107016649A (zh) | 一种基于局部低秩张量估计的视觉数据补全方法 | |
CN113255195B (zh) | 一种考虑实测分布式阻尼的接触网有限元建模方法 | |
CN104679956A (zh) | 考虑动态特性的高速压力机底座可靠性稳健设计方法 | |
CN110362786B (zh) | 一种高速铁路弓网耦合系统动力学响应分析的解耦建模方法 | |
CN108139299A (zh) | 用于确定道路或轨道车辆的模态参数并间接表征道路或轨道轮廓的方法 | |
CN108536937B (zh) | 一种虚拟环境下多分支电缆网的物理建模和布线仿真方法 | |
CN104598753A (zh) | 一种基于Brakhage v方法的桥梁移动车辆荷载识别方法 | |
CN113505534A (zh) | 一种考虑需求响应的负荷预测方法 | |
CN107121257A (zh) | 一种垭口微地貌输电导线风致振动的风洞试验方法 | |
CN116956428A (zh) | 基于物理信息神经网络的移动荷载下桥梁响应计算方法 | |
CN116049989A (zh) | 转向架数字孪生模型构建方法、系统、电子设备及介质 | |
Branlard et al. | A digital twin solution for floating offshore wind turbines validated using a full-scale prototype | |
JP2015135662A (ja) | 情報処理装置、プログラム及び方法 | |
CN111625975A (zh) | 一种弓网动力学性能确定方法及系统 | |
CN115326066A (zh) | 一种高速列车高精度定位方法及系统 | |
CN110943473A (zh) | 一种基于广域测量系统和聚类理论的发电机同调识别方法 | |
CN110110428A (zh) | 一种基于克里金模型优选与神经网络的压气机特性预测方法 | |
CN101567838A (zh) | 一种函数链神经网络的自校正方法 | |
JP6950647B2 (ja) | データ判定装置、方法、及びプログラム |
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 |