You are here
ABAQUS Element Connectivity
Hi,
I am trying to mesh the hollow sphere attachged using structured BRICK-8 elements. It turns out in order to have a structured mesh, I have to have 2 partitions such that four quadrants be generated.
What I only need from ABAQUS is the element connectivity. I do not need to solve any problem.
The issue I am facing is that each quadrant connectivity nodes' order is different. e.g. for a simple one element case with coordinates (-1,-1,-1)-(1,-1,-1)-..-(-1,1,1) you expect you number the elements in the datum plane with negative z-axis starting 1 to 4. Then you sweep the positive z-axis plane and and number them 5-8. Now if you have any other elements in the domain, the starting node in the connectivity should always be consistent by this order 1-2-3-4-5-6-7-8 however for two quadrants I have 1-2-3-4-5-6-7-8 and for the other two I have 3-4-1-2-7-8-5-6. I think this problem is due to the partitions and is inevitable.
Please let me know if you know the answer.
Cheers,
Reza
| Attachment | Size |
|---|---|
 mesh.png | 231.75 KB |
- reza.yaghmaie's blog
- Log in or register to post comments
- 12662 reads
Comments
script for element connectivity
Hello
the code below will generate a file listing (in my example)
({'connectivity': (4980, 4981, 4982, 4983), 'instanceName': None, 'label': 1, 'type': C3D4}),4980
for all elements.
The 4 numbers are the nodes forming the element identified after 'label'
I do not recall what the last number stands for, it presumably repeats the first node number
Important: node numbers are shifted by one ! A node number 0 exists.
Check the output for a few elements by picking
ABAQUS/CAE : Tools-->Query-->Element
Good luck
Frank
----------------------------------------------------------------
import part
file = open('C:/Simulations/Shallow Notch/SN08/SN08-elements.txt','w+') # path and filename you want to generate
elemArr = mdb.models['Model-1'].parts['UFTE_SN_ohne_Gewinde'].elements # pick the part you want to analyse
len(elemArr)
# print elemArr
print len(elemArr)
#a = mdb.models['Model-1'].parts['Model-1']
#print a
#allNodes = a.nodes
#print allNodes
#rootAssembly
#elemente = a.elements[1]
# for e in elemArr:
# print e.getNodes()
# print e,' : ', e.connectivity
for e in elemArr:
file.write('%s,%g\n' % (e, e.connectivity[0]))
file.close()