TemperatureFieldAnalysisandExperimentalTestofCNCLathePrecisionSpindleSystemwithoutLoadLixiuZhang1,a,QinghuaShi2,bandYuhouWu3,cShenyangJianzhuUniversity,Shenyang110168,Chinaa851617088@qq.
com,b474891749@qq.
comKeywords:mechanicalspindle;thermalanalysis;temperaturefield;thermalbalancetime.
Abstract:ThearticleisbasedonprecisionCNClathespindlesystemofthe170CP06astheresearchobject.
Thesteadyandtransienttemperaturefieldmodelofthemechanicalspindleisestablishedbyusingthefiniteelementanalysismethod.
Astheheatratefortheloadandtheconvectiveheattransfercoefficientfortheboundaryconditionsforthermalanalysis,thetemperaturefieldofspindlesystemisinvestigatedtogetthesteadystatetemperaturefielddistributionandthermalbalancetimeofthespindlesystem.
Inordertoverifythecorrectnessofthemodel,thethermalexperimentswithoutloadareexecutedtothekeypartsofthetemperatureofthespindlesystematdifferentspeeds.
Comparedwiththesimulationresults,ithasaverygoodagreementwiththereliabilityofthemodel.
Anditprovidesatheoreticalbasisforcontrollingthetemperatureriseofthespindlesystem.
IntroductionWiththerapiddevelopmentofmodernmachinerymanufacturingtechnology,therequirementsofmachine'saccuracyarealsobecominghigherandhigher.
Inthehigh-speedandprecisionmachinetools,theproblemofmachiningerrorcausedbythermaldeformationalsohasbecomeincreasinglyserious.
Alargenumberofexperimentalstudiesshowthethermalerroristhelargesterrorsourceofmachinetoolsaccountingfor40%to70%inthetotalerrorofthemachine[1].
Thespindleisoneoftheimportantcomponentsofthemachine,anditsperformanceaffectstheoveralllevelofthemachinedirectly.
Whenthemechanicalspindleisworking,theinternalheatsourceofthespindlecomesformtheheatgeneratedbythefrictionofhigh-speedrotation'sbearing.
Whenthetemperaturehasrisen,therelativespatialpositionandsizeofthespindleandotherpartsofthemachinewillbedifferentfromtheprevioustemperature.
Itwillformadifferenttemperaturefieldandproducedifferentdegreesofthermaldeformation[2].
Therefore,thethermalcharacteristicsofresearchandanalysisonthespindlesystemiscrucialtoensuretheaccuracyofthemachine.
Itisoneofthekeytechnologiestobeconsideredinhighprecisionmachinetools.
Thethermalcharacteristicsanalysisofthespindlesystemistostudythedistributionofthetemperaturefield.
Atpresent,therearemanyscholarswhostudythethermalcharacteristicsofmachinetoolspindlesystemathomeandabroad.
Forexample,JinKyungChoiofSouthKoreastudiedthermalcharacteristicsforthespindleandbearingsystemwiththefiniteelementmethod[3].
M.
H.
AttiafromCanadaanalyzedthewholetemperaturefieldofthemachineusingthefiniteelementmethod[4].
CeGuoandQinghongSunfromSoutheastUniversityhavealreadyestablishedafiniteelementmodelofthehigh-speedprecisionlathespindlecomponents'temperaturefield[5].
YouweiHongestablishedthefiniteelementmodelofthegantrymachiningcenterandsimulatedthetemperaturefielddistributionandvariation[6].
ThearticleisbasedonCNClatheprecisionspindlesystemofthe170CP06astheresearchobject.
Thetemperaturefieldmodelofthemechanicalspindleisestablishedanditisconductedthesteady-stateandtransientthermalanalysisofthespindlesystem.
Thetemperaturefieldofspindlesystemisinvestigatedtogetthetemperaturefielddistributionandthermalbalancetimeofthespindlesystem.
Thethermalexperimentswithoutloadarecarriedoutinordertoverifythereliabilityofthemodel.
CalculationofthethermalparametersTheheatsourceisthefoundationofthespindletemperaturefieldresearch.
Theheatsourceofthespindlesystemincludescuttingheatandbearingfrictionheat.
Thecuttingheatinthemachiningprocessistakenawaybycoolantliquidandswarf.
Therefore,theheatsourceofthespindlemaybemainlythebearingfrictionalheat.
CalculationofthecalorificvalueTherollingbearing'sheatismainlygeneratedbythebearingfrictionaltorque.
Theformulais41.
04710fHnM=*(1)Intheformula,Hfisthecalorificvalueofthebearing;nisthespindlespeed;M0isthebearingfrictionaltorque.
Thebearingfrictionaltorqueisthesumofthebearingrollingfriction,slidingfrictionandlubricantfrictionwhichhavegeneratedtoblockthebearingrunningtorque.
Palmgren[7]deducedtheempiricalformulawhichcalculatesthebearingfrictiontorque.
Thebearingfrictionaltorqueiscalculatedas01MMM=+(2)Intheformula,M0isrelatedtotheviscosityofthelubricant.
M1isrelatedtotheloadofthebearingregardlessofspeed.
ItonlyconsidersM0withoutload.
Palmgrenconfirmedtheexpressionusingempiricalmethods:72/3300=10()mMfnDν2000nν≥(3)730016010mMfD=*2000nν<(4)Intheformula,νisthekinematicviscosityofthelubricantinthebearingoperatingtemperature,mm2/s;0fisthecoefficientrelatedtothebearingtypesandlubrication,mDistheaveragediameterofthebearing,mm,0.
5()mDDd=+.
Whenthespindlesystemisanalyzed,thethermalloadisloadedbythewayofbearingtheheatrate.
Theformulais[9]3(/)fHqWmV=(5)Intheformula,Visthevolumeoftheheatsource;22)2/(bmDDVπ=,andDbisthediameteroftherollingelements.
Calculatedbytheequation(1)to(5),whenthespindlespeedis2000r/min,wecangettheheatratesoftheforeandrearbearingtobe342704.
4W/m3and348161.
76W/m3.
CalculationoftheheattransfercoefficientAccordingtotheheattransfertheory[10],wecanknowtheheattransferwayaretheheatconduction,theconvectionandthethermalradiation,thelossoftheradiationheatislessforthespindlesystem.
Whenweareanalyzingthetemperaturefieldofthespindlesystemonlytoconsiderthethermalconductionandtheconvectionheattransfer.
Thethermalconductivitydependsonthethermalconductivityofthecomponentsinthespindlesystem.
Theconvectiveheattransferistheconvectionofthespindlesurfaceandtheair.
Theconvectiveheattransferisformedbydrivingtherotationofthespindlearoundtheairbetweenthecaseandthespindle.
ThecoefficientoftheheattransferdependsonthecriterionofNusselt[12].
/cchNlλ=(6)Intheformula,chisthecoefficientoftheheattransfer;λisthethermalconductivityoftheair;NisthecoefficientofNusselt;clisthefeaturesize.
TheNusseltformulaoftheforcedconvectionis23130.
133RePrN=5Re4.
3100.
7Pr670)<*<<(,(7)AccordingtoReynoldsCriterion,wecanknowRe=/cculν(8)Intheformula,ReistheReynoldsnumber;PristhePrandtlnumberofthefluid;cuistheaveragevelocityoftheair;νisthecoefficientoftheair'smotionviscosity.
Whenthespindlerotatesbyhighspeed,theaveragevelocityoftheairis[13]=/60cudnπ(9)1122=nndldldldl++(10)12nllll11)Intheformula,distheaveragediameterofthespindle;nisthespeedofthespindle.
Calculatedbytheequation(6)to(11),whenthespindlespeedis2000r/min,wegeteachheattransfercoefficientswhichare26.
62/()Wmk,35.
632/()Wmk,38.
32/()Wmk,32.
52/()Wmk,35.
22/()Wmk.
FiniteelementmodelofthespindlesystemTheprecisionspindlesystemofthe170CP06ismainlycomposedofthespindle,case,beltpulley,bearingcoverandtheangularcontactballbearings,etc.
Themotordrivesbyabeltpulleymountedontheshafttoprovideadrivingforceforthespindle.
Themainparametersofthespindleis:ItsmaterialisalloySteel.
Thetotallengthis497mm.
Themaximumdiameterofthespindleismm100φ.
Theratedspeedofthespindlesystemis4500r/min.
Themaximumpoweris35kWandthemaximumtorqueis180N.
m.
Fig.
1isthestructurediagramforspindlesystem.
1—Bignut2—Beltpulley3—Leftcover4—Case5—Rightcover6—Spindle7—Forebearing8—RearbearingFig.
1StructurediagramforspindlesystemThemodelofthespindlesystemisestablishedbyusingthethree-dimensionalmappingsoftware—SoildWorks.
Anditsimplifiesappropriatelythestructurethatinfluencesanalysisresultsless.
Accordingtothespindleoftheboundaryconditionsandheattransfercharacteristics,Simplifiedsolidmodelisintroducedtothefiniteelementanalysissoftwaretobemeshingofthegrid[12].
Fig.
2isasimplifiedandmeshingmodelofthespindlesystem.
Theentiremodelhas128,143unitsand232,855nodes.
Fig.
2MeshingofthespindlesystemSimulationanalysisofthespindlesystem'stemperaturefieldThethermalanalysisisestablishedinthesoftwareoffiniteelementanalysis.
Thematerialofthespindleandbearingaredefinedasalloysteel.
Theambienttemperatureis26C°.
Accordingtotheboundaryconditionsofthespindleandthecharacteristicsoftheheattransfer,theheatgenerationrateoftheforeandrearbearingsis342704.
4W/m3and348161.
76W/m3whichareappliedtotheforeandrearbearings.
Thecoefficientsoftheheattransferconvectionareaddedtotherespectivemembersurfacespindlesystemasboundaryconditions.
Itisconductedthesteady-stateandtransientthermalanalysisofthespindlesystem.
TheanalysisofthesteadytemperaturefieldWhenthespindlespeedis2000r/min,wecangetthesimulationofthesteady-statetemperaturefieldresultinginthesteady-statetemperaturefieldcontourofthespindlesystem.
AsisshowninFig.
3.
Fig.
3Steady-statetemperaturedistributionofspindlesystemFig.
3showsthatthemaximumtemperatureofthespindlesystemappearsattheinnerringoftherearbearing.
Itstemperatureis38.
5C°.
Themaximumtemperatureriseis12.
5C°.
Theheatoftherearbearingislargerthantheforebearing.
Anditiscausedbythepoorcoolingconditionoftheinnerring.
TheanalysisofthetransienttemperaturefieldWhenthespindlesystemistheanalysisofthetransienttemperaturefieldtogetthebearingtemperaturecurveandthethermalequilibriumtimebysettingtheoperatingtimeof65min.
Fig.
4arethecontoursofthetemperaturefieldatdifferenttimepoints.
1600s2400s3200s3600sFig.
4ContoursofthetemperaturefieldatdifferenttimepointsThespindlesystemhasreachedthethermalequilibriumafterrunning60min.
Thehighesttemperaturesoftheforebearingandrearbearingare34.
19C°and35.
9C°.
ExperimentaltestTheautomatictestsystemofmechanicalspindlecharacteristicisusedforthetemperaturerisetestwhichisinthestatusoftheofmechanicalspindletorqueoutput,theloadandwithoutloadateachspeed.
Thetestsystemconsistsofthemachinebase,torqueandspeedsensors,industrialcomputer,ervomotor,multi-channeltemperaturecontrolinstrumentandsoon.
Thetestprincipleisthatthespindleconnectsthearresterdetentviacouplings,torqueandspeedsensorandthemotordrivesbyabeltdrivetomakethespindlerotatebychangingthebrakeexcitationcurrentmethodtotheloadcontrolforthespindlesystem.
Fig.
5istheoperatingprincipleoftheexperimentplatform.
Fig.
6istheexperimentdeviceofmechanicalspindlesystem.
Fig.
5OperatingprincipleoftheexperimentplatformFig.
6TheexperimentdeviceofmechanicalspindlesystemThetemperaturesensorsaredisposedwithintheouterringoftheforeandrearbearings,therearendsurfaceoftherearbearingandandoutersurfaceofthecase.
Fig.
7isthedistributionoftestpoints.
Atthespeedof2000r/minand4500r/min,thespindlesystemistestingtemperaturerisethewithoutload.
Whenthespindlesystemhasreachedtheequilibrium,wecanmeasurethetemperatureofeachofthekeyparts.
Wewilltaketestpoint2andtestpoint8forexampleandcanseethechangesinthebearingtemperaturewithtime.
Fig.
8isthetemperaturerisecurvefortestpoint2and8atthespeedof2000r/minand4500r/min.
123456789111210Fig.
7Thedistributionoftestpoints05101520253035404550556065202224262830323436T/°Ct/minSimulationvalueExperimentalvalue05101520253035404550556065202224262830323436T/°Ct/minSimulationvalueExperimentalvaluea.
Thetemperatureoftestpoint2(Speedof2000r/min)b.
Thetemperatureoftestpoint8(Speedof2000r/min)051015202530354045505560652022242628303234363840424446T/℃t/minSimulationvalueExperimentalvalue0510152025303540455055606520222426283032343638404244T/℃t/minSimulationvalueExperimentalvaluec.
Thetemperatureoftestpoint2(Speedof4500r/min)d.
Thetemperatureoftestpoint8(Speedof4500r/min)Fig.
8Temperaturerisecurvefortestpoint2and8(Speedof2000r/minand4500r/min)AscanbeseenfromFig.
8,wecanseethattherearesomeerrorscomparedwithexperimentaldataandsimulationdata.
Butthereisthesameupwardtrendbasicallythatexplainsthemodelisreliable.
Whenthespindlespeedis2000r/minandthesystemhasbeenrunningfor60mins,thesystemhasreachedsteadystatetemperature.
Andthehighesttemperaturesoftestpoint2and8are34.
19C°and35.
9C°.
Therearemoreobviouschangesofthetemperaturerisewhenthespindlesystemisrunningatthebeginning.
Thetemperaturerisechangessmallerafter55minandreachessteadystateabout57.
5min.
Whenthespindlespeedis4500r/min,thespindlesystemrisesfasterbefore52.
5minandreachesequilibriumabout55min.
Next,wewilltestthetemperaturerisewhenthespindlespeedare1000r/min,1500r/min,2500r/min,3000r/min,3500r/minand4000r/min.
Asthekeypartsofthetestpoints2and8forexample,WecanobtainthemaximumtemperatureofeachtestpointunderdifferentspeedsasshowninFig.
9.
05001000150020002500300035004000450026283032343638404244T/℃Speed(r/min)Testpoint2Testpoint8Fig.
9SpindletemperatureatdifferentspeedsAscanbeseenfromFig.
9,asthespindlespeedincreases,thetemperatureriseofthebearingsismoreandmorehigher.
Andthetestpoint2isslightlyhigherthanthetestpoint8.
Whenthespindlespeedchanges1000to3000r/min,thebearingtemperaturechangesmoreobviously.
Thechangeofthetemperatureisflatrelativelywhenthespindlespeedchanges3000to4500r/min.
ConclusionInthispaper,itisconductedthethermalanalysisofthespindlesystemcombiningthemethodofthefiniteelementanalysiswithexperimentalverificationandresultinginthefollowingconclusions:(1)Theaboveanalysisshowsthatthesimulationmodelofthemachinespindleisreliable.
Theresultoferrorsislessthan3%comparedthesimulationdatawiththeexperimentaldata.
(2)Whenthespindlespeedischanging1000to3000r/min,thetemperatureofbearingshavechangedmoreobviously.
Thechangeofthetemperatureisflatrelativelywhenspindlespeedischanging3000to4500r/min.
Itindicatesthatthespeedimpactsthetemperaturegreaterunderlowspeeds.
Andforthehighspeed,thetemperatureisimpacttothespeedrelativelysmall.
Forthespindle,whichistheratedspeedof4500r/min,thetemperaturechangesrapidlyunderitstwo-thirds'ratedspeed.
Whilethespeedismorethantwo-thirds'ratedspeed,thetemperaturechangesslowly.
(3)Thespindlesystemhasreachedthedifferentthermalequilibriumtimeatdifferentspindlespeeds.
Thehigherthespeedis,theshorterthespindlehasreachedthermalequilibriumtime.
(4)Theprecisionspindleof170CP06hasreachedthethermalequilibriumatratedspeedfor55min.
Thehighesttemperatureis45.
2C°thatoccurstothebearbearing.
References[1]D.
A.
Krulewich,Temperatureintegrationmodelandmeasurementpointselectionforthermallyinducedmachinetoolerrors,Mechantronics,1998,8:395~412.
[2]ZhenZhu,PenghaoHu,YiTao.
Theanalysisofthree-dimensionalspindletemperaturefieldandthermaldeformation[J].
Tooltechnology,2008(4):66-68.
[3]Choijin-Kyung,LeeDai-Gil.
Thermal-characteristicsofthespindlebearingsystemwithagearlocatedonthebearing[J].
[4]M.
H.
Attia,L.
Kops,ComputersimulationofnonlinearthermoPlastiebehaviorofajointinmaehinetoolstructureanditseffeetonthermaldeformation,TransaetionsoftheASME,JournalofEnginneringforIndustry,101(1979)355-361.
[5]CeGuo,QinghongSun.
Thethermal-characteristics'analysisofthehigh-speedandhigh-precisionCNClathespindlesystemandthermaldeformation.
JournalofSoutheastUniversity,2005,35(2).
[6]WenPeng,YouweiHong.
Thethermal-characteristicsnumericalanalysisoffive-axisgantrymachiningcenter'sspindlesystem[J].
ModularMachineTool&AutomaticManufacturingTechnique.
2005(08).
[7]XingjuanYing,HaolinLi.
ThefiniteelementanalysisofCNCmachinetoolspindle'ssystemthermalcharacteristics[J]tooltechnology,2010,4(1):38-40.
[8]DatongQin.
Bearing.
Beijing:ChemicalIndustryPress.
2013.
3(5).
[9]XinshengLiBobaoZhao,Yaping.
Thethermalsimulationanalysisofhigh-speedspindle[J]tooltechnology,2012(4):64-32.
[10]F.
P.
Incropera,D.
P.
DeWitt,T.
L.
Bergman.
FundamentalsofHeatandMassTransfer[M].
2012.
[11]ChongzhiGuo,LeXiao.
Asequencecouplingmethodfornumericalsimulationoftemperature[J].
ChemicalIndustryandEngineeringProgress.
2010(09).
[12]ShimingYang,WenquanTao,HeatTransfer[M]Beijing:HigherEducationPress,1998.
[13]XiangqianKong.
Thefiniteelementmethodinheattransfer.
Beijing:SciencePress,1986.
[14]CanZhou,YuhuiHe,GuilingDeng.
Thetransienttemperaturefieldanalysisandtestingofhigh-speedspindle[J].
ComputerSimulation,2012,29(2):372-377.
官方网站:点击访问星梦云活动官网活动方案:机房CPU内存硬盘带宽IP防护流量原价活动价开通方式成都电信优化线路4vCPU4G40G+50G10Mbps1个100G不限流量210元/月 99元/月点击自助购买成都电信优化线路8vCPU8G40G+100G15Mbps1个100G不限流量370元/月 160元/月点击自助购买成都电信优化线路16vCPU16G40G+100G20Mb...
bgpto怎么样?bgp.to日本机房、新加坡机房的独立服务器在搞特价促销,日本独立服务器低至6.5折优惠,新加坡独立服务器低至7.5折优惠,所有优惠都是循环的,终身不涨价。服务器不限制流量,支持升级带宽,免费支持Linux和Windows server中文版(还包括Windows 10). 特色:自动部署,无需人工干预,用户可以在后台自己重装系统、重启、关机等操作!bgpto主打日本(东京、大阪...
RAKsmart发布了新年钜惠活动,即日起到2月28日,商家每天推出限量服务器秒杀,美国服务器每月30美元起,新上了韩国服务器、GPU服务器、香港/日本/美国常规+站群服务器、1-10Gbps不限流量大带宽服务器等大量库存;VPS主机全场提供7折优惠码,同时针对部分特惠套餐无码直购每月仅1.99美元,支持使用PayPal或者支付宝等方式付款,有中英文网页及客服支持。爆款秒杀10台/天可选精品网/大...
qq空间电脑版特别关心为你推荐
文章文章toupianyuanzhu看拼音写词语http404未找到为什么网站上传,打开看不到,显示HTTP 404 - 未找到文件2019支付宝五福支付宝五福是哪五福?asp.net什么是asp.net美要求解锁iPhone苹果美版有锁和无锁的区别购物车在超市、商场中为什么需要使用购物车呢?传奇域名自己的传奇服务器怎么建设?95188是什么电话95188是什么号码我刚收到短信是什么支付宝的验证码300051三五互联170号段和三五互联什么关系
根域名服务器 香港服务器租用99idc 天津服务器租赁 cpanel paypal认证 免费ftp空间 java空间 台湾谷歌网址 免费网页空间 512mb 万网主机管理 iki lamp的音标 免费个人网页 97rb 网络速度 umax 小夜博客 mteam 香港博客 更多