水工结构的三维阶谱有限元分析.doc

水工结构的三维阶谱有限元分析 程 昭 陈胜宏 (武汉水利电力大学水电学院) 摘 要 根据p型有限单元法的阶谱特点，详细论述了三维阶谱单元法的基本的分析过程和具体的实现路径，包括基函数的构造、边界约束条件的处理、刚度矩阵和荷载列阵的形成、提高数值积分效率的途径等。引入了三维阶谱单元法的虚结点和广义结点的概念，并把三维阶谱单元法应用于水工结构计算。 关键词 有限元法，阶谱单元，三维分析，p型，水工结构。 本文于1999年5月24日收到，系国家自然科学基金资助项目. ??? 有限元法是研究工程结构问题最为广泛的数值方法。从逼近真实解的途径分类，有限元法可以分为：(1)h型，即传统的有限元法，通过减小单元尺寸h来提高有限元解的精度；(2)p型，通过增加基函数多项式的阶数p来提高有限元解的精度；(3)hp型，综合了(1)(2)两种方法。文献[1]给出了p型有限元法的理论分析，指出p型有限元的收敛速度比h型快，而对于奇点问题，至少是h型的两倍。文献[2]阐述了p型有限元法中的阶谱概念及其优点。基函数具有阶谱特性的p型有限元法，称为阶谱有限单元法或阶谱单元法。当前,p型有限元法一般都采用阶谱单元，所以在不作特别说明的情况下，p型有限元法和阶谱单元法是同一个概念。文献[3]对一维阶谱单元法作了较为详细的研究。目前对阶谱有限单元法的研究主要是针对二维问题，且偏重于理论分析，对工程计算的具体过程如何实现论述较少。比如，约束如何体现，如何提高数值积分效率等等。本文详细地给出了三维阶谱单元法的基本的分析过程和实现路径，包括基函数的构造、边界约束条件的处理、刚度矩阵和荷载列阵的形成、提高数值积分效率的途径等。引入了三维阶谱单元法中的虚结点和广义结点的概念，并应用于水工结构工程计算，取得了满意的效果。 1 三维阶谱单元的基函数 1.1 基函数的特点 以Hpi表示单元尺寸不变时pi阶阶谱单元逼近空间，阶谱的要领即低阶单元逼近空间是高阶单元逼近空间的一个子集： (1) 因此，低阶单元刚度矩阵是高阶单元刚度矩阵的子块，当为提高精度而升阶时，可继续利用已经计算出来的低阶单元刚度矩阵[3,4]。 ??? 三维阶谱单元的基函数包括点基函数、棱基函数、面基函数和体基函数。点基函数满足在本点值为1，在其余各点值为0；棱基函数满足在本棱的端点处值为0，在其余各棱上的值为0；面基函数满足在本面的棱边上的值为0，在其余各面上的值为0；体基函数满足在各面上的值为0，体基函数有时也称为内部基函数。 ??? 本文以六面体单元为例进行分析，四面体单元和楔形单元[5]分析过程类似。 1.2 基函数的形式点基函数： (2) 其中，ξ0,η0,ζ0为i点在母单元坐标系中的坐标，ξ0=(-1)i,η0=(-1)[i/2+0.5],ζ0=(-1)[i/4+0.75]。这和八结点等参单元完全相同。 ??? 棱基函数(p≥2): (3) 其中，,Lp是p阶单位Legendre多项式。 ??? 面基函数(p≥4)： ( 其中，i,j≥2,i+j=p. ??? 体基函数(p≥6): (5) 其中，i,j,k≥2,i+j+k=p. ??? 以上相当于Serendity族基函数，共有基函数的个数： (6) 1—6，再增加一个体结点作为第105个结点，相应的基函数为6B(2，2，2)；如此类推。基函数可以统一记为φi，表示第i个广义结点的基函数，其中1≤i≤fe(p). (7) 坐标的插值仍同等参单元法： (8) 于是 (9) (10) 将结构看作点数、棱数、面数和体数分别为Nn，Ne，Nf，N的一个‘大的单元’，将Nn，Ne，Nf，N分别对应3.1节中单元的点数8，棱数12，面数6和体数1，以3.1节相同的方法编排整体结点，并去掉约束所对应的行和列，能使总体刚度矩阵也成为阶谱形式。 3.3 荷载移植 荷载移植与传统的有限元法类似： (11) 需要注意的是，由于虚结点的存在，移植后的结点荷载和与原荷载不存在平衡关系，只是实结点的移植荷载和与原荷载平衡。这并不妨碍问题的解决。经3.2节方法对整体自由度进行整和后，整体荷载列阵｛F｝也是阶谱形式。 4 数值积分 ??? 当阶数升高时，基函数将是高阶多项式，为保证精度，需要更多的高斯积分点，数值积分的运算量将随着阶数升高成倍增加，对三维问题尤为突出。利用阶谱单元的特点，减少运算量的措施有： ??? (1)充分利用[B]的稀疏性。对于弹性问题，还可以利用[D]的稀疏性，将6×6的矩阵[Bi]T[D][Bj]的各个元素先算出来，再对各个元素作数值积分。 ??? (2)单元刚度矩阵的不同元素采用不同阶高斯积分来计算。实际上，如果对高阶阶谱单元，各个元素全部都用高阶高斯积分，计算