共查询到20条相似文献，搜索用时 46 毫秒
1.
Scattering of SH waves by an embedded rigid elliptic cylinder of finite length, which is partially debonded from elastic soil is studied. The debonding regions are modeled as multiple elliptic arc-shaped interface cracks with non-contacting faces. The scattered wave field is expressed as a Mathieu function expansion with unknown coefficients. The mixed boundary conditions of the problem lead to a set of singular integral equations of the first type in terms of the dislocation density functions of the cracks. A quadrature method is used to solve these integral equations numerically. The results for dynamic stress intensity factors, far-field pattern of the displacement and scattering cross sections are presented. We particularly discuss the effects of the ratio of the short radius to long radius of the cylinder. 相似文献
2.
G. S. Murty 《Pure and Applied Geophysics》1971,90(1):96-108
Summary An infinite cylinder with a small ridge on the surface parallel to the axis is assumed to be excited by an embedded periodic line source radiating compressional waves. The line source is also assumed to lie parallel to the axis, and the frequency of radiation is small enough so that the corresponding wavelength is large compared to the dimension of the ridge on the surface of the cylinder. Considering the effect of the ridge as a small perturbation, the wave equations are solved and the amplitudes of the potential fields are expressed in terms of a set of inhomogeneous algebraic equations in which the solution of the problem without the ridge supplies the source term for the first order perturbation in potential fields. 相似文献
3.
The direct finite element method is a type commonly used for nonlinear seismic soil-structure interaction(SSI) analysis. This method introduces a truncated boundary referred to as an artificial boundary meant to divide the soilstructure system into finite and infinite domains. An artificial boundary condition is used on a truncated boundary to achieve seismic input and simulate the wave radiation effect of infinite domain. When the soil layer is particularly thick, especially for a three-dimensional problem, the computational efficiency of seismic SSI analysis is very low due to the large size of the finite element model, which contains an whole thick soil layer. In this paper, an accurate and efficient scheme is developed to solve the nonlinear seismic SSI problem regarding thick soil layers. The process consists of nonlinear site response and SSI analysis. The nonlinear site response analysis is still performed for the whole thick soil layer. The artificial boundary at the bottom of the SSI analysis model is subsequently relocated upward from the bottom of the soil layer(bedrock surface) to the location nearest to the structure as possible. Finally, three types of typical sites and underground structures are adopted with seismic SSI analysis to evaluate the accuracy and efficiency of the proposed efficient analysis scheme. 相似文献
4.
Dynamic faulting under rate-dependent friction 总被引：2，自引：0，他引：2
We discuss the effects of rate-dependent friction on the propagation of seismic rupture on active faults. Several physicists using Burridge and Knopoff's box and spring model of faulting have proposed that fault complexity may arise from the spontaneous development of a self-similar stress distribution on the fault plane. If this model proves to be correct, it has important consequences for the origin of the complexity of seismic sources. In order to test these ideas on a more realistic earthquake model, we developed a new boundary integral equation method for studying rupture propagation along an antiplane fault in the presence of nonlinear rate-dependent friction. We study rupture dynamics of models with single and twin asperities. In our models, asperities are places on the fault with a higher value of prestress. Othewise all fault parameters are homogeneous. We show that for models with such asperities, a slip velocity weakening friction leads to the propagation of supersonic healing phases and to the spontaneous arrest of fracture if the prestress outside the asperities is low enough. For models with asperities, we can also observe narrow slip velocity pulses, qualitatively similar to the so-called Heaton pulses observed in some earthquake accelerograms. We also observe a complex distribution of stress after the rupture that depends on details of the initial distribution of asperities and on the details of the friction law. 相似文献
5.
通过重新编制活动断层定量参数数据库,利用中国大陆活动断裂定量资料评估了中国大陆活动断裂(段)的大震复发间隔,给出了青藏、新疆、华北3个构造区的大震复发间隔分布特征。研究结果表明,青藏高原活动断裂M7.5级以上大震复发间隔远小于华北和新疆。青藏高原从几100a到2000a,而华北和新疆则从2000a到5000a。大震复发间隔与断层滑动速率成反比,边界断裂的大震复发间隔通常比块体内部断裂短,比例关系大体为1:1.5—2.0。大震复发间隔与断裂性质相关,走滑断层的复发间隔小于逆走滑断层,逆走滑断层小于逆断层,比例关系大体为1:1.5—2.5。现有的定量数据汇编成果和大震复发间隔估计以及统计结果,既可作为地震危险性分析的重要参数,又可作为构造类比研究的参考信息。 相似文献
6.
A simple calculation procedure for estimating absolute maximum slip displacement of a freestanding rigid body placed on the ground or floor of linear/nonlinear multi‐storey building during an earthquake is developed. The proposed procedure uses the displacement induced by the horizontal sinusoidal acceleration to approximate the absolute maximum slip displacement, i.e. the basic slip displacement. The amplitude of this horizontal sinusoidal acceleration is identical to either the peak horizontal ground acceleration or peak horizontal floor response acceleration. Its period meets the predominant period of the horizontal acceleration employed. The effects of vertical acceleration are considered to reduce the friction force monotonously. The root mean square value of the vertical acceleration at the peak horizontal acceleration is used. A mathematical solution of the basic slip displacement is presented. Employing over one hundred accelerograms, the absolute maximum slip displacements are computed and compared with the corresponding basic slip displacements. Their discrepancies are modelled by the logarithmic normal distribution regardless of the analytical conditions. The modification factor to the basic slip displacement is quantified based on the probability of the non‐exceedence of a certain threshold. Therefore, the product of the modification factor and the basic slip displacement gives the design slip displacement of the body as the maximum expected value. Since the place of the body and linear/nonlinear state of building make the modification factor slightly vary, ensuring it to suit the problem is essential to secure prediction accuracy. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
7.
Based on knowledge of a commutative group calculation of the rock stiffness and on some geophysical assumptions, the simplest fractured medium may be regarded as a fracture embedded in an isotropic background medium, and the fracture interface can be simulated as a linear slip interface that satisfies non‐welded contact boundary conditions: the kinematic displacements are discontinuous across the interface, whereas the dynamic stresses are continuous across the interface. The finite‐difference method with boundary conditions explicitly imposed is advantageous for modelling wave propagation in fractured discontinuous media that are described by the elastic equation of motion and non‐welded contact boundary conditions. In this paper, finite‐difference schemes for horizontally, vertically, and orthogonally fractured media are derived when the fracture interfaces are aligned with the boundaries of the finite‐difference grid. The new finite‐difference schemes explicitly have an additional part that is different from the conventional second‐order finite‐difference scheme and that directly describes the contributions of the fracture to the wave equation of motion in the fractured medium. The numerical seismograms presented, to first order, show that the new finite‐difference scheme is accurate and stable and agrees well with the results of previously published finite‐difference schemes (the Coates and Schoenberg method). The results of the new finite‐difference schemes show how the amplitude of the reflection produced by the fracture varies with the fracture compliances. Later, comparisons with the reflection coefficients indicate that the reflection coefficients of the fracture are frequency dependent, whereas the reflection coefficients of the impedance contrast interface are frequency independent. In addition, the numerical seismograms show that the reflections of the fractured medium are equal to the reflections of the background medium plus the reflections of the fracture in the elastic fractured medium. 相似文献
8.
The induced polarization response for a 2D horizontal cylinder embedded in a half-space is calculated for a uniform electric source. Response curves, in the form of apparent charge-ability taking into account the effect of the air-earth interface, exhibit a sharp decrease in amplitude with an increase in depth of burial of the target. The resistivity contrast between the cylinder and the host plays a dominant role in determining the IP response, i.e. the amplitude decreases considerably with the increase in resistivity contrast. The decrease is due to the defocusing effect caused by the resistive cylinder. The current lines tend to deviate away from the cylindrical target. In the case of a highly conducting cylinder, apparent defocusing takes place as current lines are confined to the surface of the conducting cylinder. An increase in chargeability contrast is reflected as a steady rise in the response. The peak response at the centre is reduced by about half the magnitude when the air–earth interface is not considered. The variation of response along the profile, though noticeable, is not as high as that obtained at the centre. 相似文献
9.
Abstract Stability analysis is formulated for a two-layer fluid model in which the upper and lower layers are convectively stable and unstable, respectively. With discontinuities in viscosity and conductivity at the interface, the exchange of stability does not generally hold and overstability is possible. A detailed analytical treatment is presented for the case of small viscosity and conductivity in which viscous and conducting boundary layers are formed at the interface. The usual damping effect due to the energy dissipation by viscosity and thermal conductivity exists irrespective of whether the mode is the convection or the gravity wave, but, for larger horizontal wave lengths, the effect of the boundary layer can become more important. The jump in the thermal conductivity in the boundary layer can give rise to overstability of the gravity wave in agreement with Souffrin and Spiegel (1967). The jump in the viscosity provides a self-catalytic action for the unstable flow if the viscosity is assumed to be the nonlinear turbulent viscosity due to the motion itself. The effect, however, is not strong enough to overcome the usual viscous damping. 相似文献
10.
NI JinLong LIU JunLai TANG XiaoLing ZHAO ChunQiang & ZENG QingDong 《中国科学:地球科学(英文版)》2014,57(4):600-613
The combination of field surveys with analysis of microstructure of tectonite and Electron Backscatter Diffraction (EBSD) on quartz fabric indicated that three periods of ductile shear events developed in the Paishanlou gold deposits and the E-W and NE-striking ductile shear zones were formed during each event. The E-W-striking ductile shear zone, accompanied by compressional and dextral shear slip, was shear-cut by the NE-striking shear zones, accompanied by compressional-sinistral shear slip and sinistral-normal shear slip, successively. An E-W-striking ductile shear zone developed at a deeper tectonic level and at middle- to high-temperatures, accompanied by abundant microstructures, including microlayering between a polycrystal quartz belt and mica, and quartz deformation was depended on cylinder (10-10) 〈a〉 or 〈c〉 glide. The development of an E-W-striking shear zone can be seen as a tectonic pattern in the region of the Paishanlou gold deposits of the collision between the Mongolian tectonic belt and the North Archean Craton from Suolun to the Linxi suture zone during the Indosinian. The NE-striking ductile shear zone developed approximately 160 Ma during the early Yianshanian at middle to shallow tectonic levels and at middle- to low-temperatures, accompanied by typical microstructures, including polycrystal quartz aggregation and quartz subgrain rotation recrystallization, etc., and quartz deformation was depended on prismatic (1011) 〈a〉 glide. The last ductile shear event around the NE-striking shear zone developed at low temperatures and shallow tectonic levels, yielding to a pre-existing NE-striking shear zone, accompanied by abundant microstructures, including low-temperature quartz grain boundary migration and bulging recrystallization. The last ductile shear movement may be related to lithosphere thinning and the destruction of the North China Craton from approximately 130-120 Ma, and this shear event resulted directly in the mineralization in the Paishanlou region. 相似文献
11.
G.S. Murty 《Physics of the Earth and Planetary Interiors》1975,11(1):65-79
Assuming that the interface of two loosely bonded half spaces permits a finite amount of slip, and that a simple linear relation exists between the prevailing shearing stress and the slip, a generalised secular equation for the Stoneley mode is derived and solved numerically. The two limiting cases of smooth interface and bonded interface are shown to be special cases of this general problem.For some range of values of the elastic constants of the half spaces, unattenuated and undispersed interfacial waves can propagate along the interface only when the interface is smooth or bonded. For the same combination of elastic properties of the half spaces, the loosely bonded interface will cause the interfacial wave to be attenuated and dispersed.The usefulness of this model in relation to the problem of attenuation and dispersion of elastic body waves is briefly discussed. 相似文献
12.
V. M. Pavlov 《Izvestiya Physics of the Solid Earth》2013,49(1):24-33
A new algorithm is proposed for calculating the complete synthetic seismograms from a point source in the form of the sum of a single force and a dipole with an arbitrary seismic moment tensor in a plane layered medium composed of homogenous elastic isotropic layers. Following the idea of (Alekseev and Mikhailenko, 1978), an artificial cylindrical boundary is introduced, on which the boundary conditions are specified. For this modified problem, the exact solution (in terms of the displacements and stresses on the horizontal plane areal element) in the frequency domain is derived and substantiated. The unknown depth-dependent coefficients form the motion-stress vector, whose components satisfy the known system of ordinary differential equations. This system is solved by the method that involves the matrix impedance and propagator for the vector of motion, as previously suggested by the author in (Pavlov, 2009). In relation to the initial problem, the reflections from the artificial boundary are noise, which, to a certain degree, can be suppressed by selecting a long enough distance to this boundary and owing to the presence of a purely imaginary addition to the frequency. The algorithm is not constrained by the thickness of the layers, is applicable for any frequency range, and is suitable for computing the static offset. 相似文献
13.
An exact theoretical formulation is presented for the analysis of a thin-walled pile embedded in an elastic half-space under vertically-incident P-wave excitation. In the framework of three-dimensional elastodynamics and a shell theory, the axisymmetrical wave-scattering problem is shown to be reducible to a set of Fredholm boundary integral equations. With the incorporation of the singular characteristics of the wave-induced contact load distributions into the solution scheme, a computational boundary element method is developed for a rigorous treatment of the seismic soil-structure interaction problem. Typical results for the dynamic contact load distributions, displacements, complex-valued foundation input motion functions, and resonant pile foundation response are included for direct engineering applications. 相似文献
14.
Review on localized boundary integral equation:Discrete wavenumber method for 2D irregular layers 
下载免费PDF全文
下载免费PDF全文 The pioneer study of simulating the wave field in media with irregular interface belongs to Aki and Larner.Since that many numerical methods on the subject have been developed,such as pure numerical techniques,ray method and boundary method.The boundary method based on boundary integral equation is a semi-analytical method which is suitable to modeling wave field induced by irregular border.According to the property of the applied Green's function the boundary methods can be sorted into space domain boundary method and wavenumber domain boundary method.For both of them it is necessary to solve a large equation,which means much computation is needed.Thus,it is difficult for the boundary methods to be applied in simulating wave field with high frequency or in large range.To develop a new method with less computation is meaningful.For this purpose,localized boundary integral equation,i.e.,discrete wavenumber method is proposed.It is rooted in the Bouchon-Campillo method,an important wavenumber domain boundary method.Firstly the force on interface is separated into two parts:one is on flat part and the other on irregular part of the interface.Then Fourier transform is applied to identify their relation,the unknown distributes only on irregular part.Consequently computation efficiency is dramatically improved.Importantly its accuracy is the same as that of Bouchon-Campillo. 相似文献
15.
A numerical scheme is developed in the paper for calculating torsional, vertical, horizontal, coupling and rocking impedances in frequency domain for axial-symmetric foundations embedded in layered media. In the scheme, the whole soil domain is divided into interior and exterior domains. For the exterior domain, the analytic solutions with unknown coefficients are obtained by solving three-dimensional (3D) wave equations in cylindrical coordinates satisfying homogeneous boundary conditions. For the interior domain, the analytical solutions are also obtained by solving the same 3D wave equations satisfying the homogeneous boundary conditions and the prescribed boundary conditions. The prescribed conditions are the interaction tractions at the interfaces between embedded foundation and surrounding soil. The interaction tractions are assumed to be piecewise linear. The piecewise linear tractions at the bottom surface of foundation will be decomposed into a series of Bessel functions which can be easily fitted into the general solutions of wave equations in cylindrical coordinates. After all the analytic solutions with unknown coefficients for both interior and exterior domains are found, the variational principle is employed using the continuity conditions (both displacements and stresses) at the interfaces between interior and exterior domains, interior domain and foundation, and exterior domain and foundation to find impedance functions. 相似文献
16.
断层滑动速率,地震重复时间和平均应力降 总被引：5，自引：0，他引：5
本文综述了国内外大地震复发间隔研究的现状，系统整理了中国大陆活断层滑运速率，已发现的古地震遗迹及古地震复发间隔，结果表明，在滑动速率较在的活断层发现的古地震复发间隔较短，Ｋａｎａｍｏｒｉ和Ａｌｌｅ（１９８５）研究了具有大范围重复时间（２０年到几千年）的板内大地震的震源参数，发现复发间隔较长的大地震具有较高的平均应力降，我们在实验室内研究了应变速率不同时固体围压三轴压缩下完整花岗岩破裂后的粘滑现象。 相似文献
17.
倾斜地层地震液化和滑移的有限元分析 总被引：1，自引：0，他引：1
应用饱和多孔介质动力学分析倾斜地层的土壤动力线性反应、液化和液化滑移问题，地下水位上的地层简化为单相介质层，饱和夹砂层看作是两相介质，水是可压缩的，采用双曲线非线性本构关系，考虑了砂土的剪胀性、刚度退化、滞回特性和土水相对运动等因素。基于土力学模型，建立了适用于分析非自由场地液化的动力方程组，基于是否考虑发生渗流问题，同时建立了两种离散形式：一种是以土骨架位移和水位移为未知量的矩阵方程，另一种是以土骨架位移、水位移和孔隙水压力为未知量的矩阵方程，初步分析了适用于多孔介质波动模拟的离散模型的人工边界问题，形成的方法将有助于问题的解决。 相似文献
18.
Interplate Earthquake Fault Slip During Periodic Earthquake Cycles in a Viscoelastic Medium at a Subduction Zone 总被引：1，自引：0，他引：1
K. Hirahara 《Pure and Applied Geophysics》2002,159(10):2201-2220
--A 2-D finite-element-method (FEM) numerical experiment of earthquake cycles at a subduction zone is performed to investigate the effect of viscoelasticity of the earth on great interplate earthquake fault slip. We construct a 2-D viscoelastic FEM model of northeast Japan, which consists of an elastic upper crust and a viscoelastic mantle wedge under gravitation overlying the subducting elastic Pacific plate. Instead of the dislocation model prescribing an amount of slip on a plate interface, we define an earthquake cycle, in which the plate interface down to a depth is locked during an interseismic period and unlocked during coseismic and postseismic periods by changing the friction on the boundary with the master-slave method. This earthquake cycle with steady plate subduction is periodically repeated to calculate the resultant earthquake fault slip.¶As simulated in a previous study (Wang, 1995), the amount of fault slip at the first earthquake cycle is smaller than the total relative plate motion. This small amount of fault slip in the viscoelastic medium was considered to be one factor explaining the small seismic coupling observed at several subduction zones. Our simulation, however, shows that the fault slip grows with an increasing number of repeated earthquake cycles and reaches an amount comparable to the total relative plate motion after more than ten earthquake cycles. This new finding indicates that the viscoelasticity of the earth is not the main factor in explaining the observed small seismic coupling. In comparison with a simple one-degree-of-freedom experiment, we demonstrate that the increase of the fault slip occurs in the transient state from the relaxed initial state to the stressed equilibrium state due to the intermittent plate loading in a viscoelastic medium. 相似文献
19.
Chihiro Hashimoto Eiichi Fukuyama Mitsuhiro Matsu’ura 《Pure and Applied Geophysics》2014,171(8):1705-1728
The generation of interplate earthquakes can be regarded as a process of tectonic stress accumulation and release, driven by relative plate motion. We completed a physics-based simulation system for earthquake generation cycles at plate interfaces in the Japan region, where the Pacific plate is descending beneath the North American and Philippine Sea plates, and the Philippine Sea plate is descending beneath the North American and Eurasian plates. The system is composed of a quasi-static tectonic loading model and a dynamic rupture propagation model, developed on a realistic 3-D plate interface model. The driving force of the system is relative plate motion. In the quasi-static tectonic loading model, mechanical interaction at plate interfaces is rationally represented by the increase of tangential displacement discontinuity (fault slip) across them on the basis of dislocation theory for an elastic surface layer overlying Maxwell-type viscoelastic half-space. In the dynamic rupture propagation model, stress changes due to fault slip motion on non-planar plate interfaces are evaluated with the boundary integral equation method. The progress of seismic (dynamic) or aseismic (quasi-static) fault slip on plate interfaces is governed by a slip- and time-dependent fault constitutive law. As an example, we numerically simulated earthquake generation cycles at the source region of the 1968 Tokachi-oki earthquake on the North American-Pacific plate interface. From the numerical simulation, we can see that postseismic stress relaxation in the asthenosphere accelerates stress accumulation in the source region. When the stress state of the source region is close to a critical level, dynamic rupture is rapidly accelerated and develops over the whole source region. When the stress state is much lower than the critical level, the rupture is not accelerated. This means that the stress state realized by interseismic tectonic loading essentially controls the subsequent dynamic rupture process. 相似文献
20.
The staggered grid finite-difference method is a powerful tool in seismology and is commonly used to study earthquake source
dynamics. In the staggered grid finite-difference method stress and particle velocity components are calculated at different
grid points, and a faulting problem is a mixed boundary problem, therefore different implementations of fault boundary conditions
have been proposed. Viriuex and Madariaga (1982) chose the shear stress grid as the fault surface, however, this method has
several problems: (1) Fault slip leakage outside the fault, and (2) the stress bump beyond the crack tip caused by S waves
is not well resolved. Madariaga et al. (1998) solved the latter problem via thick fault implementation, but the former problem remains and causes a new issue;
displacement discontinuity across the slip is not well modeled because of the artificial thickness of the fault. In the present
study we improve the implementation of the fault boundary conditions in the staggered grid finite-difference method by using
a fictitious surface to satisfy the fault boundary conditions. In our implementation, velocity (or displacement) grids are
set on the fault plane, stress grids are shifted half grid spacing from the fault and stress on the fictitious surface in
the rupture zone is given such that the interpolated stress on the fault is equal to the frictional stress. Within the area
which does not rupture, stress on the fictitious surface is given a condition of no discontinuity of the velocity (or displacement).
Fault normal displacement (or velocity) is given such that the normal stress on the fault is continuous across the fault.
Artificial viscous damping is introduced on the fault to avoid vibration caused by onset of the slip. Our implementation has
five advantages over previous versions: (1) No leakage of the slip prior to rupture and (2) a zero thickness fault, (3) stress
on the fault is reliably calculated, (4) our implementation is suitable for the study of fault constitutive laws, as slip
is defined as the difference between displacement on the plane of z = + 0 and that of z = − 0, and (5) cessation of slip is achieved correctly. 相似文献