首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 234 毫秒
1.
Summary. We consider the problem of the unilateral extension of a two-dimensional anti-plane crack that initiates spontaneously at a point. The crack extends under the influence of cohesive resistance at the edge and dynamical friction along the crack walls. The stresses in the region beyond the edge of the crack are approximated so that they are exactly equal to the cohesive stresses near the edge of the crack, and are zero on the wavefront. An exact method of solving such problems is also given and can be used to determine the validity of the approximation. We find that the crack will not grow if the cohesion exceeds some critical value; this is consistent with an earlier result obtained by Knopoff, Mouton & Burridge for a similar one-dimensional model of crack propagation.  相似文献   

2.
Summary. Two-dimensional crack problems in elastic homogeneous isotropic media are considered which describe rupture over a fault surface characterized by non-uniform stress drop. Solutions can be found in which the stress field is finite at the crack tips and the rupture surface is not assigned a priori , but is part of the solution. These crack models are found to be consistent with the frictional stress threshold criterion for slip arrest over pre-existing fault surfaces. A crack is found to stop when its contribution to the stress field is opposite to the stress drop at the crack tips. The quasi-static propagation of a crack up to the arrest configuration is studied in terms of the minimum energy principle. The crack spontaneously propagates in such a way as to make the value of the stress intensity factor at one tip equal to the value at the other tip. Furthermore a tip propagating in a region with higher friction is found to move more slowly than the other tip propagating in a region with lower friction. Simple criteria for fracture arrest are derived, in terms of a properly averaged stress drop. Piecewise constant stress drop profiles are explicitly considered yielding a variety of solutions which can be applied to modelling asperities or barriers over a fault plane. The evaluation of the amount of the energy released during the quasi-static crack propagation shows that stopping phases cannot be efficiently radiated if the crack comes to rest in a low friction region.  相似文献   

3.
Summary. Following the classic work of Eshelby, the slip and stress distributions due to an elliptical plane shear crack are evaluated. The relation between average (or maximum) slip on the crack and the (constant) static stress drop, for faults of different aspect ratios, is found. The slip vector is not parallel to the applied stress but makes a small angle to it, except when the stress is applied along the major or minor axis of the ellipse. The stress -distribution around the crack shows that in addition to the expected stress concentration along the crack edge, there are broad regions of stress increase off the crack plane for circular and elliptical cracks, similar to those known to exist for in-plane but not for antiplane shear cracks. Whether the stress- intensity factor at the end of one axis is greater or less than that at the end of the other axis ( ka ≶ kb ), depends on the condition: √ b/a ≶ (1 − v ) where a and b are the semi-axes of the ellipse, ka and kb are the stress-intensity factors at the end of the a- and b -axes and v is Poisson's ratio. The total stress-intensity factor varies smoothly along the edge of the ellipse from one axis to the other and it is found that this variation is rather small.  相似文献   

4.
Summary. For a one-dimensional model for which the stress drop is a function of the relative displacement, the absence of an initial instantaneous stress drop from the static level inhibits the formation of a spontaneous, growing crack.  相似文献   

5.
Summary. This paper derives exact solutions to the equations of static plane strain by means of propagator matrices for homogeneous, gravitating and non- gravitating elastic media. These solutions are immediately verifiable and are flexible under a variety of boundary conditions. Propagator matrices are eminently suitable for computer encoding and, through their multiplication, are applicable to depth-dependent structures. Attention is focused upon the bending of floating plates which are loaded to simulate the deflection of oceanic lithosphere in the vicinity of trenches. By comparing responses computed with and without body forces, I find that gravity does not meaningfully change deflection profiles; however, it can influence important aspects of the internal stress state. Gravitational stresses are proportional to the gradient of the vertical deformation and amount to about 10 per cent of the bending stresses in these models. Propagators which include gravity are used to investigate the effect of regional horizontal stresses upon bending plates. I conclude that applied compressive forces can transport the neutral surface through 10–20 km of depth without significantly deforming the plate profile or increasing the maximum internal stress by more than 30 percent. These calculations support the contention that variable compressive stresses resulting from interplate coupling could account for observed regional differences in neutral surface height. For elastic-plastic material, the fundamental equations of motion become non-linear; however, there appears to be no a priori objection to their linearization. I speculate that the propagator formalism, when applied in an iterative approach, could be a powerful method for computing the deformation of such media.  相似文献   

6.
Summary. Most crustal earthquakes of the world are observed to occur within a seismogenic layer which extends from the Earth's surface to a depth of a few tens of kilometres at most. A model is proposed in which the shear zone along a transcurrent plate margin is represented as a viscoelastic medium with depth-dependent power-law rheology. A frictional resistance linearly increasing with depth is assumed on a vertical transcurrent fault within the shear zone. Such a model is able to reproduce a continuous transition from the brittle behaviour of the upper crust to the ductile behaviour at depth. Assuming that the shear zone is subjected to a constant strain rate from the opposite motions of the two adjacent plates, it is found that there exists a maximum depth H below which tectonic stress can never reach the frictional threshold: this may be identified as the maximum depth of earthquake nucleation. The value of H is consistent with observations for plausible values of the model parameters. The stress evolution in the shear zone is calculated in the linear approximation of the constitutive equation. A change in rigidity with depth, which is also introduced in the model, may reproduce the high vertical gradient of shear stress, which has been measured across the San Andreas fault, and the fact that most earthquakes are nucleated at some depth in the seismogenic layer. A crack which drops the ambient stress to the dynamic frictional level is then introduced in the model. To this aim, a crack solution is employed without a stress singularity at its edges, which is compatible with a frictional stress threshold criterion for fracture. A constraint on the vertical friction gradient is obtained if such cracks are assumed to be entirely confined within the seismogenic layer.  相似文献   

7.
Earthquake prediction: a new physical basis   总被引：16，自引：0，他引：16  
Summary. Subcritical crack growth in the laboratory occurs slowly but progressively in solids subjected to low stresses at low strain rates. The cracks tend to grow parallel to the maximum compressive stress so that, when this stress is aligned over a large enough region, the cracks will also be aligned and possess effective seismic anisotropy. It is suggested that such subcritical crack growth produces extensive-dilatancy anisotropy (EDA) throughout earth-quake preparation zones. This process is a possible driving mechanism for earthquake precursors observed at substantial distances from impending focal zones, and provides, in the shear-wave splitting which has been observed in several seismic regions, a possible technique for monitoring the build-up of stress before earthquakes.  相似文献   

8.
Summary. Current models of the structure of an active fault zone recognize two important subdivisions – an upper zone, extending to mid-crustal depths, in which processes associated with brittle fracture and friction dominate the fault behaviour, and a lower zone, extending into the mantle, within which stresses may be relieved by ductile flow. Anisotropy directly or indirectly induced by stress might occur throughout the fault zone, especially if caused by some form of stress-induced crack alignment. Dilatancy associated with high stresses is likely to be a very localized phenomenon in the vicinity of high strength regions (asperities), but alignments caused by subcritical crack growth at low stress and strain rate ( extensive-dilatancy anisotropy ) could give rise to anisotropy throughout the fault region.  相似文献   

9.
Summary. Dynamical rupture process on the fault is investigated in a quasi-three-dimensional faulting model with non-uniform distributions of static frictions or the fracture strength under a finite shearing pre-stress. The displacement and stress time functions on the fault are obtained by solving numerically the equations of motion with a finite stress—fracture criterion, using the finite difference method.
If static frictions are homogeneous or weakly non-uniform, the rupture propagates nearly elliptically with a velocity close to that of P waves along the direction of pre-stress and with a nearly S wave velocity in the direction perpendicular to it. The rise time of the source function and the final displacements are larger around the centre of the fault. In the case when the static frictions are heavily non-uniform and depend on the location, the rupture propagation becomes quite irregular with appreciably decreased velocities, indicating remarkable stick-slip phenomena. In some cases, there remain unruptured regions where fault slip does not take place, and high stresses remain concentrated up to the final stage. These regions could be the source of aftershocks at a next stage.
The stick—slip faulting and irregular rupture propagation radiate high-frequency seismic waves, and the near-field spectral amplitudes tend to show an inversely linear frequency dependence over high frequencies for heavily non-uniform frictional faults.  相似文献   

10.
Summary. A possible mechanism for the occurrence of slow earthquakes is investigated by calculating numerical solutions for the dynamical rupture process on a quasi-three-dimensional fault with heterogeneous frictional strengths. Experimental friction laws for the dependence of sliding frictional stress on slip velocity, which are based on the cohesive properties of fault asperities, are taken into considerations.
It is found that the applied stress does not drop very rapidly with time and the rupture velocity remarkably decreases as the dependence on slip-velocity becomes smaller. These deceleration effects for the rupture propagation are greatly enhanced with increasing heterogeneities in the distribution of frictional strength and as the initial shear stress has lower levels with respect to the average strength. For these cases, the growth of rupture is extremely slow in a nucleus region with the dimension as large as 10 times the initial rupture length, and gains a terminal velocity dependent on the above factors. The displacement-time function becomes noticeably extended in these cases, and indicates a stick—slip-like phenomena in the extended time interval for a strongly heterogeneous fault.
It seems that these results could explain the characteristic features of slow earthquakes.  相似文献   

11.
We perform analytical and numerical studies of scaling relations of earthquakes and partition of elastic strain energy between seismic and aseismic components using a thermodynamically based continuum damage model. Brittle instabilities occur in the model at critical damage level associated with loss of convexity of the strain energy function. A new procedure is developed for calculating stress drop and plastic strain in regions sustaining brittle instabilities. The formulation connects the damage rheology parameters with dynamic friction of simpler frameworks, and the plastic strain accumulation is governed by a procedure that is equivalent to Drucker–Prager plasticity. The numerical simulations use variable boundary forces proportional to the slip-deficit between the assumed far field plate motion and displacement of the boundary nodes. These boundary conditions account for the evolution of elastic properties and plastic strain in the model region. 3-D simulations of earthquakes in a model with a large strike-slip fault produce scaling relations between the scalar seismic potency, rupture area, and stress drop values that are in good agreement with observations and other theoretical studies. The area and potency of the simulated earthquakes generally follow a linear log–log relation with a slope of 2/3, and are associated with stress drop values between 1 and 10 MPa. A parameter-space study shows that the area-potency scaling is shifted to higher stress drops in simulations with parameters corresponding to lower dynamic friction, more efficient healing, and higher degree of seismic coupling.  相似文献   

12.
Summary. Anomalous stress distributions near arc—arc junctions have been found in several parts of the world. As a simple two-dimensional model, a thin elastic wedge is introduced to model the continental side of the arc—arc junction, and the stress distribution within it is calculated to interpret those anomalies. Uniform normal and/or shear stress boundary conditions are given separately to each edge to approximate the effect of the two neighbouring converging plate boundaries; Stress components are obtained by numerical integration. Among the various results, the following features are predominant and seem to be actually taking place at arc—arc junctions.
(1) The anomaly in stress pattern within the wedge is basically due to the difference in stress conditions between the two boundaries.
(2) A large stress concentration takes place when symmetrical shear stresses are applied to the two edges.  相似文献   

13.
With the purpose of improving the ice physics underpinning time–dependent glacier flowline models, three independent approaches for solving longitudinal stresses in glaciers are discussed and verified by application to Haut Glacier d'Arolla. To highlight any shortcomings, the reduced and much utilised driving stress approximation is also applied and compared. Modelled velocity patterns using the three full stress schemes exhibit consistency with one another and good coincidence with observed velocities for the 1991 summer melt season. Furthermore, these stress patterns indicate that longitudinal stresses are significant and of a similar order of magnitude as the basal shear stress components. However, the driving stress approximation yields erratic fluctuations in the stress and velocity fields which are neither realistic in terms of mass continuity nor agree with observations. Basal decoupling experiments indicate a complex relationship between basal velocity and englacial stresses with considerable dampening of any basal perturbation occurring as it is dissipated towards the surface and transferred throughout the ice mass. The driving stress approximation fails to account at all for any such coupling. Experiments to identify the length scale over which longitudinal effects operate indicate that they are significant even up to 10 ice thicknesses. The implication here is that longitudinal stresses play a significant role in determining glacier dynamics on length scales up to at least 2 km and that the predictive power of models of glacier flow based purely on the driving stress approximation is therefore subject to significant limitations. Inclusion of longitudinal stresses overcomes one of the main limitations imposed on such models and, given the potential ease of incorporation of the schemes described here, this deficiency may readily be resolved.  相似文献   

14.
Summary. The space and time characteristics of earthquake sequences, including a main shock, aftershocks and the recurrence of major shocks in a long time range, are investigated on a frictional fault model with non-uniform strengths and relaxation times, which is subjected to a time-dependent shear stress. Aftershocks with low stress drop take place successively in spaced regions so as to fill the gaps which have not yet been ruptured since the main shock, while those with high stress drop occur in and around the regions left unruptured during the main faulting. The frequency decay of aftershocks with time follows a hyperbolic law with the rates p consistent with observations. There are good linear relations in logarithmic scales for source area versus frequency and seismic moment versus frequency of the generated aftershocks. The b -value obtained in the present experiments appears slightly larger than that for observations. It was found that more heterogeneous distribution of the fault strength give smaller p and larger b -values. The recurrence of major shocks, particularly of very large shocks with high stress drop, is often preceded by a completely silent period of activity or very low activity with a small number of foreshocks. The major shocks take place successively in adjacent unruptured regions and sometimes show slow-speed migrations. These results provide explanations to various observations of earthquake sequences.  相似文献   

15.
A theoretical approach to the propagation of interacting cracks   总被引：2，自引：0，他引：2  
We propose a scheme to compute interaction effects between two randomly oriented cracks under compressive stresses and we discuss the role crack interactions play in the crack coalescence process. Stress intensity factors are computed by using an iterative technique based on the method of successive approximations. Once crack propagation occurs, curved wing cracks grow from the initial crack tips. The stress intensity factors at the wing crack tips are calculated as the sum of two terms: a component for a single wing crack subjected to both the applied stresses and the interaction effect, and a component due to the sliding of the initial crack. We have applied our procedure to various crack geometries. Our results show that interaction effects act on the crack propagation path. For cracks under tension, our approach correctly predicts the curving, hook-shaped paths of interacting cracks that have been observed in various materials. For en echelon compressive cracks, interaction effects depend on the geometry of stepping. For right-stepping cracks, no mode I crack coalescence occurs. A mixedmode propagation criterion may be introduced to check whether coalescing secondary shear fractures initiate. For left-stepping cracks, depending on whether or not there is overlapping, crack coalescence is achieved by tension wing cracks at the inner crack tips. Without overlapping, the growing wing cracks delimit a region where a tensile secondary fracture may develop and lead to coalescence. These results are consistent with previous work and show that our procedure may be now extended to a population of cracks.  相似文献   

16.
Wave speeds and attenuation of elastic waves in material containing cracks   总被引：38，自引：0，他引：38  
Summary. Expressions now exist from which may be calculated the propagation constants of elastic waves travelling through material containing a distribution of cracks. The cracks are randomly distributed in position and may be randomly orientated. The wavelengths involved are assumed to be large compared with the size of the cracks and with their separation distances so that the formulae, based on the mean taken over a statistical ensemble, may reasonably be used to predict the properties of a single sample. The results are valid only for small concentrations of cracks.
Explicit expressions, correct to lowest order in the ratio of the crack size to a wavelength, are derived here for the overall elastic parameters and the overall wave speeds and attenuation of elastic waves in cracked materials where the mean crack is circular, and the cracks are either aligned or randomly orientated. The cracks may be empty or filled with solid or fluid material. These results are achieved on the basis of simply the static solution for an ellipsoidal inclusion under stress.
The extension to different distributions of orientation or to mixtures of different types of crack is quite straightforward.  相似文献   

17.
Planetary topography can either be modelled as a load supported by the lithosphere, or as a dynamic effect due to lithospheric flexure caused by mantle convection. In both cases the response of the lithosphere to external forces can be calculated with the theory of thin elastic plates or shells. On one-plate planets the spherical geometry of the lithospheric shell plays an important role in the flexure mechanism. So far the equations governing the deformations and stresses of a spherical shell have only been derived under the assumption of a shell of constant thickness. However, local studies of gravity and topography data suggest large variations in the thickness of the lithosphere. In this paper, we obtain the scalar flexure equations governing the deformations of a thin spherical shell with variable thickness or variable Young's modulus. The resulting equations can be solved in succession, except for a system of two simultaneous equations, the solutions of which are the transverse deflection and an associated stress function. In order to include bottom loading generated by mantle convection, we extend the method of stress functions to include loads with a toroidal tangential component. We further show that toroidal tangential displacement always occurs if the shell thickness varies, even in the absence of toroidal loads. We finally prove that the degree-one harmonic components of the transverse deflection and of the toroidal tangential displacement are independent of the elastic properties of the shell and are associated with translational and rotational freedom. While being constrained by the static assumption, degree-one loads can deform the shell and generate stresses. The flexure equations for a shell of variable thickness are useful not only for the prediction of the gravity signal in local admittance studies, but also for the construction of stress maps in tectonic analysis.  相似文献   

18.
We present a theory for the radiation of high-frequency waves by earthquake faults. We model the fault as a planar region in which the stress drops to the kinematic friction during slip. This model is entirely equivalent to a shear crack. For two-dimensional fault models we show that the high frequencies originate from the stress and slip velocity concentrations in the vicinity of the fault's edges. These stress concentrations radiate when the crack expands with accelerated motion. The most efficient generation of high-frequency waves occurs when the rupture velocity changes abruptly. In this case, the displacement spectrum has an ω-2 behaviour at high frequencies. The excitation is proportional to the intensity of the stress concentration near the crack tips and to the change in the focusing factor due to rupture velocity. We extend these two-dimensional results to more general three-dimensional fault models in the case when the rupture velocity changes simultaneously on the rupture front. Results are similar to those described for two-dimensional faults. We apply the theory to the case of a circular fault that grows at constant velocity and stops suddenly. The present theory is in excellent agreement with a numerical solution of the same problem.
Our results provide upper bounds to the high-frequency radiation from more realistic models in which rupture velocity does not change suddenly. The ω-2 is the minimum possible decay at high frequencies for any crack model of the source.  相似文献   

19.
Seismic wave propagation through the earth is often strongly affected by the presence of fractures. When these fractures are filled with fluids (oil, gas, water, CO2, etc.), the type and state of the fluid (liquid or gas) can make a large difference in the response of the seismic waves. This paper summarizes recent work on methods of deconstructing the effects of fractures, and any fluids within these fractures, on seismic wave propagation as observed in reflection seismic data. One method explored here is Thomsen's weak anisotropy approximation for wave moveout (since fractures often induce elastic anisotropy due to non-uniform crack-orientation statistics). Another method makes use of some very convenient crack/fracture parameters introduced previously that permit a relatively simple deconstruction of the elastic and wave propagation behaviour in terms of a small number of crack-influence parameters (whenever this is appropriate, as is certainly the case for small crack densities). Then, the quantitative effects of fluids on these crack-influence parameters are shown to be directly related to Skempton's coefficient B of undrained poroelasticity (where B typically ranges from 0 to 1). In particular, the rigorous result obtained for the low crack density limit is that the crack-influence parameters are multiplied by a factor  (1 − B )  for undrained systems. It is also shown how fracture anisotropy affects Rayleigh wave speed, and how measured Rayleigh wave speeds can be used to infer shear wave speed of the fractured medium in some cases. Higher crack density results are also presented by incorporating recent simulation data on such cracked systems.  相似文献   

20.
In the cold regions of northern China, incidents of municipal underground gas pipeline rupture and leakage occur quite frequently, most often in winter. To prevent harm to citizen safety and property, analysis of the causes of such cracking and leakage is therefore valuable. Two incident analyses are discussed here and the reasons why most of these types of cases occur during winter are clarified. The effects of vehicle loadings above buried pipelines are calculated and compared with the effects and calculations of frost heaving forces. We demonstrate that when the soil layer above a pipeline freezes rapidly, the soil generates repeated frost heaving, which exerts heaving forces on the pipeline that can result in fatigue crack propagation and ultimate pipeline failure. Therefore, the incident induced by frost heaving is one of the primary reasons of gas pipeline failure. Based on these analyses, we present some recommendations pertaining to the proper design, construction, and management of gas pipelines.  相似文献   

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

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