fbe_1.m

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

function Be=fBe_1(ECC)
%三角单元应变矩阵函数
%%%%%%%%%%%%%%%%参数说明
%结点编号向量e单位cm
%%%%%%%%%%%%%%%%%%%%%%%%产生结点坐标
xi=ECC(1);
yi=ECC(2);
xj=ECC(3);
yj=ECC(4);
xm=ECC(5);
ym=ECC(6);
%%%%%%%%%%%%%%%%%%%%%%%控制程序
A=1/2*det([1 xi yi;1 xj yj;1 xm ym]);%单元面积cm^2
bi=-det([1 yj;1 ym]);
ci=det([1 xj;1 xm]);
bj=-det([1 ym;1 yi]);
cj=det([1 xm;1 xi]);
bm=-det([1 yi;1 yj]);
cm=det([1 xi;1 xj]);
Be=1/(2*A)*[bi 0 bj 0 bm 0;0 ci 0 cj 0 cm;ci bi cj bj cm bm]/32;