e6条件对斜齿轮动力学模态分析的影响
第32卷第1期2015年3月上海第二工业大学学报
JOURNALOFSHANGHAISECONDPOLYTECHNICUNIVERSITYVol.32No.1
Mar.2015
文章编号:1001-4543(2015)01-0046-05
基于ANSYS边界约束条件对斜齿轮动力学
模态分析的影响
刘赛,吴飞科,严剑刚,杨国策
(上海第二工业大学工程训练中心,上海201209)
摘
要:在CATIA和ANSYS软件环境下,分别建立了斜齿圆柱齿轮的三维几何模型和动力学模态分析有限元模
型,进而对一定转速、不同边界约束条件下的齿轮进行了模态分析,得到了各阶固有频率和振型。结果表明:不同的边界条件对斜齿轮模态分析影响较大;目前通常采用的边界约束条件与实际情况存在差异,固有频率和振型图均会发生变化;为提高斜齿轮模态分析的精度,应根据实际工作状态确定边界约束条件。关键词:斜齿轮;模态分析;边界约束条件;固有频率;振型中图分类号:TH132.413
文献标志码:A
0引言
齿轮作为最重要的机械传动部件之一,被广泛应用在各类工程领域中。齿轮工作时受到周期性载荷的作用,可能发生强烈的共振,而共振会引起很大的激振力,这些激振力与外界激振力的相互作用,会导致齿轮损坏甚至报废。当今工程中对齿轮动态性能的要求越来越高,对齿轮传动机构的运动平稳性、疲劳强度乃至噪音的要求愈加严格。近年来,国内外一些学者对齿轮动力学性能进行了深入的研究
[1-9]
1ANSYS模态分析的理论基础
模态分析可以确定一个结构的固有频率和
振型,它们是承受动态载荷结构设计中的重要参数。ANSYS模态分析是线性分析,任何非线性特性,即使定义也将被忽略[9]。其数学模型可以由达朗贝尔原理推出齿轮系统的运动微分方程:
¨+CX˙+KX=F(t)MX
X=[X1X2···Xn]TF(t)=[F1F2···Fn]T
式中:M,C,K分别是齿轮系统质量矩阵、阻尼矩
¨,X˙,X分别为齿轮系统振动加速阵和刚度矩阵;X
度向量、速度向量和位移向量;F(t)为齿轮所受外界激振力向量。
假设不受外力影响,则F(t)=0。在实际工程中,阻尼对结构自振频率和振型影响不大,可忽略阻尼力,得到无阻尼自由振动的运动方程:
¨+KX=F(t)MX
(2)(1)
。在齿轮传动系统中,分析其模态,有益于掌
握齿轮的振动性能,得到齿轮结构的固有频率,在设计工程中可有效避开这些频率或最大程度上减少对这些频率的激励,从而达到消除过度振动或降低噪声、提高齿轮寿命的目的[5]。然而,在目前的齿轮模态分析中,往往分析其静态模态且在齿轮轮毂内圈施加全约束,这不符合实际工作情况。因此,本文运用有限元软件ANSYS对斜齿圆柱齿轮进行不同约束条件的动态模态分析,从而分析一定转速下的固有频率、振型与边界约束的关系,力求为齿轮动态设计提供参考依据。
收稿日期:2014-12-31
通讯作者:刘赛(1983–),男,河北人,讲师,硕士,主要研究方向为机械结构有限元及优化分析、逆向工程等。电子邮箱
其对应的特征值方程为
(K−ω2
iM)Xi=F(t)
(3)
式中,ωi为第i阶模态的自振频率。对于每个自振频率,由式(3)可确定一组各节点的振幅值Xi=[Xi1Xi2···Xin]T,它们互相之间应保持固定的比值,但绝对值可任意变化,它们构成一个向量,称为特征向量,也就是工程上结构的振型。
ANSYS提供了7种模态提取方法:BlockLanc-zos(分块法),PCGLanczos(子空间法),Supernode(动态功率法),Reduced(缩减法),Unsymmetric(非对称法),Damped(阻尼法)和QRDamped(QR阻尼法)。文中使用的是第一种方法,即BlockLanczos,它适用于大型对称特征值求解问题,比PCGLanczos具有更快的收敛速度。
2斜齿圆柱齿轮实体建模
ANSYS的实体建模能力相对较差,斜齿轮在
ANSYS中直接建模有一定难度。而CATIA软件在实体建模、曲面造型方面功能强大,且与ANSYS具有良好的数据接口,因此本文是通过CATIA软件对齿轮进行参数化建模,保存为iges或model格式,然后将模型导入到ANSYS软件中。其中,斜齿轮的设计参数:模数m为2mm,齿数z为24,螺旋角为10◦,齿宽为20mm,内孔直径为20mm,键宽为6mm。斜齿圆柱齿轮实体模型如图1
所示。
图1斜齿圆柱齿轮实体模型
Fig.1Theentitymodelofhelicalcylindricalgear
3斜齿圆柱齿轮有预应力的模态分析
斜齿轮往往处于一定转速的工作状态,单纯进
行静态模态分析,不符合实际工况,而需对其进行动态模态分析。斜齿轮旋转时产生的离心应力作为预
应力,对其进行有预应力的模态分析,对齿轮动态设计具有一定的指导作用。有预应力的模态分析除了
首先要进行静力学分析把预应力施加到结构上外,有预应力模态分析的过程与普通模态分析基本一致。ANSYS模态分析包括4个主要步骤:建模,加载和求解,扩展模态,查看结果和后处理。3.1斜齿圆柱齿轮有限元模型及约束
将CATIA中建立的三维模型通过iges或model格式导入有限元分析软件ANSYS中,采用自由划分方法,选取Brick8node185单元进行网格划分,网格划分后的模型如图2所示,共划分节点8358个,单元37837个。齿轮材料选用20CrMnTi,其弹性模量E=200GPa,泊松比µ=0.3,密度ρ=7.8×103kg/m3
。
图2斜齿圆柱齿轮网格划分
Fig.2Themeshingofhelicalcylindricalgear
在典型的模态分析中唯一有效的“载荷”是零位移约束,如果在某个自由度(DOF)处指定了一个非零位移约束,则以零位移约束替代该DOF处的设置。除零位移约束之外的其他载荷,即使施加也会被忽略,因此模态分析时,不需要施加斜齿轮所受的径向力、圆周切向力、轴向推力等其他载荷[10]。要准确预测斜齿轮的动力学模态特性,除了要施加一定的转速,还必须考虑合理的边界条件的影响[11]。为了比较不同边界条件对模态分析的影响,按不同情况设置边界约束条件。
约束类型1约束斜齿轮轮毂内圈节点与键槽
两侧的所有自由度(见图3(a)),这种约束类型也是目前通常所采用的边界条件。
约束类型2考虑斜齿轮工作时键槽一侧受力的影响,因此约束斜齿轮轮毂内圈节点与键槽一侧的所有自由度(见图3(b))。
约束类型3考虑斜齿轮工作时两端面靠轴肩和套约束,与其他零件并非全部接触,因此在斜齿轮
(a) 䊝㭙㏁㾮1(b) 䊝㭙㏁㾮2(c) 䊝㭙㏁㾮
3
图3斜齿轮不同约束形式
Fig.3Thedifferentconstraintformsofhelicalcylindricalgear
轮毂内圈与键槽一侧约束x方向自由度,在斜齿轮两端面孔中心∅20∼∅30接触区间约束z方向自由度(见图3(c))。
3.2斜齿轮有预应力的模态分析加载及求解
(1)静力学分析。对斜齿轮施加角速度ω=200rad/s,将离心应力作为预应力施加到斜齿轮上,此时预应力效果选项必须打开。该命令路径为MainMenu—Solution—AnalysisType—AnalysisOp-tions。
(2)重新进入Solution,进行模态分析。同样,预应力效果选项也必须打开(PSTRES,ON)。该命令路径为MainMenu—Solution—AnalysisType—
AnalysisOptions。
3.3不同边界条件下的斜齿轮模态分析3.3.1固有频率分析
低阶频率对齿轮的影响较大,一般取5∼10阶。在不同边界约束下,斜齿轮前10阶模态的固有频率如表1和图4所示。由图4可见,随着阶数的增加,每一种约束类型的固有频率均逐渐增大;3种约束类型的计算结果有很大的出入,从整体上看,约束程度越强,其固有频率越大;约束类型1与约束类型2边界的约束条件类似,其各阶固有频率相近。
表1不同约束类型的固有频率
Tab.1Thenaturalfrequenciesofdifferentconstraintforms
约束类型123
第1阶31.90630.76611.239
第2阶37.28236.46730.191
第3阶37.30737.05537.171
第4阶40.04438.13039.325
固有频率/kHz第5阶43.60042.13542.650
第6阶43.63943.00443.078
第7阶46.36243.04145.309
第8阶47.12244.95947.987
第9阶57.30956.09155.878
第10阶58.08256.60756.783
䇱㠖㔫/kHz3最为贴近实际,应优先选用。约束类型1、2的前2阶固有频率误差很大,第1阶约为实际值的3倍,其他几阶固有频率也有不小的误差,呈现低阶不符、高阶相似的趋势。3.3.2振型分析
ㅸ㭞
选用约束类型1和约束类型3进行振型比较分析。各提取前5阶振型图进行比较。图5为约束类型1的前5阶振型图,第1阶为圆周振,第2、3阶为扭转振,第4阶为伞振,第5阶为对折振。图6为约束类型3的前5阶振型图,第1阶为圆周振,表现为
图4不同约束类型的固有频率
Fig.4Thenaturalfrequenciesofdifferentconstraintforms
在对斜齿轮进行模态分析时,很明显,约束类型
绕轴线发生摆动,轴向基本无振动;第2阶和第3阶为摆振,表现为分别绕着与轴线垂直的两条轴线转动;第4阶和第5阶为扭转振。比较两种不同边界约
束条件下的振型可知,斜齿轮的振型除了第1阶振型图均为圆周振以外,其他4阶振型图均不相同。
(a) ⭻1ㅸ(b) ⭻2ㅸ(c) ⭻3ㅸ
(d) ⭻4ㅸ(e) ⭻5
ㅸ
图5约束类型1的前5阶振型图
Fig.5Thetop5modeshapesofconstrainttype1
(a) ⭻1ㅸ(b) ⭻2ㅸ(c) ⭻3ㅸ
(d) ⭻4ㅸ(e) ⭻5
ㅸ
图6约束类型3的前5阶振型图
Fig.6Thetop5modeshapesofconstrainttype3
4结论
(1)在CATIA和ANSYS软件环境下,分别建立了斜齿圆柱齿轮的三维几何模型和动力学模态分析有限元模型,进而对一定转速下不同边界约束条件的齿轮进行了模态分析,得到了各阶固有频率和振型。
(2)不同边界约束条件对斜齿轮模态分析影响较大。目前常常采用的边界约束条件下的前10阶固有频率值大部分要高于实际值,尤其是第1、2阶固有频率误差较大,且第1阶固有频率约为实际值的3倍,呈现低阶不符、高阶相似的趋势;前5阶振型除了第1阶振型图均为圆周振以外,其他4阶振型图均发生了变化。
(3)目前通常计算的斜齿轮模态分析与实际情况存在差异,不能满足高精确性要求。为提高斜齿轮模态分析的精度,应考虑转速的影响及边界约束条件的施加。根据实际工作状态确定边界约束条件:在斜齿轮轮毂内圈与键槽一侧约束x方向自由度,在斜齿轮两端面接触区域约束z方向自由度。在此基础上,考虑转速对斜齿轮的离心作用,进行有预应力的模态分析,此方法为齿轮动态分析与设计提供了很好的参考依据。
参考文献:
[1]
MAIANMM,SILVAJMM.Theoreticalandexperimen-talmodalanalysis[M].Hertfordshire:ResearchStudiesPressLtd,1997:110-124.[2][3][4]
慕灿.斜齿圆柱齿轮模态的有限元分析[J].宜春学院学报,2013,35(12):53-55.
于英华,陈棕桢.斜齿轮的参数化建模及有限元模态分析[J].组合机床与自动化加工技术,2007(09):25-27.周建钊,徐肖攀,储伟俊.基于Pro/E与ANSYS的减速箱齿轮模态分析[J].机械制造与自动化,2013,42(5):136-139.[5][6][7][8][9]
王犇,华林.高速旋转状态下汽车弧齿锥齿轮的动力学模态分析[J].汽车工程,2011,33(5):447-451.
叶友东,周哲波.基于ANSYS直齿圆柱齿轮有限元模态分析[J].机械传动,2006,30(5):63-65.
杨伟,马星国,尤小梅.基于ANSYS的齿轮装配体模态分析[J].沈阳理工大学学报,2008,27(4):71-75.胡瑞帆,吴新跃.高速旋转状态下齿轮轴系耦合的模态分析[J].机械研究与应用,2012(4):82-89.
张洪信,王怀敏,孟祥踪,等.ANSYS基础与实例教程[M].北京:机械工业出版社,2013.
[10]高耀东,李强,张玉宝,等.ANSYS机械工程应用精华
60例[M].北京:电子工业出版社,2012.
[11]孙军,桂长林,汪景峰.边界条件处理对曲轴有限元分
析的影响研究[J].汽车工程,2005,27(6):724-726.
EffectsofBoundaryConstraintConditionsonDynamicModalAnalysisof
HelicalGearsBasedonANSYS
LIUSai,
WUFei-ke,YANJian-gang,YANGGuo-ce
(EngineerTrainingCenter,ShanghaiSecondPolytechnicUniversity,Shanghai201209,P.R.China)
Abstract:The3D(Three-Dimensional)geometricmodelandFE(FiniteElement)modelfordynamicmodalanalysisofhelicalgearsareestablishedwithsoftwareCATIAandANSYSrespectively.Thenmodelanalysesonaconstantspeedrevolvinggearindifferentboundaryconstraintconditionsareperformedandthenaturalfrequenciesandthemodeshapesareobtained.Theresultsshowthattheboundaryconstraintconditionshavesignificanteffectsonmodalanalysesofgears.Theboundaryconstraintconditionsadoptedgenerallyaredifferentfromtheactualsituation,whichwillleadtochangesinnaturalfrequenciesandmodeshapes.Toimprovetheaccuracyofmodalanalysisofhelicalgears,theboundaryconstraintsconditionsshouldbedeterminedaccordingtoactualworkingconditions.
Keywords:helicalgear;modalanalyses;boundaryconstraintconditions;naturalfrequency;modeshape