摘要：A mathematical method was developed to calculate the yield,energy spectrum and angular distribution of neutrons from D(d,n)~3He(D-D)reaction in a thick deuterium-titanium target for incident deuterons in energies lower than *** data of energy spectrum and angular distribution were applied to set up the neutron source model for the beam-shaping-assembly(BSA)design of Boron-Neutron-Capture-Therapy(BNCT)using MCNP-4C code. Three cases of D-D neutron source corresponding to incident deuteron energy of 1000,400 and 150 key were inves- *** neutron beam characteristics were compared with the model of a 2.45 MeV mono-energetic and isotropic neutron source using an example BSA designed for BNCT *** results show significant differences in the neutron beam characteristics,particularly the fast neutron component and fast neutron dose in air,between the non-isotropic neutron source model and the 2.5 MeV mono-euergetic and isotropic neutron source model.