激光激发黏弹表面波有限元数值模拟

研究黏弹性材料中激光激发的Rayleigh波的传播特征. 考虑到黏弹性材料的黏性特征,在频域内建立黏弹性材料中激光激发Rayleigh波的有限元数值模型. 在验证有限元频域数值模型正确性的基础上,模拟脉冲激光作用在黏弹性材料上激发出Rayleigh波,进而讨论激光激发的黏弹Rayleigh波的传播特征,并比较黏弹性材料与弹性材料中激光激发的Rayleigh波差异,同时分析了材料的黏性劲度参量变化对Rayleigh波特征的影响. 关键词： 表面波 激光超声 有限元方法 黏弹性     

Wave dispersion and attenuation in viscoelastic isotropic media containing multiphase flow and its application

In this paper,we introduce the complex modulus to express the viscoelasticity of a medium.According to the correspondence principle,the Biot-Squirt(BISQ)equations in the steady-state case are presented for the space-frequency domain described by solid displacements and fluid pressure in a homogeneous viscoelastic medium.The effective bulk modulus of a multiphase flow is computed by the Voigt formula,and the characteristic squirt-flow length is revised for the gas-included case.We then build a viscoelastic BISQ model containing a multiphase flow.Through using this model,wave dispersion and attenuation are studied in a medium with low porosity and low permeability.Furthermore,this model is applied to observed interwell seismic data.Analysis of these data reveals that the viscoelastic parameter tanδ is not a constant.Thus,we present a linear frequency-dependent function in the interwell seismic frequency range to express tanδ.This improves the fit between the observed data and theoretical results.     

Viscoelastic effects on asymmetric two‐dimensional vortex patterns in a strongly coupled dusty plasma

Strongly coupled dusty plasma medium is often described as a viscoelastic fluid that retains its memory. In a flowing dusty plasma medium, vortices of different sizes appear when the flow does not remain laminar. The vortices also merge to transfer energy between different scales. In the present work, we study the effect of viscoelasticity and compressibility over a localized vortex structure and multiple rotational vortices in a strongly coupled viscoelastic dusty plasma medium. In case of single rotating vortex flow, a transverse wave is generated from the localized vortex source and the evolution time of generated waves is found to be reduced due to finite viscoelasticity and compressibility of the medium. It is found that the viscoelasticity suppresses the dispersion of vorticity. In the presence of multiple vortices, we find, the vortex mergers get highly affected in the presence of memory effect of the fluid, and thus the dynamics of the medium gets completely altered compared to a non‐viscoelastic fluid. For a compressible fluid, viscoelasticity dampens the energy in the sonic waves generated in the medium. Thus a highly viscoelastic and compressible fluid, in some cases, behaves similarly to an incompressible viscoelastic fluid. The wave‐front like rings propagate in elliptical orbits keeping the footprint of the earlier position of the point‐vortex. The rings collide with each other even within the patch vortex region forming regions of high vorticity at the point of intersection and pass through each other.     

Viscoelastic BISQ Model for Low-Permeability Sandstone with Clay

A modified BISQ （Blot/Squirt） model for wave propagation in low-permeability sandstone is developed by introducing the viscoelastic mechanism of a porous skeleton into Dvorkin＇s model. The linear viscoelasticity of the Kelvin Voigt constitutive law is employed to describe the stress-strain relation of a solid frame with clay while the ultrasonic waves propagate through the fluid-saturated sandstone. The phase velocity and attenuation of two p-waves are given based on the present BISQ model. The comparisons between numerical results and experimental data indicate that our viscoelastic model is more realistic and feasible for wave propagation in the low-permeability sandstone, especially with clay, than traditional BISQ models.     

黏弹介质包裹的充液管道中低频轴对称波传播特性

研究埋地充液管道中低频轴对称波传播特性。将土壤考虑为黏弹介质,结合Kennard薄壳方程和Kelvin-Voigt线性黏弹性模型,引入土壤载荷矩阵,推导出土-管滑移情形下流体主导波和管壁压缩波的相速度表达式。通过数值模拟计算得到流体主导波和管壁压缩波的频散和衰减曲线并进行可靠性验证,分析两种波引起的管壁径向位移之比,讨论厚径比和品质因子对流体主导波传播的影响。结果表明,黏弹介质对流体主导波和管壁压缩波的相速度影响较小,但对衰减影响较大;流体主导波对管壁径向位移有较大的影响,是泄露噪声传播的主要载体;厚径比越大,流体主导波的相速度越大,衰减越小;而品质因子越大,流体主导波的频散和衰减都越小。研究结果可为埋地充液管道的泄漏检测提供一定的理论参考。      

Viscosity and Diffusion Effects at the Boundary Surface of Viscous Fluid and Thermoelastic Diffusive Solid Medium

This paper concentrates on the wave motion at the interface of viscous compressible fluid half-space and homogeneous isotropic, generalized thermoelastic diffusive half-space. The wave solutions in both the fluid and thermoelastic diffusive half-spaces have been investigated; and the complex dispersion equation of leaky Rayleigh wave motion have been derived. The phase velocity and attenuation coefficient of leaky Rayleigh waves have been computed from the complex dispersion equation by using the Muller's method. The amplitudes of displacements, temperature change and concentration have been obtained. The effects of viscosity and diffusion on phase velocity and attenuation coefficient of leaky Rayleigh waves motion for different theories of thermoelastic diffusion have been depicted graphically. The magnitude of heat and mass diffusion flux vectors for different theories of thermoelastic diffusion have also been computed and represented graphically.     

Ultrasound aided smooth dispensing for high viscoelastic epoxy in microelectronic packaging

Epoxy dispensing is one of the most critical processes in microelectronic packaging. However, due its high viscoelasticity, dispensing of epoxy is extremely difficult, and a lower viscoelasticity epoxy is desired to improve the process. In this paper, a novel method is proposed to achieve a lowered viscoelastic epoxy by using ultrasound. The viscoelasticity and molecular structures of the epoxies were compared and analyzed before and after experimentation. Different factors of the ultrasonic process, including power, processing time and ultrasonic energy, were studied in this study. It is found that elasticity is more sensitive to ultrasonic processing while viscosity is little affected. Further, large power and long processing time can minimize the viscoelasticity to ideal values. Due to the reduced loss modulus and storage modulus after ultrasonic processing, smooth dispensing is demonstrated for the processed epoxy. The subsequently color temperature experiments show that ultrasonic processing will not affect LED’s lighting. It is clear that the ultrasonic processing will have good potential to aide smooth dispensing for high viscoelastic epoxy in electronic industry.     

各向异性粘弹性孔隙介质地震波场伪谱法正演模拟

地球介质各向异性、粘弹性以及孔隙特征是地震资料描述油气储层时应综合考虑的. 建立了各向异性粘弹性孔隙介质模型，导出了各向异性粘弹性孔隙介质的弹性波波动方程，采用伪谱法正演模拟了各向异性粘弹性孔隙介质地震波，进行了其波场特征分析. 该研究有益于加深对地震波在实际地球介质中传播规律的认识.     

Extension of the mode method for viscoelastic media and focused ultrasonic beams

In the present study viscoelasticity is introduced in the mode model and the orthogonality condition is adapted for viscous media. The expansion of convergent acoustic Gaussian beams in terms of radiation modes for viscoelastic media is studied as well. The effects on the reflected and transmitted profiles of acoustic beams incident from an ideal liquid onto a viscoelastic plate are shown and physically explained. It is shown that focusing the incident beam can suppress divergence effects and gives the possibility to measure shear wave attenuation coefficients.     

Viscoelastic wave propagation in the viscoelastic single walled carbon nanotubes based on nonlocal strain gradient theory

In this paper, the viscoelastic wave propagation in an embedded viscoelastic single-walled carbon nanotube (SWCNT) is studied based on the nonlocal strain gradient theory. The characteristic equation for the viscoelastic wave in SWCNTs is derived. The emphasis is placed on the influence of the tube diameter on the viscoelastic wave dispersion. A blocking diameter is observed, above which the wave could not propagate in SWCNTs. The results show that the blocking diameter is greatly dependent on the damping coefficient, the nonlocal and the strain gradient length scale parameters, as well as the Winkler modulus of the surrounding elastic medium. These findings may provide a prospective application of SWCNTs in nanodevices and nanocomposites.     

Dispersion curves for a viscoelastic Timoshenko beam with fractional derivatives

The Kramers–Kronig dispersion relation, often used as a viscoelastic constitutive law for polymeric materials, is based on purely mathematical properties of linearity, convergence of improper integrals, and causality; thus, it may also be valid as a viscoelastic constitutive law for general structural materials. Accordingly, the motion equation of a Timoshenko beam composed of conventional elastic structural materials is extended to one composed of viscoelastic materials. From the derived governing equation, a dispersive equation is derived for a viscoelastic Timoshenko beam. By plotting phase velocity curves and group velocity curves for a beam of solid circular cross-section composed of a viscoelastic material (polyvinyl chloride foam), the influence of the fractional order of viscoelasticity is examined. As a result, it is found that, in the high frequency range, only the first mode of a Timoshenko beam converged to the propagation velocity of the Rayleigh wave, which takes account of the fractional order of viscoelasticity. In addition, the phase velocity and the group velocity were found to increase as the fractional order approaches 0, and to decrease as the fractional order approaches 1. Furthermore, the rate of velocity change becomes greater as the fractional order approaches 0, and becomes smaller as the fractional order approaches 1.     

平面波与线性粘弹海底的反射和折射

采用基于广义标准线性体Boltzmann叠加原理的固体流变模型来描述粘弹性海底底质结构,分析了沉积层和基岩的纵、横波速度频散和衰减的特点,计算了海水沉积层界面和沉积层基岩界面位移势反射和折射系数和能量反射和折射系数。结果表明,采取本模型的粘弹性海底的速度频散、衰减和已有文献的实验数据符合较好,能准确的描述海底底质的主要性质,反射和折射系数和声波的频率有关。      

Modeling of magnetoelastic shear waves due to point source in a viscoelastic crustal layer over an inhomogeneous viscoelastic half space

In the present study, propagation of magnetoelastic shear wave due to a momentary point source in a viscoelastic crustal layer over inhomogeneous viscoelastic half space has been discussed. Green’s function technique and Fourier transform along with method of successive approximation are used to find the closed-form solutions for displacement and generalized shear wave period equation. Attenuation of the resultant shear wave is computed and effects of magnetic field, width of the layer, complex wave number, viscosity, and inhomogeneity parameters are distinctly marked on dissipation curves using two-dimensional and surface plots. It is found that effect of layer’s magnetoelastic coupling parameter on attenuation pattern of shear wave is just the reverse of half space magnetoelastic coupling parameter. Similarly, internal friction of layer has somewhat different effect on shear wave angular frequency than lower half space viscosity. Certain published results are also derived as special cases to the present study.     

热粘弹波在变温非均匀合金熔体中的传播

机械波在金属凝固过程中传播的定量计算一直是一个难题,主要原因就是在这个过程中的熔体结构非常复杂.本研究考虑到熔体的变温、非均匀和粘弹性的特点,采用Kelvin粘弹性介质模型,建立了具有粘热损失特性的热粘弹性波动方程,通过隐式有限差分方法对波动方程进行求解,并以ZL203A合金熔体为研究对象,探究了热粘弹波在变温非均匀介质中的传播规律.结果表明:热粘弹波从合金熔体的低温区向高温区传播时,非均匀的温度场对波的传播有较大影响;相反,当波从合金熔体的高温区向低温区传播时,非均匀的温度场对波的传播几乎没有影响.热粘弹波在合金熔体中的衰减系数随频率的增大呈线性增大,而随温度的升高先增大后减小,在熔体的枝晶搭接温度附近达到最大值.     

Shape stability of a gas bubble in a soft solid

Predicting the onset of non-spherical oscillations of bubbles in soft matter is a fundamental cavitation problem with implications to sonoprocessing, polymeric materials synthesis, and biomedical ultrasound applications. The shape stability of a bubble in a Kelvin-Voigt viscoelastic medium with nonlinear elasticity, the simplest constitutive model for soft solids, is analytically investigated and compared to experiments. Using perturbation methods, we develop a model reducing the equations of motion to two sets of evolution equations: a Rayleigh-Plesset-type equation for the mean (volume-equivalent) bubble radius and an equation for the non-spherical mode amplitudes. Parametric instability is predicted by examining the natural frequency and the Mathieu equation for the non-spherical modes, which are obtained from our model. Our theoretical results show good agreement with published experiments of the shape oscillations of a bubble in a gelatin gel. We further examine the impact of viscoelasticity on the time evolution of non-spherical mode amplitudes. In particular, we find that viscosity increases the damping rate, thus suppressing the shape instability, while shear modulus increases the natural frequency, which changes the unstable mode. We also explain the contributions of rotational and irrotational fields to the viscoelastic stresses in the surroundings and at the bubble surface, as these contributions affect the damping rate and the unstable mode. Our analysis on the role of viscoelasticity is potentially useful to measure viscoelastic properties of soft materials by experimentally observing the shape oscillations of a bubble.     

The nonlinear propagation of acoustic waves in a viscoelastic medium containing cylindrical micropores

Based on an equivalent medium approach, this paper presents a model describing the nonlinear propagation of acoustic waves in a viscoelastic medium containing cylindrical micropores. The influences of pores' nonlinear oscillations on sound attenuation, sound dispersion and an equivalent acoustic nonlinearity parameter are discussed. The calculated results show that the attenuation increases with an increasing volume fraction of micropores. The peak of sound velocity and attenuation occurs at the resonant frequency of the micropores while the peak of the equivalent acoustic nonlinearity parameter occurs at the half of the resonant frequency of the micropores. Furthermore, multiple scattering has been taken into account, which leads to a modification to the effective wave number in the equivalent medium approach. We find that these linear and nonlinear acoustic parameters need to be corrected when the volume fraction of micropores is larger than 0.1%.     

Rayleigh wave propagation along the boundary of a non-Newtonian fluid-saturated porous medium

The Frenkel-Biot theory is used to study the reflection of elastic waves from the boundary of a non-Newtonian (Maxwell) fluid-saturated porous medium. The velocity and attenuation of a Rayleigh surface wave propagating along the boundary of the medium are determined. Two models of a fluid-saturated porous medium are used for calculation: with pore channels of a fixed diameter and with a lognormal distribution of pore channels in size. The results of calculations show that, when the fluid in the porous medium is characterized by a small Deborah number (i.e., exhibits non-Newtonian properties), the velocity of Rayleigh waves exhibits a considerable frequency dispersion. The results also suggest that, in principle, it is possible to estimate the Deborah number from the measured frequency dispersion of the Rayleigh wave velocity.     

Attenuation and dispersion of Rayleigh waves propagating on a cracked surface: an effective field approach

A new effective field approach describing the attenuation and dispersion of a Rayleigh wave propagating on a surface containing a distribution of one-dimensional, surface-breaking cracks is presented. Limited by the validity of the independent scattering approximation, the model utilizes the complex transmission coefficient of a single crack to build expressions for the phase velocity and attenuation coefficient of an effective Rayleigh wave. The model is shown to be able to accommodate the effect of compressive residual stresses that tend to close the distributed cracks at their mouth, and therefore substantially reduce the attenuation and the velocity change caused by an equivalent distribution of open cracks. Wherever possible, the predictions of the new model are compared with those of other approaches and critical remarks are offered which discuss the advantages of the effective field approach over the others. Finally, an extension of this approach to distributions of two-dimensional surface-breaking cracks is outlined.     

Phase velocities and attenuations of shear, Lamb, and Rayleigh waves in plate-like tissues submerged in a fluid (L)

In the past several decades, the fields of ultrasound and magnetic resonance elastography have shown promising results in noninvasive estimates of mechanical properties of soft tissues. These techniques often rely on measuring shear wave velocity due to an external or internal source of force and relating the velocity to viscoelasticity of the tissue. The mathematical relationship between the measured velocity and material properties of the myocardial wall, arteries, and other organs with non-negligible boundary conditions is often complicated and computationally expensive. A simple relationship between the Lamb-Rayleigh dispersion and the shear wave dispersion is derived for both the velocity and attenuation. The relationship shows that the shear wave velocity is around 20% higher than the Lamb-Rayleigh velocity and that the shear wave attenuation is about 20% lower than the Lamb-Rayleigh attenuation. Results of numerical simulations in the frequency range 0-500 Hz are presented.     

Guided waves propagating in sandwich structures made of anisotropic,viscoelastic, composite materials

The propagation of Lamb-like waves in sandwich plates made of anisotropic and viscoelastic material layers is studied. A semi-analytical model is described and used for predicting the dispersion curves (phase velocity, energy velocity, and complex wave-number) and the through-thickness distribution fields (displacement, stress, and energy flow). Guided modes propagating along a test-sandwich plate are shown to be quite different than classical Lamb modes, because this structure does not have the mirror symmetry, contrary to most of composite material plates. Moreover, the viscoelastic material properties imply complex roots of the dispersion equation to be found that lead to connections between some of the dispersion curves, meaning that some of the modes get coupled together. Gradual variation from zero to nominal values of the imaginary parts of the viscoelastic moduli shows that the mode coupling depends on the level of material viscoelasticity, except for one particular case where this phenomenon exists whether the medium is viscoelastic or not. The model is used to quantify the sensitivity of both the dispersion curves and the through-thickness mode shapes to the level of material viscoelasticity, and to physically explain the mode-coupling phenomenon. Finite element software is also used to confirm results obtained for the purely elastic structure. Finally, experiments are made using ultrasonic, air-coupled transducers for generating and detecting guided modes in the test-sandwich structure. The mode-coupling phenomenon is then confirmed, and the potential of the air-coupled system for developing single-sided, contactless, NDT applications of such structures is discussed.