哈工大电算实习报告.doc

任务 已完成如下了任务： 增加荷载类型（原来程序中只有集中力和均布力）。 计算指定截面的弯矩和剪力。 计算指定截面的弯矩影响线。 计算指定截面的剪力影响线。 计算指定支座的反力影响线。 计算均布活荷作用下的弯矩包络图。 计算均布活荷作用下的剪力包络图。 二 添加的子程序框图 1．计算指定截面的弯矩流程图 2.计算指定截面的弯矩流程图 3.支座反力计算的流程图 3.1计算左侧荷载引起的支座反力的流程图 3.2计算右侧荷载引起的支座反力的流程图 4.弯矩影响线的流程图 5.剪力影响线流程图 6.支座反力的影响线 7.弯矩包络图的流程图（注：对于弯矩的求解运用叠加原理，对每一跨荷载引起连续梁的内力进行叠加） 8.剪力包络图的流程图(此时所加的荷载与求弯矩包络图时的荷载相同,对于剪力的求解运用叠加原理，对每一跨荷载引起连续梁的内力进行叠加)： 三 程序功能 本程序可以解决任意跨数连续梁在集中力偶，集中力以及任意分布的均布荷载作用下连续梁的内力（包括剪力和弯矩）以及支座反力的求解问题；可以解决不同跨的截面惯性矩不同的问题，但要求连续梁的弹性模量相同以及任意一跨内惯性矩保持不变；也可以计算在单位荷载作用下指定截面的弯矩，剪力以及引起某支座的反力的数据并结合Excel求出对应的弯矩，剪力和支座反力影响线；还可以用于求解在均布荷载作用下各截面最大，最小的剪力与弯矩的数据并结合Excel求出弯矩和剪力包络图。 四 使用说明 4.1建立文本文档用于存储初始数据，文本文档应按顺序输入连续梁以及荷载的具体情况，应分别按以下顺序输入：支承类型，节点荷载个数，非节点荷载个数，弹性模量；杆长,惯性矩；节点荷载大小,对应未知数序号（如没有节点荷载可不用输入）；非结点荷载值,起点位置,终点位置,单元号,荷载类型号（如没有非节点荷载可不用输入）；所求弯矩对应截面所在单元以及在单元中的位置；所求剪力对应截面所在单元以及在单元中的位置；所求支座反力的节点号；弯矩影响线对应截面所在单元以及对应截面在单元中的位置；剪力影响线对应截面所在单元以及对应截面在单元中的位置；支座反力影响线对应的节点号。 具体数据如下： 支承类型1：两端铰接 支承类型2：左端固结，右端铰接 支承类型3：左端铰接，右端固结 支承类型4：两端固结 对于集中荷载（如集中力，集中力偶），起点位置,终点位置输入相同的数据（即作用点的位置），对于均匀分布荷载对应起始点以及终止点的位置。 4.2力的正负号规定： 外荷载：竖向荷载以向上为正，向下为负；弯矩以逆时针为正，顺时针为负。 杆端弯矩以顺时针为正，以逆时针为负。 各截面的内力：弯矩以下侧受拉为正，上侧受拉为负；剪力使隔离体产生顺时针旋转为正，产生逆时针旋转为负。 4.3由于数组定义数组大小的限制，本程序最多可以解决20跨连续梁的结构，对于一般情况都能满足要求。 4.4单位的规定：力的单位：KN 力偶单位：KN*m 长度单位：m 4.5荷载类型的规定：1.集中力 2.集中力偶 3.任意均布荷载 五．源程序如下 连续梁静力计算程序 PROGRAM LXL DIMENSION GC(20),GX(20),PJ(20,2),PF(10,5),DK(2,2),P(45) DIMENSION F0(2),WY(2),F(2),ZK1(45),ZK2(45),GDM(20,2) character*12 dat0,dat1 write(*,*)输入:初始数据文件名 read(*,10)dat0 write(*,*)输入:计算结果数据文件名 read(*,10)dat1 10 format(A12) open(8,file=dat0,status=old) open(9,file=dat1,status=new) write(9,11) 11 format(8X,*************连续梁内力计算*****************) READ(8,*)NE,NZ,NP,NF,E0 WRITE(9,12)NE,NZ,NP,NF,E0 12 FORMAT(1X,单元数=,6x,I3,5x,支承类型=,I3,5x,节点荷载个数=I3,/, 1 1x,非节点荷载个数=,5x,I3,弹性模量=,F12.4) NJ=NE+1 CALL SRSJ(NE,NP,NF,GC,GX,PJ,PF,N,X,NN,XX,K,M,XL,MM,XLL,KK) CALL