Numerical calculation of three-dimensional elastic wave equation of piles staggered grid finite difference using method with variable step lengths
LIU Hua-xuan1, 2, LIU Dong-jia1, LU Zhi-tang3, TAO Jun1, JIANG Jing1
1. School of Resources and Environmental Engineering, Hefei University of Technology, Hefei 230009, China; 2. Anhui Institute of Building Research and Design, Hefei 230001, China; 3. Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China
Abstract:To improve the computational efficiency of low-strain numerical simulation of piles, the variable step staggered grid finite difference method is introduced into the calculation of three-dimensional pile-soil model based on the traditional difference method. By differentiating three-dimensional elastic wave equation and introducing absorbing boundary conditions on the computational model interface, the low-strain numerical simulation of the response of piles is calculated under transient vertical excitating force. The results show that the computational accuracy is guaranteed and the computational process is more flexible and simple because of varied difference grids. The process can also reduce the memory requirements and the computational time. The absorbing boundary which absorbs the incident waves reaching the computational model interface is effective. The proposed method has achieved good effects on fitting low-strain measured curve and numerical simulation of asymmetric defective piles.
刘华瑄, 刘东甲, 卢志堂, 陶俊, 姜静. 桩基三维弹性波动方程变步长交错网格有限差分数值计算[J]. 岩土工程学报, 2014, 36(9): 1754-1760.
LIU Hua-xuan, LIU Dong-jia, LU Zhi-tang, TAO Jun, JIANG Jing. Numerical calculation of three-dimensional elastic wave equation of piles staggered grid finite difference using method with variable step lengths. Chinese J. Geot. Eng., 2014, 36(9): 1754-1760.
[1] 陈凡, 王仁军. 尺寸效应对基桩低应变完整性检测的影响[J]. 岩土工程学报, 1998, 20(5): 92-96. (CHEN Fan, WANG Ren-jun. Dimension effect on low strain integrity testing of piles[J]. Chinese Journal of Geotechnical Engineering, 1998, 20(5): 92-96. (in Chinese)) [2] CHOW Y K, PHOON K K, CHOW W F, et al. Low strain integrity testing of piles: three dimensional effects[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2003, 129(11): 1057-1062. [3] CHAI H Y, PHOON K K, ZHANG D J. Effects of the source on wave propagation in pile integrity testing[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2010, 136(9): 1200-1208. [4] AKI K, RICHARDS P G. Quantitative seismology, theory and methods[M]. San Francisco: W H Freeman & Co, 1980. [5] SCHUBERT F, PEIFFER A, KOHLER B. The elastodynamic finite integration technique for waves in cylindrical geometries[J]. Acoustical Society of America, 1998, 104(5): 2064-2614. [6] GSELL D, LEUTENEGGER T, DUAL J. Modeling three-dimensional elastic wave propagation in circular cylindrical structures using a finite-difference approach[J]. Acoustical Society of America, 2003, 116(6): 3284-3293. [7] PITARKA A. 3D elastic finite-difference modeling of seismic motion using staggered grids with nonuniform spacing[J]. Bulletin of the Seismological Society of America, 1999, 89(1): 54-68. [8] 董良国, 马在田, 曹景忠, 等. 一阶弹性波方程交错网格高阶差分解法[J]. 地球物理学报, 2000, 43(3): 411-419. (DONG Liang-guo, MA Zai-tian, CAO Jing-zhong, et al. A staggered-grid high-order difference method of one-order elastic wave equation[J]. Chinese Journal of Geophysics, 2000, 43(3): 411-419. (in Chinese)) [9] 柴华友. 桩土相互作用对基桩定量分析的影响[J]. 岩土力学, 1996, 17(4): 41-47. (CHAI Hua-you. Effect of pile-soil interaction on pile integrity quantitative analysis[J]. Rock and Soil Mechanics, 1996, 17(4): 41-47. (in Chinese)) [10] 刘东甲. 完整桩瞬态纵向振动的模拟计算[J]. 合肥工业大学学报(自然科学版), 2000, 23(5): 683-687. (LIU Dong-jia. Simulation and calculation of transient longitudinal vibration of integrate piles[J]. Journal of Hefei University of Technology (Natural Science), 2000, 23(5): 683-687. (in Chinese)) [11] 柯宅邦, 刘东甲. 低应变反射波法测桩的轴对称问题数值计算[J]. 岩土工程学报, 2006, 28(12): 2111-2115. (KE Zhai-bang, LIU Dong-jia. Numerical computation of axisymmetric problems in low strain integrity tests on piles[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(12): 2111-2115. (in Chinese)) [12] 段 军, 刘东甲, 程晓东, 等. 桩的一维纵向振动问题的交叉网格有限差分法数值计算[J]. 工程地球物理学报, 2008, 5(1): 54-59. (DUAN Jun, LIU Dong-jia, CHENG Xiao-dong, et al. Numerical calculation of one dimensional longitudinal vibration of pile by finite difference method of staggered grid[J]. Chinese Journal of Engineering Geophysics, 2008, 5(1): 54-59. (in Chinese)) [13] 卢志堂, 刘东甲, 龙丽丽, 等. 基桩低应变检测三维问题的数值计算[J]. 合肥工业大学学报(自然科学版), 2011, 34(6): 905-909. (LU Zhi-tang, LIU Dong-jia, LONG Li-li, et al. 3-D numerical computation for low strain integrity testing of piles[J]. Journal of Hefei University of Technology (Natural Science), 2011, 34(6): 905-909. (in Chinese)) [14] LU Z T, WANG Z L, LIU D J. Study on low-strain integrity testing of pipe-pile using the elastodynamic finite integration technique[J]. International Journal for Numerical and Analytical methods in Geomechanics, 2013, 37: 536-550. [15] 孙卫涛. 弹性波动方程的有限差分数值方法[M]. 北京: 清华大学出版社, 2009. (SUN Wei-tao. The finite difference numerical method of elastic wave equation[M]. Beijing: Tsinghua University Press, 2009. (in Chinese)) [16] HIGDON R L. Absorbing boundary conditions for difference approximations to the multi-dimensional wave equation[J]. Mathematics of Computation, 1986, 47: 437-459.