Huang Xiaolin, Zhao Qi, Qi Shengwen, Xia Kaiwen, Grasselli Giovanni, Chen Xuguang
Key Laboratory of Shale Gas and Geoengineering, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China.
Department of Civil Engineering, University of Toronto, Toronto, ON M5S 1A4, Canada.
Materials (Basel). 2016 Dec 27;10(1):13. doi: 10.3390/ma10010013.
This paper numerically investigates the seismic response of the filled joint under high amplitude stress waves using the combined finite-discrete element method (FDEM). A thin layer of independent polygonal particles are used to simulate the joint fillings. Each particle is meshed using the Delaunay triangulation scheme and can be crushed when the load exceeds its strength. The propagation of the 1D longitude wave through a single filled joint is studied, considering the influences of the joint thickness and the characteristics of the incident wave, such as the amplitude and frequency. The results show that the filled particles under high amplitude stress waves mainly experience three deformation stages: (i) initial compaction stage; (ii) crushing stage; and (iii) crushing and compaction stage. In the initial compaction stage and crushing and compaction stage, compaction dominates the mechanical behavior of the joint, and the particle area distribution curve varies little. In these stages, the transmission coefficient increases with the increase of the amplitude, i.e., peak particle velocity (PPV), of the incident wave. On the other hand, in the crushing stage, particle crushing plays the dominant role. The particle size distribution curve changes abruptly with the PPV due to the fragments created by the crushing process. This process consumes part of wave energy and reduces the stiffness of the filled joint. The transmission coefficient decreases with increasing PPV in this stage because of the increased amount of energy consumed by crushing. Moreover, with the increase of the frequency of the incident wave, the transmission coefficient decreases and fewer particles can be crushed. Under the same incident wave, the transmission coefficient decreases when the filled thickness increases and the filled particles become more difficult to be crushed.
本文采用组合有限离散元法(FDEM)对高幅值应力波作用下填充节理的地震响应进行了数值研究。用一层独立的多边形颗粒来模拟节理填充物。每个颗粒采用德劳内三角剖分方案进行网格划分,当荷载超过其强度时会发生破碎。研究了一维纵波通过单个填充节理的传播,考虑了节理厚度和入射波特性(如幅值和频率)的影响。结果表明,高幅值应力波作用下的填充颗粒主要经历三个变形阶段:(i)初始压实阶段;(ii)破碎阶段;(iii)破碎压实阶段。在初始压实阶段和破碎压实阶段,压实主导节理的力学行为,颗粒面积分布曲线变化不大。在这些阶段,透射系数随入射波幅值(即峰值颗粒速度,PPV)的增加而增大。另一方面,在破碎阶段,颗粒破碎起主导作用。由于破碎过程产生的碎片,颗粒尺寸分布曲线随PPV突然变化。这个过程消耗了部分波能,降低了填充节理的刚度。在这个阶段,由于破碎消耗的能量增加,透射系数随PPV的增加而减小。此外,随着入射波频率的增加,透射系数减小,可破碎的颗粒减少。在相同入射波作用下,填充厚度增加时透射系数减小,填充颗粒更难破碎。