请教关于模态应变能法计算的若干问题

Juneliu

问题一、在计算某单元的模态应变能时，MSE=Ф'KФ中的K和Ф分别是单元刚度矩阵和结构第i阶模态振型，他指的是局部坐标系下的值，还是整体坐标系下的值？
举个例子来说，那对于一个由3维beam44单元组成的系统（假设共100个自由度) ，K是否就是12X12阶的矩阵？Ф则是整体坐标系下的振型，是由该单元两端节点位移模态值构成的12X1的向量？对不？

问题二、能否在ansys中直接提取模态应变能？即在第i阶模态下，提取某单元的应变能(SENE)(从table项get)，这个值是不是等于该阶模态下的模态应变能？不知道有没有人做过？
我算了个例子，发现两者值差距比较大。不知道是不是错了？

请指教,谢谢!!

[ 本帖最后由 Juneliu 于 2009-3-24 13:36 编辑 ]