biot-cosserat连续体-离散颗粒集合体模型的非饱和土-应用力学学报.pdfVIP

第 30 卷 第 3 期 应 用 力 学 学 报 Vol.30 No.3 2013 年 6 月 CHINESE JOURNAL OF APPLIED MECHANICS Jun. 2013 文章编号：1000- 4939(2013) 03-0297-07 Biot-Cosserat 连续体-离散颗粒集合体模型的 非饱和土连接尺度方法 万柯 李锡夔 (大连理工大学工业装备结构分析国家重点实验室 116024 大连) 摘要：耦合了非饱和多孔多相介质有限元模型和颗粒介质离散元(DEM)模型，提出了以宏、细观 尺度分别耦合 Biot-Cosserat 连续体模型和离散颗粒集合体模型的连接尺度方法来分析非饱和含液 颗粒材料的力学渗流耦合问题。根据被动空气压力假定和对空间离散孔隙水质量守恒方程的约化， 从非饱和土有限元控制方程的基本未知量中消去了孔隙水压力，而将其取作有限元积分点上定义 取值的内状态变量，进而建立了节点未知量仅包含固相线位移和转角的非饱和 Cosserat 多孔连续 体约化有限元数值模型。基于连接尺度方法(BSM)宏、细观尺度数值过程的解耦计算的特点，对 宏、细观两尺度数值模型的时域积分分别采用隐式 Newmark 方法和显式中心差分法，且取不同时 间步长以提高计算效率。与全域采用DEM 的精细分析方法相比，本文 BSM 在保证计算精度的前 提下可大幅节省计算时间。在不考虑湿化效应的边坡稳定算例中，在得到类似计算精度条件下它 比全域采用DEM 节省计算时间高达 86.65%。二维边坡稳定算例结果验证了本文连接尺度方法的 有效性，以及在揭示含液颗粒结构细观破坏机理上的优点。数值算例结果显示，边坡承载能力因 降雨大幅下降约 50% ，这表明本文发展的计及伴随湿化过程的颗粒材料结构中饱和度及吸力分布 演变及其对结构破坏失效影响的非饱和颗粒材料多尺度计算模型是很有必要的。 关键词：非饱和土；多孔多相介质；多尺度；Cosserat 连续体；离散单元法 中图分类号：TU43 文献标识码：A DOI ：10.11776/cjam.30.03.B063 Biot 理论模型。基于饱和土 Biot 理论，文献[4]采用了 1 引 言 被动空气压力假定，即假定非饱和土域内孔隙气体压 力均等于周围大气压力，发展了非饱和土数值模型。 文献[1]首先建立了固液两相耦合作用的饱和土 文献[5]在 ZienKiewicz 两相模型的基础上，舍弃了被 的动力控制方程和Biot 理论框架。文献[2]从土体骨架 动空气压力，针对含不混溶流体的变形多孔多相非饱 的平衡方程和流相的动量守恒和质量守恒方程出发， 和土建立了多相流与变形多孔介质相互作用的三相 基于 Terzaghi 有效应力原理[3]，提出了饱和土的广义 数学模型和有限元数值分析模型。 基金项目：国家自然科学基金 ；国家重点基础研究发展计划资助项目(2010CB731502) 收稿日期：2012-05-25 修回日期：2012-12-31 第一作者简介：万柯，男，1982 年生，大连理工大学工程力学系，博士生；研究方向—— 固体力学、计算力学。 通讯作者：李锡夔，E-mail ：xikuili@ 298 应 用 力 学 学 报 第 30 卷 近几十年来 Cosserat 连续体理论得到了快速发 u u +u′ (1) 展和应用，文献[6]总结和完善了 Cosserat 连续体的 式中细尺度位移