流形方法的理论与应用研究

【摘要】：在有限覆盖基础上新发展的流形方法统一解决了有限元、非连续变形分析和解析法的计算问题,是一种具有广泛应用前景的数值方法。本文的主要工作是分别采用标准的矩形网格、传统的有限单元网格和离散的扇形来构成流形方法的覆盖系统,并将其应用于二维连续体的弹性静力分析和张裂纹扩展过程的数值模拟。 文中首先选取标准的矩形作为流形方法的基本覆盖,详细研究了流形方法里边界条件的施加,单元权函数的选取,连续的高精度应力场的计算,流形单元的数值积分方法以及各阶覆盖位移函数对流形方法精度的影响等关键问题。由于对各种复杂形状的结构体均可使用统一的规则矩形网格,本文方法的实施、使用及其与CAD 技术的一体化等都十分简单、方便。 针对提出的理论和方法,论文首次把面向对象的思想引入流形方法的程序设计中,将流形方法的有限覆盖系统抽象为一些独立的数据类,给出了类的描述和它们的实现方法,并用树状结构对这些数据类进行管理。对任意形状的复杂结构体,用流形方法的数据类实现了流形单元有限覆盖系统的全自动生成。 接着,文中采用围线积分法和流形方法的高阶位移函数来计算得到比较准确的混合型裂纹尖端的应力强度因子,然后基于线弹性断裂力学的裂纹扩展准则,用流形方法实现了平面Ⅰ-Ⅱ混合型张裂纹扩展的数值模拟,为数值方法研究裂纹的演化行为提供了一条简单而有效的途径。 随后本文结合流形方法的覆盖理论和滑动最小二乘法,尝试提出了一种改进的无单元方法。它摒弃了单元和网格,采用圆形或扇形的覆盖图形并直接用结构域内的离散点来简化整体近似位移函数的构造,可以看作是流形方法实现的另一种特殊形式。该方法形式简单、计算时间少,比以往的一些无单元法具有更高的效率,并且更易于编程实现。