import part

mdb.models['Model-1'].parts['Part-1'].queryGeometry()

There will be the part’s bounding box.

Hi Luis,

def Max_z(mdbname):

op=mdb.models[mdbname].rootAssembly.instances

instances=[]

for instance in op.items():

instances.append(instance)

i=0

for instance in instances:

for node in instance[1].nodes:

i=i+1

Lst_z=[0]*i

i=0

for instance in instances:

for node in instance[1].nodes:

Lst_z[i]=node.coordinates[2]

i=i+1

Max_z=max(Lst_z)

return(Max_z)

I hope it works.