平面稳态渗流的随机有限元分析程序  

Programming of Stochastic Finite Element Method in Plant Steady Seepage

作　　者：王珍[1] 刘卫群[1] 彭晓庆[1] 

机构地区：[1]中国矿业大学力学与建筑工程学院,徐州221008

出　　处：《力学季刊》2011年第2期203-210,共8页Chinese Quarterly of Mechanics

基　　金：国家自然科学基金(50774083;41074040);国家"973计划"项目子课题(2009CB219605);教育部"新世纪优秀人才支持计划"课题(NCET-07-0803)

摘　　要：岩体渗流模型中存在着许多不确定的因素,将渗透参数的随机性引入渗流有限元分析计算中会使分析结果更客观、更明确。文章首先基于小参数摄动随机有限元方法,推导出相应的渗流场随机分析中响应量(水头与水力梯度)的随机响应公式,进一步得到流速以及流量的随机响应公式。由于渗流问题与平面弹性问题模型具有较大的相似性,文章在平面弹性问题的摄动随机有限元程序基础上,依据相似关系改造其为二维稳定渗流计算的摄动随机有限元程序,并扩充改进程序功能。最终程序可计算节点水头、单元水力梯度以及节点与单元流速、单元中线流量的均值与方差,文中介绍了程序主要算法的实现过程。There are a large number of uncertain factors existing in rock seepage model,to take the randomness of seepage parameters into finite element method analysis in seepage is a useful way aiming at making analysis results more scientific and more objective.The formulas of response variables（water head and hydraulic gradient） in stochastic were deduced based on the small parameter perturbation stochastic finite element method.The random response of velocity and flux was obtained.Based on the program of perturbation stochastic FEM on the plane elasticity model,we convert it to the one on the seepage model and perfect it at the same time,due to the high similarity between them.In summary,the final program can work out the effect of standard deviation and correlation structure of the permeability on water head and hydraulic gradient and the velocity and quantities of seepage.The program was also introduced in this context.

关 键 词：渗流场 摄动随机有限元法 局部平均离散 FORTRAN程序 

分 类 号：O357[理学—流体力学]