I'm trying to get the element density from the abaqus output database. I know you can request a field output for the volume using 'EVOL', is something similar possible for the density?

I'm afraid it's not because of this: Getting element mass in Abaqus postprocessor

What would be the most efficient way to get the density? Look for every element in which section set it is?

I thought I found the solution but in the end this doesn't work either:

odb_file_path=r'your_path\file.odb'
odb = session.openOdb(name=odb_file_path)

instance = odb.rootAssembly.instances['MY_PART']
material_name = instance.elements[0].sectionCategory.name[8:-2]
density=odb.materials[material_name].density.table[0][0])

note: the 'name' attribute will give you a string like, 'solid MATERIALNAME'. So I just cut out the part of the string that gave me the real material name. So it's the sectionCategory attribute of an OdbElementObject that is the answer.

EDIT: This doesn't seem to work after all, it turns out that it gives all elements the same material name, being the name of the first material.