关于单元划分的讨论.docVIP

有限元分析中单元划分的讨论 伍妮妮 （长江大学机械工程学院） 摘要 有限元分析(FEA,Finite ElementAnalysis)为解决复杂的工程分析、计算问题提供了有效的途径,为新产品的开发提供了最佳的优化方案。一个有限元模型分析的成败关键在于有限元分析的前处理,计算以及后处理,其中前处理中单元划分尤其重要,直接决定着分析的精度与效率,对此以ANSYS为例分析探讨了在建立有限元分析模型时单元的划分,为新产品的设计提供了参考。 关键词 有限元分析;单元划分; 0引言 划分单元是建立有限元模型的一个重要环节，它要求考虑的问题较多，需要的工作量较大，所划分的网格形式对计算精度和计算规模将产生直接影响。要建造更为快速的交通工具、更大规模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备,要求工程师在设计阶段就能精确地预测出产品和工程的技术性能,需要对结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。近年来在计算机技术和数值分析方法支持下发展起来的有限元分析(FEA, Finite ElementAnalysis)方法则为解决这些复杂的工程分析计算问题提供了有效途径。FEA一般包括前处理、计算、后处理三个步骤其中前处理的好坏直接影响到后面的计算,前处理中单元的划分直接决定着分析的精度与效率,面对纷繁复杂的工程问题,工程师们建立了各种各样的单元类型,如何在解决实际工程问题中恰当的选择合适的单元划分方式,需要经验与知识的积累,本文针对常用有限元软件ANSYS与有限元理论中常用单元类型以及日常工程项目中的实际经验相结合,探讨了各单元类型的适用范围与工程问题中如何正确选择合适的单元划分,从而使在今后的有限元分析中少走弯路。[1] 1单元类型 单元类型是因研究对象的力学特性不同而定义的，单元类型的选择,跟你要解决的问题本身密切相关。在选择单元类型前,首先要对问题本身有非常明确的认识,然后,对每一种单元类型,每个节点的自由度数,特性,使用条件在ANSYS的帮助文档中都有非常详细的描述,要结合自己的问题,对照帮助文档里面的单元描述来选择恰当的单元类型。 (1)按形状分类,如图1 ~3所示。点单元:MASS;线单元: LINK、BEAM、COMBIN;面单元:PLANE、SHELL。 (2)按单元阶次分类:线性单元:对于结构分析问题,单元内的位移数值按线性变化,因而每个单元内的应力状态是保持不变的;二次单元:对于结构分析问题,单元内的位移数值按二次函数变化,因此每个单元内的应力状态是线性变化的; P单元:对于结构分析问题,单元内的位移数值按二阶到八阶函数变化,而且具有求解收敛自动控制功能,自动确定各位置上应采用的函数阶数[2]。如图4所示。 1.1点单元 几何形状为点型的结构,可用以下单元模拟,如图6所示。MASS单元主要用于动力学分析质量块结构的模拟[3]。 1.2　线单元 几何形状为线型的结构,可以用以下单元模拟,如图7、8所示。Link单元用于桁架、螺栓、螺杆等连接件的模拟。Beam单元用于梁、螺栓、螺杆、连接件等的模拟。Pipe单元用于管道、管件等结构的模拟;Combin单元用于弹簧,细长构件等的模拟。 1.3　面单元 几何形状为面型的结构,可用以下单元模拟。 (1) SHELL单元:主要用于薄板或曲面结构的模拟,壳单元分析应用的基本原则是每块面板的主尺寸不低于其厚度的10倍。 (2) PLANE单元:用于总体直角坐标系下X-Y平面内结构的应力、应变、轴对称问题,其中平面应变适用于如图9所示状态中。Z方向上的几何尺寸远远大于X和Y方向上尺寸的情况。Z方向上的应变为0,但应力不为0。所有的在和均作用在XY平面内。仅允许XY平面内的运动。在Z方向具有无线厚度。 轴对称适用于X-Y面内的横截面绕Y轴旋转360°所形成的三维实体结构,利用对称性只需建立X-Y平面内的横截面模型,不允许在-X方向有模型结构,坐标轴规定:X方向是径向, Y方向是轴向,Z方向是周向,周向位移是0,而周向应力、应变十分明显,只能考虑轴向载荷。 1.4　体单元 几何形状为体型的结构,可用SOLID单元模拟,主要用于三维实体结构的模拟。 2单元类型的选择 ANSYS的单元库提供了100多种的单元类型,单元类型选择的工作就是将单元的选择范围缩小到少数几个单元上,通常采用以下方法。 (1)设定物理场过滤菜单,将单元全集缩小到该物理场涉及的单元[4]。 (2)根据模型的几何形状选定单元的大类,如线性结构则只能用“LINK、BEAM、PIPE和COMBIN”这类单元去模拟,面结构则只能用“PLANE、SHELL”单元去模拟。 (3)根据模型结构的空间维数细化单元类别,如确定为“BEAM”单元大类后,在对话框的右栏中有“2D”“3D”的单