How do I calculate rise times for a push pull GPIO with load with hand calculations?
For this, and if you have access to it, you should better use IBIS models for the I/O cell.

The IBIS model not only provide for rise time and fall time modelling but also includes ESD diodes modelling which will make your simulations much more accurate (see Fig 9. for diagram).