第33卷第3期 天 津 理 工 大 学 学 报 JoURNAL oF TIANJIN UNIVERSITY oF TECHN0LoGY Vo1.33 No.3 2017年6月 Jun.2017 文章编号:1673—095X(2017)03—0001—05 基于ANSYS Workbench对圆柱面过盈配合 接触应力的研究 李亚洲 ,汤易升 ,陈 炜 一,张西正1,2 (1.天津理工大学机械工程学院天津市先进机电系统与智能控制重点实验室,天津300384; 2.军事医学科学院卫生装备研究所,天滓300161) 摘要:利用有限元分析软件ANSYSWorkbench,模拟空心轴与包容件的过盈配合过程,通过改变单一因素的实验方 法,分析圆柱面过盈配合中接触应力分布情况和各个因素的关系.本文分析了过盈量、包容件外径、接触宽度、轴向拉力 和压力的影响因素.对实验结果进行分析得出如下结论:轴和包容件的接触区域的最高等效应力区,受到边缘效应的影 响,最高应力集中区出现在端部边缘处.包容件在靠近边缘的中间部位出现低应力区域.空心轴表面的应力值呈现出从 一端到另端逐渐增加的现象. 文献标识码:A doi:10.39690.issn.1673—095X.2017.003.001 关键词:过盈配合;有限元;影响因素;等效应力 中图分类号:TH140 Research 0n contact stress 0f cylinder interference fit based - On the ANSYS Workbench LI Ya—zhou ,TANG Yi—sheng ,CHEN Wei 一,ZHANG Xi—zheng , (1.School of Mechanical Engineering,Tianjin key Laboratory of the Design and Intelligent Control of the Advanced Mechatronical System,Tianjin University of Technology,Tianjin 300384,China;2.Institute of Medical Equipment,Military Academy of Medical Sciences,Tianjin 30161,China) Abstract: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. Key words:interference it;ffinite element;influence factors;equivalent stress 过盈连接普遍采用的方法是热胀冷缩或强力压 入将有过盈量的两个接触面装配到一起l1 .在接触 与断裂的缺点,也解决了部分金属材料和机械结构 进行螺纹加工困难的问题,常用于薄壁件连接处(不 面上会产生很大的径向力.在过盈连接承受轴向力 和扭矩时,接触区域将产生相应的摩擦力与力矩来 抵御外部施加的力和力矩_7_ .过盈连接是一种半永 易加工螺纹),具有良好的对中性,应力分布相对比 较均匀,承载能力强,结构紧凑美观等优 [91. 现阶段对过盈配合的研究,已经在很多方面进 行了.有从过盈连接的接触面变形类型来研究的,将 过盈面的的变形分为弹性变形和塑性变形两个阶段 久性连接.它有很多优点,结构简单,生产成本低,消 除了焊接产生局部应力,零部件容易产生局部疲劳 2016-07-12. 收稿日期: 基金项目: 国家自然科学基金(11302147,51275353). 作者简介: 李亚洲(1990一通讯作者: 张西正(1956一),男,硕士研究生,E—mail:13132208359@163.coin. ),男,博士,研究员,博士生导师,E—mail:z56787@sohu.com. . 一, . 进 小同的研究 I弹 址指 过盈连接件,拆卸后, 的接触 能够恢复剑 行的直径不发 永久性 彤. 性变彤,拨珩虫 Ifif的材料受到很大的力,趟 …』 强度极限, 致接触 发,卜J,永久性变形. 种研究分类,将过 连接的研究,根据接触向figK.状 小¨,分为 锥 洲佳 …对旧柱叫的研究卞『l对比 譬 l目 誊 tl 』I L 较多. 研究领域 巾 轴和轮投酉己合 ,研究 力 法 遍采川j PI-: J 、 j材料力学的 式进行 图l模拟过盈配合过程图 Fig.1 Figure simulating interference fit 推 汁 .应川f/2统的 弹性力学,F 以斛决轴 I 位容什的应力特 … 题.似随荷 程研究的深入,人 f『j ^ 要更lJJfI血观fll准确研究接触 域的非线性 学 川题..些学者捉…J 埘过垃连接应力分伽影响fkl 豢的研究.f{1于 的弹性 J学力 法已经 通川, 尤法满址埘过㈨ 合J 域的 观的定性分析.现阶 段普遍采朋的研究 法是将数值计算的方法fjf入 Y1刈‘过 连接影I 尺1 的研究rf1.随着有 几分忻 的方法的发腮, 多 、 利用有限元分析软件,刈‘ 过 合而进仃J、 分析.能蜓JJ【l准确、直观、细致 的分析ff{过盈 的变形过 I I.本文利川仃 几 分析软件ANs wt I-kt,en(-h,通过控制 一变 的 疗法,埘过 链接的影响 素进行r研究,僻划每 个rfI-变艟 素埘接触I)‘=域心力分布的影响、 力 的影响规律. 1 过盈配合有限元模型的建立 迎 的有 摸' 材料选择为结构铡. 心辅 总长为30 I1Ill1,外 为(1)48 11|l11,内径为 40 111111;411I 企总长为 30 llllll。sic.:-小 合外径为 40 11|111,过盈 世人小取(0.042 IIIIII,0.059 111111,0.076 Illln)l』、J径为 q)32 lllll/轴和 ti企均采Jl J合金钢材,其材料常数为: 性帙毓 =2.0I×10 Pa, 松比v=0.27,材料密 度JD=7 800 kg/hi .fI!lJ 的建● 采用中 格式倒人的 力‘法, 尤 Solid 建立伞心轴与也容什I,I,J 馍, ,干IJfld fI1 I’曰j十}}』l= 人 IJ ANSYS Workbent-h软f,f: lll进仃过德配合的卡Il拟装配.模拟装配过 l 所,J . 利JIJ Solid w…・ks缱 维模型,将缱 枷勺 维模 转化为 f1 I'llJ格式X—I的文件, 入刮 __、 sYS\)'/ 1 t・h It1, 埘网格属性,材料 性进行 ,进行 格的 ̄i]/L7僻剑 1b.刚格划分 fJ ,×lf 包容什施加约 fIJf 移JJl】载.方法,对包容什端邴同 定,川时对 04ill・端施Jjli1'74移,完成的模拟装 过 .得到_『』l】网lt・所爪模拟装酣后的应 分 【殳1. 2对不同过盈配合影响因素的计算结果分析 2.1过盈量变化对等效应力的影响的计算研究 本7欠研究f 助ANs、S Wm,kb<-II(・h软什,迎过控 制单一影响[大I崇,观察结 变化的方法,研究过 li 大小,埘过 合 域的接触应 分布的影响.小}欠 实验, 40 n1111孔 轴的 合采片J TL;h0,过 { 据罔标GBl80l—l979选取H7/n7、H7/s6、I17Iu6, 种过盈配合关系,通过汁箅得到最大的过 分圳 是0.042 lllnl,0.059 1111I1,0.076 Illll1.根据控制f’『L・ 的原理,保持接触宽度,轴 力,包容什外f 人小, 郁小变,单一改变过 till}} 埘 组实验对象进 十I!拟 装 .如2罔所爪, 2a、2b、2 c-分别代表过 :l 0.042 lnlll、0.059 IIIIII、0.072 111111位容件的等效心/J分 布 ,l冬i 2 f1包 f't:lf"J心力分 仃明 的规fit , 靠近边缘的fI1部…脱r低畸 Jx:域.f纠3a、 、 分 别代表伞心轴的J、 力分m情况, 心轴外 I(ii ̄ 诎 I)<域 力分布, 现…从・端到乃一端逐f4 ̄fJff;Jlii'19 变化,最大应力集ffl划插入端. 4a所示 过 为0.042 lllnl、O.059 111111、0.072 111111的包容件,从术做 约束端到被约 端V,j、 J9- ̄qi情况,从曲线 4a【lJ‘ 知,不M过盈 卜‘包 什表面的应 分佰均发 』 突变,…班低啦 J l≮.从 条If}j线的总体趋势水仃. 接触 域的应乃 大小 过 量大小成lI 匕.随过 僦 的增加而增Jil1. 4h代表 过盈 i-,轴的 接触 域的等效 J分 圳I线.从}¨1线中_1f以 …, 轴的接触 域 ,力 ,从 端到 jl_一端不断增JJIl,从 不同过盈量X.J ̄lt刈‘比发现,轴的接触 :域 ,Jf…她 着过能量的增大Ill St . 2.2接触宽度的变化对应力集中情况影响的研究 对接触宽度的 化刈’心力集中情况的 验,采J:fJ同样的力‘法,俅持过 量,位容件外 人IJ、, 轴向力 发 变化的情 下,单 改变接触宽 的大小,建立接m虫宽度『J=30、40、50 111111的分析 J 洲.等:琏于ANSYS Workbench对网卡jm 过 配合接触鹿力的研究 ・3・ r.. ..Lr..-..Lr ..L...r..-..Lr.._..Lr.. ..州 Ⅻ州州 Ⅻ Ⅲ L 彳:彳: 一 .一 一 一 一 一m m 一 一 一 一 一 一 - , 、 、 一 .一 一 一 一 一 一- ‘ 』. ●f d II t- r 图2包容件的接触区域应力分布 Fig.2 Include a contact area stress distribution 图5不同接触宽度包容件的接触区域应力分布 型 ≤茸 蓉颦 Fig.5 DiⅢ Ⅻ ferent contact揪 stress distribution width tolerance of contact area 图3轴的接触区域应力分布 Fig.3 Axial stress distribution of contact area 图6不同接触宽度轴的接触区域应力分布 Fig.6 Diferent contact width of the axial stress distributi0n of contact area 0 5 lO l5 20 25 30 】J】f,兰茸 錾 窖一幕 接触的区域位置 姒 接 差i 一 Il 5【Il l5 【】25“I 5 4【】45 5 c】55 5【】5 10 l 2《I 5“}35 4I】45 51】55 川接触f● 图4包容件与轴的应力分布曲线 Fig.4 Tolerance and axial stress distribution curve iI l J 图7不同接触宽度对接触区应力分布影响曲线图 Fig.7 Diierent contact wifdth influence on the contact zone of stress distribution curve 此实验过 }f}为0.059 nlnl,包容件的外径qb48 inll1. 『殳l 5a、51,、5 c・表爪不 接触宽度下包容件的应 力分 情况.通过对结果的观察分析可知,包容件的 接触 域应 分市受到边缘效应的影响,存靠近边 缘处的巾 部位…现_r低应力 域.冈6a、6l 、6c为 心轴的心力分布情 ,从冈rfI可以观察到 心轴 2.3包容件外径大小对接触应力的影响研究 刈‘包容件外径大小埘配合面的应力分布影响 的研究,是在保证过盈量.接触面的基本直 ,轴向 力,接触宽度。小发乍变化的情况下.甲.一改变包容 件的外径大小,来进行模拟装配过 .实验分析巾, 建立A径接触面qb40 ilun,轴向 =0,接触宽度 的成力分 情况旱现一定的规律性,轴的应力值从 一端剑另一端逐渐增JJI1.I刳7a、7b分别代表包容什 轴,存不 宽度影响下接触Ⅸ域的应力fl}1线 .观 b=30 inin,过值最为0.056,包容件外径大小分别足 48 mm、 56 II1nl、 62 111Ill的分析模型.对。_i个模 进行模拟装配分析,分析情况如图8、 9所/l . 8a、8h、8(,分别为外 大小(1)48 nlln、q)56 II1I11、 察『幺1 7a可知,包容件的接触应力在其他影响因素不 变的情况下,随着接触宽度的增加I,包容件的等效应 力伉降低,宽度越小边缘效应越fJJ】 ,边缘处应力集 【}I越 .最大应力值越大.『冬I 71 为空心轴的应力 分布.从一端到乃一端应力值逐渐增加,宽度越大应 力值的变化速度越慢,fffJ线的越半缓,边缘处的应力 集r 寿况越不f月1 . qb62 Ill1I1包容件的接触 域应力分布 ,从分析结果 分析IIJ‘知,在接近边缘的 域m现低应力区.『车1 9a、 91)、9( 分别是外径为 48 I ̄1111、 56 Inln、 62 Inin的 包容件对应空心轴,接触 域的应力分布情况,从 ・4・ 火津理 I:大学学报 第33巷 第3期 初始 域应力值变化都比较平缓,接触巾期 域的 应力位发生_r急剧增加,包容件 径越大的对应的 轴应力急剧增加的速度越慢.相反包容件半径d,;19 等效应力值的急剧增加速度越大. 2.4轴向力对接触区应力影响的研究 研究轴向力对配合 域接触心力的影响所建 的分析模 ,控制单一变量,保持包容件外径,过盈 图8不同外径包容件的接触区域应力分布 Fig.8 Different diameter tolerance area of contact 量,接触宽度 变的情况下,单一改变轴向所 载 的力的大小乖¨力的方向,本次实验所建立的模 , stress distributi0n 过盥量为0.059 mm,包容件外径@48 inm,接触宽度 b=30 nllll,轴向力加载,分四种情况,F=±60 MPa f¨ F=±l30 MPa, 代表压 ,负号代表托力.模拟过髓 装配的分析结果如下. 1la、llI)分别代表轴向受60 MPa和l20 MPa 的压力作用下的配合[f1包容件的等效应力分布 . “ h (- 冈J】e代表空心轴在左右两端未受到力的作JI1H,t, 包容件的接触I 域的应力分布图. 12a、l2I 代表 轴 受到60 MPa和l20 MPa拉力 作用的配合中 图9不同外径轴的接触区域应力分布 Fig.9 Stress distribution of different diameter shaft contact area 包容件的等效应 分布 .冈12e代表空心轴小受 轴向 时的应力分布图.从 1Ia、lll 分析呵知.轴 向力没有影响边缘效应的¨;现,靠近边缘的rfl问 域仍然…现 低应力 ,边缘处…现应力集 fI.1 I 12a、l2b可知窄心轴的应力分布仍然保持从一端刮 另一端逐渐增加的趋势,边缘处应力值最大. r 香… 心轴的等效应力值,从一端到 一端逐渐 增加,存边缘处ff 5现最大值. j到10a所爪为小I叫外径下包容件的接触 域应 力分布情况,包容件的接触rx:域受边缘效应的影响, 最大等效 /fJ值fn现 边缘处。同时在接近边缘的 r…1 1_IJ部位仃低应力I 域fIJ现.从l lOa的『ffI线 可 以看 ,包容件的半径越大,从低应力 剑应力集巾 的边缘处这段 域变化越平缓.包容件半径越大,低 应力 所I 据的区域范同越大,低应力 域的应力 伉过渡越平缓.包容件半径越大.低应力 的最小等 效 力他和边缘处的最大等效应 值越小.罔10b 昕示为 心轴的应力分布情 .从图中可知,接触的 280 260图l 1轴向力影响下的包容件接触应力分布图 Fig.1 l Include a contact stress distribution under the influence of axial force 萋 2:40。 蠡200 蔫120 .。。80 图l0不同包容件直径下的包容件与轴的 接触区域应力曲线图 Fig.10 Diferent tolerance diameter tolerance parts in contact with the shaft stress graph area 图12不同轴向力影响下的轴接触应力分布图 Fig.12 Diferent axial force under the influence of shaft contact stress distribution 2017年6月 李亚洲,等:基于ANSYS Workbench对圆柱面过盈配合接触应力的研究 .5. 图13a、13b分别代表不同轴向力作用下,包容 件的接触区域应力分布和空心轴的接触区域应力分 布情况.由图13a所示可知,包容件在接触的开始区 接触区域越宽,包容件的应力集中现象越弱,边缘处 的最大等效应力值越小,过盈连接越可靠.包容件外 径大小对应力分布有影响,但实验结果的曲线具有 域与中部区域应力值变化都非常平缓.其中这部分 定的复杂性,总体表现是半径越小,包容件的边缘 一区域的应力值,在不同的轴向力作用下呈现一定的 规律性,120 MPa压力>60 MPa压力>0作用力> 60 MPa拉力>120 MPa拉力.图13b分析可知,轴的 接触区域应力值分布情况的总体趋势与包容件相 反,120MPa压力<60 MPa压力<0作用力<60 MPa 拉力<120 MPa拉力.分析轴向力对应力值分布影响 规律可知,轴向力对包容件应力值的影响,实质上是 轴向力影响了接触区域的过盈量的值.轴在受到压 力时,发生了弹性变形,半径变大,相当于增加了轴 的过盈量.故压力作用下应力值高于没有作用力时 的应力值,高于拉力作用下的应力值. 图13不同轴向力作用下的包容件与轴应力分布曲线 Fig.13 Diferent axial force under the action of tolerance and axial stress distribution curve 3结论 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] 张洪才.Ansysl4.0理论解析与工程应用实例[M].北京: 机械工业出版社,2012. [5] 刘 江.Ansys14.5 workbench机械仿真实例详解[M]. 北京:机械工业出版社,2014. 『6] 张杰,韩传军.轴向载荷下空心轴过盈联接的力学特 性[J].四川大学学报,2013,45(S1):53—57. [7] 李伟建,潘存云.圆柱面过盈连接的应力分析[J].机械 制造与技术,2008,27(3):313—317. [8] 张敬佩,李初晔.过盈配合产生的接触压力和拔出力计 算[JJ.机械设计与制造,2010(10):195—197. 『9]王少江,李学明,张安鹏,等.大型焊接齿轮与轴过盈配 合有限元分析[J].煤矿机械,2015,36(5):134—136. f10]滕瑞静,张余斌,周晓军,等.圆柱面过盈连接的力学特 性及设计方法[J].机械工程学报,2012,48(13):160- 166. [11]李伟建,潘存云.圆柱面过盈连接的应力分析[J].机械 科学与技术,2008,27(3):313—317. 『12]张杰,韩传军.轴向载荷下空心轴的过盈连接力学特 性『J].四川大学学报:工程科学版,2013(S1):53—57.