导热问题有限元刚度矩阵二维平面四节点

导热问题有限元刚度矩阵二维平面四节点

二维平面4节点等参有限元的热传导刚度矩阵计算很简单，这里直接给了计算结果

function KE = elementMatVec2D(D0)
KE1 = D0(1)*[2 -2 -1 1;-2 2 1 -1;-1 1 2 -2;1 -1 -2 2];
KE2 = D0(2)*[2 1 -1 -2;1 2 -2 -1;-1 -2 2 1;-2 -1 1 2];
KE = (KE1+KE2)/6;
end

---

转载请注明来源，欢迎评论区讨论