压砖机的有限元分析

XX学院设计说明书课 题： 压砖机的有限元分析 子课题: 同课题学生姓名: 专 业 学生姓名 班 级 学 号 指导教师 完成日期 第1章 绪论1.1 课题研究背景液压压砖机是陶瓷工业中用于墙地砖压制成形的必不可少的机械设备。目前我国在实际生产中采用的压砖机类型主要有钢丝缠绕型液压压砖机、开式液压压砖机、闭式液压压砖机等类型，但工作时都普遍存在压砖机变形较大，刚性较差，而压砖机的强度、刚度将直接影响到零件加工精度、压砖机导轨的磨损和模具的寿命等。因此如何优化压砖机结构，提高压砖机的静态、动态特性，同时又能降低压砖机的结构重量，对于液压压砖机的设计尤为重要。随着科学技术的发展，压砖机正向大吨位、高精度和高速度发展。要保证这些大吨位、高精度和压砖机的正常工作，首先应该在设计压力机压砖机时必须保证有足够的强度和刚度，同时考虑工作时的压力机的振动情况。目前我国压力机压砖机的设计至今大多沿用经验、类比的传统设计方法，设计出的压砖机不仅性能差，结构笨重，速度、精度提不高，而且设计周期长，制造成本高，更新换代慢，国产高档次的压力机领与国外压力机相比存在很大的差距。随着中国加入WTO，中国的制造企业的形势将变得更加严峻，并面临更为强大的竞争对手，为此，中国的压力机制造企业必须改变原有的传统设计方法，以先进的设计制造手段作为技术支撑，来提高我国压力机的设计与制造水平，在新的市场环境中积极参与竞争。随着CAD/CAM/CAE技术的日益普及和应用，有限元方法等现代结构分析方法己为工程技术设计人员广为认识和发展，在压砖机设计中得到广泛的应用，并取得了显著的技术经济效益。1.2 国内外全自动液压压砖机研究现状与发展趋势1.2.1国外全自动液压压砖机的发展世界各国生产陶瓷砖除了塑性法、注浆法成型坯体之外，只要是采用颗粒状粉料压力成型工艺的基本上都是走过手工锤打半机械化的摩擦压力锤机械式压力机摩擦液压机成型全自动液压机成型的道路。因此，当今各地企业选用的自动液压压砖机其实是实践经验总结的应用，是目前最普遍最先进的方法，但不是唯一的方法。世界上对陶瓷墙地砖制品与专用装备开发比较早又长期坚持下来的国家，首推意大利，其次是德国、美国。1983年佛山引进投产的生产线配用的是台2000型摩擦液压压砖机，此机的实质还是摩擦锤，液压机构作为操作手段。之后，引进萨克米公司的500、550、600、1000型，西蒂公司的800、1200型，德国莱斯公司的600、1100、1500型等等，在上世纪的80年代末90年代初，才真正出现以液压机构为主的液压压砖机，级别上属于中小型机。中国在1995年之后，出现生产制作500毫米500毫米至1000毫米1000毫米大规格砖的趋势，为此，娜塞堤等公司才研发了5000吨级的大型压砖机并首先进入中国市场。1999年，世界最大的7200吨压砖机进入中国市场。由此可知，意大利、德国、日本生产制造全自动液压压砖机的历史是1520年的事情。这些国家产品的共同点都是以液压机构为工作机构，引入当代的科技成果，逐步完善，最终走向机电液一体化的高技术装备。但他们又都保持着企业的特点，如在液压油路的设计、控制显示手段、机架结构、布料系统设计等。国外的压砖机企业紧跟行业终端产品而开发，紧跟社会的新科技不断提高，不轻易放弃自己的特色等做法应对中国压砖机的研制有重要的借鉴作用。1.2.2 国内全自动液压压砖机的发展我国由国家主管部门立项、自主研究取得成果，并一直影响坚持到今天的现代化全自动液压压砖机的研制应从1986年开始算起。由原国家建材局立项，咸阳陶瓷研究设计院牵头，华南理工大学、佛山陶机总厂联合研制的600型机于1989年通过鉴定，1993年之后，逐步形成批量生产。力泰公司以此机型延伸形成系列化产品，由600吨至7200吨级，直至今天许多型号仍是国内市场销售的主要机型。上世纪90年代，福建海源企业以其充沛的开发能力，引入新技术，最早突破千吨级以上压砖机的研制开发并率先通过国家级鉴定，也先后成功推出超千吨系列产品投放国内市场。1999年，广东科达机电公司一举成功推出3200吨级压砖机并通过国家级鉴定，开创中国设计生产大型压砖机的年代，而发展到现在更是出现了7800吨的压机。因此中国对现代压砖机的研制开发是前后近20年的事，形成了由“力泰”、“海源”、“科达”家能生产大型压机的企业作为龙头，另有十几家生产中小型压砖机的企业并存局面。估计不久的将来，还会出台超万吨级的压砖机，满足高效生产超大规格尺寸砖的需求1。1.3 本课题研究的主要内容与技术方案本课题基于有限元分析软件ANSYS WORKBENCH和PRO/MECHANICA平台，对WL1700全自动液压压砖机关键零、部件进行有限元结构分析。技术方案流程图，如图1.1所示。具体的方案如下：研究WL1700型全自动液压压砖机的相关技术参数、模型特点基于有限元分析软件ANSYS WORKBE和PRO/MECHANICA平台，进行有限元结构分析关键零、部件PRO/E模型有限元分析结果评价图1.1技术方案流程图（1）研究WL1700型全自动液压压砖机的相关技术参数，和模型特点，为分析提供相应的理论基础；（2）利用有限元方法对压砖机进行模态分析，分析出它的振动特性，找出振动中危险的位置；（3）存在最大静载荷时，对机架、活动横梁进行静力学分析，分析结构的应力、应变和变形位移的分布规律，根据分析结果进行强度较核；（4）基于活动横梁的模态分析，对活动横梁动态时域分析，找出应力最集中的部位应力随时间变化的规律；（5）对活动横梁进行疲劳分析，找出疲劳循环次数、疲劳破坏、安全系数、疲劳灵敏性。分析出活动横梁能不能达到疲劳强度的设计的要求；（6）在活动横梁静力学分析的基础上，对活动横梁进行拓扑优化，找出活动横梁能够去除材料的部位，为活动横梁的质量优化提供参考；（7）综合以上的分析结果，对压砖机的结构提出相关的修改方案，以使液压压砖机的结构趋于合理达到优化的目的。44第2章 压砖机有限元分析的理论基础2.1 全自动液压压砖机的工作原理、基本参数及结构特点2.1.1 全自动液压压砖机的工作原理 全自动液压压砖机的工作原理其实可以简化说明为液压机的工作原理1。液压机是利用液压压力能来传递能量，以实现各种压力加工工艺的机器。液压机根据帕斯卡原理制成。1小柱塞，在力F1=4P/d12作用下，将经过连接管道等值传至大柱塞或活塞2上，使砖坯3受到作用力F2=D12P/4，所以可得F2=F1（D2/D1），由于D2D1,则F2F1。由此可知，全自动液压压砖机就是利用在小柱塞上的较小作用力F1，可以在大柱塞上产生很大的作用力 F2 ，其工作原理图2.1所示。 图2.1 全自动液压压机工作原理图1-小柱塞 2-大柱塞 3-坯体 2.1.2 全自动液压压砖机的基本参数 基本参数是全自动液压压砖机的基本技术数据，是根据液压机的工艺用途及结构类型来确定的，它反映了它的工作能力及特点，也基本上定下了它的轮廓尺寸及本体总重。另外，基本参数也是用户选购时的主要数据1。全自动液压压砖机的基本参数有以下内容：（1）公称压力：指液压机名义上能产生的最大力量，它反映了液压机的主要工作能力；（2）最大净空距：指活动横梁停在上限位时，从工作台上表面到活动横梁下表面的距离，它反映了液压机在高度方向上工作空间的大小；（3）最大行程：指活动横梁能够移动的最大距离，应根据工件成型过程中所要求的最大工作行程来确定，它直接影响工作缸和回程缸及其柱塞的长度及整个机架的高度；（4）工作台尺寸：指工作台面上可以利用的有效尺寸，它取决于模具的平面尺寸及工艺过程的安排；（5）回程力：滑块返程时所需要的驱动力；（6）活动横梁运动速度：分为工作行程速度、空行程速度及回程速度；（7）允许最大偏心距：指工件变形阻力接近公称压力所允许的最大偏心值。本课题中分析的WL1700全自动液压压砖机的主要技术参数如表2.1表2.1 WL1700全自动液压压砖机的主要技术参数压机技术性能单位数值公称压制力kN17000模芯顶出力kN210最大行程mm140立柱间距mm1700活动横梁宽度mm680最大填料深度mm60主油缸内最大工作压力MPa33.8液压系统工作压力MPa16主电机功率kW75工作循环次数次/min18工作周期中加压次数次2或3整机质量t462.1.3 WL1700全自动液压压砖机的结构WL1700全自动液压压砖由本体(主机)及液压系统两部分组成，压砖机主机部分,如图2.2所示。它由上横梁、下横梁、四个立柱和十六个内外螺母组成一个封闭架框，框架承受全部工作载荷。工作油缸固定在上横梁上，工作油缸内装有工作柱塞，它与活动横梁连接。活动横梁以四根立柱为导向，在上、下横梁之间往复运动。上模装在活动横梁的下表面上，下模装在下横梁的工作台上。当高压液体进入工作缸后，在工作柱塞上产生很大的压力，并推动柱塞、活动横梁及上模向下运动，使坯体在上模、下模之间产生塑变形。回程油缸固定在下横梁上，其中有回程柱塞，它与活动横梁相连接。回程时,工作油缸为低压，回程油缸为高压，使回程柱塞向上运动，带动活塞横梁回到原始位置，完成一个工作循环。图2.2 WL1700全自动液压压砖本文研究的WL1700型压砖机机架采用螺栓预应力结构，设计时充分考虑到WL1700全自动液压压砖机的工况特点以及机架结构中的横梁和立柱的刚度问题，结构合理，变形小，寿命长。该机的主要特点如下：（1）台面宽，频率快，生产效率高；（2）采用强大的立柱导向，工作精度高；（3）顶出装置采用全液压控制，液压锁模，操作方便；（4）泵站采用全过滤液压装置，确保液压系统清洁度，降低压机故障率；（5）活动横梁的运动速度由比例插装阀控制，线性位移传感器检测位置，控制精度高，可在任意位置进行速度转换控制；（6）装置为独立辅机，可移开，拆装模具方便；料耙由油马达驱动，速度快，比例阀控制，可在任意位置进行速度转换控制，精度高；（7）大型可编程序控制器，友好的人机界面，丰富的信息提示，操作、调节、维护极为方便。2.2 有限元法基本原理2.2.1 有限元法基本思想有限元的基本思想是：将连续的结构离散成有限个单元，并在每一个单元中设定有限个节点，将连续体看作是只在节点处相连续的一组单元的集合体，同时选定场函数的节点值作为基本未知量，并在第一单元中假设一插值函数以求表示单元中中场函数的分布规律，进而将一个连续域中的无限自由度问题转化为离散域中的有限自由度问题6。2.2.2 有限元分析步骤有限元法分析问题的基本步骤11：(1)结构的离散化。离散化就是将要分析的结构分割成有限个单元体，并在单元体的指定点设置节点，使相邻单元的有关系数具有一定的连续性，并构成一个单元的集合体以代替原来的结构。结构离散化时，划分的单元大小和数目应根据计算精度的要求和计算机的容量来决定。(2)选择位移差值函数。为了能用节点位移表示单元体的位移、应变和应力，在分析连续体问题时，必须对单元中位移的分布作一定的假设，即假定位移是坐标的某个简单的函数。选择适当的位移函数是有限单元法分析中的关键。通常采用多项式作为位移函数。 (3)分析单元的力学特性。利用几何方程、结构方程和变分原理最终得到单位刚度矩阵。(4)集合所有单元的平衡方程，建立整体结构的平衡方程。先将各个单元刚度矩阵合成整体刚度矩阵，然后将各单元的等效节点力列阵集合成总的载荷列阵。(5)由平衡方程组求解未知节点位移和计算单元应力。2.2.3有限元法求解步骤对于不同物理性质和数学模型的问题，有限元求解法的基本步骤是相同的，只是具体公式推导和运算求解不同11。有限元求解问题的基本步骤如下：（1）问题及求解域定义。根据实际问题近似确定求解域的物理性质和几何区域。（2）求解域离散化。将求解域近似为具有不同有限大小和形状且彼此相连的有限个单元组成的离散域，习惯上称为有限元网络划分。显然单元越小(网络越细)则离散域的近似程度越好，计算结果也越精确，但计算量及误差都将增大，因此求解域的离散化是有限元法的核心技术之一。（3）确定状态变量及控制方法。一个具体的物理问题通常可以用一组包含问题状态变量边界条件的微分方程式表示，为适合有限元求解，通常将微分方程化为等价的泛函形式。（4）单元推导。对单元构造一个适合的近似解，即推导有限单元的列式，其中包括选择合理的单元坐标系，建立单元试函数，以某种方法给出单元各状态变量的离散关系，从而形成单元矩阵(结构力学中称刚度阵或柔度阵)。为保证问题求解的收敛性，单元推导有许多原则要遵循。对工程应用而言，重要的是应注意每一种单元的解题性能与约束。例如，单元形状应以规则为好，畸形时不仅精度低，而且有缺陷的危险，将导致无法求解。（5）第五步，总装求解。将单元总装形成离散域的总矩阵方程(联合方程组)反映对近似求解域的离散域的要求，即单元函数的连续性要满足一定的连续条件。总装是在相邻单元结点进行，状态变量及其导数(可能的话)连续性建立在结点处。（6）联立方程组求解和结果解释。有限元法最终导致联立方程组。联立方程组的求解可用直接法、迭代法和随机法。求解结果是单元结点处状态变量的近似值。对于计算结果的质量，将通过与设计准则提供的允许值比较来评价并确定是否需要重复计算。2.2.4 有限元法的优点有线元法的特点是适用于求解各种形式(几何上、物理上)复杂的问题，精度高，通用性强，对问题的处理既彻底又系统，适用于采用电子计算机方式。它本是线性问题的解法，但通过迭代法(如牛顿一拉裴森迭代法)也能巧妙地解决非线性问题12。其优点如下：(1)概念浅显，容易掌握。可以在不同的水平上建立起对该法的理解；可以通过非常直观的物理概念来理解；也可以建立基于严格的数学分析的理论。(2)适用性强，应用广泛，几乎适用于求解所有的连续介质和场问题。(3)采用矩阵形式表达，便于编制计算机程序，可以充分利用高速计算机所提供的方便。2.3 计算软件的选择 本课题采用ANSYS WORKBENCH和PRO/MECHANICA 有限元分析软件，对WL1700全自动液压压砖机关键零、部件进行有限元结构分析。2.3.1 ANSYS WORKBENCH软件简介ANSYS WORKBENCH软件是融合结构、热、流体、电磁、声学于一体的大型通用有限元商用分析软件，可广泛应用于核工业、铁道、石油化工、航空航天、机械制造、能源、电子、造船、汽车交通、国防工业、土木工程、生物医学、轻工、地矿、水利、日用家电等一般工业及科学研究。其分析类型有：结构静力分析、结构动力分析、结构非线性分析、动力学分析、热分析、电磁场分析、流体动力学分析、声场分析、电压分析15。ANSYS WORKBENCH软件具有以下三个方面的特点： （1）强大而广泛的分析功能：可广泛用于求解结构、热、流体、电磁、声学等多物理及多场耦合的线性、非线性问题；（2）一体化的处理技术：主要包括几何模型的建立、自动网格划分、求解、后处理、优化设计等许多功能及使用工具；（3）丰富的产品系列和完善的开放体系：不同的产品配套可应用于各种工业领域。ANSYS WORKBENCH分析过程包含3个主要的步骤：（1）创建有限元模型，并定义材料属性和划分网格；（2）施加载荷并求解；（3）查看结果，并分析结果的正确性。2.3.2 PRO/MECHANICA软件简介PRO/MECHANICA是美国PTC公司推出的可以完全实现几何建模和有限元分析的无缝集成软件。用户在PRO/ENGINEER 环境下完成零件的几何建模后，无需退出设计环境就能进行有限元分析17。（1）PRO/MECHANICA 软件包是由以下3个模块组成的。PRO/MECHANICA Structure 结构分析模块，可以进行机械零件、汽车结构、桥梁和航空结构等的结构分析和优化设计。它能够完成的分析种类有静力分析、模态分析、屈曲分析、疲劳分析、非线性大变形分析等。（2）PRO/MECHANICA Thermal 温度分析模块，可以进行零件的稳态和瞬态温度场分析。其分析数据可以返回到结构分析模块，进行灵敏度分析和优化设计。（3）PRO/MECHANICA Motion 运动分析模块，可以进行机构的运动学分析、动力学分析、三维静态分析和干涉检查等。2.4 本章小结本章主要介绍了全自动液压压砖机的工作原理，结构特点以及其基本参数，同时详细论述了有限元法基本思想、有限元分析步骤、有限元法求解步骤及有限元法的优点。为WL1700全自动液压压砖机关键零、部件有限元结构分析打下基础。第3章 全自动液压压砖机的模态分析3.1 模态分析的概述模态分析用于提取结构的固有频率和固有振型，模态分析的重要性不仅是为了避免机械零件和结构在工作时发生共振，事先算出它们的固有频率和振型，而且还是分析结构动态响应和其他动力特性的基础。在结构设计时，固有频率和振型对于动载荷下的响应是非常重要的参数，因为结构的基本模态信息能有助于得出其动力响应的特征。3.2 压砖机模态分析全自动液压压砖机在实际的工况中，当活动横梁快速返回到极限位置的时候，将会对上横梁及其它部件产生最大的冲击振动，因此也将该时刻压砖机的模型作为有限元模态分析的对象。3.2.1模型参数把压砖机的PRO/E模型直接导入到ANSYS WORKBECH中。定义压砖机模型各零件的材料属性，把立柱、大螺母、柱塞、油缸套、下横梁、上横梁定义为结构钢，活动横梁定义为铸铁。由于在ANSYS WORKBECH中，只要定义好了属于那种材料之后，材料的泊松比、弹性模量、密度的材料属性就自动设置好了，材料属性如表3.1。表3.1压砖机的材料属性材料名称泊松比弹性模量(MPa)密度（kg/mm3）结构钢0.321057.8510-6铸铁0.281.11057.210-63.2.2网格划分模态分析中，把模型划为较粗的网格，对于模型的固有频率和振型的计算精度影响不大，但是可以减少计算时间。所以在压砖机模态分析中，把模型的网格划得较粗。由于压砖机是一种框架结构，用三维实体单元来描述该结构，更能反映机架的实际情况。在ANSYS WORKBENCH软件种，三维实体单元有两种：六面体单元和四面体单元。由于六面体单元在划分是要求结构比较规则，而对于压砖机这类较复杂的结构，对其进行六面体网格的自动划分十分困难，而用四面体单元分析三维结构，单元划分比较灵活，可以逼进复杂的几何形状，因此，对机架的网格划分采用四面体单元。因此使用系统默认（四面体单元）的网格划分方式对压砖机模型进行网格划分，得到：137377个节点，91568个单元，其模型网格图，如图3.1所示。图3.1 网格3.2.3约束在实际情况中，压砖机的下横梁通过螺栓固定于地基之上。因此，对压砖机模型的下横梁的下底面添加固定约束，如图3.2所示。图3.2 约束3.2.4压砖机模态分析结果评价压砖机自由振动前六阶振型云图，如图3.3，3.4，3.5，3.6，3.7，3.8所示。压砖机模态分析结果，总结如表3.2 。图3.3一阶模态云图 图3.4 二阶模态云图图3.5 压砖机三阶模态 图3.6 压砖机四阶模态图3.7五阶模态云图 图3.8六阶模态云图表3.2 压砖机模态分析结果总结模态阶数自由振动频率（Hz）最大振幅（）模态振型描述117.020.349压砖机上部前后弯曲振动237.70.247压砖机上部左右弯曲振动352.190.364压砖机左右两端扭转振动490.10.419上横梁的上部和活动横梁下部绕其中部作前后转动振动5114.70.211上横梁和活动横梁的中部上下伸缩震动6148.540.450上横梁和活动横梁绕它们中部左右转动振动一阶振型（f1= 17.02 Hz ）为压砖机上部前后弯曲振动，该振动将会增大压砖机四个拐角处的弯矩，从而增加其应力，同时加剧导套的磨损。二阶振型（f2=37.7Hz ）为压砖机上部左右弯曲振动，该振动将增大上下横梁的弯矩，同时加剧导套的磨损，给上模及液压缸的使用寿命带来不利影响。三阶振型（f352.197Hz ）为压砖机左右两端受扭转振动，该振动将会增大上下横梁的扭矩，同时加剧立柱的磨损。四阶振型（f4=90.1Hz）为上横梁的上部和活动横梁下部绕其中部作前后转动振动，该振动将会增大四个拐角处的弯矩，同时，加剧压砖机导向部分的磨损。五阶振型 （f5=114.7Hz ）为上横梁和活动横梁上下伸缩振动，该振动也将增大机架四个拐角的应力，同时，对与压砖机上横梁相连接的上模以及与下横梁相联接的液压缸、立柱、大螺母使用寿命产生不利影响。六阶振型（f6148.547Hz ）为上横梁和活动横梁绕它们中部左右转动振动，该振动将会增大上横梁的弯矩和四个拐角的应力同时加剧立柱部分的磨损。从以上的振型分析可以看到，压砖机不仅有前后，上下及左右方向的弯曲振动，而且有扭转振动，这些振动将影响压砖机的强度和刚度，加重压砖机立柱的磨损，影响的陶瓷砖精度及模具和液压缸的使用寿命。因此，在设计过程中应适当增加局部刚度和阻尼来抑制这些振动的影响。对压砖机前六阶模态分析可得：（1）压砖机整体的刚度和质量分布较为均衡，无明显的薄弱部位和过剩部位，这有利于压砖机的动力性能；（2）压砖机的最小固有频率为17.02Hz 。（而机架的冲击频率最大为0.3Hz 即18次每分钟），远远大于其冲击频率，因而该压砖机不会发生共振综上所述， WL1700全自动液压压砖机具有合理的动态参数，符合设计要求。3.3 本章小结 本章通过对压砖机的模态分析，可知压砖机的模态特性主要取决于：（1）压砖机结构的质量分布；（2）压砖机结构的刚度及其所受的约束情况。第4章 机架的静力学分析压砖机的机架在传统的设计中，都显得笨重，而且在一些零件的关键部位出现强度不够。本章通过对机架的静力学分析，得到机架的全场应力、应变、位移的情况，并进行强度较核，找出设计的缺陷。4.1 结构静力学分析概述结构静力学分析是用于分析稳态载荷或固定不变的惯性载荷（如重力和离心力）引起的系统或部件的应力、应变、变形位移等。结构静力学分析包括线性和非线性分析11。结构静力学分析的方程：Kx=F （4.1）结构静力学分析就是按方程4.1平衡关系确定计算的，其中K为刚度系数矩阵，x为位移矢量，F为力矢量11。4.2机架静力学分析4.2.1 模型参数在机架的结构静力学分析时，认为活动横梁与立柱之间的运动不应有干涉，因此去掉活动横梁及与活动横梁连接的柱塞和模具。所以压砖机整机PRO/E模型进行相应的简化之后，直接导入到ANSYS WORKBENCH中。定义机架的材料属性，定义为上横梁、下横梁、立柱、油缸套、大螺母为结构钢。材料参数如表4.1。表4.1压砖机的材料属性材料名称泊松比弹性模量(MPa)密度（kg/mm3）结构钢0.321057.8510-64.2.2 网格划分由于机架是一种框架结构，用三维实体单元来描述机架结构，更能放映机架的实际情况。在ANSYS WORKBENCH软件种，三维实体单元有两种：六面体单元和四面体单元。由于六面体单元在划分是要求结构比较规则，而对于机架这类较复杂的结构，对其进行六面体网格的自动划分十分困难，而用四面体单元分析三维结构，单元划分比较灵活，可以逼进复杂的几何形状，因此，对机架的网格划分采用四面体单元。并且对可能出现危险部位进行网格细化，以至分析的结果能较为精确。划分网格的结果是模型中共有213752个节点，140050个单元，机架模型网格如图4.1所示。 图4.1 机架模型网格图4.2.3 约束和载荷根据实际情况，对下横梁的下表面施加固定约束。从WL1700全自动液压压砖机的主要技术参数表2.1，可知主油缸的最大工作压力为33.8MPa，作用在下横梁上表面的作用力F由工程压制力和活动等零部件的重力产生，作用力的大小为17701301N，它与下模具接触的面积为1041253mm2,所以作用于下横梁上表面的压力P的大小如下：P=F/S=177000130/1041253=17MPa (4.2)上横梁主油缸内的载荷按均布载荷作用在油缸内壁；下横梁与下模具接触表面的载荷也按均布载荷考虑1。施加约束和载荷的结果，如图4.2所示。如图4.2施加压机架模型的约束和载荷4.2.3机架静力学分析结果评价 机架静力学分析云图，如图4.3，4.4，4.5所示。分析结果总结如表4.2。 图4.3机架的应力云图 图4.4 机架的应变云图图4.5 机架位移图表4.2机架静力学分析分析结果总结机架的位置应力（MPa）应变（10-2mm/mm）位移(mm)上横梁油缸顶角处200.370.170立柱与上、下横梁的拐角处224.1530. 2600.930上横梁中部1.753上横梁使用的ZG310-570,立柱使用的材料为45钢，它们的材料主要力学性能参数如表4.3表4.3 ZG310-570和45钢的主要力学性能材料名称屈服强度（s/MPa）抗拉强度（b/ MPa）ZG310-57031050045钢353598上横梁的安全系数：n1=310/200.153=1.549,立柱的安全系数n2=353/224.153=1.575从参考文献19中查得，铸钢和45钢的完全系数都为1.21.5,所以它们的设计满足静力强度的设计要求。 4.3 本章小结最大的应力存在于上横梁、底座和立柱之间的交接的拐角处,是这些部位存在附加弯矩，故产生较大的应力；在油缸内上部角处产生较大的应力，是在高压作用下，产生应力集中现象。上横梁中部的位移变形横立柱位移变形位都会影响活动横梁的运动精度，还使得立柱与活动活动横梁之间产生相互的作用力，同时增加它们之间的磨损。在上横梁油缸顶角处的应力较大，通过适当的增大过渡圆角值，使应力值减小。通过静力强度较核发现，立柱和上横梁的安全系数要高于许用的安全系数，这可能是在传统的设计中，太保守的设计所造成的，这样的设计导致压砖机比较笨重，而且提高了生产成本。第5章 活动横梁的有限元分析与优化活动横梁是压砖机中压制陶瓷砖的关键的运动零件，它受到柱塞对它的冲击，模具对它的挤压力，并且长期工作在脉动载荷下，常发生断裂的严重事故。所以对活动横梁采用比较全面的有限元分析与优化。5.1静力学分析在活动横梁存在最大公称压制力的情况时，对活动横梁的进行静力学分析。得到活动横梁的应力、应变、位移，并且对其强度进行较核。5.1.1 模型参数把用活动横梁PRO/E模型直接导入到ANSYS WORKBECH中，定义活动横梁材料为铸铁，即材料属性如表5.1。表5.1活动横梁的材料属性材料名称泊松比弹性模量(MPa)密度（kg/mm3）铸铁0.281.11057.210-65.1.2 划分网格选用四面体单元对对活动横梁进行网格划分，并且划为较细的网格，划分网格的结果是模型中共有22580个节点，14650个单元，活动横梁模型网格图，如图5.1所示。图5.1活动横梁模型网格图5.1.3 约束和载荷活动横梁的四个力柱孔分别受到四根立柱的仅有压缩约束。仅有压缩约束是仅仅限制这个表面在约束的法向正方向的移动。在任何给定的表面可以施加法向仅有压缩的约束。解释这种约束的一种方法就是将它想象为一个“刚性”结构，它与选择的表面有相同的形状。注意到这些接触（压缩）面事先不知道。可以在一个圆柱面上模拟“扣牢的圆柱约束”。如图5.2所示，显示出了没有变形的圆柱的轮廓。有压缩力的表面阻止原始圆柱变形，而可伸长的表面自由变形。这种约束需要一个迭代来求解。由于事先不知道压缩面的行为，所以需要利用迭代求解器来判断哪个表面显示的是压缩行为。图5.2 仅有压缩约束从WL1700全自动液压压砖机的主要技术参数可知，活动横梁的上部与柱塞接触的表面受到的最大公称压制力为17000000N，而且接触的上表面面积385456mm2，所以压力为44.104MPa,活动横梁的下表面受到的来自坯料的作用力等于公称压制力与活动横梁的重力之和，其大小为17042359N，活动横梁与下模接触的面积为1041253 mm2,所以下表面的挤压力为16.356MPa。活动横梁模型约束和载荷，如图5.3所示。图5.3 活动横梁模型约束和载荷5.1.4 活动横梁静力学分析结果活动横梁静力学分析的应力、应变、位移结果，如图5.4，5.5，5.6所示，分析结果总结，如表5.2。图5.4 活动横梁的应力云图图5.5 活动横梁的应变云图图5.6 活动横梁位移图表5.2 活动横梁静力学分析结果出现的位置最大应力（MPa）最大应变（10-2mm/mm)最大位移(mm)立柱孔和中间孔处122.0760.1110.476活动横梁选用铸铁为HT350，HT350材料的主要力学性能如表5.3。对活动横梁进行静力学的强度较核：n=350/122.076=2.87。从参考文献19中查得，铸铁零件要求的安全系数为1.52.5。所以能满足活动横梁静力学强度设计的要求。表5.3 HT350材料的主要力学性能材料名称抗拉强度（MPa）抗压强度（MPa）疲劳强度（MPa）抗剪强度（MPa）HT35035011001300140180350500活动横梁的上表面的孔处和四个立柱孔处，出现最大应力集中现象。立柱孔处的应力集中主要是弯矩引起；中间孔处主要是柱塞的压力所引起。5.2 动态时域分析5.2.1 动态时域分析概述动态时域分析用于研究随时间变化的载荷作用下结构应力、应变、位移等的响应技术16。求解结构动态时域分析的方法有振型叠加法和直接积分法。前者用模态振型矩阵作为变形矩阵，将结构在物理空间上耦合的动力学方程转化为在模态空间上的一组非耦合动力学方程，然后利用单自由度振动方程的动态响应分析方法求解，最后叠加到结构在物理空间的运动响应。该法应先求出结构的固有频率和固有振型，因而适应于阻尼矩阵C可以对角化和不太复杂的激振力的情况。而后者则对动力方程直接用数值积分法求解，将时间离散化，在一个时间间隔内对位移、速度和加速度的关系做某种假设，将运动方程转化成该时刻点上位移矢量的代数方程，最后用递推公式求解。该法不须先求出结构的固有频率和振型，因而适合动载荷复杂情况10。动态时域分析是研究结构在外部激振力的作用下的各种响应，包括应力、位移等的响应。结构在外部激振的响应导致系统内部产生动态应力和动态位移，从而影响产品的使用寿命和工作性能。采用Pro/MECHANICA 的分析软件，对活动横梁进行动态时域分析，分析随时间变化的载荷下，不同时刻模型应力、应变及位移。此外，创建一个应力变化的测量值。5.2.2 模型参数在使用Pro/MECHANICA 有限元分析元件中，自定义活动材料属性，如表5.4，定义好材料属性之后，分配于整个活动横梁模型上。表5.4.活动横梁的材料属性材料名称泊松比弹性模量(MPa)密度（kg/mm3）铸铁0.281.11057.210-65.2.3 划分网格划分网格与在ANSYS WORKBENCH中类似，使用四面体实体单元进行网格划分，对于活动横梁应力集中的部位采用较细的网格。划分的网格的结果是模型中有26834个节点，16935个单元，其网格图，如图5.7所示。如图5.7活动横梁网格图5.2.4 约束和载荷添加约束条件，约束活动横梁两侧的四个立柱孔，x、z方向作为固定，y方向保持自由，x、y、z的转动，也作为固定。定义随时间变化的载荷，通过表格的方式绘制。活动横梁的上表面与柱塞接触的面积S1=385456mm2,作用力为F1；下表面受到挤压力的面积S2=1041253mm2,作用力为F2。所以，上表面的压力P1=F1/S2，下表面的压力P2=F2/S2,它们表面作用力随时间变化，如表7.1。根据表格，在软件Pro/MECHANICA 的分析中绘制的函数图，如下图5.8和图5.9所示。活动横梁模型在动态时域分析约束和载荷如图5.10所示。表5.5 活动横梁随时间变化的载荷受压力的表面活动横梁压制坯料的时间(s)00.10.40.50.80.91.3上表面S1(MPa)011.02611.02625.72725.72744.10444.104下表面S2(MPa)0.0444.1224.12209.5599.55916.35616.356图5.8 上表面随时间变化载荷的函数图图5.9下表面随时间变化载荷的函数图图5.10活动横梁模型在动态时域分析约束和载荷5.2.5 活动横梁的动态时域分析结果评价依据前面的静态应力分析的结果，在最大应力处建立测量点Measure1,即选择活动横梁中间孔的边缘一点作为测量点。动态时域分析之后，得到测量点在此时域上的应力变化如图5.11所示，其最大的应力值为131.5MPa，而静态分析的应力值122.076MPa,在动载荷下的应力值比最大的静载荷的应力值增大了9%。动态时域分析的应力云图，如图5.12，5.13，5.14所示。图5.11 测量点的应力变化图5.12 活动横梁的应力图5.13活动横梁的应变图5.14活动横梁的位移经分析发现，活动横梁的中间孔和两侧的四个立柱孔，存在显著的应力集中的现象。在活动横梁铸造时应严格铸造工艺，减少缺陷，不能出现微裂纹或微孔隙。四个立柱孔处应力集中最为显著，因此，在活动横梁的四个立柱孔处装上结构钢的导套，来缓解应力集中导致立柱孔处的破坏，并且提高活动横梁的导向性。5.3 疲劳分析5.3.1 疲劳分析概述疲劳是指结构在底于静态强度极限的载荷重复作用下出现疲劳断的现象14。疲劳破坏的主要因素有载荷的循环次数；每一个循环的应力副值和平均应力；局部应力集中14。疲劳破坏的过程是:零部件在循环载荷作用下，在局部的最高应力处，最弱的及应力最大的晶粒上形成微裂纹，然后发展成宏观裂纹，裂纹继续扩展，最终导致疲劳断裂。所以，疲劳破坏经历了裂纹形成、裂纹扩展和瞬断三个阶段。疲劳分析原理是采用了简化的弹性假设和Miner累积疲劳准则。Miner累积疲劳准则是材料在重复外载荷作用下，其内部的微缺陷和空穴不断地产生与发展，从宏观上表现为材料损伤的累积，当材料的损伤累积达到临界值时则产生失效15。5.3.2 活动横梁疲劳分析活动横梁在运动过程中，各点应力变化符合脉动应力变化规律。脉动循环变应力曲线如图5.15所示，纵坐标值1表示活动横梁的最大应力值。活动横梁的疲劳分析是建立在静力学分析的基础上的，因此，活动横梁的疲劳分析处理和其线性静力分析很相似,因此这里不做详细的介绍。不通在于要定义材料的S-N曲线，只是在求解时，选用疲劳求器，即可。图5.15 脉动循环变应力一般的SN曲线有这样的共同特点，即当破坏循环次数小于103，疲劳极限x接近于b，而且在1到103之间，变化不大；而且当循环次数大于106时，疲劳极限x趋近于-1，而在103106之间，存在着关系式，m和c为常数，对等式两边取对数后得： （5.1）因此在双对数坐标的SN曲线图和之间存在线形关系，只要知道103和106两点的值，这个直线就确定了。在N103时，近似取x=0.9b，而在N106时，近似取x=-1，由于活动横梁的材料为铸铁，型号为HT350，它的b=350MPa，-1=180MPa。所以，当N=103时，x=315 MPa；当N=106时，x=180 MPa，所以在铸铁材料HT350的SN曲线，如图5.16所示。图5.16 SN曲线活动横梁的设计寿命为循环109次，即为高周疲劳，受应力幅控制，Mean Stress Theory选用工程上较为常用的Goodman曲线，如图5.17所示。图5.17 Goodman曲线5.3.3 活动横梁疲劳分析结果（1）活动横梁的疲劳寿命分析结果，如下图5.18所示。从分析结果中得到活动横梁的疲劳寿命最小值也能达了109，到达了设计的寿命要求，因此符合活动横梁寿命的设计要求。图5.18活动横梁寿命（2）安全系数是指计算出来的各个节点寿命与指定的设计寿命之比。分析结果如图5.19所示，安全系数为1.741，从参考文献19中查得，铸铁零件要求的安全系数为1.52.5。所以活动横梁的设计符合安全可靠的要求。 图5.19 活动横梁安全系数（3）活动横梁的疲劳灵敏性是指在静态分析载荷在0与150%之间作脉动循环时，活动横梁的疲劳寿命变化情况。活动横梁载荷的脉动变化的实际变化为0到100%之间，所以通过分析结果图5.20，可知，活动横梁的疲劳灵敏性非常好。图5.20 活动横梁的疲劳灵敏性通过对活动横梁的疲劳计算，活动横梁的疲劳寿命、安全系数，都能达到活动横梁疲劳设计的要求；活动横梁的疲劳灵敏性在脉动载荷变化的情况下，灵敏性非常好。 5.4 拓扑优化5.4.1 拓扑优化概述拓扑优化是形状优化，有时也称为外型优化。外形优化的目标是寻找承受单载荷或多载荷的物体的最佳材料分配方案。这种方案在拓扑优化中表现为“最大刚度优化”设计。与传统的优化设计不同的是，拓扑优化不需要给出参数和优化变量。目标函数、状态变量和设计变量都是预定义好的。用户只需要给出结构的参数（如材料特性、模型、载荷等）和需要省去的材料百分比。拓扑优化只适用于实体。形状优化的目标函数是在满足结构约束的情况下减少结构的变形能。减少结构的变形能相当于提高结构的刚度。这种技术通过使用设计变量给每个单元赋予内部伪密度来实现16。5.4.2 ANSYNS WORKBENCH拓扑优化“Shape Finder”是一个优化问题,其结构能量在减少结构体积的基础上的最小化。它尽量得到关于体积比率的最好刚度；尽可能的找寻可以在对整体结构的强度不产生负面影响的可去除的面积。由于“Shape Finder”在载荷和约束的基础上尽量减少体积和增大刚度， 因此载荷和约束十分重要并会影响结果。“Shape Finder”一般在载荷施加的地方和约束反作用于载荷的地方保持材料不变。缩减材料不能太多，否则会导致行架式结构。“Shape Finder”基于一个单独的静力结构环境下，“Shape Finder”不能用于复合环境下。5.4.3 活动横梁的拓扑优化活动横梁的拓扑优化处理和其线性静力分析很相似,因此这里不是所有的步骤都详细介绍。只是在求解时，选用拓扑优化器即可。优化结果，如图5.20，5.21，5.22，5.23所示，图中橙色代表可去除的材料，褐色的是边缘。活动横梁的初始和最终结构的质量边缘质量详细比较如表5.1。图5.20去除10%材料 图5.21去除15%材料 图5.22 去除20%材料图 图 5.23去除25%材料表5.2初始和最终结构的质量的边缘质量详细比较优化质量百分比初始结构质量（）最终结构质量（）边缘质量（）10%45286418316720415%45286401263805220%45286384168443825%452863674347100活动横梁的上表面中间处受到柱塞对它的压力，而下表面受到挤压力，因此活动横梁上部两侧的斜角部分所受到的影响很小；四个力柱孔受到立柱的仅压缩约束，但其下边缘处受到立柱的约束影响很小,因此，活动横梁的上部两侧的斜角部分和立柱下边缘处，在合理的刚度情况下，可以去除这些部分的材料，达到优化活动横梁质量的目的。5.5 本章小结（1）通过对活动横梁的静力学分析，找出活动横梁应力集中的最大部位存在于活动横梁的四个立柱孔处和中间的孔处。（2）通过活动横梁的动态时域分析，找出活动横梁应力集中处的测量点的应力随时间变化的情况。（3）通过活动横梁的疲劳分析，分析出疲劳寿命值、疲劳不会破坏的可能性、 安全系数及疲劳灵敏性，从分析中分析结果中，可以发现活动横梁满足疲劳强度设计得要求。（4）通过对活动横梁的拓扑优化，发现活动横梁能够进行材料去除的部位，使得活动横梁在满足强度要求的情况下达到质量优化的目的。结 论本文基于有限元分析软件ANSYS WORKBENCH和PRO/MECHANICA平台，对WL1700全自动液压压砖机关键零、部件进行有限元结构分析。进行压砖机整机模态分析，得到前六阶固有频率及相应的固有振型等动态特性参数，得到以下两点结论：（1） 压砖机的最小频率为17.0159Hz ，而压砖机的冲击频率为0.3Hz，压砖机的固有频率远远大于其冲击频率，因而该压砖机不会发生共振；（2）压砖机不仅有前后，上下及左右方向的弯曲振动，而且有扭转振动，这些振动将影响压砖机的强度和刚度，具体表现为上横梁的弯曲断裂以及加重压砖机导向部分即立柱的磨损，也影响陶瓷砖的制造精度，因此，在设计过程中应增加局部刚度和阻尼来抑制这些振动的影响。进行机架和立柱的静力学分析之后，得到以下结论：（1）通过对机架的静力学分析，找到机架应力集中的位置为上横梁、底座和立柱之间的交接处的立柱处，其大小为224.153MPa，较大的出现在上横油缸顶角处，其大小为200.37MPa,强度设计满足要求，但是有些保守，造成机架较为笨重。上横梁油缸部位的变形位移最大，四根立柱也产生了较大的变形位移，向机架中间弯曲，这将影响到活动横梁的导向性。进行活动横梁比较全面的有限元分析和优化得到以下结论：（1）活动横梁应力集中的最大部位存在于活动横梁的四个立柱孔处和中间的孔处。（2）通过活动横梁的动态时域分析，找出活动横梁应力集中处的测量点的应力随时间变化的情况。提取第一步断点分析结果，应力值为131.5MPa，比静力学中的分析结果增大了9%。 （3）活动横梁的疲劳强度满足活动横梁疲劳寿命的设计要求。（4）活动横梁的上部两侧的斜角部分和立柱孔边缘处，在合理的刚度情况下，可以去除这些部分的材料，优化活动横梁质量。本文通过对全自动液压压砖机的关键零、部件进行有限元结构分析，可以得到关键零、部件精确的静动态力学参数和拓扑优化，这对于了解现有的压砖机结构特性并进而改善其性能将具有重要的意义，并且为全自动液压压砖机的整体改进设计，提供理论依据。参 考 文 献1 张柏清自动液压压砖机M南昌：江西科学出版社,200152 张柏清，戴哲敏，高文杰KD3200型全自动液压压砖机的特点及液压系统分析J 中国陶瓷工业 ，2000（2）：3135 3 彭娅清，唐华平，万海如，张红字3200吨全自动液压压砖机的有限元强度分析J 制造信息化，2005（8）：49514 张柏清，胡伟文 基于有限元方法的液压压砖机活动横梁结构参数研究J 陶瓷学报，2006（1）：29335 杨金利300kN压力机上横梁的有限元分析J 机械设计，2006（4）：24256 白葳，喻海良通用有限元分析ANSYS8.0基础教程M北京：清华大学出版社，200557 张胜民基于有限元软件ANSYS7.0的结构分析M北京：清华大学出版社，200318 张朝晖ANSYS8.0结构分析及实例解析北京：机械工业出版社M，200539 王勖成有限单元法M 北京：清华大学出版社，2002510 刘国庆，杨庆东ANSYS工程应用教程北京：.中国铁道出版社M，200311 李重伟有限元分析方法综述J 天津建设科技20O6增刊（13）12 姜勇，张波ANSYS 7.0实例精解M 北京：清华大学出版社，2003313 龙志飞，芩松有限元法新论：原理、程序、进展M 北京：中国水利水电出版社，2001914 商跃进有限元原理与ANSYS应用指南北京：清华大学出版社，2005615 龚曙光ANSYS基础应用及范例解析北京：机械工业出版社，2003116 强锋科技，李黎明ANSYS有限元使用教程北京：清华大学出版社，2006317 二代龙震工作室Pro/Mechanism/MECHANICA wildfire 2.0 机构/运动/结构/ 热力分析M北京：电子工业出版社，2006118 Zissimos P. Mourelatos. A crankshaft system model for structural dynamic anslysis of internal combustion enginesJ Cimputers and Structures 2001 79(6): 2009202719 徐锦康机械设计M北京：高等教育出版社，2004420 机械工程师手册第二版编辑委员会机械工程师手册M北京：机械工业出版社，2000621 俞新陆液压机现代设计理论M北京：机械工业出版社，19871