* Extraction of HU values is now via nodes, and then interpolated to the integration points via the element shape function, rather than directly extracting the HU values to the integration points.
* Abaqus example job files added for a shoulder model, including user subroutine
* Updated readme to render correctly on PyPi