In this paper, the numerical solutions of the nonlinear Hadamard fractional diffusion-wave model with the initial singularity are investigated. Firstly, the model is transformed into coupled equations by virtue of a symmetric fractional-order reduction method. Then the Llog,2-1σ formula on nonuniform grids is applied to approach to the time fractional derivative. In addition, the discrete fractional Grönwall inequality is used to analyze the optimal convergence of the constructed numerical scheme by the energy method. The accuracy of the theoretical analysis will be demonstrated by means of a numerical experiment at the end.