[平面刚架静力分析程序.pptVIP

平面刚架静力分析程序 刘云贺 2009年5月22日 主要变量名表 NE : 单元数 NJ： 结点数 NTYPE: 单元类型数 NDISP: 结点位移未知量总数 NW: 最大半带宽 NPJ: 结点荷载数 NPF: 非结点荷载数 IND: 非结点荷载数代码（见表3-1） BL: 单元长度 SI: 单元的 sin(a) 值 CO: 单元的 cos (a)值 JE(3,NE): JE(1,M) JE(2,M)分别是M号单元始端 i，末端j结点编号 JE(3,M) 是单元类型号 JN(3,NJ): JN(1,I) ，JN(2,I)， JN(3,I)，先是存放结点I的x方向、y方向平移和转动的约束代码，然后存放I的三个方向位移编号。 TYPEE(3,NTYPE): TYPEE(1,I) ，TYPEE(2,I) ，TYPEE(3,I)分别是第I号类型单元的弹性模量E,截面面积A和惯性矩I。 JC(6): 存放单元定位向量数组 X(NJ)，Y(NJ): 结点坐标数组 PJ(3,NPJ): 结点荷载数组 PJ(1,I) ：第 I个结点荷载作用的结点号 PJ(2,I)：第 I个结点荷载作用方向号， =1，x方向， =2， y方向， =3集中力矩 PJ(3,I)：第 I个结点荷载值，集中力沿坐标轴正方向为正，集中力矩顺时针方向为正。 PF(4,NPF): 非结点荷载数组 PF(1,I): 第I个非结点荷载作用单元号 PF(2,I): 荷载类型IND(见表3-1） PF(3,I): 荷载位置a (见表3-1） PF(4,I): 荷载数值q (见表3-1） KD(6,6): 存放局部坐标下单元刚度矩阵数组 KE(6,6)：存放整体坐标下单元刚度矩阵数组 T(6,6)： 存放单元坐标转换矩阵数组 KB(N,NW)：存放整体刚度矩阵数组 P(N)：结点综合荷载数组，解方程组后存放结点位移 F0(6)：先存放整体坐标下的单元等效结点荷载，后存放局部坐标下的杆端力 D(6)： 整体坐标系中单元杆端位移数组 C PROGRAM OF PLANE RIGID FRAME ANALYSIS PROGRAM PRF CHARACTER INDAT*10,OUTDAT*10 C WRITE(*,100) READ(*,200) INDAT WRITE(*,300) READ(*,400) OUTDAT OPEN (5,FILE=INDAT,STATUS=OLD) OPEN (6,FILE=OUTDAT,STATUS=NEW) C WRITE(6,500) INDAT,OUTDAT 100 FORMAT(4X,What is your input filename?) 200 FORMAT(A10) 300 FORMAT(4X,What is your output filename?) 400 FORMAT(A10) 500 FORMAT(4X,INPUT FILENAME:,A10/ * 4X,OUTPUT FILENAME:,A10) C I/O CONTRAL DATA,NODE,ELEMENT,LOAD C AND CONSTRAINTS MESSAGE CALL INPUT C NUMBER JOINT DISPLACEMENT CALL NJDIS C DETERMINE HALF-BAND WIDTH----NW CALL HBW C CALCULATE ASSEMBLE STIFFNESS MATRIX CALL STIFFS C CALCULATE TOTAL NODAL LOADS CALL LOADS C SOLVE SYSTEM EQUATION CALL SOLVE C CALCULATE INTERNAL FORCES AT ENDS OF ELEMENT CALL INFORCE CLOSE(5) CLOSE(6) END INPUT子程序框图 输入单元数NE，结点数NJ，单元类型数NTYPE 结点号，x坐标、y坐标、约束代码 x y z 单元号，I端，j端编号，单元类型号 单元类型号，弹性模量E、截面面积A、惯性矩I 结点荷载数NPJ，非结点荷载数NPF 结点荷载作用结点，方向号，荷载值 非结点荷载作用单元号，类型，位置a 、数值q SUBROUTINE INPUT COMMON /COM1/ NE, NJ,NTYPE,NDISP,NW COMMON /COM2