线性黏弹性球面波的特征线分析

基于ZWT黏弹性本构方程建立了体现高应变率效应的黏弹性球面波的控制方程组,包含5个偏微分方程,解5个未知量v、σr、σθ、εr和εθ。采用特征线法,问题转化为解3族特征线上的5个常微分方程,物理上图像清晰,数学上易于求解。特征线数值分析显示,黏弹性球面波的衰减和弥散效应超过线弹性球面波。球面扩散引起的环向拉应力是导致介质拉伸破坏的主因。进一步还针对强间断黏弹性球面波得出其衰减特性的解析表达式,表明这种更强的衰减特性是几何扩散效应和本构黏性效应两者共同作用的后果。     

线黏弹性球面发散应力波的频率响应特性

基于线黏弹性球面波Laplace域的理论解, 得到了不同传播距离处粒子速度、粒子位移、应力、应变等力学量的传递函数。以标准线性固体模型为例, 重点讨论了粒子速度频率响应函数的传播特征, 指出随着传播距离的增加, 粒子速度幅频响应函数的高频响应会低于低频响应, 而在理想弹性条件下, 粒子速度幅频响应函数的高频响应一直高于低频响应。以弹性半径为0.025 m的空腔爆炸为例, 采用Laplace数值逆变换方法对粒子速度波形的演化进行了分析, 给出了粒子速度强间断幅值及粒子速度峰值随传播距离变化的衰减规律曲线, 指出黏弹性介质中粒子速度幅值的衰减曲线介于理想弹性介质中粒子速度幅值衰减曲线和黏弹性介质中粒子速度强间断幅值衰减曲线之间。     

饱和黏弹性多孔介质中的平面波及能量耗散

研究了流体饱和不可压黏弹性多孔介质中的非均匀平面波及其能量流和能量耗散规律. 在流 相和固相物质微观不可压、固相骨架宏观服从积分型本构关系和小变形的假定下，利用 Helmholtz分解，得到了饱和黏弹性多孔介质中非均匀平面波的一般解以及纵波、横波相速 度和衰减率等的解析表达式，分析了平面波传播矢量和衰减矢量之间的关系. 数值结果表明 孔隙流体与固相骨架间的相互作用以及固相骨架的黏性对波的相速度、衰减率等有着显著的 影响. 同时，得到了饱和黏弹性多孔介质的能量方程，给出了能量流矢量和能量耗散率. 对 非均匀平面纵波和横波，推导了平均能量流矢量和平均能量耗散率的解析表达式.     

黏弹性固体中地下爆炸辐射地震波能量的演化

地下爆炸与介质的能量耦合和介质中的波传播机制是理解地下爆炸源物理的重要基础。为研究地下爆炸辐射地震波能量的传播衰减规律,分析了黏弹性介质中地下爆炸地震波能量的组成。基于无限介质中黏弹性球面波理论,给出了速度、位移、应力、应变等物理量Laplace域的理论解。利用Laplace数值逆求解方法,建立了黏弹性介质中地下爆炸辐射地震波场的计算方法。以干黄土作为典型黏弹性材料,计算给出了地震波能量的传播特征,分析了地下爆炸辐射能量的传播衰减规律。结果表明:（1）在黏弹性介质中,某球面处流入的能量随半径增加而逐渐降低。在理想弹性介质中,某球面处流入的能量在几倍弹性半径外即可稳定到某一定值;（2）在某一固定的有限观测区域内,当观测时间足够长时,势能和耗散能均趋于某一定值,辐射动能趋于零;（3）当有限的观测区域能容纳一个完整波长的地震波时,地震波辐射动能的稳态值随波传播距离的增大而减小,总体上可以用指数函数和幂函数进行分段拟合。     

流体饱和标准线性粘弹性多孔介质中的平面波

研究了流体饱和不可压标准线性粘弹性多孔介质中平面波的传播和反射问题．在固相骨架小变形的假定下,得到了粘弹性多孔介质中波动方程的一般解,讨论了弥散关系和波的衰减特性．结果表明：在流体饱和不可压粘弹性多孔介质中,仅存在一个耦合纵波和一个耦合横波,纵波和横波的波速、衰减率等取决于孔隙流体与固相骨架间的相互作用以及固相骨架本身的粘性．同时,研究了半空间自由边界上入射波(纵波、横波)的反射问题。得到了非均匀反射波的波速、反射系数、衰减率等的表达式及其相关的数值结果．     

黏弹性模型微分型本构方程的矩阵形式

针对多个弹性元件和黏性元件以任意连接方式组成的线性黏弹性模型,本文探究了其本构方程的通用矩阵形式表述。首先将研究问题扩展为由Maxwell基本单元构成的标准模型,然后转化为有向图,根据独立路径和闭包围的形式表征出基本应力方程和应变方程,进一步推导得到了任意线性黏弹性模型的微分型本构方程的一般矩阵形式。论文最后总结并建立了一套适合计算机编程的固定范式,利用Python编程实现了该算法、获得了一些数值计算结果。     

黏弹性模型微分型本构方程的矩阵形式

针对多个弹性元件和黏性元件以任意连接方式组成的线性黏弹性模型,本文探究了其本构方程的通用矩阵形式表述。首先将研究问题扩展为由Maxwell基本单元构成的标准模型,然后转化为有向图,根据独立路径和闭包围的形式表征出基本应力方程和应变方程,进一步推导得到了任意线性黏弹性模型的微分型本构方程的一般矩阵形式。论文最后总结并建立了一套适合计算机编程的固定范式,利用Python编程实现了该算法、获得了一些数值计算结果。     

不同检测方法对混凝土波速的影响研究

为建立超声波法和低应变反射波法检测所得混凝土基桩波速的相互转换关系,本文给出了标准线性固体模型下粘弹性介质的三维和一维纵波速度之比及其实用简化公式,并且通过实验得出,混凝土构件的超声波速与1.054倍低应变动测波速之比处于1.03~1.12,具体数值会随着混凝土龄期、弹性模量和强度等级的增加而有所减小。     

基于非局部理论的分数阶黏弹性场地地震放大效应分析

基于非局部理论和分数阶导数理论,研究上覆黏弹性场地土的地震放大效应。利用Eringen非局部理论考虑土体颗粒尺度等非局部效应的影响,通过分数阶黏弹性本构模型刻画场地土的应力应变本构关系,建立基于非局部理论的分数阶黏弹性场地土的振动微分方程;考虑分数阶导数的性质和黏弹性场地土的边界条件,得到了简谐地震波作用下黏弹性场地土的位移和剪切应力的解析解,并在频率域内给出了位移放大系数和应力放大系数的表达式;最后通过数值算例分析了非局部效应、分数阶导数的阶数和土体黏性参数等对黏弹性场地地震放大效应的影响。数值分析结果表明,在低频时位移放大系数和应力放大系数随频率变化曲线存在波动,高频时逐渐趋于稳定;非局部效应对场地土位移放大系数的影响与频率有关,对应力放大系数的影响较大,在研究场地土振动效应时有必要考虑土体非局部效应的影响;分数阶导数的阶数越小,位移放大系数和应力放大系数随频率变化曲线波动越大;场地土的力学性质对场地土的振动效应的影响较大;上覆场地土的黏性对位移放大系数的影响与频率有关,高频时,土体黏性越大,位移放大系数越大;越接近基岩,土体的应力放大系数越大,且土体深度对应力放大系数的影响越大。     

数值仿真轴向运动黏弹性梁非线性参激振动

分别通过两种直接数值方法研究速度变化的经典边界条件下轴向运动黏弹性梁参数振动的稳定性。在控制方程的推导中,采用物质导数黏弹性本构关系和只对时间取偏导数的黏弹性本构关系;分别运用有限差分法和微分求积法对两种经典边界下轴向变速运动黏弹性梁的非线性控制方程求数值解,计算得到梁中点非线性参数振动的稳定稳态响应。数值结果表明,两种黏弹性本构关系对应的稳态响应存在明显差别,同时发现两种直接数值方法的仿真结果基本吻合,证明数值仿真具有较高精度。     

接地网圆钢杆中高频纵向超声导波无损检测方法研究

针对接地网圆钢杆腐蚀检测问题,本文进行了接地网圆钢杆高频纵向超声导波无损检测方法研究。首先对埋地环境下圆钢杆中纵向导波传播特性进行理论研究,分析了不同模态导波的群速度和衰减频散特性。研究发现,高阶纵向模态导波在衰减最小和群速度最大对应的频率处,在圆钢杆中传播能力强,是适合进行地埋圆钢杆导波检测试验的频率范围。在此基础上,进行了埋地圆钢杆高频纵向超声导波无损检测试验研究。结果发现,利用优选的检测参数可以很好地实现埋地圆钢杆中腐蚀缺陷检测。研究工作为接地网运行状态评价提供了很好的技术支撑。     

水泥砂浆的一个热粘弹性率型损伤本构模型

利用SHPB实验系统及自行研制的混凝土类材料快速高温加热设备,对水泥砂浆试件进行了不同 温度(20~600℃)和3种冲击速度下的实验,得到了不同温度和冲击速度下水泥砂浆试件的应力应变关系曲 线。基于ZWT粘弹性本构模型,并且考虑高温下水泥砂浆损伤演化规律都服从Weibull分布,提出了一个水 泥砂浆的热粘弹性率型损伤本构模型。通过数据拟合,获得了本构模型的相关参数,结果表明:理论预测和实 验结果吻合良好。     

Visco-elastic effects on wave dispersion in three-phase acoustic metamaterials

This paper studies the wave attenuation performance of dissipative solid acoustic metamaterials (AMMs) with local resonators possessing subwavelength band gaps. The metamaterial is composed of dense rubber-coated inclusions of a circular shape embedded periodically in a matrix medium. Visco-elastic material losses present in a matrix and/or resonator coating are introduced by either the Kelvin–Voigt or generalized Maxwell models. Numerical solutions are obtained in the frequency domain by means of $\mathbf{k}\left(\omega \right)$-approach combined with the finite element method. Spatially attenuating waves are described by real frequencies ω and complex-valued wave vectors $\mathbf{k}$. Complete 3D band structure diagrams including complex-valued pass bands are evaluated for the undamped linear elastic and several visco-elastic AMM cases. The changes in the band diagrams due to the visco-elasticity are discussed in detail; the comparison between the two visco-elastic models representing artificial (Kelvin–Voigt model) and experimentally characterized (generalized Maxwell model) damping is performed. The interpretation of the results is facilitated by using attenuation and transmission spectra. Two mechanisms of the energy absorption, i.e. due to the resonance of the inclusions and dissipative effects in the materials, are discussed separately.It is found that the visco-elastic damping of the matrix material decreases the attenuation performance of AMMs within band gaps; however, if the matrix material is slightly damped, it can be modeled as linear elastic without the loss of accuracy given the resonator coating is dissipative. This study also demonstrates that visco-elastic losses properly introduced in the resonator coating improve the attenuation bandwidth of AMMs although the attenuation on the resonance peaks is reduced.     

TRANSVERSAL INERTIAL EFFECT ON RELAXATION/RETARDATION TIME OF CEMENT MORTAR UNDER HARMONIC WAVE

Under dynamic loading, the constitutive relation of the cement mortar will be significantly affected by the transversal inertial effect of specimens with large diameters. In this paper, one-dimensional theoretical analysis is carried out to determine the transversal inertial effect on the relaxation/retardation time of the cement mortar under the harmonic wave. Relaxation time or retardation time is obtained by means of the wave velocity, attenuation coefficient and the frequency of the harmonic wave. Thus, the transversal inertial effect on the relaxation time from Maxwell model, as well as on retardation time from Voigt model is analyzed. The results show that the transversal inertial effect may lead to the increase of the relaxation time, but induce the decrease of the retardation time. Those should be taken into account when eliminating the transversal inertial effect in applications.     

Propagation of longitudinal elastic waves in composites with a random set of spherical inclusions (effective field approach)

The work is devoted to the problem of plane monochromatic longitudinal wave propagation through a homogeneous elastic medium with a random set of spherical inclusions. The effective field method and quasicrystalline approximation are used for the calculation of the phase velocity and attenuation factor of the mean (coherent) wave field in the composite. The hypotheses of the method reduce the diffraction problem for many inclusions to a diffraction problem for one inclusion and, finally, allow for the derivation of the dispersion equation for the wave vector of the mean wave field in the composite. This dispersion equation serves for all frequencies of the incident field, properties and volume concentrations of inclusions. The long and short wave asymptotics of the solution of the dispersion equation are found in closed analytical forms. Numerical solutions of this equation are constructed in a wide region of frequencies of the incident field that covers long, middle, and short wave regions of propagating waves. The phase velocities and attenuation factors of the mean wave field are calculated for various elastic properties, density, and volume concentrations of the inclusions. Comparisons of the predictions of the method with some experimental data are presented; possible errors of the method are indicated and discussed.     

地震波倾斜入射下两相饱和土层响应数值模拟

基于Biot两相饱和多孔介质理论建立了频域内饱和土层的刚度矩阵,对从基岩倾斜入射地震波下饱和土层的波动响应进行了数值研究,入射波采取修正的Kanai-Tajimi谱。通过土层地表面处的水平以及竖直放大系数研究了在不同工况条件下土层的波动响应,研究结果表明,土层的波动响应在入射角接近临界入射角时达到峰值;土层的波动响应随着基岩埋深基以及土层阻尼比的增大显著减小;而土层与基岩的波速比对土层的波动响应的影响存在临界波速比,使土层响应达到极值。对于含有夹层的地质环境,研究了波动响应沿深度的分布,通过与均质土层下响应对比,显示非均质夹层能显著减小土层地表处的波动响应,尤其软夹层能较好地阻止入射波向上的传播。     

高应变率大变形下的聚丙烯/尼龙共混高聚物损伤型本构特性

对两种采用不同相容剂的聚丙烯（PP）和尼龙（PA）共混高聚物材料在大变形下的粘弹性力学行为进行研究,着重考察应变率效应和损伤的演化,从而分析不同的界面分子设计对共混体系材料宏观性能的作用。在准静态及冲击实验研究的基础上,基于ZWT非线性粘弹性模型,并结合了遗传算法,分别得到了能有效描述两种共混高聚物大变形阶段计及损伤的非线性粘弹性本构关系。两种材料在不同加载条件下表现出明显不一致的性能,原因在于其损伤演化的率相关性,且两种材料的大变形机制存在一定的差别,能用ZWT方程进行描述的范围也不一样。     

Time–transient response for ultrasonic guided waves propagating in damped cylinders

Herein, an enhanced spectral finite element (SFE) formulation to calculate the time–transient response in cylindrical waveguides is proposed. The original aspect over SFE-based formulations consists in the possibility to account for the effect of material absorption, i.e. guided waves attenuation, on the calculation of the time–transient response.First, the damped steady-state response is constructed by a weighted superposition of the waveguide modal properties obtained from the spectral decomposition of the governing wave equation. To this purpose an enhanced spectrally formulated finite element is developed, in which material damping is included allowing for complex stress–strain viscoelastic constitutive relations in force of the correspondence principle. Dispersive modal properties for the damped waveguide (phase velocity, energy velocity, attenuation and wavestructures) follow straightforwardly by simple formulae. Next, the frequency response of the problem is calculated by weighting the modal data and the spectrum of the applied time-dependent force via Cauchy residue theorem. Finally, the inverse Fourier transform of the frequency response leads to the time–transient response for propagative damped guided waves.The approach is not restricted to any anisotropy degree, holds for any linear viscoelastic constitutive relation that can be characterized and formulated in the frequency domain and it can be applied to SFE formulations for arbitrary cross-section waveguides. A study on guided waves propagating in a scheduled 4.in-40 ANSI steel pipe is presented, where the steel is considered first as perfectly elastic and then as an hysteretic viscoelastic medium, in order to show the effect of material absorption on the time–transient response.