基于三维动态有限元模型的轮轨磨耗数值分析.docVIP

　 90中国铁道科学第29卷 轨轮和车轮的表面硬度的测量结果见表1。轨轮旋转6×105圈相当于轨道通过1亿t轴重。文献[10]对双轮盘磨耗实验做了详细的描述,本文数值分析以此为目标,验证轮轨磨耗预测方法的有效性。 重 褪 基 彝 链 图1轮轨磨耗试验装置简图 图2轮轨磨耗试验中的轨轮廓形变化曲线 表l磨耗试验轮轨表面硬度变化情况HV 2轮轨磨耗的预测方法 对于轮缘和轨距角的接触计算,基于Hertz理论的接触算法受半空间假设和线性材料模型的限制,但有限元方法不受这2个假设的限制。用三维的弹塑接触有限元方法可以模拟车轮和轨轮的瞬态滚动接触。轮轨磨耗预测方法可以分成4个主要组成部分:基于FEM的瞬态滚动接触仿真,接触和材料的本构模型,磨耗模型和自适应网格调整技术,如图3所示。 图3轮轨磨耗的预测方法流程 2.1接触本构模型 ABAQUS软件包用于轮轨滚动接触的仿真计算。该软件包可以进行三维结构的非线性静态和瞬态分析计算,在计算过程中,可以考虑应力波和惯性影响[11]。 接触本构模型包括接触法向应力的定义和接触切向应力的定义。在某一点处2个接触面间的法向应力P可以定义成相互渗透h的函数,该模型可表示为 p:『o矗j o(1 p 2 J LlJ 1 1Kh h≥0。。式中:K为接触刚度。 接触切向应力可定义成Coulomb摩擦模型,在标准的Coulomb摩擦模型中,假定等效摩擦应力珀一￣/者+磊,当如小于临界摩擦应力zcri。时,表明没有相对运动发生,确。和P成以下正比关系: %。一肋(2式中:口为摩擦系数,它可以定义成接触压力或接触点处相对滑移速度的函数。 若r舶=“扪就会出现滑移。假如摩擦特性是各项异性的,滑移方向和摩擦力方向一致,又可表示成如下形式: 旦=旦(3 req口q 其中,　 第4期基于三维动态有限元模型的轮轨磨耗数值分析91 氓一√乞彳j二丽(4P2芦k+(卢s—ttke--dcj,q(5式中:祝为i方向的相对滑移速度;记为等效滑移式中:肚为动摩擦系数;户e为静摩擦系数;dc为速度。衰减系数。 2.2有限元模型 根据轮轨磨耗试验,用有限元法模拟轮轨滚轮 的三维瞬态滚动接触。因为轮轨的FEM模型是关于轴向对称的,轨头和车轮踏面的网格可以分别通过旋转基准平面内的2D网格生成。为了减少整个模型的单元和节点数,其他部分的网格可用别的方法生成三维网格单元。最小单元长度取2rnIn,以保证计算精度和不浪费太多的机时。轮轨滚轮用8节点三维体单元(C3D8R进行离散,C3DSR单元既适用于隐式求解,也适用于显式求解。该有限元模型总共有232041个单元和261963个节点,三维有限元模型如图4所示。 (b接触点附近的网格放大图 图4三维轮轨有限元模型 在整个轮轨滚动接触分析过程中,需要把ABAQUS的隐式求解(静态和显式求解(瞬态联合使用。首先用隐式求解(静态方法分析轮轨滚轮在轴向和径向载荷作用下的轮轨接触状态,然后把静态接触状态和初始的速度状态作为显式求解(瞬态求解的初始状态,再进行显式求解,这样可以使得瞬态计算在很短的时间内就可以进入稳定状态。 轮、轨材料本构模型为经典的双线性弹塑模本构模型,材料本构曲线如图5所示。 轮轨表面的接触形式采用自动的面面接触运算方法。在轮轨接触特性中,静摩擦、动摩擦和粘滑运动可用摩擦模型表示,在该摩擦模型中,摩擦系数是接触面的相对滑动速度的函数。ABAQUS也提供了1个指数衰减摩擦模型[15|,该摩擦模型的数学表达式为 ￡ 善 .R 趔 虚变唐 图5轮(曲线A轨(曲线B材料的非线性本构曲线 文献[13]根据Poiree在里昂铁路上车辆沿钢轨滑动的制动试验数据(试验时轮轨滑动速度高达22.0m?S-1归纳了轮轨动摩擦系数与相对滑移速度的关系曲线,如图6中的试验曲线所示。用式(5对试验曲线进行拟合可得:胁一.o.45,地=O.2339,d。一4.8950×10一,拟合结果如图6中的拟合曲线所示。 氟 惴 鼙 慢 相对滑移速度/(nun?S以 图6轮轨动摩擦系数与相对滑移速度的关系曲线 2.3磨耗模型 文献[1]一文献E53研究认为,在轮轨接触表面的磨损量与能量耗散(由摩擦力或蠕滑力和蠕滑率产生的摩擦功有关。本文之所以选用Ar—chard磨损模型,是因为它可以在摩擦学中更加普遍地用来模拟由摩擦副相对滑动产生的粘着磨损。例如,在文献[14]中它已成功地用于滚动轴承的磨耗预测,轴承的滚动接触和轮轨的滚动接触有很多相似之处。在文献[6]中,Archard的磨损模型已被用于铁道车轮踏面的磨耗预测分析。 在Archard磨损模型[15“6]中,材料的磨损量 和相对滑移距离以及正压力成正比,和材料的表面　 　 第4期 基于三维动态有限元模型的轮轨磨耗数值分析 93 Yi一∑wii(而--3If (13 j