Fully implicit, stabilised, three-field material point method for dynamic coupled problems