Manually , I do it using the option " by angle" and clicking on the face. But I am looking to automate this as I have to deal with multiple simulations
I used the following script :
rind = [[i+1 for i in range(len(elementConnectivity)) if item in elementConnectivity[i]]
for item in TOP]
rind = [item for sublist in rind for item in sublist]
rind=list(set(rind))
SFE = myPart.SetFromElementLabels('SFE', rind)
SFelems =SFE.elements
myPart.Surface(face1Elements=SFelems, name='Surf-11')
Where,
elementconnectivity is the element connectivty matrix
TOP is the list of nodes that belong to the top face of my cube
But I am not able to obtain the result I am looking for, Is there an alternative approach ?