ansys实例-正确地模拟过盈配合

合集下载

基于ANSYS的电能表过盈配合接插件稳态热分析

基于ANSYS的电能表过盈配合接插件稳态热分析

基于ANSYS的电能表过盈配合接插件稳态热分析陈盛;徐英辉;张蓬鹤;袁翔宇;井天军【摘要】为提高智能电能表现场安装效率,同时确保安装可靠性,国家电网标准Q/GDW 11008-2013《低压计量箱技术规范》明确了插拔式电能表结构及其过盈配合接插件方案,并且进一步规定了接插件温升参数.目前过盈配合接插件在初始设计阶段重点考虑结构形式,并没有统筹考虑温升因素,温升参数均采用成品实验测定.为此,建立了详细的电能表插拔式安装结构及其过盈配合接插件热传导模型,并采用ANSYS有限元仿真软件对其进行稳态热分析,为过盈配合接插件结构进行温升参数仿真验证和理论设计,避免成品实验不合格导致重新设计.【期刊名称】《电测与仪表》【年(卷),期】2018(055)020【总页数】6页(P96-101)【关键词】电能表;过盈配合接插件;温升;ANSYS;稳态热分析【作者】陈盛;徐英辉;张蓬鹤;袁翔宇;井天军【作者单位】中国电力科学研究院有限公司,北京100192;中国电力科学研究院有限公司,北京100192;中国电力科学研究院有限公司,北京100192;中国电力科学研究院有限公司,北京100192;中国农业大学,北京100083【正文语种】中文【中图分类】TM9330 引言随着我国智能电网建设进程的快速推进以及分时阶梯电价政策的逐步实施,传统电能表已无法满足电力信息采集的多样化功能需求。

为此,国家电网公司大力开展智能电能表推广应用,以此推动我国智能电网建设,同时更好地承担计量、计费等客户营销服务工作。

在此背景下,新装或者传统电能表换装智能电能表数量将持续增加。

电能表现场安装作业过程中,传统采用导线直连电能表尾端接线端子并利用螺钉压紧的方法,由于施工条件限制、作业人员能力素质等原因,往往存在以下问题[1-2]:(1)现场施工工艺导致的接线质量和安装可靠性问题:压线螺钉压力不够致使导线松动,压力过大导致压断导线;剥线长度过短,压线螺钉压在绝缘外皮导致接触不良,剥线过长则电能表尾端导线外露,造成触电事故;(2)施工人员素质能力参差不齐,现场接线可能造成接线错误,导致计量故障,引发纠纷和后续排查改线工作;(3)现场安装施工效率低,以单相电能表为例,需分别压接火线进出线和零线进出线,操作繁琐复杂,往往导致工作效率较低,减缓工程进度;(4)电能表更换作业过程中,需要对用户停电后拆除旧表并更换新表,导致客户停电。

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 normaland tangent - 法向和切向均采用拉格朗日乘子算法。

各种不同的约束算法各有其优缺点,各有各自最适用的场合,具体情况需要具体对待。

大部分情况下,默认选择KEYOPT(2)=0 就够用了。

过盈配合所致的接触分析的难点在于如何确定初始接触状态。

初始接触状态设置得不对,会导致错误的计算结果或者不准确的计算结果,下面举两个例子来说明。

ANSYS仿真计算代做:模态分析,瞬态动力学,谐响应分析和谱分析、械结构的疲劳、损伤,CFD流体;结构的强度评估和优化;企鹅:690294845例1.两个圆柱体在几何上是刚好接触,划分网格后有限元模型有间隙。

基于Ansys Workbench的行星架组件过盈配合仿真研究及结构改善

基于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 引言减速器是机械零部件中的重要成分,在我国的装备与制造业中发挥着重要作用。

Ansys_过盈配合接触设置

Ansys_过盈配合接触设置
(下面这个例子实际上是一本 ANSYS 书上的一个例子,这个例子的 PDF 版本在 网上流传甚广,但是原书上的分析结果是错误的,具体错误之处,将在后面提及)。
例 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 Workbench对圆柱面过盈配合接触应力的研究

w i d t h. t h e i n l f u e n c e f a c t o r s o f a x i a l t e n s i o n a n d s t r e s s .An a l y s i s o f e x p e r i me n t a l r e s u l t s t h e f o l l o wi n g c o n c l u s i o n s :s h a t f a n d
p r o c e s s , b y c h a n g i n g t h e s i n g l e f a c t o r e x p e r i me n t me t h o d, a n a l y s i s o f i n t e f r e r e n c e i f t o f c y l i n d e r i s c o n t a c t s t r e s s d i s t r i b u t i o n
Abs t r ac t : Us i n g f i n i t e e l e me n t a na l y s i s s o f t wa r e ANSYS Wo r k be n c h, t h e s i mu l a t i o n o f t h e ho l l o w s ha f t a nd c o n t a i ni n g a
关键词 : 过盈配合 ; 有 限元 ; 影 响 因素 ; 等效应力
中图 分 类 号 : T H 1 4 0
Re s e a r c h 0 n c o nt a c t s t r e s s 0 f c y l i n de r i n t e r f e r e n c e it f ba s e d

Ansysworkbench参数化过盈配合模拟计算确定所需过盈量校核计算

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对圆柱面过盈配合接触应力的研究李亚洲;汤易升;陈炜;张西正【摘要】利用有限元分析软件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组合模具过盈配合有限元分析
由图7可知,我们可以清楚的看到接触面上在过盈 配合中的压力变化和分布,以及最大压力发生的位置。便 于我们了解过盈配合的动态过程,这是传统计算方法所 达不到的。
利用时间历程后处理器分析100端面处某一节点的 约束反力随时间变化曲线。如图8所示。
3结语 通过采用ANSYS有限元软件的分析,该组合模具的
过盈配合过程中产生的最大应力为1240MPa,没有超出 材料的许用应力(1600MPa),表明设计没有使产品产生
基于ANSYS组合模具过盈配合有限元分析
郝兆朋,张永军,范依航 (哈尔滨理工大学机械动力工程学院,哈尔滨150080)
摘要:在细长杆冷成形自动锻压机新产品的研制过程中,为判断其组合模具结构设计是否符合要求,文中用接触有限 元法建立了内衬与外套过盈配合的有限元力学模型。针对内衬和外套的过盈配合状态。采用大型通用有限元ANSYS软
比较方案一与方案二,对于所给零件,由于两小孔比
较接近边缘,复合模冲裁零件时受到壁厚的限制,模具结 构与强度方面相对较难实现和保证,故采用级进模加工。 2模具设计计算 2.,排样、计算条料宽度及确定步距
采用单排方案,如图2所示。
查表确定搭边值,根据零件形状两零件间按矩形取 搭边值b=2.Omm,侧边取搭边值a=2.Omm。
万方数据
机械工程师2008年第5期《27
专题报道 譬墨蚕蕊圈瞄冒誓墨赣县加工菝7It与设备■●—●—■—●————■———●—■—■—■—————■—■—■■—●——■——————一
垫圈连续冲栽模具的设计
李文 (韶关机电技校,广东韶关512000)



”。
j秀一
摘要:分析了垫圈的工艺特点,介绍了垫圈的排样设计方案,以及模具结构的设计计算、模具主要零件计算、冲床参数

Ansys_过盈配合接触设置

Ansys_过盈配合接触设置
KEYOPT(2)=4 Pure Lagrange multiplier on contact normal and tangent 法向和切向均采用拉格朗日乘子算法。 各种不同的约束算法各有优缺点,各有各自的最适用的场合,具体情况具体对待。 大部分情况下,默认选择 KEYOPT(2)=0 就够用了。 过盈配合所致的接触分析的难点在于如何确定初始接触状态。初始接触状态 设置得不对,会导致错误的计算结果或者不准确的计算结果,下面举 2 个例子来 说明。 例 1.两个圆柱体在几何上是刚好接触,划分网格后有限元模型有间隙。 如图 1 所示。
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的汽轮发电机转轴的过盈配合分析黄鹏,倢任,张金华(中国长江动力集团有限公司,武汉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接触实例分析

1.实例描述一个钢销插在一个钢块中的光滑销孔中。

已知钢销的半径是0.5 units, 长是2.5 units,而钢块的宽是4 Units, 长4 Units,高为1 Units,方块中的销孔半径为0.49 units,是一个通孔。

钢块与钢销的弹性模量均为36e6,泊松比为0.3.由于钢销的直径比销孔的直径要大,所以它们之间是过盈配合。

现在要对该问题进行两个载荷步的仿真。

(1)要得到过盈配合的应力。

(2)要求当把钢销从方块中拔出时,应力,接触压力及约束力。

2.问题分析由于该问题关于两个坐标面对称,因此只需要取出四分之一进行分析即可。

进行该分析,需要两个载荷步:第一个载荷步,过盈配合。

求解没有附加位移约束的问题,钢销由于它的几何尺寸被销孔所约束,由于有过盈配合,因而产生了应力。

第二个载荷步,拔出分析。

往外拉动钢销1.7 units,对于耦合节点上使用位移条件。

打开自动时间步长以保证求解收敛。

在后处理中每10个载荷子步读一个结果。

本篇先谈第一个载荷步的计算。

下篇再谈第二个载荷步的计算。

3.读入几何体首先打开ANSYS APDL然后读入已经做好的几何体。

从【工具菜单】-->【File】-->【Read Input From】打开导入文件对话框找到ANSYS自带的文件(每个ansys都自带的)\Program Files\Ansys Inc\V145\ANSYS\data\models\block.inp 【OK】后,四分之一几何模型被导入。

4.定义单元类型只定义实体单元的类型SOLID185。

至于接触单元,将在下面使用接触向导来定义。

5.定义材料属性只有线弹性材料属性:弹性模量36E6和泊松比0.36.划分网格打开MESH TOOL,先设定关键地方的网格划分份数然后在MESH TOOL中设定对两个体均进行扫略划分,在volumeSweeping中选择pick all,按下【Sweep 】按钮,在主窗口中选择两个体,进行网格划分。

盘轴过盈配合的ANSYS分析

盘轴过盈配合的ANSYS分析
(3)建立几何模型并划分网格ANSYS提供了 强大的模型生成功能,可以方便地建立以上数学模 型。也可以用Pro/E、Solidworks、AutoCAD等软件建 立实体模型,导人到ANSYS中进行网格划分得到有 限元模型,如图2所示。
(4)创建接触对在Contact Manager(接触管理 器)中分别将圆盘的盘心面轴的外环面选中,在Coef- ficient of Friction(摩擦系数)文本框中输入“0.2”。 将Normal Penalty Stiffness(正则处罚刚度)定义为 “0.1”将Stiffness matrix(刚度矩阵)下拉框设定为 “Unsymmetric”选项。创建的接触对,如图3所示。
然后施加位移载荷(将轴沿轴向平移40ram拔出 盘孔)。选取菜单路径Utility menu I Select I Enti- ties。单击对话框中最上面下拉框中的“Node”,指定 选择对象为节点。接着在下面的下拉框中单击“By location”,指定选择方式为根据坐标值来选取。单击
万方数据
图7接触单元上的压应力分布云图
4计算结果的提取分析pJ (I)利用通用后处理器观察结果 首先将建立
的I/4模型扩展为整个盘轴结构模型,然后利用通用 后处理器来提取求解结果。依次提取过盈配合应力 分布云图如图6所示,接触单元上的压力分布云图如 图7所示。
荷子步数(Number of substeps)为1,关闭自动时间
步,其求解过程时间曲线,如图4所示。
本文读者也读过(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)

Ansys WB轴孔过盈配合仿真分析(转载)

Ansys WB轴孔过盈配合仿真分析(转载)

Ansys WB轴孔过盈配合仿真分析ANSYS有限元仿真 3月10日319本期通过一个简单实例来介绍怎么通过ANSYS Workbench有限元分析工具来模拟过盈配合。

问题描述:一个圆环与空心轴过盈连接,过盈量2mm(半径方向)。

求圆环与轴的最大变形,和最大等效应力。

1分析模型模型通过WB DM(DesignModeler)绘制。

如下图所示。

从外到内,直径依次是100、70、45mm。

拉伸厚度为18、10mm。

材料钢。

注意建模时:圆环的内径等于空心轴的外径。

2分析流程图添加了几何模型和静力学两个模块。

3接触静力分析模块Static Structural会自动将所有低于某一特定值的间隙识别为绑定接触(bonded),这需要自己手动修改接触类型为摩擦接触(Frictional);Contact和target分别为空心轴的外表面和圆环的内表面;然后添加摩擦系数这里设为0.2。

然后输入过盈量offset 2mm。

(过盈输入正值,间隙输入负值) 接触设置如下图所示。

4网格这个几何模型外形很规范可以将模型网格划分得规则一些,以便提高计算准确性。

方法:对这两零件使用扫描方式划分,发现效果并不满意,然后再次划分选择两零件的边线进行等分80份。

设置如下图。

划分效果如下:注:由于该模型单元数量少,并不会给计算机造成多大的负担,而且节点排列规则有助于提供精度。

所以这里使用六面体网格单元是合理的。

5边界条件将空心轴内表面固定约束。

6求解结果总变形云图如下等效应力云图如下通过隐藏零部件很容易查看到单个零件的力和变形。

结论如下:空心轴最大变形为0.236mm 最大应力为5836MPa。

圆环最大变形1.885mm,最大应力11478MPa。

其中圆环变形加上空心轴变形之和为2.121mm与2mm 很接近。

复合变形协调条件,仿真结果可信。

另外还计算了两种情况:①过盈量3mm,验证了圆环和轴的变形都变大并且两者之和接近3mm。

②将offset设置为-2仿真间隙配合情况,结果是变形应力均为零,符合理论情况。

基于ansys的过盈配合接触应力分析

基于ansys的过盈配合接触应力分析

基于ansys的过盈配合接触应力分析摘要介绍了基于ansys的接触分析步骤,并通过ansys软件,将对一个盘轴紧配合结构进行接触分析,来说明接触分析的有限元计算方法。

关键词ansys 过盈配合接触分析引言在工程结构中,经常会遇到大量的接触问题。

火车车轮与钢轨之间,齿轮的啮合是典型的接触问题。

接触问题是一种高度非线性行为,需要较大的计算资源,为了进行实为有效的讣算,理解问题的特性和建立合理的模型是很重要的。

接触问题存在两个较大的难点:其一,在你求解问题之前,你不知道接触区域,表面之间是接触或分开是未知的,突然变化的,这随载荷、材料、边界条件和其它因素而定;其二,大多的接触问题需要计算摩擦,有儿种摩擦和模型供你挑选,它们都是非线性的,摩擦使问题的收敛性变得困难。

本文以ansys软件为工具,以某转子中轴和盘的连接为例,分析轴和盘的过盈配合的接触应力。

1. 面面接触分析的步骤:在涉及到两个边界的接触问题中,很自然把一个边界作为“LI标”面而把另一个作为“接触”面,对刚体一柔体的接触,“目标”面总是刚性的,“接触”面总是柔性面,这种两个面合起来叫作“接触对”。

使用Targel69和Contdl71或Contal72来定义2D接触对,使用TargelYO和Contal73或Contal74来定义3D接触对,程序通过相同的实常数号来识别“接触对”。

在接触问题中,两个相互接触的物体必须满足边界不穿透的约束条件,施加边界不穿透约束的方法主要有罚函数算法和扩增的拉格朗日算法。

罚函数算法是在总势能泛函中加入惩罚项,来近似满足接触约束条件。

从物理意义上讲,罚函数法相当于在接触边界上加入线弹簧以防止接触面之间的相互渗透,而罚函数因子相当于弹簧的刚度系数。

罚函数法的优点在于不增加系统未知数总数,可保持刚度矩阵的对称性,提高了求解效率,但罚函数因子的取值对计算结果的精度影响很大,必须根据渗透情况对其进行多次调整。

扩增的拉格妙日算法是为了找到精确的拉格妙日乘子而对罚函数修正项进行反复迭代,与罚函数的方法相比,拉格朗日方法不易引起病态条件,对接触刚度的灵敏度较小,然而,在有些分析中,扩增的拉格朗日方法可能需要更多的迭代,特别是在变形后网格变得太扭曲时。

ansys实例-正确地模拟过盈配合

ansys实例-正确地模拟过盈配合

过盈配合在‎机械产品的‎装配中使用‎的相当普遍‎。

比如轴与轴‎承、轴与轴瓦、汽车的制动‎盘等,都是通过一‎定的过盈量‎来使两个装‎配部件紧密‎地连接起来‎。

下面讨论如‎何在 ANSYS‎中正确地模‎拟过盈配合‎。

过盈配合在‎有限元分析‎中是一种典‎型的非线性‎接触行为。

在有限元分‎析中设定了‎接触,从本质上来‎讲就是对相‎互接触的两‎个部件施加‎了某种约束‎,不同的接触‎算法对于接‎触约束的处‎理方法有所‎不同。

接触约束的‎理论算法的‎选择,在 ANSYS‎中是通过设‎置 conta‎c t 单元的KEOPT‎(2) 选项来实现‎的。

在 ANSYS‎中目前主要‎有5 种接触约束‎算法:KEYOP‎T(2)=0 Augme‎n ted Lagra‎n gian‎- 加强的拉格‎朗日算法,这是 ANSYS‎的缺省选择‎;KEYOP‎T(2)=1 Penal‎t y funct‎i on - 罚函算法;KEYOP‎T(2)=2 Multi‎p oint‎const‎r aint‎(MPC) - 多点约束算‎法;KEYOP‎T(2)=3 Lagra‎n ge multi‎p lier‎on conta‎c t norma‎l and penal‎t y on tange‎n t -接触法向采‎用拉格朗日‎乘子,接触切向采‎用罚函数的‎综合算法。

KEYOP‎T(2)=4 Pure Lagra‎n ge multi‎p lier‎on conta‎c t norma‎l andtange‎n t - 法向和切向‎均采用拉格‎朗日乘子算‎法。

各种不同的‎约束算法各‎有其优缺点‎,各有各自最‎适用的场合‎,具体情况需‎要具体对待‎。

大部分情况‎下,默认选择K‎E YOPT‎(2)=0 就够用了。

过盈配合所‎致的接触分‎析的难点在‎于如何确定‎初始接触状‎态。

初始接触状‎态设置得不‎对,会导致错误‎的计算结果‎或者不准确‎的计算结果‎,下面举两个‎例子来说明‎。

HyperMesh与ANSYS过盈接触联合仿真 1

HyperMesh与ANSYS过盈接触联合仿真 1
HyperMesh与ANSYS过盈接触联 合仿真(1)
内容介绍
接触非线性分析实例中的过盈接触问题仿真,本例主要讲解在HyperMesh部分的设置。
方法/步骤
1)将轴对称的模型进行四分之一对 称处理,只保留四分之一部分进行格 划分:
方法/步骤
左图表示了本步骤的内容,以图片内 容为准。
方法/步骤
左图表示了本步骤的内容,以图片内 容为准。
方法/步骤
2)点击F12,对两个物体分别进行格 剖分:
方法/步骤
左图表示了本步骤的内容,以图片内 容为准。
方法/步骤
3)定义单元和材料属性,并创建接 触对,点击Tool>contact manager>new>若采用几何穿透模拟过 盈,设置关键字如下:K5=0(不进行 自动调整初始间隙和穿透)K9=2(考 虑单元的初始接触和穿透),同时, 可以设置K10=2,允许ANSYS自动更新 接触刚度。
方 容为准。
方法/步骤
4)导出,并进入ANSYS,查看接触对定义情况。
谢谢观看

在ANSYS中怎样精确地设置过盈量

在ANSYS中怎样精确地设置过盈量

在ANSYS中怎样精确地设置过盈量下面讨论如何在ANSYS中正确地模拟过盈配合。

过盈配合在有限元分析中是一种典型的非线性接触行为。

在有限元分析中设定了接触,从本质上来讲就是对相互接触的两个部件施加了某种约束,不同的接触算法对于接触约束的处理方法有所不同。

过盈配合所致的接触分析的难点在于如何确定初始接触状态。

初始接触状态设置得不对,会导致错误的计算结果或者不准确的计算结果,下面举2个例子来说明。

例1.两个圆柱体在几何上是刚好接触,划分网格后有限元模型有间隙。

如图1所示。

图1. 两个在几何上刚好接触的圆柱体这两个圆柱体,在几何上是刚好相切的,即处于几何上刚好接触的初始状态。

划分网格后,由于在圆周上用小段直线代替了弧线,两个圆柱体之间产生了一定的间隙,两个圆柱体的有限元的初始状态不再是有接触的,此时,如果接触参数设置不当,就会因为初始约束不足,圆柱体出现刚体位移,得到错误的结果。

例2.有的人把两个接触部件的几何位置设定一定的过盈量,想用这个过盈量来模拟过盈配合,这种做法是错误的,几何上的过盈量不等于划分网格后有限元模型的实际过盈量。

下面的图2中,是一个孔类零件和一个轴类零件的截面图,轴和孔在几何位置上预设了过盈量。

(内圈的红色圆是孔边界,外圈的蓝色圆是轴边界,轴和孔在几何上是相互侵入的)。

图2:一个轴类零件和孔类零件的过盈配合的截面图在几何上,图2的轴和孔有一定的过盈配合量,其大小等于两个圆的半径之差,我们的本意是想用这个几何位置上的过盈量来模拟过盈配合。

不幸的是,两个部件划分网格之后,实际的过盈量应该为单元之间的距离,即图上中靠得比较近的两条线段之间的距离,显然,这个距离不再等于我们预先设置的过盈量了。

更何况,上面这个图还是两个部件的网格对对得比较整齐的情况,如果网格对的不整齐,过盈量就和我们预设的差的更远了。

对于过盈配合来讲,过盈量的数值变化对于过盈产生的应力的影响是很大的。

上面两种方法都不能精确地设置过盈量,下面介绍第3种方法。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

过盈配合在机械产品的装配中使用的相当普遍。

比如轴与轴承、轴与轴瓦、汽车的制动盘等,都是通过一定的过盈量来使两个装配部件紧密地连接起来。

下面讨论如何在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 normaland tangent - 法向和切向均采用拉格朗日乘子算法。

各种不同的约束算法各有其优缺点,各有各自最适用的场合,具体情况需要具体对待。

大部分情况下,默认选择KEYOPT(2)=0 就够用了。

过盈配合所致的接触分析的难点在于如何确定初始接触状态。

初始接触状态设置得不对,会导致错误的计算结果或者不准确的计算结果,下面举两个例子来说明。

ANSYS仿真计算代做:模态分析,瞬态动力学,谐响应分析和谱分析、械结构的疲劳、损伤,CFD流体;结构的强度评估和优化;企鹅:690294845例1.两个圆柱体在几何上是刚好接触,划分网格后有限元模型有间隙。

如图1 所示。

这两个圆柱体,在几何上是刚好相切的,即处于几何上刚好接触的初始状态。

划分网格后,由于在圆周上用小段直线代替了弧线,两个圆柱体之间产生了一定的间隙,两个圆柱体的有限元模型的初始状态不再是接触的。

此时,如果接触参数设置不当,就会因为初始约束不足,圆柱体出现刚体位移,得到错误的结果。

(说明:例1 本来与设置过盈量是无关的,为了说明初始接触状态的重要性顺带说说。

)例2.有的人把两个接触部件的几何位置设定一定的过盈量,想用这个过盈量来模拟过盈配合,这种做法是错误的,几何上的过盈量不等于划分网格后有限元模型的实际过盈量。

下面的图2 中,是一个孔类零件和一个轴类零件的截面图,轴和孔在几何位置上预设了过盈量。

(内圈的红色圆是孔边界,外圈的蓝色圆是轴边界,轴和孔在几何上是相互侵入的)。

在几何上,图2 的轴和孔有一定的过盈配合量,其大小等于两个圆的半径之差,我们的本意是想用这个几何位置上的过盈量来模拟过盈配合。

不幸的是,两个部件划分网格之后,实际的过盈量应该为单元之间的距离,即图中靠得比较近的两条线段之间的距离,显然,这个距离不再等于我们预先设置的过盈量了。

更何况,上面这个图还是两个部件的网格对得比较整齐的情况,如果网格对得不整齐,过盈量就和我们预设的差得更远了。

对于过盈配合来讲,过盈量的数值变化对于过盈产生的应力的影响是很大的。

在ANSYS 中,要正确的设置过盈配合,主要分3 步:第一步:设置KEYOPT(9) = 4KEYOPT(9) 的默认值为0,意思是既考虑两个接触部件由于初始几何位置造成的初始侵入量(或者间隙),同时也考虑CNOF 参数设置的偏移量。

意即接触部件的初始接触状态是由CNOF 和初始侵入量(或间隙) 共同决定的。

在这种情况下,两个接触部件的初始几何位置对初始接触状态是有影响的,这对于准确设置过盈量是很不利的。

(前面例2 已经说明了通过几何位置设置初始过盈量是不准确的)。

在设置了KEOPT(9) = 4 之后,程序在计算初始接触状态的时候就只考虑CNOF 的设置值,不考虑接触部件的几何位置造成的侵入或间隙,而且过盈量是以ramp 方式施加的。

(ramp 施加方式即逐步施加)。

第二步:通过设置Icont 实常数划分网格后,通常情况下,Target surface 和contact surface 上的单元之间会有间隙或者过盈量,如果间隙或者过盈量在Icont 设定误差范围内,间隙或者过盈量会被消除掉,程序会使contact surface 和target surface 上的单元处于刚好接触的状态。

这个值的的具体设置可以参见帮助文档,本文中设置为0.2。

第三步:通过设置实常数CNOF 来设置过盈量在第二步中,通过Icont 的设置,已经使得Contact surface 上的单元和Target surface 上面的单元处于刚好接触的位置了,此时再设置CNOF,CNOF 的值就是过盈量。

(CNOF 的本意并不是过盈量,只是在有了前面的设定后,它的值就是我们所要的过盈量,其具体含义请参考ANSYS 的帮助文档)。

下面这个例子实际上是一本ANSYS 书上的一个例子,这个例子的PDF 版本在网上流传甚广,但是原书上的分析结果是错误的,具体错误之处,将在后面提及。

例3.一个简单的轴和带孔圆盘的过盈配合的实例。

圆盘的基本尺寸为:内径Rpin = 35 mm (原书中此值为34 mm),外径Rpout = 100 mm,盘高Hp = 25mm;轴的基本尺寸为:内径Rain = 25 mm,外径Raout = 35 mm,轴长La = 150 mm。

(原书中圆盘孔内径为Rpin =34 mm,和轴在几何上形成1 mm 的过盈量。

由于结构是完全轴对称的,故可只取四分之一模型分析之。

本例分析中,取过盈量f = 0.01 mm,而且本例仅仅计算由于过盈配合所产生的应力。

按照本例各个物理量所取的单位,最终的计算结果中,应力单位应该为MPa;完整的命令流如下:Finish/clear,start/TITLE,Contact analysis with initial interference/PREP7! 带孔圆盘的基本尺寸;Rpin=35Rpout=100Hp=25! 轴的基本尺寸;Rain=25Raout=35La=150! 过盈量f;f=0.01! 实体的单元类型为带中间节点的2阶六面体单元;ET,1,solid186MP,EX,1,2.1E5 ! 弹性模量;MP,PRXY,1,0.3 ! poisson系数;! 生成带孔圆盘的1/4实体模型;CYL4,0,0,Raout,0,Rpout,90,Hp! 轴的1/4 实体模型;CYL4,0,0,Rain,0,Raout,90,La! 把轴的位置沿着轴向移动一段距离;VGEN, ,2, , , , ,-10, , ,1!!!!************对实体划分网格****************************!*LESIZE,17, , ,15, , , , ,1LESIZE,19, , ,15, , , , ,1!*LESIZE,18, , ,2, , , , ,1LESIZE,20, , ,2, , , , ,1!*LESIZE,22, , ,20, , , , ,1!*LESIZE,5, , ,10, , , , ,1LESIZE,7, , ,10, , , , ,1!*LESIZE,6, , ,8, , , , ,1LESIZE,8, , ,8, , , , ,1!*LESIZE,10, , ,3, , , , ,1!*VSWEEP,ALL!*************网格划分完毕**********************************!!/COM, CONTACT PAIR CREATION - STARTMP,MU,1,0.2MAT,1R,3REAL,3ET,2,170ET,3,174R,3,,,0.2,0.2,0.9,0 ! Icont 系数设置为0.2;! 设置过盈量为f;RMORE,,,1.0E20,f,1.0,0RMORE,0.0,0,1.0,,1.0,0RMORE,0,1.0,1.0,0.0,,1.0KEYOPT,3,4,0KEYOPT,3,5,0NROPT,UNSYMKEYOPT,3,7,0KEYOPT,3,8,0!KEYOPT(9)=4;! 不考虑初始几何位置造成的过盈或者间隙, ! 只考虑CNOF 参数设置的值,即过盈量。

KEYOPT,3,9,4KEYOPT,3,10,2KEYOPT,3,11,0KEYOPT,3,12,0KEYOPT,3,2,0KEYOPT,2,5,0! Generate the target surfaceASEL,S,,,4CM,_TARGET,AREATYPE,2NSLA,S,1ESLN,S,0ESLL,UESEL,U,ENAME,,188,189ESURFCMSEL,S,_ELEMCM! Generate the contact surfaceASEL,S,,,9CM,_CONTACT,AREATYPE,3NSLA,S,1ESLN,S,0ESURFALLSELESEL,ALLESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3/PSYMB,ESYS,1/PNUM,TYPE,1/NUM,1EPLOTESEL,ALL/COM, CONTACT PAIR CREATION - ENDFINISH/SOLFLST,2,4,5,ORDE,4FITEM,2,5FITEM,2,-6FITEM,2,11FITEM,2,-12DA,P51X,SYMM ! 设置对称约束;FLST,2,1,5,ORDE,1FITEM,2,3/GODA,P51X,ALL,!*ANTYPE,0NLGEOM,1 ! 必须打开大变形效应;NSUBST,1,0,0AUTOTS,0TIME,100/STATUS,SOLUSOLVEFINISH/POST1SET,1,LAST,1,PLNSOL,s,eqv ! 查看单元的Von Misesstress;esel,s,ename,,174PLNSOL, CONT,PRES, 0,1.0 ! 在contact 单元上查看contactpressure;SAVE计算结果:说明:修改f 值即可修改过盈量。

如果过盈量设置的过大,接触部件可能会产生塑性变形,如果要考虑塑性变形,则应该修改材料属性定义为塑性材料。

相关文档
最新文档