轮轨接触力学课件.pptx

1、温泽峰，赵鑫温泽峰，赵鑫西南交通大学西南交通大学 牵引动力国家重点实验室牵引动力国家重点实验室 只够1小时讲，再加点内容。一一 轮轨接触动力力学的研究内容与对象轮轨接触动力力学的研究内容与对象二二 轮轨接触几何关系和滚动接触蠕滑率轮轨接触几何关系和滚动接触蠕滑率三三 Hertz接触理论（法向解开创工作）接触理论（法向解开创工作）四四 Carter二维滚动接触理论（切向解开创工作）二维滚动接触理论（切向解开创工作）五五 Vermeulen-Johnson无自旋三维滚动接触理论无自旋三维滚动接触理论六六 Kalker线性蠕滑理论线性蠕滑理论七七 沈氏理论沈氏理论八八 Kalker简化理论简化理论九

2、九 Kalker三维弹性体非三维弹性体非Hertz滚动接触理论滚动接触理论十十 轮轨黏着问题研究简介轮轨黏着问题研究简介十一十一 三维弹塑性滚动接触有限元建模简介三维弹塑性滚动接触有限元建模简介十二十二 轮轨接触载荷与伤损研究简介轮轨接触载荷与伤损研究简介十三十三 快速接触算法开发快速接触算法开发十四十四 接触问题杂谈接触问题杂谈十五十五 轮轨试验台简介轮轨试验台简介NorTanGContact at one pointContact at an area due to elasticityPressure(Nor)Surface shear stress(Tan)1). 法向接触：接触斑形状

3、、大小及法向应力分布 Hertz接触理论：第一个法向接触解（1882）2). 切向接触：基于法向解，求摩擦力分布（大小、方向） Carter第一个求解切向滚动接触（1926）Concentrated contact线弹性力学，小应变假设集中力作用下应力应变场椭圆形接触斑椭球形应力分布接触斑半长、最大应力无限半空间假设二次曲面积分*1 3*1 33/(4()3/(4()am PGABbn PGAB1/32max0*313 2()22PABppPabmnG12121122GGG1()cos ()BAABP载荷，A、B几何参数基于赫兹解（二次曲面、无限半空间、线弹性材料、相对狭小的椭圆接触斑、椭球应

4、力分布等等）；二维简化（平面应变问题），库伦摩擦定律，摩擦系数恒定；方形接触斑内存在方形粘着区，应力为两个椭圆分布相减；稳态滚动。1211011zPfarfP 21*0 4EPraz2212*12111 EEE21111zfPPaaBy J.J. KalkerGeneralization of Galins TheoremMathematical Programming TheoryCattaneoMindlinDuvaut-Lions支持性、基础性理论JProf. J.J. Kalker,1933-20061）可解析的滚动接触理论Carter V-J, Strip theories；V-J把

5、Carter拓展至3D，无法考虑Spin；J 提出自旋概念，被K线性理论采用；沈氏理论，小自旋2）数值滚动接触解理论K简化理论FASTSIMK精确理论CONTACTBy J.J. Kalker自旋！自旋！Kalker工作的直接基础；可用来导出Hertz理论、Cattaneo shift解和V-J理论。By J.J. Kalker滚动接触数值解建模数学基础；控制方程和边界条件虚功原理变成等效积分弱形式最小位能原理、最小余能原理，将原问题转化为等效形式的变分问题。By J.J. Kalker推导了接触问题的变分不等式（基于虚功原理）。被Kalker用于数值求解摩擦接触问题。By J.J. Kalk

6、erBy J.J. KalkerHertz接触的no-slip问题切向力作用下的切向解旋转趋势下的切向解By J.J. KalkerV-JStrip theories；Johnson自旋理论；Kalker 线性理论；沈氏理论。a) a) VermeulenVermeulen-Johnson-Johnson理论模型（无自旋三维滚动接触理论模型）理论模型（无自旋三维滚动接触理论模型） l1958年，Johnson扩展Carter二维滚动接触理论到三维弹性球滚动接触第一个三维解；l1964年，Vermeulen和Johnson又将上述研究推广到椭圆斑的接触情形，接触区中粘着区和滑动区的划分仍按Car

7、ter的研究思路。311aTafP 122212012221212 S) C)1 ( ,),0 ( ,)Txxfpx xTabpx xx x11222222121120122222 H)() 11 ( ,)Txxaxcxfpx xTabaabfPTbbaa13131CarterV-JV-J非线性蠕滑率非线性蠕滑率/ /力定律为力定律为 *3*11 1 1iiiifPTTfP212*22*1*2*22*21*11*13 , 3fPabGCfPabGC为规格化的纵横向蠕滑率 *2*1 、为V-J蠕滑系数 *22*11 CC 、Vermeulen-JohnsonVermeulen-Johnson理论

8、与其实验值的比较理论与其实验值的比较 英国Hobbs（1967）用Kalker的蠕滑系数代替V-J的蠕滑系数*3*11 1 1iiiifPTTfP当蠕滑率很小时，上式按Taylor级数展开，取一阶线性项： *00|3dTTTfPd分量形式 *2*22*22*1*11*1133abGCfpTabGCfpTV-JV-J线性蠕滑率线性蠕滑率/ /力定律力定律l值得注意，VJ理论没有考虑滚动物体的自旋效应。事实上许多滚动接触问题的自旋运动是存在的，如轮轨滚动接触、轴承滚珠或滚子的运动等，物体间的自旋运动是不可避免的，它不仅影响接触斑上总的横向蠕滑力，而且对接触斑表面产生疲劳破坏作用。l此外，VJ理论对

9、接触区中粘滑区的划分有不完善之处。如下图有线条的阴影部分划分为滑动区是错误的，是因为在这个区域里，用VJ理论确定的滑动方向和切向力方向一致，这违背了摩擦定律。 b) Johnsonb) Johnson纯自旋理论模型纯自旋理论模型 横移越大，接触角、自旋蠕滑分量越大横移越大，接触角、自旋蠕滑分量越大小自旋小自旋大自旋大自旋l1958年，Johnson另外一篇文章讨论了纯自旋问题纵、横向蠕滑率为0，no-slip，圆形接触斑；l条形理论可以解决V-J理论粘滑区划分问题。l条形理论首先是由Hains.D.J和Ollerton.E于1963年提出来，假设仅有纵向蠕滑，且接触斑为沿横向细长的稳态滚动接触

10、情形，将椭圆分成若干个平行于方向矩形条，在每一个矩形条中，按Carter的二维求解方法，确定有关要求的参数，这些参数与无关，见下图（a）。l如果接触存在粘着区的话，只可能出在若干个矩形条域的前沿区域，则在整个接触区前沿区域上形成粘着区，如下图（b）所示。后来他们用三维光弹方法，证实了他们的正确划分。l1967年Kalker又将这种研究扩展到横向蠕滑不为零和小自旋情形。Kalker研究中发现，Hains和Ollerton对粘滑区划分还适合纵横向蠕滑率不为零且自旋蠕滑率为零的稳态滚动接触情形 。l条形理论研究中，忽略矩形条之间力学行为互相之间的影响，这种理论对于沿横向方向狭长的椭圆情形非常适合，但

11、对a接近或大于b的情形会产生较大的误差 。c ) c ) 条形理论条形理论 d) d) KalkerKalker线性蠕滑理论模型线性蠕滑理论模型 Kalker于1967年在他的博士论文中，用级数方法专门讨论了具有椭圆接触区的三维稳态滚动接触情形，研究中考虑到纵横向蠕滑率和自旋滑率对接触斑蠕滑力的影响。 Kalker J.J. On thc rolling contact of two elastic bodies in the presence of dry friction. Ph.D. Thesis. The Netherland, Delft University, 1967 33323

12、2333232222211121cCCGcMcCCGcFCGcFKalkerKalker线性蠕滑率线性蠕滑率/ /力关系力关系 cabCij：蠕滑系数 l线性模型在车辆动力学中应用极为广泛，但使用时需要特别注意，仅适合小蠕滑和小自旋情形。l这里的小蠕滑和小自旋没有明确规定的临界值。建议，当椭圆接触区处于全粘着情况下可以使用该理论模型计算轮轨蠕滑力，求得轮轨之间切向力是足够精确的，但接触区滑动或出现局部滑动情况下，使用该模型就会产生较大的误差。l因为Kalker理论研究过程中，没有考虑到当接触斑滑动或部份滑动时，Coulomb定律的限制条件，即l线性模型是很重要，后来Kalker借助于它发展了简

13、化理论。沈志云、Hedrick.J.K和Elkins.J.A利用它改进了Vermeulen.Johnson三维滚动接触非线性数学模型。 2132121222121,xxfpxxpxxpd) d) 沈氏理论（小自旋情形下三维非线性蠕滑率沈氏理论（小自旋情形下三维非线性蠕滑率/ /力计算模型力计算模型 ） l在不考虑自旋的情况下，Vermeulen-Johnson理论模型是理想的三维Hertz非线性蠕滑率/力计算模型。l通常许多滚动接触物体之间除了有刚性滑动外，也同时存在相对转动。例如，机车车辆在运动过程中，轮轨接触表面之间除了存在纵横蠕滑率外，由于车轮踏面的锥度或接触斑公法线与轮对滚动轴线不垂直

14、，则自旋效应就会产生。也就是接触表面之间产生了因滚动引起的相对转动。所以，使用V-J理论求解轮轨滚动接触问题就会产生误差。l因此，沈志云、Hedrick和Elkins对V-J非线性蠕滑定律作了改造 。Shen Z Y, Hedrick J K, Elkins J. A. A comparison of alternative creep-force models for rail vehicle dynamic analysis. Proc, 8th IAVSD Symp., Cambridge, Ma.,1984, 591605 *2*22*22*1*11*1133abGCfpTabGCfp

15、TV-JV-J线性蠕滑率线性蠕滑率/ /力定律力定律333232333232222211121cCCGcMcCCGcFCGcFKalkerKalker线性蠕滑率线性蠕滑率/ /力关系力关系 *2111 1*23222223333fPc GCfPc GCc GC 假设它们所表达的切向力相等 *2111 1222233*2233c GCfPCCcc GfP 规则化了的刚性横向蠕滑率中包含了自旋 总的正则化刚性蠕滑率 212*22*1VJ非线性模型 1 1 11*3*fPfPF2/12 22 1FFFR311 1 1RfPFfP122222311 1222223312221213133Rc GCc

16、GCc GCfPFFFfPfP代入 *2111 1222233*2233c GCfPCCcc GfP 2311 13273 13RRRRRRFFFFfPfPfPfPfPFFfPfP RRFF332211MMFTFTsss缩减因子 沈-Hedrick-Elkins的理论结果与V - J实验结果比较 沈-Hedrick-Elkins的理论与CONTACT、FASTSIM结果比较 30简化，假设，不适合应力分析；接触斑内部的情况是假设的；基于赫兹假设；伪三维问题，真正三维问题求解需要电脑。理论理论接触斑应力假设接触斑应力假设蠕滑率蠕滑率- -力关系力关系V-JStrip theoryLinear theoryShen-Hedrick-ElkinsAny questions? Thank you for your attention!Thank you for your attention!? ?