结构力学矩阵位移法上机做法详解

1、构力学矩阵位移法上机实践报告指导老师：韩衍 群班级：工程力学0901班姓名：学号：但是还是写一下吧，红字部分为解释部分，黑字为内容）矩阵位移法上机实践报告一、实践任务掌握矩阵位移法, 能够利用电算法计算各种结构的内力. 计算各种典型结构的内力 , 主要包括钢架, 连续梁和桁架。计算完成之后, 写出实践报告。二、实践要求（ 1）用程序计算给定结构的内力；（ 2）根据计算结果手工绘制出来各个结构的内力图。三、实践方法（ 1）为节点编号：选定适当的整体坐标系和局部坐标系，为各个结构的节点和杆件编号。（ 2）建立输入文件：根据已知数据，建立格式为文本文档的输入数据文件，输入文件需要严格按照格式，不能有

2、遗漏。（ 3 计算：运行程序，得到输出结果文件。（ 4）绘制内力图：根据结果绘制内力图。（以上最好重新组织一下语言写一下，当然整个报告可以去罗三那里借一份其他同学做出来的报告复印，但是那样的话，必然每个人的分数都降低，对做出来的同学不公）四、实践心得（ 1） （这些地方随便写点东西）（ 2）（ 3）（ 4） （这点是重点）pf 程序对于非节点荷载只能够处理四种情况：水平集中荷载、 竖直集中荷载、水平均布荷载和竖直均布荷载。故而在非节点处若有弯矩作用，需要在此处设定一个节点。（下面是各题的计算过程，这里重新给结构编了号，运行结果是几乎一样的。如果直接copy 下面的结果，那么可能会与别人重复，自

3、己编个再算一遍为好）工力 09 级结构力学矩阵位移法上机题2011 年 5 月（要求：根据电算结果，绘制出内力图，并附上结点位移和杆端力电算结果，该成绩作为最终考试成绩的一部分。）1.作图示刚架的Fn、Fs、M图，已知各杆截面均为矩形，柱截面宽0.5m,高 0.5m,梁截面宽 0.5m,高 0.6m,各杆 E=x 104 MPa 4 分25kN-mW：50kNlOkN/m外OkN!150kNIjOxNlOkN/rr*尸*5m5m5m -2-4对结点和单元进行编号，选定整体坐标系和局部坐标系（以下都是从左到右: 上面节点编号顺序，上边四个：1,2,3,4 中间五个：包括左边梁中非节点的弯矩作用处

4、，5,6,7,8,9F 面四个:10,11,12,13 杆件编号顺序 <号码是带圈字符 >:上边横向三根: 第二层竖直四根: 第三层横向四根，左边那根分为两根: 最下竖直四根: 则对照下边输入数据文件表，大方框中为需要在输入文件中的数据:txt文本为:* 1 a simple suppoted frame* *13 12 130e-230e-230e-225e-225e-225e-225e-230e-230e-230e-230e-225e-225e-225e-2365e5 1512233415273849566778895107118129130 1610 1620 1630 16

5、0 85 810 820 830 80 010 020 030 0101 0102 0103 0111 0112 0113 0121 0122 0123 0131 0132 025e-2133 0 42 0 0 254 150 0 06 0 -150 -507 0 0 -2572 4 -10 103 2 -250 54 3 20 88 4 -5 510 2 -150 511 4 -10 1012 3 20 8运行方法如下：打开windows的dos运行界面，就是“附件”里面的“运行”先把和上面这个名为 zs的输入文件放在c盘下（哪个盘其实都可以），如果将pf运行结果存入 c盘的文件中，则先键入

6、（回车）接着键入输入文件名然后c盘中就出现了名为1的结果记事本文件，输出结果在下面；Dos运行方法截图如下：C； Misers Xnichael>c ； xpf >>c ； 1, txtC ： 2S . tMt输入文件数据表弹性模量 杆件数 结点数 已知位移数荷载情况数365e51513121杆件始端结点号末端结点号截面面积面积矩11230e-222330e-233430e-241525e-252725e-263825e-274925e-285630e-296730e-2107830e-2118930e-21251025e-21371125e-21481225e-215913

7、25e-2结点号X坐标Y坐标10162101632016430165086587108820893081000111001220013300已知位移收位移10101010201030111110112011301212101220123013131013201330（结点他载4LJFXFY FM2002541500060-150-50700-25非结点荷载7LMLTVF DST24-101032-25054320884-55102-1505114-1010123208、输出结果文件*1 a simple suppoted frame*The Input DataThe General Info

8、rmationE NM NJ NS NLC+071513121123456789101112131415joint12345678910111213IS101102103111112The Information of Members member start end1 22 33 41 52 73 84 95 66 77 8895 107 118 129 13The Joint CoordinatesXY.000000.000000.000000.000000.000000.000000.000000The Information of SupportsVS.000000.000000.00

9、0000.000000.000000113121122123131132133ILJ2467ILM2348101112joint123.000000.000000.000000.000000.000000.000000.000000( NA= 438 )( NW= 1319 )Loading Case 1The Loadings at JointsNLJ= 4PXPYPM.0000.0000.0000.00000.0000.0000.0000The Loadings at MembersNLM= 7ITL PV DST4234243The Results of CalculationThe Joint Displacementsphi45678910111213The Terminal Forces member N(st) Q(st) M(st) N(en) Q(en) M(en)123456789101112