EAch00EAch00012EIchhh6EIchh012EIchhh6EIchh06EIchh2EIch06EIchh4EIch斜杆的单元刚度矩阵HsqrthhLLko0EAoH00EAoH00012EIoHHH6EIoHH012EIoHHH6EIoHH06EIoHH4EIoH06EIoHH2EIoHEAoH00EAoH00012EIoHHH6EIoHH012EIoHHH6EIoHH06EIoHH2EIoH06EIoHH4EIoHkbT1kb0T1总体坐标下梁的单元刚度矩阵kcT2kc0T2总体坐标下柱子的单元刚度矩阵koT3ko0T3总体坐标斜杆的单元刚度矩阵Xzeros108108Yzeros108108Zzeros108108定义108阶0矩阵K1zeros108108K2zeros108108K3zeros108108K4zeros108108K5zeros108108K6zeros108108K7zeros108108K8zeros108108K9zeros108108把梁杆单元矩阵整合到总体刚度矩阵的循环语句forii111X3ii23ii33ii23ii3kbK1K1XXzeros108108e
dforii1323Y3ii23ii33ii23ii3kbK1K1YYzeros108108e
dforii2535Z3ii23ii33ii23ii3kbK1K1ZZzeros108108e
d把柱杆单元矩阵整合到总体刚度矩阵的循环语句forjj136K23jj23jj3jj23jjkc4646e
dforjj124K33jj23jj3jj23jjkc1313e
dforjj124
fK43jj23jj3jj343jj36kc1346e
dforjj124
K53jj343jj363jj23jjkc4613e
d把斜杆杆单元矩阵整合到总体刚度矩阵的循环语句forgg31227
K63gg23gg3gg23ggko4646e
dforgg21214
K73gg23gg3gg23ggko1313e
dforgg21214
K83gg23gg3gg373gg39ko1346e
dforgg21214
K93gg373gg393gg23ggko4613e
dKK1K2K3K4K5K6K7K8K9总体刚度矩阵Pzeros1081P11hqP371hqP731hq2P751qhh12AKP结构位移B1kbA1031081B2kbA67721B3kbA31361D1zeros61D1131A70721D1461A1061081D2zeros61D2131A34361D2461A70721C1kcD1C2kcD2C3kc4646A34361M111B161M211C161M311B261M411C261M511C131M611B361M711C331M811C231
fori13m36i2
fpri
tf第d层最右侧节点的位移是d
iAm1e
dfpri
tf第1层最右侧节点左侧杆的弯矩是f
M6fpri
tf第1层最右侧节点下侧杆的弯矩是f
M7fpri
tf第1层最右侧节点上侧杆的弯矩是f
M8fpri
tf第2层最右侧节点左侧杆的弯矩是f
M3fpri
tf第2层最右侧节点下侧杆的弯矩是f
M4fpri
tf第2层最右侧节点上侧杆的弯矩是f
M5fpri
tf第3层最右侧节点左侧杆的弯矩是f
M1fpri
tf第3层最右侧节点下侧杆的弯矩是f
M2
五、试算算例
f输入数据：输入单层高h：1输入单跨度L：1输入柱子的抗弯刚度EIc：1输入柱子的抗压刚度EAc：1输入梁的抗r