CN108664707B - 基于有限元模拟的车轮轨道接触循环加卸载仿真分析方法 - Google Patents
基于有限元模拟的车轮轨道接触循环加卸载仿真分析方法 Download PDFInfo
- Publication number
- CN108664707B CN108664707B CN201810341275.XA CN201810341275A CN108664707B CN 108664707 B CN108664707 B CN 108664707B CN 201810341275 A CN201810341275 A CN 201810341275A CN 108664707 B CN108664707 B CN 108664707B
- Authority
- CN
- China
- Prior art keywords
- unloading
- loading
- model
- wheel
- contact
- 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
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]
Abstract
一种基于有限元模拟的车轮轨道接触循环加卸载仿真分析方法,主要步骤如下:(a)利用有限元软件ABAQUS,建立单轴拉压轴对称模型;(b)进行单轴拉压循环试验模拟,并根据室内试验结果和Chaboche模型,对材料非线性硬化参数进行校正;(c)建立车轮‑轨道接触模型,对材料参数和接触参数进行设置;(d)将单次循环分解为加载阶段、平稳滚动阶段和卸载阶段等3个阶段,并分别设置加载条件、约束条件、接触条件和边界条件;(e)利用restart机制,导入上一次计算结果,在此基础上对轮轨接触加‑卸载模拟,并保存卸载阶段结果,以作为下一次循环的初始条件。本发明基于非线性有限元模拟方法和接触力学理论,可对轮轨疲劳、损伤等循环加、卸载行为进行仿真分析。
Description
技术领域
本发明涉及轨道交通安全研究领域,涉及车轮-轨道接触的理论与实践计算机仿真方法。
背景技术
当前,地铁、高铁等轨道交通体系在我国得到了快速发展和大规模建设。如何评估服役期轮轨体系的安全性和可靠性,是我国轨道交通安全领域所面临的主要问题。经过参数校正后的有限元模拟技术,可以实现各种工况下的轮轨接触仿真计算,已被国内外广泛用于轮轨摩擦、断裂及损伤分析,因而可作为服役期轮轨安全评估及磨损分析的可靠预测手段。由于软、硬件计算性能及存储性能的限制,轮轨接触模拟尚存在大规模、复杂模型难以运行、循环加-卸载难以实现等缺陷。
借助于大型非线性有限元软件ABAQUS,利用其优异的接触处理性能和模型变换技术,可实现对复杂受力条件下的轮轨相互作用进行仿真分析,并大规模降低物理试验的成本,有助于推动轨道交通科技的发展。轮轨接触分析需针对车辆平稳运行状态进行评估,此前的有限元轮轨循环模拟尚未考虑到加、卸载引起的边界效应影响,因而难以消除误差。此外,尚未见轮轨循环作用实施方法的相关描述。
发明内容
本发明的目的是针对现有研究的不足,提供了一种能够消除加、卸载边界效应影响,准确描述循环荷载作用下轮轨接触的计算机仿真分析方法。
技术方案:本发明提供的一种基于有限元模拟的车轮轨道接触循环加卸载仿真分析方法,包括以下步骤:
(1)利用有限元软件ABAQUS,建立单轴拉压轴对称模型;
(2)进行单轴拉压循环试验模拟,并根据室内试验结果和Chaboche模型,对材料非线性硬化参数进行校正;
(3)建立车轮-轨道接触模型,对材料参数和接触参数进行设置;
(4)将单次循环分解为加载阶段、平稳滚动阶段和卸载阶段等3个阶段,并分别设置加载条件、约束条件、接触条件和边界条件;
(5)利用restart机制,导入上一次计算结果,在此基础上对轮轨接触加-卸载模拟,并保存卸载阶段结果,以作为下一次循环的初始条件。
所述的步骤(1)中,采用常规狗骨式(dog-bone type)光滑圆柱试样模型,其圆柱直径为14mm,圆柱长度25.4mm;
所述的步骤(2)中,对单轴循环加压载荷试验模拟,基于唯象理论和Chaboche拟合模型,利用相同条件下的室内试验结果,即应力-应变曲线校正材料的屈服强度、硬化参数力学参数;
所述的步骤(3)中,选取车轮模型和轨道模型,并设置车轮参考点及其耦合刚性面,以模拟车轴刚性体;
所述的步骤(4)中,在加载阶段,需休眠无关部件和约束关系及接触关系,仅保留本次循环所必须的车轮、轨道部件和接触关系;在平稳滚动阶段,设置所关心的监测部位和变量,以做分析所用;在卸载阶段,需保留计算结果,以做下一轮分析所用;
所述的步骤(5)中,在加载阶段Model-Edit Attributes选项中设置restart数据,在卸载阶段Output-Restart requests选项中保留计算结果。
有益效果:本发明基于车轮-轨道接触理论进行计算机仿真分析,能够消除加、卸载边界效应影响,可实现对复杂受力条件下的轮轨相互作用进行仿真分析,可解决当前轮轨接触研究难以对循环荷载作用进行描述的缺陷,适用于大规模轮轨接触疲劳、损伤等循环加、卸载行为仿真分析。
附图说明
图1为本发明的总体流程示意图。
图2为本发明的单轴拉压循环试验模拟模型图。
图3为本发明的三维轮轨接触有限元模型图。
图4为本发明的载荷幅值曲线(载荷放大系数-时间关系图)。
图5为本发明的模型变更及接触设置示例。
具体实施方式
下面结合附图和实例对本发明进行详细的描述。
本发明提供了一种基于有限元模拟的车轮-轨道接触循环加-卸载仿真分析方法,其操作流程如图1所示,具体步骤如下:
(1)利用有限元软件ABAQUS,建立狗骨式(dog-bone type)光滑圆柱试样模型,其圆柱直径为14mm,圆柱长度25.4mm,如图2(a)所示;对圆柱体内部单元的竖向应力、应变进行监测,并绘制真实应力应变曲线。进行多次加载的竖向应变云图如2(b)所示。
(2)进行单轴拉压循环试验模拟,并根据室内试验结果和Chaboche模型,对材料非线性硬化参数进行校正。循环加、卸载曲线拟合的关键在于根据唯象理论进行的棘轮效应拟合。基于Chaboche非线性随动硬化模型,首先需要确定backstress(反力)α的数量和大小,然后将其叠加,可获得循环加-卸载试验结果的棘轮效应拟合。其拟合公式如下:
其中,K为backstress的编号,εp为单轴塑性应变,CK和γK为需要校正的Chaboche参数,α为反力,αk为反力的分量,N为反力的数量。
(3)建立车轮-轨道接触模型,利用已校正好的材料参数和接触参数对车轮、轨道模型进行赋值。图3所示为直径860mm,踏面形态为LMA的车轮模型和钢轨型号为CHN60的轨道模型装配图。循环加-卸载试验需要考虑加、卸载引发的边界效应的影响,因而需要设置足够长的平稳滚动区域。本例所示的轨道模型长1.5m,由12万个C3D8单元组成;车轮由4万个C3D8单元组成。
(4)需要在车轮中心参考点上施加法向荷载P、切向荷载Q和以边界条件形式施加的转角R。单个循环P、Q、R的幅值曲线如图4所示。在0.5ms内,逐渐对车轮中心施加法向荷载至设定值;在0.5ms-3.843ms内,保持竖向荷载不变;在3.843-4.0ms内,将法向荷载逐渐减小至零。在法向加载完成后,进行切向加载和滚动加载;同时,在法向卸载完成前,进行切向卸载和滚动卸载。
(5)多次循环加载需要利用ABAQUS的model change功能和restore机制。图5所示为5个循环的轮轨接触设置示例。对于单个循环,需要休眠无关部件和约束关系、接触关系;同时激活相关部件及其约束关系、接触关系。从图中可以看到,每个循环均包含3个分析步,分别对应于加载阶段、平稳滚动阶段和卸载阶段。在加载阶段,需休眠无关部件和约束关系及接触关系,仅保留本次循环所必须的车轮、轨道部件和接触关系;在平稳滚动阶段,设置所关心的监测部位和变量,以做分析所用;在卸载阶段,需保留计算结果,以做下一轮分析所用,及在Output-Restart requests中保留计算结果。对于下一个循环的加载阶段,则需在Model-Edit Attributes中选择读取上一轮循环的结果数据作为起点进行分析。理论上,本方法可适用于无限次数车轮-轨道滚动接触模拟,并可单独提取特定循环次数的结果进行分析。
上文中,Chaboche拟合模型即Chaboche非线性随动硬化模型。材料参数包括材料硬化参数和材料的屈服强度;接触参数指轮轨接触面的参数,包括切向、法向接触面滑动设置和摩擦设置。
Claims (2)
1.一种基于有限元模拟的车轮轨道接触循环加卸载仿真分析方法,特征在于包括以下步骤:
(1)利用有限元软件ABAQUS,建立单轴拉压轴对称模型;
(2)进行单轴拉压循环试验模拟,基于唯象理论和Chaboche拟合模型,利用相同条件下的室内试验结果,即应力-应变曲线校正材料的屈服强度、硬化参数力学参数,对材料的屈服强度以及硬化参数进行校正;
(3)建立车轮-轨道接触模型,对材料参数和接触参数进行设置:选取车轮模型和轨道模型,并设置车轮参考点及其耦合刚性面,以模拟车轴刚性体;
(4)单次循环:将单次循环分解为加载阶段、平稳滚动阶段和卸载阶段3个阶段,并分别设置加载条件、约束条件、接触条件和边界条件;在加载阶段,需休眠无关部件和约束关系及接触关系,仅保留本次循环所必须的车轮、轨道部件和接触关系;在平稳滚动阶段,设置所关心的监测部位和变量,以做分析所用;在卸载阶段,需保留计算结果,以做下一轮分析所用;
(5)多次循环:利用restart机制,导入上一次计算结果,在此基础上对轮轨接触加-卸载模拟,并保存卸载阶段结果,以作为下一次循环的初始条件;在加载阶段Model-EditAttributes选项中设置restart数据,在卸载阶段Output-Restart requests选项中保留计算结果。
2.根据权利要求1所述的一种基于有限元模拟的车轮轨道接触循环加卸载仿真分析方法,其特征是,所述单轴拉压轴对称模型采用常规狗骨式dog-bone type光滑圆柱试样模型,其圆柱直径为14mm,圆柱长度25.4mm;
对圆柱体内部单元的竖向应力、应变进行监测,并绘制真实应力应变曲线;
所述步骤(2)中,循环加、卸载曲线拟合的关键在于根据唯象理论进行的棘轮效应拟合;基于Chaboche拟合模型,首先需要确定backstress反力α的数量和大小,然后将其叠加,获得循环加-卸载试验结果的棘轮效应拟合;其拟合公式如下:
其中,K为反力的编号,εp为单轴塑性应变,CK和γK为需要校正的Chaboche参数,α为反力,αk为反力的分量,N为反力的数量;
所述步骤(3)中,利用已校正好的材料的屈服强度和硬化参数对车轮模型和轨道模型进行赋值;车轮模型为直径860mm,踏面形态为LMA,轨道模型的钢轨型号为CHN60;循环加-卸载试验需要考虑加、卸载引发的边界效应的影响,因而需要设置足够长的平稳滚动区域;轨道模型长1.5m,由12万个C3D8单元组成;车轮由4万个C3D8单元组成;
所述步骤(4)中,需要在车轮中心参考点上施加法向荷载P、切向荷载Q和以边界条件形式施加的转角R;在0.5ms内,逐渐对车轮中心施加法向荷载至设定值;在0.5ms-3.843ms内,保持竖向荷载不变;在3.843-4.0ms内,将法向荷载逐渐减小至零,在法向加载完成后,进行切向加载和滚动加载;同时,在法向卸载完成前,进行切向卸载和滚动卸载;
所述步骤(5)中,多次循环加载需要利用ABAQUS的model change功能和restore机制;对于单个循环,需要休眠无关部件和约束关系、接触关系,同时激活相关部件及其约束关系、接触关系,在加载阶段,需休眠无关部件和约束关系及接触关系,仅保留本次循环所必须的车轮、轨道部件和接触关系;在平稳滚动阶段,设置所关心的监测部位和变量,以做分析所用;在卸载阶段,需保留计算结果,以做下一轮分析所用,及在Output-Restartrequests中保留计算结果;对于下一个循环的加载阶段,则需在Model-Edit Attributes中选择读取上一轮循环的结果数据作为起点进行分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810341275.XA CN108664707B (zh) | 2018-04-17 | 2018-04-17 | 基于有限元模拟的车轮轨道接触循环加卸载仿真分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810341275.XA CN108664707B (zh) | 2018-04-17 | 2018-04-17 | 基于有限元模拟的车轮轨道接触循环加卸载仿真分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108664707A CN108664707A (zh) | 2018-10-16 |
CN108664707B true CN108664707B (zh) | 2021-04-09 |
Family
ID=63783592
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810341275.XA Expired - Fee Related CN108664707B (zh) | 2018-04-17 | 2018-04-17 | 基于有限元模拟的车轮轨道接触循环加卸载仿真分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108664707B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109376503B (zh) * | 2018-12-29 | 2022-12-09 | 中南大学 | 考虑轮轨滚动接触的高速铁路路基耦合建模分析方法 |
CN110926942B (zh) * | 2019-11-19 | 2022-01-07 | 西安理工大学 | Abaqus中铁轨滚动接触疲劳裂纹的数值分析方法 |
CN111209701B (zh) * | 2020-01-02 | 2023-09-12 | 中车青岛四方机车车辆股份有限公司 | 一种轨道列车车体强度仿真方法及装置 |
CN112131678B (zh) * | 2020-08-18 | 2024-02-20 | 中铁物总运维科技有限公司 | 一种铁路钢轨打磨廓形设计方法 |
CN114254542B (zh) * | 2021-12-29 | 2023-07-28 | 武汉大学 | 球墨铸铁分段管道滑入式接口抗拉弯性能预测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102841958A (zh) * | 2012-07-10 | 2012-12-26 | 北京交通大学 | 高速铁路桥上有砟轨道结构力学分析及选型方法 |
JP2013018353A (ja) * | 2011-07-11 | 2013-01-31 | Bridgestone Corp | シミュレーション方法及びシミュレーション装置 |
CN103150458A (zh) * | 2013-04-01 | 2013-06-12 | 中南大学 | 车辆-轨道-桥梁-地基基础耦合系统及其动力分析方法 |
CN105224758A (zh) * | 2015-10-13 | 2016-01-06 | 河海大学 | 一种用于设计计算的有砟轨道高速铁路有限元建模分析方法 |
CN105677985A (zh) * | 2016-01-11 | 2016-06-15 | 北京交通大学 | 城市轨道交通板式减振轨道结构分析选型方法 |
CN107451384A (zh) * | 2016-05-31 | 2017-12-08 | 中南大学 | 一种列车‑轨道‑结构耦合系统高效动力分析方法 |
-
2018
- 2018-04-17 CN CN201810341275.XA patent/CN108664707B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013018353A (ja) * | 2011-07-11 | 2013-01-31 | Bridgestone Corp | シミュレーション方法及びシミュレーション装置 |
CN102841958A (zh) * | 2012-07-10 | 2012-12-26 | 北京交通大学 | 高速铁路桥上有砟轨道结构力学分析及选型方法 |
CN103150458A (zh) * | 2013-04-01 | 2013-06-12 | 中南大学 | 车辆-轨道-桥梁-地基基础耦合系统及其动力分析方法 |
CN105224758A (zh) * | 2015-10-13 | 2016-01-06 | 河海大学 | 一种用于设计计算的有砟轨道高速铁路有限元建模分析方法 |
CN105677985A (zh) * | 2016-01-11 | 2016-06-15 | 北京交通大学 | 城市轨道交通板式减振轨道结构分析选型方法 |
CN107451384A (zh) * | 2016-05-31 | 2017-12-08 | 中南大学 | 一种列车‑轨道‑结构耦合系统高效动力分析方法 |
Non-Patent Citations (1)
Title |
---|
基于摩擦自激振动引起高速列车车轮多边形磨耗的仿真研究;王科;《中国优秀硕士学位论文全文数据库》;20170715;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108664707A (zh) | 2018-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108664707B (zh) | 基于有限元模拟的车轮轨道接触循环加卸载仿真分析方法 | |
CN110132515B (zh) | 一种基于模型更新的时程级迭代实时混合试验方法 | |
CN106979861A (zh) | 齿轮接触疲劳全寿命评估方法及装置 | |
CN105260574A (zh) | 一种基于临界面法疲劳破坏准则的高周多轴疲劳寿命预测方法 | |
CN113868777B (zh) | 一种副车架强度台架试验载荷获取方法和装置 | |
CN104461593A (zh) | 差分升级包的制作方法及装置 | |
CN105628486A (zh) | 测量页岩岩石力学性质的方法及装置 | |
CN106326587A (zh) | 一种测力轮对轮轨力的修正方法及系统 | |
CN114330042A (zh) | 基于sn曲线的程序载荷谱编制方法、系统和存储介质 | |
Azizi et al. | Identification of nonlinear viscoelastic models of flexible polyurethane foam from uniaxial compression data | |
CN106803007B (zh) | 一种多次振动处理修复铜薄膜试件疲劳损伤的方法 | |
CN113378421B (zh) | 一种机械弹性车轮径向疲劳寿命评价及预测方法 | |
Su | A road load data processing technique for durability optimization of automotive products | |
CN104166778B (zh) | 汽车座椅导轨接触部位刚度和阻尼优化方法 | |
CN104517012B (zh) | 一种计算非线性结构的应变时间历程的方法和装置 | |
CN109877832A (zh) | 一种位置确定方法、系统及相关装置 | |
CN110955933B (zh) | 一种基于响应面法的机械结构模糊疲劳可靠度计算方法 | |
CN104834795A (zh) | 包带连接结构接触摩擦非线性特性模拟方法及系统 | |
WO2024082486A1 (zh) | 一种车辆段振动噪声仿真荷载实现方法 | |
CN110989497A (zh) | 基于迭代优化的多轴多激励振动控制组合选择方法及系统 | |
CN114722704A (zh) | 车轮磨损预测网络模型训练方法、装置 | |
CN105354374B (zh) | 连接点的刚度模拟方法及装置 | |
Wang et al. | Development in the acquisition of vehicle loads integrated with a rigid and flexible multi-body model | |
CN114564797A (zh) | 一种道岔曲梁可靠性计算方法及系统 | |
Li et al. | Parameter identification of in-plane flexible ring tyre model based on static load-deflection data: some insights |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210409 |
|
CF01 | Termination of patent right due to non-payment of annual fee |