You don’t have any subdomain information extracted into your mesh. In fact you have commented out the cell_data field that is required to get information about subdomains/boundaries which you can then pass on to MeshValueCollection/MeshFunction and so on…
See Transitioning from mesh.xml to mesh.xdmf, from dolfin-convert to meshio - #104 by dokken and the comment above it to get an idea on how to do if for your Abaqus mesh.