结构有限元程序.doc

有限单元法程序设计 姓 名：专 业：隧道工程2011级 学 号：指导教师： 目 录 第1章 平面刚架静力分析程序设计（PF.FOR）………………………1 第2章 平面三节点程序设计……………………………………17 第3章 四节点矩形薄板单元程序设计………………………………2 第1章平面刚架静力分析程序设计（PF.FOR） 用平面刚架静力分析程序分析所示结构。图示刚架材料为钢筋混凝土，杆件截面为矩形，柱截面宽0.4m，高0.4m；梁截面宽0.4m，高0.7m。材料的弹性模量E=3.55E7。 原始数据的准备与输入 （1）控制数据：对结构的结点、单元编号，并取结点1为原点，各控制数据如下表： E NM（杆件数） NJ(结点数) NS（支承条件） NLC（荷载情况数） 3.55E7 15 12 9 1 （2）结点坐标：(单位：m) 结点 1 2 3 4 5 6 x 0 0 0 0 6 6 y 0 4 7.5 11 11 7.5 7 8 9 10 11 12 6 6 12 12 12 12 4 0 11 7.5 4 0 (3)各单元始末端结点号及A、I值： 单元号 始端号 末端号 A截面积m2 I惯性矩m4 1 1 2 0.16 2.133E-3 2 2 3 0.16 2.133E-3 3 3 4 0.16 2.133E-3 4 4 5 0.28 11.43E-3 5 5 6 0.16 2.133E-3 6 6 7 0.16 2.133E-3 7 7 8 0.16 2.133E-3 8 5 9 0.28 11.43E-3 9 9 10 0.16 2.133E-3 10 10 11 0.16 2.133E-3 11 11 12 0.16 2.133E-3 12 3 6 0.28 11.43E-3 13 6 12 0.28 11.43E-3 14 2 7 0.28 11.43E-3 15 7 11 0.28 11.43E-3 （4）单元荷载 LM（受载杆件号） LT（荷载类型号） VF（荷载数值） DST（荷载终点至杆件始端距离） 1 3（水平均布荷载） 20 4 2 3 20 3.5 3 3 20 3.5 13 4（竖向均布荷载） -10 6 14 4 -10 6 将以上数据以下面的方式输入到程序中： ************************************************************* * * * shichunsheng gangjia 2009.6 * * * ************************************************************** 3.55E7 15 12 9 1 1 2 0.16 2.133E-3 2 3 0.16 2.133E-3 3 4 0.16 2.133E-3 4 5 0.28 11.43E-3 5 6 0.16 2.133E-3 6 7 0.16 2.133E-3 7 8 0.16 2.133E-3 5 9 0.28 11.43E-3 9 10 0.16 2.133E-3 10 11 0.16 2.133E-3 11 12 0.16 2.133E-3 3 6 0.28 11.43E-3 6 10 0.28 11.43E-3 2 7 0.28 11.43E-3 7 11 0.28 11.43E-3 0 0 0 4 0 7.5 0 11 6 11 6 7.5 6 4 6 0 12 11 12 7.5 12 4 12 0 11 0 12 0 13 0 81 0 82 0 83 0 121 0 122 0 123 0 0 5 1 3 20 4 2 3 20 3.5 3 3 20 3.5 13 4 -10 6 14 4 -10 6 程序运行后，输出结果为： ************************************************************* * * * shichunsheng gangjia 2009.6 * *