离散元2边界元动力耦合模型-中国科技论文在线.PDF

水 　　利 　　学 　　报 　2001 年 1 月 SHU IL I 　　XU EBAO 第 1 期 文章编号 : (2001) 离散元边界元动力耦合模型 1 1 1 金 　峰 , 贾伟伟 , 王光纶 ( 1 清华大学水利 水电工程系 , 北京 　100084) 摘 　要 : 本文提出了一种二维变形体离散元与时域边界元的耦合模型 , 这一模型可以将非连续体的模拟与无限 域的模拟统一在一个模型中 , 可用于在地震波动输入条件下 , 考虑辐射阻尼的岩体边坡或地下结构等的动力稳 定和变形分析 , 拓宽了离散元动力分析的领域. 算例分析表明本耦合分析模型具有较高的精度. 关键词 : 离散元 ; 边界元 ; 耦合 中图分类号 : O344 　　　　　　文献标识码 : A 离散元法是一种模拟离散介质的计算方法 , 自Cundall 在 70 年代提出以来 , 在岩石力学 、土力 学 、结构分析等领域的数值模拟中得到广泛应用 , 是一种新兴的非连续体分析方法[1 ,2 ,3 ] . 在动力分 析中 , 许多分析均表明辐射阻尼对分析结果影响很大 , 应该充分重视[4 ] , 动力边界元方法由于在其 基本解中包含了无限远处的辐射条件 , 在处理辐射阻尼的影响时十分方便 , 是分析地下结构动力响应 的一种有力的工具[5 ] . 动力边界元法可以分为频域边界元法[5 ] 与时域边界元法[6 ,7 ] , 前者主要适用于 线性问题. 是发展较早 、相对成熟的动力边界元方法. 时域边界元方法直接在时域内求解 , 适合解决 非线性问题 , 根据求解的问题具有非线性的特点 , 本文选择动力时域边界元方法与可变形体离散元耦 合 , 提出了一个二维时域动力边界元 ———可变形体离散元耦合模型 , 充分发挥离散元与边界元的优 点 , 将非连续体的模拟与无限介质辐射阻尼的模拟统一到一个模型中 , 为地下结构和岩质边坡的抗震 稳定分析提供了全新的手段. 1 　离散元边界元耦合模型 11 　二维可变形体离散元原理 　可变形体离散元法将模拟的 计算区域看作是若干可变形块体的组合 , 见图 1. 这些块体可 以任意平移 、旋转 , 块体之间的相互作用力 , 用法向和切向弹 簧表示 , 称为接触 , 接触力的大小由块体的相对位置决定. 每 个块体又划分为三角形差分网格以模拟变形 , 每个差分三角形 的顶点称为节点 , 其应变假定为常数 , 可 以由节点的位置确 定. 因为一旦节点位置确定 , 便可以得到差分三角形乃至整个 块体的变形和应力 , 从而得到所有的响应历程 , 因此 , 所有的 图 1 　离散体系示意 计算将围绕节点进行 , 边界条件也可以通过给定边界节点位移 或节点力来实现. 离散元计算采用显式步进的方法 , 首先将计算的过程分为若干等长的时步 , 静力分析作为动力分 收稿 日期 : 基金项 目: “九五”国家攻关项 目. 作者简介 : 金峰 ( 1966 - ) , 男 , 贵州遵义人 , 教授 、博士生导师 , 主要从事结构动力分析研究.