2D型活塞压缩机关键零部件的有限元分析论文

大庆石油学院本科生毕业论文 I 目 录 第 1 章 绪 论 . 1 1.1 课题的目的 . 1 1.2 背景 . 1 1.3 课题内容及 意义 . 3 第 2 章 活塞式压缩机简介 . 5 2.1 活塞式压缩机概述以及相关参数 . 5 2.2 活塞式压缩机分类及结构 . 7 2.3 活塞式压缩机的工作原理 . 8 2.4 活塞式压缩机在石油化工工业中的应用 . 9 第 3 章 曲轴有限元分析 . 11 3.1 对曲轴进行有限元分析的重要性 . 11 3.2 曲轴相关力学计算 . 12 3.3 曲轴的有限元分析 . 19 3.4 小结 . 23 第 4 章 连杆有限元分析 . 24 4.1 对连杆进行有限元分析的重要性 . 24 4.2 连杆相关力学计算 . 24 4.3 连杆的有限元分析 . 32 4.4 小结 . 35 第 5 章 十字头有限元分析 . 37 5.1 对十字头进行有限元分析的重要性 . 37 5.2 十字头相关力学计算 . 37 5.3 十字头的有限元分析 . 42 5.4 小结 . 46 结 论 . 48 大庆石油学院本科生毕业论文 II 参考文献 . 49 致 谢 . 51 大庆石油学院本科生毕业论文 1 第 1章 绪 论 1.1 课题的 目的 本次毕业设计的题目 是“ 2D型活塞压缩机关键零部件的有限元分析”。 课题的目 的 是 对活塞式压缩机的结构、功能以及 应用 等方面进行一些基础的了解。重点在于 运用 ANSYS软件 分析 2D型活塞式压缩机的关键零部件 ，并 通过对 ANSYS软件的学习， 掌握 该软件 完成 有限元分析的方法与步 骤， 并运用软件 对毕业设计中 具体 给定的压缩机实例进行分析，通过 完成整个 分析过程 使自己达到 对 ANSYS软件的熟练操作 ，证明 ANSYS软件可以简洁快速的完成相关课题的分析。 解决以上问题就是本文的中心目的。 1.2 背景 活塞式压缩机 是容积型压缩机中应用最广泛的一种 ，也是石油化工行业中比较重要的机械之一，各项工业流程几乎都可以看到它的身影。不论是物料的运输、加工及成型都少不了它的参与，因此对压缩机进行力学分析是一项必不可少的关键性工作。活塞式压缩机的易损件很多，其形状，受力状态又不是很规律，因此用常规的计算方法很难展开 工作，而且由于理论工况与实际工况存在着不小的差距，所以计算的结果的误差也很难让人接受。 计算的结果有时也很难让人接受，从而无法解决实际问题。 这时，有限元方法的出现让人眼前一亮。 它把分析的物体划分成有限多个 单元，再分别对这些单元进行分析，分析的结果大大减小了误差，导入相应软件进行分析也大大减轻了繁琐的计算工作量。 有限元法就是这样发展起来的。 有限元方法是随着电子计算机的发展而迅速发展起来的一种在计算数学、计算力学和计算工程科学领域最有效的现代计算方法。经过 40 多年的发展已经使各种不同的有限元方法形态相当丰富，理论 基础相当完善，并且开发了一批通用和专用有限元软件，如 ANSYS、 MSCNASTRAN、 MSCMARC、 ABAQUS和 ALGOR 等。有限元方法正被日益广泛的用于计算机辅助工程分析（ CAE） ,并与计算机辅助设计（ CAD）和计算机辅助制造（ CAM）等集合成一体为计算机整体制造 (CIM)。随着新产品设计与产品开发的需求日益增加，产品质量与成本控制越来越严格，有限元计算技术显示出极大的优越性。它可以帮助工业界降低成本，改进质量，缩短产品设计与开发的周期 。在众多可用的通用和专用有限元软件中，大庆石油学院本科生毕业论文 2 ANSYS 是最为通 用和有效的商用有限元软件之一，用户也是最广。 ANSYS 软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。由世界上最大的有限元分析软件公司之一的美国 ANSYS 开发 ,它能与多数 CAD 软件接口 ,实现数据的共享和交换。 作为一个流行的工程设计软件，具有三维立体几何构造，计算模拟，虚拟成形试样产品等功能，已被工业界公认为十分有效的计算和分析工具。它是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件 。 该项目在我国发展正处于初步阶段，相关资料并不充分 ， 在我国 该软件 大部分应用 于对静设备进行 分析。 从单纯的结构力学计算发展到求解许多物理场问题有限元分析方法最早是从结构化矩阵分析发展而来，逐步推广到板、壳和实体等连续体固体力学分析，实践证明这是一种非常有效的数值分析方法。而且从理论上也已经证明，只要用于离散求解对象的单元足够小，所得的解就可足够逼近于精确值。所以近年来有限元方法已发展到流体力学、温度场、电传导、磁场、渗流和声场等问题的求解计算，最近又发展到求解几个交叉学科的问题 。由求解线性工程问题进展到分析非线性问题随着科学技术的发展，线性理论已经远远不能满足设计的要求。例如建 筑行业中的高层建筑和大跨度悬索桥的出现，就要求考虑结构的大位移和大应变等几何非线性问题；航天和动力工程的高温部件存在热变形和热应力，也要考虑材料的非线性问题；诸如塑料、橡胶和复合材料等各种新材料的出现，仅靠线性计算理论就不足以解决遇到的问题，只有采用非线性有限元算法才能解决。众所周知，非线性的数值计算是很复杂的，它涉及到很多专门的数学问题和运算技巧，很难为一般工程技术人员所掌握。为此近年来国外一些公司花费了大量的人力和投资开发诸如 MARC、 ABQUS 和 ADINA 等专长于求解非线性问题的有限元分析软件，并广泛 应用于工程实践 。 1965 年“有限元”这个名词第一次出现 ， 到今天有限元在工程上得到广泛应用 ， 经历了 30多年的发展历史 ,理论和算法都已经日趋完善。现在从汽车到航天飞机几乎所有的设计制造都已离不开有限元分析计算 ， 其在机械制造 、 材料加工、航空航天、汽车、土木建筑、电子电器 、 国防军工 、 船舶 、 铁道 、 石化 、 能源 、科学研究等各个领域的广泛应用已使设计水平发生了质的飞跃 。 有限元方法是数值计算中的一种离散化方法。用数学术语来说 ,就是从变分原理出发 ， 通过分区插值 ， 把二次泛函 （ 能量积分 ） 的极值问题转化为一组多元线性代数方程来求解 ； 从物理和几何概念来说 ， 有限元方法是结构分析的一种计算方法 ， 是矩阵方法在结构力学和弹性力学等领域的发展和应用 ， 其基本思想是将弹性体划分为有限个单元 ， 对每个单元 ， 用有限个参数来描述它的力学特性 ， 而整个连续弹性体的力学特性可认为是这些小单元力学特性的总和 ， 从而建立起连续体的力平衡关系。这种方法常用于复杂弹性振动系统 ， 其求解方式现在有多种可利用的商业化软件 ， 这里我们应用 ANSYS软件系统 。 有限元分析步骤如下： 大庆石油学院本科生毕业论文 3 第一步 ， 问题及求解域定义 ： 根据实际问题近似确定求解域的物理性质和几何区域。 第二步 ， 求解域离散化 ： 将求解域近似为具有不同有限大小和形状且彼此相连的有限个单元组成的离散域 ， 习惯上称为有限元网络划分。显然单元越小 （ 网络越细 ） 则离散域的近似程度越好 ， 计算结果也越精确 ， 但计算量及误差都将增大 ， 因此求解域的离散化是有限元法的核心技术之一。 第三步 ， 确定状态变量及控制方法 ： 一个具体的物理问题通常可以用一组包含问题状态变量边界条件的微分方程式表示 ， 为适合有限元求解 ， 通常将微分方程化为等价的泛函形式。 第四步 ， 单元推导 ： 对单元构造一个适合的近似解 ， 即推导有限单元的列式 ，其中包括选择合理的单元坐标系 ， 建立 单元试函数 ， 以某种方法给出单元各状态变量的离散关系 ， 从而形成单元矩阵。 第五步 ， 总装求解 ： 将单元总装形成离散域的总矩阵方程 （ 联合方程组 ），反映对近似求解域的离散域的要求 ， 即单元函数的连续性要满足一定的连续条件。总装是在相邻单元结点进行 ， 状态变量及其导数连续性建立在结点处。 第六步 ， 联立方程组求解和结果解释 ： 有限元法最终导致联立方程组 ， 联立方程组的求解可用直接法、选代法和随机法。求解结果是单元结点处状态变量的近似值。对于计算结果的质量 ， 将通过与设计准则提供的允许值比较来评价 ， 并确定是否需要重复计算 。 利 用 ANSYS软件 进行模态分析的过程分为 4个步骤 ： （ 1） 建模 ；（ 2） 加载求解 ； （ 3） 扩展模态 ； （ 4） 观察结果。 1.3 课题内容 及 意义 压缩机的性能测试涉及到热动力学、流体力学、电子学、计算机科学和现代控制理论。压缩机性能测试对于研究和改善压缩机的热力性能及可靠性具有重要的意义。随着计算机技术及自动控制技术的迅猛发展 ， 工业计算机越来越广泛地被运用到各种控制系统中。利用计算机友好的人机界面 ， 实现人机对话和监控系统 ， 并可借助计算机强大的运算和管理能力 ， 对数据进行处理和保存 ， 实现压缩机性能测试。其中软件的设计是最为重要的 ， 它直接面向用户 ， 并且决定了测试系统的功能。 该软件可以将分析的物体进行有限元化，并通过设定参数和求解过程对其进行力学分析，让人可以直观、形象的了解所分析结构的受力情况。 课题的内容为 了解活塞式压缩机在石油化工工业中的 实际 应用， 同时 了解 2D型活塞式压缩机的结构及工作原理，重点完成 2D型活塞式压缩机关键零部件的有限元计算和分析 （ 应用 ANSYS软件 对压缩机重要 零部件 进行应力分析， 通过 对 相关软件的学习对所研究的零部件进行 建模、实参数的选择及确定、力学状态选择、大庆石油学院本科生毕业论文 4 载荷虚拟加载、受力加载分析，之后根据所生成的图像 观察其应力 分布情况、最危险截面位置以及零部件变形情况、位移大小等，根据这些分析的结果来确定其结构的 合理程度 并进行相关校核计算，以及 讨论 对于分析结果的解决方法），并将分析结果用图表、图像等方法具体、生动的展现出来，从而完成有限元分析过程。 论文主体 内容 分为 五 章， 第一章 绪论 是关于 应用 有限元分析 方法 的 简单介绍 ，第二 章是关于活塞式压缩机简介的内容，第 三 、四、五 章是 关于活塞 式 压缩机的曲轴、连杆 以及十字头 的相关力学计算以及有限元分析。 该毕业设计是 在理论力学与材料力学的基础上，对研究对象（压缩机传动机构）在一定工况所确定的载荷作 用下其结构的强度、刚度及屈服稳定性情况有一定的了解，并以直观的图像、表格等 形式 表现出来。 根据这些图像以及理论计算的内容相比较从而解决工程问题。 在企业中可根据所计算得到的结构上各处位移和应力的分布结果结合相应的设计标准或规范，判断所设计的结构可靠性和经济性，从而避免某些结构或零部件由于过大的应力压曲失稳而损坏、并控制机械结构整体及其零部件的刚性性能， 能够对较复杂机械零件进行静力分析，并能基本掌握机械的非线性分析和瞬态分析，能够对复杂的机械零部件进行应力和变形的分析，对简单装配件进行接触分析 ANSYS 软 件的功能如下： 基于工程学的理论以及许多数值分析的理论技术 ； 可解决大部分工程上的问题 ； 使用相当有效的解题方法 ； 以用户为导向，容易定义问题 ； 完全由问题的定义推得结果 ； 有完整且蕴含高度技巧性的图形表示能力 ； 有完整的帮助文件和完整的例题，而且经过许多验证 ； 有完整的工程及技术人员继续开发新技术 ； 每年有完整的研讨会，并有完整的电子公告栏以帮助用户解决遇到的问 题 。 大庆石油学院本科生毕业论文 5 第 2章 活塞式压缩机简介 2.1 活塞式压缩机概述 以及相关参数 从世界范围内看压缩机的发展历程和概况。活塞式压缩机的发展历史悠久，具有丰富的设计、研究、制造和运行的 经验，至今在各个领域中依然被广泛采用、发展着。然而，也必须注意到，制冷压缩机的不断进步也反映在其种类的多样性方面，活塞式以外的各类压缩机机型，如离心式、螺杆式、滚动转子式和涡旋式等均被有效地开发和利用，并各具特色，这就为我们制冷工程的业内人士在机型的选择上提供了更多的可能性。在这样的背景之下，活塞式压缩机的使用范围必然受到一定影响而出现逐渐缩小的趋势，这一趋势在大冷量范围内表现得更为显著。在中小冷量范围内，实际上还是以活塞式压缩机为主 。 活塞式压缩机 属于容积式压缩机，适用于中小输气量，排气压力可从低压直至超高压，与其它类型压缩机相比，具有 如下 一系列特点 ： 首先，运动机构的尺寸确定后，工作腔的容积变化规律也就确定了，因此，机器转速的改变对工作腔容积变化规律不发生直接影响，故机器压力与流量关系不大，工作的稳定性较好。 其次，气体的吸入与排出是靠工作腔容积变化，与气体性质关系不大，故机器适应性强并容易达到较高的压力。 再次，机器的热效率较高。 最后，容积式机器结构较复杂，尤其是往复 式压缩机易于损坏的零件多。此外，气体吸入和排出是间歇的，容易引起气柱及管道的振动 。 活塞式压缩机 在 各种场合，特别是 在 中小制冷范围内，成为应用最广、生产批量最大 的 一种机型。 本文研究的对象是 2D40-51.5/4-32-BX型原料气压缩机 。 关于活塞压缩机的曲轴、连杆及十字头的相关力学计算以及有限元分析。 选取几个有代表意义的位置进行分析。运用有限元进行力学加载分析，根据计算生成的应力云图确定其危险截面，并与理论计算进行比较，判断其危险截面所在，并根据其力学分析结果针对确定的危险截面总结相应的应对措施来减 少其失效可能，使其运行工况相对安全稳定。 以下为 原料气压缩机性能参数及主要技术指标 。 大庆石油学院本科生毕业论文 6 表 2-1 原料气压缩机性能参数及主要技术指标 项 目 单 位 参数值及技术指标 设备位号 P-1201A/B 型号名称 2D40-51.5/4-32-BX 型原料气压缩机 介质名称 原料气 组 成 (体积百分数 %) H2O NH3 H2S N2 H2 CH4 0.197 0.024 60ppm 1.319 12.35 59.37 C2H6 C3H8 iC4H10 nC4H10 iC5H12 nC5H12 6.57 7.836 4.179 5.544 0.261 2.249 C5H10 C6H14 0.004 0.097 级 数 级 级 吸气压力 MPa(G) 0.4 1.15 排气压力 MPa(G) 1.15 3.2 吸气温度 45 40 排气温度 105 105 各级安全阀开启压力 MPa(G) 1.27 3.52 排 气 量 m3/h 12768 轴 功 率 KW 1060 转向 从非驱动端看压缩机为 逆时针 冷却水总进水压力 MPa(G) 0.3 活塞杆摩擦表面温度 1100 曲轴转速 r/min 333 填料充氮压力 MPa(G) 0.1 循环油量 L/min 100 氮气消耗量 Nm3/h 6 润滑油系统供油压力 MPa(G) 0.3 0.4 润滑油系统供油温度 45 主轴承温度 65 电机定子温度 135 供气缸水温度 40 50 供填料水温度 35 机身振动 g 2 填料冷却水过滤差压 MPa(G) PreprocessorElement TypeAdd/Edit/Delete 命令，弹出Element Type对话框，单击 Add按钮，在弹出的对话框中选择 Structural及其下的Solid，在其右边的滚动栏中选择 8node 45，单击 OK按钮。 执行 Main MenuPreprocessorMaterial Props Material Model命令，弹出 Define Material Model Behavior 对话框， 在右边 的 菜 单 中 ， 连 续 双 击StructuralLinearElasticIsotropic 后又 出现一 个 Linear isotropic Properties for Material Number 1对话框，在 EX后面的输入栏中输入 2.1E+5，在 PRXY后面的输入栏中输入 0.3，单击 OK按钮，接着继续在 Define Material Model Behavior对话框中连续双击 StructuralDensity，弹出密度定义对话框，在 DENS后面的输入栏中输入7.85E-9，单击 OK按钮，关闭该对话框 。 到此 曲轴模型的 单元类型的选择及材料属性的定义工作完毕 。 3.3.4 网格划分 图 3-3 曲轴网格划分 在 这 次 设 计 中 我 采 用 的 网 格 划 分 方 法 是 自 由 网 格 划 分 。 执 行 Utility MenuSelectEntities，在弹出的对话框中， 将第一项选为 Volumes，第二项选为 By Num/Pick，下面默认为 From Full，点击 Invert，之后单击 OK按钮。之后选取整个曲轴 实 体 ， 单 击 弹 出 对 话 框 的 OK 按钮。执行 Main 大庆石油学院本科生毕业论文 21 MenuPreprocessorMeshingMeshTool，在 Smart Size这项选为 5。之后执行 Main Menu PreprocessorMeshing VolumesFree，之后选取整个曲轴实体，单击弹出对话框的 Pick All按钮，等待其运行完毕，观察确认网格已经生成后，单击工具条上的 SAVE_DB按钮，保存网格模型，操作到此为止网格划分完毕 。 3.3.5 进入求解器加载并求解 执行 Main MenuSolutionAnalysis TypeNew Analysis命令，弹出 New Analysis对话框，选择分析类型为 Static，单击 OK按钮，关闭该对话框。执行 Main MenuSolutionAnalysis TypeSoln Controls命令，弹出 Solution Controls对话框。单击 Basic选项卡，设定求解时间 1，打开自动时间步长，子步数设为 40，最大子步数为 100。之后单击 Solution Controls对话框上的 Nonlinear选项卡，在该对话框中，将 Linear search项置为 ON，同时将 DOF solution predictor项设为 On for all substp，单击 OK按钮，关闭该对话框，到此求解参数设定完毕 。 确定约束位置 ， 曲轴与连杆是圆弧接触，因此认为力是作用在一条线上。 曲轴的受力情况可以简化为其两端的轴部受到约束，两个曲拐上受到连杆的载荷先进行约束的施加 。 第一 步 ， 执行 Utility MenuSelectEntities，在弹出的对话框中，将第一项选Lines，第二项选为 By Num/Pick，下面默认为 From Full，点击 Invert，之后单击 OK按钮。 第二 步 ， 选取 实体上 所加 约束 的线，单击弹出对话框的 OK按钮。执行 Utility MenuSelectEntities，在弹出的对话框中，将第一项选 Nodes，第二项选为 Attached to，下面默认为 Lines， all，点击 Invert，之后单击 OK按钮。之后就可以进行约束的加载了 。 第三 步 ， 执行 Utility MenuPlotNodes，显示出所需加约束的几个点。再执行Main Menu SolutionDefine LoadsApplyStructuralDisplacementOn Nodes命令，选取点，单击 OK按钮，在弹出的对话框中确定约束的大小和方向，这里方向选 ALL DOF选项，大小确定为零，单击 OK按钮，约束施加完毕。 曲轴两端的约束均是这样施加的。 接下来进行力的加载操作。力的加载方法与约束加载方法大致相同 。 第四 步 ， 执行 Utility MenuSelectEntities，在弹出的对话框中，将第一项选Lines，第二项选为 By Num/Pick，下面默认为 From Full，点击 Invert，之后单击 OK按钮。之后选取 所加力的线 ，单击弹出对 话框的 OK 按钮。执行 Utility MenuSelectEntities，在弹出的对话框中，将第一项选 Nodes，第二项选为 Attached to，下面默认为 Lines， all，点击 Invert，之后单击 OK按钮。 第五 步 ， 之后执行 Utility MenuPlotNodes，显示出所需加载力的几个节点。这时再执行 Main MenuSolutionDefineLoadsApplyStructuralForce/MomentOn 大庆石油学院本科生毕业论文 22 Nodes命令，选取点，单击 OK按钮，在弹出的对话框中确定力的大小和方向，单击 OK按钮，力的加载结束。 曲轴两个曲拐处的力均是这样加载 的 。 待力和约束都施加完毕以后，就可以进行力学求解了。先执行 Utility MenuSelectEverything 命令， 选择 曲轴 的所 有分析实体 。接着执行 Main MenuSolutionSolveCurrent LS命令，弹出 Solve Current Load Step对话框，单击OK按钮开始 求解，直到弹出 Solution is done对话框，单击 Close按钮，求解过程结束。 接下来就可以查看曲轴的分析结果了。 求解结束时，执行 Utility MenuFileSave as命令，弹出 Save Database对话框，在 Save Database to输入栏中输入 Results.db，保存上述的数据结果，单击 OK按钮，关闭该对话框 。 3.3.6 查看求解结果 图 3-4 外止点 曲轴应力分布云图 图 3-5 80 时 曲轴应力分布云图 图 3-6 280 时 曲轴应力分布云图 图 3-7 内止点 曲轴应力分布云图 想 要 查 看 变 形 结 果 ， 可 执 行 Main MenuGeneral PostprocPlot ResultsDeformed Shape命令，弹出 Plot Deformed Shape对话框，在 KUND Items to be plotted中选择 Def+undef edge选项，单击 OK按钮，屏幕上自动生成曲轴变形结果图，要注意，图上的变形结果为了让人可以清晰的观察是经过夸大化了的，不要与实际变形混为一谈，但其数值结果是近似正确的 。 大庆石油学院本科生毕业论文 23 查看应力结果，可执行 Main Menu General PostprocPlot ResultsContour PlotNode Solu命令，弹出 Contour Nodal Solution Data对话框，在 Item to be contoured选项框中选择 Stressvon Mises stress，其余选 项采用默认设置，单击 OK按钮，可查看曲轴的等效应力等值线图 。 通过本章对曲轴的有限元应力分析，我们可知：在分析过程中利用 pro/e 软件建立曲轴的三维实体模型，再将其导入 ANSYS 软件进行静力分析，模型所受载荷和边界条件根据实际进行简化；最后得到曲轴的应力分布云图和变形云图；曲轴满足强度和刚度要求。由上述分析结果可知，当曲轴在上述四种拉压力作用时，应力的最大值集中在曲轴主轴、曲柄销与曲柄的连接处。 对于上述应力集中的区域，可以考虑改变主轴、曲柄销与曲柄连接处的过渡圆弧的半径来减少应力集中，提高曲轴强度 。对于连接两曲柄销间曲柄应力集中及应变较大的情况，可以考虑通过改变曲柄外形尺寸，在应变较大处增加凸台，来达到提高曲柄强度降低应变的目的。由于曲轴结构的复杂性 ， 受载的空间变化性 ， 建立在某些假设基础上的常规简化计算有较大误差 ， 很难给出较精确的静强度和疲劳强度安全系数。然而 ， 曲轴的强度直接影响着空压机的正常运转。因此 , 产品设计时 ， 对各种型号的空压机曲轴进行三维有限元计算是非常必要的 。 在曲轴调质处理时 ， 要严格控制淬火加热温度 、 保温时间及冷却速度 ， 杜绝过热现象的出现 ， 在正确回火温度下进行回火 ， 保证回火保温时间 ， 通 过正确的热处理工艺 ，以得到回火索氏体组织 ， 确保曲轴具有良好的综合机械性能。另外 ， 合理地选择加工工艺、加工精度、表面光洁度也是提高曲轴使用寿命的重要途径。 3.4 小结 本章内容为： （ 1）运用理论计算 得到进行有限元分析所需的相关力，并 算得曲轴的关键截面为 10-10 截面 ，同时 为下一步进行曲轴的有限元分析打 下了 基础。 （ 2）运用 ANSYS 软件对 曲轴进行分析 。 通过对选取的四个时刻的受力情况进行 有限元分 析， 确定危险截面为 在曲轴主轴、曲柄销与曲柄的连接处。 （ 3） 通过以上分析可知，理论计算与 ANSYS 分析结果所确定的危险截 面一致， 但由于简化条件不同算得的具体 数值有所偏差。 应用 ANSYS10.0软件进行有限元分析，可以方便地得到在工作压力下的压缩机曲轴的应力分布 ， 为压缩机曲轴强度的校核提供了一个方便可靠的计算工具 ，同时为 活塞式 压缩机曲轴的设计和改进提供了理论依据。 大庆石油学院本科生毕业论文 24 第 4章 连杆有限元分析 4.1 对连杆进行有限元分析的重要性 连杆是压缩机运动机构中的重要零件之一 ，包括大头、小头、杆体三部分，大头一端与曲柄销相连，小头端与十字头销（或活塞销）相连。 其任务是与曲轴一起将输入压缩机的旋转运动转化为活塞的往复运动。一般中、大型压缩机连杆均采用 球墨铸铁或锻钢制成。随着我国石油、化学工业的迅速发展 ， 大型合成氨等需要更多的大型高压活塞式压缩机 ， 活塞力吨位也逐渐增加 ， 因而要求运动部件强度必须充分可靠。如连杆强度不足 ， 一旦破坏 ， 危害极大 ， 若盲目增大安全系数 ， 会使机器粗 笨 ， 造成浪费。过去对连杆 进行 强度计算 ， 一般采用材料力学中所述的方法 ， 分别对连杆大、小头及杆体进行校核。对杆体进行直杆拉、压校核。将连杆大、小头作为静不定曲梁 ， 其上作用均布或余弦载荷 ， 计算几个危险截面 ， 当然这和连杆的实际情况相差很大 ， 只能粗略地进行强度校核 ， 很难详细准确地算出各部位应力及位移 ， 无助于连杆的造型设计。近年来 ， 常常采用有限元法解算大量的 这类 工程实际问题。我们在 本次 设计 中分析 2D型 活塞 式 压缩机连杆时 ， 就 采用了有限元法。现将该机连杆计算结果作一分折 ， 以求对连杆强度计算有所改进 ， 以助于设计。 4.2 连杆相关力学计算 4.2.1 连杆尺寸的确定 图 4-1 连杆结构图 大庆石油学院本科生毕业论文 25 最大活塞力 p 256KNF 曲柄半径 32 1 6 c m22SR 连杆长度 85cmL 曲柄销直径 25cmD 大头孔直径 1 26.2cmD 小头衬套内径 14.5cmd 小头衬套宽度 14.5cmb 小头衬套厚度 ( 0 . 0 6 0 . 0 7 )Sd 0 . 0 6 1 4 . 5 0 . 8 7 c mS 取 0.8cmS 小头孔直径 1 2d d S 1 1 4 . 5 2 0 . 8 1 6 . 1d 取1 16cmd 连杆中间当量直径 m (1 . 6 5 2 . 4 5 )dP m 2 . 2 2 5 . 6 1 1 . 1 3 c md 连杆中间面积 2mm 4dF 2 2m1 1 . 1 3 9 7 . 2 c m4F 大庆石油学院本科生毕业论文 26 近小头处连杆直径 m0.9dd 0 . 9 1 1 . 1 3 1 0 . 0 2 c md 取 10cmd 近大头处连杆直径 m1.1dd 1 . 1 1 1 . 1 3 1 2 . 2 4 c md 取 12cmd 连杆宽度 由图测得 11.25cmB 大头处 AA 截 面 积 Am1 . 3 8 1 . 6FF（ ） 2A 1 . 5 9 7 . 2 1 4 5 . 8 c mF 取 2A 146cmF 截面 AA 的厚度 AA FS B AA 146 1 2 . 9 8 c m1 1 . 2 5FS B 取A 13cmS 大头处截面 BB 的面积 Bm(1 .3 1 .4 )FF 2B 1 . 3 9 7 . 2 1 2 6 . 3 6 c mF 取 2B 127cmF 截面 BB 厚度 BB FS B 大庆石油学院本科生毕业论文 27 B 127 1 1 . 2 9 c m1 1 . 2 5S 取 11.3cmBS 小头处截面 CC 面积 Cm( 0 .8 1 .0 )FF 2C 0 . 8 5 9 7 . 2 8 2 . 6 2 c mF 283cmCF 截面 CC 厚度 CC FS B C 83 7 . 3 8 c m1 1 . 2 5S 取C 7.4cmS 连杆螺栓 直径 0 ( 0 . 8 5 1 . 6 )dP 0 1 . 0 5 2 5 . 6 5 . 3 1 3 c md 取0 5.4cmd 螺栓定位部分直径 20 ( 0 . 1 0 . 2 )dd 2 5 4 0 . 2 5 4 . 2 c md 两连杆螺栓间距离 0 33cml 螺栓在大头体内长度 1 ( 0 .5 5 0 .6 5 )lD 1 0 . 6 2 5 1 5 c ml 螺栓在大头盖内长度 2 ( 0 .5 0 .6 5 )lD 大庆石油学院本科生毕业论文 28 2 0 . 6 2 5 1 5 c ml 连杆螺栓个数 2Z 连杆的计算 小头衬套比压 PP db (4-1) 256000 1 2 1 7 . 6 1 3 0 0 1 5 0 01 4 . 5 1 4 . 5P 连杆的惯性半径 m4di 1 1 . 1 32 . 7 84i 柔度 850 3 0 5 . 7 62 . 7 8Li 杆体的拉压应力 p mpF 2p 256000 2 6 3 3 . 7 4 N / c m9 7 . 2 系数 41 .3 7 1 0C 查活塞式压缩机设计表 5-14，连杆材料 35 钢 惯性矩 4mx 64dJ 4 4x3 . 1 4 1 1 . 1 37 5 2 . 9 c m64J 连杆摆动平面纵弯应力 2CB xLPCJ (4-2) 大庆石油学院本科生毕业论文 29 2 - 4 2CB852 5 6 0 0 0 1 . 3 7 1 0 3 3 6 . 6 N / c m7 5 2 . 9 杆体长度 11 22D dLL 1 2 6 . 2 1 68 5 5 5 . 9 c m22L 惯性矩 4my 64dJ 4 4y3 . 1 4 1 1 . 1 37 5 2 . 9 c m64J 垂直于连杆摆动平面纵弯应力 2 1CBy4LPCJ (4-3) 2 - 4 2CB5 5 . 92 5 6 0 0 0 1 . 3 7 1 0 3 6 . 4 N / c m4 7 5 2 . 9 在连杆摆动总应力 1 p cb 21 2 6 3 3 . 7 4 3 3 6 . 6 2 9 7 0 . 3 4 N / c m 垂直于连杆摆动平面总应力 2 p C B 22 2 6 3 3 . 7 4 3 6 . 4 2 6 7 0 . 1 4 N / c m 大头盖截面 AA 的抗弯截 面模数 2AA 6BSW 2 3A1 1 . 2 5 1 2 . 9 83 1 5 . 9 c m6W 截面 AA 弯曲应力 大庆石油学院本科生毕业论文 30 BA()24DplW (4-4) 22B252 5 6 0 0 0 ( 3 3 )2 4 1 5 3 . 2 N / c m 8 0 0 0 1 2 0 0 0 N / c m4 3 1 5 . 9 小于许用应力，安全 大头盖截面 BB 的抗弯模数 2BB 6BSW 2 3B1 1 . 2 5 1 1 . 32 3 9 . 4 c m6W 大头盖截面 BB 重心到连杆螺栓线的距离 10cmb 截面 BB 与连杆螺栓线的夹角 35 截面 BB 弯曲应力 B B2PbW (4-5) 2B 2 5 6 0 0 0 1 05 3 4 6 . 7 N / c m2 2 3 9 . 4 截面 BB 拉压应力 p Bsin2P F (4-6) 2p 2 5 6 0 0 0 s i n 3 55 7 8 . 1 N / c m2 1 2 7 截面 BB 剪应力 Bcos2P F (4-7) 22 5 6 0 0 0 c o s8 1 6 . 4 N / c m2 1 2 7 截面 BB 总应力 22pB( ) 4 大庆石油学院本科生毕业论文 31 2222( 5 7 8 . 1 5 3 4 6 . 7 ) 4 8 1 6 . 46 1 4 6 . 7 N / c m 6 0 0 0 8 0 0 0 N / c m 小头处截面 CC 的抗弯模数 2CC 6BSW 2 3C1 1 . 2 5 7 . 41 0 2 . 7 c m6W 小头侧壁中心间距 1Dl d S 1 6 4 2 0 c ml 截面 CC 的弯曲应力 BC()38dplW (4-8) 22B1 4 . 52 5 6 0 0 0 ( 2 0 )3 4 7 3 6 . 1 N / c m 6 0 0 0 8 0 0 0 N / c m8 1 0 2 . 7 截面 DD 弯曲应力 B D8PlW (4-9) 2B 2 5 6 0 0 0 2 06 2 3 1 . 7 N / c m8 1 0 2 . 7 截面 DD 的拉应力 p D2PF 2p 256000 1 5 4 2 . 2 N / c m2 8 3 截面 DD 总应力 PB 221 5 4 2 . 2 6 2 3 1 . 7 7 7 7 3 . 9 N / c m 8 0 0 0 1 0 0 0 0 N / c m 经过以上分析， 连杆的小头处截面积最小， 应力最大， 因此危险截面为连杆小头处。 大庆石油学院本科生毕业论文 32 4.3 连杆的有限元分析 4.3.1 模型建 立 本文使用 pro/e软件 建立了 2D型 压缩 机的 连杆有限元三维实体模型 ， 并将其导入 ANSYS软件 中进行了 受力 分析。 4.3.2 假设条件 连杆是由连杆体、连杆盖、连杆轴瓦和连杆螺栓等零件组成的。 断裂连杆的失效形式为疲劳断裂 ， 裂纹起源于连杆螺孔齿根或齿面应力集中部位。热处理工艺不当降低了连杆的屈服强度和抗拉强度 ， 使硬度也偏低 ， 在应力集中处产生裂纹源 ， 扩展后导致连杆断裂 ， 这是连杆失效的主要原因。 本文根据连杆的结构形状特点 ， 以有限元计算的数据准备工作量、求解时间及精度等为基本尺度 ， 在连杆建模时 ， 必须做一些简化 ， 这些简化是以不 影响连杆的动静力学特性为前提的。 1 实际过程中 ， 考虑到螺栓取了足够的预紧力 ， 在有限元分析模型中 ， 连杆与连杆大头盖在有限元分析中被认为连接成一体 ； 2 为了避免局部网格过密 ， 保证网格质量 ， 节省计算时间 ， 忽略了分布在连杆上的油孔 、 定位销钉孔等造成的影响 ； 3 作用在连杆上的所有载荷均平行于连杆摆动平面 ， 且认为沿连杆厚度方向均匀分布。 4.3.3 选择单元类型及定义材料属性 执行 Main MenuPreprocessorElement TypeAdd/Edit/Delete 命令，弹出Element Type对话框，单 击 Add按钮，在弹出的对话框中选择 Structural及其下的Solid，在其右边的滚动栏中选择 8node 45，单击 OK按钮。 执行 Main MenuPreprocessorMaterial Props Material Model命令，弹出 Define Material Model Behavior 对话框，在右边的菜单中，连续双击StructuralLinearElasticIsotropic 后又 出现一 个 Linear isotropic Properties for Material Number 1对话框，在 EX后面的输入栏中输入 2.1E+5，在 PRXY后面的输入栏中输入 0.3，单击 OK按钮，接着继续在 Define Material Model Behavior对话框中连续双击 StructuralDensity，弹出密度定义对话框，在 DENS后面的输入栏中输入7.85E-9，单击 OK按钮，关闭该对话框 。 到此 连杆模型的 单元类型的选择及材料属性的定义工作完毕。 大庆石油学院本科生毕业论文 33 4.3.4 网格划分 在这次设计中我采用的网格划分方法是自由网格划分。执行 Utility MenuSelectEntities，在 弹出的对话框中，将第一项选为 Volumes，第二项选为 By Num/Pick，下面默认为 From Full，点击 Invert，之后单击 OK按钮。之后选取整个连杆 实 体 ， 单 击 弹 出 对 话 框 的 OK 按钮。执行 Main MenuPreprocessorMeshingMeshTool，在 Smart Size这项选为 5。之后执行 Main Menu PreprocessorMeshing VolumesFree，之后选取整个连杆实体，单击弹出对话框的 Pick All按钮，等待其运行完毕，观察确认网格已经生成后，单 击工具条上的 SAVE_DB按钮，保存网格模型，操作到此为止网格划分完毕。 图 4-2 连杆网格划分 4.3.5 进入求解器加载并求解 执行 Main MenuSolutionAnalysis TypeNew Analysis命令，弹出 New Analysis对话框，选择分析类型为 Static，单击 OK按钮，关闭该对话框。 执行 Main MenuSolutionAnalysis TypeSoln Controls命令，弹出 Solution Controls对话框。单击 Basic选项卡，设定求解时间 1，打开自 动时间步长，子步数设为 40，最大子步数为 100。 之后单击 Solution Controls对话框上的 Nonlinear选项卡，在该对话框中，将 Linear search项置为 ON，同时将 DOF solution predictor项设为 On for all substp，单击 OK按钮，关闭该对话框，到此求解参数设定完毕。 连杆与曲轴是圆弧接触，因此认为力 与约束都 是作用在一条线上。 连杆的受力情况可简化为二力杆的受力情况，小头一端受力，大头一端进行约束。先进行约束的施加 。 第一步， 执行 Utility MenuSelectEntities，在弹出的对话框中，将第一项选大庆石油学院本科生毕业论文 34 Lines，第二项选为 By Num/Pick，下面默认为 From Full，点击 Invert，之后单击 OK按钮。 第二步， 之后选取所加力的线，单击弹出对话框的 OK按钮。执行 Utility MenuSelectEntities，在弹出的对话框中，将第一项选 Nodes，第二项选为 Attached to，下面默认为 Lines， all，点击 Invert，之后单击 OK按钮。 之后 就 可以进行约束 的加载了 。 第三步， 执行 Utility MenuPlotNodes，显示出所需加 约束 的几个点。再执行Main Menu SolutionDefine LoadsApplyStructuralDisplacementOn Nodes命令，选取点，单击 OK按钮，在弹出的对话框中确定 约束 的大小和方向， 这里方向选 ALL DOF选项，大小确定为零， 单击 OK按钮， 约束施加完毕 。 接下来进行力的加载操作。力的加载方法与约束加载方法大致相同 。 第四步， 执行 Utility MenuSelectEntities，在弹出的对话框中，将第一项选Lines，第二项选为 By Num/Pick，下面默认为 From Full，点击 Invert，之后单击 OK按钮。之后选取 所加力的线 ，单击弹出对 话框的 OK 按钮。执行 Utility MenuSelectEntities，在弹出的对话框中，将第一项选 Nodes，第二项选为 Attached to，下面默认为 Lines， all，点击 Invert，之后单击 OK按钮。 第五步， 之后执行 Utility MenuPlotNodes，显示出所需加 载 力的几个 节 点。这时 再执行 Main MenuSolutionDefineLoadsApplyStructuralForce/MomentOn Nodes命令，选取点，单击 OK按钮，在弹出的对话框中确定力的大小和方向，单击 OK按钮， 力的加载 结束。 待力和约束都施加完毕以后，就可以进行力学求解了。先执行 Utility MenuSelectEverything 命令， 选择连杆的所 有分析实体 。接着执行 Main MenuSolutionSolveCurrent LS命令，弹出 Solve Current Load Step对话框，单击OK按钮开始求解，直到弹出 Solution is done对话框 ，单击 Close按钮，求解过程结束。 求解结束时， 执行 Utility MenuFileSave as命令，弹出 Save Database对话框，在 Save Database to输入栏中输入 Results.db，保存上述的数据结果，单击 OK按钮，关闭该对话框。 接下来就可以查看连杆的分析结果了 。 4.3.6 查看求解结果 想 要 查 看 变 形 结 果 ， 可 执 行 Main MenuGeneral PostprocPlot ResultsDeformed Shape命令，弹出 Plot Deformed Shape对话框，在 KUND Items to be plotted中选择 Def+undef edge选项，单击 OK按钮，屏幕上自动生成连杆变形结果图，要注意，图上的变形结果为了让人可以清晰的观察是经过夸大化了的，不要与实际变形混为一谈，但其数值结果是近似正确的。 查看应力结果，可执行 Main Menu General PostprocPlot ResultsContour PlotNode Solu命令，弹出 Contour Nodal Solution Data对话框，在 Item to be 大庆石油学院本科生毕业论文 35 contoured选项框中 选择 Stressvon Mises stress，其余选项采用默认设置，单击 OK按钮，可查看等效应力等值线图。 图 4-3 外止点 连杆应力分布云图 图 4-4 80 时 连杆应力分布云图 图 4-3 280时 连杆应力分布云图 图 4-4 内止点 连杆应力分布云图 通过本章对连杆的有限元应力分析，我们可知：在分析过程中利用 pro/e 软件建立连杆的 三维实体模型 ，再将其 导入 ANSYS 软件进行静力分析 ， 模型所受载荷和边界条件根据实际进行简化 ， 最后得到连杆的应力分布云图和变形云图。 在曲轴运行到 290 时应力最大。 由上述分析计算结果可知 ,当连杆在受拉状态下时 ，应力的最大值集中在连杆的小头处。同时 ， 在连杆杆身与大头连接处同样存在较大的应力。通过计算分析可以得出 ， 在设计压力下工作的压缩 机 连杆满足强度和刚度要求。对于上述应力集中的区域 ， 可以考虑在连杆小头处加工一个凸台 ， 以加强此处的强度。同时 ， 通过改变连杆杆身与大小头连接处的过渡圆弧的半径来减少应力集中 ， 提高连杆强 度。除此之外 ， 为防止连杆失效还有以下建议：严格安装工艺，螺丝紧固均匀，避免偏载，避免附加应力；机加工时保证台肩根部有足够的圆弧过度；注意控制材质的强度，使强度和塑、韧性有恰当配合。 4.4 小结 本章内容为： （ 1）运用理论计算算得连杆 有限元 分析所需的相关力。 并分析确定危险截面大庆石油学院本科生毕业论文 36 为 连杆的小头处 ，同时 为下一步进行连杆的有限元分析打 下了 基础。 （ 2）运用 ANSYS 软件对连杆进行 分析 。通过对选取的四个时刻的受力情况进行有限元分析， 应力的最大值集中在连杆的小头处。同时 ， 在连杆杆身与大头连接处同样存在较大的应力。 （ 3）通过以 上 几节的 分析可知，理论计算与 ANSYS 分析结果所确定的危险截面一致， 但由于简化条件不同 数值有所偏差。 大庆石油学院本科生毕业论文 37 第 5章 十字头有限元分析 5.1 对十字头进行有限元分析的重要性 十字头也称滑块， 作为大中型压缩机主要运动部件之一 ， 它的作用是实现运动方式的转换 ， 因而其受力较为复杂 ， 对结构设计、材质和制造过程要求较高。十字头是连接连杆和活塞杆的零件，因活塞杆与十字头销在十字头内成十字而得名 ，它 将曲柄机构的旋转运动 ， 通过连杆和十字头的运动转换变成直线运动 ， 以推动活塞的功能运动。由于十字头在这个运动过程中受力情况比连杆、活塞杆都相对复杂 ， 而且 它和滑道构成了一对运动副 ， 所以对十字头的结构、质量、材料、制造工艺等都提出了很高的要求。不仅如此 ， 十字头还是 D系列压缩机的易损件之一。通过长期的探讨和研究 ， 关于十字头的基本结构、材质、润滑系统以及运动配合间隙等 ， 都建立了行之有效的方案 ， 所以只要零件加工精度达到要求 ， 装配合格 ， 正常情况下 ， 试车前阶段不会出现故障。现在的主要毛病在于使用一段时间后 (有时甚至 1、 2 个 h)， 十字头发生事故的实例屡见不鲜。事实证明 ， 十字头失效通常是由于润滑故障引起的。如何比较准确地得到应力、变形的大小及分布和十字头的固有频率及振型 ， 对于指导十字头的设计和改进具有重要意义。 5.2 十字头相关力学计算 5.2.1 十字头尺寸确定 表 5-1 十字头尺寸确定 序号 名称 符号 单位 公式、参数的选择及计算 1 名义活塞力 PF kN P 256F 2 滑履直径 D cm 39.6D 3 滑履长度 L cm 47.2L 4 滑履宽度 B cm 25B 5 十字头销直径 d cm 14.5d 6 十字头销孔座厚度 S cm 7.6S 7 十字头体厚度 1S cm 1 4.4S 大庆石油学院本科生毕业论文 38 续表 5-1 表 5-1 十字头尺寸确定 8 与活塞杆连接处直径 0d cm 0 14.8d 9 活塞杆面积 A 2cm 70.8A 10 十字头颈部壁厚 2S cm 200 .3 6 5 .3Sd 11 十字头颈部实际直径 2d cm 2 1 22 4 2 . 2d d S 12 截面 AA 面积 A 2cm 22( 4 2 . 2 3 1 . 6 ) 6 1 44AA 13 截面 BB 面积 BA 2cm 0 . 2 5 0 . 2 5 2 5 6 6 4BAP 14 截面 BB 宽度 b cm 64 8.47.6b 5.2.2 十字头强度与比压校核 图 5-1 十字头计算图 活塞力 P 256kNF 十字头销孔座平均直径 m 14.5cmd 十字头销孔座比压 大庆石油学院本科生毕业论文 39 3Pqm2FF dS (5-1) 322q 256000 1 1 6 2 / 3 5 0 0 4 5 0 0 N / c m2 1 4 . 5 7 . 6F N c m 连杆径长比 0.188 侧向力 nPFF (5-2) n 0 . 1 8 8 2 5 6 4 8 . 1 3 k NF 十字头总重 2.9kNG 滑履比压 nq1 FGF LB (5-3) q 1 a4 8 . 1 3 2 . 94 3 P4 7 . 2 2 5F 连杆小头衬套宽度 14.5cm 十字头销支承座距离 l b S 1 4 . 5 7 . 6 2 2 . 1 c ml 十字头与活塞杆连接处活塞杆柔度 S lEA (5-4) 8S 712 0 . 8 1 02 . 1 1 0 7 0 . 8 十字头销弯曲应力 bb 44016 ()( ) 2 4PdF lbdd (5-5) 2bb 41 6 1 4 . 5 2 5 6 0 0 0 2 2 . 1 1 4 . 5( ) 3 1 7 7 N / c m1 4 . 5 2 4 十字头销为实心，所以0 0d 十字头体柔度 大庆石油学院本科生毕业论文 40 c lEA (5-6) 8c 74 7 . 2 0 . 3 7 1 02 . 1 1 0 6 1 4 外载荷系数 csc (5-7) 880 . 3 7 1 0 0 . 3 1 6( 0 . 8 0 . 3 7 ) 1 0 预紧系数 23K 取 2K 预紧力 tPF KF (5-8) t 2 2 5 6 5 1 2 k NF 螺纹拉伸力 Q t PF F F (5-9) Q 5 1 2 0 . 3 1 6 2 5 6 5 9 2 . 9 k NF 螺母面积 22f 1 0()4A d d 2 2 2f (1 3 . 4 9 . 5 ) 7 0 . 1 c m4A 螺母受拉伸应力 tbfFA 22b 512 7 3 0 3 N / c m 1 0 0 0 0 N / c m7 0 . 1 螺母台肩压碎应力 Qbc 2221()4Fdd 大庆石油学院本科生毕业论文 41 22bc 22592896 4 7 1 . 7 N / c m 7 0 0 0 N / c m( 4 2 . 2 1 3 . 4 )4 截面 AA 的面积 2 2 2A ( 4 2 . 2 3 1 . 6 ) 6 1 4 c m4A 截面 AA 拉伸应力 QPAFA 22P 592896 9 6 5 . 6 N / c m 6 0 0 0 N / c m614 截面 BB 重心 m2dbr (5-10) 1 4 . 5 8 . 41 1 . 5 c m2r 当量长度 2 2 3cmlr 截面 BB 的抗弯模量 21B 6SbW 2 3B4 . 4 8 . 4 5 1 . 7 c m6W 截面 BB 最大弯曲应力 Pmbb ()1 6 3FdlW (5-11) 22bb 2 5 6 0 0 0 1 4 . 5( 2 3 ) 5 5 8 7 N / c m 1 0 0 0 0 N / c m1 6 5 1 . 7 3 截面 CC 的面积 22C 0 0()8A D d 2 2 2C ( 2 8 2 3 ) 1 0 0 c m8A 截面 CC 到孔销中心线距离 5.5cma 图中测得 截面 CC 与活塞杆轴夹角 大庆石油学院本科生毕业论文 42 35 图中测得 截面 CC 抗弯模量 34 3C 42 3 2 3(1 ) 3 2 5 . 3 c m6 4 2 8W 截面 CC 弯曲应力 PbbC2FaW (5-12) 2bb 2 5 6 0 0 0 5 . 5 2 1 6 4 N / c m2 1 0 0 截面 CC 拉伸应力 PbCsin2F A (5-13) 2b 2 5 6 0 0 0 s i n 3 57 3 3 . 4 N / c m2 1 0 0 截面 CC 切应力 PCcos2F A (5-14) 22 5 6 0 0 0 c o s 3 51 0 4 7 . 5 N / c m2 1 0 0 截面 CC 总应力 22b b b( ) 4 (5-15) 2222( 2 1 6 4 7 3 3 . 4 ) 4 1 0 4 7 . 52 8 9 8 N / c m 6 0 0 0 N / c m 通过计算发现十字头与活塞杆连接处承受最大应力，该处为危险截面 5.3 十字头的有限元分析 5.3.1 模型建立 本文 使用 pro/e软件 建立了 2D型压缩机的十字头有限元三维实体模型 ， 并将其导入 ANSYS软件 中进行了受力分析。 大庆石油学院本科生毕业论文 43 5.3.2 假设条件 1 实际过程中 ， 考虑到螺栓取了足够的预紧力 ， 在有限元分析模型中 ， 整个十字头 在 有限元分析中被认为连接成一体 ； 2 为了避免局部网格过密 ， 保证网格质量 ， 节省计算时间 ， 忽略了分布在 十字头 上 各种孔道 等造成的影响 ； 3 作用在 十字头 上的所有载荷均平行于 十字头运动 平面 ， 且认为沿 十字头 厚度方向均匀分布。 5.3.3 选择单元类型及定义材料属性 执行 Main MenuPreprocessorElement TypeAdd/Edit/Delete 命令，弹出Element Type对话框，单击 Add按钮，在弹出的对话框中选择 Structural及其下的Solid，在其右边的滚动栏中选择 8node 45，单击 OK按钮。 执行 Main MenuPreprocessorMaterial Props Material Model命令，弹出 Define Material Model Behavior 对话框，在右边的菜单中，连续双击StructuralLinearElasticIsotropic 后又 出现一 个 Linear isotropic Properties for Material Number 1对话框，在 EX后面的输入栏中输入 2.1E+5，在 PRXY后面的输入栏中输入 0.3，单击 OK按钮，接着继续在 Define Material Model Behavior对话框中连续双击 StructuralDensity，弹出密度定义对话框，在 DENS后面的输入栏中输入7.85E-9，单击 OK按钮，关闭该对话框。 到此连杆模型的单元类型的选择及材料属性的定义工作完毕。 5.3.4 网格划分 图 5-2 十字头网格划分 大庆石油学院本科生毕业论文 44 执行 Utility MenuSelectEntities，在弹出的对话框中，将第一项选为 Volumes，第二项选为 By Num/Pick，下面默认为 From Full，点击 Invert，之后单击 OK按钮。之后选取整个连杆实体，单击弹出对话框的 OK 按钮。执行 Main MenuPreprocessorMeshingMeshTool，在 Smart Size这项选为 5。之后执行 Main Menu PreprocessorMeshing VolumesFree，之后选取整个连杆实体，单击弹出对话框的 Pick All按钮，等待其运行完毕，观察确认网格已经生成后，单击工具条上的 SAVE_DB按钮，保存网格模型，操作到此为止网格划分完毕。 5.3.5 进入求解器加载并求解 执行 Main MenuSolutionAnalysis TypeNew Analysis命令，弹出 New Analysis对话框，选择分析类型为 Static，单击 OK按钮，关闭该对话框。执行 Main MenuSolutionAnalysis TypeSoln Controls命令，弹出 Solution Controls对话框。单击 Basic选项卡 ，设定求解时间 1，打开自动时间步长，子步数设为 40，最大子步数为 100。之后单击 Solution Controls对话框上的 Nonlinear选项卡，在该对话框中，将 Linear search项置为 ON，同时将 DOF solution predictor项设为 On for all substp，单击 OK按钮，关闭该对话框，到此求解参数设定完毕。 十字头 与 连杆是圆弧接触，因此认为力与约束都是作用在一条线上。十字头的受力情况可简化 为下端约束，十字头销处受力 。先进行约束的施加 。 第一步， 执行 Utility MenuSelectEntities，在弹出的对话框中，将第一项选Lines，第二项选为 By Num/Pick，下面默认为 From Full，点击 Invert，之后单击 OK按钮。之后选取所加力的线，单击弹出对话框的 OK按钮。 第二步， 执行 Utility MenuSelectEntities，在弹出的对话框中，将第一项选Nodes，第二项选为 Attached to，下面默认为 Lines， all，点击 Invert，之后单击 OK按钮。之后就可以进行约束的加载了 。 第三步， 执行 Utility MenuPlotNodes，显示出所需加约束的几个点。再执行Main Menu SolutionDefine LoadsApplyStructuralDisplacementOn Nodes命令，选取点，单击 OK按钮，在弹出的对话框中确定约束的大小和方向，这里方向选 ALL DOF选项，大小确定为零，单击 OK按钮，约束施加完毕。接下来进行力的加载操作。 第 四 步 ， 力 的 加 载 方 法 与 约 束 加 载 方 法 大 致 相 同 ， 执 行 Utility MenuSelectEntities，在弹出的对话框中，将第一项选 Lines，第二项选为 By Num/Pick，下面默认为 From Full，点击 Invert，之后单击 OK按钮。之后选取所加力的线，单击弹出对话框的 OK按钮。执行 Utility MenuSelectEntities，在弹出的对话框中，将第一项选 Nodes，第二项选为 Attached to，下面默认为 Lines， all，点击Invert，之后单击 OK按钮。 第五步， 之后执行 Utility MenuPlotNodes，显示出所需加载力的几个节点。大庆石油学院本科生毕业论文 45 这时再执行 Main MenuSolutionDefineLoadsApplyStructuralForce/MomentOn Nodes命令，选取点，单击 OK按钮，在弹出的对话框中确定力的大小和方向，单击 OK按钮，力的加载结束。 待力和约束都施加完毕以后，就可以进行力学求解了。先执行 Utility MenuSelectEverything命令，选择 十字头 的所有分析实体。接着执行 Main MenuSolutionSolveCurrent LS命令，弹出 Solve Current Load Step对话框，单击OK按钮开始求解，直到弹出 Solution is done对 话框，单击 Close按钮，求解过程结束。求解结束时，执行 Utility MenuFileSave as命令，弹出 Save Database对话框，在 Save Database to输入栏中输入 Results.db，保存上述的数据结果，单击 OK按钮，关闭该对话框。接下来就可以查看 十字头 的分析结果了。 5.3.6 查看求解结果 图 5-3 外止点 十字头应力分布云图 图 5-4 80 时 十字头应力分布云图 图 5-5 280 时 十字头应力分布云图 图 5-6 内止点 十字头应力分布云图 想 要 查 看 变 形 结 果 ， 可 执 行 Main MenuGeneral PostprocPlot ResultsDeformed Shape命令，弹出 Plot Deformed Shape对话框，在 KUND Items to be plotted中选择 Def+undef edge选项，单击 OK按钮，屏幕上自动生成 十字头 变形结果图，要注意，图上的变形结果为了让人可以清晰的观察是经过夸 大化了的，不大庆石油学院本科生毕业论文 46 要与实际变形混为一谈，但其数值结果是近似正确的。 查看 十字头 应力结果，可执行 Main Menu General PostprocPlot ResultsContour PlotNode Solu命令，弹出 Contour Nodal Solution Data对话框，在Item to be contoured选项框中选择 Stressvon Mises stress，其余选项采用默认设置，单击 OK按钮，可查看等效应力等值线图。 通过本章对十字头的有限元应力分析，我们在在分析过程中利用 pro/e软件建立十字头的三维实体模型 ，再将其 导入 ANSYS软件进行静力分析 ， 模型所受载荷和边界条件根据实际进行简化 ； 最后得到十字头的应力分布云图和变形云图。 在曲轴 290度时应力最大。 可见十字头的小头凸台底处为危险截面。 对 以上分析结果可能引起的十字头 的各种失效因素及处理方法都有 了一定 了解。 对分析结果处失效的原因 现总结如下。 首先，十字头和滑道加工精度不够。 其次是