结构力学电算作业86
结构力学
1．杆件结构——这类结构是由杆件所组成。杆件的几何特征是横截面尺寸要比长度小得多。梁、拱、桁架、刚架是杆件结构的典型形式。
当我们拿到了一个结构对它进行几何组成分析时，我们一般可以遵循这样的步骤来进行：
首先，我们可以试着利用结构几何组成的四个规律（即三个规则或称为：“三角形规律”）来对这个结构进行分析，1、三根支杆，只看本身，四根支杆，基础算刚片；2、尽可能去掉二元体，扩大刚片，等效替换，简化结构；3、利用“两刚片规则”；4、利用“三刚片规则”（顺藤摸瓜法）；
其次，从结构的概念——建筑物和工程设施中承受、传递荷载而其骨架作用的部分成为工程结构、简称结构。我们可以看出结构在民用和军用设施中都是无处不在的，我院作为全军工程兵的最高学府，工程兵的八大专业：桥梁、道路、渡河、地雷、爆破、筑城、伪装、机械，有四大专业：桥梁、道路、爆破、筑城是直接构筑在结构力学之上的，其他的四个专业也都与结构力学密不可分。
（3） 三个虚铰在无限远处
三刚片分别用三对任意方向的平行链杆相联结，均为瞬变体系。
对于一般按照三角形规律组成的体系，以上的步骤就可以奏效，但我们知道，满足三角形规律只是结构几何不变的必要条件，即结构还可以按其它组成规律组成几何不变体系，上述方法如若不行；第二步，计算一下结构的计算自由度W。引入计算自由度W，可以根据W得到一些关于自由度S和多余约束n的定性结论。
作题练习，是学习结构力学的重要环节。不作一定数量的习题，就很难对基本概念和方法又深入的理解，也很难培养好的计算能力。但是作题也要避免各种盲目性。举例如下：
框架结构内力位移计算算例
结构力学课程大作业——多层多跨框架结构内力及位移计算班级学号姓名华中科技大学土木工程与力学学院年月日结构力学课程大作业——多层多跨框架结构内力与位移计算一、任务1、计算多层多跨框架结构在荷载作用下的弯矩和结点位移。
2、计算方法要求：（1）用迭代法、D 值法、反弯点法及求解器计算框架结构在水平荷载作用下的弯矩，并用迭代法的结果计算其结点位移。
（2）用迭代法、分层法、二次力矩分配法及求解器计算框架结构在竖向荷载作用下的弯矩，并用迭代法的结果计算其结点位移。
3、分析近似法产生误差的原因。
二、计算简图及基本数据本组计算的结构其计算简图如图1所示，基本数据如下。
混凝土弹性模量：723.010/h E kN m =⨯构件尺寸：柱：底 层：23040b h cm ⨯=⨯其它层：23030b h cm ⨯=⨯ 梁：边 梁：22560b h cm ⨯=⨯中间梁：22530b h cm ⨯=⨯ 水平荷载：'15P F kN ＝，30P F kN ＝（见图2）竖向均布恒载：17/q kN m 顶＝ 21/q kN m 其它＝（见图8） 图1各构件的线刚度：EIi L =，其中312b h I ⨯=边 梁：33410.250.6 4.51012I m -⨯==⨯F 7311 3.010 4.510225006EI i kN m L -⨯⨯⨯===⋅ 中间梁： 34420.250.3 5.6251012I m -⨯==⨯ 7422 3.010 5.6251067502.5EI i kN m L -⨯⨯⨯===⋅ 底层柱： 33440.30.4 1.61012I m -⨯==⨯ 7344 3.010 1.61096005EI i kN m L -⨯⨯⨯===⋅ 其它层柱：34430.30.3 6.751012I m -⨯==⨯ 7433 3.010 6.75106136.43.3EI i kN m L -⨯⨯⨯===⋅ 三、水平荷载作用下的计算 （一）用迭代法计算1、计算各杆的转角分配系数ikμ' 转角分配系数计算公式：()2ikikiki i i μ'=-∑结点“1”：12225000.3932(6136.422500)μ'=-=-⨯+156136.40.1072(6136.422500)μ'=-=-⨯+结点“2”：21225000.3182(67506136.422500)μ'=-=-⨯++图2232(67506136.422500)⨯++266136.40.0872(67506136.422500)μ'=-=-⨯++由于该结构是对称结构，因此结点“3”的分配系数应该等于结点“2”的，结点“4”的分配系数应该与结点“1”的相等，所以本题只需计算1、2、5、6、9、10、13、14、17、18结点的分配系数。
110～220kV同塔多回输电线路耐张塔带电检修方法研究
110～220kV同塔多回输电线路耐张塔带电检修方法研究隆晨海;邹德华;易子琦;李稳;李金亮【摘要】为确保同塔多回输电线路的供电可靠性,开展同塔多回线路带电作业研究成为当前研究的热点.本文针对目前110～220kV同塔多回输电线路耐张塔带电作业存在的技术难点,对典型110～220kV同塔多回输电线路杆塔进行结构力学分析,提出了倒序操作法更换耐张串绝缘子,研制了配套工器具,有效地解决了110～220kV同塔多回输电线路耐张塔带电检修面临的作业空间小、安全裕度低等难题,推动带电作业的发展.【期刊名称】《电气技术》【年(卷),期】2016(000)008【总页数】5页(P77-81)【关键词】输电线路;同塔多回;耐张塔;倒序操作法;带电检修【作者】隆晨海;邹德华;易子琦;李稳;李金亮【作者单位】国网湖南省电力公司带电作业中心,长沙 410007;国网湖南省电力公司带电作业中心,长沙 410007;国网湖南省电力公司带电作业中心,长沙 410007;国网湖南省电力公司带电作业中心,长沙 410007;国网湖南省电力公司带电作业中心,长沙 410007【正文语种】中文随着国民经济的发展，我国电力需求日益提高，电网改造进一步加快，输电线路走廊日趋紧张。
鉴于此，近年来同塔多回输电线路被广泛应用与输电线路建设，它可以在提高单位面积输送容量、满足电力输送需求的同时，减少了综合造价和线路走廊需求，具有明显的经济效益和社会效益。
同塔多回线路成为输电线路建设的必然趋势。
相对于常规的单回交流输电线路，同塔多回输电线路普遍存在杆塔高、杆塔型式多样、横担间距离和相间净空距离小、导体布置相对复杂等特征，在城市中间的同塔多回输电线路这一特征尤为明显，从而也导致了电场畸变严重、安全距离裕度小等一系列问题。
目前已有的带电检修作业研究成果已不能满足同塔多回输电线路带电作业发展的需求，因此，更深入的研究同塔多回输电线路带电作业技术具有重要意义。
《结构力学》第9章矩阵位移法.
对结构整体建立坐标系oxyz，则每个结点都有确定的位置坐标。
下标I表示结点编号，上标T表示矩阵转置。
结构力学
对结构所有的结点位移，统一用矢量Δ表示，称为结构整体位 移，简称结构位移或整体位移。Δ中各分量的顺序首先是结点 编号，然后是每个点本身的x,y,z顺序，即
对应结点载荷用矢量F表示，它的排序与位移排序相同
整体坐标系下单元杆端力与杆端位移间的关系—刚度方程： 简写为 其中Ke称为整体坐标系下的单元刚度矩阵。
结构力学
9.4 结构的整体刚度方程和整体刚度矩阵
上式称为结构的整体刚度方程，其中K称为结构的整体刚度 矩阵。
总体刚度矩阵是一个方阵，其阶数与结构结点位移分量总 数相同。它的分量是由单元刚度矩阵的系数叠加构成的。叠加 规律是：单元刚度矩阵的元素，按照它所处的局部行和列号， 对应单元的定位向量，在总刚度矩阵中落到新的行和列上。 总刚度矩阵的特点： (1)刚度矩阵的系数是物理量，由结构本身的长度、截面尺寸、 材料性质、连接方式等决定，与载荷、变形等量无关。 (2)总刚度系数kij表示结构沿第j个整体结点位移方向产生单位 位移Δj=1，其他所有结点位移等于0时，在第i结点位移方向所 需要施加的力(与传统位移法相同)。
结构力学
9.5 非结点荷载的等效化
计算步骤： 1. 在局部坐标系下计算单元的等效载荷 2. 将固端力转换到结构(整体)坐标系 3. 等效结点载荷FP
结构力学
9.6 计算步骤和算例
矩阵位移法的基本步骤如下：
(1)整理原始数据，对结点位移进行整体编码，得到单元定位向量等。 直接的结点载荷按它对应的结点位移编码，直接计入整体结点载荷向量 F中。 (2)单元分析，先形成局部坐标系中的单元刚度矩阵 ，用式(9-10)。 再形成整体坐标系中的单元刚度矩阵Ke，用式(9-24)。 (3)整体分析，依定位向量，将单元刚度矩阵“对号入座”集成总刚度 矩阵K。
结构力学专题习题解答_图文
圆盘转动时的惯性力偶为 平衡方程
其中 利用初始条件得
16-13试求图示梁的自振频率和主振型。梁承重可略去不计 EI=常数
P1=1
解：(1)计算自振频率 分别画出该梁在P1=1,P2=1作用 下的弯矩图M1,M2
P2=1
(2)计算主振型
16-21用振型分解法重作题16-19
解:由于此刚架振动时，各横梁不 能竖向地移动和转动而只能作水 平移动。故只有三个自由度。 (1)按刚度系数如图
(2)确定主振型
由于上式的系数行列式为0。故三个方程中只有两个是独立 的，可有三个方程中任取两个计算得
(3) 求广义质量
(4)广义荷载为
由于荷载为简谐振动, 其正则坐标幅值为
(5)求位移幅值
得 ：
， ，
试求下图楔形悬臂梁的自振频率。设梁的截面宽度b=1,截面高度为
直线变化
。
：
解 截面惯性矩 ：
单位长的质量
设其振型函数为 ：
因
，满足边界条件，
所以
如图所示为一圆轴AB,a端有一圆盘。设圆轴质量远比圆盘小 当t=0时，圆轴受有扭转变形，圆，盘具有初始角位移 和初始速度
，然后体系作自由振动，圆盘在任一时刻t的转角为 ，转动 惯量 ，试出体系自由振动的微分方程及其解答。
该刚架的极限荷载pu=32Mu/5L
θ
机构四
15-9 试用静力法求图示结构的稳定方程及临界荷载
l
l
l
解 ：
• 平衡微分方程为：
• 边界条件为 ：
• 因此得齐次方程为 ：
• 特征方程
16-9 图示悬臂梁具有一重量G=12KN的集中质量，其上受有振动荷载 其中p=5KN。若不考虑阻尼，试分别计算该梁在 振动
结构力学电算实习布置介绍
任务五：写总结，交报告
对整个电算实习过程进行总结，编写实习
报告一份。
基本要求：独立写电算实习的体会、收获、建 议、感想；用A4纸，最少写一页；机写打印或 手写都可以。
《结构力学电算实习报告》装订要求
1、实习报告统一封皮（将班级、姓名按所给空填写完整、打印）
2、《结构力学电算实习》任务书（整洁、完整，打印附上即可）
土木工程（建筑工程方向） 2010 级
结构力学 电算实习
指导教师：武胜、李国东、徐嫚
公共邮箱名： jglxds2010@ 密码： 2010jglxds
1《实习报告》统一格式（1班） 2《实习报告》统一格式（2班） 3《实习报告》统一格式（3班） 4 VB源程序文件 （缺少三个子 程序，需同学自己填充、调试） 5 VB源程序文件说明（.doc 供同学阅读、熟悉、编程使用） 6 VB中文简体安装程序
pf(0 to npf，1 to 4) ——由全部npf个非结点荷载的四种数据 组成（荷载数值、位置参数、作用杆码、荷载类型码）。 当npf=0时，可将第0行的四个元素设为0。 当npf>0时， pf(I，1)——第i号非结点荷载的数值 pf(I，2)——第i号非结点荷载的位置参数 集中荷载为荷载作用点到杆件始端距离 均布荷载为荷载长度（荷载起点规定与杆件始端 重合） pf(I，3)——第i号非结点荷载所在杆件编码。 pf(I，4)——第i号非结点荷载的类型码。 类型码1——均布荷载 类型码2——垂直集中荷载 类型码3——平行集中荷载
平面刚架电算程序数据符号说明
1. 结构基本参数
ne——单元数 nj——结点数
nz——支承数
npj——结点荷载数 npf——非结点荷载数
建筑施工高处作业安全技术规范(JGJ80-91)-推荐下载
三、钢筋横杆上杆直径不应小于 16mm，下杆直径不应小于 14mm，栏杆柱直径不应小于 18mm，采 用电焊或镀锌钢丝绑扎固定。
楼层周边及水箱与水塔周边等处，都必须设置防护栏杆。
二、头层墙高度超过 3.2m 的二层楼面周边，以及无外脚手的高度超过 3.2m 的楼层周边，必须在外围 架设安全平网一道。
三、分层施工的楼梯口和梯段边，必须安装临时护栏。顶层楼梯口应随工程结构进度安装正式防护栏 杆。
四、井架与施工用电梯和脚手架等与建筑物通道的两侧边，必须设防护栏杆。地面通道上部应装设安 全防护棚。双笼井架通道中间，应予分隔封闭。
四、钢管横杆及栏杆柱均采用 Φ48×（2.75～3.5）mm 的管材，以扣件或电焊固定。
第 2.0.8 条 因作业必需，临时拆除或变动安全防护设施时，必须经施工负责人同意，并采取相应的可 靠措施，作业后应立即恢复。
第 2.0.9 条 防护棚搭设与拆除时，应设警戒区，并应派专人监护。严禁上下同时拆除。
第 2.0.10 条 高处作业安全设施的主要受力杆件，力学计算按一般结构力学公式，强度及挠度计算按 现行有关规范进行,但钢受弯构件的强度计算不考虑塑性影响，构造上应符合现行的相应规范的要求。
第 2.0.3 条 高处作业中的安全标志、工具、仪表、电气设施和各种设备，必须在施工前加以检查，确 认其完好，方能投入使用。
第 2.0.4 条 攀登和悬空高处作业人员以及搭设高处作业安全设施的人员，必须经过专业技术培训及专 业考试合格，持证上岗，并必须定期 进行体格检查。
电算报告 (2)
H a r b i n I n s t i t u t e o f T e c h n o l o g y结构力学电算报告院系：土木工程学院班级：设计者：学号：指导教师：设计时间：哈尔滨工业大学一任务要求在读懂连续梁内力计算程序的基础上，完成如下任务：1．增加荷载类型（原来程序中只有集中力和均布力）。
2．计算指定截面的弯矩和剪力。
3．计算指定截面的弯矩影响线。
4．计算指定截面的剪力影响线。
5．计算指定支座的反力影响线。
二程序功能本程序可计算均布力、集中力、均布弯矩、集中弯矩、线性三角形荷载作用下连续梁的截面内力及影响线。
三使用说明程序运行操作步骤：按提示输入杆件与荷载的各类数据。
支承类型(1:两铰2:左刚右铰3:左铰右刚4:两刚) 节点荷载个数,非节点荷个数,弹性模量, 非结点荷载值，距离，单元号，荷载类型号(1:均布力2：集中力3：集中弯矩4：均布弯矩)，即可输出对应单元杆端力。
输入所求截面所在单元及位置，即可输出截面剪力与弯矩，并可绘制剪力与弯矩图。
程序计算过程为：初值的输入，形成总荷载矩阵，集成整体刚度矩阵，支承条件处理，方程求解--求位移输出位移，计算杆端内力，计算指定截面的剪力和弯矩，计算指定截面的剪力和弯矩影响线。
首先是初值的输入，本程序需要将所有要输入的数据依次按单元数、支承类型、节点荷载、非节点荷载、弹性模量、各单元杆长惯性矩、结点荷载大小、对应未知数序号、非结点荷载荷载大小、距离、单元号、荷载类型号、内力计算的指定截面单元号、距杆端距离、影响线指定截面单元号、距杆端距离顺序将数据被输出。
另外，在输入内力计算的指定截面单元号、距离后输入0,0表示结束，在输入影响线指定截面单元号、距杆端距离完成后也一样。
支撑类型有四种：1代表两端铰支；2代表左端固定，右端铰支；3代表右端固定，左端铰支；4代表两端固定。
非荷载类型有三种：1为均布荷载，均布荷载的分布长度；2为集中力；3为集中力偶。
方向规定：外集中力/分布力方向向上为正，剪力为绕截面顺时针方向为正，弯矩逆时针方向为正，转角位移逆时针方向为正。
结构力学电算
1-1题具体解答过程+结点,1,0,0 结点,2,3,0 结点,3,7,0结点,4,10,0单元,1,2,1,1,0,1,1,1 单元,2,3,1,1,1,1,1,1 单元,3,4,1,1,1,1,1,0 结点支承,1,2,-90,0,0 结点支承,2,1,0,0 结点支承,3,1,0,0 结点支承,4,1,0,0 单元荷载,3,-2,8,1单元荷载,2,1,8,1/2,90单元材料性质,1,3,-1,10,0,0,-1无穷大=无穷大=无穷大求解结果：（1）多余约束:2 自由度：0 结论：有多余约束的几何不变体系 （2）内力计算杆端内力值 ( 乘子 = 1) --杆端 1 杆端 2单元码 轴力 剪力 弯矩 轴力 剪力 弯矩 --1 0.00000000 -1.06666667 0.00000000 0.00000000 -1.06666667 -3.200000002 0.00000000 4.60000000 -3.20000000 0.00000000 -3.40000000 -0.800000003 0.00000000 -2.40000000 -0.80000000 0.00000000 -2.40000000 0.00000000 --弯矩图6.00=10EI EI =10EI =10 = EA EA EA 3.00 4.00 3.00-3.40-3.40剪力图（3）位移计算杆端位移值 ( 乘子 = 1)--杆端 1 杆端 2单元码 u -水平位移 v -竖直位移?-转角 u -水平位移 v -竖直位移?-转角--1 0.00000000 0.00000000 0.16000000 0.00000000 0.00000000 -0.320000002 0.00000000 0.00000000 -0.32000000 0.00000000 0.00000000 0.480000003 0.00000000 0.00000000 0.48000000 0.00000000 0.00000000 -0.84000000 2-2题具体解答过程结点,1,0,0结点,2,4,0结点,3,0,4结点,4,4,4结点,5,8,4单元,1,3,1,1,1,1,1,1单元,3,4,1,1,1,1,1,1单元,4,5,1,1,1,1,1,0单元,4,2,1,1,1,1,1,1结点支承,5,2,0,0,0结点支承,1,6,0,0,0,0结点支承,2,6,0,0,0,0结点荷载,3,-2,10结点荷载,4,-2,10单元材料性质,1,4,100,10,0,0,-1单元材料性质,4,4,100,10,0,0,-1单元材料性质,1,1,100,10,0,0,-1单元材料性质,4,4,100,10,0,0,-1单元材料性质,2,3,200,20,0,0,-1求解结果：（1）多余约束:2 自由度：0 结论：有多余约束的几何不变体 （2）内力计算杆端内力值 ( 乘子 = 1) --杆端 1 杆端 2 -单元码 轴力 剪力 弯矩 轴力 剪力 弯矩 --1 3.04780875 -1.14529645 1.45467610 3.04780875 -1.14529645 -3.126509702 -1.14529645 -3.04780875 6.87349030 -1.14529645 -3.04780875 -5.317744713 -1.75013692 -0.75639849 3.02559395 -1.75013692 -0.75639849 0.000000004 -2.29141027 -0.60484047 1.65666134 -2.29141027 -0.60484047 -0.76270054弯矩图1212=10 EI =100EA EI =20 EA =200 EI =20 EA剪力图轴力图（3）位移计算杆端位移值 ( 乘子 = 1)杆端 1 杆端 2单元码 u -水平位移 v -竖直位移 ?-转角 u -水平位移 v -竖直位移 ?-转角 -1 0.00000000 0.00000000 0.00000000 0.05790867 0.12191235 -0.33436672 2 0.05790867 0.12191235 -0.33436672 0.03500274 -0.09165641 -0.178792163 0.03500274 -0.09165641 -0.17879216 0.00000000 0.00000000 0.12376723 4 0.03500274 -0.09165641 -0.17879216 0.00000000 0.00000000 0.000000001 212( 2 )( 3 )。
《结构力学》实验课程结构数值仿真实验DOC
模型 / 材料和截面特性 /
材料
名称 ( Grade3)
设计类型 > 钢材
规范 > GB(S) ; 数据库 > Grade3
模型 / 材料和截面特性 / 截面数据 截面号 ( 1 ) ;
截面 截面形状 > 箱形截面 ;
4
《结构力学》结构仿真实验——初级课程：连续梁例题讲稿 用户：如图输入 ; 名称 > 400× 200× 12
显示 边界条件 >一般支承 ( 开)
模型 / 单元 /
单元的复制和移动
全选
形式 > 复制 ; 移动和复制 > 等间距
dx, dy, dz ( 0, 0, -5 )
; 复制次数 ( 2 )
复制节点属性 ( 开 ), 复制单元属性 ( 开)
模型 1 模型 2 模型 3
图 1.11 复制单元
11
《结构力学》结构仿真实验——初级课程：连续梁例题讲稿
( 开)
图 1.13 均布荷载引起的反力
以表格的形式查看均布荷载引起的的反力。比较外荷载总合和反力的总合来查看模型的建立和荷载 的输入是否恰当。 例题 Z轴方向荷载为 1.0 tonf/m 2×20 m× 3 = 60 tonf ，与表格中 Z轴方向的反力（ FZ）总和相等。
结果 / 分析结果表格 / 反力 荷载组合 > 均布荷载 (ST) ( 开 ) ;
结果 / 内力 /
梁单元内力图
荷载工况 / 荷载组合 > ST: 温度荷载 ; 内力 > My
显示选项 > 精确解 ; 不涂色 ; 放大 ( 2.0 )
显示类型 > 等值线 ( 开), 数值 ( 开)
数值
小数点以下位数 ( 1 ) ; 指数型 ( 关 ) ; 适用于选择确认时 ( 开)
结构力学位移法
拆成基本 结构
结构位移 加上某些条件
原结构的变形协调条件（力法基本方程）
精选可编辑ppt
7
位移法：
先求某些结点位移
解题过程：
结构内力
结构
拆成单根杆件 的组合体
加上某些条件
1.杆端位移协调条件
2.结点的平衡条件
精选可编辑ppt
8
适用范围：
力法: 超静定结构 位移法： 超静定结构，也可用于静定结构。
l
MA FB, MBFA为 杆由 端荷 弯载 矩和 ，温 称度 为变固化端引弯起矩的。
精选可编辑ppt
27
A端固定B端铰支杆的转角位移方程为
A
B
M AB 3i
A3 li
M AB
F AB
精选可编辑ppt
28
A端固定B端定向杆的转角位移方程为
MABiA MAFB MBAiA MBFA
精选可编辑ppt
根据叠加原理，共同作用等于单独作用的叠加：
R1＝R11＋R1P=0
（a)
R11为强制使结点发生转角Z1时所产生的约束反力矩。
R1P为荷载作用下所产生的约束反力矩。
精选可编辑ppt
16
为了将式(a)写成未知量Z1的显式，将R11写为
Z1=1 R11=r11Z1
R1 1r1 1Z1
r11为 单 位 转 角 （ Z1 ＝ 1 ） 产生的约束反力矩。
位移法在解题上比较规范，具有通用性，因 而计算机易于实现。
位移法可分为：手算——位移法
电算——矩阵位移法
精选可编辑ppt
5
力法与位移法最基本的区别：基本未知量不同
力法：以多余未知力基本未知量 位移法：以某些结点位移基本未知量
◆◆合集知识资料知识资料结构力学(一～十)(新版)
2024山东省安全员知识题库附答案（推荐）（单选题、多选题、判断题）一、单选题(101题)1.可能导致重大事故发生的(B )就是重大危险源。
A.隐患B.危险源C.事故隐患D.单元2.纵向水平杆（大横杆）的对接扣件应交错布置，两根相邻杆的接头，在不同步或不同跨的水平方向错开的距离应（ C ）。
A.不小于300mmB.不小于400mmC.不小于500mmD.不大于500mm3.分包单位应按规定建立安全机构并配备专职（C）。
A.材料员B.施工员C.安全员D.质检员4.当连续5天日平均气温低于（A ），一般来讲即进入冬期施工阶段。
A.5℃B.-5℃C.-8℃D.0℃5.（C）是在应急响应过程中，为最大限度地降低事故造成的损失或危害，防止事故扩大，而采取的紧急措施或行动。
A.应急预案B.应急准备C.应急响应D.应急演练6.工程式总承包单位为某建设工程项目中从事危险作业的人员办理意外伤害保险，该保险责任期限到该项目（B）之日为止。
A.提交竣工验收报告B.竣工验收合格C.竣工验收备案D.交付使用7.施工升降机是一种使用工作笼（吊笼）沿（B ）作垂直（或倾斜）运动来运送人员和物料的机械。
A.标准节B.导轨架C.导管D.通道8.根据我国《安全生产法》的规定,安全设施投入( D )。
A.必须纳入建设项目预算B.应当单独进行概算C.必须单独进行概算D.应当纳入建设项目概算9.电焊时严禁借用金属管道、金属脚手架、结构钢筋等金属物搭接代替（A）使用。
A.导线B.支撑物C.接地线D.防护物10.以下不属于施工用电组织设计内容的是（B ）。
A.现场勘探B.进行路线计算C.绘制电气平面图、立面图和接线系统图D.设计防雷装置11.易燃易爆物品应分类储藏在专用库房，并采取（B ）措施。
A.防火B.防雷C.防盗D.防汛12.实施建筑工程监理前,建设单位应当将委托的工程监理单位、监理的内容及监理权限, ( B )被监理的建筑施工企业。
结构力学电算实习报告
结构力学电算实习报告摘要:现实生活中往往存在一些很复杂的连续梁模型，它们的支撑性质和受力特点都比较复杂，结构力学电算则是将它们有规律的进行简化，最后将其所有已知的性质表达为一组数据，将这些数据输入到我们所编好的程序中，既可以得到我们所要求一些量如：杆端力，各截面内力，某界面的内力影响线等。
关键词：连续梁、电算程序、截面内力程序功能：连续梁程序可对连续梁结构进行相关的计算。
能够计算在任意形式均布荷载、跨中集中力荷载、跨中力偶荷载、节点力偶荷载组合下连续梁单元的杆端力，也可计算连续梁各截面内力，从而绘制连续梁的弯矩图。
影响线程序可以计算单位1集中力荷载作用下的指定截面中内力在单位1载荷按指定距离移动情况下的变化情况，从而提取数据绘制影响线。
使用说明：（1）连续梁内力程序在连续梁计算程序（.exe文件）所在文件夹中建立一新的txt文件，并命名（1.txt）。
在文件其中依次输入单元数、支承类型、节点荷载个数、非节点荷载个数、指定截面个数、弹性模量、杆长&惯性矩（依次按组输入）、节点荷载大小&对应节点编号（依次按组输入）、非节点荷载值&距该单元左端距离&单元号&荷载类型号（依次按组输入）、所求截面所在单元号&距该单元左端距离（依次按组输入），然后保存该文件。
运行连续梁内力程序，按提示输入起始文件即初始数据的文件（要求全名），再输入结果文件名，在相应文件夹下打开结果文件即可查看结果。
（2）影响线程序在影响线程序（.exe文件）所在文件夹中建立一新的txt文件，并命名。
在其中一次输入单元数、支承类型、要求影响线的数量、弹性模量、单位载荷每次移动的距离、杆长&惯性矩、截面所在单元号&距该单元左端距离，然后保存该文件。
运行影响线程序，输入起始文件即初始数据的文件，再输入结果文件名，在相应文件夹下打开结果文件即可查看结果。
提取期中的数据，利用excel或其它工具绘制影响线即可。
三款弹性层状体系理论电算程序应用技术评析
总第３２１期交 通 科 技ＳｅｒｉａｌＮｏ．３２１ ２０２３第６期ＴｒａｎｓｐｏｒｔａｔｉｏｎＳｃｉｅｎｃｅ＆ＴｅｃｈｎｏｌｏｇｙＮｏ．６Ｄｅｃ．２０２３ＤＯＩ１０．３９６３／ｊ．ｉｓｓｎ．１６７１ ７５７０．２０２３．０６．０１５收稿日期：２０２３ ０５ ２３第一作者：李政贤（２０００－），男，硕士生。
通信作者：蒋鑫（１９７６－），男，博士，教授。
三款弹性层状体系理论电算程序应用技术评析李政贤１，２，３ 蒋 鑫１，２，３ 沙马伍呷１，２，３ 向嫣然１，２，３ 张 免１，２，３（１．西南交通大学土木工程学院 成都 ６１００３１；２．西南交通大学道路工程四川省重点实验室 成都 ６１００３１；３．西南交通大学高速铁路线路工程教育部重点实验室 成都 ６１００３１）摘 要 基于弹性层状体系理论开发的沥青路面结构电算程序数量众多、应用广泛、各具特色。
文中选取ＥＬＳＹＭ５、ＷＥＳＬＥＡ和ＥＶＥＲＳＴＲＥＳＳ这３款具有代表性的电算程序，从结构层、荷载、计算点等方面出发，讨论它们在计算轮载作用下层状沥青路面结构附加应力时应用技术层面的差异；通过一具体算例，经与解析解对比，提出并验证适用于这３款程序的一种结构层划分方式，进而阐述ＷＥＳＬＥＡ和ＥＶＥＲＳＴＲＥＳＳ程序应用上存在的缺陷，帮助使用者合理选择电算程序，并解决其在应用这２款程序时可能的困扰。
关键词 弹性层状体系理论 电算程序 沥青路面 力学响应中图分类号 Ｕ４１２．６ Ｕ４１６．２１７ 沥青路面具有表面平整、无接缝、行车舒适、耐磨、振动小、噪声低、施工期短、养护维修简便、适宜于分期修建等优点，因而得到了广泛的应用［１］。
目前可应用于沥青路面结构力学分析的主要方法包括弹性层状体系理论、有限单元法（含轴对称有限元法、三维有限元法）、连续有限层法等［２］，这些方法各有特点。
在这些方法中，基于沥青路面各结构层水平层状分布之显著特征，弹性层状体系理论无疑是其中最为成熟、已被广泛采纳的经典理论。
桥梁结构电算--第1讲
桥梁结构电算桥梁结构计算的特点结构形式多样大型桥梁超静定次数高荷载形式复杂最终受力状态与施工方法和施工过程有关结构力学的研究内容研究结构的组成和合理形式,确定合理的计算简图研究结构内力和变形的计算方法研究结构的稳定性和动力效应结构分析的基本特点运用计算机和有限元方法进行结构内力、位移、稳定性和动力特性的研究。
方法：有限元工具：计算机本课程的基本内容1 桥梁结构受力特征及分析方法；2 重力影响的计算方法；3 活载影响的计算方法；4 其它荷载影响的计算方法；5 软件BSAS的原理和使用方法。
第一部分基本原理和方法1 1 概概述本课程的性质、特点、基本内容（1）本课程性质、特点：本课程属于专业课，旨在把学过的计算机语言、程序设计、桥梁、力学等知识结构起来，用于桥梁结构分析。
特点是既强调基本概念，又重视实际操作，基本原理与软件使用结合（结合软件“桥梁结构分析系统BSAS”教学版的使用）。
本课程的性质、特点、基本内容（2）基本内容：基本原理部分：（a）桥梁结构受力特征及分析方法；（b）重力影响的计算方法；（c）活载影响的计算方法；（d）其它荷载影响的计算方法；（e）软件BSAS的原理和使用方法。
上机操作部分（约占60%课时）主要讲解和练习软件“桥梁结构分析系统BSAS forWindows”教学版的原理和使用方法。
本课程所要求的先修课程和知识1．算法语言和程序设计（C、C++、或Fortran）；2．材料力学、结构力学、结构设计原理； 3．桥梁结构工程；4．微机操作。
第一部分基本原理和方法2 结构分析的基本方法分析方法(1)解析法建立精确的数学－物理模型，通过数学方程求解。
是一种对于模型精确求解的方法。
(2)数值法基于解析法的一种近似分析方法，包括：有限元，有限差分法，有限体积法，边界元法等有限元分析的基本概念有限元属于力学分析中的数值法，起源于航空工程中的矩阵分析，它是把一个连续的介质(或构件)看成是由有限数目的单元组成的集合体，在各单元内假定具有一定的理想化的位移和应力分布模式，各单元间通过节点相连接，并藉以实现应力的传递，各单元之间的交接面要求位移协调，通过力的平衡条件，建立一套线性方程组，求解这些方程组，便可得到各单元和结点的位移、应力。
结构力学实验报告
《结构力学》实验报告学院名称资源与安全工程学院专业班级城地专业11级 02班学生学号学生姓名郑孙斌指导老师李江腾、樊玲2013年 1月 11日目录一、实验指导书 (1)二、实验任务书 (2)三、实验名称 (2)四、实验目的 (4)五、实验步聚 (4)（1）确定单榀刚架的计算单元 (5)（2）选取平面刚架的计算简图 (5)（3）荷载计算 (5)（4）上机操作 (7)（5）输出电算实验数据 (10)六、分析实验结果 (18)七、上机操作输出图 (19)一. 《结构力学实验》指导书（一）学时与学分：实验学时：6 ；实验学分：0.5（二）实验教材：1. 《结构力学》周竞欧主编同济大学出版社2. 《结构力学实验》指导书蒋青青编3. 《PK用户手册及技术条件》中国建筑科学研究院PKPM CAD工程部（三）实验课的任务、性质与目的：实验教学是理论联系实际、培养学生实验技能和实践能力的重要教学环节。
根据城市地下空间工程专业本科培养方案和教学计划安排，结构力学是一门与工程设计密切相关的重要学科，《结构力学实验》课为三性实验（综合性、设计性、研究性）。
通过该实验课的教学，要求学生掌握结构力学的手工计算方法和电算方法，使学生熟悉结构力学电算的基本原理；正确分析结构的传力途径和准确计算作用在结构上的荷载；学会使用结构CAD系列软件PKPM中的PK计算软件来分析平面杆件结构内力和位移的实际操作步骤和方法；了解结构力学在地下空间工程中的用途，培养学生综合运用所学理论知识解决地下结构实际工程问题的动手技能和初步具备设计能力，为今后学习《混凝土结构》、《地下建筑结构》等专业课程打下一定的基础。
（四）实验项目名称：结构力学电算方法。
（五）仪器设备：在计算机房进行实验，每人一台计算机，采用中国建筑科学研究院编制的《微机结构CAD系统软件PKPM系列》中的PK平面刚架计算软件进行实验。
（六）实验方式与基本要求：1、任课教师讲清实验的基本原理、方法及要求。
【免费下载】结构力学习题难？3答案
A.全部放松； B.必须全部锁紧； C.相邻结点放松； D.相邻结点锁紧。 （ ） (D) (3 分)
203．在力矩分配法中，分配系数 AB 表示：
A.结点 A 有单位转角时，在 AB 杆 A 端产生的力矩； B.结点 A 转动时，在 AB 杆 A 端产生的力矩； C.结点 A 上作用单位外力偶时，在 AB 杆 A 端产生的力矩； D.结点 A 上作用外力偶时，在 AB 杆 A 端产生的力矩。 （ ） (C) (3 分)
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电，力根通保据过护生管高产线中工敷资艺设料高技试中术卷资，配料不置试仅技卷可术要以是求解指，决机对吊组电顶在气层进设配行备置继进不电行规保空范护载高与中带资负料荷试下卷高问总中题体资，配料而置试且时卷可，调保需控障要试各在验类最；管大对路限设习度备题内进到来行位确调。保整在机使管组其路高在敷中正设资常过料工程试况中卷下，安与要全过加，度强并工看且作护尽下关可都于能可管地以路缩正高小常中故工资障作料高；试中对卷资于连料继接试电管卷保口破护处坏进理范行高围整中，核资或对料者定试对值卷某，弯些审扁异核度常与固高校定中对盒资图位料纸置试，.卷保编工护写况层复进防杂行腐设自跨备动接与处地装理线置，弯高尤曲中其半资要径料避标试免高卷错等调误，试高要方中求案资技，料术编试交写5、卷底重电保。要气护管设设装线备备置敷4高、调动设中电试作技资气高，术料课中并3中试、件资且包卷管中料拒含试路调试绝线验敷试卷动槽方设技作、案技术，管以术来架及避等系免多统不项启必方动要式方高，案中为；资解对料决整试高套卷中启突语动然文过停电程机气中。课高因件中此中资，管料电壁试力薄卷高、电中接气资口设料不备试严进卷等行保问调护题试装，工置合作调理并试利且技用进术管行，线过要敷关求设运电技行力术高保。中护线资装缆料置敷试做设卷到原技准则术确：指灵在导活分。。线对对盒于于处调差，试动当过保不程护同中装电高置压中高回资中路料资交试料叉卷试时技卷，术调应问试采题技用，术金作是属为指隔调发板试电进人机行员一隔，变开需压处要器理在组；事在同前发一掌生线握内槽图部内 纸故，资障强料时电、，回设需路备要须制进同造行时厂外切家部断出电习具源题高高电中中源资资，料料线试试缆卷卷敷试切设验除完报从毕告而，与采要相用进关高行技中检术资查资料和料试检，卷测并主处且要理了保。解护现装场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
结构力学电算报告
1.程序功能1）求多种荷载作用下连续梁任意截面的弯矩2）求指定截面的弯矩影响线3）求弯矩包络图2.使用说明新建自定义文本文档，在第一行依次输入单元个数，支撑类型，节点荷载个数，非节点荷载个数,弹性模量。
（支撑类型说明：1.两端简支2.左端固定，右端铰支3.左端铰支，右端固定4.两端固定）接下来每行依次输入杆长，惯性矩，再接下来每行输入结点荷载大小，作用结点号。
最后每行输入非节点荷载大小，作用距离，作用单元号，荷载类型,非结点荷载类型说明：（1.集中力2.均布力3.集中力偶4.三角形分布力）运行“CHENGXU.EXE”，输入初始文件名称，根据说明按数字键后回车进行指定类型计算。
运行结果输出在“JIEGUO.TXT”中。
“WJ”表示横坐标对应的弯矩值,“WJYXX”表示横坐标对应的影响线弯矩值,“ZDWJ”表示横坐标对应的包络图最大弯矩值,“ZXWJ”表示横坐标对应的包络图最小弯矩值。
3.程序清单PROGRAM LXLDIMENSION GC(20),GX(20),PJ(20,2),PF(10,4),DK(2,2),P(45), 1F0(2),WY(2),F(2),ZK1(45),ZK2(45)character*12 dat0,dat1write(*,*)'输入:初始数据文件名'read(*,10)dat0write(*,*)'输入:计算结果数据文件名'read(*,10)dat110 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,E0WRITE(9,12)NE,NZ,NP,NF,E012 FORMAT(1X,'单元数=',6x,I3,5x,'支承类型=',I3,5x,'节点荷载个数='I3,/,1 1x,'非节点荷载个数=',5x,I3,'弹性模量=',F12.4)NJ=NE+1CALL SRSJ(NE,NP,NF,GC,GX,PJ,PF)PRINT*,'计算连续梁内力请输入1'PRINT*,'计算指定截面弯矩请输入2'PRINT*,'计算指定截面弯矩影响线请输入3'PRINT*,'计算均布活荷载作用下的包络图请输入4'READ*,IC ************************************c 计算连续梁内力C ************************************IF(I==1) THENCALL XCP(NJ,NP,NF,NE,P,PJ,PF,GC,F0,JF0)CALL JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)CALL ZCCL(NZ,NJ,ZK1,ZK2,P)CALL FCQJ(NJ,ZK1,ZK2,P)CALL SCWY(NJ,P)CALL GDL(NE,NJ,NF,E0,P,WY,DK,PF,GC,GX,F,F0,F1,F2,JF,JF0,JF1,JF2)CALL SCGDWJ(NE,F1,F2)C ************************************c 计算指定截面弯矩C ************************************ELSE IF(I==2)THENCALL XCP(NJ,NP,NF,NE,P,PJ,PF,GC,F0,FF0)CALL JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)CALL ZCCL(NZ,NJ,ZK1,ZK2,P)CALL FCQJ(NJ,ZK1,ZK2,P)CALLGDL(NE,NJ,NF,E0,P,WY,DK,PF,GC,GX,F,F0,F1,F2,FF,FF0,FF1,FF2) PRINT*,'请输入要计算内力截面位置(离连续梁最左端的距离)' READ*,WZMT=0.0DO I=1,NEML=MTMT=MT+GC(I)IF((WZ>=ML).AND.(WZ<=MT))THENWZ=WZ-MLNTWZ=IENDIFENDDOCALL WJJL(NE,F1,FF1,NF,PF,JL,WJ,WZ,NTWZ)PRINT*,' 'print*,'弯矩=',WJPRINT*,' 'call SC(NE,F1,F2,FF1,FF2,NJ,P)C ************************************c 计算指定截面弯矩影响线C ************************************ELSEIF(I==3) THENNF=1NP=0OPEN(101,FILE='ZB.TXT',STATUS='NEW')DO WHILE(WZZ<=SUM(GC)+0.03)WZZ=0.01+WZZWZ=WZZWRITE(101,*) WZENDDOCALL YXX(NE,NJ,NF,E0,P,WY,DK,PF,GC,GX,F,F0,F1,F2,FF,FF01,FF1,FF2,JL,WJ,ZK1,ZK2,NZ,PJ,NP)call SC(NE,F1,F2,FF1,FF2,NJ,P)C ************************************c 计算包络图C ************************************ELSENP=0NF=NECALL BLT(NE,NJ,NF,E0,P,WY,DK,PF,GC,GX,1F,F0,F1,F2,FF,FF0,FF1,FF2,WJ,1ZK1,ZK2,NZ,PJ,NP)ENDIFPAUSEENDC ************************************c 输入初始数据SRSJ子程序C ************************************SUBROUTINE SRSJ(NE,NP,NF,GC,GX,PJ,PF)DIMENSION GC(NE),GX(NE),PJ(NP,2),PF(NF,4)C GC(NE)：杆长 GX(NE)：惯性矩C 结点荷载:PJ(NP,1)=荷载大小 PJ(NP,2)=对应未知数序号C 非结点荷载 : PF(NE,1)=荷载大小 PF(NE,2)=距离 PF(NE,3)=单元号 PF(NE,4)=荷载类型号WRITE(9,11)11 FORMAT(1X,'杆长,惯性矩 GC(NE),GX(NE)')READ(8,*)(GC(I),GX(I),I=1,NE)write(9,12)(GC(I),GX(I),I=1,NE)12 FORMAT(10X,2F8.3)IF(NP.GT.0)THENWRITE(9,22)22 FORMAT(1X,'节点荷载大小,对应未知数序号 PJ(I,1),PJ(I,2)')READ(8,*)(PJ(I,1),PJ(I,2),I=1,NP)write(9,23)(PJ(I,1),PJ(I,2),I=1,NP)23 FORMAT(1X,2F8.3)END IFIF(NF.GT.0) THENWRITE(9,33)33 FORMAT(1X,'非结点荷载值,距离,单元号,荷载类型号') READ(8,*)((PF(I,J),J=1,4),I=1,NF)WRITE(9,34)((PF(I,J),J=1,4),I=1,NF)34 FORMAT(1X,4F8.3)END IFENDC *******************************#*************** c 计算第NHF个非结点荷截引起的等效结点荷截列阵F0 C *********************************************** SUBROUTINE DJH(NHF,NE,NF,PF,GC,F0,JF0)DIMENSION PF(NF,4),GC(NE),F0(2)DIMENSION JF0(2)！ G—荷载值,C—作用距离,NT--单元号,ID—荷载类型号 G=PF(NHF,1)C=PF(NHF,2)NT=INT(PF(NHF,3)+0.1)ID=INT(PF(NHF,4)+0.1)BL=GC(NT)D=BL-CC1=C/BLC2=D/BLGOTO(10,20,30，40),IDC 均布荷载作用下的等效荷载列阵10 F0(1)=-G*C*C*(6.0-8.0*C1+3.0*C1*C1)/12.0F0(2)=G*C*C*C*(4.0-3.0*C1)/12.0/BLJF0(1)=G*C*(1-C1*C1+C1*C1*C1/2)JF0(2)=-G*C*C1*C1*(1-C1/2)GOTO 200C 集中力作用下的等效荷截列阵20 F0(1)=-G*C*D*D/BL/BLF0(2)=G*C*C*D/BL/BLJF0(1)=G*C2*C2*(1+2*C1)JF0(2)=-G*C1*C1*(1+2*C2)GOTO 200C 集中力偶作用下的等效荷截列阵30 F0(1)=-G*C2*(2-3*C2)F0(2)=-G*C1*(2-3*C1)JF0(1)=6*G*C1*C2/BLJF0(2)=6*G*C1*C2/BLGOTO 200C 三角形分布荷载作用下的等效荷截列阵40 F0(1)=-G*C*C*(2-3*C1+1.2*C1*C1)/6F0(2)=G*C*C*C1*(1-0.8*C1)/4JF0(1)=G*C*(2-3*C1*C1+1.6*C1*C1*C1)/4JF0(2)=-G*C1*C1*C*(3-1.6*C1)/4200 RETURNENDC ******************************************c 计算第NE0个单元的单元刚度矩阵DK(2,2)C ******************************************SUBROUTINE DG(NE0,NE,E0,GC,GX,DK)DIMENSION GC(NE),GX(NE),DK(2,2)DO 15 I=1,2DO 10 J=1,2DK(I,J)=0.010 CONTINUE15 CONTINUEC DL—杆长,DI-惯性矩,S--线刚度DL=GC(NE0)DI=GX(NE0)S=E0*DI/DLDK(1,1)=4.0*SDK(1,2)=2.0*SDK(2,1)=2.0*SDK(2,2)=4.0*SENDc ********************************************************* C 集成总体刚度矩阵,ZKl,ZK2分别存主对角元素和付对角元素c ********************************************************* SUBROUTINE JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)DIMENSION ZK1(NJ),ZK2(NJ),DK(2,2),GC(NE),GX(NE)DO 10 I=1,NJZK1(I)=0.0ZK2(I)=0.010 CONTINUEDO 100 I=1,NECALL DG(I,NE,E0,GC,GX,DK)ZK1(I)=ZK1(I)+DK(1,1)ZK2(I)=ZK2(I)+DK(1,2)ZK1(I+1)=ZK1(I+1)+DK(2,2)100 CONTINUEENDC *****************************C 形成总荷裁矩阵C *****************************SUBROUTINE XCP(NJ,NP,NF,NE,P,PJ,PF,GC,F0,JF0)DIMENSION P(NJ),PJ(NP,2),PF(NF,4),GC(NE),F0(2) DIMENSION JF0(2)DO I=1,NJP(I)=0.0ENDDOIF(NP.GT.0) THENDO I=1,NPJ=INT(PJ(I,2)+0.1)P(J)=PJ(I,1)ENDDOEND IFIF(NF.GT.0) THENDO NF0=1,NFCALL DJH(NF0,NE,NF,PF,GC,F0,JF0)ND=INT(PF(NF0,3)+0.1)P(ND)=P(ND)+F0(1)P(ND+1)=P(ND+1)+F0(2)ENDDOEND IFENDc **********************************c 进行支承条件处理c **********************************SUBROUTINE ZCCL(NZ,NJ,ZK1,ZK2,P)DIMENSION ZK1(NJ),ZK2(NJ),P(NJ)GOTO(10,20,30,20),NZ10 GOTO 10020 ZK1(1)=1.0P(1)=0.0ZK2(1)=0.0IF(NZ.EQ.4) GOTO 30GOTO 10030 ZK1(NJ)=1.0ZK2(NJ-1)=0.0P(NJ)=0.0100 CONTINUEENDc *******************************c 输出位移c *******************************SUBROUTINE SCWY(NJ,P)DIMENSION P(NJ)WRITE(9,10)10 FORMAT(1X,'：：：：：：：：：位移：；：：：：：：：：')DO 100 I=1,NJWRITE(9,20)I,P(I)20 FORMAT(1X,'结点号=',I2,5X,F10.4)100 CONTINUEENDc ******************************c 计算单元杆端内力c *****************SUBROUTINE GDL(NE,NJ,NF,E0,P,WY,DK,PF,GC,GX, 1F,F0,F1,F2,JF,JF0,JF1,JF2)DIMENSION P(NJ),WY(2),DK(2,2),PF(NF,4)DIMENSION GC(NE),GX(NE)DIMENSION F(2),F0(2),F1(NE),F2(NE)DIMENSION JF(2),JF0(2),JF1(NE),JF2(NE)DO NE0=1,NECALL DG (NE0,NE,E0,GC,GX,DK)WY(1)=P(NE0)WY(2)=P(NE0+1)DO I=1,2F(I)=0.0JF(I)=0.0DO J=1,2F(I)=F(I)+DK(I,J)*WY(J)JF(I)=JF(I)+6*E0*GX(NE0)*WY(J)/(GC(NE0)*GC(NE0)) ENDDOENDDOIF(NF.GT.0) THENDO I=1,NFIF(INT(PF(I,3)+0.1).EQ.NE0) THENCALL DJH(I,NE,NF,PF,GC,F0,JF0)DO J=1,2F(J)=F(J)-F0(J)JF(J)=JF(J)+JF0(J)ENDDOEND IFENDDOEND IFF1(NE0)=F(1)F2(NE0)=F(2)ENDDOENDc *******************************c 输出杆端弯矩c *******************************SUBROUTINE SCGDWJ(NE,F1,F2)DIMENSION F1(NE),F2(NE)WRITE(9,2)2 FORMAT(1X,'.................各单元杆端内力....................')WRITE(9,150) NE0,F(1),F(2)150 FORMAT(1X,'单元号=',I2,5X,'左端弯矩=',F9.3,2X,'右端弯矩=',F9.3)200 CONTINUEWRITE(9,11)11 FORMAT('====================== 计算结束======================')ENDc *******************************c 计算指定截面弯矩c *******************************SUBROUTINE ZDWJ(NE,F1,JF1,NF,PF,WJ,WZ,NTWZ)DIMENSION F1(NE),JF1(NE)DIMENSION PF(NF,4)DO I=1,NEIF(NTWZ==I)THENWJ=-F1(I)+JF1(I)*WZENDIFENDDODO K=1,NFG=PF(K,1)C=PF(K,2)NT=INT(PF(K,3)+0.1)ID=INT(PF(K,4)+0.1)IF(NT==NTWZ)THENSELECT CASE(ID)CASE(1)IF(WZ<=C) THENWJ=WJ-G*(WZ**2)/2ELSEWJ=WJ-G*(C**2)/2-(WZ-C)*C*GENDIFCASE(2)IF(WZ<=C) THENWJ=WJELSEWJ=WJ-G*(WZ-C)ENDIFCASE(3)IF(WZ<=C) THENWJ=WJELSEWJ=WJ-GENDIFCASE(4)IF(WZ<=C) THENWJ=WJ-(WZ*WZ*WZ)*G/6/CELSEWJ=WJ-(C*C*C)*G/6/C-(WZ-C)*C*C*G/2/CENDIFENDSELECTENDIFENDDOENDC ***********************************************c 计算指定截面弯矩影响线C ***********************************************SUBROUTINE ZDYXX(NE,NJ,NF,E0,P,WY,DK,PF,GC,GX,F,F0,F1,F2,JF,JF0,JF1,JF2,TORTION,ZK1,ZK2,NZ,PJ,NP)DIMENSION P(NJ),WY(2),DK(2,2),PF(NF,4),PJ(NP,2) DIMENSION GC(NE),GX(NE)DIMENSION F(2),F0(2),F1(NE),F2(NE)DIMENSION JF(2),JF0(2),JF1(NE),JF2(NE)DIMENSION ZK1(NJ),ZK2(NJ)PRINT*,'请输入要计算内力截面位置(离连续梁最左端的距离)' READ*,WZMT=0.0DO I=1,NEML=MTMT=MT+GC(I)IF((WZ>=ML).AND.(WZ<=MT))THENWZ=WZ-MLNTWZ=IENDIFENDDOPF(1,1)=1PF(1,4)=2ME=0OPEN(12,FILE='YXXT.TXT',STATUS='NEW')CLOSE(12)OPEN(13,FILE='YXXS.TXT',STATUS='NEW')CLOSE(13)DO I=1,NEML=MEME=ME+GC(I)DO AJ=ML,ME,0.01PF(1,2)=AJ-MLPF(1,3)=ICALL XCP(NJ,NP,NF,NE,P,PJ,PF,GC,F0,FF0) CALL JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)CALL ZCCL(NZ,NJ,ZK1,ZK2,P)CALL FCQJ(NJ,ZK1,ZK2,P)CALLGDL(NE,NJ,NF,E0,P,WY,DK,PF,GC,GX,F,F0,F1,F2,FF,FF0,FF1,FF2) CALL WJJL(NE,F1,FF1,NF,PF,TORTION,WZ,NTWZ)OPEN(12,FILE='WJYXX.TXT',STATUS='OLD')WRITE(12,*) WJENDDOENDDOEND！***********************************************！计算均布荷载弯矩包络图！***********************************************SUBROUTINE BLT(NE,NJ,NF,E0,P,WY,DK,PF,GC,GX,1F,F0,F1,F2,FF,FF0,FF1,FF2,1SHEAR,TORTION,1ZK1,ZK2,NZ,PJ,NP)DIMENSION P(NJ),WY(2),DK(2,2),PF(NF,4),PJ(NP,2)DIMENSION GC(NE),GX(NE)DIMENSION F(2),F0(2),F1(NE),F2(NE)DIMENSION FF(2),FF0(2),FF1(NE),FF2(NE)DIMENSION ZK1(NJ),ZK2(NJ)OPEN(31,FILE='ZXJL.TXT',STATUS='NEW')CLOSE(31)OPEN(41,FILE='ZDWJ.TXT',STATUS='NEW')CLOSE(41)OPEN(51,FILE='ZXWJ.TXT',STATUS='NEW')CLOSE(51)OPEN(101,FILE='ZB.TXT',STATUS='NEW')PRINT*,'请输入均布荷载大小'READ*,DAXIAODO WZZ=0,SUM(GC),0.01SHEAR_MAX=-100000SHEAR_MIN=100000TORTION_MAX=-100000TORTION_MIN=100000WZ=WZZwrite(101,*) WZZMT=0.0DO I=1,NEML=MTMT=MT+GC(I)IF((WZ>=ML).AND.(WZ<=MT+0.02))THENWZ=WZ-MLNTWZ=IENDIFENDDODO J=1,2DO I=1,NEPF(I,1)=0PF(I,2)=GC(I)PF(I,3)=IPF(I,4)=1ENDDOIF(J==1) THENDO I=1,NE,2PF(I,1)=DAXIAOENDDOELSEIF(J==2) THENDO I=2,NE,2PF(I,1)=DAXIAOENDDOENDIFCALL XCP(NJ,NP,NF,NE,P,PJ,PF,GC,F0,FF0) CALL JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)CALL ZCCL(NZ,NJ,ZK1,ZK2,P)CALL FCQJ(NJ,ZK1,ZK2,P)CALLGDL(NE,NJ,NF,E0,P,WY,DK,PF,GC,GX,F,F0,F1,F2,FF,FF0,FF1,FF2) CALL WJJL(NE,F1,FF1,NF,PF,SHEAR,TORTION,WZ,NTWZ)IF(TORTION>TORTION_MAX) TORTION_MAX=TORTIONIF(TORTION<TORTION_MIN) TORTION_MIN=TORTIONIF(SHEAR>SHEAR_MAX) SHEAR_MAX=SHEARIF(SHEAR<SHEAR_MIN) SHEAR_MIN=SHEARENDDODO JJ=1,NE-1DO I=1,NEPF(I,1)=0PF(I,2)=GC(I)PF(I,3)=IPF(I,4)=1ENDDODO K=JJ,1,-2PF(K,1)=DAXIAOENDDODO K=JJ+1,NE,2PF(K,1)=DAXIAOENDDOCALL XCP(NJ,NP,NF,NE,P,PJ,PF,GC,F0,FF0)CALL JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)CALL ZCCL(NZ,NJ,ZK1,ZK2,P)CALL FCQJ(NJ,ZK1,ZK2,P)CALLGDL(NE,NJ,NF,E0,P,WY,DK,PF,GC,GX,F,F0,F1,F2,FF,FF0,FF1,FF2)CALL WJJL(NE,F1,FF1,NF,PF,SHEAR,TORTION,WZ,NTWZ)IF(TORTION>TORTION_MAX) TORTION_MAX=TORTIONIF(TORTION<TORTION_MIN) TORTION_MIN=TORTIONIF(SHEAR>SHEAR_MAX) SHEAR_MAX=SHEARIF(SHEAR<SHEAR_MIN) SHEAR_MIN=SHEARENDDOOPEN(21,FILE='ZDJL.TXT',STATUS='OLD')WRITE(21,*) SHEAR_MAXOPEN(31,FILE='ZXJL.TXT',STATUS='OLD')WRITE(31,*) SHEAR_MINOPEN(41,FILE='ZDWJ.TXT',STATUS='OLD')WRITE(41,*) TORTION_MAXOPEN(51,FILE='ZXWJ.TXT',STATUS='OLD')WRITE(51,*) TORTION_MINENDDOEND4. 算例10KN 50KN·M 80KN·M 100KN·M 40KN·M 20KN/M 12KN/M7m 8m 6m 4m(1)连续梁内力计算结果（2）.指定截面内力计算结果（3）.指定截面弯矩影响线X=13.5（4）计算弯矩包络图荷载大小为10KN0510********-1.4-1.2-1-0.8-0.6-0.4-0.200.20.4弯矩图0510********-80-60-40-2020406080弯矩包络图 最大弯矩 最小弯矩。