各位高手，小弟在进行abaqus三维单元刚度阵程序编写时，遇到个问题，请进！！谢谢！！
为了编写三维单元的刚度阵，我需要知道abaqus程序中单元的节点信息，即abaqus中单元节点顺序，以便编写单元的形函数进而编写单刚矩阵。我在读inp文件时发现这样的情况（举个例子如下）
*Element, type=C3D20R
1, 10, 11, 14, 13,1,2,5,4, 31, 30, 29, 28, 32, 33, 34,
35, 37, 36, 38, 39
2, 11, 12, 15, 14,2,3,6,5, 42, 41, 40, 30, 43, 44, 45,
33, 36, 46, 47, 38
3, 13, 14, 17, 16,4,5,8,7, 29, 50, 49, 48, 34, 51, 52,
53, 39, 38, 54, 55
4, 14, 15, 18, 17,5,6,9,8, 40, 57, 56, 50, 45, 58, 59,
51, 38, 47, 60, 54
5, 19, 20, 23, 22, 10, 11, 14, 13, 64, 63, 62, 61, 31, 30, 29,
28, 66, 65, 67, 68
6, 20, 21, 24, 23, 11, 12, 15, 14, 71, 70, 69, 63, 42, 41, 40,
30, 65, 72, 73, 67
7, 22, 23, 26, 25, 13, 14, 17, 16, 62, 76, 75, 74, 29, 50, 49,
48, 68, 67, 77, 78
8, 23, 24, 27, 26, 14, 15, 18, 17, 69, 80, 79, 76, 40, 57, 56,
50, 67, 73, 81, 77
这是一个三维20节点的单元信息 一共有8单元每个单元的节点编号都在这里 而且每个单元之间节点的排列顺序是对应的 这是不是代表在abaqus里求节点形函数时就是按照这个节点排列顺序来求的？
比方说 第一个单元的节点顺序为10, 11, 14, 13,1,2,5,4, 31, 30, 29, 28, 32, 33, 34, 35, 37, 36, 38, 39 而第一个节点（10）的坐标反映在局部坐标里是（1， -1，1） 第二个节点（11）的局部坐标是（1,1,1,）依次可写出1单元的其他节点局部坐标 而不同单元中每个节点的局部坐标都与1单元是对应的即单元1中节点10与单元2中节点11的局部坐标相同那么它们的形函数N1=1/8*(1+X)*(1-Y)*(1+Z)*(X-Y+Z-2)而1单元的第二个节点是11 其局部坐标为（1,1,1）故第二个形函数可写为N2=1/8(1+X)*(1+Y)*(1+Z)*(X+Y+Z-2) 以此类推 每个单元也是如此那么在编程的过程中是不是就按这样的单元节点顺序进行？abaqus果真是这样进行节点的计算的吗？ 我是这样理解的 但是没法求证 希望高手们能帮帮我小弟不胜感激！！！！！ 高手为啥不出现呢？{:3_64:}
页:
[1]