基于统计模型的裂隙岩体渗流场的并行数值模拟-岩石力学与工程学报.PDFVIP

第27 卷 第4 期 岩石力学与工程学报 Vol.27 No.4 2008 年 4 月 Chinese Journal of Rock Mechanics and Engineering April ，2008 基于统计模型的裂隙岩体渗流场的并行数值模拟 1 1 2 刘耀儒 ，杨 强 ，覃振朝 (1. 清华大学 水沙科学与水利水电工程国家重点实验室，北京 100084 ；2. 福建水口发电有限公司，福建 福州 350002) 摘要：为模拟天然岩体中裂隙对渗流的影响，基于裂隙的空间形态及分布的统计特性，应用裂隙网络随机模拟技 术，采用 Matlab 程序生成二维裂隙网络。由此裂隙网络生成包含裂隙的有限元计算网格。考虑到渗流场数值模拟 分析时步多、计算量大的问题，采用基于 element-by-element 策略的有限元并行计算方法进行数值模拟。该方法不 集成整体刚度矩阵，可有效节省存储量，同时可降低共轭梯度法中的迭代次数，提高收敛性。编制相应的模拟分 析并行程序，并在分布存储的并行机上实现。分别对包含 1，2 组裂隙的岩体进行渗流场的数值模拟。计算结果表 明，该方法能够较真实地反映裂隙岩体的实际状态，对于模拟岩体中的裂隙是很有效的；同时，并行计算可有效 降低数值模拟的计算时间。 关键词：数值分析；渗流；裂隙岩体；并行计算 中图分类号：O 241 文献标识码：A 文章编号：1000 –6915(2008)04 –0736 –07 PARALLEL NUMERICAL SIMULATION OF SEEPAGE IN FRACTURED ROCK MASS BASED ON STATISTIC MODEL 1 1 2 LIU Yaoru ，YANG Qiang ，QIN Zhenchao (1. State Key Laboratory of Hydroscience and Engineering ，Tsinghua University，Beijing 100084，China； 2. Fujian Shuikou Power Co.，Ltd. ，Fuzhou ，Fujian 350002，China) Abstract ：To simulate the fractures in natural rock，the fracture network has been formed with Matlab based on statistical parameters of fracture shapes and random distribution of fractures. Then finite element mesh has been generated with fracture network. Considering enormous time steps and huge workload in the numerical analysis ， the element-by-element(EBE) finite element parallel computation method is adopted for numerical simulation. This algorithm does not generate global stiffness matrix and only stores element stiffness matrices，which