In this paper, we study the nonlinear stability of travelling wave solutions with shock profile for a relaxation model with a nonconvex flux, which is proposed by Jin and Xin [<i>Comm. Pure Appl. Math.</i>, 48 (1995), pp. 555--563] to approximate an original hyperbolic system numerically under the subcharacteristic condition introduced by T. P. Liu [<i>Comm. Math. Phys.</i>, 108 (1987), pp. 153--175]. The travelling wave solutions with strong shock profile are shown to be asymptotically stable under small disturbances with integral zero using an elementary but technical energy method. Proofs involve detailed study of the error equation for disturbances using the same weight function introduced in [<i>Comm. Math. Phys.</i>, 165 (1994), pp. 83--96].