您所在位置：主题：[讨论]王瑁成FEATP程序的这个错误
whero_2000
[专家分：0] 发布于 2010-04-06 16:42:00
大家好，运用王瑁成FEATP中关于子空间迭代法子程序 SUBSPACE是有以下几行程序：
C*************** COMPUTE KI=XT* Y
C ***
DO 40 I=1,NVA
DO 40 J=1,NVA
GK(I,J) =0.0
DO 40 K= 1, NUMPT2
GK(I,J) =GK(I,J) +AA(K,I) * BB(K,J)
40 CONTINUE
其中GK是子空间刚度矩阵，它的表达式为：K=X（k+1）的转置*刚度矩阵*X（k+1）
但是上式中AA*BB中，BB代表的是上一次迭代出的特征向量矩阵，而不是刚度矩阵*X（k+1）
，所以应该改成GK(I,J) =GK(I,J) +AA(K,I) * AA(K,J)，不知小弟说的对不对？
C*************** COMPUTE KI=XT* Y
C ***
DO 40 I=1,NVA
DO 40 J=1,NVA
GK(I,J) =0.0
DO 40 K= 1, NUMPT2
GK(I,J) =GK(I,J) +AA(K,I) * BB(K,J)
40 CONTINUE
其中GK是子空间刚度矩阵，它的表达式为：K=X（k+1）的转置*刚度矩阵*X（k+1）
但是上式中AA*BB中，BB代表的是上一次迭代出的特征向量矩阵，而不是刚度矩阵*X（k+1）
，所以应该改成GK(I,J) =GK(I,J) +AA(K,I) * AA(K,J)，不知小弟说的对不对？