A spherical wave in a viscoelastic half-space

The propagation of spherical waves in an isotropie elastic medium has been studied sufficiently completely (see, e.g., [1–4]). it is proved [5, 6] that in imperfect solid media, the formation and propagation of waves similar to waves in elastic media are possible. With the use of asymptotic transform inversion methods in [7] a problem of an internal point source in a viscoelastic medium was investigated. The problem of an explosion in rocks in a half-space was considered in [8]. A numerical Laplace transform inversion, proposed by Bellman, is presented in [9] for the study of the action of an explosive pulse on the surface of a spherical cavity in a viscoelastic medium of Voigt type. In the present study we investigate the propagation of a spherical wave formed from the action of a pulsed load on the internal surface of a spherical cavity in a viscoelastic half-space. The potentials of the waves propagating in the medium are constructed in the form of series in special functions. In order to realize viscoelasticity we use a correspondence method [10]. The transform inversion is carried out by means of a representation of the potentials in integral form and subsequent use of asymptotic methods for their calculation. Thus, it becomes possible to investigate the behavior of a medium near the wave fronts. The radial stress is calculated on the surface of the cavity.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 139–146, March–April, 1976.     

The nature of the triple point singularity in the case of stationary reflection of weak shock waves

The structure of flow in the vicinity of a triple point in the problem of stationary irregular reflection of weak shock waves is numerically investigated within the framework of the Euler model, including the von Neumann paradox range. To improve the accuracy of the solution near singular points a new technology including a grid contracted toward the triple point and the discontinuity fitting is applied. It is shown that in the four-wave flow pattern the curvatures of the tangential discontinuity and the Mach front at the triple point are finite. The singularity is concentrated only in a sector between the reflected wave front and the expansion fan. When the three-wave flow pattern is realized, the curvatures of the tangential discontinuity and both wave fronts at the triple point are infinite. On the range of weak and moderate shock waves the logarithmic singularity in subsonic sectors near the triple point conserves up to transition to the regular reflection.     

Longitudinal force on an embedded filament

Summary A thin elastic filament embedded in an elastic medium is subjected to a concentrated longitudinal load. For two-dimensional geometry and a relatively stiff filament the application of the load gives rise to a system of wedge-like and cylindrical waves. The dynamic shear stresses at the interface of the filament and the matrix, and in the region of the cylindrical waves are determined by means of Fourier transform techniques and Cagniard's method [2]. At the wave fronts of the wedge-like waves the jumps in the shear stresses are computed. Along the filament, the magnitudes of propagating discontinuities decrease exponentially. Along rays normal to the wave fronts of the wedge-like waves, the magnitudes of propagating discontinuities remain unchanged.     

同震断层三维裂纹扩展波动时域超奇异积分法

应用波动时域超奇异积分法将P波、S波和磁电热弹多场耦合作用下同震断层任意形状三维裂纹扩展问题转化为求解以广义位移间断率为未知函数的超奇异积分方程组问题;定义了广义应力强度因子,得到裂纹前沿广义奇异应力增量解析表达式;应用波动时域有限部积分概念及体积力法,为超奇异积分方程组建立了数值求解方法,编制了FORTRAN程序,以三维矩形裂纹扩展问题为例,通过典型算例,研究了广义应力强度因子随裂纹位置变化规律;分析了同震断层裂纹扩展中力、磁、电场辐射规律.      

Wave propagation in a non-ideal dusty gas

In this paper we have studied the behavior of wave motion as propagating wavelets and their culmination into shock waves in a non-ideal gas with dust particles. In the absence of non-ideal effect the gas satisfies an equation of state of Mie–Gruneisen type. An expansion wave resulting from the action of receding piston is considered and the solutions to this problem showing effects of dust particles and non-idealness are obtained. The propagation of weak waves is considered and the flow variables in the region bounded by the piston and the characteristic wave front are found out. The expansive action of a receding piston undergoing an abrupt change in velocity is discussed. Cases of central expansion fan and shock fronts are studied and the solutions up to first order in the physical plane are obtained. The effects of non-idealness and dust particles are discussed in each case.     

The growth and decay of acceleration waves of arbitrary form in bkz viscoelastic fluids

The article explores the amplitude behavior of an acceleration wave of arbitrary form propagating into a particular non-linear viscoelastic fluid with memory. The media is assumed to obey the incompressible, isotropic and isothermal BKZ constitutive model. Investigation is restricted to waves propagating into regions which have been at rest in their reference configuration. Specific cases of plane, cylindrical and spherical wave fronts are examined. The results indicate that the acceleration wave amplitude (which is transverse) obeys a similar equation as found by Varley for simple materials, and hence will always decay.     

Stress intensity factor analysis of a three-dimensional interfacial corner between anisotropic bimaterials under thermal stress

A numerical method using a path-independent H-integral based on the conservation integral was developed to analyze the singular stress field of a three-dimensional interfacial corner between anisotropic bimaterials under thermal stress. In the present method, the shape of the corner front is smooth. According to the theory of linear elasticity, asymptotic stress near the tip of a sharp interfacial corner is generally singular as a result of a mismatch of the materials’ elastic constants. The eigenvalues and the eigenfunctions are obtained using the Williams eigenfunction method, which depends on the anisotropic materials’ properties and the geometry of an interfacial corner. The order of the singularity related to the eigenvalue is real, complex or power-logarithmic. The amplitudes of the singular stress terms can be calculated using the H-integral. The stress and displacement around an interfacial corner for the H-integral are obtained using finite element analysis. In this study, a proposed definition of the stress intensity factors of an interfacial corner, which includes those of an interfacial crack and a homogeneous crack, is used to evaluate the singular stress fields. Asymptotic solutions of stress and displacement around an interfacial corner front are uniquely obtained using these stress intensity factors. To prove the accuracy of the present method, several different kinds of examples are shown such as interfacial corners or cracks in three-dimensional structures.     

Singular stress fields near contact boundaries in a composite bolted joint

Singular stresses arising in the neighborhood of contact surfaces introduced in laminated orthotropic plates by mechanical joining with clamp-up were investigated by using local asymptotic solutions and full-field numerical analysis. Three-dimensional B-spline approximation of the displacements and a penalty function-based contact solution was used in the numerical analysis. Recent work has shown that fracture in bolted composite joints may initiate near the outer edge of the bolt head or washer away from the hole edge, particularly if the joint is preloaded. Material and geometric discontinuities exist in these regions, resulting in singular stress behavior. Asymptotic stress analysis was performed to obtain the power of singularity in these regions as a function of the bolt-head (washer) stiffness. Frictionless contact conditions were assumed. It was found that the characteristics of the stress singularity for such practically important combinations as titanium bolt-head and carbon fiber composite plate are similar to a crack in terms of the power of singularity and uniqueness of the singular term. Coefficients of the singular terms of the asymptotic expansion were determined by comparison with the numerical solution in the close vicinity of the singular contour. Good agreement between the asymptotic and numerical solution in the transition regions was observed.     

Dispersion of Small Amplitude Waves in a Pre-Stressed, Compressible Elastic Plate

The dispersion of harmonic waves, propagating along a principal direction in a pre-stressed, compressible elastic plate, is investigated in respect of the most general isotropic strain-energy function. Different cases, dependent on the choice of material parameters and pre-stress, are analysed. A complete long and short wave asymptotic analysis is carried out, with the approximations obtained giving phase speed (and frequency) as explicit functions of wave and mode number. Various wave fronts, both associated with the short wave limit of harmonics and arising through the combination of harmonics in a narrow wave speed region, are discussed. It is mentioned that the case of high compressibility is of particular interest. In contrast with the classical (un-strained) case, the longitudinal body wave speed may be less than the corresponding shear wave speed. In consequence, the short wave limit of all harmonics may be the appropriate longitudinal wave speed; contrasting with the classical case for which this limit is necessarily associated with a shear wave front. A further possible short wave limit is also shown to exist for which the associated wave normal has a component in the direction normal to the plate. Particularly novel numerical results are presented when the longitudinal and shear wave speeds are equal. The analysis is illustrated by numerical calculations for various strain-energy functions.     

Dynamic crack-tip fields in rate-sensitive solids

The asymptotic stress and strain fields near the tip of a crack which propagates dynamically in a rate-sensitive solid are obtained under anti-plane shear and plane strain conditions. The problem is formulated within the context of a small-strain theory for a solid whose mechanical behavior under high strain rates is described by an elastic-viscoplastic constitutive relation. It is shown that, if the stresses are singular at the crack-tip, the viscoplastic relation is equivalent asymptotically to an elastic-non-linear viscous relation. Furthermore, for a certain range of the material parameter which characterizes the rate-sensitivity of the material, the elastic strain-rates near the propagating crack tip are shown to have the same asymptotic radial dependence near the propagating crack-tip as the inelastic strain-rates. This determines the order of the stress singularity uniquely. The governing equations for anti-plane shear and plane strain are then derived. The numerical results for the stress and strain fields are presented for anti-plane shear and plane strain. For the present model, the results suggest that under small-scale yielding conditions, there exists a minimum velocity for stable steady crack propagation. The implication that a terminal velocity for a running crack may exist is also discussed.     

应力波传播与屈曲耦合情况下结构动力屈曲控制方程探讨

从应力波作用下结构动力屈曲的特点出发,指出应力波作用下结构动力屈曲与结构中应力波传播的耦合导致时间成为结构动力屈曲的参变量,从而应力波作用下结构动力屈曲问题中结构的真实运动与邻近运动是不同时刻、不同扰动区域的比较,这使得其动力屈曲控制方程的建立不宜采用传统的等时积分变分原理;以压应力波作用下弹性直杆为例,应用能量守恒原理,根据屈曲时刻结构能量的转换关系,建立了弹性压应力波作用下半无限长直杆的动力屈曲控制方程,并得到了波前附加约束条件;最后,讨论了波前附加约束条件的物理意义,指出波前附加约束条件出现的根本原因是轴向应力波的传播与屈曲不能解耦。     

Influence of singularity dominated zone for propagating cracks in finite size specimens

The singularity dominated zones for straight as well as curved cracks propagating in finite size specimens were determined experimentally by using the optical method of dynamic photoelasticity using the near-field stress equations. Experimental data was carefully analyzed using improved numerical schemes to get the complete stress field around the propagating crack. This stress field was critically examined to evaluate the size of singularity dominated zones for cracks propagating in straight as well as curved paths. For this purpose, the exact solution was compared with the singular solution using stress components σ_x, σ_y, τ_xy and the maximum shear stress τ_max as a criterion respectively. For straight cracks where the stress field is symmetric about the crack path, the singularity dominated zones can be determined by using any one of the stresses. However, for a curved crack, the zones were unsymmetric. This study shows that σ_y, the crack opening stress, yields the best result for characterizing the singularity dominated zone around a running crack tip.     

Nonlinear analysis of the flow initiated by the sudden motion of a wedge

Certain self-similar problems involving the sudden motion of a wedge which were treated in the linear approximation in [1–3] are studied by the method of matched asymptotic expansions. The nature of the wave boundary of the perturbed region is determined. Second-approximation solutions are constructed which describe flows behind weak shock fronts propagating in a stationary gas and behind fronts of weak discontinuity lines propagating by known uniform flows. A boundary-value problem is formulated whose solution describes, in first approximation, flows in the neighborhoods of points of interaction of the fronts. The existence of similarity rules of flows in these nieghborhoods is estimated. An approximate solution of the problems is given.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 37–47, May–June, 1976.     

Theoretical transient analysis and wave propagation of piezoelectric bi-materials

The transient response of piezoelectric bi-materials subjected to a dynamic anti-plane concentrated force or electric charge with perfectly bonded interface is examined in the present study. The problem is solved by using the Laplace transform method and the inverse Laplace transform is evaluated by means of Cagniard’s method. Exact transient full-field solutions of the contribution for each wave are expressed in explicit closed forms. The transient behavior of field quantities is examined in detail by numerical calculations. The existence condition of a propagating surface wave along the interface is discussed in detail. A surface wave can be guided by the interface of two semi-infinite materials in contact if one, at least, of these two materials is piezoelectric. The propagation velocity of the surface wave is explicitly expressed and is found to be less than the lower shear wave velocity of the two materials. The existence of the surface wave for piezoelectric–piezoelectric bi-materials is restricted to the situation that the shear waves of the two piezoelectric materials are very close. The possibility for the existence of the surface wave for piezoelectric–elastic bi-materials is much greater than that of the piezoelectric–piezoelectric bi-materials.     

Coarsening Fronts

We characterize the spatial spreading of the coarsening process in the Allen–Cahn equation in terms of the propagation of a nonlinear modulated front. Unstable periodic patterns of the Allen–Cahn equation are invaded by a front, propagating in an oscillatory fashion, and leaving behind the homogeneous, stable equilibrium. During one cycle of the oscillatory propagation, two layers of the periodic pattern are annihilated. Galerkin approximations and the Conley index for ill-posed spatial dynamics are used to show existence of modulated fronts for all parameter values. In the limit of small amplitude patterns or large wave speeds, we establish uniqueness and asymptotic stability of the modulated fronts. We show that the minimal speed of propagation can be characterized by a dichotomy which depends on the existence of pulled fronts. The main tools here are an Evans function type construction for the infinite-dimensional ill-posed dynamics and an analysis of the complex dispersion relation based on Sturm–Liouville theory.     

The Green’s function of the mild-slope equation:: The case of a monotonic bed profile

In the present work the Green’s function of the mild-slope and the modified mild-slope equations is studied. An effective numerical Fourier inversion scheme has been developed and applied to the construction and study of the source-generated water-wave potential over an uneven bottom profile with different depths at infinity. In this sense, the present work is a prerequisite to the study of the diffraction of water waves by localized bed irregularities superimposed over an uneven bottom. In the case of a monotonic bed profile, the main characteristics of the far-field are: (i) the formation of a shadow zone with an ever expanding width, which is located along the bottom irregularity, and (ii) in each of the two sectors not including the bottom irregularity the asymptotic behavior of the wave field approaches the form of an outgoing cylindrical wave, propagating with an amplitude of order O(R^−1/2), where R is the distance from the source, and wavelength corresponding to the sector-depth at infinity. Moreover, the weak wave system propagating in the shadow zone is of order O(R^−3/2), and along the bottom irregularity consists of the superposition of two outgoing waves with wavelengths corresponding to the two depths at infinity.     

Gaussian beam formulations and interface conditions for the one-dimensional linear Schrödinger equation

Gaussian beams are asymptotic solutions of linear wave-like equations in the high frequency regime. This paper is concerned with the beam formulations for the Schrödinger equation and the interface conditions while beams pass through a singular point of the potential function. The equations satisfied by Gaussian beams up to the fourth order are given explicitly. When a Gaussian beam arrives at a singular point of the potential, it typically splits into a reflected wave and a transmitted wave. Under suitable conditions, the reflected wave and/or the transmitted wave will maintain a beam profile. We study the interface conditions which specify the relations between the split waves and the incident Gaussian beam. Numerical tests are presented to validate the beam formulations and interface conditions.     

Wave interaction in a nonequilibrium gas flow

By employing the method of multiple time scales, we derive here the transport equations for the primary amplitudes of resonantly interacting high-frequency waves propagating into a non-equilibrium gas flow. Evolutionary behavior of non-resonant wave modes culminating into shocks or no shocks, together with their asymptotic decay behavior, is studied. Effects of non-linearity, which are noticeable over times of order O(ε^-1), are examined, and the model evolution equations for resonantly interacting multi-wave modes are derived.     

Plotnikov——Toland板模型中水弹性孤立波的迎撞

通过奇异摄动方法研究了在薄冰层覆盖的不可压缩理想流体表面上传播的两个水弹性孤立波之间的迎面碰撞.借助特殊的 Cosserat 超弹性壳 理论以及Kirchhoff--Love 板理论,冰层由 Plotnikov--Toland板模型描述.流体运动采用浅水假设和Boussinesq 近似. 应用Poincaré--Lighthill--Kuo 方法进行坐标变形,进而渐近求解控制方程及边界条件, 给出了三阶解的显式表达. 可以观察到碰撞后的孤立波不会改变它们的形状和振幅. 波浪轮廓在碰撞之前是对称的, 而在碰撞之后变成不对称的并且在波传播方向上向后倾斜. 弹性板和流体表面张力减小了波幅. 图示比 较了本文与已有结果可知线性板模型可作为本文的一个特例.      

Diagonalization method for a singularly perturbed vector Robin problem

In this paper the method and technique of the diagonalization are employed to transform a vector second-order nonlinear system into two first-order approximate diagonalized systems. The existence and the asymptotic behavior of the solutions are obtained for a vector second-order nonlinear Robin problem of singular perturbation type.