(1. State Key Laboratory of Geomechanics Geotechnical Engineering，Institute of Rock and Soil Mechanics，Chinese Academy of Sciences，Wuhan，Hubei 430071，China；2. University of Chinese Academy of Sciences，Beijing 100049，China)
Abstract：In order to avoid the outward waves being reflected back into the simulated domain at the artificial boundaries，the absorbing boundaries must be used in applying the granular discrete element method(PFC). To consider the various radii of particle elements on the artificial boundaries and the irregular boundary surfaces，the equivalent equation applicable for this discrete elements method was obtained based on the continuum viscous boundary condition. The calibration factors were introduced into the equivalent equation and a method to determine their values necessary for the optimum wave absorption was proposed. The joint network models of bedding and reverse bedding rock slopes were established with the synthetic rock mass approach (SRM). The viscous boundary conditions were built on the artificial boundaries in PFC models of rock slope based on the proposed method. At the same time，the seismic failure process of bedding and reverse bedding rock slopes with orthogonal secondary rock joints were simulated and systematic analyses of the seismic failure mechanism were conducted.
[1] LYSMER J，KUHLEMEYER R L. Finite dynamic model for infinite media[J]. Journal of the Engineering Mechanics Division，1969，95(4)：859–878.
[2] WHITE W，LEE I K，VALLIPAN S. Unified boundary for finite dynamic models[J]. Journal of the Engineering Mechanics Division，1976，103(5)：949–964.
[3] SMITH W D. A nonreflecting plane boundary for wave propagation problems[J]. Journal of Computational Physics，1974，15(4)：492–503.
[4] CERJAN C，KOSLOFF D，KOSLOFF R，et al. Absorbing boundary conditions for acoustic and elastic wave equations[J]. Bulletin of the Seismological Society of America，1977，67(6)：1 529–1 540.
[5] CLAYTON R W，ENGQUIST B. Absorbing boundary conditions for wave-equation migration[J]. Geophysics，1980，45(5)：895–904.
[6] DEEKS A J，RANDOLPH M F. Axisymmetric time-domain transmitting boundaries[J]. Journal of Engineering Mechanics，1994，120(1)：25–42.
[7] 刘晶波，吕彦东. 结构–地基动力相互作用问题分析的一种直接方法[J]. 土木工程学报，1998，31(3)：55–64.(LIU Jingbo，LU Yandong. A direct method for analysis of dynamic soil-structure interaction[J]. China Civil Engineering Journal，1998，31(3)：55–64.(in Chinese))
[8] 杜修力，赵 密，王进廷. 近场波动模拟的一种应力人工边界[J]. 力学学报，2006，38(1)：49–56.(DU Xiuli，ZHAO Mi，WANG Jingting. A stress artificial boundary in FEA for near-field wave problem[J]. Chinese Journal of Theoretical and Applied Mechanics，2006，38(1)：49–56.(in Chinese))
[9] 廖振鹏，杨柏坡，袁一凡. 暂态弹性波分析中人工边界的研究[J]. 地震工程与工程振动，1982，2(1)：1–11.(LIAO Zhenpeng，YANG Baipo，YUAN Yifan. Artificial boundary in analysis of transient elastic wave propagation[J]. Earthquake Engineering and Engineering Vibration，1982，2(1)：1–11.(in Chinese))
[10] 赵 密. 黏弹性人工边界及其与透射人工边界的比较研究[博士学位论文][D]. 北京：北京工业大学，2004.(ZHAO Mi. Study on the viscous-spring boundary and the transmitting boundary[Ph. D. Thesis][D]. Beijing：Beijing University of Technology，2004.(in Chinese))
[11] 陈万祥，郭志昆，袁正如，等. 地震分析中的人工边界及其在LS-DYNA中的实现[J]. 岩石力学与工程学报，2009，27(增2)： 3 504–3 515.(CHEN Wanxiang，GUO Zhikun，YUAN Zhengru，et al. Artificial boundary for seismic analysis and its applications to LS- DYNA[J]. Chinese Journal of Rock Mechanics and Engineering， 2009，28(Supp.2)：3 504–3 515.(in Chinese))
[12] 张 波，李术才，杨学英，等. 三维黏弹性人工边界地震波动输入方法研究[J]. 岩土力学，2009，30(3)：774–778.(ZHANG Bo，LI Shucai，YANG Xueying，et al. Research on seismic wave input with three dimensional viscoelastic artificial boundary[J]. Rock and Soil Mechanics，2009，30(3)：774–778.(in Chinese))
[13] Itasca Consulting Group. PFC2D(particle flow code in 2 dimensions) theory and background[R]. Minnesota，USA：Itasca Consulting Group，Inc.，2008.
[14] 刘新荣，傅 晏，郑颖人，等. 颗粒流细观强度参数与岩石断裂韧度之间的关系[J]. 岩石力学与工程学报，2011，30(10)：2 084– 2 089.(LIU Xinrong，FU Yan，ZHENG Yingren，et al. Relation between meso-parameters of particle flow code and fracture toughness of rock[J]. Chinese Journal of Rock Mechanics and Engineering，2011，30(10)：2 084–2 089.(in Chinese))
[15] 周 喻，吴顺川，张晓平. 岩石节理直剪试验颗粒流宏细观分析[J]. 岩石力学与工程学报，2012，31(6)：1 245–1 256.(ZHOU Yu，WU Shunchuan，ZHANG Xiaoping. Macro- and meso-analyses of rock joint direct shear test using particle flow theory[J]. Chinese Journal of Rock Mechanics and Engineering，2012，31(6)：1 245–1 256.(in Chinese))
[16] WANG C，TANNANT D D，LILLY P A. Numerical analysis of the stability of heavily jointed rock slopes using PFC2D[J]. International Journal of Rock Mechanics and Mining Sciences，2003，40(3)：415–424.
[17] SCHOLTÈS L U C，DONZÉ F V. Modeling progressive failure in fractured rock masses using a 3D discrete element method[J]. International Journal of Rock Mechanics and Mining Sciences，2012，52(6)：18–30.
[18] 李祥龙，唐辉明，胡 巍. 层面参数对顺层岩质边坡地震动力破坏过程影响研究[J]. 岩土工程学报，2014，36(3)：466–473.(LI Xianglong，TANG Huiming，HU Wei. Effect of bedding plane parameters on dynamic failure process of sliding rock slopes[J]. Chinese Journal of Geotechnical Engineering，2014，36(3)：466–473.(in Chinese))
[19] IVARS D M，PIERCE M E，DARCEL C，et al. The synthetic rock mass approach for jointed rock mass modelling[J]. International Journal of Rock Mechanics and Mining Sciences，2011，48(2)：219–244.
[20] 黄 达，岑夺丰，黄润秋. 单裂隙砂岩单轴压缩的中等应变率效应颗粒流模拟[J]. 岩土力学，2013，34(2)：535–545.(HUANG Da，CEN Duofeng，HUANG Runqiu. Influence of medium strain rate on sandstone with a single pre-crack under uniaxial compression using PFC simulation[J]. Rock and Soil Mechanics，2013，34(2)：535–545.(in Chinese))
[21] PIERCE M E，FAIRHURST C. Synthetic rock mass applications in mass mining[C]// Harmonising Rock Engineering and the Environment. Proceedings 12th ISRM International Congress. Beijing：[s. n.]，2011：9–14.
[22] THORNTON C. The conditions for failure of a face-centered cubic array of uniform rigid spheres[J]. Géotechnique，1979，29(4)：441–459.
2017, Vol. 36 
