P波作用下饱和土中深埋圆形复合式衬砌隧道动应力响应研究

采用波函数展开法,推导出了平面P波入射下饱和土中深埋圆形双层衬砌动应力集中问题的解析解,并将该解退化为弹性介质中双层衬砌对平面P波散射的稳态解,验证了计算结果的正确性。研究结果表明,饱和土中双层衬砌的内边界动应力集中系数值大于外边界;外侧衬砌及内侧衬砌弹性模量对衬砌的动应力集中系数影响较大,在保证外侧衬砌稳定性较好的条件下,内侧衬砌尽可能选择刚度较小的材料;增加外侧衬砌及内侧衬砌的厚度可明显地减小衬砌内的动应力集中系数。     

Dynamic stress intensity factors of two 3D rectangular cracks in a transversely isotropic elastic material under a time-harmonic elastic P-wave

The dynamic stress intensity factors (DSIFs) of two 3D rectangular cracks in a transversely isotropic elastic material under an incident harmonic stress wave are investigated by generalized Almansi’s theorem and the Schmidt method in the present paper. Using 2D Fourier transform and defining the jumps of displacement components across the crack surface as the unknown functions, three pairs of dual integral equations are derived. To solve the dual integral equations, the jumps of the displacement components across the crack surfaces are expanded in a series of Jacobi polynomials. Numerical examples are provided to show the effects of the geometric shape of the rectangular crack, the characteristics of the harmonic wave and the distance between two rectangular cracks on the DSIFs of the transversely isotropic elastic material.     

Dynamic investigation of a functionally graded layered structure with a crack crossing the interface

The dynamic response of a functionally graded layered structure with a crack crossing the interface is analyzed. The in-plane impact loading condition is considered. By using the Laplace and Fourier integral transforms, singular integral equation method and residue theory, the present problem is reduced to a singular integral equation in the Laplace transform domain. The influences of Young’s modulus ratio, thickness ratio, and crack length and location on the dynamic stress intensity factors (DSIFs) are investigated. Particularly, the DSIFs corresponding to different crack locations are shown in the case when the crack center moves from one layer to another layer through the interface. The peak and static values and overshoot characteristics of the DSIFs are analyzed. It is found that these values typically exhibit kinking behavior when the crack tips arrive at the interface. This study is different from previous other investigations in the following respects: (1) the dynamic response of a crack crossing the interface of a functionally graded structure is studied analytically, which has hardly been done in the past and (2) the present model can be reduced to some important problems, such as a functionally graded coating-substrate structure with a crack in the graded coating or homogeneous substrate or one intersecting the interface.     

SH波作用下地表软覆盖层中圆形夹杂的动应力分析

利用复变函数法和波函数展开法, 对地表软覆盖层中浅埋圆形夹杂在稳态SH波作用下的动应力集中问题进行研究并给出了解析解。根据SH波散射时的衰减特性, 采用了大圆弧假定的方法, 将半空间覆盖层直线边界问题转化为曲面边界问题。通过算例分析了SH波垂直入射时, 不同入射波波数和圆夹杂与半空间的波数比对圆形夹杂周边动应力集中因子的分布和动应力集中因子最大值变化的影响。算例表明, 圆形夹杂越“软”, 其波数越大, 夹杂周边的动应力集中因子越大; 入射波波数约0.35时, 夹杂周边的最大动应力集中因子达到最大值。     

半空间双相压电介质垂直边界附近圆孔对SH波的散射

利用“Green函数法”和“镜像法”对垂直边界附近含圆孔的半空间双相压电介质对SH波的散射问题进行分析,得到其稳态解。利用镜像法得到满足水平边界应力自由与电位移自由的波函数解析表达式。根据垂直边界连续性条件,利用“契合法”建立第一类Fredholm型积分方程组,得到圆孔周边的动应力集中系数与电场强度集中系数解析表达式。数值算例分析了入射波频率、入射角度、介质参数等对动应力集中系数与电场强度集中系数的影响,并与已有文献进行比较。计算表明,高频SH波垂直入射危害较大。     

动载下裂纹应力强度因子计算的改进型扩展有限元法

相较于常规扩展有限元法(extended finite element method, XFEM), 改进型扩展有限元法(improved XFEM) 解决了现有方法线性相关与总体刚度矩阵高度病态问题, 在数量级上提升了总体方程的求解效率, 克服了现有方法在动力学问题中的能量正确传递、动态应力强度因子数值震荡、精度低下问题. 本文基于改进型XFEM, 采用Newmark 隐式时间积分算法, 重点研究了动载荷作用下扩展裂纹尖端应力强度因子的求解方法, 与静力学方法相比, 增加了裂纹扩展速度项与惯性项的贡献. 通过数值算例研究了网格单元尺寸、质量矩阵、时间步长、裂尖加强区域、惯性项、扩展速度项及相互作用积分区域J-domain的网格与单元尺寸对动态应力强度因子求解精度的影响, 验证了改进型XFEM计算动态裂纹应力强度因子方法的有效性. 针对文献中具有挑战性的 "I 型半无限长裂纹先稳定后扩展"问题, 改进型XFEM给出目前为止精度最好的动态应力强度因子数值解.      

弹性波绕任意形状界面孔的散射

求解了弹性波绕任意形状界面孔的散射问题．通过入射波、反射波或折射波及孔的散射波场的叠加,得到了界面孔在SH波绕射下的总波场．总波场波函数的级数项待定系数可采用边界配点法来确定,该法不受边界正交性的限制,能够适用于任意形状的边界．最后,对界面椭圆孔进行了实例计算,得到了椭圆孔边的动应力集中系数．     

Dynamic interaction between a sub-interface crack and the interface in a bi-material: time-domain BEM analysis

Summary Transient response of a sub-interface crack in a bi-material is studied with emphasis on the dynamic interaction between the crack and the interface, by combining the traditional time-domain displacement boundary element method (BEM) and the non-hypersingular traction BEM. Computations are performed for an unbounded bi-material with a crack subjected to impact tensile loading on its faces or incident impact waves and a bounded rectangular bi-material plate under remote impact tensile loading. Numerical results of the dynamic stress intensity factors (DSIFs) and dynamic interface tractions are presented for various material combinations and crack locations. It is shown that pronounced increases in DSIFs and the interface tractions may be caused in some cases because of the dynamic interaction between the crack and the interface.This work was initialized during the second author's stay at Institute of Mechanics, TU Darmstadt, Germany under the support of the Alexander von Humboldt Foundation. Discussion on the BEM formulation with Dr. Seelig is gratefully acknowledged. The first two authors are also grateful for the partial support by the China National Natural Science Foundation under Grant No. 10025211 and the NJTU Scientific Paper Fund (PD195).     

Hopkinson拉杆的入射波波形整形实验研究

目前,分离式Hopkinson杆实验技术已经被广泛用于测试材料在10~2~10~4s~(-1)应变率范围内的动态力学特性。为了抑制入射波的高频振荡,实现恒定应变率加载,本文利用分离式Hopkinson拉杆(SHTB)实验装置,研究了加载金属短杆(2A12T4铝合金)及整形垫片(纸板、PVC软塑料及带磁性胶皮)对入射波波形的影响。实验结果表明,整形垫片降低了入射应力脉冲的高频振荡,获得了比较平滑的入射应力脉冲,延长了上升时间。同时,利用所得的波形整形结果,对2A12T4铝合金进行了拉伸应力波脉冲加载的拉伸和断裂实验测试。     

Symmetric-Galerkin boundary element analysis of the dynamic T-stress for the interaction of a crack with an auxetic inclusion

In this work, the symmetric-Galerkin boundary element method (SGBEM) for 2-D elastodynamics in the Laplace-space frequency domain (Laplace domain) is employed to study the dynamic stress intensity factors (DSIFs) and the dynamic T-stress (DTS) during the interaction between a crack and an auxetic inclusion under impact loading conditions. It is found that, while the auxeticity has virtually no effect on the DSIFs, its influence on the DTS is noticeable. This finding is particularly important as it implies the imperative need of fracture criteria based on both the DSIFs and DTS for predicting crack growth in composite materials with auxetic phases.     

Effects of electric/magnetic impact on the transient fracture of interface crack in piezoelectric-piezomagnetic sandwich structure: anti-plane case

Due to the incompatibility of the interlaminar deformations,the interface debonding or cracking usually happens in a layered magnetoelectric(ME)structure under an applied load.In this paper,the transient responses of the anti-plane interface cracks in piezoelectric(PE)-piezomagnetic(PM)sandwich structures are studied by the standard methods of the integral transform and singular integral equation.Discussion on the numerical examples indicates that the PE-PM-PE structure under electric impact is more likely to fracture than the PM-PE-PM structure under a magnetic impact.The dynamic stress intensity factors(DSIFs)are more sensitive to the variation of the active layer thickness.The effects of the material constants on the DSIFs are dependent on the roles played by PE and PM media during the deformation process.     

直立码头前船波浪力耦合计算模型

建立了外域用Boussinesq方程、内域用刚体运动方程的直立码头前二维船剖面波浪力的时 域计算耦合模型,内域与外域在交界面的匹配条件是流量连续和压力相等. 进行了相关模型 实验,并把计算结果与实验结果进行了对比. 推导了船体与水底和直立码头之间间隙内流体 运动的自振频率,研究了间隙内流体运动的共振现象.     

A fractional differential constitutive model for dynamic stress intensity factors of an anti-plane crack in viscoelastic materials

Fractional differential constitutive relationships are introduced to depict the history of dynamic stress inten- sity factors （DSIFs） for a semi-infinite crack in infinite viscoelastic material subjected to anti-plane shear impact load. The basic equations which govern the anti-plane deformation behavior are converted to a fractional wave-like equation. By utilizing Laplace and Fourier integral transforms, the fractional wave-like equation is cast into an ordinary differential equation （ODE）. The unknown function in the solution of ODE is obtained by applying Fourier transform directly to the boundary conditions of fractional wave-like equation in Laplace domain instead of solving dual integral equations. Analytical solutions of DSIFs in Laplace domain are derived by Wiener-Hopf technique and the numerical solutions of DSIFs in time domain are obtained by Talbot algorithm. The effects of four parameters α, β, b1, b2 of the fractional dif- ferential constitutive model on DSIFs are discussed. The numerical results show that the present fractional differential constitutive model can well describe the behavior of DSIFs of anti-plane fracture in viscoelastic materials, and the model is also compatible with solutions of DSIFs of anti-plane fracture in elastic materials.     

Axisymmetric longitudinal wave propagation in a finite prestretched compound circular cylinder made of incompressible materials

The propagation of an axisymmetric longitudinal wave in a finite prestrained compound (composite) cylinder is investigated using a piecewise-homogeneous body model and the three-dimensional linearized theory of wave propagation in prestressed body [13–15]. The inner and outer cylinders are assumed to be made of incompressible neo-Hookean materials. Numerical results on the influence of the prestrains in the inner and outer cylinders on wave dispersion are presented and discussed. These results are obtained for the case where the inner solid cylinder is stiffer than the outer hollow cylinder. In particular, it is established that the pretension of the cylinders increases the wave velocity     

辛解析奇异元在动荷载作用下V型切口问题中的应用

利用辛解析奇异单元,结合时域精细算法,研究了动荷载作用下的含平面V型切口问题。时域上,采用时域精细算法,并结合自适应算法控制展开项数,保证了计算精度。空间域上,切口尖端附近采用辛解析奇异单元,其余区域采用常规有限单元,避免了局部网格加密。本文使用的辛解析奇异单元不需要过渡单元和局部网格加密,且能够通过奇异单元内部的参数关系直接给出切口尖端的应力强度因子,不需要复杂的后处理过程。数值结果表明,本文方法具有良好的精度和稳定性,可以准确地计算动态应力强度因子。     

径向非均匀压电介质中圆孔对SH波的散射

利用复变函数理论对SH波作用下含圆孔径向非均匀压电介质的反平面动力特性进行了研究。压电介质的密度沿径向按幂函数形式变化,但压电参数、弹性参数、介电参数均为常数。利用变量替换法将非均匀压电介质的变系数波动方程组转化为标准的Helmholtz方程组,得到了满足边界条件的波函数解析表达式。通过数值算例分析了入射角度、入射波频率、幂次等对应力集中系数和电场强度集中系数的影响,并与已有文献进行比较。结果表明,某些参数组合下,动应力集中系数与电场强度集中系数均随幂次增大而增大。     

Scattering by penetrable acoustic targets

An acoustic target of constant density ?_t and variable index of refraction is imbedded in a surrounding acoustic fluid of constant density ?_a. A time harmonic wave propagating in the surrounding fluid is incident on the target. We consider two limiting cases of the target where the parameter ε ≡ ?_a/?_t → 0 (the nearly rigid target) or ε → ∞ (the nearly soft target). Wh en the frequency of the incident wave is bounded away from the ‘in-vacuo’ resonant frequencies of the target, the resulting scattered field is essentially the field scattered by the rigid target for ε = 0 or the soft target if ε → ∞. However, when the frequency of the incident wave is near a resonant frequency,the target oscillates and its interaction with the surrounding fluid produces peaks in the scattered field amplitude. In this paper we obtain asymptotic expansions of the solutions of the scattering problems for the nearly rigid and the nearly soft targets as ε → 0 or ε → ∞, respectively, that are uniformly valid in the incident frequency. The method of matched asymptotic expansions is used in the analysis. The outer and inner expansions correspond to the incident frequencies being far or near to the resonant frequencies, respectively. We have applied the method only to simple resonant frequencies, but it can be extended to multiple resonant frequencies. The method is applied to the incidence of a plane wave on a nearly rigid sphere of constant index of refraction. The far field expressions for the scattered fields, including the total scattering cross-sections, that are obtained from the asymptotic method and from the partial wave expansion of the solution are in close agreement for sufficiently small values of ε.     

Transient response of an interface crack between dissimilar piezoelectric layers under mechanical impacts

Using the integral transform and the Cauchy singular integral equation methods, the problem of an interface crack between two dissimilar piezoelectric layers under mechanical impacts is investigated under the permeable electrical boundary condition on the crack surface. The dynamic stress intensity factors (DSIFs) of both mode-I and II are determined. The effects of the crack configuration and the combinations of the constitutive parameters of the piezoelectric materials on the dynamic response are examined. The numerical calculation of the mode-I plane problem indicates that the DSIFs may be retarded or accelerated by specifying different combinations of material parameters. In addition, the parameters of the crack configuration, including the ratio of the crack length to the layer width and the ratio between the widths of two layers, exert a considerable influence on the DSIFs. The results seem useful for design of the piezoelectric structures and devices of high performance.     

准饱和土体中圆形衬砌对弹性波的散射

采用Vardoulakis和Beskos提出的准饱和土体的波动控制方程,根据Helmholtz矢量分解定理,得到了准饱和土中P1波(快压缩波)、P2波(慢压缩波)和S波(剪切波)的波数的势函数表达式.将准饱和土体和圆形衬砌视为各向同性的均质体,运用波函数展开法将入射波、散射波和折射波的势函数展开成Fourier-Bessel函数的级数形式,根据准饱和土体与衬砌边界处应力和位移连续及衬砌内完全自由的边界条件,得到了平面P1波入射下,准饱和土体内深埋圆形衬砌的散射系数和折射系数的理论解,通过数值计算分析了饱和度对准饱和土体和衬砌的DSCF(动应力集中因子)及准饱和土体的PPCF(孔压集中因子)的影响规律,结果表明:准饱和土体的DSCF随着饱和度的增大而减小,衬砌的DSCF基本不受饱和度的影响,而准饱和土体的PPCF则随着饱和度的增大而增大.     

Dynamic stress intensity factors for homogeneous and smoothly heterogeneous materials using the interaction integral method

Dynamic stress intensity factors (DSIFs) are important fracture parameters in understanding and predicting dynamic fracture behavior of a cracked body. To evaluate DSIFs for both homogeneous and non-homogeneous materials, the interaction integral (conservation integral) originally proposed to evaluate SIFs for a static homogeneous medium is extended to incorporate dynamic effects and material non-homogeneity, and is implemented in conjunction with the finite element method (FEM). The technique is implemented and verified using benchmark problems. Then, various homogeneous and non-homogeneous cracked bodies under dynamic loading are employed to investigate dynamic fracture behavior such as the variation of DSIFs for different material property profiles, the relation between initiation time and the domain size (for integral evaluation), and the contribution of each distinct term in the interaction integral.