Hi! When trying to write Python code to determine the external surfaces, it does not distinguish the element faces per element, as you can see on the screenshot. I was wondering if anyone knows what to do? I wrote the following script, which can be seen below:

first = len(mdb.models['Model-1'].parts['Part-1'].elements[1].getAdjacentElements())

second = len(mdb.models['Model-1'].parts['Part-1'].elements[1].getElemFaces())

if firstgetElements())

if third == 1:

   print('External surface')

   p = mdb.models['Model-1'].parts['Part-1']

   f = p.elements

   face3Elements = f[q:q+1]

   p.Surface(face3Elements=face3Elements, name='Surf-20')

else:

   print('Internal surface')

Abaqus