第１２期 机械设计与制造 ２０１１年１２月 Ｍａｃｈｉｎｅｒｙ Ｄｅｓｉｇｎ＆Ｍａｎｕｆａｃｔｕｒｅ ７３ 文章编号：１００１—３９９７（２０１１）１２—００７３—０２ 基于ＡＮＳＹＳ的主轴热变形建模与分析术 孙小帅 董辛昊 刘爱国 尹献德ｚ （’郑州大学机械工程学院，郑州４５０００１）（ 河南省特种设备安全监测研究院，郑州４５０００４） Ｍｏｄｅｌｉｎｇ ａｎｄ ａｎａｌｙｓｉｓ ｆｏｒ ｔｈｅｒｍａＩ ｄｅｆｏｒｍａｔｉｏｎ ｏｆ ｓｐｉｎｄｌｅ ｂａｓｅｄ ｏｎ ＡＮＳＹＳ ＳＵＮ Ｘｉａｏ－ｓｈｕａＰ，ＤＯＮＧ Ｘｉｎ－ｍｉｎｌ，ＬＩＵ Ａｉ－ｇｕｏ２，ＹＩＮ Ｘｉａｎ－ｄｅ２ （。Ｓｃｈｏｏｌ ｏｆ Ｍｅｃｈａｎｉｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｚｈｅｎｇ Ｚｈｏｕ Ｕｎｉｖｅｒｓｉｔｙ，Ｚｈｅｎｇｚｈｏｕ ４５０００１，Ｃｈｉｎａ） （ Ｈｅ’ｎａｎ Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｓｐｅｃｉａｌ Ｅｑｕｉｐｍｅｎｔ Ｓａｆｅｔｙ Ｉｎｓｐｅｃｔｉｏｎ ａｎｄ Ｔｅｓｔ，Ｚｈｅｎｇｚｈｏｕ ４５０００４，Ｃｈｉｎａ） 【摘要】以ＣＸ８０７５立式车铣复合加工中心主轴为研究对象，确立了热边界条件，计算了主轴的 发热量，建立了主轴的三维数字化模型，利用有限元法研究了稳态温度场分布，得出了主轴的热变形， 分析了主轴的热变形原因，奠定了主轴热变形加工误差补偿的理论基础。 关键词：热变形；主轴；稳态温度场 【Ａｂｓｔｒａｃｔ】Ｔａｋｉｎｇ ｔｈｅ ｓｐｉｎｄｌｅ ｏｆＣＸ８０７５ ｖｅｒｔｉｃａｌ ｍａｃｈｉｎｉｎｇ ｃｅｎｔｅｒ ｅｘａｍｐｌｅ，ｔｈｅ ｔｈｅｒｍａｌ ｂｏｕｎｄａｒｙ ｃｏｎｄｉｔｉｏｎｓ ｗｅｒｅ ｄｅｔｅｒｍｉｎｅｄ，ｔｈｅ ｈｅａｔｉｎｇ ｖａｌｕｅ ｗａｓ ｃａｌｃｕｌｔａｅｄ ａｎｄ ３Ｄ ｄｉｇｉｔａｌ ｍｏｄｅｌ ｏｆ ｔｈｅ ｓｐｉｄｎｌｅ ＷＣＩ８ ｅｓ－ ｔａｂｌｉｓｈｅｄ ｗｅｌ１．Ｔｈｅｎ ｔｈｅ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆｔｅｍｐｅｒｔａｕｒｅ ｉｆｅｌｄ Ｗａｓ ｓｔｕｄｉｅｄ ｂｙ ＡＮＳＹＳ ａｎｄ ｔｈｅｒｍａｌ ｄｅｆｏｒｍａｔｉｏｎ ｆｏｓｐｉｎｄｌｅ ＷａＳ ｏｂｔａｉｎｅｄ．Ａｆｔｅｒｗａｒｄｓ ｔｈｅ ｒｅａｓｏｎｓ ｏｆｔｈｅｒｍａｌ ｄｅｆｏｒｍａｔｉｏｎ ｗｅｒｅ ａｎａｌｙｚｅｄ，ｗｈｉｃｈ ｌａｉｄ ａ ｔｈｅｏｒｅｔｉ－ ｃｏｇ ｂａｓｉｓｆｏｒ ｅｒｒｏｒ ｃｏｍｐｅｎｓａｔｉｏｎ ｆｏｔｈｅｒｍａｌ ｄｅｆｏｒｍａｔｉｏｎ ｆｏｔｈｅ ｓｐｉｎｄｌｅ． Ｋｅｙ ｗｏｒｄｓ：Ｔｈｅｒｍａｌ ｄｅｆｏｒｍａｔｉｏｎ；Ｓｐｉｎｄｌｅ；Ｓｔｅａｄｙ－ｓｔａｔｅ ｔｅｍｐｅｒａｔｕｒｅ 中图分类号：ＴＨ１６文献标识码：Ａ １引言 直接耦合法和间接耦合法，直接耦合法的耦合单元需要具有热和 高速高精度加工已经成为现代工业化生产的主流，各种新 结构自由度，采用直接耦合法仅通过一次求解就能得出耦合场的 型机床不断涌现，但由于机床结构的复杂性，多种因素影响机床 分析结果，比如温度分布以及机构变形；间接耦合法又称顺序耦 的加工精度，其中热变形引起的误差是影响加工精度的重要因素 合法，是将第一次分析的结果作为第二次分析的载荷来实现两种 之一。通常情况下机床精密度越高，热误差的影响越大。“ＣＸ系列 场的耦合，比如此处的热一结构耦合分析可以将热分析得到的节 立式车铣复合加工中心”属于国家科技重大专项“高速、精密、复 点温度作为载荷施加在后续的结构分析中，实现热一结构耦合分 合数控金切机床”项目课题之一，其中子课题“ＣＸ系列立式车铣 析。分析中采用直接耦合法对电主轴进行热一结构耦合分析。 复合加工中心高速条件下机床热变形研究”是针对机床热变形和 在ＡＮＳＹＳ中采用ｓｏｌｉｄ２２７单元对电主轴划分网格，在定义 热变形补偿展开的基础性研究。 ｓｏｌｉｄ２２７单元的属性时选取Ｔｅｍｐ， ， ， 四个自由度，使该单 元可以用于直接的热一结构耦合分析。在主轴的有限元模型中简 化两端的轴承，在轴承安装部位设置结构自由度和热载荷以等效 轴承，考虑到轴承安装处是主要发热源，运动中会产生较多的热 量因此可以把此部位的网格划分的细些，远离轴承处的部位网格 可以适当稀疏一些。该有限元模型，如图１所示。 图１加工中心整机３Ｄ模型 ２主轴的有限元模型建立 Ｃｘ８０７５车铣复合加工中心铣削电主轴的最高转速≥ １２０００ｒ／ｍｉｎ；车削主轴最高转速：８００ｒ／ｍｉｎ（￣ｂ８００ｍｍ），５００ｒ／ｍｉｎ 图２ ＣＸ８０７５主轴系统的有限元模型 （￣ｂｌ０００ｍｍ）。该主轴系统主要包括主轴、轴承和轴承端盖等主要 部件。首先通过ＳｏｌｉｄＷｏｒｋｓ２００８三维实体软件建立模型，然后由 ３主轴的温度场分析 ＳｏｌｉｄＷｏｒｋｓ输出转换文件，导入ＡＮＳＹＳ分析软件。此处分析主轴 ３．１主轴发热量的确定 的热变形，为热一结构耦合场分析。热一结构耦合场分析可以采用 考虑在高转速下承受轴向载荷因素，主轴轴承采用压力角 ★来稿日期：２０１１—０２—１０★基金项目：国家科技重大专项—ｃｘ系列立式车铣复合加工中一Ｌ＂（２００９ＺＸ（Ｍ００１－０３３） ７４ 孙小帅等：基于ＡＮＳＹＳ的主轴热变形建模与分析 表１ Ｒ 与Ｃ、ｎ之间的关系 第１２期 为２５。的角接触球轴承，根据理式，轴承的发热量计算公式 为：Ｑ＝Ｉ．０４７ｘ１０－％Ｍ （１） 式中：ｒ厂—轴承转速； ｒ－轴承摩擦力矩。摩擦力矩Ｍ＝Ｍ，＋％，其 中：肘 塞度项，它与润滑剂的粘度和用量，轴承转速有关， 由Ｐａｌｍｇｒｅｎ提出的经验算法公式： Ｍ １０－７（ ３ ２．ｓ－１．ｒ．ｉｎ一 ｏ＝１０ ｆｏ ｕｎ ｄ ） ｍ ｖｎ￣＞２０００ｘｌ０ ｍ‘ｓ‘ｒ‘ｍｉｎ ｖｎ＜２０００ｘｌ０ ｍ。ｓ ａｉｒｎ 电主轴的工作情况设为环境温度２５。，转速为１２０００转，分， 此时主轴与周围空气的对流换热系数（ＣＯＮＶ）为２５０—２６０Ｗ／ ｇｏ＝１６０ｘ１０ 式中：ｄ 承中径（ｍｍ） 一与轴承类型和润滑方式有关的经 （Ｉｎ ｏｃ）。主轴材料４０Ｃｒ的导热系数选取为９１．３ Ｗ／（ｍ・Ｋ）。 将上述边界条件施加于主轴之后的有限元模型如图３所 验常数； —轴承转速（ｒ／ｍｉｎ）； 一在工作温度下润滑剂的运 动粘度（ｍ ）； ．一负荷项： ・ＪＤＩ。 （２） 式中： —确定轴承摩擦力矩的计算负荷；ｄ 承中径 一与 轴承类型和所受的负荷有关的系数： ／ 、０” ＝０・００１‘（ｅＯ／Ｃｏ ｊ （３） Ｐ０ 承的当量静载荷， 取式中的最大值ｃ 式中： —径向载荷系数； —径向载荷；Ｙｏ—轴向载荷系数； 一 轴向载荷；ｃ， —轴承的额定静载荷。当轴承在高速轻载条 件下， 占主要部分；当轴承在低速重载条件下，Ｍ．占主要 部分。 在电主轴热分析时将轴承的发热载荷加于轴承安装位置， 通过计算，前后轴承安装位置的生热率为：６ｘ１０６Ｗ／ｍ 。 ３．２边界条件的确定 热量传递的方式有对流、导热和辐射三种。由于加工中心主 轴系统的温升较小，辐射散失的热量很少，因此只考虑前两种。 热对流可分为两类：自然对流和强制对流。热对流的热量可用牛 顿冷却公式计算，即： Ｑ－－．４ａ（ ） （４） 式中： 位时间通过边界截面的热流量；Ａ—换热面积； 对 流换热系数； 一导热体壁面温度； 界流体温度。 在自然对流条件下，由经验公式求得： Ｇ，ｇ＝￣ｇｄ ｔ／ｖ （５） Ｎｕ＝ｃ（ｃｒＰ￣） （６） Ａ／ｄ （７） 式中：Ｇ —格拉晓夫准数；Ｐ，—普朗特准数；卢—体积膨胀系数； ｄ—轴外径； —流体的运动黏度；ｃ、，ｒ常数，与流动性质、 面朝向有关； —对流换热系数； ．一自然对流传热准数。 在强迫对流条件下，根据努赛尔准则： ７、７『。： 。ｌ｛ｃ．Ｒ：．Ｐｒ ｆ １Ｊ （８）８’ 传热系数： ＡＩ＿．ｊ：２（９） ‘ 式中： 一努赛尔数； 一雷诺数； 一周径；Ａ 一空气导热系数； —普朗特数（来流温度）； —普朗特数（壁温）。其中，雷诺 数Ｒ 与常数Ｃ、ｎ之间的关系，如表ｌ所示。 示，考虑主轴的支撑类型，在一个轴承位置处固定其 ， ， 自 由度，另一个轴承处固定 和 自由度。 图３加载约束的主轴有限元模型 ４主轴温度场及热变形的分析 此处对电主轴进行热稳态分析（ｓｔｅａｄｙ—ｓｔａｔｅ），对图３模型 求解后得到电主轴的温度场分布模型，如图４所示。 图４主轴稳态温度场分布图 由图可知：（１）主轴表面的温度分布很不均匀，最高温度超 过９０。，最低温度和室温持平，这样就造成了主轴热变形的分布 不均。（２）电主轴在前后轴承处温升比较高，这是因为轴承本身发 热量比较大，且都有承载。（３）前轴承内圈温升最高，这是由于前 轴径轴承靠近主轴电机，而且前轴承支撑的载荷比后支撑大，发 热量也较大，而且内圈的散热条件相对外圈不好所致。 主轴在达到热平衡后，根据温度场计算结果进行热变形分 析，得出主轴热变形，如图５所示。 图５主轴热变形图 由图５可以看出主轴在达到热平衡时，主轴前端沿轴径方 向伸长，发生巨大的热变形，发生热变形的原因主要因为靠近电 机以及前轴承重载运动产生的热量没有效导出所致。 第１２期 ２０１１年１２月 文章编号：１００１—３９９７（２０１１）１２—００７５—０３ 机械设计与制造 Ｍａｃｈｉｎｅｒｙ Ｄｅｓｉｇｎ＆Ｍａｎｕｆａｃｔｕｒｅ ７５ 单元测试模式的ＡｕｔｏＣＡＤ三维图形自动判别技术研究术 （河北科技大学机械电子工程学院，石家庄０５００１８） Ｓｔｕｄｙ ｏｎ ａｕｔｏ ｄｉｓｃｒｉｍｉｎａｔｉｏｎ ｔｅｃｈｎｏｌｏｇｙ ｏｆ ＡｕｔｏＣＡＤ ｔｈｒｅｅ—ｄｉｍｅｎｓｉｏｎａＩ ｇｒａｐｈｉｃｓ ｂａｓｅｄ ｏｎ ｕｎｉｔ ｔｅｓｔ ｍｏｄｅ ＸＵ Ｊｉｎ—ｘｉｎ，ＹＡＮＧ Ｓｏｎｇ—ｌｉｎ，ＺＨＡＮＧ Ｇａｎｇ （Ｍｅｃｈａｎｉｃａｌ ａｎｄ Ｅｌｅｃｔｒｏｎｉｃ Ｅｎｇｉｎｅｅｒｉｎｇ Ｃｏｌｌｅｇｅ ｏｆ Ｈｅｂｅｉ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｓｃｉｅｎｃｅ ａｎｄ Ｔｅｃｈｎｏｌｏｇｙ，Ｓｈｉｊｉａｚｈｕａｎｇ ０５００ １ ８，Ｃｈｉｎａ） ； 【摘要】通过对单元测试模式的ＡｕｔｏＣＡＤ三维图形自动判别技术研究，利用Ｖｉｓｕａｌ Ｂａｓｉｃ、Ｖｉｓｕ一： ｝ａｌｌｉｓｐ编程语言以及图形格式转换、图形数据信息遍历搜索方法，解决了三维ＤＷＧ图形批量自动打开｛ ｉ并完成ＤＷＧ至ＳＡＴ数据文件格式的自动转换，并对ＳＡＴ数据文件进行遍历搜索提取三维实体图形数ｉ ；据信息自动存储到Ａｃｃｅｓｓ数据库中，结合特征数据比对法实现ＣＡＤ三维图形的自动评判、结果显示、｛ ；数据输出等功能，为三维ＣＡＤ测试实现网络化、自动化、无纸化奠定重要的基础。 ｛ ｌ 关键词：测试；ＡｕｔｏＣＡＤ；三维；自动判别；ＳＡＴ；ＡＣＩＳ ２ ｉ 【Ａｂｓｔｒａｃｔ】Ｔｈｒｏｕｇｈ ｓｔｕｄｙｉｎｇ ｔｈｅ ａｕｔｏｍａｔｉｃ ｄｉｓｃｒｉｍｉｎａｔｉｏｎ ｔｅｃｈｎｏｌｏｇｙ ｏｆ ＡｕｔｏＣＡＤ ｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌ《 ；ｇｒａｐｈｉｃｓ ｂａｓｅｄ ｏｎ ｔｈｅ ｕｎｉｔ ｔｅｓｔ ｍｏｄｅ，ｔｈｅ ｐｒｏｂｌｅｍｓ ｏｆａｕｔｏｍａｔｉｃ ｂａｔｃｈ ｏｐｅｎｉｎｇ ｆｔｏｈｅ ＤＷＧ ｒａｇｐｈｉｃｓ ａｒｅ； ｝ｉ｝ｎｇ ｔｈｅ ＳＡＴ ｉｆｌｅｓ，ｅ ｔｒａｃｔｉｎｇ ｔｈｅ ｄｔａ ｉｎｆｏｒｍｔａｉ。， ｏｆ ｔｈｅ ｔｈｒｅｅ—ｄｉｍｅｓｎｉｏｎｌａ ｓ。ｌｉｄ ａｎｄ ｓｔｏｒｉｎｇ ｉｔ ， ｔｈ ｌ ｄａｔａｂａｓｅ ｗｉｔｈ ｔｈｅ ｐｒｏｒｇａｍｍｉｎｇ ｌａｎｇｕａｇｅ ＶｉｓｕａｌＢｂａｓｉｃ ａｎｄ Ｖｉｓｕａｌｌｉｓｐ ｌｎａｕａｇｇｅ ａｓ ｗｅｌｌ ａｓ ｇｒａｐｈｉｃ ｆｏｒｍａｔ ｌ ；ｃ。ｎｖｅｒｓｉ。凡ｔｅｃｈｎｏｚ。ｙｇ ｎａｄ ｅ ｔｒａｖｅｒｓａｌ ｅａｒｃｈ ｍｅｔｈ。ｄ。 ｇｒｐａｈ ｃ ｄｔａｎ ｉｎｆｏｒｍａｔｉ。ｎ．Ｔｈｕｓ ｔｈｅ ａｕｔ。，ｎａ ｃ ｅ一｛ ｝ｌ ａ ｌｕａｔｉｏｎ，，℃ｓＭ ｄｉｓｐｌａｙ ｎａｄ ｄｔａ。。“ｔｐｕｔ ｏｎ ｔｈｒｅｅ—ｄｉｍｅｓｎｉ。ｎａｌ ＣＡＤ ｇｒａｐｈｉｃｓ ｉｓ，℃ ｅ ｃ。，ｎ６ ｉｎｇ 西 ￡ ｅ ｌ ｃｏｍｐａｒｉｓｏｎ ｍｅｔｈｏｄ ｏｆｔｈｅｆｅｔａｕｒｅ ｄａｔａ，ｗｈｉｃｈ ｌａｙ ａｒｔ ｉｍｐｏｒｔａｎｔｆｏｕｎｄａｔｉｏｎｆｏｒ ｔｈｅ ｔｈｒｅｅ—ｄｉｍｅｎｓｉｏｎａｌ ＣＡＤ ｉ 《 ；ｔｅｓｔ ｔｏ ａｃｈｉｅｖｅ ｎｅｔｗｏｒｋｅｄ，ａｕｔｏｍａｔｉｏｎ ｎｄ ａｐａｐｅｒｌｅｓｓ． 中图分类号：ＴＨ１２２文献标识码：Ａ １ 已Ｉ吉 水平技术测试自动化判卷，取代人为因素多、费时效率低、工作量 随着ｃＡＤ技术的不断发展，企业加深了对自主创新的认 大、易出错的人工评判，提高了评判效率和准确性。 ｌ●。‘ —Ｋ， １ｈｈ ，— ｅｙ１ｈ ｈｗ —ｊ ｏｒ—， １ｄｓ ： Ｔ● ，ｅ —ｓｔｆ ； ｉ Ａ ｕ ｉ ｔｏ ， —ＣＡＤ； ０１Ｔｎｐｈ —ｔｒ ｅ＿—ｅ ●－ｄｉｍ －—ｅｎ—１ １ｓｉ ｏ ，—ｎ—ａ －ｌ；Ａ—０ｕ ｉｔｏ—ｍ Ｐａｔ－ｉ●ｃ｛ ｄ—ｄｉ ｓ－ｃ ｉｒ ｉｍｉｄ ＿ ｎ１ ＇ａ —ｔｉ、—ｏ ｎ—ｉ ；Ｓ１ “●—ＡＴ ；Ａ１ —，Ｃ ｔＩ ｊＳ ， —＋１、 —— ■ Ｐ＋＋１ ’ 识，提出高质量、低成本和短周期产品生产要求，推动了三维数字 ２ ＡｕｔｏＣＡＤ三维 ０试题图形数据分析 化设计制造技术在信息化和工业化的融合，这标志着企业对三维 ＣＡＤ三维上机测试以工程实际应用为原则，主要测试题型 ＣＡＤ技术人才的大量迫切需求１１１。为了保证向企业输送的三维 为主观图形题：简单三维图形绘制、三维实体编辑、精确绘制三维 ＣＡＤ人才的质量，可以通过ＣＡＤ水平技术测试准确评价这些人 图形、三维图形尺寸标注、三维实体装配、三维参数化绘图。这些 对三维ＣＡＤ技术的掌握程度。利用ｖＢ、ｖ１ｉｓｐ编程语言和Ａｃｃｅｓｓ 题包括了ＣＡＤ三维绘图的基础操作方法及技巧，考查测试者对 数据库技术开发自动化、智能化的自动评判系统，实现三维ＣＡＤ 三维ＣＡＤ的掌握程度。 ★来稿日期：２０１１－０２—１５ ★基金项目：河北省科学技术研究计划项目（０４２１３５１３４） ５结论 先用ｓｏｌｉｄｗｏｒｋ软件对主轴进行三维建模，根据主轴参数进 参考文献 ［１］李明艳，王素玉，冯明泉．高速切削温度场有限元动态仿真［Ｊ］．煤矿机 械，１９９８． 行了发热量的计算，确立了热边界条件，建立了主轴温度场分布 ［２］粱允奇．机械制造中的传热与热变形基础［Ｍ］一Ｅ京：机械工业出版社， １９８２． 张伯霖，肖曙红．高速电主轴热态特性的有限元分析［Ｊ］＿航空 及热变形分析的有限元模型，并由模型计算得出主轴稳态运行时 ［３］黄晓明，制造技术，２００３（１０）：２０—２３． 的温度分布及热变形。由分析图得出主轴高速稳态运行时前后轴 ［４］李黎明ＡＮＳＹＳ有限元分析实用教程［Ｍ Ｅ京：清华大学出版社，２００５． １５ Ｊ ＪｉｅＺｈｕ，ＪｕｎＮｉ，ＡｌｂｅｒｔＪ．Ｒｏｂｕｓｔｔｈｅｒｍａｌｅｌｌ＇ｏｒｍｏｄｅｌｉｎｇａｎｄｃｏｍｐｅｎｓａｔｉｏｎ 承处温度升高较大，且前轴承比后轴承温升大，由温度分布不均 引起的主轴热变形在前端沿轴颈方向最大。以上理论和结论为课 ｆｏｒ ＣＮＣ ｍａｃｈｉｎｅ ｔｏｏｌｓ［Ｄ］．Ｄｏｃｔｏｒａｌ ｄｉｓｓｅｒｔａｔｉｏｎ ｉｎ ｔｈｅ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｍｉｃｈｉｇａｎ，２００８． 题后续的主轴热变形综合误差补偿奠定了理论基础。 ［６］范梦吾，张伯霖，郭军，李志英．两种电主轴热态杼ｌ生的比较分析［Ｊ］．航 空制造技术，２００４（５）：８－１１．
