On the basis of Hindmarsh-Rose neuron model, dynamical robustness and the transition of firing modes of multilayer memristive neural network consisting of nonidentical neurons have been investigated in detail. The dynamic effects of memristive synapses coupling which are either cubic order flux or quadratic flux are detected. Our results suggest that for the cubic order flux-controlled memristors increasing electromagnetic induction parameter and parameter of memristor show the tendency to spoil the dynamical robustness of the multilayer memristive neural network. For quadratic flux-controlled memristors, weak electromagnetic induction can spoil the dynamical robustness while strong electromagnetic induction and parameter of memristor have little influence on dynamical robustness. We found that the ratio of inactive neurons switches the firing patterns of the active neuron among periodic bursting, chaotic bursting and spiking-like. In the case of memristive synapse coupling by cubic order flux-controlled memristors, chaotic bursting alternates with period bursting as the ratio of inactive neurons is increased, and the chaotic bursting occupy a large range of the ratio of inactive neurons. In the case of memristive synapse coupling by quadratic flux-controlled memristors, the main firing pattern of the active neuron is chaotic bursting. The distribution of interspike interval can be broadly classified into two groups: those that are long interspike intervals and those that are short interspike intervals. The dynamics of multilayer memristive neural network is also verified in the circuits built on Multisim. These results could give new mechanism explanation for aging transition by applying electromagnetic field to neural network.