General modified gravity with 21cm intensity mapping: simulations and forecast