颗粒离散元方法中接触模型的开发  

DEVELOPMENT OF CONTACT MODEL IN PARTICLE DISCRETE ELEMENT METHOD

作　　者：王涛[1,2] 吕庆[1] 李杨[1] 李宏明[1] 

机构地区：[1]武汉大学水资源与水电工程科学国家重点实验室,湖北武汉430072 [2]长江科学院水利部岩土力学与工程重点实验室,湖北武汉430010

出　　处：《岩石力学与工程学报》2009年第A02期4040-4045,共6页Chinese Journal of Rock Mechanics and Engineering

基　　金：国家自然科学基金资助项目(50879063);国家自然科学基金重大项目(90715042);长江科学院水利部岩土力学与工程重点实验室开放研究基金资助项目(G07-06)

摘　　要：在对颗粒离散元方法的理论基础进行深入研究的基础上,推导了广义Kelvin模型颗粒间接触力的理论解和数值解;采用面向对象的程序技术,通过颗粒离散元程序的接口,研制了广义Kelvin接触模型颗粒离散元计算程序;通过算例,验证了程序编制的正确性。在此基础上,以某深埋隧洞为例,计算不同接触本构模型条件下岩体的开裂破坏情况。结果表明,相对于弹性接触关系,采用广义Kelvin接触模型得出的结果与由Hoek-Brown强度准则得出的计算图更为接近。Generalized Kelvin model for calculation of contact force between particles is deduced based on the research of theoretical background of particle discrete element method.With the object-oriented programming technique and the redeveloped interface of PFC,the generalized Kelvin contact model is developed and used in a verification example.In the example,two particles are bonded together and the relaxation is calculated by PFC and theoretical method.The result shows that the contact force exponentially decreases as time increases,which coincides with the analytical solution.Based on the researches,taking a deep-buried tunnel for example,the cracking and failure of rock are analyzed under conditions of different contact constitutive models.Results show that the generalized Kelvin contact model is better than linear elastic contact model,and the former is more consistent with the calculation chart from Hoek-Brown strength criterion.

关 键 词：数值分析 接触模型 离散元方法 广义Kelvin模型 隧洞 

分 类 号：O242.1[理学—计算数学] TU45[理学—数学]