A　secondary　multidirectional　vibration　isolation　system　(MR-SMDIS)　was　designed　by　utilizing　the　multidirectional　vibration　isolation　characteristics　of　metal　rubber　(MR)　vibration　isolators.　According　to　the　spatial　distribution　angles　and　contact　form　of　the　spatial　distribution　of　metal　wires　in　MR,　the　linear　stiffness　ratio　between　axial　and　radial　directions,　damping　ratio,　and　the　ratio　of　nonlinear　stiffness　to　linear　stiffness　were　analyzed.　The　nonlinear　response　of　an　MR-SMDIS　under　random　vibrations　was　analyzed　using　the　extended　harmonic　equilibrium　method　with　alternating　frequency–time　domain　(HB–AFT).　Random　signals　were　further　created　to　analyze　the　isolation　efficiency　and　peak　amplification　of　the　MR-SMDIS.　The　particle　swarm　optimization　(PSO)　was　used　to　optimize　the　axial　and　radial　linear　stiffness,　damping,　and　nonlinear　stiffness　of　the　MR　isolator.　And　the　optimal　preparation　parameters　for　MR　were　solved　with　a　genetic　algorithm.　Furthermore,　random　vibration　experiments　of　the　MR-SMDIS　were　conducted.　The　results　showed　that　the　relative　error　in　identifying　the　dynamic　parameters　of　the　prepared　MR　isolators　using　parameter　identification　algorithms　was　within　±　13%　compared　to　the　design　parameters.　The　MR-SMDIS　met　the　design　requirements　regarding　modal　frequencies　and　peak　amplification　factors.　The　relative　error　between　the　experimental　and　numerical　simulation　results　was　within　15%.　Graphical　abstract:　(Figure　presented.)　©　The　Author(s),　under　exclusive　licence　to　Springer　Nature　B.V.　2024.