三维随机渗流摄动有限元模拟研究.pdfVIP

2016年 10月 中国矿 山工程 第 45卷 第 5期 ChinaMineEngineering ·77 · 三维随机渗流摄动有限元模拟研究 Simulationresearchofrandom seepagefieldby3一D perturbationstochasticfiniteelementmethod 宋会彬，戴海旭(中国恩菲工程技术有限公司，北京100038) 摘 要 ：在非确定性渗流研究中，本文考虑渗透 系数 的空间随机性，基于摄动随机有限元法，推导相应 的渗流场随机分析 中渗 流响应量(水头)的随机响应公式，进而实现三维稳定渗流场随机有限元分析 。通过一个算例 ，将三维摄动随机有限元法计算 的随机渗流场与蒙特卡罗法计算的渗流场进行对比，验证 自编三维随机渗流摄动有限元程序的正确性和可行性。 关键词 ：摄动随机有限元 ；蒙特卡罗法；渗流计算 Abstract：In the study ofnon-deterministic seepage，considering the spatialrandomnessofpermeability coefficient，theformulasof seepageresponsevariable(waterleve1)instochasticseepageanalysisarededucedbasedonperturbationstochasticfiniteelement method，and then the stochastic analysis of3一D steady seepage field isperformed．Finally，a numericalexample isgiven and discussed．Thecorrectnessandfeasibility oftheself—made3一D seepageperturbationstochasticfiniteelementprogram isproved by comparingtheseepagefield thatcalculatedby3-D perturbationstochasticfiniteelementmethodandMonteCarlostochasticFEM． Keywords：perturbation stochasticfiniteelement；MonteCarlomethod；seepagecalculation 1 前言 随着近几十年来计算技术及可靠性理论的发展 ，岩 土工程数值分析中不确定性的考虑越来越受到重视。岩 土渗流计算作为岩土工程数值计算分析 的一个重要分 支，存在很多不确定 因素 ，因此 ，将渗透参数的随机性引 入到渗流计算 中是十分必要的。现有的对渗流场进行随 机分析中MonteCarlo法 是 比较常用的方法 ，该方法需 对渗透参数按照一定分布形式进行大量 的随机抽样 ，通 过对确定性渗流分析的多次重复求解 ，从而得到渗流场 的统计特征 。此方法计算量 比较大，限制了其在实际工 程 中的应用 。 摄动法是用于分析非线性问题 的强有力的方法 ，由 于理论明确 、计算量小 ，摄动有限元法在力学及其他工程