有限元基础与ANSYS应用.pptVIP

南京农业大学工学院机械工程系 * * * * * * 浅谈Wachspress型多边形有限元法 姓名：贾亮景 学号：2016205018 1. Wachspress型多边形有限元法简述2. Wachspress型多边形有限元法积分方案基本理论3. 数值算例4. 结论 Wachspress型多边形有限元法简述 在二维问题中，传统有限元方法常将单元划分成常应变三结点三角形单元和双线性四结点四边形单元。多边形单元作为一种新的单元形式,对于一些复杂几何形状的区域，具有区域适应性强和前处理方便等优点。在对多晶体材料 的数值分析中，使用多边形单元能有效降低分析单元数量，大大减少了计算所需时间。在模拟裂纹扩展的问题中，应用多边形有限元只需将裂纹穿过的单元分成多边形单元。克服了传统有限元法在裂纹区域需要不断重新划分网格的缺点。 1975年，Wachspress首先提出了具有有理函数形式的凸多边形单元形函数，称为Wachspress插值法，但因为其构造方法过于复杂而没有得到重视. 从20世纪90年代，工程数值计算和计算机图形学等领域学者又对多边形单元形函数构造方法做了大量的研究,取得了许多重要的成果。提出了一些改进后的Wachspress型插值多边形形函数。 Wachspress型多边形有限元法 Wachspress插值形函数 由于多边形单元的边数是任意的，因此构造满足单元位移协调性的多项式形式的位移插值函数变得异常困难.1975年,Wachspress利用映射几何技术成功地在凸多边形上构造出具有有理函数形式的多边形单元形函数,但因为其构造过程繁琐而未能得到广泛的应用.对于具有n条边的多边形,其表达式： 式中， 为权函数； Wachspress型多边形有限元法 Wachspress型插值多边形形函数经过近十几年不断改进，已经成功地运用到多边形有限单元法中。目前比较常见的Wachspress多边形形函数有3种，图中的多边形单元权函数表达式分别为: 式中 为点 之间的距离;A表示三点围成的三角形的面积。 Wachspress型多边形有限元法 2. 积分方案 ? 多边形有限元积分方案 多边形单元由于边数不定,无法直接在整个单元上进行积分,现行处理该积分的主要方法是将多边形单元三角形化,如图所示。然后在每个三角形上进行积分,其数学公式为： 注意：现有数值积分方法在处理多项式函数积分时可获得精确值,而对非多项式函数积分理论上并不完善.对于上式,由于被积函数 中包含多边形单元形函数表达式,而多边形单元的形函数又为非多项式函数形式，因此，在理论上多边形有限元计算结果的精度要受到影响。 Wachspress型多边形有限元法 ? 三角形中的Gauss积分 针对平面三角形区域积分，对于每一个三角形区域，积分计算式可表示为： 式中， 为三角形面积坐标; 为Jacobi行列式。 为了将 Gauss 积分应用到多边形单元中，在上式的基础上通过积分的换元法对积分域进行变换,可得到 Gauss 积分式为： 式中， 为 Gauss 积分中的权函数; 为 Gauss 积分中的积分点坐标。 Wachspress型多边形有限元法 上式中积分函数关于局部坐标L1 ,L2对称， 因此在L1和L上应选择相同积分点数,但该三角 形中的 Gauss 积分点分布并不对称，因此取直角边 为单位长度的等腰直角三角形,其不同Gauss积分 点数分布如图所示 本文中所提的多边形有限元法积分点数为原多边形单元经三角形化后在每个三角形中的数值积分点个数。 三角形中Gauss积分点分布 数值算例 2. 数值算例 ?小片试验： 考虑受单位分布力作用下的单向拉伸方板,板的边长为单位1,其边界条件如图所示，上部受均布拉力 = 1 .材料常数杨氏模量E = 100 kPa, 泊松比v =0. 3，假设为平面应力状态。 数值算例 小片的网格划分如图所示，对5种不同的网格采用Wachspress形函数进行小片试验,积分方案采用Gauss积分。 数值算例 下表为Wachspres