I think this is not possible. To combine any two arrays, they must have the same dimensions. And any two dimensions in each array must be of the same size.
You can imagine (60,000, 28, 28) array as a cube. The surface looking at you has the dimension of 28 x 28. Thus, all same-size surfaces behind it are 60,000 in number. If you want to add a new entity to it, it must have the same 3-D dimension. And at least two dimensions must match those of the first cube. Otherwise, it won't get concatenated exactly.
To combine (60,000, 28, 28) with another array, the second array should have any two of 60,000, 28, 28 as its dimensions. Let's suppose, the second one has (60,000, 28, 14). Then, you can concatenate and get the result:
z = np.concatenate((array1, array2), axis=2)
z.shape
Output:
(60000, 28, 42)
Alternatively, if the second array is (30,000, 28, 28):
z = np.concatenate((array1, array2), axis=0)
z.shape
Output:
(90000, 28, 28)
np.arange((24).reshape(2,3,4). What does it mean to a a size (2,) array to that?