机床工作台的动力学分析.pdf

2008年9月第36卷第9期机床与液压MACHINETOOLHYDRAULICSSEP12008VOL136NO19收稿日期20071205作者简介张会端1972,女,博士研究生,主要研究方向为机械振动学。谭庆昌1957,男,教授,博士生导师,研究方向为机械传动与控制。电话043185094594,EMAILTANQCJLU1EDU1CN。机床工作台的动力学分析张会端,谭庆昌,裴永臣吉林大学机械科学与工程学院,吉林长春130025摘要在考虑丝杠与螺母之间的接触关系以及轴承支承刚度的情况下,利用拉格朗日方程建立了滚珠丝杠传动系统的动力学方程,用振型叠加原理和隆格库塔法进行求解,研究了滚珠丝杠传动系统的瞬态响应,按正交试验计算分析了系统参数对工作台轴向振动的影响,优化组合了系统参数,为滚珠丝杠传动系统参数的选择以及装配调整提供了理论依据,为提高工作台轴向定位精度奠定了基础。关键词机床；动力学；传动系统；参数优化中图分类号TB122文献标识码A文章编号10013881200890684THEDYNAMICANALYSISOFTHEWORKTABLEOFMACHINESZHANGHUIDUAN,TANQINGCHANG,PEIYONGCHENCOLLEGEOFMECHANICALSCIENCEANDENGINEERING,JILINUNIVERSITY,CHANGCHUNJINLIN130025,CHINAABSTRACTINTHECONDITIONOFCONSIDERINGBOTHTHEBEARINGSTIFFNESSANDTHECONTACTDEFORMATIONBETWEENTHEBALLSCREWANDTHENUT,DYNAMICMODELOFTHEBALLSCREWDRIVESYSTEMOFMACHINESWASESTABLISHEDUSINGLAGRANGEEQUATIONTHEMODELWASANALYZEDBYTHEMODESUPERPOSITIONANDRUNGE2KUTTAMETHODTOCALCULATETHETRANSIENTRESPONSEOFTHESYSTEMFROMTHEORTHOGONALEXPERIMENT,EFFECTOFTHESYSTEMPARAMETERSONAXIALVIBRATIONOFTHEMACHINEWORKTABLEWASSTUDIED,THEOPTIMISTICSYSTEMPARAMETERSFORTHEMINI2MUMVIBRATIONWEREOBTAINEDKEYWORDSMACHINE；DYNAMIC；DRIVESYSTEM；PARAMETEROPTIMIZING0引言机床传动系统一般包括伺服电机、滚珠丝杠以及随螺母做轴向往复运动的工作台。随着进给速度和加工精度的提高,工作台的运动精度和传动系统控制成为人们关注和研究的焦点14。研究传动系统的振动是利用控制技术提高工作台运动精度和定位精度的基础,特别是伺服电机越来越广泛的应用,使得传动系统瞬态响应的研究越来越重要5。滚珠丝杠传动系统的振动,可以看作是带有移动质量块的梁的振动。ARASHYAVARI等6用离散单元研究了带有移动质量块的TIMOSHENKO梁的横向振动。S1ZIBDEH等7研究了EULER2BERNOULLI梁、RAYLEIGH梁和TIMOSHENKO梁在移动的随机力作用下的横向振动。P1LEE8研究了简支条件下受移动质量块作用的TIMO2SHENKO梁的横向振动,并与在移动力作用下的横向振动做了比较。U1C1CU和C1C1CHENG9研究了由伺服电机驱动的带有移动质量块轴的扭转和弯曲振动。PAOLOGALLINA10研究了带有移动螺母的丝杠的扭转振动与纵向振动。机床工作时,工作台受力如图1所示垂直于导轨底面的主切削力PC,与吃刀方向相反的背向力PR,沿X方向的进给力PE,工作台移动时产生的摩擦力FF。主切削力与背向力由机床导轨的底面及侧面支承,因此,沿滚珠丝杠轴向受到的外力为工作台与导轨之间的摩擦力FF以及工作台受到的进给力PE,二者组成轴向力PU通过工作台作用到丝杠上。图1传动系统示意图在滚珠丝杠传动系统中,丝杠和螺母之间产生的轴向接触变形和扭转接触变形,是导致丝杠产生纵向振动和扭转振动的主要激励,引起螺母产生轴向位移,从而影响工作台的定位精度,所以丝杠的纵向振动和扭转振动是研究丝杠传动精度首先应考虑的问题1011。笔者考虑了丝杠与螺母之间的接触变形以及轴承的支承作用,利用拉格朗日方程,建立了滚珠丝杠与质量块即工作台系统的动力学方程,利用振型叠加原理和隆格库塔法求解,分析了传动系统的参数对工作台轴向振动的影响,获得了使振幅减小的优化参数组合,为滚珠丝杠传动系统的设计和装配提供了理论依据,为提高工作台定位精度奠定了基础。1滚珠丝杠传动系统的动力学模型图2滚珠丝杠传动系统的动力学模型滚珠丝杠传动系统的动力学模型如图2所示,各参数含义见附录。丝杠以角速度Ω转动,质量块沿X轴的位移ST,包括丝杠转动引起的轴向位移ΛΘD、丝杠的扭转变形引起的质量块的轴向位移ΛΦΛΘD,T忽略二阶小量、丝杠的轴向变形引起的质量块的轴向位移UΛΘD,T,以及质量块在螺母与丝杠接合处的局部振动位移QS,记为STΛΘDΛΦΛΘD,TUΛΘD,TQS1滚珠丝杠传动系统的动能表示为T12MS212∫L0ΡA9UX,T9T2DX12∫L0ΡJΩ9ΦX,T9T2DX2滚珠丝杠传动系统的势能表示为U12∫L0EA9UX,T9X2DX12∫L0GJ9ΦX,T9X2DX12KUQ2S12KΦΦΛΘD,T212KU1U0,T212KU2UL,T212KΦ1Φ0,T212KΦ2ΦL,T23系统耗能函数表示为D12CUQ2S12CΦ9ΦΛΘD,T9T212CU1DU0,TDT212CU2DUL,TDT212CΦ1ΩDΦ0,TDT212CΦ2ΩDΦL,TDT4系统外力为质量块所受摩擦力和进给力的合力PU,所做虚功表示为ΔWUPUΔS5由文献1213,设丝杠的纵向振动、扭转振动位移函数在时间和空间上可分离,并表示为UX,TUXQUT6ΦX,TΨXQΦT7将式6和7代入式2、3、4后,连同式5,根据拉格朗日方程,并利用式1,整理后得系统的动力学方程为M{Q}C{Q}K{Q}{FT}8其中{Q}SQUQΦT,{Q},{Q}分别代表加速度和速度,M、C、K分别为系统的质量矩阵、阻尼矩阵和刚度矩阵,{FT}为力向量。2丝杠的频率方程与振型函数211丝杠的纵向振动由文献1213,设丝杠纵向振型函数为UXCSINΩUAUXDCOSΩUAUX9边界条件为EADUXDXX0KU1U0010EADUXDXXLKU2UL011其中AUEΡ,C、D为常数。将式9代入以上两式,得EACΩUAUKU1D012EACΩUAUCOSΩUAULDΩUAUSINΩUAULKU2CSINΩUAULDCOSΩUAUL013令CKU1AU,由式12、13可得丝杠的纵向振动频率方程为E2A2Ω2UKU1KU2A2USINΩUAULEAΩUAUKU1KU2COSΩUAUL14把CKU1AU和式12代入式9,可得考虑轴承支承刚度的丝杠纵向振型函数UXKU1AUSINΩUAUXEAΩUCOSΩUAUX15为使振型函数唯一,把振型函数正则化,即令Δ2U∫L0ΡAU2XDX116把式15代入式16,得Δ2U1ΡAK2U1A2UE2A2Ω2U2LAU4ΩUE2A2Ω2UK2U1A2USIN2ΩUAULKU1EAA2U21COS2ΩUAUL17则两端轴承支承的丝杠纵向振动正则振型为UXΔUKU1AUSINΩUAUXEAΩUCOSΩUAUX1896第9期张会端等机床工作台的动力学分析212丝杠的扭转振动设丝杠扭转振动的振型函数为ΨXCSINΩΦAFXDCOSΩΦAFX19与纵向振动同理,可得两端轴承支承的丝杠扭转振动的频率方程与正则振型为G2J2Ω2FKF1KF2A2FSINΩFAFLGJΩFAFKF1KF2COSΩFAFL20ΨXΔFKF1AFSINΩFAFXGJΩFCOSΩFAFX21其中AFGΡΔ2F1ΡJK2F1A2FG2J2Ω2F2LAF4ΩFG2J2Ω2FK2F1A2FSIN2ΩFAFLKF1GJA2F21COS2ΩFAFL3模型验证为验证本文传动系统模型的正确性,在方程8中,令S0,ΦX,T0,KU1、KΦ1→∞,KU2KΦ20,M0,KUKΦ0,Ω0,不计阻尼,轴向力PTPSINΩT作用在杆的右端,则模型简化为文图3杆纵向振动的对比献14中一端固定、另一端自由的杆的纵向振动。为与该文献的结果对比,取P200N、Ω2Π进行计算,结果如图3所示。图中曲线吻合,说明本文建立的系统动力学方程8是正确的。4结果与分析把频率及相应的振型函数式18、21代入式8,利用振型叠加原理和隆格库塔法求解方程,分析工作台沿丝杠的运动。根据式1可知,丝杠转动引起的质量块位移S1ΛΘD为线性位移,而系统振动引起的质量块位移S2ΛΦΛΘD,TUΛΘD,TQS为波动位移,是影响定位精度的主要因素,所以下面分析传动系统参数对S2的影响。按正交试验,利用L12211正交表,确定要分析的参数及其取值,见表1。根据表1的数据,计算位移S2,分析各参数对质量块振幅的影响,结果见表2。由表可见,丝杠的导程越大,质量块S2的振幅越小；质量块的质量越小,振幅越小；质量块与丝杠连接的轴向接触刚度越大,扭转刚度越小,质量块振幅越小；丝杠两端的轴承纵向支承刚度越大,扭转刚度越小,质量块的振幅越小；而系统各连接处的阻尼和系统所受的轴向力对质量块的振幅影响不大。表1按正交试验选取的系统参数及其取值取值参数KU1KU2/NΜM1KΦ1KΦ2/NRAD1KU/NΜM1KΦ/NRAD1M/KGH/MMCU/NSM1CΦ/NSRAD1CU1CU2/NSM1CΦ1CΦ2/NSRAD1PU/N140001002000400300815050200304002600020030008005001030010040060800表2各参数对质量块振幅影响的分析结果104M参数1234567891011S2J131222179316221882131416421822198218921902193S2J221763119213631103167113431163100310931083105RJ01460140112601221136313001340102012001180112优值21211211111主次参数4537216118910注S2J1、S2J2为各参数取不同值时,质量块振幅的极差分析的均值；RJ为各参数的极差。根据减小S2的幅值,提高工作台运动精度的原则,通过正交计算分析,确定传动系统参数的优化组合KU1KU26000N/ΜM；KΦ1KΦ2100N/RAD；KU3000N/ΜM；KΦ400N/RAD；M300KG；H10MM；CU150N/RAD；CΦ50N/RAD,CU1CU2200NS/M；CΦ1CΦ230NS/RAD；PU400N。按优化参数计算质量块的振动曲线,如图4所示。质量块在最初的波动之后很快趋于平稳,且振幅较小,说明通过调整传动系统的参数,可以使工作台处于比较平稳的轴向移动状态。图4质量块的优化振动曲线07机床与液压第36卷5结论笔者利用拉格朗日方程建立了滚珠丝杠传动系统的动力学方程,用振型叠加原理和隆格库塔法求解,研究了滚珠丝杠传动系统的动态响应,分析了系统参数对质量块轴向振动的影响,并优化了系统参数组合。为了减小工作台轴向位移的波动,提高定位精度,设计或装配传动系统时,在满足其它条件的前提下,应遵循以下原则增大丝杠导程；减小工作台的质量；增大丝杠两端轴承的纵向支承刚度,减小扭转支承刚度；增大螺母与丝杠连接的纵向刚度,减小扭转刚度。参考文献【1】AMINKAMALZADEH,KAANERKORKMAZACCURATETRACKINGCONTROLLERDESIGNFORHIGH2SPEEDDRIVESJINTERNATIONALJOURNALOFMACHINETOOLSMANUFACTURE,2007,4713931400【2】DOMINICSGUECARRA,AKIRAKYUSOJINDEVELOPMENTOFANEWLAPPINGMETHODFORHIGHPRECISIONBALLSCREW2NDRE2PORTDESIGNANDEXPERIMENTSTUDYOFANAUTOMATICLAPPINGMACHINEWITHINPROCESSTORQUEMONITORINGSYSTEMJJOURNALOFTHEINTERNATIONALSOCIETIESFORPRECISIONENGI2NEERINGANDNANOTECHNOLOGY,2002,26389395【3】R1WHALLEY,M1EBRAHIMI,A1A1ABDUL2AMEERHYBRIDMODELINGOFMACHINETOOLAXISDRIVEJJOURNALOFMA2CHINETOOLANDMANUFACTURE,2005,4515601576【4】MIN2SEOKKIM,SUNG2CHONGCHUNGINTEGRATEDDESIGNMETHODOLOGYOFBALL2SCREWDRIVENSERVOMECHANISMSWITHDISCRETECONTROLLERSPART1MODELINGANDPERFORMANCEANALY2SISJMECHATRONICS,2006,16491502【5】H1DIKENDYNAMICBEHAVIOROFACOUPLEDELASTICSHAFT2ELASTICBEAMSSYSTEMJJOURNALOFSOUNDANDVIBRA2TION,2006,293115【6】ARASHYAVARI,MOSTAFANOURIDISCRETEELEMENTANALYSISOFDYNAMICRESPONSEOFTIMOSHENKOBEAMSUNDERMOVINGMASSJADVANCEINENGINEERINGSOFTWARE,2002,33143153【7】H1S1ZIBDEH,H1S1JUMADYNAMICRESPONSEOFAROTATINGBEAMSUBJECTEDTOARANDOMMOVINGLOADJJOURNALOFSOUNDANDVIBRATION,1999,2235741758【8】H1P1LEEDYNAMICRESPONSEOFATIMOSHENKOBEAMONWINKLERFOUNDATIONSUBJECTTOAMOVINGMASSJAP2PLIEDACOUSTICS,1998,553203205【9】U1C1CU,C1C1CHENGVIBRATIONANALYSISOFAHIGH2SPEEDSPINDLEUNDERTHEACTIONOFAMOVINGMASSJJOURNALOFSOUNDANDVIBRATION,2004,27811311146【10】PAOLOGALLINAVIBRATIONINSCREWJACKMECHANISMSEXP2ERIMENTALRESULTSJJOURNALOFSOUNDANDVIBRATION,2005,28210251041【11】王培功,彭伟,等数控机床进给系统的动态优化设计J机械设计与制造,20058101103【12】J1H1金斯伯格机械与结构振动理论与应用M中国宇航出版社,2005【13】张义民机械振动力学M吉林科学技术出版社,20009【14】师汉民机械振动系统分析测试建模对策下册M华中科技大学出版社,20043【15】任露泉实验优化技术M机械工业出版社,19872附录文中各符号的意义S为质量块的轴向位移；UX,T为丝杠的纵向振动位移；ΦX,T为丝杠的扭转振动位移；KU1、KU2分别为轴承1和轴承2的纵向支承刚度；KΦ1、KΦ2分别为轴承1和轴承2的扭转支承刚度；KU为螺母和丝杠接合部的纵向刚度；KΦ为螺母和丝杠接合部的扭转刚度；CU1、CU2分别为轴承1和轴承2的纵向阻尼；CΦ1、CΦ2分别为轴承1和轴承2的扭转阻尼；CU为螺母和丝杠接合部的纵向阻尼；CΦ为螺母和丝杠接合部的扭转阻尼；M为工作台的质量；PU为工作台所受轴向力；H为丝杠的导程；Ω为丝杠的角速度,Ω50ΠRAD/S；ΘD为丝杠的转角,ΘDΩT；L为丝杠的长度,L550MM；D为丝杠直径,D25MM；E为丝杠的弹性模量,E2107105MPA；G为丝杠的剪切模量,G813104MPA；Ρ为丝杠的密度,Ρ7185103KG/M3；A为丝杠的当量面积,AΠD24；J为丝杠截面的极惯性矩为,JΠD432；Λ为质量块移动速度系数为,ΛH2Π；TS为质量块沿丝杠的运动周期,TS2ΠLHΩ；ΔU、ΔF分别为纵向振型和扭转振型的正则系数。上接第67页【6】石德珂,金志浩材料力学性能M西安西安交通大学出版社,1998【7】赵升吨,柳伟,张立军,化春键棒料热应力预制裂纹的精密下料方法的数值模拟C成都第二届锻压装备与制造技术论坛论文集,20051117第9期张会端等机床工作台的动力学分析