比例边界有限元分析侧边界上施加不连续荷载的问题.PDFVIP

水 利 学 报 2017年2月 SHUILI XUEBAO 第48卷 第2期 文章编号：0559-9350（2017）02-0246-06 比例边界有限元分析侧边界上施加不连续荷载的问题 1，2 1，2 1，2 1，2 庞 林 ，林 皋 ，李建波 ，薛冰寒 （1. 大连理工大学 水利工程学院，辽宁 大连 116024； 2. 大连理工大学 海岸与近海工程国家重点实验室，辽宁 大连 116024） 摘要：比例边界有限元方法是一种半解析的数值计算方法。与边界元一样，比例边界有限元可以将求解问题的维 度降低一维。由于比例边界有限元对通过相似中心的侧面边界不进行离散，在实际计算时，作用在侧面边界上的 荷载分布需用径向坐标的幂函数进行逼近。然而当指定的侧面边界荷载不连续时，就不可以按照这种形式近似表 达。这时就需要在不连续点处划分为较多的多边形单元，从而导致前处理中的网格划分十分复杂。本文对比例边 界有限元侧面边界上作用不连续荷载的情况提出了一种简便、有效的处理方法。当侧边界荷载分布改变时无需调 整网格划分，使计算简化。本方法拓宽了比例边界有限元的应用范围，算例显示了所提方法的有效性。 关键词：比例边界有限元；不连续荷载；侧面边界；裂纹面荷载；断裂力学 中图分类号：TV31 文献标识码：A doi ：10.13243/ki.slxb 1 研究背景 对于重要水工结构的设计和校核通常需要补充有限元分析。比例边界有限元（ScaledBoundaryFi⁃ ［1］ niteElementMethod，SBFEM）是一种半解析数值方法 ，与传统的有限元方法相比，因其只需要进行 边界离散且具有径向解析的特点，能够以较小的前处理工作和计算量达到理想的精度，适合于分析 ［2-4］ ［5-6］ 水工结构与无限地基或库水的动力相互作用 以及裂纹扩展 等问题。近年来，针对其应用范围和 方法改进研究成为一个热点，应用范围扩展至渗流、液面晃动、波浪、高压电路中工频电场以及水 ［7-11］ 下结构的流固耦合等多个领域 。通过将比例边界有限元方法与有限元方法、扩展有限元方法、无 网格方法以及等几何方法［5，12-15］耦合的研究，以及处理移动相似中心和基于数码摄影图像自动划分网 格［16- 17］等问题的研究，进一步增强了比例边界有限元分析问题的能力。 工程实际中荷载形式十分复杂，在数值分析中常将其简化成常数或线性分布形式。如混凝土重 力坝设计规范要求叠加坝前水压力和泥沙压力，此时大坝的上游面受到的荷载形式为分段线性荷 载。此外，在使用比例边界有限元模拟裂纹扩展中考虑裂纹面混凝土残余的黏着力以及研究水力劈 裂问题［18-20］时，也常常遇到这种侧面边界存在不连续分布的荷载形式。然而，由于比例边界有限元只 离散外边界且通过相似中心的侧边界不离散，实际计算时，需要将侧边界上的荷载近似按幂级数展 ［21］ 开处理 。显见，这样的近似无法模拟不连续形式的荷载。类似于有限元方法，需要在荷载不连续 ［22］ 的位置人为地划分为较多的多边形超单元 。这使得前处理变得复杂且在一定程度上丧失了半解析 的优势。 本文将针对这一问题提出解决方法，比例边界有限元在每个分段区域内，位移解和内力解在径 向的形式是包含积分常数的幂函数。由于分段增加了未知积分常数的个数，因而需要补充兼容性条 收稿日期：2016-05-10；网络出版日期：2016-10-20 网络出版地址：http：///kcms/detail/11.1882.TV1558.001.html 基金项目：国家自然科学基金创新研究群体基金项目；国家自然科学基金项目（511