Considering the CP violating phases, we analyze the neutron electric dipole moment (EDM) in a CP violating supersymmetric extension of the standard model where baryon and lepton numbers are local gauge symmetries(BLMSSM). The contributions from the one loop diagrams and the Weinberg operators are taken into account. Adopting some assumptions on the relevant parameter space, we give the numerical results analysis. The numerical results for neutron EDM can reach $1.05times 10^{-25}(e.cm)$, which is about the experimental upper limit.