Latest trunk version of pymunk has two new (as of today) shapes and constraints properties on the Body class. I went with always return the shapes/constraints regardless of if they are added to the space or not.
So now you can just do:
>>> import pymunk
>>> b = pymunk.Body()
>>> b2 = pymunk.Body()
>>> j = pymunk.PivotJoint(b,b2,(0,0))
>>> s = pymunk.Circle(b,3)
>>> b.constraints
set([<pymunk.constraint.PivotJoint object at 0x02521890>])
>>> b.shapes
set([<pymunk.Circle object at 0x025218F0>])
(This will be included in pymunk 3.1.)