资源描述：
《结构力学 电算实习》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、一．完成任务：1.增加荷载类型（原来程序中只有集中力和均布力）。在新程序中增加集中力偶和均布力偶两种工况2.计算指定截面的弯矩和剪力。3.计算指定截面的弯矩影响线。4.计算指定截面的剪力影响线。5.计算指定支座的反力影响线。二．程序框图：调用求影响线子程序调用单元单刚子程序调用非节点荷载的等效荷载子程序开始变量与数组的说明部分调用输入初始数据子程序调用形成总荷载矩阵子程序调用集成整体刚度矩阵子程序调用支撑处理子程序调用求解位移子程序调用输出位移子程序调用计算内力子程序结束调用影响线子程序三程序功能原始程序利用fortran90编写，程序
2、中变量服从I-N规则。程序计算仅仅针对于连续梁结构体系，两端支撑类型为刚接或者铰接（梁中间支座均为半铰）。且默认各单元的抗弯刚度不变，抗拉刚度无穷大，不计剪切变形，弹性模量均相等。程序中的单元数、节点荷载数、非节点荷载数，可以根据需要通过改变程序的数组大小来实现。原程序可以求解梁杆端弯矩、节点位移的功能。可以在梁节点作用弯矩荷载，非节点处仅能作用集中力和均布力荷载。我在源程序的基础上增加了另外的荷载类型——在非节点出作用中弯矩荷载。增加了求解杆端剪力的功能。增加了求解任意指定截面弯矩与剪力的功能。增加了绘制指定截面弯矩影响线、剪力影响线
3、与指定支座反力影响线的功能。五．源程序修改部分的说明以及源程序码!连续梁静力计算程序PROGRAMLXLDIMENSIONGC(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)11FORMAT(1X,'单元数,支承类型,节点荷载,非节点荷载,弹性模量')open(6,file='初始数据.txt',status='old')REA
4、D(6,*)NE,NZ,NP,NF,E0NJ=NE+1!!!节点数比单元数多1!输入初始数据CALLSRSJ(NE,NP,NF,GC,GX,PJ,PF,FL)!形成总荷载矩阵PCALLXCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC)!集成整体刚度矩阵ZKl，ZK2CALLJCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)!支承条件处理CALLZCCL(NZ,NJ,ZK1,ZK2,P)!方程求解--求位移PCALLFCQJ(NJ,ZK1,ZK2,P)!输出位移CALLSCWY(NJ,P)!计算内力CALLGDL(N
5、E,NJ,NF,E0,P,WY,F,F0,DK,PF,GC,GX,Q,FJ,WX,FX,FL)!求影响线CALLYXX(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子程序!************************************SUBROUTINESRSJ(NE,NP,NF,GC,GX,PJ,PF,FL)DIMENS
6、iONGC(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)11FORMAT(1X,'杆长，惯性矩GC(NE),GX(NE)')READ(6,*)(GC(I),GX(I),I=1,NE)IF(NP.GT.0)THENWRITE(*,22)22FORMAT(1X,'荷载
7、大小，对应未知数序号PJ(I,1)，PJ(I,2)')READ(6,*)(PJ(I,1),PJ(I,2),I=1,NP)ENDIFIF(NF.GT.0)THENWRITE(*,33)33FORMAT(1X,'非结点荷载值，距离，单元号，荷载类型号')READ(6,*)((PF(I,J),J=1,4),I=1,NF)ENDIF!所求截面信息:FL(?,1)=所求截面的单元号FL(?,2)=所求截面距杆端距离WRITE(*,44)44FORMAT(1X,'所求截面的单元号，距杆端距离，输入完所求截面信息后，输入0,0结束')I=1MO=1D
8、OWHILE(MO>0)READ(6,*)FL(I,1),FL(I,2)MO=FL(I,1)!输入完所求截面信息后，输入0,0结束循环。I=I+1ENDDOEND!*********************
