I'm fairly new to ABAQUS and I am curious if there is an easy way to sort elements based upon the local value of an output variable. Specifically, I want to select elements with a certain range of density values returned by a topology optimization. However, I imagine this could also be applied to selecting only high-stress elements or similar. I want to be able to specify custom upper and lower bounds of the desired range. The end goal is to export a mesh consisting of only elements in this range.

I imagine there is a way to script this, (e.g: for element in part_1: if element.density > 0.2 and element.density < 0.6, add element to set_1), but I am still learning the ropes of how to write scripts to interact with abaqus, and any advice would be super useful.