New results for isotropic point scatterers: Foldy revisited

Since 1945, Foldy’s method has been used to predict velocity and attenuation for various types of scatterers. In this paper, it is shown that Foldy’s method also yields predictions of reflection and transmission of scalar waves by a random distribution of point or line scatterers contained in a slab. Results are given in two and three dimensions, and for normal and oblique incidences. Formulae are also obtained for the reflection and transmission of longitudinal waves by point or line scatterers distributed in an elastic (non-viscous) fluid. Energy equations are derived, and expressions are obtained for the energy dissipated in the slab on average over one period. Curves for the reflection and transmission coefficients are presented in the case of solid cylindrical bars immersed in a fluid. The results obtained in this paper are expected to be valid for a low density of scatterers. Potential applications of this work occur in ultrasonic evaluation of materials, seismic exploration and medical ultrasonics, where reflected (or backscattered) data are used to construct maps or images of the materials (metals, composites, earth subsurface, tissue). The formulae of this work are expected to provide useful tools for better and more efficient mapping or imaging.     

考虑密封的悬臂转子系统的裂纹故障分析

以双盘悬臂立式转子-轴承系统为研究对象，建立了系统运动微分方程，并用数值方法分析了在非线性密封力和非线性油膜力作用下的裂纹转子的动力学特性。分析表明，在一定深度裂纹下，转子系统响应随不同角频率比表现出复杂的非线性现象，出现了周期k运动、拟周期运动和混沌运动等多种运动形式。在一定角速度时，工作在远离临界角速度区的转子系统对裂纹非常敏感，而工作在近临界角速度区的转子系统对裂纹不是特别敏感，但是裂纹对它的运动状态影响较大。该研究结果为该类转子-轴承系统的安全运行与故障诊断提供了一定的理论参考。     

P- and S-wave velocity and anisotropy in saturated rocks with aligned cracks

Rocks in the Earth’s crust contain cracks in variable amounts and on variable scales. Depending on the in situ stress level, cracks are always distributed in a common normal direction (so-called aligned cracks). In this study, we focus on the velocity and anisotropy of perfectly aligned cracks (parallel cracks). The theoretical model is introduced to calculate the elastic properties of cracked rocks based on solid matrix properties, crack parameters, and fluid properties. The theoretical results using high- and low-frequency limits and intermediate frequency for velocities and anisotropy are compared with experimental data measured in the laboratory, where rocks are saturated with water and air. We used four synthetic rocks with different crack densities (0.0, 0.0243, 0.0486, and 0.0729) in this study. The comparison shows that the results for high, low, and intermediate frequencies provide good predictions of the shear velocities and anisotropy for saturated cracked rocks. Quantitatively, the results of the P-wave velocity and anisotropy for intermediate frequency fit the measured data better than that for the high- and low-frequency limits. Compared with water saturation, dry cracks with high compliance have a strong influence on the total compliance of rock. The water saturation results in high P-velocities and a low P-anisotropy, while dry cracks induce much lower P-velocities and a high P-anisotropy. The shear wave anisotropy is fluid-independent, theoretical results fit the measured data well in both water saturation and the dry case.     

Transient dynamic responses of a cracked solid subjected to in-plane loadings

Transient dynamic responses of an elastic cracked solid subjected to in-plane surface loadings are investigated in this study. Two vertical cracks, a surface-breaking crack and a sub-surface crack, are considered. The frequency responses of the plane strain problem are calculated by the computational mechanics combining the finite element method with the boundary integral equation. The finite element method is used for the near-field enclosing the crack, while the boundary integral equation is applied for the far-field to satisfy the Sommerfeld radiation condition. The transient responses are then obtained using fast Fourier transform. Surface displacements, crack opening displacements, and dynamic stress intensity factors are presented to show the significant effects of the cracks. The interaction between the elastic waves and the cracks as well as the mode conversion phenomena can be observed and understood through the numerical simulations.     

含裂纹梁参数的谐响应识别

主要研究含裂纹梁在简谐激励作用下的动力特性,提出一种依据幅值变化对裂纹参数进行识别的新方法。首先,在振动过程中考虑裂纹的呼吸特性,以悬臂梁为例建立含裂纹梁的二维有限元模型,指出在一般激励频率下,其对应的幅值均是明显信号,可用来描述裂纹梁的动力特性。其次,当激励频率分别取无裂纹梁一阶频率的1/4和二阶频率的1/4时,幅值随裂纹参数的变化明显不同,可依据响应幅值的变化对裂纹参数进行识别。然后,利用曲面拟合技术绘出幅值变化曲面,对未知参数的裂纹进行识别,验证了该方法的有效性,并归纳出利用幅值变化对任意裂纹参数进行识别的基本步骤。最后,针对无裂纹梁频率计算可能存在误差的情况,分析识别方法的鲁棒性,结果显示即使最大误差为10%,该方法也能对裂纹参数进行准确识别。     

Analysis of plane waves in anisotropic thermoelastic diffusive medium

This paper concentrates on the study of the propagation of harmonic plane waves in a homogeneous anisotropic thermoelastic diffusive medium in the context of different theories of thermoelastic diffusion. It is found that five types of waves propagate in an anisotropic thermoelastic diffusive medium, namely a quasi-elastodiffusive (QED-mode), two quasi-transverse (QSH-mode and QSV-mode), a quasi-mass diffusive (QMD-mode) and a quasi-thermo diffusive (QTD-mode) wave. The governing equations for homogeneous transversely isotropic diffusive medium in different theories of thermoelastic diffusion are taken as a special case. It is noticed that when plane waves propagate in one of the planes of transversely isotropic thermoelastic diffusive solid, purely quasitransverse wave mode(QSH) decouples from rest of the motion and is not affected by the thermal and diffusion vibrations. On the other hand, when plane waves propagate along the axis of solid, two quasi-transverse wave modes (QSH and QSV) decouple from the rest of the motion and are not affected by the thermal and diffusion vibrations. From the obtained results, the different characteristics of waves like phase velocity, attenuation coefficient, specific loss and penetration depth are computed numerically and presented graphically for a single crystal of magnesium. The effects of diffusion and relaxation times on phase velocity, attenuation coefficient, specific loss and penetration depth has been studied. Some particular cases are also discussed.     

Diffraction of antiplane shear waves by two cracks in the presence of a magnetic field

The problem of scattering of horizontal polarized shear waves by the two cracks in a uniform magnetostatic field is considered. The magnetic field is assumed to be parallel to the crack surfaces as well as perpendicular to the crack surfaces. The elastic medium under consideration is a homogeneous, isotropic and infinitely conducting one. The solution of the problem is reduced into a pair of triple integral equations having trigonometrical kernels. Using the finite Hilbert transform technique, solution of the pair of triple integral equations is obtained for the low frequencies. Finally, approximate formulae are derived for the stress intensity factors.     

Dynamic crack propagation in a magneto-electro-elastic solid subjected to mixed loads: Transient Mode-III problem

The transient response of a Mode-III crack propagating in a magneto-electro-elastic solid subjected to mixed loads is investigated through solving the corresponding boundary-initial-value problem in both the cracked solid region and the interior fluid region with treatment of electro-magnetically permeable and impermeable crack face conditions in a unified way. The closed-form results for the dynamic field intensity factors are used to evaluate the dynamic energy release rate through the crack-tip dynamic contour integral. The permeability of the interior fluid region relative to the cracked solid region significantly affects the magneto-electro-mechanical coupling coefficient in the Bleustein–Gulyaev wave function and, consequently, the horizontal shear surface wave speed, the dynamic field intensity factors and the dynamic energy release rate. It is revealed from dynamic fracture mechanics analysis that the dynamic energy release rate thus obtained has an odd dependence on the dynamic electric displacement intensity factor and the dynamic magnetic induction intensity factor. It is also found that the horizontal shear surface wave speed provides the limiting velocity for the propagation of a Mode-III crack in a magneto-electro-elastic solid when there is only applied traction loading.     

Why continuum damage is nonlocal: Justification by quasiperiodic microcrack array

Strain-softening damage due to distributed cracking is modeled by an elastic continuum with a quasiperiodic array of cracks of regular spacing but varying sizes. As a model for the initial stage, the cracks are penny-shaped and small compared to their spacing, and as a model for the terminal stage the uncracked ligaments between the cracks are circular and small compared to their spacing. The strain due to cracks and the compliance per crack are calculated. The cracked material is homogenized in such a manner that the macroscopic continuum strains satisfy exactly the condition of compatibility with the actual strains due to cracks, and the macroscopic continuum stress satisfies exactly the condition of work equivalence with the actual stresses in the cracked material. The results show that, contrary to the existing theories, the damage variable used in continuum damage mechanics should be nonlocal, while the elastic part of the response should be local. In particular, the nonlocal continuum damage should be considered as a function of the spatial average of the cracking strain rather than its local value. The size of the averaging region is determined by the crack spacing.     

A novel method for slant crack detection in rotors based on turning in two directions

This paper presents a novel way to detect fatigue slant cracks in rotors based on theoretical discussion. Hence, the dynamic behaviour of a Jeffcott rotor system with a mid-span slant crack under arbitrary crack orientations is studied. First, using concepts of fracture mechanics, the flexibility matrix and subsequently the system’s stiffness matrix are calculated. A symmetric relation for a global stiffness matrix is presented and proved. Next, the motion equations of the system that are obtained in four directions, two transverse, one torsional and one longitudinal, are solved using the Runge–Kutta numerical method. The characteristics of crack orientations for angles greater than 90° (transverse crack) are investigated in detail and their influence on the elements of the crack compliance matrix is presented. Also, slant crack characteristics with complementary angles are compared to each other. It is shown that the difference between cracked systems with complementary angles is only in 3rd row (3rd column) of the crack compliance matrix, and also it is shown that due to the presence of a slant crack, the system responses in forward and backward motion are different. Using the frequency responses of the shaft obtained, a technique to detect the existence of slant cracks on the shaft was proposed. This novel method is a simple way that can be used for slant crack detection in rotors.     

The dynamic behavior of two parallel symmetric cracks in functionally graded piezoelectric/piezomagnetic materials

The dynamic behavior of two parallel symmetric cracks in functionally graded piezoelectric/piezomagnetic materials subjected to harmonic antiplane shear waves is investigated using the Schmidt method. The present problem can be solved using the Fourier transform and the technique of dual integral equations, in which the unknown variables are jumps of displacements across the crack surfaces, not dislocation density functions. To solve the dual integral equations, the jumps of displacements across the crack surfaces are directly expanded as a series of Jacobi polynomials. Finally, the relations among the electric, magnetic flux, and dynamic stress fields near crack tips can be obtained. Numerical examples are provided to show the effect of the functionally graded parameter, the distance between the two parallel cracks, and the circular frequency of the incident waves upon the stress, electric displacement, and magnetic flux intensity factors at crack tips.     

Dynamic bulk and shear moduli due to grain-scale local fluid flow in fluid-saturated cracked poroelastic rocks: Theoretical model

Grain-scale local fluid flow is an important loss mechanism for attenuating waves in cracked fluid-saturated poroelastic rocks. In this study, a dynamic elastic modulus model is developed to quantify local flow effect on wave attenuation and velocity dispersion in porous isotropic rocks. The Eshelby transform technique, inclusion-based effective medium model (the Mori–Tanaka scheme), fluid dynamics and mass conservation principle are combined to analyze pore-fluid pressure relaxation and its influences on overall elastic properties. The derivation gives fully analytic, frequency-dependent effective bulk and shear moduli of a fluid-saturated porous rock. It is shown that the derived bulk and shear moduli rigorously satisfy the Biot-Gassmann relationship of poroelasticity in the low-frequency limit, while they are consistent with isolated-pore effective medium theory in the high-frequency limit. In particular, a simplified model is proposed to quantify the squirt-flow dispersion for frequencies lower than stiff-pore relaxation frequency. The main advantage of the proposed model over previous models is its ability to predict the dispersion due to squirt flow between pores and cracks with distributed aspect ratio instead of flow in a simply conceptual double-porosity structure. Independent input parameters include pore aspect ratio distribution, fluid bulk modulus and viscosity, and bulk and shear moduli of the solid grain. Physical assumptions made in this model include (1) pores are inter-connected and (2) crack thickness is smaller than the viscous skin depth. This study is restricted to linear elastic, well-consolidated granular rocks.     

Dispersion in oceanic crust during earthquake preparation

Dispersion of Stoneley waves is studied in a sedimentary layer of ocean bottom resting over basaltic solid half space. Sedimentary layer is assumed a transversely isotropic poroelastic medium. Lower-most solid half-space is assumed to be embedded with vertically aligned saturated micro-cracks and behaves transversely isotropic to wave propagation.Frequency equation is obtained in the form of determinantal equation. Role of phase angle is eliminated by expressing slowness of waves in terms of phase velocity and elastic constants. Numerical solutions for phase velocity and group velocity are obtained for a particular model. Calculations are made for different depths of ocean and sediments. Effect of thickness and density of cracks on these velocities are observed.Special cases are discussed which represent the absence of ocean and sediments, in the model considered. Changes in dispersion are discussed during the stress accumulation in an earthquake preparation region.     

Mean field attenuation and amplitude attenuation due to wave scattering

We point out the inadequacy of two widely used approaches of formulating the amplitude attenuation of seismic waves, the formulation of mean-field attenuation and that of scattering coefficient under the single scattering approximation. Using a one-dimensional layered slab, we show that the attenuation of the mean field is merely a statistical effect caused by phase interference among different realizations of the random wave ensemble, and does not represent the amplitude attenuation. We will call the attenuation coefficient of the mean field as the randomization coefficient in order to distinguish it from the amplitude attenuation coefficient.We also show that the scattering coefficient method leads to the same result as the first order approximation to the renormalized perturbation series (or the bilocal approximation to Dyson's equation) of the mean field. Therefore these two approaches are equivalent to a certain degree. This is also shown by using the one-dimensional layered slab model as an example.After pointing out the incorrectness of comparing experiments on amplitude attenuation with the mean field formulation, we suggest and discuss some methods of obtaining the mean field in experiments. For one-dimensional problems, the samples must be taken along the whole propagation path in order to use a spatial average to substitute for ensemble average. For a three-dimensional wave field, measurements over a large seismic array can be used to obtain the mean field. The data from Lasa measured by Aki are used to compare with theory; the agreement between them is good. Finally we compare the mean-field attenuation (randomization) and the amplitude attenuation using the back-halfspace-integration approximation introduced by Wu, and compare them with the measured data by Aki. The comparison shows further the inability of the mean-field formulation in dealing with the problem of amplitude attenuation.     

Damage analysis of thermopiezoelectric properties: Part I - crack tip singularities

This is Part I of the work on a two-dimensional analysis of thermal and electric fields of a thermopiezoelectric solid damaged by cracks. It deals with finding the singular crack tip behavior for the temperature, heat flow, displacements, electric potential, stresses and electric displacements. By application of Fourier transformations and the extended Stroh formalism, the problem is reduced to a pair of dual integral equations for the temperature field with the aid of an auxiliary function. The electroelastic field is governed by another pair of dual integral equations. The inverse square root singularity is found for the heat flow field while the logarithmic singularity prevailed for the electroelastic field regardless of whether the crack lies in a homogeneous piezoelectric solid or at an interface of two dissimilar piezoelectric materials. Results are given for the energy release rate and a finite length crack oriented at an arbitrarily angle with reference to the external disturbances. Part II of this paper considers the modelling of a piezoelectric material containing microcracks. A representative cracked area element is used to obtain the effective conductivity and electroelastic modulus. Numerical results are given for a peizoelectric Bati O₃ ceramic with cracks.     

裂隙挤喷流对孔隙介质排水体积模量的影响

实际岩石比如沉积形成的岩石往往是裂隙和孔隙并存的孔隙介质. 由于扁状的裂隙与近似球形或圆管形的孔隙具有不同的可压缩性,当孔隙介质受压时,液体会从易压缩的裂隙中挤出流入不易压缩的孔隙中,这种挤喷流会引起弹性模量的频散和能量的耗散. 着重研究了裂隙挤喷流和液体可压缩性对孔道变形的影响,推导出了动载荷作用下排水体积模量的表达式. 与挤喷流相关的裂隙附加柔度会引起排水体积模量随频率变化,使得孔隙介质呈现黏弹性. 频率越高,模量的实部越大,岩石抵抗变形的能力越强. 而模量的虚部体现了挤喷流对能量的耗散. 裂隙密度主要决定模量频散的幅度以及能量耗散的强度,且裂隙密度越大,模量频散幅度越大,能量耗散也越强. 裂隙的纵横比主要决定模量频散速率最快或能量耗散最强时对应的特征频率. 若孔隙介质中不含有裂隙,即裂隙密度是0时,排水体积模量退化为Biot理论中的排水体积模量.     

Numerical Modelling of Steady-State Flow in 2D Cracked Anisotropic Porous Media by Singular Integral Equations Method

The equations governing plane steady-state flow in heterogeneous porous media containing curved-line intersecting cracks (Pouya and Ghabezloo in Transp Porous Media 84:511?C532, 2010) and the potential solution obtained for these equations are considered here. The theoretical results are first completed for the mass balance at crack intersections points. Then, a numerical procedure based on a singular integral equations method is described concretely to derive this solution for cracked materials. Closed-form expressions of elementary integrals for special choice of collocation points lead to a very quick and easy numerical method. It is shown that this method can be applied efficiently to the study of the steady-state flow in cracked materials with anisotropic matrix permeability and a dense distribution of curved-line intersecting cracks. Some applications of this method to the permeability of cracked materials are given.     

混凝土面板堆石坝面板缝隙渗流计算模型研究

针对混凝土面板堆石坝面板接缝与面板裂缝的特点,基于等宽缝隙稳定流的运动规律,结合混凝土面板堆石坝的渗流机理,研究建立了面板接缝及面板裂缝的渗流计算模型,并用实例进行了验证。结果表明将等宽缝隙稳定流运动规律应用到面板接缝止水失效且面板产生裂缝情况下面板堆石坝的渗流分析中,可以获得较为准确的结果。且此方法计算过程简单,计算效率高,具有很大的实用性。     

Turbulent mass transport and attenuation in Stokes waves

The motion of turbulent Stokes waves on a finite constant depth fluid with a rough bed is considered. First and second order turbulent boundary layer equations are solved numerically for a range of roughness parameters, and from the solutions are calculated the mass transport velocity profiles and attenuation coefficients. A new mechanism of turbulent mass transport is found which predicts a reduction and reversal of drift velocity in shallow water in agreement with experimental observations under turbulent conditions. This transpires because the second order Stokes wave motion, in a turbulent boundary layer, can directly influence the mass transport velocity by mode coupling interactions between different second order Fourier modes of oscillation. It is also found that the Euler contribution due to the radiation stress of the first order motion is reduced to half of it's corresponding laminar value as a consequence of the velocity squared stress law. The attenuation is found to be of inverse algebraic type with the reciprocal wave height varying linearly with either distance or time. The severe wave height restriction applicable to the Longuet-Higgins [4] solution is shown not to apply to progressive waves on a finite constant depth of fluid. The existence of sand bars on sloping beaches exposed to turbulent waves is predicted.     

Combined analysis of fracture under stresses acting along cracks

A joint approach to the study of two non-classical fracture mechanisms, namely fracture of cracked materials with initial (residual) stresses acting along the crack planes and fracture under compression along parallel cracks, is considered in the framework of three-dimensional linearized solid mechanics. Mathematical statements of problems for pre-stressed solids that contain interacting circular cracks are given. Problems for an infinite solid containing two parallel co-axial cracks and for a space with the periodical set of co-axial parallel cracks as well as for a half-space with near-the-surface crack are solved. Several patterns of loading on the crack faces (normal loading, radial shear and torsion) are considered. The effects of initial stresses on stress intensity factors are analyzed for highly elastic materials with some types of elastic potentials. Formulation of fracture criteria accounting effect of initial (residual) stresses is given. Critical parameters of fracture of solids containing interacting cracks under compression along the cracks are calculated. The influence of geometrical parameters of the problems as well as physical and mechanical properties of materials on these critical parameters is analyzed.