Abstract: The study of secondary neutrons produced by light charged particles induced reaction is of great significance to the design and optimization of accelerator shielding. The neutron double differential yields from 33 MeV-d, 65 MeV-³He and 65 MeV-⁴He bombarding thick carbon, copper and lead targets in the directions of 0^\circ, 15^\circ, 45^\circ, 75^\circ, 135^\circ, are calculated by using Geant4 code with INCL, BIC and BERT physics models, and compared with the experimental data. The results show that, for the 33 MeV d-induced reaction, the results from the INCL model basically reproduced the experimental data of carbon and copper targets, but overestimated the neutron yields corresponding to the direct process in lead target. The results from the BIC model and Bert models failed to reproduce the broad peak due to the stripping of the projectiles. For the 65 MeV ³He-induced reaction, the results from the three models could not reproduce the neutrons in the stripping process in the forward angles. The calculated ones with these three models agreed well with the experimental results in other angles. For the 65 MeV ⁴He-induced reactions, the calculation results of the INCL model are in good agreement with the experimental ones for carbon and copper targets, but underestimated the neutrons for lead target. The calculation results of the BIC model and BERT model underestimated the experimental data for carbon target and slightly overestimated the experimental ones for lead target in large angles.