结构力学 电算实习.docVIP

一．完成任务： 增加荷载类型（原来程序中只有集中力和均布力）。 在新程序中增加集中力偶和均布力偶两种工况 计算指定截面的弯矩和剪力。 计算指定截面的弯矩影响线。 计算指定截面的剪力影响线。 计算指定支座的反力影响线。 二．程序框图： 调用求影响线子程序 调用单元单刚子程序 调用非节点荷载的等效荷载子程序 开始 变量与数组的说明部分 调用输入初始数据子程序 调用形成总荷载矩阵子程序 调用集成整体刚度矩阵子程序 调用支撑处理子程序 调用求解位移子程序 调用输出位移子程序 调用计算内力子程序 结束 调用影响线子程序 三 程序功能 原始程序利用fortran90编写，程序中变量服从I-N规则。程序计算仅仅针对于连续梁结构体系，两端支撑类型为刚接或者铰接（梁中间支座均为半铰）。且默认各单元的抗弯刚度不变，抗拉刚度无穷大，不计剪切变形，弹性模量均相等。程序中的单元数、节点荷载数、非节点荷载数，可以根据需要通过改变程序的数组大小来实现。原程序可以求解梁杆端弯矩、节点位移的功能。可以在梁节点作用弯矩荷载，非节点处仅能作用集中力和均布力荷载。我在源程序的基础上增加了另外的荷载类型——在非节点出作用中弯矩荷载。增加了求解杆端剪力的功能。增加了求解任意指定截面弯矩与剪力的功能。增加了绘制指定截面弯矩影响线、剪力影响线与指定支座反力影响线的功能。 五．源程序修改部分的说明以及源程序码 ! 连续梁静力计算程序 PROGRAM LXL DIMENSION GC(20),GX(20),PJ(20,2),PF(10,4),DK(2,2),P(45),F0(2),WY(2),F(2,100),ZK1(45),ZK2(45),Q(2),FJ(2,100),FL(100,2),FX(45),WX(45) WRITE(*,11) 11 FORMAT(1X,单元数,支承类型,节点荷载,非节点荷载,弹性模量) open(6,file=初始数据.txt,status=old) READ(6,*)NE,NZ,NP,NF,E0 NJ=NE+1 !!!节点数比单元数多1 ! 输入初始数据 CALL SRSJ(NE,NP,NF,GC,GX,PJ,PF,FL) ! 形成总荷载矩阵P CALL XCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC) ! 集成整体刚度矩阵ZKl，ZK2 CALL JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX) ! 支承条件处理 CALL ZCCL(NZ,NJ,ZK1,ZK2,P) ! 方程求解--求位移P CALL FCQJ(NJ,ZK1,ZK2,P) ! 输出位移 CALL SCWY(NJ,P) ! 计算内力 CALL GDL(NE,NJ,NF,E0,P,WY,F,F0,DK,PF,GC,GX,Q,FJ,WX,FX,FL) ! 求影响线 CALL YXX(NE,NJ,NF,E0,P,WY,F,F0,DK,PF,GC,GX,Q,FJ,WX,FX,FL,ZK1,ZK2,NZ,NP) close(6) ! 计算结束 END ! ************************************ ! 输入初始数据SRSJ子程序 ! ************************************ SUBROUTINE SRSJ(NE,NP,NF,GC,GX,PJ,PF,FL) DIMENSiON GC(NE),GX(NE),PJ(NP,2),PF(NF,4),FL(100,2) ! GC(NE)：杆长 GX(NE)：惯性矩 ! 结点荷载:PJ(NP,1)=荷载大小 PJ(NP,2)=对应未知数序号 ! 非结点荷载 : PF(NE,1)=荷载大小 PF(NE,2)=距离 PF(NE,3)=单元号 PF(NE,4)=荷载类型号 WRITE(*,11) 11 FORMAT(1X,杆长，惯性矩 GC(NE),GX(NE)) READ(6,*)(GC(I),GX(