计算结构力学

时间:2024.5.8

《计算结构力学》

实验报告

姓名: 鲁洋

学号: 08201010318

指导: 苏少卿 晏燕

安徽建筑工业学院土木工程学院

实验一: 形成单元定位向量的程序设计

一、求图示结构各单元的单元定位向量。

二、输入数据文件

8,10,10

9,2,9,4,10,6,10,8,1,3,3,5,5,7,9,10

1,1,1,0,2,1,1,10001,3,10001,0,0,4,3,3,10001,5,10001,0,0,6, 5,5,10001,7,10001,1,0,8,7,1,10001,9,0,0,10001,10,0,0,10001

三、计算结果:

平面杆系分析

单元数NE= 8 结点数NJ= 10 特殊结点数NJT= 10

单元两端结点号:

单元 起点 终点 单元 起点 终点 单元 起点 1 9 2 2 9 4 3 10 5 1 3 6 3 5 7 5

特殊结点信息(结点号,结点各自由度的约束特征数):

结点号 特征数 特征数 特征数 结点号 特征数 特征数 1 1 1 0 2 1 1 4 3 3 10001 5 10001 0 7 10001 1 0 8 7 1 10 0 0 10001

终点 单元 起点 终点 6 4 10 8 7 8 9 10 特征数 结点号 特征数 特征数 10001 3 10001 0 0 6 5 5 10001 9 0 0 特征数 0 10001 10001

未知量编号:

0 0 1 0 0 -1 -1 2 3 -1 2 -1 -1 4 5 -1 4 -1 -1 0 6 -1 0 -1 7 8 -1 9 10 -1

单元定位向量MW(6)

单元编号: 1 MW = 7 8 -1 0 0 -1

单元编号: 2 MW = 7 8 -1 -1 2 -1

单元编号: 3 MW = 9 10 -1 -1 4 -1

单元编号: 4 MW = 9 10 -1 -1 0 -1

单元编号: 5 MW = 0 0 1 -1 2 3

单元编号: 6 MW = -1 2 3 -1 4 5

单元编号: 7 MW = -1 4 5 -1 0 6

单元编号: 8 MW = 7 8 -1 9 10 -1

四、实验分析:

1.承受荷载的横杆为梁式杆件,其余皆为二力杆;

2.单元定位向量与结点坐标和杆件长度、截面性质、材料都无关。

五、实验日期:二〇一一年六月二十一日星期二

实验二: 形成单元刚度矩阵和结构刚度矩阵的程序设计

一、图示为梁桁,桁式杆的EA=390000,横梁的EI=14040,且不考虑轴向变形。

求单元的单元刚度矩阵和结构刚度矩阵。

二、输入数据文件

8,10,10,2

9,2,9,4,10,6,10,8,1,3,3,5,5,7,9,10

1,1,1,0,2,1,1,10001,3,10001,0,0,4,3,3,10001,5,10001,0,0

6,5,5,10001,7,10001,1,0,8,7,1,10001,9,0,0,10001,10,0,0,10001

1,1,1,1,2,2,2,1

390000,0

0,14040

0,0,2,2,6,6,8,8,2,6

2,2,2,2,2,2,2,2,0,0计算结果:

平面杆系分析

单元数NE= 8

结点数NJ= 10

特殊结点数NJT= 10

截面特性分组数NAI= 2

单元两端结点号:

单元 起点 终点 单元 起点 终点 单元 起点 终点

1 9 2 2 9 4 3 10 6 4 10 8 单元 起点终点

5 1 3 6 3 5 7 5 7 8 9 10

单元截面特性:

1 1 1 1 2 2 2 1

特殊结点信息(结点号,结点各自由度的约束特征数):

结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 1 1 1 0 2 1 1 10001 3 10001 0 0 4 3 3 10001 5 10001 0 0 6 5 5 10001 7 10001 1 0 8 7 1 10001 9 0 0 10001 10 0 0 10001

截面特性EA:

390000.000 .000

截面特性EI:

.000 14040.000

结点坐标:

结点 X坐标 Y坐标 结点 X坐标 Y坐标 结点 X坐标 Y坐标

1 .000 2.000 2 .000 2.000 3 2.000 2.000

4 2.000 2.000 5 6.000 2.000 6 6.000 2.000

7 8.000 2.000 8 8.000 2.000 9 2.000 .000

10 6.000 .000

未知量编号:

0 0 1 0 0 -1 -1 2 3 -1 2 -1 -1 4 5 -1

4 -1 -1 0 6 -1 0 -1 7 8 -1 9 10 -1

单元定位向量MW(6)

单元编号: 1 MW = 7 8 -1 0 0 -1

单元编号: 2 MW = 7 8 -1 -1 2 -1

单元编号: 3 MW = 9 10 -1 -1 4 -1

单元编号: 4 MW = 9 10 -1 -1 0 -1

单元编号: 5 MW = 0 0 1 -1 2 3

单元编号: 6 MW = -1 2 3 -1 4 5

单元编号: 7 MW = -1 4 5 -1 0 6

单元编号: 8 MW = 7 8 -1 9 10 -1

单元刚度矩阵:单元号 1

68942.91 -68942.91 .00 -68942.91 68942.91 .00

-68942.91 68942.91 .00 68942.91 -68942.91 .00

.00 .00 .00 .00 .00 .00

-68942.91 68942.91 .00 68942.91 -68942.91 .00

68942.91 -68942.91 .00 -68942.91 68942.91 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 2

.00 .00 .00 .00 .00 .00 .00 195000.00 .00 .00 -195000.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 -195000.00 .00 .00 195000.00 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 3

.00 .00 .00 .00 .00 .00 .00 195000.00 .00 .00 -195000.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 -195000.00 .00 .00 195000.00 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 4

68942.91 68942.91 .00 -68942.91 -68942.91 .00 68942.91 68942.91 .00 -68942.91 -68942.91 .00 .00 .00 .00 .00 .00 .00 -68942.91 -68942.91 .00 68942.91 68942.91 .00 -68942.91 -68942.91 .00 68942.91 68942.91 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 5

.00 .00 .00 .00 .00 .00 .00 21060.00 21060.00 .00 -21060.00 21060.00 .00 21060.00 28080.00 .00 -21060.00 14040.00 .00 .00 .00 .00 .00 .00 .00 -21060.00 -21060.00 .00 21060.00 -21060.00 .00 21060.00 14040.00 .00 -21060.00 28080.00

单元刚度矩阵:单元号 6

.00 .00 .00 .00 .00 .00 .00 2632.50 5265.00 .00 -2632.50 5265.00 .00 5265.00 14040.00 .00 -5265.00 7020.00 .00 .00 .00 .00 .00 .00 .00 -2632.50 -5265.00 .00 2632.50 -5265.00 .00 5265.00 7020.00 .00 -5265.00 14040.00

单元刚度矩阵:单元号 7

.00 .00 .00 .00 .00 .00 .00 21060.00 21060.00 .00 -21060.00 21060.00 .00 21060.00 28080.00 .00 -21060.00 14040.00 .00 .00 .00 .00 .00 .00 .00 -21060.00 -21060.00 .00 21060.00 -21060.00 .00 21060.00 14040.00 .00 -21060.00 28080.00

单元刚度矩阵:单元号 8

97500.00 .00 .00 -97500.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 -97500.00 .00 .00 97500.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00

结构刚度矩阵:

28080.0 -21060.0 14040.0 .0 .0 .0 .0 .0 .0 .0 -21060.0 218692.5 -15795.0 -2632.5 5265.0 .0 .0 -195000.0 .0 .0 14040.0 -15795.0 42120.0 -5265.0 7020.0 .0 .0 .0 .0 .0 .0 -2632.5 -5265.0 218692.5 15795.0 21060.0 .0 .0 .0 -195000.0 .0 5265.0 7020.0 15795.0 42120.0 14040.0 .0 .0 .0 .0 .0 .0 .0 21060.0 14040.0 28080.0 .0 .0 .0 .0 .0 .0 .0 .0 .0 .0 166442.9 -68942.9 -97500.0 .0 .0 -195000.0 .0 .0 .0 .0 -68942.9 263942.9 .0 .0 .0 .0 .0 .0 .0 .0 -97500.0 .0 166442.9 68942.9 .0 .0 .0 -195000.0 .0 .0 .0 .0 68942.9 263942.9

实验分析:

1.承受荷载的横杆为梁式杆件,其余皆为二力杆;

2.单元定位向量与结点坐标和杆件长度、截面性质、材料都无关。

三、实验日期:二〇一一年六月二十一日星期二

实验三: 形成载荷列阵的程序设计

一、图示为梁桁,桁式杆的EA=390000,横梁的EI=14040,且不考虑轴向变形。

求结构的总载荷向量。

二、输入数据文件

8,10,10,2,0,3

9,2,9,4,10,6,10,8,1,3,3,5,5,7,9,10

1,1,1,0,2,1,1,10001,3,10001,0,0,4,3,3,10001,5,10001,0,0

6,5,5,10001,7,10001,1,0,8,7,1,10001,9,0,0,10001,10,0,0,10001

1,1,1,1,2,2,2,1

390000,0

0,14040

0,0,2,2,6,6,8,8,2,6

2,2,2,2,2,2,2,2,0,0

5.1,2,-20,6.1,4,-20 7.1,2,-20

三、计算结果:

平面杆系分析

单元数NE= 8

结点数NJ= 10

特殊结点数NJT= 10

支座结点数NJZ= 4

截面特性分组数NAI= 2

单元两端结点号:

单元 起点 终点 单元 起点 终点 单元 起点 终点单元 起点终点

1 9 2 2 9 4 3 10 6 4 10 8

5 1 3 6 3 5 7 5 7 8 9 10

单元截面特性:

1 1 1 1 2 2 2 1

特殊结点信息(结点号,结点各自由度的约束特征数):

结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 1 1 1 0 2 1 1 10001 3 10001 0 0 4 3 3 10001 5 10001 0 0 6 5 5 10001 7 10001 1 0 8 7 1 10001 9 0 0 10001 10 0 0 10001

截面特性EA:

390000.000 .000

截面特性EI:

.000 14040.000

结点坐标:

结点 X坐标 Y坐标 结点 X坐标 Y坐标 结点 X坐标 Y坐标

1 .000 2.000 2 .000 2.000 3 2.000 2.000

4 2.000 2.000 5 6.000 2.000 6 6.000 2.000

7 8.000 2.000 8 8.000 2.000 9 2.000 .000

10 6.000 .000

未知量编号:

0 0 1 0 0 -1 -1 2 3 -1 2 -1 -1 4 5 -1

4 -1 -1 0 6 -1 0 -1 7 8 -1 9 10 -1

单元定位向量MW(6)

单元编号: 1 MW = 7 8 -1 0 0 -1

单元编号: 2 MW = 7 8 -1 -1 2 -1

单元编号: 3 MW = 9 10 -1 -1 4 -1

单元编号: 4 MW = 9 10 -1 -1 0 -1

单元编号: 5 MW = 0 0 1 -1 2 3

单元编号: 6 MW = -1 2 3 -1 4 5

单元编号: 7 MW = -1 4 5 -1 0 6

单元编号: 8 MW = 7 8 -1 9 10 -1

单元刚度矩阵:单元号 1

68942.91 -68942.91 .00 -68942.91 68942.91 .00

-68942.91 68942.91 .00 68942.91 -68942.91 .00

.00 .00 .00 .00 .00 .00

-68942.91 68942.91 .00 68942.91 -68942.91 .00 68942.91 -68942.91 .00 -68942.91 68942.91 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 2

.00 .00 .00 .00 .00 .00 .00 195000.00 .00 .00 -195000.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 -195000.00 .00 .00 195000.00 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 3

.00 .00 .00 .00 .00 .00 .00 195000.00 .00 .00 -195000.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 -195000.00 .00 .00 195000.00 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 4

68942.91 68942.91 .00 -68942.91 -68942.91 .00 68942.91 68942.91 .00 -68942.91 -68942.91 .00 .00 .00 .00 .00 .00 .00 -68942.91 -68942.91 .00 68942.91 68942.91 .00 -68942.91 -68942.91 .00 68942.91 68942.91 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 5

.00 .00 .00 .00 .00 .00 .00 21060.00 21060.00 .00 -21060.00 21060.00 .00 21060.00 28080.00 .00 -21060.00 14040.00 .00 .00 .00 .00 .00 .00 .00 -21060.00 -21060.00 .00 21060.00 -21060.00 .00 21060.00 14040.00 .00 -21060.00 28080.00

单元刚度矩阵:单元号 6

.00 .00 .00 .00 .00 .00 .00 2632.50 5265.00 .00 -2632.50 5265.00 .00 5265.00 14040.00 .00 -5265.00 7020.00 .00 .00 .00 .00 .00 .00 .00 -2632.50 -5265.00 .00 2632.50 -5265.00 .00 5265.00 7020.00 .00 -5265.00 14040.00

单元刚度矩阵:单元号 7

.00 .00 .00 .00 .00 .00 .00 21060.00 21060.00 .00 -21060.00 21060.00 .00 21060.00 28080.00 .00 -21060.00 14040.00 .00 .00 .00 .00 .00 .00 .00 -21060.00 -21060.00 .00 21060.00 -21060.00 .00 21060.00 14040.00 .00 -21060.00 28080.00

单元刚度矩阵:单元号 8

97500.00 .00 .00 -97500.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 -97500.00 .00 .00 97500.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00

结构刚度矩阵:

28080.0 -21060.0 14040.0 .0 .0 .0 .0 .0 .0 .0 -21060.0 218692.5 -15795.0 -2632.5 5265.0 .0 .0 -195000.0 .0 .0 14040.0 -15795.0 42120.0 -5265.0 7020.0 .0 .0 .0 .0 .0 .0 -2632.5 -5265.0 218692.5 15795.0 21060.0 .0 .0 .0 -195000.0 .0 5265.0 7020.0 15795.0 42120.0 14040.0 .0 .0 .0 .0 .0 .0 .0 21060.0 14040.0 28080.0 .0 .0 .0 .0 .0 .0 .0 .0 .0 .0 166442.9 -68942.9 -97500.0 .0 .0 -195000.0 .0 .0 .0 .0 -68942.9 263942.9 .0 .0 .0 .0 .0 .0 .0 .0 -97500.0 .0 166442.9 68942.9 .0 .0 .0 -195000.0 .0 .0 .0 .0 68942.9 263942.9

结点荷载数 单元荷载数

0 3

单元荷载:

5.100 2.000 -20.000 6.100 4.000 -20.000

7.100 2.000 -20.000

荷载向量:

-6.667 -60.000 -20.000 -60.000 20.000 6.667 .000 .000 .000 .000

四、实验分析:

1.承受荷载的横杆为梁式杆件,其余皆为二力杆;

2.单元定位向量与结点坐标和杆件长度、截面性质、材料都无关。

五、实验日期:二〇一一年六月二十一日星期二

实验四: 高斯消元法的程序设计

一、 图示为梁桁,桁式杆的EA=390000,横梁的EI=14040,且不考虑轴向变形。

求结构的结点位移。

二、 输入数据文件

8,10,10,2,0,3

9,2,9,4,10,6,10,8,1,3,3,5,5,7,9,10

1,1,1,0,2,1,1,10001,3,10001,0,0,4,3,3,10001,5,10001,0,0

6,5,5,10001,7,10001,1,0,8,7,1,10001,9,0,0,10001,10,0,0,10001

1,1,1,1,2,2,2,1

390000,0

0,14040

0,0,2,2,6,6,8,8,2,6

2,2,2,2,2,2,2,2,0,0

5.1,2,-20,6.1,4,-20 7.1,2,-20

三、计算结果:

平面杆系分析

单元数NE= 8

结点数NJ= 10

特殊结点数NJT= 10

支座结点数NJZ= 4

截面特性分组数NAI= 2

单元两端结点号:

单元 起点 终点 单元 起点 终点 单元 起点 终点单元 起点终点

1 9 2 2 9 4 3 10 6 4 10 8

5 1 3 6 3 5 7 5 7 8 9 10

单元截面特性:

1 1 1 1 2 2 2 1

特殊结点信息(结点号,结点各自由度的约束特征数):

结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 1 1 1 0 2 1 1 10001 3 10001 0 0 4 3 3 10001 5 10001 0 0 6 5 5 10001 7 10001 1 0 8 7 1 10001 9 0 0 10001 10 0 0 10001

截面特性EA:

390000.000 .000

截面特性EI:

.000 14040.000

结点坐标:

结点 X坐标 Y坐标 结点 X坐标 Y坐标 结点 X坐标 Y坐标

1 .000 2.000 2 .000 2.000 3 2.000 2.000

4 2.000 2.000 5 6.000 2.000 6 6.000 2.000

7 8.000 2.000 8 8.000 2.000 9 2.000 .000

10 6.000 .000

未知量编号:

0 0 1 0 0 -1 -1 2 3 -1 2 -1 -1 4 5 -1

4 -1 -1 0 6 -1 0 -1 7 8 -1 9 10 -1

单元定位向量MW(6)

单元编号: 1 MW = 7 8 -1 0 0 -1

单元编号: 2 MW = 7 8 -1 -1 2 -1

单元编号: 3 MW = 9 10 -1 -1 4 -1

单元编号: 4 MW = 9 10 -1 -1 0 -1

单元编号: 5 MW = 0 0 1 -1 2 3

单元编号: 6 MW = -1 2 3 -1 4 5

单元编号: 7 MW = -1 4 5 -1 0 6

单元编号: 8 MW = 7 8 -1 9 10 -1

单元刚度矩阵:单元号 1

68942.91 -68942.91 .00 -68942.91 68942.91 .00

-68942.91 68942.91 .00 68942.91 -68942.91 .00

.00 .00 .00 .00 .00 .00

-68942.91 68942.91 .00 68942.91 -68942.91 .00 68942.91 -68942.91 .00 -68942.91 68942.91 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 2

.00 .00 .00 .00 .00 .00 .00 195000.00 .00 .00 -195000.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 -195000.00 .00 .00 195000.00 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 3

.00 .00 .00 .00 .00 .00 .00 195000.00 .00 .00 -195000.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 -195000.00 .00 .00 195000.00 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 4

68942.91 68942.91 .00 -68942.91 -68942.91 .00 68942.91 68942.91 .00 -68942.91 -68942.91 .00 .00 .00 .00 .00 .00 .00 -68942.91 -68942.91 .00 68942.91 68942.91 .00 -68942.91 -68942.91 .00 68942.91 68942.91 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 5

.00 .00 .00 .00 .00 .00 .00 21060.00 21060.00 .00 -21060.00 21060.00 .00 21060.00 28080.00 .00 -21060.00 14040.00 .00 .00 .00 .00 .00 .00 .00 -21060.00 -21060.00 .00 21060.00 -21060.00 .00 21060.00 14040.00 .00 -21060.00 28080.00

单元刚度矩阵:单元号 6

.00 .00 .00 .00 .00 .00 .00 2632.50 5265.00 .00 -2632.50 5265.00 .00 5265.00 14040.00 .00 -5265.00 7020.00 .00 .00 .00 .00 .00 .00 .00 -2632.50 -5265.00 .00 2632.50 -5265.00 .00 5265.00 7020.00 .00 -5265.00 14040.00

单元刚度矩阵:单元号 7

.00 .00 .00 .00 .00 .00 .00 21060.00 21060.00 .00 -21060.00 21060.00 .00 21060.00 28080.00 .00 -21060.00 14040.00 .00 .00 .00 .00 .00 .00 .00 -21060.00 -21060.00 .00 21060.00 -21060.00 .00 21060.00 14040.00 .00 -21060.00 28080.00

单元刚度矩阵:单元号 8

97500.00 .00 .00 -97500.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 -97500.00 .00 .00 97500.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00

结构刚度矩阵:

28080.0 -21060.0 14040.0 .0 .0 .0 .0 .0 .0 .0 -21060.0 218692.5 -15795.0 -2632.5 5265.0 .0 .0 -195000.0 .0 .0 14040.0 -15795.0 42120.0 -5265.0 7020.0 .0 .0 .0 .0 .0 .0 -2632.5 -5265.0 218692.5 15795.0 21060.0 .0 .0 .0 -195000.0 .0 5265.0 7020.0 15795.0 42120.0 14040.0 .0 .0 .0 .0 .0 .0 .0 21060.0 14040.0 28080.0 .0 .0 .0 .0 .0 .0 .0 .0 .0 .0 166442.9 -68942.9 -97500.0 .0 .0 -195000.0 .0 .0 .0 .0 -68942.9 263942.9 .0 .0 .0 .0 .0 .0 .0 .0 -97500.0 .0 166442.9 68942.9 .0 .0 .0 -195000.0 .0 .0 .0 .0 68942.9 263942.9

结点荷载数 单元荷载数

0 3

单元荷载:

5.100 2.000 -20.000 6.100 4.000 -20.000

7.100 2.000 -20.000

荷载向量:

-6.667 -60.000 -20.000 -60.000 20.000 6.667 .000 .000 .000 .000

结构位移列阵:

-.00090170 -.00170855 -.00123425 -.00170855 .00123425 .00090170 -.00035385 -.00135470 .00035385 -.00135470

单元位移列阵:

M UA VA ZTA UB VB ZTB 1 -.0003539 -.0013547 .0000000 .0000000 .0000000 .0000000 2 -.0003539 -.0013547 .0000000 .0000000 -.0017085 .0000000 3 .0003539 -.0013547 .0000000 .0000000 -.0017085 .0000000 4 .0003539 -.0013547 .0000000 .0000000 .0000000 .0000000 5 .0000000 .0000000 -.0009017 .0000000 -.0017085 -.0012342 6 .0000000 -.0017085 -.0012342 .0000000 -.0017085 .0012342 7 .0000000 -.0017085 .0012342 .0000000 .0000000 .0009017 8 -.0003539 -.0013547 .0000000 .0003539 -.0013547 .0000000 四、实验分析:

1.承受荷载的横杆为梁式杆件,其余皆为二力杆;

2.单元定位向量与结点坐标和杆件长度、截面性质、材料都无关。

五、验日期:二〇一一年六月二十一日星期二

实验五: 计算杆端力和支座反力的程序设计

一、图示为梁桁,桁式杆的EA=390000,横梁的EI=14040,且不考虑轴向变形。

求结构的杆件杆端力和支座反力。

二、输入数据文件

8,10,10,2,0,3

9,2,9,4,10,6,10,8,1,3,3,5,5,7,9,10

1,1,1,0,2,1,1,10001,3,10001,0,0,4,3,3,10001,5,10001,0,0

6,5,5,10001,7,10001,1,0,8,7,1,10001,9,0,0,10001,10,0,0,10001

1,1,1,1,2,2,2,1

390000,0

0,14040

0,0,2,2,6,6,8,8,2,6

2,2,2,2,2,2,2,2,0,0

5.1,2,-20,6.1,4,-20 7.1,2,-20

三、计算结果:

平面杆系分析

单元数NE= 8

结点数NJ= 10

特殊结点数NJT= 10

支座结点数NJZ= 4

截面特性分组数NAI= 2

单元两端结点号:

单元 起点 终点 单元 起点 终点 单元 起点 终点 单元 起点 终点

1 9 2 2 9 4 3 10 6 4 10 8

5 1 3 6 3 5 7 5 7 8 9 10

单元截面特性:

1 1 1 1 2 2 2 1

特殊结点信息(结点号,结点各自由度的约束特征数):

结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 结点号 特征数 特征数

1 1 1 0 2 1 1 10001 3 10001 0 0 4 3 3 10001 5 10001 0 0 6 5 5 10001 7 10001 1 0 8 7 1 10001 9 0 0 10001 10 0 0 10001

截面特性EA:

390000.000 .000

截面特性EI:

.000 14040.000

结点坐标:

结点 X坐标 Y坐标 结点 X坐标 Y坐标 结点 X坐标 Y坐标

1 .000 2.000 2 .000 2.000 3 2.000 2.000

4 2.000 2.000 5 6.000 2.000 6 6.000 2.000

7 8.000 2.000 8 8.000 2.000 9 2.000 .000

10 6.000 .000

未知量编号:

0 0 1 0 0 -1 -1 2 3 -1 2 -1 -1 4 5 -1

4 -1 -1 0 6 -1 0 -1 7 8 -1 9 10 -1

单元定位向量MW(6)

单元编号: 1 MW = 7 8 -1 0 0 -1

单元编号: 2 MW = 7 8 -1 -1 2 -1

单元编号: 3 MW = 9 10 -1 -1 4 -1

单元编号: 4 MW = 9 10 -1 -1 0 -1

单元编号: 5 MW = 0 0 1 -1 2 3

单元编号: 6 MW = -1 2 3 -1 4 5

单元编号: 7 MW = -1 4 5 -1 0 6

单元编号: 8 MW = 7 8 -1 9 10 -1

单元刚度矩阵:单元号 1

68942.91 -68942.91 .00 -68942.91 68942.91 .00 特征数

-68942.91 68942.91 .00 68942.91 -68942.91 .00 .00 .00 .00 .00 .00 .00 -68942.91 68942.91 .00 68942.91 -68942.91 .00 68942.91 -68942.91 .00 -68942.91 68942.91 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 2

.00 .00 .00 .00 .00 .00 .00 195000.00 .00 .00 -195000.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 -195000.00 .00 .00 195000.00 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 3

.00 .00 .00 .00 .00 .00 .00 195000.00 .00 .00 -195000.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 -195000.00 .00 .00 195000.00 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 4

68942.91 68942.91 .00 -68942.91 -68942.91 .00 68942.91 68942.91 .00 -68942.91 -68942.91 .00 .00 .00 .00 .00 .00 .00 -68942.91 -68942.91 .00 68942.91 68942.91 .00 -68942.91 -68942.91 .00 68942.91 68942.91 .00 .00 .00 .00 .00 .00 .00

单元刚度矩阵:单元号 5

.00 .00 .00 .00 .00 .00 .00 21060.00 21060.00 .00 -21060.00 21060.00 .00 21060.00 28080.00 .00 -21060.00 14040.00 .00 .00 .00 .00 .00 .00 .00 -21060.00 -21060.00 .00 21060.00 -21060.00 .00 21060.00 14040.00 .00 -21060.00 28080.00

单元刚度矩阵:单元号 6

.00 .00 .00 .00 .00 .00 .00 2632.50 5265.00 .00 -2632.50 5265.00 .00 5265.00 14040.00 .00 -5265.00 7020.00 .00 .00 .00 .00 .00 .00 .00 -2632.50 -5265.00 .00 2632.50 -5265.00

.00 5265.00 7020.00 .00 -5265.00 14040.00

单元刚度矩阵:单元号 7

.00 .00 .00 .00 .00 .00 .00 21060.00 21060.00 .00 -21060.00 21060.00 .00 21060.00 28080.00 .00 -21060.00 14040.00 .00 .00 .00 .00 .00 .00 .00 -21060.00 -21060.00 .00 21060.00 -21060.00 .00 21060.00 14040.00 .00 -21060.00 28080.00

单元刚度矩阵:单元号 8

97500.00 .00 .00 -97500.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 -97500.00 .00 .00 97500.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00

结构刚度矩阵:

28080.0 -21060.0 14040.0 .0 .0 .0 .0 .0 .0 .0 -21060.0 218692.5 -15795.0 -2632.5 5265.0 .0 .0 -195000.0 .0 .0 14040.0 -15795.0 42120.0 -5265.0 7020.0 .0 .0 .0 .0 .0 .0 -2632.5 -5265.0 218692.5 15795.0 21060.0 .0 .0 .0 -195000.0 .0 5265.0 7020.0 15795.0 42120.0 14040.0 .0 .0 .0 .0 .0 .0 .0 21060.0 14040.0 28080.0 .0 .0 .0 .0 .0 .0 .0 .0 .0 .0 166442.9 -68942.9 -97500.0 .0 .0 -195000.0 .0 .0 .0 .0 -68942.9 263942.9 .0 .0 .0 .0 .0 .0 .0 .0 -97500.0 .0 166442.9 68942.9 .0 .0 .0 -195000.0 .0 .0 .0 .0 68942.9 263942.9

结点荷载数 单元荷载数

0 3

单元荷载:

5.100 2.000 -20.000 6.100 4.000 -20.000

7.100 2.000 -20.000

荷载向量:

-6.667 -60.000 -20.000 -60.000 20.000 6.667 .000 .000 .000 .000

结构位移列阵:

-.00090170 -.00170855 -.00123425 -.00170855 .00123425 .00090170 -.00035385 -.00135470 .00035385 -.00135470

单元位移列阵:

M UA VA ZTA UB VB ZTB 1 -.0003539 -.0013547 .0000000 .0000000 .0000000 .0000000 2 -.0003539 -.0013547 .0000000 .0000000 -.0017085 .0000000 3 .0003539 -.0013547 .0000000 .0000000 -.0017085 .0000000 4 .0003539 -.0013547 .0000000 .0000000 .0000000 .0000000 5 .0000000 .0000000 -.0009017 .0000000 -.0017085 -.0012342 6 .0000000 -.0017085 -.0012342 .0000000 -.0017085 .0012342 7 .0000000 -.0017085 .0012342 .0000000 .0000000 .0009017 8 -.0003539 -.0013547 .0000000 .0003539 -.0013547 .0000000

单元内力

单元 NA QA MA NB QB MB 1 -97.582 .000 .000 97.582 .000 .000 2 69.001 .000 .000 -69.001 .000 .000 3 69.001 .000 .000 -69.001 .000 .000 4 -97.582 .000 .000 97.582 .000 .000 5 .000 10.999 .000 .000 29.001 -18.002 6 .000 40.000 18.002 .000 40.000 -18.002 7 .000 29.001 18.002 .000 10.999 .000 8 -69.001 .000 .000 69.001 .000 .000

支座反力

支点 RA RB MB

1 .000 10.999 .000

2 -69.001 69.001 .000

7 .000 10.999 .000

8 .000 69.001 .000

四、实验分析:

1.承受荷载的横杆为梁式杆件,其余皆为二力杆;

2.单元定位向量与结点坐标和杆件长度、截面性质、材料都无关。

五、验日期:二〇一一年六月二十一日星期二

更多相关推荐:
结构力学感想

感悟结构力学这学期开设土木工程专业基础课结构力学,给我第一印象是:难并且复杂,但是实用。结构力学(StructuralMechanics)是固体力学的一个分支,它主要研究工程结构受力和传力的规律,以及如何进行结…

结构力学心得体会

我对结构力学的感悟结构力学是土木工程专业重要的专业基础课,其前导课程有高等数学、理论力学、材料力学等课程,后续课程又有混凝土结构、钢结构、隧道工程、桥梁工程等大量的专业课。学生学习起来感觉内容繁杂头绪众多,花费…

结构力学心得体会

本学期结构力学的课程已经接近尾声。主要是三部分内容,即渐近法、矩阵位移法和平面刚架静力分析的程序设计。通过为期八周的理论课学习和六次的上机课程设计,我收获颇丰。而对结构力学半年的学习,也让我对这门学科有了很大的…

结构力学基础总结

结构静力分析的学习要求(复习指导)结构(几何)组成分析Geometricstabilityanalysisofstructures除理解和记住各名词含义外,要熟练掌握利用基本组成规律进行体系分析。总的来说分析方…

结构力学上机心得

结构力学学习心得结构力学的学习马上就要结束了本学期学的主要是渐进法矩阵位移法和平面刚架静力分析程序设计相比上学期的画内力图和计算这学期貌似任务比较轻需要动手的不多但理解上难度较大上学期学的是比较基础的理论分析只...

结构力学第六次学习与讨论心得

结构力学第六次学习与讨论心得这周主要学习了静定结构影响线的作法与应用一开始接触影响线感觉与内力图有点类似但是二者完全不同影响线反映的是单位移动荷载作用下某量值随荷载位置变化情况而内力图则是在一固定位置荷载作用下...

结构力学学习心得体会-浅谈对称性在结构力学中的应用

结构力学学习心得体会浅谈对称性在结构力学中的应用摘要在工程实际问题中有很多结构都具有对称性我们对这些结构进行受力分析的时候常常将结构简化为杆系模型而结构力学研究的就是结构的杆系模型因此对称性在结构力学中有着广泛...

哈工大007材料学院【822材料结构与力学性能】19xx~20xx全套优级考研试题辅导笔记资料汇总

哈尔滨工业大学007材料科学工程学院080502材料学学科19xx20xx全套优秀考研笔记试题资料库20xx年9月20xx年9月更新注意全套资料均为本团队哈工大硕博导师团队6年精心积攒及整理并进行了严格精华筛选...

华科结构力学实习报告

结构力学认知实习报告实习时间20xx年5月25日学院土木工程与力学学院专业名称土木工程学号U20xx15430姓名HbFeng指导老师戴萍地点华中科技大学光谷同济医院一引言1通过结构力学的认识实习使同学们认识各...

《结构力学》自学指导书

结构力学自学指导书一课程编码及适用专业课程编码025221211学期总学时88学期面授学时24学期自学学时64适用专业适用于土木函授本科专业二课程性质结构力学是土木函授本科专业的重点技术基础课是一门既重理论又重...

华科结构力学实习报告

结构力学实习报告实习时间20xx年5月25日指导老师龙晓鸿实习地点华中科技大学光谷同济医院华中科技大学光谷同济医院实习报告一实习目的1通过切身的实践经历亲临建筑施工现场增强对建筑施工整个过程的感性认识为今后的理...

同济大牛结构力学与材料力学考研437分心得(geilivable experience)

下面转载同济大学20xx年初试第一名的考研心得因为私人信息不方便透露谢谢总分437分英语71政治77数学150材料力学与结构力学139考研准备从07年七月22日开始复习到08年一月19日20日参加考研初试三月9...

结构力学学习总结(27篇)