计算结构力学实验报告(3)

-89.61 134.41 194.15 89.61 -134.41 388.29

单元刚度矩阵：单元号 6

103.12 0.00 -206.25 -103.12 0.00 -206.25

0.00 55.00 0.00 0.00 -55.00 0.00

-206.25 0.00 550.00 206.25 0.00 275.00

-103.12 0.00 206.25 103.12 0.00 206.25

0.00 -55.00 0.00 0.00 55.00 0.00

-206.25 0.00 275.00 206.25 0.00 550.00

单元刚度矩阵：单元号 7

825.00 0.00 -825.00 -825.00 0.00 -825.00

0.00 110.00 0.00 0.00 -110.00 0.00

-825.00 0.00 1100.00 825.00 0.00 550.00

-825.00 0.00 825.00 825.00 0.00 825.00

0.00 -110.00 0.00 0.00 110.00 0.00

-825.00 0.00 550.00 825.00 0.00 1100.00

单元刚度矩阵：单元号 8

103.12 0.00 -206.25 -103.12 0.00 -206.25

0.00 55.00 0.00 0.00 -55.00 0.00

-206.25 0.00 550.00 206.25 0.00 275.00

-103.12 0.00 206.25 103.12 0.00 206.25

0.00 -55.00 0.00 0.00 55.00 0.00

-206.25 0.00 275.00 206.25 0.00 550.00

结构刚度矩阵：

1000.7 44.1 549.8 -119.3 -66.3 -185.6

-56.4 22.2 -89.6 0.0 0.0 0.0

0.0 0.0 0.0 44.1 304.2 -51.2

-66.3 -119.3 -185.6 22.2 -74.8 134.4

0.0 0.0 0.0 0.0 0.0 0.0

549.8 -51.2 1983.3 185.6 185.6 247.5

89.6 -134.4 194.1 0.0 0.0 0.0

0.0 0.0 0.0 -119.3 -66.3 185.6

500.6 -226.2 33.4 0.0 0.0 0.0

-278.2 292.4 -358.4 0.0 0.0 0.0

-66.3 -119.3 185.6 -226.2 696.2 723.2

0.0 0.0 0.0 292.4 -521.9 537.6

0.0 0.0 0.0 -185.6 -185.6 247.5

33.4 723.2 1821.6 0.0 0.0 0.0

358.4 -537.6 388.3 0.0 0.0 0.0

-56.4 22.2 89.6 0.0 0.0 0.0

955.0 667.3 -446.6 0.0 0.0 0.0

-795.5 -689.4 -742.5 22.2 -74.8 -134.4

0.0 0.0 0.0 667.3 925.3 -876.9

0.0 0.0 0.0 -689.4 -795.5 -742.5

-89.6 134.4 194.1 0.0 0.0 0.0

-446.6 -876.9 1928.2 0.0 0.0 0.0

742.5 742.5 495.0 0.0 0.0 0.0

-278.2 292.4 358.4 0.0 0.0 0.0

338.2 -292.4 358.4 -60.0 0.0 0.0

0.0 0.0 0.0 292.4 -521.9 -537.6

0.0 0.0 0.0 -292.4 790.7 -537.6

0.0 -268.8 0.0 0.0 0.0 0.0

-358.4 537.6 388.3 0.0 0.0 0.0

358.4 -537.6 776.6 0.0 0.0 0.0

0.0 0.0 0.0 0.0 0.0 0.0

-795.5 -689.4 742.5 -60.0 0.0 0.0

855.5 689.4 742.5 0.0 0.0 0.0

0.0 0.0 0.0 -689.4 -795.5 742.5

0.0 -268.8 0.0 689.4 1064.3 742.5

0.0 0.0 0.0 0.0 0.0 0.0

-742.5 -742.5 495.0 0.0 0.0 0.0

742.5 742.5 989.9

结点荷载数 单元荷载数

2 1

结点荷载:

5.100 14.142 5.200 -14.142

单元荷载:

8.100 4.000 15.000

荷载向量:

0.000 0.000 0.000 14.142 -14.142 0.000

-30.000 0.000 -20.000 0.000 0.000 0.000

0.000 0.000 0.000

四、 实验分析：

1．单元①是二力杆，是轴力杆单元；

2．20KN的集中力也可认为是单元④的单元荷载。

五、 实验日期：二〇〇八年三月十四日星期五

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

一、 图示结构竖杆的EA=220KN,EI=550KNm2,其余杆件

EA=150KN,EI=350KNm2。求结构的结点位移。

8 9

10

5 7

6

4

1 2 3

二、 输入数据文件

8 10 5 2 2 1

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

1 1 1 1 2 1 1 1 3 1 1 1 8 7 7 10001 9 10 10 10001

1 1 1 1 1 2 2 2

150 220

350 550

0 2 5 2 0 5 1.5 1.5 4 4

0 0 0 2 4 4 5 5 5 5

5.1 14.14213562 5.2 -14.14213562

8.1 4 15

三、 计算结果：

平面杆系分析

单元数NE= 8

结点数NJ= 10

特殊结点数NJT= 5

支座结点数NJZ= 3

截面特性分组数NAI= 2

单元两端结点号:

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

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

5 4 6 6 1 5 7 2 4 8 3 6

单元截面特性:

1 1 1 1 1 2 2 2

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

结点号 …… 此处隐藏：1910字，全部文档内容请下载后查看。喜欢就下载吧 ……