关于等参单元的高斯积分.pdfVIP

第6期 罗健豪等：关于等参单元的高斯积分 73 关于等参单元的高斯积分 罗健豪罗健文 (澳门大学机电工程系，澳门) 摘要讨论了等参单元的高斯积分．指出采用高阶高斯积分 算出的刚度矩阵虽更准确．但却不一定使有限元解更接近真 实解．从有限元法的本质看，稍为低阶的高斯积分有时反而 可能给出更好的解答．对该现象作出解释并提出了积分阶数 选择的建议． 图1 关键词等参单元，高斯积分，积分阶数 当边界条件取为梁的轴线在固定端处既无平移也无转角 时，本例的材料力学解为 1引 言 对于等参单元，其平面问题的单元刚度矩阵 PL。 △2面2 [Ke]=fl—fl一[B]T[。][B】·t·IJId∈d叩 444．8×2．543 i赢i石百忑厕≈ 的计算一般须应用数值积分，其中最常用的是高斯求积法． 关于高斯积分的阶数Ⅳ，文献[1]指出： “用数值积分代替 精确积分，无疑计算时会产生误差．为尽量减少这一误差， 下面分别用四节点及八节点等参单元和不同阶数的高斯 希望采用尽可能高的数值积分阶次”．文献[1】亦同时指出： 积分求解本例，检查积分阶数对解的影响． 鉴于Ⅳ很大时计算费时，建议在二维情况下，四节点单元的 对于四节点单元，网格取为图2中的A，B．结果示于表 Ⅳ取为2，八节点单元的Ⅳ取为3较好．文献[2】及其它一 1． 些教材对Ⅳ的取法也有类似的论述．文献[3，4]则只列出高 斯积分表，未有讨论Ⅳ的取法．文献[1,-一4】对阶数Ⅳ的论 述在现行的有限元教材中是颇具代表性的． A B 笔者发现学生对Ⅳ的取法，与上面的文献一样，都大致 分成两类：一部分学生对四节点及八节点单元取Ⅳ分别为 2和3，并认为得出之解比利用较大的Ⅳ算出的稍为差，但 C D 此时运算时间不致太长．另一部分学生则不论什么单元，一 律取十分大的Ⅳ，认为这样算出的[Ke】将更精确．因此有 E工j] 限元解也更好．第一类同学反映他们主要是根据书本知识． E 第二类同学则是由于所用的书上未有讨论Ⅳ的取法．因此 图2本文所用的各种网格 他们直觉地认为有限元解的精度与Ⅳ成正比．虽然两类同 学对Ⅳ的取法不同，但均认为：小的Ⅳ不如大的Ⅳ．因为 表1四节点等参单元给出之梁端挠度(mm) 大的Ⅳ可更加准确地算出[Ke】．因此在单元种类及划分法 一1．01591 (理论解为： mm) 选定后，有限元解精度将正比于高斯积分的阶数Ⅳ． 笔者认为：从数学角度看，越高阶的高斯积分当然能把 单元刚度矩阵[Ke】中的每一项算得更准确．从而使得总刚 度矩阵[困也更准确．但是，这个数学上更准确的[吲所给 出的有