Time domain method for calculating free field motion of a layered half-space subjected to obliquely incident body waves

In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complicates the computational area. In order to replace the complex frequency domain method, a time-domain method to calculate the free field motion of a layered half-space subjected to oblique incident body waves is developed in this paper. The new method decouples the equations of motion used in the finite element method and offers an interpolation formula of the free field motion. This formula is based on the fact that the apparent horizontal velocity of the free field motion is constant and can be calculated exactly. Both the theoretical analysis and numerical results demonstrate that the proposed method offers a high degree of accuracy.

Dynamic soil-tunnel interaction in layered half-space for incident P-and SV-waves

The dynamic soil-tunnel interaction is studied by the model of a rigid tunnel embedded in layered half-space, which is simplified as a single soil layer on elastic bedrock to the excitation of P- and SV-waves. The indirect boundary element method is used, combined with the Green＇ s function of distributed loads acting on inclined lines. It is shown that the dynamic characteristics of soil-tunnel interaction in layered half-space are different much from that in homoge- neous half-space, and that the mechanism of soil-tunnel interaction is also different much from that of soil-founda- tion-superstructure interaction. For oblique incidence, the tunnel response for in-plane incident SV-waves is com- pletely different from that for incident SH-waves, while the tunnel response for vertically incident SV-wave is very similar to that of vertically incident SH-wave.

2.5D scattering of incident plane SV waves by a canyon in layered half-space

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 fi.mctions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The free field displacements are added to the fictitious uniform moving loads induced displacements and the total response is obtained. Numerical calculations are performed for a canyon with homogenous and in one layer over bedrock. The effects of the thickness and stiffness of the layer on the amplification are studied and discussed.

Dynamic soil-tunnel interaction in layered half-space for incident plane SH waves

The dynamic soil-tunnel interaction is studied by indirect boundary element method （IBEM）, using the model of a rigid tunnel in layered half-space, which is simplified to a single soil layer on elastic bedrock, subjected to incident plane SH waves. The accuracy of the results is verified through comparison with the analytical solution. It is shown that soil-tunnel interaction in layered half-space is larger than that in homogeneous half-space and this interaction mechanism is essentially different from that of soil-foundation-superstructure interaction.

Diffraction of plane P waves by a canyon of arbitrary shape in poroelastic half-space (Ⅱ): Numerical results and discussion

This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffraction of waves. It is shown that the surface displacement amplitudes of the drained case are close to those of the undrained case, however, the surface displacement amplitudes of the dry case are very different from those of the saturated （either drained or undrained） cases. There are large phase shift between the dry case and the saturated cases, as well as slightly longer resultant wavelengths for the undrained case than those for the drained case and longer resultant wavelengths for the drained case than those for the dry case. For small porosity the surface displacement amplitudes for the saturated cases are almost identical to those for the dry case; while for large porosity, the effect of drainage condition becomes significant, and the surface displacement amplitudes for the undrained case are larger than those for the drained case. As the incident frequency increases, the effect of porosity becomes significant, and more significant for the undrained case than that for the drained case. As the porosity increases, the pore pressures increase significantly but their oscillations become smoother. As the incident frequency increases, the pore pressures become more complicated.

Diffraction of plane P waves by a canyon of arbitrary shape in poroelastic half-space (Ⅰ): Formulation

This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green＇s functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot＇s theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.

Diffraction of plane P waves around an alluvial valley in poroelastic half-space

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 fimctions 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 ficti- tious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, and the comparison between the degenerated solutions and available results in single-phase case. Finally, the nature of diffraction of plane P waves around an alluvial valley in poroelastic half-space is investigated in detail through nu- merical examples.

Diffraction of plane SV waves by a cavity in poroelastic half-space

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 fictitious wave sources close to the boundary of the cavity, and their magnitudes are determined by the boundary conditions. Verification of the accuracy is performed by： （1） checking the satisfaction extent of the boundary conditions, （2） comparing the degenerated solutions of a single-phased case with well- known solutions, and （3） examining the numerical stability of the solutions. The nature of diffraction of plane SV waves around a cavity in a poroelastic half-space is investigated by numerical examples.

Scattering of SH-wave by multiple circular cavities in half space

In this paper,an analytic method is developed to address steady SH-wave scattering and perform dynamic analysis of multiple circular cavities in half space.The scattered wave function used for scattering of SH-waves by multiple circular cavities,which automatically satisfies the stress-free condition at the horizontal surface,is constructed by applying the symmetry of the SH-wave scattering and the method of multi-polar coordinates system.Applying this scattered wave function and method of moving coordinates,the original problem can be transformed to the problem of SH-wave scattering by multiple circular cavities in the full space.Finally,the solution of the problem can be reduced to a series of algebraic equations and solved numerically by truncating the infinite algebraic equations to the finite ones.Numerical examples are provided for case with two cavities to show the effect of wave number,and the distances between the centers of the cavities and from the centers to the ground surface on the dynamic stress concentration around the cavity impacted by incident steady SH-wave.

Propagation of P-and S-waves in initially stressed gravitating half space

The present paper contributes in studying the phase velocities of P- and S-waves in a half space subjected to a compressive initial stress and gravity field. The density and acceleration due to gravity vary quadratically along the depth. The dispersion equation is derived in a closed form. It is shown that the phase velocities depend not only on the initial stress, gravity, and direction of propagation but also on the inhomogeneity parameter associated with the density and acceleration due to gravity. Various particular cases are obtained, and the results match with the classical results. Numerical investigations on the phase velocities of P- and S-waves against the wave number are made for various sets of values of the material parameters, and the results are illustrated graphically. The graphical user interface model is developed to generalize the effect.

Propagation of Love waves in non-homogeneous substratum over initially stressed heterogeneous half-space

The paper studies the propagation of Love waves in a non-homogeneous substratum over an initially stressed heterogeneous half-space. The dispersion equation of phase velocity is derived. The velocities of Love waves are calculated numerically as a function of kH and presented in a number of graphs, where k is the wave number, and H is the thickness of the layer. The case of Gibson＇s half-space is also considered. It is observed that the speed of Love waves is finite in the vicinity of the surface of the half-space and vanishes as the depth increases for a particular wave number. It is also observed that an increase in compressive initial stresses causes decreases of Love waves velocity for the same frequency, and the tensile initial stress of small magnitude in the half-space causes increase of the velocity.

Propagation of wave at the boundary surface of transversely isotropic thermoelastic material with voids and isotropic elastic half-space

The purpose of this research is to study the effect of voids on the surface wave propagation in a layer of a transversely isotropic thermoelastic material with voids lying over an isotropic elastic half-space. The frequency equation is derived after developing a mathematical model for welded and smooth contact boundary conditions. The dispersion curves giving the phase velocity and attenuation coefficient via wave number are plotted graphically to depict the effects of voids and anisotropy for welded contact boundary conditions. The specific loss and amplitudes of the volume fraction field, the normal stress, and the temperature change for welded contact are obtained and shown graphically for a particular model to depict the voids and anisotropy effects. Some special cases are also deduced from the present investigation.

铝合金薄板焊缝根部缺陷高阶半跨模式波定量检测

铝合金薄板焊接时会产生接近垂直的根部缺陷,且表面余高制约超声检测探头与楔块布置。结合常规半跨模式波的全聚焦方法有助于实现根部垂直缺陷轮廓表征与定量,但受余高影响,相控阵探头声场仍难以有效覆盖焊缝区域。该文考虑了薄板内的超声波多次反射,提出利用声束路径更长、声场覆盖范围更广的高阶半跨模式波进行缺陷轮廓重建。对存在和不存在余高的厚度6.0 mm铝合金板中端点深度3.0 mm、高度3.0 mm的根部裂纹,采用32阵元、中心频率5 MHz相控阵探头配合55◦横波楔块实施全矩阵捕捉和高阶半跨模式波成像。仿真和实验结果表明,缺陷轮廓可以被完整呈现,且缺陷高度定量相对误差均不超过5.0%。此外,需依据板材厚度合理选择半跨模式波阶数,板厚越薄,适用的半跨模式波阶数越高。

基础激励下的振动台变压器虚拟试验研究

针对变压器因振动台基础激励引起的强度问题,本文作者应用虚拟试验技术和振动学理论研究两种构型的变压器的动力学特性。基于基础激励条件下的虚拟试验模型,应用模态分析理论,通过分析变压器的振型和固有频率,得到其薄弱环节。采用模态叠加法和冲击动力学理论,并结合瑞利阻尼模型,得到两种构型的变压器在正弦信号激励和半正弦波冲击时的应力和变形云图。研究结果显示:两种构型的变压器的侧向均为结构的最薄弱方向。变压器的长边底座单边长度减短25mm,对其基频和正弦信号激励下共振点处的变形有重大影响,但对其高频、振型、正弦信号激励下共振点处的应力以及半正弦波冲击下的动强度影响较小。在半正弦波冲击下,高量值的冲击过载对两种构型变压器的动力学特性有重大影响。该研究结果可为变压器优化设计提供理论基础。

Experimental investigation on the optical remote sensing images of internal solitary waves with a smooth surface

The parameter inversion of internal solitary waves (ISWs) based on optical remote sensing images is a key work. A new approach is proposed and demonstrated for simulating the optical remote sensing images of ISWs with a smooth surface in the laboratory. An optical remote sensing simulation system used to detect ISWs is constructed by a two-dimensional ISW flume, a LED (light emitting diode) light source and two CCD (charge coupled device) cameras. The optical remote sensing images of the horizontal surface and ISWs propagation images of a vertical side are detected simultaneously, which aims to explore the response of optical remote sensing corresponding to ISWs with the smooth surface. The results show that during the propagation of ISWs, dark pattern images are obtained by CCD 1 camera. The characteristics of the dark patterns vary along with the incident angle of the light source. The characteristic parameters of the optical remote sensing images correspond to the wave factors of vertical profiles. The experiment also shows a positive correlation between the dark pattern width and the half wave width under different amplitudes of ISWs. The system has the advantages of clear phenomenon and high repeatability, which provides the scientific basis for quantitative investigation on imaging mechanism of ISW by optical remote sensing.

Analytical solutions for dynamic pressures of coupling fluid-solid-porous medium due to P wave incidence

Wave reflection and refraction in layered media is a topic closely related to seismology,acoustics,geophysics and earthquake engineering.Analytical solutions for wave reflection and refraction coefficients in multi-layered media subjected to P wave incidence from the elastic half-space are derived in terms of displacement potentials.The system is composed of ideal fluid,porous medium,and underlying elastic solid.By numerical examples,the effects of porous medium and the incident wave angle on the dynamic pressures of ideal fluid are analyzed.The results show that the existence of the porous medium,especially in the partially saturated case,may significantly affect the dynamic pressures of the overlying fluid.

Fast and accurate analysis of electromagnetic scattering from targets situated in dielectric half space

Electromagnetic scattering from targets situated in half space is solved by applying fast inhomogeneous plane wave algorithm combined with a tabulation and interpolation method. The integral equation is set up based on derivation of dyadic Green＇s functions in this environment. The coupling is divided into nearby region and well-separated region by grouping. The Green＇s function can be divided into two parts： primary term and reflected term. In the well-separated region, the two terms are both expressed as Sommerfeld integral, which can be accelerated by deforming integral path and taking interpolation and extrapolation. For the nearby region, the direct Sommerfeld integral makes the filling of impedance matrix time-expensive. A tabulation and interpolation method is applied to speed up this process. This infinite integral is pre-computed in sampling region, and a two-dimensional table is then set up. The impedance elements can then be obtained by interpolation. Numerical results demonstrate the accuracy and efficiency of this algorithm.

Green's function of multi-layered poroelastic half-space for models of ground vibration due to railway traffic

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.

Frequency domain fundamental solutions for a poroelastic half-space

In frequency domain, the fundamental solutions for a poroelastic half-space are re-derived in the context of Biot＇s theory. Based on Biot＇s theory, the governing field equations for the dynamic poroelasicity are established in terms of solid displacement and pore pressure. A method of potentials in cylindrical coordinate system is proposed to decouple the homogeneous Biot＇s wave equations into four scalar Helmholtz equations, and the general solutions to these scalar wave equations are obtained. After that, spectral Green＇s functions for a poroelastic full-space are found through a decomposition of solid displacement, pore pressure, and body force fields. Mirror-image technique is then applied to construct the half-space fundamental solutions.Finally, transient responses of the half-space to buried point forces are examined.

The influence of divergence angle on the deposition of neutral chromium atoms using a laser standing wave

The characteristics of neutral chromium atoms in the standing wave field are discussed. Based on a semi-classical model, the motion equation of neutral atoms in the laser standing wave field is analyzed, and the trajectories of the atoms are obtained by simulations with the different divergence angles of the atomic beam. The simulation results show that the full width at half maximum （FWHM） of the stripe is 2.75 nm and the contrast is 38.5 ： 1 when the divergence angle equals 0 mrad, the FWHM is 24.1 nm and the contrast is 6.8：1 when the divergence angle equals 0.2 mrad and the FWHMs are 58.6 and 137.8 nm, and the contrasts are 3.3 ： 1 and 1.6 ： i when the divergence angles equal 0.5 and 1.0 mrad, respectively.