实例二
图示静定桁架，共14个结点，25根杆，3个支座反力。单位荷载分别沿上弦（2、4、6、8、10、12、14各点）和下弦（1、3、5、7、9、11、13各点）移动， 计算并绘制各杆件的内力影响线。
桁架数据输入如下，结果略。
XY1:=Matrix([[0,0],[0,3],[3,0],[3,3],[6,0],[6,3],[9,0],[9,3],[12,0],[12,3],[15,0],[15,3],[18,0],[18,3]]):
NE1:=Matrix([[1,2],[1,3],[2,3],[2,4],[3,4],[3,5],[4,5],[4,6],[5,6],[5,7],[6,7],[6,8],[7,8],[7,9],[7,10],[8,10],[9,10],[9,11],[9,12],[10,12],[11,12],[11,13],[11,14],[12,14],[13,14]]): NC1:=Matrix([[1,1],[1,2],[13,2]]):
- 10 -
第六章 结论
有限元法作为结构分析的一种有效方法已经被广泛采用，而有限元法与计算机程序的结合使得大型结构分析得以实现。通过以上实例，可以看出，通过有限元的方法编程计算桁架影响线的结果准确，与传统方式比较具有相当程度的优越性。效率高，结果精度高。在处理传统方式很难或无法求解的大型、复杂结构时，通过有限元程序也能够进行计算。当然程序并不是完美的，仍有许多地方可以完善，这些将在以后的工作中完成。
参考文献
[1]科学出版社 《maple9.0符号处理及应用》 黎捷 编著 [2]高等教育出版社 《结构力学I》 龙驭求 包世华 主编
- 11 -
附录1
XY:=Matrix([[0,0],[0,6],[3,0],[3,3],[3,7],[6,0],[6,6],[6,8],[9,0],[9,9],[12,0],[12,6],[12,8],[15,0],[15,3],[15,7],[18,0],[18,6]]):
NE:=Matrix([[1,3],[1,4],[1,2],[2,5],[3,4],[4,5],[3,6],[4,6],[4,7],[5,8],[6,7],[7,8],[6,9],[7,9],[7,10],[8,10],[9,10],[9,11],[9,12],[10,12],[10,13],[11,12],[12,13],[11,14],[11,15],[12,15],[13,16],[14,15],[15,16],[14,17],[15,17],[16,18],[17,18]]): NC:=Matrix([[1,1],[1,2],[17,2]]):
P1:=Matrix([[2],[5],[8],[10],[13],[16],[18]]): P2:=Matrix([[1],[3],[6],[9],[11],[14],[17]]):
??13????14?????12?????25??????34??????00????45?????06???36????????30????46??????33??????47???????37???58???????67????60???????78????66??????2???5????68????69?????1??79????8???XY := ???90??????710????NC := ?1?2???????1?99?????10??????17?P1 := ?10???2???120???????NE := ?8????13???910???????126????????16????911???128?????
??18?????912??????150?????1012????153???????1013??????157???1112????????180???1213?????6??????18??1114?????1115?????1215?????1316??????1415??????1516????1417??????1517??????1618?????1718??? ??1???3????6???P2 := ????9??????11? ?????14????17????- 12 -
附录2
P1路径下各杆件影响线 总体内力矩阵{RR}：
?0.??0.???-1.??0.???0.???0.??0.???0.???0.???0.??0.???0.???0.???0.??0.???0.???0.???0.??0.???0.???0.???0.??0.???0.???0.???0.??0.???0.???0.??0.???0.???0.???0.??0.???1.???0.?.8333-1.1790.0.0.-1..8333-.7072-.47150..50000..3333-.3727-.23570..3333.16670.-.23570.0.0..16670.-.23570.0.0..1667-.23570.0.0..8333.1667.6667-.94300.0.0.0..66670.-.94300.0.-1..6667-.7454-.47150..6667.33330.-.47150.0.0..33330.-.47150.0.0..3333-.47150.0.0..6667.3333.5000-.70720.0.0.0..50000.-.70720.0.0..50000.-.70720.0..50000.-.70720.0.0..50000.-.70720.0.0..5000-.70720.0.0..5000.5000.3333-.47150.0.0.0..33330.-.47150.0.0..33330.-.47150..6667.6667-.7454-.47150.0.-1..66670.-.94300.0.0..6667-.94300.0.0..3333.6667.1667-.23570.0.0.0..16670.-.23570.0.0..16670.-.23570..3333.3333-.3727-.23570..50000..8333-.7072-.47150.0.-1..8333-1.1790.0.0..1667.83330.??0.???0.??0.???0.???0.??0.???0.???0.???0.??0.???0.???0.???0.??0.???0.???0.???0.??0.???0.???0.???0.??0.???0.???0.???0.??0.???0.???0.??0.???0.???0.???-1.??0.???0.???1.?? - 13 -
- 14 -