SPH方法在模拟线弹性波传播中的运用  

Application of smoothed particle hydrodynamics method to the simulations of elastic wave propagation in solid

作　　者：丁桦[1] 龙丽平[2] 伍彦峰[1] 

机构地区：[1]中国科学院力学研究所工程科学部,北京100080 [2]北京航空航天大学航空科学与工程学院,北京100083

出　　处：《计算力学学报》2005年第3期320-325,共6页Chinese Journal of Computational Mechanics

基　　金：国家973(2002CB412706);中国科学院知识创新项目滑坡工程中的关键力学问题资助项目.

摘　　要：通过对固体中波动问题的模拟建立了一种光滑粒子法的新形式,一种运用SPH的核函数的类似有限体积法的计算方法。通过对统计体积的修正以及对边界粒子的核函数修正,较好地解决了SPH方法中长期以来制约其被广泛应用的主要问题之一边界条件的表述。在此基础上成功地在光滑粒子法中实现了透射边界条件的模拟。同时利用反卷积修正使得较大粒子间距下的计算结果的精度大大提高。这种方法不但保持了SPH的简单性,而且很容易实现应力边界条件。A new version of SPH method is developed to solve the initial-boundary value problems for elastic wave propagation. In this algorithm, a very weak kernel estimate is used. In particular, the stress boundary conditions (especially the transmitting boundary condition) are invoked through residual terms in the integration by parts. The very weak kernel estimate avoids the double kernel estimate in the conventional residual method of SPH. The estimates of boundary particles are done by modified kernel function defined according to the local feature of a boundary. The final solution is sharpened by de-convolution in order to obtain reasonable results with relatively large particle distance. The conservative smoothing approach is used to stabilizing the numerical procedure.

关 键 词：SPH 弹性波传播 透射边界条件 反卷积 守恒滤波 

分 类 号：O340[理学—固体力学]