原创力文档- 1、本文档共5页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
有限元编程作业
湖南大学有限元编程作业
姓名：邹同球 班级：土硕
题目：编制可以分别考虑和不考虑梁柱轴向变形时，平面框架结构的节点定位向量自动生成程序。
要求：加输入参数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
文档评论（0）