Surface motion of a semi-elliptical hill for incident plane SH waves

A closed-form analytical solution of surface motion of a semi-elliptical cylindrical hill for incident plane SH waves is presented. Although some previous analytical work had already dealt with hill topography of semi-circular and shallow circular, our work aims at calculating surface motion of very prolate hill for high incident frequency, and explaining the special vibrating properties of very prolate hill. Accuracy of the solution is checked by boundary conditions, numerical results for surface motion of oblate and prolate hills are calculated, and some conclusions are obtained.     

Quasilinear infiltration from an elliptical cavity

We develop analytic solutions to the linearized steady-state Richards equation for head and total flowrate due to an elliptic cylinder cavity with a specified pressure head boundary condition. They are generalizations of the circular cylinder cavity solutions of Philip [Philip JR. Steady infiltration from circular cylindrical cavities. Soil Sci Soc Am J 1984;48:270–8]. The circular and strip sources are limiting cases of the elliptical cylinder solution, derived for both horizontally- and vertically-aligned ellipses. We give approximate rational polynomial expressions for total flowrate from an elliptical cylinder over a range of sizes and shapes. The exact elliptical solution is in terms of Mathieu functions, which themselves are generalizations of and computed from trigonometric and Bessel functions. The required Mathieu functions are computed from a matrix eigenvector problem, a modern approach that is straightforward to implement using available linear algebra libraries. Although less efficient and potentially less accurate than the iterative continued fraction approach, the matrix approach is simpler to understand and implement and is valid over a wider parameter range.     

Dynamic stress concentration and scattering of SH-wave by interface elliptic cylindrical cavity

In this paper, the dynamic stress concentration and scattering of SH-waves by bi-naterial structures that possess an interface elliptic cavity are investigated. First, by using the complex function method, the Green‘s function is constructed. This yields the solution of the displacement field for an elastic half space with a semi-elliptic canyon impacted by an anti-plane harmonic line source loading on the horizontal surface. Then, the problem is divided into an upper and lower half space along the horizontal interface, regarded as a harmony model. In order to satisfythe integral continuity condition,the unknown anti-plane forces are applied to the interface. The integral equations with unknown forces can be established through the continuity condition, and after transformation, the algebraic equations are solved numerically. Finally, the distribution of the dynamic stress concentration factor (DSCF) around the elliptic cavity is given and the effect of different parameters on DSCF is discussed.     

层状半空间中周期分布凸起地形对平面SH波的散射

提出了一种新的以层状半空间中周期分布斜线荷载动力格林函数为基本解的间接边界元方法,研究了周期分布凸起地形对平面SH波的散射问题.方法将散射波场分解为凸起内部散射波场和凸起外部散射波场.凸起内部散射波场通过在凸起闭合边界上施加虚拟斜线荷载产生的动力响应来模拟,而凸起外部散射波场则通过在凸起与半空间交界面上施加虚拟周期分布斜线荷载产生的动力响应来模拟.周期分布斜线荷载动力格林函数的引入,使得本文方法仅需针对一个凸起进行边界单元的离散和求解,便可完成问题的求解,避免了通过截断无限边界求解而引入的误差,方法具有较高精度的同时显著降低了求解自由度.文中通过与已有结果的比较,验证了方法的正确性,并以均匀半空间和基岩上单一土层中周期分布凸起为例进行了数值计算分析.研究表明,凸起间距对凸起地形间的动力相互作用有着显著的影响,同时层状半空间中周期分布凸起地形对SH波的散射与均匀半空间情况也有着显著的差别.

     

Scattering of shear waves by an elliptical cavity in a radially inhomogeneous isotropic medium

Complex function and general conformal mapping methods are used to investigate the scattering of elastic shear waves by an elliptical cylindrical cavity in a radially inhomogeneous medium. The conformal mappings are introduced to solve scattering by an arbitrary cavity for the Helmholtz equation with variable coefficient through the transformed standard Helmholtz equation with a circular cavity. The medium density depends on the distance from the origin with a power-law variation and the shear elastic modulus is constant. The complex-value displacements and stresses of the inhomogeneous medium are explicitly obtained and the distributions of the dynamic stress for the case of an elliptical cavity are discussed. The accuracy of the present approach is verified by comparing the present solution results with the available published data. Numerical results demonstrate that the wave number, inhomogeneous parameters and different values of aspect ratio have significant influence on the dynamic stress concentration factors around the elliptical cavity.     

线源荷载对半圆形凸起圆形夹杂附近浅埋圆孔的动力作用

采用复变函数法和Green函数法研究了在水平界面承受出平面线源荷载时弹性半空间内半圆形凸起的圆柱形弹性夹杂对浅埋圆孔的动力作用.该问题采用“分区”、“契合”思想求解.首先,将整个求解区域分割成两部分,其一为含半圆形凹陷和圆孔的弹性半空间,其二为圆柱形弹性夹杂;其次,构造满足含半圆形凹陷半空间水平界面应力自由和圆孔边界应力自由的散射波,构造满足圆形夹杂上半表面应力自由下半表面应力连续条件的驻波和散射波;最后,在两个区域的“公共边界”上实施“契合”,满足公共边界处位移和应力的连续性条件,同时满足圆孔边界应力自由的边界条件,建立起求解该问题的无穷代数方程组,并就具体算例分析讨论了浅埋圆孔边缘动应力集中系数(DSCF)的数值结果.结果表明:圆柱形弹性夹杂的“软”、“硬”对浅埋圆孔孔边动应力集中系数有不同的影响.     

P波入射时含圆柱形空腔板弹性波散射与动应力集中

采用波函数展开法给出了板的横截面处孔洞在平面P波入射下动应力集中问题的解。本文使用大圆弧假定法来满足板中的自由表面的边界条件。数值结果表明,上表面和孔洞间的距离与下表面和孔洞间的距离都对结果有影响。如果孔洞半径相对板的厚度很小,解近似等于孔洞全空间解。所给出的解是半解析解。因此能被用来验证诸如BEM,FEM,FDM等数值法。     

Scattering of SH waves by a scalene triangular hill with a shallow cavity in half-space

In this study, a theoretical approach is used to investigate the scattering problem of circular holes under a scalene triangle on the surface. The wave displacement function is obtained by solving the Helmholtz equation that meets the zero-stress boundary conditions by adopting the method of separation of variables. Based on the complex function,multi-polar coordinate method, and region-matching technique, algebraic equations are established at auxiliary boundaries and free boundaries conditions...     

Dynamic stress concentration near a fluid-filled permeable borehole induced by general modal vibrations of an internal cylindrical radiator

Stress distribution in the vicinity of a permeable cylindrical cavity surface (borehole wall) arising due to modal vibrations of an internal cylindrical radiator of infinite extent is studied. Biot phenomenological model is used to represent the behavior of sound in the fluid-saturated elastic porous medium and closed-form solution in the form of an infinite series is developed. A numerical example for the infinite cylindrical surface excited in vibrational modes of zeroth and first order while immersed in a water-filled cavity embedded within a water-saturated Ridgefield Sandstone environment is presented and several limiting cases are examined. Effects of axial and radial vibration frequencies, porosity, frame stiffness, and interface permeability condition on stress distribution at the borehole surface are presented and discussed.     

Dynamic response of twin circular tunnels due to incident SH-waves

The steady-state solution for two parallel under ground tunnels of circular cross-section subjected to incident plane harmonic SH-waves is obtained in closed form, by using the method of wave function expansion. The image technique is employed to account for the reflection of waves at the ground surface. Numerical studies are carried out to determine the influence of spacing between the tunnels on the shear stresses in concrete and steel linings.     

地下圆形衬砌洞室在平面P波入射下的动应力集中:三维级数解

利用波函数展开法,研究了地下圆形衬砌洞室在入射平面P波作用下的动应力集中问题并给出了三维级数解。研究表明:(1)衬砌刚度对环向与轴向动应力集中系数都具有重要影响,刚性衬砌、无衬砌和柔性衬砌三种情况下的动应力集中系数在空间上的分布基本相同,但刚性衬砌情况的动应力集中系数最大,无衬砌情况次之,柔性衬砌情况最小,波阵面夹角入射时1δ=45°,环向与轴向动应力集中系数最大可分别达到16.62与8.53;(2)随着入射频率的增大,环向动应力集中系数在多数情况下逐渐减小;(3)波入射角度对动应力集中系数也有很大影响,随着入射波与y轴夹角2δ的增加,环向动应力集中系数总体上逐渐减小,而轴向动应力集中系数多数情况下逐渐增大。     

Scattering of plane sh waves by a semi-elliptical canyon

In this paper we analyse the two-dimensional scattering and diffraction of plane SH waves by a semi-elliptical canyon. The exact series solution of the problem, for general angle of incidence of the plane SH waves, has been used to examine the dependence of surface amplifications inside and near the canyon. The nature of ground motion has been found to depend on two key parameters: (a) The angle of incidence. (b) The ratio of the canyon width to the wave length of incident SH waves. For short incident waves surface displacement amplitudes change rapidly from one point to another, while for the long waves and shallow canyons displacement amplitudes display only minor departure from the uniform half-space amplification of 2. For shallow canyons and long incident waves, the angle of incidence introduces only minor changes into the overall behaviour of surface amplitudes. For deep canyons and nearly grazing incidences, a prominent shadow zone is realized behind the canyon.     

Dynamic stress concentration in pre-stressed poroelastic media due to moving punch influenced by shear wave

During the occurrence of earthquake, the shear wave propagates in the rocks present inside/at the Earth’s crust. The propagation of shear wave may lead to the progression of punch present inside the rock medium. As a result of this, substantial stress accumulated at the vicinity of propagating punch inside rock medium which significantly affects the stability of various geological and human-made structure and, hence, may cause failure of structure. Therefore, the analysis of stress concentration at the vicinity of punch moving due to shear wave propagation has become prominent in the area of seismology. In the present paper, an analytical perspective has been employed to discuss the influence of velocity of moving punch associated with the propagation of shear wave on developed dynamic stress concentration (DSC) in three types of pre-stressed vertical transversely isotropic (VTI) poroelastic media viz. granite (an igneous rock); sandstone (a sedimentary rock); and marble (a metamorphic rock). The closed form expression of DSC for the force of constant intensity has been derived with the aid of Weiner-Hopf technique along with Galilean and two-sided Fourier integral transformations. The noticeable influence of different affecting parameters (viz. velocity of moving punch associated with the shear wave propagation, horizontal compressive/tensile initial stresses, vertical compressive/tensile initial stress, porosity, and anisotropy parameter) on dynamic stress concentration has also been reported. Numerical computation and graphical illustrations have been carried out for the aforementioned three different types of porous rocks to investigate the profound impact of affecting parameters on DSC. Moreover, some noteworthy peculiarities have also been derived from the obtained expression of dynamic stress concentration.     

Reflection and transmission of plane SH-waves in an initially stressed inhomogeneous anisotropic magnetoelastic medium

This study investigates the reflection and transmission of plane SH-waves in two semi-infinite anisotropic magnetoelastic media. The lower half-space is considered as initially stressed and inhomogeneous. The density of lower half-space is taken exponentially varying with depth. The solutions for half-spaces are obtained analytically. The expressions for reflection and transmission coefficient are obtained in the closed form subject to continuity conditions at the interfaces of anisotropic magnetoelastic half-spaces and the Snell’s law. It is found that these coefficients depend on the initial stress, inhomogeneity parameter, the magnetoelastic coupling parameter, and the angle at which wave crosses the magnetic field of the half-spaces. Numerical computations are performed for these coefficients for a specific model of two different anisotropic magnetoelastic half-spaces. The numerical results are illustrated by the graph of reflection and transmission coefficient versus the angle of incidence. In general, as the initial stress increases the reflection and transmission coefficient increases, the affect is more prominent for more than 10 GPa. Inhomogeneity in the density of the material also increases the reflection and transmission coefficient. The anisotropic magnetoelastic parameter and the angle at which the wave crosses the magnetic field for both the half-spaces have a quite significant effect on the reflection and transmission coefficient.     

Dynamic displacements and stresses in the vicinity of a cylindrical cavity embedded in a half-space

The two-dimensional response of a viscoelastic half-space containing a buried, unlined, infinitely long cylindrical cavity of circular cross-section subjected to harmonic plane SH, P, SV and Rayleigh waves is obtained by an indirect boundary integral method based on the two-dimensional Green's functions for a viscoelastic half-space. An extensive critical review of the existing numerical results obtained by other techniques is presented together with some new numerical results describing the motion on the ground surface and the motion and stresses on the wall of the cavity for P, SV, SH and Rayleigh waves.     

SH波入射时多个含孔半圆形凸起附近的多个浅埋圆孔的动应力分析

文中使用"契合"、复变函数法和移动坐标的方法,给出了SH波入射条件下多个含孔半圆形凸起地形附近多个浅埋圆形孔洞动力分析问题的解答.将模型分成两个区域,在区域Ⅰ和区域Ⅱ中分别构造波函数,最后在两个区域的公共边界上实施"契合",即满足两个区域的公共边界上位移、应力的连续条件.通过算例给出了浅埋圆孔动应力与入射波数、入射波角度的变化规律.     

SH波入射时半圆形凸起地形附近浅埋圆形衬砌结构的动应力分析

本文给出了地下圆形衬砌结构与地面上的半圆形凸起地形对垂直于地面入射的SH波散射问题的解答。方法是将求解区域分割成两部分。其一为包含半圆形凸起地形在内的圆形区域Ⅰ，其二为带有一个半圆形凹陷和一个圆形衬砌结构的弹性半空间Ⅱ，半圆形凹陷部分为其公共边界，在区域Ⅰ和Ⅱ中分别构造其位移解，然后再通过移动坐标，使其满足“公共边界”上的条件和地下衬砌的边界条件，建立起求解该问题的无穷代数方程组。最后，本文给出了算例，并讨论了数值结果，给出了圆形衬砌结构周边上的动应力集中系数变化规律。     

Diffraction of SH-waves by topographic features in a layered transversely isotropic half-space

The scattering of plane SH-waves by topographic features in a layered transversely isotropic (TI) half-space is investigated by using an indirect boundary element method (IBEM). Firstly, the anti-plane dynamic stiffness matrix of the layered TI half-space is established and the free fields are solved by using the direct stiffness method. Then, Green’s functions are derived for uniformly distributed loads acting on an inclined line in a layered TI half-space and the scattered fields are constructed with the deduced Green’s functions. Finally, the free fields are added to the scattered ones to obtain the global dynamic responses. The method is verified by comparing results with the published isotropic ones. Both the steady-state and transient dynamic responses are evaluated and discussed. Numerical results in the frequency domain show that surface motions for the TI media can be significantly different from those for the isotropic case, which are strongly dependent on the anisotropy property, incident angle and incident frequency. Results in the time domain show that the material anisotropy has important effects on the maximum duration and maximum amplitudes of the time histories.     

Lateral response of dams in semi-elliptical rigid canyons

An analytical closed-form solution is developed for the lateral response of earth and rockfill dams built in semi-elliptical canyons. The dam is idealized as a linearly hysteretic elastic body deforming only in shear, whereas the canyon is assumed to be rigid. The solution for the dam response is given in terms of prolate spheroidal radial and angular functions of the first kind and zero order. Results are presented for natural frequencies, modal displacement shapes, participation factors, and response to transient and steady-state harmonic base excitation for various dam length-to-height ratios. Comparisons are made of the effects of the length-to-height ratio and the canyon shape on the response of dams built in semi-elliptical and rectangular canyons. A subsequent study (Dakoulas, P. & Hsu, C.H., Response of earth dams in semi-elliptical flexible canyons to oblique SH waves, Report, Rice University, Houston, Texas, 1993) extends this model to a semi-elliptical canyon consisting of flexible elastic rock, subjected to obliquely incident harmonic SH waves.