关于节点振型与刚度矩阵的对应关系以求模态应变能

tang_jian

我想用公式计算出关于单元的一个指标值MSE=L*K*LT, 其中L代表单元的振型，K代表单元的刚度矩阵，LT代表振型的装置。遇到这样一个问题，先贴个图：

如图所示，我使用了SHELL63单元，每个节点有6个自由度，单元1的刚度矩阵K1，以及每个节点的振型（L1,L2,L3,L4,其中L1=[x1,y1,z1,rx1,ry1,rz1]，其余类似)，我都已经提取出来了，现在的问题是，用上面的那个公式计算MSE时，是不是应该有个顺序问题，是用[L1,L2,L3,L4]*K1*[L1,L2,L3,L4]T,还是用[L4,L3,L2,L1]*K1*[L4,L3,L2,L1]T ？？？就是说它的对应关系是怎样的

tang_jian

图：

                               
登录/注册后可看大图

tang_jian

都快沉底了 还没人来救我啊