首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
This paper further examines the possibility of modelling landslide as a consequence of the unstable slip in a steadily creeping slope when it is subject to perturbations, such as those induced by rainfall and earthquakes. In particular, the one-state variable friction law used in the landslide analysis by Chau is extended to a two-state variable friction law. According to this state variable friction law, the shear strength (τ) along the slip surface depends on the creeping velocity (V) as well as the two state variables (θ1 and θ2), which evolve with the ongoing slip. For translational slides, a system of three coupled non-linear first-order ordinary differential equations is formulated, and a linear stability analysis is applied to study the stability in the neighbourhood of the equilibrium solution of the system. By employing the stability classification of Reyn for three-dimensional space, it is found that equilibrium state (or critical point) of a slope may change from a ‘stable spiral’ to a ‘saddle spiral with unstable plane focus’ through a transitional state called ‘converging vortex spiral’ (i.e. bifurcation occurs), as the non-linear parameters of the slip surface evolve with its environmental changes (such as those induced by rainfall or human activities). If the one-state variable friction law is used in landslide modelling, velocity strengthening (i.e. dτss/dV > 0, where τss is the steady-state shear stress) in the laboratory always implies the stability of a creeping slope containing the same slip surface under gravitational pull. This conclusion, however, does not apply if a two-state variable friction law is employed to model the sliding along the slip surface. In particular, neither the region of stable creeping slopes in the non-linear parameter space can be inferred by that of velocity strengthening, nor the unstable region by that of velocity weakening. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

2.
The ultimate bearing capacity of a group of equally spaced multiple rough strip footings was determined due to the contribution of soil unit weight. The analysis was performed by using an upper bound theorem of limit analysis in combination with finite elements and linear programming. Along the interfaces of all the triangular elements, velocity discontinuities were considered. The value of ξγ was found to increase continuously with a decrease in S/B, where (i) ξγ is the ratio of the failure load of an interfering strip footing of a given width (B) to that of a single isolated strip footing having the same width and (ii) S is the clear spacing between any two adjacent footings. The effect of the variation of spacing on ξγ was found to be very extensive for small values of S/B; ξγ approaches infinity at S/B=0. In all the cases, the velocity discontinuities were found to exist generally in a zone only around the footing edge. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

3.
《Mathematical Geology》1995,27(6):731-748
Laser Doppler Anemometry (LDA) has proved a powerful tool for quantifying fluid turbulence and is increasingly being applied in fields such as fluvial sedimentology and geomorphology. When operated in the burst-signal processing mode, high-frequency velocity fluctuations are measured at irregular time internals. In many situations, there is a need to transform these data to obtain evenly spaced velocity values but at a lower frequency. However, clear guidelines for this type of data processing are lacking. Three steps are necessary in order to transform the original files into evenly spaced data: (1) resampling at the average sampling rate, (2) low-pass filtering with half-power frequency adjusted to the final sampling frequency, and (3) decimating at the desired frequency. The decision taken at each step will affect the resulting signal and may cause, if not assessed carefully, severe problems in the signal such as aliasing errors. This paper examines each stage of data processing and details the advantages and drawbacks of different techniques in relation to the effects on turbulence statistics (variance, instantaneous shear stress, etc.). A standard method and specific guidelines are finally proposed. An erratum to this article is available at .  相似文献   

4.
The aftershock sequence of the September 30th, 1993 Killari earthquake in the Latur district of Maharashtra state, India, recorded by 41 temporary seismograph stations are used for estimating 3-D velocity structure in the epicentral area. The local earthquake tomography (LET) method of Thurber (1983) is used. About 1500P and 1200S wave travel-times are inverted. TheP andS wave velocities as well asV P/VSratio vary more rapidly in the vertical as well as in the horizontal directions in the source region compared to the adjacent areas. The main shock hypocentre is located at the junction of a high velocity and a low velocity zone, representing a fault zone at 6–7 km depth. The estimated average errors ofP velocity andV P/VSratio are ±0.07 km/s and ±0.016, respectively. The best resolution ofP and S-wave velocities is obtained in the aftershock zone. The 3-D velocity structure and precise locations of the aftershocks suggest a ‘stationary concept’ of the Killari earthquake sequence.  相似文献   

5.
An analysis of the vane test using an Arbitrary Lagrangian–Eulerian formulation within a finite element framework is presented. This is suitable for soft clays for which the test is commonly used to measure in situ undrained shear strength. Constitutive laws are expressed in terms of shear stress–shear strain rates, and that permits the study of time effects in a natural manner. An analysis of the shear stress distributions on the failure surface according to the material model is presented. The effect of the constitutive law on the shear band amplitude and on the position of the failure surface is shown. In general, the failure surface is found at 1–1·01 times the vane radius, which is consistent with some experimental results. The problem depends on two dimensionless parameters that represent inertial and viscous forces. For usual vane tests, viscous forces are predominant, and the measured shear strength depends mainly on the angular velocity applied. That can explain some of the comparisons reported when using different vane sizes. Finally, the range of the shear strain rate applied to the soil is shown to be fundamental when comparing experimental results from vane, triaxial and viscosimeter tests. Appart from that, an experimental relation between undrained shear strength and vane angular velocity has been reproduced by this simulation. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

6.
A three-dimensional kinematic limit analysis approach based on the radial point interpolation method (RPIM) has been used to compute collapse loads for rectangular foundations. The analysis is based on the Mohr-Coulomb yield criterion and the associated flow rule. It is understood that the internal plastic power dissipation function and flow rule constraints can be expressed entirely in terms of plastic strain rates without involving stresses. The optimization problem has been solved on basis of the semidefinite programming (SDP) by using highly efficient primal-dual interior point solver MOSEK in MATLAB. The results have been presented in terms of the variation of the shape factors with changes in the aspect ratio (L/B) of the footing for different values of soil internal friction angle (ϕ). Computations have revealed that the shape factors, sc and sq, due to effects of cohesion and surcharge increase continuously with (1) decrease in L/B and (2) increase in ϕ. On the other hand, the shape factor sγ, due to the effect of soil unit weight, increases very marginally with an increase in L/B up to (1) ϕ = 25° for a rough footing and (2) ϕ = 35° for a smooth footing. Thereafter, for greater values of ϕ, the variation of sγ with L/B has been found to be quite similar to that of the factors sc and sq. The variations of (1) nodal velocity patterns, (2) plastic power dissipation, and (3) maximum plastic shear strain rates have also been examined to interpret the associated failure mechanism.  相似文献   

7.
The Ou Backbone Range strikes northwards through the central northeastern Japan arc and is bounded on both sides by the active reverse Uwandaira and Sen'ya faults. We have applied a traveltime inversion method (seismic tomography) with spatial velocity correlation to active and passive seismic data in order to investigate a three-dimensional (3-D) velocity structure. The data set contains 33,993 P- and 18,483 S-wave arrivals from 706 natural sources and 40 blasts, as well as 2803 P-wave traveltime data from 10 explosions detonated during the 1997 controlled source experiment. The traveltime inversion reveals a zone beneath the Ou Backbone Range in which P-wave velocities (VP) are approximately 6–8% lower than the average velocity at equivalent depths. The low VP and a low VP to S-wave velocity (VS) ratio (VP/VS) of about 1.65 suggest the presence of aqueous fluids in the middle crust.  相似文献   

8.
Physico-mechanical properties of rocks have great significance in all operational parts in mining activities, from exploration to final dispatch of material. Compressional wave velocity (p-wave velocity) and anisotropic behaviour of rocks are two such properties which help to understand the rock response under varying stress conditions. They also influence the breakage mechanism of rock. There are different methods to determine thep-wave velocity and anisotropyin situ and in the laboratory. These methods are cumbersome and time consuming. Fuzzy set theory, Fuzzy logic and Neural Networks techniques seem very well suited for typical geotechnical problems. In conjunction with statistics and conventional mathematical methods, hybrid methods can be developed that may prove to be a step forward in modeling geotechnical problems. Here, we have developed and compared two different models, Neuro-fuzzy systems (combination of fuzzy and artificial neural network systems) and Artificial neural network systems, for the prediction of compressional wave velocity.  相似文献   

9.
Using inversion of SV receiver functions, defined for various directions at each of the three broad-band stations located in the Baikal rift zone, detailed S velocity models of the crust and upper mantle down to 260 km have been obtained. These models reflect peculiarities of the velocity structure beneath Baikal depressions and mountains.  相似文献   

10.
Toward a model for airflow on the lee side of aeolian dunes   总被引：8，自引：0，他引：8  
The interaction between dunes and the primary wind results in a complex pattern of secondary airflow on the lee side of dunes. From 15 dunes studied during transverse flow conditions at Padre Island in Texas, White Sands in New Mexico, and the Algodones in California, distinct flow regions can generally be recognized, with the overall flow structure comparing favourably to that proposed for subaqueous bedforms. Downwind of dunes with flow separation is a back-flow eddy that extends about four dune-brink heights downwind from the brink of the dune. Beyond the separation cell, the velocity profiles can be divided into regions based upon segments separated by ‘kinks’ in the velocity profiles. The interior is an area above the dunes of relative high wind speed but low velocity gradient. Beneath the interior is the wake, which consists of two layers. The upper wake exhibits an uppermost portion where the flow decelerates while the remainder exhibits accelerating flow, so that the overall velocity gradient decreases downwind. The lower wake exhibits low velocity gradients and wind speeds that accelerate downwind at all heights, but primarily near the top of the layer, thereby causing the velocity gradient to increase downwind. At about eight dune heights downwind, the upper and lower wakes equilibrate to a single profile with the kink between them no longer apparent. The lowest recognizable region is the internal boundary layer. It is recognized by a relatively steep velocity gradient below the wake, and never exceeds a few tens of centimetres in height for our data set. Because of acceleration and increasing shear stress within this layer, interdune flats are at least potentially erosional. Overall, the wake and internal boundary layer show a downward transfer of momentum from upper regions so that the flow recovers. Where flow separation does not occur, simple flow expansion down the lee-face causes flow deceleration.  相似文献   

11.
The three dimensional scattering of near-field, from a point source, is studied for acceleration in the time domain. The perturbation method is applied to define the acceleration for the first order scattering from a weak inhomogeneity in a homogeneous surrounding. A body force, arising from the interaction between the primary waves and the inhomogeneity, acts as the source generating the scattered motion. The acceleration of scattered waves is related to the velocity and density fluctuations of the inhomogeneity. No restrictions are placed on the inhomogeneity size or locations of the source and receiver. Decoupling of scattered motion enables the identification of different phases. Integral expressions are derived for the scattering acceleration due to the incidence of near-field wave (from an impulsive point force) at a radially inhomogeneous volume element. These integrals are solved further for scattering from an inhomogeneous spherical shell. The accelerations for back scattering are obtained as a special case. These accelerations are simple analytically solvable expressions in closed form. Only spherical asymmetry ofP wave velocity inhomogeneity can affect the scatteredS acceleration. ScatteredP acceleration is affected by the gradient ofS wave velocity inhomogeneity. The back scattering of near-field from a spherical shell, is independent of radial inhomogeneity ofP wave velocity. Inhomogeneity with smoothly perturbedS wave velocity does not back-scatter any acceleration. Accelerations are computed numerically for scattering from a part of inhomogeneous spherical shell. Hypothetical models are considered to study the effects of the distances of spherical shell from source, receiver, its thickness and its position relative to the direction of impulsive force.  相似文献   

12.
The spatial distribution of soil shear-wave velocity and the fundamental period of vibration were selected as input parameters for the determination of potential seismic site effects in the Saguenay region, Canada. The methodology used in this study involved three clear steps. First, a 3D geological model of the surficial deposits was built taking into consideration the type, spatial distribution and thickness of the deposits. Second, representative average Vs values were determined for each of the major soil units. Finally, the average shear-wave velocity from the ground surface to bedrock (Vsav), the shear-wave velocity of the upper 30 m (Vs30) and the fundamental site resonance period (T0) were calculated over a regular grid for the study area. The results include the spatial distribution of the fundamental site resonance period, the average shear-wave velocity in the first 30 m of the ground and the spatial distribution of National Building Code of Canada seismic soil classes for the Saguenay region.  相似文献   

13.
This paper presents a method of calculating long-term settlement of a loaded pad on soft clay at Klang, Selangor, Malaysia where the soil model is treated as an anelastic material of viscoelastic property. Initially, an elastic shear modulus (G) value from shear wave velocity profiles of the seismic tests from spectral analysis of surface wave (SASW) and continuous surface wave (CSW) tests was obtained. A value of damping (D) at an equivalent elastic strain is then calculated from the hysteresis of the plate load tests corrected to equivalent strain using the Damping–Strain formula. The calculated elastic settlement and its equivalent damping are then used to calculate the long-term settlement by applying the generalised viscoelastic formula. Comparisons to traditional methods of settlement predictions were made and the viscoelastic formula has shown better agreement to the observed settlement. Further modification of the settlement formula is introduced to improve the settlement accuracy to 10%.  相似文献   

14.
One hundred and twenty-eight strong ground motion CALTECH (Earthquakes in the U.S.A., 1940–1971) records of five Californian earthquakes, recorded at ground level and in basements, for which the modified Mercalli macroseismic intensities are known, were statistically processed to calculate 165 wave parameters. Correlation of the peak values of particle acceleration, velocity and displacement, energy, impulse and root-mean-square amplitude of the vibrations, durations of vibrations for certain levels of the peak amplitude, and spectral energies of 10 frequency windows with macroseismic intensity are discussed from the viewpoint of classification of ground-motion records.Paper presented at the 21st General Assembly of the European Seismological Commission held in Sofia 1988.  相似文献   

15.
H.J. Melosh 《Tectonophysics》1976,35(4):363-390
This paper investigates the effect of shear heating in the asthenosphere on the thermal structure of the upper mantle. Equations describing the motion of the lithosphere over the asthenosphere in the presence of a strongly temperature-dependent stress-strain rate relation are derived and solved with the help of several approximations. These approximations are shown to be valid under conditions appropriate for the earth.Two sets of solutions are found. For one set (the “subcritical” solutions) a normal shear stress—velocity relation is found for small stresses. The velocity increases as the stress increases, reaching a maximum velocity σc for a critical stress σc. The subcritical solutions have a negligible effect on the thermal structure of the earth, even at the critical stress. The other set of solutions (the “supercritical” solutions) has the bizarre property that a decrease of applied shear stress leads to an increase of velocity. Thus, as the shear stress goes to zero, the velocity becomes infinite. At larger shear stresses the velocity decreases until it reaches σc at a stress σc (the two sets of solutions share this point in common). There are no steady solutions of any kind for shear stresses in excess of σc. We discard the supercritical solutions as candidates for the thermal structure of the earth on the basis of their instability to small perturbations of applied stress and temperature.The realm of subcritical solutions (stress less than σc, velocity less than σc) thus defines a regime of plate motion in which the thermal effects of shear heating are negligible. If the shear stresses acting on plates exceed σc, however, new physical processes must come into play to dissipate the excess heat generated. Assuming that the velocities of plates on the earth today are less than σc, relative to the deep mantle, a strict upper limit of a few tens of bars can be derived for σc, corresponding to effective viscosities of ca. 1019 poise in the asthenosphere.  相似文献   

16.
Variations in elastic wave velocity, low-frequency internal friction and acous-tic emission in granite were experimentally studied as a function of temperature.It is found that the wave velocity and Young‘s modulus tend to decrease with temperature.In association with the α-β transition of quartz a sharp internal friction peak can be recognized accompanied by a drastic drop in wave velocity and modulus and by a second peak of acoustic emission rate.  相似文献   

17.
The links between large‐scale turbulence and the suspension of sediment over alluvial bedforms have generated considerable interest in the last few decades, with past studies illustrating the origin of such turbulence and its influence on flow resistance, sediment transport and bedform morphology. In this study of turbulence and sediment suspension over large sand dunes in the Río Paraná, Argentina, time series of three‐dimensional velocity, and at‐a‐point suspended sediment concentration and particle‐size, were measured with an acoustic Doppler current profiler and laser in situ scattering transmissometer, respectively. These time series were decomposed using wavelet analysis to investigate the scales of covariation of flow velocity and suspended sediment. The analysis reveals an inverse relationship between streamwise and vertical velocities over the dune crest, where streamwise flow deceleration is linked to the vertical flux of fluid towards the water surface in the form of large turbulent fluid ejections. Regions of high suspended sediment concentration are found to correlate well with such events. The frequencies of these turbulent events have been assessed from wavelet analysis and found to concentrate in two zones that closely match predictions from empirical equations. Such a finding suggests that a combination and interaction of vortex shedding and wake flapping/changing length of the lee‐side separation zone are the principal contributors to the turbulent flow field associated with such large alluvial sand dunes. Wavelet analysis provides insight upon the temporal and spatial evolution of these coherent flow structures, including information on the topology of dune‐related turbulent flow structures. At the flow stage investigated, the turbulent flow events, and their associated high suspended sediment concentrations, are seen to grow with height above the bed until a threshold height (ca 0·45 flow depth) is reached, above which they begin to decay and dissipate.  相似文献   

18.
Abstract: A three-dimensional local-scale P-velocity model down to 25 km depth around the main shock epicenter region was constructed using 83821 event-to-receiver seismic rays from 5856 aftershocks recorded by a newly deployed temporary seismic network. Checkerboard tests show that our tomographic model has lateral and vertical resolution of ~2 km. The high-resolution P-velocity model revealed interesting structures in the seismogenic layer: (1) The Guanxian-Anxian fault, Yingxiu-Beichuan fault and Wenchuan-Maoxian fault of the Longmen Shan fault zone are well delineated by sharp upper crustal velocity changes; (2) The Pengguan massif has generally higher velocity than its surrounding areas, and may extend down to at least ~10 km from the surface; (3) A sharp lateral velocity variation beneath the Wenchuan-Maoxian fault may indicate that the Pengguan massif’s western boundary and/or the Wenchuan-Maoxian fault is vertical, and the hypocenter of the Wenchuan earthquake possibly located at the conjunction point of the NW dipping Yingxiu-Beichuan and Guanxian-Anxian faults, and vertical Wenchuan-Maoxian fault; (4) Vicinity along the Yingxiu-Beichuan fault is characterized by very low velocity and low seismicity at shallow depths, possibly due to high content of porosity and fractures; (5) Two blocks of low-velocity anomaly are respectively imaged in the hanging wall and foot wall of the Guanxian-Anxian fault with a ~7 km offset with ~5 km vertical component.  相似文献   

19.
The heterogeneous upper mantle low velocity zone   总被引：2，自引：1，他引：2  
The upper mantle low velocity zone (LVZ) is a depth interval with slightly reduced seismic velocity compared to the surrounding depth intervals. The zone is present below a relatively constant depth of 100 km in most continental parts of the world, both in cratonic areas with high average velocity and tectonically active areas with low average velocity. Evidence for the low velocity zone arises from controlled and natural source seismology, including studies of surface waves and of primary and multiple reflections of body waves from the bounding interfaces, calculations of receiver functions, and absolute velocity tomography. The available data indicates a more pronounced reduction in seismic velocity and Q-value for S-waves than P-waves as well as high electrical conductivity in the LVZ. Seismic waves are strongly scattered by the zone, which demonstrates the existence of small-scale heterogeneity. The depth to the base of the LVZ is systematically shallower in cold, stable cratonic areas than in hot, active regions of the world. Because of its global occurrence below a relative constant depth of 100 km, the LVZ cannot be explained by metamorphic or compositional variation and rheological changes. Calculated upper mantle temperatures indicate that the rocks are close to the solidus in an interval with variable thickness below 100 km depth, provided that the rocks contain water and carbon dioxide. The presence of, even small amounts of such fluids in the mantle rocks will lower the solidus by several hundred degrees and introduce a characteristic kink on the solidus curve around 80–100 km depth. The seismic velocities and Q-values are significantly reduced of rocks, which are close to the solidus or contain small amounts of partial melt. Hence, the LVZ may be explained by upper mantle temperatures being close to the solidus in a depth interval below 100 km. Assuming that the rocks contain only limited amounts of fluids, this mechanism may explain the low velocities, Q-values, and resistivity, as well as the intrinsic scattering, and the characteristic variation in thickness of the low velocity zone.  相似文献   

20.
A finite element algorithm for frictionless contact problems in a two‐phase saturated porous medium, considering finite deformation and inertia effects, has been formulated and implemented in a finite element programme. The mechanical behaviour of the saturated porous medium is predicted using mixture theory, which models the dynamic advection of fluids through a fully saturated porous solid matrix. The resulting mixed formulation predicts all field variables including the solid displacement, pore fluid pressure and Darcy velocity of the pore fluid. The contact constraints arising from the requirement for continuity of the contact traction, as well as the fluid flow across the contact interface, are enforced using a penalty approach that is regularised with an augmented Lagrangian method. The contact formulation is based on a mortar segment‐to‐segment scheme that allows the interpolation functions of the contact elements to be of order N. The main thrust of this paper is therefore how to deal with contact interfaces in problems that involve both dynamics and consolidation and possibly large deformations of porous media. The numerical algorithm is first verified using several illustrative examples. This algorithm is then employed to solve a pipe‐seabed interaction problem, involving large deformations and dynamic effects, and the results of the analysis are also compared with those obtained using a node‐to‐segment contact algorithm. The results of this study indicate that the proposed method is able to solve the highly nonlinear problem of dynamic soil–structure interaction when coupled with pore water pressures and Darcy velocity. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

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

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