Following the previous post, I have added a further new function, and in the process discovered that the py_Estress function was not handling tensile axial loads correctly. That has now been fixed, and the updated spreadsheet and Python files can be downloaded from:
The new function calculates creep and shrinkage effects on section curvature for a rectangular reinforced or prestressed concrete section, under any specified bending moment and axial load.
See the download file for full details of the required input and an example:

As before, for details of the pyxll app required to connect the Python code to Excel see Python and pyxll.