The uncertain nonlinear systems can be modeled with fuzzy equations or fuzzy differential equations (FDEs) by incorporating the fuzzy set theory. The solutions of them are applied to analyze many engineering problems. However, it is very difficult to obtain solutions of FDEs.In this paper, the solutions of FDEs are approximated by two types of Bernstein neural networks. Here, the uncertainties are in the sense of Z-numbers. Initially the FDE is transformed into four ordinary differential equations (ODEs) with Hukuhara differentiability. Then neural models are constructed with the structure of ODEs. With modified back propagation method for Z-number variables, the neural networks are trained. The theory analysis and simulation results show that these new models, Bernstein neural networks, are effective to estimate the solutions of FDEs based on Z-numbers.