结构力学上机实验报告.doc
实验报告一 平面刚架内力计算程序APF日期：2012年3月23日星期五实验地点：综合楼504实验目的：(1)分析构件刚度与外界温度对结构位移的影响，如各杆刚度改变对内力分布的影响、温度因数对内力分布的影响。(2)观察并分析刚架在静力荷载及温度作用下的内力和变形规律，包括刚度的变化，结构形式的改变，荷载的作用位置变化等因素对内力及变形的影响。对结构静力分析的矩阵位移法的计算机应用有直观的了解(3)掌握杆系结构计算的结构力学求解器的使用方法。通过实验加深对静定、超静定结构特性的认识。实验设计1： 计算图示刚架当梁柱刚度分别为、时结构的内力和位移，由此分析当刚架在水平荷载作用下横梁的水平位移与刚架梁柱比（）之间的关系。(计算时忽略轴向变形)。数据文件：(1)变量定义,EI1=1,EI2=0.2(1,5,10)结点,1,0,0 结点,2,0,4 结点,3,6,4 结点,4,6,0单元,1,2,1,1,1,1,1,1 单元,2,3,1,1,1,1,1,1 单元,3,4,1,1,1,1,1,1结点支承,1,6,0,0,0,0 结点支承,4,6,0,0,0,0 结点荷载,2,1,100,0单元材料性质,1,1,-1,EI1,0,0,-1 单元材料性质,2,2,-1,EI2,0,0,-1 单元材料性质,3,3,-1,EI1,0,0,-1(2)变量定义,EI1=5(1,0.2,0.1),EI2=1结点,1,0,0 结点,2,0,4 结点,3,6,4 结点,4,6,0单元,1,2,1,1,1,1,1,1 单元,2,3,1,1,1,1,1,1 单元,3,4,1,1,1,1,1,1结点支承,1,6,0,0,0,0 结点支承,4,6,0,0,0,0 结点荷载,2,1,100,0单元材料性质,1,1,-1,EI1,0,0,-1 单元材料性质,2,2,-1,EI2,0,0,-1 单元材料性质,3,3,-1,EI1,0,0,-1主要计算结果： 位移：(1)令I1=1时，I2=0.2,1,5,10 梁柱刚度比I2：I11：51：15：110：1横梁水平位移711.11/EI1(m)426.67/EI1(m)304.76/EI1(m)286.18/EI1(m)(2)令I2=1时，I1=5,1,0.2,0.1 梁柱刚度比I2：I11：51：15：110：1横梁水平位移142.22/EI1(m)426.67/EI1(m)1523.81/EI1(m)2861.78/EI1(m)弯矩：(1) 令I1=1时，I2=0.2,1,5,10梁柱刚度比I2：I1为1：5时的刚架弯矩图如下： 梁柱刚度比I2：I1为1：1时的刚架弯矩图如下： M图（单位：KN·m） M图（单位：KN·m）梁柱刚度比I2：I1为5：1时的刚架弯矩图如下：梁柱刚度比I2：I1为10：1时的刚架弯矩图如下： M图（单位：KN·m）M图（单位：KN·m）(2)令I2=1时，I1=5,1,0.2,0.1 梁柱刚度比I2：I1为1：5时的刚架弯矩图如下： 梁柱刚度比I2：I1为1：1时的刚架弯矩图如下：M图（单位：KN·m） M图（单位：KN·m）梁柱刚度比I2：I1为5：1时的刚架弯矩图如下：梁柱刚度比I2：I1为10：1时的刚架弯矩图如下：M图（单位：KN·m） M图（单位：KN·m） 三、结果分析及结论：无论EI1和EI2的值如何改变，只要EI2：EI1的值不改变，那么刚架的弯矩图都是相同的；且随着梁柱刚度比EI2：EI1的增大，两柱的弯矩的反弯点向下移动；横梁的弯矩的反弯点保持在中点不变；当I1=1，I2=0.2,1,5,10时，随着梁柱刚度比EI2：EI1的增大，刚架在水平荷载作用下的横梁的水平位移变小711.11/EI1(m)426.67/EI1(m) 304.76/EI1(m) 286.18/EI1(m)；当I2=1，I1=5,1,0.2,0.1时，随着梁柱刚度比EI2：EI1的增大，刚架在水平荷载作用下的横梁的水平位移变大142.22/EI1(m) 426.67/EI1(m) 1523.81/EI1(m) 2861.78/EI1(m)，且其变化的幅度远远大于当I1=1，I2=0.2,1,5,10时的幅度（因为I1=5,1,0.2,0.1是慢慢变小的）；当I1=1，I2=0.2,1,5,10或当I2=1，I1=5,1,0.2,0.1时，随着梁柱刚度比EI2：EI1的增大，梁柱交点处的梁端与柱端的弯矩逐渐变大（44.4480.0095.2497.56）（单位：KN·m）；柱底端弯矩逐渐变小（155.56120.00104.76102.44）（单位：KN·m）。结论：在水平荷载保持不变的情况下：增大梁柱的抗弯刚度比值可使梁端的弯矩值增大，同时也会使柱底端的弯矩值变小，即刚架梁柱刚度比的变化影响梁柱的内力分配；梁的水平位移与梁柱的抗弯刚度比值关系不大，而与柱的抗弯刚度的大小的关系大，柱的抗弯刚度越小，梁的水平位移越大，即横梁的水平位移主要由柱刚度控制；横梁的水平位移随结构总刚度的增加而减少；实验设计2： 求作图示刚架温度变化时的弯矩图。设N/cm2 ，各杆的截面尺寸均为50cm×40cm。试分析当改变截面尺寸（设各杆截面尺寸均为50cm×60cm，及50cm×80cm）时，对结构的内力的影响。如果截面尺寸均为50cm×40cm，温度变化为内侧降温300C，外侧升温100C，分析温度变化对内力的影响。 结点、单元编号示意图一、 数据文件：(1)第一问的数据文件：变量定义,E=1.5e7,B=0.5,H=0.4（0.6，0.8）,EI=E*B*H*H*H/12,EA=E*B*H结点,1,0,0 结点,2,0,8 结点,3,6,0 结点,4,6,8 结点,5,12,0结点,6,12,8 结点,7,18,0 结点,8,18,8 结点,9,24,0 结点,10,24,8单元,1,2,1,1,1,1,1,1 单元,2,4,1,1,1,1,1,1 单元,4,6,1,1,1,1,1,1单元,6,8,1,1,1,1,1,1 单元,8,10,1,1,1,1,1,1 单元,9,10,1,1,1,1,1,1单元,3,4,1,1,1,1,1,1 单元,5,6,1,1,1,1,1,1 单元,7,8,1,1,1,1,1,1结点支承,1,6,0,0,0,0 结点支承,3,6,0,0,0,0结点支承,5,6,0,0,0,0 结点支承,7,6,0,0,0,0 结点支承,9,6,0,0,0,0单元材料性质,1,9,EA,EI,0,0,-1 单元温度改变,1,5,-10,-40,0.00001,H单元温度改变,6,6,-10,40,0.00001,H 单元温度改变,7,9,10,0,0.00001,H(2)第二问的数据文件：变量定义,E=1.5e7,B=0.5,H=0.4,EI=E*B*H*H*H/12,EA=E*B*H结点,1,0,0 结点,2,0,8 结点,3,6,0 结点,4,6,8 结点,5,12,0结点,6,12,8 结点,7,18,0 结点,8,18,8 结点,9,24,0 结点,10,24,8单元,1,2,1,1,1,1,1,1 单元,2,4,1,1,1,1,1,1 单元,4,6,1,1,1,1,1,1单元,6,8,1,1,1,1,1,1 单元,8,10,1,1,1,1,1,1 单元,9,10,1,1,1,1,1,1单元,3,4,1,1,1,1,1,1 单元,5,6,1,1,1,1,1,1 单元,7,8,1,1,1,1,1,1结点支承,1,6,0,0,0,0 结点支承,3,6,0,0,0,0 结点支承,5,6,0,0,0,0 结点支承,7,6,0,0,0,0 结点支承,9,6,0,0,0,0单元材料性质,1,9,EA,EI,0,0,-1 单元温度改变,1,5,-10,40,0.00001,H单元温度改变,6,6,-10,-40,0.00001,H 单元温度改变,7,9,-30,0,0.00001,H二、计算结果（弯矩和轴力）：第一问的弯矩和轴力图如下： 当H0.4, B=0.5时刚架的弯矩图和轴力图： M图（单位：KN·m） N图（单位：KN） 当H0.6, B=0.5时刚架的弯矩图和轴力图： M图（单位：KN·m） N图（单位：KN） 当H0.8, B=0.5时刚架的弯矩图和轴力图： M图（单位：KN·m） N图（单位：KN）(2)第二问的弯矩和轴力图如下：当H0.4, B=0.5内侧降温300C，外侧升温100C时刚架的弯矩图和轴力图： M图（单位：KN·m） N图（单位：KN）三、结果分析及结论：（定义：对称柱与刚架对称轴重合的柱子；中柱对称柱与边柱之间的柱子）由第一问的结果可知，当刚架外侧降温300C，内侧升温100C时： 在刚架截面的宽度不变（50cm），随着高度增大（分别为40cm、60cm、80c），有弯矩的杆件的弯矩值都增大，所有杆件的轴力都增大； 对于上下表面温差不为零的杆件（即边柱和梁），温度降低的一侧，杆件受拉；温度升高的一侧，杆件受压； 而刚架内部上下表面温差为零的杆件（即两根中柱和对称柱），两中柱底端外侧受拉，顶端里侧受拉，对称柱没有弯矩； 从轴力图可知，两中柱受压，其它杆件受拉。由第二问的结果可知，当刚架的外侧升温100C，内侧降温300C时，在刚架的宽度为50cm，高度为40cm的情况下： 上下表面温差不为零的杆件（即边柱和梁），降温的一侧，杆件受拉；升温的一侧，杆件受压； 而刚架内部上下表面温差不变的杆件（即两根中柱和对称柱），两中柱底端外侧受拉，顶端里侧受拉，对称柱没有弯矩； 从轴力图可知，两根边柱和对称柱受压，中柱和梁受拉。由此可得出结论： 当杆件有温差时，弯矩图的竖矩出现在降温面的一侧，升温面产生压应力，降温面产生拉应力； 随着杆件截面高度的增大（宽度不变，即刚度增大），各有弯矩的杆件的弯矩值增大，有轴力的杆件轴力也增大，即刚度变化影响内力的变化，刚度越大，内力越大，反之，内力越小； 两端约束作用下，杆件轴心降温轴力为正，升温为负。实验报告二 平面桁架内力计算日期：2012年3月30日星期五实验地点：综合楼504实验目的：(1)考察并分析桁架的结构形式、刚度、荷载作用位置等因素对桁架内力及变形的影响及规律。（2）对比平行桁架和三角桁架的受力特点及杆件内力的变化规律；平行弦桁架和三角形桁架在相同荷载作用下内力的不同；荷载上承（荷载作用在上弦杆上）和下承（荷载作用在相应的下弦杆上）时各杆内力的变化。（3）掌握杆系结构计算的结构力学求解器的使用方法和对矩阵位移法计算位移的了解。通过实验加深对静定、超静定结构特性的认识。实验设计： 计算图示桁架的内力和位移。分析：1. 刚度对内力和位移（跨中竖向位移）的影响；2.平行弦桁架和三角形桁架在相同荷载作用下内力的不同；3.荷载上承（荷载作用在上弦杆上，如图）和下承（荷载作用在相应的下弦杆上）时各杆内力的变化。 （a） (b) 二、 数据文件：（1） 图(a)所示的桁架的数据文件如下：(荷载上承)变量定义,EA=1（2，5，10）结点,1,0,0 结点,2,0,2 结点,3,2,0 结点,4,2,2 结点,5,4,0 结点,6,4,2 结点,7,6,0 结点,8,6,2 结点,9,8,0 结点,10,8,2 结点,11,10,0 结点,12,10,2 结点,13,12,0 结点,14,12,2单元,1,2,1,1,0,1,1,0 单元,2,3,1,1,0,1,1,0 单元,3,4,1,1,0,1,1,0 单元,4,5,1,1,0,1,1,0单元,5,6,1,1,0,1,1,0 单元,6,7,1,1,0,1,1,0 单元,7,8,1,1,0,1,1,0 单元,7,10,1,1,0,1,1,0单元,9,10,1,1,0,1,1,0 单元,9,12,1,1,0,1,1,0 单元,11,12,1,1,0,1,1,0 单元,11,14,1,1,0,1,1,0单元,13,14,1,1,0,1,1,0 单元,1,3,1,1,0,1,1,0 单元,2,4,1,1,0,1,1,0 单元,3,5,1,1,0,1,1,0单元,4,6,1,1,0,1,1,0 单元,5,7,1,1,0,1,1,0 单元,6,8,1,1,0,1,1,0 单元,7,9,1,1,0,1,1,0单元,8,10,1,1,0,1,1,0 单元,9,11,1,1,0,1,1,0 单元,10,12,1,1,0,1,1,0 单元,11,13,1,1,0,1,1,0单元,12,14,1,1,0,1,1,0 结点支承,1,3,0,0,0 结点支承,13,1,0,0单元材料性质,1,25,EA,0.00000001,0,0,-1结点荷载,2,1,0.5,-90 结点荷载,4,1,1,-90 结点荷载,6,1,1,-90 结点荷载,8,1,1,-90结点荷载,10,1,1,-90 结点荷载,12,1,1,-90 结点荷载,14,1,0.5,-90(荷载下承)变量定义,EA=1结点,1,0,0 结点,2,0,2 结点,3,2,0 结点,4,2,2 结点,5,4,0 结点,6,4,2 结点,7,6,0 结点,8,6,2结点,9,8,0 结点,10,8,2 结点,11,10,0 结点,12,10,2 结点,13,12,0 结点,14,12,2单元,1,2,1,1,0,1,1,0 单元,2,3,1,1,0,1,1,0 单元,3,4,1,1,0,1,1,0 单元,4,5,1,1,0,1,1,0 单元,5,6,1,1,0,1,1,0单元,6,7,1,1,0,1,1,0 单元,7,8,1,1,0,1,1,0 单元,7,10,1,1,0,1,1,0 单元,9,10,1,1,0,1,1,0单元,9,12,1,1,0,1,1,0 单元,11,12,1,1,0,1,1,0 单元,11,14,1,1,0,1,1,0 单元,13,14,1,1,0,1,1,0单元,1,3,1,1,0,1,1,0 单元,2,4,1,1,0,1,1,0 单元,3,5,1,1,0,1,1,0 单元,4,6,1,1,0,1,1,0 单元,5,7,1,1,0,1,1,0单元,6,8,1,1,0,1,1,0 单元,7,9,1,1,0,1,1,0 单元,8,10,1,1,0,1,1,0 单元,9,11,1,1,0,1,1,0单元,10,12,1,1,0,1,1,0 单元,11,13,1,1,0,1,1,0 单元,12,14,1,1,0,1,1,0结点支承,1,3,0,0,0 结点支承,13,1,0,0 单元材料性质,1,25,EA,1,0,0,-1结点荷载,1,1,0.5,-90 结点荷载,3,1,1,-90 结点荷载,5,1,1,-90 结点荷载,7,1,1,-90 结点荷载,9,1,1,-90结点荷载,11,1,1,-90 结点荷载,13,1,0.5,-90（2）图(b)所示的桁架的数据文件如下：(荷载上承)变量定义,EA=1结点,1,0,0 结点,2,2,0 结点,3,2,2/3 结点,4,4,0 结点,5,4,4/3 结点,6,6,0结点,7,6,2 结点,8,8,0 结点,9,8,4/3 结点,10,10,0 结点,11,10,2/3 结点,12,12,0单元,1,2,1,1,0,1,1,0 单元,2,3,1,1,0,1,1,0 单元,3,4,1,1,0,1,1,0 单元,4,5,1,1,0,1,1,0单元,5,6,1,1,0,1,1,0 单元,6,7,1,1,0,1,1,0 单元,6,9,1,1,0,1,1,0 单元,8,9,1,1,0,1,1,0单元,8,11,1,1,0,1,1,0 单元,10,11,1,1,0,1,1,0 单元,10,12,1,1,0,1,1,0 单元,1,3,1,1,0,1,1,0单元,3,5,1,1,0,1,1,0 单元,5,7,1,1,0,1,1,0 单元,7,9,1,1,0,1,1,0 单元,9,11,1,1,0,1,1,0单元,11,12,1,1,0,1,1,0 单元,2,4,1,1,0,1,1,0 单元,4,6,1,1,0,1,1,0 单元,6,8,1,1,0,1,1,0单元,8,10,1,1,0,1,1,0结点支承,1,3,0,0,0 结点支承,12,1,0,0 结点荷载,1,1,0.5,-90 结点荷载,3,1,1,-90结点荷载,5,1,1,-90 结点荷载,7,1,1,-90 结点荷载,9,1,1,-90 结点荷载,11,1,1,-90结点荷载,12,1,0.5,-90 单元材料性质,1,21,EA,0.00000001,0,0,-1(荷载上承)变量定义,EA=1结点,1,0,0 结点,2,2,0 结点,3,2,2/3 结点,4,4,0 结点,5,4,4/3 结点,6,6,0 结点,7,6,2结点,8,8,0 结点,9,8,4/3 结点,10,10,0 结点,11,10,2/3 结点,12,12,0单元,1,2,1,1,0,1,1,0 单元,2,3,1,1,0,1,1,0 单元,3,4,1,1,0,1,1,0 单元,4,5,1,1,0,1,1,0 单元,5,6,1,1,0,1,1,0单元,6,7,1,1,0,1,1,0 单元,6,9,1,1,0,1,1,0 单元,8,9,1,1,0,1,1,0 单元,8,11,1,1,0,1,1,0单元,10,11,1,1,0,1,1,0 单元,10,12,1,1,0,1,1,0 单元,1,3,1,1,0,1,1,0 单元,3,5,1,1,0,1,1,0单元,5,7,1,1,0,1,1,0 单元,7,9,1,1,0,1,1,0 单元,9,11,1,1,0,1,1,0 单元,11,12,1,1,0,1,1,0单元,2,4,1,1,0,1,1,0 单元,4,6,1,1,0,1,1,0 单元,6,8,1,1,0,1,1,0 单元,8,10,1,1,0,1,1,0结点支承,1,3,0,0,0 结点支承,12,1,0,0 结点荷载,1,1,0.5,-90 结点荷载,2,1,1,-90 结点荷载,4,1,1,-90结点荷载,6,1,1,-90 结点荷载,8,1,1,-90 结点荷载,10,1,1,-90 结点荷载,12,1,0.5,-90单元材料性质,1,21,EA,1,0,0,-1二、计算结果：1. 结构各单元内力：（1）图（a）所示桁架在图中力的作用下产生的单元内力图如下： N图（荷载上承） N图（荷载下承） (2) 图（b）所示桁架在图中力的作用下产生的单元内力图如下：N图（荷载上承） N图（荷载下承）2. 不同刚度时，图（a）跨中的竖向位移。跨中下弦位置EA=1EA=2EA=5EA=10结点竖向位移-109.456/EA（m）-54.728/EA（m）-21.891/EA（m）-10.946/EA（m）三、结果分析：刚度对内力和位移（跨中竖向位移）的影响：在相同荷载的作用下，平行弦桁架随着杆件刚度的增大，桁架内力不变，跨中的竖向位移变小；平行弦桁架和三角形桁架在相同荷载作用下内力的不同： 在相同荷载的作用下，平行弦桁架的上弦杆受压，且由外向里（由边直腹杆往对称轴的方向）杆件受力是逐渐变大的，三角形桁杆的上弦杆也受压，但由外向里杆件受力是逐渐变小的； 在相同荷载的作用下，平行弦桁架的下弦杆受拉（其中下弦边杆受力为零），且由外向里杆件受力是逐渐变大的，三角形桁杆的下弦杆也受拉，但由外向里杆件受力情况是：下弦边杆与第二边杆受力相同，且比中杆受力大； 在相同荷载的作用下，平行弦桁架的直腹杆受压，且由外向里受力是逐渐变小的，三角形桁杆的直腹杆受拉，且由外向里受力是逐渐变大的（三角形桁杆的边直腹杆受力为零）； 在相同荷载的作用下，平行弦桁架的斜腹杆受拉，且由外向里受力是逐渐变小的，三角形桁杆的斜腹杆受压，且由外向里受力是逐渐变大的。荷载上承和下承时各杆内力的变化： 对于平行弦桁架而言，荷载上承与荷载下承作用下，上弦杆、下弦杆和斜腹杆的受力情况是相同的； 对于平行弦桁架而言，荷载上承作用下的直腹杆的受力大小比在相同位置荷载下承作用下的直腹杆的受力大小大1。 对于三角形桁杆而言，荷载上承与荷载下承作用下，上弦杆、下弦杆和斜腹杆的受力情况是相同的； 对于三角形桁杆而言，荷载上承作用下的直腹杆的受力大小比在相同位置荷载下承作用下的直腹杆的受力大小小1。实验报告三 平面任意杆系内力计算日期：2012年3月30日星期五实验地点：综合楼504实验目的：(1)求解平面任意杆系内力，铰接排架的剪力计算，多层框架结构的内力、动力、稳定计算。(2) 通过试验，计算当铰接排架各柱刚度比值不同时，各排架柱的柱底弯矩和剪力，从而分析刚度比对铰接排架柱内力的影响。通过计算内力，分析各柱剪力及弯矩分配、横梁水平位移的规律；学会运用结构力学求解器求解结构动力问题。实验设计1： 图示铰接排架，设各杆弹性模量为一常数，计算时取E=1，并且忽略轴向变形，计算当I1I2I3I4 分别为10.50.51、1111、1234时各排架柱的柱底弯矩和剪力。试分析各柱剪力及弯矩分配、横梁水平位移的规律。 结点、单元编号示意图一、数据文件：变量定义,E=1,h=2变量定义,I1=1,I2=0.5（1，2）,I3=0.5（1，3）,I4=1（1，4）结点,1,0,0 结点,2,2,0 结点,3,4,0 结点,4,6,0 结点,5,0,h结点,6,2,h 结点,7,4,h 结点,8,6,h单元,1,5,1,1,1,1,1,0 单元,2,6,1,1,1,1,1,0 单元,3,7,1,1,1,1,1,0 单元,4,8,1,1,1,1,1,0单元,5,6,1,1,0,1,1,0 单元,6,7,1,1,0,1,1,0 单元,7,8,1,1,0,1,1,0结点支承,1,6,0,0,0,0 结点支承,2,6,0,0,0,0 结点支承,3,6,0,0,0,0 结点支承,4,6,0,0,0,0 结点荷载,5,1,1,0变量定义,EI1=E*I1 变量定义,EI2=E*I2变量定义,EI3=E*I3 变量定义,EI4=E*I4单元材料性质,1,1,-1,EI1,0,0,-1 单元材料性质,2,2,-1,EI2,0,0,-1单元材料性质,3,3,-1,EI3,0,0,-1 单元材料性质,4,4,-1,EI4,0,0,-1单元材料性质,5,7,-1,1,0,0,-1一、 有关计算结果摘录：梁柱刚度比横梁水平位移柱1剪力及柱底弯矩柱2剪力及柱底弯矩柱3剪力及柱底弯矩柱4剪力及柱底弯矩1:0.5:0.5:10.111*h*h*h0.33/-0.33h0.17/-0.17h0.17/-0.17h0.33/-0.33h1:1:1:10.083*h*h*h0.25/-0.25h0.25/-0.25h0.25/-0.25h0.25/-0.25h1:2:3:40.033*h*h*h0.10/0.10h0.20/0.20h0.30/0.30h0.40/0.40h内力图 当排架柱刚度比为1:0.5:0.5:1时的内力图如下:N图 Q图 M图 排架柱刚度比为1:1:1:1时的内力图如下:N图 Q图 M图 当排架柱刚度比为1:2:3:4时的内力图如下: N图 Q图 M图二、 结果分析及结论：结果分析： 在相同荷载作用下，横梁的产生的轴力与柱的抗弯刚度比有关； 在相同荷载作用下，柱产生的剪力的比值与柱的抗弯刚度比值相等； 在相同荷载作用下，柱底的弯矩比值与抗弯刚度比值相等； 在相同荷载作用下，柱刚度越大，横梁的水平位移越小；结论： 在相同荷载作用下及相同柱刚度比条件下，高为h的柱的横梁产生的水平位移是高为1的柱的横梁产生的水平位移的h的三次方倍； 在相同荷载作用下及相同柱刚度比条件下，高为h的柱所产生柱底弯矩是高为1的柱所产生的柱底弯矩的h 倍； 通过结构力学求解器计算得出如上内力图。由计算结果可以看出，排架水平位移随着柱子刚度比的不断增大，水平位移呈减小趋势。在水平荷载的作用下，各柱子柱顶剪力按刚度比例分配，剪力综合与外荷载相等。由于各柱子高度相同，故各柱顶弯矩也按刚度比例分配。实验设计2： 分析图示五层钢架的内力、动力特性、稳定分析。梁截面尺寸mm，柱截面尺寸mm，材料弹性模量。一、 计算结果：数据文件：变量定义,I1=0.2*0.3*0.3*0.3/12 变量定义,I2=0.3*0.3*0.3*0.3/12 变量定义,E=2.1E11变量定义,EI1=E*I1,EI2=E*I2 变量定义,A1=0.2*0.3,A2=0.3*0.3 变量定义,EA1=E*A1,EA2=E*A2变量定义,P=7.8E3 变量定义,PA1=P*A1,PA2=P*A2结点,1,0,0 结点,2,6,0 结点,3,12,0 结点,4,0,4 结点,5,6,4 结点,6,12,4 结点,7,0,8 结点,8,6,8 结点,9,12,8 结点,10,0,12 结点,11,6,12 结点,12,12,12结点,13,0,16结点,14,6,16 结点,15,12,16 结点,16,0,20 结点,17,6,20 结点,18,12,20单元,1,4,1,1,1,1,1,1 单元,4,7,1,1,1,1,1,1 单元,7,10,1,1,1,1,1,1 单元,10,13,1,1,1,1,1,1单元,13,16,1,1,1,1,1,1 单元,2,5,1,1,1,1,1,1 单元,5,8,1,1,1,1,1,1 单元,8,11,1,1,1,1,1,1单元,11,14,1,1,1,1,1,1 单元,14,17,1,1,1,1,1,1 单元,3,6,1,1,1,1,1,1 单元,6,9,1,1,1,1,1,1单元,9,12,1,1,1,1,1,1 单元,12,15,1,1,1,1,1,1 单元,15,18,1,1,1,1,1,1 单元,4,5,1,1,1,1,1,1单元,5,6,1,1,1,1,1,1 单元,7,8,1,1,1,1,1,1 单元,8,9,1,1,1,1,1,1 单元,10,11,1,1,1,1,1,1单元,11,12,1,1,1,1,1,1 单元,13,14,1,1,1,1,1,1 单元,14,15,1,1,1,1,1,1 单元,16,17,1,1,1,1,1,1单元,17,18,1,1,1,1,1,1结点支承,1,6,0,0,0,0 结点支承,2,6,0,0,0,0 结点支承,3,6,0,0,0,0单元荷载,16,3,5000,0,1,90 单元荷载,17,3,5,0000,1,90 单元荷载,18,3,5000,0,1,90单元荷载,19,3, 5000,0,1,90 单元荷载,20,3, 5000,0,1,90 单元荷载,21,3, 5000,0,1,90单元荷载,22,3, 5000,0,1,90 单元荷载,23,3, 5000,0,1,90 单元荷载,24,3, 5000,0,1,90单元荷载,25,3, 5000,0,1,90 单元荷载,1,5,0,2000,0,1,90 单元荷载,2,5,2000,4000,0,1,90单元荷载,3,5,4000,6000,0,1,90 单元荷载,4,5,6000,8000,0,1,90 单元荷载,5,5,8000,10000,0,1,90单元材料性质,1,15,EA2,EI2,PA2,0,-1 单元材料性质,16,25,EA1,EI1,PA1,0,-1自振频率参数,5,1,0.00000005 屈曲荷载参数,5,1,0.00000005内力图、动力特性图、稳定分析： N图（单位：KN） Q图（单位：KN） M图（单位：KN·m） 第一阶频率为12.17/s的 第二阶频率为39.59/s的 第三阶频率为74.52/s的 自由振动图 自由振动图 自由振动图 第四阶频率为116.97/s的 第五阶频率为159.26/s的 第一阶临界荷载系数为372.42 自由振动图 自由振动图 的失稳模态图第二阶临界荷载系数为586.17 第三阶临界荷载系数为827.42 第四阶临界荷载系数为1114.88 的失稳模态图 的失稳模态图 的失稳模态图 第五阶临界荷载系数为1187.66 稳定分析:的失稳模态图 因为，当第 1 阶失稳荷载系数 = 372.42时，结构已失稳，所以失稳荷载值为372.42。
