首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 657 毫秒
1.
Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green’s functions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of fictitious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary...  相似文献   

2.
Diffraction of plane SV waves by a cavity in poroelastic half-space   总被引：2，自引：0，他引：2  
This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green’s functions of compressive and shear wave sources are derived based on Biot’s theory. The scattered waves are constructed using fi ctitious wave sources close to the boundary of the cavity, and their magnitudes are determined by the boundary conditions. Verifi cation of the accuracy is performed by: (1) checking the satisfaction extent of the boundary c...  相似文献   

3.
This paper presents 2.5D scattering of incident plane SH waves by a canyon in layered half-space by the indirect boundary element method (IBEM). The free field response is carried out to give the displacements and stresses on the line which forms boundary of the canyon. The fictitious uniform moving loads are applied to the same line to calculate the Green’s functions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The displacements due to the free ...  相似文献   

4.
This paper presents a 2.5D scattering of incident plane SV waves by a canyon in a layered half-space by using the indirect boundary element method (IBEM). A free field response analysis is performed to provide the displacements and stresses on the boundary of the canyon where fictitious uniform moving loads are applied to calculate the Green’s functions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The free field displacements are added to the fic...  相似文献   

5.
The indirect boundary element method （IBEM） is used to study three-dimensional scattering of obliquely incident plane SH waves by an alluvial valley embedded in a layered half-space. The free-field response of the layered half-space is calculated by the direct stiffness method, and dynamic Green＇s functions of moving distributed loads acting on inclined lines in a layered half-space are calculated to simulate the scattering wave field. The presented method yields very accurate results since the three-dimensional dynamic stiffness matrix is exact and the moving distributed loads can act directly on the valley boundary without singularity. Numerical results and analyses are performed for amplification of obliquely incident plane SH waves around an alluvial valley in a uniform half-space and in single layer over half-space. The results show that the three-dimensional responses are distinctly different from the two-dimensional responses, and the displacement amplitudes around alluvial valleys in a uniform haft-space are obviously different from those in a layered half-space.  相似文献   

6.
This paper reports a series solution of wave functions for two-dimensional scattering and diffraction of plane SH waves induced by a symmetrical V-shaped canyon with different shape ratios. A half-space with a symmetrical V-shaped canyon is divided into two sub-regions by using a circular-arc auxiliary boundary. The two sub-regions are represented by global and local cylindrical coordinate systems, respectively. In each coordinate system, the wave field satisfying the Helmholtz equation is represented by the separation of variables method, in terms of the series of both Bessel functions and Hankel functions with unknown complex coefficients. Then, the two wave fields are described in the local coordinate system using the Graf addition theorem. Finally, the unknown coefficients are sought by satisfying the continuity conditions of the auxiliary boundary. To consider the phase characteristics of the wave scattering, a parametric analysis is carried out in the time domain by assuming an incident signal of the Ricker type. Surface and subsurface transient responses demonstrate the characteristics and mechanisms of wave propagating and scattering.  相似文献   

7.
The diffraction of elastic waves by a sedimentary valley in a homogeneous elastic half-space is studied in this paper. The sediment-filled valley is composed of a fluid layer over a soft soil deposit whose characteristics may be significant and should be carefully considered when designing long span bridges with high piers. The method of analysis adopted in the paper is to decompose the problem into an interior region and an exterior region. In the exterior region, the scattered wave fields are constructed with the linear combinations of two independent sets of Lamb’s singular solutions, i.e., the integralsolutions for two concentrated surface loads in two directions; and their derivatives are used to represent the scattered wave fields. A technique is proposed to calculate the integrals in the wave-number domain based on the method of steepest descent. For the interior region, the wave fields for the fluid layer and soft soil deposit are expressed in terms of wave functions which satisfy the equation of motion. The continuity condition at the interface of the media is satisfied in the least square sense. The effects of geometric topography, soil amplification and fluid layer subject to different types of incident harmonic plane waves are analyzed and discussed..  相似文献   

8.
The indirect boundary element method is used to study the 3D dynamic response of an infinitely long alluvial valley embedded in a saturated layered half-space for obli- quely incident SV waves. A wave-number transform is first applied along the valley＇s axis to reduce a 3D problem to a 2D plane strain problem. The problem is then solved in the section perpendicular to the axis of the valley. Finally, the 3D dynamic responses of the valley are obtained by an inverse wave-number transform. The validity of the method is con- firmed by comparison with relevant results. The differences between the responses around the valley embedded in dry and in saturated poroelastic medium are studied, and the effects of drainage conditions, porosity, soil layer stiffness, and soil layer thickness on the dynamic response are dis- cussed in detail resulting in some conclusions.  相似文献   

9.
The theory of elastic wave scattering is a fundamental concept in the study of elastic dynamics and wave motion,and the wave function expansion technique has been widely used in many subjects.To supply the essential tools for solving wave scattering problems induced by an eccentric source or multi-sources as well as multi-scatters,a whole-space transform formula of cylindrical wave functions is presented and its applicability to some simple cases is demonstrated in this study.The transforms of wave functions in cylindrical coordinates can be classifi ed into two basic types: interior transform and exterior transform,and the existing Graf’s addition theorem is only suitable for the former.By performing a new replacement between the two coordinates,the exterior transform formula is fi rst deduced.It is then combined with Graf’s addition theorem to establish a whole-space transform formula.By using the whole-space transform formula,the scattering solutions by the sources outside and inside a cylindrical cavity are constructed as examples of its application.The effectiveness and advantages of the whole-space transform formula is illustrated by comparison with the approximate model based on a large cycle method.The whole-space transform formula presented herein can be used to perform the transform between two different cylindrical coordinates in the whole space.In addition,its concept and principle are universal and can be further extended to establish the coordinate transform formula of wave functions in other coordinate systems.  相似文献   

10.
Sediment movement in the wave boundary layer above a mobile sediment bed is complex.A velocity formula for the boundary layer is proposed for sheet flow induced by asymmetric waves above a mobile sediment bed.The formula consists of a free stream velocity and a defect function which contains a phase-lead,boundary layer thickness and mobile sediment bed.Phase-lag of sediment movement is considered in the formula for the mobile sediment bed.The formula needs six dependent variables about asymmetric wave and sediment characteristics.Asymmetry effects on parameters(orbital amplitude,roughness height,bed shear stress,and boundary layer thickness)are properly considered such that the formula can yield velocity differences among onshore,offshore,acceleration,and deceleration stages.The formula estimates the net boundary layer velocity resulting from the mobile sediment bed and asymmetric boundary layer thickness.In addition,a non-constant phase-lead also contributes to the net boundary layer velocity in asymmetric oscillatory sheet flow.Results of the formula are as good as that of a two-phase numerical model.Sheet flow transport induced by asymmetric waves,and the offshore net sediment transport rate with a large phase-lag under velocity-skewed waves,can be adequately estimated by the formula with a power sediment concentration function.  相似文献   

11.
Based on one type of practical Biot＇s equation and the dynamic-stiffness matrices ofa poroelastic soil layer and half-space, Green＇s functions were derived for unitformly distributed loads acting on an inclined line in a poroelastie layered site. This analysis overcomes significant problems in wave scattering due to local soil conditions and dynamic soil-structure interaction. The Green＇s functions can be reduced to the case of an elastic layered site developed by Wolf in 1985. Parametric studies are then carried out through two example problems.  相似文献   

12.
The dynamic stiffness method combined with the Fourier transform is utilized to derive the in-plane Green’s functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic(TI)half-space.The loaded layer is fixed to obtain solutions restricted in it and the corresponding reactions forces,which are then applied to the total system with the opposite sign.By adding solutions restricted in the loaded layer to solutions from the reaction forces,the global solutions in the wavenumber domain are obtained,and the dynamic Green’s functions in the space domain are recovered by the inverse Fourier transform.The presented formulations can be reduced to the isotropic case developed by Wolf(1985),and are further verified by comparisons with existing solutions in a uniform isotropic as well as a layered TI halfspace subjected to horizontally distributed loads which are special cases of the more general problem addressed.The deduced Green’s functions,in conjunction with boundary element methods,will lead to significant advances in the investigation of a variety of wave scattering,wave radiation and soil-structure interaction problems in a layered TI site.Selected numerical results are given to investigate the influence of material anisotropy,frequency of excitation,inclination angle and layered on the responses of displacement and stress,and some conclusions are drawn.  相似文献   

13.
Antiplane response of isosceles triangular hill to incident SH waves   总被引：6，自引：1，他引：6  
In this paper, antiplane response of an isosceles triangular hill to incident SH waves is studied based on the method of complex function and by using moving coordinate system. The standing wave function, which can satisfy the governing equation and boundary condition, is provided. Furthermore, numerical examples are presented; the influences of wave number and angle of the incident waves and the angle of the hill‘s peak on ground motion are discussed.  相似文献   

14.
The nonsplitting perfectly matched layer （NPML） absorbing boundary condition （ABC） was first provided by Wang and Tang （2003） for the finite-difference simulation of elastic wave propagation in solids. In this paper, the method is developed to extend the NPML to simulating elastic wave propagation in poroelastic media. Biot＇s equations are discretized and approximated to a staggered-grid by applying a fourth-order accurate central difference in space and a second-order accurate central difference in time. A cylindrical twolayer seismic model and a borehole model are chosen to validate the effectiveness of the NPML. The results show that the numerical solutions agree well with the solutions of the discrete wavenumber （DW） method.  相似文献   

15.
Active seismic sources are critical for obtaining high resolution images of the subsurface. For active imaging in urban areas, environment friendly and green seismic sources are required. In present work, we introduce a new type of green active source based on the gaseous detonation of methane and oxygen. When fired in a closed container, the chemical reaction, i.e. gaseous detonation, will produce high pressure air over 150MPa. Seismic waves are produced when high pressure air is quickly released to impact the surroundings. The first field experiment of this active source was carried out in December, 2017 in Jingdezhen, Jiangxi Province, where a series of active sources were excited to explore their potential in mine exploration. In current work, we analyzed the seismic waves recorded by near-field accelerators and a dense short-period seismic array and compared them with those from a mobile airgun source, another kind of active source by releasing high pressure air into water. The results demonstrate that it can be used for high resolution near surface imaging. Firstly, the gaseous detonation productions are harmless CO2 and water, making it a green explosive source. Secondly, the dominant seismic frequencies are 10-80Hz and a single shot can be recorded up to 15km, making it suitable for local structure investigations. Thirdly, it can be excited in vertical wells, similar to traditional powder explosive sources. It can also act as an additional on-land active source to airgun sources, which requires a suitable water body as intermediate media to generate repeating signals. Moreover, the short duration and high frequency signature of the source signals make it safe with no damage to nearby buildings. These make it convenient to excite in urban areas. As a new explosive source, the excitation equipment and conditions, such as gas ratio, sink depth and air-releasing directions, need further investigation to improve seismic wave generation efficiency.  相似文献   

16.
Based on one type of practical Biot's equation and the dynamic-stiffness matrices of a poroelastic soil layer and half-space, Green's functions were derived for uniformly distributed loads acting on an inclined line in a poroelastic layered site. This analysis overcomes significant problems in wave scattering due to local soil conditions and dynamic soil-structure interaction. The Green's functions can be reduced to the case of an elastic layered site developed by Wolf in 1985. Parametric studies are then carried out through two example problems.  相似文献   

17.
This study proposes a Green’s function,an essential representation of water-saturated ground under moving excitation,to simulate ground borne vibration from trains.First,general solutions to the governing equations of poroelastic medium are derived by means of integral transform.Secondly,the transmission and reflection matrix approach is used to formulate the relationship between displacement and stress of the stratified ground,which results in the matrix of the Green’s function.Then the Green’s function is combined into a train-track-ground model,and is verified by typical examples and a field test.Additional simulations show that the computed ground vibration attenuates faster in the immediate vicinity of the track than in the surrounding area.The wavelength of wheel-rail unevenness has a notable effect on computed displacement and pore pressure.The variation of vibration intensity with the depth of ground is significantly influenced by the layering of the strata soil.When the train speed is equal to the velocity of the Rayleigh wave,the Mach cone appears in the simulated wave field.The proposed Green’s function is an appropriate representation for a layered ground with shallow ground water table,and will be helpful to understand the dynamic responses of the ground to complicated moving excitation.  相似文献   

18.
用尾波干涉法监测介质波速变化研究进展   总被引：1，自引：1，他引：0  
Multiply scattered waves are sensitive to media changes owing to the effect of repeated sampling,superposition and amplification. Based on this characteristic,small-medium changes could be detected by using coda wave interferometry. In recent years,coda wave interferometry has been widely used in estimating velocity variation with high precision in areas such as seismology and non-destructive testing. This paper systematically presents the principle and research status of coda wave interferometry,and especial focus is placed on the research of media velocity variations by using repeating earthquakes,artificial sources,and ambient noise. Applications of coda wave interferometry can contribute to the more subtle understanding of dynamic evolution process in the medium.  相似文献   

19.
Applying a fully nonlinear numerical scheme with second-order temporal and spatial precision,nonlinear interactions of gravity waves are simulated and the matching relationships of the wavelengths and frequencies of the interacting waves are discussed.In resonant interactions,the wavelengths of the excited wave are in good agreement with the values derived from sum or difference resonant conditions,and the frequencies of the three waves also satisfy the matching condition.Since the interacting waves obey the resonant conditions,resonant interactions have a reversible feature that for a resonant wave triad,any two waves are selected to be the initial perturbations,and the third wave can then be excited through sum or difference resonant interaction.The numerical results for nonresonant triads show that in nonresonant interactions,the wave vectors tend to approximately match in a single direction,generally in the horizontal direction.The frequency of the excited wave is close to the matching value,and the degree of mismatching of frequencies may depend on the combined effect of both the wavenumber and frequency mismatches that should benefit energy exchange to the greatest extent.The matching and mismatching relationships in nonresonant interactions differ from the results of weak interaction theory that the wave vectors are required to satisfy the resonant matching condition but the frequencies are permitted to mismatch and oscillate with amplitude of half the mismatching frequency.Nonresonant excitation has an irreversible characteristic,which is different from what is found for the resonant interaction.For specified initial primary and secondary waves,it is difficult to predict the values of the mismatching wavenumber and frequency for the excited wave owing to the complexity.  相似文献   

20.
The mechanism of the effects of the upwelling mean on the ENSO event mature phase locking is ex-amined by using a mixed-mode model. The results show that the positive feedback process of the ef-fects of the seasonal variation of the upwelling mean on the Kelvin wave is the mechanism of the locking of the event mature phase to the end of the calendar year. The memory of the Rossby waves for the sign-shifting of the sea surface temperature anomaly from positive to negative 6 months before the cold peak time is the other mechanism of the locking of the La Nia event mature phase to the end of the calendar year. The results here are different from previous ones which suggest that the balance between cold and warm trends of sea surface temperature anomaly is the mechanism involved. The cold trend is caused by the upwelling Kelvin wave from upwelling Rossby wave reflected at the western boundary, excited by the westerly anomaly stress over the central Pacific and amplified by the seasonal variation of the coupled strength in its way propagating westward. The warm trend is caused by the Kelvin wave forced by the western wind stress over the middle and eastern equatorial Pacific. The cause of the differences is due to the opposite phase of the seasonal variation of the upwelling mean to that in the observation and an improper parameterization scheme for the effects of the seasonal varia-tion of the upwelling mean on the ENSO cycle in previous studies.  相似文献   

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

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