[硕士]非线性自然单元法的程序实现及其应用研究.pdf

上海交通大学硕士学位论文 中文摘要 非线性自然单元法的程序实现及其应用研究 摘 要 自然单元法是一种新兴的用于求解偏微分方程的无网格数值计算 方法，其本质是基于自然相邻插值的 Galerkin 方法。它的基本流程是利 M O 用需求解点的自然相邻结点和二阶 Voronoi 胞来构造场变量插值函数， C . 在局部的Delaunay 子域上采用标准的Galerkin 方法建立整理求解的控制 G 方程。对于本文中的偏微分方程，由于已知其变分原理，可以等效地采 N O 用变分原理进行推导。 L U 由于自然单元法采用三角形或四面体单元作为背景积分网格，在网 H 格划分以及自适应加密等方面相对四边形或六面体单元更为简单。而且 Z . 同常规有限元法相比，自然单元法的计算精度大大高于三角形或四面体 W 单元，与四边形或六面体有限元法相当。由于自然相邻插值函数在凸域 W W 边界上的相邻结点间是严格线性的，因此自然单元法在边界面的处理上 也相当简单。 本文阐述了自然单元法的基本原理与算法，包括 Delaunay 三角化和 网 Voronoi 图的概念、自然相邻插值函数及其导数的建立过程等；结合几 龙 何大变形基本理论和增量弹塑性基本理论，研究了利用自然单元法解决 筑 几何非线性以及几何材料耦合非线性问题。另外结合 Biot 固结理论，研 究了利用自然单元法解决固结问题。设计了二维非线性自然单元法算法 和二维固结自然单元法算法，对算法进行了程序实现，利用 MATLAB 编制了相应的计算程序。作为验证，本文的算例采用了数值计算软件的 经典考题：均布荷载下的悬臂梁，单级等速施加的竖向均布荷载下单层 I 上海交通大学硕士学位论文 中文摘要 均质地基压缩的问题。针对每类问题，进行了自然单元法计算、相同网 格划分的有限元法计算、解析解或高精度有限元计算。各个算例均验证 了本文自然单元算法程序的正确性，并显示了自然单元法在精度方面的 优势。 关键词：自然单元法，自然相邻插值，大变形，弹塑性，几何材料耦合 非线性，Biot 固结 M