Avisco-elastic dynamic constitutive model based on discrete element method and its application#br#
HE Xun1，HE Yong1，PAN Xuchao1，GUO Lei1，2
(1. School of Mechanical Engineering，Nanjing University of Science and Technology，Nanjing，Jiangsu 210094，China；2. Institute of Rocket Force，Beijing 100094，China)
Abstract：The visco-elastic dynamic constitutive relation is introduced into the Flat-Joint model based on the discrete element method，and a visco-elastic force-displacement formula and a rate-dependent broken criterion of contacts are derived. The rate-dependent force-displacement relationship and bond strength of contacts between particles are obtained. A subroutine(DLL) which could be called by PFC2D is generated with C++ based on the deduced formulas. By simulating the dynamic compression of rocks under different strain rates，relationships between the stress-strain curve，the quantity change curve of broken contacts and the distribution of broken strain of contacts with the strain rate are obtained，and the effects of the strain rate on the generation and distribution of cracks are analyzed. At the same time，the effects of visco-elastic parameters on the macro-parameters such as stress-strain curves and Poisson¢s rate，and deformation and failure characteristics of specimens are also analyzed. It is shown that the larger the elastic coefficient of elastic mediumis，the stronger the principal shear crack is，and the smaller the elastic coefficient of Maxwell mediumis，the more the shear cracks and dispersive cracks are. The viscous damping coefficient of Maxwell mediumhas little effect on the principle shear cracks，but the secondary cracks extending from the principal shear cracks increase with decreasing .
何 珣1，何 勇1，潘绪超1，郭 磊1，2. 基于Flat-Joint接触模型的黏弹性动态本构开发及其参数影响分析[J]. 岩石力学与工程学报, 2020, 39(5): 949-960.
HE Xun1，HE Yong1，PAN Xuchao1，GUO Lei1，2. Avisco-elastic dynamic constitutive model based on discrete element method and its application#br#. , 2020, 39(5): 949-960.
CUNDALL P A. A Discrete numerical model for granular assemble[J]. Geotechnique，1979，19(1)：47-65.
[11]
Itasca Consulting Group Inc. Particle flow code in 2 dimensions[R]. Minneapolis：Itasca Consulting Group Inc，2008.
[18]
HAO H，HAO Y，LI Z X. Advances in protectives structures research[M]. Boca Raton：CRC Press，2012：107-110.
[1]
DAI B，ZHAO G， KONIETZKY H，et al. Experimental and Numerical study on the damage evolution behavior of granitic rock during loading and unloading[J]. SCE Journal of Civil Engineering，2018，22(9)：3 278-3 291.
[4]
梁昌玉，李 晓，张 辉，等. 中低应变率范围内花岗岩单轴压缩特性的尺寸效应研究[J]. 岩石力学与工程学报，2013，32(3)：528-536.(LIANG Changyu，LI Xiao，ZHANG Hui，et al. Research on size effect of uniaxial compression properties of granite under medium and low strain rates[J]. Chinese Journal of Rock Mechanics and Engineering，2013，32(3)：528-536.(in Chinese))
[6]
XU Y， DAI F， XU N W，et al. Numerical Investigation of dynamic rock fracture toughness determination using a semi-circular bend specimen in split Hopkinson pressure bar testing[J]. Rock Mechanics and Rock Engineering，2016，49(3)：731-745.
[8]
李晓峰，李海波，刘 凯，等. 冲击荷载作用下岩石动态力学特性及破裂特征研究[J]. 岩石力学与工程学报，2017，36(10)：2 393-2 404. (LI Xiaofeng，LI Haibo，LIU Kai，et al. Dynamic properties and fracture characteristics of rocks subject to impact loading[J]. Chinese Journal of Rock Mechanics and Engineering，2017，36(10)：2 393- 2 404.(in Chinese))
[9]
黄 达，岑夺丰，黄润秋. 单裂隙砂岩单轴压缩的中等应变率效应颗粒流模拟[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))
[14]
石 崇，徐卫亚. 颗粒流数值模拟技巧与实践[M]. 北京：中国建筑工业出版社，2015：11-14. (SHI Chong，XU Weiya. Technique and practice of numerical simulation of particle flow[M]. Beijing： China Architecture and Building Press，2015：11-14.(in Chinese))
[16]
丛 宇，王在泉，郑颖人，等. 基于颗粒流原理的岩石类材细观参数的试验研究[J]. 岩土工程学报，2015，37(6)：1 031-1 040. (CONG Yu，WANG Zaiquan，ZHENG Yingren，et al. Experimental study on microscopic parameters of brittle materials based on particle flow theory[J]. Chinese Journal of Geotechnical Engineering，2015，37(6)：1 031-1 040.(in Chinese))
[2]
ZOU F，FANG Z，XIA M. Study on dynamic mechanical properties of limestone under uniaxial impact compressive loads[J]. Mathematical Problems in Engineering，2016，1：1-11.
[5]
梁书锋. 恒应变率冲击作用下花岗岩的损伤演化与本构模型研究[博士学位论文][D]. 北京：中国矿业大学(北京)，2016.(LIANG Shufeng. The study on damage evolution and constitutive model of the granite under constant strain rate impact[Ph. D thesis][D]. Beijing：China University of Mining and Technology (Beijing)，2016. (in Chinese))
[3]
卢志堂，王志亮. 中高应变率下花岗岩动力特性三轴试验研究[J]. 岩土工程学报，2016，38(6)：1 087-1 094.(LU Zhitang，WANG Zhiliang. Triaxial test study on dynamic properties of granite under medium-high strain rate loading[J]. Chinese Journal of Geotechnical Engineering，2016，38(6)：1 087-1 094.(in Chinese))
[7]
DU H，DAI F，XU Y，et al. Numerical investigation on the dynamic strength and failure behavior of rocks under hydrostatic confinement in SHPB testing[J]. International Journal of Rock Mechanics and Mining Sciences，2018，108：43-57.
[10]
张学朋，蒋宇静，王 刚，等. 基于颗粒离散元模型的不同加载速率下花岗岩数值试验研究[J]. 岩土力学，2016，39(9)：2 679-2 685. (ZHANG Xuepeng，JIANG Yujing，WANG Gang，et al. Numerical experiments on rate-dependent behaviors of granite base on particle discrete element model[J]. Rock and Soil Mechanics，2016，39(9)： 2 679-2 685.(in Chinese))
[13]
POTYONDY D O，CUNDALL P A. A Bonded-particle model for rock[J]. International Journal of Rock Mechanics and Mining Sciences，2004，41(8)：1 329-1 364.
[15]
金爱兵，刘 刚，杨振伟. 基于颗粒流程序的仿真节理岩体模型及其应用[J]. 岩土工程学报，2017，39(3)：540-546.(JIN Aibing，LIU Gang，YANG Zhenwei. Numerical model for Jointed Rock Mass in Particle Flow Code and its application[J]. Chinese Journal of Geotechnical Engineering，2017，39(3)：540-546.(in Chinese))
[17]
冯 春，李世海，刘晓宇. 基于颗粒离散元方法的连接键应变软化模型及其应用[J]. 力学学报，2016，48(1)：76-85.(FENG Chun，LI Shihai，LIU Xiaoyu. Particle-DEM based linked bar strain softening model and its application[J]. Chinese Journal of Theoretical and Applied Mechanics，2016，48(1)：76-85. (in Chinese))
2020, Vol. 39 
