模态与谐响应分析之三数控机床动态特性优化全解

1、数控机床动态特性优化% 匕发布日期：2014-04-29 来源：金属加工世界 浏览次数：66摘 要：考虑了滚动结合面和主要固定结合面的影响，对机床进行模态分析 和谐响应分析，找出了其薄弱模态；通过分析可知改变各 结合面刚度对整机动态 特性的影响，识别出对整机动态特性影响较大的结合面，并对这些结合面刚度值 进行优化，提高了整机的动态特 性。随着数控机床向高速、精密及复合化方向的发展，要求机床具有较好的抵抗 切削颤振的性能，因此，对机床的动态性能提出了越来越高的要求。 建立机床结 构准确的动力学模型，研究其动态特性，对于了解机床的薄弱环节、进行结构优 化、提高机床的加工精度和稳定性具有重要意义。机

2、床是由多个零部件组成的复杂组合结构， 各零部件之间结合部对动态性 能的影响很大。许多研究表明，机床上出现的振动问题有60%以上是源自结合面，机床的静刚度中30%50%决定于结合面的刚度特性，其阻尼值的 90% 以上来源于结合面的阻尼。因此建立考虑结合面的机床动力学模型，对于准确预测机床的动态性能具有重要意义。机床整机特性预测的分析方法主要有 3种：集中参数法、均质梁法、有限 兀法。其中前两种方法在建模时机床结构过于简化， 计算精度低；有限兀法计算 量大，但计算结果精确，目前已成为主要分析方法。在考虑结合面参数的机床整 机建模中，日本学者吉村允孝建立了双柱立式车床的分布质量梁动力学模型，识别了

3、20个结合面，其中2个为导轨结合面，其余为固定结合面。由于考虑了结 合面特性，因而计算所得结果比较接近实验值。 国内西安理工大学的黄玉美、张 广鹏等比较系统地研究了考虑结合面特性的机床动力学建模方法，他们基于结合面基础特性参数，研究了机床导轨结合部特性的建模方法， 并将其用于机床整机 特性分析的有限元模型中，由于同时考虑了固定结合面和导轨结合面的影响，分析结果与实验值较接近。中国石油大学的赵宏林等从理论上推导出各种结合部的 建模融入方法，并应用在XK712B立式镗铣床的建模计算中，得到了比较准确 的计算结果。王学林、覃文洁、张耀满等也分别建立了考虑结合面的整机动力学 模型。由于上述研究都没有考

4、虑丝杠和轴承结合面轴向刚度的影响，因此对各轴轴向振动的预测还欠精确。1有限元模型的建立1.1几何模型建立文中研究的机床为高速加工中心，主轴最大转速可达24000r/mi n ，该型机床主要由床身、墙体、工作台、十字滑座、滑枕、电主轴组成，十字滑座带动 滑枕和电主轴左右方向（x向）运动，滑枕带动电主轴可在竖直方向（y向）运 动，工作台相对于床身做前后方向（z向）运动。由于机床结构非常复杂，因此在 CAD软件中建立其几何模型，然后将其 导入Hypermesh 中进行网格划分，最后将划分好的网格模型导入ANSYA中进行前处理。在建立几何模型时，为提高网格划分质量，节省计算资源，可对结 构作以下简化：

5、（1）对于结构中比较大的孔洞，在模型中仍然保留，若孔的尺寸较小或处 于边缘的工艺小孔及螺栓孔等可在模型中略去。（2 ）结构中的倒角、倒圆在模型中均作直线化处理。（3）结构中的一些影响单元质量的局部特征可做适当简化处理。整个机床的实体部分采用实体单元 Solid45模拟，丝杠采用Beam188单 元模拟，结合面单元采用 ANSYS中的用户自定义单元Matrix27 模拟。由于B eam188 单元有6个自由度，Matrix27 单元虽然也有6个自由度，但是在处 理结合面单元时，只代入了其3个方向的移动自由度，所以为保证单元在连接 节点处自由度协调，需要在 Beam188 单元与Matrix27

6、单元连接节点处建立 约束方程。整个模型划分了 179723个实体单元，90个Beam188 单元。1.2有限元模型中结合面的等效当将相互连接的两个子结构看成是柔性连接，在受到外力激励时该结合的 两个子结构之间发生相对运动。两个子结构通过结合面连接起来，结合面的柔性 即表现为粘弹特性，故可用弹簧阻尼单元等效。结合部解析时，只能用有限个点 代替实际的结合点来描述结合部的静动态特性，将这有限个能描述结合部特性的点称为等效结合点。找到结合部的等效结合点后，用弹簧阻尼单元代替实际结合 点，则结合面即表现为粘弹特性。该机床共有两种结合面：固定结合面、滚动结合面。固定结合面根据等效 结合点概念可以由结合面的

7、面积和螺栓的型号、数目确定结合部等效结合点的数目。网格划分时在等效结合点位置处生成节点，只需要连接对应的两个节点便可以生成一个Matrix27弹簧阻尼单元。而Matrix27弹簧阻尼单元需要输入的参 数即为每个结合面单元的等效刚度阻尼值。根据前苏联学者M . M . BooapMk ob的结合面在影响锥范围内受到挤压理论，可得到式（1）K0 = n （ © + 2htan 丫 ） 2 - d2K/4（1）式中：K单位结合面刚度，为实验测定值；©螺栓帽外切圆直径；h结合面面板厚度；d螺栓孔直径；丫锥角，一般取 丫 =20 °根据结合面的结合状况（面压、粗糙度、结合面

8、介质、材料等），选择单 位结合面参数后代入式（1 ）即可得出等效参数，其中面压来自于结合面所受的 正压力，包括螺栓预紧力、重力、切削力，因此要注意各种力和力矩的等效。该 研究中采用的单位结合面数据为本课题组实验测定的基础数据。同理可以完成对阻尼的等效。滚动结合面包括导轨与滑块、丝杠副、轴承内外挡圈之间的结合面。 导轨与滑块之间结合面采用法向和侧向弹簧阻尼单元模拟，丝杠副和轴承组采用径向和轴向弹簧阻尼单元模拟。由于滑块质量较小，建模时无需考虑其影响，并 且由于导轨与床身、工作台与滑块之间的固定结合面刚度远高于导轨滑块间的滚 动结合面刚度，故可以不考虑其影响，因此导轨与滑块结合面只需在滑块位置处

9、用一组弹簧阻尼单元来模拟。导轨滑块结合面单元的刚度和阻尼数据采用本课题 组实验测试的数据带入，丝杠副和轴承组采用厂家提供的刚度阻尼数据。共建立 了 384个Matrix27弹簧阻尼单元。1.3边界条件模拟由于加工时该机床一般是直接放置于地面上，因此可以利用 ANSYS中的 单向杆单元Link10模拟与地面垂直的方向上的单向位移约束，而在水平面上， 可以认为是由机床自身重力施加的压力与地面之间形成的切向结合面，故切向约束可以参照结合面切向单元定义。1.4模态分析结果在ANSYS中利用BlockLanczos 法求解机床前30阶模态。前4阶固有 频率和振型结果如表1所示，前4阶振型如图1所示。阶次

10、固有频率/Hz振型描述133.862墙体前后摆动242.562墙体、十字滑座、 滑枕左右摆动356J81十字滑座*滑枕绕z轴旋转456.336十字滑座、滑枕绕x轴旋转1阶振型2阶抵型3阶振型4阶振型圏I 就生阶振里1.5机床结合面刚度对整机动态特性的影响为分析各主要结合面对机床动态特性的影响， 分别将各主要结合面刚度提 高20%和40%以比较整机固有频率的变化。表 2中列出了改变x向导轨法向 和侧向刚度对整机部分阶次固有频率的影响， 限于篇幅原因，其余结合面刚度对 整机固有频率的影响未在表中列出。12345X导轨法向 提高20%a 177a 1880,471320206x导轨法向 提高40%0

11、.3010.32407422,0470.355得轨侧向 提高20%0.032090356L5020.279x导轨侧向 提高40%0.0530.320.532.7250.511¥导轨法向 提高20%0.0120.0730.2210.3730.0936789x导轨法向 提高20%0,0140.5730.050753x导轨法向提高40%0.0241,0170,0791.204x导轨侧向 提高20%0.0050.6350.020.039x导轨侧向 提高40%0.0071.0930Q30.068丫导轨法向 提高20%0.0080.0570.051,076由分析可知，改变x向导轨法向刚度，对整机

12、4, 7, 9 , 28阶固有频率较 为敏感；改变x向导轨侧向刚度，对整机4，7，25阶固有频率较为敏感；改 变y向导轨法向刚度，对整机 6，14，15，18，21，25，29，30阶固有频 率较为敏感；改变y向导轨侧向刚度，对整机10阶和14阶固有频率较为敏感； 改变z向导轨法向刚度，对整机9阶和30阶固有频率较为敏感；改变y轴丝杠轴向刚度，对整机6阶固有频率较为敏感；改变z轴丝杠轴向刚度，对整机7 阶固有频率较为敏感。其余结合面刚度的改变对整机固有频率影响不大。由振型动画可发现，上述对整机固有频率较为敏感的结合面都是导致相应阶次振动主要 原因。因此优化时为重点考虑对象。2谐响应分析2.1机

13、床谐响应分析谐响应分析是用于确定线性结构在承受随时间按正弦规律变化的载荷时的稳态响应的一种技术，分析的目的是计算出结构在谐波激振力下的响应， 即位移 响应与应力响应，并得到系统的动态响应与系统激振力频率的曲线， 称为幅频曲 线。机床动态特性的评价，主要是考察其抵抗切削颤振能力，可以刀具与工件 之间的相对位移最小为评价指标，则有如下公式：R=D/F式中：R 动柔度，R值越小说明机床抵抗切削颤振能力越好；D由激励力引起的主轴上刀具安装点与工作台上相应节点之间的相对位移；F激励力。分析中对机床主轴下端安装刀具位置的节点的 3个方向分别施加幅值为5 00N的简谐力，得到电主轴和工作台之间 3个方向的相

14、对位移幅频曲线。xlO4(a) “7wis®1.60 802575125175225275频率/Hz(b) 向xIO*（c） 7.简困2相时位移幅频曲釵2.2谐响应分析结果机床3个方向的谐响应分析幅频曲线如图 2所示。由相对位移幅频曲线 可知，x方向（左右向）相对位移的最大值发生在 3阶和4阶处，2，21，22 阶处的相对位移也较大；y方向（竖直向）相对位移的最大值发生在 3阶和4 阶处，7阶和12阶处的相对位移也较大；z方向（前后向）相对位移的最大值 发生在4阶处，6，9，24，25阶处的相对位移也较大。根据振型动画可以看出，x方向3阶和4阶的振动主要是由于导轨刚度偏 低而引起的十

15、字滑座和滑枕的旋转造成；2阶振动主要是由于床身、墙体左右摆 动和x轴丝杠轴向刚度偏低造成；21阶和22阶振动主要是由于滑枕自身摆动 造成；y方向4阶振动主要由x向和y向导轨刚度的薄弱导致的十字滑座与滑 枕绕z轴旋转造成；7阶振动主要由于y轴丝杠轴向刚度薄弱环节导致的滑枕上 下运动造成；12阶振动主要由于床身和墙体之间结合面法向刚度及 y轴丝杠轴 向刚度薄弱造成；z方向4阶振动主要由x向和y向导轨刚度薄弱导致的十字滑 座与滑枕绕z轴旋转造成；6阶振动主要由于z轴丝杠轴向刚度薄弱导致的工作 台前后运动造成；9阶振动主要由于x轴上导轨、y轴导轨下端刚度偏低及滑枕 前后弯曲造成；24阶和25阶振动主要

16、由滑枕的自身摆动和电主轴在滑枕内的 旋转造成。wijs湮阶次123456789101112131415频率变化/%0.430.733.593.620.952.572.860.442.727.011.480.751.213.112.22阶次161718192021222324252627282930频率变化/%2.611.342.840.430.941.460.91.271.654.281.210.561.315.062.72频率/Hz(a) x 向0.80246802468/654432L2575125175225频率/Hz(b) y 向3结合面参数优化根据谐响应分析的曲线峰值和模态灵敏度分析结果，可以找出结构中结合 面的薄弱环节为x轴和y轴导轨的法向和侧向刚度，z轴导轨的法向刚度，y轴 和z轴丝杠的轴向刚度，为此将x轴和y轴导轨法向刚度提高40%，z轴导轨 法向刚度提高20% ,y轴和z轴丝杠轴向刚度提高30%后分别进行模态分析和 谐响应分析。表3列出了改变结合面刚度后整机的固有频率变化比，表4为