To make this, I started by using the “morph duals by truncation” function in Stella 4d (available here) on a rhombicosidodecahedron. After converting the result to the base model, I then applied the “morph duals by expansion” function.

The second image shown here was made by applying “rainbow color mode” to the polyhedron above.