If your polygon has a uniform density, its center of mass coincides with its centroid. In shapely, the centroid can be directly calculated as:
from shapely.geometry import Polygon
P = Polygon([[0, 0], [1, 0], [1, 1], [0, 1]])
print(P.centroid)
#POINT (0.5 0.5)