结构力学实验报告
(土木工程学院)
班级
姓名
学号
实验一 平面体系的机动分析
一、输入结构体系
二、结构图形生成
三、几何组成分析
1、自动求解
2、智能求解
实验二 静定平面刚架弯矩图绘制
一、输入结构体系
二、结构图形生成
三、内力计算
1、杆端内力值
2、约束反力值
3、弯矩图
实验三 静定平面桁架内力计算
一、输入结构体系
二、结构图形生成
三、内力计算
1、几何组成分析
2、内力求解
3、内力计算结果
实验四 超静定结构内力计算
一、输入结构体系
二、结构图形生成
三、内力计算
1、几何组成分析
2、位移计算
3、内力计算
实验五 影响线绘制
一、输入结构体系
二、结构图形生成
三、内力计算
1、参数设定
2、影响线求解
第二篇:计算结构力学实验报告
《计算结构力学》
计算报告
姓名: 周剑敏
学号: 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