Not sure what you want, but if you define the set on the geometric edge rather than selecting the nodes, then the set will update when you resize/remesh/etc and you can also access the set with python.

Dear Soheil Hashemi ,

1. with the following code extract the nodes and their corresponding coordinates:

Node_List = ()

for I in range ("number of nodes")

t = mdb.models['Model_Name'].rootAssembly.instances['Part_Name-1'].nodes[i]

Node_List = Node_List + ((i,t),)

2. after extracting the nodes and their corresponding coordinates, with the code below you can select the nodes on the outer edge of the geometry (Here it is a circle):

Outer_Nodes = ()

R = "radius of the circle"

Delta = 0.01

for i in range (len(Node_List )):

X = Node_List[i][1][0]

Y = Node_List[i][1][1]

if X^2 + Y^2 = R^2 - Delta:

Outer_Nodes = Outer_Nodes + (Node_List, )

3. Now you have the labels and coordinates of the nodes that are on the outer edge of the geometry. After that, with the following code you can make a path with the node labels:

Nodes1= []

for i in range (len(Outer_Nodes)):

Nodes1 = Nodes1.append[Outer_Nodes[0]]

a = mdb.models['Model_Name'].rootAssembly

n1 = a.instances['Part_Name-1'].nodes a.Set(nodes=n1[Nodes1], name='Set_Name')

Best regards,

Ali.

Yo know the diameter of your disk. if the distance between all nodes and the centre of ypur disk is equal or larger than specific value, they will be your node sets. I think this would simply solve your problem.

If I am misunderstanding, please let me know. I might be abale to help you

Hi Soheil,

a possible solution for your specific problem is to define in your variables:

R = outer radius

toll = small value