In this paper, we study the numerical solution of initial boundary-value problem for the fourth-order partial integro-differential equations with a weakly singular kernel. We use the forward Euler scheme for time discretization and the quasi-wavelet based numerical method for space discretization. Detailed discrete formulations are given to the treatment of three different boundary conditions, including clamped-type condition, simply supported-type condition and a transversely supported-type condition. Some numerical experiments are included to demonstrate the validity and applicability of the discrete technique. The comparisons of present results with analytical solutions show that the quasi-wavelet based numerical method has a distinctive local property. Especially, the method is easy to implement and produce very accurate results.