微分求积升阶谱有限元方法及其在薄板自由振动中应用.docVIP

微分求积升阶谱有限元方法及其在薄板自由振动中应用 　　摘要： 提出了一种基于Hermite插值的微分求积升阶谱有限元方法。单元在几何映射上采用了混合函数方法，而在形函数的构造上，单元边界上采用非均匀节点的Hermite插值基函数，单元内部升阶谱形函数的构造则采用雅克比正交多项式的张量积形式。将单元形函数与高斯-洛巴托积分法结合起来离散薄板的势能泛函从而得到相应的单元矩阵。提出的薄板单元在单元边界以及单元内部的节点配置完全自由，因而可以用于不同阶次的单元的连接。通过在薄板自由振动中的应用计算以及与精确解的对比，结果表明：提出的微分求积升阶谱有限元方法不仅计算精度高，而且收敛速度快，同时在阶次较高时仍然具有良好的数值稳定性。 　　关键词： 薄板自由振动； 微分求积方法； Hermite插值； 升阶谱有限元方法 　　中图分类号： O323； O241.82 文献标志码：A文章编号1004-4523（2018）02-0343-09 　　DOI：10.16385/j.cnki.issn.1004-4523.2018.02.019 　　引言 　　有限元方法自从20世纪60年代被正式提出以来，目前已经广泛应用于工程实际中，并对促进经济发展，解决工程科学问题做出了巨大贡献[1-3]。其中，h-型有限元方法以其技术简单、成熟，目前已发展成为有限元分析的主流方法。通过加密网格，该方法理论上可以得到足够高精度的计算结果。然而，在实际计算过程中，h-型有限元方法的收敛速度一般较慢，高精度的计算结果往往需要非常精细的网格才能得到，因而对CAD系统的前处理过程要求更加严格。伴随着h-型方法发展的还有p-型有限元方法，这种方法在计算过程中可以保持网格划分不变，通过提高单元的阶次便可得到高精度的计算结果。Babuka 等从理论分析和实际应用表明p-型有限元方法具有比h-型方法更快的收敛速度[3-5]。作为一种典型的p-方法，升阶谱有限元方法以其独特的优势在有限元方法发展的早期得到了广泛的研究[6-9]。其最显著的特点是低阶单元的基函数是高阶单元基函数的子集，这样在计算过程中中，低阶单元的结构矩阵是高阶单元结构矩阵的子矩阵，即所谓“嵌套”特性。进而在单元阶次升高时，前一步的计算结果仍然可以保留，从而降低重复分析的成本。基于这些优势，该方法一直被认为是未来实现自适应分析的主流技术。然而，由于历史和该方法自身的技术困难，升阶谱有限元方法的应用始终未能得到足够的推广[10]。 　　传统升阶谱有限元方法的基函数一般不具有插值特性，从而导致节点自由度不具有明确的物理意义，这样使得单元的组装以及边界条件的施加都变得困难。另外，随着单元阶次升高，高阶（几十阶或更高）多项式的计算带来的数值稳定性问题一直是一个技术难题[8]。最近，相关学者提出了微分求积升阶谱有限元方法，并将其应用于Mindlin板的自由振动问题中[11]。与传统升阶谱方法不同，该方法在单元边界上采用了非均匀节点的插值基函数，而在单元内部则仍采用升阶谱面函数。这样使得单元组装和边界条件施加的问题得到了解决。此外， 该方法采用了正交多项式的迭代公式进行计算，有效地克服了高阶多项式计算的数值动荡特性。 　　基于微分求积升阶谱有限元方法的上述优点，本文将把该方法推广到C1型单元的构造上。单元边界的形函数的构造通过一维高阶Hermite插值函数与三次Hermite插值函数混合得到，而单元内部面函数的构造则通过Jacobi正交多项式的张量积形式得到。其中，单元边界节点数目配置自由，因此可以适用于不同阶次的单元拼接以及自适应分析中局部自由度的调整。由于单元边界插值基函数的节点分布将对结构矩阵的数值特性具有重要的影响，本文提出了一类非均匀节点，即Gauss-Jacobi点，将有助改善单元的数值稳定性。此外，研究表明这类非均匀节点具有稳定的插值特性，因此可作为微分求积方法（强形式）的非均匀节点使用。计算表明，本文提出的C1型微分求?e升阶谱有限元方法具有精度高，收敛速度快等特点，在单元阶次较高时，仍然具有较好的数值稳定性。 　　第2期伍洋，等： 微分求积升阶谱有限元方法及其在薄板自由振动中的应用振 动 工 程 学 报第31卷1微分求积升阶谱薄板单元〖*2〗1.1一维Hermite插值基函数对给定定义在区间[-1，1]上的一维函数fξ，其Hermite插值近似函数为fξ≈H0ξf′-1+HN+1ξf′1+ 　　∑Nj=1Hjξfξj=∑N+1j=0Hjξj（1）式中ξj， j=1～N为插值节点，边界上的基函数为H0ξ=1-ξ22L1ξ 　　H1ξ=c1ξ+c21-ξ2L1ξ 　　HNξ=c3ξ+c41+ξ2LNξ 　　HN+1ξ=ξ2-12LNξ（2）内部节点对应的基函数为Hjξ=1-ξ21-ξ2jL