用三维映射无限元模拟半空间弹性地基.pdfVIP

维普资讯 用三维映射无限元模拟半空间弹性地基 燕 柳 斌 (广西大学土木系) 利用有限元对半空间弹性地基进行位移 与精确解及衰减无限元 的比较表明，这样 及应力分析时，往往是截取一定的深度和宽 的模拟是成功的。三维映射无限元在精度上 度区域，并在边界上施加人为的边界条件， 是可靠的，并且有很高的经济性，值得在工 将无限区域简化为有限区域。但截取多大范 程推广使用。 围才能得到较为符合实际的解答，一直是个 一 、 三维映射无限竞的 不明确的问题 。特别是在坝体应力分析中， 基本原理与公式 坝体与地基的相互作用属半空间弹性介质的 问题 。地基的弹性对坝体的作用是不能忽视 关于一维及二维映射无限元 的基 本 原 的。有关文献认为计算模型包括 1— 5倍坝 理 ，详见有关文献 、 、一 。由二维映 高的地基范围才能给出满意的结果 但这样 射无 限元 的概念很容易推广到三堆情况。本 不仅增椰了数据准备工作量 ， (尤其是对于 文采用 12节点的无限元来模拟地基远场。为 三维问题 ，数据准备工作量及计算机时的耗 便于与有限元在耦合处的变形协调，采用20 费就更为可观) ．而且往往 由于计算机容量 节点的等参单元与其耦合。图 1为一映射无 使得解题规模受到限制。 自无限元模型提出 限元的实际单元。节点 1— 8与 有 限元 耦 后 ．立即运用于这类无跟 连 续 介 质 的问题 合，节点9，10，l1，12为 无 限元 中间节 中，并取得了较高精度和较为合理 的计算结 点，其余 4个节点在无穷远处。 果 。 构造一组映射函数将图 1 (B)所示的 目前无限元按其方法分为衰减函数无限 实际无限区域单元从几何上映射为一有跟的 元和映射无限元。前者利 用Largr~ge插值函 母单元。无限元的节点设置如图 1 (b)所 数与反映位移衰减特征的衰减函数之乘积来 示，单元沿E方向为无限远。 构造形函数。但衰减函数的选用须根据不 同 整体坐标x—y— z与局部坐 标 §一 一￡ 问题来选用，靠无限元 中间节点的调整来控 的变换关系为； 制无限元位移 的衰减程度 ；而且在采用高斯 x妻MiX Y ： M Y 积分形成单刚时存在广义积分收簋性的判别 问题 。映射无限元则采用映射技术，通过一 组映射函数将无限区域从几何上映射为一有 Z兰MYi 限区域，然而对这一有限区域按普通等参元 其 中，Mi为硖射函数，取- 考虑即可求出无限元的单元剐度矩阵。 M1= (1一1)[1·11一l一 + #一}+ 专一2 本文在文献 )、 )的基础 上，将这一 +211)／C2(’1一￡)]， 模型扩展到三维情况，利用l2节点的映射无 Mt= (i一§‘)(i一11)／(i一￡)’ 限元与20节点等参有限元耦合来模拟半空间 M。= (1+})(1