基于ANSYS的齿轮动态应力强度因子的计算.docx
- 文档编号:2170037
- 上传时间:2022-10-27
- 格式:DOCX
- 页数:10
- 大小:171.80KB
基于ANSYS的齿轮动态应力强度因子的计算.docx
《基于ANSYS的齿轮动态应力强度因子的计算.docx》由会员分享,可在线阅读,更多相关《基于ANSYS的齿轮动态应力强度因子的计算.docx(10页珍藏版)》请在冰豆网上搜索。
基于ANSYS的齿轮动态应力强度因子的计算
基于ANSYS的齿轮动态应力强度因子的计算
吴军荣,干洪,何芝仙
(安徽工程大学机械与汽车工程学院先进数控和伺服驱动技术安徽省
重点实验室,安徽芜湖241000
摘要
摘要:
:
以普通渐开线齿轮为研究对象,考虑时变刚度、相对阻尼,建立齿轮振动微分方程,求取齿轮传动齿面动载荷;并基于ANSYS软件平台,建立齿根含裂纹的齿轮模型,利用APDL语言加载齿面动载荷和求解裂纹动态应力强度因子,得出其变化规律。
关键词:
动载荷啮合时变刚度应力强度因子
CalculationofDynamicStressIntensityFactorforWheelgearbyFiniteElementSoftwareANSYS
WUJunrong,GANHong,HEZhixian
(AnhuiKeyLaboratoryofAdvancedNumericalControl&ServoTechnology,SchoolofMechanical&AutomobileEngineering,AnhuiPolytechnicUniversity,Wuhu241000,China
Abstract:
Withordinaryinvolutegearforthestudy,Consideringtime-varyingstiffness、Relativedamping,Buildingupawheelgearvibrationdifferentialequationtostrikegeartoothsurface’sdynamicload;andtheestablishmentofthegeartoothrootcrackmodelwithfiniteelementsoftwareANSYS,thenusesAPDLlanguageloadedthetoothsurfacedynamicloadandsolvesthedynamicstressintensityfactorwithcrack,atlastobtainstheirchangingpatterns.
Keywords:
DynamicloadEngagingtime-varyingstiffnessStressintensityfactor
0引言
齿轮传动是一种使用十分广泛的传动形式,具有传动效率高、传动平稳、结构紧
凑等优点。
由于轮齿承受复杂的变载荷作用,齿根弯曲强度不足以导致齿轮出现裂纹。
一方面,实际工作的齿轮出现裂纹往往不容易及时发现;另一方面,出现裂纹的齿轮仍具有一定的剩余寿命。
因此,研究具有裂纹的齿轮动态应力强度因子的计算,具有十分重要的工程意义。
本文以普通渐开线齿轮为研究对象,考虑沿啮合线变化的啮合时变刚度、恒定齿轮副相对阻尼,利用数值解法求解齿轮振动微分方程,得到一个啮合周期内的齿轮振动位移离散值,根据齿轮载荷力学模型求取齿轮传动齿面动载荷沿啮合线的变化规律
以及齿轮载荷作用角;并基于ANSYS有限元软件平台,建立齿根含裂纹的齿轮模型,利用APDL语言加载齿面动载荷和求解裂纹动态应力强度因子,得出齿轮单齿啮合过程中裂纹应力强度因子动态变化规律。
1齿轮传动齿面动载荷求解1.1齿轮机构的振动模型
从理论上说,渐开线齿轮传动具有平稳性。
但实际上由于齿轮的误差及轮齿啮合刚度的变化,可激发起齿轮的振动,从而产生噪音、激发机器
的振动。
图1
齿轮振动力学模型
齿轮传动系统具有明显的质量集中的特点,所以本文采用质量集中法,建立齿轮的动力学模型。
在忽略箱体、支承轴承和传动轴变形的前提下,齿轮的振动可简化为单自由度扭转振动模型如图1所示。
在不计齿轮齿廓误差的情况下,齿轮振动的微分方程为:
L
bbbmbbFrbrrtkbrr
22221122112(((bbmbbbL
JθcrθrθbktrθrθbrF++++=̇̇̇̇(2
式中,θ为齿轮弹性变形产生的扭转角;J、LF为齿轮转动惯量和法向载荷;br、b为
齿轮的基圆半径和齿宽;(mkt、c为
齿轮时变刚度和阻尼系数。
若假设两齿轮沿啮合线的相对位移为
2
211θθbbrrx+=(3将(3、(4两式合并得e
LememmF
xmtkxmtkx=++((2̇̇̇ζ(4
式中,em为齿轮的诱导质量
(2
122212
1bberJrJbJJm+=
;ζ为齿轮副的相对阻尼系数,取ζ=0.07。
齿轮时变刚度(((0
((
(zsmsZktktTtTktktTtT++≤≤⎧=⎨
≤≤⎩(5
11
60
ZTzn=
(6
(1
szTTε=−(7
式(5-(7中ZT为一个啮合周期所用的时间、sT为一个啮合周期内双齿啮合时间、1z为齿轮1齿数、1n为齿轮1转速、(kt为齿轮单对齿刚度。
计算齿轮单对齿刚度(kt有Weber弹性力学方法、石川法、有限元
法、边界元法等[3]
本文采用石川法进行计算,计算结果如图2所示。
图2k
(t随时间的变化曲线
图3km(t随时间的变化曲线
由(4式可知,齿轮的振动微分方程为二阶变系数微分方程。
其中含一个时变系数(tkm,即齿轮啮合时变刚度,
将(tkm表示为t的函数,然后利用数值迭代法(本文使用四阶——龙格库塔法求解(4式。
在一个啮合周期内,起点的振动参数与终点的振动参数相同,所以相对位移nx和相对速度nnx
v̇=的计算值应与初始值0x和00xv̇=作比较,若其差别较大,则将nx与nv代替初值重复计算。
直到nx与0x及0v与nv的差别足够小为止。
齿轮动态计算流程图如图所示:
图5齿轮动态计算流程图
渐开线齿轮的传动参数为:
齿数
1z=32,2z=96,模数m=2mm,齿宽
b=60mm,转速1n=1000n/min,压力角
�20=α,弹性模量111010.2E×=Pa,输
入转矩1T=114.6mN⋅。
图6为一个啮合周期内两齿轮沿
啮合线的相对位移,可以看出两齿轮沿啮合线的相对位移在整个啮合周期中都有波动,双齿啮合开始后,由于时变刚度的变化和阻尼的存在,波动的幅度逐渐变小,在双、单齿交替的时刻相对位移有了较大的突变,这与实际情况相符。
图6相对位移x随啮合时间的变化规律
图7为齿轮齿面动载荷随啮合时间的变化规律。
动载荷呈波动性变化,双、单齿交替致使动载峰值大幅增加。
变换解的整体波动幅度比实际解小,而且变换解在单齿啮合区未出现载荷
波动,只出现一个载荷高峰。
图7齿面动载荷随啮合时间的变化规律
2.动态应力强度因子求解
2.1建立齿根裂纹单齿有限元模型
建立单齿模型,在齿根处开深度为5mm的裂纹。
(1根据齿轮啮合关系确定啮合起始点,建模中从啮合起始点将啮合齿廓线分成啮合区(线和非啮合区(线;如图8为齿根裂纹单齿模型图。
(2在裂纹尖点处,利用KSCON命令对裂纹尖点进行处理,划分网格时自动生成1/4节点;
(3根据已确定的起始啮合点,在啮合区(线划分99个单元,为后文
加载动态啮合力、求解齿根裂纹应力强度因子作准备;如图9所示为有限元网格模型。
(4
边界条件的处理?
图8
齿根裂纹单齿模型
图9齿根裂纹单齿有限元模型
2.2APDL语言编程处理力的加载和计算
利用APDL语言编程在齿轮啮合区动态加载不同时刻的啮合力,某一时刻的力对应着一节点,因此须将啮合区中节点进行编号,以便于动态啮合力的加载,将起始啮合点到分离点全部节点进行编号,数据存入数组;同时,前文数值计算的啮合力为矢量,加载中将力分解到X、Y两个方向上进行加载。
最后,进行求解并计算单齿啮合过程中动态裂纹应力强度因子。
如图10所示,为齿根裂纹动态应力强度因子计算流程图。
图10动态应力强度因子计算流程图
图11为齿轮啮合过程齿根裂纹动态应力强度因子变化情况。
可以看出:
(1在啮合过程中,齿根裂纹应力强度因子值随啮合点的改变而不断变化的;
(2最大值出现在单齿啮合区,与动态啮合力最大值出现在单齿啮合区相对应;
(3在单齿啮合区后的双齿啮合区段,啮合力呈减小趋势,但力臂不断增大,此时的齿根裂纹应力强度因子值不断增大;在齿轮啮合过程中,齿根裂纹应力强度因子的动态变化趋势和啮合力的变化趋势基本相
同。
图11动态应力强度因子变化情况
3.结果和讨论
(1齿轮齿面动载荷整体呈现出较大波动,在双齿与单齿交替的时刻有较大突变,作用在齿面上的载荷作用角也逐渐增大。
因此齿轮动力学效应会对齿轮裂纹动态强度因子有较大影响。
(2在齿轮啮合过程中,齿根裂纹动态强度因子随时间呈变化趋势;其最大值出现在齿轮啮合区的单齿啮合区,且整体变化趋势与动载荷变化趋势一致。
参考文献:
[1]王玉新,柳杨,王仪明.考虑啮合时变刚度和传递误差的齿轮振动分析[J].机械传动,2002,26(1:
5-8
[2]祝效华,余志祥等.ANSYS高级工程有限元分析范例精选[M].北京:
电子工业出版社,2005:
411-421.
[5]季维英,陈荣.基于ANSYS的应力强度因子计算[J].南通职业大学学报,2008,22(2,92-94.
[6]陈银强.圆柱体表面椭圆裂纹应力强度因子的有限元分析和研究[D].武汉:
武汉工程大学,2006.
[7]钟小强,蒋维,汪磊磊.变速传动下的齿轮啮合力计算仿真研究[J].机械设计,2008,25
(1):
42-44[8]卫丽,亓秀梅.齿轮传动动载荷沿啮合线分布规律的分析[J].机械工程与自动化,2006,(6):
63-66[9]ZhiwenGao,YouheZhou.Fracturebehaviorforalongcenterslantcrackedrectangularslabofsuperconductorunderelectromagneticforce[J].ModernPhysicsLettersB,2009Vol:
23No:
6pp.825-834.[10]季庆芬.断裂力学及其工程应用[M].哈尔滨:
哈尔滨工程大学出版社,2007:
67-68.97:
100
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 ANSYS 齿轮 动态 应力 强度 因子 计算
![提示](https://static.bdocx.com/images/bang_tan.gif)