大连理工大学硕士研究生有限元课程大作业参考模板.doc

PAGE PAGE 1盘轴紧配合结构分析214000000 Alex 183******** 1*********@摘要接触是一种高度非线性行为，需要较多的计算资源，为了进行有效的计算，理解问题的特性和建立合理的模型是很重要的。在本文中，将在Ansys环境中用有限元分析方法对一个盘轴紧配合结构进行接触分析。第一个载荷步分析轴和盘在过盈配合时的应力，第二个载荷步分析将该轴从盘心拔出时轴和盘的接触应力情况。关键词: 接触，盘轴，有限元方法，接触应力0.前言有限元法(Finite Element Method，FEM)，是计算力学中的一种重要的方法，它是20世纪50年代末60年代初兴起的应用数学、现代力学及计算机科学相互渗透、综合利用的边缘科学。接触是一种高度非线性行为，需要较大的计算资源，为了进行有效的计算，理解问题的特性和建立合理的模型是很重要的。本文中将通过轴盘点接触（面面接触）说明接触分析的方法。1.问题描述在旋转机械中通常会遇到轴与轴承、轴与齿轮、轴与盘连接的问题，根据各自的不同情况可能有不同的连接形式。但大多数连接形式中存在过盈配合，也就是涉及到接触问题的分析。这里我们以某转子中轴和盘的连接为例，分析轴和盘的配合应力以及将轴从盘中拔出时盘轴连接处的应力情况。本实例的轴为一等直径空心轴，盘为等厚度圆盘，其结构及尺寸如图1.1所示。由于模型和载荷都是轴对称的，可以用轴对称方法进行分析。这里为了后处理时观察结果更直观，我们采用整个模型的四分之一进行建模分析，最后将其进行扩展，来观察整个结构的变形及应力分布、变化情况。盘和轴共用同一种材料，其性质如下： 弹性模量： EX=2.5E5 泊松比： NUXY=0.35 接触摩擦系数 ： MU=0.2 图1.1 盘轴结构图1.系统建模2.1 建立几何模型并划分网格调整后的两个圆环几何体如图2.1所示。读者可以单击对话框下部的动态显示选择框，然后通过鼠标右键来调整视图位置。图 2.1 建立的盘轴几何模型对盘进行网格划分,建立的有限元模型如图2.2所示。图 2.2 盘轴结构的有限元模型2.2 创建接触对图形输出窗口中显示的接触对，如图2.3所示。图2.3 定义的接触对3.加载求解因为本实例是对盘轴结构中，将轴从盘心拔出的过程进行非线性分析。根据条件知道在这个过程中盘的外缘节点的自由度该全部约束。又由于建模时为了节约费用，根据模型结构的特点，只建立了四分之一的模型，所以需要在分析时定义对称边条，来模拟真实的情况。本实例的分析过程由两个载荷步组成，第一个载荷步为过盈分析，求解盘轴过盈安装时的应力情况。第二个载荷步为将轴从盘心拔出时的接触分析，分析在这个过程中盘心面和轴的外表面之间的接触应力。它们都属于大变形问题，属于非线性问题。在分析时需要定义一些非线性选项来帮助问题的收敛。下面进行本实例的加载和求解操作。3.1 定义轴对称边条和位移约束1．定义轴对称边条。选取菜单路径Main Menu | Solution | Define Loads | Apply | Structural | Displacement | Symmetry B. C. | On Areas，将弹出Apply SYMM on Areas (在面上施加轴对称边条) 拾取对话框。在图形输出窗口中单击选取盘和轴的四个径向截面，如图3.1所示。然后，单击拾取对话框中的按钮“OK”关闭对话框，对它们施加轴对称边条。图3.1 定义轴对称边条2. 定义盘外缘的位移约束。3. 在对话框中，单击DOFs to be constrained (约束自由度)列表框中的“All DOF”，使其高亮度显示，将其选中。其余设置保持缺省值(缺省的位移值为0)，单击“OK”按钮关闭拾取对话框，完成对位移约束的定义。3.2 定义并求解第一个载荷步1. 指定分析类型。选择菜单路径Main Menu | Analysis Type | Solution | New Analysis，将会弹出New Analysis (新分析)对话框(图略)，单击对话框中的“Static”单选按钮，指定分析类型为静力分析。2. 设定分析选项。3．单击Basic (基本)选项卡左边的Analysis Options (分析选项)下拉框中的“Large Displacement Static”将其选中，使分析中考虑大变形影响。然后在Time Control (时间控制区)中设定载荷步结束时间Time at end of load step＝100，并关掉自动时间步(Automatic time stepping为“Off”)。其余设置保持缺省，然后单击对话框中的“OK”按钮关闭对话框。4. 进行第一载荷步求解。5. 仔细阅读求解状态窗口中的信息，检查所有的选项和设