LINK1的单元质量阵怎么和我编的FORTRAN程序算的不一样啊？

kfcmmdb

RT,刚度矩阵是一样的

下面是DEBUG，输出的
    MASS MATRIX FOR ELEMENT       1
   1   0.4714045E+00  0.0000000E+00  0.2357023E+00  0.0000000E+00
   2   0.0000000E+00  0.4714045E+00  0.0000000E+00  0.2357023E+00
   3   0.2357023E+00  0.0000000E+00  0.4714045E+00  0.0000000E+00
   4   0.0000000E+00  0.2357023E+00  0.0000000E+00  0.4714045E+00

   MASS MATRIX FOR ELEMENT       2
   1   0.3333333E+00  0.0000000E+00  0.1666667E+00  0.0000000E+00
   2   0.0000000E+00  0.3333333E+00  0.0000000E+00  0.1666667E+00
   3   0.1666667E+00  0.0000000E+00  0.3333333E+00  0.0000000E+00
   4   0.0000000E+00  0.1666667E+00  0.0000000E+00  0.3333333E+00

   MASS MATRIX FOR ELEMENT      3
   1   0.4714045E+00  0.0000000E+00  0.2357023E+00  0.0000000E+00
   2   0.0000000E+00  0.4714045E+00  0.0000000E+00  0.2357023E+00
   3   0.2357023E+00  0.0000000E+00  0.4714045E+00  0.0000000E+00
   4   0.0000000E+00  0.2357023E+00  0.0000000E+00  0.4714045E+00

我自己FORTRAN程序算的
MASS MATRIX FOR ELEMENT       1
   1   0.4714045E+00  0.4714045E+00  0.2357023E+00  0.2357023E+00  
   2   0.4714045E+00  0.4714045E+00  0.2357023E+00  0.2357023E+00
   3   0.2357023E+00  0.2357023E+00  0.4714045E+00  0.4714045E+00
   4   0.2357023E+00  0.2357023E+00  0.4714045E+00  0.4714045E+00

MASS MATRIX FOR ELEMENT       2
   1   0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
   2   0.0000000E+00  0.3333333E+00  0.0000000E+00  0.1666667E+00
   3   0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
   4   0.0000000E+00  0.1666667E+00  0.0000000E+00  0.3333333E+00

MASS MATRIX FOR ELEMENT       3
   1   0.4714045E+00      －0.4714045E+00   0.2357023E+00     －0.2357023E+00
   2   －0.4714045E+00  0.4714045E+00      －0.2357023E+00   0.2357023E+00
   3   0.2357023E+00      －0.2357023E+00  0.4714045E+00      －0.4714045E+00  
   4   －0.2357023E+00  0.2357023E+00      －0.4714045E+00   0.4714045E+00


计算质量阵元素的程序块
LENTH1=SQRT((X2-X1)*(X2-X1)+(Y2-Y1)*(Y2-Y1))
C=(X2-X1)/LENTH1
S=(Y2-Y1)/LENTH1
F=BCAREA*BDENS*LENTH1/6
EM1(1,1)=2*F*C*C
EM1(2,1)=2*F*C*S
EM1(2,2)=2*F*S*S
EM1(3,1)=F*C*C
EM1(3,2)=F*C*S
EM1(3,3)=2*F*C*C
EM1(4,1)=F*C*S
EM1(4,2)=F*S*S
EM1(4,3)=2*F*S*C
EM1(4,4)=2*F*S*S

kfcmmdb

/PREP7  
N,1
N,2,-1,-1   
N,3,0,-1
N,4,1,-1
ET,1,LINK1  
MP,EX,1,1   
MP,DENS,1,1
MP,DENS,2,1
MP,EX,2,1   
R,1,1   
MAT,2   
REAL,1  
E,1,2   
E,1,4   
MAT,1   
REAL,1  
E,1,3   
FINISH  
/OUTPUT,cp,out,, ! 将输出信息送到cp.out  
/debug,-1,,,1 !
/SOLU   
SOLVE   
finish  
/OUTPUT, TERM ! 将输出信息送到output windows

[ 本帖最后由 kfcmmdb 于 2007-12-13 18:39 编辑 ]