文档介绍：作业一:有限元分析实例
实例:请对一个盘轴配合机构进行接触分析。轴为一等直径空心轴,盘为等厚度圆盘,其结构及尺寸如图所示。盘和轴为一种材料,材料参数为:弹性模量Ex=,泊松比NUXY=,摩擦系数MU=,试采用有限元计算方法分析轴和盘在过盈配合时的应力应变分布以及将轴从盘心拔出时轴和盘的接触情况。
问题分析说明
本题主要分析装配过程中结构的静态响应,所以分析步选择通用静态分析步。 由于为过盈配合,属于大变形,故应考虑几何非线性的影响。
 模型具有轴对称性,所以可以采取轴对称模型来进行分析,先建立二维模型计算,再转换为三维模型计算,这样可以节省计算时间。
 分析过程由两个载荷步组成, 第一个载荷步为过盈分析, 求解过盈安装时的情况。第二个载荷步为将轴从盘心拔出时的接触分析, 分析在这个过程中盘心面和轴的外表面之间的接触应力。它们都属于大变形问题, 属于非线性问题。在分析时需要定义一些非线性选项来帮助问题的收敛。
接触面之间有很大的相对滑动,所以模型要使用有限滑移。
模型建立的分析说明
进定义单元类型此项实例分析的问题中涉及到大变形, 故选用 So li d185 单元类型来建立本实例入部件模块,的模型。盘轴接触问题属于面面接触, 目标面和接触面都是柔性的, 将使用接触单元 T ARGET 170 和CO NTAT17 4来模拟接触面。分别创建名为为part1、part2的部件。
定义材料属性,在线性各向同性材料属性对话框中的EX (弹性模量) 文本框中输入 2 . 5E5,PRX Y (泊松比) 文本框中输入 0 . 3,并将定义的材料属性赋予给part1和part2。如下图所示。
(3)进入装配模块,创建两者间的装配关系。
(4)进入分析步模块 
定义名为step1和step2的两个分析步。
(5)进入相互作用模块,创建相互作用属性,设置摩擦系数;然后定义接触关系。如下图所示。
 
(6)进入载荷模块,创建边界条件,依次定义名为BC-2(类型为:完全固定)、BC-3(类型为:位移/转角,约束U1、UR3),分析步均为Initial 。定义名为BC-4(类型为:位移/转角,U2=,给微量位移,先接触上)的边界条件,分析步为step1,并inactive分析步 step2。最后创建名为BC-5的边界条件,分析步为step2,类型为位移/转角,U2=50。如下图所示。

(7)进入网格模块,划分网格。结果如上右图所示。
(8)进入作用模块,提交作业。
(9) 进入可视化模块,观察运行结果,下图为轴插入盘心时轴和盘的接触情况。
1)分析步刚开始时Mises应力图如下所示:
下图为应力最大时刻。
2)分析步中间某时刻,Mises应力图如下所示:
3)分析步最后的Mises应力图如下所示:
有限元计算分析说明
,在轴压入盘心的过程中,轴上的最大应力值发生的位置为轴端以及轴与盘心端面接触的位置,+04,盘心的最大应力出现在轴进入端之后,但和报告文件不符合,不知为何。
,+03,-03,,最小应变为--03。
报告文件
. In-P
***@Loc 1 ***@Loc 1
-----------------------------------------------------------
最小--03
在单元 748 913
结点 757 924
最大 -03 +03
在单元 1993 84
结点 2101 169
总 +06
下图为接触配合面的某一节点应力应变图。
下为轴和盘在过盈配合时的应变分布图。
下为轴和盘在过盈配合时的应力分布图。
对非线性有限元方法的认识
通过对非线性有限元方法与应用这门课的学****我们可以了解到,有限元法的力学基础是弹性力学,方程求解的原理是泛函极值原理(即变分原理),实现的方法是数值离散技术,最后的技术载体是有限元分析软件(例如ANSYS、ABAQUS等)。有限元方法的基本原理
是将连续的结构离散化,用有限个容易分析的单元来表示复杂的对象,单元之间通过有限个节点按一定方式相互连接,形成单元的组合体。
(确定单元形状、单元之间相互联结的点称为节点。单元上节点处的结构内力为节点力,外力(有集中力、分布力等)为节点载荷。)由于单元的数目是有限的,节点的数目也是有限的,所以称为有限元法。