无限弹性土体中洞室在反平面冲击荷载作用下的瞬态响应

为了探讨无限弹性土体内圆柱形洞室在突加反平面冲击荷载作用下的瞬态响应,利用Laplace变换及围道积分逆变换,得到土体位移和应力的一般解析表达式,并给出了数值解。在时域内分析了无限弹性土体内圆柱形孔洞在轴向荷载作用下的动力响应,并将计算结果与采用拉普拉斯数值反变换得到的结果以及静力情况下的结果作了比较。研究结果显示:波到达后,该点土体的应力和位移均瞬间增大,随后慢慢减小,并逐渐趋于静力值;波向外发散传播,并沿半径方向衰减,衰减速度比静力情况的应力衰减慢。     

三层流体中斜入射波作用下半无限板的水弹性响应

解析地研究了在三层流体中斜入射波浪与半无限弹性板的相互作用引起的波散射和板的水弹性响应. 三层流体在界面处的密度发生阶跃, 各层为一常数. 假设流体不可压缩、无黏、流体运动无旋. 在线性势流理论框架下, 使用本征函数展开法和内积式给出波板相互作用的半解析解. 根据色散关系分析, 得到了表面波模态和界面波模态入射时的临界入射角. 随着物理参数的变化, 临界角将随之发生变化. 临界角决定了当由开阔水域向板覆盖水域传播的表面波或界面波的存在性: (1)板覆盖水域入射界面上, 透射波能否存在; (2)入射界面之上界面中, 板覆盖水域中的透射波以及开阔水域中的反射波能否存在. 当下界面波入射时并且入射角足够大时, 开阔水域中的下界面波模态是整个流体域中唯一存在的模态.      

三层流体中斜入射波作用下半无限板的水弹性响应

解析地研究了在三层流体中斜入射波浪与半无限弹性板的相互作用引起的波散射和板的水弹性响应.三层流体在界面处的密度发生阶跃,各层为一常数.假设流体不可压缩、无黏、流体运动无旋.在线性势流理论框架下,使用本征函数展开法和内积式给出波板相互作用的半解析解.根据色散关系分析,得到了表面波模态和界面波模态入射时的临界入射角.随着物理参数的变化,临界角将随之发生变化.临界角决定了当由开阔水域向板覆盖水域传播的表面波或界面波的存在性:(1)板覆盖水域入射界面上,透射波能否存在;(2)入射界面之上界面中,板覆盖水域中的透射波以及开阔水域中的反射波能否存在.当下界面波入射时并且入射角足够大时,开阔水域中的下界面波模态是整个流体域中唯一存在的模态.     

梯度波阻板的地基振动控制研究

波阻板(wave impeding block,WIB)隔振体系是一种有效的振动污染治理措施,虽逐渐被应用在工程实际中,但以往的研究多集中于单相固体均质材料的情形,而对材料特性沿空间连续变化的非均匀固体材料的波阻板隔振性能的研究相对较少.基于功能梯度材料(functionally graded material,FGM)特点,本文提出了以功能梯度波阻板作为隔振屏障的一类新型的地基振动控制体系.考虑在弹性地基内部设置梯度波阻板,基于线弹性理论,利用傅里叶积分变换,根据Helmholtz矢量分解原理,建立了弹性地基在动载荷作用下的回传射线矩阵法(reverberation ray matrix method,RRMM)计算列式.假设梯度波阻板的物理力学性质沿深度方向按幂函数连续变化,采用数值傅里叶逆变换获得了弹性地基的位移和应力等物理量的数值解.通过数值算例,与单相固体均质波阻板进行了对比,并分析讨论了梯度波阻板的材料梯度因子、埋深以及梯度波阻板厚度等物理力学参数对地基隔振性能的影响规律.结果表明,梯度波阻板能有效降低振动的振幅,与单相固体均质波阻板相比,梯度波阻板具有更好的减振隔振效果.地基的位移幅值和应力幅值随着梯度因子的增大而减小.梯度波阻板的隔振效果随着波阻板厚度的增大而提高,而随着梯度波阻板埋深的增大而降低.     

ANALYSIS OF GROUND VIBRATION CONTROL BY GRADED WAVE IMPEDING BLOCK

波阻板（wave impeding block,WIB）隔振体系是一种有效的振动污染治理措施,虽逐渐被应用在工程实际中,但以往的研究多集中于单相固体均质材料的情形,而对材料特性沿空间连续变化的非均匀固体材料的波阻板隔振性能的研究相对较少.基于功能梯度材料（functionally graded material,FGM）特点,本文提出了以功能梯度波阻板作为隔振屏障的一类新型的地基振动控制体系.考虑在弹性地基内部设置梯度波阻板,基于线弹性理论,利用傅里叶积分变换,根据Helmholtz矢量分解原理,建立了弹性地基在动载荷作用下的回传射线矩阵法（reverberation ray matrix method,RRMM）计算列式.假设梯度波阻板的物理力学性质沿深度方向按幂函数连续变化,采用数值傅里叶逆变换获得了弹性地基的位移和应力等物理量的数值解.通过数值算例,与单相固体均质波阻板进行了对比,并分析讨论了梯度波阻板的材料梯度因子、埋深以及梯度波阻板厚度等物理力学参数对地基隔振性能的影响规律.结果表明,梯度波阻板能有效降低振动的振幅,与单相固体均质波阻板相比,梯度波阻板具有更好的减振隔振效果.地基的位移幅值和应力幅值随着梯度因子的增大而减小.梯度波阻板的隔振效果随着波阻板厚度的增大而提高,而随着梯度波阻板埋深的增大而降低.     

非理想界面弹性层/压电柱结构中SH波的传播特性

研究了各向同性弹性层与压电柱之间非理想连接时沿周向传播的SH波的频散特性.弹性层表面力学自由;弹性层与压电柱之间应力连续、位移间断.通过求解控制方程,将问题的解用Bessel函数表示,利用界面条件和边界条件得到频散方程,然后对其进行数值求解,分析了界面性态、材料常数和几何尺寸对SH波传播特性的影响.     

铁木辛柯梁中的卸载弯曲波及二次断裂

半无限长梁承受恒定弯矩作用后, 如果自由端的初始弯矩突然释放, 将在梁中激发出一列卸载弯曲应力波. 采用铁木辛柯梁和瑞利梁来研究突然卸载所激发出的弯曲波的传播特征. 利用拉普拉斯变换方法进行分析, 首先推导出铁木辛柯梁和瑞利梁中的卸载弯曲波的像函数解析解, 采用数值反变换方法给出了时域上波传播的响应解, 并研究了梁中各点的横向位移、弯矩和剪力随时间的变化规律. 计算结果表明: 与简化的欧拉梁不同, 旋转惯性的引入使铁木辛柯梁和瑞利梁中的弯曲波传播具有强烈的局部化效应, 特别是梁中各点经历的弯矩变化, 和其距离自由端的位置相关, 不同时刻的弯矩峰值大小不同;瑞利梁中离自由端不同距离各点的峰值弯矩先增大后降低, 最后达到一个渐近值;铁木辛柯梁中各点的峰值弯矩总体上随着时间单调增大到同一个渐近值, 该渐近值与欧拉梁中的峰值弯矩值相同, 均为1.43.切应力效应的引入进一步降低了铁木辛柯梁中卸载弯曲波的波速, 同时也使得铁木辛柯梁中弯矩峰值的最大值小于瑞利梁中的最大值. 对于脆性细长梁的纯弯曲断裂, 铁木辛柯梁可以较好地预测二次断裂的发生位置, 相应的碎片尺寸约为7倍梁横截面厚度.      

考虑热流固耦合作用的多孔介质孔隙尺度两相流动模拟

为了研究深层油气资源在岩石多孔介质内的运移过程, 使用一种基于Darcy-Brinkman-Biot的流固耦合数值方法, 结合传热模型, 完成了Duhamel-Neumann热弹性应力的计算, 实现了在孔隙模拟多孔介质内的考虑热流固耦合作用的两相流动过程. 模型通过求解Navier-Stokes方程完成对孔隙空间内多相流体的计算, 通过求解Darcy方程完成流体在岩石固体颗粒内的计算, 二者通过以动能方式耦合的形式, 计算出岩石固体颗粒质点的位移, 从而实现了流固耦合计算. 在此基础上, 加入传热模型考虑温度场对两相渗流过程的影响. 温度场通过以产生热弹性应力的形式作用于岩石固体颗粒, 总体上实现热流固耦合过程. 基于数值模型, 模拟油水两相流体在二维多孔介质模型内受热流固耦合作用的流动过程. 研究结果表明: 热应力与流固耦合作用产生的应力方向相反, 使得总应力比单独考虑流固耦合作用下的应力小; 温度的增加使得模型孔隙度增加, 但当注入温差达到150 K后, 孔隙度不再有明显增加; 温度的增加使得水相的相对渗流能力增加, 等渗点左移.      

多孔介质——自由流界面应力跳跃条件下流动特性解析解

对非对称多孔介质--自由流复合通道内多孔介质内部及多孔介质与自由流体界面处复杂质量、动量输运特性进行研究. 在多孔介质区采用Brinkman-extended Darcy模型并结合速度连续,剪切应力跳跃的界面条件对此复合通道内流体的传递现象进行求解,提出了考虑界面应力跳跃时非对称复合通道各区域流体运动速度及摩擦系数的解析式,分析了界面应力跳跃系数,达西数及无量纲多孔层偏心厚度对流体速度及摩擦系数的影响. 结果表明:改变界面性质可在一定条件下明显控制各区域流体速度分布;在达西数、多孔层偏心厚度一定情况下,界面应力系数的增大会使界面流速减小,而使流体摩擦系数增大,特别是界面应力系数小于0的情况下变化更明显,此时若不考虑界面应力系数则会造成较大误差. 当界面应力系数及多孔层偏心厚度均为较小负数值时,改变多孔层偏心厚度对界面速度的影响要大于改变界面应力系数的情况;而当界面应力系数及多孔层偏心厚度为较大正数值时,情况则相反. 较大达西数下,界面应力系数及多孔层偏心厚度对流体摩擦系数的影响均较大,继续减小达西数至一定程度时,界面应力系数对流体摩擦系数的影响可忽略不计而认为只与多孔层偏心厚度相关,且对较大多孔层偏心厚度更敏感.      

流固耦合地震波动问题的显式谱元模拟方法

流固耦合地震波动问题主要研究由流体和固体构成的复杂系统中地震波传播特性及其规律. 传统模拟方法中一般以声波方程、弹性波方程的数值解分别描述理想流体和弹性固体中的波动, 并实时地处理两种不同性质介质之间的相互耦合作用, 数值格式复杂且限制数值模拟精度与计算效率. 本文采用谱元法结合多次透射公式人工边界条件实现了一种流固耦合地震波动问题的高阶显式数值计算方法. 该方法利用了流固耦合问题统一计算框架,可将饱和多孔介质的Biot波动方程分别退化为理想流体的声波方程和弹性固体的弹性波方程. 通过P波垂直入射的水平成层理想流体-饱和多孔介质-弹性固体场地模型、P波斜入射的不规则层状界面以及任意形状界面的理想流体-饱和多孔介质-弹性固体场地模型等三个算例, 与传递函数法解析解以及集中质量有限元法计算结果进行对比分析, 证明了本文方法的正确性与有效性. 数值模拟结果表明, 本文方法相较传统有限元法可以少得多的节点数量获得更高的数值精度, 并且在较宽的频率范围内都能可靠地模拟出流固耦合系统的动力响应, 充分体现出本文方法兼顾高精度、计算效率和复杂场地建模灵活的特点.      

Propagation of waves in an electro-microstretch generalized thermoelastic semi-space

The present investigation is concerned with wave propagation in an electro-microstretch generalized thermoelastic solid half space. Two different cases have been discussed： （i） reflection of plane wave at the free surface of an electro-microstretch generalized thermoelastic solid; and （ii） propagation of Rayleigh waves in an electro-microstretch generalized thermoelastic solid half space. In case （i）, the amplitude ratios of the various reflected waves have been computed numerically and depicted graphically against angle of incidence. In case （ii）, the frequency equation is derived and dispersion curves giving phase velocity and attenuation coefficient as a function of wave number, have been plot- ted graphically for a specific model. Some special cases of interest are also deduced, for both the cases.     

高速荷载下多孔饱和地基的动力响应

研究高速荷载作用下梁与多孔饱和半空间的动力响应。由Fourier变换求解多孔饱和固体的动力基本方程，根据梁与半空间的接触条件得出多孔饱和半空间上梁的垂直位移的表达式。文中的数值算例考虑了荷载移动速度对梁的动力位移的影响，并与相应的弹性半空间问题作了对比。从算例中可以发现荷载移动速度对动力位移有很大的影响，当移动速度与半空间的表面波速相近时，地面会当产生很大的振动，同时还发现当速度大于介质的剪切波速时，多孔饱和半空间上梁的动力响应与弹性半空间上梁的动力响应有很大的差别。     

Effects of initial stress on the reflection and transmission waves at the interface between two piezoelectric half spaces

The effects of initial stress on the reflection and transmission waves at the interface between two piezoelectric half spaces are studied in this paper. First, the secular equations in the traverse isotropic piezoelectric half space are derived from the general dynamic equation with initial stress taken into consideration. Then, the interface conditions that displacement, stress, electric potential, and electric displacement are continuous across interface are required to be satisfied by three sets of coupled waves, namely, quasi-longitudinal wave, quasi-transverse wave and the electric–acoustic wave. The algebraic equations resulting from the interface conditions are solved to obtain the amplitude ratio of various waves and furthermore the energy reflection and transmission coefficients of various waves. The numerical results are shown graphically and the effects of initial stress are discussed.     

Scattering of a pulsed Rayleigh wave by a spherical cavity in an elastic half space

The time-dependent scattering by a spherical cavity in an elastic half space is considered. The incoming wave is a pulsed Rayleigh wave. The stationary part of the problem is solved by the T-matrix method, and an integration in frequency is performed with a modified gaussian weight function. The displacement components at some points on the surface of the half space are computed and shown in a number of plots.     

川藏公路地质环境与整治改建方案的思考

川藏公路由于地质环境复杂、建设标准低、后遗病害多,抗灾能力差,泥石流、滑坡、山崩、雪害、水毁等自然灾害频繁发生,公路阻车断道严重。国家投入巨资进行整治改建,并取得了明显的效果,但由于自然环境特殊、影响因素复杂,许多特大型、大型工程地质病害问题还没有可行、可靠的解决方案。本文通过分析川藏公路沿线的地质环境和灾害特点,总结历年整治改建和经验的教训,提出川藏公路建设的途径、可能达到的目标和应采用的原则。     

初应力对压电层状结构声表面波传播性能的影响

研究了压电层状结构中初应力对广义Ｒａｙｌｅｉｇｈ波传播相速度和机电耦合性能的影响,通过求解含初应力的运动微分方程,对自由界面电学开路和短路两种情况得到了相应的相速度方程。给出了具体的数值算例,所得结果对于提高和改善声表面波器件性能有参考意义。     

Squirt-Flow in Fluid-Saturated Porous Media: Propagation of Rayleigh Waves

Propagation of attenuated waves is studied in a squirt-flow model of porous solid permeated by two different pore regimes saturated with same viscous fluid. Presence of soft compliant microcracks embedded in the grains of stiff porous rock defines the double-porosity formation. Microcracks and pores respond differently to the compressional effect of a propagating wave, which induces the squirt-flow from microcracks to pores. Elastodynamics of constituent particles in porous aggregate is represented through a single-porosity formulation, which involves the frequency-dependent complex moduli. This formulation is deduced as a special case of double-porosity formation allowing the wave-induced flow of pore-fluid. This squirt-flow model of porous solid supports the attenuated propagation of two compressional waves and one shear wave. Superposition of these body waves, subject to stress-free surface, defines the propagation of Rayleigh wave. This wave is governed by a complex irrational dispersion equation, which is solved numerically after rationalising into an algebraic equation. For existence of Rayleigh wave, a complex solution of the dispersion equation should represent a leaky wave, which decays for propagation along any direction in the semi-infinite medium. A numerical example is solved to analyse the effects of squirt-flow on phase velocity, attenuation and polarisation of the Rayleigh waves, for different combinations of parameters. Numerical results suggest the existence of an additional (second) Rayleigh wave in the squirt-flow model of dissipative porous solids.     

Leaky and non-leaky behaviours of guided waves in CF/EP sandwich structures

The aim of this study is to investigate the leaky and non-leaky behaviours of guided waves, between the composite skin and the core in CF/EP sandwich structures, focusing on the fundamental symmetric like and anti-symmetric like guided wave modes and Rayleigh waves. In investigating the core effect on the guided wave propagation different types of cores are used, namely Nomex honeycomb (HRH 10 1/8-3) 10 and 20 mm in thickness and foam (Divinycell^®  PVC). The behaviour of the guided wave modes is characterised and the conversion mechanism to the Rayleigh wave is investigated. Further, leaky and non-leaky behaviours of guided waves upon interacting with debonded areas are explored, where the ability of guided waves to identify debonding of different sizes was assessed. Finite element analysis simulations are presented to support the experimental analysis, where propagation of ultrasonic waves and their interaction with debonded areas are quantitatively examined.     

The transient elastodynamic field excited by trans-Rayleigh trains

The elastic wave field due to a surface load in motion over an elastic half-space is investigated. The model serves as a canonical solution for the modelling of high speed ‘trans-Rayleigh’ trains. The analysis presented leads to closed form expressions for the particle displacement, conical waves and Rayleigh waves as separate contributions. The linearized elastodynamic equations are mapped into a proper form in order to apply the Cagniard-de Hoop technique and find closed form time domain solutions for the particle displacement in the subsonic state, transonic state and supersonic state. A special transformation is used that yields closed form space-time domain expressions for the Conical wave as well as the Rayleigh wave contributions. Attention is focussed on surface source speeds in the neighbourhood of the Rayleigh wave speed and speeds that exceed the wave speed of the shear wave. Numerical results for the conical wave field and Rayleigh wave field are presented at observation points just below the surface showing the enormous effects of the Rayleigh wave at source speeds in the near vicinity of the Rayleigh wave speed.     

Mechanism of shear attenuation near the interface under combined compression and shear impact loading

We investigate the compressional/shear coupling plastic wave propagation characteristics analytically for ideal elastic–plastic materials in both stress and particle velocity spaces, focusing on the shear wave attenuation near the interface occurring in pressure–shear plate impact experiments. The results show that the shear attenuation is strongly associated with the wave propagation characteristics of the coupling waves. In the stress space, as the shear stress increases, an adjustment of the stress components is observed and the final stress state along the wave path is a combined pure shear- and hydrostatic pressure-state. In the particle velocity space, the wave structures with different loading and maximal transverse particle velocity are obtained. The maximal transverse particle velocity varies with the longitudinal velocity and forms a boundary line. Once the loading transverse velocity exceeds this line, a transverse particle velocity discontinuity occurs at the impact interface. If the bonding strength is sufficiently high, there will be a shear band in the target in the extreme vicinity of the interface.