岩土工程中数值流形方法的应用与研究

【摘要】： 数值流形方法用于计算结构或材料的位移和变形,能对连续变形与非连续变形问题进行统一求解,对问题具有较高的求解精度,是目前岩土工程数值分析方法的一个研究热点。 本文针对数值流形方法理论体系及其在岩土工程中的应用,并结合一些学者已有的研究成果,开展了对数值流形方法部分研究工作。本论文的研究工作及取得的成果主要有以下几个方面: (1)研究了覆盖位移函数对刚度矩阵形成的影响,提出了将覆盖位移函数的基本级数函数基中x、y由材料的绝对坐标改进为以流形单元的形心坐标为原点的相对坐标的建议来改善刚度矩阵的局部性态,以利于求解计算。同时指出了选择合理的坐标原点以便充分发挥覆盖位移函数的全阶项对总平衡方程组成时刚度矩阵和荷载矩阵的贡献,使求解精度得到提高。 (2)基于广义变分原理的理论体系,利用罚函数法,推导了梁板单元分析的修正泛函;构造了梁板流形单元及相应的流形单元格式(位移函数、应变矩阵和刚度矩阵等)。 (3)针对工程中广泛应用中厚板的弯曲分析问题,以数值流形方法的基本思想为基础,建立了Winkler弹性地基上Mindlin板的数值流形方法;对工程中分析弹性地基上的中厚板弯曲问题具有一定的参考意义。 (4)利用数值流形方法对隧道工程中衬砌结构进行了尝试性的分析,提出了结合规则网格和有限网格并行、高低阶覆盖位移函数混合计算的思想,以达到前处理和程序计算、计算效率和求解精度的协调。 (5)通过对目前数值流形方法中采用的几种本构模型的研究,提出了适用于非线性分析的数值流形方法。文中本构模型采用非线性弹性模型(Duncan-Chang模型和K-G模型),为数值流形方法进行岩土工程非线性分析提供了有益的数值分析方法。 (6)在前人研究的基础上,构造了六面体有限覆盖的三维流形单元,推导了三维流形单元覆盖位移函数、应变矩阵、刚度矩阵、荷载矩阵等表达式。通过算例分析及与解析解和ANASYS计算结果的比较,表明该思路是有效可行的。 (7)自编了较大型的面向对象计算机程序,该程序采用C++语言在VC++平台中编制,能对数值流形方法的一般问题(如文中所涉及的相关问题等)进行计算,具有较好的程序前后处理功能,并通过算例验证了程序是有效和合理的。