We present a study of weak lensing shear measurements for simulated galaxy images at radio wavelengths. We construct a simulation pipeline into which we can input galaxy images of known shapelet ellipticity, and with which we then simulate observations with eMERLIN and the international LOFAR array. The simulations include the effects of the CLEAN algorithm, uv sampling, observing angle and visibility noise, and produce realistic restored images of the galaxies. We apply a shapelet-based shear measurement method to these images and test our ability to recover the true source shapelet ellipticities. We model and deconvolve the effective point spread function, and find suitable parameters for CLEAN and shapelet decomposition of galaxies. We demonstrate that ellipticities can be measured faithfully in these radio simulations, with no evidence of an additive bias and a modest (10 per cent) multiplicative bias on the ellipticity measurements. Our simulation pipeline can be used to test shear measurement procedures and systematics for the next generation of radio telescopes.