How to read multiple .odb files in Python to extract the reaction?

I have multiple .odb files for the same model. I would like to extract reaction forces at the set of nodes defined in each model. Please help me with the python script for this work.