RHS殘余力矢量RHS計(jì)算方法一:使用內(nèi)置函數(shù) 上策RHS(1:NDOFEL,1) = RHS(1:NDOFEL,1) - MATMUL(AMATRX,U(1:NDOFEL))方法二:中策DO K1=1,8 DO K2=1,8 RHS(K1,1)=RHS(K1,1)- AMATRX(K1,K2)*U(K2) END DOEND DO 方法三:下策
form and assemble stiffness matrix and internal force vector amatrx(1:16,1:16)=amatrx(1:16,1:16)+dvol* 2 ((hphi+xkap)*matmul(matmul(transpose(b),ddsdde),b)) rhs(1:16,1)=rhs(1:16,1