单元刚度矩阵，桁架

clear all
%i点坐标为（xi，yi）一样的方法假设j、m点，A为面积这里为1
syms xi yi xj yj xm ym u
b=[yj-ymym-yiyi-yj]
c=[-xj+xm-xm+xi-xi+xj]
%假设常数Et/（4*（1-u^2））为h或者为1，为方便计算为1
for （i=1:3）
    for （j=1:3）
        k1（ij）=b（i）*b（j）+（1-u）*c（i）*c（j）/2
        k2（ij）=u*b（