基于速度-状态摩擦本构定律的三维DDA方法
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
原有三维非连续变形分析(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 data measured by Beeler et al.(1994)with a 0.999-goodness of fits.Even though both data sets are mixed together,there is also a 0.997-goodness of fits.In addition,the friction healing rateβcan be evaluated by the friction parameter b when the slip rate approximates to 0.The slope of these numerical data comes close to the b value of assumption with a 5% relative error.Besides,velocity stepping tests are also modeled using the improved 3D DDA.A slider keeps a steady sliding under a loading rate of 10μm/s with a distance of 25μm,and then decreases to 1μm/s during a characteristic distance.The results show a strong velocity dependence of friction which is consistent with laboratory data.Comparison between numerical results and laboratory data shows that the 3DDDA method in combination with rata-and state-dependent friction laws is capable of simulating velocity dependence of sliding friction and time dependence of static friction,which resolves a basic problem when using 3DDDA in geodynamics research.The improved 3D DDA method still has the limitation such as unbalanced embedment due to the uncertainty of stiffness on the discontinuous interface,which leads to inexact results of contact force determination.In the future,this modification can be used in quantitative simulation of regional crustal deformation in combination with observations,such as GPS data.
引文
Beeler N M,Tullis T E,Weeks J D.1994.The roles of time anddisplacement in the evolution effect in rock friction.Geophys.Res.Lett.,21(18):1987-1990.
    Chen B,Jiang Z S,Wang S X,et al.2000.Preliminary study onblock movement and its stress field by discontinuous deformationanalysis.Crustal Deformation and Earthquake(in Chinese),20(1):38-42.
    Chen Z A,Lin B H,Bai W M.2008.3-D numerical simulation oninfluence of 1997 Mani earthquake occurrence to stability oftectonic blocks system in Qingzang and Chuandian zone usingDDA+FEM method.Chinese J.Geophys.(in Chinese),51(5):1422-1430,doi:10.3321/j.issn:0001-5733.2008.05.015.
    Dieterich J H.1978.Time-dependence friction and the mechanics ofstick-slip.Pure Appl.Geophys.,116(4-5):790-806.
    Dieterich J H.1979.Modeling of rock friction:1.Experimentalresults and constitutive equations.J.Geophys.Res.,84(B5):2161-2168.
    Dieterich J H.1981.Constitutive properties of faults with simulatedgouge.//Carter N L,Friedman M,Logan J M eds.MechanicalBehavior of Crustal Rocks:The Handin Volume.Washington,DC:Am.Geophys.Union:24:103-120.
    Dieterich J H,Kilgore B D.1994.Direct observation of frictionalcontacts:new insights for state-dependent properites.PureAppl.Geophys.,143(1-3):283-302.
    He C R.1999.Comparing two types of rate and state dependentfriction laws.Seismology and Geology(in Chinese),21(2):137-146.
    Lapusta N,Rice J R.2003.Nucleation and early seismicpropagation of small and large events in a crustal earthquakemodel.J.Geophys.Res.,108(B4),doi:10.1029/2001JB00793.
    Law H K,Lam I P.2003.Evaluation of seismic performance fortunnel retrofit project.J.Geotech.Geoenviron.Eng.,129(7):575-589.
    Liu G,Shu D Q,Jiang Q H.2003.Application of discontinuousdeformation analysis method in analyzing shaft well stabilization.Blasting(in Chinese),20(4):38-44.
    Liu L B,Linde A T,Sacks I S,et al.1996.A seismic fault slip andblock deformation in north China.Pure Appl.Geophys.,146(3-4):717-740.
    Liu Y J,Rice J R.2007.Slow slip predictions based on granite andgabbro friction data compared to GPS measurements in northernCascadia.J.Geophys.Res.,114(B9):B09407,doi:10.1029/2008JB006142.
    Marone C.1998.Laboratory-derived friction laws and theirapplication to seismic faulting.Annu.Rev.Earth Planet.Sci.,26:643-696.
    Qin X J,Zhou S Y,Zhao L Q,et al.2001.Research on horizontalcrustal movement from 1994to 1998in Jiashi and the northeastto Pamir&its relation to strong earthquake swarm on the basisof GPS data.Crustal Deformation and Earthquake(inChinese),21(4):1-7.
    Ruina A L.1983.Slip instability and state variable friction law.J.Geophys.Res.,88(B12):10359-10370.
    Scholz C H.2002.The Mechanics of Earthquakes and Faulting.2nded.New York:Cambridge University Press.
    Shi G H.1992.Discontinuous deformation analysis:A newnumerical model for the statistics and dynamics of deformableblock structures.Engng.Comp.,9(2):157-168.
    Shi G H.2001.Three-dimensional discontinuous deformationanalysis.//Elsworth D ed.Proceedings of the 38th U.S.RockMechanics Symposium.Taylor&Francis,Washington,D.C.,USA,1421-1428.
    Sitar N,MacLaughlin M M,Doolin D M.2005.Influence ofkinematics on landslide mobility and failure mode.J.Geotech.Geoenviron.Eng.,131(6):716-728.
    Wang H,Zhang G M,Wu Y,et al.2003.The deformation ofactive crustal-blocks on the Chinese Mainland and its relationwith seismic activity.Earthquake Research in China(in Chinese),19(3):243-254.
    Wu A Q,Ding X L,Chen S H,et al.2006.Researches ondeformation and failure characteristics of an undergroundpowerhowse with complicated gelolgical conditions by DDAmethod.Chinese J.Rock Mech.Engng.(in Chinese),25(1):1-8.
    Wu J H.2010.Seismic landslide simulations in discontinuousdeformation analysis.Comp.Geotech.,37(5):594-601.
    Wu J H,Chen C H.2011.Application of DDA to simulatecharacteristics of the Tsaoling landslide.Comp.Geotech.,38(5):741-750.
    Wu Y Q,Chen G Q,Jiang Z S,et al.2012.The algorithm ofsimplex integration in three-dimension and its characteristicanalysis.Int.J.Adv.Comp.Tech.,4(10):246-256.
    Zhang L,Jiang Z S,Wu Y Q.2013.Review of rate-and statedependent friction laws and their applications to seismic faulting.Progress in Geophysics(in Chinese),28(5):2352-2362,doi:10.6038/pg20130517.
    陈兵,江在森,王双绪等.2000.利用非连续变形数值方法研究块体运动及其应力场初探.地壳形变与地震,20(1):38-42.
    陈祖安,林邦慧,白武明.2008.1997年玛尼地震对青藏川滇地区构造块体系统稳定性影响的三维DDA+FEM方法数值模拟.地球物理学报,51(5):1422-1430,doi:10.3321/j.issn:0001-5733.2008.05.015.
    何昌荣.1999.两种摩擦本构关系的对比研究.地震地质,21(2):137-146.
    刘刚,舒大强,姜清辉.2003.应用DDA方法分析竖井爆振稳定问题.爆破,20(4):38-44.
    秦小军,周硕愚,赵齐乐等.2001.依据GPS数据研究伽师及帕米尔东北侧1994—1998年地壳水平运动及其与强震群的关系.地壳形变与地震,21(4):1-7.
    王辉,张国民,吴云等.2003.中国大陆活动地块变形与地震活动的关系.中国地震,19(3):243-254.
    邬爱清,丁秀丽,陈胜宏等.2006.DDA方法在复杂地质条件下地下厂房围岩变形与破坏特征分析中的应用研究.岩石力学与工程学报,25(1):1-8.
    张龙,江在森,武艳强.2013.速度-状态摩擦本构定律及其在地震断层中的应用研究进展.地球物理学进展,28(5):2352-2362,doi:10.6038/pg20130517.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心