Juneliu 发表于 2009-3-25 10:35

提取刚度矩阵的问题

论坛里有个提取刚度矩阵的贴子,其中可以用/debug命令提取某单元的单元刚度矩阵,我想知道的是此时 提取的单元刚度矩阵是整体坐标系下的,还是单元局部坐标系下的刚度矩阵,也就是说若用它集成整刚还需不需要乘上坐标转换矩阵???
另外,计算某单元的第i阶模态应变能时,是不是可以直接利用该单元刚度矩阵和第i阶整体坐标系下的各个单元节点位移模态值组成的振型来求解(而不把该该单元刚度矩阵集成整体矩阵的阶数),举个例子来说:那对于一个由3维beam44单元组成的系统(假设系统共100个自由度) ,K是否就是12X12阶的矩阵?Ф则是整体坐标系下的振型,是由该单元两端节点位移模态值构成的12X1的向量,通过MSE=Ф'KФ就能计算出第i阶单元的模态应变能?

还有,能否在ansys中直接提取模态应变能?即在第i阶模态下,提取某单元的应变能(SENE)(从table项get),这个值是不是等于该阶模态下的模态应变能?我算了个例子,发现两者值就是2倍的关系。好像是对的。
页: [1]
查看完整版本: 提取刚度矩阵的问题