FEA-07-单元构造分析报告.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 有限元建模分析的一些实际考虑 有限元建模分析的一些实际考虑 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 数值积分方法 对于平面三角形和空间三棱锥(四面体)单元以及空间三棱柱(五面体)单元, 无法使用高斯积分, 此时应考虑Hammer积分, 其基本形式为: Hammer积分 数值积分方法 数值积分方法 数值积分方法 数值积分方法 数值积分方法 数值积分方法 数值积分方法 数值积分方法 数值积分方法 数值积分方法 　　当在计算中必须进行数值积分时，如何选择数值积分的阶次将直接影响计算的精度和计算的工作量。如果选择不当，甚至会导致计算的失败。 等参单元计算中数值积分阶次的选择 * * 保证积分精度对积分阶次的选择原则 一维问题:如果插值函数N(x)中的多项式阶数为p，微分算子L中导数的阶次是m，则(对于等参元假设|J|=const，即：单元形状为矩形或平行四边形时)有限元得到的被积分函数是2(p-m)次多项式。 因为，积分函数形式是：[?mN/ ?Xm]T[?mN/ ?Xm] 为保证原积分的精度，应选择高斯积分的阶次n=p-m+1，这时可以精确积分至：2(p–m+1)-1= 2(p–m)+1 次多项式，可以达到精确积分刚度矩阵的要求。 完全积分 减缩积分 选择积分(混合积分) * * 对于二维4结点双线性单元，假设|J|=const，即：单元形状为矩形或平行四边形时。刚度矩阵的被积函数ξ和η方向的最高方次为2，所以要达到精确积分，应选择高斯积分的阶次n=p-m+1=2。即：2×2阶高斯积分。如果单元|J|≠常数，则需要选取更多的3×3阶高斯积分。 对于二维8结点单元也可作类似的分析。同样得到结论是：为精确积分单元刚度矩阵，在|J|＝常数条件下，应采用3×3阶高斯积分。如果单元|J|≠常数，则需要选取更多的4×4阶高斯积分。 完全积分(Full Integration):高斯积分阶数等于被积函数所有项次精确积分所需要阶数的积分方案。 完全积分 * * 即对上述二种单元，为保证这部分被积函数积分的精度，应由完全多项式阶数p选择高斯积分的阶次为（n=p-m+1）只需要分别采用1×1和2×2的高斯积分。 二维4结点（ p =1）和8结点单元（ p =2）中，它们的插值函数中完全多项式阶数p分别等于1和2。微分算子L中导数的阶次是m=1。由完全多项式所产生的刚度矩阵中被积函数在ξ方向和η方向的最高方次，分别为0和2（在| J |=常数的条件下为2(p - m)） 减缩积分：由完全多项式最高方次p选择高斯积分。高斯积分阶数往往低于被积函数最高项次精确积分所需要阶数的积分方案。 减缩积分 * * 即：仍按n=p - m + 1 来确定积分方案。式中p是插值函数中完全多项式的方次，m是微分算子L中导数的阶次，二维单元和三维单元分别采用n×n和n×n×n高斯积分来进行单元刚度矩阵的计算。 实际计算表明：采用缩减积分往往可以取得较完全积分更好的精度。这是因为： （1）精确积分常常是由插值函数中非完全项的最高方次要求，而决定有限元精度的是完全多项式的方次。这些非完全的最高方次项往往不能提高精度，反而可能带来不好的影响。 取较低阶的高斯积分，使积分精度正好保证完全多项式方次的要求，而不包括更高次的非完全多项式的要求，其实质是相当用一种新的插值函数替代原来的插值函数，从而一定情况下改善了单元的精度。 减缩积分 * * 　　 （2）基于最小位能原来基础上建立的位移有限元，其解答具有下限性质。即有限元的计算模型具有较实际结构偏大的整体刚度。选取缩减积分方案将使有限元计算模型的刚度有所降低，因此可能有助于提高计算精度。 另外，这种缩减积分方案对于泛函中包含罚函数的情况也常常是必须的，用以保证和罚函数相应的矩阵的奇异性，否则将可能导致完全歪曲了的结果。 注意：在采用缩减积分方案计算K时需要检查矩阵K的秩，以保证矩阵K非奇异的必要条件恒被满足。 采用精确积分方案计算矩阵K时，矩阵K非奇异的必要条件是恒被满足的。 采用减缩积分时，系统的应变能并未被精确计算。有时可能导致某些除刚体运动而外的位移模式，即：所谓零能模式（hourglassing－沙漏现象）。实质是：在减缩积分点的应变可能正好等于零。 减缩积分 * * 完全积分与单元锁定（Locking）现象 * * * * * * * * * * * * * 采用完全积分,线性平面单元要有2×2个积分