杆系结构示例.doc

一源程序中主要标识符的意义说明如下：（1）、整型变量NJ结点总数；NM单元总数；N结点位移未知量总数，即整个结构的总自由度数，亦即总刚矩阵的阶数；MB最大的半带宽；NME材料种类数；NGE截面几何特性种类数；NS受约束的结点位移数；NSJ支承结点数；NJP结点荷载个数；JN荷载作用的结点号；JD荷载作用的方向号（1水平力；2垂直力；3集中力偶）；NLM非结点荷载个数；ILT非结点荷载类型号；K单元序号；（2）、整型数组ISE（NM，2）各单元首端和未端结点号；IMG（NM，2）各单元的材料特性号和截面几何特性号；IS（NS）各约束自由度的总体编号；（3）、实型数组WE（NME，2）各种材料特性的容重和弹性模量；AI（NGE，2）各种截面几何特性的截面积和抗弯惯性矩；XY（NJ，2）各结点坐标；SD（NS）各约束自由度的约束位移值（已知的支座位移）；P（N）结构总体荷载列阵，总刚求解后为结点位移增量；TF（NM，6）各单元杆端力增量；BP（NS）各约束点的约束反力增量；DT（N）结点位移累计量；TFT（NM，6）各单元杆端力累计量；二计算示例：1 求图示刚架的内力。设各杆为矩形截面，横梁为，立柱为，立柱的均布荷载集度为，刚架跨度为12m，高度为6m，设弹性模量E=1。(图见结构力学教程(龙驭球、包世华编)P340) 1) 建立数据文件（名字为E1.dat）如下：3,4,6,2,1, NM,NJ,NS,NGE,NME1,1,2,1,1, M,(ISE(I,J),J=1,2),(IMG(I,J),J=1,2)2,2,3,1,23,3,4,1,11, 0.0, 0.0, M,(XY(I,J),J=1,2),I=1,NJ2, 0.0, 6.03,12.0, 6.04,12.0, 0.01,22.0,1.0, M,WE(I,J),J=1,2),I=1,NME1, 0.5, 4.167E-2, M,(AI(I,J),J=1,2),I=1,NGE2, 0.63,8.333E-22, NSJ1,0.0,0.0,0.0, ISJ,(SXYM(I),I=1,3)4,0.0,0.0,0.00,NJP1,NLM1,3,1.0,1.0, K,ILT,PV,DX2) 计算结果 =OUTPUT= NM= 3 NJ= 4 NS= 6 NGE= 2 NME= 1 IE IST IEN IME IGE 1 1 2 1 1 2 2 3 1 2 3 3 4 1 1 IJ X Y 1 .000 .000 2 .000 6.000 3 12.000 6.000 4 12.000 .000 NUM GRAV E 1 .2200E+02 .1000E+01 NUM A I 1 .50000E+00 .41670E-01 2 .63000E+00 .83330E-01 NUM ID 1 1 2 2 3 3 4 10 5 11 6 12 NJP= 0 NLM= 1 K ILT PV DX 1 3 .1000E+01 1.000 THE JOINT DISPLACEMENTS: JOINT u v q 1 .1764E-19 .4277E-20 -.5488E-19 2 .8471E+03 .5133E+01 -.2841E+02 3 .8235E+03 -.5133E+01 -.9648E+02 4 .1236E-19 -.4277E-20 -.4379E-19 THE TERMINAL FORCES: MEMBER N(ST) Q(ST) M(ST) N(EN) Q(EN) M(EN) 1 -.4277E+00 .4764E+01 .8488E+01 .4277E+00 .1236E+01 .2094E+01 2 .1236E+01 -.4277E+00 -.2094E+01 -.1236E+01 .4277E+00 -.3039E+01 3 .4277E+00 .1236E+01 .3039E+01 -.4277E+00 -.1236E+01 .4379E+01三练习1 试利用程序计算图示刚架的结点变形和杆端力，并作出轴力、剪力和弯矩图。