M7475B平面磨床立柱结构ANSYS有限元分析本科毕业论文

《M7475B平面磨床立柱结构ANSYS有限元分析本科毕业论文》由会员分享，可在线阅读，更多相关《M7475B平面磨床立柱结构ANSYS有限元分析本科毕业论文（47页珍藏版）》请在装配图网上搜索。

1、本科毕业论文M7475B平面磨床立柱的有限元分析 目 录前 言4第一章 绪论61.1 机床有限元分析国内外的研究现状61.1.1 国内的研究现状61.1.2 国外的研究现状71.2 本课题主要研究内容71.3 本课题的意义8第二章 M7475B平面磨床简介及立柱结构受力分析92.1 平面磨床结构简介92.2 立柱磨头受力概况102.2.1 砂轮速度的计算：102.2.2 轴向磨削力的计算：112.2.3 切向磨削力和法向磨削力的计算：112.3 磨床振源频率的确定122.4 本章小结13第三章 立柱结构有限元模型的建立133.1 PRO/E与ANSYS的连接133.2 立柱结构建模方法143.

2、3 立柱有限元模型153.4 本章小结15第四章 立柱结构有限元静力学分析164.1 有限元方法简介164.2 ANSYS软件简介164.2.1 ANSYS的组成及主要技术特点:164.2.2 ANSYS结构分析过程174.3 立柱结构的ANSYS结构刚度分析174.3.1 定义单元类型174.3.2 定义材料属性174.3.3 网格划分184.3.4 施加约束194.3.5 施加载荷194.3.6 计算结果204.4 本章小结25第五章 模态分析255.1 磨床动态特性参数265.2 模态分析的基本思想265.3 模态分析的基本理论275.4 模态分析的一般过程275.5 M7475B平面磨

3、床立柱的模态分析285.5.1 定义单元类型285.5.2 定义材料属性285.5.3 网格划分285.5.4 施加约束295.5.5 分析计算295.5.6 观察结果295.6 本章小结32第六章 立柱结构优化设计336.1 优化设计概况336.1.1 优化设计的发展336.1.2 优化设计的概念336.1.3 优化设计的经济效益336.1.4 优化设计发展方向346.2 机床结构设计准则346.3 结构改进优化356.3.1 立柱结构优化方案一356.3.2 立柱结构优化方案二396.4 本章小节43第七章 结论与展望437.1 结论437.2 论文不足之处447.3 展望44参考文献45

4、致谢46 前 言据权威部门的一项统计资料显示，目前我国乡以及乡以上独立核算的机械行业企业金属切削机床拥有量已达300万台，高居世界各国之首。此外，再加上非独立核算的和乡以下企业的拥有量，机床总数当在500万台左右，这是一种雄厚的生产资源。然而，另一方面，我国机床完好率之低也是惊人的。某机械工业集中地区有关部门的调查显示，中、小型机械企业“带病”工作的机床竟有60%。这是造成企业产品质量低、经济效益差的一个重要原因。磨床加工车间废品率上升，经常使工厂受到严重经济损失，研究分析发现，造成此种后果的重要原因一是磨床完好率差，故障频繁，精度不够，得不到及时的调整与维修，二是在机械工业以及金属切削加工技

5、术不断发展的今天，进行高速磨削提高生产效率以及加工精度已经成为磨床的发展趋势，磨床的一些主要结构部件的前几阶固有频率可能会处于工作频率范围之内，将导致磨床的共振。 “工欲善其事，必先利其器。”改变上述状况，必须改善目前企业金属切削机床基础部件的刚性，提高机床抗震性，只有这样才能从根本解决加工精度和机床寿命问题。因此利用大型有限元分析软件对磨床床身部件进行静力学分析和动力学分析以改善上述情况变得十分必要。ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件，能与多数CAD软件接口，实现数据的共享与交换。利用ANSYS软件能进行复杂的应力应变分析的巨大优势，通过大型三维

6、建模软件Pro/ENGINEER建立立柱模型，导入ANSYS软件进行分析，可以清楚并定量的表现出部件变形的情况，从而为提高部件刚性提供理论依据以及数据支持。振动现象也是机床设计中所面临的问题之一，它能造成加工误差，影响零件的加工精度。模态分析主要用于确定结构或机器部件的振动特性。建立机床基础部件三维模型， 利用大型有限元分析软件ANSYS，对机床部件进行模态分析，得出床身部件前几阶固有频率和振型，了解床身部件的各阶振动模态的特点，对于我们研究床身部件的动态特性是十分必要的，有利于机床床身系统的整体设计。提高机床零部件的前几阶固有频率是提高机床刚性、避免共振、降低振幅的有效措施。M7475B型平

7、面磨床是机械工业中广泛使用的重要金属磨削加工工具，采用砂轮端面进行磨削，磨削面积大，立柱结构受力情况极为复杂，本课题用三维CAD软件Pro/ENGINEER对M7475B型磨床立柱结构进行建模，并将三维模型导入ANSYS有限元分析软件以实现数据的共享和交换，利用其强大的分析能力对M7475B型磨床立柱结构进行应力及应变分析，探讨通过改变立柱结构以增加立柱结构刚性, 并基于ANSYS有限元方法对立柱结构进行模态分析，探讨不同结构对提高立柱前几阶固有频率的影响，通过提高前几阶固有频率提高立柱结构抗振性，对提高磨床加工质量和精度有重要意义。第一章 绪论1.1 机床有限元分析国内外的研究现状 国际上早

8、在60年代初就开始投入大量的人力和物力开发有限元分析程序，但真正的CAE软件是诞生于70年代初期，而近15年则是CAE软件商品化的发展阶段，目前流行的CAE分析软件主要有NASTRAN、 ADINA 、ANSYS、ABAQUS、MARC、MAGSOFT、COSMOS等。ANSYS软件致力于耦合场的分析计算，能够进行结构、流体、热、电磁四种场的计算，已博得了世界上数千家用户的钟爱。同时，21世纪ANSYS在机械行业的广泛应用已给传统机械工业带来了新的革命，更高的可靠性设计、更好的抗震性设计等均是ANSYS给机械工业带来的惊喜，特别是在机床结构和抗震性设计方面国内外已有大量成果斐然的研究。1.1.

9、1 国内的研究现状东南大学机械工程系的倪晓宇、易红等利用有限元法对机床床身进行静、动态分析，并使用渐进结构优化算法对床身结构进行基于基频约束和刚度约束的拓扑优化，为 ESO方法在机床大件结构拓扑优化中的应用做了有益的尝试1；内蒙古工业大学的杨明亚等建立立柱的三维有限元模型，利用大型有限元分析软件ANSYS对立柱部件进行了模态分析,得出了立柱前五阶固有频率和振型2；浙江大学现代制造工程研究所的杨晓京等基于ANSYS有限元分析软件对XK640 数控铣床的立柱进行结构优化，比较了四种结构形式的立柱的动力特性，确定了在XK640数控铣床中选用内侧加强筋结构立柱，提高了机床设计水平3；东南大学机械工程系

10、的伍建国等在对M2120A原机床床身动态测试的基础上，建立床身的有限元分析模型，并对床身进行有限元计算，找出原床身设计中的缺陷，从而对各种改进后的机床身再进行有限元分析，通过多方案的比较，得出最优设计方案4 ；东北大学机械工程与自动化学院的张耀满等在CHH6125卧式车削加工中心的研制开发过程中，在产品设计阶段对其采用有限元分析方法，对机床的原有结构进行动力学分析计算、对机床原有机构进行适当改进，提高机床的动力学性能，并对机床进行了试验，对分析结果进行验证5；西安工业大学机电工程学院的朱育权等通过ANSYS软件建立1CL50机床立柱的几何模型，采用四面体单元对立柱进行网格划分，分析了1CL50

11、立柱的一阶、二阶和三阶振型，得出一阶振动为整机摇晃和横断面内弯曲振动，二阶振型为横断面内扭曲振动加垂直方向弯曲振动，三阶振型为横断面内弯曲振动加垂直面内弯曲振动，分析指出了加工过程中应该避开的激振频率，避免一阶振动应加十字型筋板，避免二阶振动应加对角交叉筋板，避免三阶振动应加菱形筋板6。1.1.2 国外的研究现状Mohammed Alfares等研究了磨床动载情况下,对应磨削力的变化情况,以及对工件材料的影响7。N.zhang和I.Kirpitchkenko建立了磨削过程的动态模型,找出了磨削力和模型固有频率变化之间的对应关系,给出了磨削力的估算公式8。1.2 本课题主要研究内容本课题用三维C

12、AD软件Pro/ENGINEER对M7475B型磨床立柱结构进行建模，并将三维模型导入ANSYS有限元分析软件以实现数据的共享和交换，利用其强大的分析能力对M7475B型磨床立柱结构进行应力及应变分析，探讨通过改变立柱结构以增加立柱结构刚性, 并基于ANSYS有限元方法对立柱结构进行模态分析，探讨不同结构对提高立柱前几阶固有频率的影响，通过提高前几阶固有频率提高立柱结构抗振性。1.3 本课题的意义本课题利用ANSYS有限元分析软件对M7475B平面磨床立柱结构进行有限元的静力学分析和模态分析。ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件，能与多数CAD软件接

13、口，实现数据的共享与交换。利用ANSYS软件能进行复杂的应力应变分析的巨大优势，通过大型三维建模软件Pro/ENGINEER建立立柱模型，导入ANSYS软件进行分析，可以清楚并定量的表现出部件变形的情况，从而为提高部件刚性提供理论依据以及数据支持。振动现象也是机床设计中所面临的问题之一，它能造成加工误差，影响零件的加工精度。模态分析主要用于确定结构或机器部件的振动特性。建立机床基础部件三维模型， 利用大型有限元分析软件ANSYS，对机床部件进行模态分析，得出床身部件前几阶固有频率和振型，了解床身部件的各阶振动模态的特点，对于我们研究床身部件的动态特性是十分必要的，有利于机床床身系统的整体设计。

14、通过结构优化设计以提高机床零部件的前几阶固有频率是提高机床刚性、避免共振、降低振幅的有效措施。第二章 M7475B平面磨床简介及立柱结构受力分析2.1 平面磨床结构简介平面磨床是用磨削方法加工工件的机床,可以进行不同精度的各种表面的加工。同一台磨床可以一次性完成粗、精磨削,减少了工序间的搬运、装夹等辅助时间,具有较高的生产效率。平面磨床一般由床身、工作台、磨头、横向进给机构、升降机构、液压系统、冷却系统、电气系统等组成9 。 按照磨头主轴位置和结构布局的不同,主要分为卧轴矩台平面磨床 、卧轴圆台平面磨床、立轴矩台平面磨床和立轴圆台平面磨床四种类型。M7475B平面磨床属于立轴圆台平面磨床，有圆

15、形电磁工作台和立式磨头,采用砂轮端面磨削。该机床为高效率的平面磨床，主要是用于粗磨毛坯或磨削一般精度的工件，适用于成批生产。磨头的回转、机动进给和快速升降、工作台的回转和移动，都分别采用单个电机驱动。M7475B平面磨床结构如图所示。图2-1 M7475B平面磨床结构简图2.2 立柱磨头受力概况M7475B平面磨床M7475B平面磨床属于立轴圆台平面磨床，采用砂轮端面磨削。它的立柱结构主要受砂轮架的重力和端面磨削时产生的磨削力，其中端面磨削时产生的磨削力分为切向力，法向力，和纵向进给产生的轴向力。在本课题中估算砂轮架的重力G=1500N,磨削力的相关计算如下：2.2.1 砂轮速度的计算：砂轮速

16、度（m/s）砂轮直径（mm）砂轮转速（r/min）本课题中砂轮直径为450mm，为970r/min，代入数据计算得出：=22.86m/s2.2.2 轴向磨削力的计算：本课题考虑磨床加工铸铁时的受力情况,轴向磨削力的计算公式如下：10去除单位体积的磨屑所需的能，(见表1)，kgfmm 砂轮线速度，mms磨粒为圆锥形时的锥顶半角， 计算时一般取R 电磁工作台半径，mm工作台单位时间所转动的圈数，rs工件单位时间垂直进给量, mm/s表1 各种工件材料的值（GB80RA和TL80RA）工件材料花岗岩硅纯铁铸铁高速钢(kgf/mm)4305502804901800本课题中取=1mm/min,电磁工作台

17、的转速=20r/min，电磁工作台半径R=375mm，取最大值，代入公式计算得出=1330N（轴向磨削力）。2.2.3 切向磨削力和法向磨削力的计算：根据切向力估算公式，式中为砂轮的切向力, 为电机的输入功率（kw）, 为电机的传动效率, 为砂轮的转速（r/s）, 为砂轮的直径(mm)。所研究磨床的主电机功率是25kw,砂轮转速16.17 r/s,砂轮直径450mm 。将相应的数据代入公式可得砂轮的切向力=1083N,根据文献 22工程上加工铸铁材料工件时法向力和切向力之比/=0.35，求得法向力=3093N 。磨床立柱受力分析示意图2-2、2-3。 图2-2 磨床立柱受力分析示意图 图2-3

18、 磨床立柱受力分析示意图 2.3 磨床振源频率的确定磨床在工作过程中受到多种激振频率的影响,其中以主轴电机和机动进给驱动电机的影响最为明显。M7475B平面磨床所使用的主轴电机转速为970r/min, 机动进给驱动电机的转速为1410r/min,电机转动引起的受迫振动的频率可由共求出，则主轴电机引起的频率为16.17Hz,机动进给驱动电机引起的频率为23.5Hz。电机在使用过程中,由于转子绕组不对称,使得定子和转子主磁力波相互作用的径向分量引起振动,即,由成对磁极产生的电磁拉力引起振动,其频率是电机转动频率的2倍。据此,考虑磁拉力的影响,主轴电机引起的振动频率为32Hz, 机动进给驱动电机引起

19、的振动频率为47Hz。2.4 本章小结本章针对M7475B平面磨床采用立轴圆台端面磨的特点，分别近似模拟出立柱所受磨削力的示意图，并查阅相关文献计算出磨削力的大小。在这里有一点需要说明的是，在本课题中法向磨削力的方向近似与立柱结构的对称中心相平行，而在实际加工情况下法向磨削力的方向是与立柱对称中心偏离一个微小角度的。这样做是为了简化计算，并且简化对于结果影响不大。第三章 立柱结构有限元模型的建立3.1 PRO/E与ANSYS的连接ANSYS软件提供了与大多数CAD软件进行数据共享和交换的图形接口,ANSYS自带的图形接口能识别IGES、ParaSolid、CATIA、Pro/E、UG等标准的文

20、件，使用这些接口转换模型的方法很简单,只要在CAD中将建好的模型使用另存为或者导出命令,保存为ANSYS能识别的标准图形文件,通常使用的有IGES和ParaSolid文件，在ANSYS中使用File-Import导入模型,然后进行模型拓扑结构修改。对于CATIA和Pro/E等CAD软件,ANSYS能直接识别它们的文件,不需要另存其他格式的文件。使用这些图形接口虽然快速方便,但是往往会出现很多问题,甚至会发生不能识别的问题。例如IGES文件是ISO标准中规定的标准图形交换格式之一,IGES作为一种表达产品数据并将其转换成中性文件格式的行业标准,实现文件之间的交换具有很大的优势,但是ANSYS对I

21、GES的支持不够,在导入IGES文件的时候,无法识别小的几何元素,造成所生成拓扑结构不连续,无法生成实体,导入的模型只是由一些面组成，而且,ANSYS读入IGES所需的时间很长。ANSYS软件为了扩充与其他CAD软件之间的接口的功能,开发设置了ANSYS与Pro/E的连接模块,安装的时候必须选择ANSYS Connection For Pro/ENGINEER模块(代号82)，这个模块的功能就是将Pro/E里面的模型生成ANSYS的命令流文件,这个命令流文件的扩展名是anf,ANSYS读入anf文件后,即可自动将模型建好,而且这样建立的模型几乎没有误差。3.2 立柱结构建模方法利用Pro/EN

22、GINEER三维建模软件建立模型如下： 3.3 立柱有限元模型由于磨床在工作的过程中,工作状况比较复杂,受多方面因素的影响,所以进行有限元分析时,为简化计算,需要假设认定磨床材料是各向同性材料,密度均匀分布,在工作过程中始终处于弹性阶段，并且假定位移和变形都是微小的12。磨床立柱结构模型相对复杂，在导入ANSYS进行分析前需要对模型进行相应的简化，如去除相应圆倒角、凸台、螺钉孔、销孔、圆角以及退刀槽等对分析结果影响不大的细节结构13，在本课题中根据实际情况去除了大部分的圆倒角，这样可以简化电脑运算时间，并且对有限元分析结果影响甚微。有限元模型的生成如下图:3.4 本章小结本课题采用三维CAD软

23、件Pro/ENGINEER对M7475B型磨床立柱结构进行建模，并将三维模型导入ANSYS有限元分析软件以实现数据的共享和交换，生成有限元模型。第四章 立柱结构有限元静力学分析4.1 有限元方法简介有限元方法是用于求解工程中各类问题的数值方法。结构强度、刚度分析中的静力、线性或者非线性问题，热传导中稳态、瞬态或者热应力问题，以及流体力学和电磁学中的很多问题都可以用有限元方法解决。有限元方法的基本步骤如下:(1) 将实际求解范围离散化，即将求解域划分成节点和单元。(2) 选择合适的形函数，即选择一个用单元节点解描述整个单元解的连续函数。(3) 对每个单元建立单元刚度矩阵。(4) 按照一定节点编码

24、顺序，将各个单元刚度矩阵叠加以构造结构整体刚度矩阵。(5) 写出以节点自由度（DOF）为未知量的结构整体刚度方程，并将边界条件、初始条件应用到方程中。(6) 求解步骤（5）中得到的方程组，以得到节点上的自由度值。(7) 根据节点的值和形函数，得到其他的物理量。例如，应力、支座反力、弯矩图、热流量等。4.2 ANSYS软件简介ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。由世界上最大的有限元分析软件公司之一的美国ANSYS开发，它能与多数CAD软件接口，实现数据的共享和交换，如Pro/Engineer, NASTRAN, Alogor, IDEAS, Auto

25、CAD等， 是现代产品设计中的高级CAD工具之一。4.2.1 ANSYS的组成及主要技术特点:软件主要包括三个部分：前处理模块，分析计算模块和后处理模块。前处理模块提供了一个强大的实体建模及网格划分工具，用户可以方便地构造有限元模型；分析计算模块包括结构分析（可进行线性分析、非线性分析和高度非线性分析）、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析，可模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力；后处理模块可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示（可看到结构内部）等图形方式显示出来，也可将计算结果以图表、

26、曲线形式显示或输出。软件提供了100种以上的单元类型，用来模拟工程中的各种结构和材料。4.2.2 ANSYS结构分析过程有限元分析的一般流程为：（1）从三维实体建模模块进入有限元分析模块。（2）在实体上施加约束。（3）在实体上施加载荷。（4）计算（包括网格自动划分），解方程和生成应力应变结果。（5）分析计算结果，单元网格、应力或变形显示。（6）对关心的区域细化网格、重新计算。4.3 立柱结构的ANSYS结构刚度分析4.3.1 定义单元类型由于磨床立柱模型比较复杂，不宜简化为板、壳单元的有限元模型,需要采用三维实体单元对磨床结构进行网格划分。本课题中采用软件提供的solid92单元进行网格划分

27、。该单元为10 节点四面体结构线性单元,每个节点有沿X 、Y和Z方向的三个平移自由度, 并且单元有可塑性、蠕动、膨胀、应力钢化, 大变形, 和大张力的能力，如图4-1。4.3.2 定义材料属性M7475B磨床立柱采用材料HT300，材料详细信息如下表表4-1 15结构刚度分析所需材料属性信息材料名称杨氏弹性模量EX泊松比PRXYHT300 N/0.27图4-1 SOLID92 几何模型4.3.3 网格划分由于磨床立柱结构比较复杂,建模的时候出现了很多小线段,采用智能网格划分时容易在不重要的局部结构上产生过多的单元，而采用整体控制单元尺寸的自由网格划分方法,单元的数目可明显减少，这样就简化了计算

28、机进行分析计算的时间。因此,进行网格划分时,合理选择整体单元尺寸为80mm,采用自由网格划分,最后得到的平面磨床整体结构的有限元模型如图4-2所示。图4-2 立柱有限元模型网格划分图4.3.4 施加约束因立轴圆台平面磨床磨削工件分粗磨和精磨两种。粗磨时，为提高磨削效率，可将砂轮轴在垂直平面内旋转一个角度，即用倾斜的砂轮端面磨削；精磨时，必须用水平砂轮磨削。为适应砂轮主轴的旋转角度，本机床采用立柱三点调整装置，即在立柱与床身的结合处安装三套螺钉调整装置，每套相隔距离一致。所以本课题中在三个螺钉调整装置孔的圆柱面上施加约束，定义ALL DOF值为0 ，约束后的结果如图4-3所示。图4-3 施加约束

29、后立柱有限元模型4.3.5 施加载荷立柱主要承受砂轮架重力、磨削力以及重力和磨削力所带给它的弯矩，将重力以及磨削力以均布力的形式施加到立柱上，并且将弯矩简化成力偶矩，选择合适的力臂，经计算得出力偶中力的大小，同时为了避免单个力施加在单个节点上可能引起的应力集中，将力偶简化成均布力施加在多个节点上，如此对立柱进行载荷的施加，结果如图4-4所示。 图4-4 施加载荷后立柱有限元模型4.3.6 计算结果进入求解器Main Menu Solution Solve Current LS。求解完成后，进入通用后处理器观察计算结果。床身变形如图4-5、4-6、4-7、4-8所示。图4-5 X方向位移变形云图

30、图4-6 Y方向位移变形云图 图4-7 Z方向位移变形云图 图4-8 立柱总位移云图从位移云图中可以看出，立柱导轨位移最大，且导轨右半部分有明显位移，达到0.01毫米，这与磨床电磁工作台的转向相符合。 图4-9 X方向应力云图图 4-10 Y方向应力云图 图 4-11 Z方向应力云图 图 4-12 立柱总应力云图从应力云图（图4-9、4-10、4-11、4-12）来看，该型号磨床立柱受载后平均应力不高，基本不会形成应力集中，总体来说满足强度要求。4.4 本章小结本章对M7475B平面磨床的立柱结构进行了静力学分析，综合分析结果，需要注意以下几点：（1）磨床立柱结构在静载下的变形量都比较小,最大

31、变形量在0.01mm左右 ；（2）在实际应用中，M7475B平面磨床立柱的受力情况随加工工件的不同而不同，本课题只考虑了加工铸铁工件时的受力情况，对于加工其它工件时立柱受力，有待进一步研究探讨。（3）实际加工过程中，磨床立柱还会收到热应力的影响，在本课题中也作了忽略。第五章 模态分析机床的动态特性是影响机床性能的重要因素,直接影响着工件的加工精度和生产效率,随着机床加工性能的不断提高,对机床动态特性的要求也愈来愈高。所以,有必要在静态分析的基础上对磨床结构进行进一步的动态分析,包括对磨床立柱原始结构及其改进结构的模态分析。结构的固有频率是衡量结构动态特性的重要指标,磨床立柱本身的固有频率对整机

32、的动态特性有重要影响，在远离振源频率的前提下,提高立柱的固有频率可以改善整机的动态特性,从而减小磨床加工过程中的变形14。5.1 磨床动态特性参数磨床的动态性能是指机床系统在振动状态下的特性,是一项综合性的指标,主要指标有：1.固有频率和主振型,固有频率表示弹性系统每秒的振动次数（单位Hz）,可以用下式来表示：其中为k系统刚度,m为系统质量。多自由度的振动系统有多个主振动,每个主振动中,系统各个坐标之间振动位移的比值,称为主振型。2.动刚度,是指系统抵抗动载荷引起变形的能力,可以表示为：式中,F为激振力,A为共振振幅。一般来说动刚度要小于静刚度,而且随着阻尼比的变化而变化。3.当量静刚度,当系

33、统的静刚度用动刚度和阻尼比来表示时,称为当量静刚度。当量静刚度是指结构发生共振时的静刚度。可以用下式表示：由上可知,机械系统的振动特性主要取决于：(1)系统的静刚度K；(2)激振力条件,如激振频率、激振力F的大小等；(3)系统的物理特性,如质量、固有频率等。5.2 模态分析的基本思想模态分析属于结构动力学分析，它与静力分析的主要区别是动力学分析需要考虑惯性力和运动阻力的影响。5.3 模态分析的基本理论模态分析用于确定设计结构或机器部件的振动特性（固有频率和振型），即结构的固有频率和振型，它们是承受动态载荷结构设计中的重要参数。同时，也可以作为其它动力学分析问题的起点，例如瞬态动力学分析、谐响应

34、分析和谱分析，其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析，后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。ANSYS产品家族中的模态分析是一个线性分析。任何非线性特性，如塑性和接触（间隙）单元，即使定义了也将被忽略。ANSYS提供了七种模态提取方法，它们分别是子空间法、分块 Lanczos法、PowerDynamics 法、 缩减法、非对称法、阻尼法和QR阻尼法。阻尼法和QR阻尼法允许在结构中存在阻尼，本课题中采用的是子空

35、间法。5.4 模态分析的一般过程模态分析的主要步骤为：(1) 建立有限元模型模态分析所用的离散模型与结构静态分析时所用的模型可以相同，否则就应该在前处理中定义单元类型，单元实常数，材料特性和模型的几何性质。需要指出的是，建立模型时所定义的单元应该是线性单元。材料模型可以是线性的、各向同性的或各向异性的，但是必须通过弹性模量和密度或者以其他方式对材料的刚度和质量进行定义。(2) 施加载荷并进行求解这一步主要定义分析类型、指定分析设置、定义载荷和边界条件、指定加载过程设置，然后进行固有频率的有限元求解，在得到初始解后，应该对模态进行扩展以供查看。(3) 扩展模态从严格意义上来说，“扩展”这个词意味

36、着将减缩解扩展到完整的DOF集上。“减缩解”常用主DOF表达。而在模态分析中，我们用“扩展”这个词指将振型写入结果文件。也就是说，“扩展模态”不仅适用于Reduced模态提取方法得到的减缩模型，而且也适用于其他模态提取方法得到的完整振型。因此，如果想在后处理器中查看振型，必须先扩展之(也就是将振型写入结果文件)。(4) 观察结果和后处理模态分析的结果(即扩展模态处理的结果)被写入到结构分析的结果文件Jobname.RST中。分析包括：固有频率、已扩展的振型、相对应力和力分布。可以在POST1/POST1即普通后处理器中观察模态分析结果。5.5 M7475B平面磨床立柱的模态分析5.5.1 定义

37、单元类型由于磨床立柱模型比较复杂，不宜简化为板、壳单元的有限元模型,需要采用三维实体单元对磨床结构进行网格划分。本课题中采用软件提供的solid92单元进行网格划分 。该单元为10 节点四面体结构线性单元,每个节点有沿X 、Y和Z方向的三个平移自由度, 并且单元有可塑性、蠕动、膨胀、应力钢化, 大变形, 和大张力的能力。5.5.2 定义材料属性 磨床立柱采用材料HT300，材料详细信息如下表: 表5-1 模态分析所需的材料属性信息材料名称密度 Density杨氏弹性模量EX泊松比PRXYHT3007800 kg/ N/0.275.5.3 网格划分磨床立柱结构比较复杂,建模的时候出现了很多小线段

38、,采用智能网格划分时容易在不重要的局部结构上产生过多的单元，而采用整体控制单元尺寸的自由网格划分方法,单元的数目可明显减少，这样就简化了计算机进行分析计算的时间。因此,进行网格划分时,合理选择整体单元尺寸为80mm,采用自由网格划分,最后得到的平面磨床整体结构的有限元模型如图4-2所示。5.5.4 施加约束立轴圆台平面磨床磨削工件分粗磨和精磨两种。粗磨时，为提高磨削效率，可将砂轮轴在垂直平面内旋转一个角度，即用倾斜的砂轮端面磨削；精磨时，必须用水平砂轮磨削。为适应砂轮主轴的旋转角度，本机床采用立柱三点调整装置，即在立柱与床身的结合处安装三套螺钉调整装置，每套相隔距离一致。所以本课题中在三个螺钉

39、调整装置孔的圆柱面上施加约束，定义ALL DOF值为0 。5.5.5 分析计算进入求解器Main Menu Solution Solve Current LS。求解完成后，进入通用后处理器观察计算结果。分析计算磨床立柱各阶振型频率如表5-2所示： 表5-2 磨床立柱各阶振型频率振型阶数1 2345振型频率（Hz）100.74131.60251.04446.99460.80最大变形（mm）1.7691.7412.1675.982.475.5.6 观察结果磨床立柱模态分析各阶振型图如下：图5-1 1阶振型节点位移变形云图图 5-2 2阶振型节点位移变形云图图 5-3 3阶振型节点位移变形云图图 5

40、-4 4阶振型节点位移变形云图图 5-5 5阶振型节点位移变形云图由第二章磨床振源频率的确定可知主轴电机引起的振动频率为32Hz, 机动进给驱动电机引起的振动频率为47Hz。分析计算结果中1阶振型频率为100.74Hz，1阶固有振型频率已经远离驱动电机所引起的振动频率，而2到5阶振型频率更难达到，所以不容易引起磨床立柱的共振。5.6 本章小结本章节通过对磨床立柱的有限元模型进行动力学模态分析，掌握了立柱各阶振型下所引起的变形情况，通过观察振型频率和变形，了解原结构设计存在的缺陷，为后面磨床立柱的结构优化设计，提高磨床加工精度奠定了数据基础。第六章 立柱结构优化设计6.1 优化设计概况6.1.1

41、 优化设计的发展机械优化设计应用的发展历史，经历了由怀疑、提高认识到实践收效，从而引起广大工程界日益重、视的过程。在6070年代，计算机价格昂贵，优化设计在应用实践方面多数限于高等院校、研究所和少数大型企业中开展。从70年代到80年代，计算机价格大幅度下降，优化设计应用的诱人威力，市场竞争日益激化，作为产品开发和更新的第一关是如何极大地缩短设计周期、提高设计质量和降低设计成本已成为企业生存的生命线，从而引起广大企业和设计师的高度重视。用优化设计方法来改造传统设计方法已成为竞相研究和推广并可带来重大变革的发展战略，优化设计在设计领域中开拓了新的途径。6.1.2 优化设计的概念机械优化设计是将机械

42、工程的设计问题转化为最优化问题，然后选择适当的最优化方法，利用电子计算机从满足要求的可行设计方案中自动寻找实现预期目标的最优化设计方案。设计步骤：把实际机械设计问题转化为数学模型，然后根据数学模型的特性，选择某种适当的优化方法及其程序，通过电子计算机，求得最优解。6.1.3 优化设计的经济效益机械优化设计是把数学规划理论与计算方法应用于机械设计，按照预定的目标。借助于电子计算机的运算寻求最优设计方案的有关参数，从而获得好的技术经济效果:(1.)可以降低机械产品成本，提高它的性能;(2)优化设计过程中所获得的大量数据，可以帮助我们摸清各项指标的变化规律，有利于对今后设计结果作出正确的判断，从而不

43、断提高系列产品的性能;(3)用优化设计方法可合理解决多参数多目标的复杂产品设计问题。6.1.4 优化设计发展方向21世纪机械优化的对象已不是简单的机械零部件，而是复杂的零部件、整机直至系列产品和组合产品(统称为复杂系统)。主要向两个方向发展：1.对产品的性能要求由单一性要求发展到多性能和全性能要求，这将导致基于复杂系统多性能和全性能要求的优化。2.从产品的单一设计优化发展到产品的加工制造、使用、维修与管理直至回用的全寿命周期优化。6.2 机床结构设计准则为使机床具有较好的加工性能和经济效益,机床的设计一般应当满足以下的准则：具有较广的工艺范围,能满足各种零件加工的需要较高的精度,包括几何精度,

44、运动精度和传动精度等高可靠性,要求机床在规定的使用期内尽可能少发生故障17。另外,设计不能单纯从结构上考虑,有关人机关系的因素,外观、制造、装配和维修的方便性等问题也都必须予以注意，即,机床的设计必须使其结果在技术上是可以接受的 ,在外观上是令人满意的,并且在经济上是有竞争能力的。机床大件是整台机床的基础和支架,其它零部件以它们作为支承、固定和运动的基础,所以床身、立柱、拖板等大件,是决定整机性能的关键件。机床大件的功能主要有如下方面的要求18：1.静刚度要求高,在最大允许载荷时,变形量不超过规定值；在大件本体移动时,或其它部件在大件上移动时,静刚度的变化应小。2.动刚度要好,在预定的切削条件

45、下工作时,能稳定地达到预期的加工精度。3.连接刚度要好,结构要合理,便于调整和装卸。4.温度场分布合理,工作时的热变形对加工精度影响小。5.导轨面受力合理,耐磨性要好。6.结构设计合理,材质稳定,能长期保持规定精度。7.排屑容易,冷却液和润滑油的渠道畅通,而且不漏油或透油。8.重量要轻,体形要合理,便于吊装和运输。由此可见,一个好的结构不仅要有良好的静态性能,还要有良好的动态和热态性能,在此条件下还要尽量减小质量,节省材料,降低成本。6.3 结构改进优化一般机械结构系统动态优化准则为：提高各阶固有频率；各阶固有频率尽量均布；避免固有频率与外界激励频率一致引起共振；各子结构的动刚度不出现明显薄弱

46、环节。6.3.1 立柱结构优化方案一在本优化方案下，在立柱底座垂直边角出加了一边两条共四条斜筋板（见图6-1），并将改进后的立柱结构导入ANSYS软件进行模态分析，计算磨床立柱各阶振型频率如表6-1。 表6-1 磨床立柱改进型1各阶振型频率振型阶数1 2345振型频率（Hz）108.35135.54259.68448.55467.86最大变形（mm）1.7861.752.1986.3922.808图6-1 磨床立柱改进型1有限元模型表6-2 磨床立柱改进型1各阶振型频率与原结构振型频率的比较振型阶数1 2345原有结构振型频率（Hz）100.74131.60251.04446.99460.80

47、改进型1振型频率（Hz）108.35135.54259.68448.55467.86通过比较磨床立柱改进型1有限元模型各阶振型频率与原设计结构振型频率可以看出，改进型1由于在立柱底座两边加了4个斜筋板，各阶振型均有所提高，影响立柱动态特性的主要因素是立柱的第一、第二阶固有频率，其中尤其以第一阶固有频率的影响最为明显，而在本优化方案中，一阶振型频率比原设计结构一阶振型频率提高了约8Hz，表明本优化方案效果明显，有效地提高了一阶固有频率使其更加远离能引起立柱共振的驱动电机频率。本优化方案立柱结构模态分析结果如图6-2、6-3、6-4、6-5和6-6所示。图 6-2 磨床立柱改进型1一阶振型图 6-

48、3 磨床立柱改进型1二阶振型图 6-4 磨床立柱改进型1三阶振型图 6-5 磨床立柱改进型1四阶振型图 6-6 磨床立柱改进型1五阶振型6.3.2 立柱结构优化方案二在本方案中，对立柱中筋板的厚度和宽度均有所加大，其中筋板厚度加厚8mm，筋板宽度加厚20mm（如图6-7）。在PRO/E中修改原模型，优化完毕后导入ANSYS软件进行模态分析，分析计算后，在通用后处理器中观察各阶振型频率的变化，本优化方案立柱振型频率结果如表6-3所示。表6-3 磨床立柱改进型2各阶振型频率表振型阶数1 2345振型频率（Hz）99.76132.75254.71448.21467.30最大变形（mm）1.751.6

49、892.1985.8183.366本优化方案立柱结构模态分析结果如图6-8、6-9、6-10、6-11和6-12所示。图6-7 磨床立柱改进型2有限元模型图 6-8 磨床立柱改进型2一阶振型图6-9 磨床立柱改进型2二阶振型图6-10 磨床立柱改进型2三阶振型图6-11 磨床立柱改进型2四阶振型图 6-12 磨床立柱改进型2五阶振型磨床立柱改进型2各阶振型频率与原结构振型频率的比较表振型阶数1 2345原有结构振型频率（Hz）100.74131.60251.04446.99460.80改进型2振型频率（Hz）99.76132.75254.71448.21467.30从磨床立柱结构改进型2各阶振

50、型频率与原结构各阶振型频率的比较表中可以看出，除了一阶振型无明显变化外，2到5阶振型均有所提高，在改进方案二下第二到五阶固有振型频率均有所提高，但立柱一阶固有振型频率有微小降低，一阶振型最大位移变形量也有所降低，不过固有频率降低了0.97%，最大位移变形量降低了1.07%,固有频率降低幅度小于最大位移变形量降低幅度，所以该方案亦可行。6.4 本章小节本章通过分析两种不同优化方案，分别比较了两种优化方案下各阶振型频率的变化，改进型1通过在磨床立柱底座两边各加两个斜筋板，筋板厚度100mm，经导入ANSYS软件分析，结果显示各阶振型频率均有所提高。改进型2通过将立柱中的筋板加厚8mm，加宽20mm

51、，经导入ANSYS软件分析，结果显示除了1阶振型无明显变化外，2到5阶振型均有所提高，通过比较两个方案的模态分析结果，方案一效果明显，且实际应用简单，方案二通过加宽加厚筋板，虽然结果表明方案可行，不过没有方案一实际应用简单。第七章 结论与展望7.1 结论 本课题利用ANSYS有限元分析软件建立了M7475B立轴圆台平面磨床立柱的有限元模型，并对立柱结构进行了静刚度分析和动力学模态分析，特别依据模态分析的结果，对磨床立柱结构进行了简单的优化。综合以上分析结果，得出如下结论：1. 磨床立柱结构在静载下的变形量都比较小,在本课题中最大变形量在0.01mm左右 ，且立柱结构在受载情况下平均应力不高，应

52、力集中现象不明显，原立柱设计结构满足静刚度要求；2. 在立柱结构静刚度分析中，本课题考虑了磨头加工时离工作台距离最大时的极限情况，因为在此种情况下，载荷对磨床立柱的影响最大；3. 立柱结构模态分析结果表明，原立柱结构一阶固有振型频率较高，高于磨床驱动电机所引起的激振力的频率；4. 在立柱结构的优化设计中，本课题简单地研究了两种改进方案对立柱各阶固有振型频率的影响，在改进方案一下立柱各阶固有振型频率有明显提高，而最大位移变形量变化不大；在改进方案二下二到五阶固有振型频率均有所提高，但立柱一阶固有振型频率有微小降低，一阶振型最大变形量也有所降低，不过固有频率降低了0.97%，最大位移变形量降低了1

53、.07%,固有频率降低幅度小于最大位移变形量降低幅度，所以该方案亦可行。7.2 论文不足之处由于时间和作者水平等的限制，本课题的研究空间还有待扩展：1. 本课题只对磨床立柱进行了研究，对于磨床床身以及整机的有限元分析还有待研究。2. 由于作者水平有限，在对磨床立柱结构建立有限元模型以及进行受力分析时进行了近似处理。7.3 展望1. 固有频率的减小比例小于变形的减小比例,说明改进结构是可行的,此结论可以推广到类似的机械结构中。2. 采用体单元对结构进行优化设计,迭代过程中模型容易出现不能划分网格的问题,可以尝试使用壳单元模拟筋板等薄壁结构,进行优化分析,避免此类问题的出现。参考文献1倪晓宇等机床

54、床身结构的有限元分析与优化J 制造技术与机床2005(2)：47-502杨明亚等机床立柱动态特性的分析J机械制造与研究2007 (1)：29-313 杨晓京等基于ANSYS的XK640数控铣床立柱优化设J机械设与制造2006(11)：73-744 伍建国等内圆磨床床身结构的动态分析与优化设计J精密制造与自动化2002(2)：25-265 张耀满等高速机床有限元分析及其动态性能试验J组合机床与自动化加工技术2004(12) ：15-176 朱育权等1CL50型机床立柱振动模态分析J西安工业大学报2007(3) ：215-2187. Mohammed Alfares,Abdallah Elshar

55、kawy.Effect of grinding forces on the vitration of grinding machine spindle system.International Journal of Machine Tools&Manufacture,2005,40(2000):2003-2030.8. N.Zhang,I.kirpitchenko,D.k.Linux,Dynamic model of the grinding process,Journal of Sound and Vibration,280(2005):425-432.9. 孙恭寿,等.磨床精化与改造.北京

56、:机械工业出版社1986.10.赵希学.关于立轴圆台平面磨床磨削力的计算问题J.制造技术与机床1992(01):25-26.11邱言龙磨工技师手册.北京：机械工业出版社2002：201-20312赵汝嘉机械结构有限元分析西安：西安交通大学出版社199013王伟伟，翁泽宇，巫少龙大型高效数控铣床有限元模态分析.机床与液压2005(6)：2114胡如夫高精度内圆磨床结构动态优化设计研究中国机械工程2002,13(18) ：154215郑伟中机床的振动及其防治北京：科学出版社198116王月飞机床的动态特性分析河北工业科技2001,18(4) ：1117金属切削机床设计编写组金属切削机床上海：上海科学技术出版社198518叶瑞汉机床大件焊接结构设计北京机械工业出版社1986致谢