fge_1.m

来自「在有限元分析中计算结构整体刚度矩阵、单元刚度矩阵、单元结点位移、应力应变等」· M 代码 · 共 10 行

function Ge=fGe_1(n,ecn)
%三角单元结点自由度转换矩阵函数
%结构结点个数n;单元结点编号向量ecn
l=6;
I=eye(2);
Ge=zeros(l,2*n);
Ge(1:2,(2*ecn(1)-1):(2*ecn(1)))=I;
Ge(3:4,(2*ecn(2)-1):(2*ecn(2)))=I;
Ge(5:6,(2*ecn(3)-1):(2*ecn(3)))=I;