CN114781158B - 一种简化轮轨非赫兹切向接触计算的等效方法 - Google Patents
一种简化轮轨非赫兹切向接触计算的等效方法 Download PDFInfo
- Publication number
- CN114781158B CN114781158B CN202210422625.1A CN202210422625A CN114781158B CN 114781158 B CN114781158 B CN 114781158B CN 202210422625 A CN202210422625 A CN 202210422625A CN 114781158 B CN114781158 B CN 114781158B
- Authority
- CN
- China
- Prior art keywords
- contact
- formula
- wheel
- hertz
- curvature
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
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)
- Machines For Laying And Maintaining Railways (AREA)
Abstract
本发明涉及铁道工程轮轨接触计算领域,具体为一种简化轮轨非赫兹切向接触计算的等效方法,其包括以下步骤:S1、已知两个接触体的弹性模量和泊松比分别为E和v,纵向曲率接触斑纵向边界a(y),法向接触应力分布p0(y),将接触斑划分为若干沿滚动方向的条带,假设每个接触条带位于局部椭圆的中心,求解局部接触椭圆的短半轴b(y)和横向曲率B(y);S2、定义两个函数和S3、定义c=1e‑9,得到h1与h2,h1=f1(c)·f1(1‑c),h2=f2(c)·f2(1‑c);S4、如果h1<0,则a>b,根据公式得到e和b值;如果h2<0,则a<b,根据公式得到e和b值;S5、共形接触时,设定b=10a。本发明将基于赫兹假设的切向接触算法应用于非赫兹工况,从而大幅度提升非赫兹切向接触的计算效率。
Description
技术领域
本发明涉及铁道工程轮轨接触计算领域,具体是一种简化轮轨非赫兹切向接触计算的等效方法。
背景技术
轮轨系统是连接车辆与轨道的重要纽带,其中,宏观的轮轨蠕滑力是计算车辆-轨道耦合动力学的重要参量,轮轨滚动接触力学行为则为求解轮轨滚动接触疲劳和磨损问题提供必要的基础解。目前在铁路行业中得到广泛应用的接触模型主要有FASTSIM,沈氏理论与Polach方法以及CONTACT。然而,这些模型仍存在两点不足:一方面,如FASTSIM,沈氏理论和Polach方法等简化模型均是基于赫兹假设建立,无法更为准确地描述轮轨非赫兹接触斑;另一方面,如CONTACT等准确的计算模型需要大量的迭代过程,其计算效率难以满足如车辆动力学、轮轨磨损预测等大规模计算要求。
现有技术中较高效的模型为STRIPES模型,STRIPES模型采用的是局部椭圆法。STRIPES使用改进的FASTSIM来进行轮轨切向接触分析。首先将接触斑沿纵向划分为若干条带,并假设每个条带位于局部椭圆的中心,如图1所示,STRIPES模型针对每个条带计算出局部椭圆,基于局部椭圆计算每个条带上局部的Kalker线性蠕滑系数Cij和FASTSIM柔度系数,其求解局部椭圆的步骤如下:
首先通过法向接触得到每个条带对应的局部接触椭圆长半轴a(y):式中,g(y)为轮轨间隙函数,A(y)为修正后的轮轨纵向组合曲率。根据轮轨实际型面求解轮轨横向组合曲率式中,Cw和Cr分别是车轮和钢轨型面的横向曲率。
B(y)需要进行平滑和修正,平滑是由于某些位置的轮轨型面曲率半径变化剧烈,修正的原因是根据轮轨型面曲率得到的B(y)有可能为负值,但是根据Hertz方法中,B(y)必须恒为正,因此需要将B(y)中的负值替换掉。根据a(y)与B(y)求解局部接触椭圆短半轴:式中δ为刚性渗透量,ε为渗透量修正系数,依据刚性接触点处的轮轨曲率确定:式中,n和r由赫兹理论计算得到。依据每个条带的a与b值,可得到每个条带对应的Kalker线性蠕滑系数Cij和FASTSIM柔度系数,进而求解非赫兹工况下的切应力。
在STRIPES模型中,求解a(y)之前,需要对潜在接触区域的轮轨横向组合曲率B(y)进行平滑与修正。主要原因如下:首先,根据Hertz理论,B(y)应该恒为正值。对于大多数轮轨接触工况,通过轮轨型面实际曲率计算得到的B(y)为负值,此时需要把负曲率变为正值,否则不能应用于基于Hertz理论的切向接触模型。其次,一般来讲,接触斑内轮轨型面曲率半径变化剧烈,而接触斑形状在曲率半径不连续点处会出现突变。实际上,接触斑形状应该是光滑的,因此需要对其横向曲率进行人为的平滑处理。而平滑方式的选择较主观,且接触结果易受平滑方式的影响,有文献采用移动平滑方法对磨耗的LMA车轮踏面曲率进行平滑,对比了不同平滑窗口的移动平滑方法对轮轨法向力的影响,如图2所示。由此可以发现平滑方式的选择会明显改变轮轨接触力。在工程计算时,很难评判哪种平滑方式最优,导致STRIPES模型稳定性不足,计算精度依赖于人为的主观干预。因此,以往研究中通过修正轮轨型面曲率的方法的计算稳定性有待提高。
发明内容
本发明目的是针对背景技术中存在的问题,提出一种将非椭圆接触斑等效为若干局部椭圆的方法,可将基于Hertz假设建立的切向接触模型应用于非赫兹接触工况,从而大幅度提升非赫兹切向接触的计算效率的简化轮轨非赫兹切向接触计算的等效方法。
本发明的技术方案:一种简化轮轨非赫兹切向接触计算的等效方法,包括以下步骤:
S1、轮轨接触斑压应力分布p0(y)、纵向接触边界a(y)和车轮在接触点处的滚动半径R作为已知条件给定,假设轮轨接触满足弹性半空间假设,则纵向曲率在所有条带保持恒定,将接触斑划分为若干沿滚动方向的条带,假设每个接触条带位于局部椭圆的中心,为了在每一个条带上应用切向接触算法,需要求解接触斑内横向曲率分布B(y)和每个局部椭圆的横向半轴长度b(y);
S2、根据Hertz理论,如果a>b,公式一,公式二,如果a<b,公式三,公式四, 上述公式中,K(e)和E(e)分别是第一类椭圆积分和第二类椭圆积分;获得e2后,椭圆接触斑横向半轴长度b和横向组合曲率B可进行求解;但不能提前知道a与b的大小关系,发明TRIAL算法求解b和B;
S4、定义一个极小值c=1e-9,得到h1与h2,公式七,h1=f1(c)·f1(1-c),公式八,h2=f2(c)·f2(1-c);
S5、如果h1<0,则a>b,e的值可通过f1(e)=0得到,b由公式二求得;如果h2<0,则a<b,e的值可通过f2(e)=0得到,b由公式四求得;
S6、此外,存在一种特殊的情况e=1,B=0,即共形接触,由于无法通过椭圆积分求解该种接触类型,因而其不包含于S5,因此,设定b=10a。
与现有技术相比,本发明具有如下有益的技术效果:提出一种局部椭圆算法TRIAL,利用法向应力分布和接触斑边界计算得到每一横向位置处的横向曲率与每个条带的等效椭圆宽度b,可将基于赫兹假设的切向接触算法(如FASTSIM)应用于非赫兹工况,从而大幅度提升非赫兹切向接触的计算效率。且利用所提出的方法计算求得局部椭圆,能够避免由轮轨型面获得的横向组合曲率B(y)出现负值或不平滑的情况,即克服STRIPES方法的局限性。
附图说明
图1为STRIPES模型求解等效椭圆示意图;
图2为横向曲率平滑对法向力的影响的示意图;
图3为TRIAL方法流程图;
图4为TRIAL与通过轮轨型面计算得到的横向组合曲率对比示意图;
图5为旋转体示意图;
图6为双半椭圆接触的法向解示意图;
图7为切向接触解示意图。
具体实施方式
实施例一
如图1-7所示,本发明提出的一种简化轮轨非赫兹切向接触计算的等效方法,包括以下步骤:
S1、轮轨接触斑压应力分布p0(y)、纵向接触边界a(y)和车轮在接触点处的滚动半径R作为已知条件给定,假设轮轨接触满足弹性半空间假设,则纵向曲率在所有条带保持恒定,将接触斑划分为若干沿滚动方向的条带,假设每个接触条带位于局部椭圆的中心,为了在每一个条带上应用切向接触算法,需要求解接触斑内横向曲率分布B(y)和横向半轴长度b(y);
S2、假设两个接触体的材料特性相等,其弹性模量和泊松比分别为E和v,轮轨椭圆接触斑长半轴为a,最大应力为p0;根据Hertz理论,,求解局部接触椭圆的短半轴b或横向曲率B;如果a>b,公式一,公式二, 如果a<b, 公式四, 上述公式中,K(e)和E(e)分别是第一类椭圆积分和第二类椭圆积分;获得e2后,椭圆接触斑横向半轴长度b和横向组合曲率B可进行求解;但不能提前知道a与b的大小关系,发明TRIAL算法求解b和B;
S4、定义一个极小值c=1e-9,得到h1与h2,公式七,h1=f1(c)·f1(1-c),公式八,h2=f2(c)·f2(1-c);
S5、如果h1<0,则a>b,e的值可通过f1(e)=0得到,b由公式二求得;如果h2<0,则a<b,e的值可通过f2(e)=0得到,b由公式四求得;
S6、此外,存在一种特殊的情况e=1,B=0,即共形接触,由于无法通过椭圆积分求解该种接触类型,因而其不包含于S5,因此,设定b=10a。
本实施例提出一种局部椭圆算法TRIAL,利用法向应力分布和接触斑边界计算得到每一横向位置处的横向曲率,可将基于赫兹假设的切向接触算法(如FASTSIM)应用于非赫兹工况,从而大幅度提升非赫兹切向接触的计算效率。且利用所提出的方法计算求得局部椭圆,能够避免由轮轨型面获得的横向组合曲率B(y)出现负值或不平滑的情况,即克服STRIPES方法的局限性。
以LMb车轮踏面与CHN60钢轨型面为例,钢轨轨底坡设为1:40,车轮半径为460mm。图4对比了通过本方法计算得到的横向组合曲率B以及由轮轨型面得到的结果,其中轮对横移量为0与6mm时。可以发现,1、由轮轨型面得到的横向组合曲率出现了负值,导致该结果无法与基于Hertz假设的切向接触模型联合使用;2、在接触点附近,由轮轨型面获得的曲率结果存在突变,这会导致接触结果不准确。相较而言,本发明提出的TRIAL方法改善了此问题,得到的横向组合曲率变化较为平缓且数值未出现负值,不需要再进行修正即可直接应用到接触模型。
实施例二
如图1-7所示,列举一个利用本发明提出的一种简化轮轨非赫兹切向接触计算的等效方法求解局部接触椭圆进而得到切向接触解的应用案例,具体步骤如下:
S1、假设一个旋转体与平面发生接触,该旋转体的横向曲率半径在接触点附近由100mm变化为500mm,旋转体的滚动圆半径为460mm。以Kalker的CONTACT程序所得结果作为参考。定义接触坐标系,如图5所示,原点O为接触点,x轴为纵向,z轴垂直于接触平面向下,y轴通过右手法则确定,垂直于x轴与z轴,指向右侧。
S2、已知轮轨接触点(0,0)和滚动圆半径R=460mm,通过CONTACT法向算法得到轮轨接触斑每个条带对应的接触斑纵向边界a(y),以及法向应力分布p0(y),结果如图6(a)和(b)所示,纵向曲率A=1/2/R,R为接触点处的滚动圆半径。
S4、定义c=1e-9,得到h1与h2,h1=f1(c)·f1(1-c),h2=f2(c)·f2(1-c)。
S5、如果h1<0,则e2数值范围为[c,1-c],则a>b,e的值可通过f1(e)=0得到,求得e之后,通过下式求解接触斑横向半轴长度b(y)与横向组合曲率B(y),如果h2<0,则e2数值范围为[c,1-c],则a<b,e的值可通过f2(e)=0得到,求得e之后,通过下式求解接触斑横向半轴长度b(y)与横向组合曲率B(y),e=1时,不能采用椭圆积分,它存在于共形接触工况,此时B=0,这种情况不能采用基于Hertz的切向接触模型,因此,引入一个假设,b(y)=10a(y),B(y)=c。
S6、基于局部椭圆计算每个条带上局部的Kalker线性蠕滑系数Cij,Cij是关于接触斑长短半轴a,b的比值的函数。
S8、在每个条带上应用FASTSIM模型,结果如图7所示,可以发现通过TRIAL+FASTSIM得到的切应力与CONTACT结果吻合。
上面结合附图对本发明的实施方式作了详细说明,但是本发明并不限于此,在所属技术领域的技术人员所具备的知识范围内,在不脱离本发明宗旨的前提下还可以作出各种变化。
Claims (1)
1.一种简化轮轨非赫兹切向接触计算的等效方法,其特征在于,包括以下步骤:
S1、轮轨接触斑压应力分布p0(y)、纵向接触边界a(y)和车轮在接触点处的滚动半径R作为已知条件给定,假设轮轨接触满足弹性半空间假设,则纵向曲率在所有条带保持恒定,将接触斑划分为若干沿滚动方向的条带,假设每个接触条带位于局部椭圆的中心,为了在每一个条带上应用切向接触算法,需要求解接触斑内横向曲率分布B(y)和每个局部椭圆的横向半轴长度b(y);
S2、根据Hertz理论,如果a>b,公式一,公式二,如果a<b,公式三,公式四, 上述公式中,E和v分别为弹性模量和泊松比,a为轮轨椭圆接触斑长半轴,p0为最大应力,K(e)和E(e)分别是第一类椭圆积分和第二类椭圆积分;获得e2后,椭圆接触斑横向半轴长度b和横向组合曲率B可进行求解;但不能提前知道a与b的大小关系,因此发明TRIAL算法求解e、b和B;
所述的TRIAL算法包括步骤S3-S6:
S4、定义一个极小值c=1e-9,得到h1与h2,公式七,h1=f1(c)·f1(1-c),公式八,h2=f2(c)·f2(1-c);
S5、如果h1<0,则a>b,e的值可通过f1(e)=0得到,b由公式二求得;如果h2<0,则a<b,e的值可通过f2(e)=0得到,b由公式四求得;
S6、此外,存在一种特殊的情况e=1,B=0,即共形接触,由于无法通过椭圆积分求解该种接触类型,因而其不包含于S5,因此,设定b=10a。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210422625.1A CN114781158B (zh) | 2022-04-21 | 2022-04-21 | 一种简化轮轨非赫兹切向接触计算的等效方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210422625.1A CN114781158B (zh) | 2022-04-21 | 2022-04-21 | 一种简化轮轨非赫兹切向接触计算的等效方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114781158A CN114781158A (zh) | 2022-07-22 |
CN114781158B true CN114781158B (zh) | 2023-05-05 |
Family
ID=82432100
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210422625.1A Active CN114781158B (zh) | 2022-04-21 | 2022-04-21 | 一种简化轮轨非赫兹切向接触计算的等效方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114781158B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117744236A (zh) * | 2023-11-07 | 2024-03-22 | 西南交通大学 | 一种加速轮轨接触边界元计算的预估方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108984933A (zh) * | 2018-07-25 | 2018-12-11 | 太原科技大学 | 弹流润滑条件下计算滚动轴承载荷和压力的边界元法 |
CN111368444A (zh) * | 2020-03-09 | 2020-07-03 | 南京理工大学 | 一种轮轨滚动接触斑分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9861319B2 (en) * | 2015-03-23 | 2018-01-09 | University Of Kentucky Research Foundation | Noncontact three-dimensional diffuse optical imaging of deep tissue blood flow distribution |
-
2022
- 2022-04-21 CN CN202210422625.1A patent/CN114781158B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108984933A (zh) * | 2018-07-25 | 2018-12-11 | 太原科技大学 | 弹流润滑条件下计算滚动轴承载荷和压力的边界元法 |
CN111368444A (zh) * | 2020-03-09 | 2020-07-03 | 南京理工大学 | 一种轮轨滚动接触斑分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114781158A (zh) | 2022-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114781158B (zh) | 一种简化轮轨非赫兹切向接触计算的等效方法 | |
CN104420378B (zh) | 轮胎加固钢帘线及使用其的子午线轮胎 | |
EP2228556B1 (en) | Method of manufacturing a roller bearing | |
CN111368444A (zh) | 一种轮轨滚动接触斑分析方法 | |
CN104878667A (zh) | 道岔转辙区的钢轨打磨方法 | |
CN102672059A (zh) | 根据仿真冲压工件厚度确定出模具凹凸模修改型面的方法 | |
CN105195574B (zh) | 大截面尺寸构件的冷弯成形方法 | |
CN116542023B (zh) | 一种模拟轮轨共形接触的计算方法、系统及存储介质 | |
EP2979075B1 (en) | Tire uniformity improvement using estimates based on convolution/deconvolution with measured lateral force variation | |
CN107145683B (zh) | 一种UniTire轮胎模型参数的辨识方法 | |
CN110738728B (zh) | 基于线性组合过渡算法的叶片修复模型重建方法 | |
JP7091823B2 (ja) | タイヤのトレッド部の摩耗予測方法 | |
CN114925498B (zh) | 考虑曲率过渡的三维弹塑性轮轨法向载荷快速计算方法 | |
CN110929354A (zh) | 一种高铁驱动电机吊挂板簧端部平直段长度的设计方法 | |
CN112507410B (zh) | 轨道梁图纸的生成方法和生成装置 | |
Wu et al. | Non-Hertzian conformal contact at wheel/rail interface | |
JP4429841B2 (ja) | ころ軸受 | |
CN105930607B (zh) | 非端部接触式少片端部加强型主副簧各片应力的计算方法 | |
Knothe et al. | Determination of the tangential stresses and the wear for the wheel-rail rolling contact problem | |
CN111553029B (zh) | 一种棒材矫直弹复预测方法 | |
CN115408791A (zh) | 一种考虑裂纹位置的故障齿轮时变啮合刚度计算方法 | |
US20200324575A1 (en) | Non-pneumatic tire having offset spokes | |
EP2979074B1 (en) | Tire uniformity estimation using estimates based on convolution/deconvolution | |
CN115031590B (zh) | 一种承担接头载荷的加筋壁板结构稳定性优化方法 | |
CN111159954B (zh) | 用于弹性元件的自由型面网格布局及有限元分析方法、系统及介质 |
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 |