Nastran Rod 杆单元质量和刚度矩阵

标准的耦合质量和刚度矩阵使用单元线性插值可以直接得出。


使用形函数为N=[1-x/l x/l]
在与nastran的比较中发现质量矩阵存在差异，遂探讨nastran rod单元质量矩阵的构成方式。
E杨氏模量 A面积 EleL单元长度 rho密度 EleC和EleS分别是单元指向的cos和sin
I=E*A/EleL;
EleK=[I 0 -I 0;
0 0 0 0;
-I 0 I 0;
0 0 0 0];

I=rho*A*EleL;
EleM=[I/12*5 0 I/12 0;
0 I/2 0 0;
I/12 0 I/12*5 0;
0 0 0 I/2];

EleT=[EleC EleS 0 0;
-EleS EleC 0 0;
0 0 EleC EleS;
0 0 -EleS EleC];

EleK=EleT'*EleK*EleT