有限元程序设计5.docVIP

有 限 元 设 计 报 告 一. 程序功能说明 子程序说明 本子程序包含5个子程序，其主要功能是： 子程序CXCY（E，L，IE，JE）的功能 用于计算各单元的杆件的长度及在结构坐标系中的倾角α的三角函数cosα、sinα。 子程序ELEA（E，L，IE，JE，A）的功能 用于计算结构坐标系下各单元刚度矩阵中的刚度系数。 子程序TOTAL（E，LNODE，KE，KZ）的功能 用于形成半带宽的结构原始刚度矩阵。 子程序FIFX（NP，FXI，FXJ，FYI，FYJ，MMI，MMJ，WL，IE，JE）的功能 它是用来计算各单元的等效节点荷载。 子程序STIFF（E，KE）的功能 用于形成结构坐标系下的单元刚度矩阵 。 二．程序框图设计 三．程序标识符及数组说明 1、数组说明 COORD 节点坐标数组 LNODE 单元节点数组 LL 单元各杆件的长度 PP 存放杆上的荷载 BH（100，2）存放单元杆件的截面宽度（面积）和截面高度（惯性矩） RES(20，2) 存放约束的位移 JP（200，2）存放节点载荷 PP（200，4）存放杆件载荷 KE（6，6） 存放单元刚度矩阵 KZ（200，20）存放半带宽存贮的结构原始刚度矩阵 U（60）、V（60）、PHI（60）分别用于存放节点的位移 变量说明 NPOIN 节点总数 NELEM 单元总数 NJP 节点荷载总数 NPP 杆上荷载总数 NRES 最大位移约束个数 YOUNG 杨氏模量 WT 自重 NP3 位移总数 LL 单元杆长 INF 特征量 当INF＜0（整数） 输入杆件截面宽度B和截面高度H,＞0 当INF≤0（整数） 输入杆件截面面积A和截面的惯性矩I NHBW 最大的半带宽宽度 四．源程序 c 平面刚架内力计算程序 INTEGER E,NELEM,EP,S,DH,ZL,DL,H REAL LL,NI,NJ,MI,MJ,MMI,MMJ,K,JP,KE,KZ COMMON /C1/NPOIN,NELEM,YOUNG/C2/LL,CX,CY COMMON /C3/COORD(60,2),LNODE(100,2)/C4/BH(100,2) COMMON /C5/NPP,PP(200,4)/C6/NJP,JP(200,2) COMMON /C7/NRES,RES(20,2) DIMENSION WL(6),A(6),U(60),V(60),PHI(60), P(180),KE(6,6),KZ(200,20) OPEN(2,FILE=D:/N2.DAT,STATUS=NEW) READ(1,*) NPOIN,NELEM,NJP,NPP,NRES,INF,YOUNG,WT DO 20 I=1,NPOIN 20 READ(1,*) (COORD(I,J),J=1,2) DO 25 I=1,NELEM 25 READ(1,*) (LNODE(I,J),J=1,2) DO 30 I=1,NELEM 30 READ(1,*) (BH(I,J),J=1,2) DO 35 K1=1,NRES 35 READ(1,*) (RES(K1,K2),K2=1,2) IF(INF.LT.0) THEN DO 100 E=1,EO B=BH(E,1) H=BH(E,2) BH(E,1)=B*H 100 BH(E,2)=B*H**3/12 ENDIF C ***************计算半带宽******************* NP3=3*NPOIN NHBW=0 DO 132 E=1,NELEM M=ABS(LNODE(E,1)-LNODE(E,2)) IF(NHBW.LT.M) NHBW=M 132 CONTINUE NHBW=3*(NHBW+1)