1、一 源程序中主要标识符的意义说明如下：（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
2、）各种材料特性的容重和弹性模量；AI（NGE，2）各种截面几何特性的截面积和抗弯惯性矩；XY（NJ，2）各结点坐标；SD（NS）各约束自由度的约束位移值（已知的支座位移） ；P（N）结构总体荷载列阵，总刚求解后为结点位移增量；TF（NM，6）各单元杆端力增量；BP（NS）各约束点的约束反力增量；DT（N）结点位移累计量；TFT（NM，6）各单元杆端力累计量；二 计算示例：1 求图示刚架的内力。设各杆为矩形截面，横梁为 ，立柱为mhb26.15.02，立柱的均布荷载集度为 ，刚架跨度为 12m，高度为mhb5.0 mkNq/16m，设弹性模量 E=1。(图见 结构力学教程( 龙驭球、包世华编)P
3、340) 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
4、,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= 1IE IST IEN IME IGE1 1 2 1 12 2 3 1 23 3 4 1 112kN/21634IJ X Y1 .000 .0002 .000 6.0003 12.000 6.0004 12.000 .000NUM GRAV E1 .2200E+02 .1000E+01NUM A I1 .50000E+00 .41670E-012 .63000E+
5、00 .83330E-01NUM ID1 12 23 34 105 116 12NJP= 0NLM= 1K ILT PV DX1 3 .1000E+01 1.000THE JOINT DISPLACEMENTS:JOINT u v q1 .1764E-19 .4277E-20 -.5488E-192 .8471E+03 .5133E+01 -.2841E+023 .8235E+03 -.5133E+01 -.9648E+024 .1236E-19 -.4277E-20 -.4379E-19THE TERMINAL FORCES:MEMBER N(ST) Q(ST) M(ST) N(EN) Q(
6、EN) M(EN)1 -.4277E+00 .4764E+01 .8488E+01 .4277E+00 .1236E+01 .2094E+012 .1236E+01 -.4277E+00 -.2094E+01 -.1236E+01 .4277E+00 -.3039E+013 .4277E+00 .1236E+01 .3039E+01 -.4277E+00 -.1236E+01 .4379E+012431.8kN/m23三 练习 1 试利用程序计算图示刚架的结点变形和杆端力，并作出轴力、剪力和弯矩图。设各杆几何尺寸相同， ， ， ， 。要求：ml52.0A41I27/103kNE(1) 写出基本数据输入文件 LX1.dat；(2) 运行上述执行程序，并得到输出结果文件 LX1.out；(3) 根据输出结果文件绘出轴力、剪力和弯矩图。
