模态及谐响应范例--五轴加工中心耳轴式工作台动力学与仿真技巧.doc

五轴加工中心耳轴式工作台动力学分析与仿真研究 文：?|?2015年08月19日? HYPERLINK javascript:void(0) （0）?|? HYPERLINK javascript:addcollect(1727) \o 点击收藏 （0） 0 引言 随着市场对机械加工质量要求的不断提高，作为机床附件的回转工作台，其台面的加工稳定性也备受该行业设计人员的关注。遇到实际振动、工作台变形等问题，传统的方法就是点对点解决，所以很难预先发现回转工作台台面工作过程中潜在的一些刚度薄弱环节。其中文献对五轴联动机床有限元结构分析找到影响机床稳定性的最低刚度方向，文献对回转工作台台面进行了稳态分析，文献讲述了数控回转工作台的设计过程，但都不能有效发现实际存在问题的潜在规律性。为解决这一难题，本文以课题中设计的耳轴式回转工作台为研究基础，借助Pro/E软件建模和ANSYS软件的分析的强大功能，对耳轴工作台进行动力学分析与仿真，实现对回转工作台深层次的理论研究，揭示台面的固有频率、最低刚度方向、共振频率、最大变形等可能造成潜在安全隐患的工作参数，从而可避免在加工过程中出现问题，使加工平稳、零件满足精度要求。同时可为耳轴式工作台参数选择、结构优化提供了重要参考依据。 1 建立耳轴式工作台的有限元模型 1.1 工作台性能参数 本文分析所选用的耳轴工作台要实现的功能是高精度五轴联动加工，工作台所属立式加工中心以车铣加工为主，也可实现钻、镗等其他工序，可对结构复杂的零件进行高精度五面加工。工作台性能指标如表1所示。 1.2 材料属性设置 设定材料是完成有限元分析的重要环节。本文中，摆动工作台、回转工作台和机座、转子驱动轴选用灰口铸铁HT300，定子和转子选用35ww300硅钢，轴承选用Gcr15，冷却套选用铝合金AL3003，材料参数如表2所示。 1.3 对工作台划分网格 采取自由划分网格法，共得到107641个结点，66288个单元，工作台的模型如图1所示，网格划分情况如图2所示。 2 动力学分析 2.1 模态分析 用有限元方法进行动力学分析和静力学分析的本质区别在于作用在结构上的载荷以及结构产生的位移、应变、应力等，前者与时间相关，而后者时间无关。但是??有限元分析方法步骤是基本类似的，两者之间不同的是动态分析必须依据结构振动理论来建立动力学方程，在动态特性分析中不仅需形成刚度矩阵，还需形成质量矩阵和阻尼矩阵，在此基础上求解特征值问题和动力响应问题。 根据结构振动理论，动力学分析时结构的平衡方程为： 式中［M］为结构质量矩阵，［C］为结构阻尼矩阵，［K］为结构刚度矩阵，｛F(t)}为激振力。求解结构的特征值与特征向量就是求解结构的固有频率与固有振型，这是动态分析基本内容。实际经验证明，阻尼对结构的自振频率和振型的影响不大，所有略去不计。再令激振力为零，则得到系统的振动方程如下： 以上即模态分析时结构的平衡方程。 在此基础上对工作台有限元模型进行模态分析，步骤如下：打开Ansys Workbench软件，将Modal模块拖至Static Structural中的Solution选项上，会弹出如图3所示的模态分析对话框。 Max Modes to find设置为默认前六阶，然后单击Solve选项，对工作台模型进行有限元求解运算。则可得到前6阶固有频率和振型，如图4a～图4f工作台前六阶振型图所示；工作台各阶固有频率如表3工作台固有频率表所示；工作台各阶振型如表4工作台固有振型表所示。 第一阶振型是摆动工作台上下振动，主要取决于摆动工作台Y方向刚度和固有频率；第二、四阶振型是支撑部件左右扭转，变形比摆动工作台更大，对装配体动态性能有更大影响，与其抗扭转刚度相关；第三阶振型是工作台整机以YOZ平面为中心的左右俯仰和扭转，取决于支撑部件上半部分的抗扭刚度；第五阶振型为两侧支撑部件左右俯仰和扭转，与其刚度相关；第六阶振型是工作台整机的振动，涉及到工作台各处的刚度。电机是工作台主要振源，转速最高为150r/min，振动频率为2.5Hz，低于工作台各阶固有频率，说明电机回转运动不会引起共振。 2.2 谐响应分析 谐响应分析用于计算工作台承受简谐变化载荷时结构的稳态响应，目的是验证工作台能否避免共振响应，能否经受不同频率受迫振动的影响。谐响应分析仅计算稳态受迫振动，并不包括初始时刻的瞬态分析，分析中结构响应与载荷的频率都相同。 谐响应分析系统动力学响应实际上是解一个二阶常系数线性微分方程，动力学方程如下： 式中，［M］为质量矩阵，[C]为阻尼矩阵，［K］刚度矩阵，为系统中节点的位移、速度和加速度矩阵，｛F(t)｝为力矩阵，力矩阵与位移矩阵是简谐矩阵。求解该线性微分方程有两种方法：①完全法：使用完全的结构矩阵进行求解；②模态叠加法：通过模态分析得