堤防渗流场参数敏感性三维随机有限元分析.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 堤防渗流场参数敏感性三维随机有限元分析 王亚军 张我华 王沙义 浙江大学软弱土与环境土工教育部重点实验室 浙江杭州 摘要 采用随机有限元分析方法 研究了长江荆南干堤的三维各向异性非均质稳定随机渗流场中的随机参数和随 机边界条件对渗流场模拟结果的影响 并进行敏感性分析 通过统计模拟和假设检验 建立并验证了该干堤渗流 场内服从独立正态分布的渗透系数张量 讨论了堤防土层的地质结构特性 堤防渗透变形及破坏的特征 堤防渗透 参数的敏感性 研究了上 下游水头随机边界条件 垂直防渗边界条件以及减压 导渗 沟随机特性的敏感性 实现 了对渗流场的全面随机模拟和分析 得到的结论通过统计检验 并根据实测工程数据对照证明是可靠的 关键词 三维随机有限元 非均质随机渗流场 各向异性随机渗透系数 随机边界条件 随机敏感性分析 中图分类号 文献标识码 土石坝的渗流问题包含有许多随机因素 由于受这些随机因素的制约和影响 从而使土石坝渗流场 的分析具有很大程度的不确定性 造成渗流场数学模型具有不确定性的原因很多 主要有以下两个方 面原因 模型中出现的各种参数的不确定性 在模型定解条件 初始条件 边界条件和 载荷 中 也包含着不同程度的不确定性 由于模型中包含着这些不确定性 所以由模型得到的渗流场解 如水 头 压力 流速等 应当看作具有一定概率分布的随机函数的随机过程 把随机渗流场模型应用于预测 和管理土石坝的渗流稳定性时 一个很重要的问题是怎样把模型的不确定性参数定量化 以保证模型使 用的可靠性 一般来说 若渗流场定解问题中控制方程的系数或定解条件包含随机因素时 应当首先对 这些随机因素进行统计假设检验 并对这些随机参数对干堤典型段渗流场随机特性的影响的参数敏感 性进行较为深入的探讨 才能提高随机渗流场模型使用的可靠性 等 人和 等 人对此进 行了研究 国内这方面的研究工作也已开始 本文基于渗流场的随机特性 对所建立的长江荆南干 堤三维各向异性非均质稳定随机渗流场模型进行参数敏感性的分析和讨论 渗流场的随机有限元 确定性渗流场问题可以归结为求解如下泛函的极值 式中 为水头势函数 为关于 的泛函 为流场体积 为直角坐标系统 下各向异 性主方向渗透系数 对于如上问题在三维流场内进行离散即得求和泛函为 式中 为单元数 为单元体积 引入插值函数 为单元节点数 后得到传导矩阵 为 收稿日期 作者简介 王亚军 男 山西平定人 博士生 主要从事岩土工程力学的研究 其中 就随机渗流场而言 若不考虑单元的变异则随机特性主要来自于 的变异 笔者分析并验证了 中主元 呈标准正态分布 而模型中传导矩阵线性依赖于高斯随机过程 模型多种随机性的检验 土性参数的随机性 荆南干堤座落在第四系冲 湖积 层 且主要为全新统土层上 其 段截面材料的地质构造形式如图 所示 通过勘探发现 堤基土类主要为 人工填土 粉 质壤土 粉细砂 粉质黏土等 局部见有少量基岩 图 段截面材料 单位 因此 土层按其工程地质特性分为 类 对应的渗透 系数为 对于长江荆南干堤三维各向异性非均质稳定随机渗流场模型 主要的土性参数是渗透 系数张量 按土的渗透性分成两大类 其一为砂性土类 具有中等透水或中弱透水性 其二为黏性 土类 具有弱透水 微透水性 在荆南干堤段内 长江河道曲折蜿蜒 摆动频繁 并多次出现改道 沉积环境在不断变化 因而组成 堤基的土层各地差异较大 堤基地质结构复杂多变 根据堤基工程地质条件 堤基砂性土和黏性土呈组 合型式的结构 结合防渗工程特点 将堤基地质结构分为 个主要类型 单一