Normal vector of just one of the boundaries

How are you going to transform this into a numpy array?

I would suggest using a projection into a suitable function space, which you can do as:

which shows you how to get a projection of the facet normal over a certain boundary marked with a mesh tag.
You can of course adapt this to take in your stress and project in to a suitable space.