如何计算单元切线刚度矩阵

我现在在非线性有限元做结构力学分析，遇到关于计算索单元的非线性的应变矩阵，最后要得到关于单元切线刚度矩阵的计算，因为ansys和abaqus没有关于索单元的，所以这一块得自己编程计算，想请教如何进行矩阵的积分计算?

矩阵积分还是每个元素进行积分的，最常用的就是gauss积分了，把函数表达出来，写一个gauss积分的subroutine然后算就行了，基本上是计算在一些节点的值，乘上相应的权重然后相加即可。你可以wiki gauss quadrature或高斯积分

单元编号不会影响单元刚度矩阵。单元刚度矩阵里包含单元节点的坐标数据，因此节点编号会对单刚里的数值所在的位置有影响。集成总刚的时候也是按节点编号集成的。 但不管如何编号，对计算的最终结果都是一样的。