ansys导出单元刚度与自己编程的不一致?

liusing

FINISH
/CLEAR
/PREP7!　Enter　PREP7
k,1,0.,0.,
k,2,0.1,0.,
k,3,0.1,0.1,
k,4,0,0.1,
a,1,2,3,4
ET,1,PLANE42
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,1
MPDATA,PRXY,1,,.25   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,2560
allsel,all
LESIZE,all, , ,100, , , , ,1
MSHAPE,0,2D
MSHKEY,0
amesh,1
finish

/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件  
/debug,-1,,,1 ! 指定输出单元矩阵


/SOLU
ANTYPE,2
MODOPT,SUBSP,10
MXPAND,10
NSEL,S,LOC,Y,0.
D,ALL,UY,0.
NSEL,R,LOC,X,0.05
D,ALL,UX,0.
ALLSEL,ALL
SOLVE
FINISH
/POST1
SET,LIST

导出的单元刚度见图1
而计算的单刚见图2
单元为平面四节点等参元，计算时的单元形状大小物理参数完全与ansys的一致，但是刚度矩阵有所不同，这是为什么呢？另外计算应该没有问题反复核对了，计算采用曾攀的《有限元分析基础教程》中平面四节点等参元刚度矩阵的显示表达，见图3和图4.
我是用来计算特征值问题的，不知道是不是ansys对单刚矩阵是不是进行了什么处理？

liusing

问题解决！
参看帖子 A04：ABAQUS » 自定义刚度矩阵-问题已经解决！！！！！
看来还是形函数的问题！
“是ansys 输出的单元刚度矩阵出现了问题！
( K. O/ e% {& p4 w$ z  R5 b+ }3 Z; r. ^
ansys中的plane42单元采用平面应变，但是形函数查看理论手册却采用的是有额外形函数的形函数形式 ，可参见12.7.4. 2-D and Axisymmetric 4 Node Quadrilateral Solids with ESF (QM6)These shape functions are for the 2-D 4 node and axisymmetric solid elements with extra shape functions。”——引自《自定义刚度矩阵-问题已经解决！！！！！》