基于快速拉格朗日分析–并行粒子群算法的黏弹塑性参数反演及其应用  

INVERSION OF VISCOELASTO-PLASTIC PARAMETERS BASED ON FAST LAGRANGIAN ANALYSIS OF CONTINUUM-PARALLEL PARTICLE SWARM ALGORITHM AND ITS APPLICATION

作　　者：陈炳瑞[1] 冯夏庭[1] 黄书岭[1] 杨成祥[1] 张传庆[1] 

机构地区：[1]中国科学院武汉岩土力学研究所岩土力学与工程国家重点实验室,湖北武汉430071

出　　处：《岩石力学与工程学报》2007年第12期2517-2525,共9页Chinese Journal of Rock Mechanics and Engineering

基　　金：国家杰出青年科学基金资助项目(50325414);中国科学院知识创新工程青年人才领域专项前沿项目(0711031Q01);国家自然科学基金资助项目(50504004)

摘　　要：将改进粒子群算法和大型岩土数值计算软件FLAC3D有机结合,提出基于快速拉格朗日分析的全局并行改进粒子群参数反演方法。该方法既没有采用EBE方法也没有利用区域分解技术对计算区域进行处理,而是对用于参数反演的优化方法进行全局并行,这既避免了区域分解难于确定子区域边界的问题,又避免了EBE方法的大量数据交换,无疑将大大提高并行效率。基于锦屏二级水电站辅助洞BK14+599断面的算法性能分析证实了这一特性。锦屏二级水电站辅助洞黏弹塑性参数反演表明该方法对大规模岩土工程参数反演提供了一种行之有效的方法,获得了可靠的长期稳定性分析参数。利用获得的参数进行正向数值分析,计算结果表明:变形主要发生在顶拱和两侧墙靠上部位;洞周出现了一定范围的塑性区,以两侧为最,洞底次之,内壁出现了拉伸屈服;洞肩和墙角出现了应力集中现象,洞底出现底鼓现象,这些现象与现场观测到的结果一致。因此,引水隧洞在设计、施工时应加以注意,防止造成不必要的损失。An algorithm for back analysis of displacement is proposed. The method, regarding the monitoring displacement as a goal, estimates rheological mechanical parameters inversely in terms of residual sum of squares of the behavior predicted by the model and the actual behavior of rock mass under certain conditions using the improved parallel particle swarm algorithm hybridized with FLAC3D. Both EBE method and domain decomposition technology are not used to deal with the calculated zone, while optimization algorithm which designs calculating schemes is parallel in the methodology. The difficulties of sub-domains division and exchange of a great deal of data among personal computers are avoided, so efficiency and velocity must be improved well for back analysis of large-scale geotechnical engineering using numerical method. These are proven by analysis of capacity of the algorithm based on the example of section BK14＋599 at assistant tunnel B of Jinping second stage hydropower station. The viscoelasto-plastic parameters of assistant tunnel B of Jinping second stage hydropower station covered by about 1 538 m-thick rock mass, as an example, are obtained inversely using the proposed method and the obtained parameters are used for estimation of mechanical and deformation characteristics of assistant tunnel. The results indicate that. （1） the maximum deformation arises at vault and two side walls, mainly distributes along the two side walls; and （2） the phenomenon of floor heave comes up, and stress concentration comes into being at basement and spandrel of assistant tunnel B. These all agree with field investigation. So attentions must be paid to the occurrence of the above phenomena during the design and excavation of diversion tunnels.

关 键 词：数值分析 黏弹塑性 反分析 粒子群 并行算法 快速拉格朗日分析 锦屏二级水电站 

分 类 号：O246[理学—计算数学]