I am not getting the solution from the above posts you have given.
The reason why I commented out the cell_data field because there is no Cell data: gmsh:physical, gmsh:geometrical in my mesh_from_file because it is the Abaqus file, that data will be present in the GMSH mesh file.
The other thing which I am not able to follow from the given post Transitioning from mesh.xml to mesh.xdmf, from dolfin-convert to meshio - #104 by dokken
is facet_data = mesh.cell_data_dict["gmsh:physical"]["line"], because this dictionary is empty in my case.
I request you to please provide me the code if you have it in which the Abaqua mesh file is deal with and following line is worked successfully.
with XDMFFile("facet_mesh.xdmf") as infile:
infile.read(mvc, "name_to_read")
mf = cpp.mesh.MeshFunctionSizet(mesh, mvc)
ds = Measure("ds", domain=mesh, subdomain_data=mf)
So that, I will be able to apply boundary conditions.
Thanks in advance