首页/文章/ 详情

接触磨损怎么做?细说ANSYS帮助文档中的接触磨损仿真案例!

7月前浏览633


前两期文章中,本公 众号论述了一些磨损的基础知识,并讲解了ANSYS中磨损的理论知识。链接如下:

材料磨损与材料性能关系简析!减小摩擦和磨损有哪些措施?

磨损分析案例前瞻,一文速通ANSYS 接触磨损分析理论!


本文以ANSYS中的官方案例,描述磨损分析的过程。理论结合实践,更深入理解磨损分析。



--01--

建模与有关设置


①模型及模型简化

帮助文档中,官方案例是一个半径为30mm的半球形铜环在由钢制成的矩形环上旋转,该矩形环的内半径为50mm,外半径为150mm。半球形环从旋转轴的中心(位于100mm处)开始与平环相触。该半球环承受着4000N/mm2的压力载荷,并以100,000转/秒的频率进行旋转。半球环在矩形环上的滑动会导致环体磨损。


在spaceclaim中建模如下图所示。这里需要注意:ansysworkbench中二维轴对称模型软件要求对称轴必须为整体坐标系下的y轴,因此在建模时,必须以y轴为对称轴建模。

以上模型可以简化为两个环的二维轴对称模型,以减少计算量。环体由2D轴对称PLANE182单元(KEYOP(3)=1)进行网格化处理。通过用接触单元(CONTA172)和目标单元(TARGE169)覆盖两个环体表面,来模拟两者之间的无摩擦接触,如下图。


在接触磨损分析理论中提到,磨损模型的磨损量是通过公式计算得到的。磨损计算时速度可以隐式包含在磨损系数中(不显式建模),即形成对磨损系数的缩放,简化分析模型,减少计算量。


本例中,模型被简化为2D轴对称模型,半球形铜环旋转时,速度方向垂直于模型二维平面,因此速度无法被显示建模,因此通过对磨损系数的缩放实现,详见后续设置。也正因如此,可采用稳态静力学模块(static structural)进行分析。

 



② 材料设置及轴对称设置

在spaceclaim中建立模型后,在engineering data中输入铜(cu)、钢材(steel)材料数据,如下图所示:

材料参数输入完成后,双击model进入分析界面前,单击geometry,在右侧属性栏中将分析类型改为2D,如果此处不更改,软件将会将这个平面识别为壳体。如下图:


进入model界面后,首先将各材料赋予模型。如下图所示:

接着,单击geometry,在详细设置面板将2D行为改为轴对称,如下图所示:

③ 接触设置
本例中,接触采用无摩擦接触(frictionless),分析先采用非对称接触,即只能看到接触面的磨损结果。接触算法采用增广拉个朗日算法,接触探测采用节点探测(任意节点探测算法均可),如下图各所示。接触这么设置的理由详见《磨损分析案例前瞻,一文速通ANSYS 接触磨损分析理论!》,此处不赘述。

本例接触磨损分析通过插入命令流实现。因此,在接触选项下插入commands命令,进入命令编辑界面,写入如下命令流:
关于命令流,需要解释的是:
KEYOPT,CID,5,1
即单元关键字KEYOPT(5)是自动调整接触的初始闭合或接触偏移实常数的选项,如下图所示。KEYOPT(5)=1,使用自动接触面闭合间隙,即设置为刚好接触,adjust to touch(也可通过GUI设置)。

 
 

KEYOPT,CID,10,2
是刚度更新(update stiffness)选项。磨损过程中接触节点的重新定位可能会导致接触状态的变化,因此必须合理确定刚度更新选项,通常接触刚度需在每次迭代中更新(KEYOPT(10)=2),也可通过GUI设置。
 
 
 
关于命令流的其他内容,公 众号 在《磨损分析案例前瞻,一文速通ANSYS 接触磨损分析理论!》一文已有论述。这里需要明白磨损系数kcopper是考虑了滑动速度的,滑动速度值作为了磨损系数的缩放因子。即:
kcopper=10e-13*(转速*2πR)

④ 网格划分
网格划分设置如下各图。

网格划分后,可查阅网格划分质量:

网格划分后如果要查看三维模型,可以插入对称设置,如下图所示,网格划分最终效果:

⑤ 边界条件设置

本例中,钢制成的矩形环底部设置为固定约束。如下图:


半球形环顶部需要承受4000N/mm2的压力载荷。为了避免在加载时出现垂直于2D平面法线的转动,即绕着圆中心的转动和X方向的滑动,可以通过远程点和远程位移进行限制,只保持沿着Y方向的平动,如下各图设置:

半球形环顶部需要承受4000N/mm2的压力载荷。由于模型的简化,顶部荷载需转换。

本例采用远程力施加荷载,远程力施加在远程点上,相当于4000N/mm2的压力施加到线体上。可计算出的作用在远程点上的力为F=15079632N.周。
F=4000*3.14159*(130*130-70*70)
  =150796320N≈1.508e8N
因此荷载施加如下:
⑥ 分析设置
本例中,荷载施加为两个载荷步。在第一个加载步等效压力从0上升到4000N/mm2(斜坡加载),在第二个加载步骤中保持恒定在4000N/mm2。磨损在第二个加载步骤中被激活。

因此,分析设置需设置2步。设置如下:
需要注意,磨损分析时,磨损过程中接触节点的重新定位可能会导致接触状态的变化。如果磨损增量太大,所有接触单元可能会从闭合状态变为打开状态,导致刚体运动,无法收敛。为了防止这种情况,通常采用自动时间步,并使用非常小的时间增量,这样磨损增量也很小,接触状态的变化也最小化,有助于收敛。


此外,磨损分析中,存在尺寸的变化,即包括几何非线性,通常需要打开大变形开关。


磨损涉及重新定位接触面节点以模拟材料去除过程。计算时,接触单元下方的实体单元的质量会迅速恶化,通常需插入命令流,激活非线性网格自适应,以帮助收敛。本例设置如下:

关于非线性网格自适应命令流设置可参阅:《ANSYS中单元类型、材料表字段变量、非线性自适应网格定义命令详解!


--02--

结果情况与分析


① 设置为非对称接触时的结果
首先看看磨损量情况。如下图插图磨损结果:

下图为磨损引起的体积减少图。图中,每条斜线的间断处对应着网格自适应的触发,即网格重新划分时间点,可以与收敛过程对应。也可以看到,time为0~1时,即第一个载荷步,未产生磨损,与磨损命令流设置一致。总的磨损量为各斜线峰值相加之和。
磨损体积损失

求解过程中力收敛过程

下图为磨损发生前(time=1)、磨损发生后(time=4)的法向应力(Y轴)情况。可以看到磨损发生后,导致了法向应力减小了。本例计算结果与ANSYS中的帮助文档中的结果基本一致。
接着可以看看接触压力,如下图。接触压力也因为磨损减小了。

为了方便观察可以插入如下命令,查看APDL中的结果,如下图:
下图为最大接触压力的整个过程的变化情况:

② 设置为对称接触时的结果
对称接触时,需要修改接触为对称接触,并且考虑钢圈的磨损,增加钢圈磨损的命令流,如下图,重新计算。


可以看到,磨损前,法向应力和接触压力与非对称时差不多。但磨损后,因为对称接触两个面均存在磨损,法向应力和接触压力都比非对称接触时更小一点。

 

好了,以上就是本期分享的全部内容,如有疑问或意见还请留言交流。本公 众号不定期分享仿真学习知识,案例和感悟,不想错过的朋友可以关注一下本公众 号。


   


来源:薛定谔的Cube
SpaceClaimLS-DYNAWorkbench静力学非线性理论材料ANSYS
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-01-24
最近编辑:7月前
巴郡撸猫人
硕士 签名征集中
获赞 15粉丝 25文章 73课程 0
点赞
收藏
作者推荐

谈谈应力集中、应力集中系数、应力奇异! 一文吃透应力集中的应对策略

应力集中是物体中应力局部增高的现象,这是造成产品发生破坏的重要原因。有统计表明,实际发生破坏的产品零件中,80%是疲劳破坏导致的,而应力集中又是疲劳破坏的主要因素。因此应力集中问题值得重视,笔者收集整了相关应力集中知识点,今天咱们就来聊一聊应力集中的那些事儿!--01--应力集中的概念在材料或结构中,由于几何形状突变(如缺口、孔洞、裂纹、锐边等缺陷)或载荷作用导致局部应力显著高于平均应力的现象,我们称之为应力集中。这种局部应力峰值可能远超材料屈服强度的现象,是结构失效(如疲劳断裂、脆性断裂)的关键诱因。例如下图中,阶梯轴角部处的应力与非角部处的应力差异,在倒角后,角部应力减小,表明角部存在应力集中。同理,开孔部位也存在一样的情况。图1应力集中效应示意再比如,当一个集中力施加在一个构件上的时候,靠近这个力施加点的位置的应力相比其他位置的应力要大。这些都是应力集中的具体表现。通常应力在构件中传递,如同水流,我们可称之为“力流”,在力的传递路径上,空间形态如果发生改变,材料内部应力就会发生改变。因此,经常发生应力集中的部位包括:①曲率发生极大变化的地方,即几何形状突变。如开孔、开槽、折弯根部、台阶轴;图2零部件几何形状改变②材料缺陷部位,如气孔、夹杂等。这些缺陷部位不规则、不光滑,会引起局部的应力集中,造成微裂纹等。图3零部件缺陷--02--应力集中系数通常结构计算中常用应力集中系数对构件的应力状况进行评估。应力集中系数是指危险界面上的应力峰值除以危险界面的名义应力,即K值:这里名义应力指在不考虑几何不连续性(如孔、槽、带、波纹等)的情况下,在试样的有效横截面上计算得到的应力。例如,下图4中,板件原截面积A0,受单向拉伸F作用,发生变形后面积变为A,则名义应力为snom=F/A0。图4板件开孔应力集中同样如上图4,板件危险截面的峰值应力=speak;则应力集中系数=speak/snom。在实际设计中,我们应该如何获取应力集中系数呢?目前,有三种获取应力集中系数的方式:查表法、有限元法、光弹试验法(应力实测)。1)查表法查表法一般用于粗略的校核计算,也是设计中最常用的方式。工程师通过《应力集中手册》,如下图,可方便查取相应的应力集中系数。关于应力集中系数,除了工具书外,也有专业的研究类书籍,可供深入研究。图5应力集中手册如上图5,从应力集中手册的图表可以看出对于阶梯面,倒角直径应尽量大于截面直径的10%以上,才能有效减少应力集中。此外,从上图中还可以看出,应力集中的应力分布方式与构件尺寸和材质无关,它只取决于构件相关尺寸的比值,例如图5中的W/d,r/d。这对于产品设计来说无疑是个好消息,在进行产品设计的时候,无需额外进行试验就可以确定应力集中下最薄弱区域的最大应力。需要强调的是,设计时考虑应力集中系数进行设计不是我们设计的最终目的,而是为了保证设计的零部件有足够的安全储备,即安全系数。如下图6所示。当应力集中系数下降时,根据最大峰值应力计算所得的安全系数将会得到提升,零部件将更加可靠。图6应力集中系数Kt与安全系数Ktg的关系那么通过《应力集中手册》我们如何进行应力集中处应力的计算呢?我们以一个阶梯轴选材为例。首先我们可以根据相应的尺寸,查阅应力集中系数,计算最大的应力,计算过程如下:图7查表法计算峰值应力然后我们可以根据以上峰值应力,推算满足安全系数的材料强度要求,从而选择材料。2)有限元法对于一些零部件,零件异形、对精度要求较高或尺寸在手册中没有的情况,不能通过查表法获取应力集中系数时,我们可以通过有限元分析获取应力集中系数。根据应力集中系数的定义,我们只需获取到应力集中部位的峰值应力与该位置截面的平均应力,取其比值即可。如图8所示,一个平面应力状态的薄板缺口处的应力云图。图8缺口薄板应力集中云图为获取缺口处的应力集中系数,可以首先建立圆弧面路径,如下图9所示。图9建立圆弧面路径然后计算基于等效应力的线性化应力(linearizedstress),如图10所示,计算截面的平均应力(膜应力membrane)。图10计算截面的平均应力接下来,计算截面的总应力。其中,内侧总应力(inside),就是路径起点处的总应力,外侧总应力(outside),就是路径终点的总应力,取二者较大值即可。图11计算截面的总应力最后,我们就可以根据以上获得的平均应力和总应力,计算基于净截面的应力集中系数:K=2.999/1.215=2.468。3)光弹法光弹法是一种利用光学原理和弹性力学知识来测试应力的方法,主要应用于科学研究和工业领域,航空航天、军工等应用较多。是一种非侵入性的测量技术,广泛应用于各种需要精确应力分析的场合。如下图所示。图12光弹测量显示应力以上便是获取应力集中系数或应力的三种方法。从以上内容我们知道,应力集中系数反映出了缺口或圆弧等截面变化处的应力集中情况。在早期计算机还没有这么普遍被应用的年代,由于没有有效便捷的计算方法,工程师们编制了常见结构的应力集中系数手册,但是目前随着有限元等数值方法的广泛应力,我们可以直接求解得到应力集中位置的准确应力,当然应力集中系数还是可以作为结构性能设计的一个评价指标。那么我们采用有限元法计算应力或应力集中系数时,如何判定所获得的应力结果满足精度要求呢?通常我们会认为应力集中处软件默认应力红色区域完整覆盖了2层单元时,网格细化合理,计算精度满足要求。如下图13所示。这主要是因为这种配置可以确保在应力集中区域有足够的网格密度来精确地捕捉应力梯度变化,减小计算误差。有兴趣的小伙伴可以阅读《有限元分析中的“黄金准则”:为什么最大应力处需要至少2层网格?》一文,以做更深入的了解。图13应力集中处最大应力部位网格此外,在应力集中系数计算时,某些时候我们需要考虑线性化应力的路径的问题,如下图14所示。主要原因在于每个截面处都需要满足应力安全评定的要求。图14线性化应力路径--03--如何评估和应对应力集中?通过以上内容我们可以获得应力集中系数,也明白了应力集中系数与安全系数的关系,那么设计中我们如何评估和应对应力集中带来的影响呢?通常我们在考虑应力集中的影响时,需要根据零部件的材料类型、荷载类型、装配要求等来评估应力集中的影响。举个例子,如下图15所示的某孔处的应力云图。图15某板件开孔处的应力集中图中,孔部位出现了应力集中,假设其峰值应力超过了材料的屈服强度。对于塑性材料的静态受载,且为非配合区域,我们通常可以不用在意此孔的应力集中。原因在于:承受静态载荷的塑性材料,例如常用的45钢、6061铝等其应力应变曲线具有明显的屈服延展和应力强化阶段,如下图16所示。图16塑性材料应力应变曲线因此在发生屈服的部位,塑性材料具有较好的延展性,应变在应力作用下会增加,应力会因此下降下来,屈服停止,相当于此处发生了冷作硬化,弹性比例极限提高了。由于材料进入塑性,应变增大,对与有公差配合的部位,会产生公差的损失。因此一般对于静力载荷下,塑性件安全系数(材料的屈服强度或极限强度与工作应力之比)取2.5。这里安全系数计算时的工作应力,应按应力集中部位的峰值应力取值。对于脆性材料,比如铸铁,陶瓷、混凝土等,由于没有屈服后的较大的塑性变形,安全系数会取值高一些,一般安全系数取4.0。对于动载荷或交变载荷,不管是脆性材料还是塑性材料均容易发生脆性断裂,如下所示。因此,在设计中应尽量避免和减小应力集中。图17脆性断裂断面我们知道,应力集中现象是无法消除的,只能通过措施减小应力集中的影响,那么在设计时如何应对即减小应力集中呢?这里咱们介绍一些方法。1)加大倒角或渐变过渡(最通用的做法)在前面的应力集中概念中我们提到,应力在构件中传递,如同水流。“力流”在空间曲率的变化处会出现“拥堵”,即应力集中,因此缓解或减小应力集中的方法就需要在“力流”方向上削弱这种“拥堵”。设计中对于这种“拥堵”处理可以在几何上采取过渡,如锐角、直角、变截面等部位可采用圆弧过渡(如齿轮齿根R角优化),过渡尽可能地加大过渡部分的圆角,使过渡区接近于流线形,同时也要考虑到工艺性。一般可降低应力集中系数30%-50%。图18倒角对应力集中的影响2)合理优化开孔、开槽部位设计中应根据零件的受力方向和位置选择适当的开孔部位,孔一般应开在低应力区,如果必须开在高应力区,则应采取补强措施。椭圆形的长轴应与主应力方向平行,以降低应力集中系数。3)特定区域增加卸载槽,卸载孔,退刀槽在应力集中区附近的低应力部位增开缺口和圆孔这样可使应力的流线平缓,从而降低最大应力峰值,提高材料或零件的抗疲劳性能,如下图19所示。图19减小应力集中增加抗疲劳性能的一些设计4)尽量采用塑性材料采用塑性材料,即伸长率(>5%)和断面收缩率较大的材料,一方面塑性材料的塑性变形可降低材料的峰值应力,提供更高的安全系数。另一方面材料破坏前有比较明显的特征,方便识别,可为零部件的状态提供有效的判断和提前处理,如下图所示。常用的塑性材料包括:20、30、45号钢,Q235钢,5052、6061铝,SUS304、SUS316不锈钢,紫铜等。图20塑性和脆性材料的断裂特征5)合理进行表面处理表面处理可以改善基材表面平整度和粗糙度,能够有效地改善材料的抗疲劳能力,减少应力集中,从而提高零件的整体性能和可靠性。改善应力集中的表面处理技术包括:①表面强化技术:表面感应淬火、渗碳、渗氮和复合处理等,这些技术可以在表层形成高硬度和残余压应力,从而减少应力集中的影响。②薄壳淬火和喷丸强化,这些方法通过在零件表面形成硬化层和残余压应力,有效降低应力集中的危害。③表面形变强化:滚压强化和喷丸强化,这些处理使零件表面形变强化并产生残余压应力,从而降低应力集中的有害作用。图21一些表面处理后的零件--04--应力集中与应力奇异谈到应力集中我们一定会想到应力奇异问题,因为应力奇异的部位也是应力较大的部位,那么他们之间有什么区别和联系呢?应力奇异是由于几何形状的急剧变化(如尖角、裂纹尖端)或荷载、约束作用,在数值分析中应力趋向无穷大,尽管在实际应用中由于网格限制,应力值可能非常高的现象。应力集中和应力奇异的区别通过直角处的倒角和不倒角的处理解释,图22,更加直观。如下:图22应力集中和应力奇异的区别可以看出在有限元分析中,应力随着网格细化,应力一直增大,无法收敛,这种现象就是应力奇异。应力集中和应力奇异的区别在于应力是否随着网格加密收敛。即是否是应力收敛解,也称为网格无关解。应力奇异问题,其实是一种广义的应力集中。它是数值计算问题导致的。在复杂的几何形状中,数值计算问题导致了应力值异常。通常应力奇异的处理策略是通过细化网格或添加圆角等细节特征来缓解。其出现往往伴随着模型处理不当带来的人为高应力状态。因此在设计时,我们在处理模型时需要特别注意模型在点载荷、约束、凹角的处理,避免模型的过度约束和接触尖角等。关于模型的处理推荐小伙伴们去看看《使用AnsysMechanical获得精确数值结果(3)理解和处理人为高应力》这篇文章,一定会有所收获。好了,今天就和大家分享到这里,还请点赞、收藏和转发,感谢老铁支持!来源:薛定谔的Cube

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈