设整体坐标下单元刚度矩阵为，杆端力列阵为，杆端位移列阵为，杆件固端力列阵为，则有