说明：  对构件进行优化分析时，需要将其剖分成四面体网格单元，本源码可以用来求每一个小四面体的刚度矩阵
(Optimization analysis of the components need to be partitioned into tetrahedral grid cell, the source can be used to seek the stiffness matrix of each small tetrahedron)