建筑工业学院计算结构力学实验报告

《计算结构力学》计算报告姓名:学号:班级:指导:安徽建筑工业学院土木工程学院实验一:形成单元定位向量的程序设计一、求图示结构各单元的单元定位向量。二、输入数据文件6,6,43,1,4,2,4,1,3,2,3,4,5,61,5,5,10001,2,6,6,10001,3,1,1,10001,4,0,1,10001三、计算结果：平面杆系分析单元数NE=6结点数NJ=6特殊结点数NJT=4单元两端结点号:单元起点终点单元起点终点单元起点终点单元起点终点131242341432534656特殊结点信息(结点号,结点各自由度的约束特征数):结点号特征数特征数特征数结点号特征数特征数特征数结点号特征数特征数特征数15510001266100013111000140110001未知量编号:23-156-100-110-1234567单元定位向量MW(6)单元编号:1MW=00-123-1单元编号:2MW=10-156-1单元编号:3MW=10-123-1单元编号:4MW=00-156-1单元编号:5MW=00-110-1单元编号:6MW=234567四、实验分析：1．单元⑥是梁式杆，其余均为轴力杆单元；2．单元定位向量与结点坐标和杆件长度、截面性质、材料都无关。五、实验日期：二〇一二年六月二十日星期三实验二:形成单元刚度矩阵和结构刚度矩阵的程序设计一、图示桁架腹杆EA1=186KN，EI1=368KNm2;弦杆EA2=213KN，EI2=559KNm2。求各单元的单元刚度矩阵和结构刚度矩阵。二、输入数据文件6,6,4,23,1,4,2,4,1,3,2,3,4,5,61,5,5,10001,2,6,6,10001,3,1,1,10001,4,0,1,100011,1,1,1,2,2186,213368,5590,6,0,6,0,65,5,0,0,5,5三、计算结果：平面杆系分析单元数NE=6结点数NJ=6特殊结点数NJT=4截面特性分组数NAI=2单元两端结点号:单元起点终点单元起点终点单元起点终点单元起点终点131242341432534656单元截面特性:111122特殊结点信息(结点号,结点各自由度的约束特征数):结点号特征数特征数特征数结点号特征数特征数特征数结点号特征数特征数特征数15510001266100013111000140110001截面特性EA:186.000213.000截面特性EI:368.000559.000结点坐标:结点X坐标Y坐标结点X坐标Y坐标结点X坐标Y坐标1.0005.00026.0005.0003.000.00046.000.0005.0005.00066.0005.000未知量编号:23-156-100-110-1234567单元定位向量MW(6)单元编号:1MW=00-123-1单元编号:2MW=10-156-1单元编号:3MW=10-123-1单元编号:4MW=00-156-1单元编号:5MW=00-110-1单元编号:6MW=234567单元刚度矩阵：单元号135.33.00-88.32-35.33.00-88.32.0037.20.00.00-37.20.00-88.32.00294.4088.32.00147.20-35.33.0088.3235.33.0088.32.00-37.20.00.0037.20.00-88.32.00147.2088.32.00294.40单元刚度矩阵：单元号235.33.00-88.32-35.33.00-88.32.0037.20.00.00-37.20.00-88.32.00294.4088.32.00147.20-35.33.0088.3235.33.0088.32.00-37.20.00.0037.20.00-88.32.00147.2088.32.00294.40单元刚度矩阵：单元号317.85-7.15-23.17-17.857.15-23.17-7.1515.23-27.817.15-15.23-27.81-23.17-27.81188.4723.1727.8194.24-17.857.1523.1717.85-7.1523.177.15-15.2327.81-7.1515.2327.81-23.17-27.8194.2423.1727.81188.47单元刚度矩阵：单元号417.857.15-23.17-17.85-7.15-23.177.1515.2327.81-7.15-15.2327.81-23.1727.81188.4723.17-27.8194.24-17.85-7.1523.1717.857.1523.17-7.15-15.23-27.817.1515.23-27.81-23.1727.8194.2423.17-27.81188.47单元刚度矩阵：单元号535.50.00.00-35.50.00.00.0031.0693.17.00-31.0693.17.0093.17372.67.00-93.17186.33-35.50.00.0035.50.00.00.00-31.06-93.17.0031.06-93.17.0093.17186.33.00-93.17372.67单元刚度矩阵：单元号635.50.00.00-35.50.00.00.0031.0693.17.00-31.0693.17.0093.17372.67.00-93.17186.33-35.50.00.0035.50.00.00.00-31.06-93.17.0031.06-93.17.0093.17186.33.00-93.17372.67结构刚度矩阵：88.7-17.97.2.0-35.3.0.0-17.988.7-7.2.0-35.5.0.07.2-7.283.593.2.0-31.193.2.0.093.2372.7.0-93.2186.3-35.3-35.5.0.088.77.2.0.0.0-31.1-93.27.283.5-93.2.0.093.2186.3.0-93.2372.7四、实验分析：1．单元刚度矩阵和结构刚度矩阵均与杆件的材料性质有关。2．单元刚度矩阵具有对称的性质；结构刚度矩阵不具有。五、实验日期：二〇一二年六月二十日星期三实验三:形成载荷列阵的程序设计一、图示桁架腹杆EA1=186KN，EI1=368KNm2;弦杆EA2=213KN，EI2=559KNm2。求结构的总载荷向量。二、输入数据文件6,6,4,2,2,13,1,4,2,4,1,3,2,3,4,5,61,5,5,10001,2,6,6,10001,3,1,1,10001,4,0,1,100011,1,1,1,2,2186,213368,5590,6,0,6,0,65,5,0,0,5,51.1,-5,2.1,56.1,6,-3三、计算结果：平面杆系分析单元数NE=6结点数NJ=6特殊结点数NJT=4支座结点数NJZ=2截面特性分组数NAI=2单元两端结点号:单元起点终点单元起点终点单元起点终点单元起点终点131242341432534656单元截面特性:111122特殊结点信息(结点号,结点各自由度的约束特征数):结点号特征数特征数特征数结点号特征数特征数特征数结点号特征数特征数特征数15510001266100013111000140110001截面特性EA:186.000213.000截面特性EI:368.000559.000结点坐标:结点X坐标Y坐标结点X坐标Y坐标结点X坐标Y坐标1.0005.00026.0005.0003.000.00046.000.0005.0005.00066.0005.000未知量编号:23-156-100-110-1234567单元定位向量MW(6)单元编号:1MW=00-123-1单元编号:2MW=10-156-1单元编号:3MW=10-123-1单元编号:4MW=00-156-1单元编号:5MW=00-110-1单元编号:6MW=234567单元刚度矩阵：单元号135.33.00-88.32-35.33.00-88.32.0037.20.00.00-37.20.00-88.32.00294.4088.32.00147.20-35.33.0088.3235.33.0088.32.00-37.20.00.0037.20.00-88.32.00147.2088.32.00294.40单元刚度矩阵：单元号235.33.00-88.32-35.33.00-88.32.0037.20.00.00-37.20.00-88.32.00294.4088.32.00147.20-35.33.0088.3235.33.0088.32.00-37.20.00.0037.20.00-88.32.00147.2088.32.00294.40单元刚度矩阵：单元号317.85-7.15-23.17-17.857.15-23.17-7.1515.23-27.817.15-15.23-27.81-23.17-27.81188.4723.1727.8194.24-17.857.1523.1717.85-7.1523.177.15-15.2327.81-7.1515.2327.81-23.17-27.8194.2423.1727.81188.47单元刚度矩阵：单元号417.857.15-23.17-17.85-7.15-23.177.1515.2327.81-7.15-15.2327.81-23.1727.81188.4723.17-27.8194.24-17.85-7.1523.1717.857.1523.17-7.15-15.23-27.817.1515.23-27.81-23.1727.8194.2423.17-27.81188.47单元刚度矩阵：单元号535.50.00.00-35.50.00.00.0031.0693.17.00-31.0693.17.0093.17372.67.00-93.17186.33-35.50.00.0035.50.00.00.00-31.06-93.17.0031.06-93.17.0093.17186.33.00-93.17372.67单元刚度矩阵：单元号635.50.00.00-35.50.00.00.0031.0693.17.00-31.0693.17.0093.17372.67.00-93.17186.33-35.50.00.0035.50.00.00.00-31.06-93.17.0031.06-93.17.0093.17186.33.00-93.17372.67结构刚度矩阵：88.7-17.97.2.0-35.3.0.0-17.988.7-7.2.0-35.5.0.07.2-7.283.593.2.0-31.193.2.0.093.2372.7.0-93.2186.3-35.3-35.5.0.088.77.2.0.0.0-31.1-93.27.283.5-93.2.0.093.2186.3.0-93.2372.7结点荷载数单元荷载数21结点荷载:1.100-5.0002.1005.000单元荷载:6.1006.000-3.000荷载向量:.000-5.000-9.000-9.0005.000-9.0009.000四、实验分析：1．单元载荷作用在梁式杆上。2．同一结点上的桁架杆与梁式杆有不同的截面，须分开标注。五、实验日期：二〇一二年六月二十日星期三实验四:高斯消元法的程序设计一、图示桁架腹杆EA1=186KN，EI1=368KNm2;弦杆EA2=213KN，EI2=559KNm2。求结构的结点位移。二、输入数据文件6,6,4,2,2,13,1,4,2,4,1,3,2,3,4,5,61,5,5,10001,2,6,6,10001,3,1,1,10001,4,0,1,100011,1,1,1,2,2186,213368,5590,6,0,6,0,65,5,0,0,5,51.1,-5,2.1,56.1,6,-3三、计算结果：平面杆系分析单元数NE=6结点数NJ=6特殊结点数NJT=4支座结点数NJZ=2截面特性分组数NAI=2单元两端结点号:单元起点终点单元起点终点单元起点终点单元起点终点131242341432534656单元截面特性:111122特殊结点信息(结点号,结点各自由度的约束特征数):结点号特征数特征数特征数结点号特征数特征数特征数结点号特征数特征数特征数15510001266100013111000140110001截面特性EA:186.000213.000截面特性EI:368.000559.000结点坐标:结点X坐标Y坐标结点X坐标Y坐标结点X坐标Y坐标1.0005.00026.0005.0003.000.00046.000.0005.0005.00066.0005.000未知量编号:23-156-100-110-1234567单元定位向量MW(6)单元编号:1MW=00-123-1单元编号:2MW=10-156-1单元编号:3MW=10-123-1单元编号:4MW=00-156-1单元编号:5MW=00-110-1单元编号:6MW=234567单元刚度矩阵：单元号135.33.00-88.32-35.33.00-88.32.0037.20.00.00-37.20.00-88.32.00294.4088.32.00147.20-35.33.0088.3235.33.0088.32.00-37.20.00.0037.20.00-88.32.00147.2088.32.00294.40单元刚度矩阵：单元号235.33.00-88.32-35.33.00-88.32.0037.20.00.00-37.20.00-88.32.00294.4088.32.00147.20-35.33.0088.3235.33.0088.32.00-37.20.00.0037.20.00-88.32.00147.2088.32.00294.40单元刚度矩阵：单元号317.85-7.15-23.17-17.857.15-23.17-7.1515.23-27.817.15-15.23-27.81-23.17-27.81188.4723.1727.8194.24-17.857.1523.1717.85-7.1523.177.15-15.2327.81-7.1515.2327.81-23.17-27.8194.2423.1727.81188.47单元刚度矩阵：单元号417.857.15-23.17-17.85-7.15-23.177.1515.2327.81-7.15-15.2327.81-23.1727.81188.4723.17-27.8194.24-17.85-7.1523.1717.857.1523.17-7.15-15.23-27.817.1515.23-27.81-23.1727.8194.2423.17-27.81188.47单元刚度矩阵：单元号535.50.00.00-35.50.00.00.0031.0693.17.00-31.0693.17.0093.17372.67.00-93.17186.33-35.50.00.0035.50.00.00.00-31.06-93.17.0031.06-93.17.0093.17186.33.00-93.17372.67单元刚度矩阵：单元号635.50.00.00-35.50.00.00.0031.0693.17.00-31.0693.17.0093.17372.67.00-93.17186.33-35.50.00.0035.50.00.00.00-31.06-93.17.0031.06-93.17.0093.17186.33.00-93.17372.67结构刚度矩阵：88.7-17.97.2.0-35.3.0.0-17.988.7-7.2.0-35.5.0.07.2-7.283.593.2.0-31.193.2.0.093.2372.7.0-93.2186.3-35.3-35.5.0.088.77.2.0.0.0-31.1-93.27.283.5-93.2.0.093.2186.3.0-93.2372.7结点荷载数单元荷载数21结点荷载:1.100-5.0002.1005.000单元荷载:6.1006.000-3.000荷载向量:.000-5.000-9.000-9.0005.000-9.0009.000结构位移列阵:.03587846-.03529678-.18136728-.04830054.07117523-.18136728.04830054单元位移列阵:MUAVAZTAUBVBZTB1.0000000.0000000.0000000-.0352968-.1813673.00000002.0358785.0000000.0000000.0711752-.1813673.00000003.0358785.0000000.0000000-.0352968-.1813673.00000004.0000000.0000000.0000000.0711752-.1813673.00000005.0000000.0000000.0000000.0358785.0000000.00000006-.0352968-.1813673-.0483005.0711752-.1813673.0483005四、实验分析：同一结点上的桁架杆与梁式杆截面有相同的水平位移和竖直位移。五、实验日期：二〇一二年六月二十日星期三实验五:计算杆端力和支座反力的程序设计一、图示桁架腹杆EA1=186KN，EI1=368KNm2;弦杆EA2=213KN，EI2=559KNm2。求结构的杆件杆端力和支座反力。二、输入数据文件6,6,4,2,2,13,1,4,2,4,1,3,2,3,4,5,61,5,5,10001,2,6,6,10001,3,1,1,10001,4,0,1,100011,1,1,1,2,2186,213368,5590,6,0,6,0,65,5,0,0,5,51.1,-5,2.1,56.1,6,-3计算结果：平面杆系分析单元数NE=6结点数NJ=6特殊结点数NJT=4支座结点数NJZ=2截面特性分组数NAI=2单元两端结点号:单元起点终点单元起点终点单元起点终点单元起点终点131242341432534656单元截面特性:111122特殊结点信息(结点号,结点各自由度的约束特征数):结点号特征数特征数特征数结点号特征数特征数特征数结点号特征数特征数特征数15510001266100013111000140110001截面特性EA:186.000213.000截面特性EI:368.000559.000结点坐标:结点X坐标Y坐标结点X坐标Y坐标结点X坐标Y坐标1.0005.00026.0005.0003.000.00046.000.0005.0005.00066.0005.000未知量编号:23-156-100-110-1234567单元定位向量MW(6)单元编号:1MW=00-123-1单元编号:2MW=10-156-1单元编号:3MW=10-123-1单元编号:4MW=00-156-1单元编号:5MW=00-110-1单元编号:6MW=234567单元刚度矩阵：单元号135.33.00-88.32-35.33.00-88.32.0037.20.00.00-37.20.00-88.32.00294.4088.32.00147.20-35.33.0088.3235.33.0088.32.00-37.20.00.0037.20.00-88.32.00147.2088.32.00294.40单元刚度矩阵：单元号235.33.00-88.32-35.33.00-88.32.0037.20.00.00-37.20.00-88.32.00294.4088.32.00147.20-35.33.0088.3235.33.0088.32.00-37.20.00.0037.20.00-88.32.00147.2088.32.00294.40单元刚度矩阵：单元号317.85-7.15-23.17-17.857.15-23.17-7.1515.23-27.817.15-15.23-27.81-23.17-27.81188.4723.1727.8194.24-17.857.1523.1717.85-7.1523.177.15-15.2327.81-7.1515.2327.81-23.17-27.8194.2423.1727.81188.47单元刚度矩阵：单元号417.857.15-23.17-17.85-7.15-23.177.1515.2327.81-7.15-15.2327.81-23.1727.81188.4723.17-27.8194.24-17.85-7.1523.1717.857.1523.17-7.15-15.23-27.817.1515.23-27.81-23.1727.8194.2423.17-27.81188.47单元刚度矩阵：单元号535.50.00.00-35.50.00.00.0031.0693.17.00-31.0693.17.0093.17372.67.00-93.17186.33-35.50.00.0035.50.00.00.00-31.06-93.17.0031.06-93.17.0093.17186.33.00-93.17372.67单元刚度矩阵：单元号635.50.00.00-35.50.00.00.0031.0693.17.00-31.0693.17.0093.17372.67.00-93.17186.33-35.50.00.0035.50.00.00.00-31.06-93.17.0031.06-93.17.0093.17186.33.00-93.1