基于ANSYS的轮对过盈配合微动分析
基于ANSYS的轨道车辆车轮与轨道接触分析及优化设计
基于ANSYS的轨道车辆车轮与轨道接触分析及优化设计摘要：针对特殊车辆中的轨道车辆车轮与轨道相互接触的问题，结合赫兹接触计算理论[1-3]，改进计算因子。
建立三维模型，并且对模型进行有限元分析[4-6]，得到位移及应力云图。
建立数学优化模型，对于结果进行优化。
结合实例加以说明。
关键词：接触分析；优化模型；优化分析Abstract：To the wheel and the rail track vehicle special vehicle in mutual contact problems，combining with the theory of Hertz contact，improved calculation factor. The three-dimensional model is established，and the finite element analysis on the model，get the displacement and stress cloud chart. The mathematical optimization model is established to optimize the result. Combined with examples to illustrate.Keywords：Contact analysis Optimization model Optimization analysis1.概述本文主要针对轨道车辆车轮与轨道之间的接触问题建立模型，确定接触问题为刚性和柔性接触。
对于高度非线性问题，运用ANSYS分析工具，在确定接触区域及接触时间前提下，仿真分析得到位移及应力的云图，运用改进接触分析计算因子，通过建立优化模型，对结果进行优化处理。
2.构建模型简图及确定相关参数通过调节优化参数的数值后，在施加载荷相同的情况下，车轮直径一定的情况下，车轮和轨道截面的厚度在约束范围内越大，应力造成的破坏范围越小，最大应力的数值先减小，后增大，位移所显示的形变量越小。
基于Ansys Workbench的行星架组件过盈配合仿真研究及结构改善
基于Ansys Workbench 的行星架组件过盈配合仿真研究及结构改善徐 鹏1 吕小波2 胡吉全1 李 波1 杨艳芳11.武汉理工大学物流工程学院 武汉 4300632.湖北科峰传动设备有限公司 黄冈 438000摘 要：针对某典型小型行星减速器输出轴组件，对其销轴与行星架的过盈配合展开研究。
首先根据理论计算得到过盈量的初始范围；再利用Ansys Workbench 软件模拟分析设计参数对过盈装配的影响，其中重点分析了过盈量和摩擦系数对接触压力、等效应力和压装曲线的影响关系。
研究发现过盈量对过盈配合的影响最大，摩擦系数对接触压力的影响可以忽略不计，在此基础上根据仿真结果对理论计算的初始过盈量范围进行优化，得到最佳过盈量范围为0.004～0.008 mm；最后将销轴结构进行改善，结果表明所提出的阶梯销轴不会改变联接的可靠性，且能改善结构应力情况，提升压装质量，为小型行星减速器输出轴组件的设计提供参考。
关键词：行星减速器；过盈配合；仿真优化；压装质量中图分类号：TH132.46 文献标志码：A 文章编号：1001-0785（2020）14-0032-07Abstract: As for components of the output shaft of a typical small planetary reducer, the interference fit between its pin shaft and planet carrier was studied. Firstly, the initial magnitude of magnitude of interference was calculated according to theories. Then, Ansys Workbench software was used to simulate and analyze the influence of design parameters on interference fit, among which the influence of magnitude of interference and friction coefficient on contact pressure, equivalent stress and press-fitting curve was emphatically analyzed. It was found that the magnitude of interference affected the interference fit to the greatest degree, and the influence of friction coefficient on contact pressure can be ignored. Based on the above, the initial magnitude of interference calculated theoretically was optimized according to simulation results, and the optimal magnitude of interference was 0.004 to 0.008 mm; finally, the pin shaft structure is improved, and results showed that the proposed stepped pin shaft will not change the reliability of the connection, and can improve the structural stress and the quality of press fitting, which can provide a reference for the design of output shaft components of small planetary reducerKeywords: planetary reducer; interference fit; simulation optimization; press-fitting quality0 引言减速器是机械零部件中的重要成分，在我国的装备与制造业中发挥着重要作用。
过盈配合的有限元分析
过盈配合的有限元分析工程力学系张晨朝20803001过盈配合的有限元分析摘要: 在工程应用中,利用接触有限元法建立了内轴与外套过盈配合的有限元力学模型来判断结构设计是否符合要求。
针对内轴和外套的过盈配合状态，采用大型通用有限元ANSYS 软件对组合模具进行了有限元分析, 得出了内轴与外套在过盈配合状态下的应力分布规律及接触面压力分布状况, 找到了应力集中位置和大小。
结论表明结构配合尺寸设计没有使结构产生变形, 该结构完全符合产品的设计要求。
关键词: 过盈配合; ANSYSAbstract: In the project application, in order to judge whether the structural design meets the requirement, the finite element and mechanical model of the interference joint between inside lining and outside wrap is established by used contact -finite- element methods. Aimed at condition of the interference joint between inside lining and outside wrap, we carry on the finite element analysis based on ANSYS and attain the stress distribution in interference joint; the pressure distribution in contact face and the location and the size of stress concentration. It is concluded that the structure interference joint size of combined die do not make the mold have distortion and the combined die completely meets the product design requirement.Key words: interference joint; ANSYS1 引言过盈配合[1]是机械工业中一种常见的零部件组装方式,齿轮、轴承以及火车车轮等与其装配轴之间的配合大多采用过盈配合。
Ansys_过盈配合接触设置
例 3.一个简单的轴和带孔圆盘的过盈配合的实例。 圆盘的基本尺寸为： 内径 Rpin=35mm（原书中此值为 34mm）,外径 Rpout=100,盘高 Hp=25mm； 轴的基本尺寸为： 内径 Rain=25mm，外径 Raout=35mm，轴长 La=150mm。
在 ANSYS 中正确地模拟过盈配合
过盈配合在机械产品的装配中使用相当普遍,譬如轴与轴承,轴与轴瓦,汽车 的制动盘等,都是通过一定的过盈量来使两个装配部件紧密连接起来。
下面讨论如何在 ANSYS 中正确地模拟过盈配合。过盈配合在有限元分析中 是一种典型的非线性接触行为。在有限元分析中设定了接触，从本质上来讲就是 对相互接触的两个部件施加了某种约束，不同的接触算法对于接触约束的处理方 法有所不同。接触约束的理论算法的选择，在 ANSYS 中是通过设置 contact 单 元的 KEOPT(2)选项来实现的。在 ANSYS 中目前主要有 5 种接触约束算法：
图 2：一个轴类零件和孔类零件的过盈配合的截面图
在几何上，图 2 的轴和孔有一定的过盈配合量，其大小等于两个圆的半径之 差，我们的本意是想用这个几何位置上的过盈量来模拟过盈配合。
不幸的是，两个部件划分网格之后，实际的过盈量应该为单元之间的距离， 即图上中靠得比较近的两条线段之间的距离，显然，这个距离不再等于我们预先 设置的过盈量了。更何况，上面这个图还是两个部件的网格对对得比较整齐的情 况，如果网格对的不整齐，过盈量就和我们预设的差的更远了。对于过盈配合来 讲，过盈量的数值变化对于过盈产生的应力的影响是很大的。
NSLA,S,1 ESLN,S,0 ESURF ALLSEL ESEL,ALL ESEL,S,TYPE,,2 ESEL,A,TYPE,,3 ESEL,R,REAL,,3 /PSYMB,ESYS,1 /PNUM,TYPE,1 /NUM,1 EPLOT ESEL,ALL /COM, CONTACT PAIR CREATION - END
基于ANSYS Workbench对圆柱面过盈配合接触应力的研究
ｗ ｉ ｄ ｔ ｈ． ｔ ｈ ｅ ｉ ｎ ｌ ｆ ｕ ｅ ｎ ｃ ｅ ｆ ａ ｃ ｔ ｏ ｒ ｓ ｏ ｆ ａ ｘ ｉ ａ ｌ ｔ ｅ ｎ ｓ ｉ ｏ ｎ ａ ｎ ｄ ｓ ｔ ｒ ｅ ｓ ｓ ．Ａｎ ａ ｌ ｙ ｓ ｉ ｓ ｏ ｆ ｅ ｘ ｐ ｅ ｒ ｉ ｍｅ ｎ ｔ ａ ｌ ｒ ｅ ｓ ｕ ｌ ｔ ｓ ｔ ｈ ｅ ｆ ｏ ｌ ｌ ｏ ｗｉ ｎ ｇ ｃ ｏ ｎ ｃ ｌ ｕ ｓ ｉ ｏ ｎ ｓ ：ｓ ｈ ａ ｔ ｆ ａ ｎ ｄ
ｐ ｒ ｏ ｃ ｅ ｓ ｓ ， ｂ ｙ ｃ ｈ ａ ｎ ｇ ｉ ｎ ｇ ｔ ｈ ｅ ｓ ｉ ｎ ｇ ｌ ｅ ｆ ａ ｃ ｔ ｏ ｒ ｅ ｘ ｐ ｅ ｒ ｉ ｍｅ ｎ ｔ ｍｅ ｔ ｈ ｏ ｄ， ａ ｎ ａ ｌ ｙ ｓ ｉ ｓ ｏ ｆ ｉ ｎ ｔ ｅ ｆ ｒ ｅ ｒ ｅ ｎ ｃ ｅ ｉ ｆ ｔ ｏ ｆ ｃ ｙ ｌ ｉ ｎ ｄ ｅ ｒ ｉ ｓ ｃ ｏ ｎ ｔ ａ ｃ ｔ ｓ ｔ ｒ ｅ ｓ ｓ ｄ ｉ ｓ ｔ ｒ ｉ ｂ ｕ ｔ ｉ ｏ ｎ
Ａｂｓ ｔ ｒ ａｃ ｔ ： Ｕｓ ｉ ｎ ｇ ｆ ｉ ｎ ｉ ｔ ｅ ｅ ｌ ｅ ｍｅ ｎ ｔ ａ ｎａ ｌ ｙ ｓ ｉ ｓ ｓ ｏ ｆ ｔ ｗａ ｒ ｅ ＡＮＳＹＳ Ｗｏ ｒ ｋ ｂｅ ｎ ｃ ｈ， ｔ ｈ ｅ ｓ ｉ ｍｕ ｌ ａ ｔ ｉ ｏ ｎ ｏ ｆ ｔ ｈ ｅ ｈｏ ｌ ｌ ｏ ｗ ｓ ｈａ ｆ ｔ ａ ｎｄ ｃ ｏ ｎ ｔ ａ ｉ ｎｉ ｎ ｇ ａ
关键词 ： 过盈配合 ； 有 限元 ； 影 响 因素 ； 等效应力
中图 分 类 号 ： Ｔ Ｈ １ ４ ０
Ｒｅ ｓ ｅ ａ ｒ ｃ ｈ ０ ｎ ｃ ｏ ｎｔ ａ ｃ ｔ ｓ ｔ ｒ ｅ ｓ ｓ ０ ｆ ｃ ｙ ｌ ｉ ｎ ｄｅ ｒ ｉ ｎ ｔ ｅ ｒ ｆ ｅ ｒ ｅ ｎ ｃ ｅ ｉｔ ｆ ｂａ ｓ ｅ ｄ
基于ansys的过盈配合接触应力分析
基于ansys的过盈配合接触应力分析摘要介绍了基于ansys的接触分析步骤，并通过ansys软件，将对一个盘轴紧配合结构进行接触分析，来说明接触分析的有限元计算方法。
关键词ansys 过盈配合接触分析引言在工程结构中，经常会遇到大量的接触问题。
火车车轮与钢轨之间，齿轮的啮合是典型的接触问题。
接触问题是一种高度非线性行为，需要较大的计算资源，为了进行实为有效的计算，理解问题的特性和建立合理的模型是很重要的。
接触问题存在两个较大的难点：其一，在你求解问题之前，你不知道接触区域，表面之间是接触或分开是未知的，突然变化的，这随载荷、材料、边界条件和其它因素而定；其二，大多的接触问题需要计算摩擦，有几种摩擦和模型供你挑选，它们都是非线性的，摩擦使问题的收敛性变得困难。
本文以ansys软件为工具，以某转子中轴和盘的连接为例，分析轴和盘的过盈配合的接触应力。
1.面面接触分析的步骤：在涉及到两个边界的接触问题中，很自然把一个边界作为“目标”面而把另一个作为“接触”面，对刚体—柔体的接触，“目标”面总是刚性的，“接触”面总是柔性面，这种两个面合起来叫作“接触对”。
使用Targe169和Conta171或Conta172来定义2D接触对，使用Targe170和Conta173或Conta174来定义3D接触对，程序通过相同的实常数号来识别“接触对”。
在接触问题中，两个相互接触的物体必须满足边界不穿透的约束条件，施加边界不穿透约束的方法主要有罚函数算法和扩增的拉格朗日算法。
罚函数算法是在总势能泛函中加入惩罚项，来近似满足接触约束条件。
从物理意义上讲，罚函数法相当于在接触边界上加入线弹簧以防止接触面之间的相互渗透，而罚函数因子相当于弹簧的刚度系数。
罚函数法的优点在于不增加系统未知数总数，可保持刚度矩阵的对称性，提高了求解效率，但罚函数因子的取值对计算结果的精度影响很大，必须根据渗透情况对其进行多次调整。
扩增的拉格朗日算法是为了找到精确的拉格朗日乘子而对罚函数修正项进行反复迭代，与罚函数的方法相比，拉格朗日方法不易引起病态条件，对接触刚度的灵敏度较小，然而，在有些分析中，扩增的拉格朗日方法可能需要更多的迭代，特别是在变形后网格变得太扭曲时。
Ansysworkbench参数化过盈配合模拟计算确定所需过盈量校核计算
Ansysworkbench参数化过盈配合模拟计算确定所需过盈量在机构设计中，常常需要校核设计的过盈量能否满足设计需求，一些规则简单的构件能够通过手工计算校核，但一些复杂的零件会增加计算的难度和误差，这时候可以利用有限元软件进行计算确认。
1、打开软件，建立所需模型（可直接在CAD软件中建立导入）2、双击Static Structural,右键Geometry---Import Geometry---Browse导入建好的模型3、双击进入Gometry，检查模型，退出；再双击Model,进入载荷和边界条件设置材料默认为structural steel,设置网格大小，点击生成网格设置配合接触面为摩擦接触，摩擦系数设置为0.2（可根据实际材料设定摩擦系数）4、右键Frictional-1,插入commands命令右侧命令行输入keyopt,cid,9,6 （消除模型及网格划分造成的过盈量误差）5、插入contact Tool，将插入命令前后的接触信息进行对比插入前，过盈值误差3.4634e-5插入命令后，过盈值7.2172e-15,几乎可以忽略为06、增加所需的过盈量值，左键点击Frictional----offset---设置为0.02，并勾选上前面的框，出现一个P，进行参数化7、设置边界条件及载荷，端面A设置传递扭矩Moment为50N.m,端面B设置fixed，下图所示：将Moment勾选，进行参数化8、右键solution，分别插入:---Equivalent Stress---勾选Maximum---contact Tool---sliding Distance---勾选Maximum---Moment Reaction---勾选Z Axis退出，进入操作平台9、双击parameter Set左侧出现下图，分为输入和输出参数右侧出现下图列表栏根据设计需求，可分别增加参数，本例增加扭矩200N.m及过盈量0.025，进行计算，得下图结果判断依据：D栏等效应力能否满足材料强度要求；E栏两配合面相对滑动值是否过大，导致传递误差过大：F栏输出端扭矩值是否和输入端一致或者接近：。
基于ANSYS Workbench对圆柱面过盈配合接触应力的研究
基于ANSYS Workbench对圆柱面过盈配合接触应力的研究李亚洲;汤易升;陈炜;张西正【摘要】利用有限元分析软件ANSYS Workbench,模拟空心轴与包容件的过盈配合过程,通过改变单一因素的实验方法,分析圆柱面过盈配合中接触应力分布情况和各个因素的关系.本文分析了过盈量、包容件外径、接触宽度、轴向拉力和压力的影响因素.对实验结果进行分析得出如下结论:轴和包容件的接触区域的最高等效应力区,受到边缘效应的影响,最高应力集中区出现在端部边缘处.包容件在靠近边缘的中间部位出现低应力区域.空心轴表面的应力值呈现出从一端到另端逐渐增加的现象.%Using finite element analysis software ANSYS Workbench,the simulation of the hollow shaft and containing a process,by changing the single factor experiment method,analysis of interference fit of cylinder is contact stress distribution and the relationship between the factors.The interference are analyzed in this paper,the quantity,inclusive a diameter,contact width,the influence factors of axial tension and stress.Analysis of experimental results the following conclusions:shaft and tolerance of contact area of the maximum equivalent stress area,under the influence of edge effect,high stress,high in the end edge.Containing a appears near the edge of the middle low stress area.On the surface of the hollow shaft stress value present a growing phenomenon from one end to the other end.【期刊名称】《天津理工大学学报》【年(卷),期】2017(033)003【总页数】5页(P1-5)【关键词】过盈配合;有限元;影响因素;等效应力【作者】李亚洲;汤易升;陈炜;张西正【作者单位】天津理工大学机械工程学院天津市先进机电系统与智能控制重点实验室,天津300384;天津理工大学机械工程学院天津市先进机电系统与智能控制重点实验室,天津300384;天津理工大学机械工程学院天津市先进机电系统与智能控制重点实验室,天津300384;军事医学科学院卫生装备研究所,天津300161;天津理工大学机械工程学院天津市先进机电系统与智能控制重点实验室,天津300384;军事医学科学院卫生装备研究所,天津300161【正文语种】中文【中图分类】TH140过盈连接普遍采用的方法是热胀冷缩或强力压入将有过盈量的两个接触面装配到一起[1-6].在接触面上会产生很大的径向力.在过盈连接承受轴向力和扭矩时，接触区域将产生相应的摩擦力与力矩来抵御外部施加的力和力矩[7-8].过盈连接是一种半永久性连接.它有很多优点，结构简单，生产成本低，消除了焊接产生局部应力，零部件容易产生局部疲劳与断裂的缺点，也解决了部分金属材料和机械结构进行螺纹加工困难的问题，常用于薄壁件连接处（不易加工螺纹），具有良好的对中性，应力分布相对比较均匀，承载能力强，结构紧凑美观等优点[9].现阶段对过盈配合的研究，已经在很多方面进行了.有从过盈连接的接触面变形类型来研究的，将过盈面的的变形分为弹性变形和塑性变形两个阶段进行不同的研究[10].弹形是指在过盈连接件，拆卸后，原先的接触面能够恢复到原有的直径不发生永久性变形.塑性变形，接触表面的材料受到很大的力，超出其强度极限，导致接触面发生了永久性变形.另一种研究分类，将过盈连接的研究，根据接触面的形状不同，分为圆锥面和圆柱面[11]对圆柱面的研究相对比较多，主要研究领域集中在轴和轮毂配合方面，研究方法普遍采用弹性力学与材料力学的理论公式进行推导计算.应用传统的经典弹性力学，可以解决轴和包容件的应力特性问题.但随着工程研究的深入，人们需要更加直观和准确研究接触区域的非线性力学问题.一些学者提出了对过盈连接应力分布影响因素的研究.由于经典的弹性力学方法已经不再适用，无法满足对过盈配合区域的直观的定性分析.现阶段普遍采用的研究方法是将数值计算的方法引入到对过盈连接影响因素的研究中.随着有限元分析的方法的发展，很多学者利用有限元分析软件，对过盈配合面进行应力分析，能更加准确、直观、细致的分析出过盈配合的变形过程[12].本文利用有限元分析软件ANSYS Workbench，通过控制单一变量的方法，对过盈链接的影响因素进行了研究，得到每个单变量因素对接触区域应力分布的影响、应力值的影响规律.建立的有限元模型，材料选择为结构钢，空心轴总长为30 mm，外径为Φ48 mm，内径为Φ40 mm；轴套总长为Φ30 mm，基本配合外径为Φ40 mm，过盈量大小取（0.042 mm，0.059 mm，0.076 mm）内径为Φ32 mm轴和轴套均采用合金钢材，其材料常数为：弹性模量E=2.01×109Pa，泊松比ν=0.27，材料密度ρ=7 800 kg/m3.模型的建立采用中间格式倒入的方法，首先在Solid works 里建立空心轴与包容件的模型，利用中间格式导入到ANSYS Workbench软件中进行过盈配合的模拟装配.模拟装配过程如图1所示.利用Solid works建立图1a三维模型，将建立的三维模型转化为中间格式x-t的文件，导入到ANSYS Workbench中，后对网格属性，材料属性进行设置，进行网格的划分得到图1b.网格划分完成，对包容件施加约束和位移加载，方法，对包容件端部固定，同时对空心轴一端施加位移，完成的模拟装配过程.得到如图1c所示模拟装配后的应力分布图.2.1 过盈量变化对等效应力的影响的计算研究本次研究借助ANSYS Workbench软件，通过控制单一影响因素，观察结果变化的方法，研究过盈量大小，对过盈配合区域的接触应力分布的影响.本次实验，Φ40 mm孔与轴的配合采用基孔制，过盈量根据国标GB1801-1979选取H7/n7、H7/s6、H7/u6，三种过盈配合关系，通过计算得到最大的过盈量分别是0.042 mm，0.059 mm，0.076 mm.根据控制单一变量的原理，保持接触宽度，轴向力，包容件外径大小，都不变，单一改变过盈量，对三组实验对象进行模拟装配.如2图所示，图2a、2b、2c分别代表过盈量为0.042 mm、0.059 mm、0.072 mm包容件的等效应力分布图，图2中包容件的应力分布有明显的规律性，在靠近边缘的中部出现了低应力区域.图3a、3b、3c分别代表空心轴的应力分布情况，空心轴外表面接触区域应力分布，呈现出从一端到另一端逐渐增加的变化，最大应力集中到插入端.图4a所示为过盈量为0.042 mm、0.059 mm、0.072 mm的包容件，从未被约束端到被约束端的应力分布情况，从曲线图4a可知，不同过盈量下包容件表面的应力分布均发生了突变，出现低应力区，从三条曲线的总体趋势来看，接触区域的应力值大小与过盈量大小成正比，随过盈量的增加而增加.图4b代表不同过盈量下，轴的接触区域的等效应力分布曲线.从曲线中可以看出，轴的接触区域应力值，从一端到另一端不断增加，从不同过盈量的轴对比发现，轴的接触区域应力值随着过盈量的增大而升高.2.2 接触宽度的变化对应力集中情况影响的研究对接触宽度的变化对应力集中情况的研究实验，采用同样的方法，保持过盈量，包容件外径大小，轴向力F不发生变化的情况下，单一改变接触宽度的大小，建立接触宽度b=30、40、50 mm的分析模型.此实验过盈量为0.059 mm，包容件的外径Φ48 mm.图5a、5b、5c表示不同接触宽度下包容件的应力分布情况，通过对结果的观察分析可知，包容件的接触区域应力分布受到边缘效应的影响，在靠近边缘处的中间部位出现了低应力区域.图6a、6b、6c为空心轴的应力分布情况，从图中可以观察到空心轴的应力分布情况呈现一定的规律性，轴的应力值从一端到另一端逐渐增加.图7a、7b分别代表包容件与轴，在不同宽度影响下接触区域的应力曲线图.观察图7a可知，包容件的接触应力在其他影响因素不变的情况下，随着接触宽度的增加，包容件的等效应力值降低，宽度越小边缘效应越明显，边缘处应力集中越明显，最大应力值越大.图7b为空心轴的应力分布，从一端到另一端应力值逐渐增加，宽度越大应力值的变化速度越慢，曲线的越平缓，边缘处的应力集中情况越不明显. 2.3 包容件外径大小对接触应力的影响研究对包容件外径大小对配合面的应力分布影响的研究，是在保证过盈量，接触面的基本直径，轴向力，接触宽度，不发生变化的情况下，单一改变包容件的外径大小，来进行模拟装配过程.实验分析中，建立直径接触面Φ40 mm，轴向力F=0，接触宽度b=30 mm，过盈量为0.056，包容件外径大小分别是Φ48 mm、Φ56 mm、Φ62 mm的分析模型.对三个模型进行模拟装配分析，分析情况如图8、图9所示. 图8a、8b、8c分别为外径大小Φ48 mm、Φ56 mm、Φ62 mm包容件的接触区域应力分布图，从分析结果分析可知，在接近边缘的区域出现低应力区.图9a、9b、9c分别是外径为Φ48 mm、Φ56 mm、Φ62 mm的包容件对应空心轴，接触区域的应力分布情况，从图中看出空心轴的等效应力值，从一端到另一端逐渐增加，在边缘处出现最大值.图10a所示为不同外径下包容件的接触区域应力分布情况，包容件的接触区域受边缘效应的影响，最大等效应力值出现在边缘处，同时在接近边缘的中间部位有低应力区域出现.从图10a的曲线图可以看出，包容件的半径越大，从低应力区到应力集中的边缘处这段区域变化越平缓.包容件半径越大，低应力区所占据的区域范围越大，低应力区域的应力值过渡越平缓.包容件半径越大，低应力区的最小等效应力值和边缘处的最大等效应力值越小.图10b所示为空心轴的应力分布情况，从图中可知，接触的初始区域应力值变化都比较平缓，接触中期区域的应力值发生了急剧增加，包容件半径越大的对应的轴应力急剧增加的速度越慢.相反包容件半径小的等效应力值的急剧增加速度越大.2.4 轴向力对接触区应力影响的研究研究轴向力对配合区域接触应力的影响所建立的分析模型，控制单一变量，保持包容件外径，过盈量，接触宽度不变的情况下，单一改变轴向所加载的力的大小和力的方向，本次实验所建立的模型，过盈量为0.059 mm，包容件外径Φ48 mm，接触宽度b=30 mm，轴向力加载，分四种情况，F=±60 MPa和F=±130 MPa，正代表压力，负号代表拉力.模拟过盈装配的分析结果如下.图11a、11b分别代表轴向受60 MPa和120 MPa的压力作用下的配合中包容件的等效应力分布图.图11e代表空心轴在左右两端未受到力的作用时，包容件的接触区域的应力分布图.图12a、12b代表轴向受到60 MPa和120 MPa拉力力作用的配合中包容件的等效应力分布图.图12e代表空心轴不受轴向力时的应力分布图.从图11a、11b分析可知，轴向力没有影响边缘效应的出现，靠近边缘的中间区域仍然出现了低应力区，边缘处出现应力集中.图12a、12b可知空心轴的应力分布仍然保持从一端到另一端逐渐增加的趋势，边缘处应力值最大.图13a、13b分别代表不同轴向力作用下，包容件的接触区域应力分布和空心轴的接触区域应力分布情况.由图13a所示可知，包容件在接触的开始区域与中部区域应力值变化都非常平缓.其中这部分区域的应力值，在不同的轴向力作用下呈现一定的规律性，120 MPa压力＞60 MPa压力＞0作用力＞60 MPa拉力＞120 MPa拉力.图13b分析可知，轴的接触区域应力值分布情况的总体趋势与包容件相反，120MPa压力＜60 MPa压力＜0作用力＜60 MPa拉力＜120 MPa拉力.分析轴向力对应力值分布影响规律可知，轴向力对包容件应力值的影响，实质上是轴向力影响了接触区域的过盈量的值.轴在受到压力时，发生了弹性变形，半径变大，相当于增加了轴的过盈量.故压力作用下应力值高于没有作用力时的应力值，高于拉力作用下的应力值.1）圆柱面形状的过盈配合应力分布情况，受过盈量、结合面宽度、包容件直径、轴向力的影响.同时包容件的应力分布情况近似成对号函数的分布，轴的应力分布呈现近似于线性分布，包容件与轴都呈现出边缘效应，包容件的现象更明显，由边缘效应所产生的应力集中现象，应力值远远高于其他部位的应力值，这对过盈配合的可靠性产生极坏的影响.2）过盈量的大小对应力值的影响有着直接对应的关系，随着过盈量的增加，应力值增加.接触区域宽度对过盈配合区域应力值分布有一定的影响，接触区域越宽，包容件的应力集中现象越弱，边缘处的最大等效应力值越小，过盈连接越可靠.包容件外径大小对应力分布有影响，但实验结果的曲线具有一定的复杂性，总体表现是半径越小，包容件的边缘处应力集中越明显，应力值变化越快.3）轴向力对应力分布影响显著，包容件的区域应力值分布呈现出的规律为轴两端受到压力时，大于两端不受力，大于两端受到拉力作用.在承受范围内，两端力越大，这种现象越明显.轴的表面接触区域的应力值分布，在受到拉力作用时，区域的应力值高于受到压力作用，两端的力越大这种现象越明显.影响包容件应力值分布的原因，最终可以归结于过盈量的变化，轴在受到压力和拉力时，发生弹性变形改变原有过盈量，引起包容件应力分布的变化.【相关文献】［1］杨广雪，谢基龙，李强，等.过盈配合微动损伤的关键参数［J］.机械工程学报，2010，46（16）：53-59.［2］黄庆学，王建梅，静大海，等.油膜轴承锥套过盈装配过程中的压力分布及损伤［J］.机械工程学报，2006，42（10）：102-108.［3］陈连.过盈连接可靠性设计研究［J］.中国机械工程，2005，16（1）：28-30.［4］张洪才.Ansys14.0理论解析与工程应用实例［M］.北京：机械工业出版社，2012.［5］刘江.Ansys14.5 workbench机械仿真实例详解［M］.北京：机械工业出版社，2014. ［6］张杰，韩传军.轴向载荷下空心轴过盈联接的力学特性［J］.四川大学学报，2013，45（S1）：53-57.［7］李伟建，潘存云.圆柱面过盈连接的应力分析［J］.机械制造与技术，2008，27（3）：313-317.［8］张敬佩，李初晔.过盈配合产生的接触压力和拔出力计算[J].机械设计与制造，2010（10）：195-197.［9］王少江，李学明，张安鹏，等.大型焊接齿轮与轴过盈配合有限元分析［J］.煤矿机械，2015，36（5）：134-136.［10］滕瑞静，张余斌，周晓军，等.圆柱面过盈连接的力学特性及设计方法［J］.机械工程学报，2012，48（13）：160-166.［11］李伟建，潘存云.圆柱面过盈连接的应力分析［J］.机械科学与技术，2008，27（3）：313-317.［12］张杰，韩传军.轴向载荷下空心轴的过盈连接力学特性［J］.四川大学学报：工程科学版，2013（S1）：53-57.。
基于ANSYS组合模具过盈配合有限元分析
利用时间历程后处理器分析１００端面处某一节点的 约束反力随时间变化曲线。如图８所示。
３结语 通过采用ＡＮＳＹＳ有限元软件的分析，该组合模具的
过盈配合过程中产生的最大应力为１２４０ＭＰａ，没有超出 材料的许用应力（１６００ＭＰａ），表明设计没有使产品产生
基于ＡＮＳＹＳ组合模具过盈配合有限元分析
郝兆朋，张永军，范依航 （哈尔滨理工大学机械动力工程学院，哈尔滨１５００８０）
摘要：在细长杆冷成形自动锻压机新产品的研制过程中，为判断其组合模具结构设计是否符合要求，文中用接触有限 元法建立了内衬与外套过盈配合的有限元力学模型。针对内衬和外套的过盈配合状态。采用大型通用有限元ＡＮＳＹＳ软
比较方案一与方案二，对于所给零件，由于两小孔比
较接近边缘，复合模冲裁零件时受到壁厚的限制，模具结 构与强度方面相对较难实现和保证，故采用级进模加工。 ２模具设计计算 ２．，排样、计算条料宽度及确定步距
采用单排方案，如图２所示。
查表确定搭边值，根据零件形状两零件间按矩形取 搭边值ｂ＝２．Ｏｍｍ，侧边取搭边值ａ＝２．Ｏｍｍ。
万方数据
机械工程师２００８年第５期《２７
专题报道 譬墨蚕蕊圈瞄冒誓墨赣县加工菝７Ｉｔ与设备■●—●—■—●————■———●—■—■—■—————■—■—■■—●——■——————一
垫圈连续冲栽模具的设计
李文 （韶关机电技校，广东韶关５１２０００）
“
。
。
”。
ｊ秀一
摘要：分析了垫圈的工艺特点，介绍了垫圈的排样设计方案，以及模具结构的设计计算、模具主要零件计算、冲床参数
Ansys_过盈配合接触设置
Finish /clear,start
/TITLE,Contact analysis with initial interference /PREP7 !带孔圆盘的基本尺寸； Rpin=35 Rpout=100 Hp=25
!轴的基本尺寸； Rain=25 Raout=35
La=150
! 过盈量 f; f=0.01
FINISH /POST1 SET,1,LAST,1, PLNSOL,s,eqv
!查看单元的 Von Mises stress;
esel,s,ename,,174
PLNSOL, CONT,PRES, 0,1.0 !在 contact 单元上查看 contact pressure; SAVE 计算结果：
在 ANSYS 中正确地模拟过盈配合
过盈配合在机械产品的装配中使用相当普遍,譬如轴与轴承,轴与轴瓦,汽车 的制动盘等,都是通过一定的过盈量来使两个装配部件紧密连接起来。
下面讨论如何在 ANSYS 中正确地模拟过盈配合。过盈配合在有限元分析中 是一种典型的非线性接触行为。在有限元分析中设定了接触，从本质上来讲就是 对相互接触的两个部件施加了某种约束，不同的接触算法对于接触约束的处理方 法有所不同。接触约束的理论算法的选择，在 ANSYS 中是通过设置 contact 单 元的 KEOPT(2)选项来实现的。在 ANSYS 中目前主要有 5 种接触约束算法：
ANSYS中正确地模拟过盈配合
在ANSYS中正确地模拟过盈配合过盈配合在机械产品的装配中使用的相当普遍。
比如轴与轴承、轴与轴瓦、汽车的制动盘等，都是通过一定的过盈量来使两个装配部件紧密地连接起来。
下面讨论如何在 ANSYS 中正确地模拟过盈配合。
过盈配合在有限元分析中是一种典型的非线性接触行为。
在有限元分析中设定了接触，从本质上来讲就是对相互接触的两个部件施加了某种约束，不同的接触算法对于接触约束的处理方法有所不同。
接触约束的理论算法的选择，在 ANSYS 中是通过设置 contact 单元的 KEOPT(2) 选项来实现的。
在 ANSYS 中目前主要有 5 种接触约束算法：KEYOPT(2)=0 Augmented Lagrangian -加强的拉格朗日算法，这是 ANSYS 的缺省选择；KEYOPT(2)=1 Penalty function -罚函算法；KEYOPT(2)=2 Multipoint constraint (MPC) -多点约束算法；KEYOPT(2)=3 Lagrange multiplier on contact normal and penalty on tangent -接触法向采用拉格朗日乘子，接触切向采用罚函数的综合算法。
KEYOPT(2)=4 Pure Lagrange multiplier on contact normal and tangent -法向和切向均采用拉格朗日乘子算法。
各种不同的约束算法各有其优缺点，各有各自最适用的场合，具体情况需要具体对待。
大部分情况下，默认选择KEYOPT(2)=0 就够用了。
过盈配合所致的接触分析的难点在于如何确定初始接触状态。
初始接触状态设置得不对，会导致错误的计算结果或者不准确的计算结果，下面举两个例子来说明。
例1．两个圆柱体在几何上是刚好接触，划分网格后有限元模型有间隙。
如图 1 所示。
这两个圆柱体，在几何上是刚好相切的，即处于几何上刚好接触的初始状态。
过盈配合的有限元分析
过盈配合的有限元分析工程力学系张晨朝20803001过盈配合的有限元分析摘要: 在工程应用中,利用接触有限元法建立了内轴与外套过盈配合的有限元力学模型来判断结构设计是否符合要求。
针对内轴和外套的过盈配合状态，采用大型通用有限元ANSYS 软件对组合模具进行了有限元分析, 得出了内轴与外套在过盈配合状态下的应力分布规律及接触面压力分布状况, 找到了应力集中位置和大小。
结论表明结构配合尺寸设计没有使结构产生变形, 该结构完全符合产品的设计要求。
关键词: 过盈配合; ANSYSAbstract: In the project application, in order to judge whether the structural design meets the requirement, the finite element and mechanical model of the interference joint between inside lining and outside wrap is established by used contact -finite- element methods. Aimed at condition of the interference joint between inside lining and outside wrap, we carry on the finite element analysis based on ANSYS and attain the stress distribution in interference joint; the pressure distribution in contact face and the location and the size of stress concentration. It is concluded that the structure interference joint size of combined die do not make the mold have distortion and the combined die completely meets the product design requirement.Key words: interference joint; ANSYS1 引言过盈配合[1]是机械工业中一种常见的零部件组装方式,齿轮、轴承以及火车车轮等与其装配轴之间的配合大多采用过盈配合。
在ANSYS中怎样精确地设置过盈量
在ANSYS中怎样精确地设置过盈量下面讨论如何在ANSYS中正确地模拟过盈配合。
过盈配合在有限元分析中是一种典型的非线性接触行为。
在有限元分析中设定了接触，从本质上来讲就是对相互接触的两个部件施加了某种约束，不同的接触算法对于接触约束的处理方法有所不同。
过盈配合所致的接触分析的难点在于如何确定初始接触状态。
初始接触状态设置得不对，会导致错误的计算结果或者不准确的计算结果，下面举2个例子来说明。
例1．两个圆柱体在几何上是刚好接触，划分网格后有限元模型有间隙。
如图1所示。
图1. 两个在几何上刚好接触的圆柱体这两个圆柱体，在几何上是刚好相切的，即处于几何上刚好接触的初始状态。
划分网格后，由于在圆周上用小段直线代替了弧线，两个圆柱体之间产生了一定的间隙，两个圆柱体的有限元的初始状态不再是有接触的，此时，如果接触参数设置不当，就会因为初始约束不足，圆柱体出现刚体位移，得到错误的结果。
例2．有的人把两个接触部件的几何位置设定一定的过盈量，想用这个过盈量来模拟过盈配合，这种做法是错误的，几何上的过盈量不等于划分网格后有限元模型的实际过盈量。
下面的图2中，是一个孔类零件和一个轴类零件的截面图，轴和孔在几何位置上预设了过盈量。
(内圈的红色圆是孔边界，外圈的蓝色圆是轴边界，轴和孔在几何上是相互侵入的)。
图2：一个轴类零件和孔类零件的过盈配合的截面图在几何上，图2的轴和孔有一定的过盈配合量，其大小等于两个圆的半径之差，我们的本意是想用这个几何位置上的过盈量来模拟过盈配合。
不幸的是，两个部件划分网格之后，实际的过盈量应该为单元之间的距离，即图上中靠得比较近的两条线段之间的距离，显然，这个距离不再等于我们预先设置的过盈量了。
更何况，上面这个图还是两个部件的网格对对得比较整齐的情况，如果网格对的不整齐，过盈量就和我们预设的差的更远了。
对于过盈配合来讲，过盈量的数值变化对于过盈产生的应力的影响是很大的。
上面两种方法都不能精确地设置过盈量，下面介绍第3种方法。
基于ansys的汽轮发电机转轴的过盈配合分析
基于ANSYS的汽轮发电机转轴的过盈配合分析黄鹏，倢任，张金华（中国长江动力集团有限公司，武汉435200）摘要：汽轮发电机转轴在超转速的情况下，套筒可能会因为过盈量不够而与转轴分离，从而影响发电机的正常运行，所以对转轴与套筒的最小过盈量的计算分析就变得尤为重要。
文中以一台汽轮发电机转轴与励磁套筒过盈配合的结构为例，建立了其有限元计算模型并计算了最小过盈量。
将有限元仿真的结果与理论计算的结果进行对比，结果表明，该转轴与套筒过盈配合所产生过盈量能满足设计的要求，同时验证了其结构设计的合理性。
关键词：转轴；过盈配合；最小过盈量；有限元仿真中图分类号:TM311；TH124文献标志码:A文章编号：1002-2333（2020）02-0128-03 Analysis on Interference Fit of Turbo Generator Shaft Based on ANSYSHUANG Peng,REN Jie,ZHANG Jinhua（China Changjiang Energy Co.,Group,Wuhan435200,China）Abstract:At the super-rotation speed of the turbine generator shaft,the sleeve may be separated from the shaft due to insufficient interference,which affects the normal operation of the generator.Therefore,the calculation and analysis of the minimum interference of the shaft and the sleeve become especially important.Taking a structure of the turbine generator shaft and the excitation sleeve interference fit as an example,the finite element calculation model is established and the minimum interference is paring the results of the finite element simulation with the theoretical calculation results,the results show that the interference generated by the interference fit between the shaft and the sleeve can meet the design requirements,and the rationality of the structural design is verified.Keywords:shaft;interference fit;minimum interference;finite element simulation0引言转轴是汽轮发电机的关键部件之一，其主要作用是用来传递转矩。
基于ANSYS的圆弧端齿环过盈配合分析
基于ANSYS的圆弧端齿环过盈配合分析为研究过盈配合摩擦系数及过盈量对圆弧端齿环应力和变形的影响，以圆弧端齿环为研究对象，利用有限元数值模拟方法研究了圆弧端齿环的应力分布及变形情况，分析了过盈配合面的摩擦系数及过盈量对最大等效应力及最大总变形量的影响规律。
标签：圆弧端齿环；摩擦系数；过盈量Interference Fitting Analysis of Curvic Couplings Rings Based on ANSYSKuang Kaixin Zhou Zetang（AECC Hunan Aviation Powerplant Research Institute，Hunan Key Laboratory of Turbomachinery on Small and Medium Aero-Engine，Zhuzhou ，China 412002 ）Abstract：Under the interference fitting condition，in order to study the influence of friction coefficient and magnitude of interference on Curvic Couplings Rings，Curvic Couplings Rings was sat up as research object. Using finite numerical simulation method，the influence of stress distribution and deformation on Curvic Couplings Rings was carried out，and the influence law of maximum equivalent stress and maximum deformation on Curvic Couplings Rings were analyzed.Key words：Curvic Couplings Rings；Friction coefficient；Magnitude of interference圆弧端齿是端齿盘的一种特殊结构形式，轮齿在端面沿圆周方向分布并具有弧形特征，具有自动定心、承载荷能力强、定位精度高、装配简便、多次拆装不影响其原有的精度和可靠性等诸多优点。
盘轴过盈配合的ANSYS分析
（４）创建接触对在Ｃｏｎｔａｃｔ Ｍａｎａｇｅｒ（接触管理 器）中分别将圆盘的盘心面轴的外环面选中，在Ｃｏｅｆ－ ｆｉｃｉｅｎｔ ｏｆ Ｆｒｉｃｔｉｏｎ（摩擦系数）文本框中输入“０．２”。 将Ｎｏｒｍａｌ Ｐｅｎａｌｔｙ Ｓｔｉｆｆｎｅｓｓ（正则处罚刚度）定义为 “０．１”将Ｓｔｉｆｆｎｅｓｓ ｍａｔｒｉｘ（刚度矩阵）下拉框设定为 “Ｕｎｓｙｍｍｅｔｒｉｃ”选项。创建的接触对，如图３所示。
然后施加位移载荷（将轴沿轴向平移４０ｒａｍ拔出 盘孔）。选取菜单路径Ｕｔｉｌｉｔｙ ｍｅｎｕ Ｉ Ｓｅｌｅｃｔ Ｉ Ｅｎｔｉ－ ｔｉｅｓ。单击对话框中最上面下拉框中的“Ｎｏｄｅ”，指定 选择对象为节点。接着在下面的下拉框中单击“Ｂｙ ｌｏｃａｔｉｏｎ”，指定选择方式为根据坐标值来选取。单击
万方数据
图７接触单元上的压应力分布云图
４计算结果的提取分析ｐＪ （Ｉ）利用通用后处理器观察结果 首先将建立
的Ｉ／４模型扩展为整个盘轴结构模型，然后利用通用 后处理器来提取求解结果。依次提取过盈配合应力 分布云图如图６所示，接触单元上的压力分布云图如 图７所示。
荷子步数（Ｎｕｍｂｅｒ ｏｆ ｓｕｂｓｔｅｐｓ）为１，关闭自动时间
步，其求解过程时间曲线，如图４所示。
本文读者也读过(5条) 1. 安吉阁.李建宾.廖熠.AN Ji-ge.LI Jian-bin.LIAO Yi 基于ANSYS转子曲轴过盈配合有限元分析[期刊论文]-机 械工程师2010(11) 2. 杨加军.胡安.李维波.章向明.YANG Jia-jun.HU An.LI Wei-bo.ZHANG Xiang-ming 过盈配合数值分析与优化设 计[期刊论文]-海军工程大学学报2008,20(1) 3. 曹燎原 ANSYS在过盈联接中的应用[期刊论文]-现代制造工程2005(5) 4. 韦念龙.邓群.周晔.尹玉波 有限元分析软件在过盈配合联接安全核算中的应用[期刊论文]-安全与环境工程 2002,9(3) 5. 韩志峰.李云峰.HAN Zhi-feng.LI Yun-feng 基于ANSYS的电主轴轴承预紧与过盈配合分析[期刊论文]-机械 2010,37(12)
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第47卷第5期2011年3月机械工程学报JOURNAL OF MECHANICAL ENGINEERINGVol.47 No.5Mar. 2011DOI：10.3901/JME.2011.05.121基于ANSYS的轮对过盈配合微动分析*曾飞1, 2陈光雄1周仲荣1(1. 西南交通大学摩擦学研究所成都 610031；2. 中国工程物理研究院总体工程研究所绵阳 621000)摘要：铁道机车车辆在运行过程中，轮对过盈配合面边缘由于微动产生微动损伤。
微动幅值是影响轮对过盈配合面微动运动特性的重要因素之一，由于难以用仪器进行测量，利用通用有限元软件ANSYS对RD2轮对过盈配合面在210 kN轴重载荷作用下的微动情况进行了模拟分析。
以轮对受轴重载荷静止于轨道上的弯曲变形情况来表征其运行过程中某一时刻的弯曲变形情况，通过计算，获得轮对过盈配合面内、外侧区域的轮座与轮毂某接触节点副的相对滑动规律和应力分布。
结果表明，在轮对转动过程中，轮座与轮毂配合面内、外侧接触区内的节点副相对运动模式为变方向、变应力的复合微动。
在210 kN 轴重载荷作用下，轮对过盈配合面间所计算节点副的最大轴向相对位移为32 μm，最大切向相对位移为2.6 μm。
关键词：铁路轮对过盈配合复合微动有限元中图分类号：TH117Fretting Analysis of Interference Fitting of Wheel-set Based on ANSYSZENG Fei1, 2 CHEN Guangxiong1 ZHOU Zhongrong 1(1.Tribology Research Institute, Southwest Jiaotong University, Chengdu 610031;2. Institute of Systems Engineering, China Academy of Engineering Physics, Mianyang 621000)Abstract：During wheel-set of railway vehicle rolling, the edges of the interference fitting surfaces between wheel seat and wheel hub are worn due to fretting. Fretting amplitudes between wheel and axle are difficult to be measured, which is one of the important factors influencing fretting characteristics. The fretting of the interference fitting of RD2 wheel-set under a load of 210 kN is simulated and analyzed by using ANSYS. The bending deformation of the wheel-set at a moment during its rolling is characterized by that of the wheel-set under a load of 210 kN which is statically on railway. The relative sliding mode and stress distribution of a certain contact node pair on the interference fitting surfaces between the wheel seat and wheel hub are obtained. The result shows that during the wheel-set rolling the contact node pair is subjected to a composite fretting consisting of reciprocating sliding and varying stress. Under an applied load of 210 kN, the maximum relative displacement in the axial direction of the contact node pair is about 32 μm and the maximum relative displacement in the circumferential direction is about 2.6 μm.Key words：Railway wheel-set Interference fitting Composite fretting Finite element0 前言微动是两个接触面之间极小幅度(通常为微米量级)的周期性运动[1]。
在铁道机车车辆的运行过程中，车轴受到弯曲交变载荷的作用，由于轮座与轮毂结构之间的刚性不同，轮对过盈配合面边缘存在不同程度的相互微滑，即为微动。
∗国家自然科学基金(50821063)和国家重点基础研究发展计划(973计划，2007CB714700)资助项目。
20100210收到初稿，20101222收到修改稿黄梦妮等[2]为了观察轮座边缘的微动损伤情况，并避免因为退轮给轮座带来二次损伤，采用了原位剖切的方法，将一个使用了20年但未经过退轮维修的轮对轴向剖开，观察到轮座与轮毂配合面边缘约20 mm的范围为微动区，其表面分布有许多磨坑，且堆积了大量的磨屑及油垢，同时还观察到早期疲劳裂纹。
由此可见，微动使轮对配合表面产生磨损、腐蚀，并萌生疲劳裂纹[3]，随着机车车辆的运行，疲劳裂纹发生扩展，最终导致车轴断裂[4]。
轮对的微动现象十分复杂，它受轮对过盈配机械工程学报第47卷第5期122合表面间相对滑移幅值、轴应力、接触载荷、几何形状、切向力、摩擦因数等诸多因素的影响[5]，而在这些因素中，滑移幅值是影响微动运动特性的重要因素之一[6]。
张进德等[7]测量出RD2型车轴的轮毂边缘与轮座的相对滑移幅度为27~63 μm，但是轮对过盈配合面微动区中部的相对滑移幅值却难以用仪器来测量。
因此，本文利用有限元软件ANSYS对轮座与轮毂过盈配合面上内、外侧微动区内某节点副的相对位移和应力进行计算[8]，获得了该节点副的相对滑移幅值和应力变化规律。
1 计算过程1.1实体建模选取RD2型车轴为研究对象，根据铁道部相关标准尺寸进行建模。
车轴长2 146 mm，轮座直径D1=194 mm。
轮座与轮毂配合过盈量为轮座直径的0.08%～0.15%(0.155～0.291 mm)，计算选取过盈量为0.2 mm，因此轮毂内径D2=193.6 mm。
由于轮对的对称性，取1/4轮对建立模型。
1.2 划分网格选取单元类型为实体单元SOLID185，为了简化计算，将车轴与轮毂设定为相同的结构钢材料，其弹性模量E x=206 GPa，泊松比μ=0.3。
对1/4的轮毂与轮座进行网格划分，同时将轮座部位网格进行细化，共划分网格数为11 732个，见图1。
图1 1/4轮对网格图1.3 加载与求解对轮对运行过程中轮座与轮毂配合面间的微动情况进行模拟，为了简化计算，以轮对在实际工况下静止于轨道上时的应力与变形情况来表征轮对运行过程中某一时刻的应力与变形情况。
计算进行两个载荷步，第一步计算轮座与轮毂装配好后，无其他约束与载荷的情况下，过盈配合面之间的应力；第二步对车轮与轨道的接触点施加约束，同时对车轴轴颈施加轴重载荷，计算轮座与轮毂之间的应力及变形情况。
1.3.1载荷步一利用ANSYS的接触向导设置轮座与轮毂的配合面，设定二者间的摩擦因数为0.2，进行大变形非线性计算，计算出轮对过盈配合等效应力如图2，轮座与轮毂配合表面轴向等效应力分布如图3。
可以看出，轮毂内表面的等效应力大于轮座表面，最大值约为329 MPa；轮毂内表面两边缘应力大于中部应力，存在应力集中的现象；轮座表面两边缘应力小于中部应力。
图2 轮对过盈配合等效应力云图图3 轮座与轮毂配合表面轴向等效应力分布1.3.2载荷步二(1) 受力分析。
以轮对轴颈受轴重载荷，静止在轨道上时的状态来施加约束与载荷。
对车轮与轨道接触处节点施加x、y、z方向的约束，在轴颈与轴承配合的上表面中心节点施加指向地面的力52 500 N，以模拟210 kN轴重，计算后等效应力云图见图4。
由图4知，最大应力发生在轮毂内孔表面下侧外边缘处，应力值约为375 MPa，轮座下表面的应力小于上表面的应力。
2011年3月曾飞等：基于ANSYS的轮对过盈配合微动分析123图4 轮对等效应力分布图(2)微动分析。
从轮对的变形中可以分析出轮座与轮毂配合表面的微动情况。
通过计算后的变形显示，轮对在轴重载荷F的作用下，车轴中部向上弯曲，车轮绕与轨道的接触点各向两侧转过了一个微小的角度，见图5(为了便于叙述，定义轮座轴颈一侧为外侧，另一侧为内侧)。
由图5中可以看出，车轴以轮座与轮毂的配合部位为支点产生了弯曲，轮座表面因此而变形，轮毂内孔表面随轮座的变形而变形，但变形程度较轮座小。
在轮对的反复循环滚动中，轮座与轮毂之间的配合单元产生了持续的变形差与错位，这即是产生微动的主要原因。
图5 实际工况下轮对弯曲变形示意图以车轴轴线上某点为坐标原点，如图6a建立坐标系，曲线ab表示轴线在Oxz平面内的弯曲变形。
如图6b所示等角度均匀选取轮座表面上距内侧边缘约10 mm的半圆周上的17个节点，分别命为节点Z1～Z7，与该17个节点相配合的轮毂内表面的节点分别为节点G1～G17，其中，节点Z1在轮座顶部，节点Z17在轮座底部。
图7列出了节点Z1~ Z17分别相对于节点G1～G17沿轴向及圆周切向的相对位移，由图7可见，轮座节点相对于与其配合的轮毂节点的轴向位移为负，表示轮座表面节点向–z方向做相对运动；轴向相对位移在顶部最大，约22 μm，到底部逐渐减小；在1～17个节点副中，仅有第1个和第17个节点副的相对切向位移为0，它们分别位于轮座的顶部及底部，其余节点副均存在不同幅度的切向相对位移，且均偏向Oyz平面。
由此说明，在轮座与轮毂的微图6 轮座节点选择示意图图7 轮对内侧节点副相对位移滑区内，除了顶部与底部的节点沿轴线方向相互微滑，其余节点的相互滑移方向均与轴线有很小的夹角；轮对在单纯的过盈配合状态下，轮座表面没有产生弯曲变形，且各节点受对称的正压力作用，没有发生切向位移。