有限元编程作业.doc

湖南大学有限元编程作业 姓名：邹同球 班级：土硕 题目：编制可以分别考虑和不考虑梁柱轴向变形时，平面框架结构的节点定位向量自动生成程序。 要求：加输入参数NA NA=1，计梁柱轴向变形 NA=2，不计柱的轴向变形 NA=3，不计梁的轴向变形 NA=4，不计梁柱轴向变形 1 编制Fortran语言程序如下： INTEGER NTJ(90,4) COMMON/C17/JDW(150,3) CALL JDDWXL(NZ,N,NCJ,NTJ) print *,节点定位向量:节点号，x向位移，y向位移，转角向位移编号如下: WRITE(*,(4I10))(I,(JDW(I,J),J=1,3),I=1,NZ) END SUBROUTINE JDDWXL(NZ,N,NCJ,NTJ) COMMON/C17/JDW(150,3) DIMENSION NTJ(90,4) DIMENSION G(90,2) INTEGER I,J,NZ,N,NCJ,NA PRINT *,请输入轴向变形情况:1,都考虑;2,梁考虑柱不考虑;3,都不考虑;4,柱考虑粱不考虑. READ(*,*)NA PRINT *,请输入总节点数，可动节点数，从节点数: READ(*,*)NZ,N,NCJ !总节点数，可动节点数，从节点数 IF(NCJ.NE.0)THEN PRINT *,请输入从节点信息: READ(*,*)((NTJ(I,J),J=1,4),I=1,NCJ) END IF P=1 IF(NA.EQ.2)THEN P=0 ELSE IF(NA.EQ.3)THEN P=0 END IF IF(NA.GE.3)THEN PRINT *,请输入跟节点数目(跟节点即高度相同的节点数减去高度种类和从节点数): READ(*,*)X PRINT *,请输入跟节点信息: READ(*,*)((G(I,J),J=1,2),I=1,X) END IF DO 10 I=1,NZ DO 10 J=1,3 10 JDW(I,J)=0 IF(NCJ.NE.0)THEN DO 20 I=1,NCJ DO 20 J=1,3 20 JDW(NTJ(I,1),J)=NTJ(I,J+1) END IF IF(NA.GE.3)THEN IF(X.NE.0)THEN DO 80 I=1,X 80 JDW(G(I,1),1)=G(I,2) END IF END IF K=0 IF(P.EQ.0)THEN DO 40 I=1,N DO 40 J=1,3,2 IF(JDW(I,J))50,60,70 50 JDW(I,J)=0 GOTO 40 60 K=K+1 JDW(I,J)=K GOTO 40 70 JDW(I,J)=JDW(JDW(I,J),J) 40 CONTINUE ELSE IF(P.EQ.1)THEN DO 90 I=1,N DO 90 J=1,3 IF(JDW(I,J))100,110,120 100 JDW(I,J)=0 GOTO 90 110 K=K+1 JDW(I,J)=K GOTO 90 120 JDW(I,J)=JDW(JDW(I,J),J) 90 CONTINUE END IF IF(P.EQ.0)THEN DO I=1,N JDW(I,2)=0 END DO END IF END 2 给定框架调试，框架结构及节点编号如下图所示： 3 调试结果： 请输入轴向变形情况:1,都考虑;2,梁考虑柱不考虑;3,都不考虑;4,柱考虑粱不考虑. 1 请输入总节点数，可动节点数，从节点数: 17,14,7 请输入从节点信息: 2,1,1,0 4,0,-1,0 6,5,5,0 8,7,7,0 9,7,7,0 10,0,-1,0 14,13,13,0 节点定位向量:节点号，x向位移，y向位移，转角向位移编号如下: 1 1 2 3 2 1 2 4 3 5 6 7 4 8 0 9 5 10 11 12 6 10 11 13 7 14 15 16 8 14 15 17 9 14 15 18 10 19 0 20 11 21 22