三维饱和地基人工黏弹性边界.pdfVIP

第13卷第2期 中国水利水电科学研究院学报 v。1．13N。．2 2015年4月 April，2015 Journal“ChinaInstituteofWaterResourcesand Research Hydropower 文章编号：1672—3031(2015)02一0128—08 三维饱和地基人工黏弹性边界 马怀发“2，杨 茹2 (1．中国水利水电科学研究院流域水循环模拟与调控国家重点实验室，北京100038 2．中国水利水电科学研究院工程抗震中心， 北京 100048) 摘要：本文针对中低频率工程振动问题，不计流体加速度，基于球面波动方程导出了无限大三维饱和地基人工黏 弹性边界及其流量边界条件。在此基础上给出了具有人工黏弹性边界的饱和土体动力固结问题虚位移原理。根据 有限元语言编写出算法脚本文件，并基于有限元自动生成系统(FEPG)，生成了计算源代码程序。通过在无限大 地基表面施加冲击荷载和突加荷载的数值算例，验证了所建立方程和程序的正确性。本文方法可以应用于三维饱 和地基的动力固结问题数值模拟。 关键词：半无限饱和地基；三维人工黏弹性边界；虚位移原理；孑L隙水压力；有效应力 中图分类号：TU470 文献标识码：A 1 研究背景 根据地震的传播机制，地震荷载是从远域地基波及到土体结构，而后产生地震响应。地震波由 远及近，再由近及远传播。在采用有限元法进行动力分析时，一般将半无限地基划分为邻近结构的 近域地基和其外围的远域地基。近域地基及土体结构构成了地震响应分析的研究对象，而切割近域 和远域地基的边界即形成了人工边界。人工边界分为全局边界和局部边界，全局边界虽然是对无限 介质的精确模拟，但其时空偶联的频域求解方法难以在有限元中实现，非线性问题不好考虑。局部 人工边界是时空解耦的，并且易于实现、计算量小，因而得到了广泛的研究和应用。目前，用的比 较多的是人工透射边界…、黏性边界‘21和在黏性边界基础上发展的黏弹性边界¨‘71。黏性边界因其概念 清楚、易于应用的特点，一度得到了广泛的研究和应用，但其存在低频稳定性问题，而且从一维推 广到多维有很大的误差。黏弹性边界是在黏性边界的基础上增加了弹簧，即通过沿人工边界设置一 系列线性弹簧和阻尼器组成的简单力学模型来吸收射向人工边界的波动能量，从而达到消除反射的 效果，并能模拟地基的弹性恢复能力，可以很好地模拟波动在无限地基中的真实传播过程。Deeks”1 基于柱面波给出了二维黏弹性人工边界的应力公式，刘晶波∞“1等人基于球面波推导了三维黏弹性人 工边界，并给出了边界的等效刚度和阻尼系数，马怀发等…在分析黏弹性边界方法的理论基础上建立 了统一的虚位移方程，郝明辉¨1等将其应用于拱坝计算。 强震作用下土坝、土体边坡及其地基会产生液化失稳问题，因此对其进行抗震安全评价理应将 土体看作由土骨架和孔隙水两相材料组成。动力分析不仅要给出土骨架的位移时程还要计算出孔隙 水压力的变化时程。与单向固体结构的人工黏弹性边界¨。71不同，对于饱和土体的人工边界，除了要 方程，针对实际土木工程的中低频率振动并且渗透系数比较小的情况，忽略第二类压缩波，将傍轴 近似应用于两相介质，提出了饱和介质动力方程u—p形式的黏性边界。Akiyoshi等¨驯进一步给出了 收稿日期：2014一ll一24 (K儿242) 作者简介：马怀发(1962一)，男，山东枣庄人，博士，教授级高级工程师，从事计算力学、水工结构抗震及混凝土细观力学分析 hr．com 研究。E—mail：mahf@iw 通讯作者：杨茹(1989一)，女，陕西，硕士生，从事工程数值模拟。E—mail：yrl002151615@163．com 一128— 万方数据 三维饱和地基人工黏弹性边界 马怀发 杨 茹 u一似和u—u形式的饱