Abstract: In the acoustic convolution quadrature time-domain boundary element method, the calculation of convolution quadrature coefficients is computationally intensive and may not converge. To address this issue, this paper presents a division method on the basis of the direct expansion method to compute these coefficients. This method replaces the overall calculation in the direct expansion method which leads to the amplification of coefficient errors, with a distributed transfer calculation that reduces errors in multiple steps and overcomes the divergence problem. Numerical experiments on convolution coefficients for element pairs with different distances show that this method overcomes the problem of divergence in the calculation of the direct expansion method, and improves the accuracy of the acoustic convolution quadrature time-domain boundary element method. Meanwhile, the calculation time of this method is on the same order of magnitude as the direct expansion method which is much smaller than the traditional method. This promotes the application of the direct expansion method in time-domain acoustic model simulation calculations with larger sizes or smaller time steps.