计算结构力学程序设计报告.docVIP

计算结构力学程序设计报告 学 院 土木工程学院 专业班级 09城建土木x班 学 号 xxxxxxxxxxx 学生姓名 xxx 指导老师 苏少卿 提交时间 2012年06月 实验一: 形成单元定位向量的程序设计 一 求图示结构各单元的单元定位向量。 二 输入数据文件 7 5 5 3 1 4 1 4 2 5 2 1 2 3 4 4 5 1 0 0 10001 2 0 0 10001 3 1 1 10001 4 0 0 10001 5 0 1 10001 三 程序编写过程 ! 输入数据说明(数据之间建议用空格作为分隔符) ! 同一项内容的数据可以在同一行输入,也可分多行输入,建议分多行输入,便于阅读 ! 不同的数据项必须分行输入。 ! 1.单元数NE、结点数NJ、特殊结点数NJT、材料特性组数NAI ! 2.杆端结点号码JH ! 3.特殊结点信息JTX ! 4.杆件材料分组信息JMH ! 5.材料EA ! 6.材料EI ! 7.结点X坐标 ! 8.结点Y坐标 ! 允许指定多重从属关系 MODULE JGLX INTEGER NE,NJ,NJT,NAI,N,MW(6) ! 定义可分配大小的数组 INTEGER,ALLOCATABLE::JH(:,:),JW(:,:),JMH(:),JTX(:,:) DOUBLE PRECISION,ALLOCATABLE::EA(:),EI(:),X(:),Y(:),SL(:),CX(:),SY(:) DOUBLE PRECISION,ALLOCATABLE::ZK(:,:),XSA(:,:) END MODULE JGLX ! 计算单元刚度矩阵,装配结构刚度矩阵 USE JGLX CHARACTER*20 INAME,ONAME INTEGER*2 ID CALL GETARG(1,INAME,ID) IF(ID 0) THEN WRITE(*,(1X,A\)) 请指定输入数据文件： READ(*,(A)) INAME ENDIF OPEN(1, FILE = INAME,STATUS=OLD,ERR=2) ! 打开输入数据文件 CALL GETARG(2,ONAME,ID) IF(ID 0) THEN WRITE(*,(1X,A\)) 请指定输出数据文件： READ(*,(A)) ONAME ENDIF OPEN(2, FILE = ONAME,STATUS=REPLACE,ERR=2) ! 打开输出数据文件 ! 读入单元数、结点数、特殊结点数、材料特性组数 READ(1,*) NE,NJ,NJT,NAI ! 分配数组内存 ALLOCATE(JH(2,NE),JW(3,NJ),JTX(4,NJT),JMH(NE)) ALLOCATE(EA(NAI),EI(NAI),X(NJ),Y(NJ),SL(NE),CX(NE),SY(NE),XSA(NE,7)) READ(1, *) ((JH(I,J),I=1,2),J= 1,NE) ! 读入杆端结点号码 READ(1, *) ((JTX(l,J),l=1,4),J = 1,NJT) ! 读入特殊结点信息 READ(1,*) (JMH(I),I=1,NE) ! 读入杆件材料分组信息 READ(1,*) (EA(I),I=1,NAI) ! 读入材料EA READ(1,*) (EI(I),I=1,NAI) ! 读入材料EI READ(1,*) (X(I),I=1,NJ) ! 读入结点X坐标 READ(1,*) (Y(I),I=1,NJ) ! 读入结点Y坐标 CALL QJW ! 各结点求知量在结构位移列阵中的编号JW数组 WRITE(2,1) NE,NJ,NJT,NAI WRITE(2,(/1X,单元两端结点号:/4(4X,单元, 2X,起点,2X,终点))) WRITE(2,(4(4X,I4,I6,I6))) (J,(JH(I,J),I=1,2),J=1,NE) WRITE(2,(/1X,单元截面特性:/20(1X,I3))) JMH WRITE(2,(/1X,特殊结点信息