中文摘要：

原有三维非连续变形分析（DDA）方法采用常摩擦系数的Mohr-Coulomb定律作为切向破坏准则,然而当描述更大尺度构造块体的运动与变形时,常摩擦系数不再适用.速度-状态摩擦本构定律能够定量描述地震周期各阶段断层面剪应力变化,解释发震断层行为.本文将速度-状态摩擦定律与三维DDA方法相结合,首先推导了计算摩擦系数的实用公式,随后通过滑动-保持-滑动实验与速度步进实验算例对改进的三维DDA方法进行了验证.结果表明,应用速度-状态摩擦本构定律的三维DDA方法能够比较准确地模拟静摩擦的时间依赖性与动摩擦的速度依赖性,解决了将三维DDA方法在地学中应用的基本问题.

英文摘要：

The elastic medium in lithosphere should be considered as discontinuities since the complex tectonic background has produced many active faults as separation.The Mohr-Coulomb joint failure criterion with a constant friction coefficient,adopted in the original 3Ddiscontinuous deformation analysis（DDA）method,cannot meet the requirement of highly accurate calculations for motion and deformation of tectonic block systems. The rate-and state-friction laws,which are capable of reproducing virtually the entire range of observed fault behaviors,are combined into the 3D DDA method.Firstly,the formula of computing coefficient of friction on the interface with rate-and state friction laws is derived.In order to calculate the value of state variableθ,slip velocity Vand friction coefficientμin every time step,a first-order differential equation about Vis deduced.The increment of Vis calculated by the second-order Taylor series expansion in our scheme.Secondly,the evolution law is determined by the Runge-Kutta scheme with adaptive step-size control.The friction submatrix,which consists of discrete forms of V andθ,is rewritten and then combined into the 3D DDA method.Finally,on the basis of reasonable geometry and mechanical properties of the numerical model,slide-hold-slide tests and velocity stepping tests are designed to examine the accuracy of the modified method. Suits of numerical slide-hold-slide tests are performed using hold time from 1to 10000 seconds and then we take the numerical test of the 10-second hold time as an example to make a brief illustration.A rigid block moves along the base in uniform linear motion at the first five seconds,because of the equality between friction and point loading.At the fifth second,the point loading is set to zero,and the friction strength drops immediately.After ten seconds,the block is reset as the previous loading.Strength increases,reaches a peak value and returns to its previous steady-state value subsequently.The numerical results are consistent with the laboratory da