一种干滑动摩擦热-应力-磨损分步耦合的模拟方法
【专利摘要】本发明涉及一种干滑动摩擦热-应力-磨损分步耦合的模拟方法，步骤为：将摩擦热-应力-磨损耦合过程分步；对热-应力耦合分析进行初始化设置；采用瞬态热传导分析程序依次完成N个增量步的热传导分析，获得输出温度场；对应力-磨损分析进行初始化设置；设置有限元模型应力分析的单元类型，进行应力分析，获得输出接触压力场；根据温度场、接触压力场及两接触面接触节点的相对滑移速率，计算磨损量并确定其空间方向；计算接触压力偏差；根据磨损量和方向，修正接触节点的位移，更新有限元模型；判断接触压力的偏差是否超过容差，更新接触压力；判断模拟过程是否结束。本发明计算精度和计算效率较高且适用于干滑动摩擦热-应力-磨损强耦合问题的模拟。
【专利说明】一种干滑动摩擦热-应力-磨损分步耦合的模拟方法
【技术领域】
[0001] 本发明涉及一种干滑动摩擦耦合模拟方法，特别是关于一种干滑动摩擦热-应 力-磨损分步耦合的模拟方法。
【背景技术】
[0002] 干滑动摩擦界面存在摩擦热-应力-摩擦磨损耦合现象。摩擦热主要产生于两 个接触面之间的摩擦接触表层，摩擦热的热流密度与切向摩擦应力和相对滑移速度有关。 摩擦热会直接导致接触面的温度迅速升高，迅速升高的温度致使接触面发生热变形，从而 最终使得接触压力在接触面上重新不均匀的分布，在接触压力和相对滑移速度的综合作用 下，接触面逐渐磨损。由于接触面上接触压力分布不均以及接触面不同位置处相对滑移速 度存在差异，因此，接触面不同位置处的磨损量不同。这种接触面上磨损量的不均匀性导致 接触压力的再次重新分布，进而改变切向摩擦应力的分布和热流密度的大小。另一方面，接 触面温度的改变还会显著改变两接触材料的摩擦特性，比如摩擦系数和磨损系数,除此之 夕卜，还会改变材料的本构特性，比如屈服极限和硬度等。因此，接触面的摩擦热（温度场）、 应力（接触压力场）和磨损（磨损量）之间是相互耦合的。
[0003] 这种耦合行为难以求得解析解。在采用数值方法进行求解数值解时，如果采用完 全耦合分析方法，每个时间增量步都考虑摩擦热-应力-磨损三者之间的双向相互耦合效 应，则计算效率极低，无法求解较长时间的摩擦热-应力-磨损耦合问题，但如果采用顺序 耦合方法，即先分析温度场，再进行应力-磨损耦合分析，忽略应力和磨损对热传导单方向 的耦合项，则牺牲了计算精度，仅适用于摩擦热-应力-磨损的弱耦合问题。
【发明内容】
[0004] 针对上述问题，本发明的目的是提供一种计算精度和计算效率较高，且适用于强 耦合问题的干滑动摩擦热-应力-磨损分步耦合的模拟方法。
[0005] 为实现上述目的，本发明采取以下技术方案：一种干滑动摩擦热-应力-磨损分步 耦合的模拟方法，包括以下步骤：1)时域分步：将摩擦热-应力-磨损耦合过程分步，初始 化步数计数器i = 1和仿真总时间计时器t = 0,计算初始接触压力Pd = Pi ;2)进入第i步 热-应力-磨损顺序耦合，其包括以下步骤：（1)热-应力耦合分析初始化：预设第i步耦 合过程的总时间步AT，增量步的步数为N，AT为耦合过程的总时间T'与增量步的步数N 之比，初始化第i步耦合过程的当前时间增量A & = 0和增量步计数器j = 1 ;(2)热传导 分析：在已有的瞬态热传导分析程序ABAQUS中依次完成N个增量步的热传导分析，每个增 量步对应输出一个温度场，用于应力-磨损分析；（3)应力-磨损分析初始化：开始进入应 力-磨损耦合分析之前，再次初始化增量步计数器j = 1，将步骤2)结束时的状态作为本步 应力-磨损分析的初始状态，并记录初始接触压力P(l ; (4)应力分析：进入第j个增量步的 应力分析，设置有限元模型应力分析的单元类型，定义边界条件、载荷和接触条件，将热传 导分析的温度场作为热载荷施加在有限元模型中，应用瞬态热传导分析程序ABAQUS对有 限元模型进行应力分析，获得并输出接触压力场；（5)确定磨损量：根据热传导分析得到的 温度场、应力分析得到的接触压力场以及两接触面接触节点的相对滑移速率，计算两接触 面接触节点的磨损量，并且确定其空间方向；(6)接触压力偏差计算：计算当前增量步的接 触压力P」与初始接触压力p Q之间的相对偏差e (p」，pQ) ; (7)更新网格：根据两接触面接触节 点的磨损量和方向，修正该节点的位移，更新有限元模型；(8)判断相对偏差dpppj是否 超过容差：考察接触压力的相对偏差e(Pp Ptl)是否达到容许值TOL，如果偏差没有超过容许 值TOL，转到步骤（4)，进入下一个增量步的应力分析；否则，抛出异常，终止本步耦合过程， 并更新本步耦合过程的当前时间增量△、= j At，进入步骤3);如果所有增量步都没有抛 出异常，则在最后一个增量步结束后更新总时间增量A、= NAt，进入步骤3) ;3)更新总 时间t :根据第i步稱合过程的当前时间增量Ah更新总时间t = t+Ati,判断当前总时 间t和耦合过程总时间T'的大小，如果t < T'，则更新计数器i = i+1和初始接触压力 P〇=Pp进入步骤2)，开始下一步摩擦热-应力-磨损耦合分析过程，否则，模拟过程结束。 [0006] 所述步骤2)中，瞬态热传导分析的生热热流边界由摩擦功率引起，其摩擦热流q 表达式为：V AW.__'.，式中，q为接触节点的摩擦热流；U为接触节点的摩擦系数；Pci为接 触面的初始接触压力；f为两接触表面在接触节点的相对滑动速率；摩擦系数U为接触面 温度T和初始接触压力P(l等因素的函数4 y (T，P(I，…）；3>根据两构件的具体结构和 相对运动求取；摩擦热流q通过热流用户子程序DFLUX施加在其中一构件的底面上，两接 触面接触节点温度T、接触节点到构件边沿的距离「均从计算结果文件中读取，用户子程序 DFLUX包括以下步骤：①从瞬态热传导分析程序ABAQUS中跳转至用户子程序DFLUX;②初 始化常量和系数；③根据两接触面接触节点坐标计算回转半径？，根据两接触面接触节点 温度T和平均接触压力f计算摩擦系数y ;④根据回转半径？和摩擦系数y计算两接触 面接触节点处的摩擦热流孑；⑤返回瞬态热传导分析程序ABAQUS。
[0007] 所述步骤2)中磨损量的计算采用广义非线性Archard磨损公式：△ h = kAs;式中，Ah为两接触面接触节点的磨损增量；k为广义Archard磨损公式的 磨损系数；P为两接触面接触节点的接触压力；△ s为两接触表面在接触节点的相对滑动 位移增量；磨损系数k是温度T、接触压力p和接触节点的相对滑动速率#等因素的函数：
[0008] 所述步骤2)中磨损量的空间方向的确定方法如下：两接触面接触节点分为两类， 一类是位于接触面边缘的节点，简称边界节点；另一类是位于接触面内部的节点，简称非边 界节点，对于非边界节点的磨损量方向由节点法方向确定，节点法方向直接从有限元分析 结果文件中读取。
[0009] 所述步骤2)中磨损量的确定在瞬态热传导分析程序ABAQUS中通过网格移动用户 子程序UMESHM0TI0N实现，两接触面接触节点的温度T、接触压力p、节点到构件边缘的距离 I*和增量步时间长度At均从计算结果文件中读取，用户子程序UMESHM0TI0N包括以下步 骤：①初始化常量和系数；②访问计算结果文件，读取两接触面接触节点的接触压力P和温 度T ;③计算两接触面接触节点的磨损量，确定两接触面接触节点的磨损方向；④出两接触 面接触节点磨损信息，用于绘制磨损量云图。
[0010] 所述步骤2)中，接触节点位移的修正和有限元模型的更新采用任意拉格朗日欧 拉法实现；需要修正位移的接触节点为接触压力与相对切向滑移均不为零的接触节点。
[0011] 本发明由于采取以上技术方案，其具有以下优点：1、本发明由于是在顺序耦合 方法的基础上，当接触压力场的变化达到一定量级时，修正接触压力，重新计算摩擦热流， 然后继续开始顺序耦合分析过程，依次反复进行直至完成整个耦合过程，因此，这种摩擦 热-应力-磨损分步耦合的模拟方法兼顾了计算精度和计算效率，也适合于摩擦热-应 力-磨损的强耦合过程。2、本发明由于是基于现有的商用有限元软件平台，因此，设置简 单、容易上手且便于推广。3、本发明由于磨损计算的程序在商用有限元软件的应力分析算 法基础上二次开发完成，完全嵌入商用有限元软件之中，免去了数据传输和交互，因此，大 大提高了计算效率。综上所述，本方法具有非常广泛的适用范围，可以用于各类干滑动摩擦 的摩擦热-应力-磨损耦合问题的模拟。
【专利附图】
【附图说明】
[0012] 图1是本发明的总体流程示意图；
[0013] 图2是本发明的第i步摩擦热-应力-磨损耦合的计算框图；
[0014] 图3是本发明的用户子程序UMASFL的流程示意图；
[0015] 图4是本发明的用户子程序DFLUX的流程示意图；
[0016] 图5是本发明的用户子程序UMESHM0TI0N的流程示意图；
[0017] 图6是本发明的销-盘有限元模型示意图；
[0018] 图7是本发明的销接触节点示意图；
[0019] 图8是本发明中销-盘磨损试验的温度分布示意图；
[0020] 图9是本发明与现有完全耦合方法和顺序耦合方法得到的销-盘磨损试验的接触 面温度对比示意图，其中，实线表示节点N101501采用完全耦合方法得到的销-盘磨损试验 的接触面温度，虚线表示节点N101501采用顺序耦合方法得到的销-盘磨损试验的接触面 温度，点划线表示节点N101501采用分步耦合方法得到的销-盘磨损试验的接触面温度；
[0021] 图10是本发明与现有完全耦合方法和顺序耦合方法得到的销的接触压力对比示 意图，其中，实线表示节点N101501采用完全耦合方法得到的销的接触压力，虚线表示节点 N101501采用顺序耦合方法得到的销的接触压力，点划线表示节点N101501采用分步耦合 方法得到的销的接触压力；
[0022] 图11是本发明与现有完全耦合方法和顺序耦合方法得到的销的磨损量对比示 意图，其中，实线表示节点N101501采用完全耦合方法得到的销的磨损量，虚线表示节点 N101501采用顺序耦合方法得到的销的磨损量，点划线表示节点N101501采用分步耦合方 法得到的销的磨损量。
【具体实施方式】
[0023] 下面结合附图和实施例对本发明进行详细的描述。
[0024] 本发明提供一种干滑动摩擦热-应力-磨损分步耦合的模拟方法，实施本发明方 法之前，首先建立干滑动摩擦的有限元模型，有限元模型包括两个具有相对干滑动摩擦的 构件，约束不能使两构件产生相对干滑动摩擦的自由度，初始化两构件之间的作用载荷，初 始化两构件的材料参数，材料参数包括不同温度对应的比热容、导热系数、线膨胀系数、杨 氏模量、密度和泊松比。本发明方法包括以下步骤（如图1、图2所示）：
[0025] 1)时域分步：将摩擦热-应力-磨损耦合过程分步，初始化步数计数器i = 1和 仿真总时间计时器t = 0,计算初始接触压力Pd = Pi。
[0026] 2)进入第i步热-应力-磨损顺序耦合，其包括以下步骤（如图2所示）：
[0027] (1)热-应力耦合分析初始化：预设第i步耦合过程的总时间A T，增量步的步数 为N，AT为耦合过程的总时间T'与增量步的步数N之比，初始化第i步耦合过程的当前 时间增量A & = 0和增量步计数器j = 1。
[0028] (2)热传导分析：在已有的瞬态热传导分析程序ABAQUS中依次完成N个增量步的 热传导分析，每个增量步对应输出一个温度场，用于应力-磨损分析；热传导分析步骤为：
[0029] ①初始化有限元模型的温度，设置有限元模型中两构件的单元类型、生热热流边 界和散热边界；
[0030] ②对每个增量步进行瞬态热传导分析；
[0031] ③获得并输出温度场。
[0032] 如果其中一单元类型允许其对应的构件的材料脱离网格运动，那么材料的运动利 用用户子程序UMASFL实现（如图3所示），以销-盘磨损实验为例，用户子程序UMASFL包 括以下步骤：
[0033] ①从瞬态热传导分析程序ABAQUS中跳转至用户子程序UMASFL ;
[0034] ②初始化构件的常量和系数，如角速度6、材料密度P等；
[0035] ③根据销-盘表面接触节点坐标计算回转半径？，如果不是旋转件，直接读取该节 点速度即可；对于销_盘实验，知道销-盘表面接触节点坐标和旋转角速度即可计算出节点 速度；
[0036] ④根据回转半径^计算销-盘表面接触节点处的质量流速______
[0037] ⑤返回瞬态热传导分析程序ABAQUS。
[0038] 瞬态热传导分析的生热热流边界由摩擦功率引起，其摩擦热流q表达式为：
[0039] q=up{j, ( 1 )
[0040] 式中，q为接触节点的摩擦热流；U为接触节点的摩擦系数；为接触面的初始接 触压力；f为两接触表面在接触节点的相对滑动速率。
[0041] 摩擦系数y为接触面温度T和初始接触压力P(l等因素的函数：
[0042] U = U (T,p〇, --?) (2)
[0043] #根据两构件的具体结构和相对运动求取。
[0044] 摩擦热流q通过热流用户子程序DFLUX施加在其中一构件的底面上，节点温度T、 回转半径?均可以从计算结果文件中读取，以销-盘实验为例，用户子程序DFLUX包括以下 步骤（如图4所示）：
[0045] ①从瞬态热传导分析程序ABAQUS中跳转至用户子程序DFLUX ;
[0046] ②初始化常量和系数，如常量和平均接触压力f等；
[0047] ③根据销-盘表面接触节点坐标计算回转半径M艮据销_盘表面接触节点温度T 和平均接触压力]H十算摩擦系数U ;
[0048] ④根据回转半径卩和摩擦系数U计算销-盘表面接触节点处的摩擦热流f
[0049] ⑤返回瞬态热传导分析程序ABAQUS。
[0050] (3)应力-磨损分析初始化：开始进入应力-磨损耦合分析之前，再次初始化增量 步计数器j = 1，将步骤2)结束时的状态作为本步应力-磨损分析的初始状态，并记录初始 接触压力P0。
[0051] (4)应力分析：进入第j个增量步的应力分析，设置有限元模型应力分析的单元 类型，定义边界条件、载荷和接触条件，将热传导分析的温度场作为热载荷施加在有限元模 型中，应用瞬态热传导分析程序ABAQUS对有限元模型进行应力分析，获得并输出接触压力 场。
[0052] (5)确定磨损量：根据热传导分析得到的温度场、应力分析得到的接触压力场以 及两接触面接触节点的相对滑移速率，计算两接触面接触节点的磨损量，并且确定其空间 方向。其中，磨损量的计算采用广义非线性Archard磨损公式：
[0053] A h =k? p ? A s (3)
[0054] 式中，Ah为两接触面接触节点的磨损增量；k为广义Archard磨损公式的磨损 系数；P为两接触面接触节点的接触压力；△ s为两接触表面在接触节点的相对滑动位移增 量。
[0055] 磨损系数k是温度T、接触压力p和接触节点的相对滑动速率?等因素的函数：
[0056] K= k{T， p，￥,?) ( 1)
[0057] 两接触面接触节点磨损量的空间方向的确定方法如下：两接触面接触节点可以分 为两类，一类是位于接触面边缘的节点，简称边界节点；另一类是位于接触面内部的节点， 简称非边界节点。对于非边界节点的磨损量方向可以由节点法方向确定，节点法方向可以 直接从有限元分析结果文件中读取。而对于边界节点，由于不存在法方向，因此，采用边界 节点与其厚度方向对应点的连线方向作为磨损量方向。
[0058] 磨损量的确定在瞬态热传导分析程序ABAQUS中通过网格移动用户子程序 UMESHM0TI0N实现，两接触面接触节点的温度T、接触压力p、节点到构件边缘的距离I* 和增量步时间长度At均可以从计算结果文件中读取，以销-盘实验为例，用户子程序 UMESHM0TI0N包括以下步骤（如图5所示）：
[0059] ①初始化常量和系数，如常量和角速度5等；
[0060] ②访问计算结果文件，读取销-盘表面接触节点的接触压力p和温度T ;
[0061] ③计算销-盘表面接触节点的磨损量，确定销-盘表面接触节点的磨损方向；
[0062] ④输出销-盘表面接触节点磨损信息，用于绘制磨损量云图；
[0063] (6)接触压力偏差计算：计算当前增量步的接触压力Pj与初始接触压力之间的 相对偏差e (Pj, pQ)。
[0064] (7)更新网格：根据两接触面接触节点的磨损量和方向，修正接触压力与相对切 向滑移均不为零的两接触面接触节点的位移，更新有限元模型，更新网格是通过更新两 接触面接触节点位移，然后重画磨损区域网格实现，本发明采用了任意拉格朗日欧拉法 (ALE)。ALE需要设置ALE作用域和ALE约束条件，其中，ALE的作用域是所有包含需要修正 位移的接触节点的单元，ALE的约束范围为所有需要修正位移的接触节点，ALE的约束类型 为位移约束或者速度约束。
[0065] 以速度约束为例，若接触节点i的磨损速率为/丨;，其在全局坐标系下的方向向量 为| ,接触节点i的局部坐标系三个轴的归一化方向矢量分别为，则该接触节点的 ALE速度约七欠R/I以郃來系屮的分R丨"，Vli，vzi!T,山F面公式计算：
[0066]
【权利要求】
1. 一种干滑动摩擦热-应力-磨损分步耦合的模拟方法，包括以下步骤： 1) 时域分步：将摩擦热-应力-磨损耦合过程分步，初始化步数计数器i= 1和仿真 总时间计时器t= 0,计算初始接触压力Pci =Pi ; 2) 进入第i步热-应力-磨损顺序耦合，其包括以下步骤： (1) 热-应力耦合分析初始化：预设第i步耦合过程的总时间步AT，增量步的步数为 Ν，ΛΤ为耦合过程的总时间T'与增量步的步数N之比，初始化第i步耦合过程的当前时间 增量Λti = 0和增量步计数器j= 1 ; (2) 热传导分析：在已有的瞬态热传导分析程序ABAQUS中依次完成N个增量步的热传 导分析，每个增量步对应输出一个温度场，用于应力-磨损分析； (3) 应力-磨损分析初始化：开始进入应力-磨损耦合分析之前，再次初始化增量步计 数器j= 1，将步骤2)结束时的状态作为本步应力-磨损分析的初始状态，并记录初始接触 压力P0 ; (4) 应力分析：进入第j个增量步的应力分析，设置有限元模型应力分析的单元类型， 定义边界条件、载荷和接触条件，将热传导分析的温度场作为热载荷施加在有限元模型中， 应用瞬态热传导分析程序ABAQUS对有限元模型进行应力分析，获得并输出接触压力场； (5) 确定磨损量：根据热传导分析得到的温度场、应力分析得到的接触压力场以及两 接触面接触节点的相对滑移速率，计算两接触面接触节点的磨损量，并且确定其空间方 向； (6) 接触压力偏差计算：计算当前增量步的接触压力ρ」与初始接触压力Ptl之间的相对 偏差e(ρ』，p0); (7) 更新网格：根据两接触面接触节点的磨损量和方向，修正该节点的位移，更新有限 元模型； (8) 判断相对偏差e(Pj，Ptl)是否超过容差：考察接触压力的相对偏差e(Pj，Ptl)是否 达到容许值TOL，如果偏差没有超过容许值TOL，转到步骤（4)，进入下一个增量步的应力 分析；否则，抛出异常，终止本步耦合过程，并更新本步耦合过程的当前时间增量Ati = jΛt，进入步骤3);如果所有增量步都没有抛出异常，则在最后一个增量步结束后更新总 时间增量Ati =NAt，进入步骤3); 3) 更新总时间t:根据第i步稱合过程的当前时间增量Ati更新总时间t=t+Ati, 判断当前总时间t和耦合过程总时间Ti的大小，如果t<Ti，则更新计数器i=i+l和 初始接触压力Ptl =Pp进入步骤2)，开始下一步摩擦热-应力-磨损耦合分析过程，否则， 模拟过程结束。
2. 如权利要求1所述的一种干滑动摩擦热-应力-磨损分步耦合的模拟方法，其特征 在于：所述步骤2)中，瞬态热传导分析的生热热流边界由摩擦功率引起，其摩擦热流q表达 式为： <//IpJ' 式中，q为接触节点的摩擦热流；μ为接触节点的摩擦系数；Ptl为接触面的初始接触压 力；彳为两接触表面在接触节点的相对滑动速率； 摩擦系数μ为接触面温度T和初始接触压力Ptl等因素的函数： μ = μ (τ, ρ〇,···)； f根据两构件的具体结构和相对运动求取； 摩擦热流q通过热流用户子程序DFLUX施加在其中一构件的底面上，两接触面接触节 点温度T、接触节点到构件边沿的距离/均从计算结果文件中读取，用户子程序DFLUX包括 以下步骤： ① 从瞬态热传导分析程序ABAQUS中跳转至用户子程序DFLUX; ② 初始化常量和系数； ③ 根据两接触面接触节点坐标计算回转半径?,根据两接触面接触节点温度T和平均 接触压力#计算摩擦系数μ; ④ 根据回转半径F和摩擦系数μ计算两接触面接触节点处的摩擦热流f; ⑤ 返回瞬态热传导分析程序ABAQUS。
3. 如权利要求1所述的一种干滑动摩擦热-应力-磨损分步耦合的模拟方法，其特征 在于：所述步骤2)中磨损量的计算采用广义非线性Archard磨损公式： Ah=κ·ρ·As； 式中，Ah为两接触面接触节点的磨损增量；κ为广义Archard磨损公式的磨损系数；P为两接触面接触节点的接触压力；△s为两接触表面在接触节点的相对滑动位移增量； 磨损系数κ是温度T、接触压力p和接触节点的相对滑动速率#等因素的函数： /; = /·;ν
4. 如权利要求2所述的一种干滑动摩擦热-应力-磨损分步耦合的模拟方法，其特征 在于：所述步骤2)中磨损量的计算采用广义非线性Archard磨损公式： Ah=κ·ρ·As； 式中，Ah为两接触面接触节点的磨损增量；κ为广义Archard磨损公式的磨损系数；P为两接触面接触节点的接触压力；△s为两接触表面在接触节点的相对滑动位移增量； 磨损系数κ是温度T、接触压力p和接触节点的相对滑动速率等因素的函数：
5. 如权利要求1或2或3或4所述的一种干滑动摩擦热-应力-磨损分步耦合的模 拟方法，其特征在于：所述步骤2)中磨损量的空间方向的确定方法如下：两接触面接触节 点分为两类，一类是位于接触面边缘的节点，简称边界节点；另一类是位于接触面内部的节 点，简称非边界节点，对于非边界节点的磨损量方向由节点法方向确定，节点法方向直接从 有限元分析结果文件中读取。
6. 如权利要求1或2或3或4所述的一种干滑动摩擦热-应力-磨损分步耦合的模拟 方法，其特征在于：所述步骤2)中磨损量的确定在瞬态热传导分析程序ABAQUS中通过网格 移动用户子程序UMESHM0TI0N实现，两接触面接触节点的温度T、接触压力p、节点到构件边 缘的距离/"和增量步时间长度Λt均从计算结果文件中读取，用户子程序UMESHM0TI0N包括 以下步骤： ① 初始化常量和系数； ② 访问计算结果文件，读取两接触面接触节点的接触压力P和温度T; ③ 计算两接触面接触节点的磨损量，确定两接触面接触节点的磨损方向； ④ 出两接触面接触节点磨损信息，用于绘制磨损量云图。
7. 如权利要求5所述的一种干滑动摩擦热-应力-磨损分步耦合的模拟方法，其特征 在于：所述步骤2)中磨损量的确定在瞬态热传导分析程序ABAQUS中通过网格移动用户子 程序UMESHMOTION实现，两接触面接触节点的温度T、接触压力p、节点到构件边缘的距离 Γ和增量步时间长度At均从计算结果文件中读取，用户子程序UMESHMOTION包括以下步 骤： ① 初始化常量和系数； ② 访问计算结果文件，读取两接触面接触节点的接触压力P和温度T; ③ 计算两接触面接触节点的磨损量，确定两接触面接触节点的磨损方向； ④ 出两接触面接触节点磨损信息，用于绘制磨损量云图。
8. 如权利要求1或2或3或4或7所述的一种干滑动摩擦热-应力-磨损分步耦合的 模拟方法，其特征在于：所述步骤2)中，接触节点位移的修正和有限元模型的更新采用任 意拉格朗日欧拉法实现；需要修正位移的接触节点为接触压力与相对切向滑移均不为零的 接触节点。
9. 如权利要求5所述的一种干滑动摩擦热-应力-磨损分步耦合的模拟方法，其特征 在于：所述步骤2)中，接触节点位移的修正和有限元模型的更新采用任意拉格朗日欧拉法 实现；需要修正位移的接触节点为接触压力与相对切向滑移均不为零的接触节点。
10. 如权利要求6所述的一种干滑动摩擦热-应力-磨损分步耦合的模拟方法，其特征 在于：所述步骤2)中，接触节点位移的修正和有限元模型的更新采用任意拉格朗日欧拉法 实现；需要修正位移的接触节点为接触压力与相对切向滑移均不为零的接触节点。
【文档编号】G06F17/50GK104239627SQ201410458032
【公开日】2014年12月24日 申请日期:2014年9月10日 优先权日:2014年9月10日
【发明者】桂良进, 张方宇, 范子杰 申请人:清华大学