《计算结构力学》
计算报告
姓名: 周剑敏
学号: 09201010628
班级: 09土木6班
指导: 苏少卿
安徽建筑工业学院土木工程学院
实验一: 形成单元定位向量的程序设计
一、求图示结构各单元的单元定位向量。 76.008
56.006
36.004
二、输入数据文件
9 8 2
1 3 3 5 5 7 2 4 4 6 6 8 7 8 1 1 1 1 2 1 1 1
5 6 3 4 2
三、计算结果：
平面杆系分析
单元数NE= 9
结点数NJ= 8
特殊结点数NJT= 2
单元两端结点号:
单元 起点 终点 单元 起点 终点 单元 起点 终点 单元 起点 终点 1 1 3 2 3 5 3 5 7 4 2 4 5 4 6 6 6 8 7 7 8 8 5 6 9 3 4
特殊结点信息(结点号,结点各自由度的约束特征数):
结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 结点号 特征数 特征数 1 1 1 1 2 1 1 1
未知量编号:
0 0 0 0 0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
单元定位向量MW(6)
单元编号: 1 MW = 0 0 0 1 2 3
单元编号: 2 MW = 1 2 3 7 8 9
单元编号: 3 MW = 7 8 9 13 14 15
单元编号: 4 MW = 0 0 0 4 5 6
单元编号: 5 MW = 4 5 6 10 11 12
单元编号: 6 MW = 10 11 12 16 17 18
单元编号: 7 MW = 13 14 15 16 17 18
单元编号: 8 MW = 7 8 9 10 11 12
单元编号: 9 MW = 1 2 3 4 5 6
四、实验日期：二〇一二年六月二十二日星期五
特征数 3
实验二: 形成单元刚度矩阵和结构刚度矩阵的程序设计
一、图示结构竖杆的EA=468000KN,EI=14040KNm2,其余杆件EA=390000KN,EI=8125KNm2。求结构的总载荷向量。 76.008
56.006
36.004
二、输入数据文件
9 8 2 2
1 3 3 5 5 7 2 4 4 6 6 8 7 8
1 1 1 1 2 1 1 1
1 1 1 2 1 1 1 2 2
468000 390000
14040 8125
0 6 0 6 0 6 0 6
0 0 4 4 7 7 10 10
三、计算结果：
5 6 3 4 4
平面杆系分析
单元数NE= 9
结点数NJ= 8
特殊结点数NJT= 2
截面特性分组数NAI= 2
单元两端结点号:
单元 起点 终点 单元 起点 终点 单元 起点 终点 单元 起点 终点 1 1 3 2 3 5 3 5 7 4 2 4 5 4 6 6 6 8 7 7 8 8 5 6 9 3 4
单元截面特性:
1 1 1 2 1 1 1 2 2
特殊结点信息(结点号,结点各自由度的约束特征数):
结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 结点号 特征数 特征数 1 1 1 1 2 1 1 1
截面特性EA:
468000.000 390000.000
截面特性EI:
14040.000 8125.000
结点坐标:
结点 X坐标 Y坐标 结点 X坐标 Y坐标 结点 X坐标 Y坐标 1 .000 .000 2 6.000 .000 3 .000 4.000
4 6.000 4.000 5 .000 7.000 6 6.000 7.000
7 .000 10.000 8 6.000 10.000
未知量编号:
0 0 0 0 0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
单元定位向量MW(6)
单元编号: 1 MW = 0 0 0 1 2 3
单元编号: 2 MW = 1 2 3 7 8 9
单元编号: 3 MW = 7 8 9 13 14 15
单元编号: 4 MW = 0 0 0 4 5 6
单元编号: 5 MW = 4 5 6 10 11 12
单元编号: 6 MW = 10 11 12 16 17 18
单元编号: 7 MW = 13 14 15 16 17 18
单元编号: 8 MW = 7 8 9 10 11 12
单元编号: 9 MW = 1 2 3 4 5 6
单元刚度矩阵：单元号 1
2632.50 .00 -5265.00 -2632.50 .00 -5265.00
特征数 5
.00 117000.00 .00 .00 -117000.00 .00 -5265.00 .00 14040.00 5265.00 .00 7020.00 -2632.50 .00 5265.00 2632.50 .00 5265.00 .00 -117000.00 .00 .00 117000.00 .00 -5265.00 .00 7020.00 5265.00 .00 14040.00
单元刚度矩阵：单元号 2
6240.00 .00 -9360.00 -6240.00 .00 -9360.00 .00 156000.00 .00 .00 -156000.00 .00 -9360.00 .00 18720.00 -6240.00 .00 9360.00 .00 -156000.00 .00 -9360.00 .00 9360.00
单元刚度矩阵：单元号 3
6240.00 .00 -9360.00 .00 156000.00 .00 -9360.00 .00 18720.00 -6240.00 .00 9360.00 .00 -156000.00 .00 -9360.00 .00 9360.00
单元刚度矩阵：单元号 4
1523.44 .00 -3046.88 .00 97500.00 .00 -3046.88 .00 8125.00 -1523.44 .00 3046.88 .00 -97500.00 .00 -3046.88 .00 4062.50
单元刚度矩阵：单元号 5
6240.00 .00 -9360.00 .00 156000.00 .00 -9360.00 .00 18720.00 -6240.00 .00 9360.00 .00 -156000.00 .00 -9360.00 .00 9360.00
单元刚度矩阵：单元号 6
6240.00 .00 -9360.00 .00 156000.00 .00 -9360.00 .00 18720.00 -6240.00 .00 9360.00 .00 -156000.00 .00
9360.00 6240.00 .00 9360.00 -6240.00 .00 9360.00 6240.00 .00 9360.00 -1523.44 .00 3046.88 1523.44 .00 3046.88 -6240.00 .00 9360.00 6240.00 .00 9360.00 -6240.00 .00 9360.00 6240.00 .00 .00 .00 156000.00 .00 .00 -156000.00 .00 .00 156000.00 .00 .00 -97500.00 .00 .00 97500.00 .00 .00 -156000.00 .00 .00 156000.00 .00 .00 -156000.00 .00 .00 156000.00 9360.00 9360.00 .00 18720.00 -9360.00 .00 9360.00 9360.00 .00 18720.00 -3046.88 .00 4062.50 3046.88 .00 8125.00 -9360.00 .00 9360.00 9360.00 .00 18720.00 -9360.00 .00 9360.00 9360.00 .00
6
-9360.00 .00 9360.00 9360.00 .00 18720.00
单元刚度矩阵：单元号 7
78000.00 .00 .00 -78000.00 .00 .00 .00 780.00 2340.00 .00 -780.00 2340.00 .00 2340.00 9360.00 .00 -2340.00 4680.00 -78000.00 .00 .00 78000.00 .00 .00 .00 -780.00 -2340.00 .00 780.00 -2340.00 .00 2340.00 4680.00 .00 -2340.00 9360.00
单元刚度矩阵：单元号 8
65000.00 .00 .00 .00 451.39 1354.17 .00 1354.17 5416.67 -65000.00 .00 .00 .00 -451.39 -1354.17 .00 1354.17 2708.33
单元刚度矩阵：单元号 9
65000.00 .00 .00 .00 451.39 1354.17 .00 1354.17 5416.67 -65000.00 .00 .00 .00 -451.39 -1354.17 .00 1354.17 2708.33
结构刚度矩阵：
73872.5 .0 -4095.0 -6240.0 .0 -9360.0 .0 .0 .0 .0 273451.4 1354.2 .0 -156000.0 .0 .0 .0 .0 -4095.0 1354.2 38176.7 9360.0 .0 9360.0 .0 .0 .0 -65000.0 .0 .0 .0 .0 .0 .0 .0 .0 .0 -451.4 -1354.2 .0 .0 .0 .0 .0 .0 .0 1354.2 2708.3 .0 .0 .0
-65000.00 .00 .00 65000.00 .00 .00 -65000.00 .00 .00 65000.00 .00 .00 -65000.0 .0 .0 .0 .0 .0 .0 .0 .0 72763.4 -6240.0 .0 .0 .0 .0 -6313.1 9360.0 .00 -451.39 -1354.17 .00 451.39 -1354.17 .00 -451.39 -1354.17 .00 451.39 -1354.17 .0 .0 .0 -451.4 .0 .0 -1354.2 .0 .0 .0 .0 .0 253951.4 -156000.0 .0 -1354.2 .0 .00 1354.17 2708.33 .00 -1354.17 5416.67 .00 1354.17 2708.33 .00 -1354.17 5416.67 .0 .0 .0 1354.2 .0 .0 2708.3 .0 .0 -6313.1 -9360.0 .0 -1354.2 .0 .0 32261.7 9360.0
7
.0 .0 .0 .0 .0 .0 -6240.0 .0 9360.0 .0 .0 .0 77480.0 .0 .0 -65000.0 .0 .0 -6240.0 .0 -9360.0 .0 .0 .0 .0 -156000.0 .0 .0 .0 .0 .0 312451.4 1354.2 .0 -451.4 1354.2 .0 -156000.0 .0 .0 .0 .0 -9360.0 .0 9360.0 .0 .0 .0 .0 1354.2 42856.7 .0 -1354.2 2708.3 9360.0 .0 9360.0 .0 .0
.0 .0 .0 -6240.0 .0
-65000.0 .0 .0 77480.0 .0
.0 .0 .0 -6240.0 .0
.0 .0 .0 .0 -156000.0
.0 -451.4 -1354.2 .0 312451.4
.0 .0 .0 .0 -156000.0
.0 .0 .0 -9360.0 .0
.0 1354.2 2708.3 .0 -1354.2
.0 .0 .0 9360.0 .0
.0 .0 .0 .0 .0
-6240.0 .0 9360.0 .0 .0
84240.0 .0 9360.0 -78000.0 .0
.0 .0 .0 .0 .0
.0 -156000.0 .0 .0 .0
.0 156780.0 2340.0 .0 -780.0
.0 .0 .0 .0 .0
-9360.0 .0 9360.0 .0 .0
9360.0 2340.0 28080.0 .0 -2340.0
.0 .0 .0 .0 .0
.0 .0 .0 -6240.0 .0
-78000.0 .0 .0 84240.0 .0
.0 .0 .0 .0 .0
.0 .0 .0 .0 -156000.0
.0 -780.0 -2340.0 .0 156780.0
.0 .0 .0 .0 .0
.0 .0 .0 -9360.0 .0
.0 2340.0 4680.0 9360.0 -2340.0
四、实验日期：二〇一二年六月二十二日星期五
.0 9360.0 .0 -9360.0 .0 -1354.2 .0 9360.0 42856.7 9360.0 .0 .0 .0 .0 .0 2340.0 .0 .0 4680.0 .0 9360.0 9360.0 .0 .0 -2340.0 .0 9360.0 28080.0 8
实验三: 形成载荷列阵的程序设计
一、图示结构竖杆的EA=468000KN,EI=14040KNm2,其余杆件
EA=390000KN,EI=8125KNm2。求结构的总载荷向量。 76.008
56.006
36.004
二、输入数据文件
9 8 2 2 2 0 3
1 3 3 5 5 7 2 4 4 6 6 8 7 8
1 1 1 1 2 1 1 1
1 1 1 2 1 1 1 2 2
468000 390000
14040 8125
0 6 0 6 0 6 0 6
0 0 4 4 7 7 10 10
4 .1,6,-3,8.1,6,-3,9.1,6,-3
三、计算结果：
5 6 3 4 9
平面杆系分析
单元数NE= 9
结点数NJ= 8
特殊结点数NJT= 2
支座结点数NJZ= 2
截面特性分组数NAI= 2
单元两端结点号:
单元 起点 终点 单元 起点 终点 单元 起点 终点 单元 起点 终点 1 1 3 2 3 5 3 5 7 4 2 4 5 4 6 6 6 8 7 7 8 8 5 6 9 3 4
单元截面特性:
1 1 1 2 1 1 1 2 2
特殊结点信息(结点号,结点各自由度的约束特征数):
结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 结点号 特征数 特征数 1 1 1 1 2 1 1 1
截面特性EA:
468000.000 390000.000
截面特性EI:
14040.000 8125.000
结点坐标:
结点 X坐标 Y坐标 结点 X坐标 Y坐标 结点 X坐标 Y坐标 1 .000 .000 2 6.000 .000 3 .000 4.000
4 6.000 4.000 5 .000 7.000 6 6.000 7.000
7 .000 10.000 8 6.000 10.000
未知量编号:
0 0 0 0 0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
单元定位向量MW(6)
单元编号: 1 MW = 0 0 0 1 2 3
单元编号: 2 MW = 1 2 3 7 8 9
单元编号: 3 MW = 7 8 9 13 14 15
单元编号: 4 MW = 0 0 0 4 5 6
特征数 10
单元编号: 5 MW = 4 5 6 10 11 12
单元编号: 单元编号: 单元编号: 单元编号:
6 7 8 9
MW = 10 11 12 16 17 18 MW = 13 14 15 16 17 18 MW = 7 8 9 10 11 12 MW = 1 2 3 4 5 6
单元刚度矩阵：单元号 1
2632.50 .00 -5265.00 -2632.50 .00 -5265.00 .00 117000.00 .00 .00 -117000.00 .00 -5265.00 .00 14040.00 -2632.50 .00 5265.00 .00 -117000.00 .00 -5265.00 .00 7020.00
单元刚度矩阵：单元号 2
6240.00 .00 -9360.00 .00 156000.00 .00 -9360.00 .00 18720.00 -6240.00 .00 9360.00 .00 -156000.00 .00 -9360.00 .00 9360.00
单元刚度矩阵：单元号 3
6240.00 .00 -9360.00 .00 156000.00 .00 -9360.00 .00 18720.00 -6240.00 .00 9360.00 .00 -156000.00 .00 -9360.00 .00 9360.00
单元刚度矩阵：单元号 4
1523.44 .00 -3046.88 .00 97500.00 .00 -3046.88 .00 8125.00 -1523.44 .00 3046.88 .00 -97500.00 .00 -3046.88 .00 4062.50
单元刚度矩阵：单元号 5
6240.00 .00 -9360.00 .00 156000.00 .00 -9360.00 .00 18720.00 -6240.00 .00 9360.00 .00 -156000.00 .00
5265.00 2632.50 .00 5265.00 -6240.00 .00 9360.00 6240.00 .00 9360.00 -6240.00 .00 9360.00 6240.00 .00 9360.00 -1523.44 .00 3046.88 1523.44 .00 3046.88 -6240.00 .00 9360.00 6240.00 .00 .00 .00 117000.00 .00 .00 -156000.00 .00 .00 156000.00 .00 .00 -156000.00 .00 .00 156000.00 .00 .00 -97500.00 .00 .00 97500.00 .00 .00 -156000.00 .00 .00 156000.00 7020.00 5265.00 .00 14040.00 -9360.00 .00 9360.00 9360.00 .00 18720.00 -9360.00 .00 9360.00 9360.00 .00 18720.00 -3046.88 .00 4062.50 3046.88 .00 8125.00 -9360.00 .00 9360.00 9360.00 .00
11
-9360.00 .00 9360.00 9360.00 .00 18720.00
单元刚度矩阵：单元号 6
6240.00 .00 -9360.00 -6240.00 .00 -9360.00 .00 156000.00 .00 .00 -156000.00 .00 -9360.00 .00 18720.00 9360.00 .00 9360.00 -6240.00 .00 9360.00 6240.00 .00 9360.00 .00 -156000.00 .00 .00 156000.00 .00 -9360.00 .00 9360.00 9360.00 .00 18720.00
单元刚度矩阵：单元号 7
78000.00 .00 .00 .00 780.00 2340.00 .00 2340.00 9360.00 -78000.00 .00 .00 .00 -780.00 -2340.00 .00 2340.00 4680.00
单元刚度矩阵：单元号 8
65000.00 .00 .00 .00 451.39 1354.17 .00 1354.17 5416.67 -65000.00 .00 .00 .00 -451.39 -1354.17 .00 1354.17 2708.33
单元刚度矩阵：单元号 9
65000.00 .00 .00 .00 451.39 1354.17 .00 1354.17 5416.67 -65000.00 .00 .00 .00 -451.39 -1354.17 .00 1354.17 2708.33
结构刚度矩阵：
73872.5 .0 -4095.0 -6240.0 .0 -9360.0 .0 .0 .0 .0 273451.4 1354.2 .0 -156000.0 .0 .0 .0 .0 -4095.0 1354.2 38176.7 9360.0 .0 9360.0 .0 .0 .0
-78000.00 .00 .00 78000.00 .00 .00 -65000.00 .00 .00 65000.00 .00 .00 -65000.00 .00 .00 65000.00 .00 .00 -65000.0 .0 .0 .0 .0 .0 .0 .0 .0 .00 -780.00 -2340.00 .00 780.00 -2340.00 .00 -451.39 -1354.17 .00 451.39 -1354.17 .00 -451.39 -1354.17 .00 451.39 -1354.17 .0 .0 .0 -451.4 .0 .0 -1354.2 .0 .0 .00 2340.00 4680.00 .00 -2340.00 9360.00 .00 1354.17 2708.33 .00 -1354.17 5416.67 .00 1354.17 2708.33 .00 -1354.17 5416.67 .0 .0 .0 1354.2 .0 .0 2708.3 .0 .0
12
-65000.0 .0 .0 72763.4 .0 -6313.1
.0 .0 .0 -6240.0 .0 -9360.0 .0 .0 .0 .0 .0 .0 .0 -451.4 -1354.2 .0 253951.4 -1354.2 .0 .0 .0 .0 -156000.0 .0
.0 .0 .0 .0 .0 .0 .0 1354.2 2708.3 -6313.1 -1354.2 32261.7 .0 .0 .0 9360.0 .0 9360.0 .0 .0 .0 .0 .0 .0 -6240.0 77480.0 -6240.0 .0 .0 .0 -9360.0 .0 9360.0 .0 -65000.0 .0 .0 .0 .0
.0 .0 .0 .0 -6240.0 84240.0 .0 .0 .0 .0 -9360.0 9360.0 .0 .0 -78000.0 .0 .0 .0 .0 .0
.0 .0 .0 -156000.0 312451.4 -156000.0 .0 1354.2 .0 .0 .0 .0 .0 -451.4 .0 .0 1354.2 .0 .0 .0 .0 .0 -156000.0 156780.0 .0 .0 2340.0 .0 .0 .0 .0 .0 -780.0 .0 .0 9360.0 .0 -9360.0 .0 1354.2 .0 9360.0 42856.7 9360.0 .0 .0 .0 .0 -1354.2 .0 .0 2708.3 .0 .0 9360.0 9360.0 .0 .0 2340.0 .0 9360.0 28080.0 .0 .0 .0 .0 .0 -2340.0 .0 .0 .0 -65000.0 .0 .0 .0 .0 .0 .0 .0 -6240.0 77480.0 -6240.0 .0 .0 .0 -9360.0 .0 9360.0 .0 .0 -78000.0 .0 .0 .0 .0 .0 .0 .0 -6240.0 84240.0 .0 .0 .0 .0 -9360.0 .0 .0 .0 .0 -451.4 .0 .0 -1354.2 .0 .0 .0 .0 -156000.0 312451.4 -156000.0 .0 -1354.2 .0 .0 .0 .0 .0 .0 -780.0 .0 .0 -2340.0 .0 .0 .0 .0 -156000.0 156780.0 .0 .0 .0 .0 .0 .0 1354.2 .0 .0 2708.3 .0 9360.0 .0 -9360.0 .0 -1354.2 .0 9360.0
42856.7 9360.0 .0 .0 .0 .0 .0 2340.0 .0 .0 4680.0 .0 9360.0 9360.0 .0 .0 -2340.0 .0 9360.0
13
.0 2340.0 4680.0 9360.0 -2340.0 28080.0
结点荷载数 单元荷载数
2 0
结点荷载:
4.000 .100 6.000 -3.000
荷载向量:
.000 .000 .100 .000
.000 .000 -3.000 .000
.000 .000 .000 .000
四、实验日期：二〇一二年六月二十二日星期五
.000 .000 .000 .000 .000 .000 14
实验四: 高斯消元法的程序设计
一、图示结构竖杆的EA=468000KN,EI=14040KNm2,其余杆件
EA=390000KN,EI=8125KNm2。求结构的总载荷向量。 76.008
56.006
36.004
二、输入数据文件
9 8 2 2 2 0 3
1 3 3 5 5 7 2 4 4 6 6 8 7 8 5 6 3 4 1 1 1 1 2 1 1 1
1 1 1 2 1 1 1 2 2
468000 390000
14040 8125
0 6 0 6 0 6 0 6
0 0 4 4 7 7 10 10
4 .1,6,-3,8.1,6,-3,9.1,6,-3
三、计算结果：
15
平面杆系分析
单元数NE= 9
结点数NJ= 8
特殊结点数NJT= 2
支座结点数NJZ= 2
截面特性分组数NAI= 2
单元两端结点号:
单元 起点 终点 单元 起点 终点 单元 起点 终点 单元 起点 终点 1 1 3 2 3 5 3 5 7 4 2 4 5 4 6 6 6 8 7 7 8 8 5 6 9 3 4
单元截面特性:
1 1 1 2 1 1 1 2 2
特殊结点信息(结点号,结点各自由度的约束特征数):
结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 结点号 特征数 特征数 1 1 1 1 2 1 1 1
截面特性EA:
468000.000 390000.000
截面特性EI:
14040.000 8125.000
结点坐标:
结点 X坐标 Y坐标 结点 X坐标 Y坐标 结点 X坐标 Y坐标 1 .000 .000 2 6.000 .000 3 .000 4.000
4 6.000 4.000 5 .000 7.000 6 6.000 7.000
7 .000 10.000 8 6.000 10.000
未知量编号:
0 0 0 0 0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
单元定位向量MW(6)
单元编号: 1 MW = 0 0 0 1 2 3
单元编号: 2 MW = 1 2 3 7 8 9
单元编号: 3 MW = 7 8 9 13 14 15
单元编号: 4 MW = 0 0 0 4 5 6
特征数 16
单元编号: 5 MW = 4 5 6 10 11 12
单元编号: 单元编号: 单元编号: 单元编号:
6 7 8 9
MW = 10 11 12 16 17 18 MW = 13 14 15 16 17 18 MW = 7 8 9 10 11 12 MW = 1 2 3 4 5 6
单元刚度矩阵：单元号 1
2632.50 .00 -5265.00 -2632.50 .00 -5265.00 .00 117000.00 .00 .00 -117000.00 .00 -5265.00 .00 14040.00 -2632.50 .00 5265.00 .00 -117000.00 .00 -5265.00 .00 7020.00
单元刚度矩阵：单元号 2
6240.00 .00 -9360.00 .00 156000.00 .00 -9360.00 .00 18720.00 -6240.00 .00 9360.00 .00 -156000.00 .00 -9360.00 .00 9360.00
单元刚度矩阵：单元号 3
6240.00 .00 -9360.00 .00 156000.00 .00 -9360.00 .00 18720.00 -6240.00 .00 9360.00 .00 -156000.00 .00 -9360.00 .00 9360.00
单元刚度矩阵：单元号 4
1523.44 .00 -3046.88 .00 97500.00 .00 -3046.88 .00 8125.00 -1523.44 .00 3046.88 .00 -97500.00 .00 -3046.88 .00 4062.50
单元刚度矩阵：单元号 5
6240.00 .00 -9360.00 .00 156000.00 .00 -9360.00 .00 18720.00 -6240.00 .00 9360.00 .00 -156000.00 .00
5265.00 2632.50 .00 5265.00 -6240.00 .00 9360.00 6240.00 .00 9360.00 -6240.00 .00 9360.00 6240.00 .00 9360.00 -1523.44 .00 3046.88 1523.44 .00 3046.88 -6240.00 .00 9360.00 6240.00 .00 .00 .00 117000.00 .00 .00 -156000.00 .00 .00 156000.00 .00 .00 -156000.00 .00 .00 156000.00 .00 .00 -97500.00 .00 .00 97500.00 .00 .00 -156000.00 .00 .00 156000.00 7020.00 5265.00 .00 14040.00 -9360.00 .00 9360.00 9360.00 .00 18720.00 -9360.00 .00 9360.00 9360.00 .00 18720.00 -3046.88 .00 4062.50 3046.88 .00 8125.00 -9360.00 .00 9360.00 9360.00 .00
17
-9360.00 .00 9360.00 9360.00 .00 18720.00
单元刚度矩阵：单元号 6
6240.00 .00 -9360.00 -6240.00 .00 -9360.00 .00 156000.00 .00 .00 -156000.00 .00 -9360.00 .00 18720.00 9360.00 .00 9360.00 -6240.00 .00 9360.00 6240.00 .00 9360.00 .00 -156000.00 .00 .00 156000.00 .00 -9360.00 .00 9360.00 9360.00 .00 18720.00
单元刚度矩阵：单元号 7
78000.00 .00 .00 .00 780.00 2340.00 .00 2340.00 9360.00 -78000.00 .00 .00 .00 -780.00 -2340.00 .00 2340.00 4680.00
单元刚度矩阵：单元号 8
65000.00 .00 .00 .00 451.39 1354.17 .00 1354.17 5416.67 -65000.00 .00 .00 .00 -451.39 -1354.17 .00 1354.17 2708.33
单元刚度矩阵：单元号 9
65000.00 .00 .00 .00 451.39 1354.17 .00 1354.17 5416.67 -65000.00 .00 .00 .00 -451.39 -1354.17 .00 1354.17 2708.33
结构刚度矩阵：
73872.5 .0 -4095.0 -6240.0 .0 -9360.0 .0 .0 .0 .0 273451.4 1354.2 .0 -156000.0 .0 .0 .0 .0 -4095.0 1354.2 38176.7 9360.0 .0 9360.0 .0 .0 .0
-78000.00 .00 .00 78000.00 .00 .00 -65000.00 .00 .00 65000.00 .00 .00 -65000.00 .00 .00 65000.00 .00 .00 -65000.0 .0 .0 .0 .0 .0 .0 .0 .0 .00 -780.00 -2340.00 .00 780.00 -2340.00 .00 -451.39 -1354.17 .00 451.39 -1354.17 .00 -451.39 -1354.17 .00 451.39 -1354.17 .0 .0 .0 -451.4 .0 .0 -1354.2 .0 .0 .00 2340.00 4680.00 .00 -2340.00 9360.00 .00 1354.17 2708.33 .00 -1354.17 5416.67 .00 1354.17 2708.33 .00 -1354.17 5416.67 .0 .0 .0 1354.2 .0 .0 2708.3 .0 .0
18
-65000.0 .0 .0 72763.4 .0 -6313.1
.0 .0 .0 -6240.0 .0 -9360.0 .0 .0 .0 .0 .0 .0 .0 -451.4 -1354.2 .0 253951.4 -1354.2 .0 .0 .0 .0 -156000.0 .0
.0 .0 .0 .0 .0 .0 .0 1354.2 2708.3 -6313.1 -1354.2 32261.7 .0 .0 .0 9360.0 .0 9360.0 .0 .0 .0 .0 .0 .0 -6240.0 77480.0 -6240.0 .0 .0 .0 -9360.0 .0 9360.0 .0 -65000.0 .0 .0 .0 .0
.0 .0 .0 .0 -6240.0 84240.0 .0 .0 .0 .0 -9360.0 9360.0 .0 .0 -78000.0 .0 .0 .0 .0 .0
.0 .0 .0 -156000.0 312451.4 -156000.0 .0 1354.2 .0 .0 .0 .0 .0 -451.4 .0 .0 1354.2 .0 .0 .0 .0 .0 -156000.0 156780.0 .0 .0 2340.0 .0 .0 .0 .0 .0 -780.0 .0 .0 9360.0 .0 -9360.0 .0 1354.2 .0 9360.0 42856.7 9360.0 .0 .0 .0 .0 -1354.2 .0 .0 2708.3 .0 .0 9360.0 9360.0 .0 .0 2340.0 .0 9360.0 28080.0 .0 .0 .0 .0 .0 -2340.0 .0 .0 .0 -65000.0 .0 .0 .0 .0 .0 .0 .0 -6240.0 77480.0 -6240.0 .0 .0 .0 -9360.0 .0 9360.0 .0 .0 -78000.0 .0 .0 .0 .0 .0 .0 .0 -6240.0 84240.0 .0 .0 .0 .0 -9360.0 .0 .0 .0 .0 -451.4 .0 .0 -1354.2 .0 .0 .0 .0 -156000.0 312451.4 -156000.0 .0 -1354.2 .0 .0 .0 .0 .0 .0 -780.0 .0 .0 -2340.0 .0 .0 .0 .0 -156000.0 156780.0 .0 .0 .0 .0 .0 .0 1354.2 .0 .0 2708.3 .0 9360.0 .0 -9360.0 .0 -1354.2 .0 9360.0
42856.7 9360.0 .0 .0 .0 .0 .0 2340.0 .0 .0 4680.0 .0 9360.0 9360.0 .0 .0 -2340.0 .0 9360.0
19
.0 2340.0 4680.0 9360.0 -2340.0 28080.0
结点荷载数 单元荷载数
2 0
结点荷载:
4.000 .100 6.000 -3.000
荷载向量:
.000 .000 .100 .000 .000 .000 .000 .000 -3.000 .000 .000 .000 .000 .000 .000 .000 .000 .000
结构位移列阵:
.00004778 .00000394 -.00001362 .00005279 -.00000472 -.00004415
.00024659 .00000641 -.00014793 .00024633 -.00000720 -.00005586
.00044198 .00000716 -.00000955 .00043871 -.00000794 -.00004518
单元位移列阵:
M UA VA ZTA UB VB ZTB 1 .0000000 .0000000 .0000000 .0000478 .0000039 -.0000136 2 .0000478 .0000039 -.0000136 .0002466 .0000064 -.0001479 3 .0002466 .0000064 -.0001479 .0004420 .0000072 -.0000095 4 .0000000 .0000000 .0000000 .0000528 -.0000047 -.0000442 5 .0000528 -.0000047 -.0000442 .0002463 -.0000072 -.0000559 6 .0002463 -.0000072 -.0000559 .0004387 -.0000079 -.0000452 7 .0004420 .0000072 -.0000095 .0004387 -.0000079 -.0000452 8 .0002466 .0000064 -.0001479 .0002463 -.0000072 -.0000559 9 .0000478 .0000039 -.0000136 .0000528 -.0000047 -.0000442
四、实验日期：二〇一二年六月二十二日星期五
20
实验五: 计算杆端力和支座反力的程序设计
一、图示结构竖杆的EA=468000KN,EI=14040KNm2,其余杆件
EA=390000KN,EI=8125KNm2。求结构的总载荷向量。 76.008
56.006
36.004
二、输入数据文件
9 8 2 2 2 0 3
1 3 3 5 5 7 2 4 4 6 6 8 7 8 5 6 3 4 1 1 1 1 2 1 1 1
1 1 1 2 1 1 1 2 2
468000 390000
14040 8125
0 6 0 6 0 6 0 6
0 0 4 4 7 7 10 10
4 .1,6,-3,8.1,6,-3,9.1,6,-3
三、计算结果：
21
平面杆系分析
单元数NE= 9
结点数NJ= 8
特殊结点数NJT= 2
支座结点数NJZ= 2
截面特性分组数NAI= 2
单元两端结点号:
单元 起点 终点 单元 起点 终点 单元 起点 终点 单元 起点 终点 1 1 3 2 3 5 3 5 7 4 2 4 5 4 6 6 6 8 7 7 8 8 5 6 9 3 4
单元截面特性:
1 1 1 2 1 1 1 2 2
特殊结点信息(结点号,结点各自由度的约束特征数):
结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 结点号 特征数 特征数 1 1 1 1 2 1 1 1
截面特性EA:
468000.000 390000.000
截面特性EI:
14040.000 8125.000
结点坐标:
结点 X坐标 Y坐标 结点 X坐标 Y坐标 结点 X坐标 Y坐标 1 .000 .000 2 6.000 .000 3 .000 4.000
4 6.000 4.000 5 .000 7.000 6 6.000 7.000
7 .000 10.000 8 6.000 10.000
未知量编号:
0 0 0 0 0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
单元定位向量MW(6)
单元编号: 1 MW = 0 0 0 1 2 3
单元编号: 2 MW = 1 2 3 7 8 9
单元编号: 3 MW = 7 8 9 13 14 15
单元编号: 4 MW = 0 0 0 4 5 6
特征数 22
单元编号: 5 MW = 4 5 6 10 11 12
单元编号: 单元编号: 单元编号: 单元编号:
6 7 8 9
MW = 10 11 12 16 17 18 MW = 13 14 15 16 17 18 MW = 7 8 9 10 11 12 MW = 1 2 3 4 5 6
单元刚度矩阵：单元号 1
2632.50 .00 -5265.00 -2632.50 .00 -5265.00 .00 117000.00 .00 .00 -117000.00 .00 -5265.00 .00 14040.00 -2632.50 .00 5265.00 .00 -117000.00 .00 -5265.00 .00 7020.00
单元刚度矩阵：单元号 2
6240.00 .00 -9360.00 .00 156000.00 .00 -9360.00 .00 18720.00 -6240.00 .00 9360.00 .00 -156000.00 .00 -9360.00 .00 9360.00
单元刚度矩阵：单元号 3
6240.00 .00 -9360.00 .00 156000.00 .00 -9360.00 .00 18720.00 -6240.00 .00 9360.00 .00 -156000.00 .00 -9360.00 .00 9360.00
单元刚度矩阵：单元号 4
1523.44 .00 -3046.88 .00 97500.00 .00 -3046.88 .00 8125.00 -1523.44 .00 3046.88 .00 -97500.00 .00 -3046.88 .00 4062.50
单元刚度矩阵：单元号 5
6240.00 .00 -9360.00 .00 156000.00 .00 -9360.00 .00 18720.00 -6240.00 .00 9360.00 .00 -156000.00 .00
5265.00 2632.50 .00 5265.00 -6240.00 .00 9360.00 6240.00 .00 9360.00 -6240.00 .00 9360.00 6240.00 .00 9360.00 -1523.44 .00 3046.88 1523.44 .00 3046.88 -6240.00 .00 9360.00 6240.00 .00 .00 .00 117000.00 .00 .00 -156000.00 .00 .00 156000.00 .00 .00 -156000.00 .00 .00 156000.00 .00 .00 -97500.00 .00 .00 97500.00 .00 .00 -156000.00 .00 .00 156000.00 7020.00 5265.00 .00 14040.00 -9360.00 .00 9360.00 9360.00 .00 18720.00 -9360.00 .00 9360.00 9360.00 .00 18720.00 -3046.88 .00 4062.50 3046.88 .00 8125.00 -9360.00 .00 9360.00 9360.00 .00
23
-9360.00 .00 9360.00 9360.00 .00 18720.00
单元刚度矩阵：单元号 6
6240.00 .00 -9360.00 -6240.00 .00 -9360.00 .00 156000.00 .00 .00 -156000.00 .00 -9360.00 .00 18720.00 9360.00 .00 9360.00 -6240.00 .00 9360.00 6240.00 .00 9360.00 .00 -156000.00 .00 .00 156000.00 .00 -9360.00 .00 9360.00 9360.00 .00 18720.00
单元刚度矩阵：单元号 7
78000.00 .00 .00 .00 780.00 2340.00 .00 2340.00 9360.00 -78000.00 .00 .00 .00 -780.00 -2340.00 .00 2340.00 4680.00
单元刚度矩阵：单元号 8
65000.00 .00 .00 .00 451.39 1354.17 .00 1354.17 5416.67 -65000.00 .00 .00 .00 -451.39 -1354.17 .00 1354.17 2708.33
单元刚度矩阵：单元号 9
65000.00 .00 .00 .00 451.39 1354.17 .00 1354.17 5416.67 -65000.00 .00 .00 .00 -451.39 -1354.17 .00 1354.17 2708.33
结构刚度矩阵：
73872.5 .0 -4095.0 -6240.0 .0 -9360.0 .0 .0 .0 .0 273451.4 1354.2 .0 -156000.0 .0 .0 .0 .0 -4095.0 1354.2 38176.7 9360.0 .0 9360.0 .0 .0 .0
-78000.00 .00 .00 78000.00 .00 .00 -65000.00 .00 .00 65000.00 .00 .00 -65000.00 .00 .00 65000.00 .00 .00 -65000.0 .0 .0 .0 .0 .0 .0 .0 .0 .00 -780.00 -2340.00 .00 780.00 -2340.00 .00 -451.39 -1354.17 .00 451.39 -1354.17 .00 -451.39 -1354.17 .00 451.39 -1354.17 .0 .0 .0 -451.4 .0 .0 -1354.2 .0 .0 .00 2340.00 4680.00 .00 -2340.00 9360.00 .00 1354.17 2708.33 .00 -1354.17 5416.67 .00 1354.17 2708.33 .00 -1354.17 5416.67 .0 .0 .0 1354.2 .0 .0 2708.3 .0 .0
24
-65000.0 .0 .0 72763.4 .0 -6313.1
.0 .0 .0 -6240.0 .0 -9360.0 .0 .0 .0 .0 .0 .0 .0 -451.4 -1354.2 .0 253951.4 -1354.2 .0 .0 .0 .0 -156000.0 .0
.0 .0 .0 .0 .0 .0 .0 1354.2 2708.3 -6313.1 -1354.2 32261.7 .0 .0 .0 9360.0 .0 9360.0 .0 .0 .0 .0 .0 .0 -6240.0 77480.0 -6240.0 .0 .0 .0 -9360.0 .0 9360.0 .0 -65000.0 .0 .0 .0 .0
.0 .0 .0 .0 -6240.0 84240.0 .0 .0 .0 .0 -9360.0 9360.0 .0 .0 -78000.0 .0 .0 .0 .0 .0
.0 .0 .0 -156000.0 312451.4 -156000.0 .0 1354.2 .0 .0 .0 .0 .0 -451.4 .0 .0 1354.2 .0 .0 .0 .0 .0 -156000.0 156780.0 .0 .0 2340.0 .0 .0 .0 .0 .0 -780.0 .0 .0 9360.0 .0 -9360.0 .0 1354.2 .0 9360.0 42856.7 9360.0 .0 .0 .0 .0 -1354.2 .0 .0 2708.3 .0 .0 9360.0 9360.0 .0 .0 2340.0 .0 9360.0 28080.0 .0 .0 .0 .0 .0 -2340.0 .0 .0 .0 -65000.0 .0 .0 .0 .0 .0 .0 .0 -6240.0 77480.0 -6240.0 .0 .0 .0 -9360.0 .0 9360.0 .0 .0 -78000.0 .0 .0 .0 .0 .0 .0 .0 -6240.0 84240.0 .0 .0 .0 .0 -9360.0 .0 .0 .0 .0 -451.4 .0 .0 -1354.2 .0 .0 .0 .0 -156000.0 312451.4 -156000.0 .0 -1354.2 .0 .0 .0 .0 .0 .0 -780.0 .0 .0 -2340.0 .0 .0 .0 .0 -156000.0 156780.0 .0 .0 .0 .0 .0 .0 1354.2 .0 .0 2708.3 .0 9360.0 .0 -9360.0 .0 -1354.2 .0 9360.0
42856.7 9360.0 .0 .0 .0 .0 .0 2340.0 .0 .0 4680.0 .0 9360.0 9360.0 .0 .0 -2340.0 .0 9360.0
25
.0 2340.0 4680.0 9360.0 -2340.0 28080.0
结点荷载数 单元荷载数
2 0
结点荷载:
4.000 .100 6.000 -3.000
荷载向量:
.000 .000 .100 .000 .000 .000
.000 .000 -3.000 .000 .000 .000
.000 .000 .000 .000 .000 .000
结构位移列阵:
.00004778 .00000394 -.00001362 .00005279 -.00000472 -.00004415
.00024659 .00000641 -.00014793 .00024633 -.00000720 -.00005586
.00044198 .00000716 -.00000955 .00043871 -.00000794 -.00004518
单元位移列阵:
M UA VA ZTA UB VB ZTB 1 .0000000 .0000000 .0000000 .0000478 .0000039 -.0000136
2 .0000478 .0000039 -.0000136 .0002466 .0000064 -.0001479
3 .0002466 .0000064 -.0001479 .0004420 .0000072 -.0000095
4 .0000000 .0000000 .0000000 .0000528 -.0000047 -.0000442
5 .0000528 -.0000047 -.0000442 .0002463 -.0000072 -.0000559
6 .0002463 -.0000072 -.0000559 .0004387 -.0000079 -.0000452
7 .0004420 .0000072 -.0000095 .0004387 -.0000079 -.0000452
8 .0002466 .0000064 -.0001479 .0002463 -.0000072 -.0000559
9 .0000478 .0000039 -.0000136 .0000528 -.0000047 -.0000442
单元内力
单元 NA QA MA NB QB MB 1 -.460 .054 .156 .460 -.054 .060 2 -.386 -.272 .221 .386 .272 -1.036 3 -.116 -.255 -1.030 .116 .255 .265
4 .460 -.054 -.019 -.460 .054 -.198
5 .386 .272 .462 -.386 -.272 .353 6 .116 .255 .332 -.116 -.255 .432 7 .255 -.116 -.265 -.255 .116 -.432
8 .017 -.270 -.934 -.017 .270 -.685
9 -.326 -.074 -.182 .326 .074 -.264
26
支座反力
支点 RA RB MB 1 -.054 -.460 .156 2 .054 .460 -.019
四、实验日期：二〇一二年六月二十二日星期五
27
实验一金属试件打磨以及贴片和焊线实验目的1掌握金属试件打磨方法2初步掌握电阻应变片的粘贴技术3初步掌握焊接和检查实验设备和器材1电…
一实验一材料循环应力应变曲线的探究实验目的通过对材料进行循环加载并利用所获取数据绘出材料的循环应力应变曲线认识与分析材料的循环特性…
纯弯曲梁正应力电测实验一实验目的1电测法测定纯弯曲梁正应力分布规律2验证纯弯曲梁正应力计算公式3计算弹性模量E二实验装置与仪器1纯…
材料力学实验报告评分标准拉伸实验报告一、实验目的（1分）1.测定低碳钢的强度指标（s、b）和塑性指标（、）。2.测定铸铁的强度极限…
实验一拉伸实验一实验目的1测定低碳钢Q235的屈服点s强度极限b延伸率断面收缩率2测定铸铁的强度极限b3观察低碳钢拉伸过程中的各种…
一弹塑性力学发展史一弹性力学的发展近代弹性力学可认为始于柯西CauchyAL在1882年引进应变与应力的概念建立了平衡微分方程边界…
姓名图尔荪江斯拉吉学号1083310402理论力学材料力学以及结构力学的关系摘要通过学习一个学期的结构力学课程对结构力学分析及计算…
结构动力学读书报告转眼间这个学期就快结束了庆幸的是跟着张老师学习完本门课程我受益颇多在此就主要根据上课所用的河海大学出版社出版的结…
高等结构动力学实验报告学院航天学院专业固体力学姓名沈延臣学号12S11800320xx年11月11日悬臂桁架结构模态分析实验报告一…
高等岩石力学读书报告学院专业姓名学号国土资源工程学院地质工程120xx20xx71高等岩石力学读书报告岩石力学是研究岩石在外界因素…
结构力学认知实习报告实习时间20xx年5月25日学院土木工程与力学学院专业名称土木工程学号U20xx15430姓名HbFeng指导…