A Novel Frictional Contact Interaction Algorithm for FDEM and Its Application in TBM Jamming
摘要
The three-dimensional combined finite-discrete element method (3D-FDEM) is widely used to simulate the fracture process and contact interaction between different blocks of rock materials. However, the conventional FDEM contact theory only includes the normal contact force and does not characterize the frictional force during relative motion. A novel FDEM contact algorithm is proposed by inheriting the advantages of the conventional algorithm. The magnitude and the direction of frictional force are determined based on the volume of the contact area and the direction of relative velocity between contacting objects. The novel algorithm can characterize the energy attenuation during collision and slide between different blocks, which is consistent with the physical scene; a main shear penetration crack is generated on the rock indentation simulation with the novel algorithm, consistent with the experimental results. The frictional penalty value PII is calibrated: when PII = μ × PI, the simulation results are consistent with the theoretical solutions. Finally, a tunnel excavation simulation is performed, and the frictional force on the shield is calculated when TBM is restarted after a while of standstill. The longer the analysis steps of standstill, the greater the thrust force required to overcome the frictional force on the shield. When the analysis steps of standstill exceed 845,000, TBM (shield) jamming will occur. The risk of TBM jamming will increase under the geological conditions of soft rock and high geostress. The risk of shield jamming can be mitigated by reducing the time of standstill. The novel algorithm contributes to optimizing the design of tunnel excavation patterns and advance rate, guaranteeing the stability of the surrounding rockmass.