首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 921 毫秒
1.
基于速率与状态依赖性摩擦本构关系理论框架,在热水条件下研究了角闪石断层泥的摩擦滑动性质并与闪长岩的另一种主要矿物斜长石的摩擦滑动性质进行了对比.摩擦实验是在三轴实验系统上完成,有效正应力200 MPa,孔隙压力30 MPa,并将加载速率在1.22 μm/s和 0.122 μm/s之间实施了切换.结果表明角闪石的摩擦系数均值为0.70±0.01,随着温度增加没有系统性的变化,整体低于斜长石的摩擦系数(0.75±0.01);计算与实验表明,角闪石和斜长石的摩擦系数的体积分数加权平均值与闪长岩的摩擦系数基本一致;角闪石在实验温度范围内(100~614 ℃)显示速率强化(a-b>0),与斜长石在整个温度范围内的速率弱化(a-b<0)正好相反;角闪石的速率依赖性在整个实验温度范围内无系统性的变化.  相似文献   

2.
唐荣江  朱守彪 《地球物理学报》1954,63(10):3712-3726
地震是断层的摩擦失稳过程.摩擦本构关系对断层的破裂成核、破裂传播、破裂速度、能量释放、破裂终止等起着至关重要的控制作用.为了比较不同摩擦关系在断层自发破裂动力学过程中的影响,文中引入目前应用最为广泛的4种典型摩擦本构关系,它们分别是:滑移弱化摩擦关系,速率弱化摩擦关系,以及速率-状态相依摩擦关系中的老化定律和滑动定律.研究中利用有限单元方法对上述4种摩擦关系控制的断层自发破裂过程分别进行模拟计算,模拟结果显示:当模型参数相同时,不同摩擦关系模拟的破裂行为总体上具有一致性,都可以产生亚剪切破裂或超剪切破裂,并且破裂传播速度的大小与摩擦本构关系的类型无关.此外,它们之间还存在着较大的差异:(1)速率弱化摩擦关系可以模拟脉冲型破裂;而其他3个摩擦关系只能模拟裂纹性破裂.(2)不同摩擦关系模拟的超剪切破裂转换长度不同,速率-状态相关摩擦关系的老化定律相比其他摩擦关系需要更大的转换长度才能实现亚剪切破裂转变为超剪切破裂;而速率弱化的摩擦关系的超剪切转换长度可以为0,即不需要转换距离就直接产生超剪切破裂.(3)速率弱化摩擦关系模拟的破裂速度自成核后很快就达到稳定值;而其他类型摩擦关系模拟的破裂传播则要经历由缓慢破裂到逐渐加速直至达到稳定破裂的过程.值得特别指出的是,本文所使用的4种摩擦关系都不能完整地反映实际大地震破裂过程的摩擦属性,需要进一步深入研究.  相似文献   

3.
Seismicity simulation with a rate- and state-dependent friction law   总被引:3,自引:0,他引:3  
The dynamic motions and stabilities of a single-degree-of-freedom elastic system controlled by different friction laws are compared. The system consists of a sliding block connected to an elastic spring, driven at a constant velocity. The friction laws are a laboratory-inferred friction law called the rate-and-state-dependent friction law, proposed by Dieterich and Ruina, and a simple friction law described by dynamic and static frictions. We further extend the solution to a one-dimensional mass-spring model which is an analog of a fault controlled by the rate-and-state-dependent friction law. This model predicts non uniform slip and stress drop along the rupture length of a heterogeneous fault. This result is very different from some earlier modelings based on the simple friction law and a slip weakening friction law. In those earlier modelings the stress and slip functions become smoother with time along the length of the fault rupture, owing to the interactions between fault segments during slip. Because of this smoothing process the number of small events will decrease with time, and the universilly observed stationary magnitude-frequency relation cannot be explained. The interaction between a fault segment and its neighboring segments can be measured when the post-slip stress on this segment is compared with the stress on an identical segment (represented by a block in this modeling) without neighboring segments. If the post-slip stress of the former is much higher than that of the latter, strong interaction exists; if the two are close, only weak interaction exists. The interaction is determined by the relative motion between fault segments and the time duration of interaction. Our new modeling with the rate-and-state-dependent friction law appears to show no such smoothing effect and provides a physical mechanism for the roughening process in the difference between the fault strength and stress that is necessary to explain the observed stationary magnitude-frequency relation. The noninstantaneous healing predicted by the rate-and-state-dependent friction law may be repsonsible for the recurring nonuniform slip and stress drop, and may be explained by the reduction of interaction among fault segments due to the low frictional strength during the fault stopping. The very low friction during slip stopping allows much longer times than does the higher friction due to instantaneous healing for the fault segments to adjust their motions from an upper-limit slip velocity to almost rest. According to newton's second law, a process with fixed masses and constant velocity changes involves low forces and weak interactions if it is accomplished in a long time period, and vice versa. Our modeling also indicates that the existence of strong patches with higher effective stress on a fault is needed for the occurrence of major events. The creeping section of a fault, such as the one along the San Andreas fault in central California, on the other hand, can be simulated with the rate-and-state-dependent friction law by certain model parameters, which, however, must not include strong patches. In this case small earthquakes and aseismic creep relieve the accumulating strain without any large events.  相似文献   

4.
两种摩擦本构关系的对比研究   总被引:14,自引:1,他引:13       下载免费PDF全文
何昌荣 《地震地质》1999,21(2):42-146
目前有两种最常用的岩石摩擦本构关系,其主要区别在于是否与静接触时间有关。对这两种本构关系的基本性质进行了对比研究,得出如下结论:(1)在匀速滑动的稳态附近,两种本构关系趋向一致;(2)在正应力恒定条件下,两种本构关系的主要差异在于克服摩擦力所需能量的大小。与静接触时间相关的本构关系在粘滑中需消耗较大的能量来克服摩擦阻力,在粘滑的减速段,这种本构关系可达到的最低速率比另一种本构关系低10个数量级;(3)在正应力变化条件下,与静接触时间相关的本构关系基本继承了正应力恒定条件下的行为特征,而另一种本构关系在剪应力-速度相平面上却出现了一个长尾巴,使最低滑动速率降到比正应力恒定时的值低10多个数量级  相似文献   

5.
6.
Based on Zoback's method for estimating the tectonic stress magnitude and the two assumptions, we consider the conditions that three principal stresses are vertical principal stresses respectively(corresponding to three kinds of tectonic stress types). We deduced the formulae for estimating the tectonic stress magnitude by using the stress form factor and frictional strength of the fault and discussed the correlative influence of friction coefficient, pore pressure parameter and stress form factor on the stress value. When the maximum principal stress is approximately horizontal (when stress regime is strike-slip or reverse), the maximum principal stress (or the slope of stress increasing linearly with depth) is positively related with the friction coefficient and negatively related with the pore pressure coefficient. When the minimum principal stress is approximately horizontal (when stress regime is strike-slip or normal), the minimum principal stress (or the slope with depth) is negatively related to the friction coefficient, and positive to the pore pressure. Besides, these three parameters have great influence on the estimation of the tectonic stress magnitude. If the friction coefficient is too big and the pore pressure is too small, there could be a wide difference between the slope of the maximum principal stress increasing with depth and the slope of the minimum principal stress increasing with depth, which could lead to an unreasonable result. Our method is just an approximate estimation for the tectonic stress magnitude when crustal rocks have undergone brittle rupture or frictional sliding. The estimated results are not the tectonic stress magnitude when crust is in steady state.  相似文献   

7.
地震往往受控于滑动面的摩擦性质,这种摩擦性质可以由速率状态摩擦定律较好地描述.速率状态摩擦定律中的本构参数a和b与动态摩擦系数相关,从而影响着同震位移与剪切应力的时空演化.本文在前人工作的基础上,采用三维边界积分方程法模拟速率状态摩擦定律控制下均匀全空间中平面断层的自发破裂传播过程,并详细讨论了a和b对滑动速率、剪切应力和破裂传播速度的影响.数值结果表明a和b的不同取值将导致不同的破裂行为,b-a的值越大,断层越不稳定,这种不稳定性有利于裂纹的产生与扩展.但滑动速率的时空分布不只依赖b-a,而且还与a和b的具体取值有关,断层面上滑动速率峰值与剪切破裂强度均随着a的减小而增大,随着b的增大而增大.相关结果有助于加深对断层自发破裂传播的认识.  相似文献   

8.
The two one-state-variable, rate- and state-dependent friction laws, i.e., the slip and slowness laws, are compared on the basis of dynamical behavior of a one-degree-of-freedom spring-slider model through numerical simulations. Results show that two (normalized) model parameters, i.e., Δ (the normalized characteristic slip distance) and β?α (the difference in two normalized parameters of friction laws), control the solutions. From given values of Δ, β, and α, for the slowness laws, the solution exists and the unique non-zero fixed point is stable when Δ>(β?α), yet not when Δ < (β?α). For the slip law, the solution exists for large ranges of model parameters and the number and stability of the non-zero fixed points change from one case to another. Results suggest that the slip law is more appropriate for controlling earthquake dynamics than the slowness law.  相似文献   

9.
Based on dynamic rupture simulations on a planar fault in a homogeneous half-space, we investigated the nucleation processes using the time-weakening friction law. Both the characteristic time and the rupture speed in the nucleation asperity play an important role in determining rupture behaviors on a fault plane following the time-weakening friction law, with which rupture starts from a single point in the nucleation asperity and propagates at a given speed toward the boundary of the nucleation area. Rupture with a small characteristic time or a large rupture speed in the nucleation asperity propagates earlier from the hypocenter. Rupture following the slip-weakening friction law requires a smaller radius of nucleation patch to have similar rupture front contours of the time-weakening friction law. Even if the rupture velocity in the nucleation patch of the time-weakening friction law increases to infinity, the peak slip rate in the nucleation asperity is smaller than that of the slip-weakening law. The peak ground velocity distributions of ruptures following the two friction laws are also compared.  相似文献   

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

11.
In this study, we investigate what basic mechanisms a Discrete Element Model should have in order to reproduce the realistic wing crack extension, a widely observed phenomenon in uni-axial compression of brittle material with pre-existed crack. Using our Discrete Element Model—the Lattice Solid Model, we study how cracks propagate when different force-displacement laws are emplyed. Our results suggest that the basic features of crack propagation observed in laboratories cannot be reproduced under the following circumstances: 1) When only normal forces between two bonded particles exist and particle rotation is prohibited; 2) normal and shear stiffnesses are present and particle rotation is prohibited; 3) normal, shear stiffnesses and particle rotation are present and bending (rolling) stiffness is absent. Only when normal, shear and bending stiffness exist and particle rotation is permitted, is it possible to reproduce laboratory tests. We conclude that particle rotations and rolling resistance play a significant role and cannot be neglected while modeling such phenomenon. The effects of friction in the crack plane and confining pressure on extension of the cracks are also discussed.  相似文献   

12.
In this paper we consider a Linear Regression Model with a design matrix that fits the periodic structure of a time series. As a consequence, the residuals are very often autocorrelated. The main problem is that residual autocorrelation does not necessarily entail error autocorrelation. To analyse the effects of selecting different formulations to accommodate the autocorrelation in the residuals, we consider two seemingly different ways to deal with this problem: the Linear Regression Model with the error terms following an Autoregressive Stationary Process and the Partial Adjustment Model. We study the equivalence between the two formulations. We go over the problem of estimating the parameters and, especially, of making inferences in this framework. After parameter estimation, we analyse the adequacy of the models. We demonstrate that the issue of selecting the most appropriate model to capture the autocorrelation in the residuals is, in this context, a kind of an artefact since the main results concerning the fitted values and forecasting features are the same. These modelling procedures are applied to the Portuguese coastal upwelling data and we compare the estimated models.  相似文献   

13.
The freshwater budget of a tidal flat area is evaluated from long-term hydrographic time series from an observation pole positioned in a tidal channel in the Hörnum Basin (Germany). For each tidal cycle, the freshwater budget is calculated from the total imported and exported water volumes and the corresponding mean densities. The variability of the budget on a tidal scale is characterised by a period of twice the tidal period, exhibiting a minimum when the tidal flats are dry around daylight hours during the foregoing low tide, and a maximum when low tide occurs at night; enhanced evaporation on the flats at daylight hours is identified as the driving process. On the average over one year, while winter observations are missing, the freshwater budget is negative for the years 2002–2005 and positive only for 2006. The interannual mean is negative and amounts to a freshwater loss of about 2 mm day−1, although the large-scale climate in this region is humid. The results demonstrate that the bulk parametrisations for the latent and sensible heat flux between the ocean and the atmosphere must not be applied for the tidelands.  相似文献   

14.
Finite Element Analysis of a Sandwich Friction Experiment Model of Rocks   总被引:1,自引:0,他引:1  
-- Sandwich friction experiments are one of the most widely used standard methods for measuring the frictional behavior between rocks. A finite element code for modeling the nonlinear friction contact between elastoplastic bodies has been developed and extended to analyze the sandwich friction experiment model with a rate- and state-dependent friction law. The influences of prescribed slip velocity and variation of movement direction and state on the friction coefficient, the relative slip velocity, the normal contact force, the frictional force, the critical frictional force and the transition of stick-slip state between the deformable rocks are thoroughly investigated, respectively. The calculated results demonstrate the usefulness of this code for simulating the friction behavior between rocks.  相似文献   

15.
16.
震源动力学破裂过程数值模拟研究   总被引:2,自引:1,他引:1  
张丽芬  姚运生 《地震学报》2013,35(4):604-615
首先阐述了震源动力学过程研究的重要意义, 在此基础上, 研究了复杂的断层几何形态及介质模型对动力学破裂过程的影响, 并对常用的有限元方法、 离散元方法、 有限差分方法和边界积分方程方法等进行了相应介绍. 讨论了这些数值模拟方法各自的优缺点, 建议在方法的选择上应视具体问题及计算的精度而定. 最后对动力学数值模拟的关键部分, 滑动摩擦准则进行了论述. 常用的滑动摩擦准则有滑动弱化准则、 速率弱化准则和速率-状态依赖摩擦准则. 在单纯考虑某个地震的动力学破裂传播过程时, 滑动弱化准则较为常用, 其中滑动弱化距离的选取至关重要. 但若考虑整个地震循环, 速率-状态依赖摩擦准则更为合适.   相似文献   

17.
— Recent simulations using the particle dynamics method (PDM) have successfully captured many features of natural faults zones as illuminated in laboratory studies. However, 2-D simulations conducted on idealized assemblages of particles using simple elastic-frictional contact laws, yield friction values considerably lower than natural materials, and lack time- and velocity-dependent changes in strength that influence dynamic fault slip. Here, preliminary results of new PDM simulations are described, in which particle motions are restricted as a proxy for particle interlocking and out of plane contacts, and time-dependent contact healing is introduced to capture temporal strengthening of granular assemblages. Frictional strength is increased, and in the absence of interparticle rolling, can attain values observed in the laboratory. The resulting mechanical behavior is qualitatively similar to that described by empirically-based rate-state friction laws, providing new physical insight into the discrete mechanics of natural faults.  相似文献   

18.
热水条件下斜长石和辉石断层泥的摩擦滑动研究   总被引:4,自引:2,他引:2       下载免费PDF全文
罗丽  何昌荣 《地震地质》2009,31(1):84-96
在热水条件下对辉长岩的2种主要组成矿物斜长石和辉石进行了摩擦滑动实验研究。结果表明:1)斜长石在整个实验温度范围内(<607℃)显示速度弱化,而辉石除了在200℃附近表现出微弱的速度强化外,在其他温度范围也都表现为速度弱化;2)斜长石在低温时强度随温度的升高而略有增加,高温时强度则随温度的升高而降低。辉石的摩擦强度随着温度的增加没有系统性的变化,在0.74上下波动;3)在热水条件下辉长岩的摩擦性质并不是这2种主要矿物的叠加,有可能是辉长岩中的次要矿物的存在导致了辉长岩的速度强化;4)本实验所得热水条件下的实验结果与干燥条件下的结果有很大差异,再一次表明水的存在极大地影响了岩石的摩擦稳定性  相似文献   

19.
We develop two algorithms for estimating member stiffnesses and masses of a structure from measured modal response in conjunction with a finite element model of the structure. The mathematical model has known geometry and topology and parameterized constitutive properties. A few of the natural frequencies are measured and the corresponding modes are sampled at certain locations in space. The proposed algorithms are based on the concept of minimizing the sum of the squares of errors, specified as an index of discrepancy between the model and the structure, over all of the measured modes. The recursive quadratic programming method is used to solve the non-linear constrained estimation problem. Both proposed estimators can handle incompletely measured models, have robust convergence, and are amenable to modelling of complex structures. We demonstrate the use of our parameter estimation algorithm by applying it to identify the properties of a building from measured data.  相似文献   

20.
By using a single-degree-of-freedom spring-slider analog fault model, we generate a synthetic catalog of nearly 500 different seismic sequences. We explore the parameter space by assuming different values of constitutive parameters and tectonic environment. We also consider three different versions of the rate-dependent and state-dependent friction laws [the Dieterich-Ruina (DR), the Ruina-Dieterich (RD) and the Chester-Higgs (CH) models], and different approximations of the behavior of the friction at high sliding speeds, as well as the radiation damping effects. Our results indicate that for all the considered models, the recurrence time (T cycle) exhibits an inverse proportionality on the loading rate; a linear, positive dependence on the effective normal stress; and a linear, negative dependence on the characteristic distance controlling the state variable evolution. These results confirm and generalize previous studies. Remarkably, we found here that the coefficients of proportionality strongly depend on the adopted friction model, on the high speed behavior and on the reference set of parameters. Notably, we also found that the positive proportionality between T cycle and the difference ba, confirmed for DR and RD laws, does not hold in general for the CH law. Overall, we conclude that even in the simplest (and idealized) case of characteristic earthquakes considered here, in which the limiting cycle is reached by the system, and even in the framework of a very simplified fault model, the possibility to a priori predict, through an universal analytical relation, the inter-event time of an impending earthquake still remains only a dream. On the other hand, a numerical prediction of T cycle would require the exact knowledge of the rheological model (and its parameters at all times over the entire life of the fault) and the actual state of the fault, which indeed are often unknown.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号